(57)【特許請求の範囲】
【請求項１】
遠方を飛行する鳥類の飛行経路を計測する飛行経路計測システムであって，
短焦点光学系を備え、あらかじめ定めた目的とする広い計測領域を撮影する複数台の撮像手段と，
前記複数台の撮像手段のうち、鳥類を撮影した撮像手段により撮影された画像の画像データを時系列的に記録する記録手段と，
前記記録された画像から鳥類の軌跡を導出する処理手段と，
得られた鳥類の軌跡を飛行経路として表示する表示手段を備え，
前記処理手段は，
特定の画像データを取り込み，取り込んだ画像データのうち計測の対象とする鳥類の像の指定を受け付け，前記像を参照パターンとして粒子追跡法を用い，以後の時刻の画像から像を追跡して鳥類の軌跡を導出することを特徴とする飛行経路計測システム。
【請求項２】
前記処理手段は，撮影エリアの地図情報を参照し，前記表示手段において，前記飛行経路の地図上への表示を可能とする処理を行うことを特徴とする請求項１に記載の飛行経路計測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，鳥類の飛行経路を計測する飛行経路計測システムに関する。
【背景技術】
【０００２】
従来，鳥類などの遠方の飛行物体が飛行する経路を計測する場合には，専門の調査員が双眼鏡を使用するなどして計測を行っている。例えば，発電所の周囲環境への影響評価を行うときの鳥類の調査方法が非特許文献１に記載されている。ポイントセンサス法による調査であり，見通しのきく場所に定点を設定し，出現する種類等を直接観察により記録する。視野の範囲内の識別が可能な距離までを対象とし，一定時間観察を行うものである。この場合，鳥類の日々の飛行経路を計測する必要が生じることがある。このような場合，専門の調査員が計測対象エリアの所定の地点に観測小屋を建て，一定期間，飛行経路の計測を行う。
【０００３】
また，移動物体の三次元軌跡を計測する技術として，球技スポーツにおけるボールのように高速で移動する物体の移動軌跡を求める技術が特許文献１に記載されている。これは，２以上のカメラで撮像した画像から三次元の軌跡を求めるものである。
【非特許文献１】発電所に係る環境影響評価の手引き（平成１９年１月改訂 経済産業省 原子力安全・保安院）
【特許文献１】特開２００７−１１５２３６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
従来の鳥類の飛行経路計測方法では，天候や計測対象の鳥の挙動によっては，計測の労力が大きくなるという問題がある。また，計測期間が長期にわたると，人件費が高くなるという問題がある。
また，特許文献１の技術は，例えば，背景差分法を用いており，この場合は，背景が時間的に変化するときに軌跡を求めることが難しい。
【０００５】
一方，近年，ある時間間隔で画像中の各トレーサ粒子の移動を自動的に追跡する粒子追跡法が，主に流れ場を計測する方法の一つとして用いられている。
【０００６】
本発明は，主に流れ場を計測する方法として用いられている粒子追跡法を，遠方を飛行する鳥類の追跡に適用し，背景が静止していない屋外における鳥類の飛行経路を計測することを可能とし，計測にかかる労力や人件費を抑制可能とする飛行経路計測システムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
請求項１記載の発明では，遠方を飛行する鳥類の飛行経路を計測する飛行経路計測システムであって，短焦点光学系を備え、あらかじめ定めた目的とする広い計測領域を撮影する複数台の撮像手段と，前記複数台の撮像手段のうち、鳥類を撮影した撮像手段により撮影された画像の画像データを時系列的に記録する記録手段と，前記記録された画像から鳥類の軌跡を導出する処理手段と，得られた鳥類の軌跡を飛行経路として表示する表示手段を備え，前記処理手段は，特定の画像データを取り込み，取り込んだ画像データのうち計測の対象とする鳥類の像の指定を受け付け，前記像を参照パターンとして粒子追跡法を用い，以後の時刻の画像から像を追跡して鳥類の軌跡を導出することを特徴とする飛行経路計測システムを提供する。
【０００８】
請求項２記載の発明では，前記処理手段は，撮影エリアの地図情報を参照し，前記表示手段において，前記飛行経路の地図上への表示を可能とする処理を行うことを特徴とする請求項１に記載の飛行経路計測システムを提供する。
【０００９】
前記処理手段では，粒子追跡法（PTV：Particle Tracking Velocimetry）と呼ばれる技術を適用する。すなわち，遠方を飛行する鳥類を粒子とみなし，その移動を自動的に追跡する。
【発明の効果】
【００１０】
本発明の飛行経路計測システムによれば，遠方を飛行する鳥類は複数台の撮像手段によって撮影され，その撮像データは処理手段においてPTVを適用して処理されることにより，鳥類の飛行経路が計測される。従って，天候や鳥の挙動の影響を受けることなく鳥類の飛行経路を計測し，計測にかかる労力や人件費を抑制することができる。
【００１１】
また，PTVの適用により，背景が静止していない屋外において飛行する鳥類の軌跡を求めることが可能となる。
【発明を実施するための最良の形態】
【００１２】
以下に添付図面を参照しながら，本発明の好適な実施形態について説明する。
【００１３】
図１には，本発明の第一の実施形態にかかる飛行経路計測システムの全体構成図を示す。
【００１４】
本発明にかかる飛行経路計測システムは，撮像手段１０，記録手段２０，処理手段３０，及び表示手段４０から構成される。
【００１５】
撮像手段１０は記録手段２０に接続される。撮像手段１０は，あらかじめ定めた計測領域を撮影し，撮影により得られた画像データを記録手段２０に伝送する役割を担う。撮像手段１０と記録手段２０の間の伝送は有線でも無線でも良い。同一の鳥類を２以上の撮像手段１０で計測することにより，後述する通り，鳥類の三次元計測が可能である。
【００１６】
より広いエリアを撮影する為には，撮像手段１０の数を増やせば良いが，短焦点光学系を備えた撮像手段１０を用いれば，広角撮影が可能であり，より少ない撮像手段１０の数で目的とする広いエリアを撮像できる。短焦点光学系としては，焦点距離が１６ｍｍ以下のレンズを用いると良い。
【００１７】
撮像手段１０には，例えばCCD（Charge Coupled Device）撮像素子を備えたカメラ（CCDカメラ）が用いられる。CMOS（Complementary Metal Oxide Semiconductor）撮像素子を備えたカメラを用いても良い。CCDカメラは例えば，毎秒３０コマの頻度で調査エリアを撮影する。
【００１８】
記録手段２０は，撮像手段１０から伝送された画像データを時系列的に記録する役割を担う。記録手段２０は装置として独立していても良いし，撮像手段１０や後述する処理手段３０に，その一部として内蔵されていても良い。記録手段２０には，処理手段３０にて処理すべき画像データを記録可能な容量の，例えばメモリ素子を用いる。
【００１９】
記録手段２０は，処理手段３０に接続される。処理手段３０は，記録手段２０に記録された画像データを取り込む。取り込んだ画像データに対してPTVを適用する。
【００２０】
具体的には，ある画像データにおいて，計測の対象とする鳥類の像を，例えば計測者が指定する。指定は，例えば，表示手段４０に現れた十字カーソルにより選択することにより行う。選択された鳥類の像の情報は，例えば，画素が有する色情報を利用することにより，画像データから抽出される。
【００２１】
処理手段３０は指定指令を受け付け，指定された像を画像データから抽出し，参照パターンとして認識する。そして，以後の時系列の画像データから順次参照パターンと符合する像を自動追跡し，自動追跡により得られた像の位置を画像座標上に得る。
【００２２】
パターンの符合には，例えば，相互相関法が用いられる。自動追跡する鳥類が数体程度であれば，２時刻間の鳥類位置の最も近いもの同士を対応付ける方法でも追跡できる。相互相関法の代わりに，粒子マスク相関法，二値化相関法，三時刻パターンマッチング法なども適用可能である。
【００２３】
鳥類の三次元計測には，例えばステレオ法が用いられる。ステレオ法では，例えば２台のCCDカメラで撮影された画像中で同一時刻における同一の鳥類を同定できれば，三角測量の原理で鳥類の三次元位置の計測が可能である。
【００２４】
PTVの手法については，以下の非特許文献２に詳しい解説がある。
【非特許文献２】流れの可視化ハンドブック（浅沼強編，朝倉書店）
【００２５】
処理手段３０は，撮影エリアの地図情報を有し，これを参照することにより，鳥類の飛行経路を地図上へ表示する可能とする処理も行う。
【００２６】
処理手段３０は，表示手段４０に接続される。表示手段４０は，処理手段３０にて得られた鳥類の像の位置を飛行経路として表示する役割を担う。表示手段４０にて，飛行経路を地図上に表示することも可能である。
【００２７】
表示手段４０には，例えば，液晶表示装置を用いる。
【００２８】
図２は，処理手段３０におけるパターン符合方法の例の一つである相互相関法の説明図である。図２のうち，黒い丸は参照パターンとして抽出された計測対象の鳥類の像であり，斜線の丸は参照パターンが抽出された時刻から撮像間隔Δｔ経過後における鳥類の像である。参照パターンの位置情報を含む所定の数の画素群を画素群０とする。画素群０の一部を含むように，Δｔ経過後の画素群１，画素群２，・・・，画素群Ｎを選ぶ。図２の場合は，Ｎ＝４である。画素群０を基準として，画素群１，画素群２，・・・，画素群Ｎとの相関を調べる。具体的には，画素群０に含まれる参照パターンと符合するパターンを含む画素群を探索する。図２の場合，画素群０との相関が最も大きいと判断されるのは，画素群３であるから，参照パターンが抽出された時刻からΔｔ経過後の鳥類の像は画素群３にあると判断される。参照パターンを含む画素群と，探索対象の画素群との時間間隔は，Δｔの倍数であっても良い。
【００２９】
図３は，本計測システムを用いた計測処理の一例を示すフロー図である。１台の撮像手段１０により，鳥類の二次元の飛行経路表示を得る例である。
【００３０】
まず，カウンタｎに計測開始時刻に対応する値である１が代入される（ステップＳ３０１）。撮像手段１０は，あらかじめ定めた計測領域を，カウンタｎに対応する時刻tｎにて撮影する。撮影により取得された画像データは，記録手段２０に伝送され，時系列で記録される（ステップＳ３０２）。時刻tｎが計測終了時を示すｔｅｎｄ以上であれば，ステップＳ３０４に進み，ｔｅｎｄより小さければ，ステップＳ３０５に進む。ステップＳ３０５では，カウンタｎに１が加えられ，tｎには撮像間隔Δｔが加えられて撮影が継続される（ステップＳ３０３）。
【００３１】
ステップＳ３０４で，カウンタｎに，軌跡表示の開始時刻を示すｔｋに対応するｋが代入されると，時刻tｎで取得された画像が，表示手段４０に呼び出される（ステップＳ３０６）。呼び出された画像上において，例えば，計測者が，表示手段４０に現れた十字カーソル等により，追跡する鳥類を指定する。これを参照パターンとする（ステップＳ３０７）。
【００３２】
追跡する鳥類が指定されると，カウンタｎに１が加えられ，tｎには撮像間隔Δｔが加えられて（ステップＳ３０８），時刻tｎで取得された画像が呼び出される（ステップＳ３０９）。処理手段３０は，時刻tｎの画像で，時刻ｔｎ−１で特定した像を参照パターンとして追跡し，時刻tｎの画像における鳥類の像を特定する（ステップＳ３１０）。処理手段３０は，時刻tｎで特定された像の位置（ｘ，ｙ）と時刻tｎを対応づけて格納する（ステップＳ３１１）。
【００３３】
ステップＳ３１２において，ステップＳ３０３と同様，時刻tｎが，計測終了時を示すｔｅｎｄと比較され，ｔｅｎｄより小さければ，像の追跡が継続される。時刻ｔｎがｔｅｎｄ以上であれば，ステップＳ３１３に進み，カウンタｎに，上記の軌跡表示の開始時刻ｔｋに対応するｋが代入される。ここで，カウンタｎに，ｋより大きい値であるｍが代入され，時刻ｔｋ以後の任意の時刻ｔｍを軌跡表示の開始点とすることも可能である。
【００３４】
ステップＳ３１４において，時刻tｎで取得された画像が呼び出される。続いて，カウンタｎに１が加えられ，tｎには撮像間隔Δｔが加えられる（ステップＳ３１５）。時刻tｎの画像上に，ステップＳ３１１で格納された時刻tｋ，tｋ＋１，・・・，tｎにおける像の位置が重ねて表示されることにより（ステップＳ３１６），計測対象の鳥類の軌跡が表示される。時刻tｎがｔｅｎｄより小さい場合は，鳥類の軌跡の表示が継続され，ｔｅｎｄ以上の場合は，フローは終了する（ステップＳ３１７）。
【００３５】
図４は，本計測システムを用いた計測処理の次の例を示すフロー図である。１台の撮像手段１０により，鳥類の二次元の飛行経路表示を地図上に得る例である。
【００３６】
本処理フローにおけるステップＳ４０１からＳ４１１は，図３に示す処理フローにおけるステップＳ３０１からＳ３１１までと同じ為，説明を省略する。
【００３７】
ステップＳ４１２では，あらかじめ処理手段２０に格納された地図情報が参照され，画像座標上に得られた時刻ｔｎの鳥類の像の位置から，地図上における時刻tｎの像の位置（Ｘ，Ｙ）が算出され，格納される。
【００３８】
時刻tｎが，時刻終了時を示すｔｅｎｄと比較され，ｔｅｎｄより小さければ，像の追跡及び追跡された像の地図上の位置の算出が継続され，ｔｅｎｄ以上であれば，ステップＳ４１４に進む（ステップＳ４１３）。
【００３９】
ステップＳ４１４では，カウンタｎに，ステップＳ４０４で代入された，軌跡表示の開始時刻を示すｔｋに対応するｋが代入される。時刻tｎで取得された画像が呼び出され（ステップＳ４１５），軌跡を表示させたい鳥類の像を指定すると（ステップＳ４１６），表示手段４０に表示された地図上に，ステップＳ４１２で格納された時刻tｎの像の位置が表示される（ステップＳ４１７）。
【００４０】
続いて，カウンタｎに１が加えられ，tｎには撮像間隔Δｔが加えられて（ステップＳ４１８），地図上に，時刻tｎの鳥類の像の位置が重ねて表示される（ステップＳ４１９）。時刻tｎがtｅｎｄと比較され，tｅｎｄより小さい場合は，鳥類の軌跡の地図上への表示が継続され，tｅｎｄ以上の場合は，フローは終了する（ステップＳ４２０）。
【００４１】
図５は，本計測システムを用いた計測処理の次の例を示すフロー図である。１台の撮像手段１０により，鳥類の二次元の飛行経路表示を得る他の例である。
【００４２】
本処理フローにおけるステップＳ５０１からＳ５１０は，図３に示す処理フローにおけるステップＳ３０１からＳ３１０までと同じ為，説明を省略する。
【００４３】
ステップＳ５１１では，時刻tｎの画像上に，ステップＳ５１０で特定された時刻tｋ，tｋ＋１，・・・，tｎの鳥類の像の位置が重ねて表示されることにより，計測対象の鳥類の軌跡が表示される。なお，ステップＳ５１０で特定された時刻tｎの像の位置を格納しても良い。
【００４４】
時刻tｎが，時刻終了時を示すｔｅｎｄと比較され，ｔｅｎｄより小さければ，鳥類の軌跡の表示が継続され，ｔｅｎｄ以上であれば，フローは終了する（ステップＳ５１２）。
【００４５】
図６及び図７は，本計測システムを用いた計測処理の次の例を示すフロー図である。複数の撮像手段１０により，鳥類の三次元の飛行経路表示を得る例である。
【００４６】
本フロー図では，撮像手段１０の数が２の場合を示す。撮像手段１０の数がＮの場合も，Ｎ台のうち，同一の鳥類を撮影した２台を選び本フローを適用することにより，鳥類の三次元の飛行経路表示が可能である。
【００４７】
まず，カウンタｉ１，ｉ２に，それぞれ撮像手段１０を示す番号である１，２が代入される（ステップＳ６０１１，Ｓ６０１２）。
【００４８】
次に，ステップＳ６５０１，Ｓ６５０２に進み，それぞれの撮像手段１０により，あらかじめ定めた計測領域の撮影と画像データの記録が行われる。ここで，ステップＳ６３０１では，番号１の撮像手段１０（カメラ）で撮影される。ステップＳ６３０２では，番号２の撮像手段１０（カメラ）で撮影される。ステップＳ６５０１，Ｓ６５０２は，図３に示す処理フローにおけるステップＳ３０１からＳ３０５まで（Ｓ３０４を除く）と同じ為，説明を省略する。
【００４９】
ステップＳ６０５１，Ｓ６０５２で，カウンタｉ１，ｉ２に，撮像手段１０の番号を示す１，２がそれぞれ代入される。
【００５０】
次に，ステップＳ６６０１，Ｓ６６０２に進む。ステップＳ６６０１では，番号１の撮像手段１０により撮影，記録された画像データについて，処理が行われる。ステップＳ６６０２では，番号２の撮像手段１０による画像データについてステップＳ６６０１と同様の処理が行われる。
【００５１】
ステップＳ６６０１では，カウンタｎに，軌跡表示の開始時刻を示すｔｋに対応するｋが代入される（ステップＳ６０７１）。ステップＳ６０８１からＳ６１４１までは，図３に示す処理フローにおけるステップＳ３０６からＳ３１２までと同じ為，説明を省略する。ただし，ステップＳ６１３１では，撮像手段１０の番号ｉ１も対応づけて記録される。
【００５２】
ステップＳ６１５に進み，ステップＳ６０７１と同様に，カウンタｎに，軌跡表示の開始時刻ｔｋに対応するｋが代入される。ステップＳ６１３１，Ｓ６１３２で記録された時刻tｎの２つの像の位置（ｔｎ，ｘ１，ｙ１，ｉ１（＝１）），（ｔｎ，ｘ２，ｙ２，ｉ２（＝２））から，例えば，ステレオ法により像の三次元位置を取得し，時刻tｎの情報を対応づけて（ｔｎ，ｘ，ｙ，ｚ）として格納する（ステップＳ６１６）。あらかじめ格納された地図情報が参照され，地図上における時刻tｎの鳥類の像の三次元位置（Ｘ，Ｙ，Ｚ）が算出され，（ｔｎ，Ｘ，Ｙ，Ｚ）として格納される（ステップＳ６１７）。
【００５３】
時刻tｎがｔｅｎｄと比較され，ｔｅｎｄ以上であればステップＳ６１９に進み，ｔｅｎｄより小さければ，ステップＳ６２０に進む（ステップＳ６１８）。ステップＳ６２０では，カウンタｎに１が加えられ，tｎには撮像間隔Δｔが加えられて，地図上における像の三次元位置の算出及び格納が継続される。
【００５４】
ステップＳ６１９からＳ６２６（Ｓ６２０を除く）は，処理する像の位置情報が，図４に示す処理フローにおけるステップＳ４１４からＳ４２０が処理する二次元のものから，三次元のものに替わっただけであり，実質的には同じ為，説明を省略する。以上の処理により，鳥類の軌跡の地図上への三次元表示が可能となる。
【００５５】
本処理フローは，図４の処理フローで得られる鳥類の二次元の飛行経路表示を三次元で得る為のフローであるが，同様に，図３及び図５に示す処理フローで得られる二次元の飛行経路表示を三次元で得ることも可能である。
【００５６】
また，鳥類の像に関する情報として，位置と時刻が得られていることから，これから速度を求め，求めた速度を表示手段４０に表示することも可能である。
【００５７】
図８は，地図上への飛行経路の表示例を示す。等高線によって示された地形図に飛行経路が表された例である。図８のうち，細い実線は等高線，黒い四角は風力発電システム，太い実線は計測された鳥類の飛行経路を示す。図８は平面図による二次元表示の例であるが，立面図，三次元的な立体図による表示も可能である。
【産業上の利用可能性】
【００５８】
本発明の利用分野としては，遠方を飛行する鳥類の飛行経路の計測を必要としている分野がある。
【図面の簡単な説明】
【００５９】
【図１】本発明にかかる飛行経路計測システムの一実施例の全体構成図である。
【図２】処理手段における適用例の一つである相互相関法の説明図である。
【図３】本計測システムを用いた計測処理の一例を示すフロー図である。鳥類の二次元の飛行経路表示を得る例である。
【図４】本計測システムを用いた計測処理の一例を示すフロー図である。鳥類の二次元の飛行経路表示を地図上に得る例である。
【図５】本計測システムを用いた計測処理の一例を示すフロー図である。鳥類の二次元の飛行経路表示を得る他の例である。
【図６】本計測システムを用いた計測処理の一例を示すフロー図である。フローの前半を示す。鳥類の三次元の飛行経路表示を得る例である。
【図７】本計測システムを用いた計測処理の一例を示すフロー図である。フローの後半を示す。鳥類の三次元の飛行経路表示を得る例である。
【図８】表示手段に表示された飛行経路の一例である。
【図５】
0005051540.tif 000002
【図１】
0005051540.tif 000003
【図２】
0005051540.tif 000004
【図３】
0005051540.tif 000005
【図４】
0005051540.tif 000006
【図６】
0005051540.tif 000007
【図７】
0005051540.tif 000008
【図８】
0005051540.tif 000009