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 Thu 20 July 2017 15:51

Anne-Sophie_RD
Participant occasionnel
Lieu: Bagnères-de-Bigorre
Date d'inscription: 2 Apr 2012
Messages: 36

QGIS 2.18: sélection d'un ensemble de "sommets" depuis un MNT

Bonjour à tous,

Cela fait un moment que je planche sur le sujet sans trouver de solution satisfaisante, alors je fais appel à la communauté de GeoRezo pour trouver des pistes.

Pour les besoins d'une étude, je dois pré-sélectionner sur ma zone d'étude (massif des Pyrénées) des groupes de 4 "sommets" proches répondant à plusieurs critères :
1) chacun des 4 sommets doit se trouver dans un étage de végétation différent (ces étages sont définis par des limites d'altitude)
2) les 4 sommets doivent être à peu près distribués sur des intervalles d'altitude égaux
3) les 4 sommets doivent être dans une même zone climatique et sur un même type de roche mère
4) les 4 sommets ne doivent pas être soumis à l'impact d'activités humaines (pastoralisme ou tourisme)

Est défini comme "sommet" non seulement le point le plus haut d'une montagne (pic), mais aussi toute "bosse" sur une crête qui dépasse de plus de 20m les terrain environnants.

J'ai à ma disposition un MNT de résolution 25m

Mon idée est de procéder de la manière suivante : dans un premier temps extraire du MNT tous mes "sommets" et dans un second temps croiser avec les différents critères pour sélectionner des groupes de sommets qui conviennent.

Je bloque déjà sur ce premier temps. Je n'arrive pas à extraire de mon MNT les sommets. J'ai tenté d'utiliser la commande GRASS r.param.scale dans la boite à outil de QGis mais cela donne des résultats étranges (avec param=feature, j'obtiens des pixels identifiés sommets dans des zones un peu improbables en plaine). Après, je ne maitrise pas vraiment le paramétrage de cette fonction.
Peut-être faut-il combiner plusieurs analyses pour arriver à la bonne sélection ?

En bref, je rame un peu, je n'ai que peu d'expérience dans l'analyse raster. Je sens que quelque chose est possible (avec tous les outils disponibles avec GRASS et la calculatrice raster qui semble si puissante!) mais je n'arrive pas à trouver la bonne logique ou la bonne procédure pour arriver à faire ce que je veux.

Je suis donc preneuse de toute piste, tout conseil, qui me permettrait d'avancer un peu sur ma problématique.

Anne-Sophie

Hors ligne

 

#2 Thu 20 July 2017 16:21

MathieuR
Membre
Lieu: aix-en-provence
Date d'inscription: 16 Feb 2009
Messages: 1690
Site web

Re: QGIS 2.18: sélection d'un ensemble de "sommets" depuis un MNT

La méthode r.param.scale est celle que j'aurais conseillée.

Votre MNT est-il assez précis selon vos besoins, notamment par rapport à l'exigence "20 cm" de plus que les terrains environnants. D'ailleurs, l'usage du terme environnants appelle cette question : selon quelle étendue (bandwidth) considère-t-on l'environnement à considérer ?

La précision de votre MNT est également à confronter par rapport à vos couches d'occupation de sol et climatiques.

Dans l'utilisation de r.param.scale, un facteur déterminant est le paramètre size qui donne la taille de la fenêtre. Par défaut, c'est 3. Ce qui veut dire que votre fenêtre fera 3*25 = 75 m de large. Plus la taille est large, et plus les résultats seront "lissés"


geodata au cerema et petits billets en géomatique

Hors ligne

 

#3 Fri 21 July 2017 15:21

Anne-Sophie_RD
Participant occasionnel
Lieu: Bagnères-de-Bigorre
Date d'inscription: 2 Apr 2012
Messages: 36

Re: QGIS 2.18: sélection d'un ensemble de "sommets" depuis un MNT

Merci pour cette réponse.

Le MNT de résolution 25m est le plus précis à ma disposition, je vais devoir faire ce que je peux avec.

Par ailleurs, l'exigence est de 20m (pas 20cm, faute de frappe?). Le protocole que je dois suivre pour la sélection des sommets ne précise rien quant à l'étendue de l'environnement à considérer. Tout reste malheureusement assez flou et très théorique.

Concernant les autres couches avec lesquelles je souhaite croiser les sommets, je suis encore en cours de recherche mais je ne pense pas trouver des informations bien précises.

Je vais faire quelques tests avec r.param.scale en faisant varier le paramètre size pour voir ce que cela donne...

Hors ligne

 

Pied de page des forums

Powered by FluxBB