LAS の抽出 (Extract LAS) (3D Analyst)—ArcGIS Pro

サマリー

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

図

使用法

このツールは、LAS データセット、点群 .slpk ファイル、I3S 点群シーンレイヤーサービスから .las または .zlas ファイルを作成します。入力が LAS データセットの場合、LAS データセットによって参照される .las ファイルまたは .zlas ファイルごとに 1 つの出力ファイルが作成されます。入力が点群シーンレイヤーの場合、抽出されるポイントの数に応じて出力ファイルの数も異なります。各出力ファイルには約 1680 万ポイントというサイズ制限があります。この閾値を超えるファイルや、これよりポイント数が少ないファイルもあります。
分類コード、分類フラグ、リターン値の任意の組み合わせを使用して入力点群をフィルターすることで、抽出されるポイントを制限できます。ポイントフィルターを定義するには、マップまたはシーンにレイヤーを読み込んでから、[レイヤープロパティ] ダイアログボックスでフィルターを指定します。ポイントフィルターは、[LAS データセットレイヤーの作成 (Make LAS Dataset Layer)] ツールを使用して LAS データセットに適用することもできます。さらに、処理範囲と抽出境界を指定することで、入力点群のサブセットを抽出することもできます。処理範囲を抽出境界とともに指定すると、両者の共通部分を使用して出力ファイルが作成されます。処理範囲と抽出境界がどちらも定義されない場合、入力点群全体が処理されます。
入力が LAS データセットの場合、出力される LAS ファイルバージョンと出力のポイント記録形式はソースファイルと同じになります。入力が I3S 点群シーンレイヤーサービスまたは .slpk ファイルの場合、出力は LAS ファイルバージョン 1.4 で作成されますが、ポイント記録形式は、次のようにソース点群の属性によって異なります。
- 入力に RGB 値や近赤外の値が含まれていない場合は、ポイント形式 6 が使用されます。
- 入力に RGB 値が含まれているが、近赤外の値は含まれていない場合は、ポイント形式 7 が使用されます。
- 入力に RGB 値と近赤外の値が含まれている場合は、ポイント形式 8 が使用されます。
異なる LAS ファイルバージョンとポイント記録形式を使用するには、このツールからの出力ファイルを [LAS の変換 (Convert LAS)] ツールを使用して処理することで生成します。
.las ファイルの空間参照は通常、このファイル内で定義されています。この情報が見つからないか、正しくない場合は、補助的な .prj ファイルを使用して LAS ファイルの空間参照を定義できます。 .prj ファイルは .las ファイルと同じ場所に配置し、.las ファイルと同じ名前にする必要があります。存在する場合、.las ファイルの空間参照情報が .prj ファイルによってオーバーライドされます。空間参照が見つからないか、誤って定義されている .las ファイルで使用する .prj ファイルを作成するには、個別の .las または .zlas ファイルに対して [投影法の定義 (Define Projection)] ツールを使用するか、[LAS データセットの作成 (Create LAS Dataset)] ツールで [LAS ファイルの PRJ を作成] パラメーターを使用します。
出力座標系環境設定を指定することで、入力点群の水平座標系を別の空間参照に再投影することができます。また、鉛直座標系を通じて鉛直単位が定義されるようにして、z 単位変換を実施することもできます。ただし、現在、鉛直測地基準系変換はサポートされていません。水平座標系変換で元の鉛直座標系を維持できるようにする必要があります。
ポイントの空間クラスタリングに対応しないバイナリシーケンスの .las ファイルに LAS ポイントレコードを保存することは珍しいことではありません。このような分布のデータが検索される場合、LAS ポイントを表すバイナリレコードへのアクセス効率が低下する可能性があります。作成される .las ファイルのポイントを再配置することで、視覚化やその他の空間操作のデータが最適化されます。再配置オプションが有効化されると、統計情報が自動的に計算されます。 LAS ポイントを再配置しない場合、統計情報の計算を有効化または無効化することができます。統計情報の計算により、空間クエリが最適化され、クラスコードのサマリーが提供され、.las ファイルに表示される値が返されます。ただし、このツールの処理に時間がかかるようになります。作成される .las ファイルが ArcGIS で使用されない場合、ツールの処理時間を短縮するために、統計情報の計算を無効化することができます。
[抽出境界] パラメーターを使用すると、入力点群が不規則な形状にクリップされます。このパラメーターでは、比較的広い地域を含む 1 つのポリゴンフィーチャなど、隣接する分析範囲を定義するポリゴン境界をサポートしています。隣接する対象地域を表す複数のポリゴンを入力として指定できますが、境界フィーチャでは、樹木や建物など、多数の小規模なフィーチャを表すポリゴンはサポートしていません。入力点群を効率よくクリップできるように、複数のポリゴンがマージされます。このため、多数の小規模なポリゴンを含む境界フィーチャを指定すると、パフォーマンスが低下するので、実行しないでください。

パラメーター

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

派生した出力

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

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 ファイル)、または I3S 点群サービス。 I3S 点群サービスを処理するには、エクスポートプロパティが有効になっている必要があります。	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 (オプション)	入力点群をクリップする際に使用されるポリゴン境界。この境界は、一般的に 1 つの隣接ポリゴンで定義される主要な分析範囲を表す必要がありますが、空間的なギャップのある共通の分析範囲を表す場合には、複数のポリゴンで構成することもできます。	Feature Layer
process_entire_files (オプション)	処理範囲の適用方法を指定します。 PROCESS_EXTENT—対象地域と交差している LAS ポイントだけが処理されます。これがデフォルトです。 PROCESS_ENTIRE_FILES—.las ファイルの一部が対象地域と交差している場合、対象地域の外側にあるポイントも含め、そのファイル内のすべてのポイントが処理されます。	Boolean
name_suffix (オプション)	各出力 .las ファイルの名前に追加されるテキスト。各ファイルは、ソースファイルからベース名を継承し、このパラメーターで指定された接尾辞が付加されます。	String
remove_vlr (オプション)	可変長レコード (VLR) を削除するかどうかを指定します。各 .las ファイルには、それを生成したソフトウェアによって追加された VLR のセットが含まれていることがあります。これらのレコードの意味は、通常は元のソフトウェアしか認知していません。この情報を理解するアプリケーションによって出力 LAS データが処理されるのでなければ、この VLR を保持しても付加価値的な機能は何ももたらされません。 VLR を削除することで、そのサイズ合計やそれを含むファイルの数によっては、かなり大きなディスク容量を節約することができます。 MAINTAIN_VLR—処理対象の .las ファイルの可変長レコードに格納されているすべてのデータを削除せず、抽出されたファイル内に維持します。これがデフォルトです。 REMOVE_VLR—余分な可変長レコードを .las ファイルから削除します。	Boolean
rearrange_points (オプション)	点群分類の読み込みや更新のパフォーマンスを最適化するために、.las ファイルまたは .zlas ファイルのポイントを整理するかどうかを指定します。 MAINTAIN_POINTS—ポイントの順序を整理しません。 REARRANGE_POINTS—ポイントの順序は、データの読み込みを最適化する空間クラスターに整理されます。整理されたポイントにより、それ以降に点群で行われる操作のパフォーマンスが向上します。これがデフォルトです。	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

コードのサンプル

ExtractLas の例 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')

ExtractLas の例 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

サマリー

図

使用法

パラメーター

派生した出力

派生した出力

コードのサンプル

環境

ライセンス情報

関連トピック

このトピックの内容