Shallow sea water depth inversion method and system based on multi-angle optical radiation transmission model

Abstract

The invention provides a shallow sea water depth inversion method and system based on a multi-angle optical radiation transmission model, and relates to the technical field of remote sensing space information. The method comprises the following steps: on the basis of a water body biological optical model, constructing a shallow sea water depth inversion model in combination with a water body and water bottom optical radiation transmission process; and according to the obtained multi-angle satellite remote sensing observation data, calculating and obtaining shallow sea water depth information by using the shallow sea water depth inversion model. According to the invention, by fusing multi-angle satellite remote sensing observation data and a physical model and comprehensively considering water body absorption, scattering and bottom reflection characteristics, a high-robustness and high-precision shallow sea water depth inversion method and system are constructed. The method has adaptive capacity to a complex optical environment in a shallow water area, can realize large-range and high-precision water depth inversion, provides a new technical means for shallow water depth measurement, and has wide application prospects.

Description

Technical Field

[0001] The invention relates to the field of remote sensing space information technology, and in particular to a shallow sea water depth inversion method and system based on a multi-angle optical radiation transmission model. Background Art

[0002] With the continuous expansion of human activities, the importance of water depth information in shallow sea areas has become increasingly prominent in the fields of marine ecological protection, resource development, waterway maintenance, and disaster assessment. Although traditional water depth measurement technologies, such as multi-beam echo sounders and laser depth sounders, can provide high-precision water depth data, they have the limitations of expensive equipment, complex operation, low measurement efficiency, and restrictions on weather and operating conditions. It is difficult to meet the needs of real-time, fast, and economical water depth information acquisition in large areas of shallow sea. Therefore, shallow sea water depth inversion technology based on remote sensing images has become one of the research hotspots due to its advantages such as wide coverage, low cost, and high efficiency.

[0003] The core of remote sensing water depth inversion is to convert optical signals into water depth information. Remote sensing reflection signals in shallow sea areas are not only affected by the absorption and scattering characteristics of water bodies, but also by the combined effects of bottom reflection and water depth, so they are significantly complex. At present, remote sensing-based water depth inversion methods face many challenges, such as nonlinear changes in optical signals caused by complex seabed topography, optical signal confounding caused by changes in suspended matter concentration and water color, and accurate separation of the superimposed effects of bottom reflection on water optical properties in shallow water areas. These technical bottlenecks limit the applicability and scalability of existing inversion models under different water conditions.

[0004] At present, water depth inversion based on optical remote sensing data is mainly carried out at home and abroad, using multispectral images such as Landsat TM / ETM+, IKONOS, and GF series. According to the property of sunlight penetrating water bodies, theoretical models based on radiation transfer equations, empirical models based on experience and assumed mathematical models, and semi-theoretical and semi-empirical models combining theoretical models with empirical parameters have gradually been formed.

[0005] Summarizing existing research findings, most traditional optical remote sensing water depth inversion methods only use optical image data acquired from a single angle on a satellite or aerial platform. For example, only the spectral information of the zenith observation angle is used for water depth inversion. This method ignores the different scattering and reflection characteristics of light when it is transmitted at multiple angles in the water body, and cannot fully utilize the multi-angle optical signal differences generated by the interaction between the water body and the seabed. In shallow sea areas, the seabed topography is complex, and the reflection direction and intensity of light from the seabed with different slopes and roughness vary significantly at multiple angles. It is difficult for single-angle data to fully capture this information, thereby limiting the accuracy and reliability of water depth inversion, especially in shallow sea areas with large terrain undulations.

[0006] In addition, most existing methods rely on ideal environments such as high bottom reflectivity and clear water, and fail to fully consider the comprehensive interaction between light and complex factors such as water composition and seabed topography during multi-angle transmission in water, resulting in inaccurate inversion of shallow sea depth and significant restrictions on its application. In summary, existing shallow sea depth inversion technology still has many limitations. Summary of the invention

[0007] In view of the above problems, the purpose of the present invention is to provide a shallow sea water depth inversion method and system based on a multi-angle optical radiation transmission model, so as to improve the accuracy, reliability and versatility of shallow sea water depth inversion by more comprehensively considering various factors in the multi-angle optical radiation transmission process.

[0008] In order to solve the above technical problems, the present invention provides the following technical solutions:

[0009] On the one hand, a shallow sea water depth inversion method based on a multi-angle optical radiation transmission model is provided, the method comprising the following steps:

[0010] S1. Based on the water body bio-optical model and combined with the light radiation transmission process of the water body and the bottom, a shallow sea water depth inversion model is constructed;

[0011] S2. Based on the acquired multi-angle satellite remote sensing observation data, the shallow sea water depth information is calculated using the shallow sea water depth inversion model.

[0012] Optionally, in step S1, the shallow sea water depth inversion model constructed is as follows:

[0013]

[0014] Among them, H is the final calculated shallow sea depth information; K wat is the sea surface reflection signal of optically shallow water; r is the remote sensing reflectivity of optically deep water; α is the first water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body; A0 and A1 are model fitting parameters; T1 and T2 are the average diffuse scattering attenuation coefficients related to the uplink irradiance and downlink irradiance, respectively; ρ is the bottom reflectivity.

[0015] Optionally, the optically shallow water surface reflection signal K wat , calculated from the remote sensing reflectance R of the upper surface of the sea:

[0016]

[0017] Among them, ζ is the water-air interface factor; τ is the internal reflection coefficient of the water-air interface. The above two parameters are derived based on the Hydrolight water optical radiation transmission model.

[0018] Optionally, the remote sensing reflectivity r of the optical deep water is calculated as follows:

[0019]

[0020] Among them, g0 and g1 are parameters related to the phase function of each component of seawater, which are determined by the water body type and observation geometry; a and b b They are the total absorption coefficient and backscattering coefficient of water bodies, respectively, which are provided by satellite remote sensing observation data.

[0021] Optionally, the first water quality parameter α is composed of the total absorption coefficient of the water body and the backscattering coefficient b b The specific calculation is as follows:

[0022] α=a+b b (4)

[0023] Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

[0024] Optionally, the average diffuse scattering attenuation coefficients T1 and T2 associated with the uplink irradiance and the downlink irradiance are calculated as follows:

[0025]

[0026] Among them, θ s is the solar zenith angle; D0, D1, D2, and D3 are empirical parameters, which are obtained by fitting based on the measured water depth data; μ is the second water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body, and is specifically calculated as follows:

[0027]

[0028] Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

[0029] Optionally, the bottom reflectivity ρ is calculated using the Hapke model, and the formula is as follows:

[0030]

[0031] Among them, function B represents the backscattering effect, function P HG represents the average scattering phase function of soil particles, function H represents multiple scattering phenomenon, S H is the hotspot peak value of backscattering, ω is the average single scattering albedo, g is the asymmetry factor, h is the hotspot width, ω and g are calculated based on Mie scattering theory; h is determined by the half-angle width at half of the hotspot peak; φ is the scattering angle, which is determined by the solar zenith angle θ in the satellite remote sensing observation data s , observation zenith angle θ v, relative azimuth The calculation is obtained as follows:

[0032]

[0033] Function B is expressed as follows:

[0034]

[0035] Function P HG It is expressed as follows:

[0036]

[0037] The function H is expressed as follows:

[0038]

[0039] On the other hand, a shallow sea water depth inversion system based on a multi-angle optical radiation transmission model is provided, which is used to implement any of the above methods, and the system comprises:

[0040] The model building module is used to build a shallow sea water depth inversion model based on the water body bio-optical model and the light radiation transmission process of the water body and the bottom;

[0041] The calculation module is used to calculate the shallow sea water depth information using the shallow sea water depth inversion model according to the acquired multi-angle satellite remote sensing observation data.

[0042] In another aspect, an electronic device is provided, the electronic device comprising:

[0043] processor;

[0044] A memory having computer-readable instructions stored thereon, wherein when the computer-readable instructions are loaded and executed by the processor, the steps of the shallow sea water depth inversion method based on the multi-angle optical radiation transmission model as described above are implemented.

[0045] On the other hand, a computer-readable storage medium is provided, in which program code is stored. The program code can be called by a processor to execute the steps of the shallow sea water depth inversion method based on the multi-angle optical radiation transmission model as described above.

[0046] The beneficial effects brought about by the technical solution provided by the present invention include at least:

[0047] The present invention integrates multi-angle satellite remote sensing observation data and physical models, and comprehensively considers the absorption, scattering and bottom reflection characteristics of water bodies to construct a highly robust and high-precision shallow water depth inversion method and system. The method has the ability to adapt to the complex optical environment of shallow water areas, can achieve large-scale, high-precision water depth inversion, and provides a new technical means for shallow water depth measurement, with broad application prospects. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0049] Figure 1 It is a flow chart of a shallow sea water depth inversion method based on a multi-angle optical radiation transmission model provided by an embodiment of the present invention;

[0050] Figure 2 Schematic diagram of the architecture of a shallow sea water depth inversion model based on multi-angle remote sensing data provided by an embodiment of the present invention;

[0051] Figure 3 is a schematic diagram of algorithm accuracy evaluation provided by an embodiment of the present invention;

[0052] Figure 4 It is a structural schematic diagram of a shallow sea water depth inversion system based on a multi-angle optical radiation transmission model provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0053] In order to make the purpose, technical solution and advantages of the embodiment of the present invention clearer, the technical solution of the embodiment of the present invention will be clearly and completely described below in conjunction with the drawings of the embodiment of the present invention. Obviously, the described embodiment is a part of the embodiment of the present invention, not all of the embodiments. Based on the described embodiment of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0054] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "exemplary" in the present invention should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the word "exemplary" is used to present concepts in a specific way.

[0055] The embodiment of the present invention provides a shallow sea water depth inversion method based on a multi-angle optical radiation transmission model, which can be implemented by an electronic device, which can be a terminal or a server. Figure 1 As shown, the processing flow of the method may include the following steps:

[0056] S1. Based on the water body bio-optical model and combined with the light radiation transmission process of the water body and the bottom, a shallow sea water depth inversion model is constructed;

[0057] S2. Based on the acquired multi-angle satellite remote sensing observation data, the shallow sea water depth information is calculated using the shallow sea water depth inversion model.

[0058] In step S1, the shallow sea water depth inversion model constructed is as follows:

[0059]

[0060] Among them, H is the final calculated shallow sea depth information; K wat is the sea surface reflection signal of optically shallow water; r is the remote sensing reflectivity of optically deep water; α is the first water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body; A0 and A1 are model fitting parameters; T1 and T2 are the average diffuse scattering attenuation coefficients related to the uplink irradiance and downlink irradiance, respectively; ρ is the bottom reflectivity.

[0061] Figure 2 It is a schematic diagram of the architecture of the shallow water depth inversion model. The solution process of each parameter is as follows:

[0062] The optical shallow water surface reflection signal K wat , calculated from the remote sensing reflectance R of the upper surface of the sea:

[0063]

[0064] Among them, ζ is the water-air interface factor; τ is the internal reflection coefficient of the water-air interface. The above two parameters are derived based on the Hydrolight water optical radiation transmission model.

[0065] The remote sensing reflectivity r of the optical deep water is calculated as follows:

[0066]

[0067] Among them, g0 and g1 are parameters related to the phase function of each component of seawater, which are determined by the water body type and observation geometry; a and b b They are the total absorption coefficient and backscattering coefficient of water bodies, respectively, which are provided by satellite remote sensing observation data.

[0068] The first water quality parameter α is composed of the total absorption coefficient of the water body and the backscattering coefficient bb The specific calculation is as follows:

[0069] α=a+b b (4)

[0070] Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

[0071] The calculation of the average diffuse scattering attenuation coefficients T1 and T2 related to the uplink irradiance and the downlink irradiance is as follows:

[0072]

[0073]

[0074] Among them, θ s is the solar zenith angle; D0, D1, D2, and D3 are empirical parameters, which are obtained by fitting based on the measured water depth data; μ is the second water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body, and is specifically calculated as follows:

[0075]

[0076] Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

[0077] The bottom reflectivity ρ is calculated using the Hapke model, and the formula is as follows:

[0078]

[0079] Among them, function B represents the backscattering effect, function P HG represents the average scattering phase function of soil particles, function H represents multiple scattering phenomenon, S H is the hotspot peak value of backscattering, ω is the average single scattering albedo, g is the asymmetry factor, h is the hotspot width, ω and g are calculated based on Mie scattering theory; h is determined by the half-angle width at half of the hotspot peak; φ is the scattering angle, which is determined by the solar zenith angle θ in the satellite remote sensing observation data s , observation zenith angle θ v , relative azimuth The calculation is obtained as follows:

[0080]

[0081] Among them, the relative azimuth is the solar azimuth With observation azimuth The absolute value of the difference;

[0082] Function B is expressed as follows:

[0083]

[0084] Function P HG It is expressed as follows:

[0085]

[0086] The function H is expressed as follows:

[0087]

[0088] In step S2, based on the above shallow sea depth inversion model and the observation geometry information provided by multi-angle satellite remote sensing observation data, the remote sensing reflectance R of the upper surface of the sea surface layer under the observation geometry information, the total absorption coefficient a of the water body and the backscattering coefficient b are calculated. b , and the empirical parameters fitted according to the measured data, the shallow sea depth information of the corresponding water area can be calculated.

[0089] Among them, multi-angle satellite remote sensing observation data can obtain water body and bottom reflection information from different angles, and capture the details of the optical properties of the water body and the bottom reflection characteristics at the same time, thereby significantly improving the accuracy of water depth inversion. Compared with single-angle observation, multi-angle observation technology can more comprehensively characterize the light transmission law of optical shallow water and reduce the interference of environmental noise.

[0090] In an embodiment of the present invention, the radiation transmission process of light in water is utilized to construct a multi-angle optical radiation transmission model that takes into account the bottom reflection signal. Based on the model, the transmission process of light in water is analyzed to achieve accurate inversion of shallow sea water depth.

[0091] In order to verify the effect of the method of the present invention, 50 MODIS (Moderate Resolution Imaging Spectroradiometer) pixels were randomly selected along the coast of China. The relevant offshore water depth was inverted based on MODIS remote sensing observation information and compared with the ocean water depth data. The results are as follows: Figure 3 It can be seen that the prediction result of the present invention has a high accuracy.

[0092] Accordingly, an embodiment of the present invention further provides a shallow sea water depth inversion system based on a multi-angle optical radiation transmission model, such as Figure 4 As shown, the system includes:

[0093] The model building module 201 is used to build a shallow sea water depth inversion model based on a water body bio-optical model and in combination with the light radiation transmission process of the water body and the bottom of the water;

[0094] The calculation module 202 is used to calculate the shallow sea water depth information using the shallow sea water depth inversion model according to the acquired multi-angle satellite remote sensing observation data.

[0095] For ease of explanation, Figure 4 Only the main components of the system are shown. The system of this embodiment can be used to perform Figure 1 The technical solution of the method embodiment shown has similar implementation principles and technical effects, which will not be repeated here.

[0096] The present invention integrates multi-angle satellite remote sensing observation data and physical models, and comprehensively considers the absorption, scattering and bottom reflection characteristics of water bodies to construct a highly robust and high-precision shallow water depth inversion method and system. The method has the ability to adapt to the complex optical environment of shallow water areas, can achieve large-scale, high-precision water depth inversion, and provides a new technical means for shallow water depth measurement, with broad application prospects.

[0097] In an exemplary embodiment, the present invention further provides an electronic device, the electronic device comprising:

[0098] processor;

[0099] A memory having computer-readable instructions stored thereon, wherein when the computer-readable instructions are loaded and executed by the processor, the steps of the shallow sea water depth inversion method based on the multi-angle optical radiation transmission model as described above are implemented.

[0100] In an exemplary embodiment, the present invention further provides a computer-readable storage medium, wherein at least one instruction is stored in the computer-readable storage medium, and the at least one instruction is loaded and executed by a processor to implement the steps of the shallow sea water depth inversion method based on the multi-angle optical radiation transmission model as described above. For example, the computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, etc.

[0101] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of further restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

[0102] References in the specification to "one embodiment", "an embodiment", "an exemplary embodiment", "some embodiments", etc. indicate that the described embodiments may include a particular feature, structure, or characteristic, but not every embodiment may include the particular feature, structure, or characteristic. In addition, when a particular feature, structure, or characteristic is described in conjunction with an embodiment, it should be within the knowledge of a person skilled in the relevant art to implement such feature, structure, or characteristic in conjunction with other embodiments (whether or not explicitly described).

[0103] It should be understood that the term "and / or" in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. A and B can be singular or plural. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship, but it may also indicate an "and / or" relationship. Please refer to the context for specific understanding.

[0104] In the present invention, "at least one" means one or more, and "more than one" means two or more. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, at least one of a, b, or c can be represented by: a, b, c, ab, ac, bc, or abc, where a, b, c can be single or multiple.

[0105] It should be understood that in various embodiments of the present invention, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0106] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0107] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0108] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0109] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0110] The present invention covers any substitution, modification, equivalent method and scheme made on the essence and scope of the present invention. In order to make the public have a thorough understanding of the present invention, specific details are described in detail in the following preferred embodiments of the present invention, but those skilled in the art can fully understand the present invention without the description of these details. In addition, in order to avoid unnecessary confusion about the essence of the present invention, well-known methods, processes, procedures, components and circuits are not described in detail.

[0111] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A shallow sea water depth inversion method based on a multi-angle optical radiation transmission model, characterized in that: The following steps are involved: S1. Based on the water body bio-optical model and combined with the light radiation transmission process of the water body and the bottom, a shallow sea water depth inversion model is constructed; S2. Based on the acquired multi-angle satellite remote sensing observation data, the shallow sea water depth information is calculated using the shallow sea water depth inversion model.

2. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 1 is characterized in that: In step S1, the shallow sea water depth inversion model constructed is as follows: Among them, H is the final calculated shallow sea depth information; K wat is the sea surface reflection signal of optically shallow water; r is the remote sensing reflectivity of optically deep water; α is the first water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body; A0 and A1 are model fitting parameters; T1 and T2 are the average diffuse scattering attenuation coefficients related to the uplink irradiance and downlink irradiance, respectively; ρ is the bottom reflectivity.

3. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 2 is characterized in that: The optical shallow water surface reflection signal K wat , calculated from the remote sensing reflectance R of the upper surface of the sea: Among them, ζ is the water-air interface factor; τ is the internal reflection coefficient of the water-air interface. The above two parameters are derived based on the Hydrolight water optical radiation transmission model.

4. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 2 is characterized in that: The remote sensing reflectivity r of the optical deep water is calculated as follows: Among them, g0 and g1 are parameters related to the phase function of each component of seawater, which are determined by the water body type and observation geometry; a and b b They are the total absorption coefficient and backscattering coefficient of water bodies, respectively, which are provided by satellite remote sensing observation data.

5. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 3 is characterized in that: The first water quality parameter α is composed of the total absorption coefficient of the water body and the backscattering coefficient b b The specific calculation is as follows: α=a+b b (4) Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

6. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 2 is characterized in that: The calculation of the average diffuse scattering attenuation coefficients T1 and T2 related to the uplink irradiance and the downlink irradiance is as follows: Among them, θ s is the solar zenith angle; D0, D1, D2, and D3 are empirical parameters, which are obtained by fitting based on the measured water depth data; μ is the second water quality parameter, which is determined by the total absorption coefficient and backscattering coefficient of the water body, and is specifically calculated as follows: Among them, a is the total absorption coefficient of water, b is b is the backscattering coefficient.

7. The shallow sea water depth inversion method based on the multi-angle optical radiation transmission model according to claim 2 is characterized in that: The bottom reflectivity ρ is calculated using the Hapke model, and the formula is as follows: Among them, function B represents the backscattering effect, function P HG represents the average scattering phase function of soil particles, function H represents multiple scattering phenomenon, S H is the hotspot peak value of backscattering, ω is the average single scattering albedo, g is the asymmetry factor, h is the hotspot width, ω and g are calculated based on Mie scattering theory; h is determined by the half-angle width at half of the hotspot peak; φ is the scattering angle, which is determined by the solar zenith angle θ in the satellite remote sensing observation data s , observation zenith angle θ v , relative azimuth The calculation is obtained as follows: Function B is expressed as follows: Function P HG It is expressed as follows: The function H is expressed as follows:

8. A shallow sea water depth inversion system based on a multi-angle optical radiation transmission model, the system being used to implement the method according to any one of claims 1 to 7, characterized in that: The system comprises: The model building module is used to build a shallow sea water depth inversion model based on the water body bio-optical model and the light radiation transmission process of the water body and the bottom; The calculation module is used to calculate the shallow sea water depth information using the shallow sea water depth inversion model according to the acquired multi-angle satellite remote sensing observation data.

9. An electronic device, characterized in that: The electronic device comprises: processor; A memory having computer-readable instructions stored thereon, wherein when the computer-readable instructions are loaded and executed by the processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores program code, which can be called by a processor to execute the method according to any one of claims 1 to 7.

Images