CN103900989B - Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method - Google Patents
Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method Download PDFInfo
- Publication number
- CN103900989B CN103900989B CN201410160522.8A CN201410160522A CN103900989B CN 103900989 B CN103900989 B CN 103900989B CN 201410160522 A CN201410160522 A CN 201410160522A CN 103900989 B CN103900989 B CN 103900989B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- curve
- saliferous clay
- data point
- infrared band
- 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.)
- Expired - Fee Related
Links
Abstract
The invention discloses saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, comprise the steps: a) to detect saliferous clay sample by the spectral instrument with particular sample interval, thus obtaining the curve of spectrum corresponding with saliferous clay sample, the curve of spectrum has multiple data point; B) first order derivative is carried out to curve of spectrum data point, second derivative operator; C) ask radius-of-curvature at data point place, obtain curve of spectrum minimum profile curvature radius and corresponding wavelength; D) calculate the change frequency of the curve of spectrum at corresponding wavelength place, obtain the maximum change frequency f of the curve of spectrum
max; E)
be optimal spectrum resolution, thus obtain reconstructing the curve of spectrum.Method provided by the invention overcomes the shortcoming such as data redundancy in prior art, also have simultaneously can undistorted preservation complete spectrum information, meet the advantages such as saliferous clay monitoring requirements.
Description
Technical field
The present invention relates to remote sensing and spectroscopy analysis field, particularly relate to a kind of saliferous clay curve of spectrum infrared band remote-sensing inversion construction method.
Background technology
Since 21 century, the land resource of various countries is more and more deficienter, and salinization land is more and more paid attention to as potential utilized soil, and becomes the focus of research by remote sensing technology monitoring salinization soil resource.The remote sensing monitoring of saliferous clay realizes based on the saliferous clay curve of spectrum.At present, instrument ASDProFR field spectroradiometer for hyperspectral measurement can obtain the salinization soil curve of spectrum of 350-2500nm wavelength coverage, its spectral resolution is 1.4nm at visible light wave range, be 2nm at near-infrared band, AvaField-2 field spectroradiometer, its spectral response range 300-1700nm, within the scope of 300-1040nm, sampling interval is 0.59nm, and within the scope of 1040-1700nm, sampling interval is 3.8nm etc.
Based on being closely related of the saliferous clay remote sensing monitoring precision of spectral signature and the spectral resolution of sensor.Usually adopt for improving accuracy of identification the method improving sensor accuracy, this often causes the meaningless increase of the waste of spectrum resource and information processing capacity, brings unnecessary trouble to the remote sensing monitoring of saliferous clay.The optimum resolution that curve of spectrum remote-sensing inversion builds refers to and can meet the sampling interval that monitoring requirements can not cause again data redundancy, and optimum resolution can ensure the undistorted monitoring realizing sensor information.The acquisition of the undistorted saliferous clay curve of spectrum is one of important content of saliferous clay Remote Sensing Information Extraction, the method obtaining the undistorted saliferous clay curve of spectrum is remote sensing technology study frontier and focus, and therefore research and development determine that the method for curve of spectrum remote-sensing inversion structure optimum resolution has great importance.
From information-theoretical angle: determine that essence that curve of spectrum remote-sensing inversion builds optimal spectrum resolution is that the sampling interval of undistorted reconstruction information (signal) is determined.
To carry out saliferous clay remote sensing monitoring Problems existing based on spectral signature be the demand how selecting to determine remote sensing monitoring optimal spectrum resolution for current, exploitation one can either meet remote sensing monitoring demand, can not data redundancy be caused again, can also the curve of spectrum remote-sensing inversion of undistorted preservation complete information to build the defining method of optimum resolution just of crucial importance.
Summary of the invention
Because the deficiencies in the prior art, technical matters to be solved by this invention is to provide a kind of saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, its can meet saliferous clay monitoring requirements and significantly can reduce that saliferous clay spectrum redundant resource takies, undistorted preservation complete information and saliferous clay can be met detect demand.
For achieving the above object, the invention provides a kind of saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, particularly, technical scheme provided by the invention is as follows:
A kind of saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, comprises the steps:
A) spectral instrument by having particular sample interval detects saliferous clay sample, thus the curve of spectrum that the reflectivity R obtaining wavelength X and saliferous clay sample changes, the curve of spectrum has multiple data point;
B) first order derivative and second derivative is calculated respectively at multiple data point place;
C) ask the radius-of-curvature of the curve of spectrum according to first order derivative and second derivative respectively at multiple data point place, compare multiple radius-of-curvature of trying to achieve, to obtain the minimum profile curvature radius r of the curve of spectrum
minwith the corresponding wavelength λ of minimum profile curvature radius
rmin;
D) λ is calculated
rminthe change frequency at place, obtains the maximum change frequency f of the curve of spectrum
max;
E) according to Shannon's sampling theorem, saliferous clay sample light spectral curve infrared band optimal spectrum resolution is determined
utilize described optimal spectrum resolution ax λ
oRobtain reconstructing the curve of spectrum.
Further, the wavelength selection range of the curve of spectrum is more than or equal to 400nm, is less than 1650nm.
Further, first order derivative is carried out to each data point of the curve of spectrum, second derivative operator adopts with the following method:
R′(λ
i)=[R(λ
i+1)-R(λ
i-1)]/(2Δλ);
R″(λ
i)=[R(λ
i+1)-2R(λ
i)+R(λ
i-1)]/(Δλ)
2;
λ
ifor the wavelength at data point i place, R (λ
i) be the reflectivity at data point i place, R ' (λ
i) and R " (λ
i) be respectively the curve of spectrum in wavelength X
ithe single order at place and second derivative, Δ λ is λ
i-1to λ
iwavelength interval.
Further, the computing method of radius-of-curvature are:
r=|[(1+R′
2)(
3/2)]/R″|
Further, maximum frequency change calculates as follows:
1. r is worked as
minwhen>=1, f
maxbe numerically equal to r
min, optimal spectrum resolution
2. r is worked as
minwhen≤1, f
maxbe numerically equal to 1/r
min, optimal spectrum resolution ax λ
oR=1/ (2r
min).
Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method provided by the invention to overcome in prior art data redundancy and causes the shortcomings such as the meaningless increase of information processing capacity, also have simultaneously can undistorted preservation complete information, optimum efficiency can be reached namely meet the advantages such as saliferous clay monitoring requirements to distinguishing the saliferous clay of different salt content, the remote sensing resources information of saliferous clay is more optimized efficiently.
Be described further below with reference to the technique effect of accompanying drawing to method of the present invention and generation, to understand object of the present invention, characteristic sum effect fully.
Accompanying drawing explanation
Fig. 1 is the infrared band curve of spectrum of the saliferous clay sample 1 of the embodiment of the present invention.
Fig. 2 is the saliferous clay sample 1 infrared band radius-of-curvature variation diagram of the embodiment of the present invention.
Fig. 3 is the saliferous clay sample 1 infrared band curve of spectrum reconstruct front and back comparison diagram of the embodiment of the present invention.
Fig. 4 is the infrared band curve of spectrum of the saliferous clay sample 2 of the embodiment of the present invention.
Fig. 5 is the saliferous clay sample 2 infrared band radius-of-curvature variation diagram of the embodiment of the present invention.
Fig. 6 is the saliferous clay sample 2 infrared band curve of spectrum reconstruct front and back comparison diagram of the embodiment of the present invention.
Embodiment
The present embodiment provides a kind of saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, comprises the steps:
A) spectral instrument by having particular sample interval detects saliferous clay sample, thus the curve of spectrum that the reflectivity R obtaining wavelength X and saliferous clay sample changes, the curve of spectrum has multiple data point i;
B) first order derivative and second derivative is calculated respectively at multiple data point place;
C) ask the radius-of-curvature of the curve of spectrum according to first order derivative and second derivative respectively at multiple data point place, compare multiple radius-of-curvature of trying to achieve, to obtain the minimum profile curvature radius r of the curve of spectrum
minwith the corresponding wavelength λ of minimum profile curvature radius
rmin;
D) λ is calculated
rminthe change frequency at place, obtains the maximum change frequency f of the curve of spectrum
max;
E) according to Shannon's sampling theorem, saliferous clay sample light spectral curve infrared band optimal spectrum resolution is determined
utilize described optimal spectrum resolution ax λ
oRobtain reconstructing the curve of spectrum.
In this preferred embodiment, institute's usage data is that AvaField-2 Portable field EO-1 hyperion instrument obtains, and the sensing range of instrument is more than or equal to 300nm, is less than 1700nm, wherein 300-1040nm subband samples is spaced apart 0.59nm, and 1040-1700nm subband samples is spaced apart 3.8nm.In other embodiments, also can use other spectrometer, the wavelength coverage of the spectral instrument used is more than or equal to 300nm, is less than 1700nm.General instrument is detecting the jitter at band edge place, and therefore the wavelength selection range of the curve of spectrum is preferably greater than and equals 400nm, is less than 1650nm, and the infrared band data that the present embodiment adopts are 760-1650nm.
1, for the infrared band spectrum of a saliferous clay sample 1, its curve of spectrum Z1 as shown in Figure 1.
1) minimum profile curvature radius of the saliferous clay curve of spectrum is calculated, the input saliferous clay curve of spectrum, its radius-of-curvature of node-by-node algorithm, try to achieve this curve minimum profile curvature radius R
minand position:
First order derivative and second derivative adopt following formulae discovery:
R′(λ
i)=[R(λ
i+
1)-R(λ
i-1)]/(2Δλ);
R″(λ
i)=[R(λ
i+1)-2R(λ
i)+R(λ
i-1)]/(Δλ)
2;
λ
ifor the wavelength at data point i place, R (λ
i) be the reflectivity at described data point i place, R ' (λ
i) and R " (λ
i) be respectively the curve of spectrum in wavelength X
ithe single order at place and second derivative, Δ λ is λ
i-1to λ
iwavelength interval.
Calculate each point radius-of-curvature, as shown in Figure 2.Obtaining least radius position is 812.01nm, and radius-of-curvature is 28.93194;
2) the maximum change frequency of the saliferous clay curve of spectrum is obtained:
The size of radius-of-curvature determines the size of curvilinear motion frequency, the less circle meaning that adjacent three points are formed of radius-of-curvature is less, and the change frequency of curve is also larger, this curve of undistorted reconstruct, then need less sampling interval, then just contrary when radius-of-curvature becomes large.
1. r is worked as
minwhen>=1, the maximum change frequency of curve is comparatively large, and the resolution that undistorted reconstruct Curves needs is also higher, now, and f
maxbe numerically equal to r
min, therefore undistorted remote-sensing inversion builds the optimal spectrum resolution of this target signature spectral profile
2. r is worked as
minwhen≤1, the maximum change frequency of curve is less, and the resolution that undistorted reconstruct Curves needs is also lower, now f
maxbe numerically equal to 1/r
min, therefore undistorted remote-sensing inversion builds the optimal spectrum resolution ax λ of this target signature spectral profile
oR=1/ (2r
min);
3) above-mentioned computing is done to the curve of spectrum Z1 in 760-1650nm wavelength band, minimum radius-of-curvature must be had at 812.01nm place, obtain according to above-mentioned theory
4) above-mentioned Δ λ is utilized
oR=14nm reconstructs the curve of spectrum.Adopt Δ λ
oRthe curve of spectrum that=14nm paints well contains the information of original spectrum curve, and identical with primary curve, all has an absorption peak at 1415nm wave band place, and reflectivity presents the trend risen gradually with wave band change.The wave band number of the curve of spectrum will be 63 from original 709, greatly reduce data volume, ensure that the curve of spectrum does not change simultaneously.As shown in Figure 3, the Δ λ utilizing said method to try to achieve
oRdraw the curve of spectrum, finally overlap completely with original spectrum curve.
2, with the infrared band of a saliferous clay sample 2 (near infrared, far infrared, in infrared) spectrum is example, its curve of spectrum Z2 is as shown in Figure 4.
1) minimum profile curvature radius of the saliferous clay curve of spectrum is calculated, the input saliferous clay curve of spectrum, its radius-of-curvature of node-by-node algorithm, try to achieve this curve minimum profile curvature radius R
minand position:
First order derivative and second derivative adopt following formulae discovery:
R′(λ
i)=[R(λ
i+1)-R(λ
i-1)]/(2Δλ);
R″(λ
i)=[R(λ
i+1)-2R(λ
i)+R(λ
i-1)]/(Δλ)
2;
λ
ifor the wavelength at data point i place, R (λ
i) be the reflectivity at described data point i place, R ' (λ
i) and R " (λ
i) be respectively the curve of spectrum in wavelength X
ithe single order at place and second derivative, Δ λ is λ
i-1to λ
iwavelength interval.
Calculate each point radius-of-curvature, as shown in Figure 5.Obtaining least radius position is 1034.02nm, and radius-of-curvature is 0.256612;
2) the maximum change frequency of the saliferous clay curve of spectrum is obtained:
1. r is worked as
minwhen>=1, f
maxbe numerically equal to r
min;
2. r is worked as
minwhen≤1, f
maxbe numerically equal to 1/r
min.
3) above-mentioned computing is done to the curve of spectrum Z2 in 760-1650nm wavelength band, minimum radius-of-curvature must be had at 1034.02nm place, obtain according to above-mentioned theory
4) above-mentioned Δ λ is utilized
oR=1.95nm reconstructs the curve of spectrum.Adopt Δ λ
oRthe curve of spectrum that=1.95nm paints well contains the information of original spectrum curve, and identical with primary curve, all has an absorption peak at 1415nm wave band place, and reflectivity presents the trend risen gradually with wave band change.The wave band number of the curve of spectrum will be 354 from original 709, greatly reduce data volume, ensure that the curve of spectrum does not change simultaneously.As shown in Figure 6, the Δ λ utilizing said method to try to achieve
oRdraw the curve of spectrum, finally overlap completely with original spectrum curve.
More than describe preferred embodiment of the present invention in detail.Should be appreciated that the ordinary skill of this area just design according to the present invention can make many modifications and variations without the need to creative work.Therefore, all technician in the art, all should by the determined protection domain of claims under this invention's idea on the basis of existing technology by the available technical scheme of logical analysis, reasoning, or a limited experiment.
Claims (4)
1. a saliferous clay curve of spectrum infrared band remote-sensing inversion construction method, is characterized in that comprising the steps:
A) spectral instrument by having particular sample interval detects saliferous clay sample, thus the curve of spectrum that the reflectivity R obtaining wavelength X and described saliferous clay sample changes, the described curve of spectrum has multiple data point;
B) first order derivative and second derivative is calculated respectively at described multiple data point place;
C) ask the radius-of-curvature of the described curve of spectrum according to described first order derivative and second derivative respectively at described multiple data point place, compare described multiple radius-of-curvature of trying to achieve, to obtain the minimum profile curvature radius r of the curve of spectrum
minwith the corresponding wavelength λ of described minimum profile curvature radius
rmin;
D) λ is calculated as follows
rminthe change frequency at place, obtains the maximum change frequency f of the described curve of spectrum
max:
1. r is worked as
minwhen>=1, f
maxbe numerically equal to r
min;
2. r is worked as
minwhen≤1, f
maxbe numerically equal to 1/r
min;
E) according to Shannon's sampling theorem, described saliferous clay sample light spectral curve infrared band optimal spectrum resolution is determined
utilize described optimal spectrum resolution ax λ
oRobtain reconstructing the curve of spectrum.
2. saliferous clay curve of spectrum infrared band remote-sensing inversion construction method as claimed in claim 1, the determined wavelength scope of wherein said spectral instrument is more than or equal to 300nm, is less than 1700nm.
3. saliferous clay curve of spectrum infrared band remote-sensing inversion construction method as claimed in claim 1, wherein the determined wavelength selection range of the curve of spectrum is more than or equal to 400nm, is less than 1650nm.
4. saliferous clay curve of spectrum infrared band remote-sensing inversion construction method as claimed in claim 1, wherein said first order derivative and described second derivative are calculated as follows:
R′(λ
i)=[R(λ
i+1)-R(λ
i-1)]/(2Δλ);
R″(λ
i)=[R(λ
i+1)-2R(λ
i)+R(λ
i-1)]/(Δλ)
2;
λ
ifor the wavelength at described data point i place, R (λ
i) be the reflectivity at described data point i place, R ' (λ
i) and R " (λ
i) be respectively the described curve of spectrum in wavelength X
ithe single order at place and second derivative, Δ λ is λ
i-1to λ
iwavelength interval.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410160522.8A CN103900989B (en) | 2014-04-21 | 2014-04-21 | Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410160522.8A CN103900989B (en) | 2014-04-21 | 2014-04-21 | Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103900989A CN103900989A (en) | 2014-07-02 |
CN103900989B true CN103900989B (en) | 2016-03-23 |
Family
ID=50992450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410160522.8A Expired - Fee Related CN103900989B (en) | 2014-04-21 | 2014-04-21 | Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103900989B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105320833A (en) * | 2014-08-01 | 2016-02-10 | 中国土地勘测规划院 | Method for determining optimum spectral resolution established by remote sensing inversion of spectrum character curve |
CN111220163A (en) * | 2018-11-26 | 2020-06-02 | 厦门雅迅网络股份有限公司 | Curve judgment method based on curvature radius calculation, terminal device and storage medium |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7433806B2 (en) * | 2001-04-20 | 2008-10-07 | Spectral Sciences, Inc. | Band model method for modeling atmospheric propagation at arbitrarily fine spectral resolution |
AU2002344807B2 (en) * | 2001-06-18 | 2005-12-08 | Noxell Corporation | Natural tone coloured cosmetic compositions and methods for providing such compositions |
JP2004247380A (en) * | 2003-02-12 | 2004-09-02 | Mitsubishi Electric Corp | Evaluation method of piezo electric field |
WO2007050289A1 (en) * | 2005-10-21 | 2007-05-03 | Raytheon Company | Sonar system and method providing low probability of impact on marine mammals |
CN101493407B (en) * | 2008-12-11 | 2011-02-09 | 中山大学 | Air pollution information extraction method based on earth-atmosphere coupling |
CN101799871B (en) * | 2010-01-18 | 2012-02-29 | 浙江林学院 | Regression parameter transformation based method for extracting thematic information of remote sensing images |
CN101853503A (en) * | 2010-04-26 | 2010-10-06 | 华中科技大学 | Spectral line inflexion multi-scale optimizing segmentation method and application thereof |
CN102645410A (en) * | 2012-04-28 | 2012-08-22 | 中国神华能源股份有限公司 | Detection method, device and data processing device of coal gangue pollution scope |
CN103048276B (en) * | 2012-12-14 | 2014-09-03 | 北京农业信息技术研究中心 | Spectral index constructing method for detecting carbon nitrogen ratios of canopy leaves of crops |
-
2014
- 2014-04-21 CN CN201410160522.8A patent/CN103900989B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103900989A (en) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lin et al. | Cloud structure of galactic OB cluster-forming regions from combining ground-and space-based bolometric observations | |
Tsai et al. | Derivative analysis of hyperspectral data | |
CN102799782B (en) | Electrical performance prediction method for metal truss-type antenna cover in consideration of scattering and transmission | |
CN103472008B (en) | Embryo Gallus domesticus gender identification method in hatching early stage near-infrared hatching egg | |
Liu et al. | Spherical harmonic analyses of intensity mapping power spectra | |
Ebadi et al. | A review of applying second-generation wavelets for noise removal from remote sensing data | |
Bai et al. | Parameter-searched OMP method for eliminating basis mismatch in space-time spectrum estimation | |
CN102631198A (en) | Dynamic spectrum data processing method based on difference value extraction | |
Verde et al. | On minimally parametric primordial power spectrum reconstruction and the evidence for a red tilt | |
CN103900989B (en) | Saliferous clay curve of spectrum infrared band remote-sensing inversion construction method | |
CN103456015A (en) | SAR target detection method based on optimal fractional domain Gabor spectrum features | |
Zhang et al. | Quantitative estimating salt content of saline soil using laboratory hyperspectral data treated by fractional derivative | |
CN110954510A (en) | Nano plasma spectrum technology | |
CN111523587A (en) | Woody plant species spectrum identification method based on machine learning | |
Borcea et al. | Reduced order models for spectral domain inversion: embedding into the continuous problem and generation of internal data | |
Kireev et al. | Improving the accuracy and sensitivity of 13C online detection in expiratory air using the TDLAS method in the spectral range of 4860–4880 cm− 1 | |
CN103940750B (en) | Saline soil curve of spectrum visible light wave range remote-sensing inversion construction process | |
JP2021515247A (en) | Digitization technology for cultural objects, calligraphy and original products | |
CN108827462A (en) | Irradiation level probe and optical measuring system | |
Yang et al. | Application and evaluation of wavelet-based denoising method in hyperspectral imagery data | |
CN105745550B (en) | It simulates to transcriber | |
CN106323466A (en) | Leaf nitrogen content high spectral evaluation method for continuous wavelet transformation analysis | |
CN104880476A (en) | Digital multi-channel system supporting elemental qualitative and quantitative analysis | |
CN105785398A (en) | Virtually imaged phased array-based incoherent atmospheric wind lidar | |
CN105320833A (en) | Method for determining optimum spectral resolution established by remote sensing inversion of spectrum character curve |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160323 Termination date: 20190421 |