[发明专利]多频双极化林地积雪被动微波混合像元分解方法

权利要求书

1.一种多频双极化林地积雪被动微波混合像元分解方法，该方法的应用条件是冬季林地 观测地区的被动微波遥感数据和分辨率光谱遥感数据中的下垫面土地分类数据，具体包括如 下过程：

1)林地下垫面土地分类数据的获取和重新分类；

2)建立多频双极化林地积雪被动微波混合像元分解模型；

3)基于动态窗口数据选择策略的欠定性方程组求解；

所述林地下垫面土地分类数据的获取和重新分类的步骤：获取分辨率光谱遥感数据中的 土地覆盖分类产品，将分辨率光谱遥感数据中的下垫面土地分类数据进行重新分类，具体如 下：

a)将原分类数据中定义分类属于针叶植被类型的数据统一重新定义为针叶林类型B；

b)将原分类数据中定义为常绿阔叶植被、落叶阔叶植被和一年生阔叶植被类型的数据统 一重新定义为阔叶林类型N；

c)其他分类数据S₁,S₂,…S_n保持原类型不变；

所述建立多频双极化林地积雪被动微波混合像元分解模型的步骤：所述多频双极化林地 积雪被动微波混合像元分解模型为：

TKHM(x0,y0)TKVM(x0,y0)TKaHM(x0,y0)TKaVM(x0,y0)=αB(x0,y0)TKHBTKVBTKaHBTKaVB+αN(x0,y0)TKHNTKVNTKaHNTKaVN+αS1(x0,y0)TKHS1TKVS1TKaHS1TKaVS1+αS2(x0,y0)TKHS2TKVS2TKaHS2TKaVS2+...+αSn(x0,y0)TKHSnTKVSnTKaHSnTKaVSn+ϵ---(8)]]>

式(8)中T^M代表被动微波辐射计在不同频段和极化方式观测到的的林地积雪亮温值， 下脚标K和Ka分别代表被动微波K和Ka频段，下脚标H和V分别代表水平极化和垂直极化； 两种极化方式；ε代表残差量；

其中，

αC(x0,y0)=∫∫G‾(x,y)LC(x,y)dxdy,C=B,N,S1,S2,...Sn---(7)]]>

式(7)中，L(x,y)代表(x,y)位置的下垫面土地分类数据，上脚标C、B和S₁,S₂,…S_n分 别代表该位置的下垫面类型、阔叶林下垫面土地类型、针叶林下垫面土地类型和其他下垫面 土地类型；若(x,y)位置存在该下垫面，则L(x,y)＝1，否则，L(x,y)＝0；是天线增益函数， 且有

LB(x,y)+LN(x,y)+LS1(x,y)+...+LSn(x,y)=1---(5)]]>

G‾(x,y)=ln2π·exp[-ln2·(x2a2+y2b2)]---(2)]]>

式(2)中，a和b分别表示被动微波遥感数据中K和Ka频段所对应的足印；

所述基于动态窗口数据选择策略的欠定性方程组求解的步骤：使用多点观测约束求 解，当观测点数量为p时，假设在观测点周围内所有的同类下垫面有相同的亮值，则式 (8)转换为

TKHM(1)...TKHM(p)TKVM(1)...TKVM(p)TKaHM(1)...TKaHM(p)TKaVM(1)...TKaVM(p)=TKHBTKHNTKHS1...TKHSnTKVBTKVNTKVS1...TKVSnTKaHBTKaHNTKaHS1...TKaHSnTKaVBTKaVNTKaVS1...TKaVSnα1B...αpBα1N...αpNα1S...αpS......α1Sn...αpSn+E---(10)]]>

式(10)中p是观测点的个数，E代表残差矩阵；

然后，利用非负最小二乘法对其求解，

T＝(A^tA)^-1A^tY(17)

分解误差为：

E＝T-(A^tA)^-1A^tY(18)

式(17)和式(18)中，

Y=TKHM(1)...TKHM(p)TKVM(1)...TKVM(p)TKaHM(1)...TKaHM(p)TKaVM(1)...TKaVM(p)---(12)]]>

T=TKHBTKHNTKHS1...TKHSnTKVBTKVNTKVS1...TKVSnTKaHBTKaHNTKaHS1...TKaHSnTKaVBTKaVNTKaVS1...TKaVSn---(13)]]>

A=α1B...αpBα1N...αpNα1S1...αpS1......α1Sn...αpSn.---(14)]]>

约束条件：

TKVC>TKHC---(15)]]>

TKaVC>TKaHC---(16)]]>

对于公式(10)中p值选取，采用动态窗口数据选择策略：

首先，获取被测地区的多频双极化被动微波遥感数据，根据微波像元的空间相关性，以 待分解的位置(x₀,y₀)处的被动微波混合像元为中心,选取其周围8邻域窗口的被动微波混合 像元和/或周围4邻域窗口的被动微波混合像元作为方程组输入数据进行求解，所对应两种输 入数据的方案为：

(Ⅰ)8邻域窗口数据输入：将含待分解的位置(x₀,y₀)像元在内的9个像元的混合亮温 值代入公式(12)；如果计算得到的结果满足公式(15)和公式(16),获得合理求解，得 到分解后各下垫面对应的组分亮温数据T₁，得到误差数据E₁；

(Ⅱ)4邻域窗口数据输入：将含待分解的位置(x₀,y₀)像元在内的5个像元的混合亮温 值代入公式(12)；如果计算得到的结果满足公式(15)和公式(16),获得合理求解， 得到分解后各下垫面对应的组分亮温数据T₂，得到误差数据E₂。