(57)【特許請求の範囲】
【請求項１】
地表面上を飛翔する飛翔体を上方の衛星から監視する飛翔体監視装置であって、
前記飛翔体の前記衛星からの観測方位を入力として該飛翔体の推定軌道を出力する軌道推定装置と、
飛翔体のハードウェアにより定まる情報からなる複数の飛翔体モデルを記録した飛翔体データベースと、
前記飛翔体データベースに記録された前記各飛翔体モデルに対して前記推定軌道で表される飛行条件に応じた予測光源スペクトルを求める光源スペクトル算出手段と、
前記各予測光源スペクトルに対して前記推定軌道で表される観測条件に応じた予測観測スペクトルを求める観測スペクトル算出手段と、
実測の観測スペクトルと前記予測観測スペクトルを比較してスペクトル誤差を計算し、前記飛翔体データベースの中で最も可能性の高い飛翔体モデルを選択する比較装置と、
前記スペクトル誤差に従い前記飛翔体データベースの中の前記飛翔体モデルの修正量を求めて修正を行う学習部と、
を備えたことを特徴とする飛翔体監視装置。
【請求項２】
前記飛翔体モデルが飛翔体の推薬の種類、エンジン仕様を含むことを特徴とする請求項１に記載の飛翔体監視装置。
【請求項３】
前記推定軌道が、飛翔体の飛行条件である高度と速度および観測条件である位置を含むことを特徴とする請求項１または２に記載の飛翔体監視装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、衛星に搭載した光学センサで地表付近を飛翔するロケット等の飛翔体を探知、類別する飛翔体監視装置等に関する。
【背景技術】
【０００２】
従来のこの種の飛翔体監視装置では、光学センサで観測したスペクトル(背景のみ、又は背景とロケット)と、テンプレートデータベース(各種ロケットのスペクトルと背景のスペクトルから成る)との相関を計算することにより、ロケットの存在判定(探知)やロケットの機種判定(類別)を行っていた(例えば下記特許文献１参照)。
また、大気中を飛翔するロケットの赤外線強度の高度との関係が下記非特許文献１に示されている。
なお、ロケットプルーム(排気)のスペクトルのことを略してロケットのスペクトルと記すことにする。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第２５４７９４８号明細書
【非特許文献】
【０００４】
【非特許文献１】F. S. Simmons 著、”Rocket Exhaust Plume Phenomenology”、The Aerospace Press 刊、ISBN 1-884989-08-X、p.26、2000年
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来の装置では、例えば、画像センサ、観測スペクトル出力手段、テンプレートデータベース、相関計算手段、から構成される。画像センサで撮影した画像を元に、観測スペクトル出力手段により画素毎の観測スペクトルが計算される。その観測スペクトルは、相関計算手段においてテンプレートデータベースに収められた複数のスペクトルとの相関が計算され、画素中にロケットが存在するか否か、また存在する場合にはさらにロケット機種も判定されて出力される。
【０００６】
このような飛翔体監視装置にあっては、テンプレートデータベースとして各種ロケットのスペクトルと背景のスペクトルを用いているが、これらは特定の条件で観測されたスペクトルに過ぎない。一般的に、ロケットの探知追尾には赤外線の波長帯が多用されるが、ロケットを特定機種に限定しても飛行条件(高度や速度)によりロケットから放射されるスペクトルが変化し、また観測条件(ロケット軌道で変化するロケットと衛星の位置関係に依存する大気透過率等)により衛星に到達する観測スペクトルは大きく変化する。
【０００７】
また、背景も地球表面の状況(海、表土、植生等)、雲やエアロゾルの状況(高度、粒径等)、太陽光の状況(昼夜、季節等)に依存して時間的空間的に変化するものである。従って、特定の条件で作成されたテンプレートデータベースとの相関を計算する従来の飛翔体監視装置では、実際に起こりうる種々の飛行条件や観測条件での観測スペクトルに対して、ロケットの存在判定(探知)やロケットの機種判定(類別)が十分に機能しないという問題点があった。
【０００８】
この発明は、上記の課題を解決するためになされたものであり、飛翔体の飛行条件や衛星からの飛翔体観測条件が様々に変化した場合でも、より可能性の高い飛翔体モデルを高い確度で選択することが可能な飛翔体監視装置等を提供することを目的とする。
【課題を解決するための手段】
【０００９】
この発明は、地表面上を飛翔する飛翔体を上方の衛星から監視する飛翔体監視装置であって、前記飛翔体の前記衛星からの観測方位を入力として該飛翔体の推定軌道を出力する軌道推定装置と、飛翔体のハードウェアにより定まる情報からなる複数の飛翔体モデルを記録した飛翔体データベースと、前記飛翔体データベースに記録された前記各飛翔体モデルに対して前記推定軌道で表される飛行条件に応じた予測光源スペクトルを求める光源スペクトル算出手段と、前記各予測光源スペクトルに対して前記推定軌道で表される観測条件に応じた予測観測スペクトルを求める観測スペクトル算出手段と、実測の観測スペクトルと前記予測観測スペクトルを比較してスペクトル誤差を計算し、前記飛翔体データベースの中で最も可能性の高い飛翔体モデルを選択する比較装置と、前記スペクトル誤差に従い前記飛翔体データベースの中の前記飛翔体モデルの修正量を求めて修正を行う学習部と、を備えたことを特徴とする飛翔体監視装置等にある。
【発明の効果】
【００１０】
この発明では、飛翔体の飛行条件や衛星からの飛翔体観測条件が様々に変化した場合でも、より可能性の高い飛翔体モデルを高い確度で選択することが可能である。
【図面の簡単な説明】
【００１１】
【図１】この発明の実施の形態１による飛翔体監視装置の構成(類別運用時の構成)を示すブロック図である。
【図２】この発明の飛翔体監視装置の運用概念を説明するための図である。
【図３】この発明の実施の形態２による飛翔体監視訓練装置の構成(訓練運用時の構成)を示すブロック図である。
【図４】この発明の実施の形態３による飛翔体監視装置用データベース更新装置の構成(データベース更新運用時の構成)を示すブロック図である。
【図５】この発明による飛翔体監視装置の変形例の構成を示すブロック図である。
【図６】この発明による飛翔体監視装置用データベース更新装置の変形例の構成を示すブロック図である。
【発明を実施するための形態】
【００１２】
この発明に係わる飛翔体監視装置は、飛翔体のテンプレートとして飛翔体の推薬やエンジン仕様を情報とする飛翔体モデルを格納したテンプレートデータベース(飛翔体データベース)、飛行条件に依存する光源スペクトルを計算する流体／燃焼計算装置(光源スペクトル算出手段)、観測条件に依存する観測スペクトル(大気減衰率)を推定する大気減衰推定装置(観測スペクトル算出手段)、上記飛行条件や観測条件である飛翔体の高度、速度、および位置を知るための軌道推定装置、を備える。
テンプレートデータベースとして飛翔体の推薬やエンジン仕様といったハードウェアで決まる固定的な情報を用いており、飛翔体の飛行条件や観測条件の影響を受けない。飛翔体の飛行条件、観測条件は軌道推定装置で計算し、飛行条件の影響は流体／燃焼計算装置で、観測条件の影響は大気減衰推定装置で、それぞれ計算し、観測スペクトルを得る。従って、実観測時の個別の飛行条件や観測条件の変化に影響されずに飛翔体の機種判定を高い確度で行うことができる。
【００１３】
またこの発明に係わる飛翔体監視訓練装置は、背景画像のスペクトルを記録した背景データベースまたは背景画像を得るための画像センサ、仮想軌道を飛行する飛翔体に対して計算された予測観測スペクトルと背景画像を合成する模擬画像生成装置、を備える。
任意の飛翔体軌道に対して観測スペクトルを計算することができるので、実背景やデータベースの背景と合成することでオペレータ訓練用の模擬画像を生成することができる。
【００１４】
またこの発明に係わる飛翔体監視装置用データベース更新装置は、データベースの更新のために、観測スペクトルの実測値と予測値のスペクトル誤差を計算する比較装置、スペクトル誤差を用いてテンプレートデータベースの中の飛翔体モデルを修正する学習部、を備え得る。
テンプレートデータベースに固定的な情報を用い、スペクトルの実測値と予測値の比較装置と学習部を備えているので、飛行条件や観測条件に影響されることなくテンプレートデータベースの情報を効率良く修正することができる。
なお、データベースの更新のために、上記飛翔体監視装置に学習部を設け、類別(監視)モードとデータ更新モードとで切り換えるように構成してもよい。
【００１５】
以下、この発明による装置を各実施の形態に従って図面を用いて説明する。なお、各実施の形態において、同一もしくは相当部分は同一符号で示し、重複する説明は省略する。また、以下の説明では飛翔体がロケットの場合を例に挙げて説明するが、この発明はこれに限定されない。
【００１６】
実施の形態１．
図１はこの発明の実施の形態１による飛翔体監視装置の構成(類別運用時の構成)を示すブロック図である。また図２はこの発明の飛翔体監視装置の運用概念を説明するための図である。例えば飛翔体監視装置を搭載した衛星Ａの画像センサ１では、大気Ｅに覆われた地球Ｄおよび目視線Ｆ方向の地表面上を飛翔するロケットＢのプルームＣが撮影される。背景除去装置２では、画像センサ１で得られた観測画像から背景が除去され、得られた観測画像はスペクトル計算装置３に送られる。スペクトル計算装置３では、背景を除去した観測画像から観測スペクトルが算出されて探知装置４と比較装置１０へ送られる。
【００１７】
探知装置４では、ロケットＢのプルームＣが探知された場合には存在画素情報を方位計算装置５へ出力する。方位計算装置５では、存在画素情報から衛星Ａから見たロケットＢのプルームＣへの観測方位が計算されて軌道推定装置６に送られる。軌道推定装置６では、方位からロケットＢの推定軌道(飛行条件：高度，速度、観測条件：位置)が計算される。
【００１８】
テンプレートデータベースであるロケットデータベース７には、飛翔体のハードウェアにより定まる固定的な情報からなる複数のロケットモデル情報(推薬の種類、エンジン仕様等を含む)が予め記録されている。光源スペクトル算出手段である流体／燃焼計算装置８では、ロケットデータベース７の複数のロケットモデル情報(飛翔体モデル)それぞれに対して、軌道推定装置６で計算されたロケットＢの推定軌道中の飛行条件である高度や速度を用いて光源スペクトルの予測値が計算される。観測スペクトル算出手段である大気減衰推定装置９では、流体／燃焼計算装置８で求めた予測光源スペクトルと軌道推定装置６で計算された推定軌道中の観測条件である位置を元に大気減衰率が計算され、さらにこの大気減衰率を用いて観測スペクトルの予測値が計算される。
【００１９】
比較装置１０では、スペクトル計算装置３で計算された観測スペクトルの実測値と、大気減衰推定装置９を通して計算された観測スペクトルの予測値とを比較して誤差スペクトルを計算し、ロケットデータベース７に記録された複数のロケットモデル情報の中で誤差スペクトルが小さく最も可能性の高い最尤ロケットモデルが選択される。
【００２０】
ここで軌道推定装置６における軌道推定計算では例えば、衛星Ａの画像センサ１等で観測されたロケットの方向(方位角と仰角等)の時刻歴を入力とし、仮定したロケット運動モデルを用いて最小二乗法等によりロケットの位置や速度等を推定する処理が行われる。
【００２１】
ロケットデータベース７のロケットモデル情報は、流体／燃焼計算装置８で必要となる、ロケット推薬の種類やロケットの性能パラメータ(特にノズル膨張比や燃焼室圧力：ノズル形状，エンジン仕様に係る)で構成される。
【００２２】
流体／燃焼計算装置８に関し、流体/燃焼計算は例えば、数値流体力学等で扱われる対象であり、与えられたロケット推薬やノズル形状に対してナビエ・ストークスの方程式等を解いてロケット排気の速度分布，温度分布，構成分子モル比率等を得る一般的な計算手法である。構成分子のモル比率とその温度が分かれば、光源であるロケット排気が放射する赤外放射スペクトルを計算することが可能となる。そして光源スペクトル算出手段は、ロケットデータベース７の類薬やノズルの性能パラメータを入力として、流体/燃焼計算により得られる流体としてのロケット排気の組成分布や温度分布、更にその結果から計算される赤外線分光放射強度(輝度)を求めることを意味する。なお流体計算に関して、ロケットの飛行高度や飛行速度が大きく影響する。
【００２３】
大気減衰推定装置９に関し、大気減衰計算は例えば、地表面上を飛翔する飛翔体を上方の衛星から監視する衛星システムでは、地表からの赤外放射を抑制する(大気透過率の低い)特定の波長帯を利用するため、ロケット光源の赤外放射を推定するには大気減衰の影響を補正する必要がある。赤外線の大気の分光透過率に関しては代表的なMODTRAN等の様々な推算式が存在する。ロケット排気スペクトル(光源スペクトル)と、ロケットＢと衛星Ａの位置関係(観測条件)で決まる大気減衰とから、衛星Ａの画像センサ１で受けるスペクトル(観測スペクトル)を計算することができる。そして観測スペクトル算出手段は、得られた赤外画像から光源としてのロケットＢのスペクトルに対し、衛星Ａからの観測条件に依存する大気減衰を考慮して観測した場合のスペクトルを計算する。比較装置１０で推算値と実測値を比較する条件を同一(衛星センサが受ける赤外放射スペクトル)にするために必要である。
【００２４】
以上の構成により、ロケットデータベース７に記録される複数ロケットモデル情報は、ロケットＢの飛行条件(高度や速度)および衛星ＡからロケットＢのプルームＣを観測する観測条件に依存しない情報になっている。その依存しない情報を元に流体／燃焼計算装置８で予測光源スペクトルを、更に大気減衰推定装置９で予測観測スペクトルを計算するので、最も可能性の高いロケットモデルを、飛行条件や観測条件や背景に依存すること無くロケットデータベース７から選択する類別機能を高い確度で実現することができる。
【００２５】
実施の形態２．
図３はこの発明の実施の形態２による飛翔体監視訓練装置の構成(訓練運用時の構成)を示すブロック図である。図３の(ａ)は第１の実施例の構成、(ｂ)は第２の実施例の構成、をそれぞれ示し、両者の違いは後述する、背景画像を画像センサ１から実測画像を得るか、背景データベース１１から選択するかの違いである。
【００２６】
図３の(ａ)の場合、衛星Ａに搭載された画像センサ１で地球Ｄを撮影し、背景画像を得る。訓練教官により入力部１５から入力された仮想軌道(飛行条件：高度，速度、観測条件：位置を含む)、さらにロケットデータベース７に記録された複数のロケットモデル情報(推薬の種類、エンジン仕様等を含む)の中から訓練教官により入力部１５からの入力により選択されたロケットモデル情報を元に、選択されたロケットモデル情報に関して仮想軌道の高度や速度を用いて流体／燃焼計算装置８で予測光源スペクトルが計算される。大気減衰推定装置９では、流体／燃焼計算装置８で求めた予測光源スペクトルと仮想軌道の位置を元に大気減衰率が計算され、さらにこの大気減衰率を用いて観測スペクトルの予測値が計算される。模擬画像生成装置１２では、背景画像に、仮想軌道の位置と予測観測スペクトルから計算されたロケットＢの観測画像を合成し、模擬観測画像を得る。なお、図３の(ａ)の場合、画像センサ１で撮影された背景画像を使用しているのに対し、図３の(ｂ)の場合、予め撮影して背景データベース１１に記録した背景画像を使用している。
【００２７】
ここで背景データベース１１は、ロケットが存在しない時の赤外画像である。リアルタイムに撮影された画像ではなく、背景となる陸(都市部、砂漠、森林等)や水面、昼夜、四季により変化する実測画像を予め複数保存しておき、オペレータ訓練用として例えば入力部１５からの入力指定された季節や時刻や場所の模擬画像を適宜生成するためのものである。
【００２８】
また模擬画像生成装置１２は、ロケットＢが存在しない背景画像の中に、ロケットデータベース７に登録されたロケットＢが仮想軌道を飛翔したと仮定したロケット像を加算(ロケットの赤外放射のエネルギー相当分画素値をロケットが存在する画素に加算)するものである。
【００２９】
この構成によれば、訓練教官が与えた任意の仮想軌道に対して模擬画像を生成することができるので、この発明で対象とする飛翔体監視装置のオペレータを訓練する訓練機能をもたせることができる。
【００３０】
実施の形態３．
図４はこの発明の実施の形態３による飛翔体監視装置用データベース更新装置の構成(データベース更新運用時の構成)を示すブロック図である。実施の形態１の類別機能と同様に最も可能性の高い最尤ロケットモデルが選択された後、予測観測スペクトルと画像センサ１での観測から得られた観測スペクトルとを比較装置１０で比較してスペクトル誤差を計算する。そしてスペクトル誤差を減少させるように学習部１３でロケットモデル修正量を計算し、ロケットデータベース７に記録されている最も可能性の高いロケットモデルの情報を修正する。
【００３１】
ここで学習部１３は、観測スペクトルの実測値と予測値の誤差を入力として、スペクトル誤差が小さくなるようにロケットデータベース７の中の該当ロケットモデルを記述する推薬(例えば燃料と酸化剤の比率)やロケットノズル膨張比を修正するものである。例えば燃料と酸化剤の比率を変化させて予測観測スペクトルを計算し、実測値との誤差が小さくなる方向に燃料と酸化剤の比率を変更するものである。
【００３２】
この構成によれば、ロケットＢの飛行条件や衛星ＡとロケットＢの観測条件に依存せずにロケットデータベース７の中の選択されたロケットモデル情報を修正することができるので、さまざまな飛行条件や観測条件での観測結果を使って選択されたロケットモデル情報を効率よく修正するデータベース更新機能(学習機能)をもたせることができる。
【００３３】
なお、図１の飛翔体監視装置において、実施の形態１で説明した類別モードと上記説明によるデータベース更新モード(学習モード)を切り換えて実行する構成にしてもよく、その場合の構成は図１に破線で示す構成となる。類別モードではロケットデータベース７に記録された複数のロケットモデル情報の中で誤差スペクトルが小さく最も可能性の高い最尤ロケットモデルが選択され、データベース更新モードでは、スペクトル誤差を減少させるように学習部１３でロケットモデル修正量を計算し、ロケットデータベース７に記録されている最も可能性の高いロケットモデルの情報を修正する。これらはそれぞれ別々の独立したモードとして処理される。
【００３４】
また、上述の実施の形態１，３に関する図１、図４の構成では、背景除去装置２の出力する背景除去された観測画像に対してスペクトル計算装置３で観測スペクトルを求めていたが、それぞれ図５，図６に示すように探知装置４のロケットＢのプルームＣが探知された存在画素情報に対してスペクトル計算装置３で観測スペクトルを求めるようにしてもよい。これにより処理の効率が向上する。
【００３５】
また、上記各実施の形態の構成において、画像センサ１以外の構成については、衛星搭載計算機等の処理能力に応じて一部または全部を地上設置とし、地上処理に必要な情報は衛星から地上に通信手段(図示省略)で転送するようにしてもよい。
【００３６】
なお、この発明は上記各実施の形態に限定されるものではなく、各実施の形態の特徴の可能な組み合わせを全て含むことは云うまでもない。
【符号の説明】
【００３７】
１  画像センサ、２  背景除去装置、３  スペクトル計算装置、４  探知装置、５  方位計算装置、６  軌道推定装置、７  ロケットデータベース(飛翔体データベース)、８  光源スペクトル算出手段(流体／燃焼計算装置)、９  観測スペクトル算出手段(大気減衰推定装置)、１０  比較装置、１１  背景データベース、１２  模擬画像生成装置、１３  学習部、１５  入力部、Ａ  衛星、Ｂ  ロケット、Ｃ  プルーム、Ｄ  地球、Ｅ  大気、Ｆ  目視線。
【図１】
0005748532.tif 000002
【図２】
0005748532.tif 000003
【図３】
0005748532.tif 000004
【図４】
0005748532.tif 000005
【図５】
0005748532.tif 000006
【図６】
0005748532.tif 000007