2012. január 25., szerda

Nodata eltüntetése a szélekről

Amit az arcgis automatikusan felismer, ugyanaz open source szofverekkel kicsit döcögős.
Nodata érték az alábbi űrfelvételen nem lett beállítva, de valószínű, hogy a 255 255 255 RGB érték az. 

for x in *.tif;do gdal_translate -a_nodata 255 $x "${x/.tif}".tiff;done

 Terminálban egy for ciklussal a gdal_translate paranccsal végigfutok a tifeken, beállítom a 255 (255 255) rgb értéket nodatara (Feltehetőleg a felhő tartalmazhat 255 255 255 rgb pixeleket még).
Ezután virtuális rasztert készítek:

 gdalbuildvrt -hidenodata crto_hamis.vrt *.tiff

A hidenodata hatására a nodata nagy része átlátszóvá válik - persze nem mind :)
Ezután gdaltindex és mehet a map fájlba.

gdaltindex crto_hamis.shp crto_hamis.vrt

És az eredmény: 


Nincsenek megjegyzések:

Megjegyzés küldése