CN113495060A - Method and device for transmitting radiation reference between channels - Google Patents
Method and device for transmitting radiation reference between channels Download PDFInfo
- Publication number
- CN113495060A CN113495060A CN202010197250.4A CN202010197250A CN113495060A CN 113495060 A CN113495060 A CN 113495060A CN 202010197250 A CN202010197250 A CN 202010197250A CN 113495060 A CN113495060 A CN 113495060A
- Authority
- CN
- China
- Prior art keywords
- information
- infrared channel
- visible
- sea surface
- reflectivity
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000002310 reflectometry Methods 0.000 claims abstract description 164
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000002834 transmittance Methods 0.000 claims description 17
- 238000012546 transfer Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 238000004891 communication Methods 0.000 description 5
- 238000012937 correction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 239000005427 atmospheric aerosol Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4738—Diffuse reflection, e.g. also for testing fluids, fibrous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0037—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the heat emitted by liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/007—Radiation pyrometry, e.g. infrared or optical thermometry for earth observation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0074—Radiation pyrometry, e.g. infrared or optical thermometry having separate detection of emissivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
Abstract
The embodiment of the invention provides a method and a device for transmitting radiation reference among channels, wherein the method comprises the following steps: obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area; obtaining reflectivity ratio information of visible-near infrared channels and intermediate infrared channels according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information; and calculating the sea surface reflectivity information of the visible-near infrared channel according to the reflectivity ratio information of the visible-near infrared channel and the sea surface reflectivity information of the infrared channel, and determining the apparent reflectivity of the visible-near infrared channel according to the apparent radiance information of the spot-shining area. And obtaining visible-near infrared channel sea surface reflectivity information, and obtaining visible-near infrared channel apparent reflectivity according to the visible-near infrared channel sea surface reflectivity information, so as to realize radiation reference transmission among channels.
Description
Technical Field
The invention relates to the technical field of satellites, in particular to a method and a device for transmitting radiation references between channels.
Background
Satellite remote sensor in-orbit radiometric calibration can be divided into on-satellite calibration and alternative calibration. Since not all satellites are equipped with on-satellite calibration systems, on-orbit alternate calibration becomes the main means for calibrating satellite remote sensors at present. The on-track substitution calibration method of the visible light channel mainly comprises the following steps: site scaling, cross scaling and scene scaling.
The field calibration is an on-orbit alternative calibration method commonly used for visible-near infrared channels of satellite-borne remote sensors. By selecting a proper radiation correction field, the reflectivity and the atmospheric conditions of the field surface are observed before and after the satellite transit field, and the radiance at the entrance pupil of the satellite is obtained by combining the radiation transmission mode and the directional correction, so that the observation result is calibrated.
The cross calibration method takes an internationally recognized satellite with higher radiometric calibration precision as a reference load, and transmits the radiometric reference with higher calibration precision to a corresponding channel of the load to be calibrated by comparing the near-simultaneous-local observation results of the adjacent channels, so as to realize the calculation of reference transmission and calibration coefficients. The cross calibration can realize the advantages of radiation calibration among high-frequency and multiple remote sensors, is rapidly developed internationally, and is widely applied to radiation calibration of resources, environments, weather and other series satellite instruments in China.
However, in the prior art, the site calibration method needs a large amount of site-synchronous measurement data, is time-consuming and labor-consuming, has limited calibration frequency, such as 1-2 times a year, and is difficult to realize high-frequency calibration, and currently, internationally, the site calibration method has a calibration precision of only 5% at the highest in a reflected solar band, and is difficult to make a major breakthrough; the cross calibration has strict requirements on the channel arrangement, the channel spectral response function, the transit time, the geometric registration precision and the like between the reference remote sensor and the remote sensor to be calibrated, so that the cross calibration result can bring great uncertainty.
Therefore, how to more accurately and efficiently perform high-frequency secondary satellite calibration has become an urgent problem to be solved in the industry.
Disclosure of Invention
Embodiments of the present invention provide a method and an apparatus for transferring radiation reference between channels, so as to solve the technical problems mentioned in the foregoing background art, or at least partially solve the technical problems mentioned in the foregoing background art.
In a first aspect, an embodiment of the present invention provides an inter-channel radiation reference transfer method, including:
obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area;
obtaining visible-near infrared and intermediate infrared channel reflectivity ratio information through intermediate infrared channel sea surface theoretical reflectivity information and near infrared channel sea surface theoretical reflectivity information;
and calculating near-infrared channel sea surface reflectivity information according to the visible-near-infrared and intermediate-infrared channel reflectivity ratio information and the intermediate-infrared channel sea surface reflectivity information, determining the spot-shining region apparent radiance information according to the near-infrared channel sea surface reflectivity information, and determining the near-infrared channel apparent reflectivity according to the spot-shining region apparent radiance information.
Before the step of obtaining the reflectivity information of the intermediate infrared channel sea surface through the upward radiation information of the sea surface of the spot area and the reflected radiation information of the sea surface of the flare area, the method further comprises the following steps:
acquiring the brightness temperature information of the intermediate infrared channel, and combining the brightness temperature information of the intermediate infrared channel with a Planck function to obtain the radiation information emitted by the flare area;
acquiring the observation radiation information of the flare area, subtracting the observation radiation information of the flare area and the emission radiation information of the flare area to obtain the apparent reflection radiation information of the flare area, and performing atmospheric correction to obtain the sea surface reflection radiation information of the flare area.
After the step of obtaining the information of the radiation reflected by the sea surface in the focal spot region, the method further comprises the following steps:
and obtaining the upward radiation information of the sea surface in the blazing area according to the solar spectrum radiation data, the distance information of the day and the ground and the atmospheric absorption rate information.
Before the step of obtaining the visible-near infrared and intermediate infrared channel reflectivity ratio information through the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information, the method comprises the following steps:
calculating to obtain Fresnel reflectivity information of the intermediate infrared channel according to the incident angle information of the intermediate infrared channel and the refractive index of the solar surface of the flare area, and determining the theoretical reflectivity information of the sea surface of the intermediate infrared channel according to the Fresnel reflectivity information of the intermediate infrared channel;
and calculating to obtain near-infrared channel Fresnel reflectivity information according to the visible-near-infrared channel incident angle information and the flare area ocean surface refractive index, and determining near-infrared channel sea surface theoretical reflectivity information according to the near-infrared channel Fresnel reflectivity information.
The method comprises the following steps of determining apparent radiance information of a spot-shining area according to sea surface reflectivity information of a visible-near infrared channel, and specifically comprises the following steps:
according to the atmospheric parameter information and the observation set of the flare area, obtaining visible-near infrared channel atmospheric uplink radiation information, solar direction transmittance information and observation direction transmittance information by means of radiation transmission;
obtaining visible-near infrared channel solar reflection radiation information according to the visible-near infrared channel sea surface reflectivity information, the solar direction transmittance information and the observation direction transmittance information;
and obtaining apparent radiance information of the flare area according to the visible-near infrared channel water-leaving radiation information, the visible-near infrared channel solar reflection radiation information and the visible-near infrared channel atmospheric uplink radiation information.
In a second aspect, an embodiment of the present invention provides an inter-channel radiation reference transferring apparatus, including:
the first calculation module is used for obtaining the reflectivity information of the intermediate infrared channel sea surface through the upward radiation information of the sea surface in the spot-shining area and the reflected radiation information of the sea surface in the flare area;
the second calculation module is used for obtaining the reflectivity ratio information of the near-infrared channel and the intermediate-infrared channel through the intermediate-infrared channel sea surface theoretical reflectivity information and the near-infrared channel sea surface theoretical reflectivity information;
the reference transfer module is used for calculating visible-near infrared channel sea surface reflectivity information according to the near infrared channel reflectivity ratio information, the intermediate infrared channel sea surface reflectivity information, determining the blazing area apparent radiance information according to the near infrared channel sea surface reflectivity information, and determining the visible-near infrared channel apparent reflectivity according to the blazing area apparent radiance information.
In a third aspect, an embodiment of the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor executes the computer program to implement the steps of the inter-channel radiation reference transferring method according to the first aspect.
In a fourth aspect, an embodiment of the present invention provides a non-transitory computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps of the inter-channel radiation reference transferring method according to the first aspect.
According to the method and the device for transmitting the radiation standard between the channels, the sea surface emissivity information of the middle infrared channel is transmitted to the visible-near infrared channel by utilizing the fixed ratio relation of the reflectivity of the sea surface flare area observed between the channels, the sea surface reflectivity information of the near infrared channel is obtained, the apparent reflectivity of the near infrared channel is calculated according to the sea surface reflectivity information of the near infrared channel, and the radiation standard transmission between the channels is realized.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a flow chart illustrating an exemplary method for inter-channel radiation reference propagation according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an inter-channel radiation reference transferring apparatus according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The reference transfer process described in the embodiment of the present invention transfers radiometric calibration accurate reading of a channel with high calibration accuracy of an instrument to a channel to be examined, thereby implementing absolute calibration of the channel to be examined. The reference transmission process comprises reference transmission between stars, and can also be used for transmitting the reference from a channel with high calibration precision to a channel with low calibration precision in the same load.
When an un-calibrated channel and an accurate calibrated channel on a satellite remote sensing instrument observe the same target, a proper ground calibrated target is selected, the radiation of the un-calibrated channel can be calculated by using the channel with higher calibration precision, the radiation reference of the channel with higher calibration precision is transmitted to the un-calibrated channel by combining the correction of errors possibly existing in target emission and atmospheric influence, and the calibration method is an inter-channel reference transmission method.
In the observation of the cloud-free blazed region of the ocean, the reflectivities of different wavelengths are in a stable linear relationship under certain specific conditions. The sunlight reflected energy of the intermediate infrared band is usually weak, but in the solar flare area, the reflected solar reflected energy detected by the intermediate infrared band is considerable, and the intermediate infrared channel of the meteorological satellite remote sensor performs on-satellite calibration by utilizing on-satellite black bodies and cold air observation, so that the calibration precision is high. By means of the high-precision separation of the emitted radiation and the reflected radiation of the observation results of the mid-infrared channel in the solar blaze area and the characteristic that the reflectivity of different channels is in a linear ratio relation in the solar blaze area, the calibration precision of the mid-infrared channel can be transmitted to other visible light channels, and further radiation reference transmission among the channels is achieved.
Fig. 1 is a schematic flow chart of an inter-channel radiation reference transferring method according to an embodiment of the present invention, as shown in fig. 1, including:
step S1, obtaining the reflectivity information of the intermediate infrared channel sea surface through the upward radiation information of the sea surface in the spot-shining area and the reflected radiation information of the sea surface in the flare area;
step S2, obtaining the ratio information of the near-infrared channel reflectivity and the intermediate-infrared channel reflectivity through the intermediate-infrared channel sea surface theoretical reflectivity information and the near-infrared channel sea surface theoretical reflectivity information;
and step S3, calculating visible-near infrared channel sea surface reflectivity information according to the near infrared channel reflectivity ratio information and the mid infrared channel sea surface reflectivity information, determining the apparent radiance information of the spot exposure area according to the visible-near infrared channel sea surface reflectivity information, and determining the apparent reflectivity of the visible-near infrared channel according to the apparent radiance information of the spot exposure area.
Specifically, step S1 in the embodiment of the present invention is specifically to obtain radiation information Rad in the sea surface of the infrared channel in the blazed area according to the solar spectrum radiation data, the distance between the day and the ground information, and the atmospheric absorption rate informationsun_mid_calCalculating the brightness temperature of the middle infrared channel to obtain the radiation information emitted by the middle infrared channel in the flare area, and determining the difference between the observed radiation information of the middle infrared channel in the flare area and the calculated radiation information emitted by the middle infrared channel in the flare area as the radiation information reflected by the sea surface of the middle infrared channel in the flare areasun_mid_surfTherefore, the intermediate infrared channel sea surface reflectivity information is specifically expressed as:
wherein, Radsun_mid_surfFor reflecting radiation information, Rad, from the sea surface of the infrared channel in the blazed zonesun_mid_calAnd radiating information on the sea surface of the infrared channel in the blazed spot region.
Specifically, in step S2 of the embodiment of the present invention, since the solar flare area is a result of the fresnel reflection on the water surface, and the fresnel reflection depends on the incident angle and the refractive index of the water body, and the refractive index is independent with respect to the wavelength, if the refractive indexes of the water surfaces with different wavelengths are known, the theoretical reflectivity of the channel fresnel without the existence of the atmosphere on the calm sea surface can be calculated.
The Fresnel formula can be expressed as:
wherein the subscripts H and V represent the horizontal and vertical components, respectively, and n1And n2Complex refractive indices, θ, of medium 1 and medium 2, respectivelyiIs the angle of incidence. Thus the Fresnel reflectivity ρ can be expressed as:
in the case of light emitted from the atmosphere into the ocean,
if the refractive index of the ocean is known, no atmosphere exists and the sea surface is smooth, the sea surface reflectivity theoretically satisfies the Fresnel reflection theorem:
in a quiet channel without atmospheric influence, the ratio relation of the reflectivity between the channels is the ratio of the Fresnel reflectivity of the two channels, the atmospheric influence on the ocean surface is considered, the atmospheric transmittance in the solar incidence direction is obtained by utilizing a radiation transmission mode simulation, and therefore the sea surface reflectivity can be expressed as follows when the atmospheric exists:
by the formula (1) to the formula (4), the reflectivity ratio information of the near-infrared channel and the intermediate-infrared channel can be obtained.
Step S3 described in the embodiment of the present invention is specifically that, for the ocean face flare region subjected to quality control, the energy source received by the sensor is divided into four parts, that is, single or multiple scattering of atmospheric aerosol or molecules, near-infrared channel solar reflection radiation information, radiation reaching the sensor after transmission of the off-water radiation, and radiation reaching the sensor after transmission of the atmosphere after reflection of the ocean face white cap. Generally, the flare area selected by the scheme is usually selected to be an area with the wind speed of 7m/s or less, so that the radiation of the ocean surface white cap part is negligible. For off-water radiation, studies have shown that in the near infrared band, e.g. lambda > 750nm,
visible light lambda<750nm, it has been found that the ratio of TOA radiation at a certain wavelength to the water-free radiation of the non-flare region is approximately constant, and the absolute value of the water-free radiance of the flare region is constant with respect to the non-flare region, and therefore, the absolute value of the water-free radiance of the flare region is constant at wavelengths less than λ<And in the 750nm wave band, the out-of-water radiation value of the channel can be calculated by collecting the TOA radiance of a non-blazed spot region and exists as a constant term. The energy of the near-infrared channel received at the satellite entrance pupil, i.e. the flare table radiance information, can thus be expressed as:
wherein the content of the first and second substances,
for near-infrared channel atmospheric uplink radiation information, EiAs solar radiation energy information, thetasunIs the solar zenith angle at the current observation moment,
for the near infrared channel sea surface reflectivity,
is the information of the solar direction transmittance of the near infrared channel,
the information of the transmittance in the visible-near infrared channel observation direction is shown, and the information of the visible-near infrared channel water-leaving radiation is shown as C.
After the apparent radiance information of the visible-near infrared channel flare area is obtained, the apparent reflectivity of the visible-near infrared channel can be calculated by combining the solar spectrum radiation, the distance between the sun and the ground and the sun zenith angle corresponding to the observation time, and the radiation transfer among the channels is realized.
According to the embodiment of the invention, sea surface emissivity information of the middle infrared channel is transmitted to the visible-near infrared channel by using the fixed ratio relation of the reflectivity of the sea surface observed by the ocean surface flare area, so that the visible-near infrared channel sea surface reflectivity information is obtained, and the visible-near infrared channel apparent reflectivity is calculated according to the visible-near infrared channel sea surface reflection information, so that the transmission of radiation reference among the channels is realized.
On the basis of the above embodiment, before the step of obtaining the reflectivity information of the mid-infrared channel sea surface by the upward radiation information of the sea surface in the flare area and the reflected radiation information of the sea surface in the flare area, the method further includes:
acquiring the brightness temperature information of the intermediate infrared channel, and combining the brightness temperature information of the intermediate infrared channel with a Planck function to obtain the radiation information emitted by the flare area;
acquiring the observation radiation information of the flare area, and performing difference calculation on the observation radiation information of the flare area and the observation radiation information of the flare area to obtain the sea surface reflection radiation information of the flare area.
Specifically, the radiation transmission mode is utilized to simulate the brightness temperature of the intermediate infrared channel and the thermal infrared channel under different atmospheric conditions, the least square fitting is carried out on the brightness temperature of the intermediate infrared channel according to the result of the radiation transmission mode, and a statistical model of the brightness temperature of radiation emitted by the intermediate infrared channel in the ocean surface flare area is established, namely:
Tmid_IR＝f(TIR_1,TIR_2)
wherein, TIR_1And TIR_2Light temperature, T, of two thermal infrared channels, respectivelymid_IRAnd obtaining the brightness temperature information of the intermediate infrared channel for the brightness temperature of the intermediate infrared channel needing to be calculated.
The intermediate infrared channel radiation of satellite observation comprises two parts of emission radiation and reflection radiation, the intermediate infrared channel brightness temperature information obtained by the previous calculation is utilized, and the Planck function is combined to calculate the radiation information P (T) emitted by the spot-shining areamid_IR) And the difference is calculated between the observation radiation information of the flare area and the radiation information emitted by the flare area to obtain the sea surface reflection radiation information of the initial flare area, namely:
Radsun_mid＝Radmid-P(Tmid_IR)
wherein, RadmidObserving radiation information for flare areas, P (T)mid_IR) Radiation information is emitted for the flare area.
The method comprises the following steps of considering the absorption of gases such as water vapor and carbon dioxide in the atmosphere, and performing atmospheric correction on the sea surface reflection radiation information of the initial speckle region by using a radiation transmission mode to obtain the speckle region sea surface reflection radiation information, namely:
Radsun_mid_surf＝Radsun_mid·τ
on the basis of the above embodiment, after the step of obtaining the information of the reflected radiation of the sea surface in the focal spot region, the method further includes:
obtaining the upward radiation information of the sea surface in the blazing area according to the solar spectrum radiation data, the distance information of the day and the ground and the atmospheric absorption rate information: combining solar spectral radiation data, distance from the sun to the earth, simulated path atmospheric absorption, the upward radiation of solar radiation reaching the sea surface within the tunnel can be calculated:
wherein, Isun(v) is solar spectral radiation data, τ (v) is atmospheric transmittance, and φ (v) is a spectral response function of the channel. Thus, the sea surface emissivity of the mid-infrared channel can be expressed as:
on the basis of the above embodiment, before the step of obtaining the visible-near infrared and intermediate infrared channel reflectivity ratio information through the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information, the method includes:
calculating to obtain the Fresnel reflectivity information of the intermediate infrared channel according to the incident angle information of the intermediate infrared channel, and determining the sea surface theoretical reflectivity information of the intermediate infrared channel according to the Fresnel reflectivity information of the intermediate infrared channel;
and calculating to obtain visible-near infrared channel Fresnel reflectivity information according to the visible-near infrared channel incidence angle information, and determining visible-near infrared channel sea surface theoretical reflectivity information according to the visible-near infrared channel Fresnel reflectivity information.
Specifically, the solar flare area is a result of the water surface fresnel reflection, and since the fresnel reflection depends on the refractive index of the water body and the refractive index is independent with respect to the wavelength, if the water surface refractive indexes with different wavelengths are known, the channel fresnel theoretical reflectivity on a calm sea surface and without atmosphere can be calculated.
In the case of light emitted from the atmosphere into the ocean,
if the refractive index of the ocean is known, no atmosphere exists and the sea surface is smooth, the sea surface reflectivity theoretically satisfies the Fresnel reflection theorem:
in a quiet channel without atmospheric influence, the ratio relation of the reflectivity between the channels is the ratio of the Fresnel reflectivity of the two channels, the atmospheric influence on the ocean surface is considered, the atmospheric transmittance in the solar incidence direction is obtained by utilizing a radiation transmission mode simulation, and therefore the sea surface reflectivity can be expressed as follows when the atmospheric exists:
the method is used for respectively obtaining the sea surface theoretical reflectivity information of the intermediate infrared channel and the sea surface theoretical reflectivity information of the visible-near infrared channel, and then obtaining the reflectivity ratio information of the visible-near infrared channel and the intermediate infrared channel.
The visible-near infrared to mid-infrared channel reflectance ratio information can be expressed as:
where ρ isFs_2And ρFs_1Respectively representing the theoretical reflectivity information of the sea surface of the intermediate infrared channel and the theoretical reflectivity information of the sea surface of the visible-near infrared channel.
On the basis of the above embodiment, the step of determining the apparent radiance information of the flare area according to the sea surface reflectivity information of the visible-near infrared channel specifically includes:
the contents of atmospheric aerosol, water vapor and ozone in the spot-shining region can be obtained through a satellite data product, and according to the atmospheric parameter information, the solar zenith angle and the satellite zenith angle in the spot-shining region, near-infrared channel atmospheric uplink radiation information, near-infrared channel solar direction transmittance information and visible-near infrared channel observation direction transmittance information are obtained through radiation transmission;
obtaining visible-near infrared channel solar reflection radiation information according to visible-near infrared channel sea surface emissivity information, visible-near infrared channel solar direction transmittance information and visible-near infrared channel observation direction transmittance information;
and obtaining apparent radiance information of the flare area according to the visible-near infrared channel water-leaving radiation information, the visible-near infrared channel solar reflection radiation information and the visible-near infrared channel atmospheric uplink radiation information.
On the basis of the above embodiment, the step of determining the apparent reflectivity of the visible-near infrared channel according to the apparent radiance information of the blazed region specifically includes:
and calculating the solar radiation of the top of the atmospheric layer according to the solar spectrum radiation information, the distance between the day and the ground and the solar zenith angle information at the current moment and the apparent radiance information of the spot-shining area, wherein the ratio of the courseware-near infrared channel apparent radiation calculated in the previous step to the radiation value of the top of the atmospheric layer is the visible-near infrared channel apparent reflectivity.
Fig. 2 is a schematic structural diagram of an inter-channel radiation reference transferring apparatus according to an embodiment of the present invention, as shown in fig. 2, including: a first calculation module 210, a second calculation module 220, and a reference transfer module 230; the first calculation module 210 is configured to obtain reflectivity information of the intermediate infrared channel sea surface through upward radiation information of the sea surface in the flare area and reflected radiation information of the sea surface in the flare area; the second calculating module 220 is configured to obtain reflectivity ratio information of the visible-near infrared channel and the intermediate infrared channel according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information; the reference transfer module 230 is configured to calculate visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determine flare region apparent radiance information according to the visible-near infrared channel sea surface reflectivity information, and determine visible-near infrared channel apparent reflectivity according to the flare region apparent radiance information.
The apparatus provided in the embodiment of the present invention is used for executing the above method embodiments, and for details of the process and the details, reference is made to the above embodiments, which are not described herein again.
According to the embodiment of the invention, sea surface emissivity information of the middle infrared channel is transmitted to the visible-near infrared channel by utilizing the fixed ratio relation of the reflectivity of the sea surface observed by the ocean surface flare area, so that the sea surface reflectivity information of the visible-near infrared channel is obtained, and the apparent reflectivity of the visible-near infrared channel is calculated according to the sea surface reflection information of the visible-near infrared channel, so that the radiation reference transmission between the channels is realized.
Fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 3, the electronic device may include: a processor (processor)310, a communication Interface (communication Interface)320, a memory (memory)330 and a communication bus 340, wherein the processor 310, the communication Interface 320 and the memory 330 communicate with each other via the communication bus 340. The processor 310 may call logic instructions in the memory 330 to perform the following method: obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area; obtaining reflectivity ratio information of visible-near infrared channels and intermediate infrared channels according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information; and calculating visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determining according to the visible-near infrared channel sea surface reflectivity information, and determining the visible-near infrared channel apparent reflectivity according to the speckle region apparent radiance information.
In addition, the logic instructions in the memory 330 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
An embodiment of the present invention discloses a computer program product, which includes a computer program stored on a non-transitory computer readable storage medium, where the computer program includes program instructions, and when the program instructions are executed by a computer, the computer can execute the method provided by the above method embodiments, for example, the method includes: obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area; obtaining reflectivity ratio information of visible-near infrared channels and intermediate infrared channels according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information; and calculating visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determining according to the visible-near infrared channel sea surface reflectivity information, and determining the visible-near infrared channel apparent reflectivity according to the speckle region apparent radiance information.
Embodiments of the present invention provide a non-transitory computer-readable storage medium storing server instructions, where the server instructions cause a computer to execute the method provided in the foregoing embodiments, for example, the method includes: obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area; obtaining reflectivity ratio information of visible-near infrared channels and intermediate infrared channels according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information; and calculating visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determining according to the visible-near infrared channel sea surface reflectivity information, and determining the visible-near infrared channel apparent reflectivity according to the speckle region apparent radiance information.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods of the various embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (9)
1. An interchannel radiation reference transfer method, comprising:
obtaining the reflectivity information of the sea surface of the intermediate infrared channel through the upward radiation information of the sea surface of the spot-shining area and the reflected radiation information of the sea surface of the flare area;
obtaining reflectivity ratio information of visible-near infrared channels and intermediate infrared channels according to the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information;
and calculating visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determining spot-shining region apparent radiance information according to the visible-near infrared channel sea surface reflectivity information, and determining visible-near infrared channel apparent reflectivity according to the spot-shining region apparent radiance information.
2. The inter-channel radiation reference transfer method of claim 1, wherein, prior to said step of obtaining mid-infrared channel sea surface reflectivity information by radiating information up on a blazed sea surface and reflecting radiation information off a flare sea surface, said method further comprises:
acquiring middle infrared channel brightness temperature information, and combining the middle infrared channel brightness temperature information with a Planck function to obtain flare area emission radiation information;
acquiring the observation radiation information of the flare area, and calculating the difference between the observation radiation information of the flare area and the radiation information emitted by the flare area to obtain the sea surface reflection radiation information of the flare area.
3. The inter-channel radiation reference transferring method as claimed in claim 2, wherein after the step of obtaining the information of the blazed sea surface reflected radiation, the method further comprises:
and obtaining the upward radiation information of the sea surface in the blazing area according to the solar spectrum radiation data, the distance information of the day and the ground and the atmospheric absorption rate information.
4. The method for transferring radiation reference between channels according to claim 1, wherein before the step of obtaining the reflectivity ratio information of the visible-near infrared channel and the mid-infrared channel through the mid-infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information, the method comprises:
calculating to obtain Fresnel reflectivity information of the intermediate infrared channel according to the incident angle information of the intermediate infrared channel, and determining sea surface theoretical reflectivity information of the intermediate infrared channel according to the Fresnel reflectivity information of the intermediate infrared channel;
and calculating to obtain visible-near infrared channel Fresnel reflectivity information according to the visible-near infrared channel incidence angle information, and determining visible-near infrared channel sea surface theoretical reflectivity information according to the visible-near infrared channel Fresnel reflectivity information.
5. The inter-channel radiation reference transfer method according to claim 1, wherein the step of determining apparent radiance information of the flare region according to the sea surface reflectivity information of the visible-near infrared channel specifically comprises:
according to the atmospheric parameter information and the observation set of the flare area, obtaining visible-near infrared channel atmospheric uplink radiation information, visible-near infrared channel solar direction transmittance information and visible-near infrared channel observation direction transmittance information by means of radiation transmission;
obtaining visible-near infrared channel solar reflection radiation information according to the visible-near infrared channel sea surface emissivity information, the visible-near infrared channel solar direction transmittance information and the visible-near infrared channel observation direction transmittance information;
and obtaining apparent radiance information of the flare area according to the visible-near infrared channel water-leaving radiation information, the visible-near infrared channel solar reflection radiation information and the visible-near infrared channel atmospheric uplink radiation information.
6. The method for transferring radiation reference between channels according to claim 1, wherein the step of determining the apparent reflectivity of the visible-near infrared channel according to the apparent radiance information of the blazed region specifically comprises:
and calculating to obtain the visible-near infrared channel apparent reflectivity according to the solar spectrum radiation information, the day-ground distance and the sun zenith angle information at the current moment and by combining the apparent radiance information of the spot-shining area.
7. An interchannel radiation reference transfer device comprising:
the first calculation module is used for obtaining the reflectivity information of the intermediate infrared channel sea surface through the upward radiation information of the sea surface in the spot-shining area and the reflected radiation information of the sea surface in the flare area;
the second calculation module is used for obtaining reflectivity ratio information of the visible-near infrared channel and the intermediate infrared channel through the intermediate infrared channel sea surface theoretical reflectivity information and the visible-near infrared channel sea surface theoretical reflectivity information;
and the reference transmission module is used for calculating visible-near infrared channel sea surface reflectivity information according to the visible-near infrared and intermediate infrared channel reflectivity ratio information and the intermediate infrared channel sea surface reflectivity information, determining the spot-shining region apparent radiance information according to the visible-near infrared channel sea surface reflectivity information, and determining the visible-near infrared channel apparent reflectivity according to the spot-shining region apparent radiance information.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the inter-channel radiation reference transfer method of any of claims 1 to 7.
9. A non-transitory computer readable storage medium having stored thereon a computer program, which when executed by a processor, performs the steps of the inter-channel radiation reference transfer method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010197250.4A CN113495060A (en) | 2020-03-19 | 2020-03-19 | Method and device for transmitting radiation reference between channels |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010197250.4A CN113495060A (en) | 2020-03-19 | 2020-03-19 | Method and device for transmitting radiation reference between channels |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113495060A true CN113495060A (en) | 2021-10-12 |
Family
ID=77993584
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010197250.4A Pending CN113495060A (en) | 2020-03-19 | 2020-03-19 | Method and device for transmitting radiation reference between channels |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113495060A (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100108885A1 (en) * | 2007-03-27 | 2010-05-06 | Dsam London Llp | Imaging system with an optical mirror arrangement and self righting housing |
| CN103293522A (en) * | 2013-05-08 | 2013-09-11 | 中国科学院光电研究院 | Intermediate infrared two-channel remote sensing data surface temperature inversion method and device |
| CN105092055A (en) * | 2015-08-21 | 2015-11-25 | 国家卫星气象中心 | Cold cloud target-based weather satellite solar reflection band radiometric calibration method |
| CN105157827A (en) * | 2015-04-23 | 2015-12-16 | 中国科学院长春光学精密机械与物理研究所 | Integrating sphere transmission radiometer for ontrack spectrum radiation absolute calibration for remote sensor |
| CN106547840A (en) * | 2016-10-13 | 2017-03-29 | 国家卫星气象中心 | A kind of parsing of global three-dimensional atmospheric data and management method |
| CN207540669U (en) * | 2017-12-11 | 2018-06-26 | 中国人民解放军63908部队 | optical radiation measuring unit |
| CN108680534A (en) * | 2018-03-28 | 2018-10-19 | 北京大学 | Optical channel In-flight calibration verification method based on medium-wave infrared reflectivity benchmark |
| CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
| CN109883957A (en) * | 2018-12-21 | 2019-06-14 | 中国资源卫星应用中心 | Apparent reflectance model building method, system and calibrating method based on MODIS image |
-
2020
- 2020-03-19 CN CN202010197250.4A patent/CN113495060A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100108885A1 (en) * | 2007-03-27 | 2010-05-06 | Dsam London Llp | Imaging system with an optical mirror arrangement and self righting housing |
| CN103293522A (en) * | 2013-05-08 | 2013-09-11 | 中国科学院光电研究院 | Intermediate infrared two-channel remote sensing data surface temperature inversion method and device |
| CN105157827A (en) * | 2015-04-23 | 2015-12-16 | 中国科学院长春光学精密机械与物理研究所 | Integrating sphere transmission radiometer for ontrack spectrum radiation absolute calibration for remote sensor |
| CN105092055A (en) * | 2015-08-21 | 2015-11-25 | 国家卫星气象中心 | Cold cloud target-based weather satellite solar reflection band radiometric calibration method |
| CN106547840A (en) * | 2016-10-13 | 2017-03-29 | 国家卫星气象中心 | A kind of parsing of global three-dimensional atmospheric data and management method |
| CN207540669U (en) * | 2017-12-11 | 2018-06-26 | 中国人民解放军63908部队 | optical radiation measuring unit |
| CN108680534A (en) * | 2018-03-28 | 2018-10-19 | 北京大学 | Optical channel In-flight calibration verification method based on medium-wave infrared reflectivity benchmark |
| CN108896188A (en) * | 2018-06-06 | 2018-11-27 | 北京空间飞行器总体设计部 | A kind of high rail remote sensing satellite infrared channel atmospheric correction method of Shared aperture |
| CN109883957A (en) * | 2018-12-21 | 2019-06-14 | 中国资源卫星应用中心 | Apparent reflectance model building method, system and calibrating method based on MODIS image |
Non-Patent Citations (7)
| Title |
|---|
| Y. -G. QIAN 等: "Land Surface Temperature Retrieval Using Nighttime Mid-Infrared Channels Data From Airborne Hyperspectral Scanner", IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING, vol. 8, no. 3, 20 November 2014 (2014-11-20), pages 1208 - 1216, XP011576964, DOI: 10.1109/JSTARS.2014.2364853 * |
| YANG, ZD 等: "Overview of FY-3 Payload and Ground Application System", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol. 50, no. 12, 31 December 2012 (2012-12-31), pages 4846 - 4853, XP011472689, DOI: 10.1109/TGRS.2012.2197826 * |
| 孙立微: "空间高光谱遥感仪器辐射定标技术研究", 中国优秀博士学位论文全文数据库信息科技辑, no. 10, 15 October 2018 (2018-10-15), pages 140 - 59 * |
| 彭妮娜;易维宁;方勇华;: "400～1000nm波段反演气溶胶光学厚度的暗像元法", 红外与激光工程, no. 05, 25 October 2008 (2008-10-25), pages 878 - 883 * |
| 景欣 等: "基于海表耀斑中红外基准NPP/VIIRS 反射通道在轨验证", 红外与毫米波学报, vol. 36, no. 6, 31 December 2017 (2017-12-31), pages 694 - 700 * |
| 王玉鹏 等: "可在轨溯源的太阳反射波段光学遥感仪器辐射定标基准传递链路", 光学精密工程, vol. 23, no. 07, 15 July 2015 (2015-07-15), pages 1807 - 1812 * |
| 陈申玮 等: "空间辐射基准传递不确定性的光谱敏感性分析", 光学学报, vol. 38, no. 01, 20 September 2017 (2017-09-20), pages 348 - 361 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yoshida et al. | Retrieval algorithm for CO 2 and CH 4 column abundances from short-wavelength infrared spectral observations by the Greenhouse gases observing satellite | |
| Koepke | Vicarious satellite calibration in the solar spectral range by means of calculated radiances and its application to Meteosat | |
| Higurashi et al. | Development of a two-channel aerosol retrieval algorithm on a global scale using NOAA AVHRR | |
| CN102901516B (en) | A kind of multispectral image radiation correction method based on absolute radiometric calibration | |
| Dufour et al. | Spaceborne estimate of atmospheric CO 2 column by use of the differential absorption method: error analysis | |
| VERMOTE et al. | In-flight calibration of large field of view sensors at short wavelengths using Rayleigh scattering | |
| Newman et al. | Temperature and salinity dependence of sea surface emissivity in the thermal infrared | |
| CN102288956B (en) | Atmospheric correction method for multispectral data of remote sensing satellite | |
| Ahn et al. | Vicarious calibration of the geostationary ocean color imager | |
| CN102103203A (en) | Environmental satellite 1-based surface temperature single-window inversion method | |
| CN102338871A (en) | Method and device for calculating reflectivity of earth surface | |
| CN111795936A (en) | Multispectral remote sensing image atmospheric correction system and method based on lookup table and storage medium | |
| Nalli et al. | Emissivity and reflection model for calculating unpolarized isotropic water surface-leaving radiance in the infrared. I: Theoretical development and calculations | |
| US5075856A (en) | System for modelling low resolution atmospheric propagation | |
| Pan et al. | An improved spectral optimization algorithm for atmospheric correction over turbid coastal waters: A case study from the Changjiang (Yangtze) estuary and the adjacent coast | |
| Brennan et al. | Anisotropic reflectance characteristics of natural earth surfaces | |
| CN108680534B (en) | Optical channel on-orbit calibration verification method based on medium wave infrared reflectivity reference | |
| CN107656289A (en) | Spaceborne optics load absolute radiation calibration method and system based on ground spoke brightness | |
| CN108120510A (en) | A kind of in-orbit absolute radiation calibration method of optical sensor based on reflection mirror array | |
| CN106017678A (en) | Thermal infrared high spectral remote sensing data on-track spectral calibration method | |
| CN114970214A (en) | Aerosol optical thickness inversion method and device | |
| CN114218786A (en) | On-orbit polarization radiation characteristic inversion method for non-polarization satellite sensor | |
| Menzies et al. | Lidar in-space technology experiment measurements of sea surface directional reflectance and the link to surface wind speed | |
| Lin et al. | Radiative transfer simulations of the two-dimensional ocean glint reflectance and determination of the sea surface roughness | |
| CN106198460B (en) | Random fluctuation sea surface albedo full-spectrum segment calculation method applied to marine environment |
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 |