CN111595806A - Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum - Google Patents
Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum Download PDFInfo
- Publication number
- CN111595806A CN111595806A CN202010449425.6A CN202010449425A CN111595806A CN 111595806 A CN111595806 A CN 111595806A CN 202010449425 A CN202010449425 A CN 202010449425A CN 111595806 A CN111595806 A CN 111595806A
- Authority
- CN
- China
- Prior art keywords
- soil
- carbon
- model
- spectrum
- mid
- 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 120
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 93
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000001228 spectrum Methods 0.000 title claims abstract description 45
- 238000012544 monitoring process Methods 0.000 title claims abstract description 12
- 239000000126 substance Substances 0.000 claims abstract description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 238000001055 reflectance spectroscopy Methods 0.000 claims abstract description 5
- 238000005516 engineering process Methods 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims description 19
- 230000003595 spectral effect Effects 0.000 claims description 16
- 238000012795 verification Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 238000000227 grinding Methods 0.000 claims description 7
- 239000002245 particle Substances 0.000 claims description 7
- 238000004445 quantitative analysis Methods 0.000 claims description 7
- 238000013210 evaluation model Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000002203 pretreatment Methods 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 239000004677 Nylon Substances 0.000 claims description 4
- 229920001778 nylon Polymers 0.000 claims description 4
- 238000007873 sieving Methods 0.000 claims description 4
- 238000007605 air drying Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000002411 adverse Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims 4
- 238000000556 factor analysis Methods 0.000 claims 2
- 239000004575 stone Substances 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 8
- 238000002329 infrared spectrum Methods 0.000 abstract description 6
- 230000036541 health Effects 0.000 abstract description 2
- 238000013507 mapping Methods 0.000 abstract description 2
- 238000010238 partial least squares regression Methods 0.000 abstract 2
- 238000011160 research Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000009614 chemical analysis method Methods 0.000 description 3
- 238000004476 mid-IR spectroscopy Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 229910000661 Mercury cadmium telluride Inorganic materials 0.000 description 2
- MCMSPRNYOJJPIZ-UHFFFAOYSA-N cadmium;mercury;tellurium Chemical compound [Cd]=[Te]=[Hg] MCMSPRNYOJJPIZ-UHFFFAOYSA-N 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- KMUONIBRACKNSN-UHFFFAOYSA-N potassium dichromate Chemical compound [K+].[K+].[O-][Cr](=O)(=O)O[Cr]([O-])(=O)=O KMUONIBRACKNSN-UHFFFAOYSA-N 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000000985 reflectance spectrum Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- 241000219146 Gossypium Species 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000035558 fertility Effects 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000003041 laboratory chemical Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000009336 multiple cropping Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 235000013311 vegetables Nutrition 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- 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/47—Scattering, i.e. diffuse reflection
- G01N21/4738—Diffuse reflection, e.g. also for testing fluids, fibrous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
- G01N2021/3572—Preparation of samples, e.g. salt matrices
-
- 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
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/129—Using chemometrical methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a method for monitoring soil carbon components by utilizing a mid-infrared diffuse reflection spectrum technology, belonging to the technical field of soil detection. In particular to a method for rapidly detecting the content change of total carbon, organic carbon and inorganic carbon in farmland soil by using a mid-infrared diffuse reflection spectrometer. The method adopts a diffuse reflectance spectroscopy method, takes the measured value of the soil carbon by a chemical method as a basis, and combines a Partial Least Squares Regression (PLSR) algorithm to establish a mid-infrared spectrum (MIR) model of three soil carbon components of total soil carbon, organic carbon and inorganic carbon in Quzhou county, thereby establishing a method for quickly, accurately and efficiently detecting the content of the total soil carbon, the organic carbon and the inorganic carbon. The method can rapidly and accurately realize the determination of the carbon components of the three types of soil, has analysis efficiency greatly superior to that of the traditional method, provides a feasible method for carrying out real-time monitoring on large-area soil carbon content change in China, and also provides a new technical means for the work of precision agriculture, soil digital mapping, agricultural ecological environment protection, soil health evaluation and the like.
Description
Technical Field
The invention belongs to the technical field of soil detection, and particularly relates to a method for monitoring a soil carbon component by using a mid-infrared diffuse reflection spectrum. In particular to a method for rapidly detecting the content change of total carbon, organic carbon and inorganic carbon in soil of county farmland by utilizing a Fourier transform mid-infrared diffuse reflection spectrometer.
Background
Soil carbon, one of the key attributes of soil, plays an important role in soil structure formation, soil fertility preservation, global climate change regulation and land ecosystem stability maintenance, is the core of soil quality, and is the important basis of sustainable agriculture. Soil carbon, especially organic carbon, directly participates in various physical, chemical and biological processes of soil, and analysis of different soil carbon component contents has important significance for determining soil carbon sources, grasping soil quality change in time and guiding agricultural production, and is always important research content in the subject fields of agriculture, environment and the like. At present, the conventional chemical analysis method for soil carbon, such as potassium dichromate oxidation method, for determining the content of total carbon and organic carbon in soil has the defects of low determination efficiency, poor real-time performance, pollution and the like, so the conventional chemical analysis method cannot meet the requirement of rapidly and efficiently monitoring the content change of soil carbon, and the development of large-area soil carbon research and monitoring work is seriously hindered.
The mid-infrared spectroscopy based on the theory that the fundamental frequency is absorbed by the vibration in the molecule has the characteristics of strong spectral characteristics, easier information extraction, good reproducibility, simple and convenient operation, high testing speed, less sample consumption, environmental protection and suitability for batch sample determination, and the application of the mid-infrared spectroscopy in soil carbon analysis is increasing day by day, thereby gradually becoming the key point of research of scholars at home and abroad. The invention adopts a mid-infrared diffuse reflection spectrum method, and the basic principle is as follows: when light irradiates the sample and enters the sample through refraction, the light reacts with sample molecules to generate reflection, refraction, diffraction and absorption phenomena, finally, the light scattered to all directions of the space from the surface of the sample is called diffuse reflection light, and the diffuse reflection light enters the sample once and reacts with the sample molecules, so the diffuse reflection light carries the structure and composition information of the sample, and the basis of the technical work of the diffuse reflection spectrometer is provided.
At present, domestic research for quantifying soil carbon based on a mid-infrared diffuse reflection spectroscopy method is less. Therefore, the soil in Quzhou county is taken as a research object, the diffuse reflection spectroscopy is adopted, the soil carbon measurement value in the traditional method is taken as a basis, the PLSR algorithm is combined, and the MIR prediction model of total soil carbon, organic carbon and inorganic carbon in the Quzhou county is established, so that the rapid measurement of various soil carbon component indexes is realized, and the method has important significance and practical value for further developing the research of the mid-infrared spectroscopy technology on the aspect of soil carbon measurement in China, developing green agriculture, digital agriculture and the like.
Disclosure of Invention
The invention aims to provide a method for monitoring soil carbon components by using mid-infrared diffuse reflectance spectrum, which is characterized in that a mid-diffuse reflectance spectrum method is adopted, a chemical method soil carbon measurement value is taken as a basis, a PLSR algorithm is combined, and MIR prediction models of soil total carbon, organic carbon and inorganic carbon are established, so that a method for rapidly, accurately and efficiently detecting the content of the soil total carbon, the organic carbon and the inorganic carbon is established, and the method comprises the following steps:
(1) collecting and preparing farmland soil samples: the method comprises the following steps of collecting county soil samples by adopting a 1km x 1km grid survey method, coding according to a sampling sequence and recording sample information in detail, wherein the method specifically comprises the following steps: sample number, collection date, sampling depth, longitude and latitude, altitude, land utilization mode and sampler; air-drying all collected soil samples, removing sundries such as plant residues, gravels and the like, grinding, passing all the soil samples through a nylon sieve with the aperture of 2mm, and storing for later use;
(2) collecting diffuse reflection spectrograms of infrared bands in all soil samples to obtain an original spectrogram: a Tensor II Fourier transform mid-infrared spectrometer developed and produced by German Bruker company is adopted to measure the spectrum of a soil sample, and the instrument is provided with an MCT (cadmium mercury telluride) detector cooled by liquid nitrogen and used for collecting 7500-600 cm-1Resolution of 4cm-1Scanning 32 times; measuring two spectral data of each soil sample, and taking the data after arithmetic mean as the original spectral data of the soil sample; the data storage uses the OPUS format. However, for spectral analysis and modeling, only 4000-600 cm of the mid-infrared band spectrum is used-1;
(3) Representative sample screening: uniformly selecting samples in the spectral feature space by using a KS method to screen out soil samples with strong representativeness for determination of a soil carbon chemical method and MIR modeling;
(4) selecting a spectrum pretreatment method: in order to reduce the influence of the state of the spectrometer, the measurement environment and the sample grinding and sieving, 3 types of spectrum preprocessing of baseline correction, scattering correction and scale scaling are carried out on all spectrum data, and each type of spectrum preprocessing comprises a plurality of preprocessing methods. The baseline correction comprises a first derivative and a second derivative, and mainly deducts the influence of the background or drift of the instrument on the spectral signal; the scattering correction comprises Multivariate Scattering Correction (MSC) and standard normal variable transformation (SNV), and mainly eliminates the influence of scattering on the spectrum caused by uneven particle distribution and different particle sizes; the scale scaling comprises maximum and minimum normalization, constant offset elimination and straight line subtraction, and mainly eliminates adverse effects caused by overlarge scale difference. Comparing the model established after the spectrum pretreatment with the model established without the pretreatment, and selecting the optimal spectrum pretreatment method of each soil carbon component from the models;
(5) constructing a spectrum prediction model: dividing the representative samples screened in the step (3) into a correction set and a verification set according to a KS method, establishing a regression relation between the spectral signals and the component content by adopting a PLSR algorithm and taking the infrared spectrum data independent variable of the soil samples of the correction set and the soil carbon content as dependent variables, and verifying the accuracy of the model by using the independent verification set. In modeling PLSR, the factor number is selected as the most important factor. Setting the maximum factor number as 20, and determining by a leave-one-cross verification method, wherein the optimal factor number is the factor number corresponding to the minimum RMSECV value; the specific process of leave-one-out cross validation includes: selecting a sample test model from the correction set each time to predict the quality, using the rest samples to train data, repeating the process until each sample is set as test data, and calculating RMSECV;
(6) and (3) evaluating a model: r with verification setp 2And performing model evaluation on the three parameters of RMSEP and RPD. Wherein R is2Evaluation of the fitness of the model, R2The higher the value, the better the fitting effect of the model, and the stronger the capability of the model to explain the dependent variable; the stability of the RMSE evaluation model and the predictability of the RPD evaluation model are realized, and the higher the RPD value is, the lower the RMSE is, and the better the model is; in addition, further makeAnd (3) carrying out quality class classification on the model by using an RPD value: when the RPD is more than or equal to 3.0, the calibration effect is good, and the established calibration model can be used for actual detection; when 2.5<RPD<3.0, quantitative analysis is feasible, but the prediction precision is to be improved; when RPD is less than or equal to 2.5, it is difficult to use for quantitative analysis.
The invention has the beneficial effects that:
1. the method can realize rapid and accurate prediction of the content of the total carbon, organic carbon and inorganic carbon in the soil.
2. Makes up the defects of the conventional laboratory chemical analysis method.
3. The soil samples with strong representativeness are selected from a large number of samples to be used for the determination of soil carbon by a chemical method and MIR modeling, so that the model building speed is greatly increased, the storage space of a model base is reduced, and the model is convenient to update and maintain.
4. The sampling range covers the whole Quzhou county, and the county-area-range soil carbon component content change real-time monitoring can be realized.
5. The method is expected to be popularized to soil carbon analysis in northern China and even national farmlands and other ecological systems such as grasslands, forests and the like.
Drawings
FIG. 1 is a schematic diagram showing the correlation between the predicted value and the actual value of the soil total carbon spectrum model.
FIG. 2 is a schematic diagram showing the correlation between the predicted value and the actual value of the soil organic carbon spectrum model.
FIG. 3 is a schematic diagram showing the correlation between the predicted value and the true value of the soil inorganic carbon spectrum model.
Detailed Description
The invention provides a method for monitoring soil carbon components by utilizing a mid-infrared diffuse reflection spectrum technology, and the invention is further explained by combining with an embodiment.
The method selects a typical agricultural county Quzhou county in Huang-Huai-Hai plain of the food main producing area in China as a research area. Study area overview: in the middle of North China plain in Quzhou county, 5 towns and 342 administrative villages in the district, the total area of the land is 667km2Cultivated area 47782hm2(ii) a The parent soil is the flood and alluvial deposit of river, and the soil is selected from the group consisting of moisture soil,Saline soil and brown soil are the main, and the soil texture comprises sandy soil, sandy loam, light loam, medium loam and clay; the regional crops mainly comprise wheat, corn, cotton and vegetables, and have high multiple cropping index, large land investment and higher agricultural intensification degree.
(1) Soil sample collection and preparation
The sampling method comprises collecting farmland soil samples by 1km × 1km grid survey method, wherein the total area of the whole county of Quzhou is 667km2The method comprises the steps of setting 730 grids, not sampling county-area edge grids, leaving 580 grids after removal, collecting soil samples in a selected farmland by an S-shaped sampling method, collecting 5 soil samples in each farmland soil, and preparing 1 part of soil sample by a quartering method. Coding is carried out according to the sampling sequence and the sample information is recorded in detail, which specifically comprises the following steps: sample information such as sample number, acquisition time, sampling depth, longitude and latitude, altitude, land utilization mode, sampler and the like.
Preparing a soil sample: spreading the collected soil sample indoors for natural air drying, removing sundries such as plant residues, gravels and the like, grinding, passing through a nylon sieve with the aperture of 2mm, and storing for later use.
(2) Mid infrared spectral collection
Because the farmland soil sample is solid particles with irregular shapes, and the size, the shape and the uniformity of the particles can greatly influence the diffuse reflection spectrum data, the influence of the particle scattering effect on the spectrum acquisition can be reduced as much as possible by grinding and sieving. In addition, the soil water has strong absorption on the mid-infrared spectrum and is overlapped with the absorption waveband of the organic carbon, so that the prediction accuracy of the mid-infrared spectrum model of the soil carbon is obviously influenced. Therefore, further processing of the soil sample is required prior to collection of the mid-infrared spectrum. Specifically, the method comprises grinding, sieving with 0.25mm nylon sieve, drying in an oven at 105 deg.C for 2 hr, and storing in a desiccator until MIR collection is completed.
The spectrometer used was a Tensor II Fourier transform mid-infrared spectrometer developed by Bruker, Germany. The following preparations are required before spectrum acquisition: reducing the temperature and humidity of the laboratory environment, wherein the temperature is less than or equal to 25 ℃, and the relative humidity is less than or equal to 60%; the preheating of the instrument is more than or equal to 1h after starting up. Then the instrument is operatedThe stability of the measuring software OPUS 7.5, and the position of the interference peak is determined by using a check signal in an 'advanced measurement' menu; the basic state of the instrument is detected by using the running OQ test under the menu of the instrument state, and whether the instrument is normal is checked by using the running PQ test. Setting scanning parameters: spectral range of 7500-600 cm-1Resolution of 4cm-1The number of scans was 32. Soil samples were placed in instrument-fitted 24-well sample trays to ensure that the sample surface in each well was flat, with the wells 5mm in diameter and 3mm deep. The diffuse reflection standard sample is placed at the upper left corner of the sample disc, the background is subtracted from each scanning to correct the noise in the atmosphere and the instrument, 2 pieces of spectral data are collected from each soil sample, the data after arithmetic averaging is used as the original spectral data of the soil sample, and the data storage format is OPUS.
(3) Representative sample screening and soil carbon reference value determination
Representative sample selection within the study area: and screening 80 soil samples by using a K-S method for measuring reference values of a soil carbon chemical method and establishing an MIR model. The MIR spectrum is high in collection speed and small in sample amount, so that a representative sample in a research area is screened out based on the sample spectrum, the test time and the analysis cost of a large-batch sample chemical method can be reduced, more importantly, when a sample outside a model is encountered, the application range of the model can be expanded through fewer samples, and the model can be updated and maintained conveniently.
The chemical determination method of the soil carbon comprises the following steps: the contents of total carbon and organic carbon in the soil are measured by adopting an Elementar vacuum MACROcube element analyzer in Germany, and the content of inorganic carbon is obtained by adopting the subtraction method of the contents of total carbon and organic carbon.
(4) Spectrum prediction model construction
For spectral analysis and modeling, only 4000-600 cm of mid-infrared band spectrum is used-1. The use of different spectral preprocessing methods has different improvements and influences on the soil carbon prediction model by eliminating 10 kinds of light such as constant offset, subtracting a straight line, SNV, maximum and minimum normalization, MSC, first derivative, second derivative, first derivative + subtracting a straight line, first derivative + SNV, first derivative + MSCThe spectrum preprocessing method performs spectrum preprocessing on the raw spectrum data of the 80 soil samples. And comparing the model established after the spectrum pretreatment with the model established without the pretreatment, and selecting the optimal spectrum pretreatment method of each soil carbon component from the models. Sample set division is needed before modeling, and the area representative samples screened in the step (3) are divided into a correction set and a verification set, wherein the correction set is 70% and 56 samples, and the verification set is 30% and 24 samples. For the selection of calibration and validation set samples, the division was performed according to the range of reference values determined by chemical methods in combination with the K-S method. The modeling method adopts a classical method PLSR in spectral analysis, takes mid-infrared spectrum data independent variable of a soil sample of a correction set and soil carbon content as a dependent variable, establishes a regression relationship between a spectral signal and component content, and utilizes an independent verification set to test the accuracy of a model. In the modeling of PLSR, the most important factor (LV) is the selection. Selecting PLSR factor number: the maximum factor number is set to be 20, and the optimal factor number is the factor number corresponding to the minimum RMSECV value determined by a leave-one-out cross-validation method. The optimal factor numbers for each carbon component model are detailed in table 1.
(5) Model evaluation and application
R with verification setp 2And performing model evaluation on the three parameters of RMSEP and RPD. Wherein R is2Evaluation of the fitness of the model, R2The higher the value, the better the fitting effect of the model, and the stronger the capability of the model to explain the dependent variable; stability of the RMSE evaluation model and predictability of the RPD evaluation model, and the higher the RPD value, the lower the RMSE and the better the model. In addition, the RPD values are further used for quality class classification of the model: when the RPD is more than or equal to 3.0, the calibration effect is good, and the established calibration model can be used for actual detection; when 2.5<RPD<3.0, quantitative analysis is feasible, but the prediction precision is to be improved; when RPD is less than or equal to 2.5, it is difficult to use for quantitative analysis. The results of the model evaluations are detailed in Table 1.
After the model is established, the mid-infrared diffuse reflection spectrum information of soil samples with unknown carbon content and similar properties in Quzhou county can be brought into the model, and the content of the carbon component in the soil can be predicted.
TABLE 1 accuracy evaluation results of PLSR model for soil carbon component content in Quzhou county
The final quantitative model is established by combining the optimal PLSR factor number and a spectrum preprocessing method, and the result shows that: the RPD values of the model of the content of carbon components in the soil around the curve are all higher than 3.0, R2The values are all larger than 0.90, which indicates that the quantitative analysis result of the established model is more accurate. For the spectroscopic model, the stability of the model can be further analyzed by comparing the predicted and actual values of the model. A schematic diagram of the correlation between the chemical measured values of the total carbon, organic carbon and inorganic carbon of the soil and the PLSR predicted value is drawn (as shown in fig. 1, fig. 2 and fig. 3), and it can be seen from the diagram that the predicted values and measured values of the three soil carbon components have good consistency in the whole data set range, and the data points are uniformly distributed in 1: both sides of line 1 indicate higher stability of the model.
In conclusion, the model can be used for rapidly and accurately measuring the carbon components of the three types of soil, the analysis efficiency is greatly superior to that of the traditional method, a feasible method is provided for carrying out real-time monitoring on large-area soil carbon content change in China, and a new technical means is provided for the work of precision agriculture, soil digital mapping, agricultural ecological environment protection, soil health evaluation and the like.
Claims (5)
1. A method for monitoring soil carbon components by utilizing a mid-infrared diffuse reflection spectroscopy technology is characterized in that a diffuse reflection spectroscopy method is adopted, soil carbon measurement values in a traditional method are taken as the basis, a PLSR algorithm is combined, and an MIR prediction model of soil carbon is established, so that a method for rapidly, accurately and efficiently detecting the content of three soil carbon components including total soil carbon, organic carbon and inorganic carbon is established; the method comprises the following steps:
(1) collecting and preparing a soil sample: collecting soil samples in county by adopting a 1km multiplied by 1km grid survey method, coding according to a sampling sequence and recording sample information in detail; air-drying all collected soil samples, removing sundries such as plant residues and gravels, grinding, passing through a nylon sieve with the aperture of 2mm, and storing for later use;
(2) collecting diffuse reflection spectrograms of infrared bands in all soil samples to obtain an original spectrogram: the wavelength range is set to be 7500-600 cm-1Resolution of 4cm-1Scanning 32 times; however, for spectral analysis and modeling, only 4000-600 cm of the mid-infrared band spectrum is used-1;
(3) Representative sample selection: uniformly selecting samples in a spectral feature space by using a KS (Kennard-Stone) method based on Euclidean distances among variables to screen out soil samples with strong representativeness for determination of a soil carbon chemical method and MIR modeling;
(4) selecting a spectrum pretreatment method: in order to reduce the influence of the state of a spectrometer, the measurement environment and the sample grinding and sieving, 3 types of spectrum preprocessing of baseline correction, scattering correction and scale scaling are carried out on all spectrum data, and the optimal spectrum preprocessing method of each soil carbon component is selected according to the modeling effect;
(5) constructing a spectrum prediction model: the modeling method adopts a classical method PLSR in spectrum analysis, all samples are divided into a correction set and a verification set in the modeling process, the correction set is 70 percent and 56 samples, and the verification set is 30 percent and 24 samples; selecting a correction set and a verification set sample, and dividing according to a chemical reference value range by combining a KS method;
(6) and (3) evaluating a model: using the coefficient of determination (R)2) The Root Mean Square Error (RMSE) and the relative analysis error (RPD) are used as precision evaluation indexes to comprehensively judge the quality of the model; wherein R is2Evaluation of the fitness of the model, R2The higher the value, the better the fitting effect of the model, and the stronger the capability of the model to explain the dependent variable; stability of the RMSE evaluation model and predictability of the RPD evaluation model, and the higher the RPD value, the lower the RMSE and the better the model.
2. The method as claimed in claim 1, wherein the spectrometer used in step (2) is a Tensor II Fourier transform mid-infrared spectrometer developed and produced by Bruker, Germany; measuring two spectral data of each soil sample, and taking the data after arithmetic mean as the original spectral data of the soil sample; the data storage uses the OPUS format.
3. The method according to claim 1, wherein the class 3 spectrum preprocessing method in the step (4) specifically comprises: the baseline correction comprises a first derivative and a second derivative, and mainly deducts the influence of the background or drift of the instrument on the spectral signal; the scattering correction comprises Multivariate Scattering Correction (MSC) and standard normal variable transformation (SNV), and mainly eliminates the influence of scattering on the spectrum caused by uneven particle distribution and different particle sizes; the scale scaling comprises maximum and minimum normalization and constant offset elimination, and mainly eliminates adverse effects caused by overlarge scale difference; and comparing the model after the spectrum pretreatment with the model without the pretreatment, and selecting the optimal spectrum pretreatment method of each soil carbon component from the model.
4. The method according to claim 1, wherein the modeling method in step (5) employs a PLSR method based on a factor analysis theory; the method has the relationship between a spectrum information point matrix and a property concentration matrix to be measured, obtains the optimal principal component factor or latent variable by performing factor analysis on the spectrum matrix and a target component matrix, and selects the optimal number of the latent variables to establish an optimal regression model according to the accumulated contribution rate of the latent variables to the two matrixes.
5. The method of claim 1, wherein the RPD in step (6) evaluates the predictability of the model, uses the RPD values for quality classification of the model: when the RPD is more than or equal to 3.0, the calibration effect is good, and the established calibration model can be used for actual detection; when the RPD is more than 2.5 and less than 3.0, quantitative analysis is feasible, but the prediction precision is to be improved; when RPD is less than or equal to 2.5, it is difficult to use for quantitative analysis.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449425.6A CN111595806A (en) | 2020-05-25 | 2020-05-25 | Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449425.6A CN111595806A (en) | 2020-05-25 | 2020-05-25 | Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111595806A true CN111595806A (en) | 2020-08-28 |
Family
ID=72187748
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010449425.6A Pending CN111595806A (en) | 2020-05-25 | 2020-05-25 | Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111595806A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436153A (en) * | 2021-06-11 | 2021-09-24 | 中国科学院南京土壤研究所 | Method for predicting carbon components of undisturbed soil profile based on hyperspectral imaging and support vector machine technology |
CN117491343A (en) * | 2023-12-25 | 2024-02-02 | 北京市农林科学院智能装备技术研究中心 | Soil organic matter detection method, device and system |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101738383A (en) * | 2008-11-07 | 2010-06-16 | 中国农业科学院农业环境与可持续发展研究所 | Mid-infrared spectrum-based method for quickly detecting contents of heavy metal elements in soil |
CN101876632A (en) * | 2009-04-28 | 2010-11-03 | 中国科学院沈阳应用生态研究所 | Method for measuring content of carbon and nitrogen in soil |
CN102798607A (en) * | 2012-08-13 | 2012-11-28 | 浙江大学 | Method for estimating soil organic carbon content by using mid-infrared spectrum technology |
CN104596957A (en) * | 2015-01-12 | 2015-05-06 | 西安科技大学 | Estimation method for content of copper in soil on basis of visible-light near-infrared spectrum technology |
CN105699314A (en) * | 2016-02-29 | 2016-06-22 | 上海交通大学 | Method for detecting soil stable carbon isotope ratio through intermediate infrared spectrum |
CN105784628A (en) * | 2016-02-29 | 2016-07-20 | 上海交通大学 | Method for detecting chemical composition of soil organic matter with mid-infrared spectra |
CN105784629A (en) * | 2016-02-29 | 2016-07-20 | 上海交通大学 | Method for quick detection of stable carbon isotope ratio of soil by means of mid-infrared spectroscopy |
CN106644978A (en) * | 2016-11-16 | 2017-05-10 | 山东省科学院海洋仪器仪表研究所 | Detection method capable of judging precision and used for analyzing content of soil nutrients based on spectral characteristic wavelength |
CN110531054A (en) * | 2019-09-29 | 2019-12-03 | 河南省农业科学院农业经济与信息研究所 | Soil organic matter uncertainty in traffic estimating and measuring method based on Bootstrap sampling |
-
2020
- 2020-05-25 CN CN202010449425.6A patent/CN111595806A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101738383A (en) * | 2008-11-07 | 2010-06-16 | 中国农业科学院农业环境与可持续发展研究所 | Mid-infrared spectrum-based method for quickly detecting contents of heavy metal elements in soil |
CN101876632A (en) * | 2009-04-28 | 2010-11-03 | 中国科学院沈阳应用生态研究所 | Method for measuring content of carbon and nitrogen in soil |
CN102798607A (en) * | 2012-08-13 | 2012-11-28 | 浙江大学 | Method for estimating soil organic carbon content by using mid-infrared spectrum technology |
CN104596957A (en) * | 2015-01-12 | 2015-05-06 | 西安科技大学 | Estimation method for content of copper in soil on basis of visible-light near-infrared spectrum technology |
CN105699314A (en) * | 2016-02-29 | 2016-06-22 | 上海交通大学 | Method for detecting soil stable carbon isotope ratio through intermediate infrared spectrum |
CN105784628A (en) * | 2016-02-29 | 2016-07-20 | 上海交通大学 | Method for detecting chemical composition of soil organic matter with mid-infrared spectra |
CN105784629A (en) * | 2016-02-29 | 2016-07-20 | 上海交通大学 | Method for quick detection of stable carbon isotope ratio of soil by means of mid-infrared spectroscopy |
CN106644978A (en) * | 2016-11-16 | 2017-05-10 | 山东省科学院海洋仪器仪表研究所 | Detection method capable of judging precision and used for analyzing content of soil nutrients based on spectral characteristic wavelength |
CN110531054A (en) * | 2019-09-29 | 2019-12-03 | 河南省农业科学院农业经济与信息研究所 | Soil organic matter uncertainty in traffic estimating and measuring method based on Bootstrap sampling |
Non-Patent Citations (6)
Title |
---|
C.GRINAND 等: "Prediction of soil organic and inorganic carbon contents", 《EUROPEAN JOURNAL OF SOIL SCIENCE》 * |
J.A.BAKDICK 等: "Predicting contents of carbon and its component fractions", 《SOIL RESEARCH》 * |
NUWAN K. WIJEWARDANE 等: "Predicting Physical and Chemical Properties of US Soils", 《SOIL SCIENCE SOCIETY OF AMERICA JOURNAL》 * |
丁建军 等: "基于PLSR 的土壤有机碳预测模型中建模组与验证组最优数量关系", 《汉江大学学报(自然科学版)》 * |
沈掌泉 等: "基于变量选择的偏最小二乘回归法和田间行走式近红外光谱进行土壤碳含量测定研究", 《光谱学与光谱分析》 * |
第五鹏瑶 等: "光谱预处理方法选择研究", 《光谱学与光谱分析》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113436153A (en) * | 2021-06-11 | 2021-09-24 | 中国科学院南京土壤研究所 | Method for predicting carbon components of undisturbed soil profile based on hyperspectral imaging and support vector machine technology |
CN113436153B (en) * | 2021-06-11 | 2023-07-14 | 中国科学院南京土壤研究所 | Undisturbed soil profile carbon component prediction method based on hyperspectral imaging and support vector machine technology |
CN117491343A (en) * | 2023-12-25 | 2024-02-02 | 北京市农林科学院智能装备技术研究中心 | Soil organic matter detection method, device and system |
CN117491343B (en) * | 2023-12-25 | 2024-04-16 | 北京市农林科学院智能装备技术研究中心 | Soil organic matter detection method, device and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107796764B (en) | Method for constructing wheat leaf area index estimation model based on three-band vegetation index | |
CN101806730B (en) | Vinegar residue organic matrix moisture content detection method | |
CN108801934A (en) | A kind of modeling method of soil organic carbon EO-1 hyperion prediction model | |
CN114018833B (en) | Method for estimating heavy metal content of soil based on hyperspectral remote sensing technology | |
CN104897592A (en) | Monitoring method of salt ion content in saline soil based on hyperspectral technology | |
CN110531054B (en) | Soil organic carbon prediction uncertainty estimation method based on Bootstrap sampling | |
Ren et al. | Estimating senesced biomass of desert steppe in Inner Mongolia using field spectrometric data | |
CN101231274B (en) | Method for rapid measuring allantoin content in yam using near infrared spectrum | |
CN111488926B (en) | Soil organic matter determination method based on optimization model | |
CN110687053B (en) | Regional organic matter content estimation method and device based on hyperspectral image | |
CN113436153B (en) | Undisturbed soil profile carbon component prediction method based on hyperspectral imaging and support vector machine technology | |
CN105486655A (en) | Rapid detection method for organic matters in soil based on infrared spectroscopic intelligent identification model | |
CN111965140B (en) | Wavelength point recombination method based on characteristic peak | |
Gholizadeh et al. | Models for estimating the physical properties of paddy soil using visible and near infrared reflectance spectroscopy | |
CN114817825B (en) | Emission source CO based on hyperspectral remote sensing 2 Quick imaging method | |
CN113111799B (en) | Cultivated land soil fertility level monitoring method based on ensemble empirical mode decomposition | |
Lulu et al. | Spectroscopy-based soil organic matter estimation in brown forest soil areas of the Shandong Peninsula, China | |
CN111595806A (en) | Method for monitoring soil carbon component by using mid-infrared diffuse reflection spectrum | |
CN104596943A (en) | Indoor spectral layered measuring method for poisonous and harmful elements at mine reclamation area | |
CN107505179A (en) | A kind of soil pretreatment and nutrient near infrared spectrum detection method | |
CN114169165B (en) | Model method for estimating rice chlorophyll under cadmium stress by three-band vegetation index | |
Hu et al. | Soil phosphorus and potassium estimation by reflectance spectroscopy | |
CN113971989A (en) | Forest soil organic carbon content high-spectrum modeling method based on OPLS | |
Chen et al. | Estimating soil organic carbon density in Northern China’s agro-pastoral ecotone using vis-NIR spectroscopy | |
Jia et al. | Prediction of soil organic carbon contents in Tibet using a visible near-infrared spectral library |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200828 |