Extraire la zone de DEM ci-dessus où 60% de la région se trouve à l'aide d'ArcGIS Desktop?

J'ai un DEM d'une région montagneuse. Je voudrais extraire les zones du DEM au-dessus desquelles 60% de la surface totale du DEM se trouve. Donc, je voudrais essentiellement extraire les zones représentant les parties supérieures des montagnes.

J'ai entendu parler de ce que l'on appelle la ligne hypsométrique (ligne H60), mais je n'ai pu trouver aucune documentation sur la façon d'y parvenir.

Mon instinct est de calculer en quelque sorte les statistiques du DEM pour me dire où se trouve le niveau de 60%, puis de contourner le DEM à cette altitude. Je peux ensuite découper le DEM par ce contour et effectuer mes calculs sur les sous-ensembles DEM de sortie.

Comment déterminer ce niveau de 60%?

— Radar
source

Réponses:

Je crois que vous recherchez simplement la valeur du 60e centile pour votre raster. La façon la plus simple de procéder consiste probablement à utiliser la fonctionnalité d'histogramme. Cliquez avec le bouton droit sur le calque dans la table des matières, puis sélectionnez Propriétés du calque> Symbologie> Classé> Classifier. Choisissez Quantile comme méthode et utilisez 10 classes. La 6e classe commencera alors à la valeur de centile de 60%. Voir l'image ci-dessous pour un exemple.

Si vous souhaitez spécifier d'autres valeurs de pourcentage, vous pouvez appuyer sur le symbole "%" sur le côté droit de la fenêtre de classification, à côté de l'endroit où il est dit "Break Values". entrez la description de l'image ici

— dmahr
source

C'est vrai, sauf que Radar recherchera le 40e centile au lieu du 60e - seulement 40% des terres seraient au-dessus du 60e centile.

— nmpeterson

Méfiez-vous des bogues ArcGIS! Je n'ai pas testé Arc 10, mais les versions précédentes n'étaient pas en mesure de calculer correctement les statistiques. Avez-vous remarqué comment cette capture d'écran indique qu'il y a "9667 éléments dans la classe" parmi un "Count:" de 85842? C'est sensiblement plus de 10% (ce qui devrait être le cas). Il n'y a rien de mal à utiliser l'histogramme pour deviner le 40e centile (ou autre), mais vous devez ensuite vérifier cette supposition (en utilisant les calculs d'algèbre de carte) pour vous assurer qu'elle est correcte.

— whuber

Oh wow ... merci d'avoir signalé ça Whuber. Je n'ai jamais connu ce bug. Soupir ... il est difficile de croire qu'ArcMap ne peut pas faire des choses aussi simples de manière fiable.

— dmahr

Super - c'est exactement ce que je cherchais. Je m'assurerai de revérifier les valeurs.

— Radar

@dmahr Si vous souhaitez inclure "tous" les échantillons, augmentez la taille de "Sampling ..." à tout ce qui est supérieur au "count"; 90 000 devraient fonctionner pour les informations de la capture d'écran ci-dessus. Par défaut, ArcGIS inclura uniquement les 10 000 premiers enregistrements / valeurs.

— SaultDon

Une autre approche de l'extraction DEM proposée, en supposant que vous ayez Spatial Analyst, serait de:

Utilisez l'outil 'Slice' avec la méthode "EQUAL AREA" et 5 ou 10 zones.
Vous pouvez alors:

a) Reclassifiez la sortie Slice pour créer un masque pour extraire votre DEM

-- ou --

b) Utilisez-le avec l'outil 'CON' ou la calculatrice raster pour masquer votre DEM sur les 3 ou 6 premières zones (respectivement au nombre de zones de tranche définies)

De plus, Slice semble obtenir des quantiles un peu plus précis que la méthode de classification Quantile (comme indiqué par @whuber le 2011-11-28 20:29)

— Joe Champagne
source