Le portail francophone de la géomatique

J'ai regardé du coté de GDAL-OGR. Mais je ne trouve que des outils (comme gdaltindex ou gdal2tile) qui génèrent des tuiles dans un format shapefile ou kml à partir d'une série de rasters déjà existante. Or j'ai besoins de diviser mon raster en plusieurs morceaux avant de générer ces tuiles.

A quel outils précisément penses-tu si la librairie GDAL-OGR répond à mon besoin?

J'ai retrouvé dans l'arbo GDAL un script GDAL/python qui s'appelle gdal_retile.py qui fait exactement ce que tu veux. Tu peux le trouver là : http://trac.osgeo.org/gdal/browser/trun … format=raw

Par ex pour générer des tuiles de 1000x1000 dans le répertoire out_dir à partir du fichier in.tif et créer en même temps le tileindex tile.shp :
gdal_retile.py -ps 1000 1000 -tileIndex tile.shp -targetDir out_dir in.tif

Sinon, en plus rustique, tu peux regarder du côté de gdal_translate, en particulier le paramètre -srcwin qui permet d'extraire une partie du raster initial. Si tu veux découper en plein de morceaux, tu peux avoir besoin de scripter ces appels.

Sinon je m'interroge sur ton besoin intial. Quel est l'intérêt de découper un GeoTIFF pour ensuite en faire une mosaique avec gdaltindex ? Gagner en perfos ? Si ton GeoTIFF initial est tuilé, cela n'a aucun intérêt. S'il ne l'est pas, la première chose à faire est un gdal_translate -co TILED=YES in.tif out.tif et de réessayer.

Dernière modification par rouault (Wed 01 July 2009 21:08)

salut
j'ai un probleme au niveau de gdaltindex et je souhaite trouver quelqu'un qui pourrait me donnner un coup de main je m'explique j'ai exporté en .tif un dem via globalmapper le fichier tif ainsi obtenu contient plus de 200Mo de taille jusqu'à maintenant j'arrive a me debrouiller tout seul ( il s'agit d'un tif pour le maroc en wgs84) mon oups mes  problemes sont les suivants:
1- quand j'execute la commande gdalinfo maroc.tif il me signal ceci

    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["arc-second",4.848136811095357e-006],
    AUTHORITY["EPSG","4326"]]
Origin = (-68400.000000000000000,129600.000000000000000)
Pixel Size = (3.000000000000000,-3.000000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_XRESOLUTION=254
  TIFFTAG_YRESOLUTION=254
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  COMPRESSION=PACKBITS
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  -68400.000,  129600.000) (68400d 0'0.00"W,129600d 0'0.00"N)
Lower Left  (  -68400.000,   71997.000) (68400d 0'0.00"W,71997d 0'0.00"N)
Upper Right (       3.000,  129600.000) (  3d 0'0.00"E,129600d 0'0.00"N)
Lower Right (       3.000,   71997.000) (  3d 0'0.00"E,71997d 0'0.00"N)
Center      (  -34198.500,  100798.500) (34198d30'0.00"W,100798d30'0.00"N)

les coordonnée de l'extent me parraissent anormale (68400 degré bizar ou je comprend pas )y a t il des conseils ou des remèdes

2- mon 2eme probleme est que lorsque j'execute cette commande gdaltindex maroc.shp maroc.tif dans la racine ou se trouve mon maroc.tif il me genere cette erreur " erreur 4 : unable to open .shap or .SHP skipping mais quand meme il me crait un .shp et jarrive a visualiser mon raster sur cartoweb mais le probleme de coordonnées persiste . des coordonées anormales et pire encore  lorsque je fais des calculs de distances ou surface cartoweb me raconte des bétises


j'ai publié de signaler que je developpe sous windows xp sp3

veuillez accepter mon niveau tres bas dans la matière mais s'il y a des suggestions ou solutions je suis preneurs
merci

Dernière modification par josig (Sat 25 July 2009 12:32)

Bonjour,

J'étais abonné à ce sujet depuis le début.
Il se trouve que je me suis posé le même type de question et j'ai terminé justement hier.
C'était pour mes besoins propres, donc j'ai décidé après m'être consulté longuement que je divisais le raster en 16 parties, soit 4 x 4. Mais on peut le rendre paramétrable.

Dans ma version, je ne crée pas de fichier HDR et TFW, mais ce ne serait pas un problème.
Par contre, il y a un problème non négligeable et que je n'ai pas encore résolu, si il s'agit d'un fichier correspondant à une zone géographique importante (centaines de Km), la jonction entre les différents petits raster se fait mal. Il est bien évident que la solution à ce problème n'existe pas, mais ça n'empêche par de trouver une méthode pas trop bête 

A suivre.
Cordialement