JP6132978B2 - 壁背後のシーンを再構成するためのシステム - Google Patents
壁背後のシーンを再構成するためのシステム Download PDFInfo
- Publication number
- JP6132978B2 JP6132978B2 JP2016513154A JP2016513154A JP6132978B2 JP 6132978 B2 JP6132978 B2 JP 6132978B2 JP 2016513154 A JP2016513154 A JP 2016513154A JP 2016513154 A JP2016513154 A JP 2016513154A JP 6132978 B2 JP6132978 B2 JP 6132978B2
- Authority
- JP
- Japan
- Prior art keywords
- wall
- sparse
- signal
- scene
- received signal
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 claims description 27
- 238000007906 compression Methods 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 7
- 238000011069 regeneration method Methods 0.000 claims description 6
- 230000008929 regeneration Effects 0.000 claims description 3
- 238000003491 array Methods 0.000 description 17
- 238000003384 imaging method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 238000013461 design Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002040 relaxant effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004613 tight binding model Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000001080 multi-layer soft lithography Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Description
この発明は、包括的には、壁透過撮像(Through-the-wall-imaging：スルーザウォール撮像)に関し、より詳細には、圧縮検知及びＭＩＭＯアンテナアレイを用いて、壁背後のシーンを再構成することに関する。
壁透過撮像
壁透過撮像(ＴＷＩ)を用いて、包囲された構造物内部の物体を外部から検出することができる。ＴＷＩでは、送信機が電磁(ＥＭ)レーダーパルスを発し、パルスは壁を通って伝搬する。パルスは、壁の逆側で物体によって反射され、次に、発せられたパルスが重畳されたインパルス応答として受信機に戻って伝搬する。通常、送信機及び受信機はアンテナアレイを用いる。
壁透過撮像(ＴＷＩ)を用いて、包囲された構造物内部の物体を外部から検出することができる。ＴＷＩでは、送信機が電磁(ＥＭ)レーダーパルスを発し、パルスは壁を通って伝搬する。パルスは、壁の逆側で物体によって反射され、次に、発せられたパルスが重畳されたインパルス応答として受信機に戻って伝搬する。通常、送信機及び受信機はアンテナアレイを用いる。
受信信号は、壁の誘電率及び透磁率にもよるが、多くの場合、壁からの間接的な二次反射で破損し、この二次反射の結果、ノイズとして現れる画像内のゴーストアーチファクトが生じる。壁クラッター低減(wall clutter reduction：壁クラッター抑圧)技法は、ＴＷＩでのマルチパス反射から生じるアーチファクトをなくそうとする。
圧縮検知
スパースなアンダーサンプリング処理によるレーダーアレイシステムは、圧縮検知(ＣＳ：compressive sensing)並びに他のサブナイキストサンプリング方法及び捕捉方法を用いることができる。そのアンテナアレイは、従来のアンテナ構造に比べて著しく少ないアレイ素子を用いてレーダー信号を捕捉及び撮像できるようにし、それにより、アレイ実現コストを著しく削減する。
スパースなアンダーサンプリング処理によるレーダーアレイシステムは、圧縮検知(ＣＳ：compressive sensing)並びに他のサブナイキストサンプリング方法及び捕捉方法を用いることができる。そのアンテナアレイは、従来のアンテナ構造に比べて著しく少ないアレイ素子を用いてレーダー信号を捕捉及び撮像できるようにし、それにより、アレイ実現コストを著しく削減する。
スパースアレイは、アレイ処理のためのナイキスト間隔である、送信信号の半波長よりはるかに長い平均素子間間隔を有する。これは、グレーティングローブとして知られる根本的に解決不可能な曖昧性を排除する不均一な素子間隔を用いて達成される。
捕捉された画像を再生するために従来の方法が用いられてきたが、従来の方法は、それらのアレイによって示されるサイドローブが増えるという難点がある。しかしながら、スパース再生方法はサイドローブに対してロバストであるので、著しく少ない数のアレイ素子を用いて撮像できるようになる。本明細書において用いられるときに、「スパース性」は相対的な用語ではなく、むしろ、大部分が０値であり、ごく少数の非０値を有するデータを指すために用いられる専門用語である。
２０１３年７月２２日にMansourによって出願された「Method and System for Through-the-Wall Imaging using Sparse Inversion for Blind Multi-Path Elimination」と題する米国特許出願第１３／９４７，４２６号では、壁を通してパルスを送信することによって、壁背後のシーン内のターゲットが検出される。反射パルスに対応する受信信号に適用されるスパースな正規化最小二乗反転によって、一次インパルス応答が検出される。一次インパルス応答を受信信号内の類似のインパルス応答に一致させる遅延演算子も求められる。パルスが壁を通り抜けた後に、ターゲットによって反射される前のパルスの歪みも求めることができる。その歪みは、ターゲットの検出を改善し、ゴースティングアーティファクトを抑圧するために、反復プロセスにおいて用いられる。
圧縮検知(ＣＳ)及びスパースアレイ処理が、レーダー撮像システムを改善する新たな手法を提供する。
この発明の実施の形態は、壁透過撮像(ＴＷＩ)のコストを著しく削減し、その複雑さを大幅に緩和する多入力多出力(ＭＩＭＯ)レーダーアレイを使用する。
それらの実施の形態は、層状の無損失壁の存在時に、入れ子アレイ、互いに素なアレイ(co-prime array)及びランダムアレイを検討する。壁パラメーターを考慮に入れるスパース再構成問題と併せて壁パラメーター推定問題を定式化し、解くことによって、シーン再構成が実行される。
ＭＩＭＯアーキテクチャは、波形ダイバーシティに起因して低減されたアレイ利得を示し、より精細な空間分解能、より高い自由度、パラメーター識別可能性及びマルチパス除去に関して改善された性能を提供する。
本方法の説明は、撮像されるシーンがスパースであると仮定して、点広がり関数(ＰＳＦ)の特性を検討することによって、種々のスパースアレイアーキテクチャの撮像性能及び壁プロファイルをメインローブ及びサイドローブ構造の観点から解析する。ＰＳＦは、ビームパターンとしても知られており、スパース再生及び圧縮検知との関連において相互コヒーレンスに密接に関連付けられる。ＰＳＦの特徴は、スパース再構成の場合にいずれもアレイ性能に関して非常に良好な直観を与える。
そのアーキテクチャの増加したサードローブレベルが、従来の撮像方法の性能を低下させるので、対象となるシーンのスパース性を利用するために、スパース再構成が用いられる。一実施の形態では、反復ハードしきい値法(ＩＨＴ：iterative hard thresholding)を適用する。ＩＨＴは、壁背後の反射率マップを推定する、貪欲法に基づくスパース信号再生法である。
さらに、それらの実施の形態は、受信データから、壁プロファイルのパラメーターを推定する方法を提供する。そのプロファイルは、壁の誘電率及び透磁率、並びに厚さを含む。これらのパラメーターは、この発明の再構成法のための撮像演算子を開発するために用いられる。
システム設定
図１Ａ及び図１Ｂに示されるように、この発明の実施の形態は、シーン形状のいかなる予備知識も必要としない、物体５０の壁透過撮像(ＴＷＩ)のための方法及びシステムを提供する。その方法は、壁４０の背後のシーン６０を再構成することができる。
図１Ａ及び図１Ｂに示されるように、この発明の実施の形態は、シーン形状のいかなる予備知識も必要としない、物体５０の壁透過撮像(ＴＷＩ)のための方法及びシステムを提供する。その方法は、壁４０の背後のシーン６０を再構成することができる。
そのシステムは、多入力多出力(ＭＩＭＯ)アンテナアレイ１０と、送受信機２０と、プロセッサ３０とを含む。アンテナは、複数の素子１１を備える。従来技術では、素子の間隔は通常均一であり、送信信号の半波長に等しい。この発明の幾つかの実施の形態では、アンテナ素子の平均素子間間隔は不均一であり、送信信号の半波長よりはるかに大きい。
送受信機は、アンテナアレイのアンテナ素子１１のうちの幾つか又は全てを用いて、１つ又は複数のパルス１４を送信する。送信されたパルスは壁４０を通って伝搬し、壁４０の背後のシーン６０内の存在し得る物体５０によって反射される。各パルスに対応する反射信号(インパルス応答)１２は、後に説明されるように、アレイ１０の素子によって受信される。受信信号は、直接パスを介して受信された一次反射と、マルチパスによって受信された二次反射とを含む。幾つかの実施の形態では、アンテナ素子は、パルスを送信するためにのみ、又はパルスを受信するためにのみ、又はパルスの送信及び受信の両方のために用いることができることに留意されたい。
受信信号１２は、物体５０を含むシーン６０を再構成する画像７０を生成するために、方法２００によって処理される。その方法は、当該技術分野において既知であるように、バスを用いてメモリ及び入力／出力インターフェースに接続されるプロセッサ３０において実行することができる。
図１Ｃに示されるように、特に関心があるのは、再構成を混乱させる可能性がある、壁に起因する間接的な二次反射８０である。それゆえ、最初に、壁のパラメーターを推定し、そのパラメーターを用いて、そのシーンが壁を通っていかに反射されるかのモデルを構築する。その後、そのモデルを用いて、そのシーンのスパース再生を行う。
シーン再構成
図２に示されるように、壁４０の背後のシーン６０は、壁を通してシーンの中に信号１４を送信し、反射信号１２から、壁の誘電率及び透磁率のパラメーターを推定すること(２１０)によって画像７０として再構成される。壁のモデル２００が、それらのパラメーターを用いて生成される(２２０)。その後、そのモデル、スパース再生及び圧縮検知２３０を用いて、反射信号からシーンが再構成される。
図２に示されるように、壁４０の背後のシーン６０は、壁を通してシーンの中に信号１４を送信し、反射信号１２から、壁の誘電率及び透磁率のパラメーターを推定すること(２１０)によって画像７０として再構成される。壁のモデル２００が、それらのパラメーターを用いて生成される(２２０)。その後、そのモデル、スパース再生及び圧縮検知２３０を用いて、反射信号からシーンが再構成される。
ＴＷＩ信号モデル
ＭＩＭＯレーダーアレイ１０が原点１３に位置し、壁４０からのｄ０のスタンドオフ距離にある２Ｄ撮像シナリオを仮定する。Ｍｔ個の送信機(Ｔｘ)素子及びＭｒ個の受信機(Ｒｘ)素子の位置はそれぞれ、ｔｉ、ｉ＝１，…，Ｍｔ、及びｒｉ、ｉ＝１，…，Ｍｒである。
ＭＩＭＯレーダーアレイ１０が原点１３に位置し、壁４０からのｄ０のスタンドオフ距離にある２Ｄ撮像シナリオを仮定する。Ｍｔ個の送信機(Ｔｘ)素子及びＭｒ個の受信機(Ｒｘ)素子の位置はそれぞれ、ｔｉ、ｉ＝１，…，Ｍｔ、及びｒｉ、ｉ＝１，…，Ｍｒである。
点ターゲット近似を用いて、直接壁反射及び付加観測雑音の影響を除く受信散乱場は、周波数領域において以下のように書くことができる。
式(１)において、ｗ(ω)は、送信レーダー波形の周波数シグネチャを表し、ｓ(ｐ)は、ｐ(ｘ，ｙ)に位置する対象物体の反射率を表し、Ｓは、撮像される領域６０を表す。関数ｇ(ｐ１，ｐ２，ω)は、ｐ１からｐ２までの層状媒質のためのグリーン関数を表し、このグリーン関数は、壁の厚さｄ及び比誘電率εの関数である。グリーン関数は、指定された初期条件又は境界条件の場合に、或る定義域において定義される非同次微分方程式のインパルス応答である。
システムを離散化するために、Ｐ個の点からなる格子を用いて領域Ｓを分割し、
を用いてマップの複素反射率を表す。各Ｒｘ素子において、Ｎ個の周波数サンプルが得られる場合には、式(１)を離散化したバージョンは、
である。ただし、
であり、ここで、
である。行列Φは、アレイ多様体行列とも呼ばれる。
スパースアレイ設計
スパースアレイアーキテクチャ
実施の形態においてスパースアレイ設計は、それぞれＭｒ個及びＭｔ個の等間隔に離間された存在し得るＴｘアレイ素子及びＲｘアレイ素子からなる概念的な格子から開始する。この格子は、互いに素、入れ子又はランダムの各アーキテクチャに従ってサブサンプリングされ、実際のＴｘアンテナ素子又はＲｘアンテナ素子を含むように、数個の格子点のみが選択される。
スパースアレイアーキテクチャ
実施の形態においてスパースアレイ設計は、それぞれＭｒ個及びＭｔ個の等間隔に離間された存在し得るＴｘアレイ素子及びＲｘアレイ素子からなる概念的な格子から開始する。この格子は、互いに素、入れ子又はランダムの各アーキテクチャに従ってサブサンプリングされ、実際のＴｘアンテナ素子又はＲｘアンテナ素子を含むように、数個の格子点のみが選択される。
図３Ａに示されるように、互いに素なアレイが、Ｔｘアレイ及びＲｘアレイの場合にそれぞれ、一対の素数(チルダ)Ｍｔおよび(チルダ)Ｍｒによって規定される。Ｔｘアレイは、(チルダ)Ｍｒ格子単位の素子間間隔で(チルダ)Ｍｔ個の素子を含み、一方、Ｒｘアレイは、(チルダ)Ｍｔ格子単位の素子間間隔で(チルダ)Ｍｒ個の素子を含む。
図４Ａに示されるように、入れ子アレイも２つの均一な線形アレイ(ＵＬＡ)を含み、ただし、Ｔｘアレイは、１格子単位の素子間間隔で(チルダ)Ｍｔ個の素子を含み、Ｒｘアレイは、(チルダ)Ｍｔ単位の間隔で(チルダ)Ｍｒ個の素子を含む。
同じ開口を有するランダムアレイは、一様分布を用いて各格子から(チルダ)Ｍｔ個のＴｘ素子と、(チルダ)Ｍｒ個のＲｘ素子とをランダムに選択することによって設計される。
ＭＩＭＯ素子の全数(チルダ)Ｍｒ＋(チルダ)Ｍｔを与えられたとすると、自由度(チルダ)Ｍｒ(チルダ)Ｍｔを最大化することによって、最適なＭＩＭＯスパース入れ子アレイを得ることができる。互いに素なアレイの場合、付加的な素数も含まれる。
図３Ｂ及び図４Ｂはそれぞれ、(チルダ)Ｍｒ＝４、(チルダ)Ｍｔ＝５の場合の互いに素なアレイ及び入れ子アレイのためのＭＩＭＯアレイビームパターンの例を示す。図３Ｂ及び図４Ｂは、送信信号のためのビームパターン３０１と、受信信号のためのビームパターン３０２と示し、実線３０３は全生成ビームパターンを示す。
(チルダ)Ｍｒ×(チルダ)Ｍｔ個のスパースＭＩＭＯアレイは、Ｍｒ×Ｍｔ個の最大ＭＩＭＯアレイのサブサンプリングと見なすことができる。これは、サブサンプリング行列
を用いて表すことができる。最大アレイ及びスパースアレイの多様体行列を表すために、それぞれΦ及び(チルダ)Φを用いるとき、スパースアレイのための捕捉関数(２)は
のようになる。ただし、(チルダ)ｙはサブサンプリングされた受信データを表す。
アレイ設計特性
アレイ設計の特性を考えるとき、従来のアレイ技法は、アレイの点広がり関数(ＰＳＦ)、すなわち、ビームパターンに焦点を合わせる。適切には正規化されたＰＳＦは、多様体行列の列間の相互コヒーレンスに等価であり、圧縮検知(ＣＳ)捕捉システムにおいて関心がある重要な特性である。２つの列間のコヒーレンスは、列間の正規化された内積と定義される。行列のコヒーレンスは、行列内の全ての要素対の中の内積の最大絶対値と定義される。
アレイ設計の特性を考えるとき、従来のアレイ技法は、アレイの点広がり関数(ＰＳＦ)、すなわち、ビームパターンに焦点を合わせる。適切には正規化されたＰＳＦは、多様体行列の列間の相互コヒーレンスに等価であり、圧縮検知(ＣＳ)捕捉システムにおいて関心がある重要な特性である。２つの列間のコヒーレンスは、列間の正規化された内積と定義される。行列のコヒーレンスは、行列内の全ての要素対の中の内積の最大絶対値と定義される。
必須ではないが、最悪時のスパース再構成保証を与えるために、低い行列コヒーレンスμ((チルダ)Φ)で十分である。
一方、ＰＳＦによって記述されるコヒーレンス構造
は、従来の方法下で特に、分解能、雑音及び干渉ロバスト性のような、アレイの性能についてはるかに多くの情報を与え、撮像された領域内の点は、潜在的に再構成曖昧性を引き起こす可能性がある。
考慮する性能指数は、メインローブ面積(ＭＬＡ)及び最大サイドローブレベル(ＭＳＬ)である。ＭＬＡは、ＰＳＦが一定のレベル、通常−３ｄＢより高いシーン内の点の周囲の面積と定義される。ＭＳＬは、ＰＳＦがそのサイドローブにおいて、すなわち、メインローブ外で達する最も高いレベルと定義される。ＭＬＡは、そのシーン内の点の周囲の曖昧性を表すので、アレイの分解能の指標である。ＭＳＬは、その点とシーン内の他の点との最大相互コヒーレンスを測定するので、特定のシーン点の再生可能性の指標である。
シーン再構成
壁プロファイル推定
式(１)のグリーン関数ｇ(ｐ１，ｐ２，ω)を求めるために、壁内の全Ｌ個の層に関する誘電率εｌ及び厚さｄl、ｌ＝１，…，Ｌを求める必要がある。壁の形状及び反射特性はあらかじめわかっていないので、誘電率及び厚さは、受信信号１２から捕捉されたデータから推定される。
壁プロファイル推定
式(１)のグリーン関数ｇ(ｐ１，ｐ２，ω)を求めるために、壁内の全Ｌ個の層に関する誘電率εｌ及び厚さｄl、ｌ＝１，…，Ｌを求める必要がある。壁の形状及び反射特性はあらかじめわかっていないので、誘電率及び厚さは、受信信号１２から捕捉されたデータから推定される。
Ｔｘアレイ素子及びＲｘアレイ素子間の自己結合を除くとき、全てのＲｘ素子からの受信信号は、壁層ごとのマルチパス６０成分を含む。バイスタティックＴｘ及びＲｘ素子がΔ＝||ｒ−ｔ||だけ離間されるが、壁からの同じスタンドオフ距離ｄ０を共有すると仮定する。スネルの法則を用いるとき、第ｌの層からの反射は、τｌ(Δ)だけ遅れて到達し、すなわち、到達時刻(ＴＯＡ)は、
であり、ｒｉは各層内の一方向の移動距離である。
このようにして、θ＝｛ε１，…，εＬ，ｄ１，…，ｄＬ｝を用いてまとめて表される未知の壁パラメーターは、各層からの壁パラメーターを与えられたとすると、各反射の測定ＴＯＡ、τｌ(Δ)、ｌ＝１，…，Ｌと、予測されるＴＯＡ(ハット)τｌ(Δθ)との間の平均二乗誤差を最小化することによって得ることができる。
ただし、αｊは各Ｔｘ−Ｒｘ分離に割り当てられる重みであり、(チルダ)Ｍｔ(チルダ)ＭｒはＭＩＭＯレーダーアレイからのＴｘ−Ｒｘ分離の全数である。限られたＲＦ帯域幅、又は低いＳＮＲの応用例の場合、より厳密なＴＯＡ推定値を得るために、超分解能又は適応技法を適用することができる。
スパース画像再生
式(６)において、測定値(チルダ)ｙからシーン反射率(ハット)ｓを再生するために、そのシーンがスパースであり、ＣＳ技法を使用すると仮定する。
式(６)において、測定値(チルダ)ｙからシーン反射率(ハット)ｓを再生するために、そのシーンがスパースであり、ＣＳ技法を使用すると仮定する。
具体的には、一実施の形態は、スパース性制約付き最小化問題
を解くことができる。ただし、Ｋは、ｓの最大スパース性、すなわち、離散化されたシーン内の反射体の最大数である。一般に、その問題はＮＰ困難であるが、ｌ０ノルムをそのｌ１凸包に緩和することによって、又は貪欲法を用いることによって解くことができる。一実施の形態は、反復ハードしきい値法(ＩＨＴ)を用いることができ、それは、反復ｔにおけるスパース推定値
が
を用いて推定される反復法である。ただし、ηはステップサイズであり、
は、その引数のうちのＫ個の最も大きな大きさ成分のみを保存し、残りの成分を０に設定するハードしきい値演算子である。ＩＨＴは、代替案に比べて計算コストと再生能力との優れたバランスを与えるので、好ましい実施の形態である。また、ＩＨＴは、モデルに基づくＣＳを用いる信号モデルへのより高い適応性も可能にする。ＩＨＴは、各反復時にステップサイズ選択を適応させることによって、更に加速させることができる。他の実施の形態は、マッチング追跡(ＭＰ：matching pursuit)、直交マッチング追跡(ＯＭＰ：orthogonal matching pursuit)、部分空間追跡(ＳＰ：subspace pursuit)、圧縮サンプリングマッチング追跡(ＣｏＳａＭＰ：Compressive Sampling Matching Pursuit)及び近似メッセージ伝搬(ＡＭＰ：approximate message passing)のような他の方法を用いて、式(１２)を解くことができる。例えば、米国特許第７，８３４，７９５号を参照されたい。
別の実施の形態は、
を近似しようと試みる凸最適化手法を用いることができる。
上記の実施の形態の場合と同様に、その問題はＮＰ困難であり、ｌ０ノルムをそのｌ１凸包に緩和することによって解くことができる。(１４)を近似するか、又は解く方法は、数ある中でも、以下の定式化のうちの１つを用いることができる。
これらの定式化は、数ある中でも、反復ソフトしきい値法アルゴリズム(ＩＳＴＡ：iterative soft thresholding algorithm)、固定点連続法(ＦＰＣ：fixed point continuation)、スパース再構成のための勾配投影法(ＧＰＳＲ：gradient projection for sparse reconstruction)、平滑化近接勾配法(ＳＰＧ：smoothing proximal gradient)のような幾つかの方法を用いて解くことができる。米国特許第７，８３４，７９５号を参照されたい。
アレイ設計の評価
この発明のアンテナ設計の実験的な評価から、以下の結論が得られる。互いに素なアレイは、より大きなＭＩＭＯ実質開口長を与えるので、入れ子アレイより良好なクロスレンジ分解能(ＭＬＡによって測定される)を有する。一方、入れ子アレイは、より低いＭＳＬを示す。ランダムアレイは一般に、互いに素なアレイ及び入れ子アレイに比べて高いＭＳＬを生成する。アレイ形状がＭＳＬに及ぼす影響は、壁が高い比誘電率を有するほど大きく減少する。全体として、ＭＳＬは、比誘電率が増加するにつれて増加する。壁の誘電率が低いほど、複数の反射によって、レンジプロファイルが曖昧になり、結果としてＭＬＡが拡大する可能性がある。壁の誘電率が高いほど、複数の反射によって、より良好なレンジ分解性がもたらされ、ＭＬＡは小さくなるが、ＭＳＬは大きくなる。ｖ)エンドファイアアレイに近い対象点は、空気−壁界面のフレネル反射係数が増加するので、複数の壁から、より深刻な影響を受ける。
この発明のアンテナ設計の実験的な評価から、以下の結論が得られる。互いに素なアレイは、より大きなＭＩＭＯ実質開口長を与えるので、入れ子アレイより良好なクロスレンジ分解能(ＭＬＡによって測定される)を有する。一方、入れ子アレイは、より低いＭＳＬを示す。ランダムアレイは一般に、互いに素なアレイ及び入れ子アレイに比べて高いＭＳＬを生成する。アレイ形状がＭＳＬに及ぼす影響は、壁が高い比誘電率を有するほど大きく減少する。全体として、ＭＳＬは、比誘電率が増加するにつれて増加する。壁の誘電率が低いほど、複数の反射によって、レンジプロファイルが曖昧になり、結果としてＭＬＡが拡大する可能性がある。壁の誘電率が高いほど、複数の反射によって、より良好なレンジ分解性がもたらされ、ＭＬＡは小さくなるが、ＭＳＬは大きくなる。ｖ)エンドファイアアレイに近い対象点は、空気−壁界面のフレネル反射係数が増加するので、複数の壁から、より深刻な影響を受ける。
この発明の方法及びシステムは、多くの種類の分野に適用可能である。
Claims (8)
- 壁背後のシーンを再構成するためのシステムであって、
前記壁を通して前記シーンの中に信号を送信するように構成される、送信信号を送信し前記送信信号の反射信号を受信信号として受信する多入力多出力(ＭＩＭＯ)アンテナアレイと、
前記受信信号から前記壁のパラメーターを推定し、前記パラメーターを用いて前記壁の誘電率のモデルを生成し、前記モデル、スパース再生及び圧縮検知を用いて、前記受信信号から前記シーンを画像として再構成するように構成されるプロセッサと、
を備え、
前記アンテナアレイのアンテナ素子の平均素子間間隔は不均一であり、前記送信信号の半波長よりはるかに大きい、システム。 - 壁背後のシーンを再構成するためのシステムであって、
前記壁を通して前記シーンの中に信号を送信するように構成される、送信信号を送信し前記送信信号の反射信号を受信信号として受信する多入力多出力(ＭＩＭＯ)アンテナアレイと、
前記受信信号から前記壁のパラメーターを推定し、前記パラメーターを用いて前記壁の誘電率のモデルを生成し、前記モデル、スパース再生及び圧縮検知を用いて、前記受信信号から前記シーンを画像として再構成するように構成されるプロセッサと、
を備え、
前記アンテナアレイはランダムアレイである、システム。 - 壁背後のシーンを再構成するためのシステムであって、
前記壁を通して前記シーンの中に信号を送信するように構成される、送信信号を送信し前記送信信号の反射信号を受信信号として受信する多入力多出力(ＭＩＭＯ)アンテナアレイと、
前記受信信号から前記壁のパラメーターを推定し、前記パラメーターを用いて前記壁の誘電率のモデルを生成し、前記モデル、スパース再生及び圧縮検知を用いて、前記受信信号から前記シーンを画像として再構成するように構成されるプロセッサと、
を備え、
前記パラメーターは、前記受信信号の被測定到達時刻(ＴＯＡ)と予測されるＴＯＡとの間の平均二乗誤差を最小化することによって得られる、システム。 - 前記スパース再生は、貪欲スパース再生法を用いる、請求項１から３までのいずれか１項に記載のシステム。
- 前記貪欲スパース再生法は、反復ハードしきい値法アルゴリズム(ＩＨＴ)を用いる、請求項４に記載のシステム。
- 前記貪欲スパース再生法は、各反復時にステップサイズ選択を適応させることによって加速される、請求項４に記載のシステム。
- 前記スパース再生は凸スパース近似法を用いる、請求項１から３までのいずれか１項に記載のシステム。
- 前記凸スパース近似法は、各反復時にステップサイズ選択を適応させることによって加速される、請求項７に記載のシステム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US14/153,686 | 2014-01-13 | ||
| US14/153,686 US9261592B2 (en) | 2014-01-13 | 2014-01-13 | Method and system for through-the-wall imaging using compressive sensing and MIMO antenna arrays |
| PCT/JP2014/084734 WO2015105033A1 (en) | 2014-01-13 | 2014-12-19 | Method and system for reconstructing scene behind wall |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2016535243A JP2016535243A (ja) | 2016-11-10 |
| JP6132978B2 true JP6132978B2 (ja) | 2017-05-24 |
Family
ID=52394288
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016513154A Active JP6132978B2 (ja) | 2014-01-13 | 2014-12-19 | 壁背後のシーンを再構成するためのシステム |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US9261592B2 (ja) |
| EP (1) | EP3094987B1 (ja) |
| JP (1) | JP6132978B2 (ja) |
| CN (1) | CN105917249B (ja) |
| IL (1) | IL245818B (ja) |
| WO (1) | WO2015105033A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7022628B2 (ja) | 2017-03-31 | 2022-02-18 | 株式会社Soken | 内燃機関用のスパークプラグ |
Families Citing this family (47)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9948362B2 (en) * | 2015-01-12 | 2018-04-17 | Mitsubishi Electric Research Laboratories, Inc. | System and method for 3D imaging using a moving multiple-input multiple-output (MIMO) linear antenna array |
| WO2016153914A1 (en) | 2015-03-25 | 2016-09-29 | King Abdulaziz City Of Science And Technology | Apparatus and methods for synthetic aperture radar with digital beamforming |
| US10545107B2 (en) * | 2015-04-26 | 2020-01-28 | Vayyar Imaging Ltd | System, device and methods for measuring substances' dielectric properties using microwave sensors |
| CN108432049B (zh) | 2015-06-16 | 2020-12-29 | 阿卜杜拉阿齐兹国王科技城 | 有效平面相控阵列天线组件 |
| US9915730B2 (en) * | 2015-07-23 | 2018-03-13 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for through-the-wall radar imaging |
| US10290948B2 (en) | 2015-08-02 | 2019-05-14 | Vayyar Imaging Ltd | System and method for radio frequency penetration imaging of an object |
| JP6584369B2 (ja) * | 2015-08-20 | 2019-10-02 | 三菱電機株式会社 | 壁の背後に位置するシーンの雑音のない画像を求めるシステムおよび方法 |
| CN106546983B (zh) * | 2015-09-17 | 2021-11-12 | 松下电器产业株式会社 | 雷达装置 |
| US10955546B2 (en) | 2015-11-25 | 2021-03-23 | Urthecast Corp. | Synthetic aperture radar imaging apparatus and methods |
| CN105652273B (zh) * | 2016-03-17 | 2018-10-26 | 哈尔滨工程大学 | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 |
| GB2548875A (en) * | 2016-03-31 | 2017-10-04 | Plextek Services Ltd | Signal detector and method |
| US11456532B2 (en) * | 2016-05-04 | 2022-09-27 | California Institute Of Technology | Modular optical phased array |
| US10382140B2 (en) | 2016-06-07 | 2019-08-13 | California Institute Of Technology | Optical sparse phased array receiver |
| US11249369B2 (en) | 2016-10-07 | 2022-02-15 | California Institute Of Technology | Integrated optical phased arrays with optically enhanced elements |
| US10795188B2 (en) | 2016-10-07 | 2020-10-06 | California Institute Of Technology | Thermally enhanced fast optical phase shifter |
| US11879989B2 (en) | 2016-12-05 | 2024-01-23 | Echodyne Corp. | Antenna subsystem with analog beam-steering transmit array and sparse hybrid analog and digital beam-steering receive array |
| EP3545334A1 (en) * | 2016-12-05 | 2019-10-02 | Echodyne Corp | Antenna subsystem with analog beam-steering transmit array and digital beam-forming receive array |
| WO2018148758A1 (en) | 2017-02-13 | 2018-08-16 | California Institute Of Technology | Passive matrix addressing of optical phased arrays |
| CN106972882B (zh) * | 2017-03-01 | 2020-08-07 | 浙江大学 | 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法 |
| WO2018165633A1 (en) * | 2017-03-09 | 2018-09-13 | California Institute Of Technology | Co-prime optical transceiver array |
| EP3593408A4 (en) * | 2017-03-09 | 2020-12-23 | California Institute of Technology | ARRANGEMENT OF OPTICAL CO-PRIME TRANSCEIVER |
| CA3064586A1 (en) | 2017-05-23 | 2018-11-29 | King Abdullah City Of Science And Technology | Synthetic aperture radar imaging apparatus and methods for moving targets |
| US11506778B2 (en) | 2017-05-23 | 2022-11-22 | Spacealpha Insights Corp. | Synthetic aperture radar imaging apparatus and methods |
| CN107300694B (zh) * | 2017-06-14 | 2020-04-17 | 南京信息工程大学 | 一种基于电磁波透射系数的未知墙体参数估计方法 |
| CN107589407B (zh) * | 2017-08-29 | 2020-06-26 | 电子科技大学 | 多发多收穿墙雷达成像后前后墙振荡一次多径抑制方法 |
| CN107894591A (zh) * | 2017-09-30 | 2018-04-10 | 沈阳航空航天大学 | 基于压缩感知的穿墙雷达衍射层析成像方法 |
| US11525910B2 (en) | 2017-11-22 | 2022-12-13 | Spacealpha Insights Corp. | Synthetic aperture radar apparatus and methods |
| CN107957574B (zh) * | 2017-12-28 | 2021-08-20 | 桂林电子科技大学 | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 |
| CN108562897B (zh) * | 2018-01-26 | 2022-01-11 | 桂林电子科技大学 | 一种mimo穿墙雷达的结构稀疏成像方法和装置 |
| CN108717188B (zh) * | 2018-05-29 | 2022-02-08 | 成都理工大学 | 适用于mimo穿墙雷达成像的快速墙体补偿方法 |
| CN108717189B (zh) * | 2018-06-01 | 2022-07-15 | 西安电子科技大学 | 基于压缩感知理论的双基地mimo雷达成像方法 |
| CN108614265A (zh) * | 2018-06-01 | 2018-10-02 | 金陵科技学院 | 基于NSL0重构算法的THz-FMCW雷达成像算法 |
| IT201800006797A1 (it) * | 2018-06-29 | 2019-12-29 | Radar ad apertura sintetica basato a terra (gbsar) con antenne multiple in trasmissione e ricezione (mimo)e facente uso della tecnica di elaborazione detta compressive sensing (cs). | |
| CN109541585B (zh) * | 2018-10-09 | 2020-11-06 | 中国人民解放军第四军医大学 | 一种基于峰度评估的人体穿墙检测成像方法 |
| CN109799499B (zh) * | 2019-01-28 | 2023-04-28 | 沈阳航空航天大学 | 一种穿墙雷达墙体参数估计方法 |
| CN109917361B (zh) * | 2019-04-02 | 2023-04-25 | 电子科技大学 | 一种基于双基雷达的三维未知场景成像方法 |
| CN110579737B (zh) * | 2019-07-17 | 2022-10-11 | 电子科技大学 | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 |
| CN111796249B (zh) * | 2020-02-28 | 2023-04-21 | 北京理工大学 | 一种基于墙体前表面能量对消的墙体参数估计方法 |
| CN111551928B (zh) * | 2020-06-02 | 2023-09-15 | 西安石油大学 | 一种基于墙体低秩稀疏约束的穿墙雷达成像方法 |
| CN111812644B (zh) * | 2020-08-29 | 2023-05-30 | 西安电子科技大学 | 基于稀疏估计的mimo雷达成像方法 |
| WO2022130561A1 (ja) * | 2020-12-17 | 2022-06-23 | 日本電信電話株式会社 | ミリ波推定装置、方法及びプログラム |
| US11906651B2 (en) * | 2021-02-25 | 2024-02-20 | Nxp B.V. | Radar-based detection using sparse array processing |
| US11927664B2 (en) | 2021-02-25 | 2024-03-12 | Nxp B.V. | Radar-based detection using angle of arrival estimation based on sparse array processing |
| CN112904271B (zh) * | 2021-03-03 | 2023-08-04 | 西北大学 | 基于互质阵列和增广扩展阵列的四阶累积量doa估计方法 |
| CN113238221B (zh) * | 2021-05-11 | 2023-09-22 | 沈阳航空航天大学 | 基于二维最小相位相干因子的mimo穿墙雷达成像方法 |
| CN113589250B (zh) * | 2021-07-29 | 2024-02-23 | 上海无线电设备研究所 | 一种稀疏学习的相参捷变频雷达距离高分辨处理方法 |
| CN114280544B (zh) * | 2021-12-02 | 2023-06-27 | 电子科技大学 | 一种基于松弛优化的最小过渡带宽度方向图赋形方法 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6177903B1 (en) * | 1999-06-14 | 2001-01-23 | Time Domain Corporation | System and method for intrusion detection using a time domain radar array |
| US7283085B2 (en) * | 2005-03-24 | 2007-10-16 | Agilent Technologies, Inc. | System and method for efficient, high-resolution microwave imaging using complementary transmit and receive beam patterns |
| US7280068B2 (en) * | 2005-07-14 | 2007-10-09 | Agilent Technologies, Inc. | System and method for microwave imaging with suppressed sidelobes using a sparse antenna array |
| US20070139248A1 (en) * | 2005-12-16 | 2007-06-21 | Izhak Baharav | System and method for standoff microwave imaging |
| US20080079625A1 (en) * | 2006-10-03 | 2008-04-03 | William Weems | System and method for stereoscopic anomaly detection using microwave imaging |
| US7504993B2 (en) * | 2006-10-12 | 2009-03-17 | Agilent Technolgoies, Inc. | Coaxial bi-modal imaging system for combined microwave and optical imaging |
| IL186884A (en) * | 2007-10-24 | 2014-04-30 | Elta Systems Ltd | Object simulation system and method |
| WO2010051108A2 (en) | 2008-09-11 | 2010-05-06 | Lawrence Livermore National Security, Llc | Model-based tomographic reconstruction |
| US7834795B1 (en) | 2009-05-28 | 2010-11-16 | Bae Systems Information And Electronic Systems Integration Inc. | Compressive sensor array system and method |
| US8077072B2 (en) * | 2009-09-03 | 2011-12-13 | Tialinx, Inc. | Static RF imaging for inside walls of a premises |
| US9239387B2 (en) | 2010-09-30 | 2016-01-19 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for reconstructing scenes using virtual arrays of transducers and joint sparsity models |
| US20120235849A1 (en) | 2010-11-19 | 2012-09-20 | Tatoian James Z | Method and apparatus for through the wall radar imaging |
| US9329263B2 (en) * | 2011-05-23 | 2016-05-03 | The Regents Of The University Of Michigan | Imaging system and method |
-
2014
- 2014-01-13 US US14/153,686 patent/US9261592B2/en active Active
- 2014-12-19 EP EP14828557.0A patent/EP3094987B1/en active Active
- 2014-12-19 WO PCT/JP2014/084734 patent/WO2015105033A1/en active Application Filing
- 2014-12-19 JP JP2016513154A patent/JP6132978B2/ja active Active
- 2014-12-19 CN CN201480072967.9A patent/CN105917249B/zh active Active
-
2016
- 2016-05-24 IL IL245818A patent/IL245818B/en active IP Right Grant
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7022628B2 (ja) | 2017-03-31 | 2022-02-18 | 株式会社Soken | 内燃機関用のスパークプラグ |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105917249B (zh) | 2018-10-30 |
| IL245818B (en) | 2021-03-25 |
| IL245818A0 (en) | 2016-07-31 |
| WO2015105033A1 (en) | 2015-07-16 |
| CN105917249A (zh) | 2016-08-31 |
| EP3094987B1 (en) | 2018-12-19 |
| US9261592B2 (en) | 2016-02-16 |
| US20150198713A1 (en) | 2015-07-16 |
| EP3094987A1 (en) | 2016-11-23 |
| JP2016535243A (ja) | 2016-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6132978B2 (ja) | 壁背後のシーンを再構成するためのシステム | |
| Leigsnering et al. | Compressive sensing in through-the-wall radar imaging | |
| JP6095010B2 (ja) | 信号を処理する方法 | |
| US10042046B2 (en) | System and method for radar imaging using distributed arrays and compressive sensing | |
| Ender | A brief review of compressive sensing applied to radar | |
| Zhang et al. | Ultrawideband impulse radar through‐the‐wall imaging with compressive sensing | |
| US9948362B2 (en) | System and method for 3D imaging using a moving multiple-input multiple-output (MIMO) linear antenna array | |
| US9335408B2 (en) | Method and system for through-the-wall imaging using sparse inversion for blind multi-path elimination | |
| Mansour et al. | Sparse blind deconvolution for distributed radar autofocus imaging | |
| CN106680809B (zh) | 一种穿墙雷达自聚焦压缩感知成像方法 | |
| Mansour et al. | Blind multi-path elimination by sparse inversion in through-the-wall-imaging | |
| CN108562901B (zh) | 基于最大信杂噪比准则的isar高分辨成像方法 | |
| Li et al. | Sparse MIMO architectures for through-the-wall imaging | |
| Liu et al. | Sparsity-driven distributed array imaging | |
| Lodhi et al. | Coherent radar imaging using unsynchronized distributed antennas | |
| Subotic et al. | Distributed RADAR waveform design based on compressive sensing considerations | |
| Leier et al. | Aperture undersampling using compressive sensing for synthetic aperture stripmap imaging | |
| Leigsnering et al. | CS based wall ringing and reverberation mitigation for through-the-wall radar imaging | |
| Bacci et al. | Time-slotted FMCW MIMO ISAR with compressive sensing image reconstruction | |
| Fromenteze et al. | Passive UWB multiplexing device for a single-port MIMO RADAR | |
| Yoon et al. | High resolution through-the-wall radar image based on beamspace eigenstructure subspace methods | |
| Zou et al. | Sparse Reconstruction Using Improved SL0 Algorithm by Tabu Search and Application to Spotlight SAR Imaging | |
| Leier et al. | Quality assessment of synthetic aperture sonar images based on a single ping reference | |
| Leigsnering et al. | Compressive sensing for urban multipath exploitation | |
| Qiao et al. | Azimuth ambiguity in redundant sampled stripmap SAS imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170124 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170302 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20170321 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170418 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6132978 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |