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 PDF

Info

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
Application number
CN201811353279.6A
Other languages
Chinese (zh)
Inventor
赵敏杰
司福祺
江宇
汪世美
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201811353279.6A priority Critical patent/CN109238991A/en
Publication of CN109238991A publication Critical patent/CN109238991A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating 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

A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion
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=[δ12,…,δQ], wherein δ12,…,δ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=[δ12,…,δ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=[δ12,…,δQ], wherein δ12,…,δ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=[δ12,…,δ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=[δ12,…,δQ], wherein δ12,…,δ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=[δ12,…,δQ], wherein δ12,…,δ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.
CN201811353279.6A 2018-11-14 2018-11-14 A kind of big view field imaging spectrometer spectrum curvature correction method of EO-1 hyperion Pending CN109238991A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
ZHAO MIN-JIE: "Study on the Non-Uniformity Calibration of Space-Born Differential Optical Absorption Spectrometer", 《光谱学与光谱分析》 *
赵敏杰: "基于谱线匹配技术的星载成像光谱仪星上光谱定标方法研究", 《光谱学与光谱分析》 *
赵敏杰: "星载大气痕量气体差分吸收光谱仪光谱图像校正方法研究", 《光学学报》 *
赵敏杰: "星载大气痕量气体差分吸收光谱仪定标技术研究", 《学位论文知识发现系统》 *

Cited By (6)

* Cited by examiner, † Cited by third party
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