Skala i współczynnik Z nie mają wpływu na analizę cienia wzgórza w QGIS

10

Użyłem SRTM DEM, aby stworzyć odcień wzgórza w QGIS 2.2.0. Jednak powstały cień wydaje się być przesadzony w pionie. Również zmiana ustawień skali i współczynnika Z nie ma żadnego wpływu na wynikowy cień. Wypróbowałem obie opcje w QGIS, Analysis-> DEM (modele terenu) i Analiza terenu -> Hillshade. Obie metody dały takie same wyniki.

Czy to izolowany incydent, błąd, czy coś mi brakuje?

qgis scale hillshade PyMapr
źródło
4
Kiedy używasz DEM ze współrzędnymi latlonu w stopniach dziesiętnych i jednostkami pionowymi w metrach, musisz ustawić skalę (stosunek jednostek pionowych do jednostek poziomych) na 111120. Próbowałem z kilkoma DEM z Asteru, nie mam problemu. To ustawienie dotyczy nie tylko cieniowania wzgórza, ale także obliczeń nachylenia. Mam nadzieję, że to pomoże
Gerardo Jimenez
@GerardoJimenez, próbowałem i działało. Dzięki
PyMapr
@GerardoJimenez, Thanderbolt, którykolwiek z was powinien przesłać i odpowiedzieć na to pytanie, a Thunderbolt powinien je zaakceptować, aby usunąć to pytanie bez odpowiedzi.
SS_Rebelious

Odpowiedzi:

12

Odpowiedź Gerarda jest poprawna tylko na równiku. Gdaldem doc zauważa te wartości skali są za „LatLong projekcji pobliżu równika”

Aby obliczyć skalę lub współczynnik Z (odwrotność skali) dla innych szerokości geograficznych dla dem z jednostkami pionowymi w metrach, możesz użyć następujących równań:

s=111320*cos(latitude*pi/180)
z=1/(111320*cos(latitude*pi/180))

gdzie: 111320 to długość jednego stopnia na równiku w metrach (w oparciu o WGS84, możesz użyć 111111, jeśli jest to łatwiejsze do zapamiętania lub znasz historię ), a szerokość geograficzna to szerokość geograficzna w centrum DEM.

Możesz szybko obliczyć współczynnik skali / z w programie Excel / Libre Office Calc za pomocą następujących wzorów:

For -s scale parameter
=111320*COS(latitude*PI()/180)
or
=111320*COS(RADIANS(latitude))

For -z z factor parameter
=1/(111320*COS(latitude*PI()/180)
or
=1/(111320*COS(RADIANS(latitude))

Możesz nawet użyć prostej wyszukiwarki Google:

111320 * cos(latitude deg)
or
1/(111320 * cos(latitude deg))

Np. Skala i współczynnik z dla Latitude = 38,5S

W przypadku DEM z pionowymi jednostkami w stopach, wystarczy przekonwertować długość jednego stopnia na równiku w metrach na stopy, tj. 3,28 * 111320 = 365130. Zatem twoje równania stają się:

s=365130*cos(latitude*pi/180)
z=1/(365130*cos(latitude*pi/180))

Alternatywnie możesz ponownie skierować swój DEM do rzutowanego skoordynowanego systemu (sugerowanego również przez dokumentację gdaldem).

użytkownik2856
źródło
Ta odpowiedź naprawdę mi dziś pomogła, ale jestem trochę zdezorientowany terminami używanymi w QGIS. Kiedy idę do Rastra - Analiza terenu - Nachylenie, otrzymuję pole o nazwie „Współczynnik Z”, ale tutaj muszę wprowadzić tutaj wartość określaną jako „skala”, a nie odwrotną wartość, która jest tutaj nazywana „współczynnikiem Z”. Nie możesz wprowadzić bardzo małej liczby, którą otrzymuję przy obliczaniu odwrotności z powyższym obliczeniem, a wartość domyślna już mówi ci, że musi to być duża liczba, ale nadal jest trochę myląca, ponieważ myślałem, że oba określenie „skala” w porównaniu z „współczynnikiem z” nie można go używać synonimicznie.
Taya
@Taya Nie widzę pola „Z Factor” w narzędziu Raster-> Analiza-> DEM (modele terenu) z menu w QGIS 2.14. Widzę tylko pole skali. Być może był to błąd we wcześniejszych wersjach QGIS ...? Edycja: Ahh, używasz wtyczki analizy terenu , która musi być błędem we wtyczce.
user2856 18.04.2016
Sugerowałbym kilka zmian w tej odpowiedzi: 1. W przypadku większości aplikacji użyłbyś współczynnika skalowania s lub współczynnika z, a nie obu (wydaje się to bardziej logicznym wyborem i działa lepiej w QGIS). 2. Współczynnik skali północ-południe jest zawsze 111320, ale współczynnik skali wschód-zachód wynosi 111320*cos(latitude*pi/180). Na dużych szerokościach geograficznych ta różnica spowoduje znaczące problemy i przed obliczeniem prawdopodobnie trzeba będzie ponownie wykonać rzut na układ współrzędnych licznika (np. UTM). W niskich szerokościach geograficznych, możesz spróbować dzieląc różnicę: s = 111320 * 0.5 * (1 + cos(latitude*pi/180)).
Matthias Fripp,
4

Podczas pracy z DEM w współrzędnych długich lat z ich jednostkami w stopniach dziesiętnych i jednostkami pionowymi w metrach należy ustawić „skalę (stosunek jednostek pionowych do jednostek poziomych) na 111120. W

http://www.gdal.org/gdaldem.html

sugerują nawet, że to ustawienie dla DEM ze współrzędnymi lantlonowymi i jednostkami pionowymi w stopach, skala musi być ustawiona na 370400

To ustawienie wpływa również na obliczenia nachylenia.

Gerardo Jimenez
źródło