[发明专利]基于连续统去除法的土壤盐分测定方法

权利要求书

1.一种基于连续统去除法的土壤盐分测定方法，其特征在于，包括以下步骤： 

（1）土样的采集与处理 

根据代表性、典型性和集中性原则选取合适的采样区域采集土样，再将土样带回实验室，清除非土壤成分物质后自然风干，将风干后的样品经研磨过筛，以为后续进行光谱测试和盐分含量的测定做准备； 

（2）土样的含盐量测试 

将上述步骤（1）处理后的土样采用残渣烘干法对土样的含盐量进行测定； 

（3）土样的光谱测试 

将上述步骤（2）中测试后的土样放在一个能控制光照条件的暗室内，采用光谱仪进行土壤测试，其中，测试之前先进行白板校正，每个土样采集多条光谱曲线，通过算术平均后得出该土样的实际反射光谱数据； 

（4）光谱数据处理及吸收特征参数提取 

将上述步骤（3）中测试后的土样首先去除每个土样的光谱曲线噪声较大的边缘波段350～399nm、2400～2500nm；然后采用小波去噪法对土壤光谱进行滤波去噪和数据变换，数据变换包括一阶微分处理和连续统去除处理；接着通过连续统去除处理后，提取出吸收峰总面积、最大吸收深度以及面积归一化最大吸收深度； 

（5）光谱指数的处理 

经上述步骤（4）后，再根据相关性分析结果，选择连续统去除的波段构建连续统去除归一化指数、连续统去除差值指数、连续统去除比值指数；选择连续统去除一阶微分的波段构建连续统去除一阶微分归一化指数、连续统去除一阶微分差值指数、连续统去除一阶微分比值指数； 

（6）模型构建与精度评价 

经上述步骤（5）后，先采用偏最小二乘回归的方法建模，然后采用交叉验证法进行模型验证，接着对模型的精度进行合理评价，以得出模型的预测能力。 

2.如权利要求1所述的基于连续统去除法的土壤盐分测定方法，其特征在于：所述步骤（6）中模型的精度采用以下指标来评价，具体包括建模决定系数R_C²、交叉验证决定系数R_CV²、预测决定系数R_P²、相对分析误差RPD和均方根误差RMSE。 

3.如权利要求2所述的基于连续统去除法的土壤盐分测定方法，其特征在于：所述均方根误差RMSE分为建模均方根误差RMSE_C、交叉验证均方根误差RMSE_CV和预测均方根误差RMSE_P；所述建模决定系数R_C²、预测决定系数R_P²、相对分析误差RPD越大，所述建模均方根误差RMSE_C、交叉验证均方根误差RMSE_CV和预测均方根误差RMSE_P越小，说明模型建模、预测的精度越高，反之，则说明模型建模、预测的精度越低。 

4.如权利要求1所述的基于连续统去除法的土壤盐分测定方法，其特征在于，所述步骤（5）中的光谱指数是指某些特定波段的光谱值的组合； 

所述选择连续统去除的波段是选择连续统去除的673nm、755nm波段； 

所述选择连续统去除一阶微分的波段是选择连续统去除一阶微分658nm、678nm波段； 

所述连续统去除归一化指数定义为NICR,其计算公式为： 

所述连续统去除差值指数定义为DICR，其计算公式为： 

DI_CR=R₇₅₅-R₆₇₃； 

所述连续统去除比值指数定义为RICR，其计算公式为： 

所述连续统去除一阶微分归一化指数定义为NIFDCR，其计算公式为：

所述连续统去除一阶微分差值指数定义为DIFDCR，其计算公式 为： 

DI_FDCR=FR₆₇₈-FR₆₅₈； 

所述连续统去除一阶微分比值指数定义为RIFDCR，其计算公式为：

其中，上述计算公式中R₆₇₃、R₇₅₅表示673nm、755nm波段的连续统去除值；FR₆₅₈、FR₆₇₈表示658nm、678nm波段的连续统去除的一阶微分值。