A method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data

Abstract

The present invention discloses a method for inverting the vertical structure of the optical properties of the ocean subsurface based on ICESat‑2 ATL03 global positioning photon data. The method first performs quality control on the ICESat‑2 ATL03 original data set; then, the sea surface position is determined by screening high-confidence photons, and the depth of the underwater photons is further determined; then, the obtained depth data is refraction corrected. Using the processed ICESat‑2 photon data, a sliding window is used, and photons at different depths are slidingly accumulated with a set step size to obtain the depth distribution of the laser radar attenuated backscatter signal; then, the laser radar attenuated backscatter signal is deconvoluted; finally, the diffuse attenuation coefficient K of the optical properties of the ocean subsurface is obtained using a physical inversion model. d and particle backscattering coefficient b bp The present invention can realize vertical observation of the ocean subsurface even in areas where it is difficult to carry out on-site measurement, and can ensure high measurement accuracy.

Description

Technical Field

[0001] The present invention belongs to the technical field of ocean laser radar, and in particular to a method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 (Ice, Cloud and land Elevation Satellite-2) ATL03 global positioning photon data. Background Art

[0002] Observations of ocean subsurface optical parameters are crucial to revealing the dynamic changes of marine ecosystems and biogeochemical cycles. By deeply analyzing these parameters, we can infer phytoplankton biomass and particulate organic carbon content, thereby assessing marine primary productivity. This information is of great value in global ocean monitoring, especially in water quality assessment, biodiversity conservation, species distribution analysis, and monitoring of environmental issues such as algal blooms. However, obtaining the vertical structure of large-scale ocean subsurface optical properties still faces great challenges. Passive optical remote sensing technology is limited by insufficient vertical resolution, while in situ observations are difficult to achieve large-scale coverage due to high costs and human resource constraints.

[0003] Advances in lidar sensors have made it an innovative tool for observing the optical structure of water bodies, mainly including airborne, shipborne and satellite-borne lidar technologies. Airborne and shipborne lidars can efficiently detect the optical properties of the ocean subsurface and provide high-resolution data, but their high cost limits their application in large-scale and long-term monitoring. In comparison, satellite-borne lidars overcome these limitations. Although Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) extracts the optical properties of water bodies by preprocessing echo data, it faces certain challenges in processing weak signals. Several common inversion methods currently show good results in clear water bodies, but their application effects are limited in nearshore turbid waters. Summary of the invention

[0004] In order to invert the vertical profile of the ocean subsurface optical properties, the present invention provides a method for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data. The specific technical scheme is as follows:

[0005] A method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data, the method comprising the following steps:

[0006] S1: Perform quality control on the ATL03 data of ICESat-2, and only select photons with a quality label quality_ph of 0 in ATL03 of ICESat-2;

[0007] S2: Based on the ICESat-2 ATL03 photon confidence parameter signal_conf_ph, calculate the sea surface elevation h S3: Based on the ICESat-2 ATL03 parameter h_ph, obtain the photon ellipsoid height height, and subtract the sea surface elevation h from it to obtain the photon depth data z:

[0008]

[0009] S4: Refraction correction for the depth of photons:

[0010]

[0011] S5: Perform sliding accumulation with the set window and sliding step size to obtain the laser radar attenuated backscattering signal S at different depths m (z);

[0012] S6: The attenuated backscatter signal S obtained in S5 m (z) performs a deconvolution operation to obtain the correct attenuated backscatter signal S(z);

[0013] S7: Inversion of diffuse attenuation coefficient K d vertical section structure;

[0014] S7.1: Perform distance and logarithm calculations on the attenuated backscatter signal S(z) obtained in S6 to obtain S'(z); for water bodies with a depth of 3-10m, use the slope method to calculate the laser radar attenuation constant α. The specific calculation formula is as follows:

[0015]

[0016]

[0017] S7.2: Use Klett's method to obtain the vertical profile α(z) of the lidar attenuation coefficient;

[0018] S7.3: Use the vertical profile α(z) of the lidar attenuation coefficient as the diffuse attenuation coefficient K d (z) vertical section structure;

[0019] S8: Obtain particle backscattering coefficient b bp vertical structure.

[0020] Furthermore, the specific steps of obtaining the sea surface elevation position in S2 are:

[0021] The data with signal_conf_ph of 4 in the ATL03 data of ICESat-2 are taken out, the average is calculated for each group, and the nearest neighbor interpolation is used to quickly obtain the sea surface height data at different locations.

[0022] Furthermore, the S6 specifically includes the following sub-steps:

[0023] S6.1: Calculate the transient response function F(z) by the following equation:

[0024]

[0025] in,

[0026]

[0027] Where i and j represent the range unit numbers, and p represents the range unit number where the signal peak returns;

[0028] S6.2: Calculate the correct attenuated backscatter signal S(z) by the following formula:

[0029]

[0030] Among them, F -1 (z) is the inverse matrix of F(z).

[0031] Furthermore, the S7.2 specifically includes the following sub-steps:

[0032] (1) Assuming the 180° backscattering coefficient of the water body and the lidar attenuation coefficient The proportional relationship is as follows:

[0033]

[0034] Where const is a constant, k is a power law exponent, which depends on the wavelength of the lidar and the properties of the seawater, and is between 0.67 and 1;

[0035] (2) Perform distance correction on S(z) to obtain the distance-corrected attenuated backscatter signal RCS(z):

[0036]

[0037] (3) Calculate the laser radar attenuation coefficient using the following formula: :

[0038]

[0039]

[0040] Among them, zm represents the maximum effective depth of laser detection, α(z m ) is the reference z m The attenuation coefficient at is obtained by the slope method.

[0041] Furthermore, the S8 specifically includes the following sub-steps:

[0042] (1) Calculation of theoretical ocean surface backscatter :

[0043]

[0044] in, is the incident angle of the LiDAR, is the average square wave slope of the ocean surface estimated from the wind speed and is calculated as follows:

[0045]

[0046] (2) Calculate the calibration constant A:

[0047]

[0048] Where S'(0) represents the attenuated backscatter value of ICESat-2 at the sea surface after range and logarithmic correction;

[0049] (3) Finally, obtain the vertical profile b of the particle backscattering coefficient bp (z):

[0050]

[0051] Where n0 is the correction factor to compensate for the difference between the lidar angle and the solar zenith angle; β (π) is the backscattering phase function; and t is the transmittance of the sea-air interface.

[0052] Furthermore, the sizes of the window and sliding step set in S5 are: 4 km in the horizontal direction, 1 m in the vertical direction, and the vertical sliding step is 0.15 m.

[0053] A device for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data comprises one or more processors for implementing a method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data.

[0054] An electronic device, comprising:

[0055] one or more processors;

[0056] A storage device is used to store one or more programs. When the one or more programs are executed by the electronic device, the electronic device implements a method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data.

[0057] A computer-readable storage medium stores a program, which, when executed by a processor, implements a method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data.

[0058] The beneficial effects of the present invention are as follows:

[0059] Compared with passive optical remote sensing, the present invention uses the ATLAS (Advanced Topographic Laser Altimeter System) laser radar system carried by the ICESat-2 satellite to obtain the optical property profile of the ocean subsurface. This method has achieved significant improvements in observation range and data accuracy, and the results are highly consistent with existing data.

[0060] Compared with other satellite-borne lidar systems, the present invention utilizes the ATLAS lidar system carried by the ICESat-2 satellite, based on single-photon detection technology, which greatly enhances the system's sensitivity to weak signals, thereby reducing the energy demand of the laser and improving data quality and inversion accuracy, and providing more detailed vertical structure information.

[0061] In view of the problems that in-situ measured data or airborne or shipborne data often face limited coverage and difficulty in long-term monitoring, the method of the present invention has a wide detection range and supports long-term monitoring of global waters, including clear water bodies and near-shore turbid waters. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] Figure 1 It is a schematic flow chart of the inversion method of the present invention.

[0063] Figure 2 is b in S8 of the present invention bp Schematic diagram of profile inversion processing.

[0064] Figure 3 is a comparison chart of the results of using the method of the present invention and the existing method, wherein Figure 3 Figure A and Figure 3 Figure B in the figure compares the diffuse attenuation coefficient K obtained in this embodiment of the present invention. d and particle backscattering coefficient b bpResults from MODIS (moderate-resolution imaging spectroradiometer) products; Figure 3 Figure C and Figure 3 The D diagram in FIG. 1 shows the K d and b bp Vertical section structure.

[0065] Figure 4 The in-situ measurement data of a BGC_Argo float and the diffuse attenuation coefficient K inverted from ICESat-2 are shown. d and particle backscattering coefficient b bp The comparison results.

[0066] Figure 5 The in-situ measurement data of another BGC_Argo float and the diffuse attenuation coefficient K obtained by ICESat-2 inversion are shown. d and particle backscattering coefficient b bp The comparison results. DETAILED DESCRIPTION

[0067] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0068] The terms used in this application are for the purpose of describing specific embodiments only and are not intended to limit this application. The singular forms of "a", "said" and "the" used in this application and the appended claims are also intended to include plural forms unless the context clearly indicates other meanings. It should also be understood that the term "and / or" used herein refers to and includes any or all possible combinations of one or more associated listed items.

[0069] It should be understood that although the terms first, second, third, etc. may be used in the present application to describe various information, these information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of the present application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

[0070] It should also be noted that the "multiple" mentioned in this application refers to two or more than two. "And / or" describes the association relationship of the associated objects, indicating that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the previous and next associated objects are in an "or" relationship.

[0071] like Figure 1 As shown, on the one hand, the method for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data of this embodiment includes the following steps:

[0072] S1: Perform quality control on the ATL03 data of ICESat-2, and only select photons with a quality label quality_ph of 0 in ATL03 of ICESat-2;

[0073] S2: Based on the ICESat-2 ATL03 photon confidence parameter signal_conf_ph, the sea surface elevation position h is calculated.

[0074] In this embodiment, data with signal_conf_ph of 4 in ATL03 data of ICESat-2 are taken out, grouped and averaged, and nearest neighbor interpolation is used to quickly obtain sea surface height data at different locations.

[0075] S3: Based on the ATL03 parameter h_ph of ICESat-2, obtain the ellipsoid height of the photon, and subtract the sea surface elevation h from it to obtain the depth data z of the photon:

[0076]

[0077] S4: Refraction correction for the depth of photons:

[0078]

[0079] S5: Perform sliding accumulation with the set window and sliding step size to obtain the attenuated backscattering signal S at different depths m (z).

[0080] In this embodiment, the size of the window and the sliding step is set to 4 km in the horizontal direction and 1 m in the vertical direction, and the vertical sliding step is 0.15 m.

[0081] S6: The attenuated backscatter signal S obtained in S5 m (z) performs a deconvolution operation to obtain a correct attenuated backscatter signal S(z). S6 specifically includes the following sub-steps:

[0082] First, calculate the transient response function F(z) according to the following formula:

[0083]

[0084] in,

[0085]

[0086] The method uses the ICESat-2 lidar signal from 12 adjacent range cells to calculate the transient response function, starting from the cell before the signal peak return and extending to the tenth cell after the peak. j In the attenuated backscatter value at , i and j represent the range unit number, and p represents the range unit number where the signal peak returns. After deconvolution, the correct lidar attenuated backscatter signal is restored.

[0087]

[0088] where F -1 (z) is the inverse matrix of F(z).

[0089] S7: Inversion of diffuse attenuation coefficient K d Vertical section structure, including:

[0090] S7.1: Perform distance and logarithm calculations on the attenuated backscatter signal S(z) obtained in S6 to obtain S'(z): Then, for a water body with a depth of 3-10m, calculate the laser radar attenuation constant α using the slope method. The specific calculation formula is as follows:

[0091]

[0092]

[0093] S7.2: Obtain the vertical profile α(z) of the lidar attenuation coefficient using the Klett method, including:

[0094] (1) Assuming the 180° backscattering coefficient of water and the lidar attenuation coefficient The proportional relationship is as follows:

[0095]

[0096] Wherein, const is a constant; k is a power law exponent, the magnitude of which depends on the laser radar wavelength and the properties of seawater, and is between 0.67 and 1. In this embodiment, for the convenience of calculation, it is assumed that k=1.0.

[0097] (2) Perform distance correction on S(z) to obtain the distance-corrected attenuated backscatter signal RCS(z):

[0098]

[0099] (3) Calculate the laser radar attenuation coefficient using the following formula: :

[0100]

[0101] Among them, α(z m ) is the reference boundary depth z m The attenuation coefficient at can be obtained by the slope method. This method combines the forward and reverse inversion models and significantly reduces the impact of boundary values ​​and noise on the results. m Represents the maximum effective depth of laser detection. For ICESat-2, due to its low single pulse energy, its maximum effective detection depth on low-reflectivity targets such as the ocean is roughly equivalent to an optical thickness, expressed as:

[0102]

[0103] To enhance the inversion process, a reference depth rule is established to limit the reference depth to a maximum of 15 m. By adopting an optical thickness as the reference depth, the attenuation coefficient at different depths in seawater can be accurately inverted.

[0104] S7.3: Calculate the diffuse attenuation coefficient K d (z), for ICESat-2, K d (z) is approximately equal to the laser radar attenuation coefficient α(z), so the vertical profile α(z) of the laser radar attenuation coefficient is used as the diffuse attenuation coefficient K d (z).

[0105] S8: Obtain particle backscattering coefficient b bp The vertical structure, such as Figure 2 As shown, the specific steps include:

[0106] (1) First, calculate the theoretical ocean surface backscattering

[0107]

[0108] in, is the incident angle of the LiDAR, is the average square wave slope of the ocean surface estimated from the wind speed and is calculated as follows:

[0109]

[0110] (2) Then calculate the calibration constant A:

[0111]

[0112] Where S'(0) represents the attenuated backscatter value of ICESat-2 at the sea surface after range and logarithmic correction;

[0113] (3) Finally, obtain the vertical profile b of the particle backscattering coefficient bp (z):

[0114]

[0115] Where n0 is a correction factor to compensate for the difference between the laser radar angle and the solar zenith angle. In the embodiment, 1.4 is used as the compensation term. β (π) is the backscattering phase function, which is 0.32 sr -1 t is the transmittance of the sea-air interface, which is about 0.98.

[0116] On the other hand, corresponding to the aforementioned embodiment of the method for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data, the present invention also provides an embodiment of a device for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data.

[0117] The device for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data provided in an embodiment of the present invention includes one or more processors for implementing the method for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data in the above-mentioned embodiment.

[0118] The vertical structure inversion device for the optical properties of the ocean subsurface layer based on the ICESat-2 ATL03 global positioning photon data of the present invention can be applied to any device with data processing capabilities, and the any device with data processing capabilities can be a device or apparatus such as a computer. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. Taking software implementation as an example, as a device in a logical sense, it is formed by the processor of any device with data processing capabilities in which it is located to read the corresponding computer program instructions in the non-volatile memory into the memory for operation. From the hardware level, in addition to the processor, memory, network interface, and non-volatile memory, the vertical structure inversion device for the optical properties of the ocean subsurface layer based on the ICESat-2 ATL03 global positioning photon data of the present invention, any device with data processing capabilities in which the device is located in the embodiment can also include other hardware according to the actual function of the device with data processing capabilities, which will not be described in detail.

[0119] For the device embodiment, since it basically corresponds to the method embodiment, the relevant parts can refer to the partial description of the method embodiment. The device embodiment described above is only schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of the present invention. Ordinary technicians in this field can understand and implement it without paying creative work.

[0120] An embodiment of the present invention further provides a computer-readable storage medium having a program stored thereon. When the program is executed by a processor, the method for inverting the vertical structure of the optical properties of the ocean subsurface layer based on ICESat-2 ATL03 global positioning photon data in the above embodiment is implemented.

[0121] The computer-readable storage medium may be an internal storage unit of any device with data processing capability described in any of the aforementioned embodiments, such as a hard disk or a memory. The computer-readable storage medium may also be an external storage device, such as a plug-in hard disk, a smart memory card (SmartMedia card, SMC), an SD card, a flash card, etc. equipped on the device. Furthermore, the computer-readable storage medium may also include both an internal storage unit of any device with data processing capability and an external storage device. The computer-readable storage medium is used to store the computer program and other programs and data required by any device with data processing capability, and may also be used to temporarily store data that has been output or is to be output.

[0122] The effect of the method of the present invention is verified by a specific example below.

[0123] This example uses the Pacific nighttime lidar trajectory data on January 22, 2020, which comes from the ATL03 file of the ICESat-2 satellite (ATL03_20200122201210_04140606_006_01.h5). To verify the accuracy of the model results, this example uses the monthly average K d and b bp The corresponding files are AQUA_MODIS.20200101_20200131.L3m.MO.KD.Kd_490.9km.nc, AQUA_MODIS.20200101_20200131.L3m.MO.IOP.bbp_443.9km.nc, SR2902754_118.nc and SR2902754_119.nc. It should be noted that it is difficult to obtain buoy data that completely matches the time and spatial position of ICESat-2 observations. Therefore, this embodiment uses buoy data that is adjacent to the ICESat-2 observation position and close in time. In order to compare and verify with the 532 nm band data of ICESat-2, the Ed490 and bbp700 in the in-situ measured data were converted into the corresponding K values ​​of the 532 nm band after a series of conversions. d and b bp This process ensures the comparability of data from different bands, thereby improving the reliability of the verification results. Figure 3 Figure A and Figure 3 As shown in Figure B, the optical properties inverted by the method of the present invention are slightly different from the MODIS results. Figure 3 Figure C and Figure 3 Figure D in the figure shows the diffuse attenuation coefficient K of the laser scanning area inverted by the present invention. d and particle backscattering coefficient b bp Vertical cross-section structure. Figure 4 and Figure 5 The comparison results of the two in situ measurement data and the vertical profiles of the optical properties of the corresponding points of ICESat-2 are shown, indicating that the two are in good agreement with very small differences.

[0124] Those skilled in the art can understand that the above are only preferred examples of the invention and are not intended to limit the invention. Although the invention is described in detail with reference to the above examples, those skilled in the art can still modify the technical solutions recorded in the above examples or replace some of the technical features therein with equivalents. Any modification, equivalent replacement, etc. made within the spirit and principle of the invention shall be included in the protection scope of the invention.

Claims

1. A method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data, characterized in that: The method comprises the following steps: S1: Perform quality control on the ATL03 data of ICESat-2, and only select photons with a quality label of 0 in ATL03 of ICESat-2; S2: Based on the ICESat-2 ATL03 photon confidence parameter signal_conf_ph, the sea surface elevation position h is calculated; S3: Based on the ATL03 parameter h_ph of ICESat-2, obtain the ellipsoid height of the photon, and subtract the sea surface elevation position h from it to obtain the depth data z of the photon: z = height - h; S4: Refraction correction for the depth of photons: z=z / 1.333; S5: Perform sliding accumulation with the set window and sliding step size to obtain the laser radar attenuated backscattering signal S at different depths m (z); S6: The attenuated backscatter signal S obtained in S5 m (z) performs a deconvolution operation to obtain the correct attenuated backscatter signal S(z); S7: Inversion of diffuse attenuation coefficient K d vertical section structure; S7.1: Perform distance and logarithm calculations on the attenuated backscatter signal S(z) obtained in S6 to obtain S'(z); for water bodies with a depth of 3-10m, the slope method is used to calculate the laser radar attenuation constant α. The specific calculation formula is as follows: S′(z)=ln[S(z)·(H+z) 2 ]=lnAβ(π,z)-2α(z)·z; S7.2: Obtain the vertical profile α(z) of the lidar attenuation coefficient using the Klett method; S7.3: Use the vertical profile α(z) of the lidar attenuation coefficient as the diffuse attenuation coefficient K d (z) vertical section structure; S8: Obtain particle backscattering coefficient b bp The vertical structure of The S6 specifically includes the following sub-steps: S6.1: Calculate the transient response function F(z) by the following equation: in, Where i and j represent the range unit numbers, and p represents the range unit number where the signal peak returns; S6.2: Calculate the correct attenuated backscatter signal S(z) by the following formula: S(z)=F -1 (z)*S m (z); Among them, F -1 (z) is the inverse matrix of F(z); The S7.2 specifically includes the following sub-steps: (1) Assume that the 180° backscattering coefficient β(π,z) of the water body and the lidar attenuation coefficient α(z) have the following proportional relationship: β(π,z)=const×α(z) k ; Where const is a constant, k is a power law exponent, which depends on the wavelength of the lidar and the properties of the seawater, and is between 0.67 and 1; (2) Perform distance correction on S(z) to obtain the distance-corrected attenuated backscatter signal RCS(z): RCS(z)=S(z)·(H+z) 2 ; (3) The laser radar attenuation coefficient α(z) is calculated by the following formula: z m ＝1 / α； Among them, z m represents the maximum effective depth of laser detection, α(z m ) is the reference z m The attenuation coefficient at is obtained by the slope method; The S8 specifically includes the following sub-steps: (1) Calculation of theoretical ocean surface backscatter β s : Where θ is the incident angle of the laser radar, σ 2 is the average square wave slope of the ocean surface estimated from the wind speed and is calculated as follows: (2) Calculate the calibration constant A: Where S'(0) represents the attenuated backscatter value of ICESat-2 at the sea surface after range and logarithmic correction; (3) Finally, obtain the vertical profile b of the particle backscattering coefficient bp (z): Where n0 is the correction factor to compensate for the difference between the lidar angle and the solar zenith angle; β(π) is the backscattering phase function; and t is the transmittance of the sea-air interface.

2. The method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data according to claim 1, characterized in that: The specific steps of obtaining the sea surface elevation position in S2 are: The data with signal_conf_ph of 4 in the ATL03 data of ICESat-2 are taken out, the average is calculated for each group, and the nearest neighbor interpolation is used to quickly obtain the sea surface height data at different locations.

3. The method for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data according to claim 1, characterized in that: The sizes of the window and sliding step set in S5 are: 4 km in the horizontal direction, 1 m in the vertical direction, and the vertical sliding step is 0.15 m.

4. A device for inverting the vertical structure of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data, characterized in that: The method comprises one or more processors for implementing the method for inverting the vertical structure of the optical properties of the ocean subsurface layer based on ICESat-2 ATL03 global positioning photon data as described in any one of claims 1 to 3.

5. An electronic device, characterized in that: include: one or more processors; A storage device for storing one or more programs. When the one or more programs are executed by the electronic device, the electronic device implements the method for vertical structure inversion of ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data as described in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that: A program is stored thereon, and when the program is executed by a processor, the method for inverting the vertical structure of the ocean subsurface optical properties based on ICESat-2 ATL03 global positioning photon data as described in any one of claims 1 to 3 is implemented.

Images