gdal_merge оставляет "щели" между сшитыми кусками

[nextstopsun]

Всем привет!

Сшиваю кучу GeoTiff-ов утилитой gdal_merge.py. Сшивается все хорошо, вот только в выходном растре между сшитыми кусками остается зазор в 1 пиксель (см. приложенный скриншот).

Кто с этим сталкивался и как это побороть?
Буду признателен за любую информацию.

[gimran]

При сшивке добавьте два геотиффа в новый проект и сшивайте только два и т.д.
Покрайней мере мы так сшивали поканально растры, может и Вам поможет это в Вашем случае...

[nextstopsun]

К сожалению у меня 6000+ растров, добавлять в проект по 2 руками можно и устать.

Какой проект имеется в виду? Qgis?
Я сшиваю из командной строки ОС скриптом gdal_merge.py.
При сшивке 2х файлов между ними все равно остается "шов".

[Максим Дубинин]

может пару растров выложить, чтобы кто-то еще попробовал?

[Voltron]

Возможно, растры привязаны не по ценрту пикселя, а по левому верхнему углу

[ericsson]

Похожий эффект наблюдался у меня при использовании GDAL с минидрайвером WMS в качестве источника, при "больших" значениях размера тайла. Было ощущение, что где-то происходит какое-то округление.

[nextstopsun]

Возможно, растры привязаны не по ценрту пикселя, а по левому верхнему углу

А как это проверить?

может пару растров выложить, чтобы кто-то еще попробовал?

На самом деле это растры ASTER DEM на территорию РФ. Могу выложить в форум, но они довольно весомые.
Как поступить?

[Дмитрий Барышников]

ASTER DEM в ArcGIS шьются без "щелей"

[gimran]

ASTER DEM в ArcGIS шьются без "щелей"

Но вопрос-то о gdal_merge....

[Дмитрий Барышников]

Это я написал к тому, что дело точно не в ошибке данных. И привязано там все ОК.

[nextstopsun]

Сшил ArcGIS-ом пару растров - щели действительно не наблюдаю.....придется теперь питоновский скрипт сочинять, чтоб сшилось как мне надо. ArcGIS Model Builder структуру каталогов видимо не создает.
Не ожидал, что подведет gdal (

[Максим Дубинин]

сшил на пробу пару горизонтально и вертикально, щелей не вижу

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

gdal_merge.bat ASTGTM_N54E166_dem.tif ASTGTM_N54E167_dem.tif

[Voltron]

Хм... тоже прогнал пару тестов, нет щелей.
ArcGIS давно не видел, но разве он не использует тот же GDAL?

[nextstopsun]

Прогнал тесты с теми же растрами.

Понял, что проблема в другом, поэтому вопрос перефразирую:

Когда строю из исходных смежных растров модель затенения рельефа:
gdaldem hillshade ASTGTM_N54E166_dem.tif hillshade_ASTGTM_N54E166_dem.tif -z 5 -s 111120
и
gdaldem hillshade ASTGTM_N55E166_dem.tif hillshade_ASTGTM_N55E166_dem.tif -z 5 -s 111120

получаю "щели" (см. картинку).
То есть hillshade меньше исходного растра по краям.
Что можно сделать?

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

А если использовать r.shaded.relief + r.patch? Не зря же грасс висит на скриншоте Тем более, что gdaldem сделан на базе грассовских модулей (r.slope.aspect,r.shaded.relief,r.colors)
При желании можно накатать шелл-скрипт, который создавал бы растры теневой отмывки и сшивал их в нужном порядке. В GRASS я не сталкивался с "щелями" при сшивке.

А вообще, было бы недурно увидеть вывод команд:
gdalinfo ASTGTM_N54E166_dem.tif
gdalinfo hillshade_ASTGTM_N54E166_dem.tif