LAS データセット → ラスター (LAS Dataset to Raster) (変換)

概要

LAS データセットが参照する lidar ポイントに格納された高さ、強度、あるいは RGB の値を使用してラスターを作成します。

LAS データセット → ラスター (LAS Dataset to Raster)

使用法

  • LAS データセット レイヤーにより、レイヤーのフィルター設定内の分類コード、分類フラグ、リターンの組み合わせを選択することで表示および処理される LAS ポイントを制限できます。フィルターは、[レイヤー プロパティ] ダイアログ ボックスまたは [LAS データセット レイヤーの作成 (Make LAS Dataset Layer)] ツールを使用して定義できます。

  • LAS データセット レイヤーは、LAS データセットから参照される場合のあるサーフェス制御フィーチャの適用の制御にも使用できます。この制御は、LAS データセットをトライアングル サーフェスとして表示または処理する場合に適用されます。

  • 大規模なラスターをエクスポートする場合、Z 値の精度要件を整数データで表現できるようであれば、[出力データ タイプ] を整数に指定して、ディスク容量を節約することを検討してください。

  • 分析範囲の境界をクリップ制限として入力 LAS データセットの定義に含めることを強くお勧めします。その理由の 1 つは、調査する実データの範囲外で内挿が起こらないようにするためです。次に、データ範囲が正しく定義されていないときに Natural Neighbor 法を使うと著しくパフォーマンスが低下する可能性があるからです。

    メモ:

    BINNING を使うと、クリップ制限、消去制限ならびに置換制限のみが使用されます。しかしブレークラインとアンカー位置は使用されません。TIN 化のオプションを使うとすべての制限を使用できますが実行に長い時間がかかります。

構文

arcpy.conversion.LasDatasetToRaster(in_las_dataset, out_raster, {value_field}, {interpolation_type}, {data_type}, {sampling_type}, {sampling_value}, {z_factor})
パラメーター説明データ タイプ
in_las_dataset

処理対象の LAS データセット。

LAS Dataset Layer
out_raster

出力ラスターの場所と名前。ジオデータベースまたは Esri Grid などのフォルダーにラスター データセットを格納する場合、ラスター データセットの名前にファイル拡張子は付けません。ファイル拡張子を指定すると、フォルダーに格納する際のラスターの形式を定義できます。たとえば、*.tif を指定すると GeoTIFF が生成され、*.img を指定すると ERDAS IMAGINE 形式のファイルが生成されます。

ラスターを TIFF ファイルとして、またはジオデータベースに格納する場合は、ジオプロセシング環境設定を使用して、そのラスター圧縮のタイプと品質を指定できます。

Raster Dataset
value_field
(オプション)

ラスター出力を生成するために使用する LIDAR データ。

  • ELEVATIONlidar ファイルから得られる高さがラスターの作成に使用されます。これがデフォルトです。
  • INTENSITYlidar ファイルから得られる強度情報がラスターの作成に使用されます。
  • RGBLIDAR ポイントの RGB 値は、3 バンド画像の作成に使用されます。
String
interpolation_type
"BINNING {cell_assignment_type} {void_fill_method}" or "TRIANGULATION {interpolation_method} {point_thinning_type} {point_selection_method} {resolution}"
(オプション)

出力ラスターのセル値の決定に使用される内挿法。

ラスター化方法は、セルの範囲内に含まれるポイントを使用して各出力セル値を決定する [集約方法]、および LAS ポイントを含まないセルの値を決定する [穴埋め] を提供します。

集約方法

  • 平均 - セル内の全ポイントの平均値を割り当てます。これがデフォルトです。
  • 最小 - セル内にあるポイントの最小値を割り当てます。
  • 最大 - セル内にあるポイントの最大値を割り当てます。
  • IDW - 逆距離加重内挿法を使ってセル値を決定します。
  • 最近隣内挿法 - 最近隣内挿法を使ってセル値を決定します。

穴埋め

  • なし - セルに NoData が割り当てられます。
  • シンプル - NoData セルを取り囲むデータ セルの値を平均化して小さな隙間をなくします。
  • リニア - 無効な領域まで TIN 化し、その値にリニア内挿法を用いてセル値を決定します。これがデフォルトです。
  • Natural Neighbor - Natural Neighbor 内挿法を使ってセル値を決定します。

[TIN 化] 内挿法は、TIN に基づいてセル値を生成しますが、ウィンドウ サイズ メソッドを使用して LAS データのサンプリングを間引くことによって処理を高速化することもできます。

TIN 化

  • リニア - リニア内挿法を使用してセル値を決定します。
  • Natural Neighbors - Natural Neighbor 内挿法を使用してセル値を決定します。

ウィンドウ サイズ選択

  • 最大 - 各ウィンドウ サイズ内で最高値を持つポイントが維持されます。これがデフォルトです。
  • 最小 - 各ウィンドウ サイズ内で最小値を持つポイントが維持されます。
  • 平均に最も近い - ウィンドウ サイズ内で全ポイントの平均値に近い値をもつポイントが維持されます。
Interpolate
data_type
(オプション)

出力ラスターに格納される数値のタイプを指定します。

  • FLOAT出力ラスターは 32 ビット浮動小数点を使用します。-3.402823466e+38 から 3.402823466e+38 の範囲の値がサポートされます。これがデフォルトです。
  • INT出力ラスターは適切な整数のビット深度を使用します。このオプションは Z 値を最も近い整数に丸め、各ラスター セルの値に整数を書き込みます。
String
sampling_type
(オプション)

出力ラスターの解像度を定義するために、[サンプリング値] の解釈に使用する方法を指定します。

  • OBSERVATIONSLAS データセット範囲の最も長い辺を分割するセル数を定義します。
  • CELLSIZE出力ラスターのセル サイズを定義します。これがデフォルトです。
String
sampling_value
(オプション)

[サンプリング タイプ] と組み合わせて、出力ラスターの解像度を定義するために使用される値を指定します。

Double
z_factor
(オプション)

Z 値に乗算する係数。これは通常、Z リニア単位から XY リニア単位に変換する場合に使用されます。デフォルトは 1 です。この場合、標高値は変更されません。入力サーフェスの空間参照に距離単位の指定された Z 測地基準系がある場合、このパラメーターは無効になります。

Double

コードのサンプル

LasDatasetToRaster (LAS データセット → ラスター) の例 1 (Python ウィンドウ)

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

arcpy.env.workspace = 'C:/data'
arcpy.LasDatasetToRaster_3d('baltimore.lasd', 'baltimore.tif', 'INTENSITY',
                          'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                          'CELLSIZE', 10, 3.28)
LasDatasetToRaster (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)

    # 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)
    print(arcpy.GetMessages())

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

except Exception as err:
    print(err.args[0])

finally:
    arcpy.management.Delete(lasLyr)

ライセンス情報

  • Basic: 次のものが必要 Spatial Analyst または 3D Analyst
  • Standard: はい
  • Advanced: はい

関連トピック