LIDAR データの操作の開始
このチュートリアルでは、LIDAR データから情報を抽出します。 LIDAR (light detection and ranging) は、レーザー光線を使用して、地表を高密度でサンプリングし、極めて精度の高い X、Y、Z ポイント計測値を生成するリモート センシング手法です。 これらのポイント セットは点群と呼ばれます。 LIDAR システムから発射されたレーザー パルスは、植生、建物、橋など、地表に接しているオブジェクトおよび地表の上にあるオブジェクトの両方で反射します。 発射された 1 つのレーザー パルスは、建物の屋根や地表など平坦なサーフェスにぶつかると、1 回反射して戻ってきます (シングル リターン)。 また、地表に向かって進みながら、複数の反射面にぶつかる場合、複数回反射します (複数リターン)。 たとえば、樹木の場合、レーザー パルスは別の高さにある複数の枝葉で反射してから最後に地表で反射します。 これらのリターンはすべて点群にポイントとして表示され、地勢に関する優れた詳細情報を提供します。 LIDAR の詳細については、「LIDAR データとは」をご参照ください。
LIDAR データから 3D 建物の形状を抽出するには、まず、点群を分類し、対象地域の地表と建物を示す点群を識別します。 続いて、地表に分類されたポイントを使用して、地表の標高を連続ラスターで表す DEM を生成します。 次に、建物に分類されたポイントを使用して、2D 建物フットプリントを抽出します。 最後に、これらのエレメントをまとめて 3D 建物フィーチャを生成します。
注意:
このチュートリアルは、LIDAR データとは何か、LAS データセットを作成する方法、点群を分類する方法、点群から数値標高モデル (DEM) を生成する方法、2D 建物フットプリントと 3D マルチパッチ建物を抽出する方法など、LIDAR 関連の概念を伝えることに主眼を置いています。
一方で、Esri では、大規模な 3D ベースマップの生成を目指している都市とその他の政府機関を対象とした 3D Buildings ソリューションを用意しています。 このソリューションは、一連のワークフローで構成されており、3D ビルディング レイヤーの作成と管理を効率化することができます。 その用途については、ArcGIS ラボ「Working with the ArcGIS Solution for 3D Buildings」で紹介しています。
プロジェクトをダウンロードして開く
まず、チュートリアルに必要なデータをすべて含むプロジェクトをダウンロードし、ArcGIS Pro で開きます。
- Building_Extraction .zip ファイルをダウンロードします。
.zip ファイルがコンピューターにダウンロードされます。
- コンピューター上で、ダウンロードした [Building_Extraction.zip] ファイルを参照します。
注意:
お使いの Web ブラウザーによっては、ダウンロードを開始する前に、ファイルの場所を選択するよう求めるメッセージが表示される場合があります。 ほとんどのブラウザーでは、デフォルトでコンピューターのダウンロード フォルダーがダウンロード先の場所になります。
- [Building_Extraction.zip] ファイルを右クリックして、ドキュメント フォルダーなどの見つけやすい場所にファイルを展開します。
次に、ArcGIS Pro でプロジェクトを開きます。
- ArcGIS Pro を起動します。 プロンプトが表示されたら、ArcGIS アカウントの認証情報を入力します。
注意:
ArcGIS Pro へのアクセス権限または組織アカウントがない場合は、ソフトウェア アクセスのオプションをご参照ください。
- ArcGIS Pro で [別のプロジェクトを開く] をクリックします。
- [プロジェクトを開く] ウィンドウで、ダウンロードした [Building_Extraction] フォルダーを参照します。 [Building_Extraction.aprx] をクリックして選択し、[OK] をクリックします。
プロジェクトが開きます。
プロジェクトには、デンマークのコペンハーゲンにあるツボルグ港付近を中心とする 3D 対応ローカル シーンが含まれています。 このあたりは、ツボルグ醸造所のかつての工業用地に近く、多目的地区として再開発されました。 シーンには、この地区の航空写真として [Tuborg_Havn_Ortho_Photo.tif] が含まれており、チュートリアルでは参照用に使用します。
注意:
このチュートリアルで使用する [Tuborg_Havn_Ortho_Photo.tif] 画像およびすべてのデータは、デンマーク政府の Web サイト「Styrelsen for Dataforsyning og Effektivisering」にあります。
ArcGIS Pro では、ローカル シーンは、都市や都市近郊など、限られた空間範囲を持つ 3D データセットの表示に便利です。 ローカル シーンに表示されるデータは、投影座標系を使用する必要があります。 ローカル シーンの詳細と、グローバル シーンとの違いについては「シーン」ドキュメント ページをご参照ください。
- マウス ホイール ボタンを使用して拡大表示し、マップをドラッグして画面移動します。 ツボルグ港付近を調べます。
この地区には、さまざまな近代的な建物、ボートのあるマリーナ、西端および北端の小さな建物、南東側の開発中のエリアがあります。
注意:
この範囲は、チュートリアルをすみやかに進められるよう制限されています。 実際の運用環境では、より広い近郊や都市全体を処理できます。
次に、LIDAR データを準備して表示します。
LAS データセットの作成
LIDAR 点群データは通常、LAS ファイル形式で格納されています。 このチュートリアルでは、ツボルグ港地域は 2 つの LAS ファイルにまたがっています。 それらを 1 つの LAS データセットにまとめます。
注意:
対象地域に広がる複数の LAS ファイルを使用するのは、一般的なことです。 LAS データセットには、分析範囲に必要な数の LAS ファイルを追加できます。
まず、2 つの LAS 点群ファイルを参照します。
- ウィンドウ上部のリボンで [表示] タブをクリックします。 [ウィンドウ] グループで、[カタログ ウィンドウ] をクリックします。
[カタログ] ウィンドウが表示されます。
- [カタログ] ウィンドウで、[フォルダー]、[Building_Extraction]、[Tuborg_Havn_data]、および [LAS_data] を展開します。
[PUNKTSKY_1km_6181_724.las] および [PUNKTSKY_1km_6181_725.las] は、ツボルグ港近辺に広がる LAS ファイルです。
注意:
この 2 つの LAS ファイルは、デンマーク政府が管理する、デンマーク全体を含む LIDAR を生成したプロジェクトから取得したものです。
次に、LAS データセットを作成します。
- [LAS_data] フォルダーを右クリックし、[新規作成] をポイントして [LAS データセット] を選択します。
新しい LAS データセットが [カタログ] ウィンドウの [LAS_Data] フォルダーに追加され、データセットの名前が編集可能になります。
- 「Tuborg_Havn」と入力して Enter キーを押します。
拡張子が自動的に追加され、[Tuborg_Havn.lasd] という名前がリストに表示されます。
[Tuborg_Havn.lasd] は LAS データセットですが、この時点では空です。 次に、データを追加します。
- [Tuborg_Havn.lasd] を右クリックして [プロパティ] を選択します。
- [LAS データセット プロパティ] ウィンドウの [LAS ファイル] タブをクリックします。
- [ファイル] の下の [ファイルの追加] をクリックします。
- [開く] ウィンドウで、[Tuborg_Havn_data] フォルダーを参照し、[LAS_data] フォルダーをダブルクリックして開きます。
- Shift キーを押して [PUNKTSKY_1km_6181_724.las] および [PUNKTSKY_1km_6181_725.las] をクリックし、これらのファイルを両方とも選択します。 [開く] をクリックします。
2 つのファイルがリストに追加されました。 [ポイント数] の値に応じて、各 LAS ファイルには 300 万から 500 万の間のポイントが含まれています。 [ポイント間隔] の値は、ポイントとポイントの間が約 0.3 メートルであることを示しています。
- [一般] タブをクリックします。
このタブには、新しい LAS データセットの全体的な情報がまとめられています。 2 つの LAS データセットがあり、8,127,305 個の LAS ポイントが含まれていることがわかります。 データの [範囲] の値や、水平 (XY) および垂直 (Z) の単位 (メートル) も表示されています。
- [統計情報] タブをクリックします。
このタブでは、[分類コード] の下に、LAS ポイントに割り当てられた分類コードのリストが表示されます。 この時点では、ポイントは分類コードを割り当てられていないので、表示されている分類コードは [未分類] のみです。 このチュートリアルで、地表、建物、ノイズなどの新しいクラスを追加します。
- [座標系] タブをクリックします。
投影座標系は ETRS 1989 UTM Zone 32N です。 これは、元の LAS ファイルから [Tuborg_Havn.lasd] データセットに引き継がれた座標系です。 チュートリアルでは、この座標系を使用します。
- [OK] をクリックして [LAS データセット プロパティ] ウィンドウを閉じます。
LAS データセットの表示と調査
次に、LAS データセットをシーンに追加して調べます。 今のところオルソ フォトは必要ないのでオフにします。
- [コンテンツ] ウィンドウで、[Tuborg_Havn_Ortho_Photo.tif] のチェックボックスをオフにして、レイヤーを非表示にします。
- [カタログ] ウィンドウで、[Tuborg_Havn.lasd] データセットを右クリックして [現在のマップに追加] を選択します。
マップに [Tuborg_Havn.lasd] データセットが表示されます。
- [コンテンツ] ウィンドウで、[Tuborg_Havn.lasd] を展開します。
凡例では、LAS データセットは、ポイントの高さに応じてシンボル表示されています。 最も低いポイントは濃い紫、最も高いポイントは明るい赤で表示されます。 今度は点群を 3D で調べてみましょう。
- シーンの [ナビゲーター] ホイールの上にある [フル コントロールの表示] をクリックします。
[ナビゲーター] ホイールが拡張され、3D ナビゲーション機能が追加されます。
- 中央の [ナビゲーター] ホイールを使用して、シーンの傾斜や回転を行います。
ヒント:
C キーを押しながらシーンをドラッグして画面移動するか、V キーを押しながらシーンをドラッグして傾けることもできます。
- 個々の LAS ポイントを確認できるまで拡大します。
- LAS データセットがワイヤーフレームとして表示されるまで縮小します。
ワイヤーフレームは、両方のファイルの範囲と、そのファイルに含まれているポイントの最大の高さを示します
- 引き続き、点群データセットを調べます。 拡大/縮小、傾斜、回転、および画面移動して、建物、樹木、地面、その他の要素を横から見てみます。
- [クイック アクセス ツールバー] で、[保存] ボタンをクリックして、プロジェクトを保存します。
ヒント:
以前のバージョンの ArcGIS Pro でこのプロジェクトが作成されていることを知らせる警告が表示された場合は、[はい] をクリックして続行します。
Ctrl + S キーを押してプロジェクトを保存することもできます。
2 つの LAS ファイルを含む LAS データセットを作成しました。 また、そのデータセットをマップに追加して調べました。 続いて、LAS データセットを、地表や建物など、いくつかのクラスに分類します。
LAS データセットの分類
ここで、クラウド ポイントの位置や特性を解析し、所属するクラスを決定するツールを使用して、LAS データセットを地表や建物などの複数のクラスに分類します。
注意:
一部の LAS 点群は、すでに分類された状態でプロバイダーから取得されます。 このような場合は、このワークフローの分類ステップをすべてスキップしてください。
LAS 地表ポイントの分類
まず、LAS ポイント クラスが見やすくなるように、シンボルを変更します。
- [コンテンツ] ウィンドウで、[Tuborg_Havn.lasd] が選択されていることを確認します。 [Tuborg_Havn.lasd] を右クリックして、[レイヤーにズーム] を選択します。
- リボンの [LAS データセット レイヤー] タブの [描画] グループで、[シンボル] 下矢印をクリックし、[クラス] を選択します。
シンボルに各ポイントの分類が表示されるようになりました。 現在、すべて未割り当てのため、すべてグレーでシンボル表示されています。
[コンテンツ] ウィンドウには、設定可能なクラスのリストが表示されますが、そのほとんどが LAS データセットでは現在使用されていません。
次に、地表ポイントを分類します。
注意:
データを作成した航空機の元のフライト ラインが原因で、LAS ファイル同士がオーバーラップしている場所があります。 重複ポイントはノイズを生成する可能性があるため、まず [LAS オーバーラップの分類] ツールを使用して重複ポイントを特定した後、オフにします。 ただし、現在のデータにオーバーラップは存在しません。 このツールを実行する必要はありません。
- リボンの [表示] タブの [ウィンドウ] グループで、[ジオプロセシング] をクリックします。
[ジオプロセシング] ウィンドウが表示されます。 まず、[LAS の地表分類] ツールを使用します。
- [ジオプロセシング] ウィンドウで「LAS の地表分類」と入力し、結果リストの [LAS の地表分類] をクリックして開きます。
- [LAS の地表分類] ツール ウィンドの [入力 LAS データセット] で、[Tuborg_Havn.lasd] を選択します。
- その他のパラメーターはデフォルトのままにし、[実行] をクリックします。
注意:
LAS 点群内のすべての地表ポイントを特定するため、[LAS の地表分類] ツールでは、シーン全体で一貫して標高が最も低いポイント セットを検索するなどの手法を使用します。
プロセスが完了すると、LAS 地表ポイントはシーン内に茶色で表示されます。
- LAS 地表ポイントを拡大して、ポイントが地表エリアには表示されるものの、建物、植生、車両、またはボートの位置には表示されず、未分類のままである状況を確認します。
水域に位置する一部のポイントは地表に分類されていました。 ただし、今回のワークフローでは、地面と水を区別する必要はないため、問題ありません。
- LAS 点群全体が表示されるまで縮小します。
注意:
LAS ポイントを分類すると、元の LAS データセットが変更されます。 誤りがあって元に戻す必要がある場合、LAS クラス コードの変更ツールを使用できます。 たとえば、コード 2 (地表) のポイントをコード 1 (未分類) に変更できます。 別の選択肢として、元のデータセットのコピーを事前に作成し、バックアップとして使用することもできます。
DEM の生成
LAS 地表ポイントを使用して DEM を作成します。これはワークフローの後半で役立ちます。 DEM は、地表の標高を連続ラスターで表します。 [LAS データセット → ラスター] ツールを使用して DEM を生成します。 このツールは地表ポイントのみで実行する必要があるため、未割り当てのポイントはオフにします。
- リボンの [LAS データセット レイヤー] タブの [フィルター] グループで、[LAS ポイント] をクリックします。
[レイヤー プロパティ] ウィンドウの [LAS フィルター] タブに、さまざまなタイプのポイントが表示されます。 タイプはそれぞれオン/オフを切り替えることができます。 [分類コード] の下に、現在設定可能な 2 つのクラス ([1 未分類] および [2 地表]) が表示されます。
- [レイヤー プロパティ] ウィンドウの [LAS フィルター] タブの [分類コード] の下にある [未分類] クラスをオフにして、そのカテゴリに所属するポイントをオフにし、[OK] をクリックします。
未割り当てのポイント (グレー) がマップに表示されなくなります。 ここで、[LAS データセット → ラスター] ツールを実行して、DEM を生成します。
注意:
LIDAR LAS データセットを入力として使用するツールの実行中は、マップに表示されるポイントのみがプロセスで使用されます。 そのため、この場合は地表ポイントのみが使用され、他のポイントは無視されます。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。
- [LAS データセット → ラスター] ツールを検索して開きます。
- [LAS データセット → ラスター] ツール ウィンドウで、次のパラメーターを設定します。
- [入力 LAS データセット] で [Tuborg_Havn.lasd] レイヤーを選択します。
- [出力ラスター] に「Tuborg_Havn_DEM」と入力します。
- [内挿タイプ] で [TIN 化] を選択します。
- [内挿方法] で [Natural Neighbor] を選択します。
- [サンプリング値] に「0.5」と入力します。
[内挿タイプ] パラメーターの横に表示される警告は無視して構いません。
注意:
LAS 地表ポイントは地表の標高に関する多くの情報を提供していますが、各ポイント間で標高が不明な空間が存在します。 [LAS データセット → ラスター] ツールは数学モデルを使用し、既存の地表ポイントの値に基づいて、これらの不明な地域の標高を計算 (予測) します。 このプロセスを内挿と呼びます。 さまざまな内挿オプションについては、[LAS データセット → ラスター (LAS Dataset to Raster)] ツール ドキュメントをご参照ください。
ラスターはピクセルまたはセルのグリッドです。 [サンプリング値] パラメーターは、各セルのサイズを示しています。 ポイント間の間隔 (約 0.3 メートル) よりも大きい必要があります。 試行錯誤の結果、このデータセットで最良の結果が得られることが判明した 0.5 メートルを使用します (1 メートルも合理的な可能性があります)。 このチュートリアルでは、すべてのラスターがこのセル サイズを使用します。
- その他のパラメーターはすべてデフォルトのままにし、[実行] をクリックします。
プロセスが完了すると、出力 DEM ラスターが表示されます。 このラスターは、暗い色のポイントが最も低い標高、明るい色のポイントが最も高い標高を表す滑らかなサーフェスです。
- [コンテンツ] ウィンドウで、[Tuborg_Havn.lasd] レイヤーをオフにすると、DEM が見やすくなります。
ツボルグ港は複雑な都市景観であるため、DEM は通常の道路レベルを超えたさまざまなバリエーションを示しています。 これには、地下道や駐車場への入り口、建物群がそびえる高台、マリーナの水位などが含まれます。
- プロジェクトを保存します。
LAS ノイズ ポイントの分類
これで LAS 地表ポイントを特定したので、次のステップでは残りの未分類ポイントからノイズ ポイントを特定します。 ノイズ ポイントとは、異常に高いポイントまたは低いポイントのことを指し、LIDAR データのランダム エラーの結果と考えられます。
- [コンテンツ] ウィンドウで、[Tuborg_Havn_DEM] ラスターをオフにし、[Tuborg_Havn.lasd] レイヤーをオンにします。 必要に応じて、[Tuborg_Havn.lasd] レイヤーをクリックして選択します。
- リボンの [LAS データセット レイヤー] タブで、[LAS ポイント] をクリックします。 表示された [レイヤー プロパティ] ウィンドウの [分類コード] の下で、[地表] クラスをオフにし、[未分類] クラスをオンにします。 [OK] をクリックします。
未分類のポイントがシーン上にグレーで表示されます。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [LAS ノイズの分類] ツールを検索して開きます。
- [LAS ノイズの分類] ツール ウィンドウで、次のパラメーターを設定します。
- [入力 LAS データセット] で [Tuborg_Havn.lasd] レイヤーを選択します。
- [方法] で、[地表からの相対高度] を選択します。
- [高さ検出] セクションを展開して、次のパラメーターに値を設定します。
- [地表ラスター サーフェス] で、[Tuborg_Havn_DEM] を選択します。
- [最小高さ] に、「-2」と入力します。
- [最大高さ] は空白のままにします。
- [実行] をクリックします。
高ノイズを特定するため、今回は [絶対高度] メソッドを使用して、ツールを再度実行します。 標高が 42 メートル以上のポイントを高ノイズに分類します。
注意:
最大値に 42 メートルを選択しているのは、ツボルグ近郊の建物の最大標高は約 40 メートルであるためです。 これは LAS 点群を探索することで確認できます。 建物の最大標高を見つけるには、拡大して確認できる最も高いポイントをクリックします。 ポップアップ ウィンドウには、そのポイントの標高がメートル単位で表示されます。
- [LAS ノイズの分類] ツール ウィンドウの [方法] で、[絶対高度] を選択します。
- [高さ検出] セクションの [最小高さ] で、既存のテキストを削除します。 [最大高さ] に、「42」と入力します。
- [実行] をクリックします。
一部の LAS ポイントが低ノイズと高ノイズに分類されるようになりましたが、デフォルトでオフになっているため表示されません。
- リボンの [LAS データセット レイヤー] タブで、[LAS ポイント] をクリックします。 表示された [レイヤー プロパティ] ウィンドウの [分類コード] で、[7 ノイズ] および [18 高ノイズ] クラスをオンにし、その他のクラスをすべてオフにして、[OK] をクリックします。
- ツボルグ港近郊を探索し、赤色でシンボル表示された LAS ノイズ ポイントを見つけます。
たとえば、近郊の南側に高ノイズ ポイントが確認できます。 建設用クレーンに対応していると考えられます。
残りの解析を行うため、これらをオフにします。
- [LAS データセット レイヤー] タブで、[LAS ポイント] をクリックします。 表示された [レイヤー プロパティ] ウィンドウの [分類コード] で、[未分類] クラスをオンにし、その他のクラスをすべてオフにします。 [OK] をクリックします。
注意:
ここで、多くの GIS ワークフローで使用される、建物や樹木などのフィーチャを含む地表の標高を表示する数値表層モデル (DSM) を、別のラスター レイヤーで生成することもできます。 DEM を生成したときと同じように、[LAS データセット → ラスター] ツールを使用しますが、異なるポイントのサブセットを選択して処理します。 [レイヤー プロパティ] ウィンドウの、[LAS フィルター] タブの [分類コード] の下で、[ノイズ] と [高ノイズ] 以外のすべてのチェックボックスをオンにし、[リターン値] の下で [1] のチェックボックスをオンにします。 これは、最初のリターン ポイントのみを考慮することを意味します。最初のリターン ポイントは、樹木の最上部や車の屋根、裸地など、地形内のすべてのフィーチャのなかで最も高いポイントを表すためです。また、ノイズ ポイントは含めません。
[LAS データセット → ラスター] ツールでは、DEM のときと同じパラメーターを選択しますが、例外として、[内挿タイプ] で [ラスター化] オプション、[集約] で [最大]、[穴埋め] で [Natural Neighbor] を選択します。 DSM を使用した分析例は、「太陽光発電の能力の推定」をご参照ください。
建物ポイントの分類
これで LAS 地表ポイントとノイズ ポイントが分類され、表示されないようにしたため、次は建物ポイントを分類します。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [LAS の建物分類] ツールを検索して開きます。
- [LAS の建物分類] ツール ウィンドウで、次のパラメーターの値を設定します。
- [入力 LAS データセット] で [Tuborg_Havn.lasd] レイヤーを選択します。
- [屋上の高さの最小値] が「2」に設定され、単位がメートルに設定されていることを確認します。
- [最小エリア] に「10」と入力し、単位が平方メートルに設定されていることを確認します。
- [分類方法] で [積極的] を選択します。
注意:
[LAS の建物分類] ツールでは、建物ポイントを特定するために、複数の方法を組み合わせて使用しています。 ツールを実行する前に、地表とノイズ ポイントを分離しておくことが重要です。 このツールは、未分類ポイント (建物とノイズ ポイント以外のポイント) のうち、レーザー パルスごとのリターンが 1 回のみだった平坦なサーフェスを識別します。 一方、複数のリターンがあったエリアは、樹木や他の草木を表している可能性が高くなります。光が複数の高さレベルの葉に届いて反射するためです。
[屋上の高さの最小値] および [最小エリア] パラメーターは、標高が低すぎるか面積が小さすぎるサーフェスが、誤って建物に分類されないようにするために重要です。
[分類方法] では、点を控えめに建物として分類するか、積極的に分類するかを指定します。
このチュートリアルでは、建物に正しく分類されるポイントの数を最大限に高めながら、誤検知を最小限に抑えるパラメーター値を試行錯誤しながら選択しました。 引き続きご自身で試してみてください。 データの小さな分析範囲のコピーを作成して最適なパラメーターを特定することをお勧めします。
- [実行] をクリックします。
- 処理が完了したら、リボンの [LAS データセット レイヤー] タブで、[LAS ポイント] をクリックします。 [レイヤー プロパティ] ウィンドウの [分類コード] で、[6 建物] クラスをオンにし、その他のクラスをすべてオフにして、[OK] をクリックします。
建物ポイントがシーンに表示されます。
注意:
未分類のポイントがまだたくさんあります。 たとえば、[高さによる LAS の分類] ツールを使用して、[低植生]、[中植生]、および [高植生] クラスに所属するポイントを特定するなどして、引き続きポイントを詳細に分類できます。 ただし、3D 建物の抽出に植生情報は必要ないため、このチュートリアルではこれ以上の分類は実行しません。
- プロジェクトを保存します。
LIDAR 点群を、地表、低ノイズ、高ノイズ、および建物に分類しました。 地表の標高を連続ラスターで表す DEM も生成しました。 次に、LAS 建物ポイントを使用して、対象地域全体の 2D 建物フットプリントを抽出します。
建物フットプリントの抽出
次に、建物に分類された LAS ポイントを使用して、2D 建物フットプリントを抽出します。 まず、LAS 建物ポイントの場所を示すラスターを生成します。 続いて、ラスターによって建物の存在を示すポリゴンを生成します。 最後に、ポリゴンをクリーンナップして、最終的な建物フットプリント ベクター レイヤーを生成します。
建物ラスターの生成
ワークフローのこの部分は、表示により多くの処理リソースを使用する 3D 表示機能が不要なので、2D マップで行います。 シーンを 2D マップに変換します。
- 必要に応じて、リボンの [表示] タブにある [ウィンドウ] グループで、[カタログ ウィンドウ] をクリックします。
- [カタログ] ウィンドウで、[マップ] を展開します。 [Scene] を右クリックして [変換] をポイントし、[マップに変換] を選択します。
元の 3D シーンのレイヤーを含む新しい 2D マップが、[Scene_2D] という名前で表示されます。 名前を変更します。
- [カタログ] ウィンドウで、[Scene_2D] を右クリックして [名前の変更] をクリックします。 「2D Map」と入力して Enter キーを押します。
[2D Map] がリストに表示されます。
- [2D Map] タブが選択されていることを確認し、必要に応じて、LAS 建物ポイントを表示できるまで拡大します。
この時点では、マップには LAS 建物ポイントが表示されています。
作成する建物フットプリントと実際の建物の関係がわかるよう、オルソ フォトを有効にします。
- [コンテンツ] ウィンドウで、[Tuborg_Havn_Ortho_Photo.tif] イメージ レイヤーを有効にします。
- 必要に応じて、[コンテンツ] ウィンドウで [Tuborg_Havn.lasd] レイヤーを [Tuborg_Havn_Ortho_Photo.tif] の上の該当位置にドラッグします。
オルソ フォトは、実際の建物がどのように見えるかを示します。 ただし、写真の遠近法により、建物のフットプリントと下にある建物を視覚的に完全に一致させることが難しい場合があります。
次に、LAS 建物ポイントの位置に対応するラスターを作成します。 ラスターの作成には、[LAS ポイント統計をラスターに出力] ツールを使用します。このツールは、ラスターの各セルの LAS 建物ポイントを検出します。 建物ポイントがある場合、セルの値は建物のクラス コード (6) になります。建物ポイントがない場合、セルは空 (NoData) のままです。
- 必要に応じて、リボンの [表示] タブにある [ウィンドウ] グループで、[ジオプロセシング] をクリックします。 必要に応じて、[ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。
- [LAS ポイント統計をラスターに出力] ツールを検索して開きます。
- [LAS ポイント統計をラスターに出力] ツール ウィンドウで、以下のパラメーターを設定します。
- [入力 LAS データセット] で [Tuborg_Havn.lasd] レイヤーを選択します。
- [出力ラスター] で、デフォルト名を使用します。
- [方法] で [クラス コードの最頻値] を選択します。
- [サンプリング値] に「0.5」と入力します。
- [実行] をクリックします。
出力レイヤーが表示されます。
レイヤーを確認しやすくするため、LAS データセットをオフにしてラスターのシンボルを変更します。
- [コンテンツ] ウィンドウで、[Tuborg_Havn.lasd] レイヤーをオフにします。
- [tuborg_laspo] を [Buildings_Raster] に名前変更します。
- [コンテンツ] ウィンドウで、[Buildings_Raster] のシンボルを右クリックしてドロップダウン リストを展開し、[名前の表示] をオンにして [赤紫 (連続)] 配色を選択します。
ラスターがピンクになります。
- 拡大して、ラスターの建物をいくつか調べます。
建物の全体的な形状は正確ですが、建物にいくつかの穴があり、枠線が不規則に見えます。 また、ピンクの形状の一部が実際の建物に対して小さすぎるので、削除する必要があります。 これらの問題は、ワークフローの後半で解決します。ただし、最初の段階で、建物ラスターをポリゴン レイヤーに変換しておきます。
建物フィーチャクラスの生成
続いて、[ラスター → ポリゴン] ツールを使用して、建物ラスターをポリゴン レイヤーに変換します。 ラスター内の建物セルのクラスターが、新しいレイヤーでポリゴンになります。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [ラスター → ポリゴン] ツールを検索して開きます。
- [ラスター → ポリゴン] ツール ウィンドウで、次のパラメーターを設定します。
- [入力ラスター] で [Buildings_Raster] を選択します。
- [出力ポリゴン フィーチャ] で、デフォルト名を使用します。
- [ポリゴンの単純化] をオフにします。
注意:
[ポリゴン単純化] オプションを選択すると、ポリゴンの基本的なクリーンナップを行うことができます。 または、ワークフローの後半でより高度なツールを使用して、より質の高い結果を得ます。
- [実行] をクリックします。
新しいフィーチャ レイヤーが表示されます。 今は、[Buildings_Raster] レイヤーをオフにすることができます。
- [コンテンツ] ウィンドウで、[Buildings_Raster] レイヤーをオフにし、新しい建物フィーチャ レイヤーの名前を [Buildings_Raw] に変更します。
- [Buildings_Raw] のシンボルを右クリックして [フクシア ピンク] を選択します。
ヒント:
色の名前を表示するには色にカーソルを合わせます。
[Buildings_Raw] レイヤー シンボルが更新されます。
建物ポリゴンのクリーンナップ
建物フットプリントを表すポリゴンが作成されました。 ただし、クリーンナップする必要があります。
- 新しいポリゴン レイヤーを拡大して調べます。
建物ラスターと同じ問題があることがわかります。
- 一部のポリゴンは、実際の建物に対して小さすぎるので、削除する必要があります。
- ほとんどのポリゴンに穴があります。
- 建物の縁が不規則です。
クリーンナップ処理の最初の手順として、最も小さいポリゴンを除外します。 まず、削除するポリゴンの面積のカットオフ値を決定します。
- [コンテンツ] ウィンドウで、[Buildings_Raw] レイヤーをクリックして選択します。
- マップで、小さいピンクのポリゴンをいくつかクリックして、情報が表示されたポップアップ ウィンドウを開き、[Shape_Area] 属性の値に注目してサイズを把握します。 除外しない建物についても同じ作業を行います。
確認したポリゴンのサイズに基づき、除外する面積の値を 70 平方メートルにします。 次に、その値より大きいポリゴンを選択します。
- [ポップアップ] ウィンドウを閉じます。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [属性検索] ツールを検索して開きます。
- [属性検索] ツール ウィンドウの [入力行] で、[Buildings_Raw] を選択します。
- [式] の下で、「Where 句 Shape_Area が 70 以上」という式を作成します。
- [実行] をクリックします。
70 平方メートル以上のすべてのポリゴンが選択されました。
次に、[ポリゴン パートの削除] ツールを使用して、建物ポリゴンの穴を削除します。 調べた結果、50 平方メートル以下の穴を埋める必要があると結論付けることができます。
- 条件以上の大きさの建物ポリゴンを選択した状態で、[ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [ポリゴン パートの削除] ツールを検索して開きます。
- [ポリゴン パートの削除] ツール ウィンドウで、次のパラメーターを設定します。
- [入力フィーチャ] で [Buildings_Raw] を選択します。
- [出力フィーチャクラス] で、デフォルト名を使用します。
- [面積] に「50」と入力し、単位が平方メートルに設定されていることを確認します。
- [実行] をクリックします。
新しいフィーチャ レイヤーが表示されます。 新しいポリゴンのほとんどには穴がありません。また、条件以上の大きさのポリゴンにのみツールが適用されたので、小さいポリゴンはすべて除外されました。
新しいレイヤーが見やすくなるよう、[Buildings_Raw] レイヤーをオフにします。
- [コンテンツ] ウィンドウで、[Buildings_Raw] レイヤーをオフにし、新しいレイヤーの名前を [Buildings_Clean] に変更します。
続いて、建物の枠線を滑らかにします。
- [ジオプロセシング] ウィンドウの [戻る] ボタンをクリックします。 [建物フットプリントの正規化] ツールを検索して開きます。
このツールは、抽出された建物ポリゴンの以下のようなギザギザの枠線をクリーンナップし、滑らかにします。
- [建物フットプリントの正規化] ツール ウィンドウで、次のパラメーターを設定します。
- [入力フィーチャ] で [Buildings_Clean] を選択します。
- [出力フィーチャクラス] で、デフォルト名を使用します。
- [方法] が [直角] に設定されていることを確認します。
- [許容値] に「1」と入力します。
- [密度] に「1」と入力します。
- [精度] に「0.15」と入力します。
[直角] の方法は、完全な直角の建物の作成に重点を置いています。 [任意の角度] などのその他のオプションは、構造が複雑で規模の大きい建物の場合に有用ですが、小さくてシンプルな建物を完全に単純化できない可能性があります。 この演習では行いませんが、ツールを 2 回実行するというオプションがあります。それぞれ、異なる設定と異なる建物を選択して実行します。
[許容値] は、正規化されたフットプリントが、元のフィーチャの境界から外れることができる最大距離です。 通常は、元のラスター セル (この演習では 0.5 メートル) の約 2 倍のサイズを選択することをお勧めします。
[密度] は、正規化されたフィーチャが直線か曲線かを評価するのに使用されるサンプリング間隔を決定します。
[精度] は、正規化プロセスで使用されるグリッドの正確さを定義します。 0.05 ～ 0.25 の値を指定でき、試行錯誤により値を微調整することができます。
- [実行] をクリックします。
正規化された建物ポリゴンを示した新しいフィーチャ レイヤーが表示されます。
- [コンテンツ] ウィンドウで、[Buildings_Clean] レイヤーをオフにし、新しいレイヤーの名前を [Building_Footprints] に変更します。
- プロジェクトを保存します。
手作業による建物フットプリントの改善
[Building_Footprints] レイヤーは、全体的には正確なフットプリントを構築します。 ささいな欠陥はありますが、このワークフローではそれらを気にする必要はありません。 ただし、いくつかのポリゴンに、修正すべきやや大きな問題が含まれている可能性があります。 そのような問題には、手作業ですばやく解決する方法があります。
まず、いくつかのポリゴンは、建物ではなく高い場所に設置された運動場です。 レイヤーからそれらのポリゴンを削除します。
- 以下の画像でハイライト表示されている建物群を探します。
この地域の 2 つのポリゴンは、高い場所に設置された運動場です。
- リボンの [編集] タブの [選択] グループで、[選択] をクリックします。
- 運動場ポリゴンをクリックして選択します。
- リボンの [編集] タブの [フィーチャ] グループで、[削除] をクリックします。
- 必要に応じて、表示された [削除] 確認ウィンドウで [はい] をクリックします。
ポリゴンが削除されました。
- リボン上の [編集] タブにある [編集の管理] グループで、[保存] をクリックして、編集内容を保存します。 [編集の保存] ウィンドウでプロンプトが表示されたら、[はい] をクリックします。
- リボンの [マップ] タブにある [ナビゲーション] グループで、[マップ操作] をクリックして、選択モードを終了します。
次に、ポリゴンの形状を修正します。 場合によっては、以下の画像に見られるような、LIDAR データからの建物フットプリントの正常な抽出を妨げる問題があります。
この問題は主に、金網で保護されている空調機器が原因であると考えられます。 そのようなポリゴンを探して、形状を変更します。
- 以下の画像でハイライト表示されている建物を探して拡大します。
- リボンの [編集] タブの [ツール] グループで、[形状変更] をクリックします。
- マップで、ポリゴンをクリックして選択します。
- 以下の画像で印が付いている 3 つのポイントをクリックして、ポリゴンの形状を変更します。
注意:
正しく形状変更するには、ポリゴン内をクリックするか、ポリゴンのエッジを直接クリックする必要があります。 ポインターが既存のフィーチャに確実にスナップするためにスナップをオンにします。
- 最後のポイントをダブルクリックします。
ポリゴンの形状が変更されました。
ヒント:
操作を誤った場合は、Ctrl + Z キーを押して再度操作します。
- 編集内容を保存します。
- [編集] タブの [選択] グループで、[選択解除] をクリックしてポリゴンの選択を解除します。
- [マップ] タブにある [ナビゲーション] グループで、[マップ操作] をクリックして、選択モードを終了します。
- [フィーチャの修正] ウィンドウを閉じます。
引き続き他の建物の微調整もできますが、このチュートリアルでは、次の手順に進みます。
抽出された建物フットプリントの評価
建物フットプリントが正常に抽出された結果を評価するため、建物フットプリントを含む地形図ベースマップと比較します。ただし、これらのフットプリントは、分析目的では使用できません。
- [コンテンツ] ウィンドウで、[Tuborg_Havn_Ortho_Photo.tif] レイヤーをオフにし、ベースマップを表示します。
- 必要に応じて、[Building_Footprints] レイヤーをクリックして選択します。 リボンの [フィーチャ レイヤー] タブの [比較] グループで、[スワイプ] を選択します。
- 前後にスワイプして、[Building_Footprints] レイヤーをベースマップのフットプリントと比較します。
全体としては、LIDAR から取得した建物フットプリントは、ベースマップの建物フットプリントと非常に似ています。これは、ここで完了したワークフローが適切であったことを示しています。
- リボンで [マップ] タブの [ナビゲーション] グループで、[マップ操作] をクリックし、[スワイプ] ツールを非アクティブにします。
- プロジェクトを保存します。
ジオプロセシング ツールと編集ツールを使用して、建物フィーチャクラスを生成し、クリーンナップを行いました。
写実的な 3D 建物の抽出
LAS データセットから、写実的な 3D 建物を抽出します。 そのため、これまでに作成したすべてのエレメント (建物に分類された LAS ポイント、DEM ラスター、および建物フットプリント) を使用します。 抽出結果として、複雑な 3D ベクター フィーチャを格納できる形式のマルチパッチ レイヤーが得られます。
マルチパッチ レイヤーの作成
建物の 3D 表示が可能な 3D シーンでこのセクションを実行します。まず初めに、[Building_Footprints] レイヤーをコピーします。
- [コンテンツ] ウィンドウで [Building_Footprints] レイヤーを右クリックし、[コピー] をクリックします。
- [Scene] タブをクリックしてアクティブにします。
- [コンテンツ] ウィンドウで、[Scene] を右クリックして [貼り付け] をクリックします。
[Building_Footprints] レイヤーがシーンに追加されます。 続いて、3D 建物を生成します。
- [コンテンツ] ウィンドウで、[Tuborg_Havn.lasd] レイヤーがオフになっていることを確認します。
- [ジオプロセシング] ウィンドウを元に戻し、必要であれば、[戻る] ボタンをクリックします。
- [ジオプロセシング] ウィンドウで、[LAS 建物マルチパッチ] ツールを探して開きます。
- [LAS 建物マルチパッチ] ツール ウィンドウで、次のパラメーター値を設定します。
- [入力 LAS データセット] で [Tuborg_Havn.lasd] レイヤーを選択します。
- [入力フィーチャ] で [Building_Footprints] を選択します。
- [LAS 屋根ポイントの選択] で、[建物に分類されたポイント] が選択されていることを確認します。
- [地表の高さ] は、[ラスタ レイヤー] のままで [Tuborg_Havn_DEM] を選択します。
- [出力マルチパッチ フィーチャクラス] で、デフォルト名を使用します。
- [単純化許容値] に「0.5」と入力し、単位が [メートル] に設定されていることを確認します。
[単純化許容値] には、DEM ラスター セルのサイズを選択することをお勧めします。
- [実行] をクリックします。
出力レイヤーが表示されます。
- [コンテンツ] ウィンドウで、新しいレイヤーの名前を [Multipatch_3D_Buildings] に変更し、そのレイヤーとベースマップ以外のレイヤーをすべてオフにします。
マルチパッチ 3D 建物の色を変更して、外観を整えることができます。
- [コンテンツ] ウィンドウで、[Multipatch_3D_Buildings] レイヤーのシンボルを右クリックし、[サハラ サンド] 色を選択します。
建物の色が変わりました。
- 3D 建物を調べ、多くの建築の詳細部分が写実的に描かれている点に注目してください。
注意:
このワークフローでは、常にカートグラフィ面で満足のいく結果が生成されるとは限りませんが、位置の精度が確保されます。 この種の 3D レイヤーは、見通し線や可視領域など、特定の解析を実行する上で特に有用です。
一部の建物が、地表からやや浮き上がっています。 これは、地表レベルの参照として [Tuborg_Havn_DEM] レイヤーを使用して生成されたためです。 [Tuborg_Havn_DEM] は、非常に正確な LIDAR データから生成されているので、描写がやや異なり、地表を 3D でモデル化するためにシーンで使用されるデフォルトの [WorldElevation3D/Terrain3D] レイヤーより詳細にテレインが描かれています。
[Tuborg_Havn_DEM] レイヤーをシーンの標高サーフェスに追加し、シーンがこのレイヤーを使用して、ツボルグ港近辺のテレインをより正確にモデル化できるようにします。
- [コンテンツ] ウィンドウの [標高サーフェス] で [Ground] を右クリックし、[標高ソース レイヤーを追加] を選択します。
- [標高ソース レイヤーを追加] ウィンドウで [データベース] を開き、[Building_Extraction.gbd] を開いて [tuborg_lasda] をクリックします。 [OK] をクリックします。
DEM が新しい地表の標高サーフェスとして追加されたので、これを使用するようシーンが更新されます。 一部の建物の下部のわずかな隙間が解消されます。
- 新しい標高ソースの名前を [Tuborg DEM] に変更します。
- プロジェクトを保存します。
3D 建物マルチパッチ レイヤーは、さまざまな方法で使用できます。 シーンに追加して、都市計画活動をサポートすることができます。 洪水予測解析、日影影響解析、可視領域および見通し線解析などの実行に使用することもできます。
このチュートリアルで説明したワークフローでは、LIDAR を使用した作業方法、および 3D 建物抽出プロセスの多くの構成要素について理解を深めることに重点を置きました。 このワークフローの詳細と大規模な 3D ベースマップを生成する方法については、ArcGIS ラボ「Working with the ArcGIS Solution for 3D Buildings」をご参照ください。
このチュートリアルでは、LIDAR データから 3D 建物を抽出しました。 LAS データセットの作成と分類、DEM ラスターの生成、建物フットプリント ベクター レイヤーの抽出を行いました。 最後に、それらのエレメントをすべて使用して、3D 建物マルチパッチ ベクター レイヤーを作成しました。
他のチュートリアルについては、チュートリアル ギャラリーをご覧ください。
