CN116309185B - Water color satellite image land proximity effect correction method based on lookup table - Google Patents
Water color satellite image land proximity effect correction method based on lookup table Download PDFInfo
- Publication number
- CN116309185B CN116309185B CN202310530843.1A CN202310530843A CN116309185B CN 116309185 B CN116309185 B CN 116309185B CN 202310530843 A CN202310530843 A CN 202310530843A CN 116309185 B CN116309185 B CN 116309185B
- Authority
- CN
- China
- Prior art keywords
- land
- wave band
- wave
- proximity effect
- lookup table
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000000694 effects Effects 0.000 title claims abstract description 53
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012937 correction Methods 0.000 title claims abstract description 18
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 13
- 238000002310 reflectometry Methods 0.000 claims description 37
- 230000005855 radiation Effects 0.000 claims description 32
- 239000004973 liquid crystal related substance Substances 0.000 claims description 28
- 239000000443 aerosol Substances 0.000 claims description 10
- 239000013535 sea water Substances 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000005540 biological transmission Effects 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012896 Statistical algorithm Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009360 aquaculture Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G06T5/80—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention provides a water color satellite image land proximity effect correction method based on a lookup table, which comprises the following steps: s1, establishing a quantitative model of adjacent land pixels; s2, calculating land proximity effects of a plurality of wave bands under different conditions by using a Monte Carlo simulation method; s3, calculating contribution of adjacent land pixels of a plurality of wave bands under different conditions, and calculating contribution ratio of adjacent land pixels of two short wave infrared waves; s4, establishing a lookup table; s5, acquiring a satellite observation image, and calculating the contribution ratio of adjacent land pixels of a short wave infrared band SWIR1 and a short wave infrared band SWIR2 in the satellite observation image; s6, corresponding conditions of SWIR1 wave band and SWIR2 wave band in the lookup table are determined; s7, searching the land proximity effect of other wave bands under the corresponding conditions in the lookup table, and realizing the land proximity effect correction of all the wave bands. The invention saves a large amount of calculation time by calculating the lookup table in advance, ensures the authenticity of the calculated land proximity effect and realizes the land proximity effect correction in all wave bands.
Description
Technical Field
The invention relates to the technical field of satellite water color remote sensing detection, in particular to a water color satellite image land proximity effect correction method based on a lookup table.
Background
Satellite water color remote sensing is an important technical means for monitoring the ecological environment of marine and inland water bodies in a large area, however, satellite water color remote sensing products are obtained by inversion of different remote sensing algorithms, so that one key problem to be solved is to quantify the accuracy or uncertainty of satellite remote sensing data products. Satellite remote sensing product accuracy is generally evaluated by in-situ observation data far away from coastline, and an important reason for this is to avoid the influence of land pixels on satellite remote sensing images on offshore water pixels, namely the land proximity effect. The main reason for the land proximity effect is that because the land pixel radiance is much greater than the water radiance, the radiation from adjacent brighter land pixels will scatter into the field of view of the relatively darker water pixels, resulting in an abnormally increased intensity of the observed offshore water radiance on the satellite image. Therefore, the pixels of the offshore water body on the satellite water color remote sensing image are usually masked, for example, the data in 20 km from the MODIS data are affected by the land proximity effect to different degrees. The coastal offshore area is the most important aquaculture area and human activity area in China, and the long-term ecological environment monitoring of the offshore area by utilizing satellite water color remote sensing data has important significance for guaranteeing long-term sustainable development of aquatic product resources and marine ecological environment. Therefore, development of a satellite water color remote sensing data land proximity effect correction method is needed to eliminate errors caused by the land proximity effect and improve the accuracy of satellite water color remote sensing data observation in offshore areas.
The problems existing in the prior art are as follows: the estimation of the point spread function needs to solve a three-dimensional radiation transmission equation, and the three-dimensional radiation transmission equation can be accurately completed only through numerical values, so that a large amount of calculation time is needed; the statistical algorithm based on the satellite remote sensing image can only apply near infrared or short wave infrared wave bands, and the calculation result based on the point spread function method is a simulation result, and does not necessarily represent the actual satellite observation result, and lacks of reality.
Disclosure of Invention
The invention aims to provide a water-based satellite image terrestrial proximity effect correction method based on a lookup table, so as to solve the problems in the background art.
The invention is realized by the following technical scheme: a water color satellite image land proximity effect correction method based on a lookup table is characterized by comprising the following specific steps:
s1, establishing a quantitative model of adjacent land pixels;
s2, calculating land proximity effects of a plurality of wave bands under different conditions based on the land proximity pixel quantitative model by a Monte Carlo simulation method, wherein the different conditions comprise observation geometric conditions, land coverage types, reflectivity, aerosol types and aerosol optical thickness;
s3, calculating adjacent land pixel contributions of a plurality of wave bands under different conditions by using a Monte Carlo simulation method, and calculating an adjacent land pixel contribution ratio between a short wave infrared wave band SWIR1 and a short wave infrared wave band SWIR 2;
s4, establishing a lookup table, wherein the lookup table comprises land proximity effects of a plurality of wave bands under different conditions, adjacent land pixel contributions of the wave bands under different conditions and adjacent land pixel contribution ratio of SWIR1 wave band and SWIR2 wave band;
s5, acquiring a satellite observation image, and calculating the contribution ratio of adjacent land pixels of SWIR1 wave band and SWIR2 wave band in the satellite observation image;
s6, matching the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the satellite observation image with the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the lookup table, and determining the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band in the lookup table;
s7, based on the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band, the land proximity effect of the SWIR1 wave band and the SWIR2 wave band under the corresponding conditions and the land proximity effect of other wave bands under the corresponding conditions are searched in the lookup table, and land proximity effect correction of all wave bands is realized.
Further, in step S1, the established adjacent land pixel quantitative model includes a radiance image function and a land adjacent radiation function.
Further, the expression of the radiance image function is:
wherein, the liquid crystal display device comprises a liquid crystal display device,for radiance image, +.>For radiating sea surface of large air path, < >>For terrestrial or sea radiation, h is the point spread function,>and->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>Is a convolution.
Further, in the near-shore water body or the near-shore sea area, the land near-radiation function is composed of a land radiation intensity function and a sea water radiation intensity function, and the land radiation intensity function is:
=/>
wherein, the liquid crystal display device comprises a liquid crystal display device,for terrestrial radiance, +.>For spatially uniform ground reflectivity, +.>For ground reflectivity of +.>When the sun is descending irradiance, +.>The solar downlink irradiance is when the ground reflectivity is 0, namely when the ground reflectivity is completely absorbed, and S is the large balloon surface reflectivity;
the sea water radiance function:
wherein, the liquid crystal display device comprises a liquid crystal display device,is the radiance of sea water>For sea surface remote sensing reflectance->For the irradiance of the sun down the sea surface, +.>For the Fresnel reflectivity of the sea surface +.>The solar downlink irradiance when the ground reflectivity is 0, namely the solar downlink irradiance when the ground reflectivity is completely absorbed, and S is the large balloon surface reflectivity.
Further, the terrestrial proximity radiation function is:
=/> -/>
in the middle ofIs adjacent to the land and is irradiated by->For sea surface remote sensing reflectance->For a spatially uniform ground reflectivity,the Fresnel reflectivity of the sea surface is that of a big balloon surface, S is that of a big balloon surface>For radiation reflected by the sea surface->For solar downlink irradiance when the ground reflectivity is 0, i.e. when the ground is fully absorbed, h is the point spread function, +.>And->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>For convolution +.>Andas a weight function, +.>A matrix is distinguished for land and water.
Further, calculating adjacent land pixel contributions of a plurality of wave bands under different conditions by using a Monte Carlo simulation method specifically comprises the following steps: spectral dual-sphere reflectivities of the multiple bands at different times were calculated by monte carlo simulation and used to characterize the adjacent land pixel contribution.
Compared with the prior art, the invention has the following beneficial effects:
the invention provides a water color satellite image land proximity effect correction method based on a lookup table, which comprises the following steps:
(1) The satellite image is rapidly processed by calculating the lookup table in advance, so that the calculation processing time is saved by about 90%;
(2) The satellite remote sensing image is utilized to directly calculate the contribution ratio of adjacent land pixels in two short wave infrared bands, so that the authenticity of the contribution ratio of adjacent land pixels in the two bands is ensured;
(3) Adjacent land pixel contribution ratio observed in two short wave infrared bands by using satellite remote sensing imageAnd in different observation geometries, land coverage types, reflectivities calculating adjacent terrestrial pixel contribution ratio in a lookup table for aerosol type and aerosol optical thickness, etc>And matching to obtain the land proximity effect of other wave bands under corresponding conditions, and realizing the land proximity effect correction of all wave bands.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only preferred embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic diagram of a method for correcting a water-color satellite image terrestrial proximity effect based on a lookup table.
FIG. 2 is a flow chart illustrating the operation of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, exemplary embodiments according to the present invention will be described in detail with reference to the accompanying drawings. It should be apparent that the described embodiments are only some embodiments of the present invention and not all embodiments of the present invention, and it should be understood that the present invention is not limited by the example embodiments described herein. Based on the embodiments of the invention described in the present application, all other embodiments that a person skilled in the art would have without inventive effort shall fall within the scope of the invention.
In the following description, numerous specific details are set forth in order to provide a more thorough understanding of the present invention. It will be apparent, however, to one skilled in the art that the invention may be practiced without one or more of these details. In other instances, well-known features have not been described in detail in order to avoid obscuring the invention.
It should be understood that the present invention may be embodied in various forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. As used herein, the term "and/or" includes any and all combinations of the associated listed items.
In order to provide a thorough understanding of the present invention, detailed structures will be presented in the following description in order to illustrate the technical solutions presented by the present invention. Alternative embodiments of the invention are described in detail below, however, the invention may have other implementations in addition to these detailed descriptions.
Referring to fig. 1-2, a method for correcting a terrestrial proximity effect of a water color satellite image based on a lookup table, which is specifically implemented as follows:
step S1, a near land pixel quantitative model is established, namely, the total radiance, near shore water body and near land pixel quantitative relation observed by a satellite sensor are established through theoretical modeling, so that the near land effect of a plurality of wave bands is calculated.
Specifically, the established adjacent land pixel quantitative model comprises a radiance image function and a land adjacent radiation function.
Specifically, the expression of the radiance image function is:
(1)
wherein, the liquid crystal display device comprises a liquid crystal display device,for radiance image, +.>For radiating sea surface of large air path, < >>For terrestrial or sea radiation, h is the point spread function,>and->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>Is a convolution.
Further, in an offshore body of water or offshore area, L sfc Can be divided into land radiation contributionsAnd ocean radiation contribution->Ocean radiation contribution->Can be further divided into sea surface reflected radiation +.>And water-leaving radiation->Therefore, the radiance image +.>The formula becomes:
(2)
wherein, the liquid crystal display device comprises a liquid crystal display device,for radiance image, +.>For radiating sea surface of large air path, < >>Contribution to terrestrial radiation, < >>For a terrestrial and water differentiation matrix, m=0 for water pixels, m=1 for terrestrial pixels, +.>For radiation reflected by the sea surface->Is the sea water radiance, h is the point spread function,>and->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>Is a convolution.
Further, the terrestrial proximity effect may be due to terrestrial proximity radiationFor example, define the difference in radiance for the land and sea conditions, respectively +.>The method comprises the following steps:
(3)
wherein, the liquid crystal display device comprises a liquid crystal display device,is adjacent to the land and is irradiated by->Is the difference between the land radiance and the ocean radiance, +.>Contribution to terrestrial radiation, < >>For radiation reflected by the sea surface->Is the radiance of sea water>For a terrestrial and water differentiation matrix, m=0 for water pixels, m=1 for terrestrial pixels, +.>For radiation reflected by the sea surface, h is the point spread function,>and->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>Is a convolution.
Specifically, in a near-shore water body or a near-shore sea area, the land near-radiation function is composed of a land radiation brightness function and a sea water radiation brightness function, and the land radiation brightness function is as follows:
=/>(4)
wherein, the liquid crystal display device comprises a liquid crystal display device,for terrestrial radiance, +.>For spatially uniform ground reflectivity, +.>For ground reflectivity of +.>When the sun is descending irradiance, +.>The solar downlink irradiance when the ground reflectivity is 0, namely the solar downlink irradiance is completely absorbed, and S is a big balloon surfaceReflectivity;
the sea water radiance function:
(5)
wherein, the liquid crystal display device comprises a liquid crystal display device,for sea surface remote sensing reflectance->For the irradiance of the sun down the sea surface, +.>For the solar downlink irradiance at a ground reflectance of 0, i.e. at full absorption, +.>The Fresnel reflectivity of the sea surface is shown, and the S is the large balloon surface reflectivity.
Substituting equation (4) and equation (5) into equation (3) yields an expression for the terrestrial near radiation function:
=/> -/>(6)
wherein, the liquid crystal display device comprises a liquid crystal display device,
(7)
(8)
in the middle ofIs adjacent to the land and is irradiated by->For sea surface remote sensing reflectance->For a spatially uniform ground reflectivity,the Fresnel reflectivity of the sea surface is that of a big balloon surface, S is that of a big balloon surface>For radiation reflected by the sea surface->For solar downlink irradiance when the ground reflectivity is 0, i.e. when the ground is fully absorbed, h is the point spread function, +.>And->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>For convolution +.>Andfor the terrestrial and water differentiation matrix, m=0 for the water pixel and m=1 for the terrestrial pixel.
Specifically, the atmospheric radiation transmission software 6S can be adopted for simulation calculation to obtain the light source device comprising the large balloon surface albedo S, the atmospheric diffuse transmittance t, the large gas path radiance and the sea surface downward irradianceAtmospheric optical properties, inclusive.
And step S2, calculating land proximity effects of a plurality of wave bands under different conditions based on the land proximity pixel quantitative model by a Monte Carlo simulation method, wherein the different conditions comprise observation geometric conditions, land coverage types, reflectivity, aerosol types and aerosol optical thickness.
And S3, calculating adjacent land pixel contributions of a plurality of wave bands under different conditions by a Monte Carlo simulation method, and calculating adjacent land pixel contribution ratio values in a SWIR1 wave band and a SWIR2 wave band.
Specifically, the Meng Daka lock simulation algorithm is a backward Monte Carlo method (BMC).
Specifically, the adjacent land pixel contribution ratio of the two short wave infrared bands SWIR1 and SWIR2The calculation formula of (2) is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,to simulate the land proximity contribution ratio of two short wave infrared bands,for the short-wave infrared band SWIR1 on land proximity effect, < >>Is a terrestrial proximity effect at the short wave infrared band SWIR 2.
Step S4, a lookup table is established, the lookup table comprises land proximity effects of a plurality of wave bands under different conditions, contribution of adjacent land pixels of the wave bands under different conditions and contribution ratio of adjacent land pixels of two short-wave infrared wave bands SWIR1 and SWIR2, the lookup table is established in advance, images observed by satellites can be rapidly processed, a large amount of time is saved, and correction efficiency in practical application is improved.
Specifically, the calculation formula of the contribution ratio of adjacent land pixels of the SWIR1 and SWIR2 bands calculated by using the satellite observation image is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,the contribution ratio of the terrestrial proximity effect of two short wave infrared bands under satellite observation,for the short-wave infrared band SWIR1 on land proximity effect, < >>Is a terrestrial proximity effect at the short wave infrared band SWIR 2.
And S5, acquiring satellite observation images, calculating contribution ratios of adjacent land pixels of SWIR1 wave bands and SWIR2 wave bands in the satellite observation images, and directly calculating the contribution ratios of the adjacent land pixels by using the satellite-acquired observation images, so that the authenticity of the ratios can be further ensured, and the correction result is more approximate to the actual situation.
Step S6, matching the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the satellite observation image with the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the lookup table, and determining the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band in the lookup table;
and step S7, based on the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band, searching the land proximity effect of the SWIR1 wave band and the SWIR2 wave band under the corresponding conditions and the land proximity effect of other wave bands under the corresponding conditions in a lookup table, and realizing the land proximity effect correction of all wave bands.
The working principle is as follows: firstly, establishing a near land pixel quantitative model through theoretical modeling, and calculating the near land effect of a plurality of wave bands under different conditions through a Monte Carlo simulation methodThe different conditions include observation geometry, land coverage type, reflectivity, aerosol type and aerosol optical thickness, calculating adjacent land pixel contributions of multiple bands in different situations, calculating land adjacent contribution ratios in two short wave infrared bands (SWIR 1 and SWIR 2), and building a lookup table. The adjacent land pixel contribution ratio of the SWIR1 wave band and the SWIR2 wave band is calculated through satellite observation images, the adjacent land pixel contribution ratio is matched with the adjacent land pixel contribution ratio of the SWIR1 wave band and the SWIR2 wave band in the lookup table, corresponding conditions of the SWIR1 wave band and the SWIR2 wave band in the lookup table are determined, and the land proximity effect of the SWIR1 wave band and the SWIR2 wave band and the land proximity effect of other wave bands under corresponding conditions are searched in the lookup table based on the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band, so that the land proximity effect correction of all wave bands is realized.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather to enable any modification, equivalent replacement, improvement or the like to be made within the spirit and principles of the invention.
Claims (6)
1. A water color satellite image land proximity effect correction method based on a lookup table is characterized by comprising the following specific steps:
s1, establishing a quantitative model of adjacent land pixels;
s2, calculating land proximity effects of a plurality of wave bands under different conditions based on the land proximity pixel quantitative model by a Monte Carlo simulation method, wherein the different conditions comprise observation geometric conditions, land coverage types, reflectivity, aerosol types and aerosol optical thickness;
s3, calculating adjacent land pixel contributions of a plurality of wave bands under different conditions by using a Monte Carlo simulation method, and calculating an adjacent land pixel contribution ratio between a short wave infrared wave band SWIR1 and a short wave infrared wave band SWIR 2;
s4, establishing a lookup table, wherein the lookup table comprises land proximity effects of a plurality of wave bands under different conditions, adjacent land pixel contributions of the wave bands under different conditions and adjacent land pixel contribution ratio of SWIR1 wave band and SWIR2 wave band;
s5, acquiring a satellite observation image, and calculating the contribution ratio of adjacent land pixels of SWIR1 wave band and SWIR2 wave band in the satellite observation image;
s6, matching the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the satellite observation image with the contribution ratio of adjacent land pixels between the SWIR1 wave band and the SWIR2 wave band in the lookup table, and determining the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band in the lookup table;
s7, based on the corresponding conditions of the SWIR1 wave band and the SWIR2 wave band, the land proximity effect of the SWIR1 wave band and the SWIR2 wave band under the corresponding conditions and the land proximity effect of other wave bands under the corresponding conditions are searched in the lookup table, and land proximity effect correction of all wave bands is realized.
2. The method of claim 1, wherein in step S1, the established model for determining the amount of adjacent pixels includes a radiance image function and a terrestrial adjacent radiation function.
3. The method for correcting the terrestrial proximity effect of a water-based satellite image according to claim 2, wherein the radiance image function is expressed as:
wherein, the liquid crystal display device comprises a liquid crystal display device,for radiance image, +.>For radiating sea surface of large air path, < >>For terrestrial or sea radiation, h is the point spread function,>and->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>Is a convolution.
4. A method for correcting the terrestrial proximity effect of a water-based satellite image according to claim 3, wherein in a near-shore water body or a near-shore sea area, the terrestrial proximity radiation function is composed of a terrestrial radiance function and a seawater radiance function, and the terrestrial radiance function is:
= />
wherein, the liquid crystal display device comprises a liquid crystal display device,for terrestrial radiance, +.>For spatially uniform ground reflectivity, +.>For ground reflectivity of +.>When the sun is descending irradiance, +.>The solar downlink irradiance is when the ground reflectivity is 0, namely when the ground reflectivity is completely absorbed, and S is the large balloon surface reflectivity;
the sea water radiance function:
wherein, the liquid crystal display device comprises a liquid crystal display device,is the radiance of sea water>For sea surface remote sensing reflectance->On the sea surface for the sunThe irradiance of the downlink light is used to determine,the Fresnel reflectivity of the sea surface is shown, and the S is the large balloon surface reflectivity.
5. The method for correcting the terrestrial proximity effect of a water-based satellite image according to claim 4, wherein the terrestrial proximity radiation function is:
=/> -/>
in the middle ofIs adjacent to the land and is irradiated by->For sea surface remote sensing reflectance->For spatially uniform ground reflectivity, +.>The Fresnel reflectivity of the sea surface is that of a big balloon surface, S is that of a big balloon surface>For radiation reflected by the sea surface->For solar downlink irradiance when the ground reflectivity is 0, i.e. when the ground is fully absorbed, h is the point spread function, +.>And->For the target pixel position +.>And->For the pixel position around the target pixel point, +.>Is sun incident direction vector, +.>For satellite observation direction vector, < >>For convolution +.>And->As a weight function, +.>A matrix is distinguished for land and water.
6. The method for correcting the terrestrial proximity effect of the water-based satellite image according to claim 5, wherein the calculating the contribution of the adjacent terrestrial pixels of the plurality of wave bands under different conditions by the monte carlo simulation method specifically comprises: spectral dual-sphere reflectivities of the multiple bands at different times were calculated by monte carlo simulation and used to characterize the adjacent land pixel contribution.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310530843.1A CN116309185B (en) | 2023-05-12 | 2023-05-12 | Water color satellite image land proximity effect correction method based on lookup table |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310530843.1A CN116309185B (en) | 2023-05-12 | 2023-05-12 | Water color satellite image land proximity effect correction method based on lookup table |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116309185A CN116309185A (en) | 2023-06-23 |
CN116309185B true CN116309185B (en) | 2023-08-01 |
Family
ID=86788964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310530843.1A Active CN116309185B (en) | 2023-05-12 | 2023-05-12 | Water color satellite image land proximity effect correction method based on lookup table |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116309185B (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103204A (en) * | 2011-01-26 | 2011-06-22 | 环境保护部卫星环境应用中心 | Inversion method for land aerosols optical thickness based on environment satellite 1 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109635477A (en) * | 2018-12-20 | 2019-04-16 | 中国农业科学院农业资源与农业区划研究所 | A kind of thermal infrared radiation mode for considering to close on pixel effect |
CN113218874A (en) * | 2021-04-13 | 2021-08-06 | 中国科学院合肥物质科学研究院 | Method and system for obtaining surface target object reflectivity based on remote sensing image |
-
2023
- 2023-05-12 CN CN202310530843.1A patent/CN116309185B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103204A (en) * | 2011-01-26 | 2011-06-22 | 环境保护部卫星环境应用中心 | Inversion method for land aerosols optical thickness based on environment satellite 1 |
Also Published As
Publication number | Publication date |
---|---|
CN116309185A (en) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102176001B (en) | Permeable band ratio factor-based water depth inversion method | |
CN108152289B (en) | Remote sensing indirect monitoring method for total amount of extra-phoma algae in eutrophic lake | |
CN105115941B (en) | A kind of remote sensing inversion method for extracting Complex water body chlorophyll concentration distributed intelligence | |
CN102288956A (en) | Atmospheric correction method for multispectral data of remote sensing satellite | |
CN104089627A (en) | Absolute radiometric calibration method for satellite staring planar array CCD camera | |
CN108256186B (en) | Pixel-by-pixel atmospheric correction method for online calculation lookup table | |
CN107219497A (en) | Cross-radiometric calibration method and system for grand sight angular measurement sensor | |
CN103267533A (en) | Practical high spectrum remote sensing image atmosphere self-correction method and apparatus | |
CN110728642A (en) | Neural network atmosphere correction method for GOCI satellite morning and evening observation | |
Xia et al. | A bathymetry mapping approach combining log-ratio and semianalytical models using four-band multispectral imagery without ground data | |
CN110929423A (en) | Soil water content inversion method of comprehensive drought model | |
CN104318544A (en) | Method for estimating the number of light-induced trapping fishing boats based on satellite remote sensing data at night light | |
CN112414554B (en) | Sea surface salinity obtaining method, device, equipment and medium | |
CN110544236A (en) | Coral reef whitening remote sensing monitoring method based on time series satellite images | |
CN114544452B (en) | Multi-angle polarized water color remote sensor satellite atmosphere correction method | |
CN113610729A (en) | Method and system for correcting satellite-ground collaborative atmosphere of hyperspectral remote sensing image and storage medium | |
CN102073792A (en) | Method for performing inversion on optical property of aerosol of coastal zone by using MODIS image | |
CN111651707A (en) | Tidal level inversion method based on optical shallow water satellite remote sensing image | |
CN111198162B (en) | Remote sensing inversion method for urban surface reflectivity | |
CN114970214A (en) | Aerosol optical thickness inversion method and device | |
CN116309185B (en) | Water color satellite image land proximity effect correction method based on lookup table | |
CN111947628A (en) | Linear water depth inversion method based on inherent optical parameters | |
CN102878985B (en) | Water surface wave scale monitoring method based on image texture features | |
CN113340819B (en) | Water body atmosphere correction method and system based on image self statistical characteristics and storage medium | |
CN111650128B (en) | High-resolution atmospheric aerosol inversion method based on surface reflectivity library |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |