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 Wed 26 September 2012 13:15

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

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.

Code:

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 :

Code:

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):

Code:

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 :

Code:

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

Hors ligne

 

#2 Wed 26 September 2012 14:28

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

Pourquoi n'utilisez-vous pas le SRS défini par l'IGN ?
-a_srs IGNF:LAMB1

A propos, FwTools est une distribution obsolète. Utililsez plutôt une version récente de GDAL.

Hors ligne

 

#3 Wed 26 September 2012 14:50

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

IGNF:LAMB1 est un équivalent de mon fichier WKT ?

Les versions récentes ont le même problème mais ne permette pas d'utiliser toutes les définitions données en exemple sur divers forums (par exemple "+init=IGNF:LAMBE +wktext").

Hors ligne

 

#4 Wed 26 September 2012 15:10

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

C'est l'inverse : le registre IGNF n'existait pas dans les versions anciennes de GDAL.

Hors ligne

 

#5 Wed 26 September 2012 15:26

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

Avec la dernière version de GDAL et le SRS que j'ai donné en exemple :

Code:

gdal_translate --version
GDAL 1.9.1, released 2012/05/16
gdal_translate -of GTiff -a_srs "+init=IGNF:LAMBE +wktext" -a_ullr 303944.028442 253985.271362 308446.838013 251777.839355 image_00_00.png image_00_00.tif
Failed to process SRS definition: +init=IGNF:LAMBE +wktext

Dernière modification par msylvain (Wed 26 September 2012 15:28)

Hors ligne

 

#6 Wed 26 September 2012 18:06

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

Aucun problème chez moi. C'est certainement un problème d'installation. Vérifiez si la variable d'environnement GDAL_DATA pointe sur le bon répertoire.

Hors ligne

 

#7 Wed 26 September 2012 18:13

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

Code:

where gdal_translate
C:\Program Files (x86)\GDAL\gdal_translate.exe
echo %GDAL_DATA%
C:\Program Files (x86)\GDAL\gdal-data

Hors ligne

 

#8 Wed 26 September 2012 23:23

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

En fait, c'est PROJ.4 qui fournit le fichier IGNF. Vérifiez si vous avez la dernière version de PROJ.4, et si la variable d'environnemet PROJ_LIB est correcte.

Hors ligne

 

#9 Thu 27 September 2012 12:55

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

Ajouter la variable PROJ_LIB a corrigé le précèdent problème pour ma version à jour de GDAL, merci.

J'ai donc relancé la procédure mais le problème reste identique. Le calque se trouve autour de 0.0, 0.0. Si j'utilise le SRS que vous proposez, le problème persiste.

Hors ligne

 

#10 Thu 27 September 2012 13:07

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

Je ne suis pas certain d'avoir bien compris pourquoi vous voulez calibrer votre image en Lambert II , alors que vos dalles ECW sont en Lambert II étendu.
N'est-ce pas là la source de vos problèmes ?

Dernière modification par jcr83 (Thu 27 September 2012 14:38)

Hors ligne

 

#11 Thu 27 September 2012 13:54

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

Les dalles sont dans ce repère :

Code:

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]]

Et j'utilise le même repère pour placer mon calque.

Je pensais que cette projection était du Lambert II étendu.

Hors ligne

 

#12 Thu 27 September 2012 14:38

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

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

Non, c'est du Lambert II Centre (IGNF:LAMB2, ou EPSG:27562). Le Lambert II étendu a pour codes IGNF:LAMBE, ou ESPG:27572. Il y a 2000 km de différence entre les deux sur l'axe des Y.

Hors ligne

 

#13 Thu 27 September 2012 16:06

msylvain
Juste Inscrit !
Date d'inscription: 26 Sep 2012
Messages: 7

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

Histoire de lever un doute, j'ai placer mon calque en utilisant les coordonnées de la dalle elle même.

Code:

gdal_translate -of GTiff -a_srs NM.wkt -a_ullr 301900.000 269800.000 321750.000 265200.000 image_00_00.png image_00_00.tif

Le résultat est le même. Je commence à croire que IGNMap ne fonctionne pas.

Hors ligne

 

 

Pages: 1

Pied de page des forums

Powered by FluxBB