Нарезка на тайлы генштабовских карт

[edoroe]

Добрый вечер !

Есть карты генштаба, привязанные в формате OziExplorer (*.map)
Пытаюсь нарезать на тайлы четыре рядом стоящие карты.
1) конвертирую
gdalwarp -t_srs EPSG:900913 -of GTiff -r near M-36-083.map M-36-083.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -r near M-36-084.map M-36-084.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -r near M-36-095.map M-36-095.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -r near M-36-096.map M-36-096.tiff

2) объединяю
gdal_merge.py -o merged.tiff M-36-083.tiff M-36-084.tiff M-36-095.tiff M-36-096.tiff

Результат:


Можно как-то убрать поля по краям карт ?
Информация в *map файлах ози о реальном размере карты есть.

[Kandey]

Может, проще взять другой инструмент, который позволяет сделать требуемое без лишних танцев? И если не стоит задача самому порезать эти четыре листа, можно использовать украденное до нас уже готовое решение. )

[edoroe]

За ссылку - благодю !
Задача стоит именно нарезать самому.

[Kandey]

Задача стоит именно нарезать самому.

Я для этого пользуюсь GM, в нём это довольно просто, но навязывать не буду ). На этом форуме могут неправильно понять (

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

Для обрезки используйте параметр -te или связку -cutline + -crop_to_cutline. Например, для этих листов:

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

gdalwarp -t_srs EPSG:900913 -of GTiff -te 3895999.94 6388729.86 3951660.43 6446259.24 -r near M-36-083.map M-36-083.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -te 3951660.43 6388731.27 4007320.94 6446260.67 -r near M-36-084.map M-36-084.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -te 3896001.39 6331593.34 3951661.88 6388729.86 -r near M-36-095.map M-36-095.tiff
gdalwarp -t_srs EPSG:900913 -of GTiff -te 3951661.88 6331594.74 4007322.38 6388731.27 -r near M-36-096.map M-36-096.tiff

А Kandey дело говорит.

[edoroe]

Для обрезки используйте параметр -te или связку -cutline + -crop_to_cutline. Например, для этих листов:

То что нужно. Благодарствую !
Информация о координатах вырезки получается из tiff файла по gdalinfo.
А как её получить до создания tiff файла ?

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

Как вариант, скачайте нужную разграфку отсюда и снимайте координаты краёв листов в целевой СК.

[edoroe]

Как вариант, скачайте нужную разграфку отсюда и снимайте координаты краёв листов в целевой СК.

я получаю требуемые координаты с
gdalinfo M-36-083.map
но там они в другой системе координат (по моему в EPSG:4326)
соственно, вопрос сводится к
какой gdal утилитой преобразовать координаты точки с EPSG:4326 в EPSG:900913 ?

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

там они в другой системе координат (по моему в EPSG:4326)

Вряд ли. Должно быть 4284.

какой gdal утилитой преобразовать координаты точки с EPSG:4326 в EPSG:900913 ?

Попробуйте gdaltransform.

[edoroe]

Может, проще взять другой инструмент, который позволяет сделать требуемое без лишних танцев? И если не стоит задача самому порезать эти четыре листа, можно использовать украденное до нас уже готовое решение. )

Очень полезное "готовое решение".
А как с него можно сзабирать тайлы, указывая номера тайлов (API) ?
Что-то нигде не смог найти ответ на данный вопрос.

[Kandey]

[Очень полезное "готовое решение".
А как с него можно сзабирать тайлы, указывая номера тайлов (API) ?
Что-то нигде не смог найти ответ на данный вопрос.

Частичный ответ тут

[SergeyRyzhkov]

Полный ответ
http://gis-lab.info/qa/forum-rules.html
касается как самого ресурса, так и способов "кражи" данных

[edoroe]

Полный ответ
http://gis-lab.info/qa/forum-rules.html
касается как самого ресурса, так и способов "кражи" данных

Да никто не собирается ничего красть, даже то, что подкинуто намерено для провоцирования краж (90% ресурсов инет работают по такому принципу)

[trir]

Да никто не собирается ничего красть

Уже

Есть карты генштаба, привязанные в формате OziExplorer (*.map)

их из спец части выносить нельзя, даже в бумажном виде!

[edoroe]

их из спец части выносить нельзя, даже в бумажном виде!

я вижу что нельзя - к примеру, упомянутое выше "готовое решение" и ещё с десяток ресурсов где можно эти карты скачать. Очевидно срок прошел, карты ведь 42 года.