CN107505283A - 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 - Google Patents
近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 Download PDFInfo
- Publication number
- CN107505283A CN107505283A CN201710785047.7A CN201710785047A CN107505283A CN 107505283 A CN107505283 A CN 107505283A CN 201710785047 A CN201710785047 A CN 201710785047A CN 107505283 A CN107505283 A CN 107505283A
- Authority
- CN
- China
- Prior art keywords
- content
- spectral
- spectrum
- sample
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 48
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 title abstract description 31
- 230000003595 spectral effect Effects 0.000 claims abstract description 69
- 238000001228 spectrum Methods 0.000 claims abstract description 46
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000002790 cross-validation Methods 0.000 claims abstract description 6
- 238000009614 chemical analysis method Methods 0.000 claims abstract description 5
- 238000001514 detection method Methods 0.000 claims description 27
- 238000004458 analytical method Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 15
- 238000002310 reflectometry Methods 0.000 claims description 12
- 238000004611 spectroscopical analysis Methods 0.000 claims description 10
- 238000010238 partial least squares regression Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000013102 re-test Methods 0.000 claims description 5
- 238000010561 standard procedure Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims 2
- 229910002651 NO3 Inorganic materials 0.000 description 27
- 238000004587 chromatography analysis Methods 0.000 description 4
- 230000008439 repair process Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000004204 optical analysis method Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 150000003839 salts Chemical group 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- GRYLNZFGIOXLOG-UHFFFAOYSA-N Nitric acid Chemical compound O[N+]([O-])=O GRYLNZFGIOXLOG-UHFFFAOYSA-N 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical group [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000003929 acidic solution Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005251 capillar electrophoresis Methods 0.000 description 1
- 230000003197 catalytic effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002848 electrochemical method Methods 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910017604 nitric acid Inorganic materials 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种近红外光谱测定次生盐渍土壤中硝酸根离子(NO3 —)含量的方法。利用化学分析方法测定次生盐渍土样品的NO3 —含量;采集土壤样品光谱进行预处理后基于偏最小二乘法(PLSR)建立不同变换全波段光谱与NO3 —含量的关系模型,选择最佳光谱变换形式;在所选择的光谱变换下采用间隔偏最小二乘法(iPLS)对NO3 —的特征波段进行选择,基于优选出的特征波段运用PLSR回归,通过交叉验证获得精度最优的模型即为次生盐渍土壤中NO3 —的关系光谱模型。对未知浓度的次生盐渍土土样测定其近红外光谱数据并进行数据预处理后,代入关系光谱模型,可直接得到未知土样中NO3 —的含量。
Description
技术领域
本发明涉及土壤监测与光谱学分析领域,具体地涉及一种近红外光谱测定次生盐渍土壤中硝酸根离子(NO3 —)含量的方法。
背景技术
土壤的次生盐渍化又称“次生盐碱化”,是由于不合理的耕作灌溉或者人为不合理措施而引起的土壤次生盐渍化问题,其大大制约着我国农业的发展。在土壤次生盐渍化问题中硝酸盐的积聚是土壤次生盐渍化的主要特征之一。因此,实时、高效、快速、无损监测土壤的NO3 —含量以及定量分析对次生盐渍化土壤的形成和修复过程都具有重要的意义。
现有技术中对NO3 —含量的测定主要有传统化学分析方法、离子色谱法、光学分析法、电化学分析法以及其他方法等。在目前采用的诸多检测方法中,主要存在如下不足之处:1、传统的化学分析法,其是依据NO3 —在酸性溶液中的氧化还原性质对待测样品进行预先化学反应处理,然后再进行检测,适用于溶液中NO3 —的检测,并且此方法操作步骤多,分析耗时、耗力,检测结果需要滞后很久才能得到,并且在检测土壤样品时需进行预先的处理,否则存在较大的检测误差,另外,对于有颜色物质干扰时,检测结果也会产生较大误差;2、光学分析法所用仪器复杂,消耗资金较大,需要做标准曲线,并且检测也具有各自的检测限量,超过检测限量也会产生较大的误差;3、化学发光法实验方法步骤一级所用试剂较多,实验方法较为复杂,不易操作及广泛快速推广使用;4、色谱检测法中,如离子色谱法、气相色谱法等常用的色谱法,其精密度高,但所用仪器都为大型仪器,消耗原料和资金都较多。5、其他方法如毛细管电泳法、气相流动分析红外检测法、催化动力学法等也存在方法复杂,不易简便快速操作,结果精密度不高等限制并没有得到广泛的应用。
对于针对次生盐渍化土壤样品中NO3 —含量的检测使用上述现有技术中得方法均不能达到很好的效果,所以对进一步优化NO3 —含量的检测,实现实时、高效、快速、无损监测土壤NO3 —含量以及定量分析存在需求,对次生盐渍化土壤的形成和修复过程都具有重要的意义。
因此,本领域的技术人员致力于开发一种适用于次生盐渍化土壤样品中NO3 —含量的检测方法,解决上述现有技术中检测方法耗时、耗力、滞后、消耗药剂且易对土壤原样品造成破坏等不足之处。
发明内容
鉴于现有技术的上述缺陷,本发明的目的在于针对现有技术的不足和次生盐渍化土壤样品中NO3 —含量检测的需求,解决现有检测方法耗时、耗力、滞后、消耗药剂且易对土壤原样品造成破坏等问题。提供一种利用近红外光谱,主要结合地面高光谱扫描和化学计量学分析对次生盐渍土壤中NO3 —含量进行快速简单精确的检测的方法。
为实现上述目的,本发明提供了一种近红外光谱测定次生盐渍土壤中NO3 —含量的方法,包括首先通过使用高光谱地物波谱仪对多种次生盐渍土壤样品进行近红外光谱数据测试,得到次生盐渍土样品光谱波长与反射率数据,然后建立并优化光谱反射率数据与NO3 —含量的关系光谱模型;检测时将待测土壤样品的近红外光谱反射数据输入关系光谱模型中,得到该待测土壤样品的NO3 —含量数值。
进一步地,所述近红外光谱测定次生盐渍土壤中NO3 —含量的方法中,所述建立光谱反射率数据与NO3 —含量的关系光谱模型包括以下步骤:
步骤1、收集多种不同次生盐渍土壤样品,按照标准方法对土壤样品进行NO3 —含量的检测;
步骤2、将步骤1中次生盐渍土壤样品进行近红外光谱数据测试,得到次生盐渍土样品光谱波长与反射率数据;
步骤3、将步骤1和步骤2中次生盐渍土壤样品光谱数据随机对应抽取三分之二,建立光谱反射率数据与NO3 —含量的关系光谱模型;
步骤4、将步骤3得到的关系光谱模型进行精度优化和评价。
进一步地,所述近红外光谱测试的波长选取范围大于等于400nm,且小于1650nm;
进一步地,所述步骤1中,次生盐渍土壤样品数量为不小于50个,优选50~200个,最优选70~150个;
进一步地,所述步骤1中,标准方法为分析手册中检测NO3 —含量的方法,优选化学分析法;
进一步地,所述步骤2中,次生盐渍土壤样品进行近红外光谱数据测试时,每个样品重复测试后取平均值,得到次生盐渍土样品光谱波长与反射率数据;重复次数为不小于5次;
进一步地,所述步骤3中,建立光谱反射率数据与NO3 —含量的关系光谱模型是利用偏最小二乘法(PLSR)对土壤样品光谱数据的原始光谱数据及光谱变换后的数据分别建立全波段光谱与NO3 —含量的关系光谱模型;
进一步地,所述步骤3中,光谱反射率数据与NO3 —含量的关系模型为全波段光谱范围依次分为10,20和30个间隔区域,利用间隔偏最小二乘回归法(iPLS)分别将不同间隔区域的光谱作为输入变量建立NO3 —含量的PLSR关系光谱模型;
进一步地,所述步骤3中,光谱反射率数据与NO3 —含量的关系光谱模型为各个间隔区域下的最佳敏感波段数据作为自变量,土壤NO3 —含量作为因变量,运用PLSR法分别进行模型的构建,并用交叉验证得出精度最优的PLSR关系光谱模型。
本发明的有益效果在于:
1、本发明的近红外光谱测定次生盐渍土壤中NO3 —含量的方法,在关系光谱模型建立后,对待测的土壤样品进行近红外光谱反射数据检测,将数据输入模型即可快速得到土壤样品的NO3 —含量,与现有检测方法相比,大大减少了等待时间,节约了时间成本;
2、本发明的近红外光谱测定次生盐渍土壤中NO3 —含量的方法,在关系光谱模型建立后,对待测的土壤样品进行检测不需要消耗额外的试剂,节省了检测成本,操作简单,节省了人力成本;
3、本发明的近红外光谱测定次生盐渍土壤中NO3 —含量的方法,在关系光谱模型建立后,对待测的土壤样品进行检测时不需要对样品进行二次处理,避免了对土壤原样品造成破坏等;
4、本发明的近红外光谱测定次生盐渍土壤中NO3 —含量的方法,关系光谱模型的建立依据大量不同次生盐渍土壤样品作为基础支持,采用多种校正,精度评判基准,控制精度参数,确定最优波段区间,得到的模型精度高,检测结果准确可靠。
综上所述,本发明方法近红外光谱测定次生盐渍土壤中NO3 —含量的方法通过建立高精度的最优光谱波长区间与NO3 —含量的关系光谱模型,建立了一种快速、灵敏的土壤NO3 —含量分析检测方法,能快速简单精确检测出待测土壤样品中的NO3 —含量,解决现有检测方法耗时、耗力、滞后、消耗药剂且易对土壤原样品造成破坏等问题;大大提高了分析土壤NO3 —含量时的实效,对次生盐渍化土壤的形成和修复过程都具有重要的意义。用本发明方法监测次生盐渍土修复过程,具有较强的实用性。
附图说明
图1是实施例1的次生盐渍土样品光谱波长与反射率数据图;
图2是实施例1验证集土壤样本的NO3 —含量的关系光谱模型检测值与实际NO3 —含量值的散点比较图。
具体实施方式
本发明所提供的一种近红外光谱测定次生盐渍土壤中NO3 —含量的方法,包括以下步骤:
步骤1、收集50~200个不同次生盐渍土壤样品,按照常规化学分析方法对土壤样品进行NO3 —含量的检测;
步骤2、采用高光谱地物波谱仪对步骤1中的50~200个土壤样进行光谱数据测试,每个样品重复测试5次,取平均值得到次生盐渍土样品光谱波长与反射率数据;
步骤3、将步骤1中次生盐渍土壤样品光谱数据随机抽取三分之二作为建模集土壤样本,剩余为验证集土壤样本;
步骤4、利用偏最小二乘法(PLSR)对建模集土壤样本的原始光谱及光谱变换后的数据分别建立全波段光谱与步骤1中样品对应测得的NO3 —含量的关系光谱模型,并将每一种光谱变换分析对应的模型精度与原始光谱数据对应的PLSR模型精度进行比较,以模型的决定系数(R2)显著提高幅度,标准差(RMSEC)变小程度为评判基准,选择模型精度最佳的光谱变换方式;
步骤5、对步骤4中最佳光谱变换后全波段光谱范围依次分为10,20和30个间隔区域,利用间隔偏最小二乘回归法(iPLS)分别将不同间隔区域的光谱作为输入变量,与步骤1中样品对应测得的NO3 —含量一起建立PLSR关系光谱模型,通过对比各个子区间光谱所建立的模型的最小交互验证均方根误差(RMSECV),分别挑选出各间隔区域下NO3 —含量的的最优波段区间;
步骤6、将步骤5所得各个间隔区域下的最优波段数据作为自变量,将步骤1中对应样品测得的土壤NO3 —含量作为因变量,运用PLSR法分别进行模型的构建,并用交叉验证得出精度最优的关系光谱模型;
步骤7、对步骤3中的验证集土壤样本进行步骤4选择的最佳光谱变换分析,以步骤5中的最佳敏感波段数据作为变量,建立NO3 —含量的PLSR关系光谱模型,得到验证集样本NO3 —含量的预测值,对步骤6中模型进行精度评价,得到精度最优的关系光谱模型;
步骤8、将待测土壤样品进行近红外光谱反射数据检测,将数据输入上述步骤6得到的关系模型中,得到该待测土壤样品的NO3 —含量数值。
在本发明较优的实施方式中,光谱曲线的波长选取范围大于等于400nm,小于1650nm;
在本发明较优的实施方式中,建模集土壤样本采用含量梯度法进行选择,即将样本按NO3 —含量分为多个子范围,根据落入各个范围的样品数随机选取三分之二的样本作为建模集土壤样本,剩余的作为验证集土壤样本;
在本发明较优的实施方式中,步骤5中采用间隔偏最小二乘回归法选择10个间隔区域下的两个最优波段区间为775-899nm和1025-1149nm;
在本发明较优的实施方式中,步骤7中的精度评价还包括以NO3 —实测值为X轴,预测值为Y轴或预测值为X轴,实测值为Y轴,取验证集样品NO3 —的实测值与预测值散点图比较,并添加趋势线。
下面结合实施例对本发明的技术方案作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例1:近红外光谱测定次生盐渍土壤中NO3 —含量的方法,包括以下步骤:
步骤1、收集87个不同次生盐渍土壤样品,剔除异常样品10个,按照常规化学分析方法对77个土壤样品进行NO3 —含量的检测;
步骤2、采用AvaField-2便携式高光谱地物波谱仪对步骤1中的77个土壤样进行光谱数据测试,每个样品重复测试5次,取平均值得到次生盐渍土样品光谱波长与反射率数据,如图1所示;
步骤3、将步骤1中77个次生盐渍土壤样品根据样品中NO3 —含量采用梯度法随机选出51个样本作为建模集土壤样本,剩余26个样本作为验证集土壤样本;
步骤4、利用偏最小二乘法(PLSR)对建模集土壤样本的步骤2中的原始光谱及SGFD光谱变换后的数据分别建立全波段光谱与步骤1中样品对应测得的NO3 —含量的关系光谱模型,并将每一种光谱变换分析对应的模型精度与原始光谱数据对应的PLSR模型精度进行比较,测定模型的决定系数(R2)和标准差(RMSEC)选择模型精度最佳的光谱变换方式;
测定结果为SGFD变换后所建模型决定系数R2为0.933,数值最大,并且其标准差最小,低至0.109,因此选择SGFD变换进行变换分析;
步骤5、对步骤4中建模集土壤样本进行SGFD变换后全波段光谱范围依次分为10,20和30个间隔区域,利用间隔偏最小二乘回归法(iPLS)分别将不同间隔区域的光谱作为输入变量,与步骤1中对应样品测得的NO3 —含量一起建立PLSR关系光谱模型;通过对比各个子区间光谱所建立的模型的最小交互验证均方根误差(RMSECV),分别挑选出各间隔区域下NO3 —含量的的最优波段区间;
结果为:RMSECV值最小的两个最优波段区间为:
NO3 —的10个间隔区域下的最优区间为775-899nm和1025-1149nm,RMSECV值分别为0.281和0.286;
NO3 —的20个间隔区域下的最优区间为841-903nm和967-1029nm,RMSECV值分别为0.285和0.289;
NO3 —的30个间隔区域下的最优区间为820-861nm和988-1029nm,RMSECV值都为0.293;
步骤6、将步骤5所得各个间隔区域下的最优波段数据区间作为自变量,将步骤1中对应样品测得的土壤NO3 —含量作为因变量,运用PLSR法分别进行模型的构建,用交叉验证得出精度最优的关系光谱模型;
在间隔偏最小二乘法中,当全波段被分为10个间隔子区间时,将775-899nm和1025-1149nm作为变量建立的NO3 —含量反演模型精度最高,决定系数R2c达到了0.989,标准差RMSEC值为0.031,说明该方法下模型性能好;
步骤7、对步骤3中的验证集土壤样本进行SFGD变换分析,以775-899nm和1025-1149nm波段数据作为变量,建立NO3 —含量的PLSR关系光谱模型;
经检测,模型得到R2p达到0.962,与校正集样本R2c的0.989数值接近,同时RMSEP值为0.057,说明该方法建立的模型稳定并且预测性能很好;
将验证集土壤样本的NO3 —实测值作为X轴,步骤7所得模型下的NO3 —的预测值作为Y轴,作散点图并添加趋势线,如图2所示;显示趋势线斜率为0.93,说明预测值与实测值接近,说明上述方法建立的模型对次生盐渍土壤NO3 —含量进行检测得到的NO3 —与实际含量值误差较小,具有优良的检测准确性;
在实际运用中,对次生盐渍土壤NO3 —含量进行监控和检测时,只需将待检测的土壤样品进行近红外光谱反射数据检测,将数据输入上述步骤6得到的关系光谱模型中,即可快速得到该待测土壤样品的NO3 —含量数值。
综上所述,本发明技术方案的近红外光谱测定次生盐渍土壤中NO3 —含量的方法中模型精度优,具有简单快速检测、灵敏、精确等优点。
以上详细描述了本发明的具体实施例。应当理解,本领域的普通技术无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此,凡本技术领域中依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验得到的技术方案,皆应在权利要求书所确定的保护范围内。
Claims (10)
1.一种近红外光谱测定次生盐渍土壤中NO3 —含量的方法,其特征在于,首先通过使用高光谱地物波谱仪对多种次生盐渍土壤样品进行近红外光谱数据测试,得到次生盐渍土样品光谱波长与反射率数据,然后建立并优化光谱反射率数据与NO3 —含量的关系光谱模型;检测时将待测土壤样品的近红外光谱反射数据输入关系光谱模型中,得到该待测土壤样品的NO3 —含量数值。
2.根据权利要求1所述方法,其特征在于,所述建立光谱反射率数据与NO3 —含量的关系光谱模型包括以下步骤:
步骤1、收集多种不同次生盐渍土壤样品,按照标准方法对土壤样品进行NO3 —含量的检测;
步骤2、将步骤1中次生盐渍土壤样品进行近红外光谱数据测试,得到次生盐渍土样品光谱波长与反射率数据;
步骤3、将步骤1和步骤2中次生盐渍土壤样品光谱数据随机对应抽取三分之二,建立光谱反射率数据与NO3 —含量的关系光谱模型;
步骤4、将步骤3得到的关系光谱模型进行精度优化和评价;
其中,所述近红外光谱测试的波长选取范围大于等于400nm,且小于1650nm。
3.根据权利要求2所述方法,其特征在于,所述步骤1中,次生盐渍土壤样品数量为不小于50个,优选50~200个;
所述标准方法为分析手册中检测NO3 —含量的方法,优选化学分析法。
4.根据权利要求2所述方法,其特征在于,所述步骤2中,次生盐渍土壤样品进行近红外光谱数据测试时,每个样品重复测试后取平均值,得到次生盐渍土样品光谱波长与反射率数据;重复次数为不小于5次。
5.根据权利要求2所述方法,其特征在于,所述步骤3中,建立光谱反射率数据与NO3 —含量的关系模型是利用偏最小二乘法对土壤样品光谱数据的原始光谱数据及光谱变换后的数据分别建立全波段光谱与NO3 —含量的关系光谱模型。
6.根据权利要求5所述方法,其特征在于,所述光谱反射率数据与NO3 —含量的关系模型为全波段光谱范围依次分为10,20和30个间隔区域,利用间隔偏最小二乘回归法分别将不同间隔区域的光谱作为输入变量建立NO3 —含量的PLSR关系光谱模型。
7.根据权利要求6所述方法,其特征在于,所述光谱反射率数据与NO3 —含量的关系模型为各个间隔区域下的最佳敏感波段数据作为自变量,土壤NO3 —含量作为因变量,运用PLSR法分别进行模型的构建,并用交叉验证得出精度最优的PLSR关系光谱模型。
8.根据权利要求1所述方法,其特征在于,包括以下步骤:
步骤1、收集50~200个不同次生盐渍土壤样品,按照常规化学分析方法对土壤样品进行NO3 —含量的检测;
步骤2、采用高光谱地物波谱仪对步骤1中的50~200个土壤样进行光谱数据测试,每个样品重复测试5次,取平均值得到次生盐渍土样品光谱波长与反射率数据;
步骤3、将步骤1中次生盐渍土壤样品光谱数据随机抽取三分之二作为建模集土壤样本,剩余为验证集土壤样本;
步骤4、利用偏最小二乘法对建模集土壤样本的原始光谱及光谱变换后的数据分别建立全波段光谱与步骤1中样品对应测得的NO3 —含量的关系光谱模型,并将每一种光谱变换分析对应的模型精度与原始光谱数据对应的PLSR模型精度进行比较,以模型的决定系数显著提高幅度,标准差变小程度为评判基准,选择模型精度最佳的光谱变换方式;
步骤5、对步骤4中最佳光谱变换后全波段光谱范围依次分为10,20和30个间隔区域,利用间隔偏最小二乘回归法分别将不同间隔区域的光谱作为输入变量,与步骤1中样品对应测得的NO3 —含量一起建立PLSR关系光谱模型,通过对比各个子区间光谱所建立的模型的最小交互验证均方根误差,分别挑选出各间隔区域下NO3 —含量的的最优波段区间;
步骤6、将步骤5所得各个间隔区域下的最优波段数据作为自变量,将步骤1中对应样品测得的土壤NO3 —含量作为因变量,运用PLSR法分别进行模型的构建,并用交叉验证得出精度最优的关系光谱模型;
步骤7、对步骤3中的验证集土壤样本进行步骤4选择的最佳光谱变换分析,以步骤5中的最佳敏感波段数据作为变量,建立NO3 —含量的PLSR关系光谱模型,得到验证集样本NO3 —含量的预测值,对步骤6中模型进行精度评价,得到精度最优的关系光谱模型;
步骤8、将待测土壤样品进行近红外光谱反射数据检测,将数据输入上述步骤6得到的关系模型中,得到该待测土壤样品的NO3 —含量数值。
9.根据权利要求7所述方法,其特征在于,所述步骤3中,
建模集土壤样本采用含量梯度法进行选择,即将样本按NO3 —含量分为多个子范围,根据落入各个范围的样品数随机选取三分之二的样本作为建模集样本,剩余的作为验证集样本。
10.根据权利要求7所述方法,其特征在于,所述步骤5中,采用间隔偏最小二乘回归法选择10个间隔区域下的两个最优波段区间为775-899nm和1025-1149nm;所述步骤7中,精度评价还包括以NO3 —实测值为X轴,预测值为Y轴或预测值为X轴,实测值为Y轴,取验证集样品NO3 —的实测值与预测值散点图比较,并添加趋势线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710785047.7A CN107505283A (zh) | 2017-09-04 | 2017-09-04 | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710785047.7A CN107505283A (zh) | 2017-09-04 | 2017-09-04 | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107505283A true CN107505283A (zh) | 2017-12-22 |
Family
ID=60695659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710785047.7A Pending CN107505283A (zh) | 2017-09-04 | 2017-09-04 | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107505283A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110132871A (zh) * | 2019-06-06 | 2019-08-16 | 扬州大学 | 一种白菜受土壤次生盐渍化危害的判定方法 |
CN110296956A (zh) * | 2019-07-12 | 2019-10-01 | 上海交通大学 | 一种近红外光谱测定水稻秸秆发酵中有机质含量的方法 |
CN113592743A (zh) * | 2021-08-11 | 2021-11-02 | 北华航天工业学院 | 基于复数小波变换的光谱高频信息与低频信息分离与耦合方法 |
CN113674814A (zh) * | 2021-08-24 | 2021-11-19 | 佛山市海天(高明)调味食品有限公司 | 一种光谱定量分析模型的构建方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076292A (zh) * | 2012-12-27 | 2013-05-01 | 中国地质科学院矿产资源研究所 | 基于地面反射光谱的水体中稀土浓度定量评估系统和方法 |
CN104697937A (zh) * | 2015-03-18 | 2015-06-10 | 北京师范大学 | 一种土壤属性高光谱识别的技术方法 |
CN104897592A (zh) * | 2015-06-11 | 2015-09-09 | 石河子大学 | 基于高光谱技术的盐渍化土壤盐分离子含量监测方法 |
CN106290224A (zh) * | 2016-07-29 | 2017-01-04 | 北京工商大学 | 腊肉品质的检测方法 |
CN107024445A (zh) * | 2017-04-17 | 2017-08-08 | 中国科学院南京土壤研究所 | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 |
-
2017
- 2017-09-04 CN CN201710785047.7A patent/CN107505283A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076292A (zh) * | 2012-12-27 | 2013-05-01 | 中国地质科学院矿产资源研究所 | 基于地面反射光谱的水体中稀土浓度定量评估系统和方法 |
CN104697937A (zh) * | 2015-03-18 | 2015-06-10 | 北京师范大学 | 一种土壤属性高光谱识别的技术方法 |
CN104897592A (zh) * | 2015-06-11 | 2015-09-09 | 石河子大学 | 基于高光谱技术的盐渍化土壤盐分离子含量监测方法 |
CN106290224A (zh) * | 2016-07-29 | 2017-01-04 | 北京工商大学 | 腊肉品质的检测方法 |
CN107024445A (zh) * | 2017-04-17 | 2017-08-08 | 中国科学院南京土壤研究所 | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 |
Non-Patent Citations (1)
Title |
---|
朱赟等: "基于不同光谱变换的土壤盐含量光谱特征分析", 《土壤通报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110132871A (zh) * | 2019-06-06 | 2019-08-16 | 扬州大学 | 一种白菜受土壤次生盐渍化危害的判定方法 |
CN110296956A (zh) * | 2019-07-12 | 2019-10-01 | 上海交通大学 | 一种近红外光谱测定水稻秸秆发酵中有机质含量的方法 |
CN113592743A (zh) * | 2021-08-11 | 2021-11-02 | 北华航天工业学院 | 基于复数小波变换的光谱高频信息与低频信息分离与耦合方法 |
CN113592743B (zh) * | 2021-08-11 | 2024-01-23 | 北华航天工业学院 | 基于复数小波变换的光谱高频信息与低频信息分离与耦合方法 |
CN113674814A (zh) * | 2021-08-24 | 2021-11-19 | 佛山市海天(高明)调味食品有限公司 | 一种光谱定量分析模型的构建方法及装置 |
CN113674814B (zh) * | 2021-08-24 | 2023-06-23 | 佛山市海天(高明)调味食品有限公司 | 一种光谱定量分析模型的构建方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107505283A (zh) | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 | |
CN104062256B (zh) | 一种基于近红外光谱的软测量方法 | |
CN103983595B (zh) | 一种基于紫外-可见光谱处理的水质浊度解算方法 | |
CN101413885A (zh) | 一种快速定量蜂蜜品质的近红外光谱方法 | |
CN102252992A (zh) | 一种对中药提取过程进行近红外在线检测方法 | |
CN101609042A (zh) | 基于近红外光谱的手持式土壤养分无损测量系统 | |
CN103399006B (zh) | 基于颜色rgb分量的尿液分析装置及尿液半定量分析方法 | |
CN110987846B (zh) | 一种基于iPLS-PA算法的硝酸盐浓度预测方法 | |
CN104020129A (zh) | 一种基于近红外光谱结合氨基酸分析技术的工夫红茶发酵质量的判别方法 | |
CN103868880A (zh) | 基于光谱双峰指数的小麦叶片氮含量监测方法及其监测模型的构建方法 | |
CN104316489B (zh) | 一种近红外光谱检测灵芝提取物掺假的方法 | |
CN104764699B (zh) | 一种测定食用油酸价的方法 | |
CN101791331B (zh) | 一种快速测定丹参提取液中鞣质含量的方法 | |
CN105486662A (zh) | 一种基于近红外光谱技术的棉籽棉酚含量的无损测定方法 | |
CN103293118A (zh) | 一种基于近红外光谱分析技术的潲水油鉴别方法 | |
CN103411955A (zh) | 基于拉曼光谱的头孢类抗生素混合溶液的浓度检测方法 | |
CN106560699A (zh) | 用于武夷岩茶产地鉴别的联合检测方法 | |
CN201503392U (zh) | 基于近红外光谱的手持式土壤养分无损测量装置 | |
CN104807777A (zh) | 一种基于近红外光谱分析技术的槟榔水分含量快速检测方法 | |
CN101349638A (zh) | 果蔬维生素c含量的光谱快速无损检测方法 | |
CN103308475A (zh) | 一种同时测量后处理料液中Pu(Ⅳ)及HNO3含量的方法 | |
Lu et al. | Quantitative NIR spectroscopy determination of coco-peat substrate moisture content: Effect of particle size and non-uniformity | |
Wang et al. | Research on vehicle-mounted soil electrical conductivity and moisture content detection system based on current–voltage six-terminal method and spectroscopy | |
CN104865322A (zh) | 一种栀子萃取液浓缩过程快速检测方法 | |
CN103308507A (zh) | 一种快速测定聚氨酯中游离-nco的分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171222 |
|
RJ01 | Rejection of invention patent application after publication |