CN115909053A - Water transparency inversion method and system based on remote sensing image - Google Patents
Water transparency inversion method and system based on remote sensing image Download PDFInfo
- Publication number
- CN115909053A CN115909053A CN202211336003.3A CN202211336003A CN115909053A CN 115909053 A CN115909053 A CN 115909053A CN 202211336003 A CN202211336003 A CN 202211336003A CN 115909053 A CN115909053 A CN 115909053A
- Authority
- CN
- China
- Prior art keywords
- water body
- water
- remote sensing
- reflectivity
- determining
- 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
Images
Classifications
-
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a water transparency inversion method and a system based on remote sensing images, relating to the field of water color remote sensing, wherein the water transparency inversion method comprises the following steps: acquiring a water body remote sensing image; preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image; determining the remote sensing reflectivity of the water body according to the reflectivity of the remote sensing image of the water body; according to the remote sensing reflectivity of the water body, the inherent optical characteristics of the water body are obtained by adopting a QAA algorithm, and the inherent optical characteristics of the water body comprise: absorption coefficient and backscattering coefficient; determining a diffusion attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body; determining the FUI index of the water body according to the remote sensing reflectivity of the water body; and determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body. The invention can accurately monitor the transparency of the water body in a large range by an inversion mode, and meets the monitoring requirement of the transparency, namely the water quality parameter of the water body.
Description
Technical Field
The invention relates to the field of water color remote sensing, in particular to a remote sensing image-based water transparency inversion method and system.
Background
The transparency of water is an important basic physical parameter describing the optical properties of water, and is related to the concentration of chlorophyll, inorganic suspended matters and organic yellow substances in the water. In addition, the transparency of seawater is closely related to factors such as solar radiation on the surface of the water body, physical and chemical properties of the water body, meteorological conditions and the like. The transparency can visually reflect the pollution condition of the water body and plays an important role in monitoring the water quality of the water body.
The traditional method for measuring the seawater transparency is to use a Secchi disk (Secchi disk) on a ship for field measurement, but the method can only obtain the transparency of a measuring point, cannot obtain the seawater transparency distribution characteristic of large space-time distribution, and cannot meet the monitoring requirement of the water quality parameter of the water body, namely the transparency.
Therefore, how to accurately monitor the transparency of the water body in a large range becomes a technical problem which needs to be solved urgently in the field.
Disclosure of Invention
The invention aims to provide a remote sensing image-based water transparency inversion method and system, and the water transparency can be accurately monitored in a large range in an inversion mode.
In order to achieve the purpose, the invention provides the following scheme:
a remote sensing image-based water transparency inversion method comprises the following steps:
acquiring a water body remote sensing image;
preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, wherein the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
determining the remote sensing reflectivity of the water body according to the reflectivity of the remote sensing image of the water body;
obtaining the inherent optical characteristics of the water body by adopting a QAA algorithm according to the remote sensing reflectivity of the water body, wherein the inherent optical characteristics of the water body comprise: absorption coefficient and backscattering coefficient;
determining a diffusion attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, wherein the semi-analytical model is a function related to the absorption coefficient and the backscattering coefficient;
determining the FUI index of the water body according to the remote sensing reflectivity of the water body;
and determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body.
The invention also provides a water transparency inversion system based on the remote sensing image, which comprises:
the image acquisition unit is used for acquiring a water body remote sensing image;
the preprocessing unit is used for preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, and the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
the remote sensing reflectivity determining unit is used for determining the remote sensing reflectivity of the water body according to the reflectivity of the water body remote sensing image;
an inherent optical characteristic determining unit, configured to obtain an inherent optical characteristic of the water body by using a QAA algorithm according to the remote sensing reflectivity of the water body, where the inherent optical characteristic of the water body includes: absorption coefficient and backscattering coefficient;
a diffusion attenuation coefficient determination unit for determining a diffusion attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, the semi-analytical model being a function related to the absorption coefficient and the backscattering coefficient;
the FUI index determining unit is used for determining the FUI index of the water body according to the remote sensing reflectivity of the water body;
and the water transparency determining unit is used for determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a water transparency inversion method and a water transparency inversion system based on remote sensing images, wherein the water transparency inversion method comprises the following steps: acquiring a water body remote sensing image; preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, wherein the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration; determining the remote sensing reflectivity of the water body according to the reflectivity of the remote sensing image of the water body; obtaining the inherent optical characteristics of the water body by adopting a QAA algorithm according to the remote sensing reflectivity of the water body, wherein the inherent optical characteristics of the water body comprise: absorption coefficient and backscattering coefficient; determining a diffusion attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, wherein the semi-analytical model is a function related to the absorption coefficient and the backscattering coefficient; determining the FUI index of the water body according to the remote sensing reflectivity of the water body; and determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body. According to the invention, the reflectivity of the water body remote sensing image is determined according to the water body remote sensing image, the remote sensing reflectivity of the water body is determined according to the reflectivity of the water body remote sensing image, the inherent optical characteristic of the water body and the FUI index of the water body are respectively determined according to the remote sensing reflectivity of the water body, the diffusion attenuation coefficient is determined according to the inherent optical characteristic of the water body, and finally the water body transparency is determined according to the diffusion attenuation coefficient and the FUI index of the water body.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
Fig. 1 is a flowchart of a water transparency inversion method based on remote sensing images according to embodiment 1 of the present invention;
fig. 2 is a block diagram of a structure of a remote sensing image-based water transparency inversion system according to embodiment 2 of the present invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. 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 invention aims to provide a remote sensing image-based water transparency inversion method and system, and the water transparency can be accurately monitored in a large range through an inversion mode.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Example 1:
referring to fig. 1, the invention provides a water transparency inversion method based on remote sensing images, which includes the following steps:
s1: acquiring a water body remote sensing image;
s2: preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, wherein the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
radiometric calibration is a process in which when a user calculates spectral reflectivity or spectral radiance of a ground feature, or when images acquired by different sensors at different times need to be compared, luminance gray-scale values of the images must be converted into absolute radiance.
Atmospheric correction means that the total radiance of the ground target finally measured by the sensor is not a reflection of the true reflectivity of the ground, including the radiant quantity error caused by atmospheric absorption, especially scattering. Atmospheric correction is the process of inverting the real surface reflectivity of the ground object by eliminating the radiation errors caused by atmospheric influence.
Orthorectification generally refers to a process of resampling an image into an orthorectified image by selecting some ground control points on the image, and performing tilt correction and projective difference correction on the image simultaneously by using Digital Elevation Model (DEM) data in the range of the image, which has been acquired originally.
Geometric registration refers to an operation process of completely superimposing image points of the same name on the position and the orientation by geometric transformation of images (data) of the same region obtained by different remote sensor systems at different times and different wave bands.
S3: determining the remote sensing reflectivity of the water body according to the reflectivity of the remote sensing image of the water body; the specific calculation formula is as follows:
wherein R is rs And (lambda) is the remote sensing reflectivity of the water body, and R (lambda) is the reflectivity of the remote sensing image of the water body.
S4: obtaining the inherent optical characteristics of the water body by adopting a QAA algorithm according to the remote sensing reflectivity of the water body, wherein the inherent optical characteristics of the water body comprise: absorption coefficient and backscattering coefficient;
the parameters used by the QAA algorithm are shown in table 1.
TABLE 1 QAA Algorithm parameter Table
| Parameter(s) | Means of | Unit |
| R rs (λ) | Remote sensing reflectivity of water body | sr -1 |
| r rs (λ) | Remote sensing reflectivity below water surface | sr -1 |
| a(λ) | Total absorption coefficient | m -1 |
| b b (λ) | Total backscattering coefficient | m -1 |
| a w (λ) | Absorption coefficient of pure water | m -1 |
| b bw (λ) | Back scattering coefficient of pure water | m -1 |
| b bp (λ) | Backscattering coefficient of suspended particles | m -1 |
| λ 0 | Reference wavelength | nm |
| η | Power parameter | / |
| λ | Current wavelength | nm |
The remote sensing reflectivity below the surface of the water body can be calculated through the remote sensing reflectivity of the water body, and the formula is as follows:
wherein r is rs (λ) is the remote-sensing reflectance, R, below the surface of the body of water rs (lambda) is the remote sensing reflectivity of the water body, T is an intermediate parameter, and T = T _ T + /n 2 ，t - Is the radiation transmittance, t, of seawater from top to bottom + Is the radiation transmittance of seawater from bottom to top, n is the refractive index of water, gamma is the reflection coefficient of a water-gas interface, lambda is the current wavelength, Q is the ratio of upward irradiance to downward irradiance of the water surface, and in most ocean and coastal waters, the value of Q is compared with R rx (lambda) and r rs The conversion relationship between (lambda) has little influence, so according to the calculation result of Hydrolight, T ≈ 0.52 and gamma Q ≈ 1.7.
According to theoretical analysis and the numerical simulation result of the radiation transmission equation, the remote sensing reflectivity r below the surface of the water body rs (λ), is the total absorption coefficient a (λ) and the total backscattering coefficient b b (λ) is a function of. In particularThe functional relationship is shown in formulas (3) and (4):
r rs (λ)＝g 0 u(λ)+g 1 [u(λ)] 2 (3)
wherein u (λ) is the ratio of the total backscattering coefficient to the sum of the total absorption coefficient and the total backscattering coefficient, g 0 、g 1 Are all constants determined experimentally, in this example, g 0 、g 1 Values of (d) are averages of the Gordon and Lee suggested values, 0.089 and 0.125, respectively.
From the above formula, r can be selected from rs And (lambda) calculating the value of u (lambda), and obtaining the ratio of the total backscattering coefficient to the sum of the total absorption coefficient and the total backscattering coefficient. If the total absorption coefficient a (lambda) is solved, the total backscattering coefficient b can be solved b (λ) and vice versa. Wherein the backscattering coefficient can be expressed by formula (5):
b b (λ)＝b bw (λ)+b bp (λ) (5)
wherein, b b (λ) is the total backscattering coefficient, b bw (lambda) is the backscattering coefficient of pure water, b bp (λ) is the backscattering coefficient of the suspended particles, λ is the current wavelength.
From equation (3) we can obtain:
since u (λ) is the ratio of the total backscattering coefficient to the sum of the total absorption coefficient and the total backscattering coefficient, if the total absorption coefficient a (λ) is known, the total backscattering coefficient b can be solved according to the following equation b (λ)：
And if the overall backscattering coefficient b is known b (λ), the overall absorption coefficient a (λ) can also be solved:
the total absorption coefficient a (λ) for different wavelengths can be shown by equation (9):
a(λ)＝a w (λ)+Δa(λ) (9)
wherein, a w (lambda) is the absorption coefficient of pure water, lambda is the current wavelength, Δ a (lambda) is the absorption coefficient contributed by dissolved and suspended particles in water, the value of Δ a (lambda) is small when the wavelength is long, and thus a (lambda) is mainly composed of a w (lambda) especially for low or medium nutrient water bodies.
While the backscattering coefficients of suspended particles at other wavelengths can be solved from the backscattering coefficient of suspended particles at the reference wavelength, as shown in equation (10):
wherein, b bp (λ 0 ) The backscattering coefficient of the suspended particles at the reference wavelength is eta, the power parameter is obtained by calculation through the formula (11), and then backscattering coefficients of other wavelengths can be solved, so that the total backscattering coefficient of the corresponding wavelength can be solved.
S5: the inherent optical characteristics of the water body including the total absorption coefficient and the total backscattering parameters can be obtained by solving through a QAA algorithm, and then the diffusion attenuation coefficient K needs to be solved through a semi-analysis model developed by Lee et al based on a radiation transmission equation through the absorption coefficient and the backscattering parameters d (λ), the formula is as follows:
K d (λ)＝(1+0.005θ a )a(λ)+3.47×b b (λ) (12)
wherein, theta a For the solar zenith angle, the parameters in the model were obtained by a Hydrolight simulation of the average particle phase function.
S6: and (4) solving the remote sensing reflectivity of the water body according to the step (S3), and solving the FUI (Forel-UleIndex) index of the water body through the following calculation process.
(1) And integrating to calculate CIE color tristimulus values X, Y and Z of the water body. Setting the relative energy spectral distribution S (lambda) of the irradiation light source to be 1, and then substituting the remote sensing reflectivity of the water body as the object spectral reflectivity rho (lambda) into a formula (12) to match with the color
And &>
The product of the above steps is integrated within a visible light range (380 nm-700 nm) to obtain the CIE color tristimulus value of the water body. Wherein K is a regulatory factor.
(2) Since x + y + z =1, a color can be determined by using two values of x and y, and therefore, the chromaticity coordinates (x, y) of the water body can be directly calculated. And (5) substituting X, Y and Z into the formulae (15) and (16) for normalization calculation to obtain X and Y.
(3) And calculating the color angle of the water body. And (4) substituting the chromaticity coordinates (x, y) of the water body into an equation (17) to calculate the chromaticity angle alpha of the water body. Where ARCTAN2 is a bivariate arctangent function with thresholds of (0 deg., 360 deg.).
α＝ARCTAN2(x-0.333，y-0.3333) (17)
(4) And calculating the FUI index of the water body. According to the chroma angle alpha of the water body and an FUI lookup table (shown in table 2), a chroma value closest to alpha is found out from the FUI lookup table, and an FUI index corresponding to the chroma value is the FUI index of the water body.
TABLE 2 FUI lookup Table
| FUI | x | y | α | FUI | x | y | α |
| 1 | 0.191363 | 0.166919 | 40.467 | 12 | 0.402416 | 0.4811 | 205.0622 |
| 2 | 0.198954 | 0.199871 | 45.19626 | 13 | 0.416243 | 0.47368 | 210.5766 |
| 3 | 0.210015 | 0.2399 | 52.85273 | 14 | 0.431336 | 0.465513 | 216.5569 |
| 4 | 0.226522 | 0.288347 | 67.16945 | 15 | 0.445679 | 0.457605 | 222.1153 |
| 5 | 0.245871 | 0.335281 | 91.29804 | 16 | 0.460605 | 0.449426 | 227.6293 |
| 6 | 0.266229 | 0.37617 | 122.5852 | 17 | 0.475326 | 0.440985 | 232.8302 |
| 7 | 0.290789 | 0.411528 | 151.4792 | 18 | 0.488676 | 0.43285 | 237.3523 |
| 8 | 0.315369 | 0.440027 | 170.4629 | 19 | 0.503316 | 0.424618 | 241.7592 |
| 9 | 0.336658 | 0.461684 | 181.4983 | 20 | 0.515498 | 0.416136 | 245.5513 |
| 10 | 0.363277 | 0.476353 | 191.8352 | 21 | 0.528252 | 0.408319 | 248.9529 |
| 11 | 0.386188 | 0.486566 | 199.0383 |
S7: substituting the diffusion attenuation coefficient obtained by the solution in the step S5 and the FUI index of the water body obtained in the step S6 into a formula (18) to finish the transparency Z to the water body SD And (4) solving.
The accuracy of the established model is verified and evaluated by experimentally measured data. In this embodiment, the transparency is selected to be inverted by using the actual measurement data of the yellow sea area and the MODIS product data of the area to verify the accuracy of the improved model, and the actual measurement data of the partial transparency is shown in table 3. And compared to the conventional Poole-Atkins inversion algorithm proposed by Robert o.megard & Tom Berman (when parameter a equals constant 1.54, as shown in equation 19). The results of the precision alignment are shown in Table 4.
Table 3 partial transparency actual point data
On the basis of the traditional Poole-Atkins mode, the improved model introduces a FUI index to determine the parameters A in the mode, wherein RMSE and MRE are 1.8307 and 43.74 percent respectively. And the inversion accuracy of the improved model is higher than that of the traditional Poole-Atkins inversion algorithm proposed by Robert O.
TABLE 4 improved model error analysis table
In conclusion, the FUI (Forel-Ule) water color index has good research prospect and advantages in the aspect of monitoring water quality by using a remote sensing technology. The invention can realize large-range monitoring of the water transparency, and can greatly improve the calculation precision according to the introduction of the FUI index.
Example 2:
referring to fig. 2, the present invention provides a remote sensing image-based water transparency inversion system, which includes:
the image acquisition unit 1 is used for acquiring a water body remote sensing image;
the preprocessing unit 2 is configured to preprocess the water body remote sensing image to obtain a reflectivity of the water body remote sensing image, where the preprocessing includes: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
the remote sensing reflectivity determining unit 3 is used for determining the remote sensing reflectivity of the water body according to the reflectivity of the water body remote sensing image;
an inherent optical characteristic determining unit 4, configured to obtain an inherent optical characteristic of the water body by using a QAA algorithm according to the remote sensing reflectivity of the water body, where the inherent optical characteristic of the water body includes: absorption coefficient and backscattering coefficient;
a diffuse attenuation coefficient determination unit 5, configured to determine a diffuse attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, where the semi-analytical model is a function related to the absorption coefficient and the backscattering coefficient;
the FUI index determining unit 6 is used for determining the FUI index of the water body according to the remote sensing reflectivity of the water body;
and the water transparency determining unit 7 is used for determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (8)
1. A remote sensing image-based water transparency inversion method is characterized by comprising the following steps:
acquiring a water body remote sensing image;
preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, wherein the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
determining the remote sensing reflectivity of the water body according to the reflectivity of the remote sensing image of the water body;
obtaining the inherent optical characteristics of the water body by adopting a QAA algorithm according to the remote sensing reflectivity of the water body, wherein the inherent optical characteristics of the water body comprise: absorption coefficient and backscattering coefficient;
determining a diffusion attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, wherein the semi-analytical model is a function related to the absorption coefficient and the backscattering coefficient;
determining the FUI index of the water body according to the remote sensing reflectivity of the water body;
and determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body.
2. The remote-sensing-image-based water transparency inversion method according to claim 1, wherein the calculation formula for determining the remote-sensing reflectivity of the water according to the reflectivity of the remote-sensing image of the water is as follows:
wherein R is rs And (lambda) is the remote sensing reflectivity of the water body, and R (lambda) is the reflectivity of the remote sensing image of the water body.
3. The remote-sensing-image-based water transparency inversion method according to claim 1, wherein the method for obtaining the inherent optical characteristics of the water body by adopting a QAA algorithm according to the remote-sensing reflectivity of the water body specifically comprises the following steps:
the calculation formula of the absorption coefficient is as follows:
a(λ)＝a w (λ)+Δa(λ)
wherein a (lambda) is the total absorption coefficient, a w (λ) is the absorption coefficient of pure water, Δ a (λ) is the absorption coefficient contributed by dissolved and suspended particles in water, λ is the current wavelength;
the calculation formula of the backscattering coefficient is as follows:
b b (λ)＝b bw(λ) +b bp (λ)
wherein, b b (λ) is the total backscattering coefficient, b bw (lambda) is the backscattering coefficient of pure water, b bp (lambda) is the backscattering coefficient of the suspended particles, lambda is the current wavelength;
b bp (λ 0 ) Is the backscattering coefficient of the suspended particles at a reference wavelength, eta is a power parameter, and>
r rs (lambda) is the remote sensing reflectivity below the surface of the water body, device for selecting or keeping>
R rs (lambda) is the remote sensing reflectivity of the water body, T is an intermediate parameter, and T = T - t + /n 2 ，t _ Is the radiation transmittance, t, of seawater from top to bottom + Is the radiation transmittance of seawater from bottom to top, n is the refractive index of water, gamma is the reflection coefficient of a water-gas interface, Q is the ratio of upward irradiance to downward irradiance of the water surface, and lambda 0 Is the reference wavelength.
4. The remote-sensing-image-based water transparency inversion method according to claim 1, wherein the calculation formula for determining the diffuse attenuation coefficient of the semi-analytical model according to the inherent optical characteristics of the water is as follows:
K d (λ)＝(1+0.005θ a )a(λ)+3.47×b b (λ)
wherein, K d (λ) is the diffuse attenuation coefficient, θ a Is the zenith angle of the sun, a (lambda) is the total absorption coefficient, b b (λ) is the total backscattering coefficient.
5. The remote-sensing-image-based water transparency inversion method according to claim 1, wherein the determining the FUI index of the water body according to the remote-sensing reflectivity of the water body specifically comprises:
and calculating CIE color tristimulus values of the water body by integration according to the remote sensing reflectivity of the water body, wherein the calculation formula is as follows:
wherein X, Y and Z are CIE color tristimulus values of the water body, S (lambda) is relative energy spectrum distribution, rho (lambda) is object spectral reflectivity which is equal to the remote sensing reflectivity of the water body,
and &>
Is a color matching function, lambda is the current wavelength, K is an adjustment factor>
Determining the chromaticity coordinate of the water body according to the CIE color tristimulus value of the water body, wherein the calculation formula of the chromaticity coordinate of the water body is as follows:
wherein, (x, y) is the chromaticity coordinate of the water body;
determining the chromaticity angle of the water body according to the chromaticity coordinate of the water body, wherein the calculation formula of the chromaticity angle of the water body is as follows:
α＝ARCTAN2(x-0.333，y-0.3333)
wherein alpha is the chroma angle of the water body, and ARCTAN2 is a bivariate arctangent function;
and performing table lookup according to the hue angle of the water body to determine the FUI index of the water body.
6. The remote-sensing-image-based water transparency inversion method according to claim 1, wherein the calculation formula for determining the water transparency according to the diffusion attenuation coefficient and the FUI index of the water is as follows:
wherein, Z SD FUI is the transparency of water body, and FUI is the FUI index of water body, K d And (lambda) is the diffuse attenuation coefficient.
7. A remote sensing image-based water transparency inversion system is characterized by comprising:
the image acquisition unit is used for acquiring a water body remote sensing image;
the preprocessing unit is used for preprocessing the water body remote sensing image to obtain the reflectivity of the water body remote sensing image, and the preprocessing comprises the following steps: radiometric calibration, atmospheric correction, orthorectification, and geometric registration;
the remote sensing reflectivity determining unit is used for determining the remote sensing reflectivity of the water body according to the reflectivity of the water body remote sensing image;
an inherent optical characteristic determining unit, configured to obtain an inherent optical characteristic of the water body by using a QAA algorithm according to the remote sensing reflectivity of the water body, where the inherent optical characteristic of the water body includes: absorption coefficient and backscattering coefficient;
a diffuse attenuation coefficient determination unit for determining a diffuse attenuation coefficient of a semi-analytical model according to the inherent optical characteristics of the water body, the semi-analytical model being a function related to the absorption coefficient and the backscattering coefficient;
the FUI index determining unit is used for determining the FUI index of the water body according to the remote sensing reflectivity of the water body;
and the water transparency determining unit is used for determining the transparency of the water body according to the diffusion attenuation coefficient and the FUI index of the water body.
8. The remote-sensing-image-based water transparency inversion system of claim 7, wherein the FUI index determination unit comprises:
the CIE color tristimulus value calculation module is used for calculating CIE color tristimulus values of the water body by integration according to the remote sensing reflectivity of the water body;
the chromaticity coordinate determination module is used for determining the chromaticity coordinate of the water body according to the CIE color tristimulus value of the water body;
the color angle determining module is used for determining the color angle of the water body according to the chromaticity coordinate of the water body;
and the FUI index determining module is used for performing table lookup according to the color angle of the water body and determining the FUI index of the water body.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211336003.3A CN115909053A (en) | 2022-10-28 | 2022-10-28 | Water transparency inversion method and system based on remote sensing image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211336003.3A CN115909053A (en) | 2022-10-28 | 2022-10-28 | Water transparency inversion method and system based on remote sensing image |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115909053A true CN115909053A (en) | 2023-04-04 |
Family
ID=86487516
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211336003.3A Pending CN115909053A (en) | 2022-10-28 | 2022-10-28 | Water transparency inversion method and system based on remote sensing image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115909053A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113406015A (en) * | 2021-05-31 | 2021-09-17 | 内蒙古师范大学 | Transparency calculation method and system for water bodies of near-shore and inland waters |
| CN113793374A (en) * | 2021-09-01 | 2021-12-14 | 自然资源部第二海洋研究所 | Method for inverting water depth based on water quality inversion result by using improved four-waveband remote sensing image QAA algorithm |
| CN114894719A (en) * | 2022-03-16 | 2022-08-12 | 自然资源部国土卫星遥感应用中心 | Satellite water transparency inversion method based on wave band remote sensing reflectivity ratio |
-
2022
- 2022-10-28 CN CN202211336003.3A patent/CN115909053A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113406015A (en) * | 2021-05-31 | 2021-09-17 | 内蒙古师范大学 | Transparency calculation method and system for water bodies of near-shore and inland waters |
| CN113793374A (en) * | 2021-09-01 | 2021-12-14 | 自然资源部第二海洋研究所 | Method for inverting water depth based on water quality inversion result by using improved four-waveband remote sensing image QAA algorithm |
| CN114894719A (en) * | 2022-03-16 | 2022-08-12 | 自然资源部国土卫星遥感应用中心 | Satellite water transparency inversion method based on wave band remote sensing reflectivity ratio |
Non-Patent Citations (1)
| Title |
|---|
| JIE ZHAN,ET AL.: "MODIS-Based Research on Secchi Disk Depth Using an Improved Semianalytical Algorithm in the Yellow Sea", IEEE XPLORE, pages 1 - 9 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hakvoort et al. | Towards airborne remote sensing of water quality in The Netherlands—validation and error analysis | |
| Mishra et al. | Characterizing the vertical diffuse attenuation coefficient for downwelling irradiance in coastal waters: Implications for water penetration by high resolution satellite data | |
| CN112051226B (en) | Method for estimating total suspended matter concentration of offshore area based on unmanned aerial vehicle-mounted hyperspectral image | |
| CN109325973B (en) | Urban river network area water body atmosphere correction method | |
| Mao et al. | A new approach to estimate the aerosol scattering ratios for the atmospheric correction of satellite remote sensing data in coastal regions | |
| CN103499815A (en) | Method for conducting inland water body atmospheric correction based on oxygen and water-vapor absorption wavebands | |
| CN113506235B (en) | Adaptive weighted spectrum reconstruction method for resisting exposure change | |
| CN114219994A (en) | Ocean optical satellite radiometric calibration method based on air-sea cooperative observation | |
| CN111191380B (en) | Atmospheric aerosol optical thickness estimation method and device based on measurement data of foundation spectrometer | |
| Zhang et al. | A validation study of an improved SWIR iterative atmospheric correction algorithm for MODIS-Aqua measurements in Lake Taihu, China | |
| CN113639716A (en) | Depth residual shrinkage network-based water depth remote sensing inversion method | |
| CN111651707A (en) | Tidal level inversion method based on optical shallow water satellite remote sensing image | |
| CN102901563A (en) | Method and device for determining land surface emissivity of narrow band and broad band simultaneously | |
| CN113066057B (en) | Aerosol optical thickness monitoring method | |
| Popp | Correcting atmospheric masking to retrieve the spectral albedo of land surfaces from satellite measurements | |
| CN106769895A (en) | A kind of method for demarcating measurement whole atmosphere spectral transmittance | |
| CN113763272A (en) | Remote sensing inversion method for photosynthetic effective radiation attenuation coefficient of eutrophic lake | |
| Bulgarelli et al. | Seasonal impact of adjacency effects on ocean color radiometry at the AAOT validation site | |
| CN115909053A (en) | Water transparency inversion method and system based on remote sensing image | |
| Rossini et al. | Maps of sky relative radiance and luminance distributions acquired with a monochromatic CCD camera | |
| CN113029977B (en) | Automatic cross radiometric calibration method for wide-field-angle multispectral sensor | |
| CN110702228B (en) | Edge radiation correction method for aviation hyperspectral image | |
| CN114646601A (en) | Satellite ozone vertical profile inversion method and system based on multi-band coupling | |
| KR20100028337A (en) | Method for calibration of coms throught cloud using cloud target | |
| CN114674461A (en) | Method and device for determining sea surface temperature and readable storage medium |
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 |