Координатные преобразования в ArcGis

Question

Насколько корректно ведутся координатные преобразования в ArcGis? Простой пример берем набор точек в Гауссе-Крюгере(GK) и пересчитываем в UTM/WGS84 через GeoCalculator от Blue Marple или через Erdas. Загружаем эти файлы в ArcGis и видим, что разлет между точками GK и UTM составляет порядка 90м. Складывается впечатление, что координатные преобразования при переходе с одного эллипсоида на другой ведутся некорректно. Набор параметров при описании эллипсоида в ArcGIS ограниченный или я не знаю, где можно более подробно описать датумы. Может где-нибудь можно ввести новый эллипсоид , как в Erdas (c расширеным набором параметров)->н-р «Krasovsky» {
35 6378245.0 6356863.019
«Krasovsky» 0 0 0 0 0 0 0
«Pulkovo 1942» 27.0 -135.0 -84.5 0 0 -2.685867e-06 2.263e-07
}?

Answer 1

И преобразования ведутся корректно и сфероид можно описать с датумом и задать параметры перехода между сфероидами.

Только нужно не путаться с терминами.

Сфероиды и Красовского, и WGS84 (вместе с датумами) уже есть в в предопределённых системах координат

см. ArcGIS\Coordinate Systems\Geographic Coordinate Systems\Spheroid-based\Krasovsky 1942.prj и ArcGIS\Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj

А что касается следующих цифирек: 27.0 -135.0 -84.5 0 0 -2.685867e-06 2.263e-07 - то это параметры преобразования между сфероидами.

Т.к. их семь, то это или преобразования по методу или Position Vector, или Coordinate Frame - разница в знаках коэффициентов

Для преобразования из GCS_Pulkovo_1942 в WGS84 в ArcGis 9 есть 15 предопределённых наборов преобразования. Те цифры, которые Вы указали в них не значатся. Следовательно задайте их сами при проецировании.

Если Вы просто отображаете данные в ArcMap и выполняете проецирование "на лету", то это делается в свойствах фрейма данных. На закладке "Coordinate System" (Системы координат) нужно нажать кнопку Transformation (Преобразования). Затем нажать кнопку "New" (Новая) и ввести свои параметры преобразования. Необходимо учитывать, что преобразования действуют только на один сеанс работы ArcMap. см http://support.esri.com/index.cfm?fa=knowledgebase.techarticles.articleShow&d=17430

 

Answer 2

90м многовато для разницы в наборе элементов трансформирования - похоже, что трансформации датума вообще не происходит

Answer 3

Это Вам такажется.

Для того и предназначены преобразования между сфероидами, чтобы свести ошибки преобразований координат к миниуму. Точность преобразования между WGS84 и Красовского 1940 (по умолчанию - преобразование по трём параметрам) для всей России - порядка 100-200м. Не зря в поставке идёт 15 встроенных вариантов - они настроены на конкректные регионы, внутри которых точность может быть лучше (1-2 метра).

Вы же сами сдвигаете для чего-то центр сфероида на 27.0 -135.0 -84.5 метров. А ArcGIS его двигает по умолчанию в другую сторону

А чтобы у Вас были одинаковые результаты преобразований в Erdas и в ArcGIS нужно задать одинаковые методы преобразований и коэффициенты. Формулы преобразований у них одинаковые - значит координаты должны совпадать.

 

Answer 4

Я немного некорректно выразился.

90м означает, судя по величине цифры, что исходная и пересчитанная точка - в разных системах координат (например Pul и WGS84), а не, скажем разница в наборе элементов трансфомирования (Pul1 и Pul2). И Евгений получает именно то, что должно получаться, если его ГК на эллипсоиде Красовского.

Answer 5

Уважаемый Григорий, Вы отвечаете:

"Для преобразования из GCS_Pulkovo_1942 в WGS84 в ArcGis 9 есть 15 предопределённых наборов преобразования. Те цифры, которые Вы указали в них не значатся. Следовательно задайте их сами при проецировании."

Я хочу перепроецировать не на лету, а через arctoolbox - в Ердасе достаточно видоизменить файл sheroid.tab - а в аркгис, где нужно задавать параметры, чтобы выполнять перепроецирование по семи параметрам (не в arcmap)?

Answer 6

Для предопределённых способов преобразования нужный достаточно выбрать из списка в инструменте project (проецирование).

Если нужных Вам параметров нет в списке, то их придётся описывать в файле geogtran. Файл нужно создать и установить переменную окружения PEOBJEDITHOME на каталог в котором этот файл находится. Содержание файла должно выглядеть примерно так:

GEOGTRAN,208001,"WGS_1984_To_Pulkovo42_My",PE_GCS_WGS_1984, PE_GCS_PULKOVO_1942,PE_MTH_POSITION_VECTOR,PE_PAR_X_AXIS_TRANSLATION, 27.0,PE_PAR_Y_AXIS_TRANSLATION,-135.0,PE_PAR_Z_AXIS_TRANSLATION,-84.5,PE_PAR_X_AXIS_ROTATION,0.0,PE_PAR_Y_AXIS_ROTATION,0.0, PE_PAR_Z_AXIS_ROTATION,0.0,PE_PAR_SCALE_DIFFERENCE,0.0, PE_PAR_X_COORD_OF_ROTATION_ORIGIN,0.0,PE_PAR_Y_COORD_OF_ROTATION_ORIGIN, 0.0, PE_PAR_Z_COORD_OF_ROTATION_ORIGIN,0.0

Это всё одной строкой для одного метода преобразования.

Более подробно об описании своих проекций, параметров преобразования и пр. можно прочитать здесь:

http://arcsdeonline.esri.com/Support_files/capi_concepts/using_persistent_user.htm

Сразу замечу - ArcGIS после установки этой переменной будет медленнее восстанавливать списки проекций. Для коррекной работы пользовательских параметров преобразования необходим SP2. В качестве разделителя десятичных знаков в системе нужно установить точку. Ну и сами параметры нужно аккуратно устанавливать.

Answer 7

Спасибо большое за ответ. Будем пробовать.

Answer 8

2 Григорий ("Более подробно об описании своих проекций, параметров преобразования и пр. можно прочитать здесь: http://arcsdeonline.esri.com/Support_files/capi_concepts/using_persistent_user.htm")

А где сейчас можно найти эту страничку? По этому пути уже ничего нет... :((( А то я вообще не понимаю, как и куда установить эту переменную... :(((

Answer 9

Как устанавливать переменную окружения PEOBJEDITHOME см.
https://forum.esri-cis.ru/index.php?qa=25144

Но для 9.2 лучше описывать свои преобразования с помощью инструмента ArcToolBox->Управление данными->Проекции и преобразования->Создать пользовательское географическое преобразование.
http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Create_Custom_Geographic_Transformation_(Data_Management)

    

Answer 10

Спасибо, пошел изучать ))