发明名称 基于四分量光学物理模型的叶面积指数反演方法
摘要 基于四分量光学物理模型的叶面积指数反演方法,涉及叶面积指数反演方法,它为了解决现有二分量物理模型未能准确描述遥感图像植被覆盖像元的真实组成而导致物理模型模拟结果不准确的问题,它包括具体步骤如下:所述四分量光学物理模型为:<img file="dda00001816733100011.GIF" wi="1599" he="123" /><img file="dda00001816733100012.GIF" wi="923" he="93" />步骤一、由遥感图像成像条件与森林类型特点确定四分量光学物理模型中的输入参数;步骤二、设L值按0.1步长递增,根据所述L值及其对应的R<sub>4</sub>建立查找表,由此建立叶面积指数与植被冠层反射率数据一一对应关系;步骤三、采用均方根误差优化技术游历查找表,使均方根误差最小时所对应的叶面积指数即为最优叶面积指数,并将该最优叶面积指数作为反演结果。适用于准确地描述遥感图像植被覆盖。
申请公布号 CN102706293A 申请公布日期 2012.10.03
申请号 CN201210215468.3 申请日期 2012.06.27
申请人 黑龙江工程学院 发明人 王强;刘丹丹;张玉娟;庞勇;李增元;陈尔学;谭炳香
分类号 G01B11/28(2006.01)I 主分类号 G01B11/28(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 1.基于四分量光学物理模型的叶面积指数反演方法,其特征是,它包括具体步骤如下:所述四分量光学物理模型为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><msub><mi>R</mi><mrow><mo>&infin;</mo><mo>,</mo><mi>t</mi></mrow></msub><mi>K</mi><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mn>1</mn></mrow><mn>0</mn></msubsup><mrow><mo>(</mo><msup><mi>e</mi><mi>KLh</mi></msup><mo>-</mo><msup><mi>e</mi><mi>KLh</mi></msup><msup><mi>e</mi><mi>&alpha;Lh</mi></msup><mo>)</mo></mrow><mi>Ldh</mi><mo>+</mo><msub><mi>R</mi><mrow><mo>&infin;</mo><mo>,</mo><mi>&upsi;</mi></mrow></msub><mi>K</mi><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mn>1</mn></mrow><mn>0</mn></msubsup><msup><mi>e</mi><mi>KLh</mi></msup><msup><mi>e</mi><mi>&alpha;Lh</mi></msup><mi>Ldh</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mo>+</mo><msub><mi>R</mi><mi>s</mi></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;L</mi></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mi>KL</mi></mrow></msup><mo>+</mo><msub><mi>R</mi><mi>sz</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>KL</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>KL</mi></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;L</mi></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,α、K均为消光系数:K=(2/π)L[(β<sub>v</sub>-π/2)cos θ<sub>f</sub>+sinβ<sub>v</sub>tan θ<sub>v</sub>sin θ<sub>f</sub> ],α=(2/π)L[(β<sub>s</sub>-π/2)cos θ<sub>f</sub>+sinβ<sub>s</sub>tan θ<sub>v</sub>sin θ<sub>f</sub>],R<sub>4</sub>为四分量光学物理模型模拟的植被冠层反射率,R<sub>∞,υ</sub>为植被光照面反射率,R<sub>∞,υt</sub>为植被阴影面反射率,R<sub>sz</sub>为土壤阴影面反射率,R<sub>s</sub>为土壤光照面反射率,θ<sub>v</sub>表示视线天顶角,θ<sub>f</sub>表示叶子法线天顶角,θ<sub>s</sub>为太阳天顶角,L为叶面积指数,h为相对光学高度,h的取值范围从-1到0;β<sub>v</sub>=arccos(-cotθ<sub>v</sub>cot θ<sub>f</sub>),β<sub>s</sub>=arccos(-cotθ<sub>s</sub>cot θ<sub>f</sub>);步骤一、根据待处理的遥感图像数据的遥感图像成像条件与森林类型特点确定四分量光学物理模型中的输入参数:R<sub>∞,υ</sub>、R<sub>∞,υt</sub>、R<sub>sz</sub>和R<sub>s</sub>均为实地测量获得,θ<sub>v</sub>、θ<sub>f</sub>、θ<sub>s</sub>、α和K均由遥感图像成像条件确定,R<sub>4</sub>只与L的变化相关;所述遥感图像数据中包括六幅图像信息;<img file="FDA00001816732800013.GIF" wi="731" he="124" />为植被-土壤体系所构成的混合像元中光照植被所占比例,<img file="FDA00001816732800014.GIF" wi="481" he="123" />为植被-土壤体系所构成的混合像元中阴影植被所占比例,e<sup>-αL</sup>e<sup>-KL</sup>为植被-土壤体系所构成的混合像元中光照土壤阴影面所占比例,e<sup>-KL</sup>-e<sup>-KL</sup>e<sup>-αL</sup>为植被-土壤体系所构成的混合像元中土壤光照面所占比例;步骤二、设L值按0.1步长递增,根据每个L值通过四分量光学物理模型模拟对应的得到一个R<sub>4</sub>的值,根据所述L值及其对应的R<sub>4</sub>建立查找表,所述查找表由四分量光学物理模型模拟的六幅不同波段图像的植被冠层反射率,叶面积指数构成,由此建立叶面积指数与植被冠层反射率数据一一对应关系;步骤三、采用均方根误差优化技术游历查找表,使均方根误差最小时所对应的叶面积指数即为最优叶面积指数,并将该最优叶面积指数作为反演结果。
地址 150050 黑龙江省哈尔滨市道外区红旗大街999号