專利名稱：一種基于起伏地形的sar森林二次散射有效路徑計(jì)算方法
技術(shù)領(lǐng)域：
本發(fā)明涉及屬于合成孔徑雷達(dá)森林遙感技術(shù)領(lǐng)域。
背景技術(shù)：
近年來(lái)，合成孔徑雷達(dá)SAR在林業(yè)遙感中得到了廣泛的應(yīng)用，成為定量提取生物量、平均樹高、垂直結(jié)構(gòu)等森林參數(shù)的有效方式。對(duì)于三維森林場(chǎng)景，對(duì)其地表進(jìn)行模擬或者依據(jù)實(shí)測(cè)數(shù)據(jù)可以獲得地面數(shù)字高程DEM圖，由合成孔徑雷達(dá)SAR照射三維森林場(chǎng)景時(shí)，由于樹干表面通常情況下較為粗糙，因此呈現(xiàn)出較強(qiáng)的漫散射特征，鏡面方向反射不明顯；而地表在遠(yuǎn)場(chǎng)效應(yīng)、大尺度條件下呈現(xiàn)較強(qiáng)的鏡面反射特性，非鏡面反射方向強(qiáng)度明顯較弱。上述樹干雙基地散射、地表鏡面反射的散射路徑為二次散射的主要成分，此即為SAR森林二次散射有效路徑。由此可以看出，地形起伏是影響森林參數(shù)定量反演精度的重要因素之一。地形起伏使雷達(dá)局部入射角發(fā)生改變，造成森林雷達(dá)散射機(jī)理發(fā)生變化，大大增加了森林參數(shù)反演的難度。只有深入了解地形對(duì)森林雷達(dá)散射的影響，并能夠在地形起伏情況下對(duì)SAR森林二次散射進(jìn)行準(zhǔn)確計(jì)算，才能夠正確地反演森林結(jié)構(gòu)參數(shù)。對(duì)于山區(qū)及丘陵地帶，真實(shí)的地表大都存在坡度的變化，此時(shí)，隨著地表坡度的變化，SAR森林二次散射有效路徑可能不存在也可能存在多條。而目前已有的SAR森林二次散射計(jì)算方法均基于水平地表或小傾斜的粗糙度表，僅針對(duì)地表的斜率固定、坡度通常不發(fā)生變化的情況適用，在這種情況下，SAR森林二次散射有效路徑最多僅有一條。因此采用目前已有的計(jì)算方法，對(duì)于山區(qū)或丘陵地帶進(jìn)行二次散射有效路徑的計(jì)算時(shí)也僅僅考慮最多僅有一條有效路徑的情況則會(huì)引入較大的誤差，導(dǎo)致計(jì)算結(jié)果與真實(shí)自然環(huán)境的特點(diǎn)不一致，嚴(yán)重影響了計(jì)算結(jié)果在SAR森林遙感信號(hào)模擬和結(jié)構(gòu)參數(shù)反演方面的應(yīng)用。

發(fā)明內(nèi)容
有鑒于此，本發(fā)明提供了一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，打破了現(xiàn)有方法僅考慮最多僅有一條有效路徑的情況僅針對(duì)水平地表或小傾斜的粗糙度表的局限性，降低了對(duì)山區(qū)或丘陵地帶進(jìn)行二次散射有效路徑的計(jì)算誤差。為達(dá)到上述目的，本發(fā)明的技術(shù)方案為:一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，包括如下步驟:步驟一、針對(duì)任一三維場(chǎng)景，獲得其地表的地面數(shù)字高程DEM圖，所述DEM圖由一定數(shù)量的小面元組成，對(duì)于每個(gè)小面元進(jìn)行均勻采樣獲取采樣點(diǎn)，選取一定范圍塊對(duì)所有采樣點(diǎn)分別沿方位向和地距向進(jìn)行分塊，將分于同一塊的采樣點(diǎn)進(jìn)行平面擬合，得到的擬合平面記為大面元，則地表被劃分 出大面元集合；選取范圍塊的依據(jù)是:對(duì)范圍塊內(nèi)的采樣點(diǎn)進(jìn)行擬合得到的大面元為矩形并且邊長(zhǎng)大于入射波波長(zhǎng)；對(duì)于擬合得到的大面元Au，其方位向坡度角為a 和地距向坡度角為a gr_.^其中i為方位向大面元序號(hào)，j為地距向大面元序號(hào)；步驟二、由合成孔徑雷達(dá)SAR向所述三維場(chǎng)景發(fā)射電磁波令A(yù)為電磁波的入射矢
量，為不電磁波的出射矢量，$與^s平行，選擇三維場(chǎng)景中一散射點(diǎn)作為散射介電粒子p，本步驟分為以下步驟:步驟201、根據(jù)a a」，J和a gr_.j計(jì)算Ai, j的法向量盡』;步驟202、根據(jù)w和^計(jì)算Ai,」?jié)M足鏡面反射的入射矢量.SJ步驟203、經(jīng)過(guò)點(diǎn)P作方向?yàn)?的直線與面元Ay.所在的平面相交于點(diǎn)Ru，則
'—U
Ri, j即為滿足P與面元A。有效二次散射的鏡面反射點(diǎn)；步驟204、判斷若Ri,」在Ai,」內(nèi)部，則P與Ai,」存在二次散射有效路徑，否則P與Ai；J不存在二次散射有效路徑；步驟205、若P與面元A。存在二次散射有效路徑，則經(jīng)過(guò)點(diǎn)P作方向?yàn)?^;的直線
與Ai,j所在的平面相交于點(diǎn)Eu，Eu即為P與Ai,j有效二次散射的等效散射相位中心；步驟206、將步驟一得到的大面元均針對(duì)P進(jìn)行步驟201 步驟205的處理，即可得到P的所有二次散射有效路徑及相應(yīng)的等效相位中心；步驟三、三維場(chǎng)景中所有散射`點(diǎn)均進(jìn)行步驟二的處理，即可得到三維場(chǎng)景中任一散射點(diǎn)的所有二次散射有效路徑及相應(yīng)的等效相位中心。進(jìn)一步地，本方案中步驟206為:步驟2061、在步驟一劃分的大面元中選擇散射介電粒子在地表投影點(diǎn)所對(duì)應(yīng)的地距向面元；步驟2062、在步驟2061獲得的地距向面元中選擇比散射介電粒子在地表投影點(diǎn)更接近雷達(dá)一側(cè)的面元；步驟2063、在步驟2062獲得的面元中選擇散射介電粒子在地表投影點(diǎn)設(shè)定范圍內(nèi)的面元；步驟2064、將步驟2063獲得的面元進(jìn)行步驟201 步驟205的處理。進(jìn)一步地，本方案中步驟三為:步驟31、將三維森林場(chǎng)景均勻劃分為三維網(wǎng)格，所有的散射介電粒子根據(jù)空間位置分布在對(duì)應(yīng)的網(wǎng)格內(nèi)，選取網(wǎng)格的最大邊長(zhǎng)小于SAR 二維分辨率；步驟32、按照步驟二計(jì)算每個(gè)網(wǎng)格中心處散射介電粒子的二次散射有效路徑，建立網(wǎng)格中心與二次散射有效路徑的關(guān)系；步驟33、針對(duì)任一散射介電粒子Q，獲取Q所在網(wǎng)格中心的二次散射有效路徑作為Q的二次散射有效路徑。有益效果:1、本發(fā)明深入了解地形對(duì)森林雷達(dá)散射的影響，充分考慮了地表具有大尺度起伏特征時(shí)森林雷達(dá)二次散射可能不存在或存在多條二次散射有效路徑的情況，能夠有效地表征地形起伏情況下森林的雷達(dá)二次散射特征，從而打破了現(xiàn)有方法僅針對(duì)水平地表或小傾斜的粗糙度表的局限性，使得可為山區(qū)、丘陵地帶森林的SAR遙感數(shù)據(jù)模擬以及森林參數(shù)反演算法研究提供支撐。2、本發(fā)明考慮到算法中計(jì)算量較大的問(wèn)題，結(jié)合實(shí)際情況，提出了有效的快速算法，在保證算法精度的基礎(chǔ)上大大減少了計(jì)算量。


圖1為本發(fā)明實(shí)施例中機(jī)載激光雷達(dá)實(shí)測(cè)的徂徠山DEM圖；圖2為本發(fā)明實(shí)施例中步驟I)擬合得到的大尺度面元DEM圖；圖3為本發(fā)明實(shí)施例中二次散射有效路徑計(jì)算示意圖；圖4為本發(fā)明實(shí)施例中二次散射有效路徑快速計(jì)算示意；圖5為本發(fā)明實(shí)施例中計(jì)算結(jié)果示意圖，Ca)為闊葉樹的幾何示意圖，(b)為采用基于單一坡度地表模型的計(jì)算結(jié)果，(c)為采用本發(fā)明的計(jì)算方法的計(jì)算結(jié)果。
具體實(shí)施方式
下面結(jié)合附圖并舉實(shí)施例，對(duì)本發(fā)明進(jìn)行詳細(xì)描述。本實(shí)施例采用機(jī)載激光雷達(dá)獲取的山東省徂徠山地區(qū)的地面數(shù)字高程DEM圖作為起伏地形實(shí)例，如圖1所示。一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，步驟如下:步驟一、將地表地形沿方位向和地距向均勻劃分為大面元，大面元的要求為:大面兀的邊長(zhǎng)大于入射波波長(zhǎng)。其中所獲得的地表的地面數(shù)字高程DEM圖由一定數(shù)量的小面元組成，其每個(gè)面元可以看成具有一定坡度(二維)的無(wú)限大介電平面。由于真實(shí)地表通常疊加了小尺度的粗糙起伏，為了避免小尺度粗糙的影響，所以本發(fā)明使用平面擬合的方式將地表沿方位向和地距向均勻劃分出大面元，具體為:在小面元的DEM圖中，對(duì)于每個(gè)小面元進(jìn)行均勻采樣獲取采樣點(diǎn)，選取一定范圍塊，對(duì)所有采樣點(diǎn)沿方位向和地距向以該范圍塊進(jìn)行分塊，將劃分于同一塊的采樣點(diǎn)進(jìn)行擬合，得到的擬合平面記為大面元，則地表被劃分出大面元集合。選取范圍塊的依據(jù)是:對(duì)范圍塊內(nèi)的采樣點(diǎn)進(jìn)行擬合得到的大面元為矩形并且邊長(zhǎng)大于入射波波長(zhǎng)。對(duì)于擬合得到的大面元其方位向坡度角為a 和地距向坡度角為a gr_.^其中i為方位向大面元序號(hào)，j為地距向大面元序號(hào)。其中方位向?yàn)槔走_(dá)運(yùn)動(dòng)方向，地距向?yàn)殡姶挪▊鞑シ较虻乃椒至?。本步驟中擬合得到的大面元DEM如圖2所示。步驟二、使用合成孔徑雷達(dá)SAR發(fā)射電磁波對(duì)森林場(chǎng)景進(jìn)行散射模擬，在步驟I擬合得到的大平面下，二次散射有效路徑由入射角度、大尺度面元的坡度坡向和中心位置，以及散射介電粒子空間位置的共同決定。Jci和A分別表不SAR的電磁波入射矢量和出射矢量,P點(diǎn)為森林場(chǎng)景中任一散射介電粒子，則求解P與A。之間二次散射有效路徑的具體步驟為:步驟201、根據(jù)a a_.」和a gr_.」計(jì)算ku」的法向量"假設(shè)建立如圖3所示的空間坐標(biāo)系，以豎直向上的方向?yàn)閆軸，以水平面為xoy面，其中方位向?yàn)閅軸，地距向?yàn)閄軸，貝1J對(duì)于來(lái)說(shuō)，a a i；J為Ai,」面與y軸夾角，a gr i；J為A。面與x軸夾角，則根據(jù)A^j邊長(zhǎng)
與角度的關(guān)系，可建立Ay邊長(zhǎng)的矢量，兩相鄰邊長(zhǎng)矢量作叉乘即得A。的法向量^。步驟202、由法向量和電磁波出射矢量(計(jì)算Auj滿足鏡面反射的入射矢量￡ 和I為鏡面反射，關(guān)于對(duì)稱且在A ,和組成的面上。
r — l，J， J _lJ Sl，J1 -1J 1J步驟203、由于地表在遠(yuǎn)場(chǎng)效應(yīng)、大尺度條件下呈現(xiàn)較強(qiáng)的鏡面反射特性，非鏡面反射方向強(qiáng)度明顯較弱；因此對(duì)于點(diǎn)P，其與A。的有效二次散射的反射點(diǎn)即為鏡面反射點(diǎn)，求該鏡面反射點(diǎn)的具體方法為:經(jīng)過(guò)點(diǎn)P作方向?yàn)?(.的直線與面元A。所在的平面相交于點(diǎn)Ru，則Ru即為滿足散射介電粒子P與面元A。有效二次散射的鏡面反射點(diǎn)；步驟204、由于散射介電粒子P與面元Ay可能并不存在二次散射有效路徑，因此需要作出如下判斷:判斷Ru是否在面元A。內(nèi)部，若在內(nèi)部，則散射介電粒子P與面元A。存在二次散射有效路徑，否則不存在；步驟205、若存在二次散射有效路徑，則經(jīng)過(guò)點(diǎn)P作方向?yàn)榈闹本€與面元Ay所
在的平面相交于點(diǎn)Eu，則Eu即為P與Ai, j有效二次散射的等效散射相位中心。步驟206、針對(duì)步驟I擬合得到的大尺度平面的每個(gè)面元進(jìn)行上述步驟的處理，即可得到P的在起伏地表下的所有二次散射有效路徑及相應(yīng)的等效相位中心。步驟三、所述三維場(chǎng)景中所有散射點(diǎn)均進(jìn)行步驟二的處理，即可得到三維場(chǎng)景中任一散射點(diǎn)的所有二次散射有效路徑及相應(yīng)的等效相位中心。針對(duì)三維森林場(chǎng)景中所有的散射介電粒子進(jìn)行上述步驟的處理和判斷，則可以得到整個(gè)森林場(chǎng)景的二次散射有效路徑和對(duì)應(yīng)的等效散射中心；但是，當(dāng)場(chǎng)景較大時(shí)，對(duì)所有面元進(jìn)行遍歷判斷的計(jì)算量較大，因此對(duì)于步驟206本發(fā)明采用如下方法限定所處理面元的范圍:步驟2061、由于散射介電粒子與地表發(fā)生鏡面反射主要集中在沿著電磁波入射的方向上，即地距向，因此僅選擇散射介電粒子在地表投影點(diǎn)所對(duì)應(yīng)的地距向面元；在本實(shí)施例中，假設(shè)點(diǎn)P在地表投影點(diǎn)位于大面元A2，3，則選取A2，3所對(duì)應(yīng)的地距向面元，即“
△2, 3°步驟2062、由于散射介電粒子不可能遠(yuǎn)距方向的地表面元發(fā)生二次散射，因此在步驟2061選擇的地距向面元中選擇比散射介電粒子在地表投影點(diǎn)更接近雷達(dá)一側(cè)的面
元；在本實(shí)施例中，由SAR的電磁波入射矢量X可以看出，SAR應(yīng)位于點(diǎn)P的左側(cè)，則在步驟
2061所選擇的中k2A、k2,2為更接近雷達(dá)一側(cè)的面元。步驟2063、進(jìn)行小范圍搜索，由于 樹冠對(duì)電磁波產(chǎn)生一定衰減作用，因此即使距散射介電粒子較遠(yuǎn)處的面元存在二次散射有效路徑，二次散射波返回雷達(dá)的能量也相對(duì)較小，因此在步驟2062選擇的更接近雷達(dá)一側(cè)的面元中選擇散射介電粒子在地表投影點(diǎn)小范圍內(nèi)的面元；在本實(shí)施例中，可對(duì)該小范圍進(jìn)行設(shè)定，并進(jìn)行實(shí)驗(yàn)，選擇實(shí)驗(yàn)結(jié)果最佳的設(shè)定值作為該小范圍的設(shè)定值。在本實(shí)施例中假設(shè)該小范圍的設(shè)定值恰好僅包含與其最接近的A2，2面，則在本次二次散射有效路徑的計(jì)算中，僅計(jì)算P與面元A2，2的二次散射有效路徑即可，由此大大減少了計(jì)算量。在進(jìn)行二次散射有效路徑的計(jì)算時(shí)，選擇森林場(chǎng)景中的散射點(diǎn)作為散射介電粒子，則三維森林場(chǎng)景中每株樹木均由成千上萬(wàn)個(gè)的散射介電粒子組成，整個(gè)森林場(chǎng)景的散射粒子數(shù)量可能達(dá)到百萬(wàn)甚至千萬(wàn)的量級(jí)。因此在步驟三中，對(duì)于數(shù)量巨大的散射介電粒子，每個(gè)散射介電粒子均使用本實(shí)施例中提供的步驟二進(jìn)行計(jì)算，為減少計(jì)算量，本實(shí)施例提供如下方法，采用劃分三維網(wǎng)格的方式，實(shí)現(xiàn)二次散射有效路徑的快速計(jì)算，具體為:步驟31、將三維森林場(chǎng)景均勻劃分為三維正方體網(wǎng)格，所有的散射介電粒子根據(jù)空間位置分布在對(duì)應(yīng)的網(wǎng)格內(nèi)，選取網(wǎng)格的尺寸小于SAR 二維分辨率，如圖4所示。步驟32、計(jì)算每個(gè)網(wǎng)格中心處散射介電粒子的二次散射有效路徑，建立網(wǎng)格中心與二次散射有效路徑的關(guān)系。步驟33、針對(duì)任一散射介電粒子Q，獲取Q所在網(wǎng)格中心的二次散射有效路徑作為Q的二次散射有效路徑。針對(duì)三維森林場(chǎng)景中所有散射介電粒子均以其所在網(wǎng)格的網(wǎng)格中心的二次散射有效路徑作為其二次散射有效路徑，由此獲取整個(gè)三維森林場(chǎng)景的二次散射有效路徑。本發(fā)明針對(duì)在步驟I擬合得到的徂徠山區(qū)大尺度面元DEM圖上生長(zhǎng)一株闊葉樹，樹高20m，闊葉樹被離散化為1000個(gè)散射介電粒子，分別采用基于單一坡度地表模型的計(jì)算方法和本發(fā)明的計(jì)算方法計(jì)算闊葉樹的有效二次散射情況，計(jì)算結(jié)果如圖5 (a) (b) (c)所示，其中，基于單一坡度地表模型的計(jì)算方法計(jì)算得到的二次散射有效路徑數(shù)量為1000條，即每個(gè)散射介電粒子存在且僅存在I條，而且其等效散射相位中心分布在單一坡度的地表上；采用本發(fā)明的方法，二次散射有效路徑共計(jì)算得到2171條，平均每個(gè)散射介電粒子存在2條，而且其等效散射中心分布范圍較大，與地表面的最大高度差接近5m，與真實(shí)地形特性相符，真實(shí)的山地、丘陵等地形中大尺度起伏特征是普遍存在的，采用本發(fā)明的方法計(jì)算這些地區(qū)森林的SAR 二次散射有效路徑更加逼近實(shí)際情況，這說(shuō)明，按照本發(fā)明提供的技術(shù)方案進(jìn)行起伏地形下的SAR森林二次散射有效路徑計(jì)算可以達(dá)到預(yù)期目的，可為山區(qū)、丘陵地帶森林的SAR遙感數(shù)據(jù)模擬以及森林參數(shù)反演算法研究提供支撐。綜上所述，以上僅為本發(fā)明的較佳實(shí)施例而已，并非用于限定本發(fā)明的保護(hù)范圍。凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
權(quán)利要求
1.一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，其特征在于，包括如下步驟: 步驟一、針對(duì)地形存在起伏特征的三維森林場(chǎng)景，獲得其地表的數(shù)字高程DEM圖，所述DEM圖由一定數(shù)量的小面元組成，對(duì)于每個(gè)小面元進(jìn)行均勻采樣獲取采樣點(diǎn)，選取一定范圍塊對(duì)所有采樣點(diǎn)分別沿方位向和地距向進(jìn)行分塊，將分于同一塊的采樣點(diǎn)進(jìn)行擬合，得到的擬合平面記為大面元，則地表被劃分出大面元集合； 選取范圍塊的依據(jù)是:對(duì)范圍塊內(nèi)的采樣點(diǎn)進(jìn)行擬合得到的大面元為矩形并且邊長(zhǎng)大于入射波波長(zhǎng)； 對(duì)于擬合得到的大面元其方位向坡度角為Cia u和地距向坡度角為a ，其中i為方位向大面元序號(hào)，j為地距向大面元序號(hào)； 步驟二、由合成孔徑雷達(dá)SAR向所述三維森林場(chǎng)景發(fā)射電磁波，令&為電磁波的入射矢量，為不電磁波的出射矢量，(1J/f'平行且方向相反，選擇三維森林場(chǎng)景中任一散射點(diǎn)作為散射介電粒子P，本步驟分為以下步驟:步驟201、根據(jù)a a Uj和a gT Uj計(jì)算Auj的法向量// 步驟202、根據(jù)見.和I計(jì)算A。滿足鏡面反射的入射矢量& SI 步驟203、經(jīng)過(guò)點(diǎn)P作方向?yàn)镠r t ;的直線與面元Ai j所在的平面相交于點(diǎn)Ri」，則Ri jr」，j，，，即為滿足P與面元A。有效二次散射的鏡面反射點(diǎn)； 步驟204、判斷若Ru在Alj內(nèi)部，則P與A。存在二次散射有效路徑，否則P與Auj不存在二次散射有效路徑； 步驟205、若P與面元A。存在二次散射有效路徑，則經(jīng)過(guò)點(diǎn)P作方向?yàn)榈闹本€與Ai, j所在的平面相交于點(diǎn)Eu，Ei, j即為P與A。有效二次散射的等效散射相位中心； 步驟206、將步驟一得到的大面元均針對(duì)P進(jìn)行步驟201 步驟205的處理，即可得到P的所有二次散射有效路徑及相應(yīng)的等效相位中心； 步驟三、所述三維森林場(chǎng)景中所有散射點(diǎn)均進(jìn)行步驟二的處理，即可得到三維森林場(chǎng)景中任一散射點(diǎn)的所有二次散射有效路徑及相應(yīng)的等效相位中心。
2.如權(quán)利要求1所述的一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，其特征在于，所述步驟206為: 步驟2061、在步驟一劃分的大面元中選擇散射介電粒子在地表投影點(diǎn)所對(duì)應(yīng)的地距向面元； 步驟2062、在步驟2061獲得的地距向面元中選擇比散射介電粒子在地表投影點(diǎn)更接近雷達(dá)一側(cè)的面元； 步驟2063、在步驟2062獲得的面元中選擇散射介電粒子在地表投影點(diǎn)設(shè)定范圍內(nèi)的面元； 步驟2064、將步驟2063獲得的面元進(jìn)行步驟201 步驟205的處理。
3.如權(quán)利要求1所述的一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，其特征在于，所述步驟三為: 步驟31、將三維森林場(chǎng)景均勻劃分為三維網(wǎng)格，三維森林場(chǎng)景中所有的散射介電粒子根據(jù)空間位置分布在對(duì)應(yīng)的網(wǎng)格內(nèi)，選取網(wǎng)格的最大邊長(zhǎng)小于SAR 二維分辨率； 步驟32、按照步驟二計(jì)算每個(gè)網(wǎng)格中心處散射介電粒子的二次散射有效路徑，建立網(wǎng)格中心與二次散射有效路徑的關(guān)系； 步驟33、針對(duì)三維森林場(chǎng)景中任一散射介電粒子Q，獲取Q所在網(wǎng)格中心的二次散射有效路徑作為Q的二次散射 有效路徑。
全文摘要
本發(fā)明公開了一種基于起伏地形的SAR森林二次散射有效路徑計(jì)算方法，屬于合成孔徑雷達(dá)森林遙感技術(shù)領(lǐng)域。步驟如下一、將地表沿方位向和地距向均勻劃分為一定尺度的大尺度面元；二、對(duì)于大尺度面元Ai,j，計(jì)算Ai,j的法向量根據(jù)和電磁波出射矢量計(jì)算Ai,j滿足鏡面反射的入射矢量令點(diǎn)P(xP,yP,zP)為任一散射介電粒子的中心坐標(biāo)，過(guò)P作方向?yàn)榈闹本€與Ai,j所在平面相交于點(diǎn)Ri,j；判斷Ri,j若在Ai,j內(nèi)部，則P與Ai,j存在二次散射有效路徑，若存在，則經(jīng)P作方向?yàn)榈闹本€與Ai,j所在的平面相交于點(diǎn)Ei,j；以步驟二遍歷或者以一定的順序遍歷步驟一擬合得到的大尺度面元。本方法適用于地形起伏情況下的SAR森林場(chǎng)景。
文檔編號(hào)G01S13/90GK103116161SQ20131003048
公開日2013年5月22日 申請(qǐng)日期2013年1月25日 優(yōu)先權(quán)日2013年1月25日
發(fā)明者胡程, 孫晗偉, 曾濤, 龍騰, 曾大治 申請(qǐng)人:北京理工大學(xué)