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

[Александр Мурый]

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

Думаю, потому, что разработчики библиотеки PROJ.4 (которых сейчас чуть более, чем 1 человек) почему-то не следят за российскими ГОСТами. Приходится нам самим разбираться...

[Дмитрий Барышников]

PROJ.4 позволяет указать свои параметры перехода. Хоть 3 хоть 7. Реализовано в OGRSpatialReference.
http://www.gdal.org/ogr/classOGRSpatial ... 6416097156

[dab]

Хорошо, спасибо за пояснения.
Тогда возвращаюсь к своему старому вопросу viewtopic.php?f=17&t=10665&start=30#p71624

Если для МСК на базе СК-63 в MapInfo так заданы параметры

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

"МСК-56 Зона 3", 8, 1001, 7, 58.03333333333, 0, 1, 3300000, -5309414.70

Какой из вариантов описания МСК для QGIS будет верным?
1.

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

+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.70 +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0.35,0.82,-0.12 +units=m +no_defs

2.

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

+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.70 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,-0.35,-0.79,-0.22 +units=m +no_defs

3.

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

+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.70 +ellps=krass +units=m +no_defs

[Boris]

Dы будете смеяться, но правильный вам создал QGIS, как многократно сказано, параметры мапинфо - есть свои собственные "параметры мапинфо". Это не хорошо, не плохо, это просто фича такая.
Вы опишите задачу, потому как если все ваши вычисления происходят в пределах СК на одном датуме (Красовского) то как он позиционирован относительно датума WGS84, должно иметь мало влияния для PROJ4 и тех продуктов, что используют его библиотеки.
Исключение может быть только в том случае, если вы изучили коды PROJ4 и уверены, что трансформации между СК на одном датуме, он все равно всегда проделывает через отнесение точек на датум WGS84.
Иное дело, если вы собираетесь сопрягать данные, где есть и другие датумы. Тогда это вообще не тема беседы - надо выяснить локальные параметры, или выбрать любые понравившиеся. Благо для шейпа это можно проделать "легким движением руки" - исходные координаты меняться ведь не должны?

[dab]

Смеяться не буду, так как недостаточно пока знаю кугис и не могу отличить глюк от фичи

параметры мапинфо - есть свои собственные

Конечно да, но разные ГИС, каждая со своими собственными параметрами, пересчитывать координаты из географической долготы-широты в конкретную СК должны же одинаково?
Проект перевода кадастровых инженеров с мапинфо на кугис.
Задача - перевести данные из формате мапинфо (кварталы, участки, топооснова...) в "справочные" слои кугис, для использования в режиме "посмотреть-копировать-вставить" и создать новые слои кугис для работы с новыми объектами кадастрового учёта. Плюс значительное количество топосновы в формате Панорамы с типом карты "Местная система на базе СК-63".
Значения координат не должны меняться, дабы избежать расхождения по объектам, прошедшим кадастровый учёт, особенно на переходный период "mi => qgis"
Исходные данные все на одном датуме - "Пулково 1942" 8, 9, 10, 11 зоны плюс несколько региональных МСК на базе СК-63.
Новые слои планирую делать в соответствующих району работ региональных МСК на базе СК-63.
На сегодняшний момент времени, думаю что для данной задачи корректнее будет, например, для старых данных
СК-42 зона 10 использовать

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

+proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +units=m +no_defs

для новых, например МСК-56 зона 3

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

+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.70 +ellps=krass +units=m +no_defs

По такой схеме пересчёт точки с географическими координатами
N51.5 E58.5
в СК-42 зона 10 выдаёт координаты
5 708 879.67 10 604 159.84
во всех трёх ГИС - Mapinfo, QGis, Панорама 11

если вы изучили коды PROJ4 и уверены, что трансформации между СК на одном датуме, он все равно всегда проделывает через отнесение точек на датум WGS84

К сожалению, мне не хватит знаний и опыта для изучения кода PROJ4, но именно такие сомнения у меня присутствуют.

[Донецков]

во всех трёх ГИС - Mapinfo, QGis, Панорама 11

- так по идее же для MI не должно, т.к. у нее свое "понимание" Пулково-42?

[Voltron]

Смеяться не буду, так как недостаточно пока знаю кугис и не могу отличить глюк от фичи

Я вам даже больше скажу. Этот вопрос не имеет никакого отношения к QGIS. Никакого.

[dab]

Я вам даже больше скажу. Этот вопрос не имеет никакого отношения к QGIS. Никакого.

Не спорю, но если позволите, последний мой вопрос в этой теме по работе в QGIS c СК-42.
Есть точка с координатами N51.5 E58.5
В СК-42 ей соответствуют координаты X=5708879.67 Y=10604159.84
В QGis пересчёт из EPSG:4130 в EPSG:28410 даёт значения X=5708850.31 Y=10604246.05
Это вопрос к проекту QGis или к проекту PROJ4?

[Voltron]

Разумеется не к QGIS

[Александр Мурый]

Есть точка с координатами N51.5 E58.5
В СК-42 ей соответствуют координаты X=5708879.67 Y=10604159.84
В QGis пересчёт из EPSG:4130 в EPSG:28410 даёт значения X=5708850.31 Y=10604246.05

Вы, наверное, имели в виду не EPSG:4130, а EPSG:4326 (WGS84 Lat/Long)?

Проверка при помощи <cs2cs> показала, что PROJ.4 пересчитывает корректно:

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

echo "58.5 51.5" | cs2cs -f "%.2f" +proj=latlong +datum=WGS84 +to +proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +units=m +no_defs 
10604159.84	5708879.67 0.00

[dab]

Да

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

cs2cs -f "%.2f" +proj=latlong +datum=WGS84 +to +proj=tmerc +lat_0=0 +lon_0=57 +k=1 +x_0=10500000 +y_0=0 +ellps=krass +units=m +no_defs
58.5 51.5
10604159.84     5708879.67 0.00

cs2cs -f "%.2f" +proj=latlong +datum=WGS84 +to +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
58.5 51.5
10604246.05     5708850.31 21.65

Так почему в QGis по умолчанию для СК-42 зона 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

[Voltron]

Так почему в QGis по умолчанию для СК-42 зона 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

Очевидно, потому что точно такое же описание используется в proj.4

[Boris]

я как человек ленивый, совершенно запутался в обилии разных примеров с кодом.
автор вопроса декларирует работу в СК, построенных на эллипсоиде Красовского-1940, даже СК-95 не была упомянута ни разу. А правды ищет в переводе координат, которые с очевидность (для меня) заданы на "глобусе" WGS84.
Тут одно из двух, либо перестаньте использовать координаты в WGS84, либо покажите нам, что исходная точка на эллипсоиде Красовского в долготе/широте, переходит в разные точки на спроецированной СК Гаусс-Крюгер, хоть МСК, хоть СК-42, в случае, если вы задаете один и тот же способ отнесения эллипсоида Красовского к WGS84 для исходных и для конечных данных.
Все примеры подтверждают очевидное - проекция WGS84 на Красовского-1940 очень сильно зависит от параметров отнесения последнего к первому. И что? Всем это очевидно. А вопрос то в чем?

[dab]

...А вопрос то в чем?

Вопрос был такой

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

Ответа, по сути, было два

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

и

Dы будете смеяться, но правильный вам создал QGIS...

Спасибо. Пошел думать

[dab]

Возник ещё вопрос в эту же тему.
Создаю в мапинфо пользовательскую проекцию для МСК, согласно параметрам перехода СК-42 => WGS-84 взятым отсюда viewtopic.php?f=34&t=10551
То есть строка в mif файле выглядит так

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 58.03333333333, 0, 1, 3300000, -5309414.7

Открываю mif в кугис и получаю пользовательскую СК

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

+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.7 +ellps=krass +towgs84=23.57,-140.95,-79.8,-0,0.35,0.79,-0.22 +units=m +no_defs

То есть два параметра были отрицательные, а стали положительные:
... -0.35, -0.79 ... => ... 0.35,0.79 ...
Почему происходит смена знака?