[发明专利]基于室内标样高光谱特征的土壤重金属Cd含量反演建模及其光谱响应特征波段识别方法

说明书

技术领域

本发明涉及土壤重金属含量高光谱反演领域，特别涉及基于室内标样高光谱特征的土壤重金属Cd含量反演建模及其光谱响应特征波段识别的方法。

背景技术

随着城市化、工业化和农业集约化的发展，农用地土壤重金属污染物因残留时间长、隐蔽性强、迁移性小、毒性大、不易被微生物降解、可以通过食物链富集等特点，逐渐受到国内外研究者的普遍关注。然而，因大范围高密采样成本耗费巨大、常态化快速高效动态监测难以实施，大范围开展重金属污染治理至今仍然存在“重金属污染安全隐患点不清、污染空间格局不明”两大突出问题，已成为摆在土壤重金属污染全面与深度治理面前的一道“理论与方法”鸿沟，无法满足国家在广域地理空间快速高效查清重金属污染空间格局、发现流域重金属污染时空演化规律与形成机制、控制污染范围扩大和转移、合理规划农业生产、减少污染产品损害国民健康的重大需求。

与此同时，新兴的遥感探测技术具有能通过非接触的方式对目标物体进行快速、动态、大范围监测的特征，高光谱更是因其多且连续的光谱波段优势逐渐被广泛应用于土壤重金属污染监测研究中，为定量反演土壤重金属含量提供了新的契机，使得大范围、高时空分辨率、高精度快速探测土壤重金属污染靶区的含量空间分布成为可能。

近年来，国内外围绕土壤重金属含量高光谱反演开展了大量的研究工作。截止目前，大量研究已证实：当土壤重金属含量达到一定数值时，在地面可见光、近红外和热红外波段光谱测试中可发现明显的重金属光谱响应特征；在土壤重金属含量偏低时，仍可通过地面土壤有机质、铁猛氧化物和粘土矿物的光谱特征间接推算重金属含量。然而，土壤成分复杂，重金属含量甚微，有效排除土壤水分、有机质、铁锰氧化物等成分光谱干扰，建立重金属含量反演模型，准确识别微量重金属元素光谱响应特征波段，是亟待解决的关键技术难题。

发明内容

为了有效排除土壤水分、有机质、铁锰氧化物等成分光谱干扰，准确识别微量重金属元素光谱响应特征波段，本发明提供一种基于室内标样高光谱特征的土壤重金属Cd含量反演建模及其光谱响应特征波段识别的方法。

为了实现上述技术目的，本发明的技术方案是，一种基于室内标样高光谱特征的土壤重金属Cd含量反演建模及其光谱响应特征波段识别的方法，包括以下步骤：

步骤1：土壤标样制作及其准确重金属Cd含量与土壤光谱数据获取；

步骤2：土壤重金属Cd含量分析与光谱数据预处理；

步骤3：基于步骤2转换得到的光谱数据的土壤重金属Cd含量反演建模，实现光谱响应特征波段识别。

所述的方法，步骤1中所述的土壤标样的制作及准确其重金属Cd含量的获取，包括以下步骤：

步骤1)：充分考虑土壤重金属Cd污染含量实际变化范围及其光谱响应效果，人为设计不同Cd污染含量(>50等级)(单位：mg/kg)下的标准土壤样本；

步骤2)：采用格网采集方式以1m×1m分辨率采集研究区未经重金属污染的背景土壤样本若干(>50份)，留作标样添加土壤底样；

步骤3)：将步骤2)采集到的背景土壤经实验室烘干、研磨处理，准确称量得到标样添加土壤底样，每份100g；

步骤4)：取步骤3)处理后的底样一份，采用三酸消化——原子吸收分光光度法实验室化学分析测定背景土壤重金属Cd含量(单位：mg/kg)；

步骤5)：计算达到步骤1)中设定的Cd污染含量值需向100g底样中添加的1000mg/kg的Cd标准溶液，计算公式如下：

C＝100(A-S)/ρ

式中，C代表加入100g底样土壤中的Cd标准溶液量(单位：ml)，A代表步骤1)中设计的土壤标样的重金属Cd含量(单位：mg/kg)，S代表步骤4)中测得的背景土壤样本的重金属Cd含量，ρ代表Cd标准溶液的密度(单位：kg/m³)。

特别地，向100g土壤底样中人为添加Cd标准溶液后，需经风干处理后测定土壤重金属Cd含量，标样Cd含量最高不超过0.1mg，为底样的10^-4数量级，可忽略不计，因此土壤总量按100g计。

步骤6)：实验室按步骤5)中计算得到的Cd标准溶液量依次添加到步骤3)获取的土壤底样中，制得土壤标准样本；

步骤7)：采用三酸消化一一一原子吸收分光光度法实验室化学分析测定步骤6)中制得的土壤标样，准确获取土壤标样重金属Cd含量。

所述的方法，步骤2中所述的土壤污染标样光谱与背景样本光谱作比值的光谱预处理方法，包括以下步骤：