Привязка карт с поехали орг

[Moskus]

Сделал тестовый набор geotiff-файлов, чтобы хотя бы с исходными данными в тесте проблем не было. 8bit, сжаты LZW, с прозрачным фоном, лист "миллионки" N-37 и лист поменьше, входящий в его состав. Проекции для каждого: Гаусс-Крюгер 7я зона/Пулково-1942 = EPSG:28407, Меркатор/WGS84 = EPSG:3395, Equirectangular со стандартной параллелью, равной центральной для листа (54 и 55.5 градусов соответственно). <A HREF="http://www.sendspace.com/file/tmbxa4">h ... /A> <br>

[seno]

По поводу улетевших карт. Исходников не нашел, но нетрудно было воспроизвести ситуацию. Сделал три Geotiff из внимательно проверенного по привязке (и по восточному смещению в частности) листа ГШ 500м/см. Lat/Long Mercator Equirectangular c настройкам по рекомендации Moskus После конвертации в map2jnx от автора QLandkarte, только Lat/Long сохранила привязку. Две остальные карты уехали по Y. Кроме того, JNX от Equirectangular оказалась вытянута вдоль X примерно в два раза относительно вида этой проекции в GM. Mercator в JNX внешне был похож на Geotiff первоисходник. Просмотр в QLandkarteGT, прибора в руках нет. Отконвертировал в JNX набор карт от Moskus. С Equirectangular та же беда. т.е. улетели координаты по Y и вытянутость вдоль X. Удивил Mercator, его привязка (по крайней мере центр карты) осталась на месте. И вид сетки похож на исходный. С цветами конечно беда полная, из-за индексированного по цвету исходника, но в данном случае это побоку. Короче, непонятно. Какую брать проекцию исходника? Ну неправильно использовать Lat/Long в данном случае. А Equirectangular пока никак. С Mercator ввиду неоднозначности, сделаю попытку внимательней посмотреть. Может кто общается с автором и может разогнать туман? PS. Может туплю в использовании инструмента? На всякий случай - используемый конвертер http://ifolder.ru/22090885

[seno]

Metadata из GM по трем проекциям одного листа: FILENAME=D:\Map2Jnx\LL.tif DESCRIPTION=LL.tif UPPER LEFT X=33.9686962872 UPPER LEFT Y=61.8522296753 LOWER RIGHT X=34.5247869251 LOWER RIGHT Y=61.6406957775 WEST LONGITUDE=33.96869629╟ E NORTH LATITUDE=61.85222968╟ N EAST LONGITUDE=34.52478693╟ E SOUTH LATITUDE=61.64069578╟ N PROJ_DESC=Geographic (Latitude/Longitude) / WGS84 / arc degrees PROJ_DATUM=WGS84 PROJ_UNITS=arc degrees EPSG_CODE=4326 COVERED AREA=687605216 sq m NUM COLUMNS=9477 NUM ROWS=3605 NUM_BANDS=3 PIXEL WIDTH=0.0000587 arc degrees PIXEL HEIGHT=0.0000587 arc degrees PHOTOMETRIC=RGB Full-Color BIT_DEPTH=24 ROWS_PER_STRIP=1 COMPRESSION=None PIXEL_SCALE=( 5.86779e-005, 5.86779e-005, 1 ) TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( 33.9686962872, 61.8522296753, 0.0000000000 ) MODEL_TYPE=Geographic lat-long system RASTER_TYPE=Pixel is Area FILENAME=D:\Map2Jnx\Equ-lar.tif DESCRIPTION=Equ-lar.tif UPPER LEFT X=1786800.133 UPPER LEFT Y=6873816.060 LOWER RIGHT X=1816051.509 LOWER RIGHT Y=6850309.107 WEST LONGITUDE=33.96869629╟ E NORTH LATITUDE=61.85222968╟ N EAST LONGITUDE=34.52479156╟ E SOUTH LATITUDE=61.64070853╟ N PROJ_DESC=Equirectangular / WGS84 / meters PROJ_DATUM=WGS84 PROJ_UNITS=meters COVERED AREA=687610833 sq m NUM COLUMNS=4532 NUM ROWS=3642 NUM_BANDS=3 PIXEL WIDTH=6.45 meters PIXEL HEIGHT=6.45 meters PHOTOMETRIC=RGB Full-Color BIT_DEPTH=24 ROWS_PER_STRIP=1 COMPRESSION=None PIXEL_SCALE=( 6.45441, 6.45441, 1 ) TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( 1786800.133, 6873816.060, 0.000 ) MODEL_TYPE=Projection Coordinate System RASTER_TYPE=Pixel is Area FILENAME=D:\Map2Jnx\M.tif DESCRIPTION=M.tif UPPER LEFT X=6607834.777 UPPER LEFT Y=8786475.170 LOWER RIGHT X=6669739.020 LOWER RIGHT Y=8736809.719 WEST LONGITUDE=33.96869629╟ E NORTH LATITUDE=61.85222968╟ N EAST LONGITUDE=34.52479156╟ E SOUTH LATITUDE=61.64071480╟ N PROJ_DESC=Mercator / WGS84 / meters PROJ_DATUM=WGS84 PROJ_UNITS=meters COVERED AREA=687548882 sq m NUM COLUMNS=4532 NUM ROWS=3636 NUM_BANDS=3 PIXEL WIDTH=13.66 meters PIXEL HEIGHT=13.66 meters PHOTOMETRIC=RGB Full-Color BIT_DEPTH=24 ROWS_PER_STRIP=1 COMPRESSION=None PIXEL_SCALE=( 13.6594, 13.6594, 1 ) TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( 6607834.777, 8786475.170, 0.000 ) MODEL_TYPE=Projection Coordinate System RASTER_TYPE=Pixel is Area Тоже в конвертере: http://s44.radikal.ru/i106/1102/96/0ab1 ... br> <br>

[Moskus]

Я думаю что пока разработчик сам не покрутит разные примеры, все вот эти данные (они, к слову, не слишком информативны - выдача gdalinfo или listgeo дала бы больше информации, о проекции, например) гадать можно бесконечно и бессмысленно.

[AlexWhiter]

Могу связаться с автором, но сначала должен понять, что работает неправильно.

Скачал файлы Moskus'а.
В Global Mapper карта в Equirectangular сдвинута на запад относительно остальных. Что делаю не так?

[Moskus]

Гхм. Пофиксим, это я умудрился в какой-то древней версии GM вывод сделать, судя по всему. <A HREF="http://www.sendspace.com/file/iah0ok">h ... /iah0ok</A> - фикс

[seno]

<<<я умудрился в какой-то древней версии GM вывод сделать>>>
Geotiff в разных версиях GM может быть с разной степенью приближения к стандарту?
Ещё: надо экспериментировать на полноцвете. Скорей всего это неважно, но тем не менее. Вход 24 бита Geotiff будет правильно.

[AlexWhiter]

Теперь все три карты отображаются в GM одинаково.
Сконвертировал одну из карт в Lat/Lon, потом сделал JNX.

Загрузил все JNX в QLandkarte и поставил там-сям путевые точки на карте Lat/Lon.
Меркатор - выводится так же, как и Lat/Lon. Точки совпадают.
Пулково - точки и тайлы в центре карты совпадают с Lat/Lon, а по краям карта разваливается, расползается чем дальше от центра, тем сильнее.
Eq-Cyl - равномерно смещена на 10км на юго-запад.

Seno, у вас такая же картина?

[Moskus]

Ну были там какие-то старые баги. Просто у меня версий несколько, включая древние, мог промахнуться мимо шортката.
В 24 бита можно развернуть тем же GM, я это исключительно ради экономии объема. А то они и так здоровые.

[seno]

Мне кое что непонятно с самого начала. Давайте посмотрим исходники, перед тем как смотреть JNX и тем более автора мучить. Здесь таблица http://ifolder.ru/22109607 Пояснения: поставил GM12.01, открыл одну из карт от Moskus, меняя проекцию получил табличку. Смущают цифры помеченные красным.