未分類のLASファイルから地表のDEMを特定しますか？

UAVを使用して航空写真測量から作成されたRGB値を含むLAS形式のデータがあります。点群から地表DEMを抽出するソリューションを見つけようとしています。

SAGA、Fusion、MCC-LIDARを試してみましたが、LASファイルを分類する必要があるようです（当然分類されていません）。プロセスの簡単な説明で誰かが私を正しい方向に向けることができますか？

一般に、一度に約100ミルポイントを処理する必要があります（必要に応じてタイル化できます）。

分類されていない点群からのLiDAR DEMの生成：

• MCC-LIDAR-マルチスケール曲率分類（MCC）アルゴリズム。
  _{（LASバージョン1.1から1.3をサポート）}

MCC-LIDARは、フォレスト環境で離散リターンLIDARデータを処理するためのコマンドラインツールです（Evans＆Hudak、2007）。

ワークフロー：

• a）未分類の点群。
• b）地上リターンは分類されます。
• c）ベアアースDEM（ラスター）。

コードの例をさらに提供するために、架空の状況を作成してみましょう。

MCC-LIDARは次の場所にインストールされます。

C:\MCC

未分類のLiDAR点群（.lasファイル）は次の場所にあります。

C:\lidar\project\unclassified.las  

ベアアースDEMになる出力は次のとおりです。

C:\lidar\project\dem.asc  

以下の例では、MCCアルゴリズムを使用してグラウンドリターンを分類し、1メートルの解像度で地表DEMを作成します。

#MCC syntax: 
#command
#-s (spacing for scale domain)
#-t (curvature threshold)
#input_file (unclassified point cloud) 
#output_file (classified point cloud - ground -> class 2 and not ground -> class 1)
#-c (cell size of ground surface)
#output_DEM (raster surface interpolated from ground points)

C:\MCC\bin\mcc-lidar.exe -s 0.5 -t 0.07 C:\lidar\project\unclassified.las C:\lidar\project\classified.las -c 1 C:\lidar\project\dem.asc

スケール（s）と曲率しきい値（t）のパラメーターがどのように機能するかを理解するには、MCC-LiDARの実行方法と、Evans and Hudak（2007）。

手数料/ラベリングエラーを回避するために、パラメーターを調整する必要があります（ポイントが地面に属していると分類されているが、実際には植生または建物に属している場合）。例えば：

MCC-LIDARは、Thin Plate Spline（TPS）内挿法を使用して、地上点を分類し、地表DEMを生成します。

_参照：

• _{エヴァンス、ジェフリー・S。Hudak、Andrew T.2007 。森林環境で離散リターンLiDARを分類するためのマルチスケール曲率アルゴリズム。地球科学とリモートセンシングに関するIEEEトランザクション。45（4）：1029-1038。}

_{地上点分類アルゴリズムに関するその他のオプションについては、Mengほかを参照してください。（2010）：}

• _{メン、X; Currit、N .; Zhao、K.（2010）。空挺レーザー測量データ用グランドフィルタリングアルゴリズム：重要な課題のAレビュー。リモートセンシング、2（3）、833–860。doi：10.3390 / rs2030833}

私はと思いLasToolsは、あなたのニーズに合うかもしれませんが、参照LASGroundを。ライセンスは、ツールによっては少しおかしいです。ツールは、購入前にダウンロードして評価できます。また、製品は比較的安価です。

FUSION（ここでは手動）のGroundFilterコマンドで幸運がありました。私は4,000万ポイント（未分類）を処理するのに問題がなかったので、1億の問題は予期していません。

これは、Simple Morphological Filter（SMRF ）またはProgressive Morphological Filter（PMF）アルゴリズムを使用するpdal フィルターで実行できます。

クイック

pdal ground --cell_size=5 --extract input.laz out-bare-earth.laz

PMFを使用して5アースユニットセルサイズで地表圧縮LASファイルを作成します。（ドキュメント）

詳細については、ProgressiveMorphologicalFilterセグメンテーションチュートリアルを使用したグラウンドリターンの特定を参照してください。

より複雑な、SMRFの使用

次のパイプラインの例：

• SMRFフィルターを適用し、cellサイズオプションを2.0（座標系単位）と0.75のしきい値に拡大します。
• 新しく分類された地上点のみを選択します（地上2のLAS標準値です）
• 選択を非圧縮LAS出力ファイルに書き込みます（拡張子は.lazに変更して圧縮してください）。

コマンド： pdal pipeline "classify-ground-smrf.json"

JSONパラメータファイル：

{
    "pipeline": [
        "inputfile.laz",
        {
            "type":"filters.smrf",
            "cell": "2.0",
            "threshold": "0.75"
        },
        {
            "type":"filters.range",
            "limits":"Classification[2:2]"
        },
        "out/smurf_classifed.las"
    ]
}

地上のみ抽出

この例では、a）地面/非地面に分類し、b）「地面より上の高さ」属性を追加し、c）地面より上のポイント2.0（座標系単位）のみをエクスポートします。

{
    "pipeline": [
        "input.laz",
        {
            "type": "filters.assign",
            "assignment": "Classification[:]=0"
        },
        {
            "type": "filters.smrf"
        },
        {
            "type": "filters.hag"
        },
        {
            "type": "filters.range",
            "limits": "HeightAboveGround[2:]"
        },
        {
            "filename":"above-ground.laz"
        }
    ]
}

Brad Chambers、https： //lists.osgeo.org/pipermail/pdal/2017-July/001367.htmlから改作