Ölçeğin ve Z faktörünün QGIS'de yamaç analizi üzerine etkisi yoktur

10

QGIS 2.2.0'da bir tepe oluşturmak için bir SRTM DEM kullandım. Bununla birlikte, sonuçta ortaya çıkan tepe gölgesi dikey olarak abartılı gibi görünmektedir. Ayrıca, ölçek ve z faktörü ayarlarının değiştirilmesinin, elde edilen tepe gölgesi üzerinde herhangi bir etkisi yoktur. Her iki seçeneği de QGIS, Analiz-> DEM (Arazi modelleri) ve Arazi analizi -> Hillshade'de denedim. Her iki yöntem de aynı sonuçları verdi.

Bu izole bir olay mı, bir hata mı yoksa bir şey mi kaçırıyorum?

qgis scale hillshade

— PyMapr
kaynak

4

Ondalık derecelerde latlon koordinatları ve metre cinsinden dikey birimleri olan bir DEM kullanırken ölçeği (dikey birimlerin yatay birimlere oranı) 111120'ye ayarlamanız gerekir. Aster'den birkaç DEM ile denedim. Sorun değil. Bu ayar sadece tepe gölgeleme için değil aynı zamanda eğim hesaplamaları için de geçerlidir. Umarım yardımcı olur

— Gerardo Jimenez

@GerardoJimenez, denedim ve işe yaradı. Teşekkürler

— PyMapr

@GerardoJimenez, Thanderbolt ya bu soruya cevap vermeli ve bu soruyu cevaplamalısınız ve Thunderbolt bu soruyu cevapsızdan kaldırmak için kabul etmelidir.

— SS_Rebelious

12

@ Gerardo'nun yanıtı sadece ekvatorda doğrudur. Gdaldem doc olanlar ölçekli değerler içindir notları "ekvatora yakın LatLong projeksiyonlar"

Metre cinsinden dikey birimleri olan bir dem için ölçeği veya Z faktörünü (ölçeğin tersi) diğer enlemlerde hesaplamak için aşağıdaki denklemleri kullanabilirsiniz:

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

burada: 111320, ekvatorda metre cinsinden bir derecenin uzunluğudur (WGS84'e dayanarak, hatırlanması daha kolaysa veya tarihi takdir ediyorsanız 111111'i kullanabilirsiniz ) ve enlem, DEM'inizin merkezindeki enlemdir .

Excel / Libre Office Calc içindeki ölçek / z faktörünü aşağıdaki formüllerle hızlı bir şekilde hesaplayabilirsiniz:

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

Basit bir Google araması bile kullanabilirsiniz:

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

Örneğin , Latitude için ölçek ve z faktörü = 38.5S

Dikey birimleri ayak olan DEM'ler için, ekvatordaki bir derecenin uzunluğunu metre cinsinden metreye dönüştürün, yani 3.28 * 111320 = 365130. Böylece denklemleriniz:

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

Alternatif olarak, DEM'inizi yansıtılan bir koordineli sisteme (gdaldem dokümanları tarafından da önerilmektedir) yeniden yansıtabilirsiniz.

— user2856
kaynak

Bu cevap bugün bana gerçekten yardımcı oldu, ancak QGIS'te kullanılan terimler hakkında biraz kafam karıştı. Raster - Arazi Analizi - Eğim'e gittiğimde "Z faktörü" adlı alanı alıyorum, ancak buraya "ölçek" olarak adlandırılan değeri buraya "Z faktörü" olarak adlandırılan ters değeri girmem gerekiyor. Yukarıdaki hesaplama ile tersi hesaplarken elde ettiğim çok küçük sayıyı giremezsiniz ve varsayılan değer zaten çok büyük bir sayı olması gerektiğini söyler, ancak yine de biraz kafa karıştırıcıdır çünkü "ölçek" ile "z faktörü" terimini düşündüm eşanlamlı olarak kullanılamadı.

— Taya

@Taya QGIS 2.14'teki menüden Raster-> Analiz-> DEM (Arazi Modelleri) aracında "Z Faktörü" alanını görmüyorum. Yalnızca bir Ölçek alanı görüyorum. Belki de QGIS'in önceki sürümlerinde bir hataydı ...? Düzenleme: Ahh, bu eklenti bir hata olmalı arazi analizi eklentisi kullanıyorsunuz .

— user2856

Bu cevapta birkaç değişiklik öneriyorum: 1. Çoğu uygulama için, ölçek faktörünü veya z faktörünü kullanırsınız, ikisini birden değil (QGIS'de daha mantıklı bir seçim gibi görünür ve daha iyi çalışır). 2. Kuzey-güney ölçek faktörü her zaman 111320olmakla birlikte, doğu-batı ölçek faktörüdür 111320*cos(latitude*pi/180). Yüksek enlemlerde bu fark önemli sorunlar yaratacaktır ve muhtemelen hesaplamadan önce bir metre koordinat sistemine (örneğin, UTM) yeniden projeksiyon yapmanız gerekir. Düşük enlemlerde kısmında, bölme farkı deneyebilirsiniz: s = 111320 * 0.5 * (1 + cos(latitude*pi/180)).

— Matthias Fripp

4

DEM'lerle üniteleri ondalık derece ve dikey birimleri metre cinsinden uzun koordinatlarla çalışırken "ölçeği (dikey birimlerin yatay birimlere oranı) 111120'ye ayarlamanız gerekir.

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

Hatta bu ayarın lantlon koordinatlı DEM'ler ve feet cinsinden dikey birimler için ölçeğin 370400 olarak ayarlanması gerektiğini önermektedir.

Bu ayarlar eğim hesaplamalarını da etkiler.

— Gerardo Jimenez
kaynak