LAS の高さメトリクス (3D Analyst ツール)

サマリー

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

使用法

このツールは、森林樹冠をより深く理解するために、地表より上にあるポイントの統計派生物を生成します。クラスコード値 0 または 1 (未分類のポイント)、3、4、5 (植生ポイント) が割り当てられた LAS ポイントが解析されます。
処理される LAS データは、このツールを実行する前に地表ポイントが分類されている必要があります。 LAS の地表分類 ツールを使用して、この要件を満たすことができます。
このツールは、指定された統計オプションと高さパーセンタイルごとに 1 つのラスターを出力します。各ラスターの名前は、指定されたベース名と高さ指標または空間統計で始まります。平均ラスターの場合、ファイル名に mean が付加されます。標準偏差の場合、ファイル名に stdev が追加されます。中央絶対偏差の場合、ファイル名に mad が追加されます。
特に地表ポイントが十分に分類されていない場合、極端に低いポイントを除外するには 最小高さ パラメーターを使用してください。これらは、このツールの結果を歪める可能性があるためです。
有意義な結果を得るため、地表より上にあるポイントが少なすぎるエリアを除外するには、最小ポイント数 パラメーターを使用してください。

ラベル	説明	データタイプ
入力 LAS データセット	処理される LAS データセット。	LAS Dataset Layer
出力場所	出力ラスターデータセットが保存されるフォルダーまたはジオデータベース。出力場所がフォルダーの場合、生成されるラスターデータセットは TIFF 形式になります。	Workspace
出力ベース名 (オプション)	出力ラスターデータセットのベース名。	String
統計オプション (オプション)	出力ラスターの各セルのエリア内にある、地上の未分類ポイントと植生ポイントに対して計算される統計を指定します。 Mean—The average height of the LAS points. Kurtosis—The sharpness of the change in the height of the LAS points. Skewness—The direction of deviation from the nominal height of the LAS points, which indicates the level and direction of asymmetry. Standard Deviation—The variation of the height of the points. Median Absolute Deviation—The median value of the deviation from the median height.	String
高さのパーセンタイル (オプション)	セル内で指定したパーセンテージのポイントが下にある高さ。たとえば、値が 95 の場合、生成されるセル値は、地上のポイントの 95 パーセントが存在する高さを示します。	Long
最小高さ	評価対象のポイントの地表面からの最小高さ。	Linear Unit
最小ポイント数	高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定した最小数よりもポイント数が少ないセルでは、出力にデータが含まれません。	Long
セルサイズ	出力ラスターデータセットのセルサイズ。	Linear Unit
ラスターフォーマット (オプション)	出力場所がフォルダーの場合に作成されるラスター形式を指定します。 GeoTiff—Output will be created in the GeoTIFF format. This is the default. Erdas Imagine (IMG)—Output will be created in the ERDAS IMAGINE format. Esri Grid—Output will be created in the Esri Grid format.	String

派生した出力

ラベル	説明	データタイプ
出力場所	出力ラスターデータセットの場所。	Workspace
出力ラスター	出力ラスターデータセット。	Raster Dataset

arcpy.3d.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`—The average height of the LAS points. `KURTOSIS`—The sharpness of the change in the height of the LAS points. `SKEWNESS`—The direction of deviation from the nominal height of the LAS points, which indicates the level and direction of asymmetry. `STANDARD_DEVIATION`—The variation of the height of the points. `MEDIAN_ABSOLUTE_DEVIATION`—The median value of the deviation from the median height.	String
height_percentiles [height_percentiles,...] (オプション)	セル内で指定したパーセンテージのポイントが下にある高さ。たとえば、値が 95 の場合、生成されるセル値は、地上のポイントの 95 パーセントが存在する高さを示します。	Long
min_height	評価対象のポイントの地表面からの最小高さ。	Linear Unit
min_points	高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定した最小数よりもポイント数が少ないセルでは、出力にデータが含まれません。	Long
cell_size	出力ラスターデータセットのセルサイズ。	Linear Unit
raster_format (オプション)	出力場所がフォルダーの場合に作成されるラスター形式を指定します。 `TIFF`—Output will be created in the GeoTIFF format. This is the default. `IMG`—Output will be created in the ERDAS IMAGINE format. `ESRI_GRID`—Output will be created in the Esri Grid format.	String

派生した出力

名前	説明	データタイプ
derived_out_location	出力ラスターデータセットの場所。	Workspace
output_rasters	出力ラスターデータセット。	Raster Dataset

コードのサンプル

LasHeightMetrics の例 1 (Python ウィンドウ)

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

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

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

LasHeightMetrics の例 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())

環境

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

ライセンス情報

Basic: 3D Analyst が必要
Standard: 3D Analyst が必要
Advanced: 3D Analyst が必要

LAS の高さメトリクス (3D Analyst ツール)

サマリー

使用法

パラメーター

派生した出力

環境

ライセンス情報

関連トピック