多頻段典型地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?/h1>
2016-11-23 13:46:30張?jiān)?/span>吳振森張玉石
西安電子科技大學(xué)學(xué)報(bào) 2016年5期
張?jiān)?吳振森,張玉石
(1.西安電子科技大學(xué)物理與光電工程學(xué)院,陜西西安 710071; 2.中國(guó)電波傳播研究所,山東青島 266107)
多頻段典型地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?/p>
張?jiān)?,吳振森1,張玉石2
(1.西安電子科技大學(xué)物理與光電工程學(xué)院,陜西西安 710071; 2.中國(guó)電波傳播研究所,山東青島 266107)
為了發(fā)揮雷達(dá)遙感數(shù)據(jù)的應(yīng)用潛力,同時(shí)解決地表電磁散射模型輸入?yún)?shù)眾多以及地表回波預(yù)測(cè)不準(zhǔn)確等問(wèn)題,從電磁散射理論出發(fā),分別將地表隨機(jī)粗糙面散射和離散隨機(jī)介質(zhì)體散射模型簡(jiǎn)化,建立了單頻段地表散射回波與雷達(dá)參數(shù)間的經(jīng)驗(yàn)?zāi)Ｐ?并通過(guò)尋找其未知參數(shù)隨頻率的變化關(guān)系,進(jìn)一步創(chuàng)建了多頻段七參數(shù)地表后向散射經(jīng)驗(yàn)?zāi)Ｐ?并對(duì)該模型進(jìn)行了驗(yàn)證.結(jié)果表明,多頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)整體吻合良好,實(shí)現(xiàn)了對(duì)寬頻帶范圍內(nèi)任意頻率點(diǎn)地表后向散射回波的準(zhǔn)確預(yù)測(cè),提高了模型的計(jì)算效率.
多波段;電磁散射;地表;經(jīng)驗(yàn)?zāi)Ｐ?/p>
近年來(lái),隨著微波遙感技術(shù)的發(fā)展,各類典型地表的雷達(dá)散射特性研究已成為電波傳播、通信、環(huán)境監(jiān)測(cè)以及目標(biāo)識(shí)別等眾多領(lǐng)域共同的研究熱點(diǎn)[1-3].典型地表散射回波的精確測(cè)量不僅在民用方面具有重要的應(yīng)用價(jià)值,在軍事方面也具有迫切的研究需求,如低空或者超低空飛行的精密制導(dǎo)武器在識(shí)別目標(biāo)時(shí),由于多種地表雜波的影響,目標(biāo)的回波信號(hào)淹沒(méi)于地表雜波信號(hào)中,導(dǎo)致制導(dǎo)的精度降低,對(duì)國(guó)家安全產(chǎn)生巨大的威脅;或是由于植被雷達(dá)回波信號(hào)強(qiáng)于裸土等地表情況,產(chǎn)生虛警現(xiàn)象,對(duì)目標(biāo)識(shí)別產(chǎn)生誤判.因此,只有快速、準(zhǔn)確地預(yù)測(cè)地表電磁散射回波信號(hào),分析各典型地表的散射特征,才能有效濾除地表雜波的影響,對(duì)地表附近的軍事目標(biāo)進(jìn)行準(zhǔn)確識(shí)別與打擊,為國(guó)防工程建設(shè)提供重要的參考依據(jù).
過(guò)去幾十年里,國(guó)內(nèi)外開(kāi)展了L/S/X/Ku/Ka/W波段單層和多層裸土[4-6]表面后向散射實(shí)驗(yàn)測(cè)量及其相關(guān)理論研究工作.由于電磁波在不同頻段下的穿透深度不同,散射機(jī)理不同,文獻(xiàn)[7]采用積分方程模型(Integrated Equation Model,IEM)方法結(jié)合單層指數(shù)譜模型計(jì)算了8.6 GHz水泥路面的后向散射回波;文獻(xiàn)[8]采用半空間指數(shù)譜粗糙面以及粒子體散射模型計(jì)算了94 GHz水泥路面后向散射回波;同時(shí),文獻(xiàn)[9]將沙地模擬為光滑粗糙面和一維周期粗糙面,并采用輻射傳輸理論計(jì)算了沙地的面散射和體散射回波;文獻(xiàn)[10]采用離散面元的基爾霍夫近似法計(jì)算了入射波在沙丘迎風(fēng)坡、背風(fēng)坡以及兩側(cè)入射時(shí)的電磁散射回波.上述典型地表的電磁散射理論研究,能夠加深研究者對(duì)于電磁波與地物之間散射機(jī)制的理解,并對(duì)外場(chǎng)實(shí)驗(yàn)測(cè)量起到重要的指導(dǎo)作用.但隨著大區(qū)域地表變的越來(lái)越復(fù)雜,地表覆蓋物類型繁多且分布不均勻,地表參數(shù)的測(cè)量工作遇到了輸入量(地表組分、分層深度、濕度、均方根高度、相關(guān)長(zhǎng)度,植被莖、葉、桿的尺寸、分布密度及其濕度等)眾多、獲取困難以及不準(zhǔn)確性等問(wèn)題[11],導(dǎo)致理論計(jì)算結(jié)果和實(shí)測(cè)數(shù)據(jù)之間常出現(xiàn)一定的誤差,嚴(yán)重限制了雷達(dá)遙感數(shù)據(jù)對(duì)于地表參數(shù)的反演精度以及超低空目標(biāo)的探測(cè)能力.因此急需尋找地表散射回波與頻率、入射角以及散射角之間的關(guān)系,建立地表后向電磁散射經(jīng)驗(yàn)?zāi)Ｐ?從而采用簡(jiǎn)單的方法解決復(fù)雜的地表散射問(wèn)題,實(shí)現(xiàn)對(duì)地物后向散射回波的快速、精確預(yù)測(cè).
近些年來(lái),國(guó)內(nèi)外學(xué)者提出了眾多的地表散射經(jīng)驗(yàn)?zāi)Ｐ?例如適用于裸土表面的Oh模型[5]、Doubis模型[12]、Shi模型[13]、Zribi模型[14]以及Ulaby單頻模型[15]等,Nashashibi毫米波段沙地后向散射經(jīng)驗(yàn)?zāi)Ｐ蚚9],Ulaby多頻植被[16]散射經(jīng)驗(yàn)?zāi)Ｐ偷?上述經(jīng)驗(yàn)?zāi)Ｐ涂捎行p小人力、物力消耗,并提高地表散射回波預(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)的吻合度,具有重要的工程應(yīng)用價(jià)值.但已有的地表散射經(jīng)驗(yàn)?zāi)Ｐ痛蠖嗑窒抻趩晤l段,如Oh模型、Doubis模型、Shi模型、Zribi模型、Ulaby單頻模型.雖然Nashashibi沙地散射模型、Ulaby多頻植被散射經(jīng)驗(yàn)?zāi)Ｐ团c頻率相關(guān),但其只適用于特定地表,并且Nashashibi沙地散射經(jīng)驗(yàn)?zāi)Ｐ鸵蕾囉诘乇砦锢韰?shù)的測(cè)量值,具有一定的局限性.
國(guó)內(nèi)外對(duì)于多頻段任意地表的后向散射經(jīng)驗(yàn)?zāi)Ｐ脱芯枯^少.因此,為了建立具有較寬適用范圍的多頻地表經(jīng)驗(yàn)?zāi)Ｐ?且使得經(jīng)驗(yàn)?zāi)Ｐ透哂姓f(shuō)服力,筆者從地表電磁散射理論出發(fā),分別將隨機(jī)粗糙面散射和離散隨機(jī)介質(zhì)體散射模型進(jìn)行簡(jiǎn)化,建立了單頻段和多頻段0°～70°入射角范圍內(nèi)的地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?并采用L/S/X/Ku波段H H極化和VV極化裸土、沙地和水泥地的雷達(dá)后向散射實(shí)測(cè)數(shù)據(jù)對(duì)該經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了驗(yàn)證,進(jìn)一步表明該經(jīng)驗(yàn)?zāi)Ｐ涂蓪?shí)現(xiàn)對(duì)地表散射系數(shù)的準(zhǔn)確預(yù)測(cè),并可有效降低地表散射理論模型的復(fù)雜度,減小外場(chǎng)實(shí)驗(yàn)測(cè)量的工作量.
1 單頻段地表電磁散射回波與雷達(dá)參數(shù)間的關(guān)系
由于電磁波在不同頻段下的穿透深度不同,散射機(jī)理不同,因此,文中分別對(duì)低頻段和高頻段入射波頻率情況下的地表散射情況進(jìn)行了討論.當(dāng)雷達(dá)入射波頻段較低時(shí),裸土、水泥地、粗糙沙地以及稀疏植被等典型地表,其面散射在雷達(dá)總回波中占據(jù)主導(dǎo)作用.地表面散射回波常采用微擾法(Small Perturbation Method,SPM),基爾霍夫近似(Kirchhoff Approach,KA)法,改進(jìn)的積分方程模型(Advanced Integral Equation Model,AIEM)方法[17]以及幾何光學(xué)(Geometric Optics,GO)法進(jìn)行計(jì)算,其中AIEM模型在小尺度(kδ＜0.3,kl＜3,s＜0.3,k為入射波波數(shù),δ為地表均方根高度,l為相關(guān)長(zhǎng)度,s為均方根斜率)地表起伏情況下可以退化為SPM模型,而對(duì)于大起伏(kl＞6,l2＞2.76δλ,s＜0.25)地表,AIEM模型可以退化為KA模型.GO方法則適用于粗糙度kδ＞3的極粗糙表面,其后向散射回波類似于朗伯表面.由于AIEM模型的公式極其復(fù)雜,因此,文中將結(jié)合SPM方法、KA方法以及GO方法,建立任意地表粗糙度情況下單頻段雷達(dá)后向散射回波與入射波頻率和入射角之間的經(jīng)驗(yàn)?zāi)Ｐ?
小起伏粗糙面情況下,SPM方法計(jì)算得到的指數(shù)譜地表后向散射回波為
其中,αpq為極化系數(shù),p,q=h,v,θi為入射角.可以看出,在固定入射波頻率情況下,地表的結(jié)構(gòu)參數(shù)和電磁參數(shù)不發(fā)生變化,后向散射系數(shù)僅與入射角、菲涅耳反射系數(shù)以及地面功率譜相關(guān),因此,上述地表散射回波可化簡(jiǎn)為入射角的函數(shù),即
其中,pi,i=1,2,…,6,為待定參數(shù).式(2)轉(zhuǎn)化為分貝形式后,可進(jìn)一步寫(xiě)為
而對(duì)于中等起伏粗糙地表,在KA方法的適用范圍內(nèi),指數(shù)譜非相干散射系數(shù)可以表示為
其中,bi,i=1,2,…,4為待定參數(shù).同理,式(5)可以表示為分貝形式為
對(duì)于極粗糙地表面,其后向散射回波常采用GO方法進(jìn)行計(jì)算,該方法與地表譜函數(shù)無(wú)關(guān),其散射回波公式為
其中,s為均方根斜率,R(0)為入射角為0°情況下的菲涅耳反射系數(shù).在固定入射波頻率情況下,公式可進(jìn)一步化簡(jiǎn)為
其中,qi,i=1,2,…,4,參數(shù)待定,轉(zhuǎn)化為分貝形式為
可以看出,SPM方法、KA方法以及GO方法的后向散射回波簡(jiǎn)化公式的前3項(xiàng)均保持一致,其余項(xiàng)會(huì)出現(xiàn)分歧,經(jīng)筆者嘗試后,發(fā)現(xiàn)僅前3項(xiàng)公式即可實(shí)現(xiàn)低頻段任意粗糙面后向散射回波的預(yù)測(cè),即
但當(dāng)入射波頻率增高時(shí),電磁波穿透裸土、瀝青路面和沙地表面等地表,在其內(nèi)部發(fā)生體散射作用,其總散射回波為上層粗糙面散射和內(nèi)部體散射作用之和.一般情況下,其面散射可采用AIEM方法計(jì)算.而對(duì)于體散射作用,通常認(rèn)為地表內(nèi)部粒子滿足各向同性,其消光矩陣=k和相位矩陣=p分別退化為常數(shù)p1和κ,體散射[8]回波可以表示為
其中,VV表示水平方向,H H表示垂直方向.因此,在高頻段下,裸土、瀝青路面和平滑沙地等地表的總散射回波公式為
其中,σs為面散射作用,σv為體散射作用.在小入射角情況下,地表面散射占據(jù)主要作用,但當(dāng)θ＞θc時(shí)(θc與地表類型以及入射波頻率有關(guān)),體散射超過(guò)面散射而占據(jù)為主要作用.當(dāng)入射波頻率一定時(shí),地面的消光系數(shù)以及相函數(shù)比值不隨入射角發(fā)生變化,因此p1(2κ)保持不變,裸土、瀝青路面和沙地等地表的總后向散射回波可以表示為
其中,未知參數(shù)為N1,N2,N3.可以看出,在高頻入射波情況下,地表后向散射模型與低頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ皖愃?因此地表散射回波經(jīng)驗(yàn)?zāi)Ｐ涂山y(tǒng)一寫(xiě)為式(10)的形式,簡(jiǎn)單方便.
2 多頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ?/h2>
在單頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷幕A(chǔ)上,添加頻率項(xiàng)作用,建立寬頻帶范圍內(nèi)地表后向散射回波與入射角和頻率之間的經(jīng)驗(yàn)?zāi)Ｐ?從而可以通過(guò)反演獲取經(jīng)驗(yàn)?zāi)Ｐ偷奈粗獏⒘?確定多頻經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,實(shí)現(xiàn)對(duì)任意頻段下地表后向散射回波的預(yù)測(cè).
入射波為1～16 GHz情況下,裸土、粗糙沙地和水泥瀝青路面的后向散射回波隨頻率的變化關(guān)系較為類似,因此,文中僅以裸土[5]為例,通過(guò)分析其多頻段后向散射回波理論值的特點(diǎn),尋找單頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ臀粗獏?shù)Q1,Q2,Q3與頻率f的關(guān)系.如圖1所示,在1 GHz≤f≤9 GHz情況下,Q1隨入射波頻率的增大而增大,但在1 GHz≤f≤3 GHz范圍內(nèi),Q1變化緩慢.之后,當(dāng)10 GHz≤f≤16 GHz時(shí),Q1隨頻率的增大而增大.對(duì)于Q2,在1 GHz≤f≤5 GHz情況下,其隨入射波頻率的增大而增大,但之后保持不變,在10 GHz≤f≤16 GHz繼續(xù)隨入射波的增大而增大.Q3隨入射波頻率的增大先保持不變,而后增大,并在10 GHz≤f≤16 GHz情況下隨入射波的增大而快速增大.考慮到在較寬的頻段范圍內(nèi),Q1,Q2,Q3的值整體均隨頻率的增大而增大,因此首先在單頻段后向散射經(jīng)驗(yàn)?zāi)Ｐ偷纳咸砑觢g(f)項(xiàng).同時(shí)為滿足各參數(shù)的變化規(guī)律,且盡量減少未知參數(shù)個(gè)數(shù),文中將Q1設(shè)為常數(shù),Q2替換為Q2-m1f,Q3設(shè)為Q3-m2·(f-m3)2,從而建立了多頻段地表散射回波經(jīng)驗(yàn)?zāi)Ｐ?其中f為入射波頻率,m1,m2,m3,m4為新增的未知參量.
圖1 裸土地表經(jīng)驗(yàn)?zāi)Ｐ蛥?shù)Q1,Q2,Q3隨頻率f的變化規(guī)律
因此,多頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ涂梢詫?xiě)為
其中,未知參數(shù)為Q1,Q2,Q3,m1,m2,m3,m4.
3 地表后向散射特性測(cè)量實(shí)驗(yàn)
2013年中國(guó)電波傳播研究所在青島郊區(qū)采用車載散射計(jì)對(duì)L/S/X/Ku波段下裸土、水泥地和沙地后向散射回波進(jìn)行了測(cè)量.裸土實(shí)驗(yàn)區(qū)位于青島周邊的一塊空地,地面較為平坦.水泥路面選址在一處公路上,路面空曠.沙地測(cè)試區(qū)位于青島海邊,濕度較大,且沙地表面較為粗糙.實(shí)驗(yàn)測(cè)量系統(tǒng)主要由Agilent 8362B PNA矢量網(wǎng)格分析儀、功率放大器、功分器、衰減器和收發(fā)天線組成.在實(shí)驗(yàn)中,天線入射角的變化范圍為0°到66°,間隔3°,散射計(jì)的主要性能指標(biāo)如表1所示.采用龍伯球外場(chǎng)定標(biāo)后,L/S/X/Ku波段裸土、水泥地和沙地的同極化后向散射系數(shù)如圖2所示.
表1 散射計(jì)主要參數(shù)
4 反演與分析
以下將根據(jù)最小平方和原理與L/S/X/Ku波段各地表實(shí)測(cè)數(shù)據(jù),反演獲取多頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ偷母魑粗獏?shù)值,其中評(píng)價(jià)函數(shù)可以寫(xiě)為
圖2 多頻段地表后向散射回波實(shí)測(cè)數(shù)據(jù)
其中,fi為入射波頻率,θi為入射波角度,xi為經(jīng)驗(yàn)?zāi)Ｐ偷奈粗獏?shù),分別為Q1,Q2,Q3,m1,m2,m3.pp= H H或VV,σHH,σVV分別為HH極化和VV極化地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)值,HH和VV分別是HH和VV極化地表雷達(dá)后向散射回波的實(shí)測(cè)值.文中將采用遺傳算法[18]來(lái)反演地表經(jīng)驗(yàn)?zāi)Ｐ偷母鲄?shù)值,其具有群體啟發(fā)式隨機(jī)搜索特性,不容易陷入局部最優(yōu),能以最大概率尋找到全局最優(yōu)解.
文中將裸土、水泥路面以及沙地地表的后向散射實(shí)測(cè)數(shù)據(jù)分為兩部分,一部分用于反演多頻后向散射經(jīng)驗(yàn)?zāi)Ｐ蛥?shù),另外一部分用于驗(yàn)證模型的反演效果.鑒于實(shí)測(cè)數(shù)據(jù)的入射波頻段較寬,裸土表面選擇頻率為1.34 GHz,3.20 GHz和13.50 GHz實(shí)測(cè)數(shù)據(jù)用于反演,而瀝青路面和沙地表面分別選擇頻率為1.34 GHz,10.00 GHz和13.50 GHz的實(shí)測(cè)數(shù)據(jù)用于反演多頻經(jīng)驗(yàn)?zāi)Ｐ偷母鲄?shù)值.設(shè)定遺傳算法中每一代種群大小為200,最大迭代次數(shù)為300,染色體變異概率為0.03,交叉概率為0.95,并取Q1,Q2,Q3,m1,m2,m3,m4的取值范圍分別為[-50,50],[-50,50],[-50,50],[-30,30],[-30,30],[-30,30]和[-30,30].經(jīng)多次反演,評(píng)價(jià)函數(shù)值較小且趨于穩(wěn)定,在HH極化和VV極化下,裸土、水泥地和沙地各參數(shù)Q1,Q2,Q3,m1,m2,m3,m4的反演結(jié)果見(jiàn)表2所示.可以看出,所有地表的均方根誤差(Root Mean Square Error,RMSE) ERMSE≤3.99 d B,相關(guān)系數(shù)R≥0.921,同時(shí)還可以從圖3～圖5中看出,在小角度入射情況下,地表后向散射回波隨頻率的變化關(guān)系與地表類型有關(guān),但當(dāng)入射角大于某一固定角度θc時(shí),其后向散射回波均隨頻率的增大而增大.文中各頻段下地表散射回波數(shù)據(jù)的預(yù)測(cè)值均與實(shí)測(cè)數(shù)據(jù)吻合良好.
表2 多頻地表經(jīng)驗(yàn)?zāi)Ｐ蛥?shù)的反演結(jié)果
圖3 裸土后向散射模型預(yù)測(cè)結(jié)果
圖4 水泥地后向散射模型預(yù)測(cè)結(jié)果
圖5 沙地后向散射模型預(yù)測(cè)結(jié)果
分別將上述各地表參數(shù)的反演結(jié)果帶回各地表散射經(jīng)驗(yàn)?zāi)Ｐ椭?確定其經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,并預(yù)測(cè)10 GHz入射波頻率下裸土,3.2 GHz入射波頻率下水泥地和沙地的后向散射回波值,與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析,其結(jié)果分別如圖6～圖8所示.對(duì)于裸土散射情況,H H極化下模型預(yù)測(cè)結(jié)果的均方根誤差為3.70 dB,相關(guān)系數(shù)為0.989;而VV極化下,其均方根誤差為3.05 dB,相關(guān)系數(shù)為0.984.水泥路面下,HH極化模型預(yù)測(cè)結(jié)果的均方根誤差為4.29 dB,相關(guān)系數(shù)為0.961;而VV極化下,其均方根誤差為4.40 dB,相關(guān)系數(shù)為0.942.同理,對(duì)入射波頻率為3.2 GHz沙地散射回波進(jìn)行分析,H H極化下模型預(yù)測(cè)結(jié)果的均方根誤差為4.63 dB,相關(guān)系數(shù)為0.977;而VV極化下,其均方根誤差為3.89 dB,相關(guān)系數(shù)為0.975.當(dāng)?shù)乇砩⑸浠夭〝?shù)據(jù)震蕩劇烈時(shí),其預(yù)測(cè)結(jié)果會(huì)稍微差些.文中裸土表面、水泥路面以及粗糙沙地散射回波的預(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)整體吻合良好,可以滿足實(shí)際工程需求,可有效減小多頻外場(chǎng)實(shí)驗(yàn)測(cè)量的工作量.
圖6 10 GHz裸土散射回波預(yù)測(cè)結(jié)果
圖7 3.2 GHz水泥地散射回波預(yù)測(cè)結(jié)果
圖8 3.2 GHz沙地散射回波預(yù)測(cè)結(jié)果
5 結(jié)束語(yǔ)
筆者將面散射地表、面-體散射地表的電磁散射理論模型進(jìn)行簡(jiǎn)化,建立了入射波頻率為1～16 GHz的單頻段和多頻段地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?并采用L/S/X/Ku四頻段裸土、沙地和水泥地的后向散射實(shí)測(cè)數(shù)據(jù)對(duì)多頻散射經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了驗(yàn)證.結(jié)果表明,文中提出的多頻經(jīng)驗(yàn)?zāi)Ｐ蜔o(wú)需測(cè)量地表的粗糙度和介電常數(shù),只需3個(gè)頻段或4個(gè)頻段的地表雷達(dá)回波數(shù)據(jù),即可確定多頻經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,并實(shí)現(xiàn)對(duì)寬頻帶范圍內(nèi)任意頻率點(diǎn)下地表后向散射回波的預(yù)測(cè).該多頻地表電磁散射經(jīng)驗(yàn)?zāi)Ｐ秃?jiǎn)化了地表后向散射模型的復(fù)雜度,在保證地表后向散射系數(shù)預(yù)測(cè)精度的同時(shí),提高了計(jì)算效率,便于工程應(yīng)用.但值得注意的是,由于在不同頻段下,電磁波在地表的穿透深度不同,散射機(jī)制會(huì)產(chǎn)生差異,因此,更高頻段下大區(qū)域復(fù)雜地表的電磁散射經(jīng)驗(yàn)?zāi)Ｐ?將是下一步研究的重點(diǎn).
[1]SU Y,GUO Q,XUE B,et al.Spatial Distribution of Forest Aboveground Biomass in China:Estimation through Combination of Spaceborne Lidar,Optical Imagery,and Forest Inventory Data[J].Remote Sensing of Environment,2016,173:187-199.
[2]OZLEM K,GARCIA-RUBIA J M,NGHIA T,et al.Detection of Moving Human Micro-Doppler Signature in Forest Environments with Swaying Tree Components by Wind[J].Radio Science,2015,50:238-248.
[3]RASHIDI-RANJBAR E,DEH MOLLAIAN M.Target above Random Rough Surface Scattering Using a Parallelized IPO Accelerated by MLFMM[J].IEEE Geoscience and Remote Sensing Letters,2015,12:1481-1485.
[4]NASHASHIBI A,ULABY F T,SARABANDI K.Measurement and Modeling of the Millimeter-wave Backscatter Response of Soil Surfaces[J].IEEE Transactions on Geoscience and Remote Sensing,1996,34(2):561-572.
[5]OH Y,SARABANDI K,ULABY F T.An Empirical Model and an Inversion Technique for Radar Scattering from Bare Soil Surfaces[J].IEEE Transactions on Geoscience and Remote Sensing,1992,30(2):370-381.
[6]張連波,郭立新,茍雪銀,等.三層粗糙面電磁散射的矩量法研究[J].西安電子科技大學(xué)學(xué)報(bào),2013,40(6):147-154. ZHANG Lianbo,GUO Lixin,GOU Xueyin,et al.Method of Moment Investigation on Electromagnetic Scattering from the Three-layered Rough Interfaces[J].Journal of Xidian University,2013,40(6):147-154.
[7]FUNG A K.Microwave Scattering and Emission Models and Their Applications[M].Norwood:Artech House,1994.
[8]LI E S,SARABANDI K.Low Grazing Incidence Millimeter-wave Scattering Models and Measurements for Various Road Surfaces[J].IEEE Transactions on Antennas and Propagation,1999,47(5):851-861.
[9]NASHASHIBI A Y,SARABANDI K,AL-ZAID F A,et al.An Empirical Model of Volume Scattering from Dry Sandcovered Surfaces at Millimeter-wave Frequencies[J].IEEE Transactions on Geoscience and Remote Sensing,2013,51(6):3673-3682.
[10]鄭帆,郭立新,劉偉.風(fēng)積新月形沙丘的電磁散射特性分析[J].西安電子科技大學(xué)學(xué)報(bào),2013,40(6):132-139. ZHENG Fan,GUO Lixin,LIU Wei.Analysis of Electromagetic Scattering from the Wind-blown Barchan[J].Journal of Xidian University,2013,40(6):132-139.
[11]LIEVENS H,VERHOEST N E C.On the Retrieval of Soil Moisture in Wheat Fields from L-band SAR Based on Water Cloud Modeling,the IEM,and Effective Roughness Parameters[J].IEEE Geoscience and Remote Sensing Letters,2011,8(4):740-744.
[12]DUBOIS P C,van ZYL J,ENGMAN T.Measuring Soil Moisture with Imaging Radars[J].IEEE Transactions on Geoscience and Remote Sensing,1995,33(4):915-926.
[13]SHI J,WANG J,HSU A Y,et al.Estimation of Bare Soil Moisture and Surface Roughness Parameters Using L-band SAR Image Data[J].IEEE Transactions on Geoscience and Remote Sensing,1997,35(5):1254-1266.
[14]ZRIBI M,GORRAB A,BAGHDADI N.A New Soil Roughness Parameter for the Modelling of Radar Backscattering over Bare Soil[J].Remote Sensing of Environment,2014,152:62-73.
[15]ULABY F T,DONSON M C.Handbook of Radar Scattering Statistics for Terrain[M].Norwood:Artech House,1989.
[16]ULABY F T.Vegetation Clutter Model[J].IEEE Transactions on Antennas and Propagation,1980,28(4):538-545. [17]CHEN K L,CHEN K S,LI Z L,et al.Extension and Validation of an Advanced Integral Equation Model for Bistatic Scattering from Rough Surfaces[J].Progress In Electromagnetics Research,2015,152:59-76.
[18]HAFFNER S,PEREIRA L A,PEREIRA L F A.A Method for Optimization of Five-phase Induction Machines Based on Genetic Algorithms[J].Journal of Control Automation&Electrical Systems,2016,26(5):521-534.
(編輯:王 瑞)
Multi-band empirical model for backscattering coefficients prediction from typical land surfaces
ZHANG Yuanyuan1,WU Zhensen1,ZHANG Yushi2
(1.School of Physics and Optoelectronic Engineering,Xidian Univ.,Xi’an 710071,China; 2.China Research Institute of Radio Wave Propagation,Qingdao 266107,China)
Based on the electromagnetic scattering theory of the rough surface and discrete random medium,a multi-band empirical model for land backscattering prediction in 1～16 GHz is proposed to tap the potential of remote sensing in application for agricultural production and targets detection efficiently.It should be noticed that this empirical model only have seven unknowns,and that it overcomes the problems in numerous inputs and inaccuracies echoes prediction.It is verified by the measured data of soil,sand and concrete in the L/S/X/Ku band for HH and VV polarizations.The results show that the backscattering coefficients from soil,sand and concrete predicted by the empirical model are in good agreement with the measured data.Meanwhile,this multi-band empirical model can be used to forecast the backscattering echoes at any frequency in 1～16 GHz accurately.So this empirical model has a great application value for it can greatly improve the computational efficiency for the backscattering coefficients prediction,and also it can effectively reduce the workload in the field experiments.
multi-band;electromagnetic scattering;land;empirical model
TN959.71
A
1001-2400(2016)05-0190-07
10.3969/j.issn.1001-2400.2016.05.033
2015-12-01
國(guó)家自然科學(xué)基金資助項(xiàng)目(61571355)
張?jiān)?1986-),女,西安電子科技大學(xué)博士研究生,E-mail:zyy07063@163.com.
猜你喜歡
認(rèn)知能力、技術(shù)進(jìn)步與就業(yè)極化現(xiàn)代財(cái)經(jīng)-天津財(cái)經(jīng)大學(xué)學(xué)報(bào)(2022年5期)2022-06-01 06:08:32 SHPB入射波相似律與整形技術(shù)的試驗(yàn)與數(shù)值研究振動(dòng)與沖擊(2022年6期)2022-03-27 12:18:26 能鉆過(guò)柔軟沙地的蛇形機(jī)器人軍事文摘(2021年18期)2021-12-02 01:27:58 呼倫貝爾沙地實(shí)現(xiàn)良性逆轉(zhuǎn)內(nèi)蒙古林業(yè)(2021年6期)2021-06-26 10:42:32 沙地迷宮幼兒智力世界(2021年4期)2021-05-17 17:10:29 風(fēng)滾草科普童話·百科探秘(2018年7期)2018-09-13 11:12:08 瞬態(tài)激勵(lì)狀態(tài)下樁身速度以及樁身內(nèi)力計(jì)算建材與裝飾(2018年50期)2018-02-22 03:30:40 雙頻帶隔板極化器電子測(cè)試(2017年15期)2017-12-18 07:18:51 基于PWM控制的新型極化電源設(shè)計(jì)與實(shí)現(xiàn)電源技術(shù)(2015年1期)2015-08-22 11:16:18 對(duì)機(jī)械波半波損失現(xiàn)象的物理解釋電子科技(2015年11期)2015-03-06 01:32:24
張?jiān)?吳振森,張玉石
(1.西安電子科技大學(xué)物理與光電工程學(xué)院,陜西西安 710071; 2.中國(guó)電波傳播研究所,山東青島 266107)
多頻段典型地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?/p>
張?jiān)?,吳振森1,張玉石2
(1.西安電子科技大學(xué)物理與光電工程學(xué)院,陜西西安 710071; 2.中國(guó)電波傳播研究所,山東青島 266107)
為了發(fā)揮雷達(dá)遙感數(shù)據(jù)的應(yīng)用潛力,同時(shí)解決地表電磁散射模型輸入?yún)?shù)眾多以及地表回波預(yù)測(cè)不準(zhǔn)確等問(wèn)題,從電磁散射理論出發(fā),分別將地表隨機(jī)粗糙面散射和離散隨機(jī)介質(zhì)體散射模型簡(jiǎn)化,建立了單頻段地表散射回波與雷達(dá)參數(shù)間的經(jīng)驗(yàn)?zāi)Ｐ?并通過(guò)尋找其未知參數(shù)隨頻率的變化關(guān)系,進(jìn)一步創(chuàng)建了多頻段七參數(shù)地表后向散射經(jīng)驗(yàn)?zāi)Ｐ?并對(duì)該模型進(jìn)行了驗(yàn)證.結(jié)果表明,多頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)整體吻合良好,實(shí)現(xiàn)了對(duì)寬頻帶范圍內(nèi)任意頻率點(diǎn)地表后向散射回波的準(zhǔn)確預(yù)測(cè),提高了模型的計(jì)算效率.
多波段;電磁散射;地表;經(jīng)驗(yàn)?zāi)Ｐ?/p>
近年來(lái),隨著微波遙感技術(shù)的發(fā)展,各類典型地表的雷達(dá)散射特性研究已成為電波傳播、通信、環(huán)境監(jiān)測(cè)以及目標(biāo)識(shí)別等眾多領(lǐng)域共同的研究熱點(diǎn)[1-3].典型地表散射回波的精確測(cè)量不僅在民用方面具有重要的應(yīng)用價(jià)值,在軍事方面也具有迫切的研究需求,如低空或者超低空飛行的精密制導(dǎo)武器在識(shí)別目標(biāo)時(shí),由于多種地表雜波的影響,目標(biāo)的回波信號(hào)淹沒(méi)于地表雜波信號(hào)中,導(dǎo)致制導(dǎo)的精度降低,對(duì)國(guó)家安全產(chǎn)生巨大的威脅;或是由于植被雷達(dá)回波信號(hào)強(qiáng)于裸土等地表情況,產(chǎn)生虛警現(xiàn)象,對(duì)目標(biāo)識(shí)別產(chǎn)生誤判.因此,只有快速、準(zhǔn)確地預(yù)測(cè)地表電磁散射回波信號(hào),分析各典型地表的散射特征,才能有效濾除地表雜波的影響,對(duì)地表附近的軍事目標(biāo)進(jìn)行準(zhǔn)確識(shí)別與打擊,為國(guó)防工程建設(shè)提供重要的參考依據(jù).
過(guò)去幾十年里,國(guó)內(nèi)外開(kāi)展了L/S/X/Ku/Ka/W波段單層和多層裸土[4-6]表面后向散射實(shí)驗(yàn)測(cè)量及其相關(guān)理論研究工作.由于電磁波在不同頻段下的穿透深度不同,散射機(jī)理不同,文獻(xiàn)[7]采用積分方程模型(Integrated Equation Model,IEM)方法結(jié)合單層指數(shù)譜模型計(jì)算了8.6 GHz水泥路面的后向散射回波;文獻(xiàn)[8]采用半空間指數(shù)譜粗糙面以及粒子體散射模型計(jì)算了94 GHz水泥路面后向散射回波;同時(shí),文獻(xiàn)[9]將沙地模擬為光滑粗糙面和一維周期粗糙面,并采用輻射傳輸理論計(jì)算了沙地的面散射和體散射回波;文獻(xiàn)[10]采用離散面元的基爾霍夫近似法計(jì)算了入射波在沙丘迎風(fēng)坡、背風(fēng)坡以及兩側(cè)入射時(shí)的電磁散射回波.上述典型地表的電磁散射理論研究,能夠加深研究者對(duì)于電磁波與地物之間散射機(jī)制的理解,并對(duì)外場(chǎng)實(shí)驗(yàn)測(cè)量起到重要的指導(dǎo)作用.但隨著大區(qū)域地表變的越來(lái)越復(fù)雜,地表覆蓋物類型繁多且分布不均勻,地表參數(shù)的測(cè)量工作遇到了輸入量(地表組分、分層深度、濕度、均方根高度、相關(guān)長(zhǎng)度,植被莖、葉、桿的尺寸、分布密度及其濕度等)眾多、獲取困難以及不準(zhǔn)確性等問(wèn)題[11],導(dǎo)致理論計(jì)算結(jié)果和實(shí)測(cè)數(shù)據(jù)之間常出現(xiàn)一定的誤差,嚴(yán)重限制了雷達(dá)遙感數(shù)據(jù)對(duì)于地表參數(shù)的反演精度以及超低空目標(biāo)的探測(cè)能力.因此急需尋找地表散射回波與頻率、入射角以及散射角之間的關(guān)系,建立地表后向電磁散射經(jīng)驗(yàn)?zāi)Ｐ?從而采用簡(jiǎn)單的方法解決復(fù)雜的地表散射問(wèn)題,實(shí)現(xiàn)對(duì)地物后向散射回波的快速、精確預(yù)測(cè).
近些年來(lái),國(guó)內(nèi)外學(xué)者提出了眾多的地表散射經(jīng)驗(yàn)?zāi)Ｐ?例如適用于裸土表面的Oh模型[5]、Doubis模型[12]、Shi模型[13]、Zribi模型[14]以及Ulaby單頻模型[15]等,Nashashibi毫米波段沙地后向散射經(jīng)驗(yàn)?zāi)Ｐ蚚9],Ulaby多頻植被[16]散射經(jīng)驗(yàn)?zāi)Ｐ偷?上述經(jīng)驗(yàn)?zāi)Ｐ涂捎行p小人力、物力消耗,并提高地表散射回波預(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)的吻合度,具有重要的工程應(yīng)用價(jià)值.但已有的地表散射經(jīng)驗(yàn)?zāi)Ｐ痛蠖嗑窒抻趩晤l段,如Oh模型、Doubis模型、Shi模型、Zribi模型、Ulaby單頻模型.雖然Nashashibi沙地散射模型、Ulaby多頻植被散射經(jīng)驗(yàn)?zāi)Ｐ团c頻率相關(guān),但其只適用于特定地表,并且Nashashibi沙地散射經(jīng)驗(yàn)?zāi)Ｐ鸵蕾囉诘乇砦锢韰?shù)的測(cè)量值,具有一定的局限性.
國(guó)內(nèi)外對(duì)于多頻段任意地表的后向散射經(jīng)驗(yàn)?zāi)Ｐ脱芯枯^少.因此,為了建立具有較寬適用范圍的多頻地表經(jīng)驗(yàn)?zāi)Ｐ?且使得經(jīng)驗(yàn)?zāi)Ｐ透哂姓f(shuō)服力,筆者從地表電磁散射理論出發(fā),分別將隨機(jī)粗糙面散射和離散隨機(jī)介質(zhì)體散射模型進(jìn)行簡(jiǎn)化,建立了單頻段和多頻段0°～70°入射角范圍內(nèi)的地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?并采用L/S/X/Ku波段H H極化和VV極化裸土、沙地和水泥地的雷達(dá)后向散射實(shí)測(cè)數(shù)據(jù)對(duì)該經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了驗(yàn)證,進(jìn)一步表明該經(jīng)驗(yàn)?zāi)Ｐ涂蓪?shí)現(xiàn)對(duì)地表散射系數(shù)的準(zhǔn)確預(yù)測(cè),并可有效降低地表散射理論模型的復(fù)雜度,減小外場(chǎng)實(shí)驗(yàn)測(cè)量的工作量.
1 單頻段地表電磁散射回波與雷達(dá)參數(shù)間的關(guān)系
由于電磁波在不同頻段下的穿透深度不同,散射機(jī)理不同,因此,文中分別對(duì)低頻段和高頻段入射波頻率情況下的地表散射情況進(jìn)行了討論.當(dāng)雷達(dá)入射波頻段較低時(shí),裸土、水泥地、粗糙沙地以及稀疏植被等典型地表,其面散射在雷達(dá)總回波中占據(jù)主導(dǎo)作用.地表面散射回波常采用微擾法(Small Perturbation Method,SPM),基爾霍夫近似(Kirchhoff Approach,KA)法,改進(jìn)的積分方程模型(Advanced Integral Equation Model,AIEM)方法[17]以及幾何光學(xué)(Geometric Optics,GO)法進(jìn)行計(jì)算,其中AIEM模型在小尺度(kδ＜0.3,kl＜3,s＜0.3,k為入射波波數(shù),δ為地表均方根高度,l為相關(guān)長(zhǎng)度,s為均方根斜率)地表起伏情況下可以退化為SPM模型,而對(duì)于大起伏(kl＞6,l2＞2.76δλ,s＜0.25)地表,AIEM模型可以退化為KA模型.GO方法則適用于粗糙度kδ＞3的極粗糙表面,其后向散射回波類似于朗伯表面.由于AIEM模型的公式極其復(fù)雜,因此,文中將結(jié)合SPM方法、KA方法以及GO方法,建立任意地表粗糙度情況下單頻段雷達(dá)后向散射回波與入射波頻率和入射角之間的經(jīng)驗(yàn)?zāi)Ｐ?
小起伏粗糙面情況下,SPM方法計(jì)算得到的指數(shù)譜地表后向散射回波為
其中,αpq為極化系數(shù),p,q=h,v,θi為入射角.可以看出,在固定入射波頻率情況下,地表的結(jié)構(gòu)參數(shù)和電磁參數(shù)不發(fā)生變化,后向散射系數(shù)僅與入射角、菲涅耳反射系數(shù)以及地面功率譜相關(guān),因此,上述地表散射回波可化簡(jiǎn)為入射角的函數(shù),即
其中,pi,i=1,2,…,6,為待定參數(shù).式(2)轉(zhuǎn)化為分貝形式后,可進(jìn)一步寫(xiě)為
而對(duì)于中等起伏粗糙地表,在KA方法的適用范圍內(nèi),指數(shù)譜非相干散射系數(shù)可以表示為
其中,bi,i=1,2,…,4為待定參數(shù).同理,式(5)可以表示為分貝形式為
對(duì)于極粗糙地表面,其后向散射回波常采用GO方法進(jìn)行計(jì)算,該方法與地表譜函數(shù)無(wú)關(guān),其散射回波公式為
其中,s為均方根斜率,R(0)為入射角為0°情況下的菲涅耳反射系數(shù).在固定入射波頻率情況下,公式可進(jìn)一步化簡(jiǎn)為
其中,qi,i=1,2,…,4,參數(shù)待定,轉(zhuǎn)化為分貝形式為
可以看出,SPM方法、KA方法以及GO方法的后向散射回波簡(jiǎn)化公式的前3項(xiàng)均保持一致,其余項(xiàng)會(huì)出現(xiàn)分歧,經(jīng)筆者嘗試后,發(fā)現(xiàn)僅前3項(xiàng)公式即可實(shí)現(xiàn)低頻段任意粗糙面后向散射回波的預(yù)測(cè),即
但當(dāng)入射波頻率增高時(shí),電磁波穿透裸土、瀝青路面和沙地表面等地表,在其內(nèi)部發(fā)生體散射作用,其總散射回波為上層粗糙面散射和內(nèi)部體散射作用之和.一般情況下,其面散射可采用AIEM方法計(jì)算.而對(duì)于體散射作用,通常認(rèn)為地表內(nèi)部粒子滿足各向同性,其消光矩陣=k和相位矩陣=p分別退化為常數(shù)p1和κ,體散射[8]回波可以表示為
其中,VV表示水平方向,H H表示垂直方向.因此,在高頻段下,裸土、瀝青路面和平滑沙地等地表的總散射回波公式為
其中,σs為面散射作用,σv為體散射作用.在小入射角情況下,地表面散射占據(jù)主要作用,但當(dāng)θ＞θc時(shí)(θc與地表類型以及入射波頻率有關(guān)),體散射超過(guò)面散射而占據(jù)為主要作用.當(dāng)入射波頻率一定時(shí),地面的消光系數(shù)以及相函數(shù)比值不隨入射角發(fā)生變化,因此p1(2κ)保持不變,裸土、瀝青路面和沙地等地表的總后向散射回波可以表示為
其中,未知參數(shù)為N1,N2,N3.可以看出,在高頻入射波情況下,地表后向散射模型與低頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ皖愃?因此地表散射回波經(jīng)驗(yàn)?zāi)Ｐ涂山y(tǒng)一寫(xiě)為式(10)的形式,簡(jiǎn)單方便.
2 多頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ?/h2>
在單頻段地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷幕A(chǔ)上,添加頻率項(xiàng)作用,建立寬頻帶范圍內(nèi)地表后向散射回波與入射角和頻率之間的經(jīng)驗(yàn)?zāi)Ｐ?從而可以通過(guò)反演獲取經(jīng)驗(yàn)?zāi)Ｐ偷奈粗獏⒘?確定多頻經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,實(shí)現(xiàn)對(duì)任意頻段下地表后向散射回波的預(yù)測(cè).
入射波為1～16 GHz情況下,裸土、粗糙沙地和水泥瀝青路面的后向散射回波隨頻率的變化關(guān)系較為類似,因此,文中僅以裸土[5]為例,通過(guò)分析其多頻段后向散射回波理論值的特點(diǎn),尋找單頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ臀粗獏?shù)Q1,Q2,Q3與頻率f的關(guān)系.如圖1所示,在1 GHz≤f≤9 GHz情況下,Q1隨入射波頻率的增大而增大,但在1 GHz≤f≤3 GHz范圍內(nèi),Q1變化緩慢.之后,當(dāng)10 GHz≤f≤16 GHz時(shí),Q1隨頻率的增大而增大.對(duì)于Q2,在1 GHz≤f≤5 GHz情況下,其隨入射波頻率的增大而增大,但之后保持不變,在10 GHz≤f≤16 GHz繼續(xù)隨入射波的增大而增大.Q3隨入射波頻率的增大先保持不變,而后增大,并在10 GHz≤f≤16 GHz情況下隨入射波的增大而快速增大.考慮到在較寬的頻段范圍內(nèi),Q1,Q2,Q3的值整體均隨頻率的增大而增大,因此首先在單頻段后向散射經(jīng)驗(yàn)?zāi)Ｐ偷纳咸砑觢g(f)項(xiàng).同時(shí)為滿足各參數(shù)的變化規(guī)律,且盡量減少未知參數(shù)個(gè)數(shù),文中將Q1設(shè)為常數(shù),Q2替換為Q2-m1f,Q3設(shè)為Q3-m2·(f-m3)2,從而建立了多頻段地表散射回波經(jīng)驗(yàn)?zāi)Ｐ?其中f為入射波頻率,m1,m2,m3,m4為新增的未知參量.
圖1 裸土地表經(jīng)驗(yàn)?zāi)Ｐ蛥?shù)Q1,Q2,Q3隨頻率f的變化規(guī)律
因此,多頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ涂梢詫?xiě)為
其中,未知參數(shù)為Q1,Q2,Q3,m1,m2,m3,m4.
3 地表后向散射特性測(cè)量實(shí)驗(yàn)
2013年中國(guó)電波傳播研究所在青島郊區(qū)采用車載散射計(jì)對(duì)L/S/X/Ku波段下裸土、水泥地和沙地后向散射回波進(jìn)行了測(cè)量.裸土實(shí)驗(yàn)區(qū)位于青島周邊的一塊空地,地面較為平坦.水泥路面選址在一處公路上,路面空曠.沙地測(cè)試區(qū)位于青島海邊,濕度較大,且沙地表面較為粗糙.實(shí)驗(yàn)測(cè)量系統(tǒng)主要由Agilent 8362B PNA矢量網(wǎng)格分析儀、功率放大器、功分器、衰減器和收發(fā)天線組成.在實(shí)驗(yàn)中,天線入射角的變化范圍為0°到66°,間隔3°,散射計(jì)的主要性能指標(biāo)如表1所示.采用龍伯球外場(chǎng)定標(biāo)后,L/S/X/Ku波段裸土、水泥地和沙地的同極化后向散射系數(shù)如圖2所示.
表1 散射計(jì)主要參數(shù)
4 反演與分析
以下將根據(jù)最小平方和原理與L/S/X/Ku波段各地表實(shí)測(cè)數(shù)據(jù),反演獲取多頻段地表散射經(jīng)驗(yàn)?zāi)Ｐ偷母魑粗獏?shù)值,其中評(píng)價(jià)函數(shù)可以寫(xiě)為
圖2 多頻段地表后向散射回波實(shí)測(cè)數(shù)據(jù)
其中,fi為入射波頻率,θi為入射波角度,xi為經(jīng)驗(yàn)?zāi)Ｐ偷奈粗獏?shù),分別為Q1,Q2,Q3,m1,m2,m3.pp= H H或VV,σHH,σVV分別為HH極化和VV極化地表后向散射經(jīng)驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)值,HH和VV分別是HH和VV極化地表雷達(dá)后向散射回波的實(shí)測(cè)值.文中將采用遺傳算法[18]來(lái)反演地表經(jīng)驗(yàn)?zāi)Ｐ偷母鲄?shù)值,其具有群體啟發(fā)式隨機(jī)搜索特性,不容易陷入局部最優(yōu),能以最大概率尋找到全局最優(yōu)解.
文中將裸土、水泥路面以及沙地地表的后向散射實(shí)測(cè)數(shù)據(jù)分為兩部分,一部分用于反演多頻后向散射經(jīng)驗(yàn)?zāi)Ｐ蛥?shù),另外一部分用于驗(yàn)證模型的反演效果.鑒于實(shí)測(cè)數(shù)據(jù)的入射波頻段較寬,裸土表面選擇頻率為1.34 GHz,3.20 GHz和13.50 GHz實(shí)測(cè)數(shù)據(jù)用于反演,而瀝青路面和沙地表面分別選擇頻率為1.34 GHz,10.00 GHz和13.50 GHz的實(shí)測(cè)數(shù)據(jù)用于反演多頻經(jīng)驗(yàn)?zāi)Ｐ偷母鲄?shù)值.設(shè)定遺傳算法中每一代種群大小為200,最大迭代次數(shù)為300,染色體變異概率為0.03,交叉概率為0.95,并取Q1,Q2,Q3,m1,m2,m3,m4的取值范圍分別為[-50,50],[-50,50],[-50,50],[-30,30],[-30,30],[-30,30]和[-30,30].經(jīng)多次反演,評(píng)價(jià)函數(shù)值較小且趨于穩(wěn)定,在HH極化和VV極化下,裸土、水泥地和沙地各參數(shù)Q1,Q2,Q3,m1,m2,m3,m4的反演結(jié)果見(jiàn)表2所示.可以看出,所有地表的均方根誤差(Root Mean Square Error,RMSE) ERMSE≤3.99 d B,相關(guān)系數(shù)R≥0.921,同時(shí)還可以從圖3～圖5中看出,在小角度入射情況下,地表后向散射回波隨頻率的變化關(guān)系與地表類型有關(guān),但當(dāng)入射角大于某一固定角度θc時(shí),其后向散射回波均隨頻率的增大而增大.文中各頻段下地表散射回波數(shù)據(jù)的預(yù)測(cè)值均與實(shí)測(cè)數(shù)據(jù)吻合良好.
表2 多頻地表經(jīng)驗(yàn)?zāi)Ｐ蛥?shù)的反演結(jié)果
圖3 裸土后向散射模型預(yù)測(cè)結(jié)果
圖4 水泥地后向散射模型預(yù)測(cè)結(jié)果
圖5 沙地后向散射模型預(yù)測(cè)結(jié)果
分別將上述各地表參數(shù)的反演結(jié)果帶回各地表散射經(jīng)驗(yàn)?zāi)Ｐ椭?確定其經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,并預(yù)測(cè)10 GHz入射波頻率下裸土,3.2 GHz入射波頻率下水泥地和沙地的后向散射回波值,與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析,其結(jié)果分別如圖6～圖8所示.對(duì)于裸土散射情況,H H極化下模型預(yù)測(cè)結(jié)果的均方根誤差為3.70 dB,相關(guān)系數(shù)為0.989;而VV極化下,其均方根誤差為3.05 dB,相關(guān)系數(shù)為0.984.水泥路面下,HH極化模型預(yù)測(cè)結(jié)果的均方根誤差為4.29 dB,相關(guān)系數(shù)為0.961;而VV極化下,其均方根誤差為4.40 dB,相關(guān)系數(shù)為0.942.同理,對(duì)入射波頻率為3.2 GHz沙地散射回波進(jìn)行分析,H H極化下模型預(yù)測(cè)結(jié)果的均方根誤差為4.63 dB,相關(guān)系數(shù)為0.977;而VV極化下,其均方根誤差為3.89 dB,相關(guān)系數(shù)為0.975.當(dāng)?shù)乇砩⑸浠夭〝?shù)據(jù)震蕩劇烈時(shí),其預(yù)測(cè)結(jié)果會(huì)稍微差些.文中裸土表面、水泥路面以及粗糙沙地散射回波的預(yù)測(cè)值與實(shí)測(cè)數(shù)據(jù)整體吻合良好,可以滿足實(shí)際工程需求,可有效減小多頻外場(chǎng)實(shí)驗(yàn)測(cè)量的工作量.
圖6 10 GHz裸土散射回波預(yù)測(cè)結(jié)果
圖7 3.2 GHz水泥地散射回波預(yù)測(cè)結(jié)果
圖8 3.2 GHz沙地散射回波預(yù)測(cè)結(jié)果
5 結(jié)束語(yǔ)
筆者將面散射地表、面-體散射地表的電磁散射理論模型進(jìn)行簡(jiǎn)化,建立了入射波頻率為1～16 GHz的單頻段和多頻段地表后向散射回波經(jīng)驗(yàn)?zāi)Ｐ?并采用L/S/X/Ku四頻段裸土、沙地和水泥地的后向散射實(shí)測(cè)數(shù)據(jù)對(duì)多頻散射經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行了驗(yàn)證.結(jié)果表明,文中提出的多頻經(jīng)驗(yàn)?zāi)Ｐ蜔o(wú)需測(cè)量地表的粗糙度和介電常數(shù),只需3個(gè)頻段或4個(gè)頻段的地表雷達(dá)回波數(shù)據(jù),即可確定多頻經(jīng)驗(yàn)?zāi)Ｐ偷木唧w形式,并實(shí)現(xiàn)對(duì)寬頻帶范圍內(nèi)任意頻率點(diǎn)下地表后向散射回波的預(yù)測(cè).該多頻地表電磁散射經(jīng)驗(yàn)?zāi)Ｐ秃?jiǎn)化了地表后向散射模型的復(fù)雜度,在保證地表后向散射系數(shù)預(yù)測(cè)精度的同時(shí),提高了計(jì)算效率,便于工程應(yīng)用.但值得注意的是,由于在不同頻段下,電磁波在地表的穿透深度不同,散射機(jī)制會(huì)產(chǎn)生差異,因此,更高頻段下大區(qū)域復(fù)雜地表的電磁散射經(jīng)驗(yàn)?zāi)Ｐ?將是下一步研究的重點(diǎn).
[1]SU Y,GUO Q,XUE B,et al.Spatial Distribution of Forest Aboveground Biomass in China:Estimation through Combination of Spaceborne Lidar,Optical Imagery,and Forest Inventory Data[J].Remote Sensing of Environment,2016,173:187-199.
[2]OZLEM K,GARCIA-RUBIA J M,NGHIA T,et al.Detection of Moving Human Micro-Doppler Signature in Forest Environments with Swaying Tree Components by Wind[J].Radio Science,2015,50:238-248.
[3]RASHIDI-RANJBAR E,DEH MOLLAIAN M.Target above Random Rough Surface Scattering Using a Parallelized IPO Accelerated by MLFMM[J].IEEE Geoscience and Remote Sensing Letters,2015,12:1481-1485.
[4]NASHASHIBI A,ULABY F T,SARABANDI K.Measurement and Modeling of the Millimeter-wave Backscatter Response of Soil Surfaces[J].IEEE Transactions on Geoscience and Remote Sensing,1996,34(2):561-572.
[5]OH Y,SARABANDI K,ULABY F T.An Empirical Model and an Inversion Technique for Radar Scattering from Bare Soil Surfaces[J].IEEE Transactions on Geoscience and Remote Sensing,1992,30(2):370-381.
[6]張連波,郭立新,茍雪銀,等.三層粗糙面電磁散射的矩量法研究[J].西安電子科技大學(xué)學(xué)報(bào),2013,40(6):147-154. ZHANG Lianbo,GUO Lixin,GOU Xueyin,et al.Method of Moment Investigation on Electromagnetic Scattering from the Three-layered Rough Interfaces[J].Journal of Xidian University,2013,40(6):147-154.
[7]FUNG A K.Microwave Scattering and Emission Models and Their Applications[M].Norwood:Artech House,1994.
[8]LI E S,SARABANDI K.Low Grazing Incidence Millimeter-wave Scattering Models and Measurements for Various Road Surfaces[J].IEEE Transactions on Antennas and Propagation,1999,47(5):851-861.
[9]NASHASHIBI A Y,SARABANDI K,AL-ZAID F A,et al.An Empirical Model of Volume Scattering from Dry Sandcovered Surfaces at Millimeter-wave Frequencies[J].IEEE Transactions on Geoscience and Remote Sensing,2013,51(6):3673-3682.
[10]鄭帆,郭立新,劉偉.風(fēng)積新月形沙丘的電磁散射特性分析[J].西安電子科技大學(xué)學(xué)報(bào),2013,40(6):132-139. ZHENG Fan,GUO Lixin,LIU Wei.Analysis of Electromagetic Scattering from the Wind-blown Barchan[J].Journal of Xidian University,2013,40(6):132-139.
[11]LIEVENS H,VERHOEST N E C.On the Retrieval of Soil Moisture in Wheat Fields from L-band SAR Based on Water Cloud Modeling,the IEM,and Effective Roughness Parameters[J].IEEE Geoscience and Remote Sensing Letters,2011,8(4):740-744.
[12]DUBOIS P C,van ZYL J,ENGMAN T.Measuring Soil Moisture with Imaging Radars[J].IEEE Transactions on Geoscience and Remote Sensing,1995,33(4):915-926.
[13]SHI J,WANG J,HSU A Y,et al.Estimation of Bare Soil Moisture and Surface Roughness Parameters Using L-band SAR Image Data[J].IEEE Transactions on Geoscience and Remote Sensing,1997,35(5):1254-1266.
[14]ZRIBI M,GORRAB A,BAGHDADI N.A New Soil Roughness Parameter for the Modelling of Radar Backscattering over Bare Soil[J].Remote Sensing of Environment,2014,152:62-73.
[15]ULABY F T,DONSON M C.Handbook of Radar Scattering Statistics for Terrain[M].Norwood:Artech House,1989.
[16]ULABY F T.Vegetation Clutter Model[J].IEEE Transactions on Antennas and Propagation,1980,28(4):538-545. [17]CHEN K L,CHEN K S,LI Z L,et al.Extension and Validation of an Advanced Integral Equation Model for Bistatic Scattering from Rough Surfaces[J].Progress In Electromagnetics Research,2015,152:59-76.
[18]HAFFNER S,PEREIRA L A,PEREIRA L F A.A Method for Optimization of Five-phase Induction Machines Based on Genetic Algorithms[J].Journal of Control Automation&Electrical Systems,2016,26(5):521-534.
(編輯:王 瑞)
Multi-band empirical model for backscattering coefficients prediction from typical land surfaces
ZHANG Yuanyuan1,WU Zhensen1,ZHANG Yushi2
(1.School of Physics and Optoelectronic Engineering,Xidian Univ.,Xi’an 710071,China; 2.China Research Institute of Radio Wave Propagation,Qingdao 266107,China)
Based on the electromagnetic scattering theory of the rough surface and discrete random medium,a multi-band empirical model for land backscattering prediction in 1～16 GHz is proposed to tap the potential of remote sensing in application for agricultural production and targets detection efficiently.It should be noticed that this empirical model only have seven unknowns,and that it overcomes the problems in numerous inputs and inaccuracies echoes prediction.It is verified by the measured data of soil,sand and concrete in the L/S/X/Ku band for HH and VV polarizations.The results show that the backscattering coefficients from soil,sand and concrete predicted by the empirical model are in good agreement with the measured data.Meanwhile,this multi-band empirical model can be used to forecast the backscattering echoes at any frequency in 1～16 GHz accurately.So this empirical model has a great application value for it can greatly improve the computational efficiency for the backscattering coefficients prediction,and also it can effectively reduce the workload in the field experiments.
multi-band;electromagnetic scattering;land;empirical model
TN959.71
A
1001-2400(2016)05-0190-07
10.3969/j.issn.1001-2400.2016.05.033
2015-12-01
國(guó)家自然科學(xué)基金資助項(xiàng)目(61571355)
張?jiān)?1986-),女,西安電子科技大學(xué)博士研究生,E-mail:zyy07063@163.com.
猜你喜歡
現(xiàn)代財(cái)經(jīng)-天津財(cái)經(jīng)大學(xué)學(xué)報(bào)(2022年5期)2022-06-01 06:08:32
振動(dòng)與沖擊(2022年6期)2022-03-27 12:18:26
軍事文摘(2021年18期)2021-12-02 01:27:58
內(nèi)蒙古林業(yè)(2021年6期)2021-06-26 10:42:32
幼兒智力世界(2021年4期)2021-05-17 17:10:29
科普童話·百科探秘(2018年7期)2018-09-13 11:12:08
建材與裝飾(2018年50期)2018-02-22 03:30:40
電子測(cè)試(2017年15期)2017-12-18 07:18:51
電源技術(shù)(2015年1期)2015-08-22 11:16:18
電子科技(2015年11期)2015-03-06 01:32:24
