Форумы GIS-Lab.info

Уважаемые участники обсуждения !
Кто-нибудь когда-нибудь сталкивался с библиотекой PROJ4 для Delphi, выложенной на :
Proj4 for Delphi :
http://www.naadsm.org/iproj4-delphii
или, может быть, что-то про нее слышали ?
Можете ли вы высказать какое-то определенное мнение о работе и качествах этой библиотеки ?
Меня интересует вопрос о том, насколько ей можно доверять, ну и прочие самые разные ее аспекты ...
короче - хотел бы услышать просто мнения от тех, кто имел к ней хоть какое-то касательство.

Я скачал с указанного выше сайта архивы этой библиотеки, скомпоновал прект демо-приложения
этой библиотеки в Delphi ( это прошло относительно легко ) и опрбовал его на тех данных, которые обсуждались
в теме "Привязка растровой карты Арктики"
а именно - для проекции : и для точек ( градусы Lon/Lat --> метры проекции ) : Созданое приложение выдало следующее : Как видите - очень хорошее совпадение результатов.

Но, в то же время, когда я попробовал по аналогии сравнить результаты преобразований для одних и тех же проекций ( самых разных ) в этом приложении и на сайте http://cs2cs.mygeodata.eu/
то впечатления остались весьма неприятные - очень существенные расхождения.
Вот и хочу спросить - у кого какие мнения насчет этой библиотеки ...
Если кого-то интересуют скомпилированное мною упомянутое приложение или весь его исходный проект в Delphi с исходными кодами - пишите сюда, выложу здесь.

Как правило, ищут там где потерял, а не там где светлее.
Выложите примеры очень существенных расхождений.

stout писал(а):Выложите примеры очень существенных расхождений.

Вы правы, надо было это сделать с самого начала. Приношу извинения.
Результат преобразования из WGS84 в Pulkovo 1942, полученный на http://cs2cs.mygeodata.eu/ : Результат преобразования, полученный из оттранслированного мною проекта библиотеки Proj4 for Delphi,
полученного с http://www.naadsm.org/iproj4-delphii : Если в первом случае результат получен в градусах, то во втором - в метрах проекции.
Причем, как мне кажется - результат верный, поскольку при инверсном преобразовании получаем : Возникает вопрос - почему в первом случае результат в градусах, а во втором - в метрах ?
Что надо изменить в строке proj=... ( или где еще ... ), чтобы приложение выдало результат в градусах ?

Если кому интересно - выкладываю здесь прикрепленный архив с оттранслированными мною
демо-приложениями библиотеки Proj4 for Delphi. Все необходимые пояснения - в файле README.txt .

Какие ять метры - это радианы!!!

trir писал(а):Какие ять метры - это радианы!!!

Опередили на 2 минуты
Но расхождение на 7.8" для первого числа есть.

trir писал(а):Какие ять метры - это радианы!!!

Да, вы правы - и в самом деле радианы ! Прошу прощения - первоначально не обратил внимания. У меня сегодня опять очень тяжелый день ...
Как я понимаю, строка Proj= здесь ни при чем. Возможно, что само приложение результат преобразует в радианы.
Надо лезть смотреть в исходные коды проекта. Сегодня этим заниматься некогда, зделаю завтра.

stout писал(а):Но расхождение на 7.8" для первого числа есть.

ставлю на то, что где-то использовался float 4-х байтовый, ошибка округления в 6 знаке.

gamm писал(а):

stout писал(а):Но расхождение на 7.8" для первого числа есть.

ставлю на то, что где-то использовался float 4-х байтовый, ошибка округления в 6 знаке.

Вчера вечером я бегло просмотрел исходные коды проекта и протрассировал его в отладчике. Дело не в использовании типа float, а в том, что по всему пространству кода используются функции преобразования вещественных форматов, в которых для вещественных значений по умолчанию задано ограничение в 6 знаков после запятой. Попытался увеличить это число - ничего не дало. В кодах действует какая-то хитрая закавыка, с которой, к сожалению, нет возможности разбираться в ближайщие дни. Вернусь к этому вопросу несколько попозже.
Ну, а для чего в проекте установлено подобное ограничение на точность - интересный, конечно, вопрос, но, боюсь, мы здесь можем только гадать из пустого в порожнее ...
Однако, возможно, надо учитывать - проект по свободной лицензии ...

Ничего такого не заметил https://github.com/OSGeo/proj.4
Где именно происходит округление?
Везде вроде double (https://github.com/OSGeo/proj.4/blob/ma ... proj_api.h)

Дмитрий Барышников писал(а):Ничего такого не заметил https://github.com/OSGeo/proj.4
Где именно происходит округление?
Везде вроде double (https://github.com/OSGeo/proj.4/blob/ma ... proj_api.h)

Я имею в виду не ваши источники, а исходные тексты, которые приведены :

trir писал(а):здесь

Для вещественных чисел везде используется 8-байтовый тип Double, но при выводе вещественных значений используется ограничение по умолчанию - 6 цифр после запятой.
И, кстати, по поводу ваших источников :

Дмитрий Барышников писал(а): https://github.com/OSGeo/proj.4
https://github.com/OSGeo/proj.4/blob/ma ... proj_api.h

не будете ли вы так любезны хотя бы вкратце пояснить, что это за источники и каково их происхождение ?

trir писал(а):здесь

Для полной сборки проекта этого мало. Надо закачать файлы :
cxreedit-1.9.3-src.zip
general_purpose-1.6.4-src.zip
qclasses-0.1.7-src.zip
proj4_delphi-1.0.0-src.zip
с http://www.naadsm.org/opensource

не будете ли вы так любезны хотя бы вкратце пояснить, что это за источники и каково их происхождение ?

мда, это и есть Proj4, а "Proj4 for Delphi" - это его порт, причём старый

trir писал(а):мда, это и есть Proj4, а "Proj4 for Delphi" - это его порт, причём старый

коллеги, а что такое "порт для Delphi", если в этом самом Delphi можно без затей использовать dll от самого PROJ4? это какая-то специальная фича для школьников?

Дмитрий Барышников писал(а):Везде вроде double (https://github.com/OSGeo/proj.4/blob/ma ... proj_api.h)

так вы же смотрите сам PROJ4, а не его "порт для Delphi", сделанный неизвестными пионэрами ...