(57)【要約】  移動体観測方法は、所定の時間間隔で撮影された複数の時系列画像Ｆｎを取得し（ステップＳ１０）、移動体Ｐの推定運動に基づいて、時系列画像Ｆｎにおける移動体Ｐの推定位置Ｌｎを求め（ステップＳ２０）、時系列画像Ｆｎのうちの複数の画像から、移動体Ｐの推定位置Ｌｎを含む画像を抽出画像Ｑｎとして抽出し（ステップＳ３０）、移動体Ｐの基準点が一致するように複数の抽出画像Ｑｎを重ね合わせることによってテンプレート画像ＱＴを生成し（ステップＳ４０）、時系列画像Ｆｎの少なくとも１つに対してテンプレート画像ＱＴを用いたテンプレートマッチングを行うことで、時系列画像Ｆｎのそれぞれにおける移動体Ｐの位置ＬＳｎを特定する（ステップＳ５０）ことを含む。
2018230016.tif 000002
【特許請求の範囲】
【請求項１】
所定の時間間隔で撮影された複数の時系列画像を取得し、
移動体の推定運動に基づいて、前記時系列画像における前記移動体の推定位置を求め、
前記時系列画像のうちの複数の画像から、前記移動体の前記推定位置を含む画像を抽出画像として抽出し、
前記移動体の基準点が一致するように複数の前記抽出画像を重ね合わせることによってテンプレート画像を生成し、
前記時系列画像のうちの少なくとも１つの画像に対して前記テンプレート画像を用いたテンプレートマッチングを行うことで、前記時系列画像のうちの前記少なくとも１つの画像における前記移動体の位置を特定する
ことを含む移動体観測方法。
【請求項２】
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記推定位置である、
請求項１に記載の移動体観測方法。
【請求項３】
前記テンプレート画像の各画素値は、複数の前記抽出画像において対応する画素の平均画素値、中央値及び最頻値のうちの何れかに設定される、
請求項２に記載の移動体観測方法。
【請求項４】
複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、
二値化された各前記抽出画像における前記移動体の重心を算出する
ことを更に含み、
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、二値化された前記移動体の前記重心である、
請求項１に記載の移動体観測方法。
【請求項５】
複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、
二値化された各前記抽出画像における前記移動体の画素分布を所定の二方向に積算した２つの一次元分布に正規関数のフィッティングを行うことで前記移動体の最大輝度位置を算出する
ことを更に含み、
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記最大輝度位置である、
請求項１に記載の移動体観測方法。
【請求項６】
前記移動体は、地球軌道上のデブリ又は液体中の微生物である、
請求項１に記載の移動体観測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、揺らぎ等の影響を伴う移動体の観測方法に関する。
【背景技術】
【０００２】
従来、ＣＣＤカメラ等の機器が、所定の物理法則に則って（例えば一定の速度で）移動する観測対象としての物体（以下、「移動体」と称する。）を撮影するとき、例えば、移動体が遠距離にある場合、移動体の大きさが小さい場合、移動体の明るさが微弱な場合には、移動体が背景のノイズに埋もれてしまい、見つけ出すことが困難であった。特に、移動体の速度が不明なため、長期間露光を行なっても、移動体に合わせてカメラの向きを動かすことができなかった。そのため、移動体が視野内で流れてしまい、映し出すことができなかった。【０００３】
かかる問題点に対して、特許文献１は、時系列で取得した複数枚の画像を用いてノイズを除去する移動体検出方法を提案している。この方法では、一定の間隔で撮影した多数の画像を利用し、複数の連続した画像列毎に以下の処理を行っている。（ａ）ある速度で移動体が等速度運動していると仮定し、その仮定速度で画像をずらしながら重ね合わせ、同じ場所のデータに共通した僅かに明るい部分を見つけ出し、Ｓ／Ｎ比（信号雑音比）を向上させる。（ｂ）想定される移動体の速度範囲の全域に渡って仮定速度を変化させて前記（ａ）の処理を繰り返し行い、明るい点が見つかった際の速度及び位置を検出結果として取得する。（ｃ）画像列毎に前記（ａ）及び（ｂ）の処理により得られた結果から、それらが互いに矛盾しない（位置・速度が不連続でない）ものを同じ移動体であるとして統合する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−５１２１１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述した特許文献１の移動体検出方法はいわゆる重ね合わせ法を用いた手法であり、多数の画像を用いて移動体の位置を検出することができる。しかしながら、この方法は、移動体の軌跡や挙動を把握することができない。例えば、重ね合わせ法において、Ｎ枚の画像を用いた場合であっても、ｎ枚目（ｎ＜Ｎ）の画像における移動体の位置を検出することができるだけであって、残りの（Ｎ−１）枚の画像における移動体の位置を検出することができない。したがって、重ね合わせ法を用いた場合には、時間に対して非常に疎な情報しか得ることができない。【０００６】
本開示は、上述した問題点に鑑みて創案されたものであり、重ね合わせ法を用いて移動体を精度良く検出することが可能な移動体観測方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本開示の一態様は移動体観測方法であって、所定の時間間隔で撮影された複数の時系列画像を取得し、移動体の推定運動に基づいて、前記時系列画像における前記移動体の推定位置を求め、前記時系列画像のうちの複数の画像から、前記移動体の前記推定位置を含む画像を抽出画像として抽出し、前記移動体の基準点が一致するように複数の前記抽出画像を重ね合わせることによってテンプレート画像を生成し、前記時系列画像のうちの少なくとも１つの画像に対して前記テンプレート画像を用いたテンプレートマッチングを行うことで、前記時系列画像のうちの前記少なくとも１つの画像における前記移動体の位置を特定することを含むことを要旨とする。【０００８】
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記推定位置でもよい。【０００９】
前記テンプレート画像の各画素値は、複数の前記抽出画像において対応する画素の平均画素値、中央値及び最頻値のうちの何れかに設定されていてもよい。【００１０】
前記移動体観測方法は、複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、二値化された各前記抽出画像における前記移動体の重心を算出することを更に含んでもよい。この場合、前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、二値化された前記移動体の前記重心でもよい。【００１１】
前記移動体観測方法は、複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、二値化された各前記抽出画像における前記移動体の画素分布を所定の二方向に積算した２つの一次元分布に正規関数のフィッティングを行うことで前記移動体の最大輝度位置を算出すること、を更に含んでもよい。この場合、前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記最大輝度位置であってもよい。【００１２】
前記移動体は、地球軌道上のデブリ又は液体中の微生物であってもよい。
【発明の効果】
【００１３】
本開示によれば、重ね合わせ法を用いて移動体を精度良く検出することが可能な移動体観測方法を提供することができる。
【図面の簡単な説明】
【００１４】
【図１】本開示の一実施形態に係る移動体観測装置のブロック図である。
【図２】本実施形態に係る移動体観測方法のフローチャートである。
【図３】本実施形態に係る移動体観測方法の各ステップを説明するための図である。
【図４】本実施形態に係る移動体観測方法によって得られた画像の強度分布の一例であり、（ａ）は一枚の抽出画像の強度分布を示し、（ｂ）はテンプレート画像ＱＴの強度分布を示す。
【図５】図５は本実施形態の第１変形例に係るフローチャートである。
【図６】図６は本実施形態の第２変形例に係るフローチャートである。
【発明を実施するための形態】
【００１５】
以下、本開示の実施形態について図面を参照して説明する。図１は、本実施形態に係る移動体観測装置のブロック図である。この図に示すように、移動体観測装置１０は、観測部２０と、制御部３０とを備えている。観測部２０は、光学系２１と、撮像部２２とを備えている。光学系２１は、観測対象である移動体Ｐを含む光学像を形成するものであり、例えば望遠鏡や顕微鏡である。撮像部２２は光学系２１によって生成された光学像を画像データとして取得するものであり、例えばＣＣＤカメラである。制御部３０は、ＣＰＵ（演算部）３１と、記憶部３２と、入出力部（Ｉ／Ｏ）３３とを備えている。ＣＰＵ３１は、例えば、記憶部３２に記憶されたプログラムやデータ等を用いて種々の演算を実行する。記憶部３２は、本実施形態に係るプログラムや種々のデータを記憶するものであり、半導体メモリやハードディスクドライブ（ＨＤＤ）などで構成される。例えば、記憶部３２は、観測部２０から送られてきた画像データも記憶する。入出力部３３は、外部機器との間でデータや指令等の通信を行うものであり、本実施形態では、観測部２０の撮像部２２に接続している。従って、撮像部２２によって取得された画像データは、撮像部２２と入出力部３３の通信を介して記憶部３２に記憶され、ＣＰＵ３１の演算に用いられる。
【００１６】
図２は、本実施形態に係る移動体観測方法のフローチャートである。図３は本実施形態に係る移動体観測方法の各ステップを説明するための図である。これらの図に示すように、本実施形態に係る移動体観測方法は、画像取得ステップ（ステップＳ１０）と、位置推定ステップ（ステップＳ２０）と、画像抽出ステップ（ステップＳ３０）と、テンプレート生成ステップ（ステップＳ４０）と、位置特定ステップ（ステップＳ５０）とを含む。
【００１７】
観測対象である移動体Ｐは、例えば、地球軌道上のデブリである。地球軌道上には、軍事衛星、通信衛星、科学衛星、観測衛星、航行衛星等、種々の目的の人工衛星が周回している。これらの人工衛星は、故障して機能しなくなったり、役目を終えて寿命に達したりした場合には、そのまま軌道上に放置され、軌道上を移動するデブリ（宇宙ごみ、スペースデブリとも言う。）となることが多い。また、人工衛星等の打ち上げに使用したロケット等の残骸もデブリとして軌道上に放置されている。
【００１８】
現在、軌道上を周回しているデブリは、数千個以上にも及び、自然衝突により個数が増加する自己増殖の段階に突入している。このデブリは、現在使用中の人工衛星の軌道周回中や人工衛星等の打ち上げ時に衝突する可能性があることから、高精度な観測方法が求められている。
【００１９】
なお、本実施形態における観察対象は、上述のデブリに限定されない。即ち、移動体観測方法は、所定の時間間隔で取得された画像上の形状に揺らぎが生じる物体を観測対象として適用できる。従って、移動体Ｐは、例えば、薬剤などの液体中の微生物（例えば静止した微生物）であってもよい。この場合、光学系２１として顕微鏡が用いられる。なお、揺らぎの原因は、移動体が遠距離にあること、移動体の大きさが小さいこと、或いは、移動体の明るさが微弱であること、などが挙げられる。
【００２０】
画像取得ステップ（ステップＳ１０）では、所定の時間間隔で撮影された、移動体Ｐを含む複数（例えばＮ枚）の時系列画像Ｆｎ（ｎ＝１〜Ｎ）を取得する。本実施形態では、移動体観測装置１０の観測部２０が、所定の時間間隔で移動体Ｐを含む像を複数枚（例えばＮ枚）撮影し、その画像データを移動体観測装置１０の制御部３０（記憶部３２）に時系列画像Ｆ１〜ＦＮとして入力する。
【００２１】
なお、各時系列画像Ｆｎの取得後、各時系列画像Ｆｎに対して、誤検出の原因となり得る恒星等の明るい画素を除去する処理を行ってもよい。例えば、恒星の移動方向及び移動量に応じて各画像を平行移動させて重ね合わせ、重なった同一画素の画素値を差し引くことで、画像信号に含まれる恒星によるノイズ成分を除去することができる。
【００２２】
位置推定ステップ（ステップＳ２０）では、移動体Ｐの推定運動に基づいて、時系列画像Ｆｎ（ｎ＝１〜Ｎ）における移動体Ｐの推定位置Ｌｎ（ｎ＝１〜Ｎ）を求める。例えば、移動体Ｐが等速直線運動をしていると仮定すれば、その推定運動は以下の数式１のように表現することができる。
【数１】
2018230016.tif 000003
・・・（数式１）
ここで、ＰｏｓＸ（ｋ）は時系列画像Ｆｋ（フレーム番号ｋ）における移動体Ｐの推定位置ＬｋのＸ座標、ＰｏｓＹ（ｋ）は時系列画像Ｆｋにおける移動体Ｐの推定位置ＬｋのＹ座標、ＶｅｌＸは時系列画像Ｆｎ上の移動体Ｐの速度のＸ方向成分、ＶｅｌＹは時系列画像Ｆｎ上の移動体Ｐの速度のＹ方向成分、ｋはＮ以下でｎ以外の正の整数、を示している。ＰｏｓＸ（ｎ）及びＰｏｓＹ（ｎ）は、それぞれ時系列画像Ｆｎにおける移動体Ｐの推定位置ＬｎのＸ座標及びＹ座標であり、取得した時系列画像のうちの何れかから選択される。なお、推定位置Ｌｎは記憶部３２に記憶された移動体Ｐの観測履歴から参照されてもよく、移動体Ｐの推定移動速度と共に改めて算出してもよい。後者の場合、例えば、移動体Ｐの推定移動速度を（−ＶｅｌＸ，−ＶｅｌＹ）〜（ＶｅｌＸ，ＶｅｌＹ）の範囲で順次変えながら各画像を重ね合わせ、各画像上の移動体Ｐの候補の互いの一致度が高くなった時の速度を推定移動速度とし、当該候補を移動体Ｐとして設定する。
【００２３】
移動体Ｐが円運動をしていると仮定すれば、その推定運動は以下の数式２又は数式３のように表現することができる。なお、これらの円運動は、Ｘ方向又はＹ方向の何れかの方向が直線に近似できる曲率の小さい円軌道を示している。ここで、Ａは前記円軌道の中心のＸ座標、Ｂは前記円軌道の中心のＹ座標を示している。また、数式２はＶｅｌＸ≦ＶｅｌＹの場合を示し、数式３はＶｅｌＸ＞ＶｅｌＹの場合を示している。
【数２】
2018230016.tif 000004
・・・（数式２）
【数３】
2018230016.tif 000005
・・・（数式３）
【００２４】
なお、移動体Ｐの推定運動は、上述した等速直線運動や円運動に限定されるものではなく、観測対象に合わせて任意に変更することができる。
【００２５】
画像抽出ステップ（ステップＳ３０）では、時系列画像Ｆｎ（ｎ＝１〜Ｎ）から、移動体Ｐの推定位置Ｌｎを含む画像を抽出画像Ｑｎとして抽出する（切り出す）。抽出画像Ｑｎのサイズは、時系列画像Ｆｎのサイズよりも小さい限り任意である。抽出対象である時系列画像Ｆｎの枚数は任意の複数である。この枚数が増えるほど、後述のテンプレート画像ＱＴにおける移動体Ｐの形状の信頼性が高まる。その一方で、この枚数が増えるほど、制御部３０の処理速度が低下する。従って、抽出に用いられる時系列画像の枚数は、これらの事情を考慮した適当な値が望ましい。本実施形態では、説明の便宜上、取得したＮ枚の画像の全てから抽出画像Ｑｎを切り出すものとする。
【００２６】
テンプレート生成ステップ（ステップＳ４０）では、移動体Ｐの基準点が一致するように複数の抽出画像Ｑｎを重ね合わせることによってテンプレート画像ＱＴを生成する。基準点とは、例えば、抽出画像Ｑｎ上の推定位置Ｌｎである。推定位置Ｌｎは既に位置推定ステップ（ステップＳ２０）で得られているため、制御部３０は新たな演算処理を行うことなく、テンプレート画像ＱＴを生成できる。テンプレート画像ＱＴは、各抽出画像Ｑｎの画素値の平均値を算出することで生成される。例えば、抽出画像ＱｎがＭｘ×Ｍｙサイズの画像である場合、テンプレート画像ＱＴの各画素値Ｔ（ｉ，ｊ）は、次の数式４によって算出される。
【数４】
2018230016.tif 000006
・・・（数式４）
ここで、ＴＱｎ（ｉ，ｊ）は抽出画像Ｑｎ上の位置（ｉ，ｊ）における画素値である。つまり、テンプレート画像ＱＴにおける移動体Ｐ（即ち移動体ＰＴ）を示す画素の強度分布は、抽出画像Ｑ１〜ＱＮにおける移動体Ｐ（即ち移動体Ｐ１〜ＰＮ）を示す画素の強度分布の平均である。
【００２７】
図３に示すように、抽出画像Ｑ１〜ＱＮにおける移動体Ｐ１〜ＰＮの形状は、大気等に起因する揺らぎの影響を受けて絶えず変化する。テンプレート画像ＱＴは、このような形状変化を抑制した状態を示す。即ち、テンプレート画像ＱＴにおける移動体ＰＴは、抽出画像Ｑ１〜ＱＮにおける移動体Ｐ１〜ＰＮから得られた「平均的な」形状をもつ。例えば、抽出画像Ｑ１〜ＱＮに現れる個別の形状変化がテンプレート画像ＱＴでは均され、移動体ＰＴの輪郭において強度の低い「ボケ」として現れる。一方、移動体Ｐ１〜ＰＮに共通する形状はテンプレート画像ＱＴにそのまま残る。換言すれば、移動体ＰＴには、移動体Ｐ１〜ＰＮに共通する形状を示す明るい部分と、移動体Ｐ１〜ＰＮに共通しない形状を示す暗い部分とを含むことになる。これにより、移動体ＰＴが移動体Ｐの本来の形状に近い形状を有することができ、テンプレート画像ＱＴはテンプレートとしての高い信頼性をもつことになる。
【００２８】
上述の通り、テンプレート画像ＱＴの各画素値は、抽出画像Ｑ１〜ＱＮの各画素値の平均値である。従って、抽出画像Ｑ１〜ＱＮのそれぞれに現れていたランダムノイズが、テンプレート画像ＱＴでは低減され、Ｓ／Ｎ比（信号雑音比）が向上する。図４（ａ）は、一枚の抽出画像Ｑｎの強度分布を示し、図４（ｂ）は、テンプレート画像ＱＴの強度分布を示す。図４（ａ）及び図４（ｂ）において強度分布の中央に見える強度の高い部分がそれぞれ移動体Ｐｎ及び移動体ＰＴを示す。図４（ａ）に示すように、移動体Ｐｎの周囲にはランダムノイズが現れている。一方、図４（ｂ）に示すように、移動体ＰＴの周囲では図４（ａ）で見られたランダムノイズが著しく減少している。また、ランダムノイズの減少によって、移動体ＰＴの輪郭を示す強度のバラつきも移動体Ｐｎの輪郭を示す強度に比べて少なくなっている。これもテンプレートとしての高い信頼性に寄与する。
【００２９】
位置特定ステップ（ステップＳ５０）では、時系列画像Ｆ１〜ＦＮのうちの少なくとも１つの画像に対してテンプレート画像ＱＴを用いたテンプレートマッチングを行うことで、テンプレートマッチングが適用された時系列画像における移動体Ｐの位置を特定する。このテンプレートマッチングは、位置推定ステップ（ステップＳ２０）で得られた推定位置Ｌｎの周り（推定位置Ｌｎを含む所定の領域）で行われる。なお、テンプレートマッチングが適用される画像は、抽出画像Ｑｎ、時系列画像Ｆｎから新たに抽出される（切り出される）、推定位置Ｌｎを含む画像でもよい。何れの場合も、移動体Ｐ（Ｐｎ）とその周囲を示す画素値は同一なので、同一の結果が得られる。以下、説明の便宜上、テンプレートマッチングを適用する画像として時系列画像Ｆｎを用いるものとする。
【００３０】
テンプレートマッチングにおいて、テンプレート画像ＱＴと時系列画像Ｆｎの類似度は、例えば、数式５に示したゼロ平均正規化相互相関値ＲＺＮＣＣによって評価することができる。なお、数式５中における、テンプレート画像ＱＴの平均画素値は数式６によって得られ、時系列画像Ｆｎ（或いはテンプレートマッチングのために抽出された画像）の平均画素値は数式７によって得られる。
【数５】
2018230016.tif 000007
・・・（数式５）
【数６】
2018230016.tif 000008
・・・（数式６）
【数７】
2018230016.tif 000009
・・・（数式７）
ここで、Ｍｘ及びＭｙは、それぞれテンプレート画像ＱＴの幅及び高さを示す画素数である。Ｎｘ及びＮｙは、それぞれ時系列画像Ｆｎの幅及び高さを示す画素数である。（ｉ，ｊ）は各画像中の座標（ただし、０≦ｉ≦Ｍｘ−１、０≦ｊ≦Ｍｙ−１）である。Ｔ（ｉ，ｊ）はテンプレート画像ＱＴの画素値（強度、輝度値）である。Ｔ−（数式上はＴの上にマクロン（バー））はテンプレート画像ＱＴの平均画素値（平均強度、平均輝度値）である。Ｉ（ｉ，ｊ）は時系列画像Ｆｎの画素値（強度、輝度値）である。Ｉ−（数式上はＩの上にマクロン（バー））は時系列画像Ｆｎの平均画素値（平均強度、平均輝度値）である。
【００３１】
このテンプレートマッチングにより、時系列画像Ｆｎ内でテンプレート画像ＱＴと最も近似した移動体Ｐｎの位置を、時系列画像Ｆｎ内における移動体Ｐの正式な位置ＬＳｎ（ｎ＝１〜Ｎ）として特定する。換言すれば、時系列画像Ｆｎにおける移動体Ｐｎの推定位置Ｌｎを、テンプレートマッチングによって特定位置ＬＳｎに修正する。特定された移動体Ｐｎの位置ＬＳｎは、記憶部３２（図１参照）に保存される。
【００３２】
本実施形態では、テンプレート画像ＱＴがテンプレートマッチングの比較基準として用いられる。上述の通り、テンプレート画像ＱＴは、複数の時系列画像Ｆｎのうちの複数の抽出画像Ｑｎから所謂重ね合わせ法に基づいて生成される。従って、比較対象である移動体Ｐｎの形状が揺らぎ等の影響を受けて絶えず変化している場合でも、テンプレート画像ＱＴにおける移動体ＰＴの形状には、この変化が反映されにくくなる一方で、移動体Ｐｎに共通する形状が強く反映される。従って、テンプレート画像ＱＴはテンプレートとしての高い信頼性をもち、時系列画像Ｆｎにおける移動体Ｐｎの位置ＬＳｎを正確に特定することが可能になる。即ち、本実施形態によれば、重ね合わせ法を用いて移動体Ｐを精度良く検出することができる。
【００３３】
上述の通り、本実施形態によれば、時系列画像Ｆｎのそれぞれにおいて移動体Ｐｎの位置ＬＳｎが特定（修正）される。したがって、各画像における移動体Ｐの位置を繋ぎ合わせることにより、観測対象である移動体Ｐの軌跡や挙動を把握することができる。つまり、重ね合わせ法を用いた場合でも、観測時間内の移動体Ｐの位置について密な情報を得ることができる。
【００３４】
なお、テンプレート画像ＱＴの各画素値は、テンプレート画像ＱＴの生成に用いられる複数の抽出画像Ｑｎにおいて、対応する画素の強度（画素値）の中央値或いは最頻値でもよい。
【００３５】
また、本実施形態に係る移動体観測方法は、次の処理を含んでもよい。図５は本実施形態の第１変形例に係るフローチャートである。図５に示すように、移動体観測方法は、二値化ステップ（ステップＳ６０）と、重心算出ステップ（ステップＳ７０）とを含んでもよい。各ステップは、図１に示す移動体観測装置１０によって実行される。
【００３６】
二値化ステップ（ステップＳ６０）では、時系列画像Ｆｎ（ｎ＝１〜Ｎ）或いは抽出画像Ｑｎ（ｎ＝１〜Ｎ）中の各画素値が所定の閾値を用いて二値化される。閾値には、例えば観測履歴に由来するバックグランドノイズの値の所定倍の値などが適用できる。二値化ステップ（ステップＳ６０）では、二値化された抽出画像Ｑｎが得られる。そして、重心算出ステップ（ステップＳ７０）では、二値化された各時系列画像Ｆｎまたは各抽出画像Ｑｎにおける移動体Ｐｎの重心が算出される。
【００３７】
第１変形例におけるテンプレート生成ステップ（ステップＳ４０）は、重心算出ステップ（ステップＳ７０）で得られた複数の重心を、各抽出画像Ｑｎの重ね合わせにおける移動体Ｐｎの基準点として用いる。即ち、テンプレート生成ステップ（ステップＳ４０）では、重心が一致するように各抽出画像Ｑｎを重ね合わせることで、テンプレート画像ＱＴが生成される。その後、二値化された各画像を用いて位置特定ステップ（ステップＳ５０）が実行される。
【００３８】
図６は本実施形態の第２変形例に係るフローチャートである。図６に示すように、移動体観測方法は、第１変形例に係る二値化ステップ（ステップＳ６０）と、フィッティングステップ（ステップＳ８０）とを含んでもよい。第１変形例と同じく、第２変形例に係る各ステップも、図１に示す移動体観測装置１０によって実行される。
【００３９】
上述の通り、二値化ステップ（ステップＳ６０）では二値化された抽出画像Ｑｎが得られる。そして、フィッティングステップ（ステップＳ８０）では、二値化された各抽出画像Ｑｎにおける移動体Ｐｎの画素分布を所定の二方向に積算した２つの一次元分布に正規関数のフィッティングを行うことで移動体Ｐｎの最大輝度の位置（最大輝度位置）を算出する。
【００４０】
二値化ステップ（ステップＳ６０）によって、抽出画像Ｑｎの各画素は、移動体Ｐｎを示す画素値が１の画素と、移動体Ｐｎの背景を示す画素値が０の画素とに分類されている。そこで、フィッティングステップ（ステップＳ８０）では、これらの二次元分布に対して、Ｙ方向に積算したＸ方向の一次元分布を算出すると共に、Ｘ方向に積算したＹ方向の一次元分布を算出する。これらの２つの一次元分布に正規関数のフィッティングを行うことで移動体Ｐｎの最大輝度の位置（最大輝度位置）を算出する。
【００４１】
第２変形例におけるテンプレート生成ステップ（ステップＳ４０）は、フィッティングステップ（ステップＳ８０）で得られた複数の最大輝度位置を、各抽出画像Ｑｎの重ね合わせにおける移動体Ｐｎの基準点として用いる。即ち、テンプレート生成ステップ（ステップＳ４０）では、最大輝度位置が一致するように各抽出画像Ｑｎを重ね合わせることで、テンプレート画像ＱＴが生成される。その後、二値化された各画像を用いて位置特定ステップ（ステップＳ５０）が実行される。
【００４２】
なお、本開示は上述の実施形態に限定されず、特許請求の範囲の記載によって示され、さらに特許請求の範囲の記載と均等の意味および範囲内でのすべての変更を含む。
【図１】
2018230016.tif 000010
【図２】
2018230016.tif 000011
【図３】
2018230016.tif 000012
【図４】
2018230016.tif 000013
【図５】
2018230016.tif 000014
【図６】
2018230016.tif 000015
【手続補正書】
【提出日】平成30年5月21日(2018.5.21)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
所定の時間間隔で撮影された複数の時系列画像を取得し、
予め仮定した移動体の推定運動に基づいて、複数の前記時系列画像のそれぞれにおける前記移動体の推定位置を求め、
複数の前記時系列画像のうちの複数の画像から、前記移動体の前記推定位置を含む画像を抽出画像として抽出し、
前記移動体の基準点が一致するように複数の前記抽出画像を重ね合わせることによってテンプレート画像を生成し、
複数の前記時系列画像のうちの少なくとも１つの画像に対して前記テンプレート画像を用いたテンプレートマッチングを行うことで、複数の前記時系列画像のうちの前記少なくとも１つの画像における前記移動体の位置を特定する
ことを含む移動体観測方法。
【請求項２】
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記推定位置である、
請求項１に記載の移動体観測方法。
【請求項３】
前記テンプレート画像の各画素値は、複数の前記抽出画像において対応する画素の平均画素値、中央値及び最頻値のうちの何れかに設定される、
請求項２に記載の移動体観測方法。
【請求項４】
複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、
二値化された各前記抽出画像における前記移動体の重心を算出することを更に含み、
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、二値化された前記移動体の前記重心である、
請求項１に記載の移動体観測方法。
【請求項５】
複数の前記時系列画像または複数の前記抽出画像中の各画素値を二値化し、
二値化された各前記抽出画像における前記移動体の画素分布を所定の二方向に積算した２つの一次元分布に正規関数のフィッティングを行うことで前記移動体の最大輝度位置を算出する
ことを更に含み、
前記抽出画像の重ね合わせにおける前記移動体の前記基準点は、前記移動体の前記最大輝度位置である、
請求項１に記載の移動体観測方法。
【請求項６】
前記移動体は、地球軌道上のデブリ又は液体中の微生物である、
請求項１に記載の移動体観測方法。
【国際調査報告】
2018230016.tif 000016
2018230016.tif 000017
2018230016.tif 000018
2018230016.tif 000019