Suite

La fusion de GeoTiffs avec GDAL provoque des artefacts de type grille ?


Je traite les tuiles ASTER DEM, les fusionne en blocs plus grands (pour minimiser les bords visibles des tuiles) et importe les contours dans une base de données postgis. Flux de travail :

gdal_merge.py -of GTiff -o merged.tif tile1.tif tile2.tif gdal_contour -i 10 -a ele -snodata -9999 merged.tif contours.shp ogr2ogr -f PostgreSQL -nln srtm4 -append -t_srs 'EPSG:3785' contours.shp

Pour une raison quelconque, les contours "fusionnés" résultants ont des artefacts en forme de grille :

Lorsque je saute l'étape de fusion et que je traite les deux mêmes GeoTiffs séparément, il n'y a pas d'artefacts. j'ai essayé de fusionner avecgdalbuildvrtet a obtenu le même effet de grille.

Pourquoi cela se produit-il et comment puis-je éviter cela ?


Il s'avère que c'est un problème de projection. Avant de fusionner, je reprojette les tuiles de EPSG:4326 à EPSG:3785. Lorsque cela est fait après la fusion, aucun problème ne se produit.


Voir la vidéo: Convert between CSV and GeoTIFF with GDAL in Python (Octobre 2021).