KR20150074827A - System and Method for Preprocessing Satellite Image - Google Patents

Abstract

The present invention relates to a system and a method for pretreating a satellite image. The system according to the present invention comprises: a bolometeic correction unit for receiving data observed by a meteorological payload of a satellite, and performing basic bolometeic correction of a UV channel; and a geometric correction unit for performing geometric correction by using data in which basic bolometeic correction of the UV channel is performed. According to the present invention, data observed by a meteorological payload of a satellite is directly received to an image position correction device (INRSM), and basic bolometeic correction of a UV channel is performed. Then, geometric correction is performed, thereby reducing starting time of geometric correction.

Description

TECHNICAL FIELD The present invention relates to a satellite image preprocessing system,

The present invention relates to a satellite image preprocessing system and method, and more particularly, to a satellite image preprocessing system and method capable of receiving data observed by a satellite meteorological payload, performing radial correction, .

In general, satellites are classified into low orbit orbit, satellite orbiting satellite and geostationary orbit satellite depending on the orbit around the Earth. Of these, geostationary satellites are located at 36,000 km above the equator, You can observe the weather on the earth by revolving around the earth. Geostationary satellites can always observe the same place on the earth, and geostationary satellite reception antennas can always be fixed in the same direction. Thus, geostationary satellites can be used to monitor and predict the change of weather by continuously observing cloud movement, weather conditions, and so on.

The meteorological satellite system can perform meteorological observations that monitor, observe, and analyze environmental changes such as temperature, wind, and humidity based on images observed from geostationary satellites. At this time, in order to utilize the image observed from the geostationary satellite, a preprocessing process of radial and geometric calibration is required for the image, and a high rate information transmission (HRIT) To the user within the required time.

Conventionally, a full radiation correction is performed on an image observed from a geosynchronous satellite, and a complete radiation correction result is input to a geometric correction process, and the result of complete radiation correction is transmitted to a geometric correction device in a file form.

However, since the geometric correction can not be started until the complete radiation correction is completed, there is a problem that it takes a long time to generate the radiation correction result because the calculation formula for the perfect radiation correction is complicated.

As a result, since the start time for the geometric correction processing is delayed, there is a problem that the HRIT (High Rate Information Transmission) required time is not satisfied.

KR 1995-0013152 A

SUMMARY OF THE INVENTION Accordingly, an object of the present invention is to provide a method and apparatus for correcting a radiation correction of an infrared channel by directly receiving data observed by a satellite meteorological payload device in an image position correcting device (INRSM) To a satellite image preprocessing system and method.

Further, the present invention includes other objects that can be achieved from the construction of the present invention described later, in addition to the objects explicitly mentioned.

According to another aspect of the present invention, there is provided a satellite image preprocessing system including a basic copy correcting unit for receiving data observed in a satellite payload and performing basic copy correction of an infrared channel, And a geometric correction unit that performs geometric correction using the data on which the copy correction is performed.

And a complete radiopaque correction unit for receiving the data observed in the satellite payload and performing a full radiation correction of the infrared channel.

The basic copy correction unit performs basic copy correction of the infrared channel using Equation (1)

[Equation 1]

은 매 흑체관측마다 계산되는 기울기(nominal slope)로,

,

는 흑체관측 온도로부터 계산한 결과로,

,

는 위성 발사 전 결정된 계수이고,

는 흑체관측에서 측정하는 온도,

는 흑체관측값의 평균,

는 흑체관측 시 우주보기관측값의 평균,

는 지구관측시 우주보기관측값의 평균,

일 수 있다.Where R is the base radiation correction result, q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is a nominal slope calculated for every blackbody observation,

,

As a result of calculation from the blackbody observation temperature,

,

Is the determined coefficient before satellite launch,

The temperature measured at the blackbody observation,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observation,

Lt; / RTI >

Meanwhile, the satellite image preprocessing method according to an embodiment of the present invention includes the steps of receiving data observed in a satellite payload and performing basic radiance correction of the infrared channel, and using the data subjected to the basic radial correction of the infrared channel And performing a geometric correction on the image data.

Receiving the observed data from the satellite payload and performing a full radiation correction of the infrared channel.

The basic copy correction step performs basic copy correction of the infrared channel using Equation (1)

[Equation 1]

은 매 흑체관측마다 계산되는 기울기(nominal slope)로,

,

는 흑체관측 온도로부터 계산한 결과로,

,

는 위성 발사 전 결정된 계수이고,

는 흑체관측에서 측정하는 온도,

는 흑체관측값의 평균,

는 흑체관측 시 우주보기관측값의 평균,

는 지구관측시 우주보기관측값의 평균,

일 수 있다.Where R is the base radiation correction result, q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is a nominal slope calculated for every blackbody observation,

,

As a result of calculation from the blackbody observation temperature,

,

Is the determined coefficient before satellite launch,

The temperature measured at the blackbody observation,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observation,

Lt; / RTI >

A computer-readable medium according to another embodiment of the present invention records a program for causing a computer to execute any one of the above methods.

As described above, according to the satellite image preprocessing system and method according to the embodiment of the present invention, the data observed by the satellite vapor feeder is directly received by the video position correcting device (INRSM) to perform basic copy correction of the infrared channel, There is an advantage that the start time of the geometric correction can be shortened by performing the correction.

In addition, since the calculation formula for the basic copy correction is simpler, the copy correction can be performed more quickly, and even if the basic copy correction result is input as an input, the accuracy of the geometric correction accuracy can be satisfied.

Accordingly, there is an advantage that the file that has undergone the geometric correction can be transmitted within the HRIT (High Rate Information Transmission) required time.

1 is a configuration diagram of a satellite image preprocessing system according to an embodiment of the present invention.
2 is a view showing a result of tilt calculation for the SWIR detector A. FIG.
FIG. 3 is a view showing a result of tilt calculation for the WIN2 detector B. FIG.
FIG. 4 is a table in which an MI performance evaluation tool analyzes geometric correction results obtained by inputting a radiation correction result processed with a perfect copy correction formula and a basic copy correction formula.
FIGS. 5 and 6 are diagrams showing the result of analyzing the geometric correction accuracy using the landmark, using the results of the perfect copy correction formula and the basic copy correction formula, respectively.
7 is a flowchart illustrating a preprocessing process of a satellite image according to an embodiment of the present invention.

Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those skilled in the art can easily carry out the present invention.

1 is a block diagram of a satellite image preprocessing system according to an embodiment of the present invention.

As shown in FIG. 1, the satellite image preprocessing system 100 includes an image position correcting apparatus 120 and an image copy correcting apparatus 140.

An image navigation and registration software module (INRSM) 120 includes a basic copy correction unit 122 and a geometry correction unit 124.

The basic copy correction unit 122 receives the raw data observed at the meteorological payload 50 (data at the initial stage) and performs basic copy correction of the infrared channel.

Before describing the basic radiation correction, the radiation correction will first be described. The radiation correction converts the data transmitted by the weather payload 50 into a quantitative physical quantity that can be used for actual research. (Digital Count or Digital Number) that is measured by the sensor module on the ground and sent to the ground is converted into a radiance value. The meteorological payload 50 includes a sensor module for performing actual meteorological observations, a power module and a sensor module which suitably convert the power transmitted to the satellites to the meteorological payload 50 to supply them stably, And an electronic module for exchanging electric power and various electrical signals between the satellite and the satellite. Among them, the sensor module includes a scan mirror that allows the user to observe the desired spot, an aperture that collects the light incident on the scan mirror, a spectrometer that classifies the incident radiation through the telescope into wavelengths, A black body for ensuring the performance of infrared radiation energy, and a cooling system for keeping the temperature of the detector at a cryogenic temperature.

Accordingly, the meteorological payload 50 is directed to the target observation region on the ground, and the scan mirror observes the region moving in the east / west direction and the south / north direction, and the radiant energy collected through the aperture passes through the main and sub- The beam is divided into wavelengths within the spectroscope and transmitted to the detector. The radiant energy entering the telescope can be detected in each detector by the channel.

In order to perform the radiation correction on the earth region, the blackbody observations and the cosmic view observations inside the meteorological payload 50 as well as the earth should be performed in parallel. The basic radiation correction formula applying the slope calculated for each blackbody observation is [ As shown in Equation (1).

 [Equation 1]

)은 기본 복사 보정 결과로,

는 단위, q는 위상 발사 전에 결정되는 상수, X는 위성의 기상 탑재체에서 취득한 지구 관측값,

은 매 흑체관측마다 계산되는 기울기(nominal slope)이다.Here, R (

) Is the result of basic copy correction,

Where q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is the nominal slope calculated for every blackbody observation.

은 다음 [수학식 2]로,

은 다음 [수학식 3]으로 나타낼 수 있다.

Is expressed by the following equation (2)

Can be expressed by the following equation (3).

&Quot; (2) "

&Quot; (3) "

는 흑체관측 온도로부터 계산한 결과,

는 흑체관측값의 평균,

는 흑체관측 시 우주보기관측값의 평균,

는 지구관측시 우주보기관측값의 평균이다.here,

As a result of calculation from the blackbody observation temperature,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observations.

는 다음 [수학식 4]와 같다.In Equation (2)

Is expressed by the following equation (4).

&Quot; (4) "

는 위성 발사 전 결정된 계수이고,

는 흑체관측에서 측정하는 온도이다.here,

Is the determined coefficient before satellite launch,

Is the temperature measured at the blackbody observation.

The geometric correction unit 124 may perform geometric correction using the data on which the basic radiation correction of the infrared channel is performed.

Geometric calibration is to correct the distorted data according to the original latitude and longitude coordinates due to various reasons such as the position of the satellite located above the equator and the spherical earth surface, satellites thermal deformation, satellites attitude control characteristics .

In this case, the geometric correction includes a process of eliminating various errors caused by changes in the orbit and attitude of the satellite, and a process of determining the exact position of each pixel in the image and a process of determining relative positions To keep the position constant and to keep the position of the designated pixel constant in successive images taken over time.

The image radiometric calibrating module (IRCM) 140 comprises a perfect copy correcting unit 142.

The perfect copy correction unit 142 may receive the data observed in the satellite's meteorological payload 50 and perform a full radiation correction of the infrared channel.

The full radiation correction of the infrared channel is based on the basic photocopying formula and is based on five algorithms: ① Scan Mirror Emissivity Compensation (SMEC), ② Midnight Blackbody Calibration Correction (MBCC) ), ③ slope averageing (SlopeAvg), ④ detector characteristic normalization, and ⑤ 1 / f noise correction (Pink noise Compensation, 1 over f noise compensation).

The perfect radiation correction equation is expressed by the following equation (5).

&Quot; (5) "

)은 완전 복사 보정 결과로,

는 단위, q는 위상 발사 전에 결정되는 상수, X는 기상 탑재체에서 취득한 지구 관측값,

는 매 흑체관측마다 계산되는 기울기(nominal slope),

는

와 우주보기관측 평균 등을 이용하여 계산하는 값,

는 스캔미러 각도(θ)에 따른 스캔미러 방사율,

은 스캔미러의 온도값에 미리 결정된 계수를 곱해서 계산하는 값이다. Here, R (

) Is the result of full radiation correction,

Where q is a constant determined prior to phase firing, X is the Earth Observation value acquired from the meteorological payload,

Is the nominal slope calculated for each blackbody observation,

The

And the average of the observations in space,

The scan mirror emissivity according to the scan mirror angle?

Is a value calculated by multiplying the temperature value of the scan mirror by a predetermined coefficient.

는 [수학식 6]으로,

는 [수학식 7]로 나타낼 수 있다.Equation (5) is a copy correction formula when the scan mirror emissivity correction algorithm (SMEC) is applied,

(6)

Can be expressed by Equation (7).

&Quot; (6) "

&Quot; (7) "

는 흑체관측 및 스캔미러의 온도 등으로 계산되는 값, q는 위상 발사 전에 결정되는 상수,

는 흑체관측값의 평균,

는 흑체관측 시 우주보기관측값의 평균,

는 지구관측시 우주보기관측값의 평균,

는 우주보기관측(

)에서 스캔미러 방사율,

는 우주보기관측에서 스캔미러 온도값에 미리 결정된 계수를 곱해서 계산하는 결과이다.here,

Is a value calculated by the temperature of the blackbody observation and the scan mirror, q is a constant determined before the phase firing,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observation,

Observation of Universe View (

) Scan mirror emissivity,

Is the result of multiplying the scan mirror temperature value by a predetermined factor in an universe viewing observation.

를 계산하는 방법은 다음 [수학식 8]과 같다.In Equation (6)

Is calculated by the following equation (8).

&Quot; (8) "

는 위성 발사 전 결정된 상수, 나머지 세 값(

,

,

)는 매번 계산하여 적용하는데,

는 흑체관측 시 스캔미러 온도값에 미리 결정된 계수를 곱해서 계산하는 결과값이다.here,

Is the determined constant before the satellite launch, the remaining three values (

,

,

) Is calculated and applied every time,

Is a result obtained by multiplying the scan mirror temperature value by a predetermined coefficient at the time of blackbody observation.

, [수학식 6]으로 기울기인

를 각각 계산한 다음, 자정 부근의 흑체관측 보정(MBCC) 및 기울기 평균(SlopeAvg) 알고리즘을 적용하여

를 계산한다. 이후, [수학식 7]로

(

계산식에서

대신에

를 대입해서 계산)를, [수학식 5]로 R을 계산한 다음, 검출기 특성 정규화 및 1/f 잡음 보정 알고리즘을 R에 적용하면 완전 복사 보정 결과를 생성할 수 있다. In order to perform the full radiation correction, first, as shown in [Equation 8]

, ≪ / RTI > (6)

And then the blackbody observation correction (MBCC) and the slope average (SlopeAvg) algorithm around the midnight are applied

. Thereafter, in Equation (7)

(

In the equation

Instead of

, And then applying the detector characteristic normalization and 1 / f noise correction algorithm to R, a full radiation correction result can be generated.

The complete radiated and geometrically corrected data can then be provided to the user via HRIT.

As described above, instead of performing the geometric correction by receiving the complete radiation correction result as an input, the basic radiation correction function can be directly implemented and executed in the image position correcting apparatus 100 to shorten the start time of the geometric correction. That is, the image position correcting apparatus 100 can directly receive the data observed in the gas-phase payload 50, directly generate the radiation correction result, and perform the geometric correction.

In order to perform this method, we need to explain the experimental and analytical methods to verify the accuracy of the geometric correction process, do.

The satellite image preprocessing system 100 is set to apply the full radiation correction formula. In order to change it to the basic copy correction formula, the following setting should be performed in the option of the satellite image preprocessing system 100. [

A, the blackbody observation correction (MBCC) and the slope average (SlopeAvg) algorithm near the midnight, detector characterization, and noise correction calculation are not performed.

=0으로 설정한다.B, Equations (1) and (5) become equal

= 0.

=0으로 설정한다.C, [Equation 3] and [Equation 7] become equal

= 0.

=0으로 설정한다.D, Equations (4) and (8) are equal

= 0.

The four algorithms for selecting the functions of the satellite image preprocessing system 100 are deactivated in A, and the scan mirror emissivity correction algorithm is disabled as B, C, and D because there is no setting function.

The method of performing the geometric correction processing using the result of the radiation correction is summarized as follows.

The basic radiation correction formula is applied to the image preprocessing system for R & D separate from the satellite image preprocessing system currently operating in real time, and geometric correction is performed by using the generated two types of radiation correction results as input to the image position correction device The header file of the geometry correction result is stored.

We compare and analyze the header file of the geometry correction result generated by the real - time satellite image preprocessing system and the header file of the geometry correction result generated by the satellite image preprocessing system for R & D for the same observation.

Hereinafter, the geometric correction processing is performed by inputting the result of generating the copy correction result and the perfect copy correction formula generated using the basic copy correction formula, and then describe the comparative analysis of the geometric correction accuracy of the two processing results do.

The input data is a total of 33 hours of blackbody observations and earth observations between 11:15 on May 30, 2013 and 20:15 on May 31, 2013.

Each of the infrared channels (SWIR: Short Ware IR, WV: Water Vapor, WIN1: Window 1, WIN2: Window 2) of the satellite meteorological payload is composed of two detectors (detector A and detector B) 4 km x 4 km area is detected.

FIG. 2 is a view showing a slope calculation result for the SWIR detector A, and FIG. 3 is a view showing a slope calculation result for the WIN2 detector B. FIG.

, mG는

, mM은

, mC는

를 의미한다. SWIR 검출기 A와 WIN2 검출기 B 모두 다섯 개의 알고리즘이 적용되지 않은

이 나머지 기울기(

,

,

)보다 전체적으로 값이 작음을 알 수 있으며, SWIR 검출기 A의 경우, MBCC 알고리즘을 적용한 두 개의 기울기(

,

)는 자정(도 2 및 도 3에서 UTC 15시) 부근에 급격한 기울기 변화가 두드러지게 감소했음을 볼 수 있다.Referring to Figures 2 and 3,

, mG

, mM

, mC

. Both the SWIR detector A and the WIN2 detector B do not have the five algorithms applied

This remaining slope (

,

,

). In the case of the SWIR detector A, two slopes using the MBCC algorithm (

,

) Shows a sharp decrease in steep slope near the midnight (UTC 15 o'clock in Fig. 2 and Fig. 3).

The following explains the effect of the difference of the correction formula including inclination on the accuracy of geometric correction.

)과 기본 복사 보정식(

)으로 처리한 복사 보정 결과를 입력으로 수행한 기하 보정 결과(Level 1B 헤더파일)를 MI 수행 평가 툴로 분석한 표이다.FIG. 4 is a schematic view of a full-

) And the basic copy correction (

(Level 1B header file) obtained by inputting the copy correction result processed by the MI performance evaluation tool.

In Fig. 4, Navigation shows the difference between the landmark for geometric correction and the position / hardness difference of the image after the geometric correction, WithinFame is the position maintaining accuracy in one image, and Registration is the position for the geometric correction result within 15 minutes and 90 minutes And maintenance accuracy, respectively. All items are measured by east-west (EW) and inter-north (NS), and the unit is μRad.

), 대부분 완전 복사 보정식에 의한 결과(

)보다 정확도가 다소 떨어지지만, 정확도 요구 사항은 모두 만족시킴을 확인할 수 있다.In the case of the geometric correction process performed by inputting the result of the basic copy correction process (

), Mostly the result of perfect copy correction (

), But the accuracy requirements are all satisfied.

FIGS. 5 and 6 are diagrams showing the result of analyzing the geometric correction accuracy using the landmark, using the results of the perfect copy correction formula and the basic copy correction formula, respectively.

Referring to Figures 5 and 6, the straight lines in each graph represent the values of the geometric correction accuracy requirements, and it can be seen that most observations meet the requirements and the two results are similar.

The accuracy of the geometric correction results obtained by the input of the basic copy correction formula and the result of the perfect copy correction formula for the total of 33 hours observation can be confirmed to satisfy the geometric correction accuracy in both cases.

Therefore, the basic copy correction result can be used as an input to the geometric correction processing for the purpose of advancing the time period of the geometric correction processing. That is, it is possible to directly receive the data observed by the image position correcting device 120 in the gas-phase payload 50, perform the basic copy correction of the infrared channel, and then perform the geometric correction, It is possible to satisfy the HRIT transmission requirement of the meteorological payload 50 of the geostationary orbit composite satellite.

Hereinafter, the preprocessing process of the satellite image according to the embodiment of the present invention will be described in more detail.

7 is a flowchart illustrating a preprocessing process of a satellite image according to an embodiment of the present invention.

Referring to FIG. 7, data observed in the satellite payload is received (S700).

Next, the basic copy correction of the infrared channel is performed using the received data (S710).

More specifically, basic copy correction can be performed using Equation (1).

 [Equation 1]

)은 기본 복사 보정 결과로,

는 단위, q는 위상 발사 전에 결정되는 상수, X는 위성의 기상 탑재체에서 취득한 지구 관측값,

은 매 흑체관측마다 계산되는 기울기(nominal slope)이다.Here, R (

) Is the result of basic copy correction,

Where q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is the nominal slope calculated for every blackbody observation.

은 다음 [수학식 2]로,

은 다음 [수학식 3]으로 나타낼 수 있다.

Is expressed by the following equation (2)

Can be expressed by the following equation (3).

&Quot; (2) "

&Quot; (3) "

는 흑체관측 온도로부터 계산한 결과,

는 흑체관측값의 평균,

는 흑체관측 시 우주보기관측값의 평균,

는 지구관측시 우주보기관측값의 평균이다.here,

As a result of calculation from the blackbody observation temperature,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observations.

는 다음 [수학식 4]와 같다.In Equation (2)

Is expressed by the following equation (4).

&Quot; (4) "

는 위성 발사 전 결정된 계수이고,

는 흑체관측에서 측정하는 온도이다.here,

Is the determined coefficient before satellite launch,

Is the temperature measured at the blackbody observation.

Next, the geometric correction is performed using the result of performing the basic copy correction of the infrared channel (S720).

Geometric calibration is to correct the distorted data according to the original latitude and longitude coordinates due to various reasons such as the position of the satellite located above the equator and the spherical earth surface, satellites thermal deformation, satellites attitude control characteristics .

In this case, the geometric correction includes a process of eliminating various errors caused by changes in the orbit and attitude of the satellite, and a process of determining the exact position of each pixel in the image and a process of determining relative positions To keep the position constant and to keep the position of the designated pixel constant in successive images taken over time.

Embodiments of the present invention include a computer-readable medium having program instructions for performing various computer-implemented operations. This medium records a program for executing the above-described satellite image preprocessing method. The medium may include program instructions, data files, data structures, etc., alone or in combination. Examples of such media include magnetic media such as hard disks, floppy disks and magnetic tape, optical recording media such as CD and DVD, programmed instructions such as floptical disk and magneto-optical media, ROM, RAM, And a hardware device configured to store and execute the program. Or such medium may be a transmission medium, such as optical or metal lines, waveguides, etc., including a carrier wave that transmits a signal specifying a program command, data structure, or the like. Examples of program instructions include machine language code such as those produced by a compiler, as well as high-level language code that can be executed by a computer using an interpreter or the like.

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed exemplary embodiments, Of the right.

100: Satellite image preprocessing system
120: video position correcting device
122: basic copy correction unit
124: Geometrical Correction
140: image copy correction device
142:

Claims (7)

A basic copy correction unit that receives data observed from the satellite payload and performs basic copy correction of the infrared channel, and
And a geometric correction unit for performing geometric correction using data on which the basic radiation correction of the infrared channel is performed,
A satellite image preprocessing system. The method of claim 1,
Further comprising a complete radiopaque correcting unit for receiving the data observed in the satellite payload and performing a full radiation correction of the infrared channel. The method of claim 1,
Wherein the basic copy correcting unit comprises:
The basic copy correction of the infrared channel is performed using Equation (1)
[Equation 1]

Where R is the base radiation correction result, q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is a nominal slope calculated for every blackbody observation,

,

As a result of calculation from the blackbody observation temperature,

,

Is the determined coefficient before satellite launch,

The temperature measured at the blackbody observation,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observation,

Satellite image preprocessing system. Receiving the observed data from the satellite ' s meteorological payload and performing basic radiance correction of the infrared channel, and
Performing geometric correction using data on which the basic copy correction of the infrared channel is performed
The method comprising the steps of: 5. The method of claim 4,
Further comprising the step of receiving data observed in the satellite payload and performing a full radiation correction of the infrared channel. 5. The method of claim 4,
Wherein the basic copy correction step comprises:
The basic copy correction of the infrared channel is performed using Equation (1)
[Equation 1]

Where R is the base radiation correction result, q is a constant determined before the phase firing, X is the earth observation value acquired from the satellite meteorological payload,

Is a nominal slope calculated for every blackbody observation,

,

As a result of calculation from the blackbody observation temperature,

,

Is the determined coefficient before satellite launch,

The temperature measured at the blackbody observation,

Is the average of blackbody observations,

Is the average of the observations of the universe in blackbody observations,

Is the average of the observations of the universe in Earth observation,

A preprocessing method for satellite images. A computer-readable recording medium storing a program for causing a computer to execute the method according to any one of claims 4 to 6.

Images