LAS ポイント統計をラスタに出力（LAS Point Statistics As Raster） (データの管理)
サマリ
LAS データセットが参照する LAS ファイルの計測値についての統計情報を反映したセル値を持つラスタを作成します。
使用法
植生の存在を示すような高いリターン値を持つ場所を識別するには、[方法] パラメータに PREDOMINANT_LAST_RETURN の使用を検討してください。
外れ値を持つ可能性がある場所を判定するには、[メソッド] パラメータに Z_RANGE の使用を検討してください。
LAS データセット レイヤを使用すると、LAS ポイントをクラス コードやリターンでフィルタ処理できます。レイヤを作成するには、[LAS データセット レイヤの作成（Make LAS Dataset Layer）] ツールを使用するか、ArcMap または ArcScene に LAS データセットを読み込んで、レイヤ プロパティ ダイアログ ボックスからクラス コードまたはリターンを指定します。
構文
| パラメータ | 説明 | データ タイプ |
in_las_dataset |
入力 LAS データセット。 | LAS Dataset Layer |
out_raster |
出力ラスタの場所と名前。ジオデータベースまたは Esri Grid などのフォルダにラスタ データセットを格納する場合、ラスタ データセットの名前にファイル拡張子は付けません。ラスタをフォルダに格納するときにその形式を定義するためにファイル拡張子を付けることはできます。
ラスタを TIFF ファイルとして、またはジオデータベースに格納する場合は、ジオプロセシング環境設定を使用して、そのラスタ圧縮のタイプと品質を指定できます。 | Raster Dataset |
method (オプション) |
出力ラスタの各セル内にある LAS ポイントについて収集した統計情報の種類
| String |
sampling_type (オプション) |
出力ラスタの解像度を定義するために、[サンプリング値] の解釈に使用する方法を指定します。
| String |
sampling_value (オプション) |
[サンプリング タイプ] と組み合わせて、出力ラスタの解像度を定義するために使用される値を指定します。 | Double |
コードのサンプル
次のサンプルは、Python ウィンドウでこのツールを使用する方法を示しています。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img",
"INTENSITY_RANGE", "CELLSIZE", 15)
次のサンプルは、スタンドアロン 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
import exceptions, sys, traceback
# 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)
try:
arcpy.CheckOutExtension('3D')
# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
"PREDOMINANT_RETURNS_PER_PULSE",
"CELLSIZE", 15)
arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
print arcpy.GetMessages()
except:
# Get the traceback object
tb = sys.exc_info()[2]
tbinfo = traceback.format_tb(tb)[0]
# Concatenate error information into message string
pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
.format(tbinfo, str(sys.exc_info()[1]))
msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
# Return python error messages for script tool or Python Window
arcpy.AddError(pymsg)
arcpy.AddError(msgs)