(57)【特許請求の範囲】
【請求項１】
地形モデルを用いて地形の種別を分類するシステムであって、
前記地形モデルは、所定領域を平面分割して得られる多数の単位区画と、該単位区画が具備する標高値と、を含んで構成され、
単位区画に対して、平均傾斜量と傾斜標準偏差を求める地形量算出手段と、
オペレータの操作により、前記所定領域のうち地形の特徴に着目してサンプル領域を抽出するサンプル領域抽出手段と、
オペレータの操作により、前記サンプル領域の境界と交差する注目線を設定する注目線設定手段と、
前記注目線上にある単位区画の前記平均傾斜量及び前記傾斜標準偏差を、平均傾斜量軸及び傾斜標準偏差軸からなる２軸平面に配置して散布図を作成する散布図作成手段と、
前記散布図に基づいて、前記平均傾斜量及び前記傾斜標準偏差を変数とする閾値関数を設定する地形分類条件設定手段と、
前記平均傾斜量及び前記傾斜標準偏差の組み合わせが前記閾値関数の上方又は下方となる単位区画を、特徴単位区画として分類する特徴地形分類手段と、
を備え、
前記平均傾斜量は、単位区画と該単位区画の周辺にある複数の単位区画によってそれぞれ算出される周辺傾斜量に基づいて求められ、
前記傾斜標準偏差は、複数の前記周辺傾斜量に基づいて求められる、
ことを特徴とする地形分類システム。
【請求項２】
前記特徴地形分類手段によって分類された特徴単位区画が、複数集合した領域を特徴領域とし、該特徴領域の境界を遷緩線又は遷急線として設定する地形変化線設定手段を、
さらに備えたことを特徴とする請求項１記載の地形分類システム。
【請求項３】
請求項１又は請求項２記載の地形分類システムを用いて地形の種別を分類する方法であって、
前記地形量算出手段が、単位区画に対して前記平均傾斜量と前記傾斜標準偏差を求める地形量算出工程と、
オペレータが前記サンプル領域抽出手段を操作することにより、前記サンプル領域を抽出するサンプル領域抽出工程と、
オペレータが前記注目線設定手段を操作することにより、前記注目線を設定する注目線設定工程と、
前記散布図作成手段が、前記散布図を作成する散布図作成工程と、
前記地形分類条件設定手段が、前記閾値関数を設定する地形分類条件設定工程と、
前記特徴地形分類手段が、前記平均傾斜量及び前記傾斜標準偏差の組み合わせが前記閾値関数の上方又は下方となる単位区画を、前記特徴単位区画として分類する特徴地形分類工程と、を備えた、
ことを特徴とする地形分類方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本願発明は、地形モデルを用いて特徴ある地形を分類する技術に関するものであり、より具体的には、地形モデルを構成する単位区画ごとに平均傾斜量と傾斜量の標準偏差（以下、「傾斜標準偏差」という。）を求め、その平均傾斜量と傾斜標準偏差をもとに地形を分類するシステムと方法に関するものである。
【背景技術】
【０００２】
我が国は、国土に占める山地の割合が多く、しかも雨が多いことから、地すべりや斜面崩壊といった斜面災害による被害を受けやすい。特に、ひとたび地すべりや深層崩壊が生ずると、甚大な被害を伴うこともある。このような斜面災害による被害を受けないように（あるいは軽減するように）、これまで原因となる斜面で対策工が施され、災害時における避難計画が策定されてきた。
【０００３】
ところで、斜面対策工を行うにしても、住民の避難計画を策定するにしても、災害原因となる斜面を特定する必要がある。従来、地すべりが生じそうな地形（以下、「地すべり地形」という。）や、斜面崩壊が生じそうな地形（以下、「崩壊地形」という。）は、専門技術を有する技術者によって抽出されてきた。その際、膨大な範囲の地形から直ちに地すべり地形や崩壊地形を特定することは難しく、まずは地形的な特徴を糸口として候補となる地形を挙げ、その候補地形を詳細に検討することで地すべり地形や崩壊地形を抽出するのが一般的であった。
【０００４】
地すべり地形や崩壊地形の候補を挙げるためには、斜面勾配が急激に変化する地形に着目するのが有効であり、このような地形としては、斜面上方から見て緩斜面から急斜面に変化する「遷急線」と、急斜面から緩斜面に変化する「遷緩線」が知られている。またある断面に着目し、その断面のうち遷急線を構成する点を「遷急点」、遷緩線を構成する点を「遷緩点」と呼ぶこともある。なおここでは、遷急線と遷緩線の総称を「地形変化線」ということとし、遷急点と遷緩点の総称を「地形変化点」ということとする。
【０００５】
地すべり地形や崩壊地形と同様、地形変化線や地形変化点もまた専門技術者によってこれまで抽出されてきた。専門技術者が有する知識とともに過去の経験にしたがって地形変化線等を抽出するわけであるが、そのため技術者によってその結果はそれぞれ異なることもあり、つまり技術者によって地形変化線等の適否判断が分かれていた。既述のとおり地形変化線等の抽出は、地すべり地形等の候補地形を挙げるために行うものであり、地形変化線等の抽出漏れはすなわち地すべり地形等の抽出漏れを意味する。
【０００６】
このような人為的なばらつきを回避するため、言い換えれば客観的な結果を得る目的で、機械的に地形の特徴を抽出する取り組みはこれまでも行われてきた。特許文献１では、地形図に地すべりブロック（平面形状）を与えることで、その地すべり方向を自動的に求める技術について提案している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００５−１６４４２１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ところが特許文献１を含めこれまで、遷急線や遷緩線といった地形変化線（あるいは地形変化点）を機械的（自動的）に抽出する技術が提示されることはなく、関係者の間ではこのような技術が強く要望されるところであった。
【０００９】
本願発明の課題は、従来が抱える問題を解決することであり、すなわち技術者によるばらつきを回避し、つまり客観的な判断によって地形変化線を抽出するため、特徴ある地形を分類することができる地形分類システム、及び地形分類方法を提供することである。
【課題を解決するための手段】
【００１０】
本願発明は、地形モデルを構成する単位区画ごとに求められる平均傾斜量と傾斜標準偏差に基づいて特徴ある地形を分類する、という点に着目したものであり、従来にはなかった発想に基づいてなされた発明である。
【００１１】
本願発明の地形分類システムは、地形モデルを用いて地形の種別を分類するシステムであり、地形量算出手段と、サンプル領域抽出手段、注目線設定手段、散布図作成手段、地形分類条件設定手段、特徴地形分類手段を備えたものである。ここで地形モデルとは、所定領域を平面分割して得られる多数の単位区画と、これら単位区画がそれぞれ具備する標高値を含んで構成されるものである。地形量算出手段は、単位区画に対して平均傾斜量と傾斜標準偏差を求めるもので、サンプル領域抽出手段は、オペレータの操作により所定領域のうち地形の特徴に着目してサンプル領域を抽出するものである。注目線設定手段は、オペレータの操作によりサンプル領域の境界と交差する注目線を設定するもので、散布図作成手段は、サンプル領域内にある単位区画の平均傾斜量と傾斜標準偏差を、平均傾斜量軸及び傾斜標準偏差軸からなる２軸平面に配置して散布図を作成するものである。また地形分類条件設定手段は、散布図に基づいて平均傾斜量と傾斜標準偏差を変数とする閾値関数を設定するもので、特徴地形分類手段は、平均傾斜量と傾斜標準偏差の組み合わせが閾値関数の上方又は下方となる単位区画を、特徴単位区画として分類するものである。なお、平均傾斜量は、単位区画とその単位区画の周辺にある複数の単位区画によってそれぞれ算出されるもので、傾斜標準偏差は、複数の周辺傾斜量に基づいて求められるものである。
【００１２】
本願発明の地形分類システムは、地形変化線設定手段をさらに備えたものとすることもできる。地形変化線設定手段は、特徴地形分類手段によって分類された特徴単位区画が複数集合した領域を特徴領域とし、この特徴領域の境界を遷緩線又は遷急線として設定するものである。
【００１３】
本願発明の地形分類方法は、地形モデルを用いて地形の種別を分類する方法であり、地形量算出工程と、サンプル領域抽出工程、注目線設定工程、散布図作成工程、地形分類条件設定工程、特徴地形分類工程を備えたものである。地形量算出工程では、単位区画に対して平均傾斜量と傾斜標準偏差を求め、サンプル領域抽出工程では、オペレータの操作により所定領域のうち地形の特徴に着目してサンプル領域を抽出する。注目線設定工程では、オペレータの操作によりサンプル領域の境界と交差する注目線を設定し、散布図作成工程では、サンプル領域内にある単位区画の平均傾斜量と傾斜標準偏差を、平均傾斜量軸及び傾斜標準偏差軸からなる２軸平面に配置して散布図を作成する。また地形分類条件設定工程では、散布図に基づいて平均傾斜量と傾斜標準偏差を変数とする閾値関数を設定し、特徴地形分類工程では、平均傾斜量と傾斜標準偏差の組み合わせが閾値関数の上方又は下方となる単位区画を、特徴単位区画として分類する。なお、平均傾斜量は、単位区画とその単位区画の周辺にある複数の単位区画によってそれぞれ算出されるもので、傾斜標準偏差は、複数の周辺傾斜量に基づいて求められるものである。
【発明の効果】
【００１４】
本願発明の地形分類システム、及び地形分類方法には、次のような効果がある。
（１）技術者によるばらつきを回避し、客観的な判断によって、特徴ある地形を分類することができる。その結果、客観的に地形変化線を抽出することができ、ひいては地すべり地形や崩壊地形の抽出漏れを防ぐことができる。
（２）人による作業を大幅に省略できることから、作業コストを低減することができるとともに、地すべり地形等をこれまでより迅速に抽出することができる。
【図面の簡単な説明】
【００１５】
【図１】本願発明の地形分類システムを示すブロック図。
【図２】地形変化点を抽出するための条件を設定するための主な処理（工程）を示すフロー図。
【図３】条件設定手段を示すブロック図。
【図４】地形モデルと単位区画を説明するモデル図。
【図５】サンプル領域に対して設定された注目線を説明するモデル図。
【図６】（ａ）１次傾向面解析によるＸ軸方向の傾斜量の算出手法を説明するモデル図、（ｂ）は１次傾向面解析によるＹ軸方向の傾斜量の算出手法を説明するモデル図。
【図７】（ａ）は図５に示すサンプル領域に対して設定された注目線Ｌｉｎｅ１における散布図、（ｂ）は図５に示すサンプル領域に対して設定された注目線Ｌｉｎｅ２における散布図、（ｃ）は図５に示すサンプル領域に対して設定された注目線Ｌｉｎｅ３における散布図、（ｄ）は図５に示すサンプル領域に対して設定された注目線Ｌｉｎｅ４における散布図、（ｅ）は図５に示すサンプル領域に対して設定された注目線Ｌｉｎｅ５における散布図。
【図８】所定領域に対して特徴単位区画を抽出するための主な処理（工程）を示すフロー図。
【図９】特徴領域抽出手段を示すブロック図。
【発明を実施するための形態】
【００１６】
本願発明の地形分類システム、及び地形分類方法の一例を、図を参照しながら説明する。
【００１７】
１．全体概要
図１は、本願発明の地形分類システム１００を示すブロック図である。この図に示すように地形分類システム１００は、条件設定手段１１０と特徴領域抽出手段１２０を含んで構成され、プリンタやディスプレイといった出力手段１３０を含むこともできる。条件設定手段１１０は、特徴ある地形領域（以下、「特徴領域」という。）を構成する単位区画（以下、「特徴単位区画」という。）を抽出するための条件（以下、「地形分類条件」という。）を設定するもので、この地形分類条件に基づいて特徴領域抽出手段１２０が特徴単位区画を抽出する。なお本願発明の地形分類システム１００は、プログラムを実行するコンピュータを利用して構成するとよい。以下、条件設定手段１１０と特徴領域抽出手段１２０について詳しく説明する。
【００１８】
２．条件設定手段
図２は、地形分類条件を設定するための主な処理（工程）を示すフロー図であり、図３は、条件設定手段１１０を示すブロック図である。なお、図２の中央の列には実施する処理（工程）を示しており、左列にはその処理（工程）に必要な入力情報を、右列にはその処理（工程）から生まれる出力情報を示している。以下、これらの図を参照しながら条件設定手段１１０について説明する。
【００１９】
（地形モデル）
図４は、地形モデルと単位区画を説明するモデル図である。この図に示すように地形モデルＴＭは、所定領域ＡＲを平面分割した網目状のモデルであり、単位区画ＢＬは平面分割された結果得られる小領域のことである。言い換えれば、地形モデルＴＭは多数の単位区画ＢＬによって構成される。例えば、図４左側の破線領域は、２４個の単位区画ＢＬで構成されていることが分かる。なお、この図にも示すように、それぞれの単位区画ＢＬには識別番号（図ではＢＬ３２３〜ＢＬ８２４）が付与されることが多い。また地形モデルＴＭは、平面座標（あるいは緯度経度）が付与されており、通常は単位区画ＢＬの格子点に平面座標等が付与されている。
【００２０】
このような地形モデルＴＭとして代表的なのが、ＤＥＭ（Ｄｉｇｉｔａｌ  Ｅｌｅｖａｔｉｏｎ  Ｍｏｄｅｌ）やＤＳＭ（Ｄｉｇｉｔａｌ  Ｓｕｒｆａｃｅ  Ｍｏｄｅｌ）である。地表モデルとも言われるＤＥＭは、地形モデルＴＭを構成する個々の単位区画ＢＬに標高値が付与され、一方表層モデルとも言われるＤＳＭは、地形モデルＴＭを構成する個々の単位区画ＢＬに、被覆物等の標高値が付与される。これら標高値は、単位区画ＢＬのうち中心点や格子点に付与されるのが一般的である。なお図４では所定領域ＡＲが略正方格子状に区切られ、単位区画ＢＬの形状も略正方形となっているが、これに限らず、単位区画ＢＬの形状を長方形やひし形、あるいは長円形など任意の形状とすることが可能で、さらにそれぞれの単位区画ＢＬの形状や大きさを変えることもできる。
【００２１】
地形モデルＴＭは、図３に示す地形モデル記憶手段１１１に記憶され、地形モデル読出し手段１１２によって読み出される（Ｓｔｅｐ１０１：図２）。
【００２２】
（サンプル領域）
ここでサンプル領域とは、地形（特に勾配の変化）に特徴がある領域であり、操作者（オペレータ）によって選定される。サンプル領域としては、例えば、緩傾斜地や急傾斜地が代表的に選定され、さらにこれらの特徴に加え規模（面積）に応じて細分化してもよい。なおサンプル領域は、特徴領域を抽出しようとする所定領域ＡＲに対して１種類のみ抽出してもよいし、２種類以上を抽出してもよい。２種類以上のサンプル領域を抽出する場合、それぞれのサンプル領域に設定される注目線（後述する）の数だけ、地形分類条件が設定される。
【００２３】
サンプル領域は、図３に示すサンプル領域抽出手段１１３をオペレータが操作することで抽出される（Ｓｔｅｐ１０２：図２）。具体的には、表示手段１３０に表示した地形モデルＴＭを確認しながら、オペレータがマウスなどポインティングデバイスを利用して所定の領域を指定するとよい。
【００２４】
サンプル領域が抽出されると、注目線が設定される（Ｓｔｅｐ１０３：図２）。図５は、サンプル領域ＳＡに対して設定された注目線を説明するモデル図であり、この図に示すように注目線は、サンプル領域ＳＡの境界と交差するように設定される。なお図５では、１つのサンプル領域ＳＡに対しての５つの注目線（Ｌｉｎｅ１〜Ｌｉｎｅ５）が示されているが、注目線の数は任意に設定でき、またサンプル領域ＳＡの境界と注目線の交差角も任意に選択することができる。この注目線は、図３に示す注目線設定手段１１４をオペレータが操作することで設定され、具体的には、表示手段１３０に表示した地形モデルＴＭとサンプル領域ＳＡを確認しながら、オペレータがマウスなどポインティングデバイスを利用して注目線を指定するとよい。
【００２５】
（平均傾斜量）
サンプル領域ＳＡに対して注目線が設定されると、その注目線上にある単位区画ＢＬに対して平均傾斜量Ｉを求める。通常、注目線上には複数の単位区画ＢＬがあり、これらの単位区画ＢＬに対してそれぞれ平均傾斜量Ｉを求める。また、設定した注目線全てにかかる単位区画ＢＬに対して平均傾斜量Ｉを求めるとよい。なお、注目線上の単位区画ＢＬには、その注目線の一部を含む単位区画ＢＬのほか、注目線からの所定距離以内にある単位区画ＢＬを含むこともできる。
【００２６】
平均傾斜量Ｉは、１次傾向面解析による傾斜量として算出され、具体的には算出対象の単位区画ＢＬと、その周辺にある複数の単位区画ＢＬを用いて算出される。ここでは混乱を避けるため、平均傾斜量Ｉを求める対象となる単位区画ＢＬを「対象単位区画ＢＬａ」と、対象単位区画ＢＬａの周辺にある単位区画ＢＬを「周辺単位区画ＢＬｓ」ということとする。
【００２７】
図６は、１次傾向面解析による傾斜量の算出手法を説明するモデル図であり、（ａ）はＸ軸方向の傾斜量、（ｂ）はＹ軸方向の傾斜量を示している。この図では、Ｘ軸とＹ軸からなる平面上に配置された地形モデルＴＭの一部を示しており、中央の網掛けされたものを対象単位区画ＢＬａとし、その周辺にあるものを周辺単位区画ＢＬｓとしている。
【００２８】
対象単位区画ＢＬａの平均傾斜量Ｉを算出するに当たって、まず算出に用いる周辺単位区画ＢＬｓの範囲（以下、「窓領域」という。）を設定する。例えば図６では、対象単位区画ＢＬａを中心に５（Ｘ軸方向）×５（Ｙ軸方向）を窓領域とし、計２４個の周辺単位区画ＢＬｓが設定されている。なお窓領域は、あらかじめ定められた範囲（Ｎｘ×Ｎｙ）に基づいて自動生成することもできるし、その都度範囲（Ｎｘ×Ｎｙ）を入力し、その範囲に基づいて設定することもできる。
【００２９】
窓領域が設定できると、対象単位区画ＢＬａと周辺単位区画ＢＬｓを用いて「周辺傾斜量」が算出される。この周辺傾斜量の算出には、軸方向に沿った距離（以下、「軸距離」という。）が用いられる。例えば図６（ａ）では、対象単位区画ＢＬａ（中心）から周辺単位区画ＢＬｓ（中心）までの距離のＸ軸成分（Ｘ軸に投影した長さ）がＸ軸距離であり、周辺単位区画ＢＬｓ１１〜ＢＬｓ５１と周辺単位区画ＢＬｓ１５〜ＢＬｓ５５のＸ軸距離はＬｖ１であり、周辺単位区画ＢＬｓ１２〜ＢＬｓ５２と周辺単位区画ＢＬｓ１４〜ＢＬｓ５４のＸ軸距離はＬｖ２である。なお、周辺単位区画ＢＬｓ１３〜ＢＬｓ５３のＸ軸距離は０となるため、Ｘ軸方向の周辺傾斜量の算出にはこれらの周辺単位区画ＢＬｓは用いない。
【００３０】
同様に図６（ｂ）では、対象単位区画ＢＬａ（中心）から周辺単位区画ＢＬｓ（中心）までの距離のＹ軸成分（Ｙ軸に投影した長さ）がＹ軸距離であり、周辺単位区画ＢＬｓ１１〜ＢＬｓ１５と周辺単位区画ＢＬｓ５１〜ＢＬｓ５５のＹ軸距離はＬｈ１であり、周辺単位区画ＢＬｓ２１〜ＢＬｓ２５と周辺単位区画ＢＬｓ４１〜ＢＬｓ４５のＹ軸距離はＬｈ２である。なお、周辺単位区画ＢＬｓ３１〜周辺単位区画ＢＬｓ３５のＹ軸距離は０となるため、Ｙ軸方向の周辺傾斜量の算出にはこれらの周辺単位区画ＢＬｓは用いない。
【００３１】
対象単位区画ＢＬａと周辺単位区画ＢＬｓの標高値の差を、図６に示す軸距離で除したものが「成分傾斜量」として算出される。なお、同一の周辺単位区画ＢＬｓに対して、Ｘ軸距離によって算出されるＸ成分傾斜量と、Ｙ軸距離によって算出されるＹ成分傾斜量の２つの成分傾斜量が得られる（ただし、周辺単位区画ＢＬｓ１３〜ＢＬｓ５３と周辺単位区画ＢＬｓ３１〜ＢＬｓ３５は１つの成分傾斜量のみ得られる）。そこで、Ｘ成分傾斜量とＹ成分傾斜量の合成値（２乗和の平方根）を周辺単位区画ＢＬｓの周辺傾斜量とするとよい。そして、対象単位区画ＢＬａに対して求められた複数（図６では２４個）の周辺傾斜量が得られると、その対象単位区画ＢＬａの平均傾斜量Ｉを算出する（Ｓｔｅｐ１０４：図２）。この平均傾斜量Ｉは、周辺傾斜量に基づいて求められる値であり、算術平均としたり、相加平均としたり、軸距離によって重み付けした加重平均とするなど、種々の手法で算出することができる。
【００３２】
（傾斜標準偏差）
周辺傾斜量が得られると、平均傾斜量Ｉのほかに傾斜標準偏差σも算出する（Ｓｔｅｐ１０５：図２）。この傾斜標準偏差σは、対象単位区画ＢＬａに対して求められた複数（図６では２４個）の周辺傾斜量の標準偏差として計算された値である。
【００３３】
（地形量）
平均傾斜量Ｉと傾斜標準偏差σ（以下、これらを総称して「地形量」という。）は、図３に示す地形量算出手段１１５によって算出される。なお地形量は、設定された全ての注目線上にあるできるだけ多くの（全てでもよい）単位区画ＢＬに対して繰り返し算出される（図２）。
【００３４】
（散布図）
サンプル領域内にある単位区画ＢＬの地形量が得られると、傾斜標準偏差σを縦軸、平均傾斜量Ｉを横軸（あるいは、平均傾斜量Ｉを縦軸、傾斜標準偏差σを横軸）とする座標軸平面の上に、個々の単位区画ＢＬの地形量を配置した散布図を作成する（Ｓｔｅｐ１０６：図２）。いわば、平均傾斜量Ｉと傾斜標準偏差σの組み合わせを座標として、２軸座標平面上にプロットしたものが散布図である。なお、この散布図は注目線ごとに作成することができる。図７は、図５に示すサンプル領域ＳＡに対して設定された５つの注目線ごとに作成された散布図であり、（ａ）は注目線Ｌｉｎｅ１における散布図、（ｂ）は注目線Ｌｉｎｅ２における散布図、（ｃ）は注目線Ｌｉｎｅ３における散布図、（ｄ）は注目線Ｌｉｎｅ４における散布図、（ｅ）は注目線Ｌｉｎｅ５における散布図である。散布図は、図３に示す散布図作成手段１１６によって作成される。
【００３５】
（閾値関数）
散布図が得られると、閾値関数が設定される（Ｓｔｅｐ１０７：図２）。この閾値関数は、平均傾斜量Ｉと傾斜標準偏差σを変数とする関数であり、散布図上に散布された点群（つまり、地形量）に基づいて設定される。したがって、平均傾斜量Ｉと傾斜標準偏差σからなる座標軸平面の上に、閾値関数を描くと連続した直線又は曲線として表される。なお図７では、閾値関数が全て直線として表されている。この閾値関数は、例えば専門の技術者が散布図上の点群を概ね上下（あるいは左右）に２分するような境界線を定め、この境界線を閾値関数として設定することができる。
【００３６】
閾値関数は、平均傾斜量Ｉのばらつきを考慮したうえで、地形的に特徴あるサンプル領域とそうでない地形を分ける境界線といえることから、閾値関数より上方（あるいは下方）の地形量を示す単位区画が特徴単位区画と推定される。なお、閾値関数より上方とするか下方とするかは、サンプル領域の特徴（緩傾斜や急傾斜）に応じて決定される。したがって、この閾値関数が地形分類条件であり、図３に示す地形分類条件設定手段１１７によって設定され、地形分類条件記憶手段１１８に記憶される。
【００３７】
３．特徴領域抽出手段
図８は、所定領域ＡＲに対して特徴単位区画を抽出するための主な処理（工程）を示すフロー図であり、図９は、特徴領域抽出手段１２０を示すブロック図である。なお、図８の中央の列には実施する処理（工程）を示しており、左列にはその処理（工程）に必要な入力情報を、右列にはその処理（工程）から生まれる出力情報を示している。以下、これらの図を参照しながら特徴領域抽出手段１２０について説明する。
【００３８】
（地形モデルの読出し）
まず、地形モデル記憶手段１１１に記憶された地形モデルＴＭを、地形モデル読出し手段１１２によって読み出し（Ｓｔｅｐ２０１：図８）、地形量（平均傾斜量Ｉと傾斜標準偏差σ）を算出するための窓領域を設定する。なお、ここで設定する窓領域は、地形分類条件を設定する際の窓領域（図６の５×５）と同じとすることもできるし、異なる大きさの領域を設定することもできる。
【００３９】
（地形量の算出）
次に、対象領域ＡＲ（地形モデルＴＭ）を構成するそれぞれの単位区画ＢＬに対して、順に所定の処理（工程）を繰り返し行う。はじめに、単位区画ＢＬの平均傾斜量Ｉが、１次傾向面解析による傾斜量として算出される（Ｓｔｅｐ２０２：図８）。この場合も条件設定手段１１０と同様、算出対象となる単位区画ＢＬを「対象単位区画ＢＬａ」と、対象単位区画ＢＬａの窓領域内にある単位区画ＢＬを「周辺単位区画ＢＬｓ」ということとする。既述のとおり平均傾斜量Ｉは、対象単位区画ＢＬａと周辺単位区画ＢＬｓに基づいて求められた周辺傾斜量から得られ、また既述のとおりこの周辺傾斜量は、軸距離を用いて算出された成分傾斜量に基づいて求められる。
【００４０】
平均傾斜量Ｉのほか、条件設定手段１１０と同様の手法で、すなわち周辺傾斜量に基づいて対象単位区画ＢＬａの傾斜標準偏差σを求める（Ｓｔｅｐ２０３：図８）。平均傾斜量Ｉと傾斜標準偏差σは、図９に示す地形量算出手段１２１によって算出される。
【００４１】
（特徴単位区画の選別）
対象単位区画ＢＬａの地形量が得られると、地形分類条件記憶手段１１８に記憶された地形分類条件が、図９に示す地形分類条件読出し手段１２２によって読み出される（Ｓｔｅｐ２０４：図８）。そして、対象単位区画ＢＬａの地形量と、読み出した地形分類条件（すなわち、閾値関数）を照らし合わせて、その対象単位区画ＢＬａが特徴単位区画に該当するか否か選別する（Ｓｔｅｐ２０５：図８）。具体的には、対象単位区画ＢＬａの平均傾斜量Ｉと傾斜標準偏差σの組み合わせが、閾置関数よりあらかじめ定めた一方（上方又は下方）側にあれば、この対象単位区画ＢＬａを特徴単位区画と判定する。この特徴単位区画の選別は、図９に示す特徴地形分類手段１２３によって行われる。
【００４２】
（地形変化線の設定）
特徴地形分類手段１２３によって特徴単位区画が選別されると、複数の特徴単位区画が集合した領域が特徴領域として生成され（Ｓｔｅｐ２０６：図８）、その特徴領域の境界線が遷緩線又は遷急線として抽出される（Ｓｔｅｐ２０７：図８）。なお、特徴領域は図９に示す特徴領域生成手段１２４によって生成され、遷緩線又は遷急線は図９に示す地形変化線設定手段１２５によって抽出される。
【産業上の利用可能性】
【００４３】
本願発明の地形分類システム、及び地形分類方法は、種々の斜面において地すべり地形や深層崩壊地形、表層・浅層崩壊地形を抽出する際に、極めて有効な手がかりを提供することができる。本願発明を利用すれば、地すべり地形等の抽出に役立つとともに、地すべり等のメカニズムの把握につながることを考えれば、自然災害を未然に防ぎ、被害を軽減させることが可能となり、産業上利用できるうえに社会的にも貢献が期待できる発明といえる。
【符号の説明】
【００４４】
１００  本願発明の地形分類システム
１１０  （地形分類システムの）条件設定手段
１１１  （条件設定手段の）地形モデル記憶手段
１１２  （条件設定手段の）地形モデル読出し手段
１１３  （条件設定手段の）サンプル領域抽出手段
１１４  （条件設定手段の）注目線設定手段
１１５  （条件設定手段の）地形量算出手段
１１６  （条件設定手段の）散布図作成手段
１１７  （条件設定手段の）地形分類条件設定手段
１１８  （条件設定手段の）地形分類条件記憶手段
１２０  （地形分類システムの）特徴領域抽出手段
１２１  （特徴領域抽出手段の）地形量算出手段
１２２  （特徴領域抽出手段の）地形分類条件読出し手段
１２３  （特徴領域抽出手段の）特徴地形分類手段
１２４  （特徴領域抽出手段の）特徴領域生成手段
１２５  （特徴領域抽出手段の）地形変化線設定手段
１３０  （地形分類システムの）出力手段
ＡＲ    所定領域
ＢＬ    単位区画
ＳＡ    サンプル領域
ＴＭ    地形データ
【図１】
0006746851.tif 000002
【図２】
0006746851.tif 000003
【図３】
0006746851.tif 000004
【図４】
0006746851.tif 000005
【図５】
0006746851.tif 000006
【図６】
0006746851.tif 000007
【図７】
0006746851.tif 000008
【図８】
0006746851.tif 000009
【図９】
0006746851.tif 000010