JPH1054879A - 回転ストリップ開口映像測定を使用する運動物体および過渡事象検出 - Google Patents
回転ストリップ開口映像測定を使用する運動物体および過渡事象検出Info
- Publication number
- JPH1054879A JPH1054879A JP9117620A JP11762097A JPH1054879A JP H1054879 A JPH1054879 A JP H1054879A JP 9117620 A JP9117620 A JP 9117620A JP 11762097 A JP11762097 A JP 11762097A JP H1054879 A JPH1054879 A JP H1054879A
- Authority
- JP
- Japan
- Prior art keywords
- aperture
- noise
- image
- target
- difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G3/00—Observing or tracking cosmonautic vehicles
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/262—Analysis of motion using transform domain methods, e.g. Fourier domain methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/50—Systems of measurement based on relative movement of target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
- G01S17/90—Lidar systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Physics (AREA)
- Astronomy & Astrophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Radiation Pyrometers (AREA)
- Measurement Of Radiation (AREA)
- Studio Devices (AREA)
- Image Processing (AREA)
Abstract
回転ストリップ開口を使用して良好な解像度でセンサの
視野内の運動目標の検出を可能にするセンサシステムを
提供することを目的とする。 【解決手段】 光軸を中心に回動し、時間的に順次連続
する情景の映像を生成する回転ストリップ開口を具備す
る望遠鏡11と、この望遠鏡11の焦点平面に形成された映
像を検出するための2次元の検出器アレイ17と、この検
出器アレイ17に結合され、ストリップ開口が望遠鏡11の
光軸を中心に回動するとき映像情景の複数の映像フレー
ムを記録し、それらの複数の映像フレームからの映像を
合成し、合成された映像中に発見された運動目標を検出
する信号プロセッサ18とを備えていることを特徴とする
Description
定装置および放射線検出方法に関し、特に、合成された
等価全円形開口が映像を生成し、複数の回転ストリップ
開口映像からの運動目標の検出または合成された映像の
測定を行う回転ストリップ(部分的)開口映像生成用放
射線測定装置に関する。
高い解像度の映像を得るために、例えば通常のセンサ技
術では大型で、重く、展開可能な光学システムの能動制
御装置を使用している。要求される使命および1次反射
器の大きさに応じて、センサの全ての光学素子の位置、
傾斜、形状、および整列制御のために1次反射器の周期
的な位置および傾斜制御から配列することができる。本
発明と同じ解像度を有する全開口システムは比較的大き
い開口面積のために多量の集光能力を有している。しか
しながら、軌道上に数十メートルの直径の開口を配置す
るためには、本発明と競合する全開口システムは次のよ
うなことが必要になる。すなわち光学システムの直径が
発射装置よりも大きくなるとセグメント化された光学表
面および折畳まれた支持構造が必要になる。また、重量
を節減するために薄い変形可能な反射鏡が使用される場
合には複雑で潜在的に高い帯域幅の適応光学技術が必要
であり、さらに、位相アレイとして構成される場合には
複雑な位置および瞳孔整合制御が要求される。それ故、
全開口システムは本発明と比較して比較的重く、高い技
術的なリスクを有している。
スピンApと呼ばれる最も近い従来技術のシステム能力
を改善させるものである。このスピンApシステムおよ
び方法は米国特許第5243351号明細書に記載され
ている。このスピンApシステムは全円形開口に等価な
高い解像度を有する映像を合成することのできるデータ
処理方法である。等価なスピンAp合成全円形開口はス
トリップ開口形状と関連する最大相関長に等しい直径を
有する。等価全円形開口映像を正確に合成するために、
スピンApは、等価全円形開口の光学システム通過帯域
を完全に測定するために適切な時間インターバルで回転
するセンサの焦点平面検出器を使用する映像データを捕
捉する。映像は米国特許第5243351号明細書に記
載された方法により処理される。
開口運動目標検出と呼ぶ)との共通点は、全円形開口に
等価な解像度を有する映像を合成するために時間的に整
合したストリップ開口を使用することである。等価回転
開口の合成された全円形開口はストリップ開口の形状に
関係した最大の相関長に等しい直径を有する。等価全円
形開口映像を正確に合成するために、スピンApは適切
な時間インターバルで回転するセンサの焦点平面検出器
を使用する映像データを捕捉、し等価全円形開口の光学
システム通過帯域を完全に測定する。映像は米国特許第
5243351号明細書に記載された方法によって処理
される。
のの瞬間的視野中の運動目標は全体の合成された映像中
に運動による軌跡の目印を生じる。運動目標の検出は、
位置、寸法、軌跡、速度、加速度、および原点のような
目標現象を設定するために必要である。さらに、検出プ
ロセスは運動目標の静止映像合成を可能にする捕捉され
たデータにプロセッサおよび方法を適用するための前段
階である。それ故、運動目標検出プロセッサおよび方法
の階級について説明する。この階級は、空間的、時間
的、空間周波数の、時間周波数ドメイン検出プロセッサ
を含み、それはまた多重スペクトル背景排除技術を伴っ
ている。本発明のシステムは、スピンApの個々のフレ
ーム中の任意の1つにおけるスピンApの任意の対の視
野中の合成された映像の視野中の運動目標を検出するた
めに利用できる測定を使用する。同じ方法はまた過渡事
象の検出にも適用できる。
は、所定の実効合成開口寸法に対して低いペイロード重
量が得られ、一方センサの視野内の運動目標の検出を可
能にする変形されたスピンAp映像処理方法を提供する
ことである。さらに本発明の別の目的は、宇宙空間で使
用するための複雑な適応光学システムを必要とする複雑
で大型の全円形開口(または位相アレイ望遠鏡)を展開
するのに代えて、低いリスク、低いコスト、軽い重量、
および簡単な展開構造の映像処理および運動目標検出方
法を提供することである。
目的を達成するために、本発明の1形態においては、運
動目標を検出することのできる回転ストリップ放射線映
像形成および検出システムが提供される。このシステム
は時間的に連続する、または順番に並んだ映像を生成す
る回転ストリップ開口望遠鏡を備えている。回転ストリ
ップ開口望遠鏡は典型的にストリップ開口1次反射器お
よび2次反射器を備えている。2次元の検出器アレイが
設けられて望遠鏡の焦点平面に位置する映像を検出す
る。回転補償装置が使用されてアレイの検出器の積分時
間中に回転の軌跡が描かれることを阻止している。信号
プロセッサが設けられてストリップ開口が望遠鏡の光学
軸を中心に回転するとき望遠鏡によって映像が形成され
た目標情景の複数のフレームを記録し、等価全円形開口
映像を合成し、個々のストリップ開口フレーム、任意の
対のフレーム、合成された映像、または個々のフレーム
と合成された映像の組合わせにおける運動目標を検出す
る。
像を合成し、複数の回転ストリップ開口映像測定からロ
センサの視野中の運動目標を検出し、一方、個々のスト
リップ開口映像間のランダムおよび/または系統的な照
準線エラーを補償する。したがって、本発明は等価全円
形開口映像を合成するための手段を提供することによっ
て米国特許第5243351号明細書に記載されたスピ
ンApシステムに比較して改良された全体的なシステム
性能を提供し、同時に運動目標を検出する能力を与える
ことができる。
いて以下説明する。放射線測定システム、映像合成、お
よび運動目標検出処理方法の各実施形態は特定の使命の
要求に依存し、また放射線測定システム、映像合成、お
よび運動目標検出方法の技術的妥協に依存し、フレーム
間の照準線のランダムまたは合成によるドリフトの補
償、共通平面にない点の広がり機能、フィールド点と関
係のない映像の運動による映像の軌跡と、検出器の有限
の電子的帯域幅、フィールド点に依存する映像の回転運
動が補償されないことによる軌跡の補償等も必要に応じ
て選択される。
よび検出プロセスの実施形態について以下説明する。回
転ストリップ開口が望遠鏡の光学軸を中心に回転すると
き次のような動作が生じる。回転補償装置は検出器の積
分時間中逆回転し、それによって時間的に静止した映像
を与える。映像フレームは記録され保存される。回転
(情景に対して)検出器アレイ構造が選択される場合、
得られたフレームは合成された映像の基準グリッドに座
標変換され、補間される。データはフーリエ変換されて
記憶される。新しいフレームが記録され、保存される。
間での照準線エラーによる不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換はそれらの照準線
エラーに対して補正され、記録される。先行段階は各ス
トリップ開口映像フレームに対して順次反復され、また
はフレームは順次捕捉され、記憶されて、その後照準線
エラーの全体的な評価が決定され、フレームの整合に使
用される。個々のストリップ開口映像の任意の対または
対の組合わせが照準線エラーに対して補正されると、回
転開口運動目標検出プロセスが開始される。さらに、一
度十分な数(前記米国特許第5243351号明細書に
記載されたような)のストリップ開口映像測定が得られ
ると、映像合成処理が開始される。回転開口運動目標検
出方法は2回以上合成映像について行われる。
た2以上の合成映像を使用する回転開口運動目標検出方
法の実施形態は合成された映像を空間的に整合させ(測
定セット間の照準線エラーにより要求される場合)、合
成された映像のすべての可能な異なる対の差が計算さ
れ、差の映像が合計され、整合されたフィルタを通って
差および合計データが送られ、比較器で比較される。整
合されたフィルタからの出力が比較器の予め定められた
しきい値を越えたとき運動目標が検出されたことが決定
される。
2以上の個々の照準線の補正されたフレームを使用する
回転開口運動目標検出方法の別の実施形態は、システム
伝達関数(STF)等化フィルタを差を得る個々のフレ
ームに適用して等化されたフレームの任意のまたは全て
の対の差を生成し、得られた差のデータに対して共通で
ない通過帯域排除フィルタを適用する。濾波された差の
データは濾波された差の全体の合計されたスペクトルを
形成するために加算により結合される。加算されたスペ
クトルは逆変換され、それによって空間的ドメインにお
いて運動目標の目印を与え、それは雑音を最適にする整
合されたフィルタアレイおよび比較器に供給される。整
合されたフィルタアレイからの出力が比較器の予め定め
られたしきい値を越えたとき運動目標が検出されたこと
が決定される。
ルタ検出プロセッサは、形状、寸法、速度、加速度、背
景と目標のコントラスト比、信号対雑音比、目標の数、
軌道交差の数、等の潜在的な目標の物理的現象に対応す
る特性を有している。整合されたフィルタ検出プロセッ
サはまた等化、差の形成、濾波、ならびに回転開口セン
サシステムの特有の信号測定特性により伝播された測定
雑音統計の過去の知識を有している。
作は、個々のフレーム映像の共通の重なった通過帯域に
おける全ての時間的に静止している成分を除去し、また
合成映像は軌跡を描く運動目標の雑音を含む目印のみを
残す。この運動目標検出技術の性能はシナリオに依存
し、等化、差の形成、濾波動作を行われる個々のフレー
ム測定雑音、および残留フレーム対フレームによる背景
漏洩および、または合成された映像整合エラーに対する
合成された映像による背景漏洩によって劣化される。
法の別の実施形態は、任意の差を形成されたフレームに
対して最適の整合フィルタを適用し、全ての可能なフレ
ーム差に対して出力を組合わせることに対応する。
の典型であるが、時間周波数ドメインプロセスもまた使
用可能である。時間周波数ドメインプロセッサの1実施
形態は同じ情景を観察することによって得られた空間的
に整合された個々の回転ストリップ開口フレームの時間
的シーケンスを利用する。空間的に整合されたフレーム
は共通空間的周波数グリッドへ空間的にフーリエ変換さ
れる。共通空間的周波数グリッドは等価全円形開口回転
ストリップ開口合成空間周波数通過帯域と関係する。任
意のまたは全ての空間的に整合され、フーリエ変換され
た個々のフレーム映像は各空間周波数を有する全ての時
間的周波数と関連するデータセットを生成するために時
間的なフーリエ変換をされる。最良の可能な検出特性を
得るために、時間的および空間的に周波数変換されたデ
ータセットは雑音最適化された整合されたフィルタ検出
プロセッサに転送される。整合されたフィルタ検出プロ
セッサは目標現象、信号対雑音測定の統計、ならびに回
転ストリップ開口センサシステムのユニークな測定特性
についての過去の知識を有している。整合されたフィル
タアレイからの出力が予め定められたしきい値を越えた
とき運動目標が検出されたことが決定される。
整合されたフィルタは個々のフレームの変換された時間
ヒストリーにおける全ての時間に静止している成分を消
去し、それによって運動目標の時間的周波数および空間
的周波数目印のみを残す。運動目標検出技術の性能はシ
ナリオに依存し、個々のフレーム測定雑音伝播およびフ
レーム整合エラーに対する残留フレームによる背景漏洩
によって劣化される。観察パラメータが許容すれば、許
容可能な性能は最適の方法より少ないものに対して達成
され、それによって最適に整合されたフィルタアレイお
よび比較器はテンプレート整合されたフィルタおよび比
較器、ストリーク(筋)検出器および比較器、または単
なる単一画素または多重画素しきい値超過検出プロセッ
サによって置換されることができる。これにおいて画素
と言う用語は、映像の2次元のディスクリートな表示の
空間的、または空間周波数素子のことである。
よび方法の階級が与えられる。階級は空間、時間、空間
周波数、および時間周波数ドメイン検出プロセッサなら
びに多重スペクトル背景排除技術を含んでいる。適切な
処理手順および方法の選択はシナリオおよび運動目標の
実効的信号対雑音比特性に依存している。同じ方法は過
渡事象の検出にも適用することができる。
識、アレイの検出器の空間的応答特性、雑音統計、記録
された各ストリップ開口映像の時間的整合はこのセンサ
システムおよび処理方法による運動目標検出を可能にす
る。
を参照することにより本発明の種々の特徴および利点が
容易に理解されるであろう。おな、同じ参照符号は同じ
構造素子を示すものである。図面を参照すると、図1は
本発明の原理を使用して適切に修正された回転開口映像
生成放射線センサシステム10およびデータ処理方法20
(米国特許第 5,243,351号明細書に記載されたシステム
に対応している)を示す。回転開口映像生成放射線セン
サシステム10およびデータ処理方法20は、運動目標の検
出情報を提供し、等価全円形映像を生成すると同時に照
準線ジッタおよびフレーム整合エラーを除去する能力を
提供する。したがって、ここに記載されたシステム10お
よびデータ処理方法20は、米国特許第 5,243,351号明細
書の回転開口(SpinAP)システムを改良したものである。
米国特許第 5,243,351号明細書の内容全体がここにおい
て参考文献に含まれている。
は、1次反射器12aおよび2次反射器12bを含む回転ス
トリップ開口望遠鏡11を含んでいる。3次反射器(示さ
れていない)は、ある環境下で望遠鏡11において使用さ
れてもよい。システム10は、地球指向アンテナ13aを有
する静止セクション13を含む衛星の形態で示されてい
る。望遠鏡11はプラットフォーム14上に取付けられ、静
止セクション13aもまたそれに結合されている。回転開
口映像生成放射線センサシステム10は、1次反射器12a
(ストリップ開口12を含む)が望遠鏡の光学軸を中心と
して回転したときに、望遠鏡11によって映像化されたタ
ーゲット情景の多数の映像フレームを整合するように構
成されている。
装置16は、検出器アレイ17によるその感知の前に映像を
安定させ回転を除去するように構成され、望遠鏡11の光
路に沿って配置されている。回転除去装置16は、回転補
償制御装置19の制御の下に検出器アレイ17を含む検出器
の積分時間中映像を逆回転させ、それによって静止映像
を生成する。照準線安定化ミラー15は、広い帯域の照準
線エラーおよびシステム10の軌道ダイナミックスによる
照準線エラーを除去するために照準線制御システム(信
号プロセッサ18またはその他の専用制御システムにより
提供されてもよい)によって使用される。
るように構成された信号プロセッサ18に結合された望遠
鏡11の焦点平面に配置された検出器アレイ17上に投影さ
れ映像化される。個々の映像フレームは、等価全円形開
口映像を合成するために信号プロセッサ18において処理
され結合される。この強化された運動目標検出回転開口
システム10において、信号プロセッサ18(または独立し
たプロセッサ)は運動目標を検出するために本発明の方
法にしたがって修正される。以下、図2乃至5を参照し
て修正された処理方法を詳細に説明する。
テムにより使用される本発明の原理によるセンサシステ
ム10および処理方法20,20-2,20-3,20-4,の種々の異
なる実施形態を示す。観測シナリオおよびタイムライ
ン、利用可能な過去の知識の量、および処理チェインの
利用可能な計算処理に基づいて、特定の実施形態が選択
されている。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による簡単
な回転開口運動目標検出処理方法20のトップレベルの図
を示す。図2に示された検出方法20は、米国特許第 5,2
43,351号明細書にしたがって情景の合成された等価全円
形開口映像30を生成するために同じ情景を観察する2以
上の回転開口ストリップ開口映像セットを使用する。照
準線整合エラーが存在する場合、合成された映像30のそ
れぞれは米国特許第 5,243,351号明細書に記載された方
法によって空間的に最適に整合される。
くつかまたは全ての可能性のある別個の対は、差を計算
され40、続いて加算されて35、記憶され、それによって
差を計算された合成された映像の総和を生成する45。別
個の合成された映像30の数がNである場合、可能性のあ
る別個の差を計算された合成された映像対の総数は、N
(N−1)/2によって与えられる。
差を計算された合成された映像の結果的な和45は、雑音
の最適整合されたフィルタ検出プロセッサ50に送られ
る。整合されたフィルタ検出プロセッサ50は、形状、寸
法、速度、加速度、背景と目標のコントラスト比、信号
対雑音比、目標数、軌道交差数等のような潜在的な目標
の物理的現象に対応した特性を有する整合されたフィル
タ55のアレイを含んでいる。整合されたフィルタ検出プ
ロセッサ50はまた合成および差計算処理により伝播され
た測定雑音統計の過去の知識、および回転開口センサシ
ステム10の特有の信号測定特性に関する知識を含んでい
る。検出プロセッサ50の整合されたフィルタ55のそれぞ
れは、差を計算された合成された映像の総和45のあらゆ
る画素位置に適用される。整合されたフィルタ55のそれ
ぞれの出力は比較器60に伝送され、しきい値が各整合さ
れたフィルタ55の予め定められたレベルに設定される。
運動目標の検出は、整合されたフィルタ55からの出力が
比較器60の予め定められたしきい値を越えたときに決定
される。しきい値は、信号および雑音統計に関する演繹
的な想定によって、或は個々のフレームまたは回転開口
フレームのセットに基づく信号および雑音統計に関する
帰納的な想定によって予め決定されていてもよい。
ての時間的に静止した成分を消去し、合成された運動目
標の軌跡の目印だけを残す。運動目標検出処理方法20の
特性はシナリオ依存性であり、合成および差計算動作に
よって伝播された個々のフレーム測定雑音、および残留
フレーム対フレームの整合エラーによる背景漏洩によっ
て劣化される。観察パラメータが許容すれば、許容可能
な性能は最適以下の方法に対して実現可能であり、それ
によって最適な整合されたフィルタ55および比較器60は
テンプレート整合されたフィルタおよび比較器60、スト
リーク検出器および比較器60、或は単なる単一画素また
は多重画素しきい値超過検出プロセッサによって置換さ
れてもよい。
43,351号明細書に記載された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による回転
開口運動目標検出処理方法20-2のトップレベルの実施形
態を示す。図3に示された検出方法20-2は、同じ情景を
観察することによって得られ、米国特許第 5,243,351号
明細書に記載された方法を使用して生成された2以上の
回転開口ストリップ開口映像30-2を使用する。照準線整
合エラーが存在する場合、個々の回転開口映像30-2のそ
れぞれは、米国特許第 5,243,351号明細書に記載された
ように空間的に最適に整合され、共通の空間周波数座標
系およびグリッドにフーリエ変換される。典型的に、共
通の空間周波数グリッドは、等価全円形開口回転開口合
成された空間周波数通過帯域と関連している。
のフレーム映像30-2のいくつかまたは全ての可能性のあ
る異なる対は、個々のフレームのそれぞれと関連し、か
つ回転開口センサシステム10の回転の結果生じた時間的
に変化するシステム伝達関数を補償するように設計され
た空間周波数フィルタ31-2によって乗算される。システ
ム伝達関数の補償は、ターゲット全円形開口またはその
他の基準システム伝達関数に各フレームを等しくする
か、差を計算されるべきフレームのシステム伝達関数に
1つのフレームを等しくするか、或は差を計算された両
映像からシステム伝達関数をデコンボルブすることを含
むいくつかの方法の1つによって行われてもよい。理想
的なシステム伝達関数に対して、またπラジアン離れて
捕捉されたフレームだけが差を計算される場合にのみ、
システム伝達関数の等化は不要である。
の可能な等化過程の任意の1つに基づいた回転開口運動
目標検出方法(20-2)は、特定の方法の受信機動作特性
(ROC)によって特徴付けられる特有の雑音特性を有
している。信号テンプレート方法に対するROCの差
は、等化プロセスによって発生された異なる補償されな
いレベルの雑音“ブースト”または“デブースト”によ
って支配される。しかしながら、適切に設計された雑音
に対して最適整合されたフィルタ方法は、等化プロセス
において雑音特性を補償する。
れ(40-2)、2つのフレームと関連した重畳した通過帯
域の外側の非共通情報を除去するために空間周波数濾波
される(41-2)。空間周波数濾波されたフレーム差は、
前に濾波されたフレーム差に加算され(35-2)、記憶さ
れ、それによって差を計算された空間周波数濾波された
フーリエ変換された映像の総和を生成する(35-2)。異
なった個々のフレーム映像の数がNである場合、可能性
のある異なった差を計算された個々のフレーム映像対の
総数は、N(N−1)/2の大きさである。
差を計算された空間周波数濾波され変換された映像の結
果的な和35-2は、空間ドメイン映像45-2を得るために逆
フーリエ変換され(42-2)、この空間ドメイン映像は雑
音最適整合されたフィルタ検出プロセッサ50-2に伝達さ
れる。整合されたフィルタ検出プロセッサ50-2は、形
状、寸法、速度、加速度、背景と目標のコントラスト
比、信号対雑音比、目標の数、軌道交差の数等のような
潜在的な目標の物理的現象に対応した特性を有する整合
されたフィルタ55-2のアレイを含む。整合されたフィル
タ検出プロセッサ50-2はまた等化、差計算および濾波処
理により伝播された測定雑音統計の過去の知識、および
回転開口センサシステム10の特有の信号測定特性を含ん
でいる。検出プロセッサ50-2の整合された各フィルタ55
-2は、差を計算された濾波された映像の総和35-2の全て
の画素位置に適用される。整合されたフィルタ55-2のそ
れぞれの出力は比較器60-2に伝送され、しきい値が各整
合されたフィルタ55-2に対して予め定められたレベルに
設定される。運動目標の検出は、整合されたフィルタ55
-2のアレイからの任意の出力が比較器60-2の予め定めら
れたしきい値を越えたときに決定される。しきい値は、
信号および雑音統計に関する演繹的な仮定によって、或
は個々のフレームまたは回転開口フレームのセットに基
づく信号および雑音統計に関する帰納的な仮定によって
予め決定されていてもよい。
レーム映像の共通の重複通過帯域中の時間的に静止した
全ての成分を消去し、運動目標の軌跡の目印だけを残
す。運動目標検出処理方法20-2の性能はシナリオ依存性
であり、等化、差計算および濾波動作によって伝播され
た個々のフレーム測定雑音、およびフレームの整合エラ
ーに対する残留フレームによる背景漏洩によって劣化さ
れる。観察パラメータが許すならば、最適以下の方法に
対して許容可能な性能が実現でき、それによって最適な
整合されたフィルタ55-2のアレイおよび比較器60-2はテ
ンプレート整合されたフィルタおよび比較器60-2、スト
リーク検出器および比較器60-2、或は単なる単一画素ま
たは多重画素しきい値超過検出プロセッサによって置換
されてもよい。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による回転
開口運動目標検出処理方法20-3のトップレベルの実施形
態を示す。図4に示された検出方法20-3は、同じ情景を
観察することによって得られ、米国特許第 5,243,351号
明細書に記載された方法を使用して生成された個々の回
転開口ストリップ開口映像(50-3)またはサブ映像の時
間シーケンスを使用する。照準線整合エラーが存在する
場合、個々の回転開口映像50-3のそれぞれは、米国特許
第 5,243,351号明細書に記載された技術によって空間的
に最適に整合され、共通の空間周波数座標系およびグリ
ッドにフーリエ変換され(51-3)、それによってフーリ
エ変換された映像を生成する(52-3)。典型的に、共通
の空間周波数グリッドは、等価全円形開口回転開口合成
された空間周波数通過帯域と関連している。いくつかま
たは全ての可能性のある空間的に整合されフーリエ変換
された個々のフレーム映像52-3は、時間的にフーリエ変
換され(53-3)、各空間周波数と全ての時間周波数を関
連させるデータセットを生成する(54-3)。
時間周波数および空間周波数変換されたデータセット
は、雑音に最適整合されたフィルタ検出プロセッサ55-3
に伝送される。整合されたフィルタ検出プロセッサ55-3
は、形状、寸法、速度、加速度、背景と目標のコントラ
スト比、信号対雑音比、目標の数、軌道交差の数等のよ
うな潜在的な目標の物理的現象に対応した特性を有する
整合されたフィルタ56-3のアレイを含む。整合されたフ
ィルタ検出プロセッサ55-3はまたシステム伝達関数のバ
リエーションおよび回転開口センサシステム10に関連し
たその他特有の測定特性の過去の知識、並びに測定およ
び変換処理により伝播された測定信号および雑音統計を
含んでいる。検出プロセッサ55-3の各整合されたフィル
タ56-3は、時間的および空間的に変換されたデータのセ
ット54-3に適用される。整合されたフィルタ56-3のそれ
ぞれの出力は比較器60-3に伝送され、しきい値が各整合
されたフィルタ56-3に対して予め定められたレベルに設
定される。運動目標の検出は、整合されたフィルタ56-3
のアレイからの任意の出力が比較器60-3の予め定められ
たしきい値を越えたときに決定される。しきい値は、信
号および雑音統計に関する演繹的な仮定によって、およ
び、または個々のフレームまたは回転開口フレームのセ
ットに基づく信号および雑音統計に関する帰納的な仮定
によって予め決定されてもよい。
々のフレーム映像中の全ての時間的に静止した成分を無
条件に消去し、運動目標の時間および空間周波数の軌跡
の目印だけを残す。運動目標検出処理方法20-3の性能は
シナリオ依存性であり、個々のフレームの測定雑音伝
播、およびフレームの整合エラーに対する残留フレーム
による背景漏洩によって劣化される。観察パラメータが
許すならば、最適以下の方法に対して許容可能な性能が
実現可能であり、それによって最適な整合されたフィル
タ56-3のアレイおよび比較器60-3はテンプレート整合さ
れたフィルタおよび比較器60-3、または簡単な単一周波
数“画素”または多周波数“画素”しきい値超過検出プ
ロセッサによって置換されてもよい。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用した本発明の原理による回転
開口運動目標検出処理方法20-4のトップレベルの実施形
態を示す。図5の検出方法20-4は、同じ情景を観察する
ことによって得られ、米国特許第 5,243,351号明細書に
記載された方法にしたがって生成された2以上の回転開
口ストリップ開口映像30-4を使用する。照準線整合エラ
ーが存在する場合、個々の回転開口映像のそれぞれは、
米国特許第 5,243,351号明細書に記載された方法によっ
て空間的に最適に整合され、共通の空間周波数座標系お
よびグリッドにフーリエ変換される。典型的に、共通の
空間周波数グリッドは、等価全円形開口回転開口合成さ
れた空間周波数通過帯域と関連している。
のフレーム映像30-4のいくつかまたは全ての可能性のあ
る異なる対は、個々のフレームのそれぞれと関連し、か
つ回転開口センサシステム10の回転の結果生じた時間的
に変化するシステム伝達関数を補償するように設計され
た空間周波数フィルタ31-4によって乗算される。システ
ム伝達関数の補償はいくつかの方法の1つによって行わ
れることができ、そのうちの幾つかの例は、ターゲット
全円形開口またはその他の基準システム伝達関数に各フ
レームを等化するか、差を計算されるフレームのシステ
ム伝達関数に1つのフレームを等化するか、或は差を計
算された両映像30-4からシステム伝達関数をデコンボル
ブすることを含んでいる。理想的なシステム伝達関数に
対して、またπラジアン離れて捕捉されたフレームだけ
が差を計算される場合には、システム伝達関数の等化は
不要である。
(40-4)、2つのフレームと関連した重なった通過帯域
の外側の非共通情報を除去するために空間周波数濾波さ
れる(41-4)。空間周波数濾波されたフレーム差は逆フ
ーリエ変換され(42-4)、個々の対の差の映像45-4を生
成する。異なる個々のフレーム映像の数がNである場
合、可能性のある異なる差を計算された個々のフレーム
映像対の総数は、N(N−1)/2であることができ
る。
の可能な等化手順の任意の1つに基づいた回転開口運動
目標検出方法20-4は、特定の方法の受信機動作特性(R
OC)によって特徴付けられる特有の雑音性能を有して
いる。信号テンプレート方法に対するROCの差は、等
化プロセスによって発生された異なる補償されないレベ
ルの雑音“ブースト”または“デブースト”によって支
配される。しかしながら、適切に設計された雑音の最適
整合されたフィルタ方法は、等化プロセスにおいて雑音
特性を補償する。
差を計算された空間周波数濾波され変換された映像45-4
のセットは、雑音の最適整合されたフィルタ検出プロセ
ッサ50-4に伝送される。整合されたフィルタ検出プロセ
ッサ50-4は、形状、寸法、速度、加速度、背景と目標の
コントラスト比、信号対雑音比、目標の数、軌道交差の
数等のような潜在的な目標の物理的現象に対応した特性
を有する個々の対の差のフレーム整合されたフィルタ55
-4のアレイを含む。整合されたフィルタ検出プロセッサ
50-4はまた等化、差計算および濾波処理を受けて伝播さ
れた測定雑音統計の演繹的知識、および回転開口センサ
システム10の特有の信号測定特性を含んでいる。検出プ
ロセッサ50-4の個々の対の差を計算されたフレーム整合
されたフィルタ55-4のそれぞれは、適切なフレーム差に
対して全ての画素位置に適用され、続いて加算される。
各目標の物理的現象の全体的に整合されたフィルタ56-4
の出力は、比較器60-4に伝送され、しきい値が全体的に
整合されたフィルタ56-4のそれぞれの予め定められたレ
ベルに設定される。運動目標の検出は、全体的に整合さ
れたフィルタ55-4のアレイからの任意の出力が比較器60
-4の予め定められたしきい値を越えたときに決定され
る。しきい値は、信号および雑音統計に関する演繹的な
仮定によって、或は個々のフレームまたは回転開口フレ
ームのセットに基づく信号および雑音統計に関する帰納
的な仮定によって予め決定されていてもよい。
レーム映像の共通の重なった通過帯域中の時間的に静止
した全ての成分を消去し、運動目標の軌跡の目印だけを
残す。運動目標検出処理方法20-4の性能はシナリオ依存
性であり、等化、差計算および濾波動作を受けて伝播さ
れた個々のフレーム測定雑音、およびフレームの整合エ
ラーに対する残留フレームによる背景漏洩によって劣化
される。観察パラメータが許すならば、最適以下の方法
に対して許容可能な性能が達成されることができ、それ
によって最適な整合されたフィルタ56-4のアレイおよび
比較器60-4はテンプレート整合されたフィルタおよび比
較器60-4、ストリーク検出器および比較器60-4、或は単
なる単一画素または多重画素しきい値超過検出プロセッ
サによって置換されてもよい。
標検出処理方法20は以下の式により表される。米国特許
第 5,243,351号明細書の定義および用語を使用すると、
本発明において使用される強化された運動目標検出回転
開口センサ10およびデータ処理方法20,20-2,20-3,20
-4の数学的な基礎の説明が与えられる。
する。
られた位置/Rに配置された検出器の積分期間中に生成
された検出電子の平均数を表わしたものである。ここ
で、Os (/F,tj )は空間周波数/Fおよびフレー
ム時間tj における検出器出力の空間フーリエ変換であ
る。実効センサシステム伝達関数が検出器積分時間ΔT
中時間的に静止していると仮定すると、強度を検出され
た個々のフレームのフーリエ変換は次のように表わされ
る
空間的に平均された強度であり、《W0 》は検出器伝達
関数のゼロ空間周波数成分であり、AfpaSは実効焦点平
面アレイ面積であり、STF(/F,tj )は、フレー
ム時間tj に対応した全体的なセンサシステム伝達関数
であり、《Sn 》(M/F,t)は時間tにおける情景
の幾何学的な投影の正規化された空間フーリエ変換であ
る。
分を含む情景および映像生成平面への情景の幾何学的投
影の数学的表示を説明する。背景の前方(回転開口セン
サに関して)および透過散乱層の後方を運動している単
一の目標を含む情景から放出された強度は次のように表
わされる。
における運動目標の強度分布であり、B(/r,t)は
情景位置/rおよび時間tにおける背景の強度分布であ
り、H(/r,t)は情景位置/rおよび時間tに対応
した透過散乱層の強度分布であり、WT (/r,t)は
情景位置/rおよび時間tにおける目標の逆シルエット
である。目標の逆シルエットは次の式によって与えられ
る。
t (/r,t)を減算することは、情景位置/rおよび
時間tにおいて運動目標によってぼかされ背景光子を除
去することに対応する。
の情景の強度の幾何学的な投影は、次の式によって表わ
される
をフーリエ変換して、フーリエ変換類似性の定理を適用
することにより次式が得られる。
変数の変換は、
成分の対応した幾何学的投影は、
トルは、次のように表わされる。
クトル並びにスペクトルに対する運動目標の貢献によっ
て支配される2つの成分に分解されることが可能であ
り、これは次のように表わされ、
周波数M(/F)に対する正規化された背景および前景
成分であり、
おける空間周波数/Fに対する背景および前景からの正
規化された目標のバリエーションである。
化された背景および前景スペクトルを説明する。時間的
に静止した全体的な背景成分を有する観察シナリオに対
して、正規化された全体的な背景スペクトル、および運
動目標バリエーションは次のようになる。
標によってぼかされる背景光子の数が、合計背景光子数
に比較して無視できるほど僅かである場合、合計背景お
よび運動目標バリエーションの正規化されたスペクトル
はさらに簡単化されることができ、次式を生成する。
シルエットを説明する。目標空間強度分布の論証に時間
依存性の空間シフトを適用することによって目標の時間
依存性運動を説明することができる。映像生成プロセス
の開始時間をt0であると定めると、最初の目標および
目標逆シルエットは次式によって与えられる。
ーションは次式によって与えられる。
における目標のベクトル位置であり、目標軌道と呼ぶ。
トルを説明する。フーリエ変換シフト定理を運動目標お
よび運動目標逆シルエットに適用すると、以下の式によ
ってそれぞれ与えられる対応したスペクトルが生成され
る。
らの運動目標バリエーションを説明する。《T》{M
(/F),t}および《WT 》{M(/F),t}に対
する上記の式を目標スペクトルバリエーションΔ/TN
の式に代入すると、次式が得られる。
ベクトル位置は/RT (t)=/Vtとして表わされる
ことができ、ここで/Vは運動目標の速度成分である。
したがって、正規化された目標スペクトルバリエーショ
ンΔ《Tn 》は次のように表わされる。
説明する。映像化された運動目標および全体的な背景の
時間平均の結果、平均検出器測定が得られるため、シス
テム伝達関数が検出積分時間中一定であると仮定した測
定の時間平均されたフーリエスペクトルは次のようにな
る。
運動の時間平均は、
明する。回転開口合成プロセスの最も簡単なバージョン
は、無雑音測定の空間周波数冗長平均が等価全円形開口
映像を生成するようにシンセサイザを使用する。すなわ
ち、
形開口映像であり、《Λ》(/F)はシンセサイザであ
る。また《OS ave 》(/F)は個々の(無雑音)フレ
ーム測定スペクトルの空間周波数冗長平均である。
光学系伝達関数OTFと、正規化された検出器伝達関数
《WN 》との積であり、下付きの文字SおよびFは、回
転開口ストリップ開口センサおよびターゲット全円形開
口システムと関連したパラメータをそれぞれ示す。シン
セサイザの分母に現れる空間周波数冗長平均されたスト
リップ開口伝達関数は、
止した全ての背景および前景からなる情景の観察から結
果的に得られ、センサの照準線における運動目標が次式
によって与えられる。
観察している合成された映像の任意の対に対して同じな
ので、合成された映像の差は運動目標情報だけを含んで
いる。差を計算された合成された映像の無雑音空間周波
数ドメインの運動目標目印は、次の式によって与えられ
る。
は、
異なるセットの獲得に関連した時間遅延を表わす。
めに適切に修正された回転開口フレームの任意の対の差
の計算について説明する。2Nf 個の個々の回転開口ス
トリップ開口フレームに対して、2Nf (Nf −1)/
2の可能な異なった対でフレーム比較、すなわち差が存
在している。最適な検出技術は、任意の可能な対および
潜在的に非常に多数の運動目標の対となった目印の差を
計算するときに背景除去を利用する。しかしながら、時
間的に静止した開口システムに対して背景除去および運
動目標検出のためにフレームの差を計算するのとは異な
り、特定の技術および方法が、システム伝達関数の回転
特性を補償するために回転開口データセットに適用され
なければならない。回転開口データセットは、情景にお
ける時間変化と運動目標から派生した時間的変化とを区
別するように修正される。
よびフレーム毎のシステム伝達関数等化を説明する。個
々のストリップ開口雑音測定スペクトルは、無雑音スト
リップ開口検出映像スペクトル、プラスストリップ開口
の加算された雑音スペクトルとして表わされることがで
きる。
号および雑音スペクトルである。
過帯域フィルタを説明する。一般に、個々のフレーム測
定は、等価全円形開口光学系通過帯域の異なった部分に
関する情報を提供する。運動目標または過渡事象のよう
な時間的に静止していない現象を検出するには、フレー
ム間に共通した情報を識別して使用し、また回転開口セ
ンサの伝達関数の時間的に変化する特性による映像コン
トラストの差に対して調節を行う必要がある。フレーム
毎のシステム伝達関数補償は、いくつかの方法によって
行われることができ、そのうちの数例は、ターゲット全
円形開口またはその他の基準システム伝達関数にそれぞ
れのフレームを等しくするもの、差を計算されるフレー
ムのシステム伝達関数に1つのフレームを等しするこ
と、差を計算される両方の映像からシステム伝達関数を
デコンボルブすること、或は理想的なシステム伝達関数
に等しくすること等であり、またπラジアン離れて獲得
されたフレームだけが差を計算された場合にシステム伝
達関数の等化は不要である。テンプレート信号検出方法
に対して、多数の可能な等化過程の任意の1つに基づい
た回転開口運動目標検出方法は、特定の方法の受信機動
作特性(ROC)によって特徴付けられる特有の雑音特
性を有している。信号テンプレート方法に対するROC
における差は、等化プロセスによって発生された異なる
補償されないレベルの雑音“ブースト”または“デブー
スト”によって支配される。しかしながら、適切に設計
された雑音の最適整合されたフィルタ方法は、等化プロ
セスにおいて雑音特性を補償する。
の方法の任意の1つを使用することができるが、以下の
段落において説明される式は、各過程と関連した数学的
詳細の典型である。簡略化するために、1つのフレーム
が別のフレームに等しくされる過程だけを詳細に説明す
る。
ムと、比較のために使用される任意の異なる別のフレー
ムとを選択すると、フレームのスペクトルは、両方のフ
レームの通過帯域の共通の重なった領域における空間周
波数情報を提供するように濾波された空間周波数とな
る。この別のフレームはまたシステム伝達関数の差を補
償するために修正される。すなわち、
を濾波された空間周波数を示し、《mE 》(F,ti )
は比較されるべきフレームの出力を濾波され、伝達関数
等化された空間周波数を示し、《x》は光学系通過帯域
フィルタであり、
差を計算された信号および雑音スペクトルの和として表
されることができる。すなわち、
ペクトル《ΔNij》は、
クトル《ΔSij》(/F)は、
信号スペクトルは、以下の式によって背景からの正規化
された運動目標バリエーション間の差Δ《TN 》に関連
させられる:
合、濾波された差の信号は以下の関係に対応する:
の運動目標の検出を説明する。空間周波数および時間周
波数ドメイン回転開口運動目標の検出過程の簡単な数学
的記述は、ディスクリートな時間サンプリングを連続体
として近似し、個々の回転開口フレームのそれぞれに対
して空間周波数変換の連続したフーリエ変換を行うこと
によって得られる。個々の無雑音フレーム測定のフーリ
エスペクトルは以下の式によって与えられる:
エ変換は、積分を評価してコンボリューション定理を適
用することによって以下のようにして得られる。
時間フーリエ変換であり、《SN 》は正規化された幾何
学的に投影された情景スペクトルの時間および空間周波
数フーリエ変換である。
ペクトルの時間フーリエ変換を評価すると、ストリップ
開口測定の無雑音時間周波数変換は次のようになる。
は、時間および空間周波数整合されたフィルタの信号成
分を決定するために使用される。{OS }(F,f)に
対する式の最初の項は時間的に静止した結果を全て含ん
でおり、かつ全てのゼロでない空間周波数に対してゼロ
であるため、整合されたフィルタは本質的に時間的に静
止した成分を消去する。
対する簡単化された目印を説明する。回転開口運動目標
の検出過程に対する受信機動作特性の予測は、均一な強
度の背景、前景および運動目標に関連した観察シナリオ
の特殊な場合に対して簡単化されることができる。事
実、粗運動目標検出方法の設計はこのような簡単化され
た仮定を使用することができる。背景および前景の時間
的に静止しており、かつ空間的に均一な成分は次のよう
に表されることができる。
り、HSCN は情景座標における平均前景強度である。
を説明する。情景の背景および前景強度の平均値の幾何
学的な投影は、
る。
リエ変換が行われる。:
を説明する。情景スペクトルの正規化された幾何学的な
投影は、幾何学的に投影された情景スペクトルのゼロ空
間周波数成分で情景の幾何学的投影のフーリエ変換を除
算することによって生成される。背景の強度および運動
目標の強度の幾何学的に投影されたゼロ空間周波数成分
は、空間的に平均された平均値と、背景および目標の幾
何学的な投影された面積に関して以下のように表される
ことができる:
的投影面積であり、Afpaは回転開口センサの焦点平面
における背景のおおよその検出された面積であり、文字
上の符号(−)は、焦点平面アレイにおける空間平均を
示す。
周波数関係を説明する。情景座標系における目標逆シル
エットにフーリエ変換の定義を適用し、Mrscene =R
image により与えられる変数の変化によって映像平面座
標に変換すると、以下の式が得られる:
的な投影はまた、以下の式によって焦点平面における目
標の実効面積に関連されることができる:
下に説明される。情景の均一な背景成分に対して前述の
関係を適用することによって、正規化された背景および
正規化された目標の数学的記述を決定することができ
る。背景および前景の投影情景スペクトルの合計に対す
るそれらのゼロ空間周波数成分の比率は以下の通りであ
る。
と、正規化された目標との差《TN 》は、以下のように
表される。
おいて、均一な目標は、目標の逆シルエットWT と目標
の一定の強度レベルTO との積として表されることがで
きる。
され、
換は以下のように表わされる。
止した情景の、正規化された目標の背景からの変化Δ
《TN 》は、以下の式によって与えられ、
うに定義され、
成分は、以下のように表される。
リエ変換シフト定理を適用すると、目標の逆シルエット
は次式のようになる。
(M/F,t0 )は、時間=t0 における目標の正規化
された逆シルエットのフーリエ変換である。結果的に、
背景からの目標の変化の式は以下のように表される。
る、差を計算され濾波された運動目標の空間的周波数の
目印は次式によって表される。
計算された信号が以下に説明される。運動目標の差が計
算された濾波された空間周波数スペクトルを逆フーリエ
変換すると、空間ドメインの運動目標の目印は以下のよ
うに表される。
が以下に説明される。均一な目標および背景に対する位
置/Rにおける平均差信号対RMS差雑音の比率は次式
によって表される。
対の信号スペクトルが表される。均一な目標および背景
と関連した任意の対の差の合成された映像の雑音のない
空間周波数ドメインの目印は、次式によって表され、
インの運動目標は、以下のように表される。
幅広いクラスの性能基準に対して最適に動作する検出器
プロセッサが、以下、比較器と呼ばれる比較装置を有す
る一般化された整合されたフィルタのアレイに構成され
る。スピンAp運動目標検出プロセッサのための最適に
整合されたフィルタアレイの設計には、スピンApセン
サの特有の信号測定特性と、整合されたフィルタアレイ
に入力される前の所望されたデータ処理と、予測された
目標の現象学とに関する知識を組入れることが要求され
る。過去の目標知識は、本来統計的あるいは決定論的で
あり、形状、寸法、速度、加速度、背景対目標のコント
ラスト比、信号対雑音比、目標の数、軌道交差の数、雑
音の交差共分散等の物理的パラメータを含んでいる。特
定の目印に対する整合されたフィルタのアレイは、行列
法で以下のように表されることができ、 F=St CNN -1 ここにおいて、Sは過去の雑音のない信号のベクトルで
あり、CNN -1は雑音の交差共分散逆マトリックスであ
る。信号および雑音の交差共分散マトリックスの両方
が、所望された運動目標検出データ処理方法に適切に定
義される。
周波数のバリエーション、あるいはスピンApによる運
動目標の検出に適切となるように定義された別の目印を
任意に組み合わせることによって生じた信号マトリック
スは以下のように表され、
的なスピンAp運動目標目印と関連した独立した空間、
空間周波数、時間、時間周波数の信号の座標の任意の全
ての組み合わせを示している。
像、時間周波数の変化、あるいはスピンApによる運動
目標の検出に適切となるように定義された別の適切な目
印を任意に組み合わせることによって生じた雑音マトリ
ックスは以下のように表される。
音成分の合計として行列法で以下のように表される。
マトリックス成分がそれぞれ複素数の信号および雑音の
実数および虚数部分(すなわち、振幅および位相部分)
であることを認識することによって構成されることがで
きる。
ィルタの捕捉“ウィンドウ”内の任意の数の差のフレー
ムあるいは差の合成された映像の信号および雑音マトリ
ックスは、それぞれ以下のように表され、
タの捕捉ウィンドウにおける全ての成分の数を示してい
る。
されたフィルタの出力は、整合されたフィルタの雑音の
ない信号および付加的な“白色”雑音に対応し、次式に
よって表される。
特性が以下に説明される。比較器およびしきい値超過決
定方法を使用するスピンApの整合されたフィルタによ
る運動目標検出処理方法の受信機の動作特性(ROC)
は、個々のフレーム測定の近似のガウス雑音を求めるた
めに容易に決定されることができる。所定のしきい値の
スピンAp整合フィルタの出力に対する検出確率および
誤警報確率は次式によってそれぞれ決定され、
タである。
ックス成分が説明される。雑音が最適なスピンAp運動
目標検出整合フィルタの決定は、雑音の処理および測定
ステップから整合されたフィルタへの雑音の伝播に関す
る過去の知識あるいは仮定に基づく。本開示において説
明されたスピンAp目標検出方法のクラスは、個々の加
重されたフレームの差、合成されたフレームの差、ある
いは時間および空間周波数整合されたフィルタを使用す
る。適切に伝播された雑音の交差共分散マトリックス成
分の決定には、雑音の最適な整合されたフィルタを構成
することあるいは雑音最適化技術以下の性能が必要とさ
れる。以下のセクションによって、適切なスピンAp雑
音共分散マトリックス成分を決定するための例示的な技
術が提供される。はっきりとは説明されていない処理ス
テップに関連した共分散マトリックス成分は、類似した
決定の方法に従う。
が説明される。スピンAp運動目標検出方法の一実施形
態によって、等化された異なる対のシステム伝達関数
と、オーバーラップした通過帯域で濾波された個々のフ
レームとの差が計算される。対応する差が計算され、等
化され、濾波された雑音スペクトルNijは次式によって
表される。
た雑音スペクトルの逆フーリエ変換を行うことによって
以下のように得られる。
になる。
Rm )および(/Rn )における任意の2対の差が計算
され、等化され、濾波されたフレーム間の空間雑音交差
共分散は以下の次式で表され、
る。
雑音の効果が説明される。ゼロ平均で、時間的に静止
し、時間的に相関しない雑音に対して、個々のフレーム
の雑音スペクトルの交差共分散は、以下の式
れ、ここにおいて、δi,j およびδ(/F−/F)は、
それぞれクロネッカーおよびディラックのδ関数であ
る。
分散を空間的に差が計算された濾波された雑音の共分散
の式に代入すると、以下のような式が生成される。
差の画素毎の交差共分散は、フレーム指数がi=k、j
=l、およびi≠jの値をとるようにすることによって
得られる。対応する実効雑音空間的交差共分散は、以下
の式によって示される。
いて説明する。空間的白色雑音および光学系の通過帯域
のナイキスト・サンプリングに対して、個々のスピンA
pのフレームの雑音のスペクトル密度は以下のように近
似され、
実効焦点距離であり、DxおよびDy はそれぞれxおよ
びy方向の開口の寸法であり、Δxs およびΔysは空
間的サンプリングのインクレメントであり、それは、検
出器の間隔に対応する100%の充填率を有する焦点表
面検出器アレイに対するものであり、/σds 2 は個々の
検出器における雑音の平均的分散である。個々の白色雑
音のスペクトル密度の交差共分散を、差が計算され、濾
波された空間の雑音の共分散の式に代入すると、以下の
式が生成される。
ら画素の交差共分散に対応する特別な場合は、
雑音について説明する。理想的なシステム伝達関数に対
して、システム伝達関数の等化を必要としないという処
理上の利点を有するさらに簡単にする1例は、πラジア
ン離れたフレームの差を計算することに対応する。理想
的な方形ストリップ開口に対して、πラジアン離れた2
つのフレームに対するシステム伝達関数は、同じコント
ラスト性能を提供し、したがって等化フィルタは正確に
1であり、すなわち《H》(/F,tj ,ti )=1で
あり、重なった通過帯域は同一である。したがって、空
間的に白色で時間的に相関しない雑音に対して、交差共
分散雑音マトリクス成分は、
とができ、
マトリクスの逆数は、
口運動目標検出処理方法の別の実施形態は、回転開口合
成された映像の任意の異なった対の差を計算する。米国
特許第 5,243,351号明細書に記載されているように、最
も簡単な形態の合成プロセスに対する合成された雑音ス
ペクトルは、
(/F)は次式によって与えられる:
た雑音スペクトルを与えると、雑音スペクトルΔ《n
syn 》(/F)は、合成された画像の任意の対(完全に
異なったデータセットから得る必要はない)では次式の
ようになる。
めに使用されたスピンAPの個々のフレームセットを意
味する。
明する。差の合成された雑音スペクトルの逆フーリエ変
換を行うと、画素位置(/R)で空間的ドメインの合成
された雑音が発生する。
て説明する。合成された雑音の差の画素毎の交差共分散
は、スピンApシンセサイザと冗長度平均された差の雑
音スペクトルの交差共分散に関して表されることがで
き、即ち、
ようになる。
い雑音を説明する。ゼロ平均に対して時間的に静止して
いる相関されていない雑音では、個々のフレーム雑音ス
ペクトルの交差共分散は次式により雑音パワースペクト
ル密度《CNN》に関連される。
データを有するスピンApのセットの個々のフレームに
基づいて合成画像のポテンシャルを適応する。
関について説明する。異なったセットのスピンApの個
々のフレームから生成される合成された画像に対して
は、差の合成された画像についての雑音スペクトルの交
差共分散は以下の式で表される。
する画素と画素の差の合成された雑音は次式で与えられ
る。
白色の空間的な差の合成された雑音の交差共分散マトリ
ックスは以下の式で与えられる。
のスピンAPの運動目標システムの実施形態が実施され
試験された。画像データは図2の実施形態を用いて生成
され、したがって雑音のない差の合成された運動目標の
目印と、対応するテンプレート整合フィルタ出力を生成
する。
む低い信号対雑音比の差の合成された画像と、テンプレ
ート整合フィルタの出力を示し、したがって、16対1
程度の優れた出力の信号対雑音比を示す。このことはこ
の特定のクラスの運動目標に対して本発明を使用して生
成される強化された受信機の動作特性を示している。さ
らに、実験データが得られ、ここで説明したπ方向の運
動目標検出方法を証明するために使用された。
し、複数の回転ストリップ開口画像測定または合成され
た画像測定から運動目標を検出する回転ストリップ(部
分的)開口画像放射線センサを説明した。前述の実施形
態は本発明の原理の応用を表した多数の特定の実施形態
のうちの数例にすぎないことが理解されよう。明白に、
多数および変更された他の装置は本発明の技術的範囲を
逸脱することなく当業者により容易に製造されよう。
を使用する回転開口映像生成放射線センサシステムの概
略図。
た映像を使用する本発明の原理による回転開口運動目標
検出処理方法のトップレベルの実施形態の概略図。
化された個別の回転開口変換されたフレームを使用する
本発明の原理による回転開口運動目標検出処理方法のト
ップレベルの実施形態の概略図。
レーム映像の時間および空間周波数変換を使用する本発
明の原理による回転開口運動目標検出処理方法のトップ
レベルの実施形態の概略図。
化された個別の回転開口変換されたフレーム、および順
次差フレーム整合された濾波を使用した本発明の原理に
よる回転開口運動目標検出処理方法のトップレベルの実
施形態の概略図。
用する運動物体および過渡事象検出
定装置および放射線検出方法に関し、特に、合成された
等価全円形開口が映像を生成し、複数の回転ストリップ
開口映像からの運動目標の検出または合成された映像の
測定を行う回転ストリップ(部分的)開口映像生成用放
射線測定装置に関する。
高い解像度の映像を得るために、例えば通常のセンサ技
術では大型で、重く、展開可能な光学システムの能動制
御装置を使用している。要求される使命および1次反射
器の大きさに応じて、センサの全ての光学素子の位置、
傾斜、形状、および整列制御のために1次反射器の周期
的な位置および傾斜制御から配列することができる。本
発明と同じ解像度を有する全開口システムは比較的大き
い開口面積のために多量の集光能力を有している。しか
しながら、軌道上に数十メートルの直径の開口を配置す
るためには、本発明と競合する全開口システムは次のよ
うなことが必要になる。すなわち光学システムの直径が
発射装置よりも大きくなるとセグメント化された光学表
面および折畳まれた支持構造が必要になる。また、重量
を節減するために薄い変形可能な反射鏡が使用される場
合には複雑で潜在的に高い帯域幅の適応光学技術が必要
であり、さらに、位相アレイとして構成される場合には
複雑な位置および瞳孔整合制御が要求される。それ故、
全開口システムは本発明と比較して比較的重く、高い技
術的なリスクを有している。
スピンApと呼ばれる最も近い従来技術のシステム能力
を改善させるものである。このスピンApシステムおよ
び方法は米国特許第5243351号明細書に記載され
ている。このスピンApシステムは全円形開口に等価な
高い解像度を有する映像を合成することのできるデータ
処理方法である。等価なスピンAp合成全円形開口はス
トリップ開口形状と関連する最大相関長に等しい直径を
有する。等価全円形開口映像を正確に合成するために、
スピンApは、等価全円形開口の光学システム通過帯域
を完全に測定するために適切な時間インターバルで回転
するセンサの焦点平面検出器を使用する映像データを捕
捉する。映像は米国特許第5243351号明細書に記
載された方法により処理される。
開口運動目標検出と呼ぶ)との共通点は、全円形開口に
等価な解像度を有する映像を合成するために時間的に整
合したストリップ開口を使用することである。等価回転
開口の合成された全円形開口はストリップ開口の形状に
関係した最大の相関長に等しい直径を有する。等価全円
形開口映像を正確に合成するために、スピンApは適切
な時間インターバルで回転するセンサの焦点平面検出器
を使用する映像データを捕捉、し等価全円形開口の光学
システム通過帯域を完全に測定する。映像は米国特許第
5243351号明細書に記載された方法によって処理
される。
のの瞬間的視野中の運動目標は全体の合成された映像中
に運動による軌跡の目印を生じる。運動目標の検出は、
位置、寸法、軌跡、速度、加速度、および原点のような
目標現象を設定するために必要である。さらに、検出プ
ロセスは運動目標の静止映像合成を可能にする捕捉され
たデータにプロセッサおよび方法を適用するための前段
階である。それ故、運動目標検出プロセッサおよび方法
の階級について説明する。この階級は、空間的、時間
的、空間周波数の、時間周波数ドメイン検出プロセッサ
を含み、それはまた多重スペクトル背景排除技術を伴っ
ている。本発明のシステムは、スピンApの個々のフレ
ーム中の任意の1つにおけるスピンApの任意の対の視
野中の合成された映像の視野中の運動目標を検出するた
めに利用できる測定を使用する。同じ方法はまた過渡事
象の検出にも適用できる。
は、所定の実効合成開口寸法に対して低いペイロード重
量が得られ、一方センサの視野内の運動目標の検出を可
能にする変形されたスピンAp映像処理方法を提供する
ことである。さらに本発明の別の目的は、宇宙空間で使
用するための複雑な適応光学システムを必要とする複雑
で大型の全円形開口(または位相アレイ望遠鏡)を展開
するのに代えて、低いリスク、低いコスト、軽い重量、
および簡単な展開構造の映像処理および運動目標検出方
法を提供することである。
目的を達成するために、本発明の1形態においては、運
動目標を検出することのできる回転ストリップ放射線映
像形成および検出システムが提供される。このシステム
は時間的に連続する、または順番に並んだ映像を生成す
る回転ストリップ開口望遠鏡を備えている。回転ストリ
ップ開口望遠鏡は典型的にストリップ開口1次反射器お
よび2次反射器を備えている。2次元の検出器アレイが
設けられて望遠鏡の焦点平面に位置する映像を検出す
る。回転補償装置が使用されてアレイの検出器の積分時
間中に回転の軌跡が描かれることを阻止している。信号
プロセッサが設けられてストリップ開口が望遠鏡の光学
軸を中心に回転するとき望遠鏡によって映像が形成され
た目標情景の複数のフレームを記録し、等価全円形開口
映像を合成し、個々のストリップ開口フレーム、任意の
対のフレーム、合成された映像、または個々のフレーム
と合成された映像の組合わせにおける運動目標を検出す
る。
像を合成し、複数の回転ストリップ開口映像測定からロ
センサの視野中の運動目標を検出し、一方、個々のスト
リップ開口映像間のランダムおよび/または系統的な照
準線エラーを補償する。したがって、本発明は等価全円
形開口映像を合成するための手段を提供することによっ
て米国特許第5243351号明細書に記載されたスピ
ンApシステムに比較して改良された全体的なシステム
性能を提供し、同時に運動目標を検出する能力を与える
ことができる。
いて以下説明する。放射線測定システム、映像合成、お
よび運動目標検出処理方法の各実施形態は特定の使命の
要求に依存し、また放射線測定システム、映像合成、お
よび運動目標検出方法の技術的妥協に依存し、フレーム
間の照準線のランダムまたは合成によるドリフトの補
償、共通平面にない点の広がり機能、フィールド点と関
係のない映像の運動による映像の軌跡と、検出器の有限
の電子的帯域幅、フィールド点に依存する映像の回転運
動が補償されないことによる軌跡の補償等も必要に応じ
て選択される。
よび検出プロセスの実施形態について以下説明する。回
転ストリップ開口が望遠鏡の光学軸を中心に回転すると
き次のような動作が生じる。回転補償装置は検出器の積
分時間中逆回転し、それによって時間的に静止した映像
を与える。映像フレームは記録され保存される。回転
(情景に対して)検出器アレイ構造が選択される場合、
得られたフレームは合成された映像の基準グリッドに座
標変換され、補間される。データはフーリエ変換されて
記憶される。新しいフレームが記録され、保存される。
間での照準線エラーによる不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換はそれらの照準線
エラーに対して補正され、記録される。先行段階は各ス
トリップ開口映像フレームに対して順次反復され、また
はフレームは順次捕捉され、記憶されて、その後照準線
エラーの全体的な評価が決定され、フレームの整合に使
用される。個々のストリップ開口映像の任意の対または
対の組合わせが照準線エラーに対して補正されると、回
転開口運動目標検出プロセスが開始される。さらに、一
度十分な数(前記米国特許第5243351号明細書に
記載されたような)のストリップ開口映像測定が得られ
ると、映像合成処理が開始される。回転開口運動目標検
出方法は2回以上合成映像について行われる。
た2以上の合成映像を使用する回転開口運動目標検出方
法の実施形態は合成された映像を空間的に整合させ(測
定セット間の照準線エラーにより要求される場合)、合
成された映像のすべての可能な異なる対の差が計算さ
れ、差の映像が合計され、整合されたフィルタを通って
差および合計データが送られ、比較器で比較される。整
合されたフィルタからの出力が比較器の予め定められた
しきい値を越えたとき運動目標が検出されたことが決定
される。
2以上の個々の照準線の補正されたフレームを使用する
回転開口運動目標検出方法の別の実施形態は、システム
伝達関数(STF)等化フィルタを差を得る個々のフレ
ームに適用して等化されたフレームの任意のまたは全て
の対の差を生成し、得られた差のデータに対して共通で
ない通過帯域排除フィルタを適用する。濾波された差の
データは濾波された差の全体の合計されたスペクトルを
形成するために加算により結合される。加算されたスペ
クトルは逆変換され、それによって空間的ドメインにお
いて運動目標の目印を与え、それは雑音を最適にする整
合されたフィルタアレイおよび比較器に供給される。整
合されたフィルタアレイからの出力が比較器の予め定め
られたしきい値を越えたとき運動目標が検出されたこと
が決定される。
ルタ検出プロセッサは、形状、寸法、速度、加速度、背
景と目標のコントラスト比、信号対雑音比、目標の数、
軌道交差の数、等の潜在的な目標の物理的現象に対応す
る特性を有している。整合されたフィルタ検出プロセッ
サはまた等化、差の形成、濾波、ならびに回転開口セン
サシステムの特有の信号測定特性により伝播された測定
雑音統計の過去の知識を有している。
作は、個々のフレーム映像の共通の重なった通過帯域に
おける全ての時間的に静止している成分を除去し、また
合成映像は軌跡を描く運動目標の雑音を含む目印のみを
残す。この運動目標検出技術の性能はシナリオに依存
し、等化、差の形成、濾波動作を行われる個々のフレー
ム測定雑音、および残留フレーム対フレームによる背景
漏洩および、または合成された映像整合エラーに対する
合成された映像による背景漏洩によって劣化される。
法の別の実施形態は、任意の差を形成されたフレームに
対して最適の整合フィルタを適用し、全ての可能なフレ
ーム差に対して出力を組合わせることに対応する。
の典型であるが、時間周波数ドメインプロセスもまた使
用可能である。時間周波数ドメインプロセッサの1実施
形態は同じ情景を観察することによって得られた空間的
に整合された個々の回転ストリップ開口フレームの時間
的シーケンスを利用する。空間的に整合されたフレーム
は共通空間的周波数グリッドへ空間的にフーリエ変換さ
れる。共通空間的周波数グリッドは等価全円形開口回転
ストリップ開口合成空間周波数通過帯域と関係する。任
意のまたは全ての空間的に整合され、フーリエ変換され
た個々のフレーム映像は各空間周波数を有する全ての時
間的周波数と関連するデータセットを生成するために時
間的なフーリエ変換をされる。最良の可能な検出特性を
得るために、時間的および空間的に周波数変換されたデ
ータセットは雑音最適化された整合されたフィルタ検出
プロセッサに転送される。整合されたフィルタ検出プロ
セッサは目標現象、信号対雑音測定の統計、ならびに回
転ストリップ開口センサシステムのユニークな測定特性
についての過去の知識を有している。整合されたフィル
タアレイからの出力が予め定められたしきい値を越えた
とき運動目標が検出されたことが決定される。
整合されたフィルタは個々のフレームの変換された時間
ヒストリーにおける全ての時間に静止している成分を消
去し、それによって運動目標の時間的周波数および空間
的周波数目印のみを残す。運動目標検出技術の性能はシ
ナリオに依存し、個々のフレーム測定雑音伝播およびフ
レーム整合エラーに対する残留フレームによる背景漏洩
によって劣化される。観察パラメータが許容すれば、許
容可能な性能は最適の方法より少ないものに対して達成
され、それによって最適に整合されたフィルタアレイお
よび比較器はテンプレート整合されたフィルタおよび比
較器、ストリーク(筋)検出器および比較器、または単
なる単一画素または多重画素しきい値超過検出プロセッ
サによって置換されることができる。これにおいて画素
と言う用語は、映像の2次元のディスクリートな表示の
空間的、または空間周波数素子のことである。
よび方法の階級が与えられる。階級は空間、時間、空間
周波数、および時間周波数ドメイン検出プロセッサなら
びに多重スペクトル背景排除技術を含んでいる。適切な
処理手順および方法の選択はシナリオおよび運動目標の
実効的信号対雑音比特性に依存している。同じ方法は過
渡事象の検出にも適用することができる。
識、アレイの検出器の空間的応答特性、雑音統計、記録
された各ストリップ開口映像の時間的整合はこのセンサ
システムおよび処理方法による運動目標検出を可能にす
る。
を参照することにより本発明の種々の特徴および利点が
容易に理解されるであろう。おな、同じ参照符号は同じ
構造素子を示すものである。図面を参照すると、図1は
本発明の原理を使用して適切に修正された回転開口映像
生成放射線センサシステム10およびデータ処理方法20
(米国特許第 5,243,351号明細書に記載されたシステム
に対応している)を示す。回転開口映像生成放射線セン
サシステム10およびデータ処理方法20は、運動目標の検
出情報を提供し、等価全円形映像を生成すると同時に照
準線ジッタおよびフレーム整合エラーを除去する能力を
提供する。したがって、ここに記載されたシステム10お
よびデータ処理方法20は、米国特許第 5,243,351号明細
書の回転開口(SpinAP)システムを改良したものである。
米国特許第 5,243,351号明細書の内容全体がここにおい
て参考文献に含まれている。
は、1次反射器12aおよび2次反射器12bを含む回転ス
トリップ開口望遠鏡11を含んでいる。3次反射器(示さ
れていない)は、ある環境下で望遠鏡11において使用さ
れてもよい。システム10は、地球指向アンテナ13aを有
する静止セクション13を含む衛星の形態で示されてい
る。望遠鏡11はプラットフォーム14上に取付けられ、静
止セクション13aもまたそれに結合されている。回転開
口映像生成放射線センサシステム10は、1次反射器12a
(ストリップ開口12を含む)が望遠鏡の光学軸を中心と
して回転したときに、望遠鏡11によって映像化されたタ
ーゲット情景の多数の映像フレームを整合するように構
成されている。
装置16は、検出器アレイ17によるその感知の前に映像を
安定させ回転を除去するように構成され、望遠鏡11の光
路に沿って配置されている。回転除去装置16は、回転補
償制御装置19の制御の下に検出器アレイ17を含む検出器
の積分時間中映像を逆回転させ、それによって静止映像
を生成する。照準線安定化ミラー15は、広い帯域の照準
線エラーおよびシステム10の軌道ダイナミックスによる
照準線エラーを除去するために照準線制御システム(信
号プロセッサ18またはその他の専用制御システムにより
提供されてもよい)によって使用される。
るように構成された信号プロセッサ18に結合された望遠
鏡11の焦点平面に配置された検出器アレイ17上に投影さ
れ映像化される。個々の映像フレームは、等価全円形開
口映像を合成するために信号プロセッサ18において処理
され結合される。この強化された運動目標検出回転開口
システム10において、信号プロセッサ18(または独立し
たプロセッサ)は運動目標を検出するために本発明の方
法にしたがって修正される。以下、図2乃至5を参照し
て修正された処理方法を詳細に説明する。
テムにより使用される本発明の原理によるセンサシステ
ム10および処理方法20,20-2,20-3,20-4,の種々の異
なる実施形態を示す。観測シナリオおよびタイムライ
ン、利用可能な過去の知識の量、および処理チェインの
利用可能な計算処理に基づいて、特定の実施形態が選択
されている。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による簡単
な回転開口運動目標検出処理方法20のトップレベルの図
を示す。図2に示された検出方法20は、米国特許第 5,2
43,351号明細書にしたがって情景の合成された等価全円
形開口映像30を生成するために同じ情景を観察する2以
上の回転開口ストリップ開口映像セットを使用する。照
準線整合エラーが存在する場合、合成された映像30のそ
れぞれは米国特許第 5,243,351号明細書に記載された方
法によって空間的に最適に整合される。
くつかまたは全ての可能性のある別個の対は、差を計算
され40、続いて加算されて35、記憶され、それによって
差を計算された合成された映像の総和を生成する45。別
個の合成された映像30の数がNである場合、可能性のあ
る別個の差を計算された合成された映像対の総数は、N
(N−1)/2によって与えられる。
差を計算された合成された映像の結果的な和45は、雑音
の最適整合されたフィルタ検出プロセッサ50に送られ
る。整合されたフィルタ検出プロセッサ50は、形状、寸
法、速度、加速度、背景と目標のコントラスト比、信号
対雑音比、目標数、軌道交差数等のような潜在的な目標
の物理的現象に対応した特性を有する整合されたフィル
タ55のアレイを含んでいる。整合されたフィルタ検出プ
ロセッサ50はまた合成および差計算処理により伝播され
た測定雑音統計の過去の知識、および回転開口センサシ
ステム10の特有の信号測定特性に関する知識を含んでい
る。検出プロセッサ50の整合されたフィルタ55のそれぞ
れは、差を計算された合成された映像の総和45のあらゆ
る画素位置に適用される。整合されたフィルタ55のそれ
ぞれの出力は比較器60に伝送され、しきい値が各整合さ
れたフィルタ55の予め定められたレベルに設定される。
運動目標の検出は、整合されたフィルタ55からの出力が
比較器60の予め定められたしきい値を越えたときに決定
される。しきい値は、信号および雑音統計に関する演繹
的な想定によって、或は個々のフレームまたは回転開口
フレームのセットに基づく信号および雑音統計に関する
帰納的な想定によって予め決定されていてもよい。
ての時間的に静止した成分を消去し、合成された運動目
標の軌跡の目印だけを残す。運動目標検出処理方法20の
特性はシナリオ依存性であり、合成および差計算動作に
よって伝播された個々のフレーム測定雑音、および残留
フレーム対フレームの整合エラーによる背景漏洩によっ
て劣化される。観察パラメータが許容すれば、許容可能
な性能は最適以下の方法に対して実現可能であり、それ
によって最適な整合されたフィルタ55および比較器60は
テンプレート整合されたフィルタおよび比較器60、スト
リーク検出器および比較器60、或は単なる単一画素また
は多重画素しきい値超過検出プロセッサによって置換さ
れてもよい。
43,351号明細書に記載された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による回転
開口運動目標検出処理方法20-2のトップレベルの実施形
態を示す。図3に示された検出方法20-2は、同じ情景を
観察することによって得られ、米国特許第 5,243,351号
明細書に記載された方法を使用して生成された2以上の
回転開口ストリップ開口映像30-2を使用する。照準線整
合エラーが存在する場合、個々の回転開口映像30-2のそ
れぞれは、米国特許第 5,243,351号明細書に記載された
ように空間的に最適に整合され、共通の空間周波数座標
系およびグリッドにフーリエ変換される。典型的に、共
通の空間周波数グリッドは、等価全円形開口回転開口合
成された空間周波数通過帯域と関連している。
のフレーム映像30-2のいくつかまたは全ての可能性のあ
る異なる対は、個々のフレームのそれぞれと関連し、か
つ回転開口センサシステム10の回転の結果生じた時間的
に変化するシステム伝達関数を補償するように設計され
た空間周波数フィルタ31-2によって乗算される。システ
ム伝達関数の補償は、ターゲット全円形開口またはその
他の基準システム伝達関数に各フレームを等しくする
か、差を計算されるべきフレームのシステム伝達関数に
1つのフレームを等しくするか、或は差を計算された両
映像からシステム伝達関数をデコンボルブすることを含
むいくつかの方法の1つによって行われてもよい。理想
的なシステム伝達関数に対して、またπラジアン離れて
捕捉されたフレームだけが差を計算される場合にのみ、
システム伝達関数の等化は不要である。
の可能な等化過程の任意の1つに基づいた回転開口運動
目標検出方法(20-2)は、特定の方法の受信機動作特性
(ROC)によって特徴付けられる特有の雑音特性を有
している。信号テンプレート方法に対するROCの差
は、等化プロセスによって発生された異なる補償されな
いレベルの雑音“ブースト”または“デブースト”によ
って支配される。しかしながら、適切に設計された雑音
に対して最適整合されたフィルタ方法は、等化プロセス
において雑音特性を補償する。
れ(40-2)、2つのフレームと関連した重畳した通過帯
域の外側の非共通情報を除去するために空間周波数濾波
される(41-2)。空間周波数濾波されたフレーム差は、
前に濾波されたフレーム差に加算され(35-2)、記憶さ
れ、それによって差を計算された空間周波数濾波された
フーリエ変換された映像の総和を生成する(35-2)。異
なった個々のフレーム映像の数がNである場合、可能性
のある異なった差を計算された個々のフレーム映像対の
総数は、N(N−1)/2の大きさである。
差を計算された空間周波数濾波され変換された映像の結
果的な和35-2は、空間ドメイン映像45-2を得るために逆
フーリエ変換され(42-2)、この空間ドメイン映像は雑
音最適整合されたフィルタ検出プロセッサ50-2に伝達さ
れる。整合されたフィルタ検出プロセッサ50-2は、形
状、寸法、速度、加速度、背景と目標のコントラスト
比、信号対雑音比、目標の数、軌道交差の数等のような
潜在的な目標の物理的現象に対応した特性を有する整合
されたフィルタ55-2のアレイを含む。整合されたフィル
タ検出プロセッサ50-2はまた等化、差計算および濾波処
理により伝播された測定雑音統計の過去の知識、および
回転開口センサシステム10の特有の信号測定特性を含ん
でいる。検出プロセッサ50-2の整合された各フィルタ55
-2は、差を計算された濾波された映像の総和35-2の全て
の画素位置に適用される。整合されたフィルタ55-2のそ
れぞれの出力は比較器60-2に伝送され、しきい値が各整
合されたフィルタ55-2に対して予め定められたレベルに
設定される。運動目標の検出は、整合されたフィルタ55
-2のアレイからの任意の出力が比較器60-2の予め定めら
れたしきい値を越えたときに決定される。しきい値は、
信号および雑音統計に関する演繹的な仮定によって、或
は個々のフレームまたは回転開口フレームのセットに基
づく信号および雑音統計に関する帰納的な仮定によって
予め決定されていてもよい。
レーム映像の共通の重複通過帯域中の時間的に静止した
全ての成分を消去し、運動目標の軌跡の目印だけを残
す。運動目標検出処理方法20-2の性能はシナリオ依存性
であり、等化、差計算および濾波動作によって伝播され
た個々のフレーム測定雑音、およびフレームの整合エラ
ーに対する残留フレームによる背景漏洩によって劣化さ
れる。観察パラメータが許すならば、最適以下の方法に
対して許容可能な性能が実現でき、それによって最適な
整合されたフィルタ55-2のアレイおよび比較器60-2はテ
ンプレート整合されたフィルタおよび比較器60-2、スト
リーク検出器および比較器60-2、或は単なる単一画素ま
たは多重画素しきい値超過検出プロセッサによって置換
されてもよい。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用する本発明の原理による回転
開口運動目標検出処理方法20-3のトップレベルの実施形
態を示す。図4に示された検出方法20-3は、同じ情景を
観察することによって得られ、米国特許第 5,243,351号
明細書に記載された方法を使用して生成された個々の回
転開口ストリップ開口映像(50-3)またはサブ映像の時
間シーケンスを使用する。照準線整合エラーが存在する
場合、個々の回転開口映像50-3のそれぞれは、米国特許
第 5,243,351号明細書に記載された技術によって空間的
に最適に整合され、共通の空間周波数座標系およびグリ
ッドにフーリエ変換され(51-3)、それによってフーリ
エ変換された映像を生成する(52-3)。典型的に、共通
の空間周波数グリッドは、等価全円形開口回転開口合成
された空間周波数通過帯域と関連している。いくつかま
たは全ての可能性のある空間的に整合されフーリエ変換
された個々のフレーム映像52-3は、時間的にフーリエ変
換され(53-3)、各空間周波数と全ての時間周波数を関
連させるデータセットを生成する(54-3)。
時間周波数および空間周波数変換されたデータセット
は、雑音に最適整合されたフィルタ検出プロセッサ55-3
に伝送される。整合されたフィルタ検出プロセッサ55-3
は、形状、寸法、速度、加速度、背景と目標のコントラ
スト比、信号対雑音比、目標の数、軌道交差の数等のよ
うな潜在的な目標の物理的現象に対応した特性を有する
整合されたフィルタ56-3のアレイを含む。整合されたフ
ィルタ検出プロセッサ55-3はまたシステム伝達関数のバ
リエーションおよび回転開口センサシステム10に関連し
たその他特有の測定特性の過去の知識、並びに測定およ
び変換処理により伝播された測定信号および雑音統計を
含んでいる。検出プロセッサ55-3の各整合されたフィル
タ56-3は、時間的および空間的に変換されたデータのセ
ット54-3に適用される。整合されたフィルタ56-3のそれ
ぞれの出力は比較器60-3に伝送され、しきい値が各整合
されたフィルタ56-3に対して予め定められたレベルに設
定される。運動目標の検出は、整合されたフィルタ56-3
のアレイからの任意の出力が比較器60-3の予め定められ
たしきい値を越えたときに決定される。しきい値は、信
号および雑音統計に関する演繹的な仮定によって、およ
び、または個々のフレームまたは回転開口フレームのセ
ットに基づく信号および雑音統計に関する帰納的な仮定
によって予め決定されてもよい。
々のフレーム映像中の全ての時間的に静止した成分を無
条件に消去し、運動目標の時間および空間周波数の軌跡
の目印だけを残す。運動目標検出処理方法20-3の性能は
シナリオ依存性であり、個々のフレームの測定雑音伝
播、およびフレームの整合エラーに対する残留フレーム
による背景漏洩によって劣化される。観察パラメータが
許すならば、最適以下の方法に対して許容可能な性能が
実現可能であり、それによって最適な整合されたフィル
タ56-3のアレイおよび比較器60-3はテンプレート整合さ
れたフィルタおよび比較器60-3、または簡単な単一周波
数“画素”または多周波数“画素”しきい値超過検出プ
ロセッサによって置換されてもよい。
記載された適切に修正された回転開口放射線センサシス
テムおよび処理方法を使用した本発明の原理による回転
開口運動目標検出処理方法20-4のトップレベルの実施形
態を示す。図5の検出方法20-4は、同じ情景を観察する
ことによって得られ、米国特許第 5,243,351号明細書に
記載された方法にしたがって生成された2以上の回転開
口ストリップ開口映像30-4を使用する。照準線整合エラ
ーが存在する場合、個々の回転開口映像のそれぞれは、
米国特許第 5,243,351号明細書に記載された方法によっ
て空間的に最適に整合され、共通の空間周波数座標系お
よびグリッドにフーリエ変換される。典型的に、共通の
空間周波数グリッドは、等価全円形開口回転開口合成さ
れた空間周波数通過帯域と関連している。
のフレーム映像30-4のいくつかまたは全ての可能性のあ
る異なる対は、個々のフレームのそれぞれと関連し、か
つ回転開口センサシステム10の回転の結果生じた時間的
に変化するシステム伝達関数を補償するように設計され
た空間周波数フィルタ31-4によって乗算される。システ
ム伝達関数の補償はいくつかの方法の1つによって行わ
れることができ、そのうちの幾つかの例は、ターゲット
全円形開口またはその他の基準システム伝達関数に各フ
レームを等化するか、差を計算されるフレームのシステ
ム伝達関数に1つのフレームを等化するか、或は差を計
算された両映像30-4からシステム伝達関数をデコンボル
ブすることを含んでいる。理想的なシステム伝達関数に
対して、またπラジアン離れて捕捉されたフレームだけ
が差を計算される場合には、システム伝達関数の等化は
不要である。
(40-4)、2つのフレームと関連した重なった通過帯域
の外側の非共通情報を除去するために空間周波数濾波さ
れる(41-4)。空間周波数濾波されたフレーム差は逆フ
ーリエ変換され(42-4)、個々の対の差の映像45-4を生
成する。異なる個々のフレーム映像の数がNである場
合、可能性のある異なる差を計算された個々のフレーム
映像対の総数は、N(N−1)/2であることができ
る。
の可能な等化手順の任意の1つに基づいた回転開口運動
目標検出方法20-4は、特定の方法の受信機動作特性(R
OC)によって特徴付けられる特有の雑音性能を有して
いる。信号テンプレート方法に対するROCの差は、等
化プロセスによって発生された異なる補償されないレベ
ルの雑音“ブースト”または“デブースト”によって支
配される。しかしながら、適切に設計された雑音の最適
整合されたフィルタ方法は、等化プロセスにおいて雑音
特性を補償する。
差を計算された空間周波数濾波され変換された映像45-4
のセットは、雑音の最適整合されたフィルタ検出プロセ
ッサ50-4に伝送される。整合されたフィルタ検出プロセ
ッサ50-4は、形状、寸法、速度、加速度、背景と目標の
コントラスト比、信号対雑音比、目標の数、軌道交差の
数等のような潜在的な目標の物理的現象に対応した特性
を有する個々の対の差のフレーム整合されたフィルタ55
-4のアレイを含む。整合されたフィルタ検出プロセッサ
50-4はまた等化、差計算および濾波処理を受けて伝播さ
れた測定雑音統計の演繹的知識、および回転開口センサ
システム10の特有の信号測定特性を含んでいる。検出プ
ロセッサ50-4の個々の対の差を計算されたフレーム整合
されたフィルタ55-4のそれぞれは、適切なフレーム差に
対して全ての画素位置に適用され、続いて加算される。
各目標の物理的現象の全体的に整合されたフィルタ56-4
の出力は、比較器60-4に伝送され、しきい値が全体的に
整合されたフィルタ56-4のそれぞれの予め定められたレ
ベルに設定される。運動目標の検出は、全体的に整合さ
れたフィルタ55-4のアレイからの任意の出力が比較器60
-4の予め定められたしきい値を越えたときに決定され
る。しきい値は、信号および雑音統計に関する演繹的な
仮定によって、或は個々のフレームまたは回転開口フレ
ームのセットに基づく信号および雑音統計に関する帰納
的な仮定によって予め決定されていてもよい。
レーム映像の共通の重なった通過帯域中の時間的に静止
した全ての成分を消去し、運動目標の軌跡の目印だけを
残す。運動目標検出処理方法20-4の性能はシナリオ依存
性であり、等化、差計算および濾波動作を受けて伝播さ
れた個々のフレーム測定雑音、およびフレームの整合エ
ラーに対する残留フレームによる背景漏洩によって劣化
される。観察パラメータが許すならば、最適以下の方法
に対して許容可能な性能が達成されることができ、それ
によって最適な整合されたフィルタ56-4のアレイおよび
比較器60-4はテンプレート整合されたフィルタおよび比
較器60-4、ストリーク検出器および比較器60-4、或は単
なる単一画素または多重画素しきい値超過検出プロセッ
サによって置換されてもよい。
標検出処理方法20は以下の式により表される。米国特許
第 5,243,351号明細書の定義および用語を使用すると、
本発明において使用される強化された運動目標検出回転
開口センサ10およびデータ処理方法20,20-2,20-3,20
-4の数学的な基礎の説明が与えられる。
する。
中に得られた位置/Rに配置された検出器の積分期間中
に生成された検出電子の平均数を表わしたものである。
ここで、Os (/F,tj )は空間周波数/Fおよびフ
レーム時間tj における検出器出力の空間フーリエ変換
である。実効センサシステム伝達関数が検出器積分時間
ΔT中時間的に静止していると仮定すると、強度を検出
された個々のフレームのフーリエ変換は次のように表わ
される
たって空間的に平均された強度であり、《W0 》は検出
器伝達関数のゼロ空間周波数成分であり、AfpaSは実効
焦点平面アレイ面積であり、STF(/F,tj )は、
フレーム時間tj に対応した全体的なセンサシステム伝
達関数であり、《Sn 》(M/F,t)は時間tにおけ
る情景の幾何学的な投影の正規化された空間フーリエ変
換である。
分を含む情景および映像生成平面への情景の幾何学的投
影の数学的表示を説明する。背景の前方(回転開口セン
サに関して)および透過散乱層の後方を運動している単
一の目標を含む情景から放出された強度は次のように表
わされる。
時間tにおける運動目標の強度分布であり、B(/r,
t)は情景位置/rおよび時間tにおける背景の強度分
布であり、H(/r,t)は情景位置/rおよび時間t
に対応した透過散乱層の強度分布であり、WT (/r,
t)は情景位置/rおよび時間tにおける目標の逆シル
エットである。目標の逆シルエットは次の式によって与
えられる。
Wt (/r,t)を減算することは、情景位置/rおよ
び時間tにおいて運動目標によってぼかされ背景光子を
除去することに対応する。
の情景の強度の幾何学的な投影は、次の式によって表わ
される
な投影をフーリエ変換して、フーリエ変換類似性の定理
を適用することにより次式が得られる。
および変数の変換は、
成分の対応した幾何学的投影は、
スペクトルは、次のように表わされる。
前景スペクトル並びにスペクトルに対する運動目標の貢
献によって支配される2つの成分に分解されることが可
能であり、これは次のように表わされ、
ける空間周波数M(/F)に対する正規化された背景お
よび前景成分であり、
時間tにおける空間周波数/Fに対する背景および前景
からの正規化された目標のバリエーションである。
化された背景および前景スペクトルを説明する。時間的
に静止した全体的な背景成分を有する観察シナリオに対
して、正規化された全体的な背景スペクトル、および運
動目標バリエーションは次のようになる。
る。運動目標によってぼかされる背景光子の数が、合計
背景光子数に比較して無視できるほど僅かである場合、
合計背景および運動目標バリエーションの正規化された
スペクトルはさらに簡単化されることができ、次式を生
成する。
目標の逆シルエットを説明する。目標空間強度分布の論
証に時間依存性の空間シフトを適用することによって目
標の時間依存性運動を説明することができる。映像生成
プロセスの開始時間をt0であると定めると、最初の目
標および目標逆シルエットは次式によって与えられる。
間バリエーションは次式によって与えられる。
情景座標における目標のベクトル位置であり、目標軌道
と呼ぶ。
トルを説明する。フーリエ変換シフト定理を運動目標お
よび運動目標逆シルエットに適用すると、以下の式によ
ってそれぞれ与えられる対応したスペクトルが生成され
る。
た背景からの運動目標バリエーションを説明する。
《T》{M(/F),t}および《WT 》{M(/
F),t}に対する上記の式を目標スペクトルバリエー
ションΔ/TN の式に代入すると、次式が得られる。
て、目標のベクトル位置は/RT (t)=/Vtとして
表わされることができ、ここで/Vは運動目標の速度成
分である。したがって、正規化された目標スペクトルバ
リエーションΔ《Tn 》は次のように表わされる。
クトルを説明する。映像化された運動目標および全体的
な背景の時間平均の結果、平均検出器測定が得られるた
め、システム伝達関数が検出積分時間中一定であると仮
定した測定の時間平均されたフーリエスペクトルは次の
ようになる。
定速度の運動の時間平均は、
セスを説明する。回転開口合成プロセスの最も簡単なバ
ージョンは、無雑音測定の空間周波数冗長平均が等価全
円形開口映像を生成するようにシンセサイザを使用す
る。すなわち、
べき全円形開口映像であり、《Λ》(/F)はシンセサ
イザである。また《OS ave 》(/F)は個々の(無雑
音)フレーム測定スペクトルの空間周波数冗長平均であ
る。
し、それは光学系伝達関数OTFと、正規化された検出
器伝達関数《WN 》との積であり、下付きの文字Sおよ
びFは、回転開口ストリップ開口センサおよびターゲッ
ト全円形開口システムと関連したパラメータをそれぞれ
示す。シンセサイザの分母に現れる空間周波数冗長平均
されたストリップ開口伝達関数は、
間的に静止した全ての背景および前景からなる情景の観
察から結果的に得られ、センサの照準線における運動目
標が次式によって与えられる。
じ情景を観察している合成された映像の任意の対に対し
て同じなので、合成された映像の差は運動目標情報だけ
を含んでいる。差を計算された合成された映像の無雑音
空間周波数ドメインの運動目標目印は、次の式によって
与えられる。
空間目印は、
レームの異なるセットの獲得に関連した時間遅延を表わ
す。
めに適切に修正された回転開口フレームの任意の対の差
の計算について説明する。2Nf 個の個々の回転開口ス
トリップ開口フレームに対して、2Nf (Nf −1)/
2の可能な異なった対でフレーム比較、すなわち差が存
在している。最適な検出技術は、任意の可能な対および
潜在的に非常に多数の運動目標の対となった目印の差を
計算するときに背景除去を利用する。しかしながら、時
間的に静止した開口システムに対して背景除去および運
動目標検出のためにフレームの差を計算するのとは異な
り、特定の技術および方法が、システム伝達関数の回転
特性を補償するために回転開口データセットに適用され
なければならない。回転開口データセットは、情景にお
ける時間変化と運動目標から派生した時間的変化とを区
別するように修正される。
よびフレーム毎のシステム伝達関数等化を説明する。個
々のストリップ開口雑音測定スペクトルは、無雑音スト
リップ開口検出映像スペクトル、プラスストリップ開口
の加算された雑音スペクトルとして表わされることがで
きる。
無雑音信号および雑音スペクトルである。
過帯域フィルタを説明する。一般に、個々のフレーム測
定は、等価全円形開口光学系通過帯域の異なった部分に
関する情報を提供する。運動目標または過渡事象のよう
な時間的に静止していない現象を検出するには、フレー
ム間に共通した情報を識別して使用し、また回転開口セ
ンサの伝達関数の時間的に変化する特性による映像コン
トラストの差に対して調節を行う必要がある。フレーム
毎のシステム伝達関数補償は、いくつかの方法によって
行われることができ、そのうちの数例は、ターゲット全
円形開口またはその他の基準システム伝達関数にそれぞ
れのフレームを等しくするもの、差を計算されるフレー
ムのシステム伝達関数に1つのフレームを等しするこ
と、差を計算される両方の映像からシステム伝達関数を
デコンボルブすること、或は理想的なシステム伝達関数
に等しくすること等であり、またπラジアン離れて獲得
されたフレームだけが差を計算された場合にシステム伝
達関数の等化は不要である。テンプレート信号検出方法
に対して、多数の可能な等化過程の任意の1つに基づい
た回転開口運動目標検出方法は、特定の方法の受信機動
作特性(ROC)によって特徴付けられる特有の雑音特
性を有している。信号テンプレート方法に対するROC
における差は、等化プロセスによって発生された異なる
補償されないレベルの雑音“ブースト”または“デブー
スト”によって支配される。しかしながら、適切に設計
された雑音の最適整合されたフィルタ方法は、等化プロ
セスにおいて雑音特性を補償する。
の方法の任意の1つを使用することができるが、以下の
段落において説明される式は、各過程と関連した数学的
詳細の典型である。簡略化するために、1つのフレーム
が別のフレームに等しくされる過程だけを詳細に説明す
る。
ムと、比較のために使用される任意の異なる別のフレー
ムとを選択すると、フレームのスペクトルは、両方のフ
レームの通過帯域の共通の重なった領域における空間周
波数情報を提供するように濾波された空間周波数とな
る。この別のフレームはまたシステム伝達関数の差を補
償するために修正される。すなわち、
ムの出力を濾波された空間周波数を示し、《mE 》
(F,ti )は比較されるべきフレームの出力を濾波さ
れ、伝達関数等化された空間周波数を示し、《x》は光
学系通過帯域フィルタであり、
あり、
た差は、差を計算された信号および雑音スペクトルの和
として表されることができる。すなわち、
た雑音スペクトル《ΔNij》は、
信号スペクトル《ΔSij》(/F)は、
結果的な信号スペクトルは、以下の式によって背景から
の正規化された運動目標バリエーション間の差Δ
《TN 》に関連させられる:
ている場合、濾波された差の信号は以下の関係に対応す
る:
ドメインの運動目標の検出を説明する。空間周波数およ
び時間周波数ドメイン回転開口運動目標の検出過程の簡
単な数学的記述は、ディスクリートな時間サンプリング
を連続体として近似し、個々の回転開口フレームのそれ
ぞれに対して空間周波数変換の連続したフーリエ変換を
行うことによって得られる。個々の無雑音フレーム測定
のフーリエスペクトルは以下の式によって与えられる:
間フーリエ変換は、積分を評価してコンボリューション
定理を適用することによって以下のようにして得られ
る。
達関数の時間フーリエ変換であり、《SN 》は正規化さ
れた幾何学的に投影された情景スペクトルの時間および
空間周波数フーリエ変換である。
ペクトルの時間フーリエ変換を評価すると、ストリップ
開口測定の無雑音時間周波数変換は次のようになる。
の等式は、時間および空間周波数整合されたフィルタの
信号成分を決定するために使用される。{OS }(F,
f)に対する式の最初の項は時間的に静止した結果を全
て含んでおり、かつ全てのゼロでない空間周波数に対し
てゼロであるため、整合されたフィルタは本質的に時間
的に静止した成分を消去する。
対する簡単化された目印を説明する。回転開口運動目標
の検出過程に対する受信機動作特性の予測は、均一な強
度の背景、前景および運動目標に関連した観察シナリオ
の特殊な場合に対して簡単化されることができる。事
実、粗運動目標検出方法の設計はこのような簡単化され
た仮定を使用することができる。背景および前景の時間
的に静止しており、かつ空間的に均一な成分は次のよう
に表されることができる。
強度であり、HSCN は情景座標における平均前景強度で
ある。
を説明する。情景の背景および前景強度の平均値の幾何
学的な投影は、
ようになる。
うなフーリエ変換が行われる。:
ている。
を説明する。情景スペクトルの正規化された幾何学的な
投影は、幾何学的に投影された情景スペクトルのゼロ空
間周波数成分で情景の幾何学的投影のフーリエ変換を除
算することによって生成される。背景の強度および運動
目標の強度の幾何学的に投影されたゼロ空間周波数成分
は、空間的に平均された平均値と、背景および目標の幾
何学的な投影された面積に関して以下のように表される
ことができる:
な幾何学的投影面積であり、Afpaは回転開口センサの
焦点平面における背景のおおよその検出された面積であ
り、文字上の符号(−)は、焦点平面アレイにおける空
間平均を示す。
周波数関係を説明する。情景座標系における目標逆シル
エットにフーリエ変換の定義を適用し、Mrscene =R
image により与えられる変数の変化によって映像平面座
標に変換すると、以下の式が得られる:
の幾何学的な投影はまた、以下の式によって焦点平面に
おける目標の実効面積に関連されることができる:
トルが以下に説明される。情景の均一な背景成分に対し
て前述の関係を適用することによって、正規化された背
景および正規化された目標の数学的記述を決定すること
ができる。背景および前景の投影情景スペクトルの合計
に対するそれらのゼロ空間周波数成分の比率は以下の通
りである。
な背景と、正規化された目標との差《TN 》は、以下の
ように表される。
の座標において、均一な目標は、目標の逆シルエットW
T と目標の一定の強度レベルTO との積として表される
ことができる。
され、
された変換は以下のように表わされる。
間的に静止した情景の、正規化された目標の背景からの
変化Δ《TN 》は、以下の式によって与えられ、
以下のように定義され、
間周波数成分は、以下のように表される。
めにフーリエ変換シフト定理を適用すると、目標の逆シ
ルエットは次式のようになる。
《WTN》(M/F,t0 )は、時間=t0 における目標
の正規化された逆シルエットのフーリエ変換である。結
果的に、背景からの目標の変化の式は以下のように表さ
れる。
に対する、差を計算され濾波された運動目標の空間的周
波数の目印は次式によって表される。
ンの差が計算された信号が以下に説明される。運動目標
の差が計算された濾波された空間周波数スペクトルを逆
フーリエ変換すると、空間ドメインの運動目標の目印は
以下のように表される。
音の比率が以下に説明される。均一な目標および背景に
対する位置/Rにおける平均差信号対RMS差雑音の比
率は次式によって表される。
成された対の信号スペクトルが表される。均一な目標お
よび背景と関連した任意の対の差の合成された映像の雑
音のない空間周波数ドメインの目印は、次式によって表
され、
空間ドメインの運動目標は、以下のように表される。
される。幅広いクラスの性能基準に対して最適に動作す
る検出器プロセッサが、以下、比較器と呼ばれる比較装
置を有する一般化された整合されたフィルタのアレイに
構成される。スピンAp運動目標検出プロセッサのため
の最適に整合されたフィルタアレイの設計には、スピン
Apセンサの特有の信号測定特性と、整合されたフィル
タアレイに入力される前の所望されたデータ処理と、予
測された目標の現象学とに関する知識を組入れることが
要求される。過去の目標知識は、本来統計的あるいは決
定論的であり、形状、寸法、速度、加速度、背景対目標
のコントラスト比、信号対雑音比、目標の数、軌道交差
の数、雑音の交差共分散等の物理的パラメータを含んで
いる。特定の目印に対する整合されたフィルタのアレイ
は、行列法で以下のように表されることができ、 F=St CNN -1 ここにおいて、Sは過去の雑音のない信号のベクトルで
あり、CNN -1は雑音の交差共分散逆マトリックスであ
る。信号および雑音の交差共分散マトリックスの両方
が、所望された運動目標検出データ処理方法に適切に定
義される。
周波数のバリエーション、あるいはスピンApによる運
動目標の検出に適切となるように定義された別の目印を
任意に組み合わせることによって生じた信号マトリック
スは以下のように表され、
は決定論的なスピンAp運動目標目印と関連した独立し
た空間、空間周波数、時間、時間周波数の信号の座標の
任意の全ての組み合わせを示している。
像、時間周波数の変化、あるいはスピンApによる運動
目標の検出に適切となるように定義された別の適切な目
印を任意に組み合わせることによって生じた雑音マトリ
ックスは以下のように表される。
および雑音成分の合計として行列法で以下のように表さ
れる。
び雑音のマトリックス成分がそれぞれ複素数の信号およ
び雑音の実数および虚数部分(すなわち、振幅および位
相部分)であることを認識することによって構成される
ことができる。
ィルタの捕捉“ウィンドウ”内の任意の数の差のフレー
ムあるいは差の合成された映像の信号および雑音マトリ
ックスは、それぞれ以下のように表され、
たフィルタの捕捉ウィンドウにおける全ての成分の数を
示している。
されたフィルタの出力は、整合されたフィルタの雑音の
ない信号および付加的な“白色”雑音に対応し、次式に
よって表される。
特性が以下に説明される。比較器およびしきい値超過決
定方法を使用するスピンApの整合されたフィルタによ
る運動目標検出処理方法の受信機の動作特性(ROC)
は、個々のフレーム測定の近似のガウス雑音を求めるた
めに容易に決定されることができる。所定のしきい値の
スピンAp整合フィルタの出力に対する検出確率および
誤警報確率は次式によってそれぞれ決定され、
パラメータである。
ックス成分が説明される。雑音が最適なスピンAp運動
目標検出整合フィルタの決定は、雑音の処理および測定
ステップから整合されたフィルタへの雑音の伝播に関す
る過去の知識あるいは仮定に基づく。本開示において説
明されたスピンAp目標検出方法のクラスは、個々の加
重されたフレームの差、合成されたフレームの差、ある
いは時間および空間周波数整合されたフィルタを使用す
る。適切に伝播された雑音の交差共分散マトリックス成
分の決定には、雑音の最適な整合されたフィルタを構成
することあるいは雑音最適化技術以下の性能が必要とさ
れる。以下のセクションによって、適切なスピンAp雑
音共分散マトリックス成分を決定するための例示的な技
術が提供される。はっきりとは説明されていない処理ス
テップに関連した共分散マトリックス成分は、類似した
決定の方法に従う。
が説明される。スピンAp運動目標検出方法の一実施形
態によって、等化された異なる対のシステム伝達関数
と、オーバーラップした通過帯域で濾波された個々のフ
レームとの差が計算される。対応する差が計算され、等
化され、濾波された雑音スペクトルNijは次式によって
表される。
処理された雑音スペクトルの逆フーリエ変換を行うこと
によって以下のように得られる。
下のようになる。
位置(/Rm )および(/Rn )における任意の2対の
差が計算され、等化され、濾波されたフレーム間の空間
雑音交差共分散は以下の次式で表され、
表している。
雑音の効果が説明される。ゼロ平均で、時間的に静止
し、時間的に相関しない雑音に対して、個々のフレーム
の雑音スペクトルの交差共分散は、以下の式
付けられ、ここにおいて、δi,j およびδ(/F−/
F)は、それぞれクロネッカーおよびディラックのδ関
数である。
分散を空間的に差が計算された濾波された雑音の共分散
の式に代入すると、以下のような式が生成される。
レームの差の画素毎の交差共分散は、フレーム指数がi
=k、j=l、およびi≠jの値をとるようにすること
によって得られる。対応する実効雑音空間的交差共分散
は、以下の式によって示される。
ングについて説明する。空間的白色雑音および光学系の
通過帯域のナイキスト・サンプリングに対して、個々の
スピンApのフレームの雑音のスペクトル密度は以下の
ように近似され、
ステムの実効焦点距離であり、DxおよびDy はそれぞ
れxおよびy方向の開口の寸法であり、Δxs およびΔ
ysは空間的サンプリングのインクレメントであり、そ
れは、検出器の間隔に対応する100%の充填率を有す
る焦点表面検出器アレイに対するものであり、/σds 2
は個々の検出器における雑音の平均的分散である。個々
の白色雑音のスペクトル密度の交差共分散を、差が計算
され、濾波された空間の雑音の共分散の式に代入する
と、以下の式が生成される。
る画素から画素の交差共分散に対応する特別な場合は、
よび白色雑音について説明する。理想的なシステム伝達
関数に対して、システム伝達関数の等化を必要としない
という処理上の利点を有するさらに簡単にする1例は、
πラジアン離れたフレームの差を計算することに対応す
る。理想的な方形ストリップ開口に対して、πラジアン
離れた2つのフレームに対するシステム伝達関数は、同
じコントラスト性能を提供し、したがって等化フィルタ
は正確に1であり、すなわち《H》(/F,tj ,
ti )=1であり、重なった通過帯域は同一である。し
たがって、空間的に白色で時間的に相関しない雑音に対
して、交差共分散雑音マトリクス成分は、
されることができ、
故、共分散マトリクスの逆数は、
る。回転開口運動目標検出処理方法の別の実施形態は、
回転開口合成された映像の任意の異なった対の差を計算
する。米国特許第 5,243,351号明細書に記載されている
ように、最も簡単な形態の合成プロセスに対する合成さ
れた雑音スペクトルは、
S ave 》(/F)は次式によって与えられる:
合成された雑音スペクトルを与えると、雑音スペクトル
Δ《nsyn 》(/F)は、合成された画像の任意の対
(完全に異なったデータセットから得る必要はない)で
は次式のようになる。
成するために使用されたスピンAPの個々のフレームセ
ットを意味する。
明する。差の合成された雑音スペクトルの逆フーリエ変
換を行うと、画素位置(/R)で空間的ドメインの合成
された雑音が発生する。
散について説明する。合成された雑音の差の画素毎の交
差共分散は、スピンApシンセサイザと冗長度平均され
た差の雑音スペクトルの交差共分散に関して表されるこ
とができ、即ち、
は以下のようになる。
れていない雑音を説明する。ゼロ平均に対して時間的に
静止している相関されていない雑音では、個々のフレー
ム雑音スペクトルの交差共分散は次式により雑音パワー
スペクトル密度《CNN》に関連される。
れ、共通のデータを有するスピンApのセットの個々の
フレームに基づいて合成画像のポテンシャルを適応す
る。
関について説明する。異なったセットのスピンApの個
々のフレームから生成される合成された画像に対して
は、差の合成された画像についての雑音スペクトルの交
差共分散は以下の式で表される。
雑音に対する画素と画素の差の合成された雑音は次式で
与えられる。
空間的な白色の空間的な差の合成された雑音の交差共分
散マトリックスは以下の式で与えられる。
た本発明のスピンAPの運動目標システムの実施形態が
実施され試験された。画像データは図2の実施形態を用
いて生成され、したがって雑音のない差の合成された運
動目標の目印と、対応するテンプレート整合フィルタ出
力を生成する。
む低い信号対雑音比の差の合成された画像と、テンプレ
ート整合フィルタの出力を示し、したがって、16対1
程度の優れた出力の信号対雑音比を示す。このことはこ
の特定のクラスの運動目標に対して本発明を使用して生
成される強化された受信機の動作特性を示している。さ
らに、実験データが得られ、ここで説明したπ方向の運
動目標検出方法を証明するために使用された。
し、複数の回転ストリップ開口画像測定または合成され
た画像測定から運動目標を検出する回転ストリップ(部
分的)開口画像放射線センサを説明した。前述の実施形
態は本発明の原理の応用を表した多数の特定の実施形態
のうちの数例にすぎないことが理解されよう。明白に、
多数および変更された他の装置は本発明の技術的範囲を
逸脱することなく当業者により容易に製造されよう。
を使用する回転開口映像生成放射線センサシステムの概
略図。
た映像を使用する本発明の原理による回転開口運動目標
検出処理方法のトップレベルの実施形態の概略図。
化された個別の回転開口変換されたフレームを使用する
本発明の原理による回転開口運動目標検出処理方法のト
ップレベルの実施形態の概略図。
レーム映像の時間および空間周波数変換を使用する本発
明の原理による回転開口運動目標検出処理方法のトップ
レベルの実施形態の概略図。
化された個別の回転開口変換されたフレーム、および順
次差フレーム整合された濾波を使用した本発明の原理に
よる回転開口運動目標検出処理方法のトップレベルの実
施形態の概略図。
Claims (3)
- 【請求項1】 光軸を中心に回動し、時間的に順次連続
する情景の映像を生成する回転ストリップ開口を具備す
る望遠鏡と、 この望遠鏡の焦点平面に形成された映像を検出するため
の2次元の検出器アレイと、 この検出器アレイに結合され、ストリップ開口が望遠鏡
の光軸を中心に回動するとき映像情景の複数の映像フレ
ームを記録し、それらの複数の映像フレームからの映像
を合成し、合成された映像中に発見された運動目標を検
出する信号プロセッサとを具備していることを特徴とす
る回転ストリップ開口放射線センサシステム。 - 【請求項2】 信号プロセッサは、同じ情景を観察する
ことによって得られた2以上の合成された映像を処理
し、 合成された映像の全ての異なる対の差を生成し、 差の映像を合計し、 整合されたフィルタアレイと比較器とにより差を生成さ
れ、合計されたデータを処理して整合されたフィルタア
レイの出力が比較器の予め定められたしきい値を越えた
とき運動目標の検出を決定する処理装置によって合成さ
れた映像中に発見された運動目標を検出する請求項1記
載のシステム。 - 【請求項3】 信号プロセッサは、同じ情景を観察する
ことによって得られた2以上の個々の照準線補正フレー
ムを処理し、 システム伝達関数等化フィルタを使用して個々のフレー
ムを濾波し、 等化されたフレームの異なる対の差を生成し、 共通でない通過帯域排除フィルタを使用して差のデータ
を濾波し、 濾波された差のデータを合計して濾波された差の全体の
合計されたスペクトルを形成し、 合計されたスペクトルを逆変換して空間的ドメインにお
いて運動目標の目印を与え、 整合されたフィルタアレイと比較器とを通過することに
よって得られた合計されたデータを処理し、 整合されたフィルタアレイの出力が比較器の予め定めら
れたしきい値を越えたとき運動目標の検出を決定する請
求項1記載のシステム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US633980 | 1996-04-19 | ||
| US08/633,980 US6192322B1 (en) | 1996-04-19 | 1996-04-19 | Moving object and transient event detection using rotation strip aperture image measurements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1054879A true JPH1054879A (ja) | 1998-02-24 |
| JP2889210B2 JP2889210B2 (ja) | 1999-05-10 |
Family
ID=24541959
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP9117620A Expired - Lifetime JP2889210B2 (ja) | 1996-04-19 | 1997-04-21 | 回転ストリップ開口映像測定を使用する運動物体および過渡事象検出 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6192322B1 (ja) |
| EP (1) | EP0802426B1 (ja) |
| JP (1) | JP2889210B2 (ja) |
| DE (1) | DE69708195T2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018036260A (ja) * | 2016-08-24 | 2018-03-08 | エクセリタス テクノロジーズ シンガポール プライヴェート リミテッド | モデル化された背景差分による赤外線在感知 |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6142681A (en) | 1999-02-22 | 2000-11-07 | Vialogy Corporation | Method and apparatus for interpreting hybridized bioelectronic DNA microarray patterns using self-scaling convergent reverberant dynamics |
| US6136541A (en) | 1999-02-22 | 2000-10-24 | Vialogy Corporation | Method and apparatus for analyzing hybridized biochip patterns using resonance interactions employing quantum expressor functions |
| US7215795B2 (en) * | 2000-09-28 | 2007-05-08 | Hitachi Kokusai Electric Inc. | Intruding object detecting method and intruding object monitoring apparatus employing the method |
| JP4337266B2 (ja) * | 2001-01-31 | 2009-09-30 | コニカミノルタセンシング株式会社 | 3次元測定方法および3次元測定システム |
| US6707418B2 (en) * | 2002-07-12 | 2004-03-16 | Itt Manufacturing Enterprises, Inc. | Correlating single continuous wave radar space object observation for initial orbit determination |
| ES2392512T3 (es) * | 2004-01-27 | 2012-12-11 | Saab Ab | Método para la reducción del desenfoque angular en imágenes de radar |
| EP1792263A2 (en) | 2004-09-02 | 2007-06-06 | Vialogy Corporation | Detecting events of interest using quantum resonance interferometry |
| US8022864B2 (en) * | 2007-11-08 | 2011-09-20 | The United States Of America As Represented By The Secretary Of Commerce | Detection of transient signals in doppler spectra |
| CN101261319B (zh) * | 2008-04-18 | 2010-12-01 | 华中科技大学 | 一种综合孔径辐射计成像校正方法 |
| SG172972A1 (en) * | 2009-01-28 | 2011-08-29 | Bae Systems Plc | Detecting potential changed objects in images |
| US8254650B2 (en) * | 2009-02-27 | 2012-08-28 | General Electric Company | System and method for contrast enhancement of time-resolved fluorescence images |
| US8325981B2 (en) * | 2009-04-21 | 2012-12-04 | Nec Soft, Ltd. | Human tracking apparatus, human tracking method, and human tracking processing program |
| CN102062572B (zh) * | 2009-11-13 | 2012-07-11 | 中国科学院西安光学精密机械研究所 | 基于jtc的高精度光电混合像移测量装置及其方法 |
| US8760340B2 (en) * | 2010-07-08 | 2014-06-24 | Raytheon Company | Processing radar return signals to detect targets |
| RU2461067C1 (ru) * | 2011-05-04 | 2012-09-10 | Государственное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ГОУ ВПО ТГТУ | Способ обнаружения движущихся объектов |
| US9228895B2 (en) * | 2012-05-25 | 2016-01-05 | Raytheon Company | Imaging system with multiple focal plane array sensors |
| US9430846B2 (en) | 2013-04-19 | 2016-08-30 | Ge Aviation Systems Llc | Method of tracking objects using hyperspectral imagery |
| KR102081087B1 (ko) | 2013-06-17 | 2020-02-25 | 삼성전자주식회사 | 동기적 영상과 비동기적 영상을 위한 영상 정합 장치 및 이미지 센서 |
| CN103684605B (zh) * | 2013-11-29 | 2016-05-04 | 西安空间无线电技术研究所 | 一种基于微波光子的综合孔径辐射计多路信号传输系统 |
| JP2015186170A (ja) * | 2014-03-26 | 2015-10-22 | ソニー株式会社 | 画像処理装置および画像処理方法 |
| DE102015120967A1 (de) | 2015-12-02 | 2017-06-08 | Carl Zeiss Ag | Verfahren und Vorrichtung zur Bildkorrektur |
| CN106156769A (zh) * | 2016-07-14 | 2016-11-23 | 上海摩象网络科技有限公司 | 基于特征频率光学系统的多目标识别技术 |
| DE102017112484A1 (de) | 2017-06-07 | 2018-12-13 | Carl Zeiss Ag | Verfahren und Vorrichtung zur Bildkorrektur |
| CN107817053B (zh) * | 2017-09-28 | 2019-08-23 | 中国科学院长春光学精密机械与物理研究所 | 一种大口径红外辐射测量系统定标方法及装置 |
| AU2023397402A1 (en) * | 2022-12-14 | 2025-07-17 | Western Sydney University | Systems for imaging a target object moving relative to background features |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0181935B1 (en) * | 1984-05-07 | 1989-07-12 | Hughes Aircraft Company | Microwave radiometer using fanbeam inversion |
| US4724439A (en) * | 1984-05-07 | 1988-02-09 | Hughes Aircraft Company | Microwave radiometer using fanbeam inversion |
| US4739332A (en) * | 1984-09-17 | 1988-04-19 | Hughes Aircraft Company | Fanbeam inversion radar |
| US5072226A (en) * | 1990-06-07 | 1991-12-10 | Hughes Aircraft Company | Radiometer system incorporating a cylindrical parabolic reflector and minimum redundancy array feed |
| US5113194A (en) * | 1990-08-29 | 1992-05-12 | Hughes Aircraft Company | Acceleration compensation by matched filtering |
| FR2670978A1 (fr) * | 1990-12-21 | 1992-06-26 | Philips Electronique Lab | Procede d'extraction de mouvement comprenant la formation d'images de differences et un filtrage tridimensionnel. |
| EP0497586A3 (en) * | 1991-01-31 | 1994-05-18 | Sony Corp | Motion detection circuit |
| US5243351A (en) * | 1992-06-25 | 1993-09-07 | Hughes Aircraft Company | Full aperture image synthesis using rotating strip aperture image measurements |
| US5350911A (en) * | 1993-04-09 | 1994-09-27 | Hughes Aircraft Company | Wavefront error estimation derived from observation of arbitrary unknown extended scenes |
-
1996
- 1996-04-19 US US08/633,980 patent/US6192322B1/en not_active Expired - Lifetime
-
1997
- 1997-04-19 EP EP97106511A patent/EP0802426B1/en not_active Expired - Lifetime
- 1997-04-19 DE DE69708195T patent/DE69708195T2/de not_active Expired - Lifetime
- 1997-04-21 JP JP9117620A patent/JP2889210B2/ja not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018036260A (ja) * | 2016-08-24 | 2018-03-08 | エクセリタス テクノロジーズ シンガポール プライヴェート リミテッド | モデル化された背景差分による赤外線在感知 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2889210B2 (ja) | 1999-05-10 |
| EP0802426A3 (en) | 1998-01-07 |
| US6192322B1 (en) | 2001-02-20 |
| DE69708195D1 (de) | 2001-12-20 |
| EP0802426A2 (en) | 1997-10-22 |
| EP0802426B1 (en) | 2001-11-14 |
| DE69708195T2 (de) | 2002-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2889210B2 (ja) | 回転ストリップ開口映像測定を使用する運動物体および過渡事象検出 | |
| US5243351A (en) | Full aperture image synthesis using rotating strip aperture image measurements | |
| JP2705911B2 (ja) | クラッター測定値と信号対クラッター比を使用したイメージベースの検出および追跡システム並びに処理方法 | |
| JP2547948B2 (ja) | ミサイル監視方法及び装置 | |
| Smith et al. | Understanding image quality losses due to smear in high-resolution remote sensing imaging systems | |
| US20100067822A1 (en) | System and method of super-resolution imaging from a sequence of translated and rotated low-resolution images | |
| Roggemann | Optical performance of fully and partially compensated adaptive optics systems using least-squares and minimum variance phase reconstructors | |
| US20110267486A1 (en) | Range measurement using multiple coded apertures | |
| CN114170192B (zh) | 星载光学相机焦面震颤精细检测方法 | |
| EP0802425A2 (en) | Full aperture scene synthesis using rotating strip aperture image measurements | |
| AU2020245706A1 (en) | Method and system to detect orbiting objects | |
| CN112927177A (zh) | 微小卫星平台上的多帧图像快速重建方法 | |
| Zhaoxiang et al. | Attitude jitter compensation for remote sensing images using convolutional neural network | |
| Sutton et al. | Correcting spacecraft jitter in HiRISE images | |
| CN120405701B (zh) | 一种空间目标多维度探测装置和方法 | |
| Baranova et al. | The resident space object detection method based on the connection between the Fourier domain image of the video data difference frame and the orbital velocity projection | |
| Prabu et al. | Demonstration of orbit determination for LEO objects using the Murchison Widefield Array | |
| US8351738B2 (en) | Method of estimating at least one deformation of the wave front of an optical system or of an object observed by the optical system and associated device | |
| Lambert et al. | Superresolution in image restoration of wide-area images viewed through atmospheric turbulence | |
| Hilkert et al. | Specifications for image stabilization systems | |
| Eineder et al. | A flexible system for the generation of interferometric SAR products | |
| Wang et al. | DeepPSF: On-Orbit Point Spread Function Estimation of Space Camera with Deep Learning | |
| CN119803450B (zh) | 一种基于空间滤波测速的目标初定轨方法及装置 | |
| Muruganandan et al. | Improving the resolution of low earth orbit objects by multi-exposure imaging and deconvolution | |
| Bayrak | Video stabilization: digital and mechanical approaches |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080219 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090219 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100219 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100219 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110219 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120219 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120219 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130219 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130219 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140219 Year of fee payment: 15 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |