Логика открытия слоёв в 28410 в QGIS и в MapInfo

Вопросы по свободной ГИС QGIS. Сообщения об ошибках, предложения по улучшению, локализация.
Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Логика открытия слоёв в 28410 в QGIS и в MapInfo

Сообщение dab » 02 окт 2012, 08:46

Донецков писал(а):...В QGIS если открыть таблицу MI, которой задана проекция, отсутствующая в системах координат QGIS, то создается "безымянная" СК с исходными параметрами, которую только остается переименовать и сохранить в пользовательских СК...
В продолжении этого разговора - возник вопрос при работе с СК на базе Пулково 1942.
Например, последовательность действий(MapInfo 10.5, QGIS 1.8.0-Lisboa)
1. в мапинфо создаю новую таблицу t1.tab в EPSG:28410

Код: Выделить всё

GK Зона 10 (Пулково 1942)\p28410", 8, 1001, 7, 57, 0, 1, 10500000, 0
2. создаю полигон с координатами

Код: Выделить всё

Region  1
  4
10580000 5780000
10600000 5720000
10540000 5740000
10580000 5780000
    Pen (1,2,0) 
    Brush (1,0,16777215)
    Center 10570000 5750000
3. открываю t1.tab в кугис, и, да - создаётся пользовательская СК

Код: Выделить всё

USER:100000 - * Создана система координат (+proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +units=m +no_defs)
хотя в кугис EPSG:28410

Код: Выделить всё

Pulkovo 1942 / Gauss-Kruger zone 10 EPSG:28410 (+proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs)
4. соответственно, далее при работе с указанным полигоном в кугис:
если пользоваться USER:100000

Код: Выделить всё

"features": [
{ "type": "Feature", "id": 1, "properties": { "id": 0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 10580000.0, 5780000.0 ], [ 10600000.0, 5720000.0 ], [ 10540000.0, 5740000.0 ], [ 10580000.0, 5780000.0 ] ] ] } }
]
если пользоваться Pulkovo 1942 / Gauss-Kruger zone 10 EPSG:28410

Код: Выделить всё

"features": [
{ "type": "Feature", "id": 1, "properties": { "id": 0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 10580003.795934485, 5779989.9949698839 ], [ 10600003.48970212, 5719989.8854421275 ], [ 10540003.585843131, 5739990.1952549759 ], [ 10580003.795934485, 5779989.9949698839 ] ] ] } }
]
Координаты различаются на метры :(
Стоит задача взять накопленные кадастровые геоданные из мапинфо и в дальнейшем работать с ними в кугис.
Как корректно провести этот процесс, с точки зрения кугис, что бы не тащить тяжелое наследство СК "USER:100000 ..." и в то же время не получить расхождения в координатах по участкам, прошедшим кадастровый учёт?
Вложения
t1.zip
(2.06 КБ) 318 скачиваний
Последний раз редактировалось dab 11 окт 2012, 07:40, всего редактировалось 2 раза.

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение Максим Дубинин » 02 окт 2012, 09:46

приложите созданный вами tab посмотреть
пристегивайтесь, турбулентность прямо по курсу

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение dab » 02 окт 2012, 12:00

Максим Дубинин писал(а):приложите созданный вами tab посмотреть
Приложил к своему предыдущему посту.

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Параметры региональных МСК для MapInfo

Сообщение Александр Мурый » 02 окт 2012, 12:43

dab писал(а): Стоит задача взять накопленные кадастровые геоданные из мапинфо и в дальнейшем работать с ними в кугис.
Как корректно провести этот процесс, с точки зрения кугис, что бы не тащить тяжелое наследство СК "USER:100000 ..." и в то же время не получить расхождения в координатах по участкам, прошедшим кадастровый учёт?
Ну да, у MapInfo и PROJ.4 разные параметры "towgs84". Создайте в QGIS пользовательскую проекцию с параметрами:

Код: Выделить всё

+proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +units=m +no_defs)
и пользуйтесь ей для всех данных. В чём проблема-то?
Редактор материалов, модератор форума

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение dab » 02 окт 2012, 13:23

amuriy писал(а):...В чём проблема-то?
Теоретически - разные параметры означают разные математические модели реального мира. Так какая модель более точная и "правильная"?
Практически - ок, сделаю пользовательскую, перенесу накопленные в мапинфо данные, отказавшись от стандартной проекции кугис, и всё будет хорошо, пока не появятся, например, сторонние геоданные со стандарной проекцией кугис. Перепроецировать? Куда? В любом случае получу расхождения с координатами по документам.
Просто не представляю как это можно будет "состыковать" с необходимой для кадастрового инженера точностью.

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Параметры региональных МСК для MapInfo

Сообщение Александр Мурый » 02 окт 2012, 13:30

dab писал(а): Так какая модель более точная и "правильная"?
Зависит от задачи, разве нет? Ещё ж есть ГОСТы 2001/2008 — какие параметры перехода для СК "самые идеальные"?
dab писал(а): Просто не представляю как это можно будет "состыковать" с необходимой для кадастрового инженера точностью.
А как это обычно делается (не в QGIS)?
Редактор материалов, модератор форума

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение dab » 02 окт 2012, 13:53

amuriy писал(а):А как это обычно делается (не в QGIS)?
Мне понятна процедура открытия, например, шейп файла в мапинфо, когда через диалоговое окно система спрашивает о проекции открываемых данных. Я указываю проекцию из готового списка. Если открываю геоданные с нестандартной проекцией - тогда я должен предварительно её параметры прописать в гис.
По такому принципу работает импорт-экспорт через обменные форматы в MapInfo, RSI ENVI, Панорама.
Кугис ни о чём не спрашивая делает некую пользовательскую проекцию, а когда указываешь стандартную - получаешь изменение координат :shock:
Опыта работы с кугис у меня конечно мало, но хочется понять принципы. :)

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

Re: Параметры региональных МСК для MapInfo

Сообщение Донецков » 02 окт 2012, 14:57

Столкнулся с этим же... Для себя решил, в MI в prj добавляю новую СК на базе ГОСТ от 2008 г., с помощью утилиты ReCalc (или через Mif) меняю СК-42 по умолчанию в СК-42 по ГОСТ 2008 (значения координат сохранятся без изменений), в QGIS добавляю пользовательскую СК также с параметрами ГОСТ, открываю tab/mif в QGIS с указанием созданной СК, получаю координаты идентичные с MI...
Последний раз редактировалось Донецков 13 мар 2013, 14:26, всего редактировалось 1 раз.

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Параметры региональных МСК для MapInfo

Сообщение Александр Мурый » 02 окт 2012, 14:58

dab писал(а): Мне понятна процедура открытия, например, шейп файла в мапинфо, когда через диалоговое окно система спрашивает о проекции открываемых данных. Я указываю проекцию из готового списка. Если открываю геоданные с нестандартной проекцией - тогда я должен предварительно её параметры прописать в гис.
По такому принципу работает импорт-экспорт через обменные форматы в MapInfo, RSI ENVI, Панорама.
Очень рад за MapInfo, RSI ENVI, Панораму и прочих, но не будем сводить дискуссию к вопросам религии :)
dab писал(а): Кугис ни о чём не спрашивая делает некую пользовательскую проекцию, а когда указываешь стандартную - получаешь изменение координат :shock:
ИМХО, логично: если в самих данных прописана СК/проекция, почему бы ГИС не открыть их в этой СК/проекции, даже если она не прописана в самой ГИС? От перемены мест слагаемых ...
Редактор материалов, модератор форума

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение dab » 02 окт 2012, 15:14

Донецков писал(а):... в MI в prj добавляю новую СК на базе ГОСТ от 2008 г...
Через 6 параметров аффинного преобразования?

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Параметры региональных МСК для MapInfo

Сообщение dab » 03 окт 2012, 07:14

amuriy писал(а):Очень рад за MapInfo, RSI ENVI, Панораму и прочих, но не будем сводить дискуссию к вопросам религии :)
Нет-нет, я не религиозный человек :) Вы же спросили как это делается не в QGis - я и ответил из собственного опыта.
amuriy писал(а):ИМХО, логично: если в самих данных прописана СК/проекция, почему бы ГИС не открыть их в этой СК/проекции, даже если она не прописана в самой ГИС? От перемены мест слагаемых ...
Ну да, "перемена слагаемых" - проверка на логику в виде простого эксперимента по схеме "туда-обратно"
1. В кугис создаю шейп в СК "EPSG:28410 - Pulkovo 1942 / Gauss-Kruger zone 10" под названием 28410.shp и создаю точку с координатами (10540000, 5740000)
2. из кугис экспортирую в mif под названием mi.mif

Код: Выделить всё

Version 300
Charset "Neutral"
Delimiter ","
CoordSys Earth Projection 8, 1001, "m", 57, 0, 1, 10500000, 0
Columns 1
  id Char(3)
Data

Point 10540000 5740000
    Symbol (35,0,12)
3. открываю в кугис mi.mif (!) и автоматически получаю для этого слоя "USER:100000 - * Создана система координат ..." с координатами в EPSG:28410 (10540003.585843131, 5739990.1952549759)
Где логика?

ИМХО, если в мапинфо параметры EPSG:28410 не соответствуют ГОСТ, то пусть эта проблема и остаётся в рамках мапинфо. При импорте файлов мапинфо в кугис не нужно создавать пользовательскую СК, а считать описание

Код: Выделить всё

"GK Зона 10 (Пулково 1942)\p28410", 8, 1001, 7, 57, 0, 1, 10500000, 0
тождественным

Код: Выделить всё

Pulkovo 1942 / Gauss-Kruger zone 10 EPSG:28410 (+proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs)
Вложения
qgis-mi-qgis.zip
(4.57 КБ) 315 скачиваний

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

Re: Параметры региональных МСК для MapInfo

Сообщение Донецков » 03 окт 2012, 08:13

Через 6 параметров аффинного преобразования?
, нет в prj добавляются параметры СК на основании ГОСТа без использования аффинных преобразований, а по типу:
"GK зона 6 (Пулково 1942–WGS)", 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 33, 0, 1, 6500000, 0
в частности параметры СК с использованием аффинных преобразований, которыее можно создать в MI можно ли аналогично воспроизвести QGIS?

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 792
Ваше звание: званий не имею
Откуда: Москва

Re: Параметры региональных МСК для MapInfo

Сообщение Александр Мурый » 03 окт 2012, 12:07

dab, может быть, это поведение связано с ошибками QGIS или OGR в работе с данными в форматах MapInfo. Надо поискать в баг-трекере. Если это ошибки, то это сюда.
Редактор материалов, модератор форума

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Логика открытия слоёв в 28410 в QGIS и в MapInfo

Сообщение dab » 03 окт 2012, 15:13

amuriy, понятно, спасибо за ссылки.
Свою практическую задачу по переносу данных из MI в QGis решил.
Воспользовался утилитой MI EasyLoader и перенёс данные из MapInfo в PostGIS.
В QGis открыл и сразу получил слой в EPSG:28410 с сохранением значений координат и площадей :D

Аватара пользователя
dab
Гуру
Сообщения: 671
Зарегистрирован: 16 дек 2011, 20:02
Репутация: 170
Ваше звание: Гуру
Откуда: Москва
Контактная информация:

Re: Логика открытия слоёв в 28410 в QGIS и в MapInfo

Сообщение dab » 04 окт 2012, 09:50

А почему параметры QGis для EPSG:28410 соответствуют ГОСТ 51794-2001, а не ГОСТ 51794-2008 ?
ГОСТ 51794-2001

Код: Выделить всё

23.92, -141.27, -80.9, 0, 0.35, 0.82, -0.12
ГОСТ 51794-2008

Код: Выделить всё

23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22
Вопрос возник после топика viewtopic.php?f=34&t=10551

Ответить

Вернуться в «QGIS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 9 гостей