Шкала та коефіцієнт Z не впливають на аналіз схилу схилів у QGIS

10

Я використовував SRTM DEM для створення схилу схилу в QGIS 2.2.0. Однак результуюча схила гори здається вертикально перебільшеною. Крім того, зміна параметрів масштабу та z-коефіцієнта не впливає на результуючу схилу схилу. Я спробував обидва варіанти в QGIS, Analysis-> DEM (Terrain models) та Terrain Analysis -> Hillshade. Обидва методи дали однакові результати.

Це поодинокий випадок, помилка чи я щось пропускаю?

qgis scale hillshade

— PyMapr
джерело

4

Якщо ви використовуєте DEM з координатами латлону в десяткових градусах і вертикальними одиницями в метрах, ви повинні встановити масштаб (відношення вертикальних одиниць до горизонтальних одиниць) до 111120. Я намагався з парою DEM від Aster, у мене немає проблеми. Ця настройка застосовується не тільки для схилів, але і для розрахунків схилів. Сподіваюсь, що це допомагає

— Херардо Хіменес

@GerardoJimenez, я спробував це, і це спрацювало. Спасибі

— PyMapr

@GerardoJimenez, Thanderbolt будь-хто з вас має подати та відповісти на це питання, а Thunderbolt повинен прийняти його, щоб видалити це питання без відповіді.

— SS_Rebelious

12

@ Херардо в відповідь правильний тільки на екваторі. Документ gdaldem відзначає, що значення масштабу призначені для "LatLong-проекцій поблизу екватора"

Для обчислення шкали або коефіцієнта Z (оберненого масштабу) на інших широтах для dem з вертикальними одиницями в метрах можна використовувати такі рівняння:

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

де: 111320 - довжина одного градуса на екваторі в метрах (виходячи з WGS84, ви можете використовувати 111111, якщо це легше запам'ятати, або якщо ви розумієте історію ), а широта - широта в центрі вашої DEM.

Ви можете швидко обчислити коефіцієнт масштабу / z у Excel / Libre Office Calc за допомогою таких формул:

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))

Ви навіть можете використовувати простий пошук у Google:

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

Напр. Масштаб і z коефіцієнт z для Широти = 38,5S

Для DEM з вертикальними одиницями в стопах просто перетворіть довжину одного градуса на екваторі в метрах в фути, тобто 3,28 * 111320 = 365130. Отже, ваші рівняння стають:

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

Крім того, ви можете повторно відхилити DEM в проектовану координовану систему (також запропоновану документами gdaldem).

— user2856
джерело

Ця відповідь мені справді допомогла сьогодні, але я трохи розгублений щодо термінів, які використовуються в QGIS. Коли я переходжу до Растрового - Аналіз місцевості - Нахил, я отримую поле під назвою "фактор Z", але тут я повинен ввести значення, яке тут називається "шкала", а не зворотне значення, яке тут називається "коефіцієнт Z". Ви не можете ввести дуже маленьке число, яке я отримую при обчисленні зворотного з вищевказаним обчисленням, і значення за замовчуванням вже говорить вам, що це повинно бути велика кількість, але його все ще трохи заплутано, тому що я подумав, що і термін "шкала" проти "z фактор" не вдалося використати синонімічно.

— Тая

@Taya Я не бачу поля "Фактор Z" в інструменті Raster-> Analysis-> DEM (Terrain Models) у меню, наведеному в QGIS 2.14. Я бачу лише поле масштабу. Можливо, це була помилка в попередніх версіях QGIS ...? Редагувати: А-а, ви використовуєте плагін для аналізу місцевості , це повинен бути помилка в плагіні.

— користувач2856

Я б запропонував кілька змін у цій відповіді: 1. Для більшості застосувань ви б використовували коефіцієнт масштабу s або z-фактор, а не обидва (s здається більш логічним вибором і краще працює в QGIS). 2. Коефіцієнт шкали північ - Південь завжди є 111320, але коефіцієнт масштабу Схід - Захід є 111320*cos(latitude*pi/180). На великих широтах ця різниця створить значні проблеми, і вам, ймовірно, доведеться перепроектувати систему координат лічильника (наприклад, UTM) перед розрахунком. У низьких широтах, ви можете спробувати розщеплення різницю: s = 111320 * 0.5 * (1 + cos(latitude*pi/180)).

— Маттіас Фріпп

4

Під час роботи з DEM в lat довгих координатах з їх одиницями в десяткових градусах і вертикальних одиницях в метрах ви повинні встановити "шкала (відношення вертикальних одиниць до горизонтальних одиниць) до 111120. В

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

вони навіть припускають, що для цього параметра для DEM з координатами ланцюга та вертикальними одиницями в футах шкала повинна бути встановлена на 370400

Ці параметри впливають також на розрахунки нахилу.

— Херардо Хіменес
джерело