LAS ポイント統計をラスターに出力 (LAS Point Statistics As Raster) (データ管理)

サマリー

LAS ポイントに関する統計情報を反映するセル値を持つラスターを作成します。

使用法

  • LAS データセット レイヤーのポイント フィルターにより、分類コード、分類フラグ、および戻り値の組み合わせを使用して、このツールで処理されるポイントをフィルター処理できます。 フィルターは、[レイヤー プロパティ] ダイアログ ボックスで、または [LAS データセット レイヤーの作成 (Make LAS Dataset Layer)] ツールを使用して定義できます。 LAS ポイント フィルターは、すべての最後のリターン ポイントを評価するパルス数の方法を除くすべての統計出力で得られた結果に影響します。

  • [方法] パラメーターの [ポイント数] オプションは、LAS データセットのポイント密度を評価します。 このタイプのデータは、他のツールを使用して LAS ファイルを処理するときにパフォーマンスが低下する可能性があるため、高密度ポイントの不規則なホット スポットを検出する場合にこの情報が役立ちます。 高密度のホット スポットが発生した場合は、[LAS の間引き (Thin LAS)] ツールを使用して、より一貫した空間分布を持つ新しい点群を生成することができます。

  • [方法] パラメーターの [クラス コードの最頻値] を使用すると、特にそれらのオブジェクトが水平方向に空間的に分離されている場合に特定のクラス コードで定義された一意のオブジェクトの数を割り出せます。 そのためには、最初にポイント フィルターを適用して、特定のクラス コードを分離します。 次に、ポイントを表すセルの接続性が同じオブジェクトに属し、異なるオブジェクトのポイント間のギャップが維持されるように、十分な大きさのセル サイズを使用して統計ラスターを作成します。 たとえば、街灯に属するポイントの接続に適したセル サイズを特定するには、ポイントの全体的な間隔と隣接する街灯間の距離を考慮します。 結果として得られたラスターを使用して、各オブジェクトのフットプリントを評価することができます。 一意のオブジェクトの数を得るには、[ラスター → ポリゴン (Raster to Polygon)] ツールを使用してラスターをポリゴン フィーチャに変換するか、[リージョン グループ (Region Group)] ツールを使用してセルの各クラスターに個別値を割り当てます。

パラメーター

ラベル説明データ タイプ
入力 LAS データセット

処理対象の LAS データセット。

LAS Dataset Layer
出力ラスター

出力ラスターの場所と名前。 ジオデータベースまたは Esri Grid などのフォルダーにラスター データセットを格納する場合、ラスター データセットの名前にファイル拡張子は付けません。 ファイル拡張子を指定すると、フォルダーに格納する際のラスターの形式を定義できます。たとえば、*.tif を指定すると GeoTIFF が生成され、*.img を指定すると ERDAS IMAGINE 形式のファイルが生成されます。

ラスターを TIFF ファイルとして、またはジオデータベースに格納する場合は、ジオプロセシング環境設定を使用して、そのラスター圧縮のタイプと品質を指定できます。

Raster Dataset
方法
(オプション)

出力ラスターの各セル内の LAS ポイントに関して収集される統計情報のタイプを指定します。

  • パルス数ラスト リターン ポイントの数が収集されます。
  • ポイント数すべてのリターンのポイント数が収集されます。
  • ラスト リターンの最頻値ラスト リターン値の最頻値が収集されます。
  • クラス コードの最頻値クラス コードの最頻値が収集されます。
  • 強度値の範囲強度値の範囲が収集されます。
  • 標高値の範囲標高値の範囲が収集されます。
String
サンプリング タイプ
(オプション)

出力ラスターの解像度を定義するために、[サンプリング値] パラメーター値の解釈に使用する方法を指定します。

  • 観測点LAS データセット範囲の最も長い辺を分割するセル数を使用します。
  • セル サイズ出力ラスターのセル サイズを使用します。 これがデフォルトです。
String
サンプリング値
(オプション)

[サンプリング タイプ] パラメーターと組み合わせて、出力ラスターの解像度を定義するために使用される値。

Double

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
名前説明データ タイプ
in_las_dataset

処理対象の LAS データセット。

LAS Dataset Layer
out_raster

出力ラスターの場所と名前。 ジオデータベースまたは Esri Grid などのフォルダーにラスター データセットを格納する場合、ラスター データセットの名前にファイル拡張子は付けません。 ファイル拡張子を指定すると、フォルダーに格納する際のラスターの形式を定義できます。たとえば、*.tif を指定すると GeoTIFF が生成され、*.img を指定すると ERDAS IMAGINE 形式のファイルが生成されます。

ラスターを TIFF ファイルとして、またはジオデータベースに格納する場合は、ジオプロセシング環境設定を使用して、そのラスター圧縮のタイプと品質を指定できます。

Raster Dataset
method
(オプション)

出力ラスターの各セル内の LAS ポイントに関して収集される統計情報のタイプを指定します。

  • PULSE_COUNTラスト リターン ポイントの数が収集されます。
  • POINT_COUNTすべてのリターンのポイント数が収集されます。
  • PREDOMINANT_LAST_RETURNラスト リターン値の最頻値が収集されます。
  • PREDOMINANT_CLASSクラス コードの最頻値が収集されます。
  • INTENSITY_RANGE強度値の範囲が収集されます。
  • Z_RANGE標高値の範囲が収集されます。
String
sampling_type
(オプション)

出力ラスターの解像度を定義するために、[サンプリング値] パラメーター値の解釈に使用する方法を指定します。

  • OBSERVATIONSLAS データセット範囲の最も長い辺を分割するセル数を使用します。
  • CELLSIZE出力ラスターのセル サイズを使用します。 これがデフォルトです。
String
sampling_value
(オプション)

[サンプリング タイプ] パラメーターと組み合わせて、出力ラスターの解像度を定義するために使用される値。

Double

コードのサンプル

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

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

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
LasPointStatsAsRaster の例 2 (スタンドアロン スクリプト)

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

'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for
             each pulse in a given set of LAS files.
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)

# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
                                       "PREDOMINANT_RETURNS_PER_PULSE",
                                       "CELLSIZE", 15)

ライセンス情報

  • Basic: 次のものが必要 3D Analyst または Spatial Analyst
  • Standard: Yes
  • Advanced: Yes

関連トピック