楊 科， 廖桂生， 徐 青

(西安電子科技大學(xué) 雷達(dá)信號處理國家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710071)

筆者根據(jù)ROSAR系統(tǒng)能夠快速重訪同一場景的特點(diǎn)，并結(jié)合直升機(jī)等旋翼飛行器易于懸停的特性，建立了ROSAR 的干涉成像模型．在該模型下，成像平臺首先保持在某一高度上完成場景的一次觀測，然后上升至另一高度完成場景的二次觀測，經(jīng)過兩次或者多次觀測可以獲得多組ROSAR數(shù)據(jù)，通過對這些數(shù)據(jù)進(jìn)行干涉處理獲得場景的高程信息．在實(shí)際的數(shù)據(jù)采集過程中，由于受強(qiáng)風(fēng)、對流等天氣影響，同一系統(tǒng)平臺在不同高度上存在旋轉(zhuǎn)軸偏移問題．軸偏移將會引起天線與場景之間的距離變化，而在ROSAR模式中，由于受到物理結(jié)構(gòu)的限制，不同斜距對應(yīng)的合成孔徑大小存在差異．因此，軸偏移會導(dǎo)致同一場景目標(biāo)的方位帶寬有所變化，同時軸偏移也會引起干涉相位發(fā)生變化，使得高程信息無法正確獲取．為了使不同高度的轉(zhuǎn)軸位置相互重合，同時保障轉(zhuǎn)軸位置和場景位置的相對關(guān)系，筆者采用轉(zhuǎn)軸和場景整體平移的方法，將軸偏移等效為場景位置移動，進(jìn)而分析場景位置變化對方位帶寬和干涉相位的影響，獲取斜距差和偏移量的關(guān)系式，最后通過構(gòu)建補(bǔ)償函數(shù)來消除軸偏移的影響．

圖1 ROSAR干涉模型示意圖

1 ROSAR干涉模型

如圖1所示，雷達(dá)平臺距地面的高度為H，平臺上的天線以角速度ω作勻速圓周運(yùn)動，并指向背離轉(zhuǎn)軸方向，對周圍場景進(jìn)行圓環(huán)形照射．定義天線到轉(zhuǎn)軸中心的垂直距離(旋轉(zhuǎn)半徑)為ra，俯仰角為θ，俯仰向波束寬度為ε，方位角為0°(天線照射方向垂直于天線速度方向)，方位向波束寬度為η．考慮到天線作圓周運(yùn)動，為了便于分析和計(jì)算，文中采用圓柱坐標(biāo)系(r，φ，z)．不失一般性，假設(shè)天線旋轉(zhuǎn)至x軸平行位置時，旋轉(zhuǎn)角定義為0°．旋轉(zhuǎn)至任意位置Q時，旋轉(zhuǎn)角為ωta，則天線位置為(ra，ωta，H)．假設(shè)點(diǎn)目標(biāo)P的位置為(r0，φ0，z0)，r0表示點(diǎn)目標(biāo)的地面距離(地距)，φ0表示點(diǎn)目標(biāo)的方位向位置，z0表示點(diǎn)目標(biāo)的高度，則天線至點(diǎn)目標(biāo)的斜距表達(dá)式為

(1)

假設(shè)雷達(dá)發(fā)射線性調(diào)頻脈沖信號p(t)=rect(t/Tp) exp(j2πfct+jπγt2)，Tp為發(fā)射脈沖寬度，fc為發(fā)射信號載頻，γ為信號調(diào)頻率．點(diǎn)目標(biāo)回波經(jīng)過基頻變換，在距離快時間-方位慢時間域 (tr-ta域)可表示為

(2)

其中，λ是發(fā)射信號波長．

利用已有的二維ROSAR成像算法[2，7-8]，對原始回波數(shù)據(jù)S0進(jìn)行二維匹配成像，可得單幅ROSAR圖像為

(3)

其中，A1表示二維成像后的幅度，Br表示發(fā)射信號帶寬，Ba表示方位向帶寬，α表示天線至點(diǎn)目標(biāo)的最短斜距，α= ((ra-r0)2+(H-z0)2)1/2．

當(dāng)平臺垂直上升至另一高度H+B⊥時，對場景進(jìn)行第2次觀測，如圖1所示，可利用二維成像算法獲得第二幅ROSAR圖像．通過對兩幅圖像進(jìn)行干涉處理可獲得場景的高程信息，即能實(shí)現(xiàn)ROSAR干涉成像的目的．

初中數(shù)學(xué)學(xué)科涉及的知識范圍廣，難度大，平時還有許多考試，需要時間進(jìn)行復(fù)習(xí)，因此，對教師的教學(xué)進(jìn)度有一定要求。這使得教師的教學(xué)過程往往以提高課堂教學(xué)效率為目標(biāo)，在培養(yǎng)學(xué)生獨(dú)立思考能力方面就會流于表面，只注重形式化的過場，而沒有把培養(yǎng)學(xué)生獨(dú)立思考能力的教學(xué)活動真正落到實(shí)處，學(xué)生無法理解獨(dú)立思考的本質(zhì)，也就無法提高自身的獨(dú)立思考能力。

2 軸偏移影響及補(bǔ)償方法

在實(shí)際情況中，由于受到周圍環(huán)境影響，不同高度的旋轉(zhuǎn)軸中心可能存在偏移，如圖1中的虛線所示，兩軸之間的偏移量為xd．為了方便討論，圖中假設(shè)轉(zhuǎn)軸偏移方向與x軸平行，這一假設(shè)不影響模型實(shí)質(zhì)描述．要使第2次數(shù)據(jù)錄取坐標(biāo)與第1次的相同，需對第2個轉(zhuǎn)軸中心進(jìn)行 (-xd，0，0)的位置校正．由于轉(zhuǎn)軸位置和場景位置是相對的，為了滿足轉(zhuǎn)軸與場景的位置相對不變，需同時對照射場景進(jìn)行整體移動．原場景內(nèi)的點(diǎn)目標(biāo)P位置為(r0，φ0，z0)，假設(shè)移動后的點(diǎn)目標(biāo)P′位置為(roffset，φoffset，z0)，P與P′的關(guān)系為

整理式(4)和式(5)可得新的位置參數(shù)為

在新的坐標(biāo)位置下，式(3)中的斜距位置α、方位向位置φ0以及相位大小 -4πα/λ發(fā)生改變，其中，新的斜距大小α′= ((ra-roffset)2+ (H-z0)2)1/2．同時，由于ROSAR的特殊幾何結(jié)構(gòu)，在新的坐標(biāo)中方位向帶寬Ba也有相應(yīng)的變化，導(dǎo)致干涉圖像之間的相關(guān)性下降．以下具體分析轉(zhuǎn)軸偏移對方位帶寬和干涉相位的影響．

2.1 對方位帶寬的影響

在傳統(tǒng)的線性SAR成像系統(tǒng)中，多普勒調(diào)頻率隨斜距的增加而減少，但是合成孔徑時間隨斜距的增加而增加，因此，不同斜距的方位帶寬是相同的[9-10]．而在ROSAR系統(tǒng)中，由于物理結(jié)構(gòu)的限制，合成孔徑時間不能隨斜距的增長而線性增長，進(jìn)而導(dǎo)致不同斜距的方位帶寬不同．以下將基于ROSAR的方位向調(diào)頻率和合成孔徑時間的計(jì)算，對方位帶寬進(jìn)行推導(dǎo)，進(jìn)而分析轉(zhuǎn)軸偏移對方位帶寬的影響．

利用匹配函數(shù)對原始回波數(shù)據(jù)S0進(jìn)行距離向匹配，可得

(8)

其中，A2表示距離向匹配后的幅度．

對式(8)中的相位進(jìn)行二次求導(dǎo)，可得方位向調(diào)頻率為

(9)

圖2 ROSAR合成孔徑示意圖

(10)

其中，η′為方位向波束寬度η在水平面上的投影．對于高度位置為z0的點(diǎn)目標(biāo)，假設(shè)俯仰角為θz0，根據(jù)幾何關(guān)系可得 cosθz0= (r0-ra)/α，從而可得η′和η的關(guān)系表達(dá)式為

tan(η′/2)=tan(η/2)/cosθz0．

(11)

文獻(xiàn)[11]給出了二維場景的合成孔徑長度表達(dá)式，在不考慮平臺高度H和點(diǎn)目標(biāo)高度z0的情況下，筆者將三維場景簡化為二維場景，令η′=η，式(10)簡化為S=ηra(1-ra/r0)，此時該式與文獻(xiàn)[11]中的表達(dá)形式相一致．

假設(shè)T為方位向有效合成孔徑時間，則時間T與孔徑S的關(guān)系可表示為

ωT=S/ra．

(12)

根據(jù)式(9)的方位向調(diào)頻率Ka及式(10)～式(12)所得的合成時間T，可得方位向帶寬為

Ba=|KaT|=2βωη′(1-ra/r0)/(λα) ．

(13)

由式(13)可知，方位帶寬Ba隨著r0的變化而變化，而轉(zhuǎn)軸偏移可以等效于r0的變化，進(jìn)而引起B(yǎng)a的變化．根據(jù)文獻(xiàn)[12-13]可知，方位帶寬的差異將會導(dǎo)致信號去相關(guān)，進(jìn)而使得干涉相位的等效噪聲增加．如果轉(zhuǎn)軸偏移引起的方位帶寬變化較大，在預(yù)處理的過程中需要對原圖像進(jìn)行方位預(yù)濾波，從而保證最終的干涉成像結(jié)果．

2.2 對干涉相位的影響及補(bǔ)償方法

轉(zhuǎn)軸偏移對圖像相位也會產(chǎn)生影響，即式(3)中相位 exp(-j4πα/λ) 變?yōu)?exp(-j4πα′/λ)．以下分析兩者相位差 exp(-j4πΔα/λ) 與轉(zhuǎn)軸偏移量xd的關(guān)系表達(dá)式，其中斜距差表示為 Δα=α′-α．

將變化后的地距roffset表示為r0+Δr0，并對斜距α′進(jìn)行泰勒級數(shù)展開，可得

(14)

式(14)中二次及二次以上高次項(xiàng)為極小量，在只保留一次項(xiàng)的情況下可得斜距差Δα= (r0-ra)Δr0/α．而地距差可表示為 Δr0=roffset-r0，重新改寫式(6)，并對其進(jìn)行泰勒級數(shù)展開，可得

(15)

忽略式(15)中的二次及二次以上極小量，可得Δr0=-cosφ0xd．結(jié)合式(14)所得結(jié)果，以Δr0為中間變量，可得斜距差Δα與轉(zhuǎn)軸偏移量xd的關(guān)系式為

Δα=[-(r0-ra)/α] cosφ0xd．

(16)

由以上分析可知，旋轉(zhuǎn)中心軸的偏移導(dǎo)致圖像增加了一個附加相位 exp(-j4πΔα/λ)，其中Δα如式(16)所示．通過干涉手段獲取場景的高程信息，對干涉相位的精確性有極高的要求．為了正確恢復(fù)場景高程信息，需要對上述附加相位進(jìn)行補(bǔ)償．補(bǔ)償函數(shù)為

(17)

對兩幅ROSAR圖像進(jìn)行配準(zhǔn)處理，并通過復(fù)共軛相乘獲取原始干涉相位．該相位包含了軸偏移所引起的附加相位．在去平地相位之前，利用式(17)所示函數(shù)與原始相位進(jìn)行相乘，可實(shí)現(xiàn)附加相位校正．接著通過去平地相位、相位解纏繞及高度換算等處理，最終可獲得ROSAR場景高程信息．

3 仿真實(shí)驗(yàn)

ROSAR系統(tǒng)的主要仿真參數(shù)如下: 平臺高度 1 000 m，旋轉(zhuǎn)半徑 1.5 m，地距變化范圍 1 100 m～ 2 300 m，俯仰角60°，信號發(fā)射波長 0.03 m．圖3給出了旋轉(zhuǎn)軸偏移對方位帶寬的影響．通常情況下旋轉(zhuǎn)軸受環(huán)境影響存在數(shù)米到數(shù)十米的偏移量，圖3考慮極大偏移情況，選取偏移量為 100 m．圖3(a)以場景目標(biāo)的方位位置φ0為參變量，給了軸偏移前后的方位帶寬差ΔBa隨地距r0的變化曲線．由圖可知，軸偏移會引起方位帶寬發(fā)生變化(即存在方位帶寬差)，軸偏移方向上 (φ0=0°，實(shí)線所示)的目標(biāo)帶寬差變化較大，遠(yuǎn)離軸偏移方向 (φ0=90°，虛線所示)的目標(biāo)帶寬差變化較小，但是對于各方位位置上的目標(biāo)來說，方位帶寬變化都為較小量．圖3(b)給出了某一點(diǎn)目標(biāo)在軸偏移前后的方位帶寬截面圖，點(diǎn)目標(biāo)的具體位置為方位向φ0=0°、地距向r0=1 733 m，由圖可知，兩者的方位帶寬基本上是一致的．綜上所述，軸偏移會導(dǎo)致方位帶寬發(fā)生變化，但是該變化量較小，可忽略其影響．

圖3 方位帶寬影響分析

圖4給出了旋轉(zhuǎn)軸偏移所引起的干涉相位影響．圖4(a)為不同平臺高度的兩幅ROSAR圖像干涉相位圖，兩幅圖像在產(chǎn)生干涉相位前已經(jīng)過圖像配準(zhǔn)處理，平臺1高度為 1 000 m，平臺2高度為 1 001 m，假設(shè)平臺1的軸位置為參考位置，平臺2相對平臺1水平偏移 2 m，其它參數(shù)如上述所示．平臺2的軸位置偏移導(dǎo)致干涉相位中存在一個附加相位，圖4(b)給出了上述偏移所產(chǎn)生的附加相位圖，該相位將會影響最終的高程信息獲取，為此需要對這一附加相位進(jìn)行補(bǔ)償，補(bǔ)償結(jié)果如4(c)所示．圖4(c)中包含平地相位，對其進(jìn)行去平地相位處理，并進(jìn)行相位解纏繞及高度換算，最終可得到如圖4(d)所示的場景高程圖．

圖4 干涉成像結(jié)果

4 結(jié) 束 語

筆者構(gòu)建了ROSAR干涉成像模型，分析了不同高度系統(tǒng)平臺的軸偏移問題．通過分析發(fā)現(xiàn)，軸偏移會影響方位帶寬和干涉相位．實(shí)驗(yàn)仿真表明，軸偏移導(dǎo)致ROSAR方位帶寬產(chǎn)生變化，但是該變化量較小，在干涉處理過程中可以忽略上述影響．軸偏移對干涉相位產(chǎn)生影響，使得圖像有一個附加相位，該相位影響較大，仿真實(shí)驗(yàn)在去平地相位之前對該相位進(jìn)行了校正，獲得較好的補(bǔ)償效果．最后，經(jīng)過解纏繞及高度換算，獲得場景的三維高程圖．

[1] Jeon M, Kim Y S. Migration Technique for Rotor Synthetic Aperture Radar [J]. Electronics Letters, 1997, 33(7): 630-631.

[2] Lee H, Cho S J, Kim K E. A Ground-based Arc-scanning Synthetic Aperture Radar (arc-SAR) System and Focusing Algorithms [C]//IEEE International Geoscience and Remote Sensing Symposium. Piscataway: IEEE, 2010: 3490-3493.

[3] Ali F, Urban A, Vossiek M. A Short Range Synthetic Aperture Imaging Radar with Rotating Antenna [J]. International Journal of Electronics and Telecommunications, 2011, 57(1): 97-102.

[4] Zhang Jun. The Static Small Object Detection Based on Ground-based Arc SAR [C]//International Conference on Microwave and Millimeter Wave Technology. Piscataway: IEEE, 2012: 1489-1492.

[5] Ghasr M T, Pommerenke D, Case J T, et al. Rapid Rotary Scanner and Portable Coherent Wideband Q-band Transceiver for High-resolution Millimeter-wave Imaging Applications [J]. IEEE Transactions on Instrumentation and Measurement, 2011, 60(1): 186-197.

[6] 廖軼, 周松, 邢孟道, 等. 一種基于級數(shù)反演的機(jī)載圓跡環(huán)掃SAR成像算法 [J]. 電子與信息學(xué)報(bào), 2012, 34(11): 2587-2593.

Liao Yi, Zhou Song, Xing Mengdao, et al. An Imaging Algorithm for Airborne Circular Scanning SAR Bbased on the Method of Series Reversion [J]. Journal of Electronics & Information Technology, 2012, 34(11): 2587-2593.

[7] Bara M, Sagués L, Paniagua F, et al. High-speed Focusing Algorithm for Circular Synthetic Aperture Radar (C-SAR) [J]. Electronics Letters, 2000, 36(9): 828-830.

[8] 孫兵, 周蔭清, 陳杰, 等. 廣域觀測圓軌跡環(huán)掃SAR成像模式研究 [J]. 電子與信息學(xué)報(bào), 2008, 30(12): 2805-2808.

Sun Bing, Zhou Yinqing, Chen Jie, et al. Operation Mode of Circular Trace Scanning SAR for Wide Observation [J]. Journal of Electronics & Information Technology, 2008, 30(12): 2805-2808.

[9] 侯育星, 張磊, 徐剛, 等. 脈沖相位編碼信號大測繪帶合成孔徑雷達(dá)成像 [J]. 西安電子科技大學(xué)學(xué)報(bào), 2012, 39(5): 54-60.

Hou Yuxing, Zhang Lei, Xu Gang, et al. Wide Swath SAR Imaging with Pulse-phase-codes [J]. Journal of Xidian University, 2012, 39(5): 54-60.

[10] 梁毅, 王虹現(xiàn), 邢孟道, 等. 同航線雙基調(diào)頻連續(xù)波SAR改進(jìn)距離徙動算法 [J]. 西安電子科技大學(xué)學(xué)報(bào), 2011, 38(1): 71-79.

Liang Yi, Wang Hongxian, Xing Mengdao, et al. Modified RMA for the Tandem Bistatic FMCW SAR with Frequency Non-linearity Correction [J]. Journal of Xidian University, 2011, 38(1): 71-79.

[11] Kreitmair-Steck W, Braun G. Heliradar Technology for Helicopter Allweather Operations [C]//Proceedings of SPIE-Radar Sensor Technology Ⅱ. Bellingham: SPIE, 1997: 28-38.

[12] 郭交, 李真芳, 劉艷陽, 等.一種基于相干性本質(zhì)的InSAR方位向預(yù)濾波方法 [J]. 電子學(xué)報(bào), 2012, 40(3): 417-421.

Guo Jiao, Li Zhenfang, Liu Yanyang, et al. An Improved Method for InSAR Azimuth Prefiltering Based on Coherence Principle [J]. Acta Electronica Sinica, 2012, 40(3): 417-421.

[13] Wang R, Deng Y K, Loffeld O, et al. Focusing Bistatic SAR Data in Double Sliding Spotlight Mode with TerraSAR-X and PAMIR Based on Azimuth Chirp Filtering [C]//9the European Conference on Synthetic Aperture Radar. Piscataway: IEEE, 2012: 772-775.