CN111415392B - Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging - Google Patents
Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging Download PDFInfo
- Publication number
- CN111415392B CN111415392B CN202010176900.7A CN202010176900A CN111415392B CN 111415392 B CN111415392 B CN 111415392B CN 202010176900 A CN202010176900 A CN 202010176900A CN 111415392 B CN111415392 B CN 111415392B
- Authority
- CN
- China
- Prior art keywords
- bell
- template
- satellite
- relative radiation
- imaging
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Abstract
The invention discloses a video satellite in-orbit relative radiation calibration method imaged by a Bell template push-broom mode, which comprises the following steps: s1, reconstructing a Bell template; s2, imaging original data based on a video satellite push-broom mode; s3, determining the relative radiation calibration reference standard of each wave band of the Bell template; s4, mapping the gray distribution of each probe element of the Bell template to the relative radiation calibration reference standard of each wave band. S5, finishing the relative radiation calibration work of the video satellite under the Bell template load push-broom imaging mode through the steps S1-S4, acquiring the relative radiation correction parameters of the video for each wave band of the Bell template, and providing input for the radiation correction of the follow-up Bell template original image. S6, application of the Bell template relative to the radiation correction parameters. The method is suitable for in-orbit relative radiation calibration of the video satellite imaged in the Bell template push-broom mode, and has the advantages of high frequency, high precision, convenience and low cost.
Description
Technical Field
The invention relates to an in-orbit relative radiation calibration method for a video satellite carrying a Bell template type load, in particular to an in-orbit relative radiation calibration method for a video satellite imaged in a Bell template push-broom mode.
Background
Under the promotion of the important special construction of national science and technology such as 'high resolution earth observation system', the space resolution, time resolution and data quality of the remote sensing satellite in China are greatly improved through innovations in the aspects of platform sensor development, multi-star networking, ground data processing and the like, and information service and decision support are provided for important fields such as modern agriculture in China, disaster prevention and reduction, resource environment, public safety and the like. With the penetration of various remote sensing applications, application demands develop from regular static census to real-time dynamic monitoring, and the continuous monitoring of global hot spot areas and targets by satellites to acquire dynamic information has become an urgent requirement. The video satellite can obtain time sequence images of targets within a certain time range, has the capability of continuously monitoring moving targets, and becomes a hot spot developed in recent years.
A plurality of parallel linear array sensors covering different spectral ranges are paved on a focal plane through a conventional linear array push-broom type optical remote sensing satellite sensor, the same target is imaged instantaneously in sequence, information on different wave bands of the target is obtained, and a color image of the target is formed through wave band combination. For a video satellite carrying an area array load, if a similar method is adopted, a multichannel focal plane needs to be paved on a focal plane of a satellite sensor, when the area array size is 4096×3072, at least 3 area array detectors with the size of 4096×3072 need to be paved on the focal plane, and three area arrays are required to have the same light path condition. The sensor has the advantages of large volume and weight, high development cost, high production process difficulty, and large data size, and is unfavorable for the application of space video to earth observation. In order to reduce the size, weight and cost of the sensor, reduce the difficulty of load development and production process, and reduce the pressure of on-board solid storage and data transmission, the video satellite load generally adopts a combination mode of a detector and an optical filter. A color filter array (Color Filter Array, CFA) is arranged in front of the detector, and color information screening is completed through the color filter array, wherein the color filter array usually adopts a Bayer template (also called a Bayer template), for example, jilin-one video satellite in China adopts an RGGB type Bayer template (fig. 1 (a)) and a Zhuhai-one video satellite adopts a GBRG type Bayer template (fig. 1 (b)). When the Bayer template is imaged, each detection element position only acquires one color, so that the weight and the volume of the sensor are effectively reduced, two thirds of data volume is reduced, and the complete red, green and blue (RGB) color band information is reconstructed through a Bayer interpolation algorithm.
In order to solve the problems of single imaging mode and narrow surface coverage of the traditional video satellite, most of the domestic and foreign video satellites have multiple imaging modes, such as staring imaging (shown in figure 2) to obtain surface dynamic information, pushing scanning imaging (shown in figure 3) to expand the single imaging range of the video satellite and obtain long-strip and large-range imaging data, and the two modes can be flexibly adjusted according to actual on-orbit application requirements.
In-orbit relative radiation calibration is a key technology for guaranteeing the radiation quality of satellites, and is an indispensable ring in a ground processing system after satellites are in orbit. However, the satellite is affected by the transmission vibration, the physical environment such as warm and space environment after the satellite is in orbit, and the attenuation problem of time variation of the satellite sensor probe element, so that the satellite sensor in-orbit response state has periodic complex variation, the laboratory calibration result before the video satellite transmission cannot be applied in orbit for a long time, the quality of satellite imaging data is directly reduced, and in-orbit radiometric calibration is urgently needed to improve the image quality. There have been related studies on in-orbit relative radiometric calibration for video area array gaze patterns, but less for video area array push-broom patterns. And the on-orbit uniform field calibration is different from video area array imaging, and the on-orbit uniform field calibration is also different from the statistical calibration or yaw calibration of a linear array push-broom satellite, and the area array Bell module push-broom imaging mode brings challenges to the on-orbit radiation calibration of the video satellite. Therefore, the on-orbit radiation calibration technology of the video satellite area array sensor is researched, particularly the on-orbit radiation calibration technology of the video satellite imaged in the area array Bell template push-broom mode is adopted, the satellite radiation quality is improved, and the method has very important significance for guaranteeing the application effect of the video satellite in the dynamic observation field.
Disclosure of Invention
The invention provides a novel on-orbit relative radiation calibration method aiming at a video satellite imaged by a Bell module push-broom mode. The original Bell module data is converted into the linear array probe array by reconstructing the original Bell module data, and the video satellite on-orbit relative radiation calibration is realized based on the principle that the gray level distribution of each probe element image of the Bell template is consistent under the condition of massive samples.
The invention aims to provide an on-orbit relative radiation calibration method suitable for a video satellite push-broom mode imaging condition carrying area array Bell load, which specifically comprises the following steps:
Step 2, based on the original data imaged by the video satellite push-broom mode, counting the gray distribution P of each probe element of the reconstructed Bell area array j (k) Where j is E [0, N total ]，k∈[0,2 n -1]N is the number of quantization bits of the video satellite;
step 3, determining the relative radiation calibration reference P of each wave band of the Bell template G (k)、P B (k)、P R (k) Wherein R, G and B respectively represent three wave bands of red, green and blue;
step 4, mapping the gray distribution of each probe element of the Bell template to the relative radiation calibration reference standard of each wave band, wherein the obtained mapping relation is the relative radiation correction parameter;
and 5, completing the relative radiation calibration work under the video satellite Bell template load push-broom imaging mode through the steps 1-4, and obtaining the relative radiation correction parameters of each wave band of the Bell template.
Further, the imaging probe number N in the step 1 total Twice the number M of single-line detectors of the bell template, namely:
N total ＝2×M (1)。
further, the specific implementation manner of the step 2 is as follows,
2.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i ；
2.2, determining the initial position of the Bell module, wherein the initial position judgment rule is required to be confirmed according to a satellite-ground interface file provided by a satellite development party;
if the initial position of the ith behavior bell template is changed to step 2.3, otherwise, executing i=i+1, and changing to step 2.1;
2.3, judging whether the ith row of Bell template data is valid or not, wherein the judging rule is as follows:
|LineCounter i -LineCounter i+1 |＝1 (2)
|LineCounter i-1 -LineCounter i |＝1 (3)
if the ith row of the Bayer template data is valid, turning to the step 2.4, otherwise, discarding the ith row of the data, and turning to the step 2.1 to process the next row of the data;
2.4, reading the ith row of the Bayer template image data, performing the following statistics on the imaging pixel number and the total pixel number of the jth probe element,
PixelNums j (k)＝PixelNums j (k)+1 (4)
TotalPixelNums j ＝TotalPixelNums j +1 (5)
PixelNums in j (k) Imaging the number of pixels with the gray value equal to k of the jth probe element after reconstructing the planar array Bayer template, and defaulting to 0 initially; totalPixelNums j The total number of imaging pixels of the j-th probe element after the planar array Bell template is reconstructed is initially defaulted to 0;
2.5, repeating the step 2.4 until all the probe elements of the Bayer template are processed after the reconstruction is completed, and converting the step 2.1 to process the next row of data;
2.6, repeating the steps 2.1-2.5 until all lines of the original data imaged in the video satellite push-broom mode are processed, and turning to the step 2.7;
2.7, calculating the gray distribution P of the j-th probe element of the reconstructed Bayer planar array based on the following formula j (k)，
And 2.8, repeating the step 2.7 to obtain the gray distribution of all the probe elements of the reconstructed Bell area array, and ending.
Further, the line count calculation formula of the OVS-3 satellite data is as follows,
LineCounter i ＝LCByte[0]×16777216+LCByte[1]×65536LCByte[2]×256+LCByte[3] (7)
where LCByte is the binary data stored representing the row count.
Further, the judgment rule of the OVS-3 satellite is as follows:
(LineCounter i ％2)＝0 (8)
the above equation represents a remainder for 2, requiring the remainder of linccount divided by 2 to be 0.
Further, in step 3, P G (k)、P B (k)、P R (k) The calculation formula of (2) is as follows:
wherein R, G, B respectively represent three wave bands of red, green and blue.
Further, in the step 4, the mapping rule of mapping the gray distribution of each probe element of the bell template to the radiometric calibration reference standard of each wave band is as follows,
4.1, mapping K to K according to the following rule for the j-th probe element and the K-th gray scale;
4.1.1, determining the gray levels x and y, where x, y E [0,2 n -1]N is the quantization bit number of the video satellite, so that the gray level distribution P of the kth gray level of the jth probe element j (k) Satisfies the following formula:
P R (k-x)≤P j (k)≤P R (k+y) (12)
4.1.2, determining the K value after the K-level gray scale mapping according to the following formula,
4.2, traversing all gray levels k in sequence to obtain a mapping table Histo of the jth probe element in the video satellite quantization range j I.e. inputting a K-level gray scale to obtain a new gray scale K, as shown in the following formula:
K＝Histo j (k) (14)
the mapping table Histo is the j-th probe element relative radiation correction parameter of the video satellite area array Bell template;
4.3, repeating the steps 4.1-4.2 to obtain the relative radiation correction parameters of all the probe elements of the R wave band of the video satellite area array Bayer template;
and 4.4, repeating the steps 4.1-4.3 to obtain the relative radiation correction parameters of all the probe elements of the three wave bands of the video satellite area array bell module R, G, B.
Further, the method also comprises a step 6 of applying the relative radiation correction parameters of the Bell templates, wherein the specific implementation mode is as follows,
6.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i ；
6.2, determining the starting position of the Bell module; the initial position judgment rule needs to be confirmed according to a satellite-ground interface file provided by a satellite development party;
if the initial position of the behavior Bell template is the initial position, turning to the step 6.3, otherwise turning to the step 6.1;
6.3, reading the original image data of the ith row and storing the original image data to DN j The imaging gray value of the j-th probe element of the Bayer template after reconstruction;
6.4, correcting parameter Histo by using relative radiation of each probe element of the reconstructed Bell template j Completing the relative radiation correction of all the probe element imaging gray values of the Bell template, and setting the imaging gray value DN of the jth Bell template probe element j The gray value corrected by the relative radiation is DN j corr The following steps are:
DN j corr ＝Histo j (DN j ) (15)
6.5, repeating the steps 6.3-6.4 to sequentially finish the relative radiation correction of the imaging gray level of all the probe elements of the Bell template;
and 6.6, repeating the steps 6.1-6.5 to finish the relative radiation correction of all the line image data of the original data.
Compared with the prior art, the invention has the following characteristics and beneficial effects:
(1) An on-orbit relative radiation calibration method for video Bell module push-broom mode imaging is provided.
(2) On-orbit targeting does not require imaging of the surface uniform field.
(3) The satellite is not required to have on-board calibration processing capability.
(4) The video Bayer template is directly calibrated, so that the randomness influence of different Bayer interpolation algorithms on the relative radiometric calibration accuracy is eliminated, and the accuracy of the relative radiometric calibration algorithm is improved.
(5) The method can be widely applied to the in-orbit relative radiation calibration of the video satellite imaged in the push-broom mode of the area array sensor.
Drawings
Fig. 1 is an optical video satellite bell template (a) "gilin number one" video star bell template (b) "zhuhai number one" video star bell template).
Fig. 2 is an optical video satellite gaze imaging schematic.
Fig. 3 is a schematic illustration of optical video satellite push-broom imaging.
Fig. 4 is a video satellite area array bell 2-level integral push broom schematic (GBRG bell template).
Fig. 5 is a schematic view of reconstruction of a video satellite area array bell module (GBRG bell template).
Detailed Description
The technical scheme of the invention is described in detail below with reference to the accompanying drawings.
The invention provides a video satellite in-orbit relative radiation calibration method imaged by a Bell template push-broom mode, which specifically comprises the following steps:
s1, reconstructing a Bell template, wherein the method is concretely realized as follows;
the push-broom mode of the video satellite with the planar array Bell load is different from the push-broom mode of the conventional linear array satellite. The push-broom mode of the planar array bell load video satellite is to realize the function of multi-level integration similar to the linear array satellite load by selecting a plurality of rows of bell templates in the middle of the planar array bell load, as shown in fig. 4, which is a push-broom schematic of 2-level integration of the video star GBRG bell template.
When the video satellite performs on-orbit push-broom imaging in the form of an area array Bayer template, the minimum imaging unit GBRG of the Bayer template is equivalent to 4 imaging probe elements of linear array push-broom load. Thus, scallop is providedThe minimum imaging unit GBRG of the template is converted into 4 imaging probe elements, namely each imaging detector row (two rows of planar array Bell loads) of the Bell template is converted into one row (conversion in the Bell first-order integral mode is shown in fig. 5), and the number N of the imaging probe elements after conversion is calculated total Twice the number M of single-line detectors of the bell template, namely:
N total ＝2×M (16)
s2, based on original data imaged by a video satellite push-broom mode, counting gray distribution P of each probe element of the reconstructed Bell area array j (k) Where j is E [0, N total ]，k∈[0,2 n -1]N is the number of quantization bits for the video satellite. The method further comprises the following steps:
2.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i . Different satellite data have different calculation modes and are determined according to satellite-ground interface files provided by satellite development parties. The following formula is the line count calculation formula of the pearl sea No. 03 group video satellite OVS-3 satellite data.
LineCounter i ＝LCByte[0]×16777216+LCByte[1]×65536LCByte[2]×256+LCByte[3] (17)
Where LCByte is the binary data stored representing the row count.
2.2, determining the starting position of the Bell module. The starting position judgment rule needs to be confirmed according to the satellite-ground interface file provided by the satellite development party. The judgment rule of the pearl sea No. 03 group video satellite OVS-3 satellite is as follows:
(LineCounter i ％2)＝0 (18)
the above equation represents a remainder for 2, requiring the remainder of linccount divided by 2 to be 0; if the starting position of the ith behavior bell template is changed to step 2.3, otherwise, i=i+1 is executed, and step 2.1 is changed.
And 2.3, judging whether the ith row of the Bell template data is valid. The judgment rule is as follows:
|LineCounter i -LineCounter i+1 |＝1 (19)
|LineCounter i-1 -LineCounter i |＝1 (20)
if the ith row of the Bayer template data is valid, turning to the step 2.4, otherwise, discarding the ith row of the data, and turning to the step 2.1 to process the next row of the data.
And 2.4, reading the ith row of the Bayer template image data, and executing the following formula to count the imaging pixel number and the total pixel number of the jth probe element.
PixelNums j (k)＝PixelNums j (k)+1 (21)
TotalPixelNums j ＝TotalPixelNums j +1 (22)
PixelNums in j (k) Imaging the number of pixels with the gray value equal to k of the jth probe element after reconstructing the planar array Bayer template, and defaulting to 0 initially; totalPixelNums j And (3) imaging the total number of pixels for the j-th probe element after reconstructing the planar array Bell template, and initially defaulting to 0.
And 2.5, repeating the step 2.4 until all the probe elements of the Bayer template are processed after the reconstruction is completed, and converting the step 2.1 to process the next row of data.
And 2.6, repeating the steps 2.1-2.5 until all lines of the original data imaged by the video satellite push-broom mode are processed, and turning to the step 2.7. The number of the original data counted here takes the push-broom imaging orbit number of the video satellite as a unit, and at least 5 orbits of data are needed to participate in counting under the condition of the same integral number.
2.7, calculating the gray distribution P of the j-th probe element of the reconstructed Bayer planar array based on the following formula j (k)，
And 2.8, repeating the step 2.7 to obtain the gray distribution of all the probe elements of the reconstructed Bell area array, and ending.
S3, determining that the relative radiometric calibration reference standard of the three wave bands of the Bell template R, G, B is P respectively G (k)、P B (k)、P R (k) The calculation formula is as follows:
wherein R, G, B respectively represent three wave bands of red, green and blue;
s4, mapping the gray distribution of each probe element of the Bell template to the relative radiometric calibration reference standard of each wave band, wherein the mapping rule is as follows (taking the Bell template R wave band as an example):
4.1, for the j-th probe element, the K-th gray scale, K is mapped to K according to the following rule.
4.1.1, the gray levels x and y (x, y. Epsilon. [0,2 ] n -1]N is the number of quantization bits of the video satellite) such that the gray distribution P of the kth gray level of the jth probe element j (k) Satisfies the following formula:
P R (k-x)≤P j (k)≤P R (k+y) (27)
4.1.2, determining the K value after the K-level gray scale mapping according to the following formula,
4.2, traversing all gray levels k in sequence to obtain a mapping table Histo of the jth probe element in the video satellite quantization range j I.e. inputting a K-level gray scale to obtain a new gray scale K, as shown in the following formula:
K＝Histo j (k) (29)
the mapping table Histo is the j-th probe element relative radiation correction parameter of the video satellite area array Bell template.
And 4.3, repeating the steps 4.1-4.2 to obtain the relative radiation correction parameters of all the probe elements of the R wave band of the video satellite area array Bayer template.
And 4.4, repeating the steps 4.1-4.3 to obtain the relative radiation correction parameters of all the probe elements of the three wave bands of the video satellite area array bell module R, G, B.
S5, finishing the relative radiation calibration work of the video satellite under the Bayer template load push-broom imaging mode through the steps S1-S4, acquiring the relative radiation correction parameters of each wave band of the Bayer template, and providing input for the radiation correction of the follow-up Bayer template original image.
S6, application of the Bell template relative radiation correction parameters further comprises:
6.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i . Different satellite data have different calculation modes and are determined according to satellite-ground interface files provided by satellite development parties. The following formula is a line count calculation formula for OVS-3 satellite data.
LineCounter i ＝LCByte[0]×16777216+LCByte[1]×65536LCByte[2]×256+LCByte[3] (30)
Where LCByte is the binary data stored representing the row count.
And 6.2, determining the starting position of the Bell module. The starting position judgment rule needs to be confirmed according to the satellite-ground interface file provided by the satellite development party. The judgment rule of the OVS-3 satellite is as follows:
(LineCounter i ％2)＝0 (31)
if the starting position of the behavior bell template is changed to step 6.3, otherwise, the starting position of the behavior bell template is changed to step 6.1.
6.3, reading the original image data of the ith row and storing the original image data to DN j And reconstructing an imaging gray value of the j-th probe element of the back Bayer template.
6.4, correcting parameter Histo by using relative radiation of each probe element of the reconstructed Bell template j And (5) finishing the relative radiation correction of all the imaging gray values of the probe elements of the Bell template. Setting the imaging gray level DN of the jth Bell template probe element j The gray value corrected by the relative radiation is DN j corr The following steps are:
DN j corr ＝Histo j (DN j ) (32)
and 6.5, repeating the steps 6.3-6.4 to sequentially finish the relative radiation correction of the imaging gray level of all the probe elements of the Bell template.
And 6.6, repeating the steps 6.1-6.5 to finish the relative radiation correction of all the line image data of the original data.
The specific embodiments described herein are offered by way of example only to illustrate the spirit of the invention. Those skilled in the art may make various modifications or additions to the described embodiments or substitutions thereof without departing from the spirit of the invention or exceeding the scope of the invention as defined in the accompanying claims.
Claims (7)
1. The video satellite on-orbit relative radiation calibration method imaged in a Bell template push-broom mode is characterized by comprising the following steps of:
step 1, reconstructing a Bell template aiming at a Bell load video satellite to obtain the number N of imaging probe elements total ；
The mode of reconstructing the bell template is as follows: converting a minimum imaging unit GBRG of the Bell template into 4 imaging probe elements, namely converting each imaging detector row of the Bell template into one row;
step 2, based on the original data imaged by the video satellite push-broom mode, counting the gray distribution P of each probe element of the reconstructed Bell area array j (k) Where j is E [0, N total ]，k∈[0,2 n -1]N is the number of quantization bits of the video satellite;
the specific implementation of step 2 is as follows,
2.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i ；
2.2, determining the initial position of the Bell module, wherein the initial position judgment rule is required to be confirmed according to a satellite-ground interface file provided by a satellite development party;
if the initial position of the ith behavior bell template is changed to step 2.3, otherwise, executing i=i+1, and changing to step 2.1;
2.3, judging whether the ith row of Bell template data is valid or not, wherein the judging rule is as follows:
|LineCounter i -LineCounter i+1 |＝1 (1)
|LineCounter i-1 -LineCounter i |＝1 (2)
if the ith row of the Bayer template data is valid, turning to the step 2.4, otherwise, discarding the ith row of the data, and turning to the step 2.1 to process the next row of the data;
2.4, reading the ith row of the Bayer template image data, performing the following statistics on the imaging pixel number and the total pixel number of the jth probe element,
PixelNums j (k)＝PixelNums j (k)+1 (3)
TotalPixelNums j ＝TotalPixelNums j +1 (4)
PixelNums in j (k) Imaging the number of pixels with the gray value equal to k of the jth probe element after reconstructing the planar array Bayer template, and defaulting to 0 initially; totalPixelNums j The total number of imaging pixels of the j-th probe element after the planar array Bell template is reconstructed is initially defaulted to 0;
2.5, repeating the step 2.4 until all the probe elements of the Bayer template are processed after the reconstruction is completed, and converting the step 2.1 to process the next row of data;
2.6, repeating the steps 2.1-2.5 until all lines of the original data imaged in the video satellite push-broom mode are processed, and turning to the step 2.7;
2.7, calculating the gray distribution P of the j-th probe element of the reconstructed Bayer planar array based on the following formula j (k)，
2.8, repeating the step 2.7 to obtain the gray distribution of all the probe elements of the reconstructed Bell area array, and ending;
step 3, determining the relative radiation calibration reference P of each wave band of the Bell template G (k)、P B (k)、P R (k) Wherein R, G and B respectively represent three wave bands of red, green and blue;
step 4, mapping the gray distribution of each probe element of the Bell template to the relative radiation calibration reference standard of each wave band, wherein the obtained mapping relation is the relative radiation correction parameter;
and 5, completing the relative radiation calibration work under the video satellite Bell template load push-broom imaging mode through the steps 1-4, and obtaining the relative radiation correction parameters of each wave band of the Bell template.
2. The method for calibrating the in-orbit relative radiation of a video satellite imaged in a bayer pattern push-broom mode according to claim 1, wherein the method comprises the following steps: number of imaging probe elements N in step 1 total Twice the number M of single-line detectors of the bell template, namely:
N total ＝2×M (6)。
3. the method for calibrating the in-orbit relative radiation of a video satellite imaged in a bayer pattern push-broom mode according to claim 1, wherein the method comprises the following steps: the line count calculation formula of the pearl sea No. 03 group video satellite OVS-3 satellite data is as follows,
LineCounter i ＝LCByte[0]×16777216+LCByte[1]×65536LCByte[2]×256+LCByte[3] (7)
where LCByte is a binary number set storing a count representing a row.
4. The method for calibrating the in-orbit relative radiation of a video satellite imaged in a bayer pattern push-broom mode according to claim 1, wherein the method comprises the following steps: the judgment rule of the pearl sea No. 03 group video satellite OVS-3 satellite is as follows:
(LineCounter i ％2)＝0 (8)
the above equation represents a remainder for 2, requiring the remainder of linccount divided by 2 to be 0.
5. The method for calibrating the in-orbit relative radiation of a video satellite imaged in a bayer pattern push-broom mode according to claim 1, wherein the method comprises the following steps: p in step 3 G (k)、P B (k)、P R (k) The calculation formula of (2) is as follows:
wherein R, G, B respectively represent three wave bands of red, green and blue.
6. The method for calibrating the in-orbit relative radiation of a video satellite imaged in a bayer pattern push-broom mode according to claim 5, wherein the method comprises the following steps: in step 4, the mapping rule of mapping the gray distribution of each probe element of the bell template to the radiometric calibration reference standard of each wave band is as follows,
4.1, mapping K to K according to the following rule for the j-th probe element and the K-th gray scale;
4.1.1, determining the gray levels x and y, where x, y E [0,2 n -1]N is the quantization bit number of the video satellite, so that the gray level distribution P of the kth gray level of the jth probe element j (k) Satisfies the following formula:
P R (k-x)≤P j (k)≤P R (k+y) (12)
4.1.2, determining the K value after the K-level gray scale mapping according to the following formula,
4.2, traversing all gray levels k in sequence to obtain a mapping table Histo of the jth probe element in the video satellite quantization range j I.e. inputting a K-level gray scale to obtain a new gray scale K, as shown in the following formula:
K＝Histo j (k) (14)
the mapping table Histo is the j-th probe element relative radiation correction parameter of the video satellite area array Bell template;
4.3, repeating the steps 4.1-4.2 to obtain the relative radiation correction parameters of all the probe elements of the R wave band of the video satellite area array Bayer template;
and 4.4, repeating the steps 4.1-4.3 to obtain the relative radiation correction parameters of all the probe elements of the three wave bands of the video satellite area array bell module R, G, B.
7. The method for in-orbit relative radiometric calibration of a video satellite imaged in a bayer template push-broom mode according to claim 6, wherein the method comprises the steps of: the method also comprises a step 6 of applying the relative radiation correction parameters of the Bell templates, and the specific implementation mode is as follows,
6.1, calculating the line count LineCountr of the ith line of the original data and the Bayer template data i ；
6.2, determining the starting position of the Bell module; the initial position judgment rule needs to be confirmed according to a satellite-ground interface file provided by a satellite development party;
if the initial position of the behavior Bell template is the initial position, turning to the step 6.3, otherwise turning to the step 6.1;
6.3, reading the original image data of the ith row and storing the original image data to DN j The imaging gray value of the j-th probe element of the Bayer template after reconstruction;
6.4, correcting parameter Histo by using relative radiation of each probe element of the reconstructed Bell template j Completing the relative radiation correction of all the probe element imaging gray values of the Bell template, and setting the imaging gray value DN of the jth Bell template probe element j The gray value corrected by the relative radiation is DN j corr The following steps are:
DN j corr ＝Histo j (DN j ) (15)
6.5, repeating the steps 6.3-6.4 to sequentially finish the relative radiation correction of the imaging gray level of all the probe elements of the Bell template;
and 6.6, repeating the steps 6.1-6.5 to finish the relative radiation correction of all the line image data of the original data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010176900.7A CN111415392B (en) | 2020-03-13 | 2020-03-13 | Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010176900.7A CN111415392B (en) | 2020-03-13 | 2020-03-13 | Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111415392A CN111415392A (en) | 2020-07-14 |
| CN111415392B true CN111415392B (en) | 2023-04-21 |
Family
ID=71493038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010176900.7A Active CN111415392B (en) | 2020-03-13 | 2020-03-13 | Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111415392B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107036629A (en) * | 2017-04-20 | 2017-08-11 | 武汉大学 | The in-orbit relative radiometric calibration method and system of video satellite |
| CN107465900A (en) * | 2017-08-22 | 2017-12-12 | 长光卫星技术有限公司 | A kind of preprocess method suitable for large area array sub-meter grade remote sensing video image unit |
| CN108174127A (en) * | 2018-01-30 | 2018-06-15 | 中国科学院长春光学精密机械与物理研究所 | Relative radiometric correction methods of the face battle array CMOS under global shutter working method |
| CN109521415A (en) * | 2018-12-19 | 2019-03-26 | 上海同繁勘测工程科技有限公司 | Radiant correction apparatus and system |
| CN109712089A (en) * | 2018-12-14 | 2019-05-03 | 航天恒星科技有限公司 | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite |
| CN110501025A (en) * | 2019-07-26 | 2019-11-26 | 中国科学院长春光学精密机械与物理研究所 | A kind of In-flight calibration method of Rolling shutter face battle array CMOS camera |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10598550B2 (en) * | 2017-02-17 | 2020-03-24 | Flir Systems, Inc. | Radiometric correction and alignment techniques for thermal imager with non-contact temperature sensor |
| US10692178B2 (en) * | 2018-02-22 | 2020-06-23 | Raytheon Company | Spatial radiometric correction of an optical system having a color filter mosaic |
-
2020
- 2020-03-13 CN CN202010176900.7A patent/CN111415392B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107036629A (en) * | 2017-04-20 | 2017-08-11 | 武汉大学 | The in-orbit relative radiometric calibration method and system of video satellite |
| CN107465900A (en) * | 2017-08-22 | 2017-12-12 | 长光卫星技术有限公司 | A kind of preprocess method suitable for large area array sub-meter grade remote sensing video image unit |
| CN108174127A (en) * | 2018-01-30 | 2018-06-15 | 中国科学院长春光学精密机械与物理研究所 | Relative radiometric correction methods of the face battle array CMOS under global shutter working method |
| CN109712089A (en) * | 2018-12-14 | 2019-05-03 | 航天恒星科技有限公司 | Method suitable for the infrared shortwave load relative detector calibration of export-oriented remote sensing satellite |
| CN109521415A (en) * | 2018-12-19 | 2019-03-26 | 上海同繁勘测工程科技有限公司 | Radiant correction apparatus and system |
| CN110501025A (en) * | 2019-07-26 | 2019-11-26 | 中国科学院长春光学精密机械与物理研究所 | A kind of In-flight calibration method of Rolling shutter face battle array CMOS camera |
Non-Patent Citations (4)
| Title |
|---|
| A.s. Arya 等.Mars Colour Camera: The payload characterization/calibration and data analysis from Earth imaging phase. 《Current Science》.2015,第109卷(第6期),全文. * |
| Guo Zhang 等.On-orbit relative radiometric calibration of optical video satellites without uniform calibration sites.《International Journal of Remote Sensing》.2019,第40卷(第4期),全文. * |
| 张过 等.遥感25号无场化相对辐射定标.《测绘学报》.2017,第46卷(第8期),全文. * |
| 王东升 等.基于梯度局部一致性的自适应Bayer模式彩色图像恢复算法.《南华大学学报(自然科学版)》.2019,第33卷(第2期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111415392A (en) | 2020-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8089533B2 (en) | Fixed pattern noise removal circuit, fixed pattern noise removal method, program, and image pickup apparatus | |
| US9467628B2 (en) | High dynamic range image sensor | |
| CN102300057B (en) | Method for correcting response inconsistency of linear array CCD (Charge Coupled Device) image elements | |
| US20070019085A1 (en) | Method and apparatus providing imager noise reduction using parallel input arithmetic mean modules | |
| CN100432697C (en) | Linear array scan image non-uniform correction method based on statistical restraint model | |
| CN109459135A (en) | A kind of CCD imaging spectrometer image bearing calibration | |
| CN112529807B (en) | Relative radiation correction method and device for satellite image | |
| CN110009688A (en) | A kind of infrared remote sensing image relative radiometric calibration method, system and remote sensing platform | |
| US7796153B1 (en) | Equalization system and method for an imaging sensor | |
| CN111415392B (en) | Video satellite in-orbit relative radiation calibration method based on Bell template push-broom imaging | |
| CN107220947B (en) | Remote sensing image relative radiation correction method and system | |
| CN114125431A (en) | Non-uniformity calibration correction method for static track optical large-area array camera | |
| CN112085684A (en) | Method and device for fusing remote sensing images | |
| CN112291446A (en) | Non-uniformity correction method for large-area array CMOS image sensor | |
| GB2584743A (en) | Method and device for processing an image signal of an image sensor for a vehicle | |
| US20220155454A1 (en) | Analysis portion, time-of-flight imaging device and method | |
| EP2517172B1 (en) | Filter setup learning for binary sensor | |
| CN111751006B (en) | Data calibration method, target radiation characteristic analysis method and device | |
| US6324308B1 (en) | Non-uniformity correction method and apparatus for imaging systems | |
| KR20100054929A (en) | Satelite image correction method and apparatus using neural network | |
| CN112070663B (en) | Multi-detector spliced camera radiation calibration combined data processing method | |
| CN113347376B (en) | Method for compensating adjacent pixel crosstalk of image sensor | |
| CN115830146B (en) | On-orbit relative radiation calibration and correction method for aerospace optical remote sensing camera | |
| US7495806B2 (en) | System and method for compensating for noise in a captured image | |
| Guo et al. | Study on the relative radiometric correction of CBERS satellite CCD image |
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 |