Преобразование из WGS84 в Местную Московскую СК

[gobsy]

Добрый день, господа форумчане.
Помогите разобраться с преобразованием из WGS84 в Московскую систему координат.
исходные данные см. прикрепленные архивы.
Задача: преобразовать текущий шейп из wgs84 в местную московскую СК.
Использую для преобразования ogr2ogr.
последовательность шагов:
1) произвожу 7-ми параметрическую трансформацию из WGS84 в Пулково "Pulkovo_1942_To_WGS_1984_12" параметры взяты из geographic_transformations.pdf

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

ogr2ogr -s_srs "+proj=latlong +ellps=WGS84 +datum=WGS84" -t_srs "+proj=tmerc +ellps=krass +towgs84=21.587,-97.541,-60.925,-1.01378,-0.58117,-0.2348,-4.6121" -f "ESRI Shapefile" 42_7_param.shp reserve_point.shp

2) собственно преобразую из полученного shape-file в Московскую.

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

ogr2ogr  -t_srs "+proj=tmerc +ellps=krass +x_0=-19.1 +y_0=4.9 +lon_0=37.4998 +k_0=1.0 +lat_0=55.667 +units=m" -f "ESRI Shapefile" 42_CS_param.shp 42_7_param.shp

насколько я понимаю в результате преобразования должен создаться набор файлов в нужной мне проекции однако 42_CS_param.prj сгенерированный огром не совпадает с исходным Мoscow_cs.prj. и выглядит следующим образом:

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

PROJCS["Transverse_Mercator",GEOGCS["GCS_Krassovsky, 1942",DATUM["D_unknown",SPHEROID["krass",6378245,298.3]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",55.667],PARAMETER["central_meridian",37.4998],PARAMETER["scale_factor",1],PARAMETER["false_easting",-19.1],PARAMETER["false_northing",4.9],UNIT["Meter",1]] 

разница в том что тут DATUM = "D_unknown" и UNIT["Degree",0.017453292519943295] а должен быть
DATUM = "D_Pulkovo_1942" и UNIT["Degree",0.0174532925199433].

результат преобразования сравниваю c вариантом сделанным в ArcGis ошибка примерно 50 метров.

Вопрос: Что я делаю не так ?? и Если есть возможность помогите с преобразованием.

[Boris]

До того как углубляться в вычисления и параметры, необходимо уточнить о какой именно Московской СК идет речь. Поскольку уже много раз обсуждалось, что СК города Москвы, т.е. наборы геоданных, которые ведет Мосгеотрест, ведутся в местной московской СК, основанной на эллипсоиде Бесселя-1841 года, с центром где-то в районе МИГАИКа.

[gobsy]

Все чем я могу ответить на этот вопрос, это привести содержание того проекционного файла который у меня есть.

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

PROJCS["MoscowCS",
GEOGCS["GCS_Pulkovo_1942",DATUM["D_Pulkovo_1942",SPHEROID["Krasovsky_1940",6378245.0,298.3]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
	PROJECTION["Gauss_Kruger"],
		PARAMETER["False_Easting",-19.1],
		PARAMETER["False_Northing",4.9],
		PARAMETER["Central_Meridian",37.4998],
		PARAMETER["Scale_Factor",1.0],
		PARAMETER["Latitude_Of_Origin",55.667],
	UNIT["Meter",1.0]]

Собственно он лежит в Moskow_CS.rar
больше добавить ничего не смогу

[Boris]

тогда поподробнее вот эту строку

результат преобразования сравниваю c вариантом сделанным в ArcGis ошибка примерно 50 метров.

что вы делаете с OGR вы расписали, но что делаете в ArcGIS - нет.

[gobsy]

Виноват. Исправляюсь.
Преобразование тоже проводится в два этапа.
Этап №1
Сначала из WGS84 в Pulkovo_1942 по 12 зоне.
Открываем ArcGis идем в ArcToolbox выбираем там "Data Management Tools" раскрываем его, выбираем "Projections and Transformations" раскрываем его, выбираем "Feature", выбираем "Project".У этот путь выглядит как на прикрепленном рисунке №1.
Далее в появившемся окне в поле "Input Dataset or Feature Class" вставляем reserve_point.shp
В поле "Output Dataset or Feature Class" заполняем ручками название преобразованного слоя,можно не заполнять ArcGis сам туда что-нибудь сгенерирует.
В поле "Output Coordinate System" выбираем Pulkovo_1942.prj
В поле "Geographic Transformation" выбираем Pulkovo_1942_To_WGS_1984_12. В итоге должно получиться как на прикрепленном рисунке №2.
Нажимаем "OK" и получаем результат.

Этап №2.
Из Pulkovo_1984 в Московскую СК.
Открываем, если уже успели закрыть, ArcToolbox выбираем там "Data Management Tools" раскрываем его, выбираем "Projections and Transformations" раскрываем его, выбираем "Feature", выбираем "Project".У этот путь выглядит как на прикрепленном рисунке №1.
Далее в появившемся окне в поле "Input Dataset or Feature Class" вставляем только что сгенерированный reserve_point_Project2.shp
В поле "Output Dataset or Feature Class" заполняем ручками название преобразованного слоя,можно не заполнять ArcGis сам туда что-нибудь сгенерирует.
В поле "Output Coordinate System" выбираем Moscow_CS.prj.
В итоге должно получиться как на прикрепленном рисунке №3.

Результат сравнения приведен на Рисунке №4. Красный трек - преобразование в ArcGis, синий трек - преобразование в ogr2ogr.

[gobsy]

Сравнение полученных треков
на рисунке №5 представлена информация преобразовании сделанном в ogr2ogr
на рисунке №6 представлена информация преобразовании сделанном в ArcGis

[Boris]

Надо подумать, но мой опыт перевода координат с использованием ArcGIS и OGR (Proj4) говорит от том, что детали трансформации, задаются не только и не столько исходной и конечной СК, но и параметрами "Geographic Transformation". Вот что именно они задают и как их можно извлечь - советую спросить в ветке по АркГИС.
Я пробовал повторить трансформацию из WGS в локальную СК на основе СК-42, выполняющуюся в АркГИС со строкой "Geographic Transformation", в других продуктах - результат не удалось свести.

[gobsy]

Спасибо. Если у меня получится решить эту задачу, обязательно отпишусь.

[svt]

Здравствуйте, проблема решена была на прошлой неделе, но руки не доходили сюда написать.
Приемлемое решение нашли через конвертацию огром 7-и параметрическим методом со значениями
dx dy dz rx ry rz ds
40.187 -110.741 -123.39 1.41378 -0.78147 0.7348 0.6121

Отклонения точек в сравнении с WGS84 от 6 мм до 5 см на карте.

[Boris]

Отклонения точек в сравнении с WGS84 от 6 мм до 5 см на карте.

На карте или на местности? 5 см что-то многовато для карты.

[svt]

Для нашей задачи мы сошли погрешность приемлемой. К тому же сроки поджимали.
P.S. 5 см и "до 5 см" несколько разные понятия. В среднем получается от 1 до 2 см, что вполне нормально. В некоторых случаях, даже меньше 1 см и составляет порядка 6-7 мм.

[Boris]

а масштаб то какой?

[svt]

а масштаб то какой?

На этот вопрос, я Вам не смогу ответить. Все отображается в электронном виде, поэтому масштаб варьируемый. Все тестовые образцы при сравнении WGS 84 И Moscow_CS показали результаты точку в точку... А это местность примерно 200 на 200 м на листе формата А3. Спасибо за участие.

[Донецков]

PRJ файл для Москвы несколько другой, чем указанный:
PROJCS["Moskow",GEOGCS["MSK-77 Moskow",DATUM["<D_Moskow>",SPHEROID["<Bessel_1841>",6377397.0,299.15]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",37.5],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",55.66666666666666],UNIT["Meter",1.0]]

[Boris]

PRJ файл для Москвы несколько другой, чем указанный:

А точность какая? А положение эллипсоида Бесселя (?) по отношению к WGS84? А чем измеряли?