[Gfoss] gdal_merge: come si usa?

[Paolo Craveri]

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
-------------- parte successiva --------------
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
-------------- parte successiva --------------
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