CN108375553A - A kind of soil tenor detection method - Google Patents
A kind of soil tenor detection method Download PDFInfo
- Publication number
- CN108375553A CN108375553A CN201810130555.6A CN201810130555A CN108375553A CN 108375553 A CN108375553 A CN 108375553A CN 201810130555 A CN201810130555 A CN 201810130555A CN 108375553 A CN108375553 A CN 108375553A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- vector
- soil
- signal
- sample
- 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
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
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
The present invention provides a kind of soil tenor detection methods, which employs spectral analysis techniques, and provide the specific method of soil sampling, and a kind of data processing method of synthesis is proposed according to the characteristics of spectrum analysis, so that the prediction effect not only for the soil organism and total nitrogen content is preferable, but also the monitoring suitable for soil metal component content, also basicly stable good for the prediction effect of the metal components content such as phosphorus, potassium, calcium, magnesium.
Description
Technical field
The invention belongs to field of spectral analysis technology, and in particular to a kind of using spectral analysis technique soil tenor
Detection method.
Background technology
It is especially main to two kinds of soil with organic matter and full nitrogen at present at home and abroad about in soil constituent study on monitoring
Component content research in, it is most of to have obtained preferable achievement in research.On this basis, Many researchers also use light
Some metal element contents in spectral analysis technology quantitative analysis soil, it is pre- although also obtaining certain achievement for monitoring
Effect is surveyed but not as good as apparent to the prediction effect of organic matter and total nitrogen content.
That attempt spectral technique being introduced into Soil K+adsorption for the first time is American scientist Browsers and Hanks, 1965
Find that the organic matter in soil has relevant feature line mark near infrared spectrum region and organic matter year.Then in 1970, with
Afterwards in 1970, Salisbury and Hunt are equally found that certain minerals in soil near infrared spectrum there is also phases
The feature line mark of pass.The achievement of the above several researchers has started the new method detected for soil constituent, while in theory
On demonstrate using spectral analysis technique in soil Cucumber or ingredient carry out qualitative or quantitative analysis possibility.
At home and abroad in the result of study of scholar it is recognised that at present for soil constituent detection in, organic matter and
The prediction effect of full nitrogen is generally forged a good relationship, and, that has not consistent in each document report for the prediction effect of metallic element
The prediction effect of other scholar is preferable, and some is then not fully up to expectations, and at every in report, the processing for soil sampling object
It is not quite similar, some researchers the operations such as can be dried, be ground up, sieved to collected sample early period in experiment, experiment
As a result show that the prediction effect of spectrum analysis can be improved by carrying out processing appropriate to soil sampling object.Report proves sample
The physical properties such as water content, particle diameter precision of prediction can all be had an impact, and testing conditions such as sample whether rotate,
The factors such as detection height, angle also can all have an impact prediction effect.But the experiment conducted in having been reported that is adopted close red
External spectrum instrument carries out, and tests and carry out in laboratory environments mostly, rarely has the scene of field progress in the wild is practical to measure
Report.
When analyzing soil constituent using spectral analysis technique, the performance of spectral instrument has a fixing for prediction effect
It rings, and different prediction model method for building up has also embodied extremely important effect during prediction.Even if adopting sometimes
With the spectral instrument of same model, the method that prediction model is established is different, and obtained detection result is also different, wherein mainly adopting
With the methods of principle component regression, offset minimum binary and artificial neural network.Also have many scholars for a variety of data at
Reason method is compared, and as a result shows that each method has the characteristics that itself, but which kind of specifically used data processing method into
The testing result that row modeling can most be had, does not obtain an exact conclusion in current achievement in research.
Li Jie of China Agricultural Machinery Science Research Academy in 2012 et al. is complete to the soil sampling object near 72 Beijing
Nitrogen, full potassium, the content of organic matter and pH value are analyzed, and are modeled by Partial Least Squares, prediction result and actual value
With good consistency, the coefficient of determination reaches as high as 0.9554.What of Zhejiang University brave et al. has studied variable grain diameter
With different modeling methods for the image of soil organic matter content near-infrared spectral analytical method, sample is divided into particle in research
A diameter of 0.169-2mm and be less than two sample sets of 0.169mm, establish three kinds of prediction models, the results showed that work as particle diameter
When between 0.169-2mm, prediction related coefficient is all higher than 0.84, and predicted root mean square error is respectively less than 0.20；And when particle is straight
When diameter is less than 0.169mm, prediction related coefficient is respectively less than 0.71, and predicted root mean square error is all higher than 0.23.
From the point of view of development present situation, spectra methods is generally defined as a kind of completely new soil constituent inspection at present
Quick, the loseless method surveyed, and new development trend from now on, for solving the traditional chemical in soil constituent detection at present
The efficiency of detection method is low, it is of high cost, can not the drawbacks such as Fields detection.According to current achievement in research, spectral analysis technique for
The prediction effect of the soil organism and total nitrogen content is generally preferable, and the related coefficient of prediction result and chemical detection can generally reach
0.9 or more is even higher.The size of diameter of soil pellet has a certain impact for prediction effect, but less, soil moisture content
Also prediction result can be produced bigger effect, influence of the data processing method used by prediction mode to prediction result is little.
And generally it is not so good as the prediction effect of organic matter and total nitrogen content for the prediction effect of the metal components content such as phosphorus, potassium, calcium, magnesium
It is good, but can also reach 0.8 in some studies, what effect was relatively good in individual researchs can reach 0.9 or more, this just needs building
Suitable data processing method is selected when vertical prediction model.
Invention content
In view of the above analysis, the main purpose of the present invention is to provide a kind of to be also applied for suitable for soil in water quality
The detection method of tenor, which employs spectral analysis techniques, and provide the specific method of soil sampling, and according to
The characteristics of spectrum analysis, proposes a kind of data processing method of synthesis so that not only for the soil organism and total nitrogen content
Prediction effect is preferable, but also the monitoring suitable for soil metal component content, for the metal components content such as phosphorus, potassium, calcium, magnesium
Prediction effect it is also basicly stable good.
The purpose of the present invention is what is be achieved through the following technical solutions.
A kind of detection method being also applied for the tenor in water quality suitable for soil, includes the following steps：
Soil sampling is carried out to selection area；Soil sampling object is handled；Spectral scan is carried out to soil sampling object；
Spectral signal is handled.
It is further, described that carry out soil sampling to selection area include difference meter by variance yields and average value and true value
Number of samples P, and the natural number that P is >=19 are calculated, according to the size of sampling area, using snakelike sampling, sampling depth choosing
It is selected as 30CM, is sampled according to calculated number of samples, when sampling, a sampling section is first produced, is being parallel to section
It is sampled, adopts the earthwork to perpendicular to ground, the soil after sampling is mixed, then mixing thickness soil sampling object is done
It is broken, removal stone, plant big root system, be sufficiently mixed and transform into quadrangle, divide diagonal line and be divided into four parts, it is therein right to take
Two parts of angle, after being repeated a number of times, by the Weight control of soil sampling object needs weight.
Further, it is described to soil sampling object carry out processing include, by sample contents through hydrochloric acid, nitric acid, perchloric acid carry out
Processing.
Further, soil sampling object is positioned in sample container, gentle agitation makes its upper surface flat distribution, then
Sample is positioned on quartz window, in measurement process, soil sampling object can be rotated with sample container instrument to obtain repeatedly
The rotary speed of the averaged spectrum obtained after rotation sweep analysis, sample container is 4cm/s, and soil sampling object can quilt in 40s
Scanning 64 this, spectrum after arithmetic mean is as a sampling spectrum.Each sample is positioned in different sample containers, measures two
Secondary, the spectroscopic data after being averaged again is used for signal processing.
Further, further include absorption spectrum aligning step, including：Reference signal is acquired, measured signal is acquired, is obtained former
Beginning absorption spectrum determines correction reference data wavelength, utilizes correction reference data wavelength calibration absorption spectrum.
Further, described to carry out processing to spectral signal include spectrum sample at equal intervals, carries out two-dimensional wavelet transformation,
Then matrix of wavelet coefficients stationary window transversal scanning is carried out, the estimation of nonstationary noise variance evaluation vector sum signal variance is obtained
Vector adjusts window width, then carries out rescan point by point, obtains accurate noise variance evaluation vector sum precise signal
Variance evaluation vector, then shrinks denoising by Bayes's threshold value, finally carries out 2-d wavelet contravariant and exchanges central row vector work for
For output valve.
Further, further include being removed background interference and noise reduction process to spectral signal.
Further, it after the spectral signal that obtains that treated, is compared with the spectral signal threshold value of setting, when beyond threshold
It when value, alarms, reaches tenor testing goal.
Technical scheme of the present invention has the following advantages：
A kind of soil tenor detection method of the present invention, and synthesis proposes soil sampling and subsequent spectral signal
Processing method, multi task process especially has been carried out to spectral signal, and interfere it and noise has carried out fusion treatment, made
The spectral signal value that must be obtained disclosure satisfy that prediction effect, and confidence level is higher, overcome in the prior art to soil tenor
Prediction effect is not good, cannot achieve the defect of soil tenor monitoring.
Specific implementation mode
The soil tenor detection method of the present invention determines number of samples before specifically including sampling, is by variance yields
With the mathematic interpolation number of samples P of average value and true value, and P be >=19 natural number adopted according to the size of sampling area
With snakelike sampling, sampling depth is selected as 30CM, is sampled according to calculated number of samples, when sampling, first produces one
Section is sampled, is sampled being parallel to section, is adopted the earthwork to perpendicular to ground, the soil after sampling is mixed, then
Mixing thickness soil sampling object is broken into pieces, removal stone, plant big root system, be sufficiently mixed and transform into quadrangle, divide pair
Linea angulata is divided into four parts, diagonal two parts therein is taken, after being repeated a number of times, by the Weight control of soil sampling object in needs
Weight.
Then to handling sample contents through hydrochloric acid, nitric acid, perchloric acid.
Further, soil sampling object is positioned in sample container, gentle agitation makes its upper surface flat distribution, then
Sample is positioned on quartz window, in measurement process, soil sampling object can be rotated with sample container instrument to obtain repeatedly
The rotary speed of the averaged spectrum obtained after rotation sweep analysis, sample container is 4cm/s, and soil sampling object can quilt in 40s
Scanning 64 this, spectrum after arithmetic mean is as a sampling spectrum.Each sample is positioned in different sample containers, measures two
Secondary, the spectroscopic data after being averaged again is used for signal processing.
Further, further include absorption spectrum aligning step, including：Reference signal is acquired, measured signal is acquired, is obtained former
Beginning absorption spectrum determines correction reference data wavelength, utilizes correction reference data wavelength calibration absorption spectrum.
The correction reference data wavelength by being calculated as follows：
Spectrum matrix is denoted as X (M*K), wherein i-th of wavelength variable represents i-th of column vector in spectrum matrix X, remembers N=
Min (m-1, k), the n=1 before carrying out first time interative computation, the arbitrary a certain row for selecting spectrum matrix, are denoted as variable xk(0)；
The column vector location sets not being selected into are denoted as T,
Calculate separately xk(n-1)To remaining the projection of following variables, it is denoted as：
Note k (n-1)=arg (max (| | Pj||),j∈T)；N=n+1 enables x if n=Nk(n-1)As the row being selected into
Variable, returns to above-mentioned column vector position and calculates step and calculated again, using multiple linear regression method to the variables set that newly selects
Close { xk(0),xk(1),…,xk(n-1)Assessed, obtain correction reference data wavelength value.
Further, after obtaining spectral signal, described to carry out processing to spectral signal include spectrum sample at equal intervals, into
Row two-dimensional wavelet transformation, then carry out matrix of wavelet coefficients stationary window transversal scanning, obtain nonstationary noise variance evaluation to
Amount and signal variance estimate vector, adjust window width, then carry out rescan point by point, obtain accurate noise variance and estimate
Vector sum precise signal variance evaluation vector is counted, denoising is then shunk by Bayes's threshold value, finally carries out 2-d wavelet contravariant
Center row vector is exchanged for as signal value output.
After the spectral scan output valve for obtaining soil sampling object, have compared with high s/n ratio to remove background interference and obtain
Soil spectrum data are needed with improving the robustness and predictive ability of the soil constituent prediction model that the later stage is established to data
It is handled as follows：
Calculate the spectrum average of sample contents：
To unknown sample spectrum x (1 × m), pass through the following formula spectrum x that obtains that treatedcentered：
Then spectrum standard deviation is calculated:
Spectrum sample value is handled as follows：
Entire spectrum sample signal value is handled as follows, value X after being handled：
Wherein, n is sample number, is counted for wavelength, and the spectrum that Xn corresponds to the sample contents as n-th sample is equal
Value；
The spectral signal of sampling has a degree of inhomogeneities, and it is that significant non-stationary property is made an uproar to cause spectral noise
Sound, spectral noise matrix are the algebraical sums of the notable nonstationary noise on short-term stationarity noise and spectrum axis on a timeline.
Generally use wavelet transformation carries out denoising, however the Wavelet-denoising Method of single threshold value be difficult to cope with it is notable on spectrum axis
Nonstationary noise is easily lost the signal detail of low-noise area.To containing nonstationary noise image B, M layers points of small echo is carried out
Solution.The window W that an area is m*n is arranged to slide point by point in matrix of wavelet coefficients, the denoising of window center point wavelet coefficient
Threshold value only has the wavelet coefficient in window and is calculated, and the Spatial distributions estimation of noise-removed threshold value may be implemented.The area of W, by passing through
Value setting is tested, it is too small, then it can not utilize statistical law Accurate Prediction noise variance；It is too big, then it is unable to the fast of response noises variance
Speed variation.Therefore, it when dealing with objects different, needs to continuously attempt to its area, noise variance changes small image, selection
Larger window, conversely, then selection is compared with wicket.
But in this present embodiment, as a result of the soil collection and processing mode of the present embodiment, and being adopted to soil
For sample object, spectral noise concentrates on specific region, and noise variance is big in region and variation is strong, and rest part noise variance
Relatively small and steady, fixed window area is easy that threshold value estimation is made to generate error.Therefore, applicant proposed a kind of windows
The point-by-point adjustment mode of width dynamic carries out noise reduction process to spectral signal matrix.
The Noise Variance Estimation of window is fixed first, defines stationary window Wm×n, width m is set as empirical value, height
N is coefficient matrix height, by window width diagonally high frequency coefficient matrix center row vectorOn slide laterally point by point,
WhereinThe noise variance vector that central row is estimated using wavelet coefficient in window, according to the noise
The change rate of variance vectors, dynamic adjust window width, and specially according to noise, approaches uniformity is distributed in all wavelet coefficients
Property and wavelet basis linear phase property, the noise variance vector of high-rise wavelet coefficient can be by the noise side of last layer
Two neighboring variance is averaging to obtain in difference vector, signal variance vector is estimated, evaluation method and above-mentioned noise
Variance vectors evaluation method is similar, but since signal variance concentrates in the wavelet coefficient of part, has on each wavelet decomposition layer larger
Difference needs to carry out layering estimation, is slid laterally respectively in the center row vector of all directions coefficient matrix of the layer calculated,
Using all wavelet coefficients respectively decomposed in this layer of window on direction, the signal variance vector of this layer is estimated, according to signal variance
The change rate dynamic setting of vector, resets serial ports width, and setting signal window is adjusted threshold value, estimated again using new window
Every layer of signal variance vector, then the threshold vector of every layer of matrix of wavelet coefficients is acquired, for high frequency transverse direction coefficient matrix, longitudinal direction
Every a line of high frequency coefficient matrix, diagonal high frequency coefficient matrix, carries out contraction noise reduction using threshold vector, will be small after noise reduction
Spectrum matrix is changed in wave system number contravariant, takes the center row vector of spectrum matrix as the spectral signal after noise reduction.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
All any modification, equivalent and improvement etc., should all be included in the protection scope of the present invention made by within refreshing and principle.
Claims (2)
1. a kind of soil tenor detection method comprising carry out soil sampling to selection area；Soil sampling object is carried out
Processing；Spectral scan is carried out to soil sampling object；Spectral signal is handled；By treated, signal is carried out with given threshold
Compare, and alarms the signal value more than threshold value.
2. a kind of soil tenor detection method as claimed in claim 1, specifically includes and determines number of samples P before sampling,
It is the mathematic interpolation number of samples by variance yields and average value and true value, and P is >=19 natural number；According to the face of sampling area
Product size, using snakelike sampling, sampling depth is selected as 30CM, is sampled according to calculated number of samples, when sampling, first
One sampling section of extraction, is sampled being parallel to section, adopts the earthwork to perpendicular to ground, the soil after sampling is blended in
Together, then by mixing thickness soil sampling object it is crushed, removal stone, root system of plant mix and transform into quadrangle, divide diagonal line
Be divided into four parts, take diagonal two parts therein, after being repeated a number of times, by the Weight control of soil sampling object needs weight
Amount, then to handling sample contents through hydrochloric acid, nitric acid, perchloric acid；
Soil sampling object is positioned in sample container, shake is distributed with making its upper surface substantially flat, then places sample
In on quartz window, in measurement process, soil sampling object can be rotated with sample container instrument to obtain multiple rotary scanning point
The rotary speed of the averaged spectrum obtained after analysis, sample container is 4cm/s, and soil sampling object can be scanned 64 times in 40s, is calculated
Spectrum after number is average is as a sampling spectrum；Each sample is positioned in different sample containers, is measured twice, average again
Spectroscopic data afterwards is used for signal processing；
The wherein described signal processing further includes：Reference signal is acquired, measured signal is acquired, original absorbance spectrum is obtained, determines school
Positive reference data wavelength utilizes correction reference data wavelength calibration absorption spectrum；The correction reference data wavelength passes through as follows
It is calculated：
Spectrum matrix is denoted as X (M*K), wherein i-th of wavelength variable represents i-th of column vector in spectrum matrix X, remembers N=min
(m-1, k), the n=1 before carrying out first time interative computation, the arbitrary a certain row for selecting spectrum matrix are denoted as variable xk(0)；Not having
There are the column vector location sets being selected into be denoted as T,
Calculate separately xk(n-1)To remaining the projection of following variables, it is denoted as：
Note k (n-1)=arg (max (| | Pj||),j∈T)；N=n+1 enables x if n=Nk(n-1)As the row variable being selected into,
It returns to above-mentioned column vector position and calculates step and calculated again, using multiple linear regression method to the variables collection that newly selects
{xk(0),xk(1),…,xk(n-1)Assessed, obtain correction reference data wavelength value；
After the spectral scan output valve for obtaining soil sampling object, in order to remove the soil spectrum that background interference obtains more high s/n ratio
Data need the data to P sample to improve the robustness and predictive ability of the soil constituent prediction model that the later stage is established
It is handled as follows：
Calculate the spectrum average of sample contents：
To unknown sample spectrum x (1 × m), pass through the following formula spectrum x that obtains that treatedcentered：
Then spectrum standard deviation is calculated:
Spectrum sample value is handled as follows：
Entire spectrum sample signal value is handled as follows, value X after being handled：
Wherein, n is sample number, is counted for wavelength, and Xn corresponds to the spectrum average of the sample contents as n-th sample；
After obtaining spectral signal, the spectral signal processing further includes that spectrum sample, carries out two-dimensional wavelet transformation, connect at equal intervals
Carry out matrix of wavelet coefficients stationary window transversal scanning, obtain nonstationary noise variance evaluation vector sum signal variance estimate to
Amount, adjusts window width, then carry out rescan point by point, obtains accurate noise variance evaluation vector sum precise signal side
Then poor estimate vector shrinks denoising by Bayes's threshold value, finally carry out 2-d wavelet contravariant and exchange central row vector conduct for
Signal value output；
It is described to include to the point-by-point adjustment of window width progress：The Noise Variance Estimation of window is fixed first, defines fixed window
Mouth Wm×n, width m is set as empirical value, and height n is coefficient matrix height, by window width diagonally high frequency coefficient matrix center
Row vectorOn slide laterally point by point, whereinCentral row is estimated using wavelet coefficient in window
Noise variance vector, according to the change rate of noise variance vector, dynamic adjustment window width specially exists according to noise
The linear phase property of approaches uniformity is distributed in all wavelet coefficients property and wavelet basis, the noise of high-rise wavelet coefficient
Variance vectors can carry out signal variance vector by being averaging to obtain by two neighboring variance in the noise variance vector of last layer
Estimation, evaluation method is similar with above-mentioned noise variance vector evaluation method, but since signal variance concentrates on part wavelet systems
In number, there is larger difference on each wavelet decomposition layer, need to carry out layering estimation, in all directions coefficient matrix of the layer calculated
It is slid laterally respectively in the row vector of center, using all wavelet coefficients respectively decomposed in this layer of window on direction, estimates this layer
Signal variance vector, is set dynamically according to the change rate of signal variance vector, resets serial ports width, setting signal window tune
Whole threshold value, every layer of signal variance vector is estimated using new window again, then acquires the threshold vector of every layer of matrix of wavelet coefficients,
For every a line of high frequency transverse direction coefficient matrix, longitudinal high frequency coefficient matrix, diagonal high frequency coefficient matrix, threshold vector is used
Contraction noise reduction is carried out, spectrum matrix is changed into the wavelet coefficient contravariant after noise reduction, takes the center row vector of spectrum matrix as drop
Spectral signal after making an uproar.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810130555.6A CN108375553A (en) | 2018-02-08 | 2018-02-08 | A kind of soil tenor detection method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810130555.6A CN108375553A (en) | 2018-02-08 | 2018-02-08 | A kind of soil tenor detection method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108375553A true CN108375553A (en) | 2018-08-07 |
Family
ID=63017453
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810130555.6A Pending CN108375553A (en) | 2018-02-08 | 2018-02-08 | A kind of soil tenor detection method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108375553A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115508292A (en) * | 2022-09-20 | 2022-12-23 | 中国科学院南京土壤研究所 | Soil profile nitrogen content high spectrum detection and visualization method based on machine learning |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107462542A (en) * | 2017-09-19 | 2017-12-12 | 黑龙江省能源环境研究院 | A kind of heavy metal-polluted soil detection warning device and its application method |
-
2018
- 2018-02-08 CN CN201810130555.6A patent/CN108375553A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107462542A (en) * | 2017-09-19 | 2017-12-12 | 黑龙江省能源环境研究院 | A kind of heavy metal-polluted soil detection warning device and its application method |
Non-Patent Citations (3)
| Title |
|---|
| 吴德操: "面向地表水环境的分布式紫外-可见光谱水质在线检测监测系统关键技术研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 尚廷义 等: "《基于光谱技术的寒地水稻稻瘟病害分析及机理研究》", 30 June 2016, 哈尔滨工程大学出版社 * |
| 李颉: "土壤养分的近红外光谱快速分析系统开发与实验研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115508292A (en) * | 2022-09-20 | 2022-12-23 | 中国科学院南京土壤研究所 | Soil profile nitrogen content high spectrum detection and visualization method based on machine learning |
| CN115508292B (en) * | 2022-09-20 | 2024-05-28 | 中国科学院南京土壤研究所 | Machine learning-based hyperspectral detection and visualization method for nitrogen content of soil profile |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nawar et al. | Optimal sample selection for measurement of soil organic carbon using on-line vis-NIR spectroscopy | |
| Wang et al. | Modeling salinity effects on soil reflectance under various moisture conditions and its inverse application: A laboratory experiment | |
| Monakhova et al. | Independent components in spectroscopic analysis of complex mixtures | |
| CN110531054B (en) | Soil organic carbon prediction uncertainty estimation method based on Bootstrap sampling | |
| CN105628581B (en) | A kind of tight sandstone reservoir based on hyperspectral technique is appeared porosity characterizing method | |
| CN106645037A (en) | Method for detecting heavy metal content of coal gangue filling reclamation reconstruction soil based on high spectrum technology | |
| CN110455722A (en) | Rubber tree blade phosphorus content EO-1 hyperion inversion method and system | |
| CN112363161B (en) | Vegetation vertical structure and under-forest topography inversion method and device based on scattering mechanism decomposition | |
| Zhang et al. | Quantitative estimating salt content of saline soil using laboratory hyperspectral data treated by fractional derivative | |
| CN108344691A (en) | Soil contamination grade detecting system | |
| Xiong et al. | The total P estimation with hyper-spectrum–a novel insight into different P fractions | |
| Šestak et al. | Visible and near infrared reflectance spectroscopyfor field-scale assessment of stagnosols properties. | |
| Zhu et al. | A physics-based model for crystal orientation dictionary indexing by directional reflectance microscopy | |
| CN108375553A (en) | A kind of soil tenor detection method | |
| CN108362654A (en) | Solid waste prevents monitoring system | |
| CN108387539A (en) | Oilfield drilling waste disposes early warning equipment | |
| CN108344704A (en) | Consumer waste incineration handles monitoring device | |
| CN108362655A (en) | A kind of monitoring water quality on line system of more administrative grade linkages | |
| Zhuo et al. | Quantitative retrieving of soil organic matter using field spectrometer and hyperspectral remote sensing | |
| Liu et al. | Detection of Apple Taste Information Using Model Based on Hyperspectral Imaging and Electronic Tongue Data. | |
| CN108344705A (en) | Rubbish from cooking disposes monitoring method | |
| Zeng et al. | Hyperspectral reflectance models for soil salt content by filtering methods and waveband selection | |
| Kononov et al. | Model-associated forest parameter retrieval using VHF SAR data at the individual tree level | |
| CN108344692A (en) | Using the COD online automatic monitors of technology of Internet of things | |
| CN108362656A (en) | A kind of PH on-line computing models of networking |
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: 20180807 |
|
| RJ01 | Rejection of invention patent application after publication |