#1 Tue 20 September 2016 07:50
- Thy06
- Juste Inscrit !
- Date d'inscription: 19 Sep 2016
- Messages: 4
entête données GRID ASCII
Bonsoir,
J'ai récupéré les données BD alti (pas 75 m) du site de l'IGN ( http://professionnels.ign.fr/bdalti )
Je ne suis pas spécialiste de la géomatique. Je souhaite utiliser ces données pour simuler un relief, une montagne... en arrière plan des perspectives que je réalise en 3D pour l'architecture. Cela permet de placer le projet dans son environnement.
Les données des fichiers .asc me semblent simples à comprendre, excepté les valeurs xllcorner et yllcorner.
J'ai cherché partout, je n'ai toujours pas compris à quoi ça correspond !
Je crois que ce sont des valeurs en mètres, mais à partir de quel repère ?
Ça aurait pu être le méridien de Greenwich et l'équateur, ou le pôle nord...
Mais j'ai testé un fichier, j'ai fait des mesures sur Google earth, ça ne correspond pas à ces repères.
Merci
Hors ligne
#2 Tue 20 September 2016 09:01
Re: entête données GRID ASCII
Bonjour,
Ce sont les coordonnées d'origine de votre fichier, un des coins (je ne sais plus lequel...).
XLLCENTER or XLLCORNER = X coordinate of the origin (by center or lower left corner of the cell).
YLLCENTER or YLLCORNER = Y coordinate of the origin (by center or lower left corner of the cell).
http://resources.esri.com/help/9.3/arcg … format.htm.
En effet, ce sont surement des données en mètres, leur origine (0;0) va dépendre de la projection utilisée.
Quelles sont ces valeurs?
Avez-vous d'autres fichiers de même nom et d'extension différentes? (.xml, .prj...)
Martin
Hors ligne
#3 Tue 20 September 2016 11:17
- fbecir
- Participant assidu
- Lieu: Saint-Mandé
- Date d'inscription: 16 Sep 2008
- Messages: 519
Re: entête données GRID ASCII
Bonjour
La BDALTI à 75 m est livrée dans la projection officielle de la France métropolitaine, c'est-à-dire le Lambert 93 (c'est indiqué sur le site).
Les XLLCORNER et YLLCORNER correspondent aux coordonnées du coin Sud-Ouest de l'emprise couverte par le MNT.
Cordialement
Hors ligne
#4 Tue 20 September 2016 11:34
- Thy06
- Juste Inscrit !
- Date d'inscription: 19 Sep 2016
- Messages: 4
Re: entête données GRID ASCII
Bonjour, merci pour votre réponse.
@nitram
Ci-dessous l'entête d'un fichier (format .asc) correspondant à ma région (sud est), parmi tous les fichiers téléchargés (toute la France) :
ncols 1000
nrows 1000
xllcorner 974962.500000000000
yllcorner 6225037.500000000000
cellsize 75.000000000000
------
J'ai recherché parmi les fichiers joints, j'ai trouvé un fichier "dalles.prj"
Contenu (je l'ai reformaté avec retours lignes) :
PROJCS["RGF93_Lambert_93"
GEOGCS["GCS_RGF_1993"
DATUM["D_RGF_1993"
SPHEROID[« GRS_1980 »,6378137.0,298.257222101]]
PRIMEM[« Greenwich",0.0]
UNIT[« Degree",0.0174532925199433]]
PROJECTION["Lambert_Conformal_Conic"]
PARAMETER[« False_Easting",700000.0]
PARAMETER[« False_Northing »,6600000.0]
PARAMETER[« Central_Meridian",3.0]
PARAMETER[« Standard_Parallel_1",44.0]
PARAMETER[« Standard_Parallel_2",49.0]
PARAMETER[« Latitude_Of_Origin",46.5]
UNIT[« Meter",1.0]]
------
J'ai aussi extrait des informations des différents autres fichiers (pdf, html) joints. voir image jointe (infosBDalti01.jpg)
----
Je joins une 2ème image, correspondant à la zone couverte (France entière et même au-delà), décrite (sauf erreur de ma part) selon "l'étendue géographique" (données présentes tout en bas de la 1ère image).
J'ai supposé, que les valeurs xllcorner/yllcorner (tout en haut) sont relatives à cette zone. En x (975 km), ça pourrait fonctionner. En y (6225 km) ça n'a rien à voir.
... donc je patauge...
@fbecir
d'accord, mais ces données sont relatives à quel point sur la planète ?
dans l'exemple tout en haut, 975 km et 6225 km par rapport à où ?
Dernière modification par Thy06 (Tue 20 September 2016 11:40)
Hors ligne
#5 Tue 20 September 2016 22:04
- Jérémie N.
- Juste Inscrit !
- Date d'inscription: 10 Apr 2014
- Messages: 3
Re: entête données GRID ASCII
Bonjour
Votre fichier grid ascii a pour emprise un carré de 1000*75 m (ncols ou nrows * cellsize) et donc le coin bas gauche (sud ouest) se situe aux coordonées indiqué par xllcorner et yllcorner (xll ou yll pour x lower left ou y lower left)
Les jpgs dans votre message représente l'emprise géographique de la BD ALTI et non pas uniquement celle de votre fichier, il est donc normal que l'emprise soit un petit plus étendu.
Ensuite tout est dit dans votre fichier dalles.prj.
Ces coordonnées sont relatives à un point situé à l'intersection d'un méridien et d'un parallèle.
Ils sont indiqué pour le méridien par la ligne : PARAMETER[« Central_Meridian",3.0] donc soit 3° à l'est du méridien d'origine qui lui est indiqué par la ligne PRIMEM[« Greenwich",0.0] soit le méridien de Greenwich.
Pour le parallèle il faut regarder la ligne PARAMETER[« Latitude_Of_Origin",46.5], qui est donc un parallèle, le 46.5° nord.
Donc en résumé l'origine se situe au croisement du méridien 3°E et du parallèle 46.5°N soit quelque part dans l'Allier et comme toute origine de repère ce point aurait du avoir comme coordonnées planes (0,0).
Sauf que pour éviter de trainer des coordonnées négatives ce qui aurait vite été pénible, cette origine a été translatée de 700000m vers l'est (indiqué par le paramètre PARAMETER[« False_Easting",700000.0] ) et de 6600000m vers le nord (PARAMETER[« False_Northing »,6600000.0])
Le (0,0) du Lambert 93 se situe maintenant quelque part dans l'atlantique sud mais ce n'est pas grave et ça n'a aucune importance.
J'espère avoir été clair.
Cordialement.
Hors ligne
#6 Wed 21 September 2016 01:12
- Thy06
- Juste Inscrit !
- Date d'inscription: 19 Sep 2016
- Messages: 4
Re: entête données GRID ASCII
@ Jérémie N.
Bonsoir,
Je crois que vous avez la réponse.
J'ai fait des mesures sur Google Earth, en allant d'abord vers l'est et le nord à partir du point d'origine (0,0 centre de la France, (3.0 , 46.5)).
J'ai trouvé ça logique puisque vous dites translation de 700 km vers l'est et 6600 vers le nord. Mais je me suis retrouvé au dessus du pôle nord vers le détroit de Bering. Et en fait c'est à l'opposé (atlantique sud) qu'il fallait aller, comme vous l'indiquez aussi. Et là, mes valeurs xllcorner et yllcorner semblent correspondre.
Je vais retester par calculs dès que j'aurai le temps.
Je trouve ce repère dans l'Atlantique un peu bizarre. Mais peut-être qu'il sert aussi pour tous les autres territoires français d'outre-mer, d'où son emplacement ?
Merci grandement
Cordialement
@ autres contributeurs, merci aussi
Hors ligne
#7 Wed 21 September 2016 09:30
- fbecir
- Participant assidu
- Lieu: Saint-Mandé
- Date d'inscription: 16 Sep 2008
- Messages: 519
Re: entête données GRID ASCII
Bonjour Thy06
Je crois qu'il faudrait avant toute chose que vous lisiez cet article :
https://fr.wikipedia.org/wiki/Projection_cartographique
Si vos données ne sont pas dans la même projection, vous ne pouvez pas les superposer. Une simple translation ou une transformation un peu plus subtile (affinité, homographie) ne vous tirera pas d'affaire non plus ...
Les coordonnées (0 ; 0) dans un système de projection comme le Lambert 93 ont été fixées non pas sur un point remarquable mais simplement
définies de telle manière que les coordonnées dans la zone de définition de la projection soient simple d'utilisation.
Dans le cas du Lambert 93, on voulait que toutes les coordonnées en X sur la France métropolitaine soient positives et que les Y soient dans les 6 / 7 millions. Pourquoi cela ? Simplement parce que les systèmes de projection qui étaient utilisés avant (Lambert 2 étendu ou Lambert carto) avaient des Y dans les 1 à 4 millions. Donc pour éviter les confusions le nouveau système est dans les 6 millions ...
Pour en savoir plus, vous pouvez aussi lire ces pages : http://geodesie.ign.fr/
Bon courage
Hors ligne
#8 Wed 21 September 2016 12:27
- Thy06
- Juste Inscrit !
- Date d'inscription: 19 Sep 2016
- Messages: 4
Re: entête données GRID ASCII
Bonjour fbecir
Je crois qu'il faudrait avant toute chose que vous lisiez cet article :
https://fr.wikipedia.org/wiki/Projection_cartographique
vous me prenez pour un scientifique là ;-)
Mais c'est très intéressant.
Utilisation finale de ces données :
Je n'utilise pas de logiciel spécifique pour gérer/visualiser ce type de fichiers (.asc), auquel cas je n'aurais pas à me poser de questions.
Je vais traiter les fichiers .asc dans une base de données (FileMaker), pour associer les Z aux X et Y et exporter les données dans un autre format de fichier, spécifique au logiciel architecture que j'utilise (format GDL sur Archicad pour obtenir un objet 3D). Et donc pour placer correctement l'objet 3D correspondant à la zone de terrain que j'obtiendrai, il faudra que je sache comment le placer en XYZ dans ma scène par rapport à mon projet archi. D'où la nécessité de comprendre à quoi correspondent les valeurs xllcorner yllcorner.
Si vos données ne sont pas dans la même projection, vous ne pouvez pas les superposer. Une simple translation ou une transformation un peu plus subtile (affinité, homographie) ne vous tirera pas d'affaire non plus ...
Qu'entendez-vous par là ? La proposition de Jérémie N. n'est pas suffisamment précise ?
Comme dit plus haut, j'ai testé la solution, de manière approximative et ça fonctionne. Je retesterai par calculs précis dès que possible.
Cordialement,
Hors ligne