JPH1063821A - 回転ストリップ開口映像測定を使用する超解像度全円形開口情景合成 - Google Patents
回転ストリップ開口映像測定を使用する超解像度全円形開口情景合成Info
- Publication number
- JPH1063821A JPH1063821A JP9117621A JP11762197A JPH1063821A JP H1063821 A JPH1063821 A JP H1063821A JP 9117621 A JP9117621 A JP 9117621A JP 11762197 A JP11762197 A JP 11762197A JP H1063821 A JPH1063821 A JP H1063821A
- Authority
- JP
- Japan
- Prior art keywords
- scene
- super
- aperture
- resolution
- equation
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0871—Complete apparatus or systems; circuits, e.g. receivers or amplifiers
-
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/78—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic or electromagnetic waves, or particle emission, not having a directional significance, are being received using electromagnetic waves other than radio waves
- G01S3/782—Systems for determining direction or deviation from predetermined direction
- G01S3/785—Systems for determining direction or deviation from predetermined direction using adjustment of orientation of directivity characteristics of a detector or detector system to give a desired condition of signal derived from that detector or detector system
- G01S3/786—Systems for determining direction or deviation from predetermined direction using adjustment of orientation of directivity characteristics of a detector or detector system to give a desired condition of signal derived from that detector or detector system the desired condition being maintained automatically
- G01S3/7864—T.V. type tracking systems
- G01S3/7865—T.V. type tracking systems using correlation of the live video image with a stored image
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【課題】 本発明は、大きい円形開口に匹敵する高い解
像度のが得られる回転ストリップ開口放射線センサシス
テムおよび測定方法を得ることを目的とする。 【解決手段】 光軸を中心に回動して情景の時間的に順
次の映像を生成する回転ストリップ開口12を備えた望遠
鏡11と、この望遠鏡の焦点平面に配置され2次元の検出
器アレイ17と、検出器アレイ17に結合され、ストリップ
開口12が望遠鏡の光軸を中心に回動するときストリップ
開口映像の複数の映像フレームを記録し、回転するスト
リップ開口の空間周波数カットオフよりも大きい空間周
波数に対する記録された映像から観察された情景の評価
を合成する評価プロセッサ23を備えている信号プロセッ
サ20とを備えていることを特徴とする。
像度のが得られる回転ストリップ開口放射線センサシス
テムおよび測定方法を得ることを目的とする。 【解決手段】 光軸を中心に回動して情景の時間的に順
次の映像を生成する回転ストリップ開口12を備えた望遠
鏡11と、この望遠鏡の焦点平面に配置され2次元の検出
器アレイ17と、検出器アレイ17に結合され、ストリップ
開口12が望遠鏡の光軸を中心に回動するときストリップ
開口映像の複数の映像フレームを記録し、回転するスト
リップ開口の空間周波数カットオフよりも大きい空間周
波数に対する記録された映像から観察された情景の評価
を合成する評価プロセッサ23を備えている信号プロセッ
サ20とを備えていることを特徴とする。
Description
【0001】
【発明の属する技術分野】本発明は、回転開口放射線測
定装置および放射線検出方法に関し、特に、合成された
超解像度情景が複数の回転ストリップ開口映像の測定か
ら評価される回転ストリップ(部分的)開口映像生成用
放射線測定装置および方法に関する。
定装置および放射線検出方法に関し、特に、合成された
超解像度情景が複数の回転ストリップ開口映像の測定か
ら評価される回転ストリップ(部分的)開口映像生成用
放射線測定装置および方法に関する。
【0002】
【従来の技術】宇宙空間で使用するプラットホームから
高解像度の映像を得るために、例えば通常のセンサ技術
では大型で、重く、展開可能な光学システムの能動制御
装置を使用している。要求される使命および1次反射鏡
の大きさに応じて、センサを構成する全ての光学素子の
位置、傾斜、形状、および整列の制御のために1次反射
鏡セグメントの周期的な位置および傾斜の制御が行われ
る。
高解像度の映像を得るために、例えば通常のセンサ技術
では大型で、重く、展開可能な光学システムの能動制御
装置を使用している。要求される使命および1次反射鏡
の大きさに応じて、センサを構成する全ての光学素子の
位置、傾斜、形状、および整列の制御のために1次反射
鏡セグメントの周期的な位置および傾斜の制御が行われ
る。
【0003】本発明と同じ解像度を有する全開口システ
ムは比較的大きい開口面積のために多量の集光能力を有
している。しかしながら、軌道上に数十メートルの直径
の開口を配置するためには、本発明と競合する全開口シ
ステムは次のようなことが必要になる。すなわち光学シ
ステムの直径が宇宙船の発射装置の直径よりも大きくな
るとセグメント化された光学表面および折畳まれた支持
構造が必要になる。また、重量を節減するために薄い変
形可能な反射鏡が使用される場合には複雑で潜在的に高
い帯域幅の適応光学技術が必要である。さらに、位相ア
レイとして構成される場合には複雑な位置および瞳孔整
合制御が要求される。それ故、全開口システムは本発明
と比較して比較的重く、コストが高く、また高い技術的
リスクを有している。
ムは比較的大きい開口面積のために多量の集光能力を有
している。しかしながら、軌道上に数十メートルの直径
の開口を配置するためには、本発明と競合する全開口シ
ステムは次のようなことが必要になる。すなわち光学シ
ステムの直径が宇宙船の発射装置の直径よりも大きくな
るとセグメント化された光学表面および折畳まれた支持
構造が必要になる。また、重量を節減するために薄い変
形可能な反射鏡が使用される場合には複雑で潜在的に高
い帯域幅の適応光学技術が必要である。さらに、位相ア
レイとして構成される場合には複雑な位置および瞳孔整
合制御が要求される。それ故、全開口システムは本発明
と比較して比較的重く、コストが高く、また高い技術的
リスクを有している。
【0004】本発明に関係する従来技術は米国特許第5
243351号明細書に記載されており、スピンApシ
ステムと呼ばれている。スピンApシステムと本発明の
方法(超スピンApシステムと呼ぶ)との共通性は、映
像の合成または情景の評価のために時間的に整合したス
トリップ開口を使用する点にある。
243351号明細書に記載されており、スピンApシ
ステムと呼ばれている。スピンApシステムと本発明の
方法(超スピンApシステムと呼ぶ)との共通性は、映
像の合成または情景の評価のために時間的に整合したス
トリップ開口を使用する点にある。
【0005】
【発明が解決しようとする課題】それ故、本発明の目的
は、合成された超解像度情景が複数の回転ストリップ開
口映像測定から評価される回転ストリップ開口映像放射
線測定装置および測定方法を提供することである。
は、合成された超解像度情景が複数の回転ストリップ開
口映像測定から評価される回転ストリップ開口映像放射
線測定装置および測定方法を提供することである。
【0006】
【課題を解決するための手段】上記の、およびその他の
目的を達成するために、本発明の1形態においては、合
成された超解像度情景が複数の回転ストリップ開口映像
測定から評価される回転ストリップ開口映像放射線測定
装置および測定方法が提供される。ここで使用される超
解像度とは開口を限定している空間周波数カットオフよ
り大きい空間周波数に対して情景情報を正確に評価する
能力を言う。特にここに記載したセンサシステムおよび
データ処理方法に対して、超解像度とはストリップ開口
システムの光学的カットオフを越えて評価された情景の
情報内容を拡張すること、および/または回転する開口
の形状と関係する最大の相関長に等しい半径を有する等
価全円形開口の光学的カットオフを越えて拡張すること
を意味する。表現“ストリップ開口”は米国特許第52
43351号明細書に記載されたような一般的な開口形
状を意味するものである。
目的を達成するために、本発明の1形態においては、合
成された超解像度情景が複数の回転ストリップ開口映像
測定から評価される回転ストリップ開口映像放射線測定
装置および測定方法が提供される。ここで使用される超
解像度とは開口を限定している空間周波数カットオフよ
り大きい空間周波数に対して情景情報を正確に評価する
能力を言う。特にここに記載したセンサシステムおよび
データ処理方法に対して、超解像度とはストリップ開口
システムの光学的カットオフを越えて評価された情景の
情報内容を拡張すること、および/または回転する開口
の形状と関係する最大の相関長に等しい半径を有する等
価全円形開口の光学的カットオフを越えて拡張すること
を意味する。表現“ストリップ開口”は米国特許第52
43351号明細書に記載されたような一般的な開口形
状を意味するものである。
【0007】本発明のシステムは時間的に連続する、或
いは順次の映像を生成する回転ストリップ開口望遠鏡を
備えている。回転ストリップ開口望遠鏡は典型的に回転
ストリップ開口1次反射器および2次反射器を備えてい
る。2次元の検出器アレイが設けられて望遠鏡の焦点平
面に位置する映像を検出する。回転補償装置が設けられ
て、アレイの検出器の積分時間中の回転により軌跡が描
かれるのを阻止するために使用されている。信号プロセ
ッサが設けられて、望遠鏡の光学軸を中心としてストリ
ップ開口が回転するとき望遠鏡により映像の形成された
目標の情景の複数の映像フレームを記録し、回転するス
トリップ開口の空間周波数カットオフより大きい空間周
波数、および/または回転開口システムの全円形開口の
光学的カットオフを越える記録された映像から観察され
た情景の超解像度評価を合成する。
いは順次の映像を生成する回転ストリップ開口望遠鏡を
備えている。回転ストリップ開口望遠鏡は典型的に回転
ストリップ開口1次反射器および2次反射器を備えてい
る。2次元の検出器アレイが設けられて望遠鏡の焦点平
面に位置する映像を検出する。回転補償装置が設けられ
て、アレイの検出器の積分時間中の回転により軌跡が描
かれるのを阻止するために使用されている。信号プロセ
ッサが設けられて、望遠鏡の光学軸を中心としてストリ
ップ開口が回転するとき望遠鏡により映像の形成された
目標の情景の複数の映像フレームを記録し、回転するス
トリップ開口の空間周波数カットオフより大きい空間周
波数、および/または回転開口システムの全円形開口の
光学的カットオフを越える記録された映像から観察され
た情景の超解像度評価を合成する。
【0008】したがって、本発明は、複数の回転ストリ
ップ開口映像測定からの超解像度放射線測定情景評価を
合成する回転ストリップ(部分的)開口映像放射線測定
装置を提供する。本発明はしたがって、同じ重量の従来
のスピンApシステムと比較したとき、或いは同じ重量
の従来の全開口システムと比較したとき、改善された高
い解像度の映像を提供する。
ップ開口映像測定からの超解像度放射線測定情景評価を
合成する回転ストリップ(部分的)開口映像放射線測定
装置を提供する。本発明はしたがって、同じ重量の従来
のスピンApシステムと比較したとき、或いは同じ重量
の従来の全開口システムと比較したとき、改善された高
い解像度の映像を提供する。
【0009】上述の説明を要約すると、本発明のセンサ
およびプロセッサによって行われる合成処理の1実施形
態は以下のとおりである。望遠鏡の光学軸を中心として
ストリップ開口が回転するとき、次のような動作が行わ
れる。回転補償装置は検出器の積分時間中逆回転してそ
れによって時間的に静止している映像が与えられる。映
像フレームは記録され、保存される。回転する(情景に
対して)検出器アレイ構造が選択されるならば、捕捉さ
れたフレームは座標変換され、(超解像度ではない)合
成された映像の基準グリッドに補間される。データはフ
ーリエ変換され、記憶される。新しいフレームは記録お
よび保存される。ランダムな照準線エラーにより記録さ
れたデータのフレーム毎の不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換は照準線エラーに
対して補正され、記憶される。上述のステップは各スト
リップ開口映像フレームに対して順次反復され、或いは
フレームは順次捕捉され、記憶され、その後グローバル
な照準線の評価がフレームを整合するために決定され
る。
およびプロセッサによって行われる合成処理の1実施形
態は以下のとおりである。望遠鏡の光学軸を中心として
ストリップ開口が回転するとき、次のような動作が行わ
れる。回転補償装置は検出器の積分時間中逆回転してそ
れによって時間的に静止している映像が与えられる。映
像フレームは記録され、保存される。回転する(情景に
対して)検出器アレイ構造が選択されるならば、捕捉さ
れたフレームは座標変換され、(超解像度ではない)合
成された映像の基準グリッドに補間される。データはフ
ーリエ変換され、記憶される。新しいフレームは記録お
よび保存される。ランダムな照準線エラーにより記録さ
れたデータのフレーム毎の不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換は照準線エラーに
対して補正され、記憶される。上述のステップは各スト
リップ開口映像フレームに対して順次反復され、或いは
フレームは順次捕捉され、記憶され、その後グローバル
な照準線の評価がフレームを整合するために決定され
る。
【0010】個々のストリップ開口フレームが整合され
たならば、超解像度合成処理が監視される。超解像度合
成処理の1実施形態では以下のようなステップが行われ
る。通過帯域拡張の所望の量は合成された超解像度映像
に対するサンプル間隔を決定するために使用される。座
標系および超解像度映像と関連する空間グリッドは高い
解像度または超解像度グリッドと呼ばれる。超解像度合
成映像の最初の評価および高解像度グリッドにおける超
解像度映像の平均値が決定される。超解像度映像に対す
る最初の評価はデータに接近するスピンAp処理を適用
し、低雑音サンプルを得るためにローパス空間周波数フ
ィルタ処理を行い、全開口伝達関数に対して補償し、画
素複製により超解像度グリッドに評価を拡張することに
よって得られることができる。超解像度映像の平均値に
対する最初の評価は同様の方法で得ることができる。
たならば、超解像度合成処理が監視される。超解像度合
成処理の1実施形態では以下のようなステップが行われ
る。通過帯域拡張の所望の量は合成された超解像度映像
に対するサンプル間隔を決定するために使用される。座
標系および超解像度映像と関連する空間グリッドは高い
解像度または超解像度グリッドと呼ばれる。超解像度合
成映像の最初の評価および高解像度グリッドにおける超
解像度映像の平均値が決定される。超解像度映像に対す
る最初の評価はデータに接近するスピンAp処理を適用
し、低雑音サンプルを得るためにローパス空間周波数フ
ィルタ処理を行い、全開口伝達関数に対して補償し、画
素複製により超解像度グリッドに評価を拡張することに
よって得られることができる。超解像度映像の平均値に
対する最初の評価は同様の方法で得ることができる。
【0011】一度開始点が設定されると非線形超解像度
および相互作用超解像度技術が現在の評価の非線形関数
と超解像度情景の現在の評価との積を形成することによ
り超解像度情景の次の評価が形成される。非線形関数は
測定されたフレームの数ニヨッテスケールされた複数の
回転測定の指数積である。指数の論点は測定、ストリッ
プ開口センサの応答関数、および現在の超解像度情景評
価の関数である。各回転位置に関連する指数の論点は補
助関数により特定の回転位置に対する反射されたストリ
ップ開口の応答関数のコンボリューションを実行するこ
とによって得られる。補助関数発生の1素子は、特定の
回転位置に対するストリップ開口の空間応答関数を有す
る超解像度情景の現在の評価のコンボリューションによ
り画素測定を分割することからなる。
および相互作用超解像度技術が現在の評価の非線形関数
と超解像度情景の現在の評価との積を形成することによ
り超解像度情景の次の評価が形成される。非線形関数は
測定されたフレームの数ニヨッテスケールされた複数の
回転測定の指数積である。指数の論点は測定、ストリッ
プ開口センサの応答関数、および現在の超解像度情景評
価の関数である。各回転位置に関連する指数の論点は補
助関数により特定の回転位置に対する反射されたストリ
ップ開口の応答関数のコンボリューションを実行するこ
とによって得られる。補助関数発生の1素子は、特定の
回転位置に対するストリップ開口の空間応答関数を有す
る超解像度情景の現在の評価のコンボリューションによ
り画素測定を分割することからなる。
【0012】本発明の放射線センサシステムおよび映像
合成処理方法は特定の使命の要求および技術的妥協に依
存するから、放射線センサシステムおよび映像合成処理
方法は、フレーム間の照準線のランダムおよび/または
系統的なドリフトを補償する手段、およびイソプラナテ
ィック光学系の点の広がる作用、映像の運動および検出
器の有限の電子的帯域幅によるフィールド点と関係のな
い映像の軌跡、映像の補償されない回転運動により生じ
たフィールド点と関係する映像の軌跡のような、過去お
よび後の知識によるエラー源を補償する手段を具備する
こともできる。情景構造に関する仮定を実効するための
マルコフ(Markov )ランダムフィールドの使用、全エ
ネルギ測定基準の使用、別の条件確率密度関数の使用の
ような過去および後の知識を伴うことにより、前述の説
明した処理技術の率直な変形が得られる。本発明のシス
テムおよび方法において使用される補償技術の性能は過
去および後の知識の正確度および実効信号対雑音比に依
存する。
合成処理方法は特定の使命の要求および技術的妥協に依
存するから、放射線センサシステムおよび映像合成処理
方法は、フレーム間の照準線のランダムおよび/または
系統的なドリフトを補償する手段、およびイソプラナテ
ィック光学系の点の広がる作用、映像の運動および検出
器の有限の電子的帯域幅によるフィールド点と関係のな
い映像の軌跡、映像の補償されない回転運動により生じ
たフィールド点と関係する映像の軌跡のような、過去お
よび後の知識によるエラー源を補償する手段を具備する
こともできる。情景構造に関する仮定を実効するための
マルコフ(Markov )ランダムフィールドの使用、全エ
ネルギ測定基準の使用、別の条件確率密度関数の使用の
ような過去および後の知識を伴うことにより、前述の説
明した処理技術の率直な変形が得られる。本発明のシス
テムおよび方法において使用される補償技術の性能は過
去および後の知識の正確度および実効信号対雑音比に依
存する。
【0013】本発明のシステムおよび方法は、回転スト
リップ開口を備えた光学系から超解像度映像を提供する
ことのできる回転開口形式の改良された宇宙船搭載光学
映像センサおよび処理構造を提供する。本発明の映像測
定および合成手順は、リスクが低く、コストが安く、軽
量で、製造が簡単であり、また宇宙船搭載映像装置用と
して複雑で大きい開口の適応光学系を展開するのに比較
して簡単な展開構造を与える利点がある。
リップ開口を備えた光学系から超解像度映像を提供する
ことのできる回転開口形式の改良された宇宙船搭載光学
映像センサおよび処理構造を提供する。本発明の映像測
定および合成手順は、リスクが低く、コストが安く、軽
量で、製造が簡単であり、また宇宙船搭載映像装置用と
して複雑で大きい開口の適応光学系を展開するのに比較
して簡単な展開構造を与える利点がある。
【0014】本発明のシステムおよび方法は、特に望遠
鏡の重量がペイロードの重量の大きい割合を占めている
場合に衛星軌道から高い解像度の映像を提供する。衛星
の発射のコストは主として軌道中の重量の関数である。
動作システムの状態に応じて、所定の許容可能な望遠鏡
の重量に対して、本発明は従来のスピンApシステムよ
り50%まで解像度を増加させることができる。本発明
は高解像度(大きい開口と等価)の宇宙船搭載観察装置
としてそのようなシステムで使用されることができる。
鏡の重量がペイロードの重量の大きい割合を占めている
場合に衛星軌道から高い解像度の映像を提供する。衛星
の発射のコストは主として軌道中の重量の関数である。
動作システムの状態に応じて、所定の許容可能な望遠鏡
の重量に対して、本発明は従来のスピンApシステムよ
り50%まで解像度を増加させることができる。本発明
は高解像度(大きい開口と等価)の宇宙船搭載観察装置
としてそのようなシステムで使用されることができる。
【0015】本発明の超スピンApシステムは利用でき
る測定の使用により全開口映像を生成し、ストリップ開
口と関連した最大の相関長と関連した空間周波数カット
オフを越えた、すなわち、等価なスピンApシステム
(帯域制限されている)の合成全開口空間周波数カット
オフを越えた空間周波数カットオフを有する。本発明の
超スピンApシステムは最大の相関長を有するスピンA
pシステムから発生された映像を表示する。それ故所定
の実効合成開口寸法に対する低いペイロード重量が得ら
れる変形されたスピンAp映像処理方法が得られる利点
がある。さらに、リスクが低く、コストが安く、軽量
で、製造が簡単であり、また宇宙船搭載映像装置用とし
て複雑で大きい全円形開口の適応光学系を展開するもの
(または位相アレイ)のに比較して簡単な展開構造を与
える利点がある。さらに、超解像度評価技術はストリッ
プ開口の通過帯域の外側の情景情報の評価を可能にする
から、超スピンAp処理手順はより少ない回転サンプル
によるスピンApと同じ合成映像品質を提供することが
できる。
る測定の使用により全開口映像を生成し、ストリップ開
口と関連した最大の相関長と関連した空間周波数カット
オフを越えた、すなわち、等価なスピンApシステム
(帯域制限されている)の合成全開口空間周波数カット
オフを越えた空間周波数カットオフを有する。本発明の
超スピンApシステムは最大の相関長を有するスピンA
pシステムから発生された映像を表示する。それ故所定
の実効合成開口寸法に対する低いペイロード重量が得ら
れる変形されたスピンAp映像処理方法が得られる利点
がある。さらに、リスクが低く、コストが安く、軽量
で、製造が簡単であり、また宇宙船搭載映像装置用とし
て複雑で大きい全円形開口の適応光学系を展開するもの
(または位相アレイ)のに比較して簡単な展開構造を与
える利点がある。さらに、超解像度評価技術はストリッ
プ開口の通過帯域の外側の情景情報の評価を可能にする
から、超スピンAp処理手順はより少ない回転サンプル
によるスピンApと同じ合成映像品質を提供することが
できる。
【0016】本発明の種々の特徴および利点は、添付図
面に関連した以下の詳細な説明に関して一層容易に理解
され、添付図面において、同一の参照番号は同一の構造
素子を示している。
面に関連した以下の詳細な説明に関して一層容易に理解
され、添付図面において、同一の参照番号は同一の構造
素子を示している。
【0017】
【発明の実施の形態】米国特許第5,243,351 号明細書の
図1に対応する図1を参照すると、本発明の原理によっ
て変更された回転開口映像生成放射線センサシステム10
の一例が示されている。米国特許第5,243,351 号の内容
が本明細書において参照文献とされている。回転開口映
像生成放射線センサシステム10は、超解像度の情景の評
価を合成し、他方で、照準線ジッタを取除き、同じ重量
の通常の光学システムおよび同じ重量の回転開口システ
ムとの比較において改良された高解像度の映像を供給す
るように構成されている。回転開口映像生成放射線セン
サシステム10は、1次反射器12a および2次反射器12b
を含む回転ストリップ開口望遠鏡11を具備している。所
定の環境の下では3次反射器(図示されていない)が望
遠鏡11において使用されることがある。
図1に対応する図1を参照すると、本発明の原理によっ
て変更された回転開口映像生成放射線センサシステム10
の一例が示されている。米国特許第5,243,351 号の内容
が本明細書において参照文献とされている。回転開口映
像生成放射線センサシステム10は、超解像度の情景の評
価を合成し、他方で、照準線ジッタを取除き、同じ重量
の通常の光学システムおよび同じ重量の回転開口システ
ムとの比較において改良された高解像度の映像を供給す
るように構成されている。回転開口映像生成放射線セン
サシステム10は、1次反射器12a および2次反射器12b
を含む回転ストリップ開口望遠鏡11を具備している。所
定の環境の下では3次反射器(図示されていない)が望
遠鏡11において使用されることがある。
【0018】本発明を開示する目的のために、システム
10は、地球指向アンテナ13a を有する静止セクション13
を具備している衛星の形態で示されている。望遠鏡11
は、プラットフォーム14上に配置され、そこには静止セ
クション13も結合されている。回転開口映像生成放射線
センサシステム10は、(ストリップ開口12を有してい
る)1次反射器12a が望遠鏡11の光学軸の周囲を回転す
る際に望遠鏡11によって映された複数のターゲット情景
の映像フレームを記録するように構成されている。照準
線安定反射器15および映像回転除去装置16は、検出器ア
レイ17によって感知される前に映像の安定化および回転
除去を行うように構成された望遠鏡11の光路に沿って配
置されている。回転除去装置16は、回転補償制御装置19
の制御の下で、検出器アレイ17を構成している検出器の
積分時間中に映像を逆に回転させ、それによって静止映
像を供給する。照準線安定反射器15は、システム10の軌
道力学のために、照準線エラーの大きい帯域幅ならびに
照準線エラーを取除くために(信号プロセッサあるいは
別の専用制御システム等によって供給された)照準線制
御システムによって使用される。ターゲット情景は、映
像フレームを処理するように調整された信号プロセッサ
20に結合された望遠鏡11の焦点平面に位置された検出器
アレイ17上に映像が形成される。その代りに、信号プロ
セッサ20が専用の基礎プロセッサを具備していることも
ある。望遠鏡11、検出器アレイ17、および関連したハー
ドウェアは、本明細書においてセンサシステム21と呼ば
れる。センサシステム21によって生成された個々の映像
フレームは、信号プロセッサ20において処理および結合
され、それによって、本発明のシステム10および方法30
に従って超解像度情景評価を合成し、それは特に図2乃
至11を参照にして説明されている。
10は、地球指向アンテナ13a を有する静止セクション13
を具備している衛星の形態で示されている。望遠鏡11
は、プラットフォーム14上に配置され、そこには静止セ
クション13も結合されている。回転開口映像生成放射線
センサシステム10は、(ストリップ開口12を有してい
る)1次反射器12a が望遠鏡11の光学軸の周囲を回転す
る際に望遠鏡11によって映された複数のターゲット情景
の映像フレームを記録するように構成されている。照準
線安定反射器15および映像回転除去装置16は、検出器ア
レイ17によって感知される前に映像の安定化および回転
除去を行うように構成された望遠鏡11の光路に沿って配
置されている。回転除去装置16は、回転補償制御装置19
の制御の下で、検出器アレイ17を構成している検出器の
積分時間中に映像を逆に回転させ、それによって静止映
像を供給する。照準線安定反射器15は、システム10の軌
道力学のために、照準線エラーの大きい帯域幅ならびに
照準線エラーを取除くために(信号プロセッサあるいは
別の専用制御システム等によって供給された)照準線制
御システムによって使用される。ターゲット情景は、映
像フレームを処理するように調整された信号プロセッサ
20に結合された望遠鏡11の焦点平面に位置された検出器
アレイ17上に映像が形成される。その代りに、信号プロ
セッサ20が専用の基礎プロセッサを具備していることも
ある。望遠鏡11、検出器アレイ17、および関連したハー
ドウェアは、本明細書においてセンサシステム21と呼ば
れる。センサシステム21によって生成された個々の映像
フレームは、信号プロセッサ20において処理および結合
され、それによって、本発明のシステム10および方法30
に従って超解像度情景評価を合成し、それは特に図2乃
至11を参照にして説明されている。
【0019】本発明による超スピンApシステム10およ
び方法30の特定の実施形態が以下に説明される。特に、
図2乃至11において、本発明の原理にしたがった超ス
ピンAp処理方法30の種々の実施形態が示されており、
それらは、本発明の原理を実行するために適切に変更さ
れた図1に示された回転開口放射線センサシステムおよ
び方法によって使用される。観測シナリオおよび時間
線、使用できる従来技術の知識の量、および処理の流れ
における使用可能な計算スループットに基づいて特定の
実施形態が選択される。
び方法30の特定の実施形態が以下に説明される。特に、
図2乃至11において、本発明の原理にしたがった超ス
ピンAp処理方法30の種々の実施形態が示されており、
それらは、本発明の原理を実行するために適切に変更さ
れた図1に示された回転開口放射線センサシステムおよ
び方法によって使用される。観測シナリオおよび時間
線、使用できる従来技術の知識の量、および処理の流れ
における使用可能な計算スループットに基づいて特定の
実施形態が選択される。
【0020】図2において、本発明によるシステム10お
よび方法30の第1の実施形態が示されており、そこにお
いて、最高レベルの超スピンAp映像合成方法30が信号
プロセッサ20において実行される。回転開口センサ21
は、一組のストリップ開口映像22を捕捉および生成し、
それは超スピンApプロセッサ23に移送される。超スピ
ンApプロセッサ23は、超解像度情景評価24を合成す
る。
よび方法30の第1の実施形態が示されており、そこにお
いて、最高レベルの超スピンAp映像合成方法30が信号
プロセッサ20において実行される。回転開口センサ21
は、一組のストリップ開口映像22を捕捉および生成し、
それは超スピンApプロセッサ23に移送される。超スピ
ンApプロセッサ23は、超解像度情景評価24を合成す
る。
【0021】図3を参照すると、第2のシステム10およ
び方法30が示されている。回転開口センサ21は一組のス
トリップ開口映像22を捕捉および生成し、それらはフレ
ーム整合エラープロセッサ25に移される。フレーム整合
エラープロセッサ25は、その組の個々のフレーム間の残
りの整合エラーを決定する。例えば、フレームの整合は
以下のように達成される。回転ストリップ開口が望遠鏡
11の光学軸の周囲を回転しているとき、以下のように動
作する。回転除去装置16は、検出器アレイ17の積分時間
中、逆に回転し、それによって時間的に静止映像を供給
する。映像フレームは記録され、蓄えられる。(情景に
関して)回転する検出器アレイ構造が選択された場合、
捕捉されたフレームは、調和して変換され、(超解像度
でない)合成映像の基準グリッドに補間される。データ
はフーリエ変換され、記憶される。新しいフレームが記
録され、記憶される。ランダムな照準線エラーのために
生じた記録されたデータのフレーム毎の不整合の評価が
得られる。ストリップ開口映像、すなわちフーリエ変換
は、それらの照準線エラーを補正され、記憶される。前
述のステップは、各ストリップ開口映像フレーム毎に連
続的に反復され、そうでなければ、フレームは順次連続
的に捕捉および記憶され、その後、フレームを整合させ
るために照準線の全体的な評価が決定される。整合プロ
セッサ25は、一組の補正されたフレーム、あるいは一組
の補正されていないフレームおよびフレーム間の整合エ
ラー評価、もしくはその両方を超解像度プロセッサ23に
転送する。補正されたフレーム、あるいは補正されてい
ないフレームおよびそれに対応する一組の整合エラーを
使用して、超解像度評価プロセッサ23は、超解像度情景
評価24を生成する。
び方法30が示されている。回転開口センサ21は一組のス
トリップ開口映像22を捕捉および生成し、それらはフレ
ーム整合エラープロセッサ25に移される。フレーム整合
エラープロセッサ25は、その組の個々のフレーム間の残
りの整合エラーを決定する。例えば、フレームの整合は
以下のように達成される。回転ストリップ開口が望遠鏡
11の光学軸の周囲を回転しているとき、以下のように動
作する。回転除去装置16は、検出器アレイ17の積分時間
中、逆に回転し、それによって時間的に静止映像を供給
する。映像フレームは記録され、蓄えられる。(情景に
関して)回転する検出器アレイ構造が選択された場合、
捕捉されたフレームは、調和して変換され、(超解像度
でない)合成映像の基準グリッドに補間される。データ
はフーリエ変換され、記憶される。新しいフレームが記
録され、記憶される。ランダムな照準線エラーのために
生じた記録されたデータのフレーム毎の不整合の評価が
得られる。ストリップ開口映像、すなわちフーリエ変換
は、それらの照準線エラーを補正され、記憶される。前
述のステップは、各ストリップ開口映像フレーム毎に連
続的に反復され、そうでなければ、フレームは順次連続
的に捕捉および記憶され、その後、フレームを整合させ
るために照準線の全体的な評価が決定される。整合プロ
セッサ25は、一組の補正されたフレーム、あるいは一組
の補正されていないフレームおよびフレーム間の整合エ
ラー評価、もしくはその両方を超解像度プロセッサ23に
転送する。補正されたフレーム、あるいは補正されてい
ないフレームおよびそれに対応する一組の整合エラーを
使用して、超解像度評価プロセッサ23は、超解像度情景
評価24を生成する。
【0022】図4を参照すると、第3のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に送られる。フレーム整
合エラープロセッサ25は、その組の個々のフレームの間
の残りの整合エラーを決定する。整合プロセッサ25は、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を回転開口帯域限定映像合成プロセッサ
26および超解像度評価プロセッサ23に送る。回転開口帯
域限定映像合成プロセッサ26は、米国特許第5,243,351
号明細書に記載された回転開口処理によって全円形開口
映像を表す映像を生成する。一組の補正されたフレー
ム、あるいは一組の補正されていないフレームおよびフ
レーム間の整合エラー評価、もしくはその両方、ならび
に回転開口帯域限定映像を使用して、超解像度評価プロ
セッサは超解像度情景24の評価を生成する。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に送られる。フレーム整
合エラープロセッサ25は、その組の個々のフレームの間
の残りの整合エラーを決定する。整合プロセッサ25は、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を回転開口帯域限定映像合成プロセッサ
26および超解像度評価プロセッサ23に送る。回転開口帯
域限定映像合成プロセッサ26は、米国特許第5,243,351
号明細書に記載された回転開口処理によって全円形開口
映像を表す映像を生成する。一組の補正されたフレー
ム、あるいは一組の補正されていないフレームおよびフ
レーム間の整合エラー評価、もしくはその両方、ならび
に回転開口帯域限定映像を使用して、超解像度評価プロ
セッサは超解像度情景24の評価を生成する。
【0023】図5を参照すると、第4のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および発生し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に転送する。さ
らに、正確に整合されたフレーム、あるいは一組の補正
されていないフレームおよびフレーム間の整合エラー評
価が回転開口帯域限定映像合成プロセッサ26に送られ
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、映像合成デー
タベース27に記憶された情景あるいは物体の構造の過去
の知識を組入れる。補正されたフレーム、あるいは一組
の補正されていないフレームおよび対応する一組の整合
エラー評価、ならびに回転開口帯域限定映像を使用し
て、超解像度評価プロセッサは超解像度情景評価24を生
成し、一方で、データベース27に記憶された情景あるい
は物体の構造の過去の知識を組入れる。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および発生し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に転送する。さ
らに、正確に整合されたフレーム、あるいは一組の補正
されていないフレームおよびフレーム間の整合エラー評
価が回転開口帯域限定映像合成プロセッサ26に送られ
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、映像合成デー
タベース27に記憶された情景あるいは物体の構造の過去
の知識を組入れる。補正されたフレーム、あるいは一組
の補正されていないフレームおよび対応する一組の整合
エラー評価、ならびに回転開口帯域限定映像を使用し
て、超解像度評価プロセッサは超解像度情景評価24を生
成し、一方で、データベース27に記憶された情景あるい
は物体の構造の過去の知識を組入れる。
【0024】図6を参照すると、第5のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に送る。さら
に、正確に整合されたフレーム、あるいは一組の補正さ
れていないフレームおよびフレーム間の整合エラー評価
が回転開口帯域限定映像合成プロセッサ26に転送され
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、データベース
27に記憶された情景あるいは物体の構造の過去の知識を
組入れる。帯域限定された回転開口映像は、訓練アルゴ
リズムプロセッサ28に移送される。帯域限定された回転
開口映像から得られた情報に基づいて、訓練プロセッサ
は、情景の内容および構造に関する情報を抽出し、それ
はデータベース27に含まれた過去の情報を変更するため
に使用される。補正されたフレーム、あるいは補正され
ていないフレームおよび対応する一組の整合エラー評
価、ならびに回転開口帯域限定映像を使用して、超解像
度評価プロセッサ23は超解像度情景評価24を生成し、一
方で、データベース27に記憶され、訓練アルゴリズムプ
ロセッサ28によって変更された情景あるいは物体の構造
の過去の知識を組入れる。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に送る。さら
に、正確に整合されたフレーム、あるいは一組の補正さ
れていないフレームおよびフレーム間の整合エラー評価
が回転開口帯域限定映像合成プロセッサ26に転送され
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、データベース
27に記憶された情景あるいは物体の構造の過去の知識を
組入れる。帯域限定された回転開口映像は、訓練アルゴ
リズムプロセッサ28に移送される。帯域限定された回転
開口映像から得られた情報に基づいて、訓練プロセッサ
は、情景の内容および構造に関する情報を抽出し、それ
はデータベース27に含まれた過去の情報を変更するため
に使用される。補正されたフレーム、あるいは補正され
ていないフレームおよび対応する一組の整合エラー評
価、ならびに回転開口帯域限定映像を使用して、超解像
度評価プロセッサ23は超解像度情景評価24を生成し、一
方で、データベース27に記憶され、訓練アルゴリズムプ
ロセッサ28によって変更された情景あるいは物体の構造
の過去の知識を組入れる。
【0025】超スピンAp処理の実現方法が以下に説明
される。図7および8において、図2乃至6のシステム
10および方法30において使用される回転開口の超解像度
処理方法23の実施形態が示されている。図7および8に
示された方法23において、回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それはフレ
ーム整合エラープロセッサ25に送られる。フレーム整合
エラープロッサ25は、その一組の個々のフレーム間の残
りの整合エラーを決定する。整合プロセッサは、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に転送する。
される。図7および8において、図2乃至6のシステム
10および方法30において使用される回転開口の超解像度
処理方法23の実施形態が示されている。図7および8に
示された方法23において、回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それはフレ
ーム整合エラープロセッサ25に送られる。フレーム整合
エラープロッサ25は、その一組の個々のフレーム間の残
りの整合エラーを決定する。整合プロセッサは、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に転送する。
【0026】さらに、正確に整合されたフレームおよび
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に移される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた全円形開口映像を表す映像を生成する。帯域が限定
された回転開口映像は、超解像度情景の第1の評価33を
供給するために変更される。
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に移される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた全円形開口映像を表す映像を生成する。帯域が限定
された回転開口映像は、超解像度情景の第1の評価33を
供給するために変更される。
【0027】現在の情景評価34の観察から得られる各フ
レームの測定結果を決定するために現在の超解像度情景
評価33およびシステム応答関数37の過去の情報が使用さ
れる。予測された測定結果34は、現在の超解像度情景評
価33を個々のフレーム37のシステム応答関数でコンボリ
ューション処理することによって決定される。処理の流
れにおける次のステップ39および40では、測定結果22を
予測された測定結果34で除算し、そこから1を減算し、
その後、各フレームの反射された応答関数38でコンボリ
ューション処理する。結果として生じたコンボリューシ
ョン40は、各フレーム41毎に合計され、その後、ステッ
プ42においてスケーリング係数43で乗算される。第1の
超解像度評価における終わりの方のステップにおいて、
ステップ39乃至43の出力が指数関数化され(43)、指数
関数化された結果が現在の超解像度情景評価34で乗算さ
れる(ステップ45)。ステップ21乃至45の結果、回転開
口超解像度情景の第1の真の評価46が得られる。その
後、回転開口超解像度情景の第1の真の評価46が現在の
情景評価33として使用され、このプロセスは所望された
回数だけ反復される。
レームの測定結果を決定するために現在の超解像度情景
評価33およびシステム応答関数37の過去の情報が使用さ
れる。予測された測定結果34は、現在の超解像度情景評
価33を個々のフレーム37のシステム応答関数でコンボリ
ューション処理することによって決定される。処理の流
れにおける次のステップ39および40では、測定結果22を
予測された測定結果34で除算し、そこから1を減算し、
その後、各フレームの反射された応答関数38でコンボリ
ューション処理する。結果として生じたコンボリューシ
ョン40は、各フレーム41毎に合計され、その後、ステッ
プ42においてスケーリング係数43で乗算される。第1の
超解像度評価における終わりの方のステップにおいて、
ステップ39乃至43の出力が指数関数化され(43)、指数
関数化された結果が現在の超解像度情景評価34で乗算さ
れる(ステップ45)。ステップ21乃至45の結果、回転開
口超解像度情景の第1の真の評価46が得られる。その
後、回転開口超解像度情景の第1の真の評価46が現在の
情景評価33として使用され、このプロセスは所望された
回数だけ反復される。
【0028】図9および10において、図2乃至6のシ
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。この方
法23において、図9および10において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に移送される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に転送する。
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。この方
法23において、図9および10において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に移送される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に転送する。
【0029】また、正確に整合されたフレームおよびフ
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に送られる。回転開口帯域限定映像合成プロ
セッサ26は、米国特許第5,243,351 号明細書に記載され
た手段によって全円形開口映像を表す映像を生成する。
帯域限定された回転開口映像は、超解像度情景サンプル
手段の第1の評価33を供給するように変更される。
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に送られる。回転開口帯域限定映像合成プロ
セッサ26は、米国特許第5,243,351 号明細書に記載され
た手段によって全円形開口映像を表す映像を生成する。
帯域限定された回転開口映像は、超解像度情景サンプル
手段の第1の評価33を供給するように変更される。
【0030】現在の情景サンプル平均評価34の観察によ
り得られる各フレームの測定を決定するために現在の超
解像度情景サンプル平均評価33およびシステム応答関数
37の過去の情報が使用される。予測された測定結果34
は、現在の超解像度情景サンプル平均評価33を個々のフ
レーム37のシステム応答関数でコンボリューション処理
することによって決定される。処理の流れにおける次の
ステップ39および40では、測定結果22を予測された測定
結果34で除算し、そこから1を減算し、その後、各フレ
ームの反射された応答関数38でコンボリューションする
(ステップ40)。結果として生じたコンボリューション
40は、ステップ41において各フレーム毎に合計され、そ
の後、ステップ42においてスケーリング係数43で乗算さ
れる。第1の超解像度評価における最終のステップにお
いて、ステップ39乃至43の出力が指数関数化(ステップ
44)され、ステップ45において、指数関数化された結果
が現在の超解像度情景サンプル平均評価33と乗算され
る。ステップ21乃至46の結果、回転開口超解像度情景の
第1の真の評価46が得られる。その後、回転開口超解像
度情景の第1の真の評価46が現在の情景サンプル平均評
価33を更新するために使用され、このプロセスは所望さ
れた回数だけ反復される。
り得られる各フレームの測定を決定するために現在の超
解像度情景サンプル平均評価33およびシステム応答関数
37の過去の情報が使用される。予測された測定結果34
は、現在の超解像度情景サンプル平均評価33を個々のフ
レーム37のシステム応答関数でコンボリューション処理
することによって決定される。処理の流れにおける次の
ステップ39および40では、測定結果22を予測された測定
結果34で除算し、そこから1を減算し、その後、各フレ
ームの反射された応答関数38でコンボリューションする
(ステップ40)。結果として生じたコンボリューション
40は、ステップ41において各フレーム毎に合計され、そ
の後、ステップ42においてスケーリング係数43で乗算さ
れる。第1の超解像度評価における最終のステップにお
いて、ステップ39乃至43の出力が指数関数化(ステップ
44)され、ステップ45において、指数関数化された結果
が現在の超解像度情景サンプル平均評価33と乗算され
る。ステップ21乃至46の結果、回転開口超解像度情景の
第1の真の評価46が得られる。その後、回転開口超解像
度情景の第1の真の評価46が現在の情景サンプル平均評
価33を更新するために使用され、このプロセスは所望さ
れた回数だけ反復される。
【0031】図11および12において、図2乃至6の
システム10および方法30において使用される回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図11および12において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に送られる。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。整合されたフレーム29から、超ス
ピンApプロセッサはステップ29a において全ての検出
器およびフレームの測定結果の合計を決定し、全ての検
出器の測定結果と全てのフレームの測定結果の合計は1
より小さい。その合計は、正規化された補正項39a の成
分となる。
システム10および方法30において使用される回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図11および12において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に送られる。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。整合されたフレーム29から、超ス
ピンApプロセッサはステップ29a において全ての検出
器およびフレームの測定結果の合計を決定し、全ての検
出器の測定結果と全てのフレームの測定結果の合計は1
より小さい。その合計は、正規化された補正項39a の成
分となる。
【0032】また、正確に整合されたフレームおよびフ
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に供給される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた手段によって全円形開口映像を表す映像を生成す
る。帯域限定された回転開口映像は、超解像度情景の第
1の評価33を供給するように変更される。
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に供給される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた手段によって全円形開口映像を表す映像を生成す
る。帯域限定された回転開口映像は、超解像度情景の第
1の評価33を供給するように変更される。
【0033】現在の情景評価34の観察から得られる各フ
レームの測定を決定するために現在の超解像度情景評価
33およびシステム応答関数37の過去の情報が使用され
る。予測された測定結果34は、現在の超解像度情景評価
33を個々のフレーム37のシステム応答関数37でコンボリ
ューション処理することによって決定される。予測され
た測定結果34はまた、現在の情景評価33の観察から得ら
れる全てのフレームおよび画素の合計35を決定するため
に使用される。予測された測定結果34は、ステップ36に
おいて1より小さい全てのフレームにわたる全ての画素
および全ての画素にわたって合計され、それによって、
正規化された補正項39a の成分が生成される。
レームの測定を決定するために現在の超解像度情景評価
33およびシステム応答関数37の過去の情報が使用され
る。予測された測定結果34は、現在の超解像度情景評価
33を個々のフレーム37のシステム応答関数37でコンボリ
ューション処理することによって決定される。予測され
た測定結果34はまた、現在の情景評価33の観察から得ら
れる全てのフレームおよび画素の合計35を決定するため
に使用される。予測された測定結果34は、ステップ36に
おいて1より小さい全てのフレームにわたる全ての画素
および全ての画素にわたって合計され、それによって、
正規化された補正項39a の成分が生成される。
【0034】処理の流れにおける次のステップ40では、
測定結果22を予測された測定結果34で除算し、1を減算
した後、それは各フレームの反射された応答関数38でコ
ンボリューション処理される。結果として生じたコンボ
リューション40は、各フレーム41に対して合計され(ス
テップ41)、その後、ステップ29a,35,36 の結果を有す
る正規化された補正項がこの合計から減算される。この
結果は、ステップ42においてスケーリング係数43で乗算
される。第1の超解像度評価における最後の方のステッ
プ44において、ステップ39乃至43の出力が指数関数化さ
れ、ステップ45において、指数関数化された結果が現在
の超解像度情景評価33と乗算される。ステップ21乃至45
の結果が、回転開口超解像度情景の第1の真の評価46と
なる。その後、回転開口超解像度情景の第1の真の評価
46が現在の情景評価33として使用され、このプロセスは
所望された回数だけ反復される。
測定結果22を予測された測定結果34で除算し、1を減算
した後、それは各フレームの反射された応答関数38でコ
ンボリューション処理される。結果として生じたコンボ
リューション40は、各フレーム41に対して合計され(ス
テップ41)、その後、ステップ29a,35,36 の結果を有す
る正規化された補正項がこの合計から減算される。この
結果は、ステップ42においてスケーリング係数43で乗算
される。第1の超解像度評価における最後の方のステッ
プ44において、ステップ39乃至43の出力が指数関数化さ
れ、ステップ45において、指数関数化された結果が現在
の超解像度情景評価33と乗算される。ステップ21乃至45
の結果が、回転開口超解像度情景の第1の真の評価46と
なる。その後、回転開口超解像度情景の第1の真の評価
46が現在の情景評価33として使用され、このプロセスは
所望された回数だけ反復される。
【0035】図13および14において、図2乃至6の
システム10および方法30において使用された回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図13および14において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。フレーム整合エラープロセッサ25は、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に供給する。
システム10および方法30において使用された回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図13および14において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。フレーム整合エラープロセッサ25は、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に供給する。
【0036】さらに、正確に整合されたフレームおよび
フレーム間の整合エラー評価は、回転開口帯域限定合成
映像プロセッサ26に供給される。回転開口帯域限定映像
合成プロセッサ26は、米国特許第5,243,351 号明細書に
記載された全円形開口映像を表す映像を生成する。帯域
が限定された回転開口映像は、超解像度情景アンサンブ
ル平均の第1の評価33を供給するために変更される。
フレーム間の整合エラー評価は、回転開口帯域限定合成
映像プロセッサ26に供給される。回転開口帯域限定映像
合成プロセッサ26は、米国特許第5,243,351 号明細書に
記載された全円形開口映像を表す映像を生成する。帯域
が限定された回転開口映像は、超解像度情景アンサンブ
ル平均の第1の評価33を供給するために変更される。
【0037】現在の情景アンサンブル平均評価34の観察
から得られる各フレームの測定を決定するように現在の
超解像度情景アンサンブル平均評価33およびシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
個々のフレームのシステム応答関数37でコンボリューシ
ョンすることによって決定される。処理の流れにおける
次のステップ39では、測定結果22を予測された測定結果
34で除算し、その後、ステップ40において各フレームの
反射された応答関数38でコンボリューション処理する。
結果として生じたコンボリューション40は、ステップ41
において各フレームに対して合計され、その後、ステッ
プ42においてスケーリング係数43と乗算される。第1の
超解像度評価における終段のステップ45において、ステ
ップ39乃至43の出力が現在の超解像度情景アンサンブル
平均評価34と乗算される。ステップ21乃至45の結果が、
回転開口超解像度情景アンサンブル平均の第1の真の評
価46となる。その後、回転開口超解像度情景の第1の真
の評価46が現在の情景アンサンブル平均評価33として使
用され、このプロセスは所望された回数だけ反復され
る。
から得られる各フレームの測定を決定するように現在の
超解像度情景アンサンブル平均評価33およびシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
個々のフレームのシステム応答関数37でコンボリューシ
ョンすることによって決定される。処理の流れにおける
次のステップ39では、測定結果22を予測された測定結果
34で除算し、その後、ステップ40において各フレームの
反射された応答関数38でコンボリューション処理する。
結果として生じたコンボリューション40は、ステップ41
において各フレームに対して合計され、その後、ステッ
プ42においてスケーリング係数43と乗算される。第1の
超解像度評価における終段のステップ45において、ステ
ップ39乃至43の出力が現在の超解像度情景アンサンブル
平均評価34と乗算される。ステップ21乃至45の結果が、
回転開口超解像度情景アンサンブル平均の第1の真の評
価46となる。その後、回転開口超解像度情景の第1の真
の評価46が現在の情景アンサンブル平均評価33として使
用され、このプロセスは所望された回数だけ反復され
る。
【0038】図15および16おいて、図2乃至6のシ
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。方法23
において、図15および16において示されているよう
に、回転開口センサ21は一組のストリップ開口映像22を
捕捉および生成し、それはフレーム整合エラープロセッ
サ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。方法23
において、図15および16において示されているよう
に、回転開口センサ21は一組のストリップ開口映像22を
捕捉および生成し、それはフレーム整合エラープロセッ
サ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。
【0039】さらに、正確に整合されたフレームおよび
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に転送される。回転開口帯域限定映像合成
プロセッサ26は、米国特許第5,243,351 号明細書に記載
されたように全円形開口映像を表す映像を生成する。帯
域が限定された回転開口映像は、超解像度情景アンサン
ブル平均の第1の評価33を供給するために変更される。
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に転送される。回転開口帯域限定映像合成
プロセッサ26は、米国特許第5,243,351 号明細書に記載
されたように全円形開口映像を表す映像を生成する。帯
域が限定された回転開口映像は、超解像度情景アンサン
ブル平均の第1の評価33を供給するために変更される。
【0040】現在の情景アンサンブル平均評価34の観察
から得られる第1のフレームに対する測定結果を決定す
るように現在の超解像度情景アンサンブル平均評価33お
よび第1の回転開口フレーム測定に関連したシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
第1のフレーム37のシステム応答関数でコンボリューシ
ョン処理することによって決定される。処理の流れにお
ける次のステップ39では、第1の回転開口の測定結果22
を予測された第1の測定結果34で除算し、その後、それ
をステップ40において第1のフレームの反射された応答
関数38でコンボリューション処理する。ステップ39b に
おいて、結果として生じたコンボリューションから1を
減算し、その結果は、ステップ41a において情景アンサ
ンブル平均の現在の評価34で乗算される。この結果は、
ステップ42においてスケーリング係数43と乗算される。
第1の超解像度評価における終段のステップ44a におい
て、ステップ39乃至43の出力が現在の超解像度情景アン
サンブル平均評価34に加算される。ステップ21乃至45の
結果が、回転開口超解像度情景アンサンブル平均の第1
の真の評価46となる。その後、回転開口超解像度情景の
第1の評価46が現在の情景アンサンブル平均評価33とし
て使用され、このプロセスは所望された回数だけ反復さ
れる。
から得られる第1のフレームに対する測定結果を決定す
るように現在の超解像度情景アンサンブル平均評価33お
よび第1の回転開口フレーム測定に関連したシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
第1のフレーム37のシステム応答関数でコンボリューシ
ョン処理することによって決定される。処理の流れにお
ける次のステップ39では、第1の回転開口の測定結果22
を予測された第1の測定結果34で除算し、その後、それ
をステップ40において第1のフレームの反射された応答
関数38でコンボリューション処理する。ステップ39b に
おいて、結果として生じたコンボリューションから1を
減算し、その結果は、ステップ41a において情景アンサ
ンブル平均の現在の評価34で乗算される。この結果は、
ステップ42においてスケーリング係数43と乗算される。
第1の超解像度評価における終段のステップ44a におい
て、ステップ39乃至43の出力が現在の超解像度情景アン
サンブル平均評価34に加算される。ステップ21乃至45の
結果が、回転開口超解像度情景アンサンブル平均の第1
の真の評価46となる。その後、回転開口超解像度情景の
第1の評価46が現在の情景アンサンブル平均評価33とし
て使用され、このプロセスは所望された回数だけ反復さ
れる。
【0041】超解像度の回転開口の式を以下に示す。S
m により示される時間的に静止した情景を実現するた
め、回転開口センサの焦点平面アレイからの雑音のない
検出器出力Gm は次式のように表されることができる。
m により示される時間的に静止した情景を実現するた
め、回転開口センサの焦点平面アレイからの雑音のない
検出器出力Gm は次式のように表されることができる。
【0042】
【数1】 ★下付け文字mはストリップ開口のm番目の回転方向を
示しており、その期間にm番目の画像フレームが捕捉さ
れ、Hm は全体のシステム応答関数である。Smはm番
目の測定に関する情景の実体を表している。これは雑音
のない回転開口画像捕捉の最も一般的なモデルである。
示しており、その期間にm番目の画像フレームが捕捉さ
れ、Hm は全体のシステム応答関数である。Smはm番
目の測定に関する情景の実体を表している。これは雑音
のない回転開口画像捕捉の最も一般的なモデルである。
【0043】システム応答関数による情景のコンボリュ
ーションのディスクリートな形態は次式により表される
ことができる。
ーションのディスクリートな形態は次式により表される
ことができる。
【0044】
【数2】 ★ここでjは位置/rjからの検出器出力を示してお
り、kは焦点平面の空間的位置/rkへ幾何学的に投影
されたディスクリートな情景素子を示している。
り、kは焦点平面の空間的位置/rkへ幾何学的に投影
されたディスクリートな情景素子を示している。
【0045】システム応答関数Hm はm番目の画像捕捉
用の光学システム点分布関数PSFm のコンボリューシ
ョンであり、個々の検出器応答関数Wm はm番目のフレ
ームの光学システム点分布関数PSFm のコンボリュー
ションであり、次式のように表される。
用の光学システム点分布関数PSFm のコンボリューシ
ョンであり、個々の検出器応答関数Wm はm番目のフレ
ームの光学システム点分布関数PSFm のコンボリュー
ションであり、次式のように表される。
【0046】
【数3】 ★Gm はGm により示される列マトリックスにより表さ
れることができ、各エレメントは焦点平面アレイの各検
出器測定に対応する。Hm はHm により示されている2
次元マトリックスにより表され、エレメントは電子光セ
ンサのシステム応答関数に対応する。同様に、情景はS
m により示される列マトックスとして表され、超解像度
のグリッドへ投影されるm番目の情景実体のサンプルに
対応する素子Sx を有する。
れることができ、各エレメントは焦点平面アレイの各検
出器測定に対応する。Hm はHm により示されている2
次元マトリックスにより表され、エレメントは電子光セ
ンサのシステム応答関数に対応する。同様に、情景はS
m により示される列マトックスとして表され、超解像度
のグリッドへ投影されるm番目の情景実体のサンプルに
対応する素子Sx を有する。
【0047】超スピン開口プロセス処理の1実施形態
は、最大事後確率評価と通常呼ばれているバイエル(Ba
yes )評価方法を使用する。最大事後評価は測定を行う
目標物体の条件付き確率を最大にする未知の目標物体の
値を求める。超スピン開口に対する最大事後方法は1組
の回転開口センサ測定を行う超解像度の情景の確率を最
大にする超解像度の情景の評価を模索する。超スピン開
口のこの実施形態では、式2で与えられる画像モデルは
次式により特定化される。
は、最大事後確率評価と通常呼ばれているバイエル(Ba
yes )評価方法を使用する。最大事後評価は測定を行う
目標物体の条件付き確率を最大にする未知の目標物体の
値を求める。超スピン開口に対する最大事後方法は1組
の回転開口センサ測定を行う超解像度の情景の確率を最
大にする超解像度の情景の評価を模索する。超スピン開
口のこの実施形態では、式2で与えられる画像モデルは
次式により特定化される。
【0048】
【数4】 ★ここで情景Sは回転開口測定間隔期間の単一の情景実
体であると仮定されている。
体であると仮定されている。
【0049】バイエルの法則を回転開口測定の全ての組
に適用すると、Sの測定値を与える条件付き確率の式が
導かれる。
に適用すると、Sの測定値を与える条件付き確率の式が
導かれる。
【0050】
【数5】 ★Gは全測定マトリックスであり、これは個々のフレー
ム測定マトリックスからなるブロック列のマトリックス
である。
ム測定マトリックスからなるブロック列のマトリックス
である。
【0051】
【数6】 ★ここでnf は捕捉されるフレームの総数である。Sは
超解像度グリッドに投影される情景サンプルのマトリッ
クスである。
超解像度グリッドに投影される情景サンプルのマトリッ
クスである。
【0052】情景放射または各画素素子に対する光子の
反射が統計的に独立していると仮定すると、1次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
反射が統計的に独立していると仮定すると、1次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
【0053】
【数7】 ★ソース光子変動により左右される情景統計では、情景
画素素子の確率密度関数はポアッソン(Poisson )の分
布式であると仮定され、
画素素子の確率密度関数はポアッソン(Poisson )の分
布式であると仮定され、
【0054】
【数8】 ★ここでkは情景の投影された画素素子を示す指数であ
る。
る。
【0055】統計的に独立したフレーム測定では、情景
が与えられた全測定ベクトルのジョイントの確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
が与えられた全測定ベクトルのジョイントの確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
【0056】
【数9】 ★同様に、統計的に独立した検出器測定では、情景が与
えられたフレーム測定の発生についての条件付き確率は
情景上で条件付けされた個々の画素確率密度の積として
表されることができる。
えられたフレーム測定の発生についての条件付き確率は
情景上で条件付けされた個々の画素確率密度の積として
表されることができる。
【0057】
【数10】 ★全体の測定の条件付き確率密度はそれ故、情景上で条
件付けされた画素およびフレーム測定確率の積により与
えられる。
件付けされた画素およびフレーム測定確率の積により与
えられる。
【0058】
【数11】 ★個々のベクトル測定がポアッソンプロセスであると仮
定すると、個々の測定の条件付き確率は次式のように与
えられることができる。
定すると、個々の測定の条件付き確率は次式のように与
えられることができる。
【0059】
【数12】 ★ここでΣk Hm (j−k)S(k)は個々の検出器測
定の平均値である。条件付き確率密度を最大にする情景
を発見することは確率密度の対数を最大にする情景を発
見することに等しいので、次式を最大にするSに対する
解が求められる。
定の平均値である。条件付き確率密度を最大にする情景
を発見することは確率密度の対数を最大にする情景を発
見することに等しいので、次式を最大にするSに対する
解が求められる。
【0060】
【数13】 ★積の対数が対数の和であることを認識すると、Sが与
えられた測定マトリックスGの確率は次式のようにな
る。
えられた測定マトリックスGの確率は次式のようにな
る。
【0061】
【数14】 ★
【0062】
【数15】 ★
【0063】(13)を(12)に代入すると、以下の
スターリング(Stirling)の概算式 ln[x!]=xln[x]−x (14) をS(k)!に適用すると次式が得られる。
スターリング(Stirling)の概算式 ln[x!]=xln[x]−x (14) をS(k)!に適用すると次式が得られる。
【0064】
【数16】 ★条件付き確率を最大にするS(l)値を決定するた
め、ln[p(S|G)]の導関数がゼロになる値S
(l)が必要とされる。ln[p(S|G)]の導関数
は次式により与えられる。
め、ln[p(S|G)]の導関数がゼロになる値S
(l)が必要とされる。ln[p(S|G)]の導関数
は次式により与えられる。
【0065】
【数17】 ★導関数をゼロに設定すると、Sを解く非線形の式が得
られる。
られる。
【0066】
【数18】 ★この式はHm を指数関数化し係数化した後、次式のよ
うになる。
うになる。
【0067】
【数19】 ★ピカール(Picard)法として知られる反復プロセスは
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用される。反復プロセスは以下の式により示され
る。
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用される。反復プロセスは以下の式により示され
る。
【0068】
【数20】 ★ここで上付き文字[n]と[n+1]はそれぞれ現在
および次の情景評価を表している。
および次の情景評価を表している。
【0069】超解像度の回転開口が元来の測定よりも高
い空間周波数に対応するグリッド点を評価するので、反
復プロセスは、超解像度の座標グリッドと、測定(また
は中間基準)グリッドとの間で再度マッピングすること
を含んでいる。反復的アルゴリズムの1実施形態はブロ
ック平均により超解像度グリッドのダウンサンプリング
を行い、複製により測定グリッドのアップサンプリング
を行うことによって、測定座標グリッドと超解像度の座
標グリッドとの間のサンプル間隔の差に適応する。代わ
りに、アップサンプリングとダウンサンプリングはより
高い次数の補間技術を使用して実行されることができ
る。反復的解は次式により与えられる。
い空間周波数に対応するグリッド点を評価するので、反
復プロセスは、超解像度の座標グリッドと、測定(また
は中間基準)グリッドとの間で再度マッピングすること
を含んでいる。反復的アルゴリズムの1実施形態はブロ
ック平均により超解像度グリッドのダウンサンプリング
を行い、複製により測定グリッドのアップサンプリング
を行うことによって、測定座標グリッドと超解像度の座
標グリッドとの間のサンプル間隔の差に適応する。代わ
りに、アップサンプリングとダウンサンプリングはより
高い次数の補間技術を使用して実行されることができ
る。反復的解は次式により与えられる。
【0070】
【数21】 ★上向きの矢印でアップサンプリングを、下向きの矢印
でダウンサンプリングを示す。ここで指数関数の和は次
式を使用して簡単にされる。
でダウンサンプリングを示す。ここで指数関数の和は次
式を使用して簡単にされる。
【0071】
【数22】 ★現在の評価値/S(l)[n] の平均値は現在の評価値
S(l)[n] に置き換えられる。スケール係数はコンバ
ージェンスを制御するため指数関数に含まれてもよい。
S(l)[n] に置き換えられる。スケール係数はコンバ
ージェンスを制御するため指数関数に含まれてもよい。
【0072】先に説明した超スピン開口の実施形態は情
景Sを評価するため最大事後評価方法を使用する。別の
実施形態は情景を表すためにサンプリング統計を使用す
る。式(1)と(2)のSm とGm はランダムフィール
ド(ランダム変数の二次元アレイ)であるので、これら
のランダムフィールドの平均値を解析することが可能で
ある。測定全体にわたるサンプル平均を示すため< >
を使用すると、コンボリューションの平均値は次式のよ
うに表されることができる。
景Sを評価するため最大事後評価方法を使用する。別の
実施形態は情景を表すためにサンプリング統計を使用す
る。式(1)と(2)のSm とGm はランダムフィール
ド(ランダム変数の二次元アレイ)であるので、これら
のランダムフィールドの平均値を解析することが可能で
ある。測定全体にわたるサンプル平均を示すため< >
を使用すると、コンボリューションの平均値は次式のよ
うに表されることができる。
【0073】
【数23】 ★ここで情景実体全体のサンプル平均は以下のように表
される。
される。
【0074】
【数24】 ★ここでnf は測定された回転開口のフレームの数であ
る。先の実施形態のように反復的評価プロセスを得るた
めに最大事後方程式は公式化され、ピカール技術が適用
される。
る。先の実施形態のように反復的評価プロセスを得るた
めに最大事後方程式は公式化され、ピカール技術が適用
される。
【0075】この実施形態ではGは式6で定められ、S
は同様に情景実体で限定される。
は同様に情景実体で限定される。
【0076】情景放射または各画素素子に対する光子の
反射が統計的に独立し、実体Sm が測定フレームにわた
って統計的に独立しているものと仮定すると、1次元情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができる。即ち、
反射が統計的に独立し、実体Sm が測定フレームにわた
って統計的に独立しているものと仮定すると、1次元情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができる。即ち、
【0077】
【数25】 ★ソース光子の変動により左右される情景統計では、情
景画素素子の確率密度関数は次式のポアッソン(Poisso
n )の分布式であると仮定され、
景画素素子の確率密度関数は次式のポアッソン(Poisso
n )の分布式であると仮定され、
【0078】
【数26】 ★kは情景の投影された画素素子を示す指数である。全
体の平均値/S(k)はサンプル平均式3aにより近似
され、次式が得られる。
体の平均値/S(k)はサンプル平均式3aにより近似
され、次式が得られる。
【0079】
【数27】 ★ダミー指数pは特定のフレームmからのフレームにわ
たるカウントを弁別するために使用されている。
たるカウントを弁別するために使用されている。
【0080】統計と独立したフレーム測定では、情景を
与えられた全体の測定ベクトルのジョイント確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
与えられた全体の測定ベクトルのジョイント確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
【0081】
【数28】 ★同様に、統計と独立したベクトル測定では、情景を与
えられたフレーム測定の発生の条件付き確率は情景につ
いて条件付けされた個々の画素確率密度の積として表さ
れることができる。
えられたフレーム測定の発生の条件付き確率は情景につ
いて条件付けされた個々の画素確率密度の積として表さ
れることができる。
【0082】
【数29】 ★それ故、全体の測定の条件付き確率密度は情景につい
て条件付けされた画素およびフレーム測定確率の積によ
り与えられる。
て条件付けされた画素およびフレーム測定確率の積によ
り与えられる。
【0083】
【数30】 ★個々の検出器の測定がポアッソンプロセスであると仮
定すると、個々の測定の条件付き確率は次式により表さ
れることができる。
定すると、個々の測定の条件付き確率は次式により表さ
れることができる。
【0084】
【数31】 ★ここで/G(j)は個々の検出器測定の全体的平均値
である。全体的平均値は式3のサンプル平均値<G
m (j)>により評価され、以下の測定確率密度が導か
れる。
である。全体的平均値は式3のサンプル平均値<G
m (j)>により評価され、以下の測定確率密度が導か
れる。
【0085】
【数32】 ★条件付き確率密度を最大にする情景を発見することは
確率密度の対数を最大にする情景を発見することに等し
いので、次式を最大にするSの解が求められる。
確率密度の対数を最大にする情景を発見することに等し
いので、次式を最大にするSの解が求められる。
【0086】
【数33】 ★
【0087】積の対数が対数の和であることを認識する
と、Sが与えられる測定マトリックスGの確率は次式の
ように表されることができる。
と、Sが与えられる測定マトリックスGの確率は次式の
ように表されることができる。
【0088】
【数34】 ★
【0089】式(25)を式(24)に代入すると、以
下の式が導かれる。
下の式が導かれる。
【0090】
【数35】 ★
【0091】以下のスターリングの概算式 ln[x!]=xln[x]−x (26) をSm (k)!に適用すると次式が得られる。
【0092】
【数36】 ★
【0093】反復的な手順が得られ、ここでq番目の実
体Sq が評価され、この評価値は<S>の評価値を更新
するために使用される。条件付き確率を最大にする値S
q (l)を決定するため、ln[p(S|G)]の導関
数がゼロになる値Sq (l)が必要とされる。ln[p
(S|G)]の導関数は次式により与えられる。
体Sq が評価され、この評価値は<S>の評価値を更新
するために使用される。条件付き確率を最大にする値S
q (l)を決定するため、ln[p(S|G)]の導関
数がゼロになる値Sq (l)が必要とされる。ln[p
(S|G)]の導関数は次式により与えられる。
【0094】
【数37】 ★
【0095】導関数をゼロに設定すると、Sq を解く次
の非線形の式が得られる。
の非線形の式が得られる。
【0096】
【数38】 ★この式はHm を指数関数化し係数化した後、次式のよ
うになる。
うになる。
【0097】
【数39】 ★
【0098】ピカール法として知られる反復プロセスは
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用されることができる。反復プロセスは以下の式
により示される。
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用されることができる。反復プロセスは以下の式
により示される。
【0099】
【数40】 ★上向きの矢印はアップサンプリングを、下向きの矢印
はダウンサンプリングを示す。ここで上付き文字[n]
と[n+1]はそれぞれ現在および次の情景評価を表し
ている。先の実施形態のように、上下の矢印で示される
アップサンプリングとダウンサンプリングは測定座標グ
リッドと超解像度の座標グリッドとの間のサンプル間隔
の差に適応するために使用されることができる。
はダウンサンプリングを示す。ここで上付き文字[n]
と[n+1]はそれぞれ現在および次の情景評価を表し
ている。先の実施形態のように、上下の矢印で示される
アップサンプリングとダウンサンプリングは測定座標グ
リッドと超解像度の座標グリッドとの間のサンプル間隔
の差に適応するために使用されることができる。
【0100】全体の測定に適応する単一の確率密度を使
用する最大事後方法を以下説明する。先に示したアルゴ
リズム手順は基本的に決定しなければならない2つの未
知数と、超解像度の情景画素と超解像度の情景画素平均
値を有する。第1の組の超解像度情景評価に関するアル
ゴリズムの実行は情景画素平均値を評価するための付属
プロセスを伴う。第2の組の超解像度情景評価式に関す
るアルゴリズムの実行は各情景画素の全体平均を概算す
るためにサンプル統計を使用する。
用する最大事後方法を以下説明する。先に示したアルゴ
リズム手順は基本的に決定しなければならない2つの未
知数と、超解像度の情景画素と超解像度の情景画素平均
値を有する。第1の組の超解像度情景評価に関するアル
ゴリズムの実行は情景画素平均値を評価するための付属
プロセスを伴う。第2の組の超解像度情景評価式に関す
るアルゴリズムの実行は各情景画素の全体平均を概算す
るためにサンプル統計を使用する。
【0101】代わりの方法は定数である独立したランダ
ム変数の総数を認識し、即ち全ての合計された測定の総
数または合計された部分的測定の総数は個々の測定で条
件付けされたランダム変数である。
ム変数の総数を認識し、即ち全ての合計された測定の総
数または合計された部分的測定の総数は個々の測定で条
件付けされたランダム変数である。
【0102】それ故、合計された測定の総数と1より小
さい1組の測定が与えられた超解像度の情景の条件付き
確率を最大にするバイエル(Bayes )法則が求められ
る。代わりのプロセスの体系はここで説明した第1の実
施形態と同一である。
さい1組の測定が与えられた超解像度の情景の条件付き
確率を最大にするバイエル(Bayes )法則が求められ
る。代わりのプロセスの体系はここで説明した第1の実
施形態と同一である。
【0103】光子変動により左右される情景および測定
統計では、個々の情景および測定素子の確率密度関数は
ポアッソンであると仮定され、即ち次式のようになる。
統計では、個々の情景および測定素子の確率密度関数は
ポアッソンであると仮定され、即ち次式のようになる。
【0104】
【数41】 ★ここで/S(k↑)は情景のk↑番目の統計された
“画素”の平均値であり、Gm (k↓)はフレームmか
らのj番目の検出器測定である。
“画素”の平均値であり、Gm (k↓)はフレームmか
らのj番目の検出器測定である。
【0105】平均値/X1 と個々のポアッソン分布を有
するn個の統計的に独立したポアッソンランダム変数X
1 …Xn では以下のようになる。
するn個の統計的に独立したポアッソンランダム変数X
1 …Xn では以下のようになる。
【0106】
【数42】 ★ポアッソンランダム変数の合計は次式のようになり、
【0107】
【数43】 ★この式はF. A. Haight氏による文献Handbook of the
Poisson Distributionの71頁で与えられる分布を有する
ポアッソンランダム変数である。
Poisson Distributionの71頁で与えられる分布を有する
ポアッソンランダム変数である。
【0108】
【数44】 ★ここで、
【0109】
【数45】 ★
【0110】平均値/X1 …/Xn を有するn個の統計
的に独立したポアッソンランダム変数X1 …Xn では、
個々のポアッソン分布p(Xk )と、X1 ,…Xn-2 ,
Xn- 1 の条件付き確率分布と、以下の合計されたランダ
ム変数の総数、即ち、
的に独立したポアッソンランダム変数X1 …Xn では、
個々のポアッソン分布p(Xk )と、X1 ,…Xn-2 ,
Xn- 1 の条件付き確率分布と、以下の合計されたランダ
ム変数の総数、即ち、
【0111】
【数46】 ★は多項式であり、F. A. Haight氏による文献Handbook
of the Poisson Distributionの71頁により与えられて
いる。
of the Poisson Distributionの71頁により与えられて
いる。
【0112】
【数47】 ★
【0113】式(28)を1組の超スピン開口測定に適
用することは、条件付き確率に明白に存在しない測定を
選択することに依存する。測定のナンバリングは任意で
あるので、合計された測定の総数が与えられる条件付き
確率は次式のように表される
用することは、条件付き確率に明白に存在しない測定を
選択することに依存する。測定のナンバリングは任意で
あるので、合計された測定の総数が与えられる条件付き
確率は次式のように表される
【0114】。
【数48】 ★ここで、Nが単位フレーム当たりの検出器数−1であ
るならば、M=フレーム数であり、Nが単位フレーム当
たりの検出器数であるならば、M=フレーム数−1であ
り、次式が成り立つ。
るならば、M=フレーム数であり、Nが単位フレーム当
たりの検出器数であるならば、M=フレーム数−1であ
り、次式が成り立つ。
【0115】
【数49】 ★
【0116】データの任意のサブセットに関しては、次
式が成り立つ。
式が成り立つ。
【0117】
【数50】 ★独立したランダム変数の総数が一定であり、即ち合計
された全ての測定の総数または合計された測定部分の総
数が個々の測定で条件付けされたランダム変数であるこ
とを認識することは、超解像度情景の評価を発見するた
め別の変形された自己一貫最大事後方法を設定するため
の手段を与える。全体の測定と1より小さい1組の測定
が与えられる超解像度の情景のバイエルの評価規則は次
式のようになる。
された全ての測定の総数または合計された測定部分の総
数が個々の測定で条件付けされたランダム変数であるこ
とを認識することは、超解像度情景の評価を発見するた
め別の変形された自己一貫最大事後方法を設定するため
の手段を与える。全体の測定と1より小さい1組の測定
が与えられる超解像度の情景のバイエルの評価規則は次
式のようになる。
【0118】
【数51】 ★ここでG_は測定の1より小さい総セットに対応する
ブロックマトリックスであり、Sは超解像度グリッドに
投影される情景サンプルのマトリックスである。マトリ
ックスG_の素子は以下の組に対応する。
ブロックマトリックスであり、Sは超解像度グリッドに
投影される情景サンプルのマトリックスである。マトリ
ックスG_の素子は以下の組に対応する。
【0119】
【数52】 ★ここで、N=単位フレーム当たりの検出器数−1であ
るならば、Mはフレーム数であり、またはN=単位フレ
ーム当たりの検出器数であるならば、Mはフレーム数−
1である。
るならば、Mはフレーム数であり、またはN=単位フレ
ーム当たりの検出器数であるならば、Mはフレーム数−
1である。
【0120】自己一貫の必要性は、個々の検出器測定の
平均値を仮定することにより再度取入れられ、全体の回
転開口測定の和の平均値は雑音のない組の測定に対応
し、即ち次式のようになる。
平均値を仮定することにより再度取入れられ、全体の回
転開口測定の和の平均値は雑音のない組の測定に対応
し、即ち次式のようになる。
【0121】
【数53】 ★
【0122】↑と↓は超解像度および測定座標の空間グ
リッドを示している。
リッドを示している。
【0123】画像センサ測定が情景から伝播される光子
の検出に対応するので、情景の評価値は正でなければな
らない。それ故、付加的な自己一貫の要求は、 S≧0 (31)であり、これは超解像度アルゴ
リズムの実施形態により不明瞭に満足される。
の検出に対応するので、情景の評価値は正でなければな
らない。それ故、付加的な自己一貫の要求は、 S≧0 (31)であり、これは超解像度アルゴ
リズムの実施形態により不明瞭に満足される。
【0124】確率密度の最大値の決定はln[p(S|
G_,Gtot )]の最大値の決定と等価である。
G_,Gtot )]の最大値の決定と等価である。
【0125】個々の測定は統計的に独立したポアッソン
ランダム変数であると仮定され、合計された全体の測定
の全てまたはその一部の合計値も同様に個々の測定の平
均値はポアッソンランダム変数であると仮定される。さ
らに、ポアッソンランダム変数の平均値は雑音のない測
定であると仮定され、即ち合計された総回転開口測定の
平均値は以下のようになり、
ランダム変数であると仮定され、合計された全体の測定
の全てまたはその一部の合計値も同様に個々の測定の平
均値はポアッソンランダム変数であると仮定される。さ
らに、ポアッソンランダム変数の平均値は雑音のない測
定であると仮定され、即ち合計された総回転開口測定の
平均値は以下のようになり、
【0126】
【数54】 ★個々の検出器測定の平均値は次式のようになる。
【0127】
【数55】 ★
【0128】それ故、合計された全体の測定が与えられ
ている測定の確率の多項分布に平均値を代入すると、情
景の出現で条件付けられた対応する確率が得られ、即
ち、
ている測定の確率の多項分布に平均値を代入すると、情
景の出現で条件付けられた対応する確率が得られ、即
ち、
【0129】
【数56】 ★積の対数は対数の合計であり、情景で条件付けされた
測定の対数は次式のようになる。 ln[p(S|G_,Gtot )]=ln[p(G_,G
tot IS)]+ln[p(S)]−ln[p(G_,G
tot )] スターリング近似を超解像度情景の投影された画素の統
計に適用し、積の対数が対数の合計であることを認識
し、適切な確率分布をバイエルの定理に代入すると次式
が得られる。
測定の対数は次式のようになる。 ln[p(S|G_,Gtot )]=ln[p(G_,G
tot IS)]+ln[p(S)]−ln[p(G_,G
tot )] スターリング近似を超解像度情景の投影された画素の統
計に適用し、積の対数が対数の合計であることを認識
し、適切な確率分布をバイエルの定理に代入すると次式
が得られる。
【0130】
【数57】 ★超スピン開口の最大事後の実施形態の最大が条件付き
確率を最大にする超解像度の情景を発見することに依存
するので、条件付き確率の対数の微分から生じる非線形
の部分的微分方程式を解く情景評価が求められ、即ち、
次式が得られる。
確率を最大にする超解像度の情景を発見することに依存
するので、条件付き確率の対数の微分から生じる非線形
の部分的微分方程式を解く情景評価が求められ、即ち、
次式が得られる。
【0131】
【数58】 ★
【0132】先の方法に関しては、超解像度情景の解
は、情景の超解像度画素に関する条件付き確率の導関数
がゼロになることを必要とすることによって生成される
1組の非線形の結合式を満足しなければならない。この
セクションで説明している実施形態の1組の連立方程式
を以下に示す。
は、情景の超解像度画素に関する条件付き確率の導関数
がゼロになることを必要とすることによって生成される
1組の非線形の結合式を満足しなければならない。この
セクションで説明している実施形態の1組の連立方程式
を以下に示す。
【0133】
【数59】 ★
【0134】1組の連立方程式の数値解は指数関数とピ
カールの反復的方法を適用することにより実行されるこ
とができる。
カールの反復的方法を適用することにより実行されるこ
とができる。
【0135】合計された測定の超解像度の実施形態を説
明する。この超解像度の回転開口実施形態では付加的な
超解像度式の組が回転開口フレーム中の測定された総光
子の数に基づいて設定される。評価された情景値が、測
定された回転開口フレームが与えられた情景値に対する
条件付き確率密度の極値であるのと同様に、全体の測定
の発生が与えられた情景値に対する条件付き確率密度の
極値であることを必要とすることにより超解像度の情景
値は決定される。2つの組の式は未知の平均画素値を消
去するために使用されることができ、超解像度の情景画
素値用の1組の非線形の式を残す。1組の非線形の式は
ピカールの方法を用いて解かれる。この代わりの超解像
度手段の定義に導かれる手順の詳細について以下説明す
る。
明する。この超解像度の回転開口実施形態では付加的な
超解像度式の組が回転開口フレーム中の測定された総光
子の数に基づいて設定される。評価された情景値が、測
定された回転開口フレームが与えられた情景値に対する
条件付き確率密度の極値であるのと同様に、全体の測定
の発生が与えられた情景値に対する条件付き確率密度の
極値であることを必要とすることにより超解像度の情景
値は決定される。2つの組の式は未知の平均画素値を消
去するために使用されることができ、超解像度の情景画
素値用の1組の非線形の式を残す。1組の非線形の式は
ピカールの方法を用いて解かれる。この代わりの超解像
度手段の定義に導かれる手順の詳細について以下説明す
る。
【0136】雑音のない測定を実施する別の方法を説明
する。回転開口の超解像度アルゴリズムの導出の第1の
バージョンに関しては、雑音のない画像測定はセンサシ
ステム応答関数Hm のコンボリューションとして表さ
れ、観察された情景Sまたは次式を有する。
する。回転開口の超解像度アルゴリズムの導出の第1の
バージョンに関しては、雑音のない画像測定はセンサシ
ステム応答関数Hm のコンボリューションとして表さ
れ、観察された情景Sまたは次式を有する。
【0137】
【数60】 ★m番目の雑音のない画像測定のディスクリートな形態
と、全てのフレームに対する全ての検出器の雑音のない
全体の測定は次式で表される。
と、全てのフレームに対する全ての検出器の雑音のない
全体の測定は次式で表される。
【0138】
【数61】 ★
【0139】先の式では、jは検出器位置ベクトル/r
j に関連する一次元指数であり、kは検出器位置ベクト
ル/rk に関連する一次元指数である。
j に関連する一次元指数であり、kは検出器位置ベクト
ル/rk に関連する一次元指数である。
【0140】イソプラネティック光学システムと空間的
に均一な応答焦点平面アレイ検出器に対して、m番目の
フレームシステム応答関数は個々の検出器応答関数Wm
を有する光学システムの点分布関数PSFm のコンボリ
ューションとして表され、即ち次式のようになる。
に均一な応答焦点平面アレイ検出器に対して、m番目の
フレームシステム応答関数は個々の検出器応答関数Wm
を有する光学システムの点分布関数PSFm のコンボリ
ューションとして表され、即ち次式のようになる。
【0141】
【数62】 ★
【0142】変更された最大事後方法は、超解像度の情
景の評価値を発見するため2つの異なったバイエルの確
率規則を使用する。2つのバイエル規則は次式であり、
景の評価値を発見するため2つの異なったバイエルの確
率規則を使用する。2つのバイエル規則は次式であり、
【0143】
【数63】 ★この式は1組の回転開口測定全体が与えられる超解像
度情景が出現する確率と、全ての回転開口測定の和が与
えられる超解像度情景が出現する確率にそれぞれ対応す
る。先に定めたように、Gは全体の測定マトリックスで
あり、これは個々のフレーム測定マトリックスからなる
ブロックマトリックスであり、Sは超解像度座標グリッ
ドに投影される個々の情景素子からなる列マトリックス
であり、即ち、
度情景が出現する確率と、全ての回転開口測定の和が与
えられる超解像度情景が出現する確率にそれぞれ対応す
る。先に定めたように、Gは全体の測定マトリックスで
あり、これは個々のフレーム測定マトリックスからなる
ブロックマトリックスであり、Sは超解像度座標グリッ
ドに投影される個々の情景素子からなる列マトリックス
であり、即ち、
【0144】
【数64】 ★
【0145】ここでnf は測定フレームの総数であり、
npixel は超解像度情景の評価の超解像度点の総数であ
る。
npixel は超解像度情景の評価の超解像度点の総数であ
る。
【0146】情景が与えられた全体の測定に対する条件
付き確率を有する自己一貫の超解像度方法は、情景Sに
対する評価値を求め、これは回転開口測定の完全な組が
与えられる情景の条件付き確率を最大にし、全ての回転
開口測定の和が与えられる情景の条件付き確率を最大に
し、即ち、同時にp(S|G)とp(S|Gtot )を最
大にするSの解が求められる。一般的な分布では、両者
の条件を満たす極値の解は可能ではないが、以下のセク
ションで説明する手順および統計的な必要条件は非線形
の1組の方程式に関する解を許容する。
付き確率を有する自己一貫の超解像度方法は、情景Sに
対する評価値を求め、これは回転開口測定の完全な組が
与えられる情景の条件付き確率を最大にし、全ての回転
開口測定の和が与えられる情景の条件付き確率を最大に
し、即ち、同時にp(S|G)とp(S|Gtot )を最
大にするSの解が求められる。一般的な分布では、両者
の条件を満たす極値の解は可能ではないが、以下のセク
ションで説明する手順および統計的な必要条件は非線形
の1組の方程式に関する解を許容する。
【0147】自己一貫の要求は個々の検出器測定の平均
値を仮定することにより組入れられ、回転開口測定の和
の合計の平均値は雑音のない組の測定に対応し、即ち、
値を仮定することにより組入れられ、回転開口測定の和
の合計の平均値は雑音のない組の測定に対応し、即ち、
【0148】
【数65】 ★
【0149】↑と↓は超解像度と“測定”座標空間グリ
ッドを示す。
ッドを示す。
【0150】画像センサ測定は情景から伝播する光子の
検出に対応し、情景の評価値は常に正である。それ故、
付加的な自己一貫の要求は、 S≧0, (34)であり、これは前述の超解像
度アルゴリズムにより暗黙的に満たされる。
検出に対応し、情景の評価値は常に正である。それ故、
付加的な自己一貫の要求は、 S≧0, (34)であり、これは前述の超解像
度アルゴリズムにより暗黙的に満たされる。
【0151】対数は単調な関数であるので、確率密度の
最大値の決定はln[p(SIG)]とln[p(SI
Gtot )]の確率密度の最大値を決定することに等し
い。対数をバイエルの定理に適用すると、次式の最大値
が求められる。
最大値の決定はln[p(SIG)]とln[p(SI
Gtot )]の確率密度の最大値を決定することに等し
い。対数をバイエルの定理に適用すると、次式の最大値
が求められる。
【0152】
【数66】 ★
【0153】統計的に独立したフレーム測定に対して、
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
【0154】
【数67】 ★
【0155】同様に、統計的に独立した検出器測定に対
して、情景が与えられるフレーム測定の発生の条件付き
確率は情景で条件付けされる個々の画素確率密度の積と
して表されることができる。
して、情景が与えられるフレーム測定の発生の条件付き
確率は情景で条件付けされる個々の画素確率密度の積と
して表されることができる。
【0156】
【数68】 ★
【0157】全体の測定の条件付き確率密度はそれ故、
情景について条件付けされる画素とフレーム測定確率の
積により与えられる。
情景について条件付けされる画素とフレーム測定確率の
積により与えられる。
【0158】
【数69】 ★
【0159】情景放射または各画素素子に対する光子の
反射が統計的に独立していると仮定すると、一次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
反射が統計的に独立していると仮定すると、一次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
【0160】
【数70】 ★
【0161】ソースの光子変動により左右される情景統
計では、情景画素素子の確率密度関数は次式のように仮
定される。
計では、情景画素素子の確率密度関数は次式のように仮
定される。
【0162】
【数71】 ★ここで/S(k↑)はk↓番目の投影された“画素”
の平均値である。
の平均値である。
【0163】個々の検出器測定もポアッソンランダムプ
ロセスであると仮定すると、個々の検出器測定確率密度
関数は次式のようになることが示唆される。
ロセスであると仮定すると、個々の検出器測定確率密度
関数は次式のようになることが示唆される。
【0164】
【数72】 ★m番目のフレーム測定j↓番目の“画素”の平均値を
置換すると、
置換すると、
【0165】
【数73】 ★この式は情景を与える検出器測定の確率の式を導出す
る。
る。
【0166】
【数74】 ★
【0167】統計的に独立したポアッソンランダム変数
の和もポアッソンランダム変数であるので、合計された
全ての回転開口検出器確率の全体はポアッソンランダム
変数であり、次式により与えられる。
の和もポアッソンランダム変数であるので、合計された
全ての回転開口検出器確率の全体はポアッソンランダム
変数であり、次式により与えられる。
【0168】
【数75】 ★合計された回転開口測定の全体の平均値は次式のよう
になり、
になり、
【0169】
【数76】 ★これをポアッソン確率密度に代入すると、所定のSが
発生するGtot の条件付き確率が得られる。
発生するGtot の条件付き確率が得られる。
【0170】
【数77】 ★
【0171】積の対数は対数の和であるので、情景で条
件付けされた測定の対数は以下のようになる。
件付けされた測定の対数は以下のようになる。
【0172】
【数78】 ★
【0173】同様に、情景確率密度の対数は以下のよう
になる。
になる。
【0174】
【数79】 ★
【0175】先の式を情景が与えられる1組の情景測定
に対する条件付き確率に代入すると次式が得られる。
に対する条件付き確率に代入すると次式が得られる。
【0176】
【数80】 ★
【0177】p(S|Gtot )の対数の決定は情景が与
えられる合計された全体の測定の条件付き確率の対数を
決定することに依存する。
えられる合計された全体の測定の条件付き確率の対数を
決定することに依存する。
【0178】
【数81】 ★
【0179】条件付き確率の関係に対してバイエルの法
則へ代入すると、次式が得られる。
則へ代入すると、次式が得られる。
【0180】
【数82】 ★
【0181】約50を越える光子をカウントした場合、
スターリング近似である。
スターリング近似である。
【0182】ln[x!]=xln[x]−xの式は計
算手順を簡単にするために使用される。問題の2つの条
件付き確率に適用されるスターリング概算は次式を導
く。
算手順を簡単にするために使用される。問題の2つの条
件付き確率に適用されるスターリング概算は次式を導
く。
【0183】
【数83】 ★
【0184】所望の情景画素素子S(l↑)に関して適
切な確率の対数を微分することにより、確率を最大にす
る情景の解を要求する極値方程式が得られる。微分によ
って次式が得られる。
切な確率の対数を微分することにより、確率を最大にす
る情景の解を要求する極値方程式が得られる。微分によ
って次式が得られる。
【0185】
【数84】 ★
【0186】クロネッカー(Kronecker )のデルタ関数
指標積の合計を計算することにより式を簡単にすると、
指標積の合計を計算することにより式を簡単にすると、
【0187】
【数85】 ★
【0188】測定が与えられる情景確率の対数の導関数
の設定に関連する1組の式を以下示す。
の設定に関連する1組の式を以下示す。
【0189】
【数86】 ★
【0190】S(l↑)を解くことによって、1組の結
合された非線形の式が得られる。
合された非線形の式が得られる。
【0191】
【数87】 ★
【0192】指数関数化によって次式が得られる。
【0193】
【数88】 ★
【0194】合計された全体の測定が与えられる情景確
率の対数の導関数を設定することに関する1組の式を以
下示す。
率の対数の導関数を設定することに関する1組の式を以
下示す。
【0195】
【数89】 ★
【0196】S(l↑)の対数を解くと、次の1組の結
合された非線形の式が得られる。
合された非線形の式が得られる。
【0197】
【数90】 ★
【0198】指数関数化によって次式が得られる。
【0199】
【数91】 ★
【0200】2つの組の式(35)、(36)はS(l
↑)または/S(l↑)の評価プロセスを公式化するた
めに同時に解かれることができる。
↑)または/S(l↑)の評価プロセスを公式化するた
めに同時に解かれることができる。
【0201】以上説明した超スピン開口の実施形態は最
大事後評価技術を適用する。超スピン開口プロセスに対
する別の方法は最大尤度評価と通常呼ばれるバイエルの
評価方法を使用する。最大尤度評価は目標が与えられる
測定の条件付き確率を最大にする未知の目標の値を求め
る。即ち、目標の最大尤度評価は1組の観察された回転
開口センサ測定全体に導かれる可能性が最も高い目標で
ある。
大事後評価技術を適用する。超スピン開口プロセスに対
する別の方法は最大尤度評価と通常呼ばれるバイエルの
評価方法を使用する。最大尤度評価は目標が与えられる
測定の条件付き確率を最大にする未知の目標の値を求め
る。即ち、目標の最大尤度評価は1組の観察された回転
開口センサ測定全体に導かれる可能性が最も高い目標で
ある。
【0202】前述したように、Gは全測定マトリックス
であり、これは個々のフレーム測定マトリックスからな
るブロック列マトリックスであり、以下に示す。
であり、これは個々のフレーム測定マトリックスからな
るブロック列マトリックスであり、以下に示す。
【0203】
【数92】 ★ここでnf は捕捉されるフレームの総数である。Sは
ランダム情景のnf 実体に関するブロック列マトリック
スである。
ランダム情景のnf 実体に関するブロック列マトリック
スである。
【0204】
【数93】 ★
【0205】情景Sがランダムフィールドであるとする
と、情景平均値の最大尤度評価が計算される。最大尤度
評価は情景平均/Sの値であり、これは尤度関数P(G
|/S)を最大にする。
と、情景平均値の最大尤度評価が計算される。最大尤度
評価は情景平均/Sの値であり、これは尤度関数P(G
|/S)を最大にする。
【0206】統計的に独立したフレーム測定に対して、
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
【0207】
【数94】 ★
【0208】同様に、統計的に独立した検出器測定で
は、情景が与えられるフレーム発生の条件付き確率は情
景について条件付けされた個々の画素確率密度の積とし
て表されることができる。
は、情景が与えられるフレーム発生の条件付き確率は情
景について条件付けされた個々の画素確率密度の積とし
て表されることができる。
【0209】
【数95】 ★
【0210】全体の測定条件付き確率密度はそれ故、情
景について条件付けされた画素とフレーム測定確率の積
により与えられる。
景について条件付けされた画素とフレーム測定確率の積
により与えられる。
【0211】
【数96】 ★
【0212】個々の検出器測定がポアッソンプロセスで
あると仮定すると、個々の測定の条件付き確率は次式の
ように表される。
あると仮定すると、個々の測定の条件付き確率は次式の
ように表される。
【0213】
【数97】 ★
【0214】条件付き確率密度を最大にする情景を発見
することは、確率密度の対数を最大にする情景を発見す
ることに等しいので、次式を最大にする/Sの解が求め
られる。
することは、確率密度の対数を最大にする情景を発見す
ることに等しいので、次式を最大にする/Sの解が求め
られる。
【0215】
【数98】 ★積の対数が対数の和であることを認識すると、/Sが
与えられる測定マトリックスGの確率は以下のようにな
る。
与えられる測定マトリックスGの確率は以下のようにな
る。
【0216】
【数99】 ★
【0217】式(43)を式(42)に代入すると、次
式が得られる。
式が得られる。
【0218】
【数100】 ★条件付き確率を最大にする/S(l)値を決定するた
め、ln[p(G|/S)]の導関数がゼロになる値/
S(l)が要求される。ln[p(G|/S)]の導関
数は次式により与えられる。
め、ln[p(G|/S)]の導関数がゼロになる値/
S(l)が要求される。ln[p(G|/S)]の導関
数は次式により与えられる。
【0219】
【数101】 ★
【0220】導関数をゼロに設定すると、/Sを解く非
線形の式が得られる。
線形の式が得られる。
【0221】
【数102】 ★
【0222】光学システムの特性を保存する光子によ
り、ほとんどの構成例では、式の左側はnf に合計され
る。ピカールの技術を適用するため、式(44)の各側
は/S(l)により乗算され、評価は次式のように反復
して形成される。
り、ほとんどの構成例では、式の左側はnf に合計され
る。ピカールの技術を適用するため、式(44)の各側
は/S(l)により乗算され、評価は次式のように反復
して形成される。
【0223】
【数103】 ★ここで上付き文字nは反復数を示している。
【0224】超スピン開口の最大尤度評価を計算するこ
の方法は予期最大アルゴリズムである。予期最大アルゴ
リズムは情景平均の最大尤度評価へ集中する。
の方法は予期最大アルゴリズムである。予期最大アルゴ
リズムは情景平均の最大尤度評価へ集中する。
【0225】以上説明したアルゴリズム手順の各反復で
は、乗算の補正係数はnf 回転開口測定の全体の組から
計算される。反復予期最大アルゴリズムは最大尤度評価
を計算するために使用されてもよい。反復アルゴリズム
は個々に各フレームに基づく評価値を更新し、それによ
って最大尤度評価を決定する計算上の必要性をなくす。
は、乗算の補正係数はnf 回転開口測定の全体の組から
計算される。反復予期最大アルゴリズムは最大尤度評価
を計算するために使用されてもよい。反復アルゴリズム
は個々に各フレームに基づく評価値を更新し、それによ
って最大尤度評価を決定する計算上の必要性をなくす。
【0226】予期される最大アルゴリズムの反復形態
は、超スピン開口に適用するためTitterington氏による
文献Journal of the Royal Statistical Society-B、46
巻、No.2、257 〜267 頁、1984年により展開された方法
に発展される。以下の反復的予期最大アルゴリズムはn
=0,1,…,、m=0,1,…,Nf −1、l=1,
2,…,Npixel で得られる。
は、超スピン開口に適用するためTitterington氏による
文献Journal of the Royal Statistical Society-B、46
巻、No.2、257 〜267 頁、1984年により展開された方法
に発展される。以下の反復的予期最大アルゴリズムはn
=0,1,…,、m=0,1,…,Nf −1、l=1,
2,…,Npixel で得られる。
【0227】
【数104】 ★
【0228】この反復的な予期される最大アルゴリズム
は先に説明した方法を用いて初期化されるか、情景の平
均評価の各素子を回転開口フレームの光子の平均数に基
づいて一定値に初期化される。
は先に説明した方法を用いて初期化されるか、情景の平
均評価の各素子を回転開口フレームの光子の平均数に基
づいて一定値に初期化される。
【0229】
【数105】 ★
【0230】前述の超スピン開口の実施形態は未知の情
景を評価するための2つの広い解決方法、即ち最大事後
評価と、最大尤度評価を含む。最大事後評価のカテゴリ
ー内には4つの実施形態が開示されており、それらは、
情景の現在の評価に等しい情景平均を評価することによ
り情景平均/S(l)が特別な方法で評価される最大事
後方法と、サンプル統計を用いて情景平均/S(l)が
評価され、情景Sq (l)の各実体を評価するための反
復アルゴリズムを生じ、それに続いてサンプル平均の再
評価が行われる最大事後方法と、単一の確率密度関数が
測定と1組の測定された回転開口フレームの総光子とで
公式化される最大事後方法と、別々の確率密度関数が1
組の測定された回転開口フレームの総光子で公式化さ
れ、それに続いて2つの組の非線形の式を同時的に解く
最大事後方法とを含んでいる。
景を評価するための2つの広い解決方法、即ち最大事後
評価と、最大尤度評価を含む。最大事後評価のカテゴリ
ー内には4つの実施形態が開示されており、それらは、
情景の現在の評価に等しい情景平均を評価することによ
り情景平均/S(l)が特別な方法で評価される最大事
後方法と、サンプル統計を用いて情景平均/S(l)が
評価され、情景Sq (l)の各実体を評価するための反
復アルゴリズムを生じ、それに続いてサンプル平均の再
評価が行われる最大事後方法と、単一の確率密度関数が
測定と1組の測定された回転開口フレームの総光子とで
公式化される最大事後方法と、別々の確率密度関数が1
組の測定された回転開口フレームの総光子で公式化さ
れ、それに続いて2つの組の非線形の式を同時的に解く
最大事後方法とを含んでいる。
【0231】最大尤度評価のカテゴリー内では2つの実
施形態が説明されている。第1の実施形態は情景の全体
平均/S(l)の評価を形成するために1組の測定され
た回転開口フレームを使用し、ここでは回転開口フレー
ムの組全体は評価のそれぞれの更新で使用される。第2
の最大尤度実施例では、反復的アルゴリズムが使用さ
れ、ここでは情景の全体平均は、一度に1つの回転開口
フレーム測定を用いて更新される。各最大尤度実施形態
は予期最大手順であり、それ故ピカール反復が最大尤度
評価へ集中することを保証する。
施形態が説明されている。第1の実施形態は情景の全体
平均/S(l)の評価を形成するために1組の測定され
た回転開口フレームを使用し、ここでは回転開口フレー
ムの組全体は評価のそれぞれの更新で使用される。第2
の最大尤度実施例では、反復的アルゴリズムが使用さ
れ、ここでは情景の全体平均は、一度に1つの回転開口
フレーム測定を用いて更新される。各最大尤度実施形態
は予期最大手順であり、それ故ピカール反復が最大尤度
評価へ集中することを保証する。
【0232】説明を完全にするため、図17は放射状の
試験パターンの超解像度の画像を生成するために使用さ
れる超スピン開口合成方法30の実施結果を示している。
図18は超スピン開口合成動作の例を示している。動作
の定量的測定は目標の既知の基準パワースペクトルと超
解像度画像のパワースペクトルとの間の相関係数を決定
することにより得られる。図19は真の情景が水平バー
ターゲットである超スピン開口合成の例を示している。
動作の定量的測定は物体の既知の基準パワースペクトル
と超解像度画像のパワースペクトルとの間の相関係数を
決定することにより発見される。
試験パターンの超解像度の画像を生成するために使用さ
れる超スピン開口合成方法30の実施結果を示している。
図18は超スピン開口合成動作の例を示している。動作
の定量的測定は目標の既知の基準パワースペクトルと超
解像度画像のパワースペクトルとの間の相関係数を決定
することにより得られる。図19は真の情景が水平バー
ターゲットである超スピン開口合成の例を示している。
動作の定量的測定は物体の既知の基準パワースペクトル
と超解像度画像のパワースペクトルとの間の相関係数を
決定することにより発見される。
【0233】以上、回転ストリップ(部分的)開口画像
放射線センサ装置と、複数の回転ストリップ開口画像測
定からの超解像度情景評価を合成する方法を説明した。
所望の実施形態は本発明の原理の応用を表した多数の特
定の実施形態のうちの数例にすぎないことが理解されよ
う。明らかに、多数のおよび他の装置が本発明の技術範
囲を逸脱することなく当業者により容易に行われること
ができる。
放射線センサ装置と、複数の回転ストリップ開口画像測
定からの超解像度情景評価を合成する方法を説明した。
所望の実施形態は本発明の原理の応用を表した多数の特
定の実施形態のうちの数例にすぎないことが理解されよ
う。明らかに、多数のおよび他の装置が本発明の技術範
囲を逸脱することなく当業者により容易に行われること
ができる。
【図1】本発明の原理による回転開口画像放射線センサ
システムの概略図。
システムの概略図。
【図2】本発明の原理によるシステムおよび方法の第1
の実施形態の概略図。
の実施形態の概略図。
【図3】本発明の第2の実施形態の概略図。
【図4】本発明の第3の実施形態の概略図。
【図5】本発明の第4の実施形態の概略図。
【図6】本発明の第5の実施形態の概略図。
【図7】図2乃至6の本発明に関連する超解像度処理方
法の実施形態の概略図。
法の実施形態の概略図。
【図8】図2乃至6の本発明に関連する超解像度処理方
法の実施形態の概略図。
法の実施形態の概略図。
【図9】図2乃至6の本発明に関連する超解像度処理方
法の別の実施形態の概略図。
法の別の実施形態の概略図。
【図10】図2乃至6の本発明に関連する超解像度処理
方法の別の実施形態の概略図。
方法の別の実施形態の概略図。
【図11】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図12】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図13】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図14】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図15】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図16】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図17】試験パターンに関する本発明を用いた試験結
果の図。
果の図。
【図18】別の試験パターンに関する本発明を用いた試
験結果の図。
験結果の図。
【図19】さらに別の試験パターンに関する本発明を用
いた試験結果の図。
いた試験結果の図。
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成9年4月24日
【手続補正1】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正内容】
【書類名】 明細書
【発明の名称】 回転ストリップ開口映像測定を使
用する超解像度全円形開口情景合成
用する超解像度全円形開口情景合成
【特許請求の範囲】
【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、回転開口放射線測
定装置および放射線検出方法に関し、特に、合成された
超解像度情景が複数の回転ストリップ開口映像の測定か
ら評価される回転ストリップ(部分的)開口映像生成用
放射線測定装置および方法に関する。
定装置および放射線検出方法に関し、特に、合成された
超解像度情景が複数の回転ストリップ開口映像の測定か
ら評価される回転ストリップ(部分的)開口映像生成用
放射線測定装置および方法に関する。
【0002】
【従来の技術】宇宙空間で使用するプラットホームから
高解像度の映像を得るために、例えば通常のセンサ技術
では大型で、重く、展開可能な光学システムの能動制御
装置を使用している。要求される使命および1次反射鏡
の大きさに応じて、センサを構成する全ての光学素子の
位置、傾斜、形状、および整列の制御のために1次反射
鏡セグメントの周期的な位置および傾斜の制御が行われ
る。
高解像度の映像を得るために、例えば通常のセンサ技術
では大型で、重く、展開可能な光学システムの能動制御
装置を使用している。要求される使命および1次反射鏡
の大きさに応じて、センサを構成する全ての光学素子の
位置、傾斜、形状、および整列の制御のために1次反射
鏡セグメントの周期的な位置および傾斜の制御が行われ
る。
【0003】本発明と同じ解像度を有する全開口システ
ムは比較的大きい開口面積のために多量の集光能力を有
している。しかしながら、軌道上に数十メートルの直径
の開口を配置するためには、本発明と競合する全開口シ
ステムは次のようなことが必要になる。すなわち光学シ
ステムの直径が宇宙船の発射装置の直径よりも大きくな
るとセグメント化された光学表面および折畳まれた支持
構造が必要になる。また、重量を節減するために薄い変
形可能な反射鏡が使用される場合には複雑で潜在的に高
い帯域幅の適応光学技術が必要である。さらに、位相ア
レイとして構成される場合には複雑な位置および瞳孔整
合制御が要求される。それ故、全開口システムは本発明
と比較して比較的重く、コストが高く、また高い技術的
リスクを有している。
ムは比較的大きい開口面積のために多量の集光能力を有
している。しかしながら、軌道上に数十メートルの直径
の開口を配置するためには、本発明と競合する全開口シ
ステムは次のようなことが必要になる。すなわち光学シ
ステムの直径が宇宙船の発射装置の直径よりも大きくな
るとセグメント化された光学表面および折畳まれた支持
構造が必要になる。また、重量を節減するために薄い変
形可能な反射鏡が使用される場合には複雑で潜在的に高
い帯域幅の適応光学技術が必要である。さらに、位相ア
レイとして構成される場合には複雑な位置および瞳孔整
合制御が要求される。それ故、全開口システムは本発明
と比較して比較的重く、コストが高く、また高い技術的
リスクを有している。
【0004】本発明に関係する従来技術は米国特許第5
243351号明細書に記載されており、スピンApシ
ステムと呼ばれている。スピンApシステムと本発明の
方法(超スピンApシステムと呼ぶ)との共通性は、映
像の合成または情景の評価のために時間的に整合したス
トリップ開口を使用する点にある。
243351号明細書に記載されており、スピンApシ
ステムと呼ばれている。スピンApシステムと本発明の
方法(超スピンApシステムと呼ぶ)との共通性は、映
像の合成または情景の評価のために時間的に整合したス
トリップ開口を使用する点にある。
【0005】
【発明が解決しようとする課題】それ故、本発明の目的
は、合成された超解像度情景が複数の回転ストリップ開
口映像測定から評価される回転ストリップ開口映像放射
線測定装置および測定方法を提供することである。
は、合成された超解像度情景が複数の回転ストリップ開
口映像測定から評価される回転ストリップ開口映像放射
線測定装置および測定方法を提供することである。
【0006】
【課題を解決するための手段】上記の、およびその他の
目的を達成するために、本発明の1形態においては、合
成された超解像度情景が複数の回転ストリップ開口映像
測定から評価される回転ストリップ開口映像放射線測定
装置および測定方法が提供される。ここで使用される超
解像度とは開口を限定している空間周波数カットオフよ
り大きい空間周波数に対して情景情報を正確に評価する
能力を言う。特にここに記載したセンサシステムおよび
データ処理方法に対して、超解像度とはストリップ開口
システムの光学的カットオフを越えて評価された情景の
情報内容を拡張すること、および/または回転する開口
の形状と関係する最大の相関長に等しい半径を有する等
価全円形開口の光学的カットオフを越えて拡張すること
を意味する。表現“ストリップ開口”は米国特許第52
43351号明細書に記載されたような一般的な開口形
状を意味するものである。
目的を達成するために、本発明の1形態においては、合
成された超解像度情景が複数の回転ストリップ開口映像
測定から評価される回転ストリップ開口映像放射線測定
装置および測定方法が提供される。ここで使用される超
解像度とは開口を限定している空間周波数カットオフよ
り大きい空間周波数に対して情景情報を正確に評価する
能力を言う。特にここに記載したセンサシステムおよび
データ処理方法に対して、超解像度とはストリップ開口
システムの光学的カットオフを越えて評価された情景の
情報内容を拡張すること、および/または回転する開口
の形状と関係する最大の相関長に等しい半径を有する等
価全円形開口の光学的カットオフを越えて拡張すること
を意味する。表現“ストリップ開口”は米国特許第52
43351号明細書に記載されたような一般的な開口形
状を意味するものである。
【0007】本発明のシステムは時間的に連続する、或
いは順次の映像を生成する回転ストリップ開口望遠鏡を
備えている。回転ストリップ開口望遠鏡は典型的に回転
ストリップ開口1次反射器および2次反射器を備えてい
る。2次元の検出器アレイが設けられて望遠鏡の焦点平
面に位置する映像を検出する。回転補償装置が設けられ
て、アレイの検出器の積分時間中の回転により軌跡が描
かれるのを阻止するために使用されている。信号プロセ
ッサが設けられて、望遠鏡の光学軸を中心としてストリ
ップ開口が回転するとき望遠鏡により映像の形成された
目標の情景の複数の映像フレームを記録し、回転するス
トリップ開口の空間周波数カットオフより大きい空間周
波数、および/または回転開口システムの全円形開口の
光学的カットオフを越える記録された映像から観察され
た情景の超解像度評価を合成する。
いは順次の映像を生成する回転ストリップ開口望遠鏡を
備えている。回転ストリップ開口望遠鏡は典型的に回転
ストリップ開口1次反射器および2次反射器を備えてい
る。2次元の検出器アレイが設けられて望遠鏡の焦点平
面に位置する映像を検出する。回転補償装置が設けられ
て、アレイの検出器の積分時間中の回転により軌跡が描
かれるのを阻止するために使用されている。信号プロセ
ッサが設けられて、望遠鏡の光学軸を中心としてストリ
ップ開口が回転するとき望遠鏡により映像の形成された
目標の情景の複数の映像フレームを記録し、回転するス
トリップ開口の空間周波数カットオフより大きい空間周
波数、および/または回転開口システムの全円形開口の
光学的カットオフを越える記録された映像から観察され
た情景の超解像度評価を合成する。
【0008】したがって、本発明は、複数の回転ストリ
ップ開口映像測定からの超解像度放射線測定情景評価を
合成する回転ストリップ(部分的)開口映像放射線測定
装置を提供する。本発明はしたがって、同じ重量の従来
のスピンApシステムと比較したとき、或いは同じ重量
の従来の全開口システムと比較したとき、改善された高
い解像度の映像を提供する。
ップ開口映像測定からの超解像度放射線測定情景評価を
合成する回転ストリップ(部分的)開口映像放射線測定
装置を提供する。本発明はしたがって、同じ重量の従来
のスピンApシステムと比較したとき、或いは同じ重量
の従来の全開口システムと比較したとき、改善された高
い解像度の映像を提供する。
【0009】上述の説明を要約すると、本発明のセンサ
およびプロセッサによって行われる合成処理の1実施形
態は以下のとおりである。望遠鏡の光学軸を中心として
ストリップ開口が回転するとき、次のような動作が行わ
れる。回転補償装置は検出器の積分時間中逆回転してそ
れによって時間的に静止している映像が与えられる。映
像フレームは記録され、保存される。回転する(情景に
対して)検出器アレイ構造が選択されるならば、捕捉さ
れたフレームは座標変換され、(超解像度ではない)合
成された映像の基準グリッドに補間される。データはフ
ーリエ変換され、記憶される。新しいフレームは記録お
よび保存される。ランダムな照準線エラーにより記録さ
れたデータのフレーム毎の不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換は照準線エラーに
対して補正され、記憶される。上述のステップは各スト
リップ開口映像フレームに対して順次反復され、或いは
フレームは順次捕捉され、記憶され、その後グローバル
な照準線の評価がフレームを整合するために決定され
る。
およびプロセッサによって行われる合成処理の1実施形
態は以下のとおりである。望遠鏡の光学軸を中心として
ストリップ開口が回転するとき、次のような動作が行わ
れる。回転補償装置は検出器の積分時間中逆回転してそ
れによって時間的に静止している映像が与えられる。映
像フレームは記録され、保存される。回転する(情景に
対して)検出器アレイ構造が選択されるならば、捕捉さ
れたフレームは座標変換され、(超解像度ではない)合
成された映像の基準グリッドに補間される。データはフ
ーリエ変換され、記憶される。新しいフレームは記録お
よび保存される。ランダムな照準線エラーにより記録さ
れたデータのフレーム毎の不整合の評価が得られる。ス
トリップ開口映像またはフーリエ変換は照準線エラーに
対して補正され、記憶される。上述のステップは各スト
リップ開口映像フレームに対して順次反復され、或いは
フレームは順次捕捉され、記憶され、その後グローバル
な照準線の評価がフレームを整合するために決定され
る。
【0010】個々のストリップ開口フレームが整合され
たならば、超解像度合成処理が監視される。超解像度合
成処理の1実施形態では以下のようなステップが行われ
る。通過帯域拡張の所望の量は合成された超解像度映像
に対するサンプル間隔を決定するために使用される。座
標系および超解像度映像と関連する空間グリッドは高い
解像度または超解像度グリッドと呼ばれる。超解像度合
成映像の最初の評価および高解像度グリッドにおける超
解像度映像の平均値が決定される。超解像度映像に対す
る最初の評価はデータに接近するスピンAp処理を適用
し、低雑音サンプルを得るためにローパス空間周波数フ
ィルタ処理を行い、全開口伝達関数に対して補償し、画
素複製により超解像度グリッドに評価を拡張することに
よって得られることができる。超解像度映像の平均値に
対する最初の評価は同様の方法で得ることができる。
たならば、超解像度合成処理が監視される。超解像度合
成処理の1実施形態では以下のようなステップが行われ
る。通過帯域拡張の所望の量は合成された超解像度映像
に対するサンプル間隔を決定するために使用される。座
標系および超解像度映像と関連する空間グリッドは高い
解像度または超解像度グリッドと呼ばれる。超解像度合
成映像の最初の評価および高解像度グリッドにおける超
解像度映像の平均値が決定される。超解像度映像に対す
る最初の評価はデータに接近するスピンAp処理を適用
し、低雑音サンプルを得るためにローパス空間周波数フ
ィルタ処理を行い、全開口伝達関数に対して補償し、画
素複製により超解像度グリッドに評価を拡張することに
よって得られることができる。超解像度映像の平均値に
対する最初の評価は同様の方法で得ることができる。
【0011】一度開始点が設定されると非線形超解像度
および相互作用超解像度技術が現在の評価の非線形関数
と超解像度情景の現在の評価との積を形成することによ
り超解像度情景の次の評価が形成される。非線形関数は
測定されたフレームの数ニヨッテスケールされた複数の
回転測定の指数積である。指数の論点は測定、ストリッ
プ開口センサの応答関数、および現在の超解像度情景評
価の関数である。各回転位置に関連する指数の論点は補
助関数により特定の回転位置に対する反射されたストリ
ップ開口の応答関数のコンボリューションを実行するこ
とによって得られる。補助関数発生の1素子は、特定の
回転位置に対するストリップ開口の空間応答関数を有す
る超解像度情景の現在の評価のコンボリューションによ
り画素測定を分割することからなる。
および相互作用超解像度技術が現在の評価の非線形関数
と超解像度情景の現在の評価との積を形成することによ
り超解像度情景の次の評価が形成される。非線形関数は
測定されたフレームの数ニヨッテスケールされた複数の
回転測定の指数積である。指数の論点は測定、ストリッ
プ開口センサの応答関数、および現在の超解像度情景評
価の関数である。各回転位置に関連する指数の論点は補
助関数により特定の回転位置に対する反射されたストリ
ップ開口の応答関数のコンボリューションを実行するこ
とによって得られる。補助関数発生の1素子は、特定の
回転位置に対するストリップ開口の空間応答関数を有す
る超解像度情景の現在の評価のコンボリューションによ
り画素測定を分割することからなる。
【0012】本発明の放射線センサシステムおよび映像
合成処理方法は特定の使命の要求および技術的妥協に依
存するから、放射線センサシステムおよび映像合成処理
方法は、フレーム間の照準線のランダムおよび/または
系統的なドリフトを補償する手段、およびイソプラナテ
ィック光学系の点の広がる作用、映像の運動および検出
器の有限の電子的帯域幅によるフィールド点と関係のな
い映像の軌跡、映像の補償されない回転運動により生じ
たフィールド点と関係する映像の軌跡のような、過去お
よび後の知識によるエラー源を補償する手段を具備する
こともできる。情景構造に関する仮定を実効するための
マルコフ(Markov )ランダムフィールドの使用、全エ
ネルギ測定基準の使用、別の条件確率密度関数の使用の
ような過去および後の知識を伴うことにより、前述の説
明した処理技術の率直な変形が得られる。本発明のシス
テムおよび方法において使用される補償技術の性能は過
去および後の知識の正確度および実効信号対雑音比に依
存する。
合成処理方法は特定の使命の要求および技術的妥協に依
存するから、放射線センサシステムおよび映像合成処理
方法は、フレーム間の照準線のランダムおよび/または
系統的なドリフトを補償する手段、およびイソプラナテ
ィック光学系の点の広がる作用、映像の運動および検出
器の有限の電子的帯域幅によるフィールド点と関係のな
い映像の軌跡、映像の補償されない回転運動により生じ
たフィールド点と関係する映像の軌跡のような、過去お
よび後の知識によるエラー源を補償する手段を具備する
こともできる。情景構造に関する仮定を実効するための
マルコフ(Markov )ランダムフィールドの使用、全エ
ネルギ測定基準の使用、別の条件確率密度関数の使用の
ような過去および後の知識を伴うことにより、前述の説
明した処理技術の率直な変形が得られる。本発明のシス
テムおよび方法において使用される補償技術の性能は過
去および後の知識の正確度および実効信号対雑音比に依
存する。
【0013】本発明のシステムおよび方法は、回転スト
リップ開口を備えた光学系から超解像度映像を提供する
ことのできる回転開口形式の改良された宇宙船搭載光学
映像センサおよび処理構造を提供する。本発明の映像測
定および合成手順は、リスクが低く、コストが安く、軽
量で、製造が簡単であり、また宇宙船搭載映像装置用と
して複雑で大きい開口の適応光学系を展開するのに比較
して簡単な展開構造を与える利点がある。
リップ開口を備えた光学系から超解像度映像を提供する
ことのできる回転開口形式の改良された宇宙船搭載光学
映像センサおよび処理構造を提供する。本発明の映像測
定および合成手順は、リスクが低く、コストが安く、軽
量で、製造が簡単であり、また宇宙船搭載映像装置用と
して複雑で大きい開口の適応光学系を展開するのに比較
して簡単な展開構造を与える利点がある。
【0014】本発明のシステムおよび方法は、特に望遠
鏡の重量がペイロードの重量の大きい割合を占めている
場合に衛星軌道から高い解像度の映像を提供する。衛星
の発射のコストは主として軌道中の重量の関数である。
動作システムの状態に応じて、所定の許容可能な望遠鏡
の重量に対して、本発明は従来のスピンApシステムよ
り50%まで解像度を増加させることができる。本発明
は高解像度(大きい開口と等価)の宇宙船搭載観察装置
としてそのようなシステムで使用されることができる。
鏡の重量がペイロードの重量の大きい割合を占めている
場合に衛星軌道から高い解像度の映像を提供する。衛星
の発射のコストは主として軌道中の重量の関数である。
動作システムの状態に応じて、所定の許容可能な望遠鏡
の重量に対して、本発明は従来のスピンApシステムよ
り50%まで解像度を増加させることができる。本発明
は高解像度(大きい開口と等価)の宇宙船搭載観察装置
としてそのようなシステムで使用されることができる。
【0015】本発明の超スピンApシステムは利用でき
る測定の使用により全開口映像を生成し、ストリップ開
口と関連した最大の相関長と関連した空間周波数カット
オフを越えた、すなわち、等価なスピンApシステム
(帯域制限されている)の合成全開口空間周波数カット
オフを越えた空間周波数カットオフを有する。本発明の
超スピンApシステムは最大の相関長を有するスピンA
pシステムから発生された映像を表示する。それ故所定
の実効合成開口寸法に対する低いペイロード重量が得ら
れる変形されたスピンAp映像処理方法が得られる利点
がある。さらに、リスクが低く、コストが安く、軽量
で、製造が簡単であり、また宇宙船搭載映像装置用とし
て複雑で大きい全円形開口の適応光学系を展開するもの
(または位相アレイ)のに比較して簡単な展開構造を与
える利点がある。さらに、超解像度評価技術はストリッ
プ開口の通過帯域の外側の情景情報の評価を可能にする
から、超スピンAp処理手順はより少ない回転サンプル
によるスピンApと同じ合成映像品質を提供することが
できる。
る測定の使用により全開口映像を生成し、ストリップ開
口と関連した最大の相関長と関連した空間周波数カット
オフを越えた、すなわち、等価なスピンApシステム
(帯域制限されている)の合成全開口空間周波数カット
オフを越えた空間周波数カットオフを有する。本発明の
超スピンApシステムは最大の相関長を有するスピンA
pシステムから発生された映像を表示する。それ故所定
の実効合成開口寸法に対する低いペイロード重量が得ら
れる変形されたスピンAp映像処理方法が得られる利点
がある。さらに、リスクが低く、コストが安く、軽量
で、製造が簡単であり、また宇宙船搭載映像装置用とし
て複雑で大きい全円形開口の適応光学系を展開するもの
(または位相アレイ)のに比較して簡単な展開構造を与
える利点がある。さらに、超解像度評価技術はストリッ
プ開口の通過帯域の外側の情景情報の評価を可能にする
から、超スピンAp処理手順はより少ない回転サンプル
によるスピンApと同じ合成映像品質を提供することが
できる。
【0016】本発明の種々の特徴および利点は、添付図
面に関連した以下の詳細な説明に関して一層容易に理解
され、添付図面において、同一の参照番号は同一の構造
素子を示している。
面に関連した以下の詳細な説明に関して一層容易に理解
され、添付図面において、同一の参照番号は同一の構造
素子を示している。
【0017】
【発明の実施の形態】米国特許第5,243,351 号明細書の
図1に対応する図1を参照すると、本発明の原理によっ
て変更された回転開口映像生成放射線センサシステム10
の一例が示されている。米国特許第5,243,351 号の内容
が本明細書において参照文献とされている。回転開口映
像生成放射線センサシステム10は、超解像度の情景の評
価を合成し、他方で、照準線ジッタを取除き、同じ重量
の通常の光学システムおよび同じ重量の回転開口システ
ムとの比較において改良された高解像度の映像を供給す
るように構成されている。回転開口映像生成放射線セン
サシステム10は、1次反射器12a および2次反射器12b
を含む回転ストリップ開口望遠鏡11を具備している。所
定の環境の下では3次反射器(図示されていない)が望
遠鏡11において使用されることがある。
図1に対応する図1を参照すると、本発明の原理によっ
て変更された回転開口映像生成放射線センサシステム10
の一例が示されている。米国特許第5,243,351 号の内容
が本明細書において参照文献とされている。回転開口映
像生成放射線センサシステム10は、超解像度の情景の評
価を合成し、他方で、照準線ジッタを取除き、同じ重量
の通常の光学システムおよび同じ重量の回転開口システ
ムとの比較において改良された高解像度の映像を供給す
るように構成されている。回転開口映像生成放射線セン
サシステム10は、1次反射器12a および2次反射器12b
を含む回転ストリップ開口望遠鏡11を具備している。所
定の環境の下では3次反射器(図示されていない)が望
遠鏡11において使用されることがある。
【0018】本発明を開示する目的のために、システム
10は、地球指向アンテナ13a を有する静止セクション13
を具備している衛星の形態で示されている。望遠鏡11
は、プラットフォーム14上に配置され、そこには静止セ
クション13も結合されている。回転開口映像生成放射線
センサシステム10は、(ストリップ開口12を有してい
る)1次反射器12a が望遠鏡11の光学軸の周囲を回転す
る際に望遠鏡11によって映された複数のターゲット情景
の映像フレームを記録するように構成されている。照準
線安定反射器15および映像回転除去装置16は、検出器ア
レイ17によって感知される前に映像の安定化および回転
除去を行うように構成された望遠鏡11の光路に沿って配
置されている。回転除去装置16は、回転補償制御装置19
の制御の下で、検出器アレイ17を構成している検出器の
積分時間中に映像を逆に回転させ、それによって静止映
像を供給する。照準線安定反射器15は、システム10の軌
道力学のために、照準線エラーの大きい帯域幅ならびに
照準線エラーを取除くために(信号プロセッサあるいは
別の専用制御システム等によって供給された)照準線制
御システムによって使用される。ターゲット情景は、映
像フレームを処理するように調整された信号プロセッサ
20に結合された望遠鏡11の焦点平面に位置された検出器
アレイ17上に映像が形成される。その代りに、信号プロ
セッサ20が専用の基礎プロセッサを具備していることも
ある。望遠鏡11、検出器アレイ17、および関連したハー
ドウェアは、本明細書においてセンサシステム21と呼ば
れる。センサシステム21によって生成された個々の映像
フレームは、信号プロセッサ20において処理および結合
され、それによって、本発明のシステム10および方法30
に従って超解像度情景評価を合成し、それは特に図2乃
至11を参照にして説明されている。
10は、地球指向アンテナ13a を有する静止セクション13
を具備している衛星の形態で示されている。望遠鏡11
は、プラットフォーム14上に配置され、そこには静止セ
クション13も結合されている。回転開口映像生成放射線
センサシステム10は、(ストリップ開口12を有してい
る)1次反射器12a が望遠鏡11の光学軸の周囲を回転す
る際に望遠鏡11によって映された複数のターゲット情景
の映像フレームを記録するように構成されている。照準
線安定反射器15および映像回転除去装置16は、検出器ア
レイ17によって感知される前に映像の安定化および回転
除去を行うように構成された望遠鏡11の光路に沿って配
置されている。回転除去装置16は、回転補償制御装置19
の制御の下で、検出器アレイ17を構成している検出器の
積分時間中に映像を逆に回転させ、それによって静止映
像を供給する。照準線安定反射器15は、システム10の軌
道力学のために、照準線エラーの大きい帯域幅ならびに
照準線エラーを取除くために(信号プロセッサあるいは
別の専用制御システム等によって供給された)照準線制
御システムによって使用される。ターゲット情景は、映
像フレームを処理するように調整された信号プロセッサ
20に結合された望遠鏡11の焦点平面に位置された検出器
アレイ17上に映像が形成される。その代りに、信号プロ
セッサ20が専用の基礎プロセッサを具備していることも
ある。望遠鏡11、検出器アレイ17、および関連したハー
ドウェアは、本明細書においてセンサシステム21と呼ば
れる。センサシステム21によって生成された個々の映像
フレームは、信号プロセッサ20において処理および結合
され、それによって、本発明のシステム10および方法30
に従って超解像度情景評価を合成し、それは特に図2乃
至11を参照にして説明されている。
【0019】本発明による超スピンApシステム10およ
び方法30の特定の実施形態が以下に説明される。特に、
図2乃至11において、本発明の原理にしたがった超ス
ピンAp処理方法30の種々の実施形態が示されており、
それらは、本発明の原理を実行するために適切に変更さ
れた図1に示された回転開口放射線センサシステムおよ
び方法によって使用される。観測シナリオおよび時間
線、使用できる従来技術の知識の量、および処理の流れ
における使用可能な計算スループットに基づいて特定の
実施形態が選択される。
び方法30の特定の実施形態が以下に説明される。特に、
図2乃至11において、本発明の原理にしたがった超ス
ピンAp処理方法30の種々の実施形態が示されており、
それらは、本発明の原理を実行するために適切に変更さ
れた図1に示された回転開口放射線センサシステムおよ
び方法によって使用される。観測シナリオおよび時間
線、使用できる従来技術の知識の量、および処理の流れ
における使用可能な計算スループットに基づいて特定の
実施形態が選択される。
【0020】図2において、本発明によるシステム10お
よび方法30の第1の実施形態が示されており、そこにお
いて、最高レベルの超スピンAp映像合成方法30が信号
プロセッサ20において実行される。回転開口センサ21
は、一組のストリップ開口映像22を捕捉および生成し、
それは超スピンApプロセッサ23に移送される。超スピ
ンApプロセッサ23は、超解像度情景評価24を合成す
る。
よび方法30の第1の実施形態が示されており、そこにお
いて、最高レベルの超スピンAp映像合成方法30が信号
プロセッサ20において実行される。回転開口センサ21
は、一組のストリップ開口映像22を捕捉および生成し、
それは超スピンApプロセッサ23に移送される。超スピ
ンApプロセッサ23は、超解像度情景評価24を合成す
る。
【0021】図3を参照すると、第2のシステム10およ
び方法30が示されている。回転開口センサ21は一組のス
トリップ開口映像22を捕捉および生成し、それらはフレ
ーム整合エラープロセッサ25に移される。フレーム整合
エラープロセッサ25は、その組の個々のフレーム間の残
りの整合エラーを決定する。例えば、フレームの整合は
以下のように達成される。回転ストリップ開口が望遠鏡
11の光学軸の周囲を回転しているとき、以下のように動
作する。回転除去装置16は、検出器アレイ17の積分時間
中、逆に回転し、それによって時間的に静止映像を供給
する。映像フレームは記録され、蓄えられる。(情景に
関して)回転する検出器アレイ構造が選択された場合、
捕捉されたフレームは、調和して変換され、(超解像度
でない)合成映像の基準グリッドに補間される。データ
はフーリエ変換され、記憶される。新しいフレームが記
録され、記憶される。ランダムな照準線エラーのために
生じた記録されたデータのフレーム毎の不整合の評価が
得られる。ストリップ開口映像、すなわちフーリエ変換
は、それらの照準線エラーを補正され、記憶される。前
述のステップは、各ストリップ開口映像フレーム毎に連
続的に反復され、そうでなければ、フレームは順次連続
的に捕捉および記憶され、その後、フレームを整合させ
るために照準線の全体的な評価が決定される。整合プロ
セッサ25は、一組の補正されたフレーム、あるいは一組
の補正されていないフレームおよびフレーム間の整合エ
ラー評価、もしくはその両方を超解像度プロセッサ23に
転送する。補正されたフレーム、あるいは補正されてい
ないフレームおよびそれに対応する一組の整合エラーを
使用して、超解像度評価プロセッサ23は、超解像度情景
評価24を生成する。
び方法30が示されている。回転開口センサ21は一組のス
トリップ開口映像22を捕捉および生成し、それらはフレ
ーム整合エラープロセッサ25に移される。フレーム整合
エラープロセッサ25は、その組の個々のフレーム間の残
りの整合エラーを決定する。例えば、フレームの整合は
以下のように達成される。回転ストリップ開口が望遠鏡
11の光学軸の周囲を回転しているとき、以下のように動
作する。回転除去装置16は、検出器アレイ17の積分時間
中、逆に回転し、それによって時間的に静止映像を供給
する。映像フレームは記録され、蓄えられる。(情景に
関して)回転する検出器アレイ構造が選択された場合、
捕捉されたフレームは、調和して変換され、(超解像度
でない)合成映像の基準グリッドに補間される。データ
はフーリエ変換され、記憶される。新しいフレームが記
録され、記憶される。ランダムな照準線エラーのために
生じた記録されたデータのフレーム毎の不整合の評価が
得られる。ストリップ開口映像、すなわちフーリエ変換
は、それらの照準線エラーを補正され、記憶される。前
述のステップは、各ストリップ開口映像フレーム毎に連
続的に反復され、そうでなければ、フレームは順次連続
的に捕捉および記憶され、その後、フレームを整合させ
るために照準線の全体的な評価が決定される。整合プロ
セッサ25は、一組の補正されたフレーム、あるいは一組
の補正されていないフレームおよびフレーム間の整合エ
ラー評価、もしくはその両方を超解像度プロセッサ23に
転送する。補正されたフレーム、あるいは補正されてい
ないフレームおよびそれに対応する一組の整合エラーを
使用して、超解像度評価プロセッサ23は、超解像度情景
評価24を生成する。
【0022】図4を参照すると、第3のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に送られる。フレーム整
合エラープロセッサ25は、その組の個々のフレームの間
の残りの整合エラーを決定する。整合プロセッサ25は、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を回転開口帯域限定映像合成プロセッサ
26および超解像度評価プロセッサ23に送る。回転開口帯
域限定映像合成プロセッサ26は、米国特許第5,243,351
号明細書に記載された回転開口処理によって全円形開口
映像を表す映像を生成する。一組の補正されたフレー
ム、あるいは一組の補正されていないフレームおよびフ
レーム間の整合エラー評価、もしくはその両方、ならび
に回転開口帯域限定映像を使用して、超解像度評価プロ
セッサは超解像度情景24の評価を生成する。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に送られる。フレーム整
合エラープロセッサ25は、その組の個々のフレームの間
の残りの整合エラーを決定する。整合プロセッサ25は、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を回転開口帯域限定映像合成プロセッサ
26および超解像度評価プロセッサ23に送る。回転開口帯
域限定映像合成プロセッサ26は、米国特許第5,243,351
号明細書に記載された回転開口処理によって全円形開口
映像を表す映像を生成する。一組の補正されたフレー
ム、あるいは一組の補正されていないフレームおよびフ
レーム間の整合エラー評価、もしくはその両方、ならび
に回転開口帯域限定映像を使用して、超解像度評価プロ
セッサは超解像度情景24の評価を生成する。
【0023】図5を参照すると、第4のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および発生し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に転送する。さ
らに、正確に整合されたフレーム、あるいは一組の補正
されていないフレームおよびフレーム間の整合エラー評
価が回転開口帯域限定映像合成プロセッサ26に送られ
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、映像合成デー
タベース27に記憶された情景あるいは物体の構造の過去
の知識を組入れる。補正されたフレーム、あるいは一組
の補正されていないフレームおよび対応する一組の整合
エラー評価、ならびに回転開口帯域限定映像を使用し
て、超解像度評価プロセッサは超解像度情景評価24を生
成し、一方で、データベース27に記憶された情景あるい
は物体の構造の過去の知識を組入れる。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および発生し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に転送する。さ
らに、正確に整合されたフレーム、あるいは一組の補正
されていないフレームおよびフレーム間の整合エラー評
価が回転開口帯域限定映像合成プロセッサ26に送られ
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、映像合成デー
タベース27に記憶された情景あるいは物体の構造の過去
の知識を組入れる。補正されたフレーム、あるいは一組
の補正されていないフレームおよび対応する一組の整合
エラー評価、ならびに回転開口帯域限定映像を使用し
て、超解像度評価プロセッサは超解像度情景評価24を生
成し、一方で、データベース27に記憶された情景あるい
は物体の構造の過去の知識を組入れる。
【0024】図6を参照すると、第5のシステム10およ
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に送る。さら
に、正確に整合されたフレーム、あるいは一組の補正さ
れていないフレームおよびフレーム間の整合エラー評価
が回転開口帯域限定映像合成プロセッサ26に転送され
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、データベース
27に記憶された情景あるいは物体の構造の過去の知識を
組入れる。帯域限定された回転開口映像は、訓練アルゴ
リズムプロセッサ28に移送される。帯域限定された回転
開口映像から得られた情報に基づいて、訓練プロセッサ
は、情景の内容および構造に関する情報を抽出し、それ
はデータベース27に含まれた過去の情報を変更するため
に使用される。補正されたフレーム、あるいは補正され
ていないフレームおよび対応する一組の整合エラー評
価、ならびに回転開口帯域限定映像を使用して、超解像
度評価プロセッサ23は超解像度情景評価24を生成し、一
方で、データベース27に記憶され、訓練アルゴリズムプ
ロセッサ28によって変更された情景あるいは物体の構造
の過去の知識を組入れる。
び方法30が示されている。回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それらはフ
レーム整合エラープロセッサ25に転送される。フレーム
整合エラープロセッサ25は、その組の個々のフレームの
間の残りの整合エラーを決定する。整合プロセッサは、
一組の補正されたフレーム、あるいは一組の補正されて
いないフレームおよびフレーム間の整合エラー評価、も
しくはその両方を超解像度プロセッサ23に送る。さら
に、正確に整合されたフレーム、あるいは一組の補正さ
れていないフレームおよびフレーム間の整合エラー評価
が回転開口帯域限定映像合成プロセッサ26に転送され
る。回転開口帯域限定映像合成プロセッサ26は、米国特
許第5,243,351 号明細書に記載された手段によって全円
形開口映像を表す映像を生成し、一方で、データベース
27に記憶された情景あるいは物体の構造の過去の知識を
組入れる。帯域限定された回転開口映像は、訓練アルゴ
リズムプロセッサ28に移送される。帯域限定された回転
開口映像から得られた情報に基づいて、訓練プロセッサ
は、情景の内容および構造に関する情報を抽出し、それ
はデータベース27に含まれた過去の情報を変更するため
に使用される。補正されたフレーム、あるいは補正され
ていないフレームおよび対応する一組の整合エラー評
価、ならびに回転開口帯域限定映像を使用して、超解像
度評価プロセッサ23は超解像度情景評価24を生成し、一
方で、データベース27に記憶され、訓練アルゴリズムプ
ロセッサ28によって変更された情景あるいは物体の構造
の過去の知識を組入れる。
【0025】超スピンAp処理の実現方法が以下に説明
される。図7および8において、図2乃至6のシステム
10および方法30において使用される回転開口の超解像度
処理方法23の実施形態が示されている。図7および8に
示された方法23において、回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それはフレ
ーム整合エラープロセッサ25に送られる。フレーム整合
エラープロッサ25は、その一組の個々のフレーム間の残
りの整合エラーを決定する。整合プロセッサは、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に転送する。
される。図7および8において、図2乃至6のシステム
10および方法30において使用される回転開口の超解像度
処理方法23の実施形態が示されている。図7および8に
示された方法23において、回転開口センサ21は、一組の
ストリップ開口映像22を捕捉および生成し、それはフレ
ーム整合エラープロセッサ25に送られる。フレーム整合
エラープロッサ25は、その一組の個々のフレーム間の残
りの整合エラーを決定する。整合プロセッサは、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に転送する。
【0026】さらに、正確に整合されたフレームおよび
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に移される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた全円形開口映像を表す映像を生成する。帯域が限定
された回転開口映像は、超解像度情景の第1の評価33を
供給するために変更される。
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に移される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた全円形開口映像を表す映像を生成する。帯域が限定
された回転開口映像は、超解像度情景の第1の評価33を
供給するために変更される。
【0027】現在の情景評価34の観察から得られる各フ
レームの測定結果を決定するために現在の超解像度情景
評価33およびシステム応答関数37の過去の情報が使用さ
れる。予測された測定結果34は、現在の超解像度情景評
価33を個々のフレーム37のシステム応答関数でコンボリ
ューション処理することによって決定される。処理の流
れにおける次のステップ39および40では、測定結果22を
予測された測定結果34で除算し、そこから1を減算し、
その後、各フレームの反射された応答関数38でコンボリ
ューション処理する。結果として生じたコンボリューシ
ョン40は、各フレーム41毎に合計され、その後、ステッ
プ42においてスケーリング係数43で乗算される。第1の
超解像度評価における終わりの方のステップにおいて、
ステップ39乃至43の出力が指数関数化され(43)、指数
関数化された結果が現在の超解像度情景評価34で乗算さ
れる(ステップ45)。ステップ21乃至45の結果、回転開
口超解像度情景の第1の真の評価46が得られる。その
後、回転開口超解像度情景の第1の真の評価46が現在の
情景評価33として使用され、このプロセスは所望された
回数だけ反復される。
レームの測定結果を決定するために現在の超解像度情景
評価33およびシステム応答関数37の過去の情報が使用さ
れる。予測された測定結果34は、現在の超解像度情景評
価33を個々のフレーム37のシステム応答関数でコンボリ
ューション処理することによって決定される。処理の流
れにおける次のステップ39および40では、測定結果22を
予測された測定結果34で除算し、そこから1を減算し、
その後、各フレームの反射された応答関数38でコンボリ
ューション処理する。結果として生じたコンボリューシ
ョン40は、各フレーム41毎に合計され、その後、ステッ
プ42においてスケーリング係数43で乗算される。第1の
超解像度評価における終わりの方のステップにおいて、
ステップ39乃至43の出力が指数関数化され(43)、指数
関数化された結果が現在の超解像度情景評価34で乗算さ
れる(ステップ45)。ステップ21乃至45の結果、回転開
口超解像度情景の第1の真の評価46が得られる。その
後、回転開口超解像度情景の第1の真の評価46が現在の
情景評価33として使用され、このプロセスは所望された
回数だけ反復される。
【0028】図9および10において、図2乃至6のシ
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。この方
法23において、図9および10において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に移送される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に転送する。
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。この方
法23において、図9および10において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に移送される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に転送する。
【0029】また、正確に整合されたフレームおよびフ
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に送られる。回転開口帯域限定映像合成プロ
セッサ26は、米国特許第5,243,351 号明細書に記載され
た手段によって全円形開口映像を表す映像を生成する。
帯域限定された回転開口映像は、超解像度情景サンプル
手段の第1の評価33を供給するように変更される。
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に送られる。回転開口帯域限定映像合成プロ
セッサ26は、米国特許第5,243,351 号明細書に記載され
た手段によって全円形開口映像を表す映像を生成する。
帯域限定された回転開口映像は、超解像度情景サンプル
手段の第1の評価33を供給するように変更される。
【0030】現在の情景サンプル平均評価34の観察によ
り得られる各フレームの測定を決定するために現在の超
解像度情景サンプル平均評価33およびシステム応答関数
37の過去の情報が使用される。予測された測定結果34
は、現在の超解像度情景サンプル平均評価33を個々のフ
レーム37のシステム応答関数でコンボリューション処理
することによって決定される。処理の流れにおける次の
ステップ39および40では、測定結果22を予測された測定
結果34で除算し、そこから1を減算し、その後、各フレ
ームの反射された応答関数38でコンボリューションする
(ステップ40)。結果として生じたコンボリューション
40は、ステップ41において各フレーム毎に合計され、そ
の後、ステップ42においてスケーリング係数43で乗算さ
れる。第1の超解像度評価における最終のステップにお
いて、ステップ39乃至43の出力が指数関数化(ステップ
44)され、ステップ45において、指数関数化された結果
が現在の超解像度情景サンプル平均評価33と乗算され
る。ステップ21乃至46の結果、回転開口超解像度情景の
第1の真の評価46が得られる。その後、回転開口超解像
度情景の第1の真の評価46が現在の情景サンプル平均評
価33を更新するために使用され、このプロセスは所望さ
れた回数だけ反復される。
り得られる各フレームの測定を決定するために現在の超
解像度情景サンプル平均評価33およびシステム応答関数
37の過去の情報が使用される。予測された測定結果34
は、現在の超解像度情景サンプル平均評価33を個々のフ
レーム37のシステム応答関数でコンボリューション処理
することによって決定される。処理の流れにおける次の
ステップ39および40では、測定結果22を予測された測定
結果34で除算し、そこから1を減算し、その後、各フレ
ームの反射された応答関数38でコンボリューションする
(ステップ40)。結果として生じたコンボリューション
40は、ステップ41において各フレーム毎に合計され、そ
の後、ステップ42においてスケーリング係数43で乗算さ
れる。第1の超解像度評価における最終のステップにお
いて、ステップ39乃至43の出力が指数関数化(ステップ
44)され、ステップ45において、指数関数化された結果
が現在の超解像度情景サンプル平均評価33と乗算され
る。ステップ21乃至46の結果、回転開口超解像度情景の
第1の真の評価46が得られる。その後、回転開口超解像
度情景の第1の真の評価46が現在の情景サンプル平均評
価33を更新するために使用され、このプロセスは所望さ
れた回数だけ反復される。
【0031】図11および12において、図2乃至6の
システム10および方法30において使用される回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図11および12において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に送られる。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。整合されたフレーム29から、超ス
ピンApプロセッサはステップ29a において全ての検出
器およびフレームの測定結果の合計を決定し、全ての検
出器の測定結果と全てのフレームの測定結果の合計は1
より小さい。その合計は、正規化された補正項39a の成
分となる。
システム10および方法30において使用される回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図11および12において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に送られる。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。整合されたフレーム29から、超ス
ピンApプロセッサはステップ29a において全ての検出
器およびフレームの測定結果の合計を決定し、全ての検
出器の測定結果と全てのフレームの測定結果の合計は1
より小さい。その合計は、正規化された補正項39a の成
分となる。
【0032】また、正確に整合されたフレームおよびフ
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に供給される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた手段によって全円形開口映像を表す映像を生成す
る。帯域限定された回転開口映像は、超解像度情景の第
1の評価33を供給するように変更される。
レーム間のエラー評価は、回転開口帯域限定映像合成プ
ロセッサ26に供給される。回転開口帯域限定映像合成プ
ロセッサ26は、米国特許第5,243,351 号明細書に記載さ
れた手段によって全円形開口映像を表す映像を生成す
る。帯域限定された回転開口映像は、超解像度情景の第
1の評価33を供給するように変更される。
【0033】現在の情景評価34の観察から得られる各フ
レームの測定を決定するために現在の超解像度情景評価
33およびシステム応答関数37の過去の情報が使用され
る。予測された測定結果34は、現在の超解像度情景評価
33を個々のフレーム37のシステム応答関数37でコンボリ
ューション処理することによって決定される。予測され
た測定結果34はまた、現在の情景評価33の観察から得ら
れる全てのフレームおよび画素の合計35を決定するため
に使用される。予測された測定結果34は、ステップ36に
おいて1より小さい全てのフレームにわたる全ての画素
および全ての画素にわたって合計され、それによって、
正規化された補正項39a の成分が生成される。
レームの測定を決定するために現在の超解像度情景評価
33およびシステム応答関数37の過去の情報が使用され
る。予測された測定結果34は、現在の超解像度情景評価
33を個々のフレーム37のシステム応答関数37でコンボリ
ューション処理することによって決定される。予測され
た測定結果34はまた、現在の情景評価33の観察から得ら
れる全てのフレームおよび画素の合計35を決定するため
に使用される。予測された測定結果34は、ステップ36に
おいて1より小さい全てのフレームにわたる全ての画素
および全ての画素にわたって合計され、それによって、
正規化された補正項39a の成分が生成される。
【0034】処理の流れにおける次のステップ40では、
測定結果22を予測された測定結果34で除算し、1を減算
した後、それは各フレームの反射された応答関数38でコ
ンボリューション処理される。結果として生じたコンボ
リューション40は、各フレーム41に対して合計され(ス
テップ41)、その後、ステップ29a,35,36 の結果を有す
る正規化された補正項がこの合計から減算される。この
結果は、ステップ42においてスケーリング係数43で乗算
される。第1の超解像度評価における最後の方のステッ
プ44において、ステップ39乃至43の出力が指数関数化さ
れ、ステップ45において、指数関数化された結果が現在
の超解像度情景評価33と乗算される。ステップ21乃至45
の結果が、回転開口超解像度情景の第1の真の評価46と
なる。その後、回転開口超解像度情景の第1の真の評価
46が現在の情景評価33として使用され、このプロセスは
所望された回数だけ反復される。
測定結果22を予測された測定結果34で除算し、1を減算
した後、それは各フレームの反射された応答関数38でコ
ンボリューション処理される。結果として生じたコンボ
リューション40は、各フレーム41に対して合計され(ス
テップ41)、その後、ステップ29a,35,36 の結果を有す
る正規化された補正項がこの合計から減算される。この
結果は、ステップ42においてスケーリング係数43で乗算
される。第1の超解像度評価における最後の方のステッ
プ44において、ステップ39乃至43の出力が指数関数化さ
れ、ステップ45において、指数関数化された結果が現在
の超解像度情景評価33と乗算される。ステップ21乃至45
の結果が、回転開口超解像度情景の第1の真の評価46と
なる。その後、回転開口超解像度情景の第1の真の評価
46が現在の情景評価33として使用され、このプロセスは
所望された回数だけ反復される。
【0035】図13および14において、図2乃至6の
システム10および方法30において使用された回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図13および14において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。フレーム整合エラープロセッサ25は、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に供給する。
システム10および方法30において使用された回転開口超
解像度処理方法23の別の実施形態が示されている。方法
23において、図13および14において示されているよ
うに、回転開口センサ21は一組のストリップ開口映像22
を捕捉および生成し、それはフレーム整合エラープロセ
ッサ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。フレーム整合エラープロセッサ25は、一組の
補正されたフレーム29およびフレーム間の整合エラー評
価を超解像度プロセッサ23に供給する。
【0036】さらに、正確に整合されたフレームおよび
フレーム間の整合エラー評価は、回転開口帯域限定合成
映像プロセッサ26に供給される。回転開口帯域限定映像
合成プロセッサ26は、米国特許第5,243,351 号明細書に
記載された全円形開口映像を表す映像を生成する。帯域
が限定された回転開口映像は、超解像度情景アンサンブ
ル平均の第1の評価33を供給するために変更される。
フレーム間の整合エラー評価は、回転開口帯域限定合成
映像プロセッサ26に供給される。回転開口帯域限定映像
合成プロセッサ26は、米国特許第5,243,351 号明細書に
記載された全円形開口映像を表す映像を生成する。帯域
が限定された回転開口映像は、超解像度情景アンサンブ
ル平均の第1の評価33を供給するために変更される。
【0037】現在の情景アンサンブル平均評価34の観察
から得られる各フレームの測定を決定するように現在の
超解像度情景アンサンブル平均評価33およびシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
個々のフレームのシステム応答関数37でコンボリューシ
ョンすることによって決定される。処理の流れにおける
次のステップ39では、測定結果22を予測された測定結果
34で除算し、その後、ステップ40において各フレームの
反射された応答関数38でコンボリューション処理する。
結果として生じたコンボリューション40は、ステップ41
において各フレームに対して合計され、その後、ステッ
プ42においてスケーリング係数43と乗算される。第1の
超解像度評価における終段のステップ45において、ステ
ップ39乃至43の出力が現在の超解像度情景アンサンブル
平均評価34と乗算される。ステップ21乃至45の結果が、
回転開口超解像度情景アンサンブル平均の第1の真の評
価46となる。その後、回転開口超解像度情景の第1の真
の評価46が現在の情景アンサンブル平均評価33として使
用され、このプロセスは所望された回数だけ反復され
る。
から得られる各フレームの測定を決定するように現在の
超解像度情景アンサンブル平均評価33およびシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
個々のフレームのシステム応答関数37でコンボリューシ
ョンすることによって決定される。処理の流れにおける
次のステップ39では、測定結果22を予測された測定結果
34で除算し、その後、ステップ40において各フレームの
反射された応答関数38でコンボリューション処理する。
結果として生じたコンボリューション40は、ステップ41
において各フレームに対して合計され、その後、ステッ
プ42においてスケーリング係数43と乗算される。第1の
超解像度評価における終段のステップ45において、ステ
ップ39乃至43の出力が現在の超解像度情景アンサンブル
平均評価34と乗算される。ステップ21乃至45の結果が、
回転開口超解像度情景アンサンブル平均の第1の真の評
価46となる。その後、回転開口超解像度情景の第1の真
の評価46が現在の情景アンサンブル平均評価33として使
用され、このプロセスは所望された回数だけ反復され
る。
【0038】図15および16おいて、図2乃至6のシ
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。方法23
において、図15および16において示されているよう
に、回転開口センサ21は一組のストリップ開口映像22を
捕捉および生成し、それはフレーム整合エラープロセッ
サ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。
ステム10および方法30において使用される回転開口超解
像度処理方法23の別の実施形態が示されている。方法23
において、図15および16において示されているよう
に、回転開口センサ21は一組のストリップ開口映像22を
捕捉および生成し、それはフレーム整合エラープロセッ
サ25に供給される。フレーム整合エラープロセッサ25
は、その一組の個々のフレーム間の残りの整合エラーを
決定する。整合プロセッサは、一組の補正されたフレー
ム29およびフレーム間の整合エラー評価を超解像度プロ
セッサ23に供給する。
【0039】さらに、正確に整合されたフレームおよび
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に転送される。回転開口帯域限定映像合成
プロセッサ26は、米国特許第5,243,351 号明細書に記載
されたように全円形開口映像を表す映像を生成する。帯
域が限定された回転開口映像は、超解像度情景アンサン
ブル平均の第1の評価33を供給するために変更される。
フレーム間のエラー評価は、回転開口帯域限定映像合成
プロセッサ26に転送される。回転開口帯域限定映像合成
プロセッサ26は、米国特許第5,243,351 号明細書に記載
されたように全円形開口映像を表す映像を生成する。帯
域が限定された回転開口映像は、超解像度情景アンサン
ブル平均の第1の評価33を供給するために変更される。
【0040】現在の情景アンサンブル平均評価34の観察
から得られる第1のフレームに対する測定結果を決定す
るように現在の超解像度情景アンサンブル平均評価33お
よび第1の回転開口フレーム測定に関連したシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
第1のフレーム37のシステム応答関数でコンボリューシ
ョン処理することによって決定される。処理の流れにお
ける次のステップ39では、第1の回転開口の測定結果22
を予測された第1の測定結果34で除算し、その後、それ
をステップ40において第1のフレームの反射された応答
関数38でコンボリューション処理する。ステップ39b に
おいて、結果として生じたコンボリューションから1を
減算し、その結果は、ステップ41a において情景アンサ
ンブル平均の現在の評価34で乗算される。この結果は、
ステップ42においてスケーリング係数43と乗算される。
第1の超解像度評価における終段のステップ44a におい
て、ステップ39乃至43の出力が現在の超解像度情景アン
サンブル平均評価34に加算される。ステップ21乃至45の
結果が、回転開口超解像度情景アンサンブル平均の第1
の真の評価46となる。その後、回転開口超解像度情景の
第1の評価46が現在の情景アンサンブル平均評価33とし
て使用され、このプロセスは所望された回数だけ反復さ
れる。
から得られる第1のフレームに対する測定結果を決定す
るように現在の超解像度情景アンサンブル平均評価33お
よび第1の回転開口フレーム測定に関連したシステム応
答関数37の過去の情報が使用される。予測された測定結
果34は、現在の超解像度情景アンサンブル平均評価33を
第1のフレーム37のシステム応答関数でコンボリューシ
ョン処理することによって決定される。処理の流れにお
ける次のステップ39では、第1の回転開口の測定結果22
を予測された第1の測定結果34で除算し、その後、それ
をステップ40において第1のフレームの反射された応答
関数38でコンボリューション処理する。ステップ39b に
おいて、結果として生じたコンボリューションから1を
減算し、その結果は、ステップ41a において情景アンサ
ンブル平均の現在の評価34で乗算される。この結果は、
ステップ42においてスケーリング係数43と乗算される。
第1の超解像度評価における終段のステップ44a におい
て、ステップ39乃至43の出力が現在の超解像度情景アン
サンブル平均評価34に加算される。ステップ21乃至45の
結果が、回転開口超解像度情景アンサンブル平均の第1
の真の評価46となる。その後、回転開口超解像度情景の
第1の評価46が現在の情景アンサンブル平均評価33とし
て使用され、このプロセスは所望された回数だけ反復さ
れる。
【0041】超解像度の回転開口の式を以下に示す。S
m により示される時間的に静止した情景を実現するた
め、回転開口センサの焦点平面アレイからの雑音のない
検出器出力Gm は次式のように表されることができる。
m により示される時間的に静止した情景を実現するた
め、回転開口センサの焦点平面アレイからの雑音のない
検出器出力Gm は次式のように表されることができる。
【0042】
【数1】下付け文字mはストリップ開口のm番目の回転
方向を示しており、その期間にm番目の画像フレームが
捕捉され、Hm は全体のシステム応答関数である。Sm
はm番目の測定に関する情景の実体を表している。これ
は雑音のない回転開口画像捕捉の最も一般的なモデルで
ある。
方向を示しており、その期間にm番目の画像フレームが
捕捉され、Hm は全体のシステム応答関数である。Sm
はm番目の測定に関する情景の実体を表している。これ
は雑音のない回転開口画像捕捉の最も一般的なモデルで
ある。
【0043】システム応答関数による情景のコンボリュ
ーションのディスクリートな形態は次式により表される
ことができる。
ーションのディスクリートな形態は次式により表される
ことができる。
【0044】
【数2】ここでjは位置/rjからの検出器出力を示し
ており、kは焦点平面の空間的位置/rkへ幾何学的に
投影されたディスクリートな情景素子を示している。
ており、kは焦点平面の空間的位置/rkへ幾何学的に
投影されたディスクリートな情景素子を示している。
【0045】システム応答関数Hm はm番目の画像捕捉
用の光学システム点分布関数PSFm のコンボリューシ
ョンであり、個々の検出器応答関数Wm はm番目のフレ
ームの光学システム点分布関数PSFm のコンボリュー
ションであり、次式のように表される。
用の光学システム点分布関数PSFm のコンボリューシ
ョンであり、個々の検出器応答関数Wm はm番目のフレ
ームの光学システム点分布関数PSFm のコンボリュー
ションであり、次式のように表される。
【0046】
【数3】Gm はGm により示される列マトリックスによ
り表されることができ、各エレメントは焦点平面アレイ
の各検出器測定に対応する。Hm はHm により示されて
いる2次元マトリックスにより表され、エレメントは電
子光センサのシステム応答関数に対応する。同様に、情
景はSm により示される列マトックスとして表され、超
解像度のグリッドへ投影されるm番目の情景実体のサン
プルに対応する素子Sx を有する。
り表されることができ、各エレメントは焦点平面アレイ
の各検出器測定に対応する。Hm はHm により示されて
いる2次元マトリックスにより表され、エレメントは電
子光センサのシステム応答関数に対応する。同様に、情
景はSm により示される列マトックスとして表され、超
解像度のグリッドへ投影されるm番目の情景実体のサン
プルに対応する素子Sx を有する。
【0047】超スピン開口プロセス処理の1実施形態
は、最大事後確率評価と通常呼ばれているバイエル(Ba
yes )評価方法を使用する。最大事後評価は測定を行う
目標物体の条件付き確率を最大にする未知の目標物体の
値を求める。超スピン開口に対する最大事後方法は1組
の回転開口センサ測定を行う超解像度の情景の確率を最
大にする超解像度の情景の評価を模索する。超スピン開
口のこの実施形態では、式2で与えられる画像モデルは
次式により特定化される。
は、最大事後確率評価と通常呼ばれているバイエル(Ba
yes )評価方法を使用する。最大事後評価は測定を行う
目標物体の条件付き確率を最大にする未知の目標物体の
値を求める。超スピン開口に対する最大事後方法は1組
の回転開口センサ測定を行う超解像度の情景の確率を最
大にする超解像度の情景の評価を模索する。超スピン開
口のこの実施形態では、式2で与えられる画像モデルは
次式により特定化される。
【0048】
【数4】ここで情景Sは回転開口測定間隔期間の単一の
情景実体であると仮定されている。
情景実体であると仮定されている。
【0049】バイエルの法則を回転開口測定の全ての組
に適用すると、Sの測定値を与える条件付き確率の式が
導かれる。
に適用すると、Sの測定値を与える条件付き確率の式が
導かれる。
【0050】
【数5】Gは全測定マトリックスであり、これは個々の
フレーム測定マトリックスからなるブロック列のマトリ
ックスである。
フレーム測定マトリックスからなるブロック列のマトリ
ックスである。
【0051】
【数6】ここでnf は捕捉されるフレームの総数であ
る。Sは超解像度グリッドに投影される情景サンプルの
マトリックスである。
る。Sは超解像度グリッドに投影される情景サンプルの
マトリックスである。
【0052】情景放射または各画素素子に対する光子の
反射が統計的に独立していると仮定すると、1次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
反射が統計的に独立していると仮定すると、1次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
【0053】
【数7】ソース光子変動により左右される情景統計で
は、情景画素素子の確率密度関数はポアッソン(Poisso
n )の分布式であると仮定され、
は、情景画素素子の確率密度関数はポアッソン(Poisso
n )の分布式であると仮定され、
【0054】
【数8】ここでkは情景の投影された画素素子を示す指
数である。
数である。
【0055】統計的に独立したフレーム測定では、情景
が与えられた全測定ベクトルのジョイントの確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
が与えられた全測定ベクトルのジョイントの確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
【0056】
【数9】同様に、統計的に独立した検出器測定では、情
景が与えられたフレーム測定の発生についての条件付き
確率は情景上で条件付けされた個々の画素確率密度の積
として表されることができる。
景が与えられたフレーム測定の発生についての条件付き
確率は情景上で条件付けされた個々の画素確率密度の積
として表されることができる。
【0057】
【数10】全体の測定の条件付き確率密度はそれ故、情
景上で条件付けされた画素およびフレーム測定確率の積
により与えられる。
景上で条件付けされた画素およびフレーム測定確率の積
により与えられる。
【0058】
【数11】個々のベクトル測定がポアッソンプロセスで
あると仮定すると、個々の測定の条件付き確率は次式の
ように与えられることができる。
あると仮定すると、個々の測定の条件付き確率は次式の
ように与えられることができる。
【0059】
【数12】ここでΣk Hm (j−k)S(k)は個々の
検出器測定の平均値である。条件付き確率密度を最大に
する情景を発見することは確率密度の対数を最大にする
情景を発見することに等しいので、次式を最大にするS
に対する解が求められる。
検出器測定の平均値である。条件付き確率密度を最大に
する情景を発見することは確率密度の対数を最大にする
情景を発見することに等しいので、次式を最大にするS
に対する解が求められる。
【0060】
【数13】積の対数が対数の和であることを認識する
と、Sが与えられた測定マトリックスGの確率は次式の
ようになる。
と、Sが与えられた測定マトリックスGの確率は次式の
ようになる。
【0061】
【数14】
【0062】
【数15】
【0063】(13)を(12)に代入すると、以下の
スターリング(Stirling)の概算式 ln[x!]=xln[x]−x (14) をS(k)!に適用すると次式が得られる。
スターリング(Stirling)の概算式 ln[x!]=xln[x]−x (14) をS(k)!に適用すると次式が得られる。
【0064】
【数16】条件付き確率を最大にするS(l)値を決定
するため、ln[p(S|G)]の導関数がゼロになる
値S(l)が必要とされる。ln[p(S|G)]の導
関数は次式により与えられる。
するため、ln[p(S|G)]の導関数がゼロになる
値S(l)が必要とされる。ln[p(S|G)]の導
関数は次式により与えられる。
【0065】
【数17】導関数をゼロに設定すると、Sを解く非線形
の式が得られる。
の式が得られる。
【0066】
【数18】この式はHm を指数関数化し係数化した後、
次式のようになる。
次式のようになる。
【0067】
【数19】ピカール(Picard)法として知られる反復プ
ロセスは非線形の式の最も可能性のある未知の情景の解
を解くために使用される。反復プロセスは以下の式によ
り示される。
ロセスは非線形の式の最も可能性のある未知の情景の解
を解くために使用される。反復プロセスは以下の式によ
り示される。
【0068】
【数20】ここで上付き文字[n]と[n+1]はそれ
ぞれ現在および次の情景評価を表している。
ぞれ現在および次の情景評価を表している。
【0069】超解像度の回転開口が元来の測定よりも高
い空間周波数に対応するグリッド点を評価するので、反
復プロセスは、超解像度の座標グリッドと、測定(また
は中間基準)グリッドとの間で再度マッピングすること
を含んでいる。反復的アルゴリズムの1実施形態はブロ
ック平均により超解像度グリッドのダウンサンプリング
を行い、複製により測定グリッドのアップサンプリング
を行うことによって、測定座標グリッドと超解像度の座
標グリッドとの間のサンプル間隔の差に適応する。代わ
りに、アップサンプリングとダウンサンプリングはより
高い次数の補間技術を使用して実行されることができ
る。反復的解は次式により与えられる。
い空間周波数に対応するグリッド点を評価するので、反
復プロセスは、超解像度の座標グリッドと、測定(また
は中間基準)グリッドとの間で再度マッピングすること
を含んでいる。反復的アルゴリズムの1実施形態はブロ
ック平均により超解像度グリッドのダウンサンプリング
を行い、複製により測定グリッドのアップサンプリング
を行うことによって、測定座標グリッドと超解像度の座
標グリッドとの間のサンプル間隔の差に適応する。代わ
りに、アップサンプリングとダウンサンプリングはより
高い次数の補間技術を使用して実行されることができ
る。反復的解は次式により与えられる。
【0070】
【数21】上向きの矢印でアップサンプリングを、下向
きの矢印でダウンサンプリングを示す。ここで指数関数
の和は次式を使用して簡単にされる。
きの矢印でダウンサンプリングを示す。ここで指数関数
の和は次式を使用して簡単にされる。
【0071】
【数22】現在の評価値/S(l)[n] の平均値は現在
の評価値S(l)[n] に置き換えられる。スケール係数
はコンバージェンスを制御するため指数関数に含まれて
もよい。
の評価値S(l)[n] に置き換えられる。スケール係数
はコンバージェンスを制御するため指数関数に含まれて
もよい。
【0072】先に説明した超スピン開口の実施形態は情
景Sを評価するため最大事後評価方法を使用する。別の
実施形態は情景を表すためにサンプリング統計を使用す
る。式(1)と(2)のSm とGm はランダムフィール
ド(ランダム変数の二次元アレイ)であるので、これら
のランダムフィールドの平均値を解析することが可能で
ある。測定全体にわたるサンプル平均を示すため< >
を使用すると、コンボリューションの平均値は次式のよ
うに表されることができる。
景Sを評価するため最大事後評価方法を使用する。別の
実施形態は情景を表すためにサンプリング統計を使用す
る。式(1)と(2)のSm とGm はランダムフィール
ド(ランダム変数の二次元アレイ)であるので、これら
のランダムフィールドの平均値を解析することが可能で
ある。測定全体にわたるサンプル平均を示すため< >
を使用すると、コンボリューションの平均値は次式のよ
うに表されることができる。
【0073】
【数23】ここで情景実体全体のサンプル平均は以下の
ように表される。
ように表される。
【0074】
【数24】ここでnf は測定された回転開口のフレーム
の数である。先の実施形態のように反復的評価プロセス
を得るために最大事後方程式は公式化され、ピカール技
術が適用される。
の数である。先の実施形態のように反復的評価プロセス
を得るために最大事後方程式は公式化され、ピカール技
術が適用される。
【0075】この実施形態ではGは式6で定められ、S
は同様に情景実体で限定される。
は同様に情景実体で限定される。
【0076】情景放射または各画素素子に対する光子の
反射が統計的に独立し、実体Sm が測定フレームにわた
って統計的に独立しているものと仮定すると、1次元情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができる。即ち、
反射が統計的に独立し、実体Sm が測定フレームにわた
って統計的に独立しているものと仮定すると、1次元情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができる。即ち、
【0077】
【数25】ソース光子の変動により左右される情景統計
では、情景画素素子の確率密度関数は次式のポアッソン
(Poisson )の分布式であると仮定され、
では、情景画素素子の確率密度関数は次式のポアッソン
(Poisson )の分布式であると仮定され、
【0078】
【数26】kは情景の投影された画素素子を示す指数で
ある。全体の平均値/S(k)はサンプル平均式3aに
より近似され、次式が得られる。
ある。全体の平均値/S(k)はサンプル平均式3aに
より近似され、次式が得られる。
【0079】
【数27】ダミー指数pは特定のフレームmからのフレ
ームにわたるカウントを弁別するために使用されてい
る。
ームにわたるカウントを弁別するために使用されてい
る。
【0080】統計と独立したフレーム測定では、情景を
与えられた全体の測定ベクトルのジョイント確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
与えられた全体の測定ベクトルのジョイント確率密度は
個々のm番目のフレームジョイント確率密度の積として
表されることができる。
【0081】
【数28】同様に、統計と独立したベクトル測定では、
情景を与えられたフレーム測定の発生の条件付き確率は
情景について条件付けされた個々の画素確率密度の積と
して表されることができる。
情景を与えられたフレーム測定の発生の条件付き確率は
情景について条件付けされた個々の画素確率密度の積と
して表されることができる。
【0082】
【数29】それ故、全体の測定の条件付き確率密度は情
景について条件付けされた画素およびフレーム測定確率
の積により与えられる。
景について条件付けされた画素およびフレーム測定確率
の積により与えられる。
【0083】
【数30】個々の検出器の測定がポアッソンプロセスで
あると仮定すると、個々の測定の条件付き確率は次式に
より表されることができる。
あると仮定すると、個々の測定の条件付き確率は次式に
より表されることができる。
【0084】
【数31】ここで/G(j)は個々の検出器測定の全体
的平均値である。全体的平均値は式3のサンプル平均値
<Gm (j)>により評価され、以下の測定確率密度が
導かれる。
的平均値である。全体的平均値は式3のサンプル平均値
<Gm (j)>により評価され、以下の測定確率密度が
導かれる。
【0085】
【数32】条件付き確率密度を最大にする情景を発見す
ることは確率密度の対数を最大にする情景を発見するこ
とに等しいので、次式を最大にするSの解が求められ
る。
ることは確率密度の対数を最大にする情景を発見するこ
とに等しいので、次式を最大にするSの解が求められ
る。
【0086】
【数33】
【0087】積の対数が対数の和であることを認識する
と、Sが与えられる測定マトリックスGの確率は次式の
ように表されることができる。
と、Sが与えられる測定マトリックスGの確率は次式の
ように表されることができる。
【0088】
【数34】
【0089】式(25)を式(24)に代入すると、以
下の式が導かれる。
下の式が導かれる。
【0090】
【数35】
【0091】以下のスターリングの概算式 ln[x!]=xln[x]−x (26) をSm (k)!に適用すると次式が得られる。
【0092】
【数36】
【0093】反復的な手順が得られ、ここでq番目の実
体Sq が評価され、この評価値は<S>の評価値を更新
するために使用される。条件付き確率を最大にする値S
q (l)を決定するため、ln[p(S|G)]の導関
数がゼロになる値Sq (l)が必要とされる。ln[p
(S|G)]の導関数は次式により与えられる。
体Sq が評価され、この評価値は<S>の評価値を更新
するために使用される。条件付き確率を最大にする値S
q (l)を決定するため、ln[p(S|G)]の導関
数がゼロになる値Sq (l)が必要とされる。ln[p
(S|G)]の導関数は次式により与えられる。
【0094】
【数37】
【0095】導関数をゼロに設定すると、Sq を解く次
の非線形の式が得られる。
の非線形の式が得られる。
【0096】
【数38】この式はHm を指数関数化し係数化した後、
次式のようになる。
次式のようになる。
【0097】
【数39】
【0098】ピカール法として知られる反復プロセスは
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用されることができる。反復プロセスは以下の式
により示される。
非線形の式の最も可能性のある未知の情景の解を解くた
めに使用されることができる。反復プロセスは以下の式
により示される。
【0099】
【数40】上向きの矢印はアップサンプリングを、下向
きの矢印はダウンサンプリングを示す。ここで上付き文
字[n]と[n+1]はそれぞれ現在および次の情景評
価を表している。先の実施形態のように、上下の矢印で
示されるアップサンプリングとダウンサンプリングは測
定座標グリッドと超解像度の座標グリッドとの間のサン
プル間隔の差に適応するために使用されることができ
る。
きの矢印はダウンサンプリングを示す。ここで上付き文
字[n]と[n+1]はそれぞれ現在および次の情景評
価を表している。先の実施形態のように、上下の矢印で
示されるアップサンプリングとダウンサンプリングは測
定座標グリッドと超解像度の座標グリッドとの間のサン
プル間隔の差に適応するために使用されることができ
る。
【0100】全体の測定に適応する単一の確率密度を使
用する最大事後方法を以下説明する。先に示したアルゴ
リズム手順は基本的に決定しなければならない2つの未
知数と、超解像度の情景画素と超解像度の情景画素平均
値を有する。第1の組の超解像度情景評価に関するアル
ゴリズムの実行は情景画素平均値を評価するための付属
プロセスを伴う。第2の組の超解像度情景評価式に関す
るアルゴリズムの実行は各情景画素の全体平均を概算す
るためにサンプル統計を使用する。
用する最大事後方法を以下説明する。先に示したアルゴ
リズム手順は基本的に決定しなければならない2つの未
知数と、超解像度の情景画素と超解像度の情景画素平均
値を有する。第1の組の超解像度情景評価に関するアル
ゴリズムの実行は情景画素平均値を評価するための付属
プロセスを伴う。第2の組の超解像度情景評価式に関す
るアルゴリズムの実行は各情景画素の全体平均を概算す
るためにサンプル統計を使用する。
【0101】代わりの方法は定数である独立したランダ
ム変数の総数を認識し、即ち全ての合計された測定の総
数または合計された部分的測定の総数は個々の測定で条
件付けされたランダム変数である。
ム変数の総数を認識し、即ち全ての合計された測定の総
数または合計された部分的測定の総数は個々の測定で条
件付けされたランダム変数である。
【0102】それ故、合計された測定の総数と1より小
さい1組の測定が与えられた超解像度の情景の条件付き
確率を最大にするバイエル(Bayes )法則が求められ
る。代わりのプロセスの体系はここで説明した第1の実
施形態と同一である。
さい1組の測定が与えられた超解像度の情景の条件付き
確率を最大にするバイエル(Bayes )法則が求められ
る。代わりのプロセスの体系はここで説明した第1の実
施形態と同一である。
【0103】光子変動により左右される情景および測定
統計では、個々の情景および測定素子の確率密度関数は
ポアッソンであると仮定され、即ち次式のようになる。
統計では、個々の情景および測定素子の確率密度関数は
ポアッソンであると仮定され、即ち次式のようになる。
【0104】
【数41】ここで/S(k↑)は情景のk↑番目の統計
された“画素”の平均値であり、Gm (k↓)はフレー
ムmからのj番目の検出器測定である。
された“画素”の平均値であり、Gm (k↓)はフレー
ムmからのj番目の検出器測定である。
【0105】平均値/X1 と個々のポアッソン分布を有
するn個の統計的に独立したポアッソンランダム変数X
1 …Xn では以下のようになる。
するn個の統計的に独立したポアッソンランダム変数X
1 …Xn では以下のようになる。
【0106】
【数42】ポアッソンランダム変数の合計は次式のよう
になり、
になり、
【0107】
【数43】この式はF. A. Haight氏による文献Handbook
of the Poisson Distributionの71頁で与えられる分布
を有するポアッソンランダム変数である。
of the Poisson Distributionの71頁で与えられる分布
を有するポアッソンランダム変数である。
【0108】
【数44】ここで、
【0109】
【数45】
【0110】平均値/X1 …/Xn を有するn個の統計
的に独立したポアッソンランダム変数X1 …Xn では、
個々のポアッソン分布p(Xk )と、X1 ,…Xn-2 ,
Xn-1 の条件付き確率分布と、以下の合計されたランダ
ム変数の総数、即ち、
的に独立したポアッソンランダム変数X1 …Xn では、
個々のポアッソン分布p(Xk )と、X1 ,…Xn-2 ,
Xn-1 の条件付き確率分布と、以下の合計されたランダ
ム変数の総数、即ち、
【0111】
【数46】は多項式であり、F. A. Haight氏による文献
Handbook of the Poisson Distributionの71頁により与
えられている。
Handbook of the Poisson Distributionの71頁により与
えられている。
【0112】
【数47】
【0113】式(28)を1組の超スピン開口測定に適
用することは、条件付き確率に明白に存在しない測定を
選択することに依存する。測定のナンバリングは任意で
あるので、合計された測定の総数が与えられる条件付き
確率は次式のように表される
用することは、条件付き確率に明白に存在しない測定を
選択することに依存する。測定のナンバリングは任意で
あるので、合計された測定の総数が与えられる条件付き
確率は次式のように表される
【0114】。
【数48】ここで、Nが単位フレーム当たりの検出器数
−1であるならば、M=フレーム数であり、Nが単位フ
レーム当たりの検出器数であるならば、M=フレーム数
−1であり、次式が成り立つ。
−1であるならば、M=フレーム数であり、Nが単位フ
レーム当たりの検出器数であるならば、M=フレーム数
−1であり、次式が成り立つ。
【0115】
【数49】
【0116】データの任意のサブセットに関しては、次
式が成り立つ。
式が成り立つ。
【0117】
【数50】独立したランダム変数の総数が一定であり、
即ち合計された全ての測定の総数または合計された測定
部分の総数が個々の測定で条件付けされたランダム変数
であることを認識することは、超解像度情景の評価を発
見するため別の変形された自己一貫最大事後方法を設定
するための手段を与える。全体の測定と1より小さい1
組の測定が与えられる超解像度の情景のバイエルの評価
規則は次式のようになる。
即ち合計された全ての測定の総数または合計された測定
部分の総数が個々の測定で条件付けされたランダム変数
であることを認識することは、超解像度情景の評価を発
見するため別の変形された自己一貫最大事後方法を設定
するための手段を与える。全体の測定と1より小さい1
組の測定が与えられる超解像度の情景のバイエルの評価
規則は次式のようになる。
【0118】
【数51】ここでG_は測定の1より小さい総セットに
対応するブロックマトリックスであり、Sは超解像度グ
リッドに投影される情景サンプルのマトリックスであ
る。マトリックスG_の素子は以下の組に対応する。
対応するブロックマトリックスであり、Sは超解像度グ
リッドに投影される情景サンプルのマトリックスであ
る。マトリックスG_の素子は以下の組に対応する。
【0119】
【数52】ここで、N=単位フレーム当たりの検出器数
−1であるならば、Mはフレーム数であり、またはN=
単位フレーム当たりの検出器数であるならば、Mはフレ
ーム数−1である。
−1であるならば、Mはフレーム数であり、またはN=
単位フレーム当たりの検出器数であるならば、Mはフレ
ーム数−1である。
【0120】自己一貫の必要性は、個々の検出器測定の
平均値を仮定することにより再度取入れられ、全体の回
転開口測定の和の平均値は雑音のない組の測定に対応
し、即ち次式のようになる。
平均値を仮定することにより再度取入れられ、全体の回
転開口測定の和の平均値は雑音のない組の測定に対応
し、即ち次式のようになる。
【0121】
【数53】
【0122】↑と↓は超解像度および測定座標の空間グ
リッドを示している。
リッドを示している。
【0123】画像センサ測定が情景から伝播される光子
の検出に対応するので、情景の評価値は正でなければな
らない。それ故、付加的な自己一貫の要求は、 S≧0 (31)であり、これは超解像度アルゴ
リズムの実施形態により不明瞭に満足される。
の検出に対応するので、情景の評価値は正でなければな
らない。それ故、付加的な自己一貫の要求は、 S≧0 (31)であり、これは超解像度アルゴ
リズムの実施形態により不明瞭に満足される。
【0124】確率密度の最大値の決定はln[p(S|
G_,Gtot )]の最大値の決定と等価である。
G_,Gtot )]の最大値の決定と等価である。
【0125】個々の測定は統計的に独立したポアッソン
ランダム変数であると仮定され、合計された全体の測定
の全てまたはその一部の合計値も同様に個々の測定の平
均値はポアッソンランダム変数であると仮定される。さ
らに、ポアッソンランダム変数の平均値は雑音のない測
定であると仮定され、即ち合計された総回転開口測定の
平均値は以下のようになり、
ランダム変数であると仮定され、合計された全体の測定
の全てまたはその一部の合計値も同様に個々の測定の平
均値はポアッソンランダム変数であると仮定される。さ
らに、ポアッソンランダム変数の平均値は雑音のない測
定であると仮定され、即ち合計された総回転開口測定の
平均値は以下のようになり、
【0126】
【数54】個々の検出器測定の平均値は次式のようにな
る。
る。
【0127】
【数55】
【0128】それ故、合計された全体の測定が与えられ
ている測定の確率の多項分布に平均値を代入すると、情
景の出現で条件付けられた対応する確率が得られ、即
ち、
ている測定の確率の多項分布に平均値を代入すると、情
景の出現で条件付けられた対応する確率が得られ、即
ち、
【0129】
【数56】積の対数は対数の合計であり、情景で条件付
けされた測定の対数は次式のようになる。 ln[p(S|G_,Gtot )]=ln[p(G_,G
tot IS)]+ln[p(S)]−ln[p(G_,G
tot )] スターリング近似を超解像度情景の投影された画素の統
計に適用し、積の対数が対数の合計であることを認識
し、適切な確率分布をバイエルの定理に代入すると次式
が得られる。
けされた測定の対数は次式のようになる。 ln[p(S|G_,Gtot )]=ln[p(G_,G
tot IS)]+ln[p(S)]−ln[p(G_,G
tot )] スターリング近似を超解像度情景の投影された画素の統
計に適用し、積の対数が対数の合計であることを認識
し、適切な確率分布をバイエルの定理に代入すると次式
が得られる。
【0130】
【数57】超スピン開口の最大事後の実施形態の最大が
条件付き確率を最大にする超解像度の情景を発見するこ
とに依存するので、条件付き確率の対数の微分から生じ
る非線形の部分的微分方程式を解く情景評価が求めら
れ、即ち、次式が得られる。
条件付き確率を最大にする超解像度の情景を発見するこ
とに依存するので、条件付き確率の対数の微分から生じ
る非線形の部分的微分方程式を解く情景評価が求めら
れ、即ち、次式が得られる。
【0131】
【数58】
【0132】先の方法に関しては、超解像度情景の解
は、情景の超解像度画素に関する条件付き確率の導関数
がゼロになることを必要とすることによって生成される
1組の非線形の結合式を満足しなければならない。この
セクションで説明している実施形態の1組の連立方程式
を以下に示す。
は、情景の超解像度画素に関する条件付き確率の導関数
がゼロになることを必要とすることによって生成される
1組の非線形の結合式を満足しなければならない。この
セクションで説明している実施形態の1組の連立方程式
を以下に示す。
【0133】
【数59】
【0134】1組の連立方程式の数値解は指数関数とピ
カールの反復的方法を適用することにより実行されるこ
とができる。
カールの反復的方法を適用することにより実行されるこ
とができる。
【0135】合計された測定の超解像度の実施形態を説
明する。この超解像度の回転開口実施形態では付加的な
超解像度式の組が回転開口フレーム中の測定された総光
子の数に基づいて設定される。評価された情景値が、測
定された回転開口フレームが与えられた情景値に対する
条件付き確率密度の極値であるのと同様に、全体の測定
の発生が与えられた情景値に対する条件付き確率密度の
極値であることを必要とすることにより超解像度の情景
値は決定される。2つの組の式は未知の平均画素値を消
去するために使用されることができ、超解像度の情景画
素値用の1組の非線形の式を残す。1組の非線形の式は
ピカールの方法を用いて解かれる。この代わりの超解像
度手段の定義に導かれる手順の詳細について以下説明す
る。
明する。この超解像度の回転開口実施形態では付加的な
超解像度式の組が回転開口フレーム中の測定された総光
子の数に基づいて設定される。評価された情景値が、測
定された回転開口フレームが与えられた情景値に対する
条件付き確率密度の極値であるのと同様に、全体の測定
の発生が与えられた情景値に対する条件付き確率密度の
極値であることを必要とすることにより超解像度の情景
値は決定される。2つの組の式は未知の平均画素値を消
去するために使用されることができ、超解像度の情景画
素値用の1組の非線形の式を残す。1組の非線形の式は
ピカールの方法を用いて解かれる。この代わりの超解像
度手段の定義に導かれる手順の詳細について以下説明す
る。
【0136】雑音のない測定を実施する別の方法を説明
する。回転開口の超解像度アルゴリズムの導出の第1の
バージョンに関しては、雑音のない画像測定はセンサシ
ステム応答関数Hm のコンボリューションとして表さ
れ、観察された情景Sまたは次式を有する。
する。回転開口の超解像度アルゴリズムの導出の第1の
バージョンに関しては、雑音のない画像測定はセンサシ
ステム応答関数Hm のコンボリューションとして表さ
れ、観察された情景Sまたは次式を有する。
【0137】
【数60】m番目の雑音のない画像測定のディスクリー
トな形態と、全てのフレームに対する全ての検出器の雑
音のない全体の測定は次式で表される。
トな形態と、全てのフレームに対する全ての検出器の雑
音のない全体の測定は次式で表される。
【0138】
【数61】
【0139】先の式では、jは検出器位置ベクトル/r
j に関連する一次元指数であり、kは検出器位置ベクト
ル/rk に関連する一次元指数である。
j に関連する一次元指数であり、kは検出器位置ベクト
ル/rk に関連する一次元指数である。
【0140】イソプラネティック光学システムと空間的
に均一な応答焦点平面アレイ検出器に対して、m番目の
フレームシステム応答関数は個々の検出器応答関数Wm
を有する光学システムの点分布関数PSFm のコンボリ
ューションとして表され、即ち次式のようになる。
に均一な応答焦点平面アレイ検出器に対して、m番目の
フレームシステム応答関数は個々の検出器応答関数Wm
を有する光学システムの点分布関数PSFm のコンボリ
ューションとして表され、即ち次式のようになる。
【0141】
【数62】
【0142】変更された最大事後方法は、超解像度の情
景の評価値を発見するため2つの異なったバイエルの確
率規則を使用する。2つのバイエル規則は次式であり、
景の評価値を発見するため2つの異なったバイエルの確
率規則を使用する。2つのバイエル規則は次式であり、
【0143】
【数63】この式は1組の回転開口測定全体が与えられ
る超解像度情景が出現する確率と、全ての回転開口測定
の和が与えられる超解像度情景が出現する確率にそれぞ
れ対応する。先に定めたように、Gは全体の測定マトリ
ックスであり、これは個々のフレーム測定マトリックス
からなるブロックマトリックスであり、Sは超解像度座
標グリッドに投影される個々の情景素子からなる列マト
リックスであり、即ち、
る超解像度情景が出現する確率と、全ての回転開口測定
の和が与えられる超解像度情景が出現する確率にそれぞ
れ対応する。先に定めたように、Gは全体の測定マトリ
ックスであり、これは個々のフレーム測定マトリックス
からなるブロックマトリックスであり、Sは超解像度座
標グリッドに投影される個々の情景素子からなる列マト
リックスであり、即ち、
【0144】
【数64】
【0145】ここでnf は測定フレームの総数であり、
npixel は超解像度情景の評価の超解像度点の総数であ
る。
npixel は超解像度情景の評価の超解像度点の総数であ
る。
【0146】情景が与えられた全体の測定に対する条件
付き確率を有する自己一貫の超解像度方法は、情景Sに
対する評価値を求め、これは回転開口測定の完全な組が
与えられる情景の条件付き確率を最大にし、全ての回転
開口測定の和が与えられる情景の条件付き確率を最大に
し、即ち、同時にp(S|G)とp(S|Gtot )を最
大にするSの解が求められる。一般的な分布では、両者
の条件を満たす極値の解は可能ではないが、以下のセク
ションで説明する手順および統計的な必要条件は非線形
の1組の方程式に関する解を許容する。
付き確率を有する自己一貫の超解像度方法は、情景Sに
対する評価値を求め、これは回転開口測定の完全な組が
与えられる情景の条件付き確率を最大にし、全ての回転
開口測定の和が与えられる情景の条件付き確率を最大に
し、即ち、同時にp(S|G)とp(S|Gtot )を最
大にするSの解が求められる。一般的な分布では、両者
の条件を満たす極値の解は可能ではないが、以下のセク
ションで説明する手順および統計的な必要条件は非線形
の1組の方程式に関する解を許容する。
【0147】自己一貫の要求は個々の検出器測定の平均
値を仮定することにより組入れられ、回転開口測定の和
の合計の平均値は雑音のない組の測定に対応し、即ち、
値を仮定することにより組入れられ、回転開口測定の和
の合計の平均値は雑音のない組の測定に対応し、即ち、
【0148】
【数65】
【0149】↑と↓は超解像度と“測定”座標空間グリ
ッドを示す。
ッドを示す。
【0150】画像センサ測定は情景から伝播する光子の
検出に対応し、情景の評価値は常に正である。それ故、
付加的な自己一貫の要求は、 S≧0, (34)であり、これは前述の超解像
度アルゴリズムにより暗黙的に満たされる。
検出に対応し、情景の評価値は常に正である。それ故、
付加的な自己一貫の要求は、 S≧0, (34)であり、これは前述の超解像
度アルゴリズムにより暗黙的に満たされる。
【0151】対数は単調な関数であるので、確率密度の
最大値の決定はln[p(SIG)]とln[p(SI
Gtot )]の確率密度の最大値を決定することに等し
い。対数をバイエルの定理に適用すると、次式の最大値
が求められる。
最大値の決定はln[p(SIG)]とln[p(SI
Gtot )]の確率密度の最大値を決定することに等し
い。対数をバイエルの定理に適用すると、次式の最大値
が求められる。
【0152】
【数66】
【0153】統計的に独立したフレーム測定に対して、
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
【0154】
【数67】
【0155】同様に、統計的に独立した検出器測定に対
して、情景が与えられるフレーム測定の発生の条件付き
確率は情景で条件付けされる個々の画素確率密度の積と
して表されることができる。
して、情景が与えられるフレーム測定の発生の条件付き
確率は情景で条件付けされる個々の画素確率密度の積と
して表されることができる。
【0156】
【数68】
【0157】全体の測定の条件付き確率密度はそれ故、
情景について条件付けされる画素とフレーム測定確率の
積により与えられる。
情景について条件付けされる画素とフレーム測定確率の
積により与えられる。
【0158】
【数69】
【0159】情景放射または各画素素子に対する光子の
反射が統計的に独立していると仮定すると、一次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
反射が統計的に独立していると仮定すると、一次元の情
景マトリックスの確率密度は各素子に関連する確率密度
関数の積として表されることができ、即ち、
【0160】
【数70】
【0161】ソースの光子変動により左右される情景統
計では、情景画素素子の確率密度関数は次式のように仮
定される。
計では、情景画素素子の確率密度関数は次式のように仮
定される。
【0162】
【数71】ここで/S(k↑)はk↓番目の投影された
“画素”の平均値である。
“画素”の平均値である。
【0163】個々の検出器測定もポアッソンランダムプ
ロセスであると仮定すると、個々の検出器測定確率密度
関数は次式のようになることが示唆される。
ロセスであると仮定すると、個々の検出器測定確率密度
関数は次式のようになることが示唆される。
【0164】
【数72】m番目のフレーム測定j↓番目の“画素”の
平均値を置換すると、
平均値を置換すると、
【0165】
【数73】この式は情景を与える検出器測定の確率の式
を導出する。
を導出する。
【0166】
【数74】
【0167】統計的に独立したポアッソンランダム変数
の和もポアッソンランダム変数であるので、合計された
全ての回転開口検出器確率の全体はポアッソンランダム
変数であり、次式により与えられる。
の和もポアッソンランダム変数であるので、合計された
全ての回転開口検出器確率の全体はポアッソンランダム
変数であり、次式により与えられる。
【0168】
【数75】合計された回転開口測定の全体の平均値は次
式のようになり、
式のようになり、
【0169】
【数76】これをポアッソン確率密度に代入すると、所
定のSが発生するGtot の条件付き確率が得られる。
定のSが発生するGtot の条件付き確率が得られる。
【0170】
【数77】
【0171】積の対数は対数の和であるので、情景で条
件付けされた測定の対数は以下のようになる。
件付けされた測定の対数は以下のようになる。
【0172】
【数78】
【0173】同様に、情景確率密度の対数は以下のよう
になる。
になる。
【0174】
【数79】
【0175】先の式を情景が与えられる1組の情景測定
に対する条件付き確率に代入すると次式が得られる。
に対する条件付き確率に代入すると次式が得られる。
【0176】
【数80】
【0177】p(S|Gtot )の対数の決定は情景が与
えられる合計された全体の測定の条件付き確率の対数を
決定することに依存する。
えられる合計された全体の測定の条件付き確率の対数を
決定することに依存する。
【0178】
【数81】
【0179】条件付き確率の関係に対してバイエルの法
則へ代入すると、次式が得られる。
則へ代入すると、次式が得られる。
【0180】
【数82】
【0181】約50を越える光子をカウントした場合、
スターリング近似である。
スターリング近似である。
【0182】ln[x!]=xln[x]−xの式は計
算手順を簡単にするために使用される。問題の2つの条
件付き確率に適用されるスターリング概算は次式を導
く。
算手順を簡単にするために使用される。問題の2つの条
件付き確率に適用されるスターリング概算は次式を導
く。
【0183】
【数83】
【0184】所望の情景画素素子S(l↑)に関して適
切な確率の対数を微分することにより、確率を最大にす
る情景の解を要求する極値方程式が得られる。微分によ
って次式が得られる。
切な確率の対数を微分することにより、確率を最大にす
る情景の解を要求する極値方程式が得られる。微分によ
って次式が得られる。
【0185】
【数84】
【0186】クロネッカー(Kronecker )のデルタ関数
指標積の合計を計算することにより式を簡単にすると、
指標積の合計を計算することにより式を簡単にすると、
【0187】
【数85】
【0188】測定が与えられる情景確率の対数の導関数
の設定に関連する1組の式を以下示す。
の設定に関連する1組の式を以下示す。
【0189】
【数86】
【0190】S(l↑)を解くことによって、1組の結
合された非線形の式が得られる。
合された非線形の式が得られる。
【0191】
【数87】
【0192】指数関数化によって次式が得られる。
【0193】
【数88】
【0194】合計された全体の測定が与えられる情景確
率の対数の導関数を設定することに関する1組の式を以
下示す。
率の対数の導関数を設定することに関する1組の式を以
下示す。
【0195】
【数89】
【0196】S(l↑)の対数を解くと、次の1組の結
合された非線形の式が得られる。
合された非線形の式が得られる。
【0197】
【数90】
【0198】指数関数化によって次式が得られる。
【0199】
【数91】
【0200】2つの組の式(35)、(36)はS(l
↑)または/S(l↑)の評価プロセスを公式化するた
めに同時に解かれることができる。
↑)または/S(l↑)の評価プロセスを公式化するた
めに同時に解かれることができる。
【0201】以上説明した超スピン開口の実施形態は最
大事後評価技術を適用する。超スピン開口プロセスに対
する別の方法は最大尤度評価と通常呼ばれるバイエルの
評価方法を使用する。最大尤度評価は目標が与えられる
測定の条件付き確率を最大にする未知の目標の値を求め
る。即ち、目標の最大尤度評価は1組の観察された回転
開口センサ測定全体に導かれる可能性が最も高い目標で
ある。
大事後評価技術を適用する。超スピン開口プロセスに対
する別の方法は最大尤度評価と通常呼ばれるバイエルの
評価方法を使用する。最大尤度評価は目標が与えられる
測定の条件付き確率を最大にする未知の目標の値を求め
る。即ち、目標の最大尤度評価は1組の観察された回転
開口センサ測定全体に導かれる可能性が最も高い目標で
ある。
【0202】前述したように、Gは全測定マトリックス
であり、これは個々のフレーム測定マトリックスからな
るブロック列マトリックスであり、以下に示す。
であり、これは個々のフレーム測定マトリックスからな
るブロック列マトリックスであり、以下に示す。
【0203】
【数92】ここでnf は捕捉されるフレームの総数であ
る。Sはランダム情景のnf 実体に関するブロック列マ
トリックスである。
る。Sはランダム情景のnf 実体に関するブロック列マ
トリックスである。
【0204】
【数93】
【0205】情景Sがランダムフィールドであるとする
と、情景平均値の最大尤度評価が計算される。最大尤度
評価は情景平均/Sの値であり、これは尤度関数P(G
|/S)を最大にする。
と、情景平均値の最大尤度評価が計算される。最大尤度
評価は情景平均/Sの値であり、これは尤度関数P(G
|/S)を最大にする。
【0206】統計的に独立したフレーム測定に対して、
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
情景が与えられる全体の測定ベクトルのジョイント確率
密度は個々のm番目のフレームジョイント確率密度の積
として表されることができる。
【0207】
【数94】
【0208】同様に、統計的に独立した検出器測定で
は、情景が与えられるフレーム発生の条件付き確率は情
景について条件付けされた個々の画素確率密度の積とし
て表されることができる。
は、情景が与えられるフレーム発生の条件付き確率は情
景について条件付けされた個々の画素確率密度の積とし
て表されることができる。
【0209】
【数95】
【0210】全体の測定条件付き確率密度はそれ故、情
景について条件付けされた画素とフレーム測定確率の積
により与えられる。
景について条件付けされた画素とフレーム測定確率の積
により与えられる。
【0211】
【数96】
【0212】個々の検出器測定がポアッソンプロセスで
あると仮定すると、個々の測定の条件付き確率は次式の
ように表される。
あると仮定すると、個々の測定の条件付き確率は次式の
ように表される。
【0213】
【数97】
【0214】条件付き確率密度を最大にする情景を発見
することは、確率密度の対数を最大にする情景を発見す
ることに等しいので、次式を最大にする/Sの解が求め
られる。
することは、確率密度の対数を最大にする情景を発見す
ることに等しいので、次式を最大にする/Sの解が求め
られる。
【0215】
【数98】積の対数が対数の和であることを認識する
と、/Sが与えられる測定マトリックスGの確率は以下
のようになる。
と、/Sが与えられる測定マトリックスGの確率は以下
のようになる。
【0216】
【数99】
【0217】式(43)を式(42)に代入すると、次
式が得られる。
式が得られる。
【0218】
【数100】条件付き確率を最大にする/S(l)値を
決定するため、ln[p(G|/S)]の導関数がゼロ
になる値/S(l)が要求される。ln[p(G|/
S)]の導関数は次式により与えられる。
決定するため、ln[p(G|/S)]の導関数がゼロ
になる値/S(l)が要求される。ln[p(G|/
S)]の導関数は次式により与えられる。
【0219】
【数101】
【0220】導関数をゼロに設定すると、/Sを解く非
線形の式が得られる。
線形の式が得られる。
【0221】
【数102】
【0222】光学システムの特性を保存する光子によ
り、ほとんどの構成例では、式の左側はnf に合計され
る。ピカールの技術を適用するため、式(44)の各側
は/S(l)により乗算され、評価は次式のように反復
して形成される。
り、ほとんどの構成例では、式の左側はnf に合計され
る。ピカールの技術を適用するため、式(44)の各側
は/S(l)により乗算され、評価は次式のように反復
して形成される。
【0223】
【数103】ここで上付き文字nは反復数を示してい
る。
る。
【0224】超スピン開口の最大尤度評価を計算するこ
の方法は予期最大アルゴリズムである。予期最大アルゴ
リズムは情景平均の最大尤度評価へ集中する。
の方法は予期最大アルゴリズムである。予期最大アルゴ
リズムは情景平均の最大尤度評価へ集中する。
【0225】以上説明したアルゴリズム手順の各反復で
は、乗算の補正係数はnf 回転開口測定の全体の組から
計算される。反復予期最大アルゴリズムは最大尤度評価
を計算するために使用されてもよい。反復アルゴリズム
は個々に各フレームに基づく評価値を更新し、それによ
って最大尤度評価を決定する計算上の必要性をなくす。
は、乗算の補正係数はnf 回転開口測定の全体の組から
計算される。反復予期最大アルゴリズムは最大尤度評価
を計算するために使用されてもよい。反復アルゴリズム
は個々に各フレームに基づく評価値を更新し、それによ
って最大尤度評価を決定する計算上の必要性をなくす。
【0226】予期される最大アルゴリズムの反復形態
は、超スピン開口に適用するためTitterington氏による
文献Journal of the Royal Statistical Society-B、46
巻、No.2、257 〜267 頁、1984年により展開された方法
に発展される。以下の反復的予期最大アルゴリズムはn
=0,1,…,、m=0,1,…,Nf −1、l=1,
2,…,Npixel で得られる。
は、超スピン開口に適用するためTitterington氏による
文献Journal of the Royal Statistical Society-B、46
巻、No.2、257 〜267 頁、1984年により展開された方法
に発展される。以下の反復的予期最大アルゴリズムはn
=0,1,…,、m=0,1,…,Nf −1、l=1,
2,…,Npixel で得られる。
【0227】
【数104】
【0228】この反復的な予期される最大アルゴリズム
は先に説明した方法を用いて初期化されるか、情景の平
均評価の各素子を回転開口フレームの光子の平均数に基
づいて一定値に初期化される。
は先に説明した方法を用いて初期化されるか、情景の平
均評価の各素子を回転開口フレームの光子の平均数に基
づいて一定値に初期化される。
【0229】
【数105】
【0230】前述の超スピン開口の実施形態は未知の情
景を評価するための2つの広い解決方法、即ち最大事後
評価と、最大尤度評価を含む。最大事後評価のカテゴリ
ー内には4つの実施形態が開示されており、それらは、
情景の現在の評価に等しい情景平均を評価することによ
り情景平均/S(l)が特別な方法で評価される最大事
後方法と、サンプル統計を用いて情景平均/S(l)が
評価され、情景Sq (l)の各実体を評価するための反
復アルゴリズムを生じ、それに続いてサンプル平均の再
評価が行われる最大事後方法と、単一の確率密度関数が
測定と1組の測定された回転開口フレームの総光子とで
公式化される最大事後方法と、別々の確率密度関数が1
組の測定された回転開口フレームの総光子で公式化さ
れ、それに続いて2つの組の非線形の式を同時的に解く
最大事後方法とを含んでいる。
景を評価するための2つの広い解決方法、即ち最大事後
評価と、最大尤度評価を含む。最大事後評価のカテゴリ
ー内には4つの実施形態が開示されており、それらは、
情景の現在の評価に等しい情景平均を評価することによ
り情景平均/S(l)が特別な方法で評価される最大事
後方法と、サンプル統計を用いて情景平均/S(l)が
評価され、情景Sq (l)の各実体を評価するための反
復アルゴリズムを生じ、それに続いてサンプル平均の再
評価が行われる最大事後方法と、単一の確率密度関数が
測定と1組の測定された回転開口フレームの総光子とで
公式化される最大事後方法と、別々の確率密度関数が1
組の測定された回転開口フレームの総光子で公式化さ
れ、それに続いて2つの組の非線形の式を同時的に解く
最大事後方法とを含んでいる。
【0231】最大尤度評価のカテゴリー内では2つの実
施形態が説明されている。第1の実施形態は情景の全体
平均/S(l)の評価を形成するために1組の測定され
た回転開口フレームを使用し、ここでは回転開口フレー
ムの組全体は評価のそれぞれの更新で使用される。第2
の最大尤度実施例では、反復的アルゴリズムが使用さ
れ、ここでは情景の全体平均は、一度に1つの回転開口
フレーム測定を用いて更新される。各最大尤度実施形態
は予期最大手順であり、それ故ピカール反復が最大尤度
評価へ集中することを保証する。
施形態が説明されている。第1の実施形態は情景の全体
平均/S(l)の評価を形成するために1組の測定され
た回転開口フレームを使用し、ここでは回転開口フレー
ムの組全体は評価のそれぞれの更新で使用される。第2
の最大尤度実施例では、反復的アルゴリズムが使用さ
れ、ここでは情景の全体平均は、一度に1つの回転開口
フレーム測定を用いて更新される。各最大尤度実施形態
は予期最大手順であり、それ故ピカール反復が最大尤度
評価へ集中することを保証する。
【0232】説明を完全にするため、図17は放射状の
試験パターンの超解像度の画像を生成するために使用さ
れる超スピン開口合成方法30の実施結果を示している。
図18は超スピン開口合成動作の例を示している。動作
の定量的測定は目標の既知の基準パワースペクトルと超
解像度画像のパワースペクトルとの間の相関係数を決定
することにより得られる。図19は真の情景が水平バー
ターゲットである超スピン開口合成の例を示している。
動作の定量的測定は物体の既知の基準パワースペクトル
と超解像度画像のパワースペクトルとの間の相関係数を
決定することにより発見される。
試験パターンの超解像度の画像を生成するために使用さ
れる超スピン開口合成方法30の実施結果を示している。
図18は超スピン開口合成動作の例を示している。動作
の定量的測定は目標の既知の基準パワースペクトルと超
解像度画像のパワースペクトルとの間の相関係数を決定
することにより得られる。図19は真の情景が水平バー
ターゲットである超スピン開口合成の例を示している。
動作の定量的測定は物体の既知の基準パワースペクトル
と超解像度画像のパワースペクトルとの間の相関係数を
決定することにより発見される。
【0233】以上、回転ストリップ(部分的)開口画像
放射線センサ装置と、複数の回転ストリップ開口画像測
定からの超解像度情景評価を合成する方法を説明した。
所望の実施形態は本発明の原理の応用を表した多数の特
定の実施形態のうちの数例にすぎないことが理解されよ
う。明らかに、多数のおよび他の装置が本発明の技術範
囲を逸脱することなく当業者により容易に行われること
ができる。
放射線センサ装置と、複数の回転ストリップ開口画像測
定からの超解像度情景評価を合成する方法を説明した。
所望の実施形態は本発明の原理の応用を表した多数の特
定の実施形態のうちの数例にすぎないことが理解されよ
う。明らかに、多数のおよび他の装置が本発明の技術範
囲を逸脱することなく当業者により容易に行われること
ができる。
【図面の簡単な説明】
【図1】本発明の原理による回転開口画像放射線センサ
システムの概略図。
システムの概略図。
【図2】本発明の原理によるシステムおよび方法の第1
の実施形態の概略図。
の実施形態の概略図。
【図3】本発明の第2の実施形態の概略図。
【図4】本発明の第3の実施形態の概略図。
【図5】本発明の第4の実施形態の概略図。
【図6】本発明の第5の実施形態の概略図。
【図7】図2乃至6の本発明に関連する超解像度処理方
法の実施形態の概略図。
法の実施形態の概略図。
【図8】図2乃至6の本発明に関連する超解像度処理方
法の実施形態の概略図。
法の実施形態の概略図。
【図9】図2乃至6の本発明に関連する超解像度処理方
法の別の実施形態の概略図。
法の別の実施形態の概略図。
【図10】図2乃至6の本発明に関連する超解像度処理
方法の別の実施形態の概略図。
方法の別の実施形態の概略図。
【図11】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図12】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図13】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図14】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図15】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図16】図2乃至6の本発明に関連する別の超解像度
処理方法の別の実施形態の概略図。
処理方法の別の実施形態の概略図。
【図17】試験パターンに関する本発明を用いた試験結
果の図。
果の図。
【図18】別の試験パターンに関する本発明を用いた試
験結果の図。
験結果の図。
【図19】さらに別の試験パターンに関する本発明を用
いた試験結果の図。
いた試験結果の図。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ジェラルド・エル・ラファネリ アメリカ合衆国、カリフォルニア州 92708、ファウンテイン・バレー、マウン ト・アコマ・サークル 16645 (72)発明者 フィリップ・ジェイ・セメンティリ アメリカ合衆国、アリゾナ州 85718、タ クソン、イー・カレ・エレーナ 1201 (72)発明者 スーザン・ビー・マウント アメリカ合衆国、カリフォルニア州 90503、トランス、エメラルド・ストリー ト・ナンバー75 3922 (72)発明者 アルバート・エム・ビスビー アメリカ合衆国、アリゾナ州 85749、タ クソン、エヌ・ロングホーン・ドライブ 3101 (72)発明者 ジェイムズ・エフ・モンゴメリー アメリカ合衆国、カリフォルニア州 90254、ハーモサ・ビーチ、マンハッタ ン・アベニュー 1022 (72)発明者 スティーブン・ケー・ジョンソン アメリカ合衆国、バージニア州 22124、 ハーンドン、ヒッコリー・ヒルズ・ドライ ブ 3412
Claims (3)
- 【請求項1】 予め定められた空間周波数カットオフを
有し、光軸を中心に回動して情景の時間的に順次連続す
る映像を生成する回転ストリップ開口を具備する望遠鏡
と、 この望遠鏡の焦点平面に配置され1組のストリップ開口
映像を生成する2次元の検出器アレイと、 この検出器アレイに結合され、ストリップ開口が望遠鏡
の光軸を中心に回動するとき1組のストリップ開口映像
からなる映像情景の複数の映像フレームを記録し、回転
するストリップ開口の空間周波数カットオフよりも大き
い空間周波数に対する記録された映像から観察された情
景の評価を合成する評価プロセッサを備えている信号プ
ロセッサとを具備していることを特徴とする回転ストリ
ップ開口放射線センサシステム。 - 【請求項2】 信号プロセッサは、 情景に対応する過去の知識を記憶する映像合成データベ
ースと、 帯域制限された映像合成プロセッサと映像合成データベ
ースとの間に結合されて映像から情景内容および構造情
報を抽出し、そのデータベースに含まれている過去の情
報を変更する訓練アルゴリズムプロセッサとを具備し、
前記情報は信号プロセッサによって処理されて超解像度
情景の評価を生成する請求項1記載のシステム。 - 【請求項3】 光軸を中心に回動する回転ストリップ開
口望遠鏡と、この望遠鏡の焦点平面に位置された映像を
検出する2次元の検出器アレイと、検出器アレイの検出
器の積分時間中に静止映像を生成する回転補償手段と、
望遠鏡の光軸を中心に回動するとき望遠鏡によって形成
された情景の複数の映像フレームを記憶し、記憶された
映像フレームから全円形開口映像を合成する信号プロセ
ッサとを具備している映像システムを使用して映像を合
成する映像合成方法において、 1組のストリップ開口映像を発生し、 現在の評価の非線形関数と超解像度情景の現在の評価と
の積を形成することにより超解像度の全円形開口映像ス
ペクトルを合成することを特徴とする映像合成方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US635073 | 1996-04-19 | ||
| US08/635,073 US6155704A (en) | 1996-04-19 | 1996-04-19 | Super-resolved full aperture scene synthesis using rotating strip aperture image measurements |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH1063821A true JPH1063821A (ja) | 1998-03-06 |
Family
ID=24546332
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP9117621A Pending JPH1063821A (ja) | 1996-04-19 | 1997-04-21 | 回転ストリップ開口映像測定を使用する超解像度全円形開口情景合成 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6155704A (ja) |
| EP (1) | EP0802425B1 (ja) |
| JP (1) | JPH1063821A (ja) |
| DE (1) | DE69700718T2 (ja) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6466894B2 (en) | 1998-06-18 | 2002-10-15 | Nec Corporation | Device, method, and medium for predicting a probability of an occurrence of a data |
| US20030068097A1 (en) * | 2001-06-15 | 2003-04-10 | Massachusetts Institute Of Technology | Adaptive mean estimation and normalization of data |
| US7123780B2 (en) * | 2001-12-11 | 2006-10-17 | Sony Corporation | Resolution enhancement for images stored in a database |
| US6603424B1 (en) * | 2002-07-31 | 2003-08-05 | The Boeing Company | System, method and computer program product for reducing errors in synthetic aperture radar signals |
| ES2392512T3 (es) * | 2004-01-27 | 2012-12-11 | Saab Ab | Método para la reducción del desenfoque angular en imágenes de radar |
| US7064702B1 (en) | 2005-03-01 | 2006-06-20 | The Boeing Company | System, method and computer program product for reducing quadratic phase errors in synthetic aperture radar signals |
| JP4814840B2 (ja) * | 2007-05-23 | 2011-11-16 | オリンパス株式会社 | 画像処理装置又は画像処理プログラム |
| CN101261319B (zh) * | 2008-04-18 | 2010-12-01 | 华中科技大学 | 一种综合孔径辐射计成像校正方法 |
| JP2012061141A (ja) * | 2010-09-16 | 2012-03-29 | Canon Inc | 被検体情報取得装置及び被検体情報取得方法 |
| US8362947B2 (en) * | 2011-03-29 | 2013-01-29 | Victor Gorelik | Method for obtaining object-plane field from its two images |
| US20140027577A1 (en) * | 2012-07-26 | 2014-01-30 | Alliant Techsystems Inc. | Imaging satellite system and method |
| US10154193B1 (en) | 2016-03-08 | 2018-12-11 | National Technology & Engineering Solutions Of Sandia, Llc | Noncircular aperture imaging system |
| CN108512608B (zh) * | 2018-04-09 | 2023-11-28 | 中国科学技术大学 | 一种基于多孔径合成的量子通信接收装置 |
| CN108896991B (zh) * | 2018-04-26 | 2021-04-13 | 西安空间无线电技术研究所 | 一种基于数据融合的星载分布式综合孔径微波辐射计系统 |
| DE102020107965B3 (de) * | 2020-03-23 | 2021-09-09 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung eingetragener Verein | Verfahren zur optischen Bestimmung einer Intensitätsverteilung |
| CN119179066B (zh) * | 2024-11-25 | 2025-02-07 | 中国人民解放军国防科技大学 | 多孔径成像系统协同定位点标校方法及装置 |
Family Cites Families (17)
| 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 |
| US5160934A (en) * | 1984-08-27 | 1992-11-03 | The United States Of America As Represented By The Secretary Of The Navy | Cross-switched MICRAD seeker |
| US4739332A (en) * | 1984-09-17 | 1988-04-19 | Hughes Aircraft Company | Fanbeam inversion radar |
| US4825215A (en) * | 1986-07-03 | 1989-04-25 | Hughes Aircraft Company | Radiometric imager having a frequency-dispersive linear array antenna |
| US4876602A (en) * | 1988-05-02 | 1989-10-24 | Hughes Aircraft Company | Electronic focus correction by signal convolution |
| DE3922428A1 (de) * | 1989-07-07 | 1991-01-17 | Deutsche Forsch Luft Raumfahrt | Verfahren zur extraktion von bewegungsfehlern eines ein kohaerentes abbildungsradarsystem mitfuehrenden traegers aus radar-rohdaten und einrichtung zur durchfuehrung des verfahrens |
| DE4007611C1 (ja) * | 1990-03-09 | 1991-05-16 | Deutsche Forschungsanstalt Fuer Luft- Und Raumfahrt Ev, 5300 Bonn, De | |
| US5161204A (en) * | 1990-06-04 | 1992-11-03 | Neuristics, Inc. | Apparatus for generating a feature matrix based on normalized out-class and in-class variation matrices |
| US5072226A (en) * | 1990-06-07 | 1991-12-10 | Hughes Aircraft Company | Radiometer system incorporating a cylindrical parabolic reflector and minimum redundancy array feed |
| US5122805A (en) * | 1991-02-06 | 1992-06-16 | Radian Corporation | Radio acoustic sounding system for remotely determining atmospheric temperature profiles |
| US5128530A (en) * | 1991-05-28 | 1992-07-07 | Hughes Aircraft Company | Piston error estimation method for segmented aperture optical systems while observing arbitrary unknown extended scenes |
| US5164730A (en) * | 1991-10-28 | 1992-11-17 | Hughes Aircraft Company | Method and apparatus for determining a cross-range scale factor in inverse synthetic aperture radar systems |
| US5122803A (en) * | 1991-11-06 | 1992-06-16 | The United States Of America As Represented By The Secretary Of The Army | Moving target imaging synthetic aperture radar |
| 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 |
| EP0634669B1 (de) * | 1993-07-15 | 1997-09-24 | Daimler-Benz Aerospace Aktiengesellschaft | Verfahren zur Klassifikation eines Gegenstandes und Verwendung des Verfahrens |
-
1996
- 1996-04-19 US US08/635,073 patent/US6155704A/en not_active Expired - Lifetime
-
1997
- 1997-04-19 DE DE69700718T patent/DE69700718T2/de not_active Expired - Lifetime
- 1997-04-19 EP EP97106509A patent/EP0802425B1/en not_active Expired - Lifetime
- 1997-04-21 JP JP9117621A patent/JPH1063821A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| EP0802425A2 (en) | 1997-10-22 |
| EP0802425B1 (en) | 1999-11-03 |
| DE69700718T2 (de) | 2000-06-21 |
| DE69700718D1 (de) | 1999-12-09 |
| US6155704A (en) | 2000-12-05 |
| EP0802425A3 (en) | 1997-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5243351A (en) | Full aperture image synthesis using rotating strip aperture image measurements | |
| JP6570991B2 (ja) | 大口径望遠鏡における、非同一平面上の(anisoplanatic)画像形成のためのレンズレット、ビームウォーク(beamwalk)、及びチルトの多様化 | |
| US7933464B2 (en) | Scene-based non-uniformity correction and enhancement method using super-resolution | |
| JP2889210B2 (ja) | 回転ストリップ開口映像測定を使用する運動物体および過渡事象検出 | |
| EP0802425B1 (en) | Full aperture scene synthesis using rotating strip aperture image measurements | |
| EP0590846B1 (en) | Image processor for target detection and tracking | |
| Hunziker et al. | PCA-based approach for subtracting thermal background emission in high-contrast imaging data | |
| US20050045801A1 (en) | State space wavefront reconstructor for an adaptive optics control | |
| CN112261315B (zh) | 基于相机阵列孔径合成的高分辨率计算成像系统及方法 | |
| US20030031382A1 (en) | Process for the stabilization of the images of a scene correcting offsets in grey levels, detection of mobile objects and harmonization of two snapshot capturing apparatuses based on the stabilization of the images | |
| US7523078B2 (en) | Bayesian approach for sensor super-resolution | |
| US20240394895A1 (en) | Initial orbit determination using angular velocity and angular acceleration measurements | |
| CN114037636B (zh) | 自适应光学系统校正图像的多帧盲复原方法 | |
| Keshk et al. | Performance evaluation of quality measurement for super-resolution satellite images | |
| US11875483B2 (en) | Method and device for removing remanence in an infrared image of a changing scene | |
| CN117761716B (zh) | 融合CR和GNSS的InSAR滑坡形变时序监测方法和存储介质 | |
| Schanz et al. | 3D particle position determination and correction at high particle densities | |
| 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 | |
| Kamalabadi | Multidimensional image reconstruction in astronomy | |
| CN119146932B (zh) | 一种针对甚长基线干涉阵列数据的天体成像方法及装置 | |
| CN119738049B (zh) | 一种针对大口径红外相机的星地联合辐射定标方法、装置 | |
| Řeřábek et al. | Processing of the astronomical image data obtained from UWFC optical systems | |
| Cooley et al. | Is a Picture Worth 1000 Words: Characterization Potential of Commercial Non-Earth Imaging of Satellites | |
| Duan et al. | Saturation energy restoration technology for infrared point targets | |
| Qiu et al. | Design and Verification of Super-Resolution Reconstruction Algorithm for Remote Sensing |