基于分形结构的植被高阶相干散射模型研究-饶丽婷.pdf
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《基于分形结构的植被高阶相干散射模型研究-饶丽婷.pdf》由会员分享,可在线阅读,更多相关《基于分形结构的植被高阶相干散射模型研究-饶丽婷.pdf(7页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第38卷第10期 电 子 与 信 息 学 报 Vol. 38No.10 2016年10月 Journal of Electronics & Information Technology . Oct. 2016 基于分形结构的植被高阶相干散射模型研究 饶丽婷张晓娟*王友成方广有(中国科学院电磁辐射与探测技术重点实验室 北京 100190) (中国科学院大学 北京 100049) 摘 要:电磁波低频入射情况下,植被中散射体独立不相关的假设无效,此时应考虑散射体之间的相干效应以及近场互作用。该文提出一种基于分形结构的植被高阶相干散射模型,该模型利用分形理论生成近乎真实植被的3维几何结构,根据每个散射
2、体的空间位置信息考虑了相干效应,应用互易定理计算了相邻散射体间高阶互作用,结合非相干的分层模型中后向散射机制划分方式,给出了各项散射机制的表达式。与机载合成孔径雷达实验数据对比,验证了模型的准确性。在针叶林仿真参数下,分析了各项散射机制对总散射效应的贡献与入射频率、角度、植被结构的关系,结果表明,低频入射条件下,稀疏植被散射模型可进一步简化从而应用于参数反演中。 关键词:植被;分形结构;相干效应;近场互作用;后向散射 中图分类号:TP722.6 文献标识码:A 文章编号:1009-5896(2016)10-2502-07 DOI: 10.11999/JEIT160095 High-order
3、Coherent Scattering Model for Vegetation with Fractal Structures RAO LitingZHANG XiaojuanWANG YouchengFANG Guangyou(Key Laboratory of Electromagnetic Radiation and Sensing Technology, Chinese Academy of Sciences, Beijing 100190, China) (University of Chinese Academy of Sciences, Beijing 100049, Chin
4、a) Abstract: At low frequency, the assumption of independent scattering of the scatterers in vegetation medium is no longer valid. The coherent effect and near field interactions should be considered. In this paper, a high-order coherent scattering model for vegetation with fractal structure is pres
5、ented. The fractal theory is employed to generate a realistic 3-D spatial structure of vegetation. The near field interaction between scatterers is formulated using an efficient algorithm based on the reciprocity theorem. For the coherent effect, every scatterer with a deterministic location is take
6、n into account. The main scattering mechanisms are defined in the way of layered vegetation model, allowing better understanding of microwave interaction with trunk-crown structure. Good agreements are obtained from the comparisons of the theoretical predictions with the multifrequency and multipola
7、rization measurement results of boreal forest. Through an extensive ground truth, theoretical analysis of the contribution of the scattering mechanisms for various frequencies, incident angles and vegetation structures is carried out. It is found that under specified conditions the vegetation scatte
8、ring model can be simplified according to the main contribution scattering mechanism which can be applied to the inversion issue. Key words: Vegetation; Fractal structures; Coherent effect; Near field interaction; Backscattering 1 引言地物电磁散射模型是电磁学领域中一个热门课题,在雷达探测、无线电传播与通信以及农业科学等方面有着广泛的应用1 6。通常情况下,从微波遥感
9、观测数据中获得地物的特性参数需要首先建立收稿日期:2016-01-21;改回日期:2016-07-01;网络出版:2016-09-01 *通信作者:张晓娟 基金项目:国家自然科学基金( 61172017) Foundation Item: The National Natural Science Foundation of China (61172017) 一个相应的散射模型,然后应用散射模型反演未知参数7 14,其中,针对地表植被的散射模型,主要是以辐射传输理论和修正的Born近似方法为基础的离散植被模型12 16。离散植被模型将植被模拟成由随机分布的树干、树枝、树叶、果实等离散散射体组成
10、的介质,这些离散散射体可由形状规则的几何体比如有限长圆柱、小圆盘、椭球体等来描述17。 早期的植被散射模型,在建模过程中忽略植被散射体相互作用,如文献12提出的MIMICS模型以及在此基础上的改进模型13,15。对于低频入射情万方数据第10期 饶丽婷等: 基于分形结构的植被高阶相干散射模型研究 2503 况下,植被间距与波长尺寸相当或者更小,此时植被属于电稠密植被,将植被散射体当作是独立不相关的假设不再有效,需考虑植散射体之间相干效应以及近场互作用。文献18提出了基于矢量辐射传输理论的相位幅度修正模型,该方法通过天线阵列来考虑散射体之间的相干效应和近场相互作用,只能修正相同形状和大小的散射体,
11、对于不同种类散射体之间的相干效应及近场互作用无法考虑。文献19和文献20相继采用分形树研究植被的相干散射模型,这些模型中仅考虑了相干效应,散射机制按照体散射和面散射的角度区分。 因此,针对电稠密植被,本文建立了一种植被的高阶相干散射模型,该模型采用分形结构实现了植被的几何建模。结合非相干的分层模型中散射机制划分方式,考虑了冠层内散射体直接后向散射、冠层散射体地面一次作用后向散射、树干地面一次作用后向散射、地面直接后向散射等4项散射机制,各项散射机制中根据散射体的几何位置信息准确考虑了相干效应及近场互作用即高阶散射。与合成孔径雷达实测数据对比验证了模型的准确性,而后根据针叶林仿真参数,讨论并分析
12、了入射角、频率以及植被结构对各项散射分量的影响。 2 模型原理 2.1 植被几何结构 由于植被的形状复杂并具有随机性,对其进行精确的几何建模相当困难,在电磁散射模型中,通常将植被离散成由有限数目的不同形状规则体构成,其中树干被模拟成近乎垂直的有限长圆柱体,树枝被模拟成不同取向及大小的圆柱体,对于阔叶林,树叶被模拟成小圆盘,对于针叶林,树叶被模拟成针状体。然后采用一定的几何组合方式将离散规则体构成植被,其中,一种应用广泛的方式是采用分层建模12 13,15,植被分成冠层、树干层,在相应的层内离散规则体根据密度均匀分布,但是这种几何建模方式忽略了植被结构中各个部分之间的相互干涉,在低频入射时各部分
13、干涉会给总散射带来较大影响,由几何结构不同产生的相干效应影响很大。 分形结构考虑了自相似的多尺度几何构造,可产生与实际树木近乎一致的结构。L -系统是丹麦植物学家Lindenmayer首次提出的描述树木生长的数学模型21,它用形式语言来描述植物的拓扑结构,经过有限次数的规则替换,初始符号被扩展成一个字符序列,通过对字符序列进行数学构图解释,从而产生复杂的植物形态。本文采用L-系统方法随机生成3维植被,计算机实现分两个过程完成,方法描述如下 第1个过程:根据起始符号元和替换规则产生目标字符串。 (1)声明并设置公理、产生式规则; (2)声明并设置起始点、初始角、迭代步长以及迭代深度等参数; (3
14、)根据产生式重写规则,将初始形式的每个字符依次替换成新的字符形式; (4)反复步骤(3 ),直至到达迭代深度为止,最终生成目标字符串。 第2个过程:解释目标字符串并绘图。 (1)读取目标字符串中的每个字符; (2)根据读到的每个字符执行相应的绘图操作。 2.2 各项散射机制 图1是全局坐标系下植被散射模型示意图。假设入射到植被上的电磁波是沿着ik方向传播的平面波,表示为 0j0eikiikrEE(1) 式中,0k为数将所构造的分形树离散成具有特定位置、取向、几何尺寸的散射体,散射体包含表示树干和树枝的圆柱体、表示叶子的针状体或圆盘。为了使各项散射机制易于分析,借用分层散射模型中散射机制的划分方
15、式15,如图1所示。在远区的总后向散射场可以表示为 0j011eikNMb bg tg g is nn mnm krE S S S SEr(2) 其中,N是冠层中总散射体个数,nS表示冠层中第n个散射体的散射矩阵,M是总树干个数,mS表示第m个树干的散射矩阵。式(2)中包含4项散射机制:冠层内散射体直接后向散射bnS,冠层内散射体地面一次作用后向散射bgnS,树干地面一次作用后向散射tgmS,地面直接后向散射gS。植被散射模型中包含的散射机制示意图如图2所示。 图1 全局坐标系下植被散射模型示意图 万方数据2504 电 子 与 信 息 学 报 第38卷 图2各项散射机制示意图 以上各项散射机制
16、具体可以表示为 1j0e,di sn nn s i nS TS k k T (3) 23j0j0e,bg i r sn n n gi i nir sn n s gs nS T S k k T RTT RT S k k T (4) 23j0j0e,tg i r sm m m gi i mir sm m s gs mS T S k k T RTT RT S k k T (5) 4j0e,gig siS TS k k T (6) 其中,iT表示入射方向传播矩阵,sT表示散射方向传播矩阵,rT表示反射方向传播矩阵,R表示地面反射矩阵,0S表示不经历衰减的散射矩阵,为散射体在全局坐标系中相移,各项散射机
17、制中相移分别为 10nisk k kr (7) 20 02nnb i s g gikk k kr r n nk (8) 30 02nnb i s g gs k kr r n nk (9) 20 02mmt i s g gikk k kr r n nk (10) 30 02mmt i s g gs k kr r n nk (11) 40gisk k kr (12) 平面电磁波在植被介质中因为散射和吸收会衰减,根据Foldy近似16,将这种衰减效应采用传输矩阵(, )sTk表示为 e 0,0evvhhiM siM ssTk (13) 其中, 00(2 / ) , , ,pq pq i i s sM
18、 kn f , 0k是入射波自由空间中的波数,0n表示散射体的体密度,,pq分别表示入射波和散射波的极化方式,(, ,pq i if ,)ss是散射体的系统平均前向散射矩阵,s表示电磁波穿透植被层的经历路径的长度。 考虑植被分布在垂直方向上的不均匀性,计算传输矩阵时,将植被分成M层,每层厚度为md,1, 2, ,mM ,电磁波沿着传播方向在第m层中经历路径为mL,如图3所示,假定分层边界之间无反射或绕射,在第m层中第n个散射体经历的传播矩阵分别为 图3 电磁波在植被中传播模型 11i ii i in m mn m m M MLL LTT T T (14) 11s ss s sn m mn m
19、m M MLL LTT T T (15) 1 1 11gi gi gi gi gin m mn m mLL LTT T T (16) 1 1 11gs gs gs gs gsn m mn m mLL LTT T T (17) 其中,mnL表示散射体n沿着传播方向在第m层中经历路径。 本文中地面设定为小尺度起伏的粗糙面,用一阶小扰动模型来计算地面的散射矩阵。小扰动模型要求粗糙面满足条件:01kH , 3201kHL , 2 / 0.3HL,其中H , L为粗糙面的均方根高度和相关长度,一阶小扰动模型在粗糙面散射特性研究中非常具有代表性,在很多正向和逆向问题中被采用18。散射机制bgnS和tgmS
20、表达式中包含地面反射矩阵,该矩阵通过斯奈尔反射系数乘以粗糙度修正系数求得22 2022 20exp( 2 cos ) 0=0 exp( 2 cos )vihirHrHkRk(18) 其中vr , hr分别为垂直极化和水平极化的斯奈尔反射系数,H为描述地面粗糙度的均方根高度。 各项体散射机制中单个散射体的后向散射包含入射波对散射体的一次作用以及周围散射体对该散射体的作用,不同类型散射体的散射特性计算方法不同,根据文献1 7中方法可计算垂直有限长圆柱体、圆盘的散射矩阵,由于生成的树枝的位置和取向角是随机的,需进行欧拉角和极化基变换,求得任意取向散射体的散射矩阵。 对于面积为A的植被场景,归一化的散
21、射系数为 2114,NMb bg tg gsi n n mnmA kk S S S S (19) 2.3 两相邻散射体间高阶散射 对于电稠密植被,需要考虑散射体间近场互作用,下面应用互易定理求解两散射体间二阶散射场。万方数据第10期 饶丽婷等: 基于分形结构的植被高阶相干散射模型研究 2505 假设在空间中有两相邻电介质目标,首先考虑第1种情况如图4(a)所示,在不考虑目标2存在时入射平面波在目标1上感应了极化体电流1J,体电流1J在空间激发电场11E,电场11E激发了目标2的散射场22E,此时将电场11E和散射场22E之和记为1E。然后考虑第2种情况如图4(b)所示,在远场观测点P点处放置单
22、位点电流源 0eJ prr,单位点电流源极化方向p为散射方向矢量sk的水平极化矢量sh或者垂直极化矢量sv,在只有目标2存在下单位电流源在空间激发电场2eE。 应用Carson互易原理22,可得 11 12deVv pE J E (20) 同理,可得 22 21deVv pE J E (21) 其中,2E为由目标2的体电流2J在空间激发的电场21E和目标1对21E的散射场12E之和,1eE为只存在目标1时单位点电流eJ在空间激发的电场。式(19)中由单位点电流源在空间激发的电场2eE可以分解为点电流源直接激发的电场edE和目标2对edE的散射场2esE,因此可以将式(20)写为 11111 2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 结构 植被 相干 散射 模型 研究 饶丽婷
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内