CN109712089A - Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite - Google Patents
Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite Download PDFInfo
- Publication number
- CN109712089A CN109712089A CN201811535873.7A CN201811535873A CN109712089A CN 109712089 A CN109712089 A CN 109712089A CN 201811535873 A CN201811535873 A CN 201811535873A CN 109712089 A CN109712089 A CN 109712089A
- Authority
- CN
- China
- Prior art keywords
- calibration
- light intensity
- pixel
- data
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Image Processing (AREA)
Abstract
The invention discloses a kind of methods suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, by utilizing 90 ° of calibration factures of yaw on in-orbit star, parameter needed for generating infrared remote sensing image, relative detector calibration coefficient data file including different gains, different spectral coverage, establish the quantitative relationship between the realistic objective radiance value in the numerical quantization value pixel corresponding with the detector of each probe unit output signal, each response nonuniformity visited between member is eliminated, to ensure the reliability and accuracy of remote sensing data application.
Description
Technical field
The invention belongs to Remote Sensing Image Processing Technology field, more particularly to one kind are infrared short suitable for export-oriented remote sensing satellite
The method of wave load relative detector calibration.
Background technique
For remote sensing satellite infrared payload mounted, transmitting front and back all must output and acquisition to instrument data into
The comprehensive radiation calibration of row, to ensure that data collected are applied to greatest extent.During in satellier injection and in orbit,
In view of being mingled with noise energy, such as background, atmosphere and instrument itself transmitting in the received infrared emanation energy of infrared payload
Radiation, these can be affected to the measurement of Target Infrared Radiation characteristic, so as to cause infrared payload performance often with
There are certain deviations for test result before transmitting, therefore, for infrared payload data carry out relative detector calibration, be at present urgently
Major issue to be solved, this is most important to the quality of Satellite Product.
Although having carried out the radiation calibration research in building for the satellite infrared load of China's independent development, due to
Remote sensing strength such as is tested less input, experimental site type is not complete, subject crossing is inadequate, basic research is weak at the factors by China, sternly
The original development with full quantification of theory and technology of system is affected again.In terms of radiating radiometric calibration site, tentatively build up
Dunhuang land bracketing field and Qinghai Lake water surface bracketing field, and respectively for China FY-1C, FY-2B meteorological satellite,
The satellites such as FY-1D, HY-1 seasat and CBERS-1 landsat have carried out place external calibration method test, but due to marching into the arena
Number and weather condition limitation, are only limitted to the correction work of part remote sensor, and bearing calibration used is on the low side, main at present
Based on traditional wavelet analysis method, this method is mainly using the advantage on processing jump signal, and limitation is bigger,
It is far from satisfying the needs of practical application, the production of the pretreatment work and standardized product of satellite data of direct image.Mirror
In this, this patent proposes a kind of method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, eliminates
Each response nonuniformity visited between member solves the infrared shortwave load image of export-oriented remote sensing satellite and sweeps present on direction pushing away
Strip phenomenon.
Summary of the invention
Technical problem solved by the present invention is having overcome the deficiencies of the prior art and provide a kind of suitable for export-oriented remote sensing
The method of satellite infrared shortwave load relative detector calibration, this method can correct the factor for influencing infrared signature, solve
Strip phenomenon on infrared remote sensing image.
The object of the invention is achieved by the following technical programs: one kind being suitable for the infrared shortwave of export-oriented remote sensing satellite and carries
The method of lotus relative detector calibration, described method includes following steps: (1) the in-orbit yawed flight of satellite is preset, it is former according to acquisition
90 ° of calibrations of primary yaw then are completed under substar imaging time of integration setting, focal plane detector acquires original calibration data F；
(2) according to the imaging sequence of the pixel number N of the focal plane detector of satellite and preset pixel, different pixels are according in step (1)
Original calibration data F rearranged to obtain the data value on several rows, wherein data value be different pixels to same
The response of atural object imaging；(3) translation processing is carried out to the data value on several rows in step (2), obtains the mean value of every a line,
Compare the mean value of every a line and the preset threshold of corresponding row, it, should if the mean value of the row is greater than the preset threshold of the row
Row is rejected, and effective calibration data G is finally obtained；(4) it according to effective calibration data G of step (3), constructs standard and adds up histogram
Scheme and standard accumulative histogram is normalized to obtain normalization accumulative histogram, and then obtains the every of focal plane detector
The accumulative histogram of a detection member；(5) the accumulative histogram combination focal plane detector of each detection member for obtaining step (4)
Laboratory relative detector calibration model constructs relative detector calibration look-up table；(6) according to the relative detector calibration in step (5)
Look-up table carries out relative detector calibration to original image to be measured.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (1),
The acquisition principle are as follows: a) tonal range being imaged covers camera dynamic range；B) unsaturated scenery in imaging；C) known
The quantization digit n of line number L, focal plane detector total pixel number N and image data is imaged, then L > n+100*2n.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (3),
Translation processing is carried out to the data value on several rows in step (2) to obtain by following formula:
G (i, j)=F (i+j-1, j)；
Wherein, i is focal plane detector pixel line number, and j is focal plane detector pixel columns, and i=1,2,3 ..., L, L are total
Line number；J=1,2,3 ..., N/L, G (i, j) are response numerical value of the pixel to imaging of effective calibration data, and F (i+j-1, j) is
Response numerical value of the pixel of original calibration data to imaging.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (3),
The mean value of every a line is obtained by following formula:
Wherein, MiFor the mean value of the i-th row, work as Mi>2nThis journey data are saturation data when -1, then reject this journey.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (3),
Preset threshold is 2n- 1, wherein n is the quantization digit of image data.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (4),
Standard accumulative histogram is obtained by following formula:
Wherein, pixel scale number q1,q2,...,qk∈ [0,1], k=1,2,3 ..., K, k are the histogram for detecting pixel
Group column, K are detection pixel set of histograms column sum, [DNmin,DNmax] it is MiImages light intensity value range；NbIt is total after rejecting
Line number；For the images light intensity value of focal plane detector laboratory relative radiation；NykFor MiMiddle images light intensity value is less thanPixel
Number.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (4),
Normalization accumulative histogram is obtained by following formula:
Wherein, h[j][X]For the normalized images light intensity value of j-th of detection member, NxIt is j-th in effective calibration data G
Detect the number of pixels that first images light intensity value is less than X；NbFor total line number after rejecting；X is variable, value range [DNmin,
DNmax]。
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (5),
Relative detector calibration look-up table is obtained by following steps:
For j-th of detection member:
As the images light intensity value x ∈ [DN of original imageminj,DNmaxj] when, use radiant correction look-up tableIt looks for
To the corresponding images light intensity value Y of X, Y is assigned to the corresponding position images light intensity value after correction and carries out relative detector calibration；Its
In, k=1,2,3 ..., K, k are the histogram sequence number for detecting pixel, and K is the detection total sequence number of pixel histogram,To have
The accumulative histogram of each detection member of calibration data is imitated,For the images light intensity value of focal plane detector laboratory relative radiation,
[DNmin,DNmax] it is images light intensity value range；
As the images light intensity value x ∈ [0, DN of original imageminj] or x ∈ [DNmaxj,2nWhen -1], searched using radiant correction
TableFind the corresponding images light intensity value Y ' of Xr, Y 'rThe images light intensity value of corresponding position after being assigned to correction carries out
Relative detector calibration；Wherein,It in radiation flux is Φ for j-th of detection pixelrWhen the radiation flux that exports, Y 'rFor experiment
Images light intensity value of the room relative detector calibration model in r detection member, r=1,2,3 ..., DNminjOr r=DNmaxj,
DNmaxj+1,DNmaxj+2,...,2n- 1, j=1,2,3 ..., N/L, wherein L is imaging line number.
It is above-mentioned suitable for the method for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, in step (6),
Carrying out relative detector calibration to original image to be measured according to the relative detector calibration look-up table in step (5) includes: to utilize radiation
Correction look-up table finds the images light intensity value of corresponding original image to be measured, establishes each numerical quantization for detecting first output signal
The quantitative relationship being worth between the realistic objective radiance value in pixel corresponding with the detector eliminates each response visited between member
Inconsistency.
Compared with prior art, the present invention has the following advantages:
The present invention passes through using 90 ° of calibration factures are yawed on in-orbit star, and parameter needed for generating infrared remote sensing image is wrapped
The relative detector calibration coefficient data file for including different gains, different spectral coverage, establishes the number of each probe unit output signal
The quantitative relationship being worth between the realistic objective radiance value in quantized value pixel corresponding with the detector is eliminated between each spy member
Response nonuniformity, to ensure the reliability and accuracy of remote sensing data application.
Detailed description of the invention
By reading the following detailed description of the preferred embodiment, various other advantages and benefits are common for this field
Technical staff will become clear.The drawings are only for the purpose of illustrating a preferred embodiment, and is not considered as to the present invention
Limitation.And throughout the drawings, the same reference numbers will be used to refer to the same parts.In the accompanying drawings:
Fig. 1 is 90 ° of calibration process flow diagrams of yaw；
Fig. 2 is 90 ° of calibration mode schematic diagrames of yaw；
Fig. 3 is 90 ° of calibration time diagrams of yaw；
Fig. 4 (a) is the schematic diagram of original calibration data；Fig. 4 (b) is the schematic diagram of the calibration data after extracting；
Fig. 5 (a) is the schematic diagram of data before radiant correction；Fig. 5 (b) is the schematic diagram of data after radiant correction.
Specific embodiment
Exemplary embodiments of the present disclosure are described in more detail below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
It is limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure
It is fully disclosed to those skilled in the art.It should be noted that in the absence of conflict, embodiment in the present invention and
Feature in embodiment can be combined with each other.The present invention will be described in detail below with reference to the accompanying drawings and embodiments.
A kind of method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite is present embodiments provided,
Steps are as follows:
(1) the in-orbit yawed flight of satellite is set, completes one under substar imaging time of integration setting according to acquisition principle
90 ° of calibrations of secondary yaw, focal plane detector acquire original calibration data F.Acquire principle are as follows: a) tonal range being imaged covers camera
Dynamic range；B) unsaturated scenery in imaging；C) line number L is imaged, it is assumed that always first number is N, image data to focal plane detector
Quantization digit be n, then L > n+100*2n；
(2) according to the imaging sequence of the pixel number N of the focal plane detector of satellite and preset pixel, different pixels are according to step
Suddenly the original calibration data F in (1) is rearranged to obtain the data value on several rows, wherein, data value is different pictures
The response that same atural object is imaged in member；
(3) translation processing is carried out for the data value on several rows in step (2), obtains the mean value of every a line, compares
The preset threshold of the mean value of every a line and corresponding row picks the row if the mean value of the row is greater than the preset threshold of the row
It removes, and finally obtains effective calibration data G；
(4) the effective calibration data G extracted according to step (3), calculates the mean value of each row of data, and building standard adds up histogram
Scheme and standard accumulative histogram is normalized to obtain to normalization accumulative histogram, and then obtains focal plane detector
The accumulative histogram of each detection member, the specific method is as follows:
A. for the effective calibration data G extracted, the mean value M of each row of data is calculatedi, i=1,2,3 ..., L, L are head office
Number；
B. standard accumulative histogram is constructed, wherein pixel scale number q1,q2,...,qk∈ [0,1], construction method is as follows:
Wherein, group column k=1,2,3 ..., K, k are the histogram sequence number for detecting pixel, and K is that detection pixel histogram is total
Sequence number, [DNmin,DNmax] it is MiImages light intensity value range；NbFor total line number after rejecting；For focal plane detector experiment
The images light intensity value of room relative radiation；NykFor MiMiddle images light intensity value is less thanNumber of pixels.
C. normalized accumulative histogram is calculated for j-th of detection member:
Wherein, h[j][X]For the normalized images light intensity value of j-th of detection member, NxIt is j-th in effective calibration data G
Detect the number of pixels that first images light intensity value is less than X；NbFor total line number after rejecting；X is variable, value range [DNmin,
DNmax]。
D. the accumulative histogram of each detection member of effective calibration dataSo thatI.e. j-th spy
It surveys in member and is less thanPixel sum account for all pixel sums ratio be qk:
Wherein, group column k=1,2,3 ..., K,For h[j][X]Inverse function, j=1,2,3 ..., N/L, wherein L be at
As line number.
(5) the laboratory relative detector calibration model for the accumulative histogram combination focal plane detector for obtaining step (4), structure
Relative detector calibration look-up table is built, the specific method is as follows:
For j-th of detection member:
A. as DN value x ∈ [DNmin,DNmax] when, byCorresponding relationship establish this detection member radiant correction
Look-up table:
Wherein, group column k=1,2,3 ..., K,For the accumulative histogram of each detection member of effective calibration data,
For the images light intensity value of focal plane detector laboratory relative radiation,For h[j][X]Inverse function, j=1,2,3 ..., N/L,
Middle L is imaging line number；The two is established into corresponding relationship.
B. as DN value x ∈ [0, DNminj] and as DN value x ∈ [DNmaxj,2nWhen -1], building [0, DNminj] and [DNmaxj,2n-
1] the radiant correction look-up table in range:
1) the relative detector calibration model in laboratory is constructed
Wherein, N is detection pixel sum,It in radiation flux is Φ for j-th of detection pixelrWhen the radiation that exports it is logical
Amount, YrIt in radiation flux is Φ for j-th of detection memberrOutput after relative correction.
2) as DN value x ∈ [DNminj, DNmaxj] when, byCorresponding relationship and laboratory relative detector calibration mould
TypeCount the calibration model variable quantity Z after entering the orbitk:
Wherein, group column k=1,2,3 ..., K, Y 'kFor the output valve after relative correction,For focal plane detector laboratory
The images light intensity value of relative radiation, j-th of detection member, j=1,2,3 ..., N/L, wherein L is imaging line number.
It enables:ByIt is defeated after radiation flux is relative correction to obtain j-th of detection member at this time
It is denoted as Y ' outk, and then obtain
3) [0, DN is constructedminj] and [DNmaxj,2n- 1] look-up table in range:
Y′r=Yr·Z
Wherein, X is variable, value range [DNmin,DNmax],It in radiation flux is Φ for j-th of detection pixelrWhen
The radiation flux of output, YrIt in radiation flux is Φ for j-th of detection memberrOutput after relative correction, ZkFor calibration model variation
Amount, Y 'kFor the output valve after relative correction,For the images light intensity value of focal plane detector laboratory relative radiation, Z is after entering the orbit
Radiant correction model, Y 'rImages light intensity value for laboratory relative detector calibration model in r detection member, r=1,2,
3,...,DNminjOr r=DNmaxj,DNmaxj+1,DNmaxj+2,...,2n- 1, j=1,2,3 ..., N/L, wherein L is imaging
Line number.
(6) look-up table for acquiring step (5) carries out relative detector calibration to original image to be measured.Specific method is to utilize
Radiant correction look-up table finds the images light intensity value of corresponding original image to be measured, establishes each numerical value for detecting first output signal
The quantitative relationship between realistic objective radiance value in quantized value pixel corresponding with the detector is eliminated between each spy member
Response nonuniformity.
Strip phenomenon, this hair present on direction are swept in order to solve the infrared shortwave load image of export-oriented remote sensing satellite pushing away
Bright proposition is using 90 ° of scaling methods of yaw on in-orbit star, parameter needed for generating infrared remote sensing image, including different gains, difference
The relative detector calibration coefficient data file of spectral coverage establishes the numerical quantization value and the detection of each probe unit output signal
Device corresponds to the quantitative relationship between the realistic objective radiance value in pixel, eliminates each response nonuniformity visited between member,
To ensure the reliability and accuracy of remote sensing data application.
As shown in Figure 1, main contents are as follows for 90 ° of calibration process flow diagrams of yaw:
The present invention proposes that using 90 ° of calibration factures of yaw, main flow mainly includes six steps, the first step, setting
The in-orbit yawed flight of satellite completes 90 ° of calibrations of primary yaw under substar imaging time of integration setting, selects original scalar number
According to；Second step obtains the imaging sequence of pixel according to the pixel number of focal plane detector, different pixels to same scenery at
As data permutation, finally data on a same row are the data that same atural object is imaged in different pixels；Third step, for
Original calibration data carry out translation processing, calculate the mean value of every a line, the saturation data in image data after rejecting translation, and
Finally obtain effective calibration data；4th step calculates the mean value of each row of data, building standard accumulative histogram and detection member meter
Normalization accumulative histogram is calculated, and then obtains the accumulative histogram of each detection member；5th step, Binding experiment room relative radiation school
Positive model constructs relative detector calibration look-up table；6th step carries out relative detector calibration to original image using look-up table.
As shown in Fig. 2, main contents are as follows for 90 ° of calibration mode schematic diagrames of yaw:
It when ± 90 ° of flight of the in-orbit yaw of satellite is arranged in the present invention, realizes that the scanning direction of camera is parallel with linear array, owns
Detection member obtains identical spoke brightness by the same target in ground, in this way can one virtual uniform field of in-orbit formation with
Improve radiant correction precision.
As shown in figure 3, main contents are as follows for 90 ° of calibration time diagrams of yaw:
The present invention selects original calibration data F, does translation processing to original calibration data, extracts effective calibration data.Root
According to satellite under substar push-scanning image mode, time of integration tintCalculation formula are as follows:
Wherein, GSD is ground sampling resolution, Vsat_earthFor satellite ground flying speed, H is orbit altitude, dpixelFor
Camera detector pixel dimension, f are camera focus, and r is earth radius, and μ is Gravitational coefficient of the Earth.Timing is calibrated according to Fig. 3 to show
Shown in intention, completion 90 ° of calibrations of primary yaw, the shortest time theoretically needed are under substar imaging time of integration setting
tint× N, wherein N is total pixel number.
It is main interior for the calibration data (right side) behind original calibration data (left side) and extraction as shown in Fig. 4 (a) and Fig. 4 (b)
Hold as follows:
The present invention carries out translation processing for original calibration data F, calculates the mean value of every a line, the image after rejecting translation
Saturation data in data, and effective calibration data are finally obtained, it is embodied as follows:
Pixel arrangement mode is as follows:
A. following translation processing is carried out for original calibration data F:
G (i, j)=F (i+j-1, j)
Wherein, i is focal plane detector pixel line number, and j is focal plane detector pixel columns, and i=1,2,3 ..., L, L are total
Line number；J=1,2,3 ..., N/L, G (i, j) are response numerical value of effective calibration data pixel to imaging, and F (i+j-1, j) is original
Response numerical value of the beginning calibration data pixel to imaging.
B. the saturation data in image data after rejecting translation, and finally obtain effective calibration data G；It calculates each
Capable mean value:
Work as Mi>2nThis journey data are saturation data when -1, then reject this journey.
As shown in Fig. 5 (a) and Fig. 5 (b), data (right side), main contents after data (left side) and radiant correction before radiant correction
It is as follows:
The present invention finds corresponding original image images light intensity value using radiant correction look-up table, pair being assigned to after correcting
It answers the images light intensity value of position to carry out radiant correction, establishes the numerical quantization value and the detection of each probe unit output signal
Device corresponds to the quantitative relationship between the realistic objective radiance value in pixel, eliminates each response nonuniformity visited between member.
Specific step is as follows:
For j-th of detection member:
A. as the images light intensity value x ∈ [DN of original imageminj,DNmaxj] when, use radiant correction look-up table
The corresponding images light intensity value Y of X is found, Y is assigned to the corresponding position images light intensity value after correction and carries out relative detector calibration.Its
In, k=1,2,3 ..., K, k are the histogram sequence number for detecting pixel, and K is the detection total sequence number of pixel histogram,To have
The accumulative histogram of each detection member of calibration data is imitated,For the images light intensity value of focal plane detector laboratory relative radiation,
[DNmin,DNmax] it is images light intensity value range.
B. as the images light intensity value x ∈ [0, DN of original imageminj] or x ∈ [DNmaxj,2nWhen -1], looked into using radiant correction
Look for tableFind the corresponding images light intensity value Y ' of Xr, Y 'rBe assigned to correction after corresponding position images light intensity value into
Row relative detector calibration.Wherein,It in radiation flux is Φ for j-th of detection pixelrWhen the radiation flux that exports, Y 'rFor reality
Room relative detector calibration model is tested in the images light intensity value of r detection member, r=1,2,3 ..., DNminjOr r=DNmaxj,
DNmaxj+1,DNmaxj+2,...,2n- 1, j=1,2,3 ..., N/L, wherein L is imaging line number.
Embodiment
In order to verify the performance of infrared shortwave load relative radiometric correction method, using on in-orbit star yaw 90 ° of scaling methods,
Parameter needed for generating infrared remote sensing image establishes the numerical quantization value and the detector pair of each probe unit output signal
The quantitative relationship between the realistic objective radiance value in pixel is answered, each response nonuniformity visited between member is eliminated, passes through
The image of comparison before and after the processing, develops practicability of the invention.
According to Fig. 4, calibration data (right side) after comparing original calibration data (left side) and extracting, can show will be original fixed
Mark data carry out translation processing, calculate the mean value of every a line, after the saturation data in image data after rejecting translation, can obtain
To effective calibration data.
According to Fig. 5, corresponding original image images light intensity value is found using radiant correction look-up table, after being assigned to correction
The DN value of corresponding position carries out radiant correction, establishes the numerical quantization value and the detector pair of each probe unit output signal
Answer the quantitative relationship between the realistic objective radiance value in pixel.By comparison diagram 5, the striped of the image after radiant correction
It has substantially eliminated.
It can be found by above-mentioned verification result, the method for the present invention can be realized effective opposite spoke for infrared remote sensing image
Correction process is penetrated, may be implemented to handle the relative detector calibration of infrared remote sensing image using the method for the present invention.The invention is one
The method that kind is suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite.
Embodiment described above is the present invention more preferably specific embodiment, and those skilled in the art is in this hair
The usual variations and alternatives carried out in bright technical proposal scope should be all included within the scope of the present invention.
Claims (9)
1. a kind of method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite, which is characterized in that described
Method includes the following steps:
(1) the in-orbit yawed flight of satellite is preset, is completed under substar imaging time of integration setting according to acquisition principle primary inclined
It navigates 90 ° and calibrates, focal plane detector acquires original calibration data F；
(2) according to the imaging sequence of the pixel number N of the focal plane detector of satellite and preset pixel, different pixels are according to step
(1) the original calibration data F in is rearranged to obtain the data value on several rows, wherein data value is different pixels
Response to the imaging of same atural object；
(3) translation processing is carried out to the data value on several rows in step (2), obtains the mean value of every a line, more every a line
Mean value and corresponding row preset threshold, if the mean value of the row be greater than the row preset threshold, which is rejected, finally
Obtain effective calibration data G；
(4) it according to effective calibration data G of step (3), constructs standard accumulative histogram and standard accumulative histogram is returned
One change handles to obtain normalization accumulative histogram, and then obtains the accumulative histogram of each detection member of focal plane detector；
(5) the laboratory relative radiation school of the accumulative histogram combination focal plane detector for each detection member for obtaining step (4)
Positive model constructs relative detector calibration look-up table；
(6) relative detector calibration is carried out to original image to be measured according to the relative detector calibration look-up table in step (5).
2. the method according to claim 1 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: in step (1), the acquisition principle are as follows: a) tonal range being imaged covers camera dynamic range；B) it is imaged
In unsaturated scenery；C) the quantization digit n of line number L, focal plane detector total pixel number N and image data are imaged known to, then L
>n+100*2n。
3. the method according to claim 1 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: carrying out translation processing in step (3) to the data value on several rows in step (2) and passing through following formula
It obtains:
G (i, j)=F (i+j-1, j)；
Wherein, i is focal plane detector pixel line number, and j is focal plane detector pixel columns, and i=1,2,3 ..., L, L are head office
Number；J=1,2,3 ..., N/L, G (i, j) are response numerical value of the pixel to imaging of effective calibration data, and F (i+j-1, j) is original
Response numerical value of the pixel of beginning calibration data to imaging.
4. the method according to claim 3 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: the mean value of every a line is obtained by following formula in step (3):
Wherein, MiFor the mean value of the i-th row, work as Mi>2nThis journey data are saturation data when -1, then reject this journey.
5. the method according to claim 4 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: in step (3), preset threshold 2n- 1, wherein n is the quantization digit of image data.
6. the method according to claim 4 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: standard accumulative histogram is obtained by following formula in step (4):
Yqk=DNmin+k-1；
Wherein, pixel scale number q1,q2,...,qk∈ [0,1], k=1,2,3 ..., K, k are the set of histograms column for detecting pixel,
K is detection pixel set of histograms column sum, [DNmin,DNmax] it is MiImages light intensity value range；NbFor total line number after rejecting；
YqkFor the images light intensity value of focal plane detector laboratory relative radiation；NykFor MiMiddle images light intensity value is less than YqkNumber of pixels.
7. the method according to claim 6 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: normalization accumulative histogram is obtained by following formula in step (4):
Wherein, h[j][X]For the normalized images light intensity value of j-th of detection member, NxFor j-th of detection member in effective calibration data G
Images light intensity value is less than the number of pixels of X；NbFor total line number after rejecting；X is variable, value range [DNmin,DNmax]。
8. the method according to claim 7 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: relative detector calibration look-up table is obtained by following steps in step (5):
For j-th of detection member:
As the images light intensity value x ∈ [DN of original imageminj,DNmaxj] when, use radiant correction look-up tableFind X
Y is assigned to the corresponding position images light intensity value after correction and carries out relative detector calibration by corresponding images light intensity value Y；Wherein, k
=1,2,3 ..., K, k be detect pixel histogram sequence number, K be detection the total sequence number of pixel histogram,It is effectively fixed
The accumulative histogram of each detection member of data is marked,For the images light intensity value of focal plane detector laboratory relative radiation,
[DNmin,DNmax] it is images light intensity value range；
As the images light intensity value x ∈ [0, DN of original imageminj] or x ∈ [DNmaxj,2nWhen -1], radiant correction look-up table is usedFind the corresponding images light intensity value Y of Xr', YrThe images light intensity value of ' the corresponding position being assigned to after correction carries out phase
To radiant correction；Wherein,It in radiation flux is Φ for j-th of detection pixelrWhen the radiation flux that exports, Y 'rFor laboratory
Images light intensity value of the relative detector calibration model in r detection member, r=1,2,3 ..., DNminjOr r=DNmaxj,DNmaxj+
1,DNmaxj+2,...,2n- 1, j=1,2,3 ..., N/L, wherein L is imaging line number.
9. the method according to claim 1 suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite,
It is characterized by: carrying out phase to original image to be measured according to the relative detector calibration look-up table in step (5) in step (6)
Include: the images light intensity value for finding corresponding original image to be measured using radiant correction look-up table to radiant correction, establishes each
Detect quantifying between the realistic objective radiance value in the numerical quantization value pixel corresponding with the detector of first output signal
Relationship eliminates each response nonuniformity visited between member.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811535873.7A CN109712089A (en) | 2018-12-14 | 2018-12-14 | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811535873.7A CN109712089A (en) | 2018-12-14 | 2018-12-14 | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109712089A true CN109712089A (en) | 2019-05-03 |
Family
ID=66256564
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811535873.7A Pending CN109712089A (en) | 2018-12-14 | 2018-12-14 | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109712089A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110782429A (en) * | 2019-09-26 | 2020-02-11 | 长沙天仪空间科技研究院有限公司 | Imaging quality evaluation method based on satellite-borne remote sensing camera |
| CN111257238A (en) * | 2020-02-05 | 2020-06-09 | 国家卫星气象中心(国家空间天气监测预警中心) | Detection element relative calibration method based on satellite-borne solar diffusion plate |
| CN111415392A (en) * | 2020-03-13 | 2020-07-14 | 湖北师范大学 | Video satellite in-orbit relative radiation calibration method adopting Bell template push-broom mode imaging |
| CN113177512A (en) * | 2021-05-20 | 2021-07-27 | 国家卫星气象中心(国家空间天气监测预警中心) | Matching threshold analysis method for intersatellite cross radiometric calibration |
| CN113570523A (en) * | 2021-07-31 | 2021-10-29 | 中国人民解放军61646部队 | Method for automatic generation and extrapolation of relative radiation correction coefficients for optical images |
| CN113984673A (en) * | 2021-10-25 | 2022-01-28 | 中国科学院合肥物质科学研究院 | Satellite-borne polarization scanning remote sensor data pre-correction method and device |
| CN114125431A (en) * | 2021-11-22 | 2022-03-01 | 北京市遥感信息研究所 | Non-uniformity calibration correction method for static track optical large-area array camera |
| CN116539542A (en) * | 2023-06-20 | 2023-08-04 | 北京朗森基科技发展有限公司 | Ball-borne imaging Fourier transform spectrometer for atmospheric trace gas analysis |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901516A (en) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | Multispectral image radiation correction method based on absolute radiometric calibration |
| CN103778637A (en) * | 2014-01-24 | 2014-05-07 | 北京空间机电研究所 | Side-slither radiometric calibration processing method based on histogram equalization |
| CN106871925A (en) * | 2016-12-29 | 2017-06-20 | 中国资源卫星应用中心 | A kind of remote sensing satellite relative radiometric calibration processing method of in-orbit comprehensive dynamic adjustment |
-
2018
- 2018-12-14 CN CN201811535873.7A patent/CN109712089A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901516A (en) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | Multispectral image radiation correction method based on absolute radiometric calibration |
| CN103778637A (en) * | 2014-01-24 | 2014-05-07 | 北京空间机电研究所 | Side-slither radiometric calibration processing method based on histogram equalization |
| CN106871925A (en) * | 2016-12-29 | 2017-06-20 | 中国资源卫星应用中心 | A kind of remote sensing satellite relative radiometric calibration processing method of in-orbit comprehensive dynamic adjustment |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110782429A (en) * | 2019-09-26 | 2020-02-11 | 长沙天仪空间科技研究院有限公司 | Imaging quality evaluation method based on satellite-borne remote sensing camera |
| CN114708211B (en) * | 2019-09-26 | 2024-06-11 | 长沙天仪空间科技研究院有限公司 | Optimization processing method and system for satellite remote sensing imaging |
| CN110782429B (en) * | 2019-09-26 | 2022-04-15 | 长沙天仪空间科技研究院有限公司 | Imaging quality evaluation method based on satellite-borne remote sensing camera |
| CN114708211A (en) * | 2019-09-26 | 2022-07-05 | 长沙天仪空间科技研究院有限公司 | Optimization processing method and system for satellite remote sensing imaging |
| CN111257238B (en) * | 2020-02-05 | 2023-03-14 | 国家卫星气象中心(国家空间天气监测预警中心) | Detection element relative calibration method based on satellite-borne solar diffusion plate |
| CN111257238A (en) * | 2020-02-05 | 2020-06-09 | 国家卫星气象中心(国家空间天气监测预警中心) | Detection element relative calibration method based on satellite-borne solar diffusion plate |
| CN111415392A (en) * | 2020-03-13 | 2020-07-14 | 湖北师范大学 | Video satellite in-orbit relative radiation calibration method adopting Bell template push-broom mode imaging |
| CN111415392B (en) * | 2020-03-13 | 2023-04-21 | 湖北师范大学 | Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging |
| CN113177512A (en) * | 2021-05-20 | 2021-07-27 | 国家卫星气象中心(国家空间天气监测预警中心) | Matching threshold analysis method for intersatellite cross radiometric calibration |
| CN113177512B (en) * | 2021-05-20 | 2024-04-05 | 国家卫星气象中心(国家空间天气监测预警中心) | Matching threshold analysis method for cross radiation calibration between stars |
| CN113570523A (en) * | 2021-07-31 | 2021-10-29 | 中国人民解放军61646部队 | Method for automatic generation and extrapolation of relative radiation correction coefficients for optical images |
| CN113570523B (en) * | 2021-07-31 | 2022-01-28 | 中国人民解放军61646部队 | Method for automatic generation and extrapolation of relative radiation correction coefficients for optical images |
| CN113984673A (en) * | 2021-10-25 | 2022-01-28 | 中国科学院合肥物质科学研究院 | Satellite-borne polarization scanning remote sensor data pre-correction method and device |
| CN114125431A (en) * | 2021-11-22 | 2022-03-01 | 北京市遥感信息研究所 | Non-uniformity calibration correction method for static track optical large-area array camera |
| CN114125431B (en) * | 2021-11-22 | 2023-06-23 | 北京市遥感信息研究所 | Non-uniformity calibration correction method for stationary track optical large area array camera |
| CN116539542A (en) * | 2023-06-20 | 2023-08-04 | 北京朗森基科技发展有限公司 | Ball-borne imaging Fourier transform spectrometer for atmospheric trace gas analysis |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109712089A (en) | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite | |
| Suzuki et al. | Initial inflight calibration for Hayabusa2 optical navigation camera (ONC) for science observations of asteroid Ryugu | |
| CN102279393B (en) | Cross radiometric calibration method of hyper-spectral sensor based on multi-spectral sensor | |
| Brunner et al. | Toward more precise photometric redshifts: Calibration via CCD photometry | |
| CN102901516A (en) | Multispectral image radiation correction method based on absolute radiometric calibration | |
| CN112131746B (en) | Chlorophyll a concentration inversion method and system | |
| CN111415309B (en) | High-resolution remote sensing image atmosphere correction method based on minimum reflectivity method | |
| Yang et al. | A total sky cloud detection method using real clear sky background | |
| CN103761704B (en) | Image generating methods based on infrared remote sensing data and system | |
| CN103778637B (en) | A kind of in-orbit 90 ° of radiation calibration processing methods based on histogram balance | |
| Ren et al. | A method and results of color calibration for the Chang'e-3 terrain camera and panoramic camera | |
| Omari et al. | Retrieval of forest canopy parameters by inversion of the PROFLAIR leaf-canopy reflectance model using the LUT approach | |
| CN106940219A (en) | A kind of spectral response acquisition methods of broadband satellite remote sensor in orbit | |
| Brunner et al. | A probabilistic quantification of galaxy cluster membership | |
| Markelin | Radiometric calibration, validation and correction of multispectral photogrammetric imagery | |
| CN103808413A (en) | Method and device for temperature-emissivity separation based on noise separation | |
| CN104111080B (en) | Modulation transfer function (MTF) in-orbit measuring method for gazing satellite area array CCD camera | |
| CN109029735B (en) | Earth surface emissivity calculation method and device | |
| CN115015147A (en) | High-spatial-resolution high-spectral thermal infrared remote sensing image simulation method | |
| CN102855625A (en) | Method of obtaining rice canopy cluster indexes by using hemispherical photography method | |
| CN113570523A (en) | Method for automatic generation and extrapolation of relative radiation correction coefficients for optical images | |
| CN108154478A (en) | A kind of remote sensing image processing method | |
| Kubik et al. | Amethist: A method for equalization thanks to histograms | |
| Lobo et al. | A photometric catalogue of the Coma cluster core | |
| Ostrikov et al. | Calibration of hyperspectral data aviation mode according with accompanying ground-based measurements of standard surfaces of observed scenes |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190503 |