banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

Pages: 1

 

#1 Thu 31 January 2013 10:31

aquila96
Participant occasionnel
Lieu: Turku
Date d'inscription: 11 Dec 2009
Messages: 21

Problème de projection: décalage images Landsat5 track 190 et 191

Bonjour,

J'ai récemment téléchargé une série temporelle d'images Landsat 5 d'une région finlandaise à partir du site L'Agence Spatiale Européenne.
Lesdites images sont normalement en UTM34N (epsg 32634) ce que me confirme GDAL. Or je me suis rendu comtpe d'un petit souci en les traitant: celles provenant du Track 191 de Landsat sont à leur place, ce que mon confirme le shape de ma zone d'étude (en rouge à droite sur l'image). Par contre, toutes les images provenant du track 190 sont environ 300 km à l'ouest de leur position réelle et je n'arrive pas à comprendre d'où vient l'erreur (si c'est moi qui en fait une).

Voici les sorties de GDAL pour une image du track 191

Code:

C:\Users\alevil>gdalinfo E:/Boulot/Spatial/Landsat/Kauhava/ESA/LS05_TM__SYC_1P_1
9860428T091217_19860428T091232_KSE_11474_0000.TIFF/L5191000_00019860428_B60.TIF
Driver: GTiff/GeoTIFF
Files: E:/Boulot/Spatial/Landsat/Kauhava/ESA/LS05_TM__SYC_1P_19860428T091217_198
60428T091232_KSE_11474_0000.TIFF/L5191000_00019860428_B60.TIF
Size is 6920, 5960
Coordinate System is:
PROJCS["WGS 84 / UTM zone 34N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",21],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32634"]]
GeoTransform =
  505814.9545012147, 28.71599757058706, -8.682826931710533
  7125649.357998786, -8.682826931710533, -28.71599757058706
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  505814.955, 7125649.358) ( 21d 7'11.96"E, 64d15'24.91"N)
Lower Left  (  454065.306, 6954502.012) ( 20d 6'6.37"E, 62d43'5.18"N)
Upper Right (  704529.658, 7065564.196) ( 25d 8'0.18"E, 63d39'30.21"N)
Lower Right (  652780.009, 6894416.850) ( 23d55'55.62"E, 62d 9'2.22"N)
Center      (  579297.482, 7010033.104) ( 22d34'37.23"E, 63d12'38.32"N)
Band 1 Block=6920x64 Type=Byte, ColorInterp=Gray

et pour une provenant du track 190

Code:

C:\Users\alevil>gdalinfo  E:/Boulot/Spatial/Landsat/Kauhava/ESA/LS05_TM__SYC_1P_
19880528T091130_19880528T091145_KSE_22556_0000.TIFF/L5190000_00019880528_B60.TIF

Driver: GTiff/GeoTIFF
Files: E:/Boulot/Spatial/Landsat/Kauhava/ESA/LS05_TM__SYC_1P_19880528T091130_198
80528T091145_KSE_22556_0000.TIFF/L5190000_00019880528_B60.TIF
Size is 3500, 3044
Coordinate System is:
PROJCS["WGS 84 / UTM zone 34N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",21],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32634"]]
GeoTransform =
  289183.6034257925, 28.04314841492462, -10.65747751480202
  7063164.021574208, -10.65747751480202, -28.04314841492462
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  289183.603, 7063164.022) ( 16d44'35.54"E, 63d37'59.56"N)
Lower Left  (  256742.242, 6977800.678) ( 16d13'8.21"E, 62d50'56.30"N)
Upper Right (  387334.623, 7025862.850) ( 18d44'55.94"E, 63d20'37.33"N)
Lower Right (  354893.261, 6940499.506) ( 18d10'34.84"E, 62d34'0.49"N)
Center      (  322038.432, 7001831.764) ( 17d28'22.83"E, 63d 6'6.91"N)
Band 1 Block=3500x64 Type=Byte, ColorInterp=Gray

J'ai contacté l'ESA, parce que j'ai un doute sur le fait que ça vienne de mon côté (ils ont apparemment des soucis avec le logiciel qui traite leurs images) mais en attendant, si l'un(e) d'entre vous a une suggestion sur l'origine de ce décalage, et surtout comment le corriger, je suis preneur. Et puis étant loin d'être infaillible, je serais content d'avoir un avis extérieur.

Par avance merci

Alex

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.
Alexandre Villers, Postdoc researcher
Spatial ecologly and population dynamics
Section Ecology, University of Turku, 20014 Turku

Hors ligne

 

 

Pages: 1

Pied de page des forums

Powered by FluxBB