Suite

La mosaïque de geotiffs avec gdal crée des lacunes dues aux couches alpha


J'ai une série de géotiffs que j'essaie de mosaïquer ensemble, mais lorsqu'ils sont combinés, leurs couches alpha individuelles chevauchent des sections des autres images et créent des lacunes dans la mosaïque finale :

Mon approche est d'abord de générer un .vrt à partir de tous les geotiffs :

gdalbuildvrt -hidenodata mosaic.vrt geo_pict20140910_131*

Et puis pour convertir ce .vrt en geotiff avec gdal_translate :

gdal_translate mosaïque.vrt mosaïque.tif

Existe-t-il un moyen de mosaïquer uniquement les parties non alpha de ces géotiffs pour éviter ces lacunes ?

Voici un lien vers les geotiffs originaux : https://www.dropbox.com/s/udjq9bter9fmgxz/gtiffs.zip?dl=0

Voici un lien vers les mosaïques .vrt et .tif finales : https://www.dropbox.com/s/oov9vatz69f31zg/mosaic.zip?dl=0

Voici le gdalinfo pour la mosaïque finale :

Pilote : GTiff/GeoTIFF Fichiers : mosaic.tif La taille est 4782, 2837 Le système de coordonnées est : GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG"," 7030"]], AUTORITE["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degré",0.0174532925199433], AUTORITE["EPSG","4326"]] Origine = (- 121.460884208999985,38.583432355500015) Taille du pixel = (0.000000591354474, -0.000000591354474) Métadonnées : AREA_OR_POINT=Métadonnées de la structure de l'image de zone : INTERLEAVE=PIXEL Coordonnées des coins : en haut à gauche (-121.4608842, 38.5834324) (121d27'39.18"W, 38d35') En bas Gauche (-121.4608842, 38.5817547) (121d27'39.18"W, 38d34'54.32"N) En haut à droite (-121.4580564, 38.5834324) (121d27'29.00"W, 38d35' 0.36"N) En bas à droite (-121.4580564, 38.5817547) ( 121d27'29.00"W, 38d34'54.32"N) Centre (-121.4594703, 38.5825935) (121d27'34.09"W, 38d34'57.34"N) Band 1 Block=4782x1 Type=Byte, ColorInterp=Red Mask Flags: PER_DATASET ALPHA Band 2 Block=4782x1 Type=Byte, ColorInterp=Green Mask Flags: PER_DATASET ALPHA Band 3 Block=4782x1 Type=Byte, ColorInterp=Blue Mask Flags : PER_DATASET ALPHA Band 4 Block=4782x1 Type=Byte, ColorInterp=Alpha

J'ai trouvé une solution rapide - en remplaçant "-hidenodata" par "-srcnodata 0" dans la version .vrt :

gdalbuildvrt -srcnodata 0 mosaic.vrt geo_pict20140910_131*