Le portail francophone de la géomatique

Merci de la réactivité cependant je ne trouve pas la fonction courbure. J'utilise la version 6.4 de Grass.

Bonjour,

c'est possible avec le module grass r.slope.aspect (menu Raster > Analyse de terrain > Pente et exposition; http://grass.osgeo.org/grass64/manuals/ … pect.html) qui permet à priori de faire exactement les mêmes calculs que le module QGIS cité ...
il y a également les mêmes fonctionnalités avec r.param.scale (menu Raster > Analyse de terrain > Paramètres de terrain;  http://grass.osgeo.org/grass64/manuals/ … ale.html), qui utilise des algorithmes un peu différents ...

Sylvain
PS: si certaines fonctionnalités sont difficiles à trouver, n'hésitez pas à me proposer des noms de menus différents, ça sera pris en compte pour les prochaines versions de GRASS

Bonjour,

est-ce que ça ne serait pas un problème lié à vos paramètres de région ? Elle doit être fixée sur l'emprise et la résolution de votre mnt, sinon vous risquez d'avoir des temps de calculs démesurés si la maille est trop petite, et/ou des données incohérentes ...

en diminuant la résolution du raster ou bien en se concentrant sur une étendue plus petite, les temps de calcul seront forcément plus courts.

dans ton exemple, tu utilises r.slope.aspect alors qu'il me semble que ton besoin nécessite plutôt d'utiliser les courbures (curvature).

Pour r.param.scale, les par. sont à sélectionner parmi ceux-là:
# profc or profile curvature (curvature intersecting with the plane defined by Z axis and maximum gradient direction). Positive values describe convex profile curvature, negative values concave profile curvature.
# planc or plan curvature (horizontal curvature, intersecting with the XY plane)
# longc or longitudinal curvature (profile curvature intersecting with the plane defined by the surface normal and maximum gradient direction)
# crosc or cross-sectional curvature (tangential curvature intersecting with the plane defined by the surface normal and a tangent to the contour - perpendicular to maximum gradient direction)
# maxic or maximum curvature (can be in any direction)
# minic or minimum curvature (in direction perpendicular to the direction of of maximum curvature)

=> voir les 2 liens vers les pages de manuel que j'ai donné plus haut:
- pour r.slope.aspect : pcurv, tcurv, ou les paramètres de dérivées partielles
- pour r.param.scale : comme l'a dit baobazz, profc, planc, ...


Sylvain

Plutôt que pentes et expositions, le calcul des bassins versants ne donnerait-il pas, à leurs jonctions, les crêtes ?

ben oui (...)

sous grass: r.watershed basin

Puis conversion des polygones en lignes

Pour caractériser la concavité de la crête, il faudra utiliser la courbure. Avoir les deux: lignes + raster dessous avec valeurs de courbure peut être intéressant

Dernière modification par baobazz (Wed 14 November 2012 14:26)

En l’occurrence, si on fait les choses bien, on commence plutôt par l'exutoire du bassin versant avant tout, vu que c'est ce qui va définir ledit bassin... Les lignes de crêtes viennent en conséquence et  peuvent être automatiquement déduites de cet exutoire, enfin sauf si on commence à étudier les bassins de versant c'est un peu plus ardu.

Je parlais niveau méthodologique - à moins qu'on ne cherche à déterminer des BV par pixels contributifs minimum ou autre
Ceci dit, je n'utilise pas les pentes/ruptures de pentes directement, mais plutôt les directions et les accumulations de flux... qui sont issus de la direction de la pente entre deux pixels, certes...