(57)【要約】  【課題】夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて鳥類が飛行している場合でも、その鳥類の数をより正確に測定することができる鳥類の観測方法を提供する。  【解決手段】静止画像抽出ステップと、ノイズ除去ステップ及び面積・重心算定ステップと、移動軌跡補間ステップと、移動軌跡生成ステップと、個体数算定ステップ、を経て鳥類の数を算定するようにした鳥類の観測方法。【選択図】図３
2016090409.tif 000002
【特許請求の範囲】
【請求項１】
船舶用レーダーの空中線部を鉛直面に沿って回転させ、該空中線部から発射されて鳥類に当たって反射された電波を捕捉することにより、該鳥類を観測するようにした鳥類の観測方法において、
前記レーダーによる１回転分の走査画面を静止画像として抽出する静止画像抽出ステップと、
前記抽出された静止画像を二値化してノイズを除去するノイズ除去ステップと、
前記ノイズを除去した前記静止画像上のレーダー反応をポリゴン化して面積を算定するとともに前記ポリゴンの重心を算定する面積・重心算定ステップと、
前記レーダーが１回転する毎に前記静止画像抽出ステップから前記ポリゴンの面積・重心算定ステップまでの各処理を行い、前記各ポリゴンにおける重心の移動軌跡を描く移動軌跡生成ステップと、
その際に前記移動軌跡生成ステップで描かれた前記ポリゴンの重心の移動軌跡が途切れている場合に、その前のｔ−１回転目の移動ベクトルを延長して移動軌跡を補間する移動軌跡補間ステップと、
前記移動軌跡生成ステップで描かれた移動軌跡の数に１レーダー反応当たりの平均的な個体数を乗じ、前記鳥類の個体数を算定する個体数算定ステップと、
を経て前記鳥類の数を算定する、ことを特徴とする鳥類の観測方法。
【請求項２】
前記静止画像抽出ステップは、前記レーダーの１回転分の走査画面をＲＧＢカラー化された静止画像として抽出し、その色情報を参照してレーダー反応以外を静止画像から除去し、そのＲＧＢカラー化された静止画像をグレースケール化した静止画像に変換する、ことを特徴とする請求項１に記載の鳥類の観測方法。
【請求項３】
前記ノイズ除去ステップは、前記グレースケール化した静止画像を低、高それぞれ異なる閾値で各々二値化し、前記モルフォロジー演算によりそれぞれ低閾値でのノイズ除去と高閾値でのノイズ除去を行う、ことを特徴とする請求項１に記載の鳥類の観測方法。
【請求項４】
前記面積・重心算定ステップは、前記低閾値で二値化されたポリゴンを対象に、該ポリゴンの面積ａＬが所定の面積Ａｍｉｎ、Ａｍａｘに対して、Ａｍｉｎ≦面積ａＬ≦Ａｍａｘを満たすポリゴンの重心を抽出し、さらに所定の面積Ａに対してａＬ≧Ａを満たす場合は該低閾値で二値化されたポリゴンの重心を前記高閾値で二値化されたポリゴンから抽出した重心に補正する、ことを特徴とする請求項１に記載の鳥類の観測方法。
【請求項５】
前記移動軌跡補間ステップは、前記移動軌跡生成ステップで描かれた前記ポリゴンの重心の移動軌跡がｔ回転目で途切れている場合に、その前のｔ−１回転目の移動ベクトルを延長して移動軌跡を補間する、ことを特徴とする請求項１に記載の鳥類の観測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は鳥類の観測方法に関するものであり、特に、レーダーの走査画面に現れる鳥類の個体数や群れ数、及び飛行軌跡の観測を可能にし得る鳥類の観測方法に関するものである。
【背景技術】
【０００２】
鳥類は、個体の生息や種の維持のため、巣やねぐらと採餌地との間を往来するばかりでなく、繁殖地と越冬地との間の長い距離を定期的に渡る等、飛翔行動に大きな生理生態的特徴を持つ。このため、環境アセスメント等において鳥類の生息保全を検討する際には、飛翔行動状況の把握が必要になるケースが多い。その理由として、送電線等の高くて視認性が悪い構造物や発電用風車の設置、飛行場の建設等が、バードストライクを引き起こしたり、渡りを阻害したりすることがあり、そのような鳥類への影響の検討に際しては、事業計画空間と鳥類飛翔空間との重なりの有無や程度を確認することが求められる。
【０００３】
例えば飛行場の建設予定地周辺等の比較的狭い空間における鳥類の飛翔行動を観測する場合、従来は、測風経緯儀又はレーザー測遠システムによる観測が行われていた。測風経緯儀とは、目標の方位角及び鉛直角を２地点から同時に計測し、その交点を目標の位置として求めるものである。また、レーザー測遠システムは、１地点から目標にレーザー光を当て、目標の方位角、鉛直角、及び、斜距離を計測し、目標の位置を求めるものである。しかしながら、測風経緯儀は、２地点からの視準線が交差した場合にのみ観測値が得られ、レーザー測遠システムは、目標である鳥にレーザー光が当たった場合にのみ観測値が得られる。したがって、いずれも小形で速く飛んだり複雑な軌跡で飛ぶものの観測は困難で、観測できない個体が多数生じたり、鳥が見えない夜間や高高度域の観測はできないという問題があった。
【０００４】
一方、大陸間の渡りなど広範囲の鳥類の飛翔行動を研究するために、気象用ドップラーレーダーを用いて、レーダーから発射され鳥類に当たって跳ね返ってきた電波を捉えることにより、目視観察など通常の方法では到底把握することができない広い範囲や、夜間における鳥類の動きを観測する手法もある。しかしながら、気象用ドップラーレーダーは、観測範囲は広いものの狭い範囲の分解能が低く、また大型で運搬が困難であるとともに非常に高価であるため、必要な場所で手軽に観測に用いることができないという問題があった。
【０００５】
そこで、船舶用レーダーを使用し、その船舶用レーダーの空中線部を鉛直面に沿って回転させ、該空中線部から発射されて鳥類等の飛翔性動物に当たって反射された電波を捕捉することにより、該飛翔性動物を観測するようにした鳥類等の飛翔性動物の観測方法も知られている（例えば、特許文献１参照）。この観測方法では、船舶用レーダーをその空中線部を鉛直面に沿って回転させて用いるので、夜間や高高度域であっても鳥類等の飛翔性動物を高精度で観測することができる。また、船舶用レーダーは比較的小型で安価であるので、必要な場所に運搬して手軽に観測を行うことができるなどの効果がある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−９６１０３号公報。
【非特許文献１】Sethi, K.I. and R. Jain. Finding trajectories of feature points in a monocular image sequence. IEEE Trans. Pattern Analysis and Machine Intelligence, 9, 56-73, 1987.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしなから、特許文献１で知られる船舶用レーダーを使用した鳥類等の飛翔性動物の観測方法は、鳥類が群れをなして飛行している場合や飛行軌跡が複雑な場合、夜間や高高度にあってはその鳥類の数をより正確に測定することが難しいという問題点があった。
【０００８】
そこで、夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて鳥類が飛行している場合でも、その鳥類の数をより正確に測定することができる鳥類の観測方法を提供するために解決すべき技術的課題が生じてくるのであり、本発明はこの課題を解決することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は上記目的を達成するために提案されたものであり、請求項１に記載の発明は、船舶用レーダーの空中線部を鉛直面に沿って回転させ、該空中線部から発射されて鳥類に当たって反射された電波を捕捉することにより、該鳥類を観測するようにした鳥類の観測方法において、前記レーダーによる１回転分の走査画面を静止画像として抽出する静止画像抽出ステップと、前記抽出された静止画像を二値化してノイズを除去するノイズ除去ステップと、前記ノイズを除去した前記静止画像上のレーダー反応をポリゴン化して面積を算定するとともに前記ポリゴンの重心を算定する面積・重心算定ステップと、前記レーダーが１回転する毎に前記静止画像抽出ステップから前記ポリゴンの面積・重心算定ステップまでの各処理を行い、前記各ポリゴンにおける重心の移動軌跡を描く移動軌跡生成ステップと、その際に前記移動軌跡生成ステップで描かれた前記ポリゴンの重心の移動軌跡がｔ回転目で途切れている場合に、その前のｔ−１回転目の移動ベクトルを延長して移動軌跡を補間する移動軌跡補間ステップと、前記移動軌跡生成ステップで描かれた移動軌跡の数に１レーダー反応当たりの平均的な個体数を乗じ、前記鳥類の個体数を算定する個体数算定ステップと、を経て前記鳥類の数を算定する、ことを特徴とする。
【００１０】
この方法によれば、比較的安価で小型の船舶用レーダーを必要な場所に運搬して、船舶用レーダーの空中線部を鉛直面に沿って回転させることにより、一回転で全方位を走査することができる。また、その際、走査画面より抽出した静止画像のノイズを除去した後、レーダー反応をポリゴン(多角形)化し、該ポリゴンの面積と重心を算定して、空中線部が一回転する度に変わる該重心の位置を滑らかに結ぶことにより該ポリゴンの移動軌跡を描くことができる。また、その移動軌跡が途切れている場合には、その前の移動ベクトルを延長して移動軌跡を補間し、該移動軌跡の断片化を軽減しポリゴンの重心の位置を更に滑らかに結び、精度を高めることができる。そして、ここで描かれた移動軌跡の数に１レーダー反応当たりの平均的な個体数を乗じ、鳥類の個体数を算定することにより、例え、夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて鳥類が飛行している場合であっても、その鳥類の数をより正確に測定することができる。なお、ここでの鳥類とは、小形から大形の鳥、コウモリ類、及び小形の鳥以上の表面積を持つ飛翔性動物を含むものである。
【００１１】
請求項２に記載の発明は、請求項１記載の方法において、前記静止画像抽出ステップは、前記レーダーの１回転分の走査画面をＲＧＢカラー化された静止画像として抽出し、その色情報を参照してレーダー反応以外を静止画像から除去し、そのＲＧＢカラー化された静止画像をグレースケール化した静止画像に変換する、ことを特徴とする。
【００１２】
この方法によれば、レーダーの１回転分の走査画面をＲＧＢカラー化された静止画像として抽出し、その色情報を参照してレーダー反応以外を静止画像から除去し、そのＲＧＢカラー化された静止画像をグレースケール化し静止画像に変換することにより、レーダー反応のみに注目できる。そして、精度を更に向上させ、その鳥類の数をより正確に測定することができる。
【００１３】
請求項３に記載の発明は、請求項１に記載の方法において、前記ノイズ除去ステップは、前記グレースケール化した静止画像を低、高それぞれ異なる閾値で各々二値化し、前記モルフォロジー演算によりそれぞれ低閾値でのノイズ除去と高閾値でのノイズ除去を行う、ことを特徴とする。
【００１４】
この方法によれば、グレースケール化した静止画像を、低、高それぞれ異なる閾値で各々二値化し、更に前記モルフォロジー演算によりそれぞれ低閾値でのノイズ除去と高閾値でのノイズ除去を行うことにより、レーダー反応の強弱を反映させてポリゴンを識別し易くする。そして、精度を更に向上させて、その鳥類の数をより正確に測定することができる。
【００１５】
請求項４に記載の発明は、請求項１に記載の方法において、前記移動軌跡生成ステップは、前記低閾値で二値化されたポリゴンを対象に、該ポリゴンの面積ａＬが所定の面積Ａｍｉｎ、Ａｍａｘに対して、Ａｍｉｎ≦面積ａＬ≦Ａｍａｘを満たすポリゴンの重心を抽出し、さらに所定の面積Ａに対してａＬ≧Ａを満たす場合は該低閾値で二値化されたポリゴンの重心を前記高閾値で二値化されたポリゴンから抽出した重心に補正する、ことを特徴とする。
【００１６】
この方法によれば、低閾値で二値化されたポリゴンの重心を高閾値で二値化されたポリゴンから抽出した重心に補正することにより、精度が高まり、その鳥類の数を更に正確に測定することができる。
【００１７】
請求項５に記載の発明は、請求項１に記載の方法において、前記移動軌跡補間ステップは、前記移動軌跡生成ステップで描かれた前記ポリゴンの重心の移動軌跡がｔ回転目で途切れている場合に、その前のｔ−１回転目に得られた移動ベクトルを延長して移動軌跡を補間する、ことを特徴とする。
【００１８】
この方法によれば、移動軌跡を描く際、ポリゴンの重心の移動軌跡が途切れている場合に、その前の移動ベクトルを延長して補間し、その移動軌跡の断片化を軽減しポリゴンの重心の位置を滑らかに結ぶことにより、精度を更に高めて、その鳥類の数をより正確に測定することができる。
【発明の効果】
【００１９】
本発明による鳥類の観測方法によれば、比較的安価で小型の、船舶用レーダーを必要な場所に運搬して、その運搬した場所で、船舶用レーダーの空中線部を鉛直面に沿って回転させることにより、夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて飛行している鳥類の数をより正確に測定することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の一実施形態に係る鳥類の観測方法を実施するための観測装置のブロック構成図である。
【図２】本発明の一実施形態に係る鳥類の観測方法を概念的に説明するための図である。
【図３】本発明の一実施形態に係る鳥類の観測方法における処理手順の一例を示す解析フローである。
【図４】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、指示部の解析画面のハードコピーの一例を示す図である。
【図５】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、（ａ）は上記ハードコピーのレーダー反応のみを抽出した一例を示す図、（ｂ）はグレースケール化した一例を示す図である。
【図６】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、（ａ）はレーダー反応を高閾値で二値化した後における静止画像の一例を示す図、（ｂ）は前記静止画像をモルフォロジー演算によって微小なノイズを除去した静止画像の一例を示す図である。
【図７】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、（ａ）はレーダー反応を低閾値で二値化した静止画像の一例を示す図、（ｂ）は前記静止画像をモルフォロジー演算によって微小なノイズを除去した静止画像の一例を示す図である。
【図８】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、（ａ）は静止画像上のレーダー反応をポリゴン化した一例を示す図、（ｂ）はポリゴンの重心を抽出する一例を示す図である。
【図９】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、ポリゴンの軌跡を作成するハードコピーの一例を示す図である。
【図１０】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、（ａ）は低閾値で二値化したポリゴンの一例、（ｂ）は高閾値で二値化したポリゴンの一例を各々示す図である。
【図１１】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、移動軌跡を作成する一例を説明する図である。
【図１２】本発明の一実施形態に係る鳥類の観測方法の解析イメージ図で、途切れた移動軌跡を補間する一例を説明する図である。
【発明を実施するための形態】
【００２１】
本発明は夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて鳥類が飛行している場合でも、その鳥類の数をより正確に測定することができる鳥類の観測方法を提供するという目的を達成するために、船舶用レーダーの空中線部を鉛直面に沿って回転させ、該空中線部から発射されて鳥類に当たって反射された電波を捕捉することにより、該鳥類を観測するようにした鳥類の観測方法において、前記レーダーによる１回転分の走査画面を静止画像として抽出する静止画像抽出ステップと、前記抽出された静止画像を二値化してノイズを除去するノイズ除去ステップと、前記ノイズを除去した前記静止画像をポリゴン化して面積を算定するとともに前記ポリゴンの重心を算定する面積・重心算定ステップと、前記レーダーが１回転する度に前記静止画像抽出ステップから前記ポリゴンの面積・重心算定ステップまでの各処理を行い、前記各ポリゴンにおける重心の移動軌跡を描く移動軌跡生成ステップと、その際に前記移動軌跡生成ステップで描かれた前記ポリゴンの重心の移動軌跡がｔ回転目で途切れている場合に、その前のｔ−１回転目の移動ベクトルを延長して移動軌跡を補間する移動軌跡補間ステップと、前記移動軌跡生成ステップで描かれた移動軌跡の数に１レーダー反応当たりの平均的な個体数を乗じ、前記鳥類の個体数を算定する個体数算定ステップと、を経て前記鳥類の数を算定する、ようにして実現した。
【００２２】
以下、本発明の実施形態による鳥類の観測方法の好適な一実施例を、図面を参照しながら詳細に説明する。
【実施例】
【００２３】
図１及び図２は本発明に係る鳥類の観測方法を概念的に説明するための図で、図１はその観測装置のブロック構成図、図２はその観測説明図である。
【００２４】
図１において、観測装置１１は、鳥類、すなわち小形から大形の鳥、コウモリ類、及び小形の鳥以上の表面積を持つ飛翔性動物(以下、これらを総称して「鳥類」という)を観測するためのものであり、船舶用レーダー１２と、ＧＰＳ受信機１３と、軌跡解析用プロッター１４と、プリンター１５と、録画装置１６と、電源１７とを備えている。
【００２５】
一般に船舶用レーダー１２とは、衝突防止及び自船位置の測定のために船舶に設置され、他船、島、陸岸、航路標識等の物標を検出するためのパルスレーダーであり、電波法で定められた、３ＧＨｚ帯〜９ＧＨｚ帯内の指定周波数帯を用いたパルスレーダーを使用する。
【００２６】
さらに、前記船舶用レーダー１２は、空中線部２１と指示部２２とを備えている。空中線部２１は、長い箱型の形状をなし、中央部に接続されたモーターにより回転しながら、長手方向に沿った縁部からパルス波を発射するのと同時に、飛翔する鳥類からの反射波を受信し、受信信号として指示部２２に出力する。指示部２２は、空中線部２１から入力された受信信号を反射波の処理回路により処理して物標位置を演算し、図２に示すように、空中線部２１からの電波で捕捉された物標の位置を表示部２２ａに表示する。
【００２７】
前記ＧＰＳ受信機１３は、受信したＧＰＳ衛星からの信号に基づいて自機位置の経度緯度、標高を演算し、その情報を軌跡解析用プロッター１４に出力するものである。
【００２８】
前記軌跡解析用プロッター１４は、例えばパーソナルコンピュータであって、情報を決められた手順に従って処理する解析アルゴリズムを実施するためのプログラムが内蔵され、またプログラムで処理されたデータを読み書き可能に保持しておくためのメモリ２３を有している。その軌跡解析用プロッター１４は、船舶用レーダー１２から入力された情報と、ＧＰＳ受信機１３から入力された情報とに基づいて、捕捉された鳥類の位置及び個体数、群れ数、飛行軌跡の空間座標を演算して解析するとともに、その解析結果を出力して判りやすく表示部２２ａに表示するものである。
【００２９】
前記プリンター１５は、軌跡解析用プロッター１４からのデータに基づいて指示部２２の表示部２２ａに表示される画像等を印刷するためのカラープリンターである。
【００３０】
前記録画装置１６は、軌跡解析用プロッター１４からのデータに基づいて指示部２２の表示部２２ａに表示される画像等を録画するためのレコーダーである。
【００３１】
前記電源１７は、上記の各機器の電源となるもので、発電機と該発電機で生成された電力を蓄積しておくためのバッテリーを備えている。
【００３２】
また、観測装置１１は、図２に示すように、空中線部２１を鉛直面に沿うように支持する空中線部架台１８を備えている。
【００３３】
図３は、軌跡解析用プロッター１４が、船舶用レーダー１２から入力された情報とＧＰＳ受信機１３から入力された情報とに基づき、捕捉された鳥類の位置及び個体数、群れ数、飛行軌跡の空間座標をメモリ２３内のプログラムに従って演算して解析する手順の一例を説明する解析フローである。また、図４〜図９は指示部２２の解析画面のハードコピーの一例を示すものである。
【００３４】
そこで、図１及び図２示した鳥類の観測装置１１を用いて捕捉された鳥類の位置及び個体数、群れ数、飛行軌跡の空間座標を、メモリ２３内のプログラムに従って解析する方法の一例を、図１〜図９と共に以下説明する。
【００３５】
まず、観測装置１１の概略動作を説明すると、本実施形態の観測方法では、観測装置１１の船舶用レーダー１２を、観測対象範囲、例えば風力発電機を設置する範囲やその近傍に設置して行う。通常、船舶用レーダー１２は、空中線部２１を水平面に沿って回転させ、船舶用レーダー１２から物標までの距離と方位角を計測し、物標の平面的な位置を観測する。しかし、本実施形態の観測方法では、空中線部２１を、空中線架台１８に支持させて、水平面に対して垂直になるように立ててから、図２の一点鎖線矢印に示すように回転させる。すなわち、空中線部２１を鉛直面に沿って回転させ、空中線部２１から所定のビーム幅で電波を発射する。
【００３６】
すると、空中線部２１は、その空中線部２１の設置地点を中心とする一定幅の探査断面を空域に連続的に作り出し、その断面内を飛翔する鳥類等からの反射した電波を受信して、指示部２２に表示する。
【００３７】
そして、ここでの船舶用レーダー１２は、空中線部２１を鉛直面に沿って回転させることにより、一回転で全方位を操作し、空中線部２１の設定地点を中心とする一定幅の空域に探査断面を連続的に作り出す。その際に指示部２２の走査画面（解析画面)上に現われるレーダー反応を重心で代表し、回転の度に変わるその位置を滑らかに結ぶことによって軌跡を作成する。また、レーダー反応が微弱で一時的に位置を見失う場合に備えた、軌跡の直前の状況から位置を予測し補間する方法と、レーダー反応の強弱を考慮した重心の補正方法を組み合わせる。なお、滑らかさの評価には次式（ａ）に示すパスコヒーレンス関数を用いる（例えば、非特許文献１参照）。この関数は重心間が滑らかに繋がるほど０に近づく。角度項、速度項の重みは状況に応じて調節する。
【００３８】
【数１】
2016090409.tif 000003
【００３９】
更に、上記観測装置１１を使用して鳥類の数を算定する解析アルゴリズムを、図３に示す解析フロー(Ｓ１０１〜Ｓ３１６)と共に、ステップ(１)〜(１４)の順に従って、図４〜図９に示す解析イメージ図を参照しつつ、また、必要に応じて図１０〜図１２の模式図を加えて説明する。
【００４０】
ここでの解析アルゴリズムは、大別すると、静止画像抽出ステップ(１)、ノイズ除去ステップ及び面積・重心算定ステップ(２)〜(５)、移動軌跡補間ステップ(６)、移動軌跡生成ステップ(７)〜(１３)、個体数算定ステップ(１４)からなる。
【００４１】
［静止画像抽出ステップ］
ステップ（１）：船舶用レーダー１２の空中線部２１を鉛直面に沿って回転させ、該空中線部２１から発射されて鳥類に当たって反射された電波を捕捉することにより、ｔ＋１回転目の走査画面を静止画像として抽出する(Ｓ１０１)。図４は、その抽出された静止画像、すなわち指示部２２の解析画面のハードコピーの一例を示す。この静止画像は、ＲＧＢカラー化されて抽出される。
【００４２】
ステップ（２）：上記静止画像の色情報を参照し、図５(ａ)に示すようにレーダー反応以外を静止画像から除去し(Ｓ１０２)、更に図５(ｂ)に示すようにそのＲＧＢカラー化された静止画像をグレースケール化した静止画像に変換する(Ｓ１０３)。
【００４３】
［ノイズ除去ステップ及び面積・重心算定ステップ］
ステップ（３）：上記グレースケール化した静止画像を図６(ａ)に示すように高閾値ＫＨにて二値化し、次いで図６(ｂ)に示すようにモルフォロジー演算(収縮Ｎ回・膨張Ｎ回)によって微小なノイズを除去する。また、その後、図８(ａ)に示すように上記静止画像をポリゴン(多角形)化し、図８（ｂ）に示すように、そのポリゴンＰＨの重心ＣＨを抽出し、そのポリゴンＰＨの重心ＣＨをデータ化する(Ｓ１０９〜Ｓ１１５)。
【００４４】
ステップ（４）：上記ステップ（３）と同時に、グレースケール化した静止画像を図７（ａ）に示すように低閾値ＫＬにて二値化し、次いで図７（ｂ）に示すようにモルフォロジー演算(収縮Ｎ回・膨張Ｎ回)によって微小なノイズを除去する。また、その後、ステップ（３）での高閾値ＫＨの場合と同様に、図８（ａ）に示すようにポリゴン化し、そのポリゴンＰＬの面積ａＬを算定する (Ｓ１０４〜Ｓ１０８)。
【００４５】
すなわち、グレースケール化した静止画像を上記ステップ（３）と（４）において、高閾値ＫＨと低閾値ＫＬでそれぞれ二値化すると、図１０の（ａ）に低閾値で二値化されたポリゴンＰＬ、（ｂ）に高閾値で二値化されたポリゴンＰＨをそれぞれ示すように、グレースケール化した静止画像にレーダー反応の強弱を反映させて、ポリゴンを細かく識別し、これによって精度の向上を図ることができる。図１０の（ａ）では低閾値で二値化された１個のポリゴンＰＬが抽出されている。しかし、低閾値で二値化された段階では単一のポリゴンとして対象とできるかどうかが不明確であるが、更に高閾値で二値化することにより、図１１の（ｂ）に示すように２個のポリゴンＰＨに変化し対象とするポリゴンが明確になる。
【００４６】
ステップ（５）：ステップ(４)のポリゴンＰＬを対象に面積ａＬが所定の面積Ａｍｉｎ、Ａｍａｘに対してＡｍｉｎ≦ａＬ≦Ａｍａｘを満たすポリゴンＰＬの重心ＣＬを抽出する。そのうち、所定の面積Ａに対してａＬ≧Ａも満たす場合は、ステップ(３)のポリゴンＰＨから抽出した重心ＣＨに、ポリゴンの重心Ｃを補正する(Ｓ１１７〜Ｓ１２１)。
【００４７】
［移動軌跡補間ステップ］
ステップ（６）：図１２に示すように、対象ｉ(本実施例では各ポリゴン)、すなわち各ポリゴンの軌跡Ｍ(ｉ)がｔ−１回転目の位置Ｐｉ，ｔ−１で途切れている場合、Ｐｉ，ｔ−１に移動ベクトルＰｉ，ｔ−２Ｐｉ，ｔ−１を適用し、図１２に示すようにポリゴンの位置Ｐｉ，ｔ−１の次の位置Ｐｉ，ｔを予測して補間する(Ｓ２０１〜Ｓ２０４)。
【００４８】
［移動軌跡生成ステップ］
ステップ（７）：そして、ステップ（５）の重心Ｃを中心とする半径Ｒの円内に、ｔ回転目での対象ｉの位置Ｐｉ，ｔがある軌跡Ｍ（ｉ）を探索する(Ｓ２０５、Ｓ２０６)。
【００４９】
ステップ（８）：ステップ（７）の軌跡Ｍ（ｉ）のＰｉ，ｔを中心とする半径Ｒの円内に位置する重心Ｃを重心Ｃ’として探索する(Ｓ２０７、Ｓ２０８)。
【００５０】
ステップ（９）：各軌跡について、Ｐｉ，ｔ−１、Ｐｉ，ｔ、Ｃ’によるパスコヒーレンス関数が最も小さくなる重心Ｃ’を選定する(Ｓ２０９〜Ｓ２１２)。
【００５１】
ステップ（１０）：重心Ｃ’として重心Ｃが重複せず選定された場合は、重心Ｃを対象ｉの軌跡Ｍ（ｉ）にＰｉ，ｔ＋１として追加する。重複した場合は、パスコヒーレンス関数が最も小さい軌跡にＰｉ，ｔ＋１として追加し、いずれの軌跡にもＰｉ，ｔ＋１として追加されなかった重心Ｃは新たな軌跡の始点に設定する(Ｓ３０１〜Ｓ３０５)。
【００５２】
ステップ（１１）：ステップ（６）にて補間された対象ｉの軌跡Ｍ（ｉ）にいずれの重心ＣもＰｉ，ｔ−１として追加されなかった場合は補間したＰｉ，ｔを削除する（Ｓ３０６〜Ｓ３０９）。
【００５３】
ステップ（１２）：走査が継続する場合は次の回転時に進み（Ｓ３１１）、終了する場合はステップ（１３）へ移行する（Ｓ３１０）。
【００５４】
ステップ（１３）：軌跡Ｍ（ｉ）の長さｌが所定の長さＬに対して小さい、ｌ＜Ｌを満たす軌跡を削除する(Ｓ３１２〜Ｓ３１６)。
【００５５】
［個体数算定ステップ］
ステップ（１４）：回転毎の重心及び軌跡の数に１レーダー反応当りの平均的な個体数を乗じ、瞬間個体数及び時間個体数を算定する。
【００５６】
ただし、ｔ＝０(１回転目)では抽出した重心Ｃを新たな軌跡の始点に設定するのみとし、ｔ＝１(２回転目)では予測位置による補間は適用せず、パスコヒーレンス関数の代わりに式（ｂ）に示すノルム関数を用い、滑らかさではなく近さで評価する。
【００５７】
【数２】
2016090409.tif 000004
【００５８】
以上説明したように、本実施例によれば、比較的安価で小型の船舶用レーダー１２を必要な場所に運搬して、船舶用レーダー１２の空中線部２１を鉛直面に沿って回転させ、前記静止画像抽出ステップ（１）、ノイズ除去ステップ及び面積・重心算定ステップ（２）〜（５）、移動軌跡補間ステップ（６）、移動軌跡生成ステップ（７）〜（１３）、個体数算定ステップ（１４）、からなる解析アルゴリズムを実施することにより、鳥類が夜間や高高度で群れをなし、また複雑な飛行軌跡を描いて飛行している場合でも、その鳥類の数をより正確に測定することができる。
【００５９】
なお、本発明は、本発明の精神を逸脱しない限り種々の改変を為すことができ、そして、本発明が該改変されたものに及ぶことは当然である。
【産業上の利用可能性】
【００６０】
以上説明したように、本発明は鳥類の観測以外に、類似する飛翔性動物等の観測にも応用できる。
【符号の説明】
【００６１】
１１      観測装置
１２      船舶用レーダー
１３      ＧＰＳ受信器
１４      軌跡解析用プロッター
１５      プリンター
１６      録画装置
１７      電源
１８      空中線部架台
２１      空中線部
２２      指示部
２２ａ    表示部
２３      メモリ
ｔ        走査回数
ＫＬ      低閾値
ＫＨ      高閾値
Ｎ        収縮・膨張回数
ＰＨ      高閾値で得られたポリゴン
ＣＨ      高閾値で得られたポリゴンＰＨの重心
ＰＬ      低閾値で得られたポリゴン
ＣＬ      低閾値で得られたポリゴンＰＬの重心
ａＬ      ポリゴンＰＬの面積
Ａｍｉｎ、Ａｍａｘ、Ａ    面積
Ｃ、Ｃ’  ポリゴンの重心
Ｐｉ      軌跡ｉの構成点
Ｍ(ｉ)    軌跡ｉ
ｌ        軌跡の長さ
Ｌ        長さ
【図１】
2016090409.tif 000005
【図２】
2016090409.tif 000006
【図３】
2016090409.tif 000007
【図８】
2016090409.tif 000008
【図１２】
2016090409.tif 000009
【図４】
2016090409.tif 000010
【図５】
2016090409.tif 000011
【図６】
2016090409.tif 000012
【図７】
2016090409.tif 000013
【図９】
2016090409.tif 000014
【図１０】
2016090409.tif 000015
【図１１】
2016090409.tif 000016