① 元データ
- 国土地理院の標高データ(GSI DEM)を利用します。
- 指定した座標の「標高(m)」を取得しています。
② 傾斜の計算方法(どう出してるか)
- 各セルの中心から、東・西・南・北の4点の標高を取得します。
- その差分から、横方向と縦方向の傾きを近似し、合成して傾斜角(度)を算出します。
👉 イメージ
東 − 西 → 横方向の傾き(dz/dx)
北 − 南 → 縦方向の傾き(dz/dy)
dz/dx ≈ (E − W) / cellM
dz/dy ≈ (N − S) / cellM
傾斜角(度) = atan( sqrt( (dz/dx)^2 + (dz/dy)^2 ) ) × 180/π
③ 表示
- 算出した傾斜角(度)を、必要に応じて %(勾配)に変換します。
- 色で可視化して、マップ上に表示します。
勾配%(percent grade) = tan(傾斜角[deg]) × 100
※計算はセルごとの近似です(東西南北の4点を使う有限差分)。標高が取得できない点がある場合は N/A として扱います。