Пересчет XYZ -> XYZ (PROJ4/CS2CS)

[juffin_h]

До этого момента с proj4и cs2cs не работал. Не могу уговорить cs2cs пересчитать координаты XYZ в XYZ с разными датумами. Всякие варианты работают. BL-BL, BL-xy, BL-XYZ и .т.д. Но координаты XYZ - XYZ никак. Или не могу уговаривать, или …? Подскажите кто может.

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

Как-то так:

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

cs2cs +proj=cart +datum=WGS84 +to +proj=cart +towgs84=0,0,-1 xyz_0.tsv > xyz_1.tsv

[juffin_h]

Как-то так:
cs2cs +proj=cart +datum=WGS84 +to +proj=cart +towgs84=0,0,-1 xyz_0.tsv > xyz_1.tsv

Я тоже думал, что как то так. Но не хочет. Пишет:

cannot initialize transformation
cause: (null)
program abnormally terminated

пробовал так

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

cs2cs +proj=cart +towgs84=0,0,0 +to +proj=cart +towgs84=0,0,-10

Не ругается, но вход и выход полностью совпадают при любых значениях +towgs84

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

У меня почему-то хочет:

[gamm]

под уиндовсом - не хочет. В latlon тоже не хочет

Буфер обмена-1.png (10.42 КБ) 9535 просмотров

[juffin_h]

У меня под win такая картина
Печально.

[gamm]

версия другая, ща попробую накатить ... не получилось. Насколько я понимаю, cs2cs есть только в QGIS, в GDAL я его не нашел, но у последнего QGIS тот же 5.2.0. Или нужно ставить отдельно PROJ?

[juffin_h]

PROJ входит в GDAL. B cs2cs там же. По крайней мере я GDAL поставил и там все было.

[gamm]

значит, у меня какой-то огрызок стоит. Поставить полный OSGeo не получилось, валится. Может, Игорь Белов найдет, почему не работает ...

[juffin_h]

С cs2cs не получилось, решил попробовать с proj через pipeline. И на тебе. Вообще не понял, что proj выдает. cct при этом, выглядит адекватно

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

C:\Program Files\GDAL>proj +proj=cart
17 45 130
4320193.64      1320815.75 130

C:\Program Files\GDAL>cct +proj=cart
17 45 130 2019.0
 4320281.5502   1320842.6224  4487440.3326     2019.0000