Wyodrębniasz obszar DEM powyżej, gdzie 60% regionu leży za pomocą ArcGIS Desktop?

9

Mam DEM regionu górskiego. Chciałbym wyodrębnić obszary DEM, powyżej których leży 60% całkowitego obszaru DEM. Zasadniczo chciałbym wydobyć obszary reprezentujące górne partie gór.

Słyszałem o tym jako o linii hipsometrycznej (linia H60), ale nie byłem w stanie znaleźć żadnej dokumentacji, jak to osiągnąć.

Moim instynktem jest jakoś obliczyć statystyki DEM, aby powiedzieć mi, gdzie jest poziom 60%, a następnie obrysować DEM na tej wysokości. Następnie mogę przyciąć DEM za pomocą tego konturu i wykonać moje obliczenia na wyjściowych podzbiorach DEM.

Jak ustalić ten poziom 60%?

Radar
źródło

Odpowiedzi:

7

Uważam, że po prostu szukasz wartości 60 percentyla dla swojego rastra. Najłatwiejszym sposobem, aby to zrobić, jest prawdopodobnie funkcja histogramu. Kliknij warstwę prawym przyciskiem myszy w Spisie treści, a następnie przejdź do Właściwości warstwy> Symbolika> Sklasyfikowane> Klasyfikuj. Wybierz Quantile jako metodę i użyj 10 klas. Szósta klasa rozpocznie wówczas od wartości 60% percentyla. Zobacz przykład poniżej.

Jeśli chcesz określić inne wartości procentowe, możesz nacisnąć symbol „%” po prawej stronie okna Klasyfikacji, obok napisu „Wartości podziału”. wprowadź opis zdjęcia tutaj

dmahr
źródło
3
Zgadza się, z tym wyjątkiem, że Radar będzie szukał 40. percentyla zamiast 60. - tylko 40% powierzchni będzie powyżej 60. percentyla.
nmpeterson
6
Strzeż się błędów ArcGIS! Nie testowałem Arc 10, ale poprzednie wersje nie były w stanie poprawnie obliczyć statystyk. Czy zauważyłeś, jak na tym zrzucie ekranu jest „9667 elementów w klasie” wśród „Count:” z 85842? To znacznie więcej niż 10% (co powinno być). Nie ma nic złego w używaniu histogramu do zgadywania na 40 percentylu (lub cokolwiek innego), ale powinieneś następnie sprawdzić to przypuszczenie (używając obliczeń Algebry Map), aby upewnić się, że jest poprawne.
whuber
Och wow ... dzięki za wskazanie tego Whubera. Nigdy nie wiedziałem o tym błędzie. Westchnienie ... trudno uwierzyć, że ArcMap nie może niezawodnie robić tak prostych rzeczy.
dmahr
Świetnie - właśnie tego szukałem. Dopilnuję podwójnego sprawdzenia wartości.
Radar
4
@dmahr Jeśli chcesz uwzględnić wszystkie próbki, zwiększ rozmiar „Próbkowanie ...” do dowolnej wartości powyżej „liczby”; 90 000 powinno działać dla informacji z powyższego zrzutu ekranu. Domyślnie ArcGIS będzie zawierać tylko pierwsze 10 000 rekordów / wartości.
SaultDon
3

Alternatywnym podejściem do proponowanej ekstrakcji DEM, przy założeniu, że masz Spatial Analyst, byłoby:

  1. Użyj narzędzia „Plasterek” z metodą „RÓWNYM OBSZAREM” i 5 lub 10 strefami.
  2. Następnie możesz:

    a) Ponownie sklasyfikuj wynik Slice, aby utworzyć maskę do wyodrębnienia DEM

    - lub -

    b) Użyj go w narzędziu „CON” lub kalkulatorze rastrowym, aby zamaskować swój DEM w górnych 3 lub 6 strefach (odpowiednio do liczby zdefiniowanych stref plasterków)

Ponadto Slice wydaje się uzyskiwać nieco dokładniejsze kwantyle niż metoda klasyfikacji kwantyli (jak wskazał @whuber 28.11.2011 20:29)

Joe Champagne
źródło