CN109238991A - A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion - Google Patents
A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion Download PDFInfo
- Publication number
- CN109238991A CN109238991A CN201811353279.6A CN201811353279A CN109238991A CN 109238991 A CN109238991 A CN 109238991A CN 201811353279 A CN201811353279 A CN 201811353279A CN 109238991 A CN109238991 A CN 109238991A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- spectral
- imaging spectrometer
- visual field
- matching
- 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)
- Spectrometry And Color Measurement (AREA)
Abstract
The invention discloses a kind of big view field imaging spectrometer spectrum curvature correction methods of EO-1 hyperion, imaging spectrometer passes through the diffusion light of the sun spectrum that diffusing reflection plate obtains full filed first, then select a telluric lines as matching spectral line in central vision spectrum, and Spectral matching is carried out to obtain the pixel position of the telluric lines at different visual fields in imaging spectrometer direction in space, and then obtain spatial field of view spectrum bending value.The curved correction of imaging spectrometer spectrum is finally realized based on spectrum bending value.The present invention realizes the spectrum curvature correction based on telluric lines, eliminates the spectrum flexural distortion due to caused by big visual field, and method is reliable and simple and easy, is conducive to the spectral quality for improving imaging spectrometer.
Description
Technical field
The invention belongs to the big view field imaging spectrometer spectrum pictures of EO-1 hyperion to correct field, and in particular to a kind of EO-1 hyperion is big
View field imaging spectrometer spatial field of view spectrum curvature correction method.
Background technique
The spaceborne big view field imaging spectrometer of EO-1 hyperion is that pushing away on a kind of nadir direction sweeps type imaging spectrometer, using convex surface
Grating beam splitting, investigative range 310-406nm, spectral resolution < 1nm pass through detection earth atmosphere or earth surface reflection, scattering light
The distribution and variation for carrying out polluted gas in inverting atmosphere have the characteristics that high-resolution, big field angle.It is formed wearing rail direction
It is capable of forming 114 ° big visual field (corresponding ground 2600km when in orbit), detector uses area array CCD, can record quilt simultaneously
Survey the spectral information and spatial information of object.Since to wear rail orientation detection visual field larger for imaging spectrometer, so that measured target list
When colour images in the focal plane of imaging spectrometer, it may appear that the non-linear phenomena in space generates spectrum bending, spectrum bending
Presence can reduce spectrum picture quality, influence the inversion accuracy of spectroscopic data.Quickly to obtain the big view field imaging of EO-1 hyperion
The spectrum of spectrometer direction in space is bent, the curved correction of high-precision completion spectrum, it is necessary to propose a kind of big view of EO-1 hyperion
Field imaging spectrometer direction in space spectrum curvature correction method.
Robert A.Neville has studied the spectrum curved detection of imaging spectrometer (referring to Robert A, et at present
al.Neville.Detection of spectral line curvature in imaging spectrometer
Data.SPIE, 5093 (2003): 144-154.), the imaging spectrometer spectral resolution being directed to is lower, studies in text
AVIRIS and SFSI load spectral resolution is 10nm, and lower resolution ratio can not differentiate sunny fraunhofer line, so only
Wider infrared waveband absorbing peak, such as O can be chosen2(1268nm)、H2O(1470nm)、CO2(2007nm), not to ultraviolet band
It is studied, is not also provided the curved bearing calibration of spectrum, Bo-Cai Gao is using Spectral Matching Technique to imaging spectrometer data
The research in terms of spectral calibration has been carried out (referring to Bo-Cai Gao, et al.Refinement of wavelength
calibration of hyperspectral imaging data using a spectrum-matching
Technique.Remote Sensing of Environment, 2004 (90): 424-433), the Absorption Line chosen is
O at 760nm2Absorption peak, standard spectrum and matching spectrum between matching result determined by standard deviation, this determination method
It is low compared to Pearson correlation coefficient method precision, ultraviolet band is not studied, does not also provide the curved correction side of spectrum
Method.Jiankang Zhou research spectrum bending measuring method be based on tungsten halogen lamp-holimium oxide diffusing reflection plate (referring to
Jiankang Zhou,et al.Smile effect detection for dispersive by hypersepctral
imager based on the doped reflectance panel.Proc.Of SPIE.2012(8557):85571T-
1--85571T-8), i.e., the Absorption Line that the smooth spectral line that tungsten halogen lamp issues generates after the reflection of holimium oxide diffusing reflection plate is for examining
The spectrum bending for surveying imaging spectrometer, does not provide the curved bearing calibration of spectrum yet.
Research curved to imaging spectrometer spatial field of view spectrum is concentrated mainly on spectrally resolved lower (10nm), nothing at present
Method tells fraunhofer line, and Absorption Line is chosen at wider absorption peak (such as H of infrared band2O、O3Deng), it studies also at present not
The curved correction of spectrum is studied.For EO-1 hyperion, (spectral resolution < 1nm, high spectral resolution can be differentiated the present invention
More fraunhofer line out) big view field imaging spectrometer, Spectral matching is carried out using fraunhofer line in ultraviolet band and is determined
Imaging spectrometer spatial field of view spectrum bending value, and have studied the curved bearing calibration of spatial field of view spectrum.
Summary of the invention
Technology of the invention solves the problems, such as: overcoming the deficiencies of the prior art and provide a kind of big view field imaging spectrum of EO-1 hyperion
Instrument spectrum curvature correction method obtains the spectrogram of high quality to correct the spectrum flexural distortion due to caused by the big visual field in space
Picture improves spectra inversion precision.
In order to achieve the above object, the technical scheme adopted by the invention is as follows: imaging spectrometer is obtained by diffusing reflection plate
Full filed diffusion light of the sun, intensity is not preferably and by other Absorption Lines influenced in chosen spectrum signal, by this Absorption Line in sky
Between visual field carry out Spectral matching, according to matching result it is optimal obtain Absorption Line pixel position, the spectrum of full filed is calculated
Bending value completes spectrum curvature correction according to correcting algorithm based on this value.
Concrete methods of realizing is as follows:
A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion, the diffusion light of the sun light based on full filed
Spectrum chooses telluric lines as matching spectral line, carries out Spectral matching in imaging spectrometer direction in space to obtain spectrum bending
Value realizes the curved correction of imaging spectrometer spectrum, eliminates the spectrum flexural distortion due to caused by big visual field.
If full filed diffusion light of the sun spectrogram picture is Si,j, wherein i be spectrum channel number, value range 1~P of i ∈ by
As spectrometer detection spectral region determines that j is spatial field of view, and value range 1~Q of j ∈ is true by imaging spectrometer detection viewing field
It is fixed;
The first step, full filed diffusion light of the sun spectrogram picture obtain;
Imaging spectrometer observes diffusion light of the sun, guarantees that diffusion light of the sun observes visual field full of imaging spectrometer, obtains complete
Visual field diffusion light of the sun spectrogram picture is Si,j;
Second step, the selection of telluric lines;
From full filed diffusion light of the sun spectrogram as Si,jThe middle spectroscopic data S for extracting central visioni,center, with central vision
Spectroscopic data Si,centerFor reference spectra, from the spectroscopic data S of central visioni,centerMiddle Selection Center visual field Atmospheric Absorption
Line SΔλ,center, wherein Δ λ is the corresponding spectrum channel range of Absorption Line, the telluric lines that other visual field l, 1~Q of l ∈ choose
For SΔλ,l;
Third step, the setting of Spectral matching step-length;
Based on imaging spectrometer spectral calibration precision, central vision telluric lines S is setΔλ,centerWith other visual field l, l
∈ 1~Q telluric lines SΔλ,lCarry out matched step-length δ;
4th step, Spectral matching;
By central vision telluric lines SΔλ,centerAccording to matching step-length δ and other visual field l, l ∈ 1~Q telluric lines
SΔλ,lMatched, the matching result under each step-length determined according to Pearson correlation coefficient method, obtain matching result it is optimal when
Corresponding step-length, is visual field l, the spectral shift value δ under 1~Q of l ∈l;
5th step, spatial field of view spectrum bending value determine
According to the 4th step Spectral matching method to full filed diffusion light of the sun spectrogram as Si,jLight is carried out in spatial field of view direction
Spectrum matching, obtains the spectral shift value δ under each visual fieldj, j 1~Q of ∈, based on spectral shift value δjIt is curved to obtain spatial field of view spectrum
Bent value Cj;
6th step, spectrum curvature correction
Calibration model is established, according to spatial field of view spectrum bending value CjTo full filed diffusion light of the sun spectrogram as Si,jIt carries out
Spectrum curvature correction is finally completed complete visual field spectrum curvature correction.
The first step, full filed diffusion light of the sun spectrogram is as Si,jIt obtains, is implemented as follows:
(1) select fair weather to guarantee that it is less that diffusion light of the sun is affected by atmospheric effects, direct sunlight is incident on bright
To form uniform area source on the aluminium diffusing reflection plate of Bert, aluminium diffusing reflection plate and imaging spectrometer position angle are adjusted, is protected
Imaging spectrometer spatial field of view can be full of by demonstrate,proving the area source;
(2) after imaging spectrometer spatial field of view is full of, the time of integration, gain parameter are set, guarantee that the sun observed dissipates
Light signal strength reaches saturation value 80% is penetrated, spectroscopic data is acquired, obtains full filed diffusion light of the sun spectrogram as Si,jIt is expressed as
Si,j=[Si,1,Si,2,…,Si,Q], wherein Si,1,Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, i is spectrum channel
Number 1~P of i ∈;
The second step, central vision telluric lines SΔλ,centerWith other visual field telluric lines SΔλ,lSelection, tool
Body is accomplished by
(1) from full filed diffusion light of the sun spectrogram as Si,jMiddle extraction central vision spectroscopic data Si,center;
(2) to center visual field spectrum data Si,centerIn the telluric lines that contains analyzed, if Absorption Line signal is strong
It spends lower, gives up, give up if being influenced by other Absorption Lines, it is higher and do not influenced by other Absorption Lines to choose signal strength
Absorption Line of the Absorption Line as Spectral matching;
(3) after determining matching Absorption Line, it is based on imaging spectrometer spectrum channel number, determines telluric lines spectral region
Δ λ obtains central vision telluric lines SΔλ,center;
(4) visual field l, l ∈ 1~Q telluric lines S are obtained according to step (2), (3)Δλ,l;
The third step, the setting of Spectral matching step-length δ, is implemented as follows:
It is set as spectrometer detection spectral region being λ, spectrum channel bandwidth isP is the light of imaging spectrometer detection
Port number is composed, if the spectral calibration precision of imaging spectrometer is σ, then Spectral matching step-length
4th step, Spectral matching are implemented as follows:
(1) central vision telluric lines SΔλ,centerAs reference spectral line and visual field l, l ∈ (1-Q) telluric lines
SΔλ,lSpectral matching is carried out according to matching step-length δ;
(2) result of Spectral matching is determined using Pearson correlation coefficient method, central vision telluric lines and view
Telluric lines S under field l, 1~Q of l ∈Δλ,jRelated coefficientWherein α=∑ SΔλ,center·
SΔλ,l, β=∑ SΔλ,center·∑SΔλ,l, κ=∑ (SΔλ,center)2-(∑SΔλ,center)2, η=∑ (SΔλ,l)2-(∑SΔλ,l)2,N
For sampling number, rPearsonThe range of absolute value is 0~1, and value is bigger, and expression correlation is stronger, and spectral line matching result is more excellent;
(3) according to the matching result under each step-length, obtain matching result it is optimal when corresponding matching step-length δlFor visual field
Spectral shift value under 1~Q of l, l ∈.
5th step, spatial field of view spectrum bending value determines, is implemented as follows:
Spectral matching is iterated processing in imaging spectrometer spatial field of view, obtains the spectral shift value under each visual field
δj, and then determine spatial field of view spectrum bending value Cj, it is written as Cj=[δ1,δ2,…,δQ], wherein δ1,δ2,…,δQFor 1~Q of j ∈
Under spectral shift value, be visual field function.
6th step, the 6th step, spectrum curvature correction are implemented as follows:
(1) the spectral shift value δ under 1~Q of visual field l, l ∈lBalance correction, translation direction are carried out to the spectrum under this visual field
It is determined by spectrum bending direction, for detecting band edge, translates corresponding spectral response value and be set as 0;
(2) according to the spectral correction method in (1), it is based on spectrum amount of bow Cj=[δ1,δ2,…,δQ], wherein δ1,
δ2,…,δQFor the spectral shift value under 1~Q of j ∈ to full filed diffusion light of the sun spectrogram as Si,j=[Si,1,Si,2,…,Si,Q] into
Row corrects, wherein Si,1,Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, 1~P of i ∈ is spectrum channel number, final complete
At the curved correction of spatial field of view spectrum.
The advantages of the present invention over the prior art are that:
(1) present invention can complete the spectrum bending inspection of high spectral resolution (spectral resolution < 1nm) imaging spectrometer
It surveys, chooses H with current research2O、O3It is compared etc. wider Absorption Line (line width reaches 15nm), from center visual field spectrum data
Si,centerThe Absorption Line extracted in spectrum is that the fraunhofer line (line width < 1nm) of fine structure carries out Spectral matching, be can be improved
The computational accuracy of Pearson correlation coefficient.
(2) compared with existing Spectral matching, fraunhofer line and higher spectral calibration precision using fine structure,
Smaller matching step-length δ can be set, improve Spectral matching precision.
(3) existing spectrum is bent the technology of determination and does not study spectrum curvature correction, and the present invention is obtaining spectrum
On the basis of curved, spectrum curvature correction is increased, high-precision spectrum curvature correction, energy are completed by balance correction method
It enough further increases spectrum picture quality, improve data inversion precision.
Detailed description of the invention
Fig. 1 is big view field imaging spectrometer spatial spectral curvature correction method flow;
Fig. 2 is the full filed diffusion light of the sun spectrogram applied to this example, and spectrogram horizontal direction is spectrum dimension, spectrum dimension
Corresponding incident light dispersion direction, vertical direction is space dimension, and space dimension corresponds to the observation visual field direction of imaging spectrometer;
The fraunhofer Absorption Line for Spectral matching that Fig. 3 chooses from Fig. 2;
Fig. 4 is the central vision and peripheral field spectrum comparing result applied to the embodiment of the present invention, since spectrum is bent
Presence to deviate between two spectral lines;
Fig. 5 be applied to the spectrum comparing result after the spectrum curvature correction of the embodiment of the present invention, after spectrum curvature correction,
So that the Absorption Line of two spectrum is overlapped.
Specific embodiment
With reference to the accompanying drawing and specific embodiment further illustrates the present invention.
As shown in fig. 1~fig. 5, in the present invention, spectrum curvature correction method is to be directed to the big view field imaging spectrometer of EO-1 hyperion,
It requires the spectral resolution of imaging spectrometer higher, the telluric lines in diffusion light of the sun can be told, based on full view
Field diffusion light of the sun spectrogram picture realizes the spectrum curvature correction on direction in space.
If full filed diffusion light of the sun spectrogram picture is Si,j, wherein i be spectrum channel number, value range 1~P of i ∈ by
As spectrometer detection spectral region determines that j is spatial field of view, and value range 1~Q of j ∈ is true by imaging spectrometer detection viewing field
It is fixed.
1, full filed diffusion light of the sun spectrogram is as Si,jIt obtains, is embodied as follows:
(1) select fair weather to guarantee that it is less that diffusion light of the sun is affected by atmospheric effects, direct sunlight is incident on bright
To form uniform area source on the aluminium diffusing reflection plate of Bert, aluminium diffusing reflection plate and imaging spectrometer position angle are adjusted, is protected
Card area source can be full of imaging spectrometer spatial field of view.
(2) after imaging spectrometer spatial field of view is full of, the time of integration, gain parameter are set, guarantee that the sun observed dissipates
Light signal strength reaches saturation value 80% is penetrated, spectroscopic data is acquired, obtains full filed diffusion light of the sun spectrogram as Si,j, such as Fig. 2
The full filed diffusion light of the sun spectrogram of imaging spectrometer acquisition is shown as Si,j=[Si,1,Si,2,…,Si,Q], wherein Si,1,
Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, i is 1~P of spectrum channel number i ∈.Fig. 2 horizontal direction is spectrum dimension,
Spectrum ties up corresponding incident light dispersion direction, and vertical direction is space dimension, and space dimension corresponds to 114 ° of spatial field of view of imaging spectrometer.
Black line in spectrum picture is telluric lines, and the curved presence of spectrum is so that black line bends in vertical direction.
2, central vision telluric lines SΔλ,centerWith other visual field telluric lines SΔλ,lSelection, specific implementation such as
Under:
Because 344.25nm Absorption Line signal strength is higher and is not influenced by other spectral lines, as shown in figure 3, selecting this Absorption Line
As Spectral matching spectral line SΔλ,center, the range delta lambda of selection is 661-681 pixel, and the pixel range of selection can be retouched completely
State the line style of entire Absorption Line.
3, the setting of Spectral matching step-length δ is embodied as follows:
(1) imaging spectrometer detecting light spectrum range lambda=96nm, spectrum channel number P=1024, spectral bandwidth are
Spectral calibration precision is 0.05nm, the matching step-length of settingI.e. matching step-length is 0.5 pixel.
(2) central vision telluric lines SΔλ,centerAs reference spectral line and visual field l, l ∈ (1-Q) telluric lines
SΔλ,lSpectral matching is carried out according to matching step-length δ;
(3) according to the matching result under each step-length, obtain matching result it is optimal when corresponding matching step-length δlFor visual field
Spectral shift value under 1~Q of l, l ∈;
4, the acquisition of spectrum bending value is embodied as follows:
It is iterated processing in 114 ° of spatial field of view of imaging spectrometer according to step 3, the spectrum obtained under each visual field is inclined
Shifting value δj, and then determine spatial field of view spectrum bending value Cj, it is written as Cj=[δ1,δ2,…,δQ], wherein δ1,δ2,…,δQFor j ∈
Spectral shift value under 1~Q is the function of visual field.As Fig. 4 shows the spectrum of imaging spectrometer central vision and peripheral field
Bending value is 8 pixels.
5, spectrum curvature correction is embodied as follows:
(1) the spectral shift value δ according to visual field l, under 1~Q of l ∈lBalance correction, translation are carried out to the spectrum under this visual field
Direction is determined by spectrum bending direction, for detecting band edge, is translated corresponding spectral response value and is set as 0;
(2) according to the spectral correction method in (1), it is based on spectrum amount of bow Cj=[δ1,δ2,…,δQ], wherein δ1,
δ2,…,δQFor the spectral shift value under 1~Q of j ∈ to 114 ° of full filed diffusion light of the sun spectrograms of imaging spectrometer as Si,j=
[Si,1,Si,2,…,Si,Q] be corrected, wherein Si,1,Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, 1~P of i ∈
For spectrum channel number, it is finally completed the curved correction of spatial field of view spectrum.Fig. 5 is the diffusion light of the sun after spectrum curvature correction
It composes, after spectrum curvature correction, so that the Absorption Line of two spectrum is overlapped.
Claims (7)
1. a kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion, it is characterised in that the following steps are included:
If full filed diffusion light of the sun spectrogram picture is Si,j, wherein i is spectrum channel number, and 1~P of value range i ∈ is by imaging
Spectrometer detecting light spectrum range determines that j is spatial field of view, and value range 1~Q of j ∈ is determined by imaging spectrometer detection viewing field.
The first step, full filed diffusion light of the sun spectrogram picture obtain;
Imaging spectrometer observes diffusion light of the sun, guarantees that diffusion light of the sun observes visual field full of imaging spectrometer, obtains full filed
Diffusion light of the sun spectrogram picture is Si,j,
Second step, the selection of telluric lines;
From full filed diffusion light of the sun spectrogram as Si,jThe middle spectroscopic data S for extracting central visioni,center, with the light of central vision
Modal data Si,centerFor reference spectra, from the spectroscopic data S of central visioni,centerMiddle Selection Center visual field telluric lines
SΔλ,center, wherein Δ λ is the corresponding spectrum channel range of Absorption Line, and the telluric lines that other visual field l, 1~Q of l ∈ choose is
SΔλ,l;
Third step, the setting of Spectral matching step-length;
Central vision telluric lines S is arranged in spectral calibration precision based on imaging spectrometerΔλ,centerWith other visual field l, l ∈
1~Q telluric lines SΔλ,lCarry out matched step-length δ;
4th step, Spectral matching;
By central vision telluric lines SΔλ,centerAccording to matching step-length δ and other visual field l, l ∈ 1~Q telluric lines SΔλ,l
It is matched, the matching result under each step-length is determined according to Pearson correlation coefficient method, obtain corresponding to when matching result is optimal
Step-length, be visual field l, the spectral shift value δ under 1~Q of l ∈l;
5th step, spatial field of view spectrum bending value determine;
According to the 4th step Spectral matching method to full filed diffusion light of the sun spectrogram as Si,jSpectrum is carried out in spatial field of view direction
Match, obtains the spectral shift value δ under each visual fieldj, j 1~Q of ∈, based on spectral shift value δjObtain spatial field of view spectrum bending value
Cj;
6th step, spectrum curvature correction
Calibration model is established, according to spatial field of view spectrum bending value CjTo full filed diffusion light of the sun spectrogram as Si,jIt is curved to carry out spectrum
Qu Jiaozheng is finally completed complete visual field spectrum curvature correction.
2. according to the method described in claim 1, it is characterized by: the first step, full filed diffusion light of the sun spectrogram is as Si,j
Acquisition is implemented as follows:
(11) select fair weather to guarantee that it is less that diffusion light of the sun is affected by atmospheric effects, direct sunlight is incident on Lambertian
Property aluminium diffusing reflection plate on to form uniform area source, adjust aluminium diffusing reflection plate and imaging spectrometer position angle, guarantee institute
Imaging spectrometer spatial field of view can be full of by stating area source;
(12) after imaging spectrometer spatial field of view is full of, the time of integration, gain parameter is set, guarantee the diffusion light of the sun observed
Signal strength reaches the 80% of saturation value, acquires spectroscopic data, obtains full filed diffusion light of the sun spectrogram as Si,jIt is expressed as Si,j
=[Si,1,Si,2,…,Si,Q], wherein Si,1,Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, i is spectrum channel number i
1~P of ∈.
3. according to the method described in claim 1, it is characterized by: second step, central vision telluric lines SΔλ,centerWith it
Its visual field telluric lines SΔλ,lSelection, be implemented as follows:
(21) from full filed diffusion light of the sun spectrogram as Si,jMiddle extraction central vision spectroscopic data Si,center;
(22) to center visual field spectrum data Si,centerIn the telluric lines that contains analyzed, if Absorption Line signal strength compared with
It is low, give up, give up if being influenced by other Absorption Lines, chooses that signal strength is higher and suction that do not influenced by other Absorption Lines
Absorption Line of the take-up as Spectral matching;
(23) after determining matching Absorption Line, it is based on imaging spectrometer spectrum channel number, that is, spectrum sample interval, determines big aspiration
Take-up spectral region Δ λ, obtains central vision telluric lines SΔλ,center;
(24) visual field l, l ∈ 1~Q telluric lines S are obtained according to step (22), (23)Δλ,l。
4. according to the method described in claim 1, it is characterized by: third step, the setting of Spectral matching step-length δ, specific implementation
It is as follows:
It is set as spectrometer detection spectral region being λ, spectrum channel bandwidth isP is that the spectrum of imaging spectrometer detection is logical
Road number, if the spectral calibration precision of imaging spectrometer is σ, then Spectral matching step-length
5. according to the method described in claim 1, it is characterized by: the 4th step, Spectral matching are implemented as follows:
(41) central vision telluric lines SΔλ,centerAs reference spectral line and visual field l, l ∈ (1-Q) telluric lines SΔλ,lIt presses
Spectral matching is carried out according to matching step-length δ;
(42) result of Spectral matching is determined using Pearson correlation coefficient method, central vision telluric lines SΔλ,center
With visual field l, telluric lines S under 1~Q of l ∈Δλ,jRelated coefficientWherein α=∑ SΔλ,center·SΔλ,l,
β=∑ SΔλ,center·∑SΔλ,l, κ=∑ (SΔλ,center)2-(∑SΔλ,center)2, η=∑ (SΔλ,l)2-(∑SΔλ,l)2, N is
Sampling number, rPearsonThe range of absolute value is 0~1, and value is bigger, and expression correlation is stronger, and spectral line matching result is more excellent;
(43) according to the matching result under each step-length, obtain matching result it is optimal when corresponding matching step-length δlFor visual field l, l ∈
Spectral shift value under 1~Q.
6. according to the method described in claim 1, it is characterized by: the 5th step, spatial field of view spectrum bending value determine, specifically in fact
Now as follows: Spectral matching is iterated processing in imaging spectrometer spatial field of view, obtains the spectral shift value δ under each visual fieldj,
And then determine spatial field of view spectrum bending value Cj, it is written as Cj=[δ1,δ2,…,δQ], wherein δ1,δ2,…,δQFor under 1~Q of j ∈
Spectral shift value, be visual field function.
7. according to the method described in claim 1, it is characterized by: the 6th step, spectrum curvature correction, specific implementation is such as
Under:
(61) the spectral shift value δ under 1~Q of visual field l, l ∈lBalance correction is carried out to the spectrum under this visual field, translation direction is by light
It composes bending direction to determine, for detecting band edge, translates corresponding spectral response value and be set as 0;
(62) according to the spectral correction method in (61), it is based on spectrum amount of bow Cj=[δ1,δ2,…,δQ], wherein δ1,δ2,…,δQ
For the spectral shift value under 1~Q of j ∈ to full filed diffusion light of the sun spectrogram as Si,j=[Si,1,Si,2,…,Si,Q] be corrected,
Wherein Si,1,Si,2,…,Si,QFor the spectrum under 1~Q of spatial field of view j ∈, 1~P of i ∈ is spectrum channel number, is finally completed space
The curved correction of visual field spectrum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811353279.6A CN109238991A (en) | 2018-11-14 | 2018-11-14 | A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811353279.6A CN109238991A (en) | 2018-11-14 | 2018-11-14 | A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109238991A true CN109238991A (en) | 2019-01-18 |
Family
ID=65074650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811353279.6A Pending CN109238991A (en) | 2018-11-14 | 2018-11-14 | A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109238991A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220557A (en) * | 2019-11-21 | 2020-06-02 | 中国科学院合肥物质科学研究院 | Full-field solar reference spectrum acquisition method of hyperspectral large-field imaging spectrometer |
CN112284537A (en) * | 2020-09-15 | 2021-01-29 | 中国科学院上海技术物理研究所 | Push-broom hyperspectral imager spectrum bending and keystone distortion detection method |
CN116593001A (en) * | 2023-07-18 | 2023-08-15 | 长春理工大学 | Imaging spectrometer optical system spectral line bending correction method based on macro file |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120098924A1 (en) * | 2008-11-26 | 2012-04-26 | Honeywell International Inc. | Signal spectra detection system |
CN104316187A (en) * | 2014-11-08 | 2015-01-28 | 中国科学院合肥物质科学研究院 | Spectrum calibration device for large-visual-field hyper-spectrum imaging differential absorption spectrograph |
CN105259127A (en) * | 2015-11-27 | 2016-01-20 | 成都理工大学 | Vegetation index calculation method based on comprehensive absorption capacity |
CN105424186A (en) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | Spectrum calibration and correction method of light field imaging spectrometer |
CN107407629A (en) * | 2015-03-05 | 2017-11-28 | 生物辐射实验室股份有限公司 | The Spectral matching of optimization and display |
-
2018
- 2018-11-14 CN CN201811353279.6A patent/CN109238991A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120098924A1 (en) * | 2008-11-26 | 2012-04-26 | Honeywell International Inc. | Signal spectra detection system |
CN104316187A (en) * | 2014-11-08 | 2015-01-28 | 中国科学院合肥物质科学研究院 | Spectrum calibration device for large-visual-field hyper-spectrum imaging differential absorption spectrograph |
CN107407629A (en) * | 2015-03-05 | 2017-11-28 | 生物辐射实验室股份有限公司 | The Spectral matching of optimization and display |
CN105424186A (en) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | Spectrum calibration and correction method of light field imaging spectrometer |
CN105259127A (en) * | 2015-11-27 | 2016-01-20 | 成都理工大学 | Vegetation index calculation method based on comprehensive absorption capacity |
Non-Patent Citations (4)
Title |
---|
ZHAO MIN-JIE: "Study on the Non-Uniformity Calibration of Space-Born Differential Optical Absorption Spectrometer", 《光谱学与光谱分析》 * |
赵敏杰: "基于谱线匹配技术的星载成像光谱仪星上光谱定标方法研究", 《光谱学与光谱分析》 * |
赵敏杰: "星载大气痕量气体差分吸收光谱仪光谱图像校正方法研究", 《光学学报》 * |
赵敏杰: "星载大气痕量气体差分吸收光谱仪定标技术研究", 《学位论文知识发现系统》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220557A (en) * | 2019-11-21 | 2020-06-02 | 中国科学院合肥物质科学研究院 | Full-field solar reference spectrum acquisition method of hyperspectral large-field imaging spectrometer |
CN111220557B (en) * | 2019-11-21 | 2022-08-23 | 中国科学院合肥物质科学研究院 | Full-field solar reference spectrum acquisition method of hyperspectral large-field imaging spectrometer |
CN112284537A (en) * | 2020-09-15 | 2021-01-29 | 中国科学院上海技术物理研究所 | Push-broom hyperspectral imager spectrum bending and keystone distortion detection method |
CN112284537B (en) * | 2020-09-15 | 2022-07-08 | 中国科学院上海技术物理研究所 | Push-broom hyperspectral imager spectrum bending and keystone distortion detection method |
CN116593001A (en) * | 2023-07-18 | 2023-08-15 | 长春理工大学 | Imaging spectrometer optical system spectral line bending correction method based on macro file |
CN116593001B (en) * | 2023-07-18 | 2024-03-12 | 长春理工大学 | Imaging spectrometer optical system spectral line bending correction method based on macro file |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Suzuki et al. | Initial inflight calibration for Hayabusa2 optical navigation camera (ONC) for science observations of asteroid Ryugu | |
CN109238991A (en) | A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion | |
Marleau et al. | Extragalactic source counts at 24 microns in the Spitzer First Look survey | |
Guo et al. | Atmospheric correction comparison of SPOT-5 image based on model FLAASH and model QUAC | |
Montanes-Rodriguez et al. | Globally integrated measurements of the Earth’s visible spectral albedo | |
CN111145351B (en) | Minnarert terrain correction model optimization method considering ground feature types | |
CN109974854A (en) | A kind of radiation correction method of frame width formula FPI high spectrum image | |
CN105548032A (en) | Compact high-resolution wide-view-field spectral imaging system | |
CN109300133A (en) | A kind of city river network Clean water withdraw method | |
CN113790798B (en) | Seamless spectral imaging device, system and method for dynamic point target tracking measurement | |
Schläpfer et al. | Correction of shadowing in imaging spectroscopy data by quantification of the proportion of diffuse illumination | |
Delaney et al. | Visible and infrared reflectance imaging spectroscopy of paintings: pigment mapping and improved infrared reflectography | |
CN110702228B (en) | Edge radiation correction method for aviation hyperspectral image | |
CN104568148A (en) | Optical system for atmosphere CO2 hyper-spectrum imaging spectrometer | |
Kneubühler et al. | Geometric and radiometric pre-processing of CHRIS/PROBA data over mountainous terrain | |
Davis et al. | Calibration, characterization, and first results with the Ocean PHILLS hyperspectral imager | |
CN109781260A (en) | The fast illuminated polarization spectrum imaging detection device of ultra-compact and detection method | |
CN111220557B (en) | Full-field solar reference spectrum acquisition method of hyperspectral large-field imaging spectrometer | |
Fang et al. | Comparision of eight topographic correction algorithms applied to landsat-8 oli imagery based on the dem | |
Green et al. | Determination of in-flight AVIRIS spectral, radiometric, spatial and signal-to-noise characteristics using atmospheric and surface measurements from the vicinity of the rare-earth-bearing carbonatite at Mountain Pass, California | |
CN111735538B (en) | Airborne area array staring type hyperspectral image illumination correction method | |
Wiehr | Magnetic field strength and inclination in the penumbral fine-structure | |
Zhao et al. | Improvement of dark object method in atmospheric correction of hyperspectral remotely sensed data | |
Dubs et al. | Calibration of meteor spectra | |
Hiriart et al. | Molecular hydrogen kinematics in Cepheus A |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190118 |