LAS の高さメトリクス (LAS Height Metrics) (3D Analyst)

サマリー

LAS データで捕捉した植生ポイントの標高計測値の分布に関する統計情報を計算します。

使用法

• このツールは、森林樹冠の理解を深めることができる、地表より上のポイントの統計的派生物を生成します。未分類ポイントを表すクラス コード値 0 または 1 と、植生ポイントを表すクラス コード値 3、4、5 が割り当てられた LAS ポイントが解析されます。

• このツールを実行する前に、処理される LAS データの地表ポイントを分類する必要があります。この要件を満たすために、[LAS の地表分類 (Classify LAS Ground)] ツールを使用できます。

• このツールは、指定された統計情報オプションと高さのパーセンタイルごとに 1 つのラスターを出力します。各ラスターの名前は、指定されたベース名で始まり、高さメトリクスまたは空間統計が付きます。平均値ラスターの場合、ファイル名に「mean」という単語が追加されます。標準偏差の場合、ファイル名に「stdev」が追加されます。平均絶対偏差の場合、ファイル名に「mad」が追加されます。

• [最小高さ] パラメーターを使用すると、特に地表ポイントが十分に分類されていない場合に、このツールの結果に歪みを生じさせる可能性がある非常に低いポイントを除外できます。

• [最小ポイント数] パラメーターを使用すると、地表より上のポイントが少なすぎるエリアを除外して、意味のある結果を得ることができます。

パラメーター

ダイアログ

ラベル	説明	データ タイプ
入力 LAS データセット	処理される LAS データセット。	LAS Dataset Layer
出力場所	出力ラスター データセットが保存されるフォルダーまたはジオデータベース 出力場所がフォルダーの場合、生成されるラスター データセットは TIFF 形式になります。	Workspace
出力ベース名 (オプション)	出力ラスター データセットのベース名	String
統計情報オプション (オプション)	出力ラスターの各セルのエリア内にある、地上の未分類の植生ポイントに対して計算された統計情報を指定します。 MEAN—LAS ポイントの平均高さ。 KURTOSIS—LAS ポイントの高さの変化の大きさ。 SKEWNESS—LAS ポイントの標準高さからの歪みの方向。非対称の度合いと方向を示します。 STANDARD_DEVIATION—ポイントの高さの変量。 MEDIAN_ABSOLUTE_DEVIATION—中央高さからの偏差値の中央値。	String
高さのパーセンタイル (オプション)	セル内で指定したパーセンテージのポイントが下にある高さ。たとえば、値が 95 の場合、生成されるセル値は、地上のポイントの 95 パーセントが存在する高さを示します。	Long
最小高さ	評価対象のポイントの地表面からの最小高さ。	Linear Unit
最小ポイント数	高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定された最小数よりもポイント数が少ないセルには、出力にデータが含まれません。	Long
セル サイズ	出力ラスター データセットのセル サイズ。	Linear Unit
ラスター フォーマット (オプション)	出力場所がフォルダーの場合に作成されるラスター形式を指定します。 TIFF—出力が GeoTIFF 形式で作成されます。これがデフォルトです。 IMG—出力が ERDAS IMAGINE 形式で作成されます。 ESRI_GRID—出力が Esri Grid 形式で作成されます。	String

派生した出力

ラベル	説明	データ タイプ
出力場所	出力ラスター データセットの場所。	ワークスペース
出力ラスター	出力ラスター データセット	ラスター データセット

Python

arcpy.ddd.LasHeightMetrics(in_las_dataset, out_location, {base_name}, {statistics}, {height_percentiles}, min_height, min_points, cell_size, {raster_format})

名前	説明	データ タイプ
in_las_dataset	処理される LAS データセット。	LAS Dataset Layer
out_location	出力ラスター データセットが保存されるフォルダーまたはジオデータベース 出力場所がフォルダーの場合、生成されるラスター データセットは TIFF 形式になります。	Workspace
base_name (オプション)	出力ラスター データセットのベース名	String
statistics [statistics,...] (オプション)	出力ラスターの各セルのエリア内にある、地上の未分類の植生ポイントに対して計算された統計情報を指定します。 MEAN—LAS ポイントの平均高さ。 KURTOSIS—LAS ポイントの高さの変化の大きさ。 SKEWNESS—LAS ポイントの標準高さからの歪みの方向。非対称の度合いと方向を示します。 STANDARD_DEVIATION—ポイントの高さの変量。 MEDIAN_ABSOLUTE_DEVIATION—中央高さからの偏差値の中央値。	String
height_percentiles [height_percentiles,...] (オプション)	セル内で指定したパーセンテージのポイントが下にある高さ。たとえば、値が 95 の場合、生成されるセル値は、地上のポイントの 95 パーセントが存在する高さを示します。	Long
min_height	評価対象のポイントの地表面からの最小高さ。	Linear Unit
min_points	高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定された最小数よりもポイント数が少ないセルには、出力にデータが含まれません。	Long
cell_size	出力ラスター データセットのセル サイズ。	Linear Unit
raster_format (オプション)	出力場所がフォルダーの場合に作成されるラスター形式を指定します。 TIFF—出力が GeoTIFF 形式で作成されます。これがデフォルトです。 IMG—出力が ERDAS IMAGINE 形式で作成されます。 ESRI_GRID—出力が Esri Grid 形式で作成されます。	String

派生した出力

名前	説明	データ タイプ
derived_out_location	出力ラスター データセットの場所。	ワークスペース
output_rasters	出力ラスター データセット	ラスター データセット

コードのサンプル

LasHeightMetrics (LAS の高さメトリクス) の例 1 (Python ウィンドウ)

次のサンプルは、Python ウィンドウでこのツールを使用する方法を示しています。

arcpy.env.workspace = 'C:/data'

arcpy.ddd.ClassifyLasBuilding('Highland.lasd', minHeight='9 feet', 
                              minArea='30 Square Feet', compute_stats=True)

LasHeightMetrics (LAS の高さメトリクス) の例 2 (スタンドアロン スクリプト)

次のサンプルは、スタンドアロン Python スクリプトでこのツールを使用する方法を示しています。

'''****************************************************************************
       Name: Tile & Classify LAS Files
Description: Creates & classifies tiled LAS files.
****************************************************************************'''
# Import system modules
import arcpy
import tempfile
import math

in_las = arcpy.GetParameterAsText(1) # The LAS files that need to be tiled
out_folder = arcpy.GetParameterAsText(2) # folder for LAS files
basename = arcpy.GetParameterAsText(3) # basename for output files
out_lasd = arcpy.GetParameterAsText(4) # output LAS dataset


try:
    # Create temp LAS dataset to reference LAS files that will be tiled
    temp_lasd = arcpy.CreateUniqueName('temp.lasd', tempfile.gettempdir())
    arcpy.management.CreateLasDataset(in_las, temp_lasd)
    arcpy.ddd.TileLas(temp_lasd, out_folder, basename, out_lasd, las_version=1.4, 
                      point_format=7, file_size=300)
    arcpy.management.Delete(temp_lasd)
    arcpy.ddd.ClassifyLasGround(out_lasd, method='AGGRESSIVE')
    arcpy.ddd.ClassifyLasBuilding(out_lasd, min_height='3 Meters', min_area='4 Meters')
    arcpy.ddd.ClassifyLasByHeight(out_lasd, height_classification=[(3, 6), (4,20), (5,70)],
                                  noise='All Noise', compute_stats='COMPUTE_STATS')

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

環境

現在のワークスペース, 範囲, 出力座標系, 地理座標系変換