[fwtools 2.4.7] Créer un calque sur une orthophoto

[msylvain]

Bonjour,

Je possède une orthophoto sous la formes de dalles ECW (Lambert II étendu) et je dois placer un calque d'information dessus puis visualiser le résultat dans la dernière version de IGNMap (04092012).

J'essaye de créer un calque à partir d'une image PNG à partir du WKT de ces dalles. Je connais l'emplacement de ce calque dans le repère de l'orthophoto.

gdal_translate -of GTiff -a_srs NM.wkt -a_ullr 303466.09375 253533.796875 308741.3125 251253.90625 image_00_00.png image_00_00.tif

Avec ce fichier NM.wkt :

PROJCS["LM1FRE2D",
  GEOGCS["N.T.F.",
      DATUM["NTF",
          SPHEROID["CLA80IGN",6378249.2,293.4660213]],
      PRIMEM["Greenwich",0],
      UNIT["degree",0.0174532925199433]],
  PROJECTION["Lambert_Conformal_Conic_1SP"],
  PARAMETER["latitude_of_origin",46.8],
  PARAMETER["central_meridian",2.337229166666664],
  PARAMETER["scale_factor",0.99987742],
  PARAMETER["false_easting",600000],
  PARAMETER["false_northing",200000],
  UNIT["Meter",1]]

Pour information, gdalinfo me donne (similaire pour toutes les dalles):

Driver: ECW/ERMapper Compressed Wavelets
Files: NM_01.ecw
Size is 198500, 46000
Coordinate System is:
PROJCS["LM1FRE2D",
    GEOGCS["N.T.F.",
        DATUM["NTF",
            SPHEROID["CLA80IGN",6378249.2,293.4660213]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",2.337229166666664],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",200000],
    UNIT["Meter",1]]
Origin = (301899.999999999770000,269799.999999999830000)
Pixel Size = (0.100000000000000,-0.100000000000000)
Corner Coordinates:
Upper Left  (  301900.000,  269800.000) (  1d36'38.32"W, 47d21'38.02"N)
Lower Left  (  301900.000,  265200.000) (  1d36'27.32"W, 47d19'9.24"N)
Upper Right (  321750.000,  269800.000) (  1d20'53.38"W, 47d22'9.22"N)
Lower Right (  321750.000,  265200.000) (  1d20'43.11"W, 47d19'40.42"N)
Center      (  311825.000,  267500.000) (  1d28'40.57"W, 47d20'39.49"N)
Band 1 Block=198500x1 Type=Byte, ColorInterp=Red
  Overviews: arbitrary
Band 2 Block=198500x1 Type=Byte, ColorInterp=Green
  Overviews: arbitrary
Band 3 Block=198500x1 Type=Byte, ColorInterp=Blue
  Overviews: arbitrary

J'obtiens ces informations pour le calque :

Driver: GTiff/GeoTIFF
Files: image_00_00.tif
Size is 2048, 1004
Coordinate System is:
PROJCS["LM1FRE2D",
    GEOGCS["N.T.F.",
        DATUM["NTF",
            SPHEROID["CLA80IGN",6378249.2,293.4660213]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",2.337229166666664],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",200000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (303466.093750000000000,253533.796875000000000)
Pixel Size = (2.575790405273438,-2.270807395418327)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL


Corner Coordinates:
Upper Left  (  303466.094,  253533.797) (  1d34'45.17"W, 47d12'54.45"N)
Lower Left  (  303466.094,  251253.906) (  1d34'39.77"W, 47d11'40.71"N)
Upper Right (  308741.313,  253533.797) (  1d30'34.75"W, 47d13'2.88"N)
Lower Right (  308741.313,  251253.906) (  1d30'29.45"W, 47d11'49.14"N)
Center      (  306103.703,  252393.852) (  1d32'37.29"W, 47d12'21.81"N)

Band 1 Block=2048x1 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA 

Band 2 Block=2048x1 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA 

Band 3 Block=2048x1 Type=Byte, ColorInterp=Blue
  Mask Flags: PER_DATASET ALPHA 
Band 4 Block=2048x1 Type=Byte, ColorInterp=Alpha

Ensuite, j'ouvre mes dalles et mon calque dans IGNmap, mais le calque se retrouve placer autour de 0.0 0.0.

Je suis preneur de toutes les idées qui vous passent par la tête.

Merci d'avance