LAS の抽出 (Extract LAS) (3D Analyst)

サマリー

LAS データセットまたは点群シーン レイヤーの点群データから新しい LAS ファイルを作成します。

図

使用法

• このツールは、入力 LAS データセットまたは点群シーン レイヤー パッケージ (*.slpk ファイル) から新しい *.las ファイルまたは *.zlas ファイルを作成します。 入力が LAS データセットの場合、LAS データセットによって参照される *.las ファイルまたは *.zlas ファイルごとに 1 つの出力ファイルが作成されます。 入力が点群シーン レイヤーの場合、抽出されるポイントの数に応じて出力ファイルの数も異なります。 各出力ファイルには約 1670 万ポイントというサイズ制限があります。 複数のファイルが作成された場合、出力されるファイルに重なる範囲が含まれる場合があります。これが原因で、特定の状況において不十分なデータ処理につながることがあります。 [LAS タイルの作成 (Tile LAS)] ツールを使用して、範囲が重なっていない特定のタイル スキーマに、点群シーン レイヤーから抽出されたファイルを再編成することを検討してください。

• 処理範囲、抽出境界、またはその両方を指定することで、入力点群のサブセットを抽出できます。 処理範囲と抽出境界を指定すると、両者の交差部分によって、抽出される *.las ファイルのカバレッジが定義されます。 処理範囲と抽出境界のどちらも定義していない場合、入力点群の範囲全体が出力 *.las ファイルに抽出されます。 また、ポイント フィルターを入力点群に適用して、分類コード、分類フラグ、リターン値の任意の組み合わせに基づいて、抽出されるポイントを制限できます。 アクティブなマップまたはシーンにレイヤーを読み込み、[レイヤー プロパティ] ダイアログ ボックスに希望のフィルターを指定することでポイント フィルターを定義します。 [LAS データセット レイヤーの作成 (Make LAS Dataset Layer)] ツールを使用して、LAS データセットのポイント フィルターを定義することもできます。

• [出力座標系] 環境を設定することで、このツールを使用して別の水平座標系と鉛直座標系に入力点群の空間参照を再投影できます。 *.las ファイルの空間参照は、ファイルのヘッダーか、WKT (既知の文字列) 規則を使用して補正 *.prj ファイルで定義できます。 *.prj ファイルは *.las ファイルと同じ場所に配置し、*.las ファイルと同じ名前にする必要があります。 *.prj ファイルが存在する場合は、*.las ファイルのヘッダーの空間参照情報はオーバーライドされます。 空間参照が存在しないか誤って定義されている *.las ファイルに対する *.prj ファイルを作成するには、個別の *.las または *.zlas ファイルに対して [投影法の定義 (Define Projection)] ツールを使用するか、[LAS データセットの作成 (Create LAS Dataset)] ツールを [PRJ] オプションを有効にした状態で使用します。

• ポイントの空間クラスタリングに対応しないバイナリ シーケンスの *.las ファイルに LAS ポイント レコードを保存することは珍しいことではありません。 このような分布のデータが検索される場合、LAS ポイントを表すバイナリ レコードへのアクセス効率が低下する可能性があります。 作成される *.las ファイルのポイントを再配置することで、視覚化やその他の空間操作のデータが最適化されます。 再配置オプションが有効化されると、統計情報が自動的に計算されます。 LAS ポイントの再配置をしないことを選択した場合は、統計情報の計算を有効化または無効化することができます。 統計情報の計算により、空間検索が最適化され、クラス コードのサマリーが提供され、*.las ファイルに表示される値が返されます。 ただし、このツールの処理に時間がかかるようになります。 作成される *.las ファイルが ArcGIS で使用されない場合、ツールの処理時間を短縮するために、統計情報の計算を無効化することができます。

パラメーター

ダイアログ

ラベル	説明	データ タイプ
入力点群データ	処理する LAS データセットまたは点群シーン レイヤー パッケージ (*.slpk ファイル)。	LAS Dataset Layer; Scene Layer; File
ターゲット フォルダー	出力 *.las ファイルが書き込まれる既存のフォルダー。	Folder
処理範囲 (オプション)	評価されるデータの範囲。 [デフォルト] - 範囲はすべての入力の最大範囲に基づきます。これがデフォルトです。 [入力データのすべての領域] - この範囲は、すべての入力の最大範囲に基づきます。 [入力データの共通領域] - この範囲は、すべての入力に共通する最小領域に基づきます。 [現在の表示範囲] - 範囲は、表示範囲と同じになります。アクティブなマップが存在しない場合、このオプションは使用できません。 [以下の指定に一致] - この範囲は、指定された最小および最大範囲値に基づきます。 [参照] - 範囲は、既存のデータセットに基づきます。	Extent
抽出境界 (オプション)	クリップされる *.las ファイルのエリアを定義するポリゴン境界。	Feature Layer
範囲と重なる LAS ファイル全体を処理 (オプション)	*.las ファイルの処理方法の決定に対象地域をどのように使用するかを指定します。 対象地域は、[処理範囲] パラメーター値、[処理境界] パラメーター値、またはこれら 2 つの組み合わせによって定義されます。 オフ - 対象地域と交差している LAS ポイントだけが処理されます。 これがデフォルトです。 オン - *.las ファイルの一部が対象地域と交差している場合、対象地域の外側にあるポイントも含め、その *.las ファイル内のすべてのポイントが処理されます。	Boolean
出力ファイル名の接尾辞 (オプション)	各出力 *.las ファイルの名前に追加されるテキスト。 各ファイルは、ソース ファイルからベース名を継承し、このパラメーターで指定された接尾辞が付加されます。	String
可変長レコードの削除 (オプション)	可変長レコード (VLR) を削除するかどうかを指定します。 各 *.las ファイルには、それを生成したソフトウェアによって追加された VLR のセットが含まれている可能性があります。 これらのレコードの意味は、通常は元のソフトウェアしか認知していません。 この情報を理解するアプリケーションによって出力 LAS データが処理されるのでなければ、この VLR を保持しても付加価値的な機能は何ももたらされません。 VLR を削除することで、そのサイズ合計やそれを含むファイルの数によっては、かなり大きなディスク容量を節約することができます。 オフ - 入力 *.las ファイル内の可変長レコードを削除せず、出力 *.las ファイル内に維持します。 これがデフォルトです。 オン - 入力 *.las ファイル内の可変長レコードを出力 *.las ファイルから削除します。	Boolean
ポイントの整理 (オプション)	*.las ファイルのポイントを整理するかどうかを指定します。 オフ - *.las ファイルのポイントの順序を整理しません。 オン - *.las ファイルのポイントを整理します。 これがデフォルトです。	Boolean
統計情報の計算 (オプション)	LAS データセットで参照されている *.las ファイルの統計情報を計算するかどうかを指定します。 統計情報を計算することで、*.las ファイルごとの空間インデックスが提供され、解析と表示のパフォーマンスが向上します。 また統計によって、分類コードやリターン情報などの LAS 属性の表示が *.las ファイルに存在する値に制限されるので、フィルタリングとシンボルのエクスペリエンスも強化されます。 オン - 統計情報を計算します。 これがデフォルトです。 オフ - 統計情報を計算しません。	Boolean
出力 LAS データセット (オプション)	新しく作成された *.las ファイルを参照する出力 LAS データセット。	LAS Dataset
圧縮 (オプション)	出力 *.las ファイルを圧縮形式にするか、標準の LAS 形式にするかを指定します。 入力と同様—入力と同じ圧縮になります。 このオプションは、入力が LAS データセットである場合にのみ使用できます。その場合はこれがデフォルトです。 圧縮なし—出力は標準の LAS 形式 (**.las) になります。 これは、入力が点群シーン レイヤーである場合のデフォルトです。 zLAS 圧縮—出力 *.las ファイルは zLAS 形式で圧縮されます。	String

派生した出力

ラベル	説明	データ タイプ
出力フォルダー	*.las ファイルが書き込まれるフォルダー。	Folder

Python

arcpy.ddd.ExtractLas(in_las_dataset, target_folder, {extent}, {boundary}, {process_entire_files}, {name_suffix}, {remove_vlr}, {rearrange_points}, {compute_stats}, {out_las_dataset}, {compression})

名前	説明	データ タイプ
in_las_dataset	処理する LAS データセットまたは点群シーン レイヤー パッケージ (*.slpk ファイル)。	LAS Dataset Layer; Scene Layer; File
target_folder	出力 *.las ファイルが書き込まれる既存のフォルダー。	Folder
extent (オプション)	評価されるデータの範囲。 MAXOF - すべての入力の最大範囲が使用されます。 MINOF - すべての入力に共通する最小領域が使用されます。 DISPLAY - 範囲は、表示範囲と同じになります。 [レイヤー名] - 指定したレイヤーの範囲が使用されます。 Extent オブジェクト - 指定したオブジェクトの範囲が使用されます。 [座標のスペース区切りの文字列] - 指定した文字列の範囲が使用されます。座標は、x-min、y-min、x-max、y-max の順序で表されます。	Extent
boundary (オプション)	クリップされる *.las ファイルのエリアを定義するポリゴン境界。	Feature Layer
process_entire_files (オプション)	処理範囲の適用方法を指定します。 PROCESS_EXTENT—対象地域と交差している LAS ポイントだけが処理されます。 これがデフォルトです。 PROCESS_ENTIRE_FILES—*.las ファイルの一部が対象地域と交差している場合、対象地域の外側にあるポイントも含め、その *.las ファイル内のすべてのポイントが処理されます。	Boolean
name_suffix (オプション)	各出力 *.las ファイルの名前に追加されるテキスト。 各ファイルは、ソース ファイルからベース名を継承し、このパラメーターで指定された接尾辞が付加されます。	String
remove_vlr (オプション)	可変長レコード (VLR) を削除するかどうかを指定します。 各 *.las ファイルには、それを生成したソフトウェアによって追加された VLR のセットが含まれている可能性があります。 これらのレコードの意味は、通常は元のソフトウェアしか認知していません。 この情報を理解するアプリケーションによって出力 LAS データが処理されるのでなければ、この VLR を保持しても付加価値的な機能は何ももたらされません。 VLR を削除することで、そのサイズ合計やそれを含むファイルの数によっては、かなり大きなディスク容量を節約することができます。 MAINTAIN_VLR—処理対象の *.las ファイルの可変長レコードに格納されているすべてのデータを削除せず、抽出された *.las ファイル内に維持します。 これがデフォルトです。 REMOVE_VLR—余分な可変長レコードを *.las ファイルから削除します。	Boolean
rearrange_points (オプション)	*.las ファイルのポイントを整理するかどうかを指定します。 MAINTAIN_POINTS—*.las ファイルのポイントの順序を整理しません。 REARRANGE_POINTS—*.las ファイルのポイントを整理します。 これがデフォルトです。	Boolean
compute_stats (オプション)	LAS データセットで参照されている *.las ファイルの統計情報を計算するかどうかを指定します。 統計情報を計算することで、*.las ファイルごとの空間インデックスが提供され、解析と表示のパフォーマンスが向上します。 また統計によって、分類コードやリターン情報などの LAS 属性の表示が *.las ファイルに存在する値に制限されるので、フィルタリングとシンボルのエクスペリエンスも強化されます。 COMPUTE_STATS—統計情報を計算します。 これがデフォルトです。 NO_COMPUTE_STATS—統計情報を計算しません。	Boolean
out_las_dataset (オプション)	新しく作成された *.las ファイルを参照する出力 LAS データセット。	LAS Dataset
compression (オプション)	出力 *.las ファイルを圧縮形式にするか、標準の LAS 形式にするかを指定します。 SAME_AS_INPUT—入力と同じ圧縮になります。 このオプションは、入力が LAS データセットである場合にのみ使用できます。その場合はこれがデフォルトです。 NO_COMPRESSION—出力は標準の LAS 形式 (**.las) になります。 これは、入力が点群シーン レイヤーである場合のデフォルトです。 ZLAS—出力 *.las ファイルは zLAS 形式で圧縮されます。	String

派生した出力

名前	説明	データ タイプ
out_folder	*.las ファイルが書き込まれるフォルダー。	Folder

コードのサンプル

LAS の抽出の例 1 (Python ウィンドウ)

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

import arcpy
from arcpy import env

env.workspace = 'C:/data'
arcpy.ddd.ExtractLas('test.lasd', 'c:/lidar/subset', boundary='study_area.shp',
                     name_suffix='subset', remove_vlr=True, 
                     rearrange_points='REARRANGE_POINTS', 
                     out_las_dataset='extracted_lidar.lasd')

LAS の抽出の例 2 (スタンドアロン スクリプト)

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

'''****************************************************************************
Name: Split Large LAS File
Description: Divides a large LAS file whose point distribution covers the full
             XY extent of the data into smaller files to optimize performance
             when reading lidar data.
****************************************************************************'''
# Import system modules
import arcpy
import tempfile
import math

in_las_file = arcpy.GetParameterAsText(0)
tile_width = arcpy.GetParameter(1) # double in LAS file's XY linear unit
tile_height = arcpy.GetParameter(2) # double in LAS file's XY linear unit
out_folder = arcpy.GetParameterAsText(3) # folder for LAS files
out_name_suffix = arcpy.GetParameterAsText(4) # basename for output files
out_lasd = arcpy.GetParameterAsText(5) # output LAS dataset


try:
    temp_lasd = arcpy.CreateUniqueName('temp.lasd', tempfile.gettempdir())
    arcpy.management.CreateLasDataset(in_las_file, temp_lasd, 
                                      compute_stats='COMPUTE_STATS')
    desc = arcpy.Describe(temp_lasd)
    total_columns = int(math.ceil(desc.extent.width/tile_width))
    total_rows = int(math.ceil(desc.extent.height/tile_height))
    digits = int(math.log10(max(cols, rows))) + 1
    for row in range(1, total_rows+1):
        yMin = desc.extent.YMin + tile_height*(row-1)
        yMax = desc.extent.YMin + tile_height*(row)
        for col in range (1, total_columns+1):
            xMin = desc.extent.XMin + tile_width*(col-1)
            xMax = desc.extent.XMax + tile_width*(col)
            name_suffix = '_{0}_{1}x{2}'.format(out_name_suffix, 
                                                str(row).zfill(digits), 
                                                str(col).zfill(digits))
            arcpy.ddd.ExtractLas(temp_lasd, out_folder, 
                                 arcpy.Extent(xMin, yMin, xMax, yMax),
                                 name_suffix=name_suffix, 
                                 rearrange_points='REARRANGE_POINTS',
                                 compute_stats='COMPUTE_STATS')
    arcpy.env.workspace = out_folder
    arcpy.management.CreateLasDataset(arcpy.ListFiles('*{0}*.las'.format(out_name_suffix)),
                                      out_lasd, compute_stats='COMPUTE_STATS',
                                      relative_paths='RELATIVE_PATHS')
except arcpy.ExecuteError:
    print(arcpy.GetMessages())

環境

現在のワークスペース, 範囲, 出力座標系, 地理座標系変換

ライセンス情報

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