J'ai utilisé un SRTM DEM pour créer un ombrage dans QGIS 2.2.0. Cependant, l'ombrage résultant semble être exagéré verticalement. De plus, la modification des paramètres d'échelle et de facteur z n'a aucun effet sur l'ombrage résultant. J'ai essayé les deux options dans QGIS, Analysis-> DEM (Modèles de terrain) et Analyse de terrain -> Hillshade. Les deux méthodes ont donné les mêmes résultats.
S'agit-il d'un incident isolé, d'un bug ou est-ce que je manque quelque chose?
Réponses:
@ La réponse de Gerardo n'est correcte qu'à l'équateur. Le doc gdaldem note que ces valeurs d'échelle sont pour les "projections LatLong près de l'équateur"
Pour calculer l'échelle ou le facteur Z (inverse de l'échelle) à d'autres latitudes pour un dem avec des unités verticales en mètres, vous pouvez utiliser les équations suivantes:
où: 111320 est la longueur d'un degré à l'équateur en mètres (basé sur WGS84, vous pouvez utiliser 111111 si c'est plus facile à retenir, ou si vous avez une appréciation de l'histoire ) et la latitude est la latitude au centre de votre DEM.
Vous pouvez calculer rapidement le facteur échelle / z dans Excel / Libre Office Calc avec les formules suivantes:
Vous pouvez même utiliser une simple recherche Google:
Par exemple, échelle et facteur z pour Latitude = 38,5S
Pour les DEM avec des unités verticales en pieds, il suffit de convertir la longueur d'un degré à l'équateur en mètres en pieds, soit 3,28 * 111320 = 365130. Vos équations deviennent donc:
Alternativement, vous pouvez reprojeter votre DEM à un système coordonné projeté (également suggéré par les documents gdaldem).
la source
111320
, mais le facteur d'échelle est-ouest l'est111320*cos(latitude*pi/180)
. Aux hautes latitudes, cette différence créera des problèmes importants et vous devrez probablement reprojeter dans un système de coordonnées métriques (par exemple, UTM) avant de calculer. Aux basses latitudes, vous pouvez essayer de diviser la différence:s = 111320 * 0.5 * (1 + cos(latitude*pi/180))
.Lorsque vous travaillez avec des DEM en coordonnées longues lat avec leurs unités en degrés décimaux et les unités verticales en mètres, vous devez définir "l'échelle (rapport des unités verticales aux unités horizontales) à 111120. Dans
http://www.gdal.org/gdaldem.html
ils suggèrent même que ce paramètre pour les DEM avec les coordonnées lantlon et les unités verticales en pieds, l'échelle doit être réglée sur 370400
Ce paramètre affecte également les calculs de pente.
la source