发明名称 基于ASAR和Hyperion数据反演植被覆盖下土壤水分的方法及系统
摘要 本发明提供基于ASAR和Hyperion数据反演植被覆盖下土壤水分的方法及系统,具体包括:获取实验区的Hyperion高光谱图像数据,提取实验区中的植被区域;构建基于比值指数及经验系数的植被含水量模型;结合地面实测的植被含水量及Hyperion数据计算的植被比值指数,进行回归计算,获得所述经验系数值;匹配高光谱数据所需的ASAR雷达数据;将匹配后的高光谱数据、ASAR雷达数据和所述植被含水量模型、经验系数值输入水云模型,反演计算获得植被覆盖下土壤含水量。本发明的方法及系统既有明确的物理基础,又具有输入参数少、简单、灵活、易于操作等优点,使植被覆盖下土壤含水量测量更为准确。
申请公布号 CN103196862B 申请公布日期 2015.01.21
申请号 CN201310058961.3 申请日期 2013.02.25
申请人 北京师范大学 发明人 姚云军;程洁;刘素红;张晓通;张楠楠
分类号 G01N21/25(2006.01)I;G01N22/04(2006.01)I 主分类号 G01N21/25(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 王莹
主权项 一种基于ASAR和Hyperion数据反演植被覆盖下土壤水分的方法,其特征在于,该方法包括:S1、获取实验区的Hyperion高光谱图像数据,提取实验区中的植被区域;S2、构建基于比值指数及经验系数的植被含水量模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>SR</mi><mo>=</mo><mfrac><msub><mi>R</mi><mi>SWIR</mi></msub><msub><mi>R</mi><mi>NIR</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000570518740000011.GIF" wi="324" he="181" /></maths>m<sub>v</sub>=a<sub>0</sub>SR+b<sub>0</sub>;其中,SR为植被水分比值指数,m<sub>v</sub>为植被含水量,R<sub>SWIR</sub>和R<sub>NIR</sub>分别为短波红外波段反射率和近红外波段反射率,a<sub>0</sub>和b<sub>0</sub>分别为经验系数;S3、结合地面实测的m<sub>v</sub>及Hyperion数据计算的SR值进行回归计算,获得所述经验系数a<sub>0</sub>和b<sub>0</sub>的值;S4、匹配所述Hyperion高光谱数据所需的ASAR雷达数据;S5、将所述匹配后的Hyperion高光谱数据、ASAR雷达数据和所述植被含水量模型、所述a<sub>0</sub>和b<sub>0</sub>的值输入水云模型,所述水云模型为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><mn>0</mn></msup><mo>=</mo><msubsup><mi>&sigma;</mi><mi>veg</mi><mn>0</mn></msubsup><mo>+</mo><msup><mi>&tau;</mi><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>soil</mi><mn>0</mn></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000570518740000016.GIF" wi="454" he="121" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&sigma;</mi><mi>veg</mi><mn>0</mn></msubsup><mo>=</mo><mi>A</mi><mo>&CenterDot;</mo><msub><mi>m</mi><mi>v</mi></msub><mo>&CenterDot;</mo><mi>cos</mi><mi>&theta;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&tau;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000570518740000012.GIF" wi="653" he="116" /></maths>τ<sup>2</sup>=exp(‑2Bm<sub>v</sub>secθ);其中,σ<sup>0</sup>是植被覆盖地表下总的后向散射系数,<img file="FDA0000570518740000013.GIF" wi="106" he="96" />为植被层的后向散射系数,<img file="FDA0000570518740000014.GIF" wi="93" he="76" />是直接地表后向散射系数,τ<sup>2</sup>是农作物的双程衰减因子,m<sub>v</sub>是植被含水量,A和B分别为依赖于植被类型的参数,θ是雷达入射角;根据<img file="FDA0000570518740000015.GIF" wi="96" he="78" />与植被覆盖下土壤含水量m<sub>s</sub>的关系:σ<sup>0</sup><sub>soil</sub>=a<sub>0</sub>ln(m<sub>s</sub>)+b<sub>0</sub>反演计算获得植被覆盖下土壤含水量m<sub>s</sub>。
地址 100875 北京市海淀区新街口外大街19号