陥没のない DEM の作成
流向の過程を表す入力データとしては、陥没のないデジタル標高モデル（DEM）を使用することを推奨します。シンク（窪地）があると、流向ラスタにエラーが含まれることがあります。場合によっては、データ内のシンクが正当な場合もあります。地形の形状を十分理解して、地表の実際のシンクを表すフィーチャと、データ上の単なるエラーであるフィーチャを見分けることが重要です。
シンクの抽出
シンクは、[窪地の抽出（Sink）] ツールを使用して見つけることができます。このツールには、[流向ラスタの作成（Flow Direction）] ツールで作成した方向ラスタが必要です。ツールの結果は、データに存在するシンクを示すラスタです。結果によって、シンクを埋めたり、出力データを使用して平滑化の制限値を決定したりできます。シンクは、[サーフェスの平滑化 (Fill)] ツールを使用して埋めることができます。[窪地の抽出（Sink）] ツールからの出力データを使用して平滑化の制限値を決定する方法については、このトピック内の「シンクの深さの検出」（後述）をご参照ください。
 |
シンクのプロファイル表示 |
シンクの平滑化
[サーフェスの平滑化 (Fill)] ツールは、前述の水文解析ツールをはじめとする ArcGIS Spatial Analyst のさまざまなツールを使用して、陥没のない DEM を作成します。このツールには、入力サーフェス、平滑化の制限値、および出力ラスタが必要です。シンクを平滑化するときには、流出点、集水域の境界に沿った最小の標高まで平滑化されます。
陥没のない DEM を作成する際にシンクを特定して除去する操作は、反復的なプロセスです。シンクが平滑化されると、平滑化された領域の境界でさらに新しいシンクが作成されることがあり、このシンクも平滑化する必要があります。データの容量が大きい DEM やシンクを多数持つ DEM では、この処理に数分から数時間掛かることがあります。
シンクの深さの検出
1 つのシンク、またはシンクのグループの深さを知っておくと便利です。この情報があれば、[サーフェスの平滑化 (Fill)] ツールの適切な Z 制限値の決定、データに存在するエラーの種類の理解、およびシンクが正当な地表面の形状であるかどうかの判断が可能になります。以下のステップで、シンクの深さを検出する一般的なプロセスの概要を示します。
- [サーフェスの平滑化 (Fill)] ツールを実行してシンクの深さの値を持つラスタを作成し、ラスタ内のシンクを見つけます。
- [集水域ラスタの作成（Watershed）] ツールで、流出点の入力データとして標高ラスタの流向、および [サーフェスの平滑化 (Fill)] ツールの出力データを使用して、各シンクに影響するエリアのラスタを作成します。
- Python の [マップ代数演算] で、[最小値] オプションを指定して [ゾーン統計（Zonal Statistics）] ツールを使用し、各シンクの集水域での最小標高値を持つラスタを作成します。
sink_min = ZonalStatistics(sink_areas, "Value", elev_ras, "Minimum")
sink_areas 入力は、[集水域ラスタの作成（Watershed）] ツールの出力です。 - [ゾーン フィル（Zonal Fill）] ツールで、集水域ごとの境界ごとに最も低い標高値を持つラスタを作成します（これは、周縁の高さにまで平滑化されると、流れが流域から流出する標高に相当します）。
sink_max = ZonalFill(sink_areas, elev_ras)
- [Minus] ツールを使用して、最大値から最小値を減算して、深さを再び調べます。
sink_depth = Minus(sink_max, sink_min)
上記の手順は、[窪地の抽出（Sink）]、[集水域ラスタの作成（Watershed）]、[ゾーン統計（Zonal Statistics）]、[ゾーン フィル（Zonal Fill）]、および [Minus] のジオプロセシング ツールを使用しても実行できます。