Annonce
#1 Sun 28 May 2017 06:21
- grizzli
- Juste Inscrit !
- Date d'inscription: 27 May 2017
- Messages: 4
Obtenir une matrice python d'élévations
Bonjour,
Je suis débutant sur QGis que j'ai installé dans le but d’extraire une matrice d'élévations.
J'entend par là un tableau de dimension 2 qui serait version discrétisée d'une portion de MNT,
avec une progression de l’abscisse x selon les lignes, et de l'ordonnée y selon les colonnes.
Pouvez vous m'indiquer comment faire ?
(j'ai mentionné QGis mais je me satisferait de tout logiciel gratuit ou avec version d'essai)
(exemple: --> x
[[57,58,59,75,70,63] y
,[55,56,60,70,69,55] ^
,[45,48,60,65,70,69] i
,[37,45,53,60,68,65]
,[35,46,55,63,75,58]
,[40,48,56,76,72,45]] où le prélèvement aurait été fait, par exemple, sur une portion carrée de 500m en 500m selon x et y )
Je vous remercie.
Hors ligne
#2 Mon 29 May 2017 11:11
- Pascal Boulerie
- Participant assidu
- Lieu: France
- Date d'inscription: 12 Sep 2005
- Messages: 2879
- Site web
Re: Obtenir une matrice python d'élévations
Définissez déjà ce que vous entendez ou comprenez par le terme spécifique de MNT .
par exemple d'après :
Modèle numérique de terrain
https://fr.wikipedia.org/wiki/Mod%C3%A8 … de_terrain
Sinon, le manuel QGIS donne-t-il un aperçu des formats de données gérés ? Et des fonctions de traitement du relief ?
Par "débutant", vous avez combien de temps pour arriver à maitriser ces questions ? (par exemple dans le cadre d'un stage étudiant)
« L'État est désormais quasi déliquescent. » (José Cohen-Aknine, ingénieur X-Ponts, IGPEF, dans Déliquescence et renaissance de l'État.)
Hors ligne
#3 Mon 29 May 2017 23:46
- grizzli
- Juste Inscrit !
- Date d'inscription: 27 May 2017
- Messages: 4
Re: Obtenir une matrice python d'élévations
Bonsoir,
Tout d'abord, merci d'avoir répondu.
Par "débutant", vous avez combien de temps pour arriver à maîtriser ces questions ? (par exemple dans le cadre d'un stage étudiant)
Je suis élève en prépa et cette démarche est pour une partie d'exposé (il me reste deux semaines environ).
Je débute réellement et le logiciel QGis me semble assez complexe. Je sais qu'il comporte une interface Python et j'ai installé le module en question.
Mais le guide d'utilisation est assez vaste et mon problème de matrice n'est évidemment pas mentionné. Par ailleurs, si j'ai bien compris, l'utilisation de Python se fait en travaillant sur des calques, ce que je ne maîtrise pas du tout. Le résultat est que je demeure bloqué avec ce problème (même s'il n'est probablement pas très compliqué).
Il serait très long de tout m'expliquer en détail. Mais à la limite, ce n'est pas nécessaire : une démarche que je pourrais suivre pas à pas me conviendrait parfaitement.
PS : Pour moi, MNT signifiait carte avec élévations, et je pensais par exemple au format .dem, disponible sur : https://earthexplorer.usgs.gov/
Merci.
Hors ligne
#4 Tue 30 May 2017 10:01
- Pascal Boulerie
- Participant assidu
- Lieu: France
- Date d'inscription: 12 Sep 2005
- Messages: 2879
- Site web
Re: Obtenir une matrice python d'élévations
Donc c'est pour un exposé de " TIPE ". Si c'est pour passer les concours 2017, il ne reste effectivement que peu de temps... ça aurait été mieux de commencer plus tôt...
MNT est le terme français correspondant à DEM pour Digital elevation model :
https://en.wikipedia.org/wiki/Digital_elevation_model
Moi, je ne suis pas allé voir ce qu'il y a en catalogue sur le site USGS, mais c'est surtout dans le cas de figure de passer d'un réseau irrégulier de triangles vers une matrice que je m'inquiétais de la difficulté (TIN : Triangulated irregular network)
https://en.wikipedia.org/wiki/Digital_e … pes_of_DEM
https://fr.wikipedia.org/wiki/Mod%C3%A8 … pes_de_MNT
https://en.wikipedia.org/wiki/Triangula … ar_network
Dans QGIS, il y a une interface pour traiter les images avec le logiciel GRASS, peut-être que les MNT sont à traiter dans GRASS (mais ça fait longtemps ou pratiquement jamais que j'ai pratiqué avec GRASS, donc je ne saurais vous dire quelles sont les fonctions à utiliser).
Pour la documentation :
http://docs.qgis.org/2.14/fr/docs/user_manual/
Analyse Raster
http://docs.qgis.org/2.14/fr/docs/user_ … =elevation
Intégration du SIG GRASS
http://docs.qgis.org/2.14/fr/docs/user_ … ation.html
Vous n'avez besoin que d'extraire une zone rectangulaire à partir d'une zone plus grande, en conservant tous les points initiaux, et la précision de ces données initiales,
pas de ré-échantillonner en calculant de nouvelles altitudes moyennes pour obtenir un MNT moins détaillé ?
« L'État est désormais quasi déliquescent. » (José Cohen-Aknine, ingénieur X-Ponts, IGPEF, dans Déliquescence et renaissance de l'État.)
Hors ligne
#5 Tue 30 May 2017 13:26
- olivier.eckmann
- Participant assidu
- Date d'inscription: 14 Dec 2007
- Messages: 351
Re: Obtenir une matrice python d'élévations
Bonjour,
avec IGNMap tu peux charger des MNT de types ASC (Ascii Grid Esri) ou TIF (par exemple la BDAlti au pas de 75m est gratuite et disponible au format ASC en Lambert93), puis tu peux reprojeter (éventuellement avec le même système de projection) en limitant la zone (rectangulaire) et en spécifiant le nouveau pas de la grille, ainsi que le format ASC ou XYZ
Olivier
Hors ligne
#6 Wed 31 May 2017 00:11
- grizzli
- Juste Inscrit !
- Date d'inscription: 27 May 2017
- Messages: 4
Re: Obtenir une matrice python d'élévations
Bonsoir,
Merci pour vos réponses.
Oui, effectivement, importer ces cartes me serait utile pour mon TIPE.
J'ai essayé d'une part de télécharger IGNMap, et je suis parvenu à exporter depuis des .asc des formats jpeg, ce qui me servirait de dernier recours (avec la matrice de l'image) si je ne parvenais pas à un résultat plus précis et plus rapide à mettre en oeuvre.
Ce que je recherche est en effet une extraction de données les plus précises possible, quitte à moyenner ensuite sur Python, (ce qui ne pose plus problème).
J'ai rapidement regardé les guides d'utilisation et je prendrais le temps d'essayer le module Grass, mais cette extension n'a pas l'air simple non plus. Pensez vous que je devrais reporter ce sujet dans la partie QGis du forum ?
Merci.
Hors ligne
#7 Wed 31 May 2017 09:45
- Pascal Boulerie
- Participant assidu
- Lieu: France
- Date d'inscription: 12 Sep 2005
- Messages: 2879
- Site web
Re: Obtenir une matrice python d'élévations
Que voulez-vous dire par "les plus précises possible" ??? (parce que ça peut vouloir dire plein de choses, par exemple que vous voulez faire un relevé à l'échelle 1/100 d'un chantier archéologique par exemple...)
Quel est le contenu d'un format ASC ?? (spécifications détaillées du format)
des points X Y Z (en 3D), je suppose...
X1 Y1 Z1
X2 Y2 Z2
(ce n'est pas détaillé dans la page https://www.fileext.com/file-extension/asc.html )
> Ascii Grid Esri
Si c'est un format établi par la société ESRI , il doit être documenté quelque part. Peut-être dans un autre message :
https://georezo.net/forum/search.php
Pourquoi le convertir en JPEG ?
« L'État est désormais quasi déliquescent. » (José Cohen-Aknine, ingénieur X-Ponts, IGPEF, dans Déliquescence et renaissance de l'État.)
Hors ligne
#8 Wed 31 May 2017 10:41
Re: Obtenir une matrice python d'élévations
Bonjour,
Le format ASCII ESRI est décrit ici : http://desktop.arcgis.com/fr/arcmap/lat … format.htm
A+
Franck
Hors ligne
#9 Thu 01 June 2017 09:41
- olivier.eckmann
- Participant assidu
- Date d'inscription: 14 Dec 2007
- Messages: 351
Re: Obtenir une matrice python d'élévations
Bonjour,
si c'est pour utiliser sans interface, regarde du côte de GDAL_TRANSLATE : http://www.gdal.org/gdal_translate.html
Olivier
Hors ligne
#10 Sat 03 June 2017 17:05
- grizzli
- Juste Inscrit !
- Date d'inscription: 27 May 2017
- Messages: 4
Re: Obtenir une matrice python d'élévations
Bonjour,
Je vais voir si je peux utiliser cette méthode.
Sinon, dans tout les cas, les images elles mêmes donnent de bons résultats.
Merci beaucoup à tous.
Hors ligne
