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