Страница 1 из 3

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

Добавлено: 02 окт 2012, 08:46
dab
Донецков писал(а):...В 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 ..." и в то же время не получить расхождения в координатах по участкам, прошедшим кадастровый учёт?

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

Добавлено: 02 окт 2012, 09:46
Максим Дубинин
приложите созданный вами tab посмотреть

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

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

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)
и пользуйтесь ей для всех данных. В чём проблема-то?

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

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

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

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

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

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

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

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

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

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

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

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

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

Добавлено: 03 окт 2012, 07:14
dab
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)

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?

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

Добавлено: 03 окт 2012, 12:07
Александр Мурый
dab, может быть, это поведение связано с ошибками QGIS или OGR в работе с данными в форматах MapInfo. Надо поискать в баг-трекере. Если это ошибки, то это сюда.

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

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

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

Добавлено: 04 окт 2012, 09:50
dab
А почему параметры 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