|

»

gdal_merge: come si usa?

View: New views

4 Messages — Rating Filter: Alert me

	gdal_merge: come si usa? by pcraveri :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message Buongiorno a tutti ho un gran numero di raster che devo mosaicare 4 a 4; i raster hanno le caratteristiche contenute nel file in allegato (gdalinfo_input.txt). Ho effettuato alcune prove su due di questi files ma non riesco ad ottenere il risultato voluto: gdal_merge.py -v -o merge.tif -co "COMPRESS=CCITTFAX4" -co NBITS=1 S191090lzw.tif S191100lzw.tif il risultato che desidero è un raster ad una sola banda (0=bianco o, meglio, trasparente; 1=nero) come (credo!) i files di input. Il risultato appare codificato al contrario (0=nero, 1=biano) ed inoltre la zona di adiacenza tra i due raster di input appare completamente nera (gdalinfo_output.txt). Uso GDAL 1.6.2, released 2009/07/31, da pacchetto per ubuntu. Dove sbaglio? ciao a tutti grazie -- -- Paolo C. Lat. 44° 39' 11.08'' N Long. 7° 23' 25.26'' E Driver: GTiff/GeoTIFF Files: S191090lzw.tif Size is 8760, 5326 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.2572235629972, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (7.330839642506918,44.700193642631803) Pixel Size = (0.000009832827403,-0.000009832827403) Metadata: AREA_OR_POINT=Area Image Structure Metadata: COMPRESSION=CCITTFAX4 INTERLEAVE=BAND Corner Coordinates: Upper Left ( 7.3308396, 44.7001936) ( 7d19'51.02"E, 44d42'0.70"N) Lower Left ( 7.3308396, 44.6478240) ( 7d19'51.02"E, 44d38'52.17"N) Upper Right ( 7.4169752, 44.7001936) ( 7d25'1.11"E, 44d42'0.70"N) Lower Right ( 7.4169752, 44.6478240) ( 7d25'1.11"E, 44d38'52.17"N) Center ( 7.3739074, 44.6740088) ( 7d22'26.07"E, 44d40'26.43"N) Band 1 Block=8760x7 Type=Byte, ColorInterp=Palette Image Structure Metadata: NBITS=1 Color Table (RGB with 2 entries) 0: 255,255,255,255 1: 0,0,0,255 Driver: GTiff/GeoTIFF Files: merge.tif Size is 17226, 5326 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.2572235629972, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (7.330839642506918,44.700193642631803) Pixel Size = (0.000009832827403,-0.000009832827403) Metadata: AREA_OR_POINT=Area Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=BAND Corner Coordinates: Upper Left ( 7.3308396, 44.7001936) ( 7d19'51.02"E, 44d42'0.70"N) Lower Left ( 7.3308396, 44.6478240) ( 7d19'51.02"E, 44d38'52.17"N) Upper Right ( 7.5002199, 44.7001936) ( 7d30'0.79"E, 44d42'0.70"N) Lower Right ( 7.5002199, 44.6478240) ( 7d30'0.79"E, 44d38'52.17"N) Center ( 7.4155298, 44.6740088) ( 7d24'55.91"E, 44d40'26.43"N) Band 1 Block=17226x3 Type=Byte, ColorInterp=Palette Image Structure Metadata: NBITS=1 Color Table (RGB with 2 entries) 0: 0,0,0,255 1: 255,255,255,255 _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... http://lists.faunalia.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.
	Re: gdal_merge: come si usa? by Antonio Falciano :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message Paolo Craveri ha scritto: > Buongiorno a tutti > > ho un gran numero di raster che devo mosaicare 4 a 4; i raster hanno > le caratteristiche contenute nel file in allegato > (gdalinfo_input.txt). > > Ho effettuato alcune prove su due di questi files ma non riesco ad > ottenere il risultato voluto: > > gdal_merge.py -v -o merge.tif -co "COMPRESS=CCITTFAX4" -co NBITS=1 > S191090lzw.tif S191100lzw.tif > > il risultato che desidero Ãš un raster ad una sola banda (0=bianco o, > meglio, trasparente; 1=nero) come (credo!) i files di input. > > Il risultato appare codificato al contrario (0=nero, 1=biano) ed > inoltre la zona di adiacenza tra i due raster di input appare > completamente nera (gdalinfo_output.txt). > > Uso GDAL 1.6.2, released 2009/07/31, da pacchetto per ubuntu. > > Dove sbaglio? Prova ad impostare il valore di nodata con "-n 0" come in: http://trac.osgeo.org/gdal/wiki/UserDocs/GdalMerge#Specifyoverlapprecedence ciao Antonio -- Antonio Falciano http://www.linkedin.com/in/antoniofalciano _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... http://lists.faunalia.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.
	Re: gdal_merge: come si usa? by pcraveri :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message ciao Il 03 ottobre 2009 17.26, Antonio Falciano <afalciano@...> ha scritto: > Paolo Craveri ha scritto: >> Buongiorno a tutti >> >> ho un gran numero di raster che devo mosaicare 4 a 4; i raster hanno >> le caratteristiche contenute nel file in allegato >> (gdalinfo_input.txt). >> >> Ho effettuato alcune prove su due di questi files ma non riesco ad >> ottenere il risultato voluto: >> >> gdal_merge.py -v -o merge.tif -co "COMPRESS=CCITTFAX4" -co NBITS=1 >> S191090lzw.tif S191100lzw.tif >> >> il risultato che desidero Ãš un raster ad una sola banda (0=bianco o, >> meglio, trasparente; 1=nero) come (credo!) i files di input. >> >> Il risultato appare codificato al contrario (0=nero, 1=biano) ed >> inoltre la zona di adiacenza tra i due raster di input appare >> completamente nera (gdalinfo_output.txt). >> >> Uso GDAL 1.6.2, released 2009/07/31, da pacchetto per ubuntu. >> >> Dove sbaglio? > > Prova ad impostare il valore di nodata con "-n 0" come in: > http://trac.osgeo.org/gdal/wiki/UserDocs/GdalMerge#Specifyoverlapprecedence > > ciao > Antonio > > -- > Antonio Falciano > http://www.linkedin.com/in/antoniofalciano > grazie per la risposta: gdal_merge.py -v -n 0 -o merge.tif -co "COMPRESS=CCITTFAX4" -co NBITS=1 S191090lzw.tif S191100lzw.tif risolve il problema della zona di adiacenza tra i due raster che, con '-o 0' risulta ora visibile. Seguento il link che mi hai indicato ho trovato spunto per questo: gdalwarp -srcnodata 0 -dstnodata 0 -r cubic S191090lzw.tif S191100lzw.tif mergewarp.tif che fornisce il risultato atteso, fatta salva la dimensione del file di output (non compresso). Mi sembra che per risolvere quest'ultimo problema dovrei fa seguire gdal_translate per ottenere un .tif compresso (-co "COMPRESS=CCITTFAX4"). Mi chiedo se è possibile fare tutto in un'unica operazione. grazie, ciao -- -- Paolo C. Lat. 44° 39' 11.08'' N Long. 7° 23' 25.26'' E _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... http://lists.faunalia.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.
	Re: gdal_merge: come si usa? by Antonio Falciano :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message Paolo Craveri ha scritto: > gdal_merge.py -v -n 0 -o merge.tif -co "COMPRESS=CCITTFAX4" -co > NBITS=1 S191090lzw.tif S191100lzw.tif > > risolve il problema della zona di adiacenza tra i due raster che, con > '-o 0' risulta ora visibile. > > Seguento il link che mi hai indicato ho trovato spunto per questo: > > gdalwarp -srcnodata 0 -dstnodata 0 -r cubic S191090lzw.tif > S191100lzw.tif mergewarp.tif > > che fornisce il risultato atteso, fatta salva la dimensione del file > di output (non compresso). Mi sembra che per risolvere quest'ultimo > problema dovrei fa seguire gdal_translate per ottenere un .tif > compresso (-co "COMPRESS=CCITTFAX4"). > > Mi chiedo se è possibile fare tutto in un'unica operazione. Credo di si... prova ad eliminare le virgolette dopo il primo "-co": gdal_merge.py -v -n 0 -o merge.tif -co COMPRESS=CCITTFAX4 -co NBITS=1 S191090lzw.tif S191100lzw.tif dovendo essere -co NAME=VALUE come da: http://www.gdal.org/gdal_merge.html buon weekend Antonio -- Antonio Falciano http://www.linkedin.com/in/antoniofalciano _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... http://lists.faunalia.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.