[发明专利]采用合成孔径雷达技术测量土壤湿度的方法

权利要求书

1.一种采用合成孔径雷达技术测量土壤湿度的方法，其特征在于，包括：

步骤A：利用两个不同时刻t₁和t₂重复观测获取的观测区域的极化合成孔径雷达复图像S₁(x，y)和S₂(x，y)，计算得到观测区域内每一个像素点的归一化干涉复相干系数γ₁₂(x，y)；

步骤B：对于观测区域内的每一个像素点，计算其归一化干涉复相干系数γ₁₂(x，y)的幅度A₁₂(x，y)和相位值φ₁₂(x，y)；

步骤C，在观测区域中，排除归一化干涉复相干系数γ₁₂(x，y)的幅度A₁₂(x，y)小于a的像素组成的区域，获取有效观测区域，以及

步骤D：对有效观测区域内的每一个元素，利用地表相干散射模型，由前一个时刻的土壤湿度值Mv₁计算后一个时刻下的土壤湿度值Mv₂，从而得到有效观测区域内的土壤湿度分布。

2.根据权利要求1所述的方法，其特征在于，所述步骤D中，对有效观测区域内的一个元素，由前一个时刻的土壤湿度值Mv₁计算后一个时刻下的土壤湿度值Mv₂包括：

子步骤D1：由前一个时刻的土壤湿度值Mv₁计算其复介电常数ε_r1；

子步骤D2：由前一时刻土壤的复介电常数ε_r1计算前一时刻对应的极化反射系数α₁；

子步骤D3：若φ₁₂(x，y)＜0，则设定后一个时刻的土壤湿度初始值为Mv₂＝Mv₁-b％，若φ₁₂(x，y)＞0，则设定后一个时刻的土壤湿度初始值为Mv₂＝Mv₁+b％；

子步骤D4：由后一时刻的土壤湿度值Mv₂计算其复介电常数ε_r2；

子步骤D5：由后一时刻土壤的复介电常数ε_r2计算后一时刻对应的极化反射系数α₂；

子步骤D6：由前一时刻对应的极化反射系数α₁和后一时刻对应的极化反射系数α₂，计算归一化干涉复相干系数的估计

子步骤D7，计算归一化干涉复相干系数估计的相位值进而计算相位值与φ₁₂(x，y)的相位差绝对值Δφ；以及

子步骤D8，令Mv₁＝Mv₂，重复执行步骤D3～D7，查找使相位差绝对值Δφ最小的Mv₂取值，即为该像素点的后一个时刻的土壤湿度值。