[发明专利]一种基于DAC的土壤重金属含量估测方法

权利要求书

1.一种基于DAC的土壤重金属含量估测方法，其特征在于，包括以下步骤：

步骤1)数据获取：获取待测区域地面调查得到的土壤重金属采样点数据、MODIS卫星遥感影像数据、ASTERGDEM卫星遥感影像数据以及待测区域地表温度数据；所述的MODIS卫星遥感影像为MOD13A2数据产品；

步骤2)数据预处理：将步骤1)将ASTER GDEM卫星遥感影像数据进行聚合计算得到DEM数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算得到植被指数数据；将反演而来的地表温度数据聚合计算到空间分辨率为1km；将土壤重金属数据整理入库后提取出Hg含量数据，并进行异常值的剔除处理；

步骤3)根据陆表环境变量与Hg元素含量之间的相关关系，进行分区建模：将步骤2)处理后的Hg含量数据作为因变量，以植被指数数据、DEM数据、地表温度作为自变量，将该数据集进行分类建模预测各个采样点的Hg含量。

2.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的分区建模过程中所采用的多元回归模型形式为：

其中，N表示参数估算模型中自变量个数；a_n表示第n个环境变量的系数；a₀表示模型参数的常数项系数；y_n表示Hg含量预测值；x_n表示第n个环境变量；

a₀和a_n的计算公式如下：

其中:k代表土壤重金属采样点的个数；x_in代表第i个土壤重金属采样点的第n个环境变量的值，y_i代表的是第i个土壤重金属采样点的Hg含量，代表所有土壤重金属采样点中第n个环境变量因子的均值，代表所有土壤地面调查样点中Hg含量的均值。

3.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的DEM数据的空间分辨率为1km。

4.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的植被指数数据的空间分辨率为1km。

5.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的分区建模过程中所采用的多元回归模型具体如下：

(1)当dem≤25.02

Y_Hg＝0.07-0.0036X_lstd-2.2e-005×X_dem-0.03×X_ndvi+0.001X_lon

(2)当lat>30.13且dem≤25.02

Y_Hg＝X_ndvi X_dem＝-3.66-0.002553×X_dem+0.13×X_lat

(3)当lat≤30.13且dem≤25.02

Y_Hg＝-3.43-0.004382×X_dem+0.089×X_lat+0.038X_lstd+0.013X_lstn+0.08X_ndvi

其中Y_Hg是土壤重金属Hg含量，X_dem代表的是空间分辨率为1km的数字高程模型的栅格值，X_ndvi代表的是空间分辨率为1km的植被指数栅格值，X_lat代表的是采样点的纬度，X_lon代表的是采样点的经度，X_lstd代表的是白天的地面温度，X_lstn代表的是晚上的地面温度。