LAS データセット レイヤーの作成 (Make LAS Dataset Layer)

概要

LAS ポイントにフィルターを適用し、サーフェス制御フィーチャの適用をコントロールできる LAS データセット レイヤーを作成します。

使用法

• LAS データセット レイヤーは、LIDAR データをフィルタリングした後で解析を実行するために使用できます。たとえば、地表面の測量値を示すクラス コード 2 を対象としてフィルタリングした後、[LAS データセット → ラスター (LAS Dataset To Raster)] ツールの入力として結果レイヤーを使用することで、デジタル標高モデル (DEM) を簡単に作成できます。

  メモ:

  結果レイヤーは、[レイヤー ファイルの保存 (Save To Layer File)] ツールを使用してレイヤー ファイルとして保存することができます。

構文

arcpy.management.MakeLasDatasetLayer(in_las_dataset, out_layer, {class_code}, {return_values}, {no_flag}, {synthetic}, {keypoint}, {withheld}, {surface_constraints})

パラメーター	説明	データ タイプ
in_las_dataset	処理対象の LAS データセット。	LAS Dataset Layer
out_layer	作成される LAS データセット レイヤーの名前。バックスラッシュまたはスラッシュを使用して、グループ レイヤーを示すことができます。	LAS Dataset Layer
class_code [class_code,...] (オプション)	分類コードで LAS ポイントをフィルタリングできます。有効な値の範囲は、LAS データセットが参照する LAS ファイルのバージョンがサポートしているクラス コードによって変わります。デフォルトでは、すべてのクラス コードが選択されます。 0 —分類方法で処理しません 1 —分類方法で処理しますが、クラス コードは決定できません。 2 —地表計測値 3 —地域で標高が低いと見なされる植生 4 —地域で標高が中程度と見なされる植生 5 —地域で標高が高いと見なされる植生 6 —屋根および壁を持つ構造 7 —地面に近い誤りのデータまたは不要なデータ 8 —以降の使用のための予約、LAS 1.1 - 1.3 のモデル キー ポイントに使用されます。 9 —水域 10 —列車が使用する鉄道の線路 11 —路面 12 —以降の使用のための予約、LAS 1.1 - 1.3 のオーバーラップ ポイントに使用されます。 13 —電線の周囲の標識 14 —電線 15 —頭上の電線を支えるために使用する鉄塔 16 —電気回路を結合する機械部品 17 —橋の表面。 18 —地面から遠い誤りのデータまたは不要なデータ 19 - 63 —ASPRS 指定用の予約済みクラス コード。 64 - 255 —ユーザー定義可能なクラス コード。	String
return_values [return_values,...] (オプション)	LAS ポイントのフィルタリングに使用するリターン値。値が指定されていない場合は、すべてのリターンが使用されます。 Last Return —ラスト リターン First of Many —多数の中の最初 Last of Many —多数の中の最後 Single Return —シングル リターン 1 —1 番目のリターン 2 —2 番目のリターン 3 —3 番目のリターン 4 —4 番目のリターン 5 —5 番目のリターン 6 —6 番目のリターン 7 —7 番目のリターン 8 —8 番目のリターン 9 —9 番目のリターン 10 —10 番目のリターン 11 —11 番目のリターン 12 —12 番目のリターン 13 —13 番目のリターン 14 —14 番目のリターン 15 —15 番目のリターン	String
no_flag (オプション)	分類フラグが割り当てられていないデータ ポイントを表示および解析のために有効化するかどうかを指定します。 INCLUDE_UNFLAGGED —フラグが割り当てられていないポイントを表示します。これがデフォルトです。 EXCLUDE_UNFLAGGED —フラグが割り当てられていないポイントを表示しません。	Boolean
synthetic (オプション)	Synthetic のフラグが設定されたデータ ポイント、つまり LIDAR 以外のデータ ソースからのポイントを表示および解析のために有効化するかどうかを指定します。 INCLUDE_SYNTHETIC —Synthetic ポイントを表示します。これがデフォルトです。 EXCLUDE_SYNTHETIC —Synthetic ポイントを表示しません。	Boolean
keypoint (オプション)	モデル キー ポイントのフラグが設定されたデータ ポイント、つまり間引きの対象外の重要な測定値を表示および解析のために有効化するかどうかを指定します。 INCLUDE_KEYPOINT —モデル キー ポイントを表示します。これがデフォルトです。 EXCLUDE_KEYPOINT —モデル キー ポイントを表示しません。	Boolean
withheld (オプション)	Withheld のフラグが設定されたデータ ポイント (通常は、不要なノイズ測定値を示す) を表示および解析のために有効化するかどうかを指定します。 EXCLUDE_WITHHELD —Withheld ポイントを表示しません。これがデフォルトです。 INCLUDE_WITHHELD —Withheld ポイントを表示します。	Boolean
surface_constraints [surface_constraints,...] (オプション)	レイヤー内で有効化されるサーフェス制御フィーチャの名前。すべてのサーフェス制御がデフォルトで有効となります。	String

コードのサンプル

MakeLasDatasetLayer (LAS データセット レイヤーの作成) の例 1 (Python ウィンドウ)

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

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.MakeLasDatasetLayer_management('Baltimore.lasd', 'Baltimore Layer', 
                                     2, 'LAST', 'INCLUDE_UNFLAGGED', 
                                     'EXCLUDE_SYNTHETIC', 'INCLUDE_KEYPOINT',
                                     'EXCLUDE_WITHHELD')

MakeLasDatasetLayer (LAS データセット レイヤーの作成) の例 2 (スタンドアロン スクリプト)

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

'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
             ground measurements from LAS files to a raster using a 
             LAS dataset. This sample is designed to be used as a script
             tool.
*********************************************************************'''
# Import system modules
import arcpy

try:
    # Set Local Variables
    inLas = arcpy.GetParameterAsText(0)
    recursion = arcpy.GetParameterAsText(1)
    surfCons = arcpy.GetParameterAsText(2)
    classCode = arcpy.GetParameterAsText(3)
    returnValue = arcpy.GetParameterAsText(4)
    spatialRef = arcpy.GetParameterAsText(5)
    lasD = arcpy.GetParameterAsText(6)
    outRaster = arcpy.GetParameterAsText(7)
    cellSize = arcpy.GetParameter(8)
    zFactor = arcpy.GetParameter(9)
    if arcpy.ProductInfo == 'ArcView':
        arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
    # Execute MakeLasDatasetLayer
    lasLyr = arcpy.CreateUniqueName('Baltimore')
    arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
    # Execute LasDatasetToRaster
    arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
                              'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                              'CELLSIZE', cellSize, zFactor)
    arcpy.GetMessages()
    
except arcpy.ExecuteError:
    print(arcpy.GetMessages())
    
except Exception as err:
    print(err.args[0])
    
finally:
    arcpy.management.Delete(lasLyr)

環境

• 現在のワークスペース

ライセンス情報

• Basic: はい
• Standard: はい
• Advanced: はい