proj4, lonlat -> XY

[drug]

На форуме нашел как осуществлять преобразование из географической СК в плоскую прямоугольную с помощью proj4:

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

 echo 19.7255859375 50.17515625 | cs2cs -f "%.5f" +init=epsg:4284 +to +init=epsg:28404
4408962.13881   5561205.19867 0.00000

И далее обнаружил, что если изменить географическую СК, то результат остается тем же:

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

 echo 19.7255859375 50.17515625 | cs2cs -f "%.5f" +init=epsg:4326 +to +init=epsg:28404
4408962.13881   5561205.19867 0.00000

Т.е. при пересчете:
а) из географических Гаусса-Крюгера в плоские Гаусса-Крюгера
б) из географических WGS84 в плоские Гаусса-Крюгера
получается один и тот же результат. Но ведь у Крассовского и WGS84 разные характеристики. Или в данном случае за датум результата всегда берется датум исходной СК? Можете прояснить немного?

[trir]

почитай внимательно описание epsg:28404

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

PROJCS["Pulkovo 1942 / Gauss-Kruger zone 4",
GEOGCS["Pulkovo 1942",
DATUM["Pulkovo_1942",
SPHEROID["Krassowsky 1940",6378245,298.3,
AUTHORITY["EPSG","7024"]],
AUTHORITY["EPSG","6284"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4284"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",21],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",4500000],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","28404"],AXIS["Y",EAST],AXIS["X",NORTH]]

[drug]

Почитал, но ничего нового не обнаружил.
Я исхожу из того, что при пересчете из географических в плоские прямоугольные используются параметры эллипсоида, в частности его оси. Я предпологал, что proj4 в начале переводит из географических WGS84, далее географические на Крассовском и уже потом плоские. А здесь получается, что при данном преобразовании географические координаты всегда берутся на Крассовском? (т.е. результируещем эллипсоиде)

З.Ы. Кажется догадываюсь - epsg:28404 это именно преобразование и есть?

[trir]

Да... был не прав, действительно странно...

[gamm]

Да... был не прав, действительно странно...

ничего странного, это прописано в доках, что при некоторых вызовах используется WGS84, что бы вы ни написали в +proj строке (прямо или косвенно). Внутри библиотеки PROJ4 нужно вызывать разные функции, в rgdal та же особенность (там два варианта вызова, кошерный - spTransform(), некошерный - project(), для командной строки тоже что-то подобное наверняка есть). Скорее всего, это "тяжкое наследие прошлого" тащится ...

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

Забавно. У меня иначе:

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

$ echo 19.7255859375 50.17515625 | cs2cs -f "%.5f" +init=epsg:4284 +to +init=epsg:28404
4408962.13881	5561205.19867 0.00000
$ echo 19.7255859375 50.17515625 | cs2cs -f "%.5f" +init=epsg:4326 +to +init=epsg:28404
4409090.39139	5561243.00529 -30.54567

[Александр Мурый]

У меня результат как у ErnieBoyd. cs2cs версии 4.8.0. Как именно работает cs2cs и proj.4, можно понять из исходного кода (благо что он открыт). Если уметь бегло читать на Си, конечно

[drug]

Проблема была в том, что я собрал proj4 из свежих сходников, однако бинарники из-за динамической линковки тянули proj.so из предустановленного пакета версии 4.7.1 (мы используем астра линукс, которые на дебиане7 основан со всеми вытекающими). Удалил старый пакет и make install решил мою проблему. За подсказку про версию спасибо, Александр, а вот про беглое чтение сишного кода это вы не к месту сказали, оно тут никаким боком.

[Александр Мурый]

а вот про беглое чтение сишного кода это вы не к месту сказали, оно тут никаким боком.

Не соглашусь с вами, т.к. именно чтение исходного кода позволяет лучше всего разобраться с работе программы. Если было не к месту, то пардоньте.

[drug]

Да я уже кучу кода прочитал и кучу написал - в начале не смог до конца протестить свою реализацию ГОСТа из-за отсутствия эталонных данных. Решил заюзать proj через gdal в качестве эталона, а у него тоже результаты неверные оказались. Все это заставило меня начать немного нервничать, т.к. время идет, а результата нет. Мне нужны всего лишь эталонные данные, чтобы убедится в правильности расчетов, а вы меня послали код читать, вот я в эмоциях и прокомментировал. Так что это вы меня пардоньте, если что.