Spatial Analyst のライセンスで利用可能。
3D Analyst のライセンスで利用可能。
概要
ラスター サーフェスの曲率を計算します。オプションとして断面曲率と平面曲率も計算できます。
使用法
主たる出力は、そのセルと周辺の 8 つのセルを対象として計算したセル単位でのサーフェスの曲率です。曲率は、サーフェスの二次導関数で、傾斜角の傾斜角です。オプションの出力曲率は 2 種類あります。断面曲率は最大傾斜角の方向で、平面曲率は最大傾斜角の方向に対して鉛直です。
正の曲率は、サーフェスがそのセルで上方向に凸状であることを示します。負の曲率は、サーフェスがそのセルで上方向に凹状であることを示します。ゼロ (0) という値はサーフェスが平らであることを示します。
断面出力の場合、負の値は、サーフェスがそのセルで上方向に凸状であることを示します。正の断面曲率は、サーフェスがそのセルで上方向に凹状であることを示します。ゼロ (0) という値はサーフェスが平らであることを示します。
平面出力の場合、正の値は、サーフェスがそのセルで上方向に凸状であることを示します。負の平面曲率は、サーフェスがそのセルで上方向に凹状であることを示します。ゼロ (0) という値はサーフェスが平らであることを示します。
曲率出力ラスターの単位は、オプションの出力断面曲率ラスターと出力平面曲率ラスターの単位も、Z 単位の 1/100 です。3 つの出力ラスターすべてに合理的に期待される値は、丘陵地 (中程度の起伏) で -0.5 〜 0.5、急峻で険しい山岳地帯 (極端な起伏) で -4 〜 4 の範囲になります。この範囲を超えるラスター サーフェスを使用することもできます。
入力ラスターをリサンプリングする必要がある場合は、共一次内挿法を使用します。入力ラスをリサンプリングする例としては、出力座標系、範囲、またはセル サイズが入力と異なる場合があります。
構文
arcpy.ddd.Curvature(in_raster, out_curvature_raster, {z_factor}, {out_profile_curve_raster}, {out_plan_curve_raster})| パラメーター | 説明 | データ タイプ |
in_raster | 入力サーフェス ラスター。 | Raster Layer |
out_curvature_raster | 出力曲率ラスター。 浮動小数点タイプになります。 | Raster Dataset |
z_factor (オプション) | サーフェス Z 単位あたりの地表 X,Y 単位の数。 z 単位の測定基準が入力サーフェスの xy 単位と異なる場合、z 係数を使用して z 測定単位を調整します。最終的な出力サーフェスを計算する際、入力サーフェスの Z 値に Z 係数が乗算されます。 X,Y 単位と Z 単位の測定基準が同じ場合、Z 係数は 1 これがデフォルトです。 X,Y 単位と Z 単位の測定基準が異なる場合は、Z 係数に適切な値を設定する必要があります。Z 係数が適切でないと、正しい結果が得られません。たとえば、Z 単位がフィートであり、X,Y 単位がメートルの場合、Z 係数 0.3048 を使用して Z 単位をフィートからメートルへ変換します (1 フィート = 0.3048 メートル)。 | Double |
out_profile_curve_raster (オプション) | 出力断面曲率ラスター データセットです。 これは、傾斜角の方向におけるサーフェスの曲率です。 浮動小数点タイプになります。 | Raster Dataset |
out_plan_curve_raster (オプション) | 出力平面曲率ラスター データセットです。 これは、傾斜角の方向に鉛直なサーフェスの曲率です。 浮動小数点タイプになります。 | Raster Dataset |
コードのサンプル
Curvature (曲率) の例 1 (Python ウィンドウ)
この例では、入力サーフェス ラスターから曲率ラスターを作成し、Z 係数も適用します。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Curvature_3d("elevation", "c:/output/outcurv01", 1.094)
Curvature (曲率) の例 2 (スタンドアロン スクリプト)
この例では、入力サーフェス ラスターから曲率ラスターを作成し、Z 係数も適用します。
# Name: Curvature_3d_Ex_02.py
# Description: Calculates the curvature of a raster surface,
# optionally including profile and plan curvature.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inRaster = "elevation"
outRaster = "C:/output/outcurv02"
zFactor = 1.094
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute Curvature
arcpy.Curvature_3d(inRaster, outRaster, 1.094)
環境
ライセンス情報
- Basic: 次のものが必要 3D Analyst または Spatial Analyst
- Standard: 次のものが必要 3D Analyst または Spatial Analyst
- Advanced: 次のものが必要 3D Analyst または Spatial Analyst