发明名称 一种基于起伏地形的SAR森林二次散射有效路径计算方法
摘要 本发明公开了一种基于起伏地形的SAR森林二次散射有效路径计算方法,属于合成孔径雷达森林遥感技术领域。步骤如下:一、将地表沿方位向和地距向均匀划分为一定尺度的大尺度面元;二、对于大尺度面元A<sub>i,j</sub>,计算A<sub>i,j</sub>的法向量<img file="DDA00002779512400011.GIF" wi="98" he="105" />根据<img file="DDA00002779512400012.GIF" wi="74" he="85" />和电磁波出射矢量<img file="DDA00002779512400013.GIF" wi="40" he="92" />计算A<sub>i,j</sub>满足镜面反射的入射矢量<img file="DDA00002779512400014.GIF" wi="141" he="112" />令点P(x<sub>P</sub>,y<sub>P</sub>,z<sub>P</sub>)为任一散射介电粒子的中心坐标,过P作方向为<img file="DDA00002779512400015.GIF" wi="122" he="96" />的直线与A<sub>i,j</sub>所在平面相交于点R<sub>i,j</sub>;判断R<sub>i,j</sub>若在A<sub>i,j</sub>内部,则P与A<sub>i,j</sub>存在二次散射有效路径,若存在,则经P作方向为<img file="DDA00002779512400016.GIF" wi="82" he="85" />的直线与A<sub>i,j</sub>所在的平面相交于点E<sub>i,j</sub>;以步骤二遍历或者以一定的顺序遍历步骤一拟合得到的大尺度面元。本方法适用于地形起伏情况下的SAR森林场景。
申请公布号 CN103116161A 申请公布日期 2013.05.22
申请号 CN201310030483.5 申请日期 2013.01.25
申请人 北京理工大学 发明人 胡程;孙晗伟;曾涛;龙腾;曾大治
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 付雷杰;高燕燕
主权项 1.一种基于起伏地形的SAR森林二次散射有效路径计算方法,其特征在于,包括如下步骤:步骤一、针对地形存在起伏特征的三维森林场景,获得其地表的数字高程DEM图,所述DEM图由一定数量的小面元组成,对于每个小面元进行均匀采样获取采样点,选取一定范围块对所有采样点分别沿方位向和地距向进行分块,将分于同一块的采样点进行拟合,得到的拟合平面记为大面元,则地表被划分出大面元集合;选取范围块的依据是:对范围块内的采样点进行拟合得到的大面元为矩形并且边长大于入射波波长;对于拟合得到的大面元A<sub>i,j</sub>,其方位向坡度角为α<sub>a_i,j</sub>和地距向坡度角为α<sub>gr_i,j</sub>,其中i为方位向大面元序号,j为地距向大面元序号;步骤二、由合成孔径雷达SAR向所述三维森林场景发射电磁波,令<img file="FDA00002779512100011.GIF" wi="48" he="96" />为电磁波的入射矢量,<img file="FDA00002779512100012.GIF" wi="51" he="88" />为示电磁波的出射矢量,<img file="FDA00002779512100013.GIF" wi="44" he="84" />与<img file="FDA00002779512100014.GIF" wi="48" he="91" />平行且方向相反,选择三维森林场景中任一散射点作为散射介电粒子P,本步骤分为以下步骤:步骤201、根据α<sub>a_i,j</sub>和α<sub>gr_i,j</sub>计算A<sub>i,j</sub>的法向量<img file="FDA00002779512100015.GIF" wi="99" he="85" />步骤202、根据<img file="FDA00002779512100016.GIF" wi="72" he="85" />和<img file="FDA00002779512100017.GIF" wi="48" he="92" />计算A<sub>i,j</sub>满足镜面反射的入射矢量<img file="FDA00002779512100018.GIF" wi="140" he="118" />步骤203、经过点P作方向为<img file="FDA00002779512100019.GIF" wi="117" he="100" />的直线与面元A<sub>i,j</sub>所在的平面相交于点R<sub>i,j</sub>,则R<sub>i,j</sub>即为满足P与面元A<sub>i,j</sub>有效二次散射的镜面反射点;步骤204、判断若R<sub>i,j</sub>在A<sub>i,j</sub>内部,则P与A<sub>i,j</sub>存在二次散射有效路径,否则P与A<sub>i,j</sub>不存在二次散射有效路径;步骤205、若P与面元A<sub>i,j</sub>存在二次散射有效路径,则经过点P作方向为<img file="FDA000027795121000110.GIF" wi="83" he="83" />的直线与A<sub>i,j</sub>所在的平面相交于点E<sub>i,j</sub>,E<sub>i,j</sub>即为P与A<sub>i,j</sub>有效二次散射的等效散射相位中心;步骤206、将步骤一得到的大面元均针对P进行步骤201~步骤205的处理,即可得到P的所有二次散射有效路径及相应的等效相位中心;步骤三、所述三维森林场景中所有散射点均进行步骤二的处理,即可得到三维森林场景中任一散射点的所有二次散射有效路径及相应的等效相位中心。
地址 100081 北京市海淀区中关村南大街5号