サマリー
LAS データで捕捉した植生ポイントの標高計測値の分布に関する統計情報を計算します。
使用法
森林面積から、LIDAR コレクションのバイオマスを計算するときに、このツールの使用を検討してください。
このツールでは地盤分類ポイントが必要で、クラス コード値が 0、1、3、4、5 の LAS ポイントを処理します。
構文
LasHeightMetrics_3d (in_las_dataset, out_location, {base_name}, {statistics}, {height_percentiles}, min_height, min_points, cell_size)| パラメーター | 説明 | データ タイプ |
in_las_dataset | 処理対象の LAS データセット。 | LAS Dataset Layer |
out_location | 出力ラスター データセットが保存されるフォルダーまたはジオデータベース | Workspace |
base_name (オプション) | 出力ラスター データセットのベース名 | String |
statistics [statistics,...] (オプション) | 出力ラスターの各セルのエリア内にある、地上の未分類の植生ポイントに対して計算された統計情報
| String |
height_percentiles [height_percentiles,...] (オプション) | 地表面からの高さのパーセンタイル。これに対してヒストグラムを計算します。 | Long |
min_height | 評価対象のポイントの地表面からの最小高さ | Linear Unit |
min_points | LAS 高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定された最小数よりもポイント数が少ないセルには、出力にデータが含まれません。 | Long |
cell_size | 出力ラスター データセットのセル サイズ | Linear Unit |
派生した出力
| 名前 | 説明 | データ タイプ |
| derived_out_location | 出力ラスター データセットの場所 | ワークスペース |
| output_rasters | 出力ラスター データセット |
コードのサンプル
LasHeightMetrics (LAS の高さメトリクス) の例 1 (Python ウィンドウ)
次のサンプルは、Python ウィンドウでこのツールを使用する方法を示しています。
arcpy.env.workspace = 'C:/data'
arcpy.ClassifyLasBuilding_3d('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())
環境
ライセンス情報
- ArcGIS Desktop Basic: 次のものが必要 3D Analyst
- ArcGIS Desktop Standard: 次のものが必要 3D Analyst
- ArcGIS Desktop Advanced: 次のものが必要 3D Analyst