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
統計情報オプション
(オプション)

出力ラスターの各セルのエリア内にある、地上の未分類の植生ポイントに対して計算された統計情報を指定します。

  • 平均値LAS ポイントの平均高さ。
  • 尖度LAS ポイントの高さの変化の大きさ。
  • 歪度LAS ポイントの標準高さからの歪みの方向。非対称の度合いと方向を示します。
  • 標準偏差ポイントの高さの変量。
  • 中央絶対偏差中央高さからの偏差値の中央値。
String
高さのパーセンタイル
(オプション)

セル内で指定したパーセンテージのポイントが下にある高さ。たとえば、値が 95 の場合、生成されるセル値は、地上のポイントの 95 パーセントが存在する高さを示します。

Long
最小高さ

評価対象のポイントの地表面からの最小高さ。

Linear Unit
最小ポイント数

高さメトリクスを計算するために、指定のセルに存在する必要がある最小ポイント数。指定された最小数よりもポイント数が少ないセルには、出力にデータが含まれません。

Long
セル サイズ

出力ラスター データセットのセル サイズ。

Linear Unit
ラスター フォーマット
(オプション)

出力場所がフォルダーの場合に作成されるラスター形式を指定します。

  • GeoTiff出力が GeoTIFF 形式で作成されます。これがデフォルトです。
  • Erdas Imagine (IMG)出力が ERDAS IMAGINE 形式で作成されます。
  • Esri Grid出力が Esri Grid 形式で作成されます。
String

派生した出力

ラベル説明データ タイプ
出力場所

出力ラスター データセットの場所。

ワークスペース
出力ラスター

出力ラスター データセット

ラスター データセット

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,...]
(オプション)

出力ラスターの各セルのエリア内にある、地上の未分類の植生ポイントに対して計算された統計情報を指定します。

  • MEANLAS ポイントの平均高さ。
  • KURTOSISLAS ポイントの高さの変化の大きさ。
  • SKEWNESSLAS ポイントの標準高さからの歪みの方向。非対称の度合いと方向を示します。
  • 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())

ライセンス情報

  • Basic: 次のものが必要 3D Analyst
  • Standard: 次のものが必要 3D Analyst
  • Advanced: 次のものが必要 3D Analyst

関連トピック