Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
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

Printemps des cartes 2024

#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

nitram
Participant assidu
Lieu: Brest
Date d'inscription: 19 May 2009
Messages: 989
Site web

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

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)


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

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

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

fbecir a écrit:

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.

fbecir a écrit:

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

 

Pied de page des forums

Powered by FluxBB