汪 玲 鄧 昊 王 潔

(南京航空航天大學雷達成像與微波光子技術教育部重點實驗室， 南京，211106)

適用于寬帶和窄帶照射源的無源合成孔徑雷達成像方法

汪 玲 鄧 昊 王 潔

(南京航空航天大學雷達成像與微波光子技術教育部重點實驗室， 南京，211106)

無源合成孔徑雷達成像利用外部照射源和運動接收機對感興趣場景進行成像，具有良好的電子對抗性，并可降低系統(tǒng)成本，利于系統(tǒng)小型化，緩解頻帶擁擠。本文給出基于逆散射理論和微局部分析的無源合成孔徑成像統(tǒng)一理論，適用于沿任意軌跡運動的機載接收機和靜止或運動的照射源。考慮發(fā)射兩種發(fā)射波形：窄帶連續(xù)波(Continuous-wave, CW)和寬帶脈沖波形，相應給出兩種新的無源合成孔徑成像模式：針對寬帶脈沖波形的合成孔徑無源成像(Synthetic aperture hitchhiker, SAH)和針對單頻或超窄帶連續(xù)波的多普勒合成孔徑無源成像(Doppler synthetic aperture hitchhiker, DSAH)。首先建立回波相關信號模型，該模型可消除與發(fā)射機相關的相位項,然后利用濾波反投影理論和微局部分析對相關回波數據進行反演，對應SAH和DSAH這兩種無源模式，分別將相關回波數據反投影到無源距離和無源多普勒等值線上,給出詳細的分辨率分析。最后通過仿真實驗對成像方法進行驗證。

合成孔徑雷達;無源雷達;成像;機會照射源;窄帶

引 言

隨著廣播站、移動通訊基站和導航衛(wèi)星等機會照射源數目的迅速增長，以及低成本接收機可快速部署的應用優(yōu)勢，近年來無源成像成為雷達領域的研究熱點之一[1-25]?，F有的大多數無源探測方法都集中在固定接收機配置下的目標檢測和空中目標成像[1-18]，較少涉及地面場景的無源合成孔徑成像。而且，大部分無源雷達需要設置一個輔助接收通道接收來自照射源的直達波信號，從中獲取發(fā)射波形，與回波信號相關完成匹配濾波，這要求接收天線具有良好的方向性和預知照射源位置信息。這種接收信號處理無疑增加了雷達系統(tǒng)的復雜性。研究無需外部照射源信息的無源合成孔徑雷達成像方法將使無源探測擺脫對發(fā)射源信息的依賴，并解決地面大場景的無源成像，可更好地發(fā)揮無源探測模式的優(yōu)勢。

文獻[19-25]針對無源合成孔徑成像相繼開展了研究，給出了照射源信息已知情況下的無源成像方法[19,20]，并對基于回波相關模型、無需照射源信息的無源成像方法進行了初步研究[21-25]。本文在逆散射理論和微局部分析的基礎上給出無源合成孔徑成像統(tǒng)一理論，并進行分辨率分析。該成像理論將基于反投影的圖像重建方法與統(tǒng)計波束形成方法、模糊理論聯(lián)系起來[26-28]。與現有的無源雷達檢測方法相比，該無源成像理論的優(yōu)勢主要體現在[21,22]：(1)可以工作在發(fā)射源波形、位置信息完全未知的情況下，即適用于非合作發(fā)射源，不再依賴參考通道獲取發(fā)射波形信息；(2)適用于靜止或運動的機會照射源；(3)適用于一個或多個機載接收機；(4)適用于任意載機飛行軌跡。該無源合成孔徑成像理論也適用于聲學、地球物理學和醫(yī)學領域的無源合成孔徑成像。

文中假設場景中多架接收機沿著任意航跡運動，并考慮兩類照射源波形：窄帶或超窄帶連續(xù)波和寬帶脈沖波形?；谡瓗нB續(xù)波(Continuous-wave,CW)的高多普勒分辨率特性，利用多普勒分辨進行成像的無源SAR稱為多普勒合成孔徑無源成像(Doppler synthetic aperture hitchhiker, DSAH)[21]，而利用寬帶波形的高距離分辨率進行成像的無源SAR稱為合成孔徑無源成像(Synthetic aperture hitchhiker, SAH)[22]。首先將不同接收機或接收機不同位置的接收回波信號進行相關，構建成像數據模型，該相關回波模型可去除發(fā)射機相關的相位項，因此成像不需要獲得發(fā)射機的位置。DSAH模型不依據“停-走-?！苯?，而是基于快時間多普勒效應(Temporal Doppler)或脈內多普勒效應，而SAH模型依據“停-走-停”近似，本質上仍然基于空間多普勒(Spatial Doppler)或慢時間多普勒。忽略快時間多普勒，DSAH模型退變?yōu)镾AH模型。對DSAH和SAH成像數據模型的高頻分析表明，相關回波模型是將場景輻射率分別投影到無源等多普勒線和無源等距離線上。采用濾波反投影方法對相關信號進行成像處理，實現場景輻射率的重建，成像分辨率的分析結果與模糊理論相一致[26,27]。仿真成像驗證了DSAH和SAH兩種無源SAR成像的可行性。

1 回波模型

假設接收機對和發(fā)射機的位置分別為T和R，則接收信號可以表示為[29]

(1)

(2)

圖1 無源SAR成像場景示意圖 Fig．1 Illustration of the imaging geometry

(3)

(4)

將接收信號ei和ej進行相關,有

(5)

假設發(fā)射源為非合作，即發(fā)射機的位置T和Jtr未知，并假設Jtr和場景散射率函數ρ統(tǒng)計獨立，則相關回波的期望E[cij]可以表示為

(6)

式中Cρ和CJtr分別為Jtr和ρ的相關函數，即

(7)

(8)

ARij為接收機天線波束圖的乘積

(9)

Gij為

(10)

對Jtr和ρ進行非相干近似[30]，使Jtr和ρ滿足以下等式

(11)

(12)

式中：Rρ是場景中z處的電磁輻射平均功率[30]；RT是接收機位置T處發(fā)射的電磁輻射到目標點z處的平均功率。因此，Rρ其實就是場景輻射率，而RT為發(fā)射機的發(fā)射功率。

將式(11)和式(12)代入式(6)，得到

(13)

1.1DSAH成像模型

(14)

式中

(15)

(16)

(17)

(18)

(19)

當接收機的速度遠小于光速時，將式(20)括號里的等式左右兩邊同時乘以ω0，并進行整理后得到

(20)

圖2 兩個接收平臺沿同一圓形軌跡飛行，某一時刻對應的無源多普勒等值線Fig．2 DSAH iso-Doppler contours for two receivers traversing a circular flight trajectory

1.2SAH成像模型

(21)

為利用SAH發(fā)射源波形的高距離分辨率，按照式(22)構建SAH的前向成像數據模型,有

(22)

式(22)對應ei和ej在快時間和慢時間兩維時間域上進行相關，也稱為空時相關，其中τ對應快時間延遲。在SAH成像中，s′，s的采樣率通常與脈沖重復頻率一致。

(23)

式中

(24)

(25)

式(24)中

(26)

(27)

定義式(27)對應的曲線Hij(s′,s,τ)為無源距離等值線。圖3給出了平坦地形假設下，接收機繞圓形軌跡飛行某一孔徑采樣時刻的無源距離等值線圖。黑色和白色箭頭表示接收載機某個采樣時刻的位置。

圖3 兩個接收平臺沿同一圓形軌跡飛行，某一時刻對應的無源距離等值線Fig．3 SAH iso-range contours for two receivers traversing a circular flight trajectory

2 成像方法

2.1DSAH濾波反投影算法

DSAH成像可描述為[21]

(28)

將式(14)代入式(29)，經過化簡可得

(30)

(31)

式中Dψ(z)為地面投影算子，完成空間任一三維矢量向地表的投影，其形式為

(32)

將式(39)和式(40)代入式(37)，得到

(33)

引入如下變量變換，有

(34)

式(33)轉換為

(35)

(36)

選擇如下的濾波器使PSF盡可能接近Dirac-delta函數，即有

(37)

2.2SAH濾波反投影算法

SAH成像可描述為[22]

(38)

(39)

將式(23)代入式(39)，利用駐留相位定理，近似得到

(40)

(41)

將式(41)代入式(40)，并進行如下的變量代換

(42)

(43)

(44)

選擇如下的濾波器使PSF盡可能接近Dirac-delta函數,有

(45)

通過式(45)選擇合適的濾波器，可以在正確的位置和方向重建場景輻射率的可見邊緣，而且確保幅度不失真。

3 分辨率分析

將式(37)和式(45)分別代入式(35)和式(43)，可以得到

(46)

(47)

(48)

(49)

式中：Lφ為窗函數長度或一個孔徑采樣時刻參與成像處理的信號長度，Bω為發(fā)射信號帶寬。

另外，在DSAH或SAH成像中，孔徑采樣樣本和成像使用時間窗長度的增加，也會增大數據采集空間，提高成像分辨率。

4 仿真驗證

仿真場景為[0,22] km×[0,22] km的平坦地勢，離散成128×128個像素點，[0,0,0] km和[22,22,0]km分別對應像素點(1,1)和(128,128)。假設發(fā)射機和接收機的天線均為各向同性。仿真使用兩部機載接收機和單個固定的發(fā)射機，發(fā)射機的位置為y0=(0,0,6.5) km。兩接收機沿同一圓軌跡飛行，圓軌跡方程為

(50)

根據非相干近似，采用以下多點目標模型

(51)

(52)

圖4 DSAH和SAH成像中所用場景Fig．4 Scenes used in DSAH and SAH simulations

圖5 成像幾何關系示意圖 Fig．5 Illustration of the imaging geometry

4.1 DSAH成像仿真驗證

根據式(5)生成數據，并選擇式(5)中的窗函數為漢寧窗，長度為0.085 3 s,使用載頻為800 MHz的單頻CW波。對于合作發(fā)射機，重建圖像如圖6(a)所示，可以看出，使用DSAH法可以很好地重建圖像。圖6(b)為使用非合作發(fā)射機得到的重建圖像。由于發(fā)射機的位置未知，所接收的信號未補償與發(fā)射機相關的幾何擴展因子,因此接近發(fā)射機的目標比那些遠離發(fā)射機的更亮。

圖6 DSAH 成像結果Fig．6 DSAH imaging results

4.2 SAH成像仿真驗證

根據式(22)生成數據，和DSAH仿真一樣選擇式(22)中的窗函數為漢寧窗。載頻為1 GHz、帶寬為30 MHz的發(fā)射波形。

圖7(a)和圖7(b)分別給出了合作和非合作發(fā)射機情況下的重建圖像。在這兩種情況下，圖像均得到很好的重建，且接近發(fā)射機的目標更亮。

圖7 SAH 成像結果Fig．7 SAH imaging results

5 結束語

本文給出了無源合成孔徑成像統(tǒng)一理論。首先建立了回波相關信號模型，對不同位置接收機接收的加窗、尺度變化和平移后的回波信號進行相關。相關模型分為DSAH和SAH兩種，在“停-走-?！奔僭O下，DSAH模型轉化為SAH模型。相關過程從DSAH和SAH相關模型中去除了與發(fā)射源信息相關的相位項，無需發(fā)射機的位置信息即可進行反投影。采用濾波反投影進行成像，微局部分析表明DSAH和SAH成像分別將相關信號反投影到無源多普勒等值線和無源距離等值線上。分辨率分析表明，DSAH重建圖像的分辨率主要由參與成像處理的信號長度或窗函數長度以及發(fā)射源載頻決定，SAH重建圖像的分辨率主要由發(fā)射波形帶寬決定。這個結論與CW或超窄帶波形、寬帶波形的模糊理論相一致。仿真成像驗證了DSAH和SAH這兩種無源SAR成像方法。

[1] Baker C J, Griffiths H D, Papoutsis I. Passive coherent location radar systems. Part 2: Waveform properties[J]. IEE Proceedings Radar Sonar and Navigation, 2005, 152(3):160-168.

[2] Chetty K, Woodbridge K, Guo H, et al. Passive bistatic WiMAX radar for marine surveillance[C]// IEEE Radar Conference. Washington, D C: IEEE, 2010:188-193.

[3] Christiansen J M, Olsen K E. Range and Doppler walk in DVB-T based passive bistatic radar[C]// IEEE Radar Conference. Washington, D C: IEEE, 2010:620-626.

[4] Coleman C, Yardley H. Passive bistatic radar based on target illuminations by digital audio broadcasting[J]. IET Radar Sonar & Navigation, 2008, 2(5):366-375.

[5] Falcone P, Colone F, Bongioanni C, et al. Experimental results for OFDM WiFi-based passive bistatic radar[C]// IEEE Radar Conference. Washington, D C: IEEE, 2010:516-521.

[6] Garnier J, Papanicolaou G. Passive imaging using cross correlations of ambient noise signals[C]// IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP). Aruba, Dutch Antilles: IEEE, 2009:221-224.

[7] Garnier J, Papanicolaou G. Passive sensor imaging using cross correlations of noisy signals in a scattering medium[J]. Siam Journal on Imaging Sciences, 2010, 2(2):396-437.

[8] Garnier J, Papanicolaou G. Resolution analysis for imaging with noise[J]. Inverse Problems, 2010, 26(26):074001.

[9] Garnier J, Solna K. Passive imaging and detection in cluttered media[C]// IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP). Aruba, Dutch Antilles: IEEE, 2009:225-228.

[10]Griffiths H D, Baker C J. Passive coherent location radar systems. Part 1: Performance prediction [J]. IEE Proceedings Radar, Sonar and Navigation, 2005, 152(3):153-159.

[11]Guo H, Woodbridge K, Baker C J. Evaluation of WiFi beacon transmissions for wireless based passive radar[C]// IEEE Radar Conference. Rome: IEEE, 2008:1-6.

[12]Harms H A, Davis L M, Palmer J. Understanding the signal structure in DVB-T signals for passive radar detection[C]// IEEE Radar Conference. Washington, D C: IEEE, 2010: 532-537.

[13]Homer J, Kubik K, Mojarrabi B, et al. Passive bistatic radar sensing with LEOS based transmitters[C]// IEEE International Symposium on Geoscience and Remote Sensing Symposium. Toronto, Ontario, Canada: IEEE, 2002: 438-440.

[14]Howland P E, Maksimiuk D, Reitsma G. FM radio based bistatic radar [J]. IEE Proceedings Radar Sonar and Navigation, 2005, 152(3):107-115.

[15]O'Hagan D W, Baker C J. Passive bistatic radar (PBR) using FM radio illuminators of opportunity[C]// New Trends for Environmental Monitoring Using Passive Systems. Hyeres, French Riviera: IEEE, 2008:1-6.

[16]Poullin D. Passive detection using digital broadcasters (DAB, DVB) with COFDM modulation [J]. IEE Proceedings Radar, Sonar and Navigation, 2005, 152(3):143-152.

[17]Tan D K P, Sun H, Lu Y, et al. Passive radar using global system for mobile communication signal: Theory, implementation and measurements [J]. IEE Proc Radar, Sonar and Navigation,2005, 152(3):116-123.

[18]Wang L, Son I Y, Yazici B. Passive imaging using distributed apertures in multiple-scattering environments [J]. Inverse Problems, 2010, 26(6):065002.

[19]Antoniou M, Cherniakov M, Hu C. Space-surface bistatic SAR image formation algorithms [J]. IEEE Transactions on Geoscience & Remote Sensing, 2009, 47(6):1827-1843.

[20]Li G, Xu J, Peng Y N, et al. Bistatic linear antenna array SAR for moving target detection, location, and imaging with two passive airborne radars[J]. IEEE Transactions on Geoscience & Remote Sensing, 2007, 45(3):554-565.

[21]Wang L, Yarman C E, Yazc B. Doppler-hitchhiker: A novel passive synthetic aperture radar using ultra-narrowband sources of opportunity[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10):3521-3537.

[22]Yarman C E, Yazici B. Synthetic aperture hitchhiker imaging [J]. IEEE Trans on Image Processing, 2008, 17(11):2156-2173.

[23]汪玲，莊旭昇，龐鑠. 適用于任意接收機航跡的無源WiMAX-SAR成像方法[J].數據采集與處理，2013,28(3):404-409.

Wang Ling, Zhuang Xusheng, Pang Shuo. Passive WiMAX-SAR imaging method applicable to arbitrary flight trajectories of receiver[J]. Journal of Data Acquisition and Processing, 2013,28(3):404-409.

[24]張璇，汪玲. 一種基于回波相關的無源合成孔徑雷達成像方法[J].電子與信息學報，2012,34(6):1511-1515.

Zhang Xuan, Wang Ling. An echo-correlation based passive imaging method for synthetic aperture radar [J]. Journal of Electronics & Information Technology, 2012,34(6):1511-1515.

[25]張璇，汪玲. 適用于非合作照射源的合成孔徑雷達無源成像方法[J].航空學報，2012,34(6):1397-1404.

Zhang Xuan, Wang Ling. A passive synthetic aperture radar imaging method for non-cooperative transmitters [J]. Acta Aeronautica et Astronautica Sinica, 2012,34(6):1397-1404.

[26]Levanon N, Mozeson E. Radar signals [M]. New Jersey:John Wiley & Sons, Inc,2004.

[27]斯科尼克.雷達系統(tǒng)導論[M].北京：電子工業(yè)出版社，2007.

Skolnik M I. Introduction to radar systems [M]. Beijing:Publishing House of Electronics Industry, 2007.

[28]Veen B D V, Buckley K M. Beamforming: A versatile approach to spatial filtering [J]. IEEE ASSP Magazine, 2002, 5(2):4-24.

[29]Nolan C J, Cheney M. Synthetic aperture inversion for arbitrary flight paths and non-flat topography [J]. IEEE Transactions on Image Processing Publication of the IEEE Signal Processing Society, 2003, 12(9):1035-1043.

[30]Barrett H H, Myers K J, Rathee S. Foundations of image science [J]. Journal of Electronic Imaging, 2005, 14(4):029901.

[31]Yazici B, Cheney M, Yarman C E. Synthetic-aperture inversion in the presence of noise and clutter [J]. Inverse Problems, 2006, 22(5):1705-1729.

[32]Nolan C J, Cheney M. Synthetic aperture inversion [J]. Inverse Problems, 2002, 18(18):221.

[33]Quinto E T. Singularities of the X-ray transform and limited data tomography in r2 and r3 [J]. Siam Journal on Mathematical Analysis, 1993, 24(5):1215-1225.

Passive Synthetic Aperture Radar Imaging Methods for Wideband and Narrowband Illuminators of Opportunity

Wang Ling, Deng Hao, Wang Jie

(Key Laboratory of Radar Imaging and Microwave Photonics, Ministry of Education, Nanjing University of Aeronautics and Astronautics, Nanjing, 211106, China)

Passive synthetic aperture radar imaging uses illuminators of opportunity and moving receivers to image the scene of interest and has good electronic antagonism. It can also reduce system cost, benefit to system miniaturization and ease frequency band congestion. In this paper, we present a unified theory for passive synthetic aperture radar imaging based on inverse scattering and microlocal analysis. It is suitable for airborne receivers flying along arbitrary flight trajectories and static or moving illuminators of opportunity. Two types of waveforms are considered: narrowband continuous-wave (CW) waveforms and wideband pulsed waveforms. Our theory results in two novel synthetic aperture imaging modalities: Synthetic aperture hitchhiker (SAH) that uses wideband pulsed waveforms and Doppler synthetic aperture Hitchhiker (DSAH) that uses single-frequency or ultra-narrowband CW waveforms. First, we develop measurement models which are able to remove the transmitter-related terms in the phase of the correlated measurements. We use filtered backprojection (FBP) theory and microlocal analysis to develop approximate inversion formulas for SAH and DSAH. The inversion formulas involve backprojection of the correlated measurements onto iso-passive-range and iso-passive-Doppler contours for SAH and DSAH imaging, respectively. Detailed resolution analysis is given. Finally, we present numerical simulations to demonstrate our theoretical results.

synthetic aperture radar; passive radar; imaging; illuminators of opportunity; narrow band

南京航空航天大學研究生創(chuàng)新基地(實驗室)開放基金(kfjj20160402)資助項目; 中央高?；究蒲袠I(yè)務費專項基金資助項目；航空科學基金(20132052035)資助項目。

2017-05-03；

2017-06-10

TN951

A

汪玲(1977-)，女，教授、博士生導師，研究方向:雷達成像、合成孔徑雷達(SAR)成像、逆合成孔徑雷達(ISAR)成像以及無源雷達成像，E-mail：tulip_wling@nuaa.edu.cn。

鄧昊(1992-)，女，碩士研究生，研究方向：通信與信息系統(tǒng)和合成孔徑雷達成像。

王潔(1993-)，女，碩士研究生，研究方向：電子與通信工程和無源雷達成像。