(57)【特許請求の範囲】
【請求項１】
複数コマの観測画像から前記観測画像上を移動していく移動体を検出する方法において、少なくとも３コマ以上の前記観測画像の中で前記移動体の動きを仮定し、その仮定された動きに合わせて前記観測画像の一部を切り取り、それら複数の切り取り画像において同じピクセル位置を占める各画素について特異値を除いた一定範囲内の画素値に基づいて求められた中央値によって中央値画像を作成し、仮定した動きの前記移動体の像のみを残存させた画像を得るものであって、
前記一定範囲は、前記切り取り画像における背景の中央値を中心とし前記背景の雑音レベルに基づいて定められる第１範囲内の前記画素値を持つ複数の前記画素から中央値と標準偏差とを求め、前記中央値を中心として前記標準偏差に基づいて定められる前記第１範囲よりも狭い幅狭範囲として、又は前記幅狭範囲内の前記画素値を持つ複数の前記画素に対して前記幅狭範囲を求める手法を適用する操作を繰り返して得られる一層狭い幅狭範囲として得られる第２範囲であることを特徴とする移動体検出方法。
【請求項２】
前記複数の中央値画像の平均値画像を作成し、仮定した動きの前記移動体の像のみを残存・強調させた画像を得ることから成る請求項１に記載の移動体検出方法。
【請求項３】
前記移動体の仮定される前記動きは、使用している観測システムで検出が可能な任意の移動方向及び速度を持つ動きであることを特徴とする請求項１に記載の移動体検出方法。
【請求項４】
前記観測領域は天球領域であり、前記移動体は人工天体、スペースデブリ、小惑星、彗星等の移動天体であることを特徴とする請求項１に記載の移動体検出方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、移動体（例えば、人工天体、スペースデブリ、小惑星、彗星等の移動天体）の検出方法に関し、特に、小型で暗い人工衛星等の人工天体の監視や軌道決定、また、運用中の人工衛星等に衝突した場合に産業に多大な被害をもたらすスペースデブリの発見や軌道決定、さらには、人類の産業活動や存続そのものに大きな影響をもたらす地球衝突型の小惑星や彗星の早期発見や軌道決定のための移動体検出方法に関する。
【０００２】
【従来の技術】
図８は、従来の移動天体の検出方法を説明する図である。図８中の観測画像（１−１）及び（１−２）は、同じ天球領域を一定時間間隔で、（１−１）、（１−２）の順で撮像した画像である。画像（１−１）中の１〜９、及び画像（１−２）中の１' 〜９' は恒星の画像であり、恒星１と恒星１' 、恒星２と恒星２' 、・・・、恒星９と恒星９' が、それぞれ同じ恒星を表している。画像（１−１）中の１０、及び画像（１−２）中の１０' は、移動天体を表しており、撮像間隔の間に天球上を移動していることがわかる。
【０００３】
従来の移動天体の検出方法は、この撮像して得られた２コマの画像を解析して、それぞれの画像に写っている天体を見付け出し、カタログ化する（図８中のカタログ（１−１Ｄ）、カタログ（１−２Ｄ）参照）。カタログには、それぞれ天体１〜１０、及び天体１' 〜１０' の位置、明るさ（Ｍａｇ）等の情報が書き込まれている。この２つのカタログを比較して、２コマの画像上で同じ位置に存在しない天体の対を検出することにより、移動天体１０（１０’）を知ることができる。
【０００４】
ところで、近年、ＣＣＤは大面積化し、一回の撮像で観測できる天空領域は飛躍的に広がっている。その結果、移動天体を捉えられる確率も上昇しているが、移動天体検出にとっては雑音となる恒星の数も増大し、その数は数万に達する。この膨大な雑音の中から移動天体を検出するには、かなりの困難を伴う。また、移動天体は画像上を移動してしまうため、暗い移動天体を検出するために、露出時間を長くして、天体からの光を画像上の定位置に蓄積させるという恒星の観測に用いられる技術を利用しても、充分な成果を得ることができない。有効な露出時間は、移動天体が点像に留まっている時間となり、限界等級は、観測に用いられる望遠鏡の口径と、ＣＣＤの量子効率で決まってしまう。
【０００５】
本出願人は、既に、複数コマの画像を利用することにより、利用した観測システムで検出可能な、あらゆる移動方向、移動速度の移動体について、その検出の妨げとなる大量の明るい光源の像を除去し、１コマの観測画像では検出が不可能な程度に暗い移動体の検出を可能にする方法を提案している（特２０００−３３４３６４）。この方法によれば、複数コマの画像から切り取った切り取り画像の各画素値の中央値を算出することで、恒星のような明るい光源の影響を除去すると共に背景雑音を抑えることで、一見、検出が不可能な程度に暗い移動体をも検出することが可能になった。
【０００６】
【発明が解決しようとする課題】
ところで、明るい光源（恒星）は、画像上、その中心部や周辺部である広がりを持っており、また暗い恒星の画像や画素値を返さない画素の欠損も存在しているので、それらの星々の像は切り取り画像の各画素値の中央値の算出に少なからず影響がある。そこで、それらの星々の影響を少なくして暗い移動体をより的確に且つ効率的に検出する点で、なお改善の余地がある。
【０００７】
この発明の目的は、こうした改善の余地に応えるため、複数コマの画像について各画素の中央値を用いることによって、利用した観測システム（例えば、望遠鏡及びＣＣＤカメラ）で検出可能な、あらゆる移動方向、移動速度の移動体を検出する移動体検出方法において、中央値の算出の際に移動体検出の妨げとなる大量の明るい光源像や暗い恒星又は画素値を返さない欠損した画素の影響を除去して、１コマの観測画像では検出不可能な程度に暗い移動体を効率的に検出することを可能する移動体検出方法を提供することである。
【０００８】
【課題を解決するための手段】
本発明は、上記課題を解決するため、複数コマの観測画像から前記観測画像上を移動していく移動体検出方法において、少なくとも３コマ以上の前記観測画像の中で前記移動体の動きを仮定し、その仮定された動きに合わせて前記観測画像の一部を切り取り、それら複数の切り取り画像において同じピクセル位置を占める各画素について特異値を除いた一定範囲内の画素値に基づいて求められた中央値によって中央値画像を作成し、仮定した動きの前記移動体の像のみを残存させた画像を得る手法を採用している。
【０００９】
本発明では、複数の観測画像からの移動体の動きに合わせた切り取り画像を用いて、移動体検出にとって雑音となって検出の妨げとなる大量の明るい光源（天体の場合には恒星の中心部やその周辺部を含む）の像や暗い画像又は画素値を返さない欠損した画素を中央値の算出から予め除去することで、特異値を示す画素値の影響を軽減した中央値画像を作成する手法を採用したことにより、１コマの観測画像では捉えることのできない暗い移動体の的確且つ効率的な検出が可能となること、また、こうした特異値を除く画像処理を施すことにより、移動体の検出限界等級を下げることが可能となると共に誤検出の割合も大幅に軽減することが可能となることをベースとしている。
【００１０】
本発明の移動体検出方法においては、前記一定範囲は、前記切り取り画像における背景の中央値を中心とし前記背景の雑音レベルに基づいて定められる第１範囲設定することにより、明るい星の中央部や画素値を返さない画素の影響を軽減することをベースとしている。
【００１１】
また、前記一定範囲は、前記第１範囲内の前記画素値を持つ複数の前記画素から中央値と標準偏差とを求め、前記中央値を中心として前記標準偏差に基づいて定められる前記第１範囲よりも狭い幅狭範囲として、又は前記幅狭範囲内の前記画素値を持つ複数の前記画素に対して前記幅狭範囲を求める手法を適用する操作を繰り返して得られる一層狭い幅狭範囲として得られる第２範囲とすることができる。第２範囲を設定することにより、明るい星の周辺部や暗い恒星の影響を軽減することができる。
【００１２】
この移動体検出方法において、前記複数の中央値画像の平均値画像を作成し、仮定した動きの前記移動体の像のみを残存・強調させた画像を得ることができる。即ち、移動体からの光量に対する雑音の比を大幅に抑えて、１コマの観測画像では検出が不可能であった暗い移動体を抽出することができる。
【００１３】
この移動体検出方法において、前記移動体の仮定される前記動きは、使用している観測システムで検出が可能な任意の移動方向及び速度を持つ動きであるとすることができ、そうすることにより、観測画像上を移動していく多数の暗い移動体を検出することができる。
【００１４】
この移動体検出方法において、前記観測領域を天球領域とすることで、前記移動体として人工天体、スペースデブリ、小惑星、彗星等の移動天体を検出することができる。天球領域の観測画像上には、明るい光源としての恒星の像が大量に写っているが、この検出方法を適用することにより、移動天体の検出には大きな妨害となる大量の恒星像を完全に除去し、移動天体の画像のみを残したり、或いは強調して残すことができる。
【００１５】
更に、多数の各移動天体に対して本発明の処理を実施することにより、着目天体だけでなく該着目天体に近い動きをしている他の天体の検出が可能であり、その画像も平均値処理を施すことで明るく強調することができ、これによって得られた多数の画像を総合することにより、いままでの検出方法では捉えることができなかった観測システムの検出限界等級以下の暗い移動天体を、自動的に捉えることができる。
【００１６】
このことは、運用中の人工衛星や現在建設中の国際宇宙ステーションに深刻な影響を及ぼし、さらに現在、検出が困難とされている径が数ｃｍ〜数１０ｃｍサイズのスペースデブリの発見、軌道決定に大きく貢献する。また、人類の活動や存続を左右する、径が数１００ｍ〜ｌｋｍサイズの地球衝突型の小惑星や彗星の早期発見を可能にする。
【００１７】
この移動体検出方法において、前記観測画像は、観測システムとして、近年特に発達が目ざましいＣＣＤカメラによって、同じ観測領域を撮像した画像とすることができる。
【００１８】
【発明の実施の形態】
以下、添付した図面に基づいて、この発明による移動体検出方法の実施例を説明する。図に示す実施例は、この発明で検出すべき移動体を移動天体とした例である。図１は本発明に係わる中央値画像の作成により恒星像が除去されることを説明するための基本概念図であり、図２は本発明に係わる平均値画像ができるまでの処理手順を示した図であり、図３は本発明において走査すべき移動天体の速度ベクトル領域を示した図であり、図４は移動天体の速度ベクトルの相違による切り取り画像の大きさ及び移動方向の違いを示す図であり、そして、図５は本発明による移動天体検出方法を適用して得られた画像の一実施例を示す図である。
【００１９】
本発明は、撮像された複数コマの観測画像を利用することにより、利用した観測システム（例えば、望遠鏡及びＣＣＤカメラ）で検出可能な、あらゆる移動方向、移動速度の移動体について、その検出の妨げとなる明るい光源が大量に存在していてもその光源の像を除去し、更に１コマの観測画像では検出が不可能な程度に暗い移動体の検出を可能にするものであるが、以下に示す実施例では、本発明による移動体の検出を、観測領域を天球領域として暗い移動天体を検出する場合に適用し、そのときの処理手順を説明する。この処理手順は、主に、次の３つに分けられる。
（１）あらかじめ移動天体の動作を予測して、複数コマの観測画像の中央値画像を作成し、移動天体の検出の妨げとなる恒星像を除去する、
（２）前記（１）の処理手順で作成された中央値画像を複数コマ作成し、それらの中央値画像の平均値画像を作成し、１コマの観測画像では検出が不可能である暗い移動天体をも抽出する、
（３）前記（１）及び（２）の処理手順を、さまざまな移動天体の動作に合わせて実行する。
【００２０】
以上の処理手順（１）〜（３）を、図面を参照して説明する。
上記処理手順（１）は、本発明の一実施形態に係り、複数コマの観測画像の中央値画像を作成し、移動天体の検出の妨げとなる恒星像を除去するものである。中央値画像の作成により恒星像が除去される処理手順（１）においては、図１に示されているように、まず、同一領域の観測画像である１−１，１−２，１−３が、時間軸に示すように、この番号順に、一定時間間隔で撮像される。各画像１−１〜１−３において、☆印は恒星を表し、●印は移動天体を表す。移動天体の動きは、予め予測したものである。
次に、一定時間間隔で撮像された複数コマの画像１−１〜１−３から、それぞれ移動天体の動きに一致した領域の画像（図中、点線で囲まれた領域）２−１，２−２，２−３が切り取られる。切り取り画像２−１〜２−３では、移動天体はすべての画像で同一位置に存在しているが恒星の位置は各画像で変化している。
【００２１】
次に、切り取った全ての画像２−１〜２−３の同一ピクセルについて、中央値を取ることによって、中央値画像３−１が作成される。平均値は全ピクセルの値を均等に考慮した値であるのに対し、中央値は、切り取った全ての画像の着目ピクセルについて取り込まれた値を値順に並べてその中央に位置するデータとして採用された値であり、特異値( 即ち、明るい恒星に対応した値) を示したピクセルの値とは関係しない値である。即ち、中央値を取ることにより、特異値の影響を少なくする処理( 即ち、恒星像の影響を少なくする処理) が行われると共に、背景雑音を大幅に抑えることができる（背景雑音は、使用した画像数の２分の１乗で小さくなる）。
切り取った画像のピクセル配置を表した図１のピクセル配置図において、１升目が、１つのピクセルを示す。例えば、各画像２−１〜２−３中の、ピクセル配置における位置Ｆにあるピクセルと、位置Ｍにあるピクセルとに注目する。簡単のため、天体が存在しているピクセルは、「１」というピクセル値を示し、存在していないピクセルは「０」というピクセル値を示すとする。位置Ｆに相当するピクセルは、切り取った画像２−１，２−２，２−３のうち、２−２に恒星が存在している。つまり、それぞれ位置Ｆのピクセル値は「０」、「１」、「０」であり、これらの中央値は「０」である。
【００２２】
また、位置Ｍに相当するピクセルは、全てに移動天体が存在している。つまり、位置Ｍのピクセル値は、「１」、「１」、「１」であり、これらの中央値は「１」である。中央値を取ることにより、切り取った画像上を移動していく恒星の影響を除去し、移動天体のみを残すことが可能になる。
本実施形態では、例として、中央値画像３−１を作成するために、３コマの観測画像１−１〜１−３を用いたが、利用するコマ数及び何コマの画像を利用するかは雑音となる恒星の数による。多くの恒星が画像に含まれている場合は、その影響をなくすため、より多くの画像を用いて中央値画像を作成する必要がある。
【００２３】
処理手順（２）は、本発明の次なるステップの実施形態に係り、処理手順（１）により得られた中央値画像を複数コマ用いて、平均値画像を作成し、１コマの観測画像では検出が不可能である暗い移動天体を抽出するものである。
図２は、平均値画像ができるまでの処理手順（２）を説明する図である。前記（１）の処理手順で、中央値を取ることにより、恒星の影響を少なくするという効果の他に、恒星以外の夜空の明るさによる背景雑音の影響を抑え、暗い移動天体を検出できるという効果がある。ここで、中央値を取ったときの背景雑音の値（中央値画像の背景雑音）をσｍとすると、σｍは、式（１）のように表される。
【数１】
0003671220.tif 000002
但し、σｉは、１コマの画像の背景雑音、係数１．２は統計的に得られた値、ｎは、中央値画像を作成するのに用いた観測画像のコマ数である。
【００２４】
しかし、平均値は、中央値と比較して、背景雑音の抑制がより効果的である。平均値をとったときの背景雑音の値（平均値画像の背景雑音）をσａとすると、σａは、式（２）のように表される。
【数２】
0003671220.tif 000003
上記式（１）、（２）からわかるように、平均値をとれば、中央値の１. ２倍暗い移動天体まで検出可能であることがわかる。
本実施形態は、さらに検出限界等級を下げるため、前記実施形態に係る処理手順（１）より、それぞれ複数コマの画像の中央値画像３−１，３−２，・・・，３−ｋを作成して恒星の影響を除去し、その後、より暗い移動天体を検出するために、作成された複数コマの中央値画像３−１〜３−ｋの平均値画像４−１を作成するものである。即ち、平均値画像４−１では、移動天体が存在しないところではたとえ一部の画像が明るくても平均的には暗くなるのに対して、移動天体が存在するところでは、常に存在する画像を平均することになるので、移動天体が存在しないところと比較して相対的に明るく強調されることになる。
【００２５】
本実施形態において暗い移動天体の検出を更に効率良く行うため、先の中央値画像を求める際の各画素についての中央値の計算においては、特異値が存在する明るい星の中央部やその周辺部、暗い恒星、及び画素の欠損に起因して値を返さない画素の影響を軽減することが好ましい。こうした悪影響を与える画素を可能な限り取り除く処理が、図６に示すフローチャートに従って行われる。この処理においては、先ず、切り取り画像の背景の中央値が求められる（ステップ１（「Ｓ１」と略す。以下、同様））。切り取られたすべての画像の同一画素（同じピクセル位置の画素）について、画素値ＩｉとＳ１で求めた背景の中央値Ｃ１との差の絶対値が、背景雑音のレベル（標準偏差σｍ）のある倍数よりも小さいか否かが判定され、その範囲（第１範囲Ｗ１）内の値を持つ画素のみを採用し、それ以外の値を持つ画素が無視される（Ｓ２）。σｍは、切り取りコマ数ｎに基づいて上記の式１から求められる。この倍数は、経験的にまた演算の簡素化のため、例えば、８倍程度とするのが好ましい。Ｓ１及びＳ２の処理によって、中央値の算出に際して、明るい星の中央部や画素の欠損に起因して値を返さない画素の値が予め取り除かれ、それらの異常な値の影響を軽減することができる。
【００２６】
更に、残された画素について中央値Ｃ２と標準偏差σとが計算され（Ｓ３）、画素値Ｉｊと中央値Ｃ２との差の絶対値が標準偏差σのある倍数（例えば、３倍、先の場合よりも小さくする必要がある）よりも小さいか否かが判定され、その範囲（第２範囲Ｗ２）内の値を持つ画素のみを採用し、それ以外の値を持つ画素が無視される（Ｓ４）。Ｓ３とＳ４の操作は、すべての画素値が、中央値を中心にして前後にそれぞれ標準偏差のある倍数の広がりを持つ第２範囲Ｗ２内に入るまで繰り返される。Ｓ３及びＳ４の処理を繰り返すことにより、中央値の算出に際して、明るい星の周辺部や暗い恒星の影響を軽減することができる。最後に得られた中央値が、中央値画像におけるそのピクセル位置の画素の代表値とされる（Ｓ５）。切り取り画像上のすべての画素について、上記の二段階による画素値の選択処理を行うことにより、悪影響を取り除いた中央値画像を得ることが可能になる。
【００２７】
この二段階による画素値の選択処理による画素値の取捨の概要が、図７に図示されている。図７は、例として切り取り画像数（ｎ）を１１とし、ある特定のピクセル位置にある画素の画素値を大きさの順に並べたものである。Ｓ１で中央値Ｃ１が特定（画素Ｉ６）され、Ｓ２で中央値Ｃ１を中心とする上下８σｍの第１範囲Ｗ１内の画素値Ｉｉ（Ｉ２〜Ｉ８）のみが採用される。この処理により、明るい星の中心部のような極端に明るい画素値や画素の値が返らない画素値が取り除かれる。次に、Ｓ３において、採用された画素値についての中央値Ｃ２（Ｉ５）と標準偏差σを求め、Ｓ４において中央値Ｃ２の上下３σの第２範囲Ｗ２内の画素値Ｉｉ（Ｉ３〜Ｉ５）のみが採用され、明るい星の周辺部や暗い恒星の影響が取り除かれる。Ｓ５において、このピクセル位置における中央値（この場合Ｃ２）として画素値Ｉ４が選択される。こうした特異値を除く画像処理を施すことにより、実際の画像解析において、移動天体の検出限界等級を更に下げることが可能となり、また、恒星等の影響による誤検出の割合が大幅に軽減された。
【００２８】
上記処理手順（３）は、本発明のさらに次のステップの実施形態に係り、前記処理手順（１）、（２）の作業を、移動天体の考えられるあらゆる速度ベクトルＶに関して実行するものである。即ち、上記の処理は、移動天体が存在するとして行われる手順であるので、移動天体の動きが捉えられない場合でも、移動天体の考えられるあらゆる速度ベクトルＶに関して画像の切り出しと中央値処理、或いは更に平均値処理を施すことによって、その切り出しのベクトル、即ち、切り取った画像の移動方向と距離が、たまたま移動天体の速度ベクトルＶ、即ち、移動天体の移動方向と距離に一致したときには、移動天体の存在を浮かび上がらせることができる。
図３は、処理手順（３）において実行される走査すべき移動天体の速度ベクトル領域を示した図である。太い矢印が移動天体の画像上の速度ベクトルＶであって、実線で囲まれた灰色部分がパラメータ領域であり、Ｘ軸及びＹ軸は移動天体の画像上の速度ベクトルＶの２つの成分を表している。速度成分の最大値は、ＣＣＤ画像の大きさと、利用する画像のコマ数によって決定する。
また、図３中で原点付近の領域は、移動量が小さいため、恒星の除去が効果的に行われない領域である。すなわち、移動天体が恒星の動きに近い動きをしている移動天体に着目して本発明の処理手順（１）を実行しても、その場合には恒星の除去は効果的でない。
図３のパラメータ領域にはいる任意の速度ベクトルＶに対し、切り取る画像の大きさは一意的に決定される。
【００２９】
図４は、移動天体の速度ベクトルの相違による切り取る画像の大きさ及び移動方向の違いを示す図である。実線で囲まれた四角は、ＣＣＤカメラによる観測画像の大きさをあらわし、破線で囲まれた四角は、切り取られた画像の大きさをあらわす。また、太い矢印は、速度ベクトルＶを示す。
図４（Ａ）で示されているように、速度ベクトルＶの絶対値が小さい場合、切り取る画像の面積は大きく、図４（Ｂ）で示されるように、速度ベクトルＶの絶対値が大きい場合、その面積は小さくなる。図４（Ａ）、図４（Ｂ）より、移動量の少ない移動天体の方が、検出が可能となる有効面積が広いことがわかる。
多数の各移動天体に対してこの処理を実施することにより、着目天体だけでなく該着目天体に近い動きをしている他の天体の検出が可能であり、その画像も平均値処理を施すことで明るく強調することができる。これによって得られた多数の画像を総合することにより、従来手法では得られなかった暗い移動天体の把握が可能となる。
【００３０】
実施例
図５は、本発明の移動天体検出方法を適用して得られた画像の一実施例を示す図である。
静止軌道上の衛星やデブリを捉えるため、口径５０ｃｍ、Ｆ２. ０の広視野望遠鏡に、３ｃｍ×３ｃｍのイメージエリアをもつＣＣＤカメラを設置し、静止軌道領域を恒星追尾モードで露出時間２秒、撮影間隔１３秒で、３０コマの連続撮像を行った。静止軌道上の天体の動作は、恒星追尾の場合、１秒間で約１５″東に移動する。そのため、この動きに合わせた領域を各画像から切り取り、まず、１０コマの中央値画像を３コマ作成し、その３コマの中央値画像の平均画像を作成した。
図５（Ａ）は、１コマの観測画像から静止軌道上の天体の動作に合わせて切り取った画像である。画像には、移動天体の検出には雑音となる恒星が多数写っている。図５（Ｂ）は、３０コマの観測画像を用いて本発明の手法によって作成された最終画像である。
画像のほぼ中央に、静止軌道上の天体が明確に捉えられている。雑音となる恒星は、ほぼ完璧に除去されており、また検出された静止軌道上の天体の輝度も、１コマの観測画像と比べて明るくなっていることがわかる。このことは、１コマの観測画像では捉えることのできない暗い移動天体も、本発明によれば不必要な明るい存在である恒星をまず除去し、必要ながら比較的暗い存在である移動天体を平均処理によって強調して検出できることを示している。
【００３１】
上記の移動体検出方法において、各撮影画像から直ちに画像の切り取りについて説明したが、それに先立って、観測によって得られた数１０枚〜数１００枚の画像について、ＣＣＤカメラ固有の電子雑音を除くため、すべての画像からダークフレームを差し引き、更にＣＣＤ画素間の感度の相違を補正するため、すべての画像をフラットフィールド画像で割るという前処理を施すのが好ましい。
また、上記の移動体検出方法において、画像の切り取り、中央値画像の作成及び平均値画像の作成による画像合成作業によって直ちに移動天体を検出したとしていたが、画像合成作業によって得られた移動天体はその候補であると位置付けて、検出された移動天体候補についてその移動量を微小に変化させて小領域の画像の切り取りを行い、その小領域の切り取り画像について、上記した画像合成作業を適用し、移動量を変化させたときの移動天体候補が示す最大値調査を行うのが好ましい。移動天体候補が示す明るさの値が最大となる移動量を真の移動量として、移動天体を特定することができる。特定された移動天体は、背景の星との位置関係及び輝度関係から、観測開始時及び終了時の正確な座標及び明るさを知ることができる。
【００３２】
以上のように、この発明による移動体検出方法を、添付した図面に基づいて移動天体の検出に適用した例で説明したが、様々な応用が可能である。即ち、運用中の人工衛星に衝突する危険性のあるスペースデブリの発見及び軌道決定、また将来的に増加すると思われるが現状の技術では光学的に検出が困難な静止軌道付近の小型衛星の軌道決定、人類の存続そのものに大きな影響をもたらす地球衝突型の小惑星や彗星の早期発見及び軌道決定等に応用することができる。更に、科学的分野への応用として太陽系の外縁の存在するエッジワースカイパーベルト天体の発見や小惑星帯に存在する数百メートルサイズの小惑星の発見、未知彗星の発見等が挙げられる。これら科学への応用は、太陽系やその生成についての新たな知見を得ることに大きく貢献するものと期待される。この発明による移動体検出方法は、その外にも、本発明の精神を逸脱しない限りにおいて、例えば、夜間において人や物の移動を検出するように、明るい光源が多数存在する中を移動する暗い移動体を検出するのにも適用することができること、及び細部を変更して実施することは、当業者にとって明らかである。
【００３３】
【発明の効果】
この発明による移動体検出方法は、上記のように複数コマの観測画像から前記観測画像上を移動していく移動体検出方法において、少なくとも３コマ以上の観測画像の中で仮定された移動体の動きに合わせて観測画像の一部を切り取り、それら複数の切り取り画像において同じピクセル位置を占める各画素について特異値を除いた一定範囲内の画素値に基づいて求められた中央値によって中央値画像を作成し、仮定した動きの前記移動体の像のみを残存させた画像を得ることから成っている構成されているので、移動体検出にとって雑音となる大量の明るい光源（天体の場合には恒星）の像や、暗い画像又は画素値をさない欠損画素の影響が効率的に除去され、１コマの観測画像では検出が不可能な程度に暗い移動体の効率的な検出をすることができる。本発明を暗い移動天体の検出方法として適用した場合には、移動天体の検出限界等級を更に下げることができ、一層暗い移動天体の検出が可能となると共に、恒星等の影響による誤検出の割合も大幅に軽減することができる。
【図面の簡単な説明】
【図１】本発明に係わる中央値画像の作成により恒星像が除去されることを説明するための基本概念図である。
【図２】本発明に係わる平均値画像ができるまでの処理手順を示す図である。
【図３】本発明において走査すべき移動天体の速度ベクトル領域を示した図である。
【図４】移動天体の速度ベクトルの相違による切り取り画像の大きさ及び移動方向の違いを示す図である。
【図５】本発明の移動天体検出方法を適用して得られた画像の一実施例を示す図である。
【図６】この発明による移動天体検出方法において、中央値画像を求める際の特異値を示す画素値を取り除く処理を示すフローチャートである。
【図７】この発明による移動天体検出方法において、第１範囲及び第２範囲による画素値の取捨選択の概要を示す図である。
【図８】従来の移動天体の検出方法を説明する図である。
【符号の説明】
１−１・・・・１−ｎ  観測画像          ☆    恒星
２−１・・・          切り取り画像      ●    移動体
３−１・・・          中央値画像        Ｖ    速度ベクトル
４−１・・・          平均値画像
Ｆ，Ｍ  ピクセルの位置          ｎ  切り取りコマ数
σｍ    中央値画像の背景雑音    σａ平均値画像の背景雑音
Ｉｉ    画素値                  Ｃ１，Ｃ２  中央値
Ｗ１  第１範囲                  Ｗ２  第２範囲
【図１】
0003671220.tif 000004
【図２】
0003671220.tif 000005
【図３】
0003671220.tif 000006
【図４】
0003671220.tif 000007
【図５】
0003671220.tif 000008
【図６】
0003671220.tif 000009
【図７】
0003671220.tif 000010
【図８】
0003671220.tif 000011