CN108254074B - 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
- CN108254074B CN108254074B CN201711284886.7A CN201711284886A CN108254074B CN 108254074 B CN108254074 B CN 108254074B CN 201711284886 A CN201711284886 A CN 201711284886A CN 108254074 B CN108254074 B CN 108254074B
- 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.)
- Active
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 68
- 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 7
- 238000012360 testing method Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 9
- 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
- 238000000605 extraction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 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
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification 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, it include: 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 is 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 is 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, fixed more particularly, to a kind of in-orbit spectrum of high-spectrum remote-sensing instrument
Mark method.
Background technique
It is distant to satellite with the continuous improvement of application demand, the especially continuous improvement of Atmospheric components quantitative detection demand
The requirement of the spectral resolution of sense is also higher and higher.Therefore, Modern Satellite remote sensing instrument development direction first is that have superelevation
Spectral resolution.
The carbon satellite of China's successful launch, by weak CO2 absorption band (1.6 μm), strong CO2 absorption band (2.06 μm) and
The observation spectrum of O2-A absorption band (0.76 μm) carries out high-precision CO2 concentration quantitative 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 central wavelength position in channel.This drift influences less for fat pipe instrument, but for
EO-1 hyperion instrument can but be generated and be seriously affected.The variation of channel spectrum center caused by small drift, 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
Technical point is wanted, is the important prerequisite of accurate radiation calibration.
Summary of the invention
In order to overcome the above problem or at least be partially solved the above problem, 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 and stability of spectral calibration.
The present invention provides a kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument, comprising: and S1 is distant based on target EO-1 hyperion
The in-orbit sun observation data and ground calibration coefficient for feeling instrument, the in-orbit sun for calculating the target high-spectrum remote-sensing instrument are 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 central wavelength position;S3 is 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, the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument is calculated;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, in the step of S3 the selection principle of sun standard spectrum include: the sun standard spectrum resolution
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: choosing 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 further specially that husband Lang He takes Absorption Line.
Wherein, the step of S1 is further specifically included: extract the target high-spectrum remote-sensing instrument in orbit in
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, intended using given
Formula is closed, in-orbit each channel center's wavelength location of sun observation spectrum is calculated.
Wherein, the step of S2 further comprises: being accurately positioned velocity information based on the satellite, utilizes Doppler's frequency
Characteristic formula is moved, the in-orbit doppler drift of the target high-spectrum remote-sensing instrument is calculated;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: during extraction in-orbit each channel of sun observation spectrum is effective
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, it is calculated by the following formula the corresponding active centre observing frequency of the active centre wavelength location:
In formula, vobsIndicate effective observation spectrum frequency, vsunIndicate sun standard spectrum frequency, VrelIt indicates along EO-1 hyperion
Speed of the remote sensing instrument direction of visual lines satellite with respect to the sun, c expression vacuum light speed.
Wherein, in-orbit sun observation data described in step S1 are further specially 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, to realize target EO-1 hyperion
The in-orbit spectral calibration of remote sensing instrument can make full use of the stability feature of solar spectrum, fixed so as to effectively improve spectrum
Target accuracy and stability.
Detailed description of the invention
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 process 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 process 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 of embodiment of the present invention high-spectrum remote-sensing instrument.
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 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
Every other embodiment obtained, shall fall within the protection scope of the present invention under the premise of creative work.
As one embodiment of the embodiment of the present invention, it is fixed that the present embodiment provides a kind of in-orbit spectrum of high-spectrum remote-sensing instrument
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 Fig. 1, comprising:
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 based on each channel active centre wavelength location of the in-orbit sun observation spectrum and sun standard spectrum standard is inhaled
Take-up position calculates the in-orbit spectral calibration error of the target high-spectrum remote-sensing instrument;
S4, is based on the in-orbit spectral calibration error, and 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 in ground experiment room it is to be understood that before satellite launch
Analogue scaling test, 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 having obtained above-mentioned ground calibration coefficient and in-orbit sun observation data, according to ground calibration coefficient
With in-orbit sun observation data, each channel center's wavelength in the in-orbit sun observation spectrum of target high-spectrum remote-sensing instrument is calculated
Position νorbit。
Wherein optional, it is that a kind of in-orbit sun of the embodiment of the present invention is seen that the S1's, which is further processed step with reference to Fig. 2,
Survey the calculating process flow chart of each channel center's wavelength location of spectrum, comprising:
S11 extracts the target high-spectrum remote-sensing instrument and is seen 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 is based on the in-orbit sun observation data and the ground calibration coefficient, using given fitting formula, calculates
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 light and near infrared light spectral coverage, satellite completes the earth observation of the area Yang Zhao
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 that satellite has in-orbit sun observation, 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 obtaining in-orbit sun observation data, the bracketing for carrying out ground experiment room is also needed, 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 being directed to target high-spectrum remote-sensing instrument.
Finally, calculating each channel center's wavelength location of in-orbit sun observation spectrum on the basis of above-mentioned steps.It is herein
Based on the spectral calibration coefficient that Laboratory spectral calibration obtains, the central wavelength in each channel on each spectral coverage is calculated.For
The EO-1 hyperion CO2 survey meter carried on carbon satellite, is calculated using following formula:
In formula, λ indicates pixel central wavelength (nm) match value, ciIndicate dispersion curve coefficient of polynomial fitting, channel
It indicates 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 further specially that sun observation continues respectively
The average value of frame data.
It is to be understood that multiframe, which generally will be continued, to sun observation can be used multiframe too for the influence for eliminating random noise
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,Indicate the average value of multiframe sun observation data, RiIndicate the i-th frame sun observation data, n indicates the sun
It observes and continues frame number.
Step S2 is it is to be understood that satellite speed of service in space can introduce and high-spectrum remote-sensing instrument up to 7km per second
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
50% or more, 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, the frequency drift is corrected.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, it is that a kind of in-orbit sun of the embodiment of the present invention is seen that the S2's, which is further processed step with reference to Fig. 3,
Survey the acquisition process flow chart of each channel active centre wavelength location of spectrum, comprising:
S21 is accurately positioned velocity information based on the satellite and calculates the target using Doppler shift characteristics formula
The in-orbit doppler drift of high-spectrum remote-sensing instrument.
It is to be understood that realize the removal of drift value, drift value need to be first solved, then in the theory that above-mentioned calculating obtains
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, it is calculated by the following formula the corresponding active centre observation frequency of the active centre wavelength location
Rate:
In formula, vobsIndicate effective observation spectrum frequency, vsunIndicate sun standard spectrum frequency, VrelIt indicates along EO-1 hyperion
Speed of the remote sensing instrument direction of visual lines satellite with respect to the sun, c expression vacuum light speed.
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, also to can be used as target EO-1 hyperion distant for above formula
The calculation formula for any 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 repaired
Each channel active centre wavelength location of in-orbit sun observation spectrum is just being obtained afterwards.
Step S3 is it is to be understood that obtain in-orbit each channel active centre of sun observation spectrum calculating according to above-mentioned steps
After wavelength location, to realize to the accurate in-orbit spectral calibration of target high-spectrum remote-sensing instrument, in-orbit sun observation need to be established
The incidence relation of each channel active centre wavelength location of spectrum and stable sun standard spectrum standard absorption line position, 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 deviation of central wavelength position and the benchmark, the i.e. in-orbit spectral calibration error of target high-spectrum remote-sensing instrument.That is, passing 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: the sun standard spectrum in the step of S3
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
When sun standard spectrum, select high resolution in the sun standard spectrum of the resolution ratio of target high-spectrum remote-sensing instrument.It is wherein high
Magnitude out can be set, spectrum of the spectral resolution for such as sun standard spectrum being selected to have than 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 selecting high-resolution sun standard spectrum, the Absorption Line type for determining the standard spectrum is needed, also 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: choosing absorption depth and reaches
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 further specially that husband Lang He takes Absorption Line.
It is to be understood that in visible light, near infrared light wave band there is Absorption Line abundant, but not institute in solar spectrum
There is Absorption Line to all can serve as the standard of 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 close apart, since EO-1 hyperion instrumental resolution is limited, it is difficult to clearly distinguish.Cause
This, needs to consider selection Absorption Line, i.e. absorption depth and absorption line width at least from the two indexs.
When specific choice, setting absorbs the pre-set level of depth and the predetermined width of Absorption Line according to actual needs.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 when selecting the Absorption Line of sun standard spectrum:
Absorption Line has the absorption depth of setting rank, to guarantee 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 guaranteeing that Absorption Line is isolated.
There is husband Lang He more abundant in visible light and near infrared light spectral coverage, solar spectrum and takes Absorption Line information, husband
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, implements at one
The husband Lang He chosen in sun standard spectrum in example 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 range, 8 standard sun Absorption Lines can be extracted.
Wherein optional, it is a kind of high-spectrum remote-sensing of the embodiment of the present invention that the S3's, which is further processed step with reference to Fig. 4,
The calculating process flow chart of the in-orbit spectral calibration error of instrument, comprising:
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, central wavelength position be it is determining, in order to ask
Take observation data active centre wavelength location with respect to the deviation of the standard, need to correspond to the position of the standard absorption line, 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 central wavelength position of extraction and the progress of the central 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 central wavelength position of positive standard spectrum Absorption Line is matched.
S32 calculates each corresponding active centre wavelength location with respect to institute at the central wavelength position respectively to match
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 sets after being matched with the central wavelength position of sun standard spectrum Absorption Line, at the Absorption Line to match, with Absorption Line
Central wavelength is the wave length shift that standard determines EO-1 hyperion instrument, i.e., in-orbit spectral calibration error.Following formula meter can specifically be passed through
It calculates:
Δ v=vo-vf;
In formula, Δ v indicates the in-orbit spectral calibration error of high-spectrum remote-sensing instrument, voIndicate 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, to realize 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 the above embodiments, it is fixed that the present embodiment provides the in-orbit spectrum of high-spectrum remote-sensing instrument as shown in Figure 5
Mark method, Fig. 5 are the flow chart of the in-orbit spectrum calibration method of another kind of embodiment of the present invention high-spectrum remote-sensing instrument.
Observation when the present embodiment extracts satellite high-spectrum remote-sensing instrument in orbit first by diffusing reflection plate to the sun
Data, the calibration coefficient that Binding experiment room provides calculate preliminary in-orbit sun observation spectrum channel central 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 central wavelength position and deducts, obtain effective in-orbit sun observation spectrum channel central 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 and sun standard spectrum theoretical absorption line carry out location matches, and the spectrum for calculating target high-spectrum remote-sensing instrument is fixed
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 be 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 sun standard spectrum husband Lang He
Line (Fraunhofer lines) carries out calibration calculating, can make full use of sun original spectrum high resolution, feature in spectrum
The stable feature of position of spectral line accurately calculates in-orbit spectral calibration error, and monitors the stability of in-orbit spectral calibration.Effectively solution
Certainly high-spectrum remote-sensing instrument frequency can change with the variation of environment temperature, cause spectral calibration to evaluate unstable problem, light
Spectrum calibration precision is higher, 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
Present invention has been described in detail with reference to the aforementioned embodiments, and those skilled in the art is it is understood that it still can be right
Technical solution documented by foregoing embodiments is modified or equivalent replacement of some of the technical features;And this
It modifies or replaces, the spirit and model of technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (7)
1. a kind of in-orbit spectrum calibration method of high-spectrum remote-sensing instrument characterized by comprising
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 central wavelength position carries out frequency displacement amendment, obtains each channel active centre wavelength location of in-orbit sun observation spectrum;
S3 is 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 is based on the in-orbit spectral calibration error, carries out the in-orbit spectrum precise calibration of the target high-spectrum remote-sensing instrument;
Wherein, the step of S1 is further specifically included:
Extract the in-orbit sun observation data that the target high-spectrum remote-sensing instrument is obtained in orbit by diffusing reflection plate;
By the ground calibration test to the target high-spectrum remote-sensing instrument, the ground calibration coefficient is calculated;
It is calculated described in-orbit based on the in-orbit sun observation data and the ground calibration coefficient using given fitting formula
Each channel center's wavelength location of sun observation spectrum;
The step of S2, further comprises:
Velocity information is accurately positioned based on the satellite, and it is distant to calculate the target EO-1 hyperion using Doppler shift characteristics formula
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;
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 central 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.
2. the method according to claim 1, wherein in the step of S3 sun standard spectrum selection principle
It include: the high resolution of the sun standard spectrum 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: to choose to absorb depth and reach the width of setting rank and both wings to reach the Absorption Line of predetermined width and inhale as the standard
Take-up.
5. according to the method described in claim 4, it is characterized in that, the standard absorption line is further specially that husband Lang He takes suction
Take-up.
6. the method according to claim 1, wherein being calculated by the following formula the active centre wavelength location
Corresponding active centre observing frequency:
In formula, vobsIndicate effective observation spectrum frequency, vsunIndicate sun standard spectrum frequency, VrelIt indicates along high-spectrum remote-sensing
Speed of the instrument direction of visual lines satellite with respect to the sun, c expression vacuum light speed.
7. any method in -6 according to claim 1, which is characterized in that in-orbit sun observation data described in step S1
It is further specially the average value that sun observation continues 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 CN108254074A (en) | 2018-07-06 |
CN108254074B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN111006762B (en) * | 2019-12-10 | 2021-09-24 | 东华理工大学 | Space remote sensing instrument on-satellite calibration method based on LED piecewise fitting solar spectrum |
CN111624164B (en) * | 2020-03-19 | 2023-06-02 | 自然资源部第一海洋研究所 | Ocean spilled oil type hyperspectral identification method based on spectral standard deviation feature analysis |
CN111721734B (en) * | 2020-06-29 | 2022-12-27 | 中国科学院合肥物质科学研究院 | On-orbit spectrum calibration method for infrared very-high spectral resolution detector for high-resolution five-number satellite |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8730047B2 (en) * | 2004-05-24 | 2014-05-20 | Trutouch Technologies, Inc. | System for noninvasive determination of analytes in tissue |
CN105865624B (en) * | 2016-03-22 | 2018-02-09 | 清华大学 | The spectrum extracting method and device of hyper-spectral data gathering system |
CN106017678B (en) * | 2016-06-12 | 2017-10-13 | 中国科学院上海技术物理研究所 | A kind of in-orbit spectrum calibration method of thermal infrared high-spectrum remote sensing data |
CN107133976B (en) * | 2017-04-24 | 2020-04-14 | 浙江大学 | Method and device for acquiring three-dimensional hyperspectral information |
-
2017
- 2017-12-07 CN CN201711284886.7A patent/CN108254074B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN108254074A (en) | 2018-07-06 |
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 | |
Asmar et al. | Spacecraft Doppler tracking: Noise budget and accuracy achievable in precision radio science observations | |
Mears et al. | Construction of the Remote Sensing Systems V3. 2 atmospheric temperature records from the MSU and AMSU microwave sounders | |
Bignall et al. | Rapid interstellar scintillation of PKS 1257–326: two-station pattern time delays and constraints on scattering and microarcsecond source structure | |
CN105092055B (en) | Meteorological satellite sun reflected waveband Calibration Method based on cold cloud target | |
CN102998665A (en) | Target radar cross section measuring and calibrating processing method | |
CN104614333B (en) | Measuring method of atmospheric water vapor molecule absorption coefficient profile | |
CN106054147A (en) | Microwave radiometer brightness temperature non-linear correction method and system | |
Dyrud et al. | Ionospheric measurement with GPS: Receiver techniques and methods | |
Corstanje et al. | The effect of the atmospheric refractive index on the radio signal of extensive air showers | |
CN103644907B (en) | A kind of pulsar angle-measurement system based on two satellite platform and method | |
Resch et al. | Radiometric correction of atmospheric path length fluctuations in interferometric experiments | |
Dempsey et al. | Line profile asymmetries in chromospherically active stars | |
Inai et al. | Altitude misestimation caused by the Vaisala RS80 pressure bias and its impact on meteorological profiles | |
Santana et al. | Development of procedures for calibration of meteorological sensors. Case study: calibration of a tipping-bucket rain gauge and data-logger set | |
Gaté et al. | Computing the electric field from extensive air showers using a realistic description of the atmosphere | |
CN110174652A (en) | A kind of calibrating method and device of the in-orbit emissivity of satellite-borne microwave imager antenna | |
Kalinnikov et al. | Determination of troposphere characteristics using signals of satellite navigation systems | |
CN102519879B (en) | Total column ozone detection method | |
Zeng et al. | Intercalibration of FY-3C MWRI brightness temperature against GMI measurements based on ocean microwave radiative transfer model | |
Corbella et al. | Some results on SMOS-MIRAS calibration and imaging | |
Travouillon et al. | Using a sodar to measure optical turbulence and wind speed for the thirty meter telescope site testing. Part II: Comparison with independent instruments | |
Gao et al. | Method study and uncertainty analysis of calibration coefficients validation based on the Inner Mongolia test site | |
Zakharov et al. | Estimation of the relative humidity in the troposphere using GNSS signals |
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 |