发明名称 基于五种地物分类信息的积雪被动微波混合像元分解方法
摘要 本发明属于遥感图像处理技术领域,提供了一种基于五种地物分类信息的积雪被动微波混合像元分解方法。本发明根据影响积雪辐射特性的典型五种下垫面类型,通过获得观测地区中较高空间分辨率的地物分类数据,利用积雪在不同下垫面微波辐射之间差异特性,选择微波天线增益函数和采样率,建立积雪微波混合像元分解模型,采用具有约束条件的最小二乘法迭代计算求解欠定性方程组,实现积雪被动微波混合像元分解。本发明可以有效地解决积雪被动微波混合像元问题,改进积雪参数反演精度,在气候和水文研究以及积雪灾害评估等领域中具有重要应用价值。
申请公布号 CN102608592B 申请公布日期 2013.06.12
申请号 CN201210096736.4 申请日期 2012.04.05
申请人 吉林大学 发明人 顾玲嘉;任瑞治;张爽;王昊丰;孙健
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 长春吉大专利代理有限责任公司 22201 代理人 王恩远
主权项 1.一种基于五种地物分类信息的积雪被动微波混合像元分解方法,该方法的应用条件是冬季积雪被动微波遥感数据,其特征在于,方法包括如下过程:1)确定观测地区冬季积雪情况下的五种的地物分类数据,2)确定微波天线增益函数和采样率,3)建立积雪被动微波混合像元分解模型,4)积雪被动微波混合像元分解模型求解;所述的确定观测地区冬季积雪情况下的五种的地物分类数据,是将已有的观测地区的地物分类数据产品重新分类为水体、草地、林地、农田和裸土五种下垫面信息,确定观测地区冬季积雪情况下的五种地物分类数据的结果L;所述的确定微波天线增益函数和采样率,是根据被动微波辐射计的特性,选择被动微波数据中不同频率对应的椭圆足印a和b,设计微波天线增益函数<img file="FDA00002962794500011.GIF" wi="77" he="76" />水平采样率Δx和垂直采样率Δy,结合地物分类数据的结果L,利用公式(12)求得系数α:<maths num="0001"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mo>-</mo><mi>a</mi></mrow><mi>a</mi></munderover><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mo>-</mo><mi>bf</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow><mrow><mi>bf</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></munderover><mover><mi>G</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>&Delta;y</mi><mo>)</mo></mrow><mi>&Delta;x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mn>1</mn><mo>-</mo><mfrac><msup><mi>x</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mn>2</mn></msup></mfrac></msqrt><mo>;</mo></mrow></math>]]></maths>α是天线增益函数<img file="FDA00002962794500014.GIF" wi="62" he="77" />和对应足印内地物分类数据的结果L的卷积,(x,y)是被动微波数据的位置;所述的建立积雪被动微波混合像元分解模型,是将地物分类数据和微波天线增益函数及采样率的计算结果代入积雪被动微波混合像元分解模型,确定积雪被动微波混合像元分解模型求解方法;具体的积雪被动微波混合像元分解模型为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>T</mi><mi>B</mi></msub><mo>=</mo><mo>&Integral;</mo><mo>&Integral;</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>{</mo><msub><mi>L</mi><mi>land</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>land</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>water</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>water</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>forest</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>forest</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>+</mo><msub><mi>L</mi><mi>grass</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>grass</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>crop</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>crop</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>}</mo><mi>dxdy</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mtext>16</mtext><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mo>=</mo><msub><mi>&alpha;</mi><mi>land</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>land</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&alpha;</mi><mi>water</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>water</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&alpha;</mi><mi>forest</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>forest</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mo>+</mo><msub><mi>&alpha;</mi><mi>grass</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>grass</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&alpha;</mi><mi>crop</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>crop</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中:<maths num="0007"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mo>&Integral;</mo><mo>&Integral;</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>L</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>dxdy</mi><mo>,</mo><mi>c</mi><mo>=</mo><mi>land</mi><mo>,</mo><mi>water</mi><mo>,</mo><mi>forest</mi><mo>,</mo><mi>grass</mi><mo>,</mo><mi>crop</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(16)中,T<sub>B</sub>代表积雪混合像元亮温值,(x,y)代表混合像元空间位置,即,被动微波数据的位置;T<sub>land</sub>、T<sub>water</sub>、T<sub>forest</sub>、T<sub>grass</sub>和T<sub>crop</sub>分别代表下垫面为裸土、水体、林地、草地和农田的亮温值,L<sub>land</sub>、L<sub>water</sub>、L<sub>forest</sub>、L<sub>grass</sub>和L<sub>crop</sub>分别代表与混合像元空间位置(x,y)匹配的分类数据中裸土、水体、林地、草地和农田像元的比例;公式(17)的结果由公式(12)计算得到;所述的积雪被动微波混合像元分解模型求解,是最终将积雪被动微波混合像元数据分解成五种类型的被动微波组分亮温数据;具体过程是利用建立的积雪被动微波混合像元分解模型,结合被动微波遥感数据,建立联立方程组,对m×n范围内的被动微波混合像元进行分解,通过构建方程组和约束性最小二乘法求解m×n范围内被动微波混合像元分解后的五种地物的组分亮温T<sub>c</sub>。
地址 130012 吉林省长春市前进大街2699号