Le portail francophone de la géomatique

benoitlarcher · Wed 03 June 2009 19:16

Bonjour à tous,

Dans la série "j'ai du mal avec les projections", j'ai récupéré hier un fichier de limites d'une aire protégée (polygone, je pense... à moins que ce ne soit une polyligne. Comment peut-on en être sûr? le vérifier?) venant de ArcView, qui était originellement en projection Lambert costarricaine (référencée par rapport à je en sais quel point au Honduras. Connaissez-vous ce système de projection apparemment valable pour toute l'Amérique Centrale?). Le fichier a été transformé grâce une extension de ArcView en un fichier projeté en UTM 1983 (si j'ai bien suivi l'affaire), puis m'a été transmis comme tel en fichier shape.

Quand je l'ouvre avec mapinfo en shp, le polygone s'affiche normalement (pièce jointe 1), mais quand j'ouvre en même temps la table contenant l'image satellite de la zone (projetée elle en Long/lat WGS84), le polygone disparaît. J'enregistre donc une copie de la table contenant le polygone en le projetant en Long/Lat WGS84 pour que les 2 se superposent (l'image sat et le polygone), mais là, tout disparait de nouveau. Je fais donc "afficher toutes les couches" dans le menu qui s'ouvre en clic-droit, et là, je ne vois plus rien, car si l'image sat est elle bien placée, le polygone (qui l'était (bien placé) lorsqu'il était encore sous ArcView) ne l'est plus, et de plus, il a subi à se moment la une symétrie et une rotation (cf pièce jointe 2, vue après avoir fait "afficher la couche Carpintera", celle contenant ce foutu polygone ;-)). En gros, il faudrait je pense recaler ce polygone (comme si son géoréférencement avait disparu lors de sa conversion en *.shp, ou en *.tab...???) dans les bonnes coordonnées, mais comment faire cela? ou pensez-vous à une autre solution?

Merci beaucoup pour votre aide, et désolé de ne pas être plus clair dans mes explications...
Bonne soirée/ journée

Benoît

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

benoitlarcher · Wed 03 June 2009 19:21

et voici le deuxième fichier joint que je n'ai pas réussi à attacher au 1er post...

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Daniel · Wed 03 June 2009 20:49

Bonjour Benoît
C'est le gars du service de soir!!!!
En théorie, les objets dans Mapinfo "devraient" se superposer, lorsqu'il y a deux projections différentes, à la volée. Par exemple, une couche de polygones en UTM 18 devrait se superposée par dessus une image en Lat/Long sans problème. Le logiciel est en mesure de comprendre celà. Cependant, il peut y avoir certain problème. Selon mon expérience, la majorité des problèmes avec des images, c'est que celles-ci ne sont pas bien géoréférencées.
À ta place, je m'assurerais que l'image est bien dans le système que tu as mentionné.
Pour ton système de projection costarricaine, le seule moyen de asvoir est d'ouvrir dans Bloc Note le fichier MAPINFOW.PRJ et avec l'aide de Mapinfo, trouver la projection en regardant les paramètres. Ça ne devrait pas être trop long car les projections sont placées par groupe.
Ton fichier shape de départ, lorsque tu l'importes dans Mapinfo, utilises-tu le traducteur universelle? Si oui, enlève le crochet pour lire la projection dans le ficheir et force la projection en UTM Bande???, NAD 83.

Reviens-nous là-dessus.

Daniel

benoitlarcher · Wed 03 June 2009 23:55

Bien l'bonsoir Daniel !

En fait, quand j'ouvre le fichier d'origine ArcView (donc au format shape) directement dans MI en faisant "Fichier, ouvrir" et en choisissant "Fichier Shape .shp", le code que cela me donne est :

!table
!version 700
!charset WindowsLatin1

Definition Table
Type SHAPEFILE Charset "WindowsLatin1"
Fields 6
Id Decimal (11, 0) ;
Area Decimal (20, 3) ;
Perimeter Decimal (20, 3) ;
Hectares Decimal (20, 3) ;
Nombre Char (36) ;
Código Char (16) ;
begin_metadata
"\IsReadOnly" = "FALSE"
"\Shapefile" = ""
"\Shapefile\PersistentCache" = "FALSE"
"\Spatial Reference" = ""
"\Spatial Reference\Geographic" = ""
"\Spatial Reference\Geographic\Projection" = ""
"\Spatial Reference\Geographic\Projection\Clause" = "CoordSys Earth Projection 8, 33, ""m"", -93, 0, 0.9996, 500000, 0 Bounds (-7745844.29605, -9997964.94315) (8745844.29605, 9997964.94315)"
"\DefaultStyles" = ""
"\DefaultStyles\Pen" = ""
"\DefaultStyles\Pen\LineWidth" = "1"
"\DefaultStyles\Pen\LineStyle" = "0"
"\DefaultStyles\Pen\Color" = "0"
"\DefaultStyles\Pen\Pattern" = "2"
"\DefaultStyles\Brush" = ""
"\DefaultStyles\Brush\Pattern" = "2"
"\DefaultStyles\Brush\Forecolor" = "16777215"
"\DefaultStyles\Brush\Backcolor" = "16777215"
end_metadata

La ligne de référencement spatial ne correspond à rien dans le fichier MAPINFOW.PRJ... Bizarre, non?
Après, si j'utilise le traducteur universel, je n'ai que ça en ouvrat le .tab avec l'éditeur de texte :

!table
!version 300
!charset WindowsLatin1

Definition Table
Type NATIVE Charset "WindowsLatin1"
Fields 6
ID Decimal (11, 0) ;
AREA Decimal (20, 3) ;
PERIMETER Decimal (20, 3) ;
HECTARES Decimal (20, 3) ;
NOMBRE Char (36) ;
CóDIGO Char (16) ;

Après l'avoir ouvert de manière direct, j'ai essayé de la projeter dans 3 systèmes différents :
- Long/Lat WGS8A (qui est normalement le système de projection de mon image satellite, provenant de GoogleMaps)
- Long/Lat Nad83
- Long/Lat Nad27 for central america
A chaque fois, en ouvrant avec l'éditeur de texte,; j'obtient exactement la même chose, à l'exception du fait que pour la ligne Type SHAPEFILE Charset "WindowsLatin1", cela change en Type NATIVE Charset "WindowsLatin1"

J'avoue que je suis complètement perdu, mais si c'est inexplicable, je vais essayer de trouver une solution tout seul...
merci de ne pas vous casser la tête... je m'en voudrais dans ce cas!

Benoît

PS : pour info, voici le code de mon image, à la base un fichier .bmp, importé dans MI, calé...

!table
!version 300
!charset WindowsLatin1

Definition Table
File "imsatzpcc.bmp"
Type "RASTER"
(9.884969,-83.936127) (2396,567) Label "Pista Cartago",
(9.900313,-84.002142) (855,183) Label "San Diego",
(9.883525,-84.035532) (84,585) Label "Patarra",
(9.871843,-83.993969) (1051,858) Label "Coris"
CoordSys Earth Projection 8, 33, "m", -93, 0, 0.9996, 500000, 0
Units "degree"

...puis projeté en Long/Lat WGS84 :

!table
!version 300
!charset WindowsLatin1

Definition Table
File "imsatzpcc.bmp"
Type "RASTER"
(9.884069,-83.936127) (2394,569) Label "Pista Cartago",
(9.900313,-84.002142) (856,181) Label "San Diego",
(9.883525,-84.03532) (84,585) Label "Patarra",
(9.871843,-83.993969) (1049,857) Label "Coris"
CoordSys Earth Projection 1, 104
Units "degree"

...puis en UTM NAD 83 Zone 15 (projection dans laquelle était le fichier shape au départ quand on me l'a donné...j'ai tenté de faire cette conversion pour voir si les deux fichiers (image sat et shpe) pouvaient se superposer de manière correcte) :

!table
!version 300
!charset WindowsLatin1

Definition Table
File "imsatzpcc.bmp"
Type "RASTER"
(9.884969,-83.936127) (2396,567) Label "Pista Cartago",
(9.900313,-84.002142) (855,183) Label "San Diego",
(9.883525,-84.035532) (84,585) Label "Patarra",
(9.871843,-83.993969) (1051,858) Label "Coris"
CoordSys Earth Projection 8, 33, "m", -93, 0, 0.9996, 500000, 0
Units "degree"

benoitlarcher · Thu 04 June 2009 00:07

Effectivement, quand j'ouvre le polygone contenu dans le fichier shape, puis l'image satellite ensuite, alors le polygone disparai, et/ou se déforme quand je fais "afficher l'ensemble des couches"... d'où un problème de géoréférencement dans l'image (que j'ai pourtant fait à X reprise...).
En fait, ce que je ne comprend pas, c'est s'il faut, quand on cale l'image avec des points donnés, choisir comme projection celle d'origine de l'image (qui doit être Long/lat WGS84 vu que l'image vient de GoogleMaps) ou celle que l'on souhaite donner à l'image (genre, UTM NAD 83 zone15 ou autre...) ???

Thanks again.

Benoît

Daniel · Thu 04 June 2009 13:47

Bonjour Benoit
Pour ton image, il faut toujours donner la projection d'origine. Pour changer la projection d'une image, il faut un autre logiciel que Mapinfo.

Daniel
P.S. Je regarde les détails de tes .TAB

NicolasGeorge · Thu 04 June 2009 14:50

benoitlarcher a écrit:

En fait, ce que je ne comprend pas, c'est s'il faut, quand on cale l'image avec des points donnés, choisir comme projection celle d'origine de l'image (qui doit être Long/lat WGS84 vu que l'image vient de GoogleMaps) ou celle que l'on souhaite donner à l'image (genre, UTM NAD 83 zone15 ou autre...) ???

Quand on fait le calage à la main, il vaut mieux utiliser une projection qui colle avec les axes X et Y de l'image, mais ce n'est pas non plus obligatoire (surtout depuis que MI8.5 sait reprojeter les images pour l'affichage).

Ici, les coordonnées des points de calage de l'image sont en exprimés en degrés de latitude/longitude (il s'agit bien du Costa Rica, et d'une zone entre San José et Cartago?), il faut donc rester dans une projection de type latitude/longitude, et surtout ne pas leur associer une projection UTM en mètres. Erreur possible N°1.

Après, pour voir où est l'erreur entre entre le shape (qui semble correct ouvert directement par MapInfo) et votre image de calage, utilisez une référence extérieure comme les données échantillons MapInfo (helas assez peu précises) ou une vue satellite en accès WMS (JPL ou autre).

benoitlarcher · Thu 18 June 2009 00:41

Bonsoir à tous,

erreur résolue cf http://georezo.net/forum/viewtopic.php? … 07#p141607
je m'étais stupidement gouré en inversant les X et les Y à cause du fait que l'on parle latitude/longitude (du moins, c'est ce qui me paraissait à la base normal!)
Merci à tous pour vos explications.
Bonne journée

Benoît

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Wed 03 June 2009 19:16

Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#2 Wed 03 June 2009 19:21

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#3 Wed 03 June 2009 20:49

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#4 Wed 03 June 2009 23:55

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#5 Thu 04 June 2009 00:07

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#6 Thu 04 June 2009 13:47

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#7 Thu 04 June 2009 14:50

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

#8 Thu 18 June 2009 00:41

Re: Pbl de transformation de projection UTM 1983 (?) vers Long/Lat WGS84

Pied de page des forums