LAS データセットの統計情報（LAS Dataset Statistics） (データの管理)
サマリ
LAS データセットの統計情報を計算または更新し、オプションの統計情報のレポートを作成します。
使用法
統計情報は、ファイルの小さな領域を漸次評価するラスタ化という方法で、各 LAS ファイルの平均ポイント間隔を計算し、推定を決定します。
統計情報は、LAS データセット レイヤのフィルタ オプションを有効にして、利用可能なクラス コードを自動的に表示し、参照される LAS ファイルにある値を返します。フィルタ オプションは、ArcMap および ArcScene の [レイヤ プロパティ] ダイアログ ボックスで指定できます。
-
LAS ファイルに関する次の情報は、LAS データセットに格納され、ArcCatalog の [LAS データセット プロパティ] ダイアログ ボックスで確認することができます。
- 標高と強度値の範囲
- クラスコードあたりのデータポイントの数と返す値
- RGB 値の存在
- LAS バージョンとポイントの形式 ID
- プロジェクト ID とファイル ソース ID
- LIDAR ファイルの生成に使用されるソフトウェアの名前
オプションの統計情報のレポート ファイルは、LAS データセットのすべてのファイル、または LAS ファイルごとの LAS プロパティの概要を、任意の数のアプリケーションにインポート可能なテキスト形式で提供します。
構文
| パラメータ | 説明 | データ タイプ |
in_las_dataset |
入力 LAS データセット。 | LAS Dataset Layer |
calculation_type | 統計情報の計算をすべての LIDAR ファイルに対して行うか、あるいは統計情報を持たない LIDAR ファイルに対してのみ行うかを指定します。
| Boolean |
out_file (オプション) |
LAS データセットの統計情報のサマリを含む出力テキスト ファイル | Text File |
summary_level (オプション) |
レポートに含まれるサマリのタイプを指定します。
| String |
delimiter (オプション) |
テキスト ファイルで使用されるフィールドの区切り記号。
| String |
decimal_separator (オプション) | 整数部と小数部を区別するためにテキスト ファイルで使用される小数記号。
| String |
コードのサンプル
次のサンプルは、Python ウィンドウでこのツールを使用する方法を示しています。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.LASDatasetStatistics_3d("test.lasd", "NO_FORCE", "LAS_FILE",
"DECIMAL_POINT", "SPACE", "LAS_summary.txt")
次のサンプルは、スタンドアロン Python スクリプトでこのツールを使用する方法を示しています。
'''*********************************************************************
Name: Modify Files in LAS Dataset& Calculate Stats for LASD
Description: Adds files & surface constraints to a LAS dataset, then
calculates statistics and generates report.
*********************************************************************'''
# Import system modules
import arcpy
from arcpy import env
import exceptions, sys, traceback
try:
# Script variables
env.workspace = 'C:/data'
lasd = 'sample.lasd'
oldLas = ['2006', '2007/file2.las']
newLas = ['2007_updates_1', '2007_updates_2']
oldSurfaceConstraints = ['boundary.shp', 'streams.shp']
newSurfaceConstraints = [['sample.gdb/boundary', '<None>',
'Soft_Clip']
['sample.gdb/streams', 'Shape.Z',
'Hard_Line']]
arcpy.management.RemoveFilesFromLasDataset(lasd, oldLas,
oldSurfaceConstraints)
arcpy.management.AddFilesToLasDataset(lasd, newLas, 'RECURSION',
newSurfaceConstraints)
arcpy.management.LasDatasetStatistics(lasd, "UPDATED_FILES",
"lasd_stats.txt",
"LAS_FILE", "DECIMAL_POINT",
"SPACE", "LAS_summary.txt")
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)