Здравствуйте! Возникла проблема при пересчете в ГИС «Карта 2011» из системы Nord Sahara 1959 в WGS84. Параметры для проецирования используются следующие:
ΔX = -209,362198 m Rot X = 0,00461215 sec ΔY = -87,816200 m Rot Y = 3,47842207 sec ΔZ = 404,619830 m Rot Z = 0,58048472 sec Масштабный коэффициент: -1,4547220 ppm
Более подробно здесь В чем проблема и как ее решить?
Николай написал: Здравствуйте! Возникла проблема при пересчете в ГИС «Карта 2011» из системы Nord Sahara 1959 в WGS84. Параметры для проецирования используются следующие:
ΔX = -209,362198 m Rot X = 0,00461215 sec ΔY = -87,816200 m Rot Y = 3,47842207 sec ΔZ = 404,619830 m Rot Z = 0,58048472 sec Масштабный коэффициент: -1,4547220 ppm
Более подробно здесь В чем проблема и как ее решить?
Просьба написать каким образом получены эталонные координаты в WGS84 и Nord Sahara 1959.
Если посмотреть базу кодов EPSG версии 9, то для системы Nord Sahara 1959 в геодезической системе EPSG:4307 имеется 9 вариантов датума. Ваш вариант версия 9 "INCT-Dza" дает те значения, что вычислены в ГИС Карта 2011 и показаны в Вашем вложении. И это правильные значения. Вариант датума версия 8 "Tot-Dza Ahnet" дает те значения, что показаны в ArcGIS и они тоже правильные, но другие. Датум версии 8 имеет значения DX=-168.52, DY=-72.05, DZ=304.3. Остальные параметры равны 0.
Oleg Belenkov написал: Если посмотреть базу кодов EPSG версии 9, то для системы Nord Sahara 1959 в геодезической системе EPSG:4307 имеется 9 вариантов датума. Ваш вариант версия 9 "INCT-Dza" дает те значения, что вычислены в ГИС Карта 2011 и показаны в Вашем вложении. И это правильные значения. Вариант датума версия 8 "Tot-Dza Ahnet" дает те значения, что показаны в ArcGIS и они тоже правильные, но другие. Датум версии 8 имеет значения DX=-168.52, DY=-72.05, DZ=304.3. Остальные параметры равны 0.
Параметры проецирования: https://developers.arcgis.com/java/10-2/guide/datum-transformations.htm[TABLE][TR][TD]5660[/TD][TD]Nord_Sahara_1959_To_WGS_1984_9[/TD][TD]GEOGTRAN["Nord_Sahara_1959_To_WGS_1984_9",GEOGCS["GCS_Nord_Sahara_1959",DATUM["D_Nord_Sahara_1959",SPHEROID["Clarke_1880_RGS",6378249.145,293.465]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],METHOD["Position_Vector"],PARAMETER["X_Axis_Translation",-209.3622],PARAMETER["Y_Axis_Translation",-87.8162],PARAMETER["Z_Axis_Translation",404.6198],PARAMETER["X_Axis_Rotation",0.0046],PARAMETER["Y_Axis_Rotation",3.4784],PARAMETER["Z_Axis_Rotation",0.5805],PARAMETER["Scale_Difference",-1.4547]][/TD][/TR][/TABLE] Т.е,, все исходные данные идентичны.
Есть 2 способа пересчета по Гельмерту: "Position Vector transformation" EPSG:1033 и "Coordinate Frame Rotation" EPSG:1032. В базе данных кодов EPSG встречаются обе формы записи. Отличие заключается в том, что меняются знаки угловых поправок при вводе в программу вычислений. В ГИС Панорама реализован второй метод "Coordinate Frame Rotation" EPSG:1032. Для кода EPSG:4307 поправки заданы по варианту EPSG:1033. Поэтому при создании карты нужно поменять знаки угловых поправок и все сойдется.
Oleg Belenkov написал: Есть 2 способа пересчета по Гельмерту: "Position Vector transformation" EPSG:1033 и "Coordinate Frame Rotation" EPSG:1032. В базе данных кодов EPSG встречаются обе формы записи. Отличие заключается в том, что меняются знаки угловых поправок при вводе в программу вычислений. В ГИС Панорама реализован второй метод "Coordinate Frame Rotation" EPSG:1032. Для кода EPSG:4307 поправки заданы по варианту EPSG:1033. Поэтому при создании карты нужно поменять знаки угловых поправок и все сойдется.