Как происходят преобразования в ГИС карта 2011?

[evgen_sobolev]

Здравствуйте! Использую преобразование координат через ГИС карту 2011 для тестирования своего ПО. У меня преобразование координат происходят по средства WKT через GDAL. Никто не знает как это делается в ГИС карта 2011? А то расхождение при тестирование в 1 минуту меня смущают(

[Boris]

телепатический сеанс не получился. матрица возможных вариантов вопроса имеет размерность 20х20х20x...x.... .
для начала начните с очевидного: исходная СК, конечная СК, описание трансформации в ГИС, параметры трансформации, если они вам известны.
предположить, что там математика разная, стоит в самую последнюю очередь. а вот параметры поискать в первую. особенно в GDAL.

[evgen_sobolev]

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

[Boris]

Я уверен, что ГИС-карта не использует библиотеки GDAL, но вероятность того, что она использует сильно различающиеся способы пересчета - минимальна. Единственное отличие может быть в реализации алгоритмов расчета рядов, но и оно не может быть вопиющим - в 1 минуту.
Сперва надо найти параметры исходных преобразований - самих проекций и их эллипсоидов. GDAL может брать эти параметры из параметров по умолчанию, в нескольких местах самого пакета. Поскольку GDAL использует в свою очередь PROJ4 для географических расчетов, то лучше сразу начинать проверку с PROJ4. В утилите командной строки cs2cs есть параметр -v, который выведет все использованные параметры исходной и начальной проекции.

[Boris]

из цилиндрической Миллера на Кларк 1880 в WGS 84.

для примера пересчитаем точку на эллипсоиде Кларк 1880 с координатами (50;50) в точку на WGD84 через проекцию Miller Cylindrical (+proj=mill)

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

echo 50 50 | cs2cs +proj=latlong +ellps=clrk80 +to +proj=mill  +ellps=clrk80  +units=m +no_defs | cs2cs -v +proj=mill +ellps=clrk80 +units=m +no_defs +to +proj=latlong +ellps=WGS84 +datum=WGS84 +units=m +no_defs

дает

# ---- From Coordinate System ----
#Miller Cylindrical
# Cyl, Sph
# +proj=mill +ellps=clrk80 +units=m +no_defs
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
#
# +proj=latlong +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0
#--- following specified but NOT used
# +ellps=WGS84
50dE 50dN 0.000

что означает, что стандартно переход выполнялся с расчетом того, что центры эллипсоидов CLarck-1880 и WGS84 - совпадают, и угловые координаты (50;50) НЕ изменились! Это умолчания PROJ4, умолчания GDAL могут быть иными, а условия перехода между эллипсоидами в Панораме - третьими.

Вот сколько датумов на Clarke 1880 знает http://www.globalmapper.com/helpv13/datum_list.htm Global Mapper