CN108254074A - A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument - Google Patents
A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument Download PDFInfo
- Publication number
- CN108254074A CN108254074A CN201711284886.7A CN201711284886A CN108254074A CN 108254074 A CN108254074 A CN 108254074A CN 201711284886 A CN201711284886 A CN 201711284886A CN 108254074 A CN108254074 A CN 108254074A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- orbit
- sun
- standard
- sensing instrument
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 215
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000010521 absorption reaction Methods 0.000 claims abstract description 67
- 230000003595 spectral effect Effects 0.000 claims abstract description 51
- 238000006073 displacement reaction Methods 0.000 claims abstract description 4
- 239000000284 extract Substances 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 238000004164 analytical calibration Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J2003/2866—Markers; Calibrating of scan
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The present invention provides a kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument, including:S1, in-orbit sun observation data and ground calibration coefficient based on target high-spectrum remote-sensing instrument calculate each channel center's wavelength location of in-orbit sun observation spectrum of the target high-spectrum remote-sensing instrument;S2 is accurately positioned velocity information based on satellite, according to Doppler shift characteristics, carries out frequency displacement amendment to in-orbit each channel center's wavelength location of sun observation spectrum, obtains each channel active centre wavelength location of in-orbit sun observation spectrum;S3 based on each channel active centre wavelength location of the in-orbit sun observation spectrum and sun standard spectrum standard absorption line position, calculates the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument;S4 based on the in-orbit spectral calibration error, carries out the in-orbit spectrum precise calibration of the target high-spectrum remote-sensing instrument.The present invention can effectively improve the accuracy and stability of spectral calibration.
Description
Technical field
The present invention relates to remote sensing calibration technology fields, determine more particularly, to a kind of in-orbit spectrum of high-spectrum remote-sensing instrument
Mark method.
Background technology
It is distant to satellite with the continuous improvement of the continuous improvement of application demand, especially Atmospheric components quantitative detection demand
The requirement of the spectral resolution of sense is also higher and higher.Therefore, one of direction of Modern Satellite remote sensing instrument development is with superelevation
Spectral resolution.
The carbon satellite of China's successful launch, by weak CO2 absorption bands (1.6 μm), strong CO2 absorption bands (2.06 μm) and
The observation spectrums of O2-A absorption bands (0.76 μm) carries out high-precision CO2 concentration quantitatives inverting, highest resolution up to 0.03nm,
In the range of about 15nm, there are more than 1000 a channels.
EO-1 hyperion instrument in orbit after, due to being influenced by variation of ambient temperature and the effect of various minimal stress,
Small drift can occur for the centre wavelength position of channel.This drift influences less for fat pipe instrument, but for
EO-1 hyperion instrument can but be generated and be seriously affected.Channel spectrum center caused by small drift changes, can be to quantitative inversion
Product Precision generation seriously affect, therefore, the in-orbit spectral calibration of instrument become high-spectrum remote sensing data pretreatment work weight
Technology point is wanted, is the important prerequisite of accurate radiation calibration.
Invention content
In order to overcome the above problem or solve the above problems at least partly, the present invention provides a kind of high-spectrum remote-sensing instrument
The in-orbit spectrum calibration method of device, to effectively improve the accuracy of spectral calibration and stability.
The present invention provides a kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument, including:S1 is distant based on target EO-1 hyperion
Feel the in-orbit sun observation data of instrument and ground calibration coefficient, the in-orbit sun for calculating the target high-spectrum remote-sensing instrument is seen
Survey each channel center's wavelength location of spectrum;S2 is accurately positioned velocity information based on satellite, according to Doppler shift characteristics, to institute
It states each channel center's wavelength location of in-orbit sun observation spectrum and carries out frequency displacement amendment, obtaining in-orbit each channel of sun observation spectrum has
Imitate centre wavelength position;S3, based on each channel active centre wavelength location of the in-orbit sun observation spectrum and sun standard light
Standard absorption line position is composed, calculates the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument;S4, based on described in-orbit
Spectral calibration error carries out the in-orbit spectrum precise calibration of the target high-spectrum remote-sensing instrument.
Wherein, the selection principle of sun standard spectrum includes in the step of S3:The resolution of the sun standard spectrum
Rate is higher than the resolution ratio of the target high-spectrum remote-sensing instrument.
Wherein, the sun standard spectrum is further specifically chosen for Kurucz R.L. data.
Wherein, the selecting step of standard absorption line described in step S3 further comprises:It chooses absorption depth and reaches setting
The width of rank and both wings reaches the Absorption Line of predetermined width as the standard absorption line.
Wherein, the standard absorption line is specially further that the bright standing grain of husband takes Absorption Line.
Wherein, the step of S1 is further specifically included:In extracting the target high-spectrum remote-sensing instrument in orbit
The in-orbit sun observation data obtained by diffusing reflection plate;It is tried by the ground calibration to the target high-spectrum remote-sensing instrument
It tests, calculates the ground calibration coefficient;Based on the in-orbit sun observation data and the ground calibration coefficient, given plan is utilized
Formula is closed, calculates in-orbit each channel center's wavelength location of sun observation spectrum.
Wherein, the step of S2 further comprises:Velocity information is accurately positioned based on the satellite, utilizes Doppler's frequency
Characteristic formula is moved, calculates the in-orbit doppler drift of the target high-spectrum remote-sensing instrument;The in-orbit sun is eliminated to see
The in-orbit doppler drift component in each channel center's wavelength location of spectrum is surveyed, it is each to obtain the in-orbit sun observation spectrum
Channel active centre wavelength location.
Wherein, the step of S3 is further specifically included:Extract in-orbit each channel of sun observation spectrum it is effective in
Active centre wavelength location corresponding with the sun standard spectrum standard absorption line position in heart wavelength location, by the correspondence
Active centre wavelength location matched with the sun standard spectrum standard absorption line position;In the middle cardiac wave respectively to match
At long position, each relatively described sun standard spectrum standard absorption line position of the corresponding active centre wavelength location is calculated
Wave length shift.
Wherein, the corresponding active centre observing frequency of the active centre wavelength location is calculated by the following formula:
In formula, vobsRepresent effective observation spectrum frequency, vsunRepresent sun standard spectrum frequency, VrelIt represents along EO-1 hyperion
Remote sensing instrument direction of visual lines satellite is with respect to the speed of the sun, c expression vacuum light speeds.
Wherein, in-orbit sun observation data described in step S1 are specially further that sun observation continues the flat of each frame data
Mean value.
A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument provided by the invention is made by using sun standard spectrum
For standard, deviation between the in-orbit observation spectrum of target high-spectrum remote-sensing instrument and standard spectrum is calculated, so as to fulfill target EO-1 hyperion
The in-orbit spectral calibration of remote sensing instrument can make full use of the stability feature of solar spectrum, determine so as to effectively improve spectrum
Target accuracy and stability.
Description of the drawings
Fig. 1 is a kind of flow chart of the in-orbit spectrum calibration method of high-spectrum remote-sensing instrument of the embodiment of the present invention;
Fig. 2 is a kind of calculating process flow of in-orbit each channel center's wavelength location of sun observation spectrum of the embodiment of the present invention
Figure;
Fig. 3 is a kind of acquisition process of in-orbit each channel active centre wavelength location of sun observation spectrum of the embodiment of the present invention
Flow chart;
Fig. 4 is a kind of calculating process flow of the in-orbit spectral calibration error of high-spectrum remote-sensing instrument of the embodiment of the present invention
Figure;
Fig. 5 is the flow chart of the in-orbit spectrum calibration method of another kind high-spectrum remote-sensing instrument of the embodiment of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached in the embodiment of the present invention
Figure, is clearly and completely described the technical solution in the present invention, it is clear that described embodiment is one of the present invention
Divide embodiment, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not making
All other embodiments obtained under the premise of creative work, shall fall within the protection scope of the present invention.
As one embodiment of the embodiment of the present invention, the present embodiment provides a kind of in-orbit spectrum of high-spectrum remote-sensing instrument to determine
Mark method is a kind of flow chart of the in-orbit spectrum calibration method of high-spectrum remote-sensing instrument of the embodiment of the present invention with reference to figure 1, including:
S1, in-orbit sun observation data and ground calibration coefficient based on target high-spectrum remote-sensing instrument, calculates the mesh
Each channel center's wavelength location of in-orbit sun observation spectrum of absolute altitude spectral remote sensing instrument;
S2 is accurately positioned velocity information based on satellite, according to Doppler shift characteristics, to the in-orbit sun observation spectrum
Each channel center's wavelength location carries out frequency displacement amendment, obtains each channel active centre wavelength location of in-orbit sun observation spectrum;
S3 is inhaled based on each channel active centre wavelength location of the in-orbit sun observation spectrum and sun standard spectrum standard
Take-up position calculates the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument;
S4, based on the in-orbit spectral calibration error, the in-orbit spectrum for carrying out the target high-spectrum remote-sensing instrument is accurate
Calibration.
Step S1 carries out target high-spectrum remote-sensing instrument it is to be understood that before satellite launch in ground experiment room
Analogue scaling is tested, and obtains ground calibration coefficient.After satellite entrance is normal in orbit, target high-spectrum remote-sensing is utilized
Instrument measures sufficient amount of practical sun observation data, i.e., in-orbit sun observation data.
On the basis of above-mentioned ground calibration coefficient and in-orbit sun observation data has been obtained, according to ground calibration coefficient
With in-orbit sun observation data, each channel center wavelength in the in-orbit sun observation spectrum of target high-spectrum remote-sensing instrument is calculated
Position νorbit。
Wherein optional, the S1's is further processed step with reference to figure 2, is that a kind of in-orbit sun of the embodiment of the present invention is seen
The calculating process flow chart of each channel center's wavelength location of spectrum is surveyed, including:
S11 is seen in extracting the target high-spectrum remote-sensing instrument in orbit by the in-orbit sun that diffusing reflection plate obtains
Measured data;
S12 is tested by the ground calibration to the target high-spectrum remote-sensing instrument, calculates the ground calibration coefficient;
S13 based on the in-orbit sun observation data and the ground calibration coefficient, using given fitting formula, is calculated
Each channel center's wavelength location of in-orbit sun observation spectrum.
It is to be understood that the present embodiment is directed to visible ray and near infrared light spectral coverage, satellite completes Yang Zhao area's earth observations
Afterwards, target high-spectrum remote-sensing instrument can carry out sun observation by self-contained diffusing reflection plate, obtain in-orbit sun observation number
According to.Since the sun is highly stable reliable target source, the present embodiment is using sun standard spectrum as spectral calibration source.
Under the precondition for having in-orbit sun observation in satellite, obtained by target high-spectrum remote-sensing instrument by diffusing reflection plate
The in-orbit sun observation data taken, and the in-orbit sun observation number is extracted by target high-spectrum remote-sensing instrument when being calculated
According to.
On the basis of in-orbit sun observation data are obtained, also need to carry out the bracketing of ground experiment room, to establish
Rail observes the incidence relation between data and terrestrial reference.Specifically before satellite launch, satellite is carried out by laboratory normal
The ground rating test of rule calculates the ground calibration coefficient for target high-spectrum remote-sensing instrument.
Finally, on the basis of above-mentioned steps, each channel center's wavelength location of in-orbit sun observation spectrum is calculated.It is herein
Based on the spectral calibration coefficient that Laboratory spectral calibration obtains, the centre wavelength of each channel on each spectral coverage is calculated.For
The EO-1 hyperion CO2 survey meters carried on carbon satellite, are calculated using following formula:
In formula, λ represents pixel centre wavelength (nm) match value, ciRepresent dispersion curve coefficient of polynomial fitting, channel
It represents probe access number, is 1-1242 for 0.76 mu m waveband, be 1-500 for 1.61 mu m wavebands and 2.06 mu m wavebands.
In one embodiment, in-orbit sun observation data described in step S1 are specially further that sun observation continues respectively
The average value of frame data.
It is to be understood that generally multiframe will be continued to sun observation, to eliminate the influence of random noise, multiframe can be used too
The average value of sun observation data carries out in-orbit spectral calibration.The specific average value that multiframe sun observation data are calculated using following formula:
In formula,Represent the average value of multiframe sun observation data, RiRepresent the i-th frame sun observation data, n represents the sun
Observation continues frame number.
Step S2 is it is to be understood that satellite speed of service in space up to 7km per second, can introduce and high-spectrum remote-sensing instrument
The frequency drift of device spectral resolution (0.03nm) quite magnitude.This is equivalent to channel wavelength and drifts about for EO-1 hyperion instrument
More than 50%, radiation calibration precision and basis weight products inverting can be generated and seriously affected.In addition, since satellite is in operational process
In, the variation of the environment such as temperature can also make the above-mentioned each channel center position being calculated that small drift can occur, this is right
It is influenced in basis weight products inversion accuracy very big.
Therefore, it is necessary to the basis of each channel center's wavelength location of in-orbit sun observation spectrum obtained is calculated in above-mentioned steps
On, the spectral position drift after in orbit is calculated, corrects the frequency drift.It is specific to consider Doppler shift characteristics,
Velocity information during being accurately positioned according to satellite calculates drift value, and according to the in-orbit of drift value corrected Calculation acquisition
Each channel center's wavelength location of sun observation spectrum, revised position are each channel active centre of in-orbit sun observation spectrum
Wavelength location.
Wherein optional, the S2's is further processed step with reference to figure 3, is that a kind of in-orbit sun of the embodiment of the present invention is seen
The acquisition process flow chart of each channel active centre wavelength location of spectrum is surveyed, including:
S21 is accurately positioned velocity information based on the satellite, using Doppler shift characteristics formula, calculates the target
The in-orbit doppler drift of high-spectrum remote-sensing instrument.
It is to be understood that realize the removal of drift value, need to first solve in drift value, then the theory obtained in above-mentioned calculating
It is accordingly corrected in heart wavelength location.Therefore, the velocity information during this step is accurately positioned according to satellite, utilizes physics
Doppler characterization formula in calculates the in-orbit doppler drift of target high-spectrum remote-sensing instrument.
It is wherein optional, the corresponding active centre observation frequency of the active centre wavelength location is calculated by the following formula
Rate:
In formula, vobsRepresent effective observation spectrum frequency, vsunRepresent sun standard spectrum frequency, VrelIt represents along EO-1 hyperion
Remote sensing instrument direction of visual lines satellite is with respect to the speed of the sun, c expression vacuum light speeds.
It is to be understood that relative to standard solar spectrum, after considering Doppler frequency shift, target high-spectrum remote-sensing Instrument observation
To the corresponding observing frequency of active centre wavelength location can be calculated by above formula.In addition, above formula can also be used as target EO-1 hyperion distant
The calculation formula of arbitrary effective frequency that sense Instrument observation arrives.During the step calculates, wavelength and frequency unit need to be by changing accordingly
It calculates, is converted to uniform units level.
S22 eliminates the in-orbit doppler drift in in-orbit each channel center's wavelength location of sun observation spectrum
Component obtains in-orbit each channel active centre wavelength location of sun observation spectrum.
It is to be understood that according to above-described embodiment obtain each channel center's wavelength location of in-orbit sun observation spectrum due to
It is influenced by Doppler effect and environment temperature etc., contains Doppler drift amount.This step calculates in above-mentioned steps and obtains the drift
It moves on the basis of component, each channel center's wavelength location of in-orbit sun observation spectrum is modified according to the drift component, is repaiied
Each channel active centre wavelength location of in-orbit sun observation spectrum is just being obtained afterwards.
Step S3 being calculated according to above-mentioned steps it is to be understood that obtain in-orbit each channel active centre of sun observation spectrum
After wavelength location, to realize the accurate in-orbit spectral calibration to target high-spectrum remote-sensing instrument, in-orbit sun observation need to be established
Each channel active centre wavelength location of spectrum and the incidence relation for the sun standard spectrum standard absorption line position stablized, so as to protect
Demonstrate,prove the stability of spectral calibration.
It is effective to calculate in-orbit each channel of sun observation spectrum on the basis of sun standard spectrum standard absorption line position for this step
The in-orbit spectral calibration error of centre wavelength position and the deviation of the benchmark, i.e. target high-spectrum remote-sensing instrument.That is, pass through ratio
Line position and high-spectrum remote-sensing Instrument observation position of spectral line are absorbed compared with sun standard spectrum, obtains the difference between the two, as in-orbit light
Spectrum calibration error.Therefore, it is necessary first to choose high-resolution sun standard spectrum.
Wherein optionally, the selection principle of sun standard spectrum includes in the step of S3:The sun standard spectrum
High resolution in the resolution ratio of the target high-spectrum remote-sensing instrument.
In one embodiment, the sun standard spectrum is further specifically chosen for Kurucz R.L. data.
It is to be understood that in order to meet the requirement of the spectral calibration precision of target high-spectrum remote-sensing instrument, in selection reference
During sun standard spectrum, high resolution is selected in the sun standard spectrum of the resolution ratio of target high-spectrum remote-sensing instrument.It is wherein high
The magnitude gone out can be set, and such as select the spectral resolution that sun standard spectrum has than the spectrum of target high-spectrum remote-sensing instrument
Resolution ratio is higher by 1-2 magnitude.For example, by using the Kurucz of high-spectrum remote-sensing, R.L. spectrum.Its spectral resolution reaches as high as
0.001nm, than carbon satellite remote sensing instrument highest spectral resolution 0.03nm high an order of magnitude.Kurucz is to sun Absorption Line
It portrays fine enough.
After high-resolution sun standard spectrum is selected, also need to determine the Absorption Line type of the standard spectrum, so as to clear
Each channel center's wavelength location of in-orbit sun observation spectrum of clear identification target high-spectrum remote-sensing instrument.
Wherein optional, the selecting step of standard absorption line described in step S3 further comprises:Absorption depth is chosen to reach
Width to setting rank and both wings reaches the Absorption Line of predetermined width as the standard absorption line.
In one embodiment, the standard absorption line is specially further that the bright standing grain of husband takes Absorption Line.
It is to be understood that in visible ray, near infrared light wave band, there is abundant Absorption Line, but not institute in solar spectrum
There is the standard that Absorption Line can serve as spectral calibration.Because some are absorbed, line absorption is weaker, and EO-1 hyperion instrument can observe
Signal it is also weaker.In addition, some Absorption Lines are too near apart, since EO-1 hyperion instrumental resolution is limited, it is difficult to clearly distinguish.Cause
This, needs to consider selection Absorption Line at least from the two indexs, that is, absorbs depth and absorption line width.
During specific choice, setting according to actual needs absorbs the pre-set level of depth and the predetermined width of Absorption Line.Choosing
The width for selecting absorption depth arrival pre-set level and Absorption Line both wings reaches the Absorption Line of predetermined width as standard absorption line.I.e.
Meet following impose a condition simultaneously in the Absorption Line for selecting sun standard spectrum:
Absorption Line has the absorption depth of setting rank, to ensure under the resolution ratio of target high-spectrum remote-sensing instrument, mesh
Absolute altitude spectral remote sensing instrument can be realized to its clear identification;
Absorption Line both wings have preset width, are not overlapped with ensureing that Absorption Line isolates.
Exist in visible ray and near infrared light spectral coverage, solar spectrum and take Absorption Line information, husband compared with the bright standing grain of the husband of horn of plenty
Lang He, which takes, absorbs that line absorption depth is strong, and adjacent absorbent line is easy to be identified by high-spectrum remote-sensing instrument at a distance of moderate, in an implementation
The bright standing grain of husband chosen in example in sun standard spectrum takes Absorption Line as standard absorption line.
According to above-mentioned established standards, in the spectral region of 758-778nm, 10 standard sun can be extracted and absorbed
Line;In the spectral region of 1594-1624nm, 8 standard sun Absorption Lines can be extracted;In the spectrum of 2041-2081nm
In the range of, 8 standard sun Absorption Lines can be extracted.
Wherein optional, the S3's is further processed step with reference to figure 4, is a kind of high-spectrum remote-sensing of the embodiment of the present invention
The calculating process flow chart of the in-orbit spectral calibration error of instrument, including:
S31, extract in in-orbit each channel active centre wavelength location of sun observation spectrum with the sun standard light
The corresponding active centre wavelength location of standard absorption line position is composed, by the corresponding active centre wavelength location and the sun
Standard spectrum standard absorption line position is matched.
It is to be understood that for the standard absorption line of sun standard spectrum, centre wavelength position is to determine, in order to ask
Observation data active centre wavelength location is taken to need to correspond to the position of the standard absorption line with respect to the deviation of the standard, respectively from
Corresponding active centre wavelength location is extracted in all each channel active centre wavelength locations of in-orbit sun observation spectrum.
Then by the correspondence centre wavelength position of extraction and the progress of the centre wavelength position of sun standard spectrum Absorption Line
Match.I.e. by revised in-orbit each channel active centre wavelength location of sun observation spectrum doppler shifted in step S2 and too
The centre wavelength position of positive standard spectrum Absorption Line is matched.
S32 at the centre wavelength position respectively to match, calculates each corresponding active centre wavelength location with respect to institute
State the wave length shift of sun standard spectrum standard absorption line position.
It is to be understood that by each channel active centre wavelength position of in-orbit sun observation spectrum after Doppler frequency shift modification
It puts after being matched with the centre wavelength position of sun standard spectrum Absorption Line, at the Absorption Line to match, with Absorption Line
Centre wavelength determines the wave length shift of EO-1 hyperion instrument, i.e., in-orbit spectral calibration error for standard.Following formula meter can specifically be passed through
It calculates:
Δ v=vo-vf;
In formula, Δ v represents the in-orbit spectral calibration error of high-spectrum remote-sensing instrument, voRepresent that high-spectrum remote-sensing instrument is in-orbit
The sun absorbs line position, v in observation spectrumfLine position is absorbed in representation theory solar spectrum.
Step S4 is it is to be understood that after above-mentioned steps calculate and obtain in-orbit spectral calibration error, according to the in-orbit light
Spectrum calibration error is modified the in-orbit spectral calibration of target high-spectrum remote-sensing instrument, realizes in-orbit spectrum precise calibration.
A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument provided in an embodiment of the present invention, by using sun standard
Spectrum calculates deviation between the in-orbit observation spectrum of target high-spectrum remote-sensing instrument and standard spectrum, so as to fulfill target as standard
The in-orbit spectral calibration of high-spectrum remote-sensing instrument can make full use of the stability feature of solar spectrum, so as to effectively improve
The accuracy and stability of spectral calibration.
On the basis of above-described embodiment, the present embodiment provides the in-orbit spectrum of high-spectrum remote-sensing instrument as shown in Figure 5 to determine
Mark method, Fig. 5 are the flow chart of the in-orbit spectrum calibration method of another kind high-spectrum remote-sensing instrument of the embodiment of the present invention.
By observation of the diffusing reflection plate to the sun when the present embodiment extracts satellite high-spectrum remote-sensing instrument in orbit first
Data, the calibration coefficient that Binding experiment room provides calculate preliminary in-orbit sun observation spectrum channel centre wavelength position;Then
On the basis of satellite is accurately positioned the velocity information provided, calculates doppler drift and see it from the preliminary in-orbit sun
It surveys in spectrum channel centre wavelength position and deducts, obtain effective in-orbit sun observation spectrum channel centre wavelength position.
Subsequently, Absorption Line selection gist is taken based on high-resolution theoretical sun standard spectrum and Fu Lang standing grain, selected
Available Absorption Line, and extract the center of these Absorption Lines;It finally will be in effective in-orbit sun observation spectrum channel
Heart wavelength location carries out location matches with sun standard spectrum theoretical absorption line, and the spectrum for calculating target high-spectrum remote-sensing instrument is determined
Mark error.On this basis, the in-orbit spectral calibration of high spectral resolution remote sensing instrument can be effectively realized.
It should be understood that the precondition of the present embodiment the method is target high-spectrum remote-sensing instrument have it is in-orbit too
Positive observing capacity.
The in-orbit spectrum calibration method of high-spectrum remote-sensing instrument of the present embodiment takes absorption based on the bright standing grain of sun standard spectrum husband
Line (Fraunhofer lines) carries out calibration calculating, can make full use of feature in sun original spectrum high resolution, spectrum
The characteristics of position of spectral line is stablized, accurately calculates in-orbit spectral calibration error, and monitor the stability of in-orbit spectral calibration.Effectively solution
Certainly high-spectrum remote-sensing instrument frequency can change with the variation of environment temperature, and spectral calibration is caused to evaluate the problem of unstable, light
Spectrum calibration precision higher, it is more stable.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
The present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that:It still can be right
Technical solution recorded in foregoing embodiments modifies or carries out equivalent replacement to which part technical characteristic;And this
A little modifications are replaced, the spirit and model of various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (10)
1. a kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument, which is characterized in that including:
It is high to calculate the target for S1, in-orbit sun observation data and ground calibration coefficient based on target high-spectrum remote-sensing instrument
Each channel center's wavelength location of in-orbit sun observation spectrum of spectral remote sensing instrument;
S2 is accurately positioned velocity information based on satellite, each to the in-orbit sun observation spectrum logical according to Doppler shift characteristics
Road centre wavelength position carries out frequency displacement amendment, obtains each channel active centre wavelength location of in-orbit sun observation spectrum;
S3, based on each channel active centre wavelength location of the in-orbit sun observation spectrum and sun standard spectrum standard absorption line
Position calculates the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument;
S4 based on the in-orbit spectral calibration error, carries out the in-orbit spectrum precise calibration of the target high-spectrum remote-sensing instrument.
2. according to the method described in claim 1, it is characterized in that, in the step of the S3 sun standard spectrum selection principle
Including:The high resolution of the sun standard spectrum is in the resolution ratio of the target high-spectrum remote-sensing instrument.
3. according to the method described in claim 2, it is characterized in that, the sun standard spectrum is further specifically chosen for
Kurucz R.L. data.
4. according to the method described in claim 2, it is characterized in that, the selecting step of standard absorption line described in step S3 is into one
Step includes:It chooses to absorb depth and reach the width of setting rank and both wings and reaches the Absorption Line of predetermined width and inhaled as the standard
Take-up.
5. according to the method described in claim 4, it is characterized in that, the standard absorption line, which is specially further the bright standing grain of husband, takes suction
Take-up.
6. according to the method described in claim 1, it is characterized in that, the step of S1 further specifically include:
The in-orbit sun observation data obtained in extracting the target high-spectrum remote-sensing instrument in orbit by diffusing reflection plate;
It is tested by the ground calibration to the target high-spectrum remote-sensing instrument, calculates the ground calibration coefficient;
Based on the in-orbit sun observation data and the ground calibration coefficient, using given fitting formula, calculate described in-orbit
Each channel center's wavelength location of sun observation spectrum.
7. according to the method described in claim 1, it is characterized in that, the step of S2 further comprise:
Velocity information is accurately positioned based on the satellite, using Doppler shift characteristics formula, it is distant to calculate the target EO-1 hyperion
Feel the in-orbit doppler drift of instrument;
The in-orbit doppler drift component in each channel center's wavelength location of in-orbit sun observation spectrum is eliminated, is obtained
Each channel active centre wavelength location of in-orbit sun observation spectrum.
8. the method according to the description of claim 7 is characterized in that the step of S3, further specifically includes:
It extracts in each channel active centre wavelength location of in-orbit sun observation spectrum and is inhaled with the sun standard spectrum standard
The corresponding active centre wavelength location in take-up position, by the corresponding active centre wavelength location and the sun standard spectrum
Standard absorption line position is matched;
At the centre wavelength position respectively to match, each relatively described sun mark of corresponding active centre wavelength location is calculated
The wave length shift of quasi-optical spectrum standard absorption line position.
9. the method according to the description of claim 7 is characterized in that the active centre wavelength location is calculated by the following formula
Corresponding active centre observing frequency:
In formula, vobsRepresent effective observation spectrum frequency, vsunRepresent sun standard spectrum frequency, VrelIt represents along high-spectrum remote-sensing
Instrument direction of visual lines satellite is with respect to the speed of the sun, c expression vacuum light speeds.
10. according to the method any in claim 1-9, which is characterized in that in-orbit sun observation number described in step S1
According to being specially further that sun observation continues the average value of each frame data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711284886.7A CN108254074B (en) | 2017-12-07 | 2017-12-07 | A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711284886.7A CN108254074B (en) | 2017-12-07 | 2017-12-07 | A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108254074A true CN108254074A (en) | 2018-07-06 |
CN108254074B CN108254074B (en) | 2019-11-26 |
Family
ID=62721014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711284886.7A Active CN108254074B (en) | 2017-12-07 | 2017-12-07 | A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108254074B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109374533A (en) * | 2018-12-18 | 2019-02-22 | 成都信息工程大学 | A kind of satellite short-wave infrared EO-1 hyperion observation data spectral drift evaluation method |
CN109781258A (en) * | 2018-12-29 | 2019-05-21 | 华中科技大学 | A kind of method that map association obtains sea-surface target infrared spectroscopy |
CN111006762A (en) * | 2019-12-10 | 2020-04-14 | 东华理工大学 | Space remote sensing instrument on-satellite calibration method based on LED piecewise fitting solar spectrum |
CN111624164A (en) * | 2020-03-19 | 2020-09-04 | 自然资源部第一海洋研究所 | Spectral standard deviation characteristic analysis-based hyperspectral identification method for marine spilled oil species |
CN111721734A (en) * | 2020-06-29 | 2020-09-29 | 中国科学院合肥物质科学研究院 | On-orbit spectrum calibration method for infrared very high spectral resolution detector for high-resolution five-number satellite |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120197096A1 (en) * | 2004-05-24 | 2012-08-02 | Trent Ridder | System for Noninvasive Determination of Analytes in Tissue |
CN105865624A (en) * | 2016-03-22 | 2016-08-17 | 清华大学 | Spectrum extraction method and device for hyperspectral collection system |
CN106017678A (en) * | 2016-06-12 | 2016-10-12 | 中国科学院上海技术物理研究所 | Thermal infrared high spectral remote sensing data on-track spectral calibration method |
CN107133976A (en) * | 2017-04-24 | 2017-09-05 | 浙江大学 | A kind of method and apparatus for obtaining three-dimensional hyperspectral information |
-
2017
- 2017-12-07 CN CN201711284886.7A patent/CN108254074B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120197096A1 (en) * | 2004-05-24 | 2012-08-02 | Trent Ridder | System for Noninvasive Determination of Analytes in Tissue |
CN105865624A (en) * | 2016-03-22 | 2016-08-17 | 清华大学 | Spectrum extraction method and device for hyperspectral collection system |
CN106017678A (en) * | 2016-06-12 | 2016-10-12 | 中国科学院上海技术物理研究所 | Thermal infrared high spectral remote sensing data on-track spectral calibration method |
CN107133976A (en) * | 2017-04-24 | 2017-09-05 | 浙江大学 | A kind of method and apparatus for obtaining three-dimensional hyperspectral information |
Non-Patent Citations (1)
Title |
---|
撖芃芃: "一种高光谱成像光谱仪光谱定标方法", 《长春工业大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109374533A (en) * | 2018-12-18 | 2019-02-22 | 成都信息工程大学 | A kind of satellite short-wave infrared EO-1 hyperion observation data spectral drift evaluation method |
CN109374533B (en) * | 2018-12-18 | 2021-04-09 | 成都信息工程大学 | Spectral drift estimation method for satellite short-wave infrared hyperspectral observation data |
CN109781258A (en) * | 2018-12-29 | 2019-05-21 | 华中科技大学 | A kind of method that map association obtains sea-surface target infrared spectroscopy |
CN111006762A (en) * | 2019-12-10 | 2020-04-14 | 东华理工大学 | Space remote sensing instrument on-satellite calibration method based on LED piecewise fitting solar spectrum |
CN111006762B (en) * | 2019-12-10 | 2021-09-24 | 东华理工大学 | Space remote sensing instrument on-satellite calibration method based on LED piecewise fitting solar spectrum |
CN111624164A (en) * | 2020-03-19 | 2020-09-04 | 自然资源部第一海洋研究所 | Spectral standard deviation characteristic analysis-based hyperspectral identification method for marine spilled oil species |
CN111624164B (en) * | 2020-03-19 | 2023-06-02 | 自然资源部第一海洋研究所 | Ocean spilled oil type hyperspectral identification method based on spectral standard deviation feature analysis |
CN111721734A (en) * | 2020-06-29 | 2020-09-29 | 中国科学院合肥物质科学研究院 | On-orbit spectrum calibration method for infrared very high spectral resolution detector for high-resolution five-number satellite |
Also Published As
Publication number | Publication date |
---|---|
CN108254074B (en) | 2019-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108254074B (en) | A kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument | |
Lebach et al. | Measurement of the solar gravitational deflection of radio waves using very-long-baseline interferometry | |
Mears et al. | Construction of the Remote Sensing Systems V3. 2 atmospheric temperature records from the MSU and AMSU microwave sounders | |
Pieters et al. | The Moon mineralogy mapper (M³) on chandrayaan-1 | |
Bignall et al. | Rapid interstellar scintillation of PKS 1257–326: two-station pattern time delays and constraints on scattering and microarcsecond source structure | |
CN101598798B (en) | System for rebuilding spectrum of high spectrum intervention data and method thereof | |
Trott et al. | Spectral calibration requirements of radio interferometers for epoch of reionisation science with the SKA | |
CN103364766B (en) | External calibration method of satellite-borne interferometric synthetic aperture radar (InSAR) system | |
CN105092055B (en) | Meteorological satellite sun reflected waveband Calibration Method based on cold cloud target | |
CN106054147A (en) | Microwave radiometer brightness temperature non-linear correction method and system | |
CN103616685B (en) | Based on the ISAR image geometry calibrating method of characteristics of image | |
Bendig et al. | Solar-induced chlorophyll fluorescence measured from an unmanned aircraft system: sensor etaloning and platform motion correction | |
Panciera et al. | A proposed extension to the soil moisture and ocean salinity level 2 algorithm for mixed forest and moderate vegetation pixels | |
Gao et al. | HJ-1A HSI on-orbit radiometric calibration and validation research | |
Mateus et al. | Mapping precipitable water vapor time series from Sentinel-1 interferometric SAR | |
CN103644907B (en) | A kind of pulsar angle-measurement system based on two satellite platform and method | |
CN110388986A (en) | Land surface temperature inversion method based on TASI data | |
Hosseini et al. | Combination of optical and SAR sensors for monitoring biomass over corn fields | |
CN110174652A (en) | A kind of calibrating method and device of the in-orbit emissivity of satellite-borne microwave imager antenna | |
Tian et al. | Analysis of vegetation index NDVI anisotropy to improve the accuracy of the GOES-R green vegetation fraction product | |
CN107478904A (en) | Atmospheric phase disturbance modification method based on global position system difference | |
Grody et al. | Synoptic-scale studies using the Nimbus 6 scanning microwave spectrometer | |
CN102519879B (en) | Total column ozone detection method | |
Berg et al. | Evaluation of monthly rainfall estimates derived from the special sensor microwave/imager (SSM/I) over the tropical Pacific | |
Zeng et al. | Intercalibration of FY-3C MWRI brightness temperature against GMI measurements based on ocean microwave radiative transfer model |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230717 Address after: 100081 No. 46 South Main Street, Haidian District, Beijing, Zhongguancun Patentee after: NATIONAL SATELLITE METEOROLOGICAL CENTER (NATIONAL CENTER FOR SPACE WEATHER) Address before: National Satellite Meteorological Center, No.46, Zhongguancun South Street, Haidian District, Beijing 100081 Patentee before: Bi Yanmeng |
|
TR01 | Transfer of patent right |