Le portail francophone de la géomatique

Essaye de regarder du côté des options -srcnodata et -dstnodata de gdalwarp.

Maintenant, j'arrive à avoir un <NoDataValue> dans le fichier résultat de gdalwarp. Mais ça ne marche toujours pas...

Pour l'instant je fais :
- "-srcnodata 0" dans gdalbuildvrt
- "-dstnodata 0" dans gdalwarp
- "-a_nodata 0" dans gdal_translate

J'ai essayé aussi

- rien dans gdalbuildvrt
- "-srcnodata 0" dans gdalwarp
- "-a_nodata 0" dans gdal_translate

Et d'autres "combinaisons"...mais rien ne marche.

Ted,

Projeter chaque image ça n'est pas (dans mon cas) la meilleure solution car j'ai notamment un jeu de données où il y a 54 dalles.

Mais j'ai essayé pour un autre jeu de données où il n'y en a que deux.
Donc :
- j'ai reprojeté les deux dalles
- je les ai assemblées
- j'ai ajouté les paramètres dans le fichier vrt en sortie du gdalbuildvrt
- j'ai sorti le fichier ecw avec gdal_translate

Mais le résultat ne va pas, déjà le fond n'est toujours pas noir, et en plus les deux dalles ne se superposent pas.
Et s'il faut utiliser les options (srcnodata ...etc.) je ne sais pas comment faire pour que dans ce cas-là ça fonctionne.

Dernière modification par Sam_Dalembert (Wed 03 February 2010 14:17)

Bonjour,

Si votre format final est l'ECW, vous n'aurez jamais un contour noir net puisque la valeur peut varier de 0 à 3-4, il vous faudrait donc masquer non seulement  "0 0 0", mais ce qui est proche "0 1 1", "0 1 0"...
Par contre, en utilisant le jp2ecw, on arrive a avoir une valeur uniforme des bordures noires et alors on peut les masquer.

Mais peut-être n'ai je pas bien compris votre manip.

Bon je vais continuer là, pour peut-être pas créer un topic pour chaque question...


Concernant gdalinfo, il y a une petite subtilité que je ne comprends pas.
Exemple : je pars d'un fichier en L2, je reprojette virtuellement en CC48, et j'extrais (avec -projwin). Donc ça me met le fichier ecw, et un fichier xml.
Quand je fais un gdalinfo du fichier, ça me renseigne bien les infos de projection, mais les "Corner Coordinates" sont en "pixels".

PROJCS["RGF93 / CC48",
    GEOGCS["RGF93",
        DATUM["Reseau_Geodesique_Francais_1993",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6171"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4171"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",47.25],
    PARAMETER["standard_parallel_2",48.75],
    PARAMETER["latitude_of_origin",48],
    PARAMETER["central_meridian",3],
    PARAMETER["false_easting",1700000],
    PARAMETER["false_northing",7200000],
    AUTHORITY["EPSG","3948"],
    AXIS["X",EAST],
    AXIS["Y",NORTH]]
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,16173.0)
Upper Right (19860.0,    0.0)
Lower Right (19860.0,16173.0)
Center      ( 9930.0, 8086.5)

Quand j'enlève le fichier xml, il me met les bonnes coordonnées mais ne me donne plus les infos de projection ! Et en plus (plutôt une bonne chose), la taille du pixel apparaît alors qu'elle n'y est pas quand je laisse le fichier xml !

Coordinate System is `'
Origin = (1624451.667468465400000,7195024.645112288200000)
Pixel Size = (0.249962856981897,-0.249962856981897)
Corner Coordinates:
Upper Left  ( 1624451.667, 7195024.645)
Lower Left  ( 1624451.667, 7190981.996)
Upper Right ( 1629415.930, 7195024.645)
Lower Right ( 1629415.930, 7190981.996)
Center      ( 1626933.799, 7193003.320)

Je pensais que ça pouvait venir du gdal_translate, car je ne mets pas les options PROJ et DATUM, mais même en les mettant, ça me donne la même chose... Une idée ?

Dernière modification par Sam_Dalembert (Tue 09 February 2010 16:37)