Suite

Comment obtenir une courbe hypsométrique à partir des données d'altitude ?

Comment obtenir une courbe hypsométrique à partir des données d'altitude ?


J'ai mes données d'altitude mais je ne peux pas en faire une table attributaire et je ne sais donc pas comment obtenir ces données dans une table afin de créer une courbe hypsométrique.


Le raster est stocké au format à virgule flottante sous forme de valeurs individuelles dans chaque cellule ; il n'a donc pas de table attributaire. Une solution efficace consiste à discrétiser les valeurs et à les convertir au format entier. Cela crée une table attributaire qui, une fois tracée, produira un histogramme. Ses sommes cumulées sont la courbe hypsométrique.

L'apparence de la courbe hypsométrique (sous forme CDF) pour cet échantillon DEM (une partie de l'Illinois disponible sur http://exampledata.wolfram.com/ArcGRID.zip) varie avec la finesse de la discrétisation de la valeur. Les images suivantes utilisent des bacs de 5, 1 et 1/5 mètres, respectivement. Comme vous pouvez le voir, on peut obtenir une courbe précise en utilisant des bacs moyennement grossiers.

Discrétiser une grille implique deux opérations qui peuvent être effectuées en une seule étape : diviser par la largeur du bac et tronquer (ou arrondir, si vous préférez). Par exemple, les altitudes en mètres peuvent être discrétisées en incréments de 0,1 mètre via une division par 0,1 suivie d'une troncature, comme dans

Int( [élévation] / 0,1 )

Cette syntaxe, ou quelque chose d'assez proche, serait utilisée dans presque toutes les versions d'ArcGIS (et dans de nombreux autres SIG raster également).


J'utilise normalement R pour cela. Vous pouvez lire les ensembles de données raster GDAL en utilisantrgdal(à partir du CRAN), puis construisez une distribution cumulative empirique pour les valeurs d'altitude en utilisantecdf(intégré).

Par exemple, j'ai mon DEM dans un fichier GeoTIFF, en utilisant R :

library(rgdal) topo <- readGDAL("mytoporaster.tif") plot(ecdf(topo$band1), main="Courbe hypsométrique", xlab="Élévation (m)")


Voir aussi le tutoriel à Carleton.


Je suis tombé sur cet article qui parle d'une extension ArcGIS 9.x.

CalHypso : Une extension ArcGIS pour calculer les courbes hypsométriques et leurs moments statistiques. Applications à l'analyse des bassins versants dans le sud-est de l'Espagne

L'extension avec des exemples de données est fournie à cet emplacement de téléchargement.

Pour mémoire l'extension fonctionne sous ArcGIS 10.4.1.


Voir la vidéo: la courbe hypsometrique et la classe des altitudes by LAZRAK