У меня есть класс пространственных объектов (графства Южная Каролина, поэтому довольно большой географический район) в NAD83 SC State Plane. Его необходимо преобразовать во вторую проекцию (NAD83 UTM 17), а затем преобразовать обратно в оригинал. Я буду использовать инструмент проекта Esri для достижения этой цели.
Может ли это двойное преобразование вызвать сдвиг в расположении координат многоугольников и на сколько - сантиметры, метры, километры?
arcgis-desktop
coordinate-system
united-states
вересковые
источник
источник
Ответы:
Я не знаю, какой проекционный движок использует ArcGis, но очень интересный вопрос и для proj.4. Поэтому я попытаюсь протестировать проекционный механизм proj.4 в среде GNU-R. Я использую углы NAD 83 - UTM 17 и EPSG 26917, рекурсивно повторяю их 10000 и 1000000 раз и вычисляю разницу с начальными значениями.
Вот результаты:
Кажется, что ошибка «перепроецирования» в пределах сантиметрового диапазона для 10000 петель.
И вырастет до ошибки в диапазоне метров, если вы выполните цикл 1000000 раз.
Вот сценарий.
Дальнейшие тесты в среде статистики должны быть простыми. Сценарии и объяснение кода для среды linux доступны по адресу github.com/bigopensky .
источник
WGS84 Bounds: -84.0000, 24.0000, -78.0000, 83.0000
представляют собой интересную область. Я сделал ошибку?Esri имеет собственный проекционный движок.
Большинство методов проекций и географических / базовых преобразований хорошо работают при использовании в соответствующей области интересов. Если вы выходите слишком далеко за пределы зоны UTM, поперечный Меркатор не всегда точно «инвертирует» (конвертирует в широту-долготу). Проекции, используемые для всего мира, могут иметь некоторые проблемы на полюсах или вокруг них, или на меридиане +/- 180, или на «анти-меридиане» (меридиане, который находится напротив центра проецируемой системы координат).
Я пробежал 4 очка, которые падают за пределы Южной Каролины через проекционный двигатель Esri. Для стресс-теста в 1k, 10k или 1M баллов мне придется что-то кодировать, поскольку мой существующий аналогичный тест просто выполняет «круговую передачу» - от географического до прогнозируемого. 32133 - НАД 1983 года, государственный самолет Южной Каролины (в метрах). 26917 - NAD 1983 UTM, зона 17 север.
Итак, вы можете видеть, что у нас было две точки, которые вернулись в 10e-09.
Обработка в ArcGIS усложняется тем, что существует пространственная привязка. Пространственная привязка включает систему координат, а также некоторые значения для хранения и анализа. По умолчанию системы координат, использующие метры, хранятся с точностью до одной десятой миллиметра, 0,0001.
Раскрытие: я работаю на Esri.
источник
Я думаю, что это тот случай, когда вам нужно протестировать предложенный рабочий процесс с некоторыми функциями контрольных точек, к которым легко добавить поля координат XY.
Сравните значения XY ваших начальных точек с теми, которые вы проецировали / трансформировали (хотя и много раз), и вы будете количественно определять разницу.
источник