LAS のカラー化 (Colorize LAS) (3D Analyst)

サマリー

正射投影画像から LAS ポイントに色と近赤外の値を適用します。

図

使用法

• RGB 情報を使用して LAS ポイントを表示すると、没入感のある写実的な表示によって独特の表示体験を提供できます。LAS ポイントに重ね合わせた実際の画像のビジュアル フィードバックにより、ポイント クラウドの不連続特性に関する洞察を得ることができます。これは、データの対話的分類、新しいフィーチャのデジタイズ、3D 距離計測値を取得する際のリファレンス ポイントの確立などに役立ちます。

• LAS ポイントに色を適用する際に使用する画像は、取得されるフィーチャとベスト マッチするよう、LIDAR 測量のときと同じものが最適です。それがない場合、建設などのイベントや落葉による季節の違いによる差異を最小にするために、そのスキャン日付にできるだけ近い画像を使用します。

• ソース画像が多くのタイル化画像で構成されている場合、これらのタイルをモザイク データセットに読み込んで、単一データセットとして参照し、LAS データをカラー化するための画像入力として使用できるようにすることを検討してください。モザイク データセットの作成の詳細

• LAS ポイントの近赤外値の保存に対応しているのは、ポイント記録形式 8 の LAS ファイル バージョン 1.4 のみです。

• 処理対象範囲が定義されている場合は、処理範囲と重なる LAS ファイル全体がカラー化されます。入力 LAS ファイルのサブセットのみが必要な場合、[LAS の抽出 (Extract LAS)] ツールを使用して、サブセットをクリップし、このツールの入力として結果ファイルを使用することを検討します。

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

パラメーター

ダイアログ

ラベル	説明	データ タイプ
入力 LAS データセット	処理される LAS データセット。	LAS Dataset Layer
入力イメージ	LAS ポイントに色を適用するために使用される画像。	Mosaic Layer; Raster Layer
バンド割り当て	入力画像からのバンド。出力 LAS ポイントに関連付けられたカラー チャンネルに割り当てられます。	Value Table
ターゲット フォルダー	出力 .las ファイルが書き込まれる既存のフォルダー。	Folder
出力ファイル名の接尾辞 (オプション)	各出力 .las ファイルの名前に追加されるテキスト。 各ファイルは、ソース ファイルからベース名を継承し、このパラメーターで指定された接尾辞が付加されます。	String
LAS ファイル バージョン (オプション)	作成される出力ファイルの LAS バージョン。 1.2—LAS ファイル バージョン 1.2 が作成されます。 1.3—LAS ファイル バージョン 1.3 が作成されます。 1.4—LAS ファイル バージョン 1.4 が作成されます。これがデフォルトです。	String
ポイント形式 (オプション)	出力 LAS ファイルのポイント レコード形式。 2—ポイント記録形式 2。 3—ポイント記録形式 3 は GPS 時間の保存に対応しています。 7—ポイント記録形式 7。これがデフォルト値で、LAS バージョン 1.4 でのみ使用できます。 8—ポイント記録形式 8 は、近赤外値の保存に対応しています。これは LAS バージョン 1.4 でのみ使用できます。	Long
圧縮 (オプション)	出力 .las ファイルを圧縮形式にするか、標準の LAS 形式にするかを指定します。 NO_COMPRESSION—出力は標準の LAS 形式 (*.las ファイル) になります。 これがデフォルトです。 ZLAS—出力 .las ファイルは zLAS 形式で圧縮されます。	String
ポイントの整理 (オプション)	.las ファイルのポイントを整理するかどうかを指定します。 オフ - .las ファイルのポイントの順序を整理しません。 オン - .las ファイルのポイントを整理します。 これがデフォルトです。	Boolean
統計情報の計算 (オプション)	LAS データセットで参照されている .las ファイルの統計情報を計算するかどうかを指定します。 統計情報を計算することで、.las ファイルごとの空間インデックスが提供され、解析と表示のパフォーマンスが向上します。 また統計によって、分類コードやリターン情報などの LAS 属性の表示が .las ファイルに存在する値に制限されるので、フィルターとシンボルのエクスペリエンスも強化されます。 オン - 統計情報を計算します。 これがデフォルトです。 オフ - 統計情報を計算しません。	Boolean
出力 LAS データセット (オプション)	新しく作成された .las ファイルを参照する出力 LAS データセット。	LAS Dataset

派生した出力

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

Python

arcpy.ddd.ColorizeLas(in_las_dataset, in_image, bands, target_folder, {name_suffix}, {las_version}, {point_format}, {compression}, {rearrange_points}, {compute_stats}, {out_las_dataset})

名前	説明	データ タイプ
in_las_dataset	処理される LAS データセット。	LAS Dataset Layer
in_image	LAS ポイントに色を適用するために使用される画像。	Mosaic Layer; Raster Layer
bands [bands,...]	入力画像からのバンド。出力 LAS ポイントに関連付けられたカラー チャンネルに割り当てられます。	Value Table
target_folder	出力 .las ファイルが書き込まれる既存のフォルダー。	Folder
name_suffix (オプション)	各出力 .las ファイルの名前に追加されるテキスト。 各ファイルは、ソース ファイルからベース名を継承し、このパラメーターで指定された接尾辞が付加されます。	String
las_version (オプション)	作成される出力ファイルの LAS バージョン。 1.2—LAS ファイル バージョン 1.2 が作成されます。 1.3—LAS ファイル バージョン 1.3 が作成されます。 1.4—LAS ファイル バージョン 1.4 が作成されます。これがデフォルトです。	String
point_format (オプション)	出力 LAS ファイルのポイント レコード形式。 2—ポイント記録形式 2。 3—ポイント記録形式 3 は GPS 時間の保存に対応しています。 7—ポイント記録形式 7。これがデフォルト値で、LAS バージョン 1.4 でのみ使用できます。 8—ポイント記録形式 8 は、近赤外値の保存に対応しています。これは LAS バージョン 1.4 でのみ使用できます。	Long
compression (オプション)	出力 .las ファイルを圧縮形式にするか、標準の LAS 形式にするかを指定します。 NO_COMPRESSION—出力は標準の LAS 形式 (*.las ファイル) になります。 これがデフォルトです。 ZLAS—出力 .las ファイルは zLAS 形式で圧縮されます。	String
rearrange_points (オプション)	.las ファイルのポイントを整理するかどうかを指定します。 NO_REARRANGE_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

派生した出力

名前	説明	データ タイプ
output_folder	出力 LAS ファイルが書き込まれるフォルダー。	フォルダー

コードのサンプル

ColonizeLas (LAS のカラー化) の例 1 (Python ウィンドウ)

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

arcpy.env.workspace = 'C:/data'
arcpy.ddd.ColorizeLas('2014_lidar_survey.lasd', '2014_CIR.tif', 
                      'RED Band_1; GREEN Band_2; BLUE Band_3', 
                      'las/rgb', '_rgb', 1.3, 3, 'ZLAS', 
                      'REARRANGE_POINTS')

ColonizeLas (LAS のカラー化) の例 2 (スタンドアロン スクリプト)

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

'''****************************************************************************
       Name: Classify Lidar & Extract Building Footprints
Description: Extract footprint from lidar points classified as buildings, 
             regularize its geometry, and calculate the building height.

****************************************************************************'''
import arcpy

lasd = arcpy.GetParameterAsText(0)
dem = arcpy.GetParameterAsText(1)
footprint = arcpy.GetParameterAsText(2)

try:
    desc = arcpy.Describe(lasd)
    if desc.spatialReference.linearUnitName in ['Foot_US', 'Foot']:
        unit = 'Feet'
    else:
        unit = 'Meters'
    ptSpacing = desc.pointSpacing * 2.25
    sampling = '{0} {1}'.format(ptSpacing, unit)
    # Classify overlap points
    arcpy.ddd.ClassifyLASOverlap(lasd, sampling)
    # Classify ground points
    arcpy.ddd.ClassifyLasGround(lasd)
    # Filter for ground points
    arcpy.management.MakeLasDatasetLayer(lasd, 'ground', class_code=[2])
    # Generate DEM
    arcpy.conversion.LasDatasetToRaster('ground', dem, 'ELEVATION', 
                                        'BINNING NEAREST NATURAL_NEIGHBOR', 
                                        sampling_type='CELLSIZE', 
                                        sampling_value=desc.pointSpacing)
    # Classify noise points
    arcpy.ddd.ClassifyLasNoise(lasd, method='ISOLATION', edit_las='CLASSIFY', 
                               withheld='WITHHELD', ground=dem, 
                               low_z='-2 feet', high_z='300 feet', 
                               max_neighbors=ptSpacing, step_width=ptSpacing, 
                               step_height='10 feet')
    # Classify buildings
    arcpy.ddd.ClassifyLasBuilding(lasd, '7.5 feet', '80 Square Feet')
    #Classify vegetation
    arcpy.ddd.ClassifyLasByHeight(lasd, 'GROUND', [8, 20, 55], 
                                  compute_stats='COMPUTE_STATS')
    # Filter LAS dataset for building points
    lasd_layer = 'building points'
    arcpy.management.MakeLasDatasetLayer(lasd, lasd_layer, class_code=[6])
    # Export raster from lidar using only building points
    temp_raster = 'in_memory/bldg_raster'
    arcpy.management.LasPointStatsAsRaster(lasd_layer, temp_raster,
                                           'PREDOMINANT_CLASS', 'CELLSIZE', 2.5)
    # Convert building raster to polygon
    temp_footprint = 'in_memory/footprint'
    arcpy.conversion.RasterToPolygon(temp_raster, temp_footprint)
    # Regularize building footprints
    arcpy.ddd.RegularizeBuildingFootprint(temp_footprint, footprint, 
                                          method='RIGHT_ANGLES')

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

環境

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