スケールとZファクターは、QGISの陰影起伏分析に影響を与えません

10

SRTM DEMを使用して、QGIS 2.2.0で陰影起伏を作成しました。ただし、結果の陰影起伏は垂直方向に誇張されているようです。また、スケールとZファクターの設定を変更しても、結果の陰影起伏に影響はありません。私は、QGISの両方のオプションを試してみました。分析-> DEM（地形モデル）と地形分析->陰影起伏。どちらの方法でも同じ結果が得られました。

これは孤立したインシデント、バグですか、それとも何か不足していますか？

qgis scale hillshade

— PyMapr
ソース

4

latlon座標が10進度で垂直単位がメートルであるDEMを使用している場合は、スケール（垂直単位と水平単位の比率）を111120に設定する必要があります。AsterのDEMをいくつか試しましたが、問題はありませんでした。この設定は、陰影起伏だけでなく、勾配の計算にも適用されます。それがお役に立てば幸い

— ヘラルドヒメネス

@GerardoJimenez、私は試してみましたが、うまくいきました。ありがとう

— PyMapr 2014年

@ GerardoJimenez、Thanderboltのどちらかがこの質問に対して送信および回答し、Thunderboltがそれを受け入れて、この質問を未回答から削除する必要があります。

— SS_Rebelious 2014年

12

@ジェラルドの答えは赤道でのみ正しいです。gdaldemドキュメントは、これらのスケール値をためているノート「赤道近くLatLongブログ予測」

メートル単位の垂直単位のあるデムの他の緯度でのスケールまたはZ係数（スケールの逆）を計算するには、次の方程式を使用できます。

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

ここで：111320は赤道での1度の長さ（メートル単位）（WGS84に基づく、覚えやすいか、履歴を理解している場合は111111を使用できます）、緯度はDEMの中心の緯度です。

Excel / Libre Office Calcで次の数式を使用して、scale / z係数をすばやく計算できます。

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

たとえば、緯度= 38.5SのスケールとZ係数

垂直単位がフィートのDEMの場合、赤道での1度の長さをメートル単位でフィートに変換します。つまり、3.28 * 111320 = 365130です。したがって、方程式は次のようになります。

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

または、DEMを投影座標系に再投影することもできます（gdaldem docsでも提案されています）。

— ユーザー2856
ソース

この答えは今日私に本当に役立ちましたが、QGISで使用されている用語について少し混乱しています。Raster-Terrain Analysis-Slopeに移動すると、「Zファクター」というフィールドが表示されますが、ここでは「Zファクター」と呼ばれる逆の値ではなく、「スケール」と呼ばれる値をここに入力する必要があります。上記の計算で逆数を計算するときに取得する非常に小さい数値を入力することはできません。デフォルト値は、それが大きい数値でなければならないことをすでに示していますが、「スケール」と「Zファクター」の両方の用語を考えたため、まだ少し混乱しています。同義語として使用できませんでした。

— Taya

@Taya QGIS 2.14のメニューから、Raster-> Analysis-> DEM（Terrain Models）ツールに「Z Factor」フィールドが表示されません。スケールフィールドのみが表示されます。おそらくそれはQGISの以前のバージョンのバグでした...？編集：ああ、あなたは地形解析プラグインを使用していますが、それはプラグインのバグであるに違いありません。

— user2856

この回答にはいくつか変更を加えることをお勧めします。1.ほとんどのアプリケーションでは、スケールファクターまたはzファクターのいずれかを使用します。2.南北の倍率は常に111320ですが、東西の倍率は111320*cos(latitude*pi/180)です。高緯度では、この違いにより重大な問題が発生します。計算する前に、メートル座標系（UTMなど）に再投影する必要があります。低緯度では、違いを分割してみてくださいs = 111320 * 0.5 * (1 + cos(latitude*pi/180))。

— Matthias Fripp

4

緯度と経度の座標でDEMを扱い、単位は10進度、垂直単位はメートルです。「スケール（垂直単位と水平単位の比率）」を111120に設定する必要があります。

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

彼らはさらに、ラントロン座標と垂直単位がフィートのDEMのこの設定では、スケールを370400に設定する必要があることを示唆しています。

この設定は、勾配の計算にも影響します。

— ヘラルドヒメネス
ソース