窪地の抽出の詳細

Spatial Analyst のライセンスで利用可能。

窪地とは、流向として、流向ラスターの 8 つの有効な値のいずれも割り当てることができない 1 つのセルまたは空間的に接続している一連のセルのことです。 これは、すべての隣接セルが処理セルよりも高い場合、または 2 つのセルが相互に流れ込んでループが形成されている場合に発生します。

窪地は未定義の流向があると見なされ、可能な方向すべての合計値が割り当てられます。 たとえば、最も傾斜が急な下りで、右 (1) と左 (16) の流向が同じ場合、そのセルの流向の値は 17 になります。

流向の正確な表現を作成し、累積流量を求めるには、窪地のないデータセットを使用するのが最適です。 すべての窪地を削除する処理が行われたデジタル標高モデル (DEM) は、窪地のない DEM と呼ばれます。

標高データにおける窪地は、ほとんどの場合データのエラーによるものです。 これらのエラーは、サンプリング効果や、標高を整数値に丸めたことによって発生することがよくあります。 本来、氷河またはカルスト地域を除き、セル サイズが 10 メートル以上の標高データで窪地が発生することはまれ (Mark 1988) なので、一般にエラーと見なすことができます。 多くの場合、セル サイズが大きくなると、データセット内の窪地の数も増えます。

窪地のない DEM を作成する際の窪地の特定と削除は、反復的なプロセスです。 窪地を平滑化すると、平滑化された領域の境界に新しい窪地が生まれ、それを平滑化する必要が生じることがあります。 [ゾーン フィル (Zonal Fill)] ツールに関するトピックは、窪地のない DEM の作成方法を説明する上で役立ちます。

窪地または窪地のグループの深さを知ると役に立つことがあります。 この情報があれば、[サーフェスの平滑化 (Fill)] ツールの [Z 制限] の設定に適した値を決定したり、データに存在するエラーの種類を理解したりすることが可能になります。

窪地の断面表示
窪地の断面表示

ジオプロセシング ツールを使用して窪地を平滑化する手順は以下のとおりです。

  1. [窪地の抽出 (Sink)] ツールを使用して、深さでコード化された窪地のラスターを作成します。

    [入力流向ラスター]: flowdir

    [出力ラスター]: sinks

  2. [集水域ラスターの作成 (Watershed)] ツールを使用して、各窪地の流出寄与域のラスターを作成します。

    [入力流向ラスター]: flowdir

    [流出点データとして使用する入力ラスター、またはフィーチャ]: sinks

    [出力ラスター]: sink_areas

  3. [ゾーン統計 (Zonal Statistics)] ツールで最小統計を使用し、各窪地の集水域での最小標高値を持つラスターを作成します。

    [入力ラスター、またはフィーチャ ゾーン データ]: sink_areas

    [ゾーン フィールド]: Value

    [入力値ラスター]: elevation

    [出力ラスター]: sink_min

    [統計の種類]: 最小

  4. [ゾーン フィル (Zonal Fill)] ツールを使用し、各窪地の集水域での最大標高値を持つラスターを作成します。

    [入力ゾーン ラスター]: sink_areas

    [入力加重ラスター]: elevation

    [出力ラスター]: sink_max

  5. [Minus] ツールを使用して、最大値から最小値を減算して、深さを求めます。

    [入力ラスター 1]: sink_max

    [入力ラスター 2]: sink_min

    [出力ラスター]: sink_depth

これは ArcPy でマップ代数演算を使用するのと同じ手順です。

sinks = Sink(flowdir)
sink_areas = Watershed(flowdir, sinks)
sink_min = ZonalStatistics(sink_areas, "Value", elevation, "Minimum")
sink_max = ZonalFill(sink_areas, elevation)
sink_depth = Minus(sink_max, sink_min)

参考情報

Mark, D. M. 1988. "Network Models in Geomorphology." Modelling Geomorphological Systems, ed. M. G. Anderson. New York: John Wiley. 73-97.

関連トピック


このトピックの内容
  1. 参考情報