Тривиальная конвертация DATUM-а в WGS84

[Longy]

Здравствуйте!

О DATUM-е исходной с/к известно, что lat/lon-ы в нем связаны с lat/lon-ми WGS84 оффсетом (lat_wgs84=lat+ofs_lat, lon_wgs84=lon+ofs_lon). Можно ли через API GDAL-а сконвертировать точки из исходной с/к в с/к с WGS84 datum?

Возможно, есть способ пересчитать эти офсеты в 7 параметров TOWGS84? Или можно задать сетку трансформации из 4-х точек на весь мир для исходного датума, учитывающую эти смещения?

Буду рад любым идеям.

Спасибо!

[Игорь Белов]

Очевидно, имеется в виду ограниченная территория, иначе представленное преобразование широты вырождается вблизи полюсов. Самое простое и правильное решение — использовать сдвиг как есть:

φ_wgs84 = φ + ∆φ
λ_wgs84 = λ + ∆λ

Если же непременно хочется заморочиться, можно придумать over 9000 способов представления этого локального преобразования каким-нибудь глобальным, только следует помнить, что точность преобразования падает по мере удаления от территории. Вот первый пришедший на ум вариант семи параметров:

∆X = 0
∆Y = 0
∆Z = −M₀ cos φ₀ ∆φ
ωx = 0
ωy = 0
ωz = −∆λ
m = 0

Здесь φ₀ — широта в середине территории, M₀ — радиус кривизны меридиана; все угловые величины выражены для простоты в радианах. Для GDAL и других программ, использующих PROJ.4, меняем знак разворота: ωz = +∆λ.