(57)【要約】【課題】簡単な構成で高精度に管網シミュレーションを行うことが可能な管網シミュレーション装置を提供する。【解決手段】地図の画像を取り込む画像入力手段１１と、取り込んだ画像から輪郭を抽出するエッジ検出手段１２と、取り込んだ画像を予め設定する大きさの領域で区切る画像分割手段１３と、前記領域ごとのエッジの量を数えるエッジ量カウント手段１４と、前記エッジ量カウント手段１４によって求めたエッジ量に基づいて人口密度の推定値を算出する人口密度推定手段１５と、前記人口密度推定手段１５によって求めた人口密度の推定値を用いて上下水道の設備の状態変化を求める管網シミュレーション手段１６とを備えるようにした。【選択図】図１
2010020707.tif 000002
【特許請求の範囲】
【請求項１】
地図の画像を読み込む画像入力手段と、
画像処理により前記画像入力手段によって取り込んだ入力画像上の線の輪郭を抽出するエッジ検出手段と、
取り込んだ画像を予め設定する大きさの領域に分割する画像分割手段と、
前記画像分割手段によって区切られた領域ごとのエッジの量を求めるエッジ量カウント手段と、
前記エッジ量カウント手段によって求めたエッジの量に基づいて人口密度の推定値を算出する人口密度推定手段と
前記人口密度推定手段によって求めた人口密度の推定値を用いて上下水道の設備の状態変化を求める管網シミュレーション手段と
を備えることを特徴とする管網シミュレーション装置。
【請求項２】
地図上の二地点以上の位置における人口密度の真値に基づいて人口密度の推定に用いる係数を演算する係数演算手段
を備えることを特徴とする請求項１記載の管網シミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、管網シミュレーション装置に関し、特に画像処理により検出した地図のエッジ情報を利用してなる管網シミュレーション装置に関する。
【背景技術】
【０００２】
従来、地理情報システムを利用して配電線網や、ガス・水道管網を管理するシステムがある。地図情報システムは、地図に関する属性情報（構造物、人口密度、土地利用などの情報）を持ち、これら属性情報を加工して都市計画や土地管理・マーケティング等を支援するための情報を得ることができる。例えば、下水道分野においては、人口密度や土地利用状態を考慮して下水道の管網シミュレーションを行うことにより、近年の都市化の影響によって懸念される降雨時の急激な雨水増量による浸水等の状態を予測したり把握・検討したりすることができる。
【０００３】
また、特許文献１、非特許文献１、非特許文献２には、排水面積に対する不浸透域の割合である不浸透面積率を求める方法として以下の方法が開示されている。
１）サンプリング法：航空写真あるいは土地利用図に適当なメッシュを切り、乱数でいくつかのメッシュを選び出し、選んだメッシュ内の表面工種が浸透面か不浸透面かを読み取り、その割合から全体の不浸透面積率の推進値とする方法である。
２）画像処理法：航空写真あるいは土地利用図に対し、色抽出処理によって緑地や土部分を抽出し、不浸透面積率を求める方法である。
【０００４】
【特許文献１】特開２００３−３２２５２３号公報
【非特許文献１】渡辺、川本著、「合流式下水道の改善対策に関する研究調査」、下水道協会誌、1987年9月、第24巻、第280号、p.35-53
【非特許文献２】「下水道雨水流出解析に使用する等到達時間域不浸透面積率のカラー航空写真の色抽出処理に基づく計算」、平成14年度電気学会産業応用部門全国大会講演論文集［Ｉ］、2002年、p.541-546
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上述したような管網シミュレーションに人口密度の推定値を用いる場合、この人口密度の推定値を求める手段として、上述したサンプリング法、画像処理法を利用することが考えられる。しかしながら、上述した従来の手法を利用して人口密度の推定値を求める場合には次のような問題があった。
【０００６】
１）サンプリング法
（イ）作業量が膨大で調査に時間がかかる。
（ロ）作業者によって調査結果にバラツキが生じる場合がある。
【０００７】
２）航空写真に対する画像処理法
（イ）夕日や朝日などの太陽光やカメラ感度などの撮影条件により写真の色状態が一定ではないため、同じ地区の航空写真であっても撮影日時の異なる航空写真に対してそのまま色抽出すると結果が異なる可能性がある。
（ロ）航空写真を貼り合せて広域の画像にする場合、撮影時刻、天候、太陽光、カメラ感度などの撮影条件が異なることにより写真の境界線で色と明暗が不連続になり境界線にノイズができる可能性がある。
【０００８】
このような問題により、管網シミュレーションを煩雑にし、またシミュレーションの精度を低下させるおそれがあった。
【０００９】
このようなことから本発明は、簡単な構成で高精度に管網シミュレーションを行うことが可能な管網シミュレーション装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記の課題を解決するための第１の発明に係る管網シミュレーション装置は、地図の画像を読み込む画像入力手段と、画像処理により前記画像入力手段によって取り込んだ入力画像上の線の輪郭を抽出するエッジ検出手段と、取り込んだ画像を予め設定する大きさの領域に分割する画像分割手段と、前記画像分割手段によって区切られた領域ごとのエッジの量を求めるエッジ量カウント手段と、前記エッジ量カウント手段によって求めたエッジの量に基づいて人口密度の推定値を算出する人口密度推定手段と前記人口密度推定手段によって求めた人口密度の推定値を用いて上下水道の設備の状態変化を求める管網シミュレーション手段とを備えることを特徴とする。
【００１１】
また、第２の発明に係る管網シミュレーション装置は、第１の発明に係る管網シミュレーション装置において、地図上の二地点以上の位置における人口密度の真値に基づいて人口密度の推定に用いる係数を演算する係数演算手段を備えることを特徴とする請求項１記載の画像処理による人口密度推定装置。
【発明の効果】
【００１２】
上述した第１の発明に係る管網シミュレーション装置によれば、通常の人口密度は広域をひとまとめにした人口密度で表されるのに対し、地図上の局所的で多地点の人口密度がわかることにより、等高線表示のような局所的な人口密度分布が容易に作成でき、より高精度に管網シミュレーションを行うことが可能となる。
【００１３】
また、航空写真を用いる場合、撮影時刻、天候、太陽光、カメラ感度などの撮影条件により写真にバラツキが生じるおそれがあり、画像処理の結果にバラツキが生じる可能性があるのに対し、本発明では描写にバラツキのない地図を用いるため画像処理による人口密度の推定結果にバラツキが生じるおそれがないという利点がある。
【００１４】
さらに、航空写真を用いる場合、複数の航空写真をタイル状に繋ぎ合わせて全体の写真を構成する必要があるため、繋ぎ合わせる際に境界部分が不連続になりこれがノイズになるおそれがあるのに対し、本発明では地図を利用するためこのような境界線によるノイズが生じるおそれがなく人口密度の推定をより確実に行うことができるという利点がある。
【００１５】
さらに加えて、地図から得たエッジ量を管網シミュレータで解析することにより、次の利点がある。
【００１６】
上水道分野では、人口密度から配水状態や配水バランスのシミュレーションができるため、配水地やポンプやバルブ等設備の状態変化の把握ができる。さらに設備の運用方法の検討や課題の対策ができる利点がある。
【００１７】
下水道分野では、土地利用状態や人口密度から下水道の管網シミュレーションができるため、近年の都市化による雨水の浸水状態の把握や対策ができる利点がある。
【００１８】
また、第２の発明に係る管網シミュレーション装置によれば、二地点以上の人口密度の真値をもとに算出した係数を用いて人口密度の推定値を求めるので、多地点の人口密度を推定するための係数としてより信頼性が高い係数を得ることができ、管網視乳レーションをより高精度に行うことができる。
【発明を実施するための最良の形態】
【００１９】
本発明の実施形態の詳細を以下の実施例において説明する。
【実施例１】
【００２０】
図１乃至図５を用いて本発明の第１の実施例を説明する。図１は本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置の構成を表すブロック図、図２は本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置における処理の流れを示すフローチャート、図３は画像を複数の領域に分割する例を示す説明図、図４は画像中の輪郭を抽出した例を示す説明図、図５は本実施例により得られる人口密度を等高線表示した例を示す説明図である。
【００２１】
図１に示すように、本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置は、画像入力手段１１、エッジ検出手段１２、画像分割手段１３、エッジ量カウント手段１４、人口密度推定手段１５、及び、管網シミュレート手段１６を備え、一般的な地図を利用してこの地図上に表示された家屋やビル等といった建造物を画像処理によりエッジとして抽出し、得られたエッジの量に基づいて推定した人口密度を利用して管網シミュレーションを行うものである。なお、人口密度の推定は、図１中破線で囲んだ部分で行われる。
【００２２】
画像入力手段１１は国土地理院等が発行する属性情報のない一般的な地図を例えばスキャナー装置等により入力画像として本装置のメモリへ読み込む手段である。以下、例として、図３に模式的に示した地図から人口密度を推定する場合について説明する。
【００２３】
エッジ検出手段１２は、画像入力手段１１によって読み込んだ例えば図３に示すような入力画像に対し、ｓｏｂｅｌフィルタや微分フィルタにより入力画像中に表示されている線、即ち、家屋やビル等といった建造物などの輪郭を抽出し、これをエッジとして検出する手段である。該エッジ検出手段１２では、例えば図４に示すように、検出したエッジ部分を白色、それ以外の部分を黒色として出力する。
【００２４】
画像分割手段１３はエッジ検出手段１２によってエッジを検出された入力画像を予め設定された単位面積に対応するメッシュで縦横に区切り、複数に分割する手段である。例えば、入力画像を横方向に３１分割、縦方向に２３分割した場合は、図４に示すような領域に区切られる。
【００２５】
エッジ量カウント手段１４は画像分割手段１３によって区切られた領域、例えば図４に示すそれぞれの位置（ｉ，ｊ）（但し、１≦ｉ≦３１、１≦ｊ≦２３）の領域ごとにエッジの量をカウントし、各領域のエッジ量ｅ（ｉ，ｊ）を求める手段である。該エッジ量カウント手段１４では、例えば、図４に示すそれぞれの領域ごとに輝度を計測し、この輝度を各々の位置（ｉ，ｊ）におけるエッジ量として出力する。
【００２６】
人口密度推定手段１５は下記（１）式に示すようにエッジ量カウント手段１４によってそれぞれの領域ごとに計測されたエッジ量ｅ（ｉ，ｊ）に係数ｋを乗じて、画像分割手段１３によって区切られた領域ごとの人口密度の推定値ｄ（ｉ，ｊ）を算出する手段である。なお、係数ｋは予め実験的に求めた値とする。
【００２７】
ｄ（ｉ，ｊ）＝ｋ・ｅ（ｉ，ｊ）      ・・・（１）
【００２８】
管網シミュレート手段１６は所定の地点のエッジ量に基づき上下水道の設備の状態変化を求める手段である。
【００２９】
次に、図２を用いて本実施例における管網シミュレートの流れの一例を説明する。図２に示すように、本実施例においてはまず画像入力手段１１において人口密度推定の対象とする領域の地図を入力画像として本装置のメモリに読み込み（ステップＳ１）、エッジ検出手段１２によって画像中に表示された家屋やビル等をエッジとして検出する（ステップＳ２）。
【００３０】
続いて、画像分割手段１３によって画像を予め設定する単位面積に区切り（ステップＳ３）、エッジ量カウント手段１４でそれぞれの領域ごとのエッジ量ｅ（ｉ，ｊ）をカウントし（ステップＳ４）、人口密度推定手段１５において上述した（１）式の演算を行ってそれぞれの領域ごとの人口密度の推定値ｄ（ｉ，ｊ）を算出する（ステップＳ５）。そして、それぞれの領域ごとの人口密度の推定値ｄ（ｉ，ｊ）を算出したら、この推定値ｄ（ｉ，ｊ）を利用して、既知の手法により管網シミュレーションを行う（ステップＳ６）。
【００３１】
その後、画像の入力が終了したか否かの判定を行い（ステップＳ７）、画像の入力が終了していない場合（Ｎｏ）はステップＳ１の処理に戻りステップＳ１〜ステップＳ６の処理を繰り返す。一方、ステップＳ７における判定の結果、画像の入力が終了したと判定された場合（Ｙｅｓ）は人口密度推定処理を終了する。
【００３２】
以上に説明した本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置によれば、通常の人口密度が広域をひとまとめにして表されるのに対し、地図上の局所的で多地点の人口密度がわかる利点がある。これにより、例えば図５に示すような局所的な人口密度分布が容易に作成できる利点がある。なお、図５は人口密度が多いと推定される部分が濃色となるように表した例である。
【００３３】
また、航空写真を用いる場合、撮影時刻、天候、太陽光、カメラ感度などの撮影条件により写真にバラツキがあり、そのまま画像処理すると結果にバラツキが生じることが考えられる。これに対し本実施例に係る管網シミュレーション装置によれば、地図にバラツキがないため画像処理の結果にバラツキが生じることが無い。
【００３４】
さらに、航空写真を用いる場合、航空写真をタイル状に繋ぎ合わせて全体の写真を構成する必要があるため、繋ぎ合わせる際に境界線の部分が不連続になりこれがノイズになるおそれがあるのに対し、本実施例では地図を利用することにより境界線によるノイズが生じるおそれがない。
【００３５】
また、地図から得たエッジ量を管網シミュレータで解析することにより、次の利点がある。
【００３６】
上水道分野では、人口密度から配水状態や配水バランスのシミュレーションができるため、配水地やポンプやバルブ等設備の状態変化の把握ができる。さらに設備の運用方法の検討や課題の対策ができる利点がある。
【００３７】
下水道分野では、土地利用状態や人口密度から下水道の管網シミュレーションができるため、近年の都市化による雨水の浸水状態の把握や対策ができる利点がある。
【実施例２】
【００３８】
図６乃至図８を用いて本発明の第２の実施例を説明する。図６は本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置に利用する人口密度の推定値を求める装置の構成を表すブロック図、図７はエッジ量と人口密度真値との関係の一例を表すグラフ、図８はエッジ量と人口密度真値との関係の一例を表す他のグラフである。
【００３９】
図６に示すように、本実施例において人口密度を推定する装置は、図１に示し上述した実施例１の地図のエッジ情報を利用してなる管網シミュレーション装置に係数演算手段１６を追加した例である。その他の構成は図１に示し上述した実施例１の構成と概ね同様であり、以下、同様の作用効果を奏する部分には同一の符号を付し、重複する説明は省略する。
【００４０】
本実施例において、係数演算手段１６は地図上の二箇所以上の地点の人口密度の真値に基づきエッジ量カウント手段１４において得たエッジ量ｅを人口密度の推定値ｄに変換するための係数ｋを算出する手段である。例えば、図７に示すように、ある地点Ａにおけるエッジ量がｅa、人口密度の真値がｄaであり、他の地点Ｂにおけるエッジ量がｅb、人口密度の真値がｄbである場合、係数ｋは以下の（２）式に示す値として算出される。
【００４１】
【数１】
2010020707.tif 000003
【００４２】
そして、人口密度推定手段１５によって（２）式に示す係数ｋを用い、以下の（３）式に示す演算を行うことにより、人口密度の推定値ｄを求めることができる。
【００４３】
ｄ−ｄa＝ｋ・（ｅ−ｅa）    …（３）
【００４４】
以上に説明したように、本実施例に係る地図のエッジ情報を利用してなる管網シミュレーション装置によれば、二地点以上の人口密度の真値をもとに算出した係数を用いて人口密度の推定値を求めるので、多地点の人口密度を推定するための係数としてより信頼性が高い係数を得ることができるため、人口密度の推定を高効率に行うことができる。
【００４５】
なお、本実施例では係数ｋを二地点の人口密度の真値を用いて求める例を示したが、例えば図８に示すように、ある地点Ａ，Ｂ，Ｃにおけるそれぞれのエッジ量ｅa，ｅb，ｅcに対して、人口密度の真値がｄa，ｄb，ｄcである場合は、複数の点ごとに（３）式を立ててエッジ量カウント手段１４において得たエッジ量ｅを人口密度の推定値ｄに変換するようにしてもよい。
【００４６】
即ち、エッジ量カウント手段１４において得られたエッジ量ｅと地点Ａ，Ｂ，Ｃにおけるそれぞれのエッジ量ｅa，ｅb，ｅcとを比較し、例えばエッジ量カウント手段１４において得られたエッジ量ｅがｅb＜ｅ＜ｅcを満たす場合は、このエッジ量ｅを人口密度の推定値ｄに変換するための係数ｋとして、地点Ｂ，Ｃにおけるそれぞれのエッジ量ｅb，ｅc及び人口密度の真値ｄb，ｄcから得られる係数ｋ＝ｋ’を用いて上記（３）式からこの地点の人口密度の推定値ｄを求めるようにしても良い。
【産業上の利用可能性】
【００４７】
本発明は、地図のエッジ情報を利用してなる管網シミュレーション装置に適用して好適なものである。
【図面の簡単な説明】
【００４８】
【図１】本発明の実施例１に係る地図のエッジ情報を利用してなる管網シミュレーション装置の構成を表すブロック図である。
【図２】本発明の実施例１における管網シミュレーションの流れを示すフローチャートである。
【図３】本発明の実施例１において画像を分割する例を示す説明図である。
【図４】本発明の実施例１において画像上の輪郭を抽出する例を示す説明図である。
【図５】人口密度推定の結果の一例を示す説明図である。
【図６】本発明の実施例２に係る地図のエッジ情報を利用してなる管網シミュレーション装置の構成を表すブロック図である。
【図７】地図から抽出したエッジ量と人口密度との関係の一例を示すグラフである。
【図８】地図から抽出したエッジ量と人口密度との関係の一例を示す他のグラフである。
【符号の説明】
【００４９】
１  管網シミュレーション装置
１１  画像処理手段
１２  エッジ検出手段
１３  メッシュ分割手段
１４  エッジ量カウント手段
１５  人口密度推定手段
１６  管網シミュレート手段
【図１】
2010020707.tif 000004
【図２】
2010020707.tif 000005
【図３】
2010020707.tif 000006
【図６】
2010020707.tif 000007
【図７】
2010020707.tif 000008
【図８】
2010020707.tif 000009
【図４】
2010020707.tif 000010
【図５】
2010020707.tif 000011