This application is based upon and claims priority to Chinese Patent Application No. 202211204456.0, filed on Sep. 29, 2022, the entire contents of which are incorporated herein by reference.
The present invention belongs to the technical field of water body transmittance data identification, and especially relates to the method and application of flash imaging lidar to measure the attenuation coefficient of water bodies.
Water transmittance is a common index to evaluate the water environment. An important parameter to quantify the water transmittance is the attenuation coefficient of the water body, and at the same time, the attenuation coefficient of the water body is an important parameter for underwater detection by LIDAR systems, which characterizes the basic properties of the water body and determines the spectral characteristics of the light field transmitted in the water. Therefore, it is particularly important to measure the attenuation coefficient of the water body accurately and effectively in real time.
The directly method to measure the attenuation coefficient of water bodies is Beer-Lambert Law. Other methods to measure the attenuation coefficient of water bodies such as by measuring its Raman scattering, through the stimulated Brillouin scattering, and measuring its apparent optical properties. Although the Beer-Lambert Law is more accurate within a certain range, it requires the sampling of test water which makes it difficult to be performed remotely.
Through the above analysis, the problems and defects of the prior technology are:
To overcome the problems in related technologies, the disclosed embodiments of the present invention provide methods and applications of flash imaging lidar for measuring the attenuation coefficient of water bodies.
The technical solution is as follows: The flash imaging lidar uses the AFD method to measure the backscattering intensity values along the front and back edges of the range-gate in the water body, and the model of the relationship between the attenuation coefficient of the water body and the backscattering intensity is combined to obtain the attenuation coefficient of a location in the water body. The following steps are included:
In step 1, The equation for the total backscattering intensity in the water column is:
Where AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient, e≈2.7182.
In step 2, the horizontal distance of the target from the LIDAR system is solved to obtain the differential equation as:
Where dpb is the backscattering intensity, AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient.
In step 3, The relationship between the attenuation coefficient of the water body and the backscattering intensity is modeled as:
Where dpb is the backscattering intensity, AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient.
In one embodiment, the backscattering intensity values in a water body at the location of the front and back edges of the range-gate measured using the AFD method include:
The Δτgate is the width of the range-gate, Δτstop is the delay step size of the flash lidar system, the imaging distance is large than R1 which indicates F=1, the intensity image received by the ICCD is the backscattering intensity image of the full laser beam in the water body:
By assuming that the adjacent frames of the backscattering intensity images received by the lidar are kth frame and k+1th frame, the corresponded intensity difference image ΔIk can be acquired. The pixel intensity values of both front and back edge parts in ΔIk as the backscattering intensity at the corresponding distance. Assuming the backscattering intensities at the back and front edges of the range-gate are dPr
dP
r
=ΣΔI
(i,j)
k
, dP
r
=ΣΔI
(m,n)
k (9)
In an instance, obtaining the attenuation coefficient at a location in a water body includes: According to the Eq.(9) and the Eq.(8) we can get:
where r1=v(Δτ+kΔτstep), r2=v[Δτ+(k+1)Δτstep], v is the light speed in the water, Δτ is the initial delay of the lidar system. C is the attenuation coefficient at a certain location in the water body.
A further object of the present invention is to provide a computer device. The computer device which is include a memory and a processor. Memory stores computer program. When the computer program executed by the processor, causes processor to perform a method for implementing the flash imaging lidar measurement of the attenuation coefficient of a water body as claimed in any one of claims 1-7.
A further object of the present invention is to provide a computer readable storage media. The computer-readable storage media stores program, when the program executed by a processor, causes the processor to perform a method of implementing any one of claims 1-7 for flash imaging lidar measurement of attenuation coefficients in a water body.
Another purpose of the invention is to provide an information data processing terminal which, when implemented on an electronic device, provides a user input interface to implement the flash imaging lidar method for measuring backward attenuation coefficient of water bodies.
Another object of the present invention is to provide a method for the application of said flash imaging lidar to measure the attenuation coefficient of a water body in liquid measurements.
Combining all of the above technical solutions, the advantages and positive effects of the present invention are: Based on the backscattering model of water bodies, this invention realizes the real-time remote sensing measurement of the attenuation coefficient of water bodies using the adjacent frame difference method. The relationship between the attenuation coefficient and the backscattering intensity of the water body was first established based on a model of the total backscattering intensity in the water body. Subsequently, the adjacent frame difference method proposed according to the present invention can measure the backscattering intensity at a location in the water body in real time, and thus obtain the attenuation coefficient of the water body. The experimental and theoretical analyses fully illustrate the feasibility of real-time remote sensing measurement of the attenuation coefficient of water bodies using the adjacent frame difference method, which has a strong practical value.
The accompanying figures herein are incorporated into and form part of the specification, illustrate embodiments consistent with the present disclosure, and are used in conjunction with the specification to explain the principles of the present disclosure.
In order to make the above-mentioned objects, features and advantages of the present invention more obvious and understandable, the following detailed description of specific embodiments of the present invention is made in conjunction with the accompanying figures. Many specific details are set forth in the following description to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways than those described herein, and similar improvements can be made by those skilled in the art without contradicting the content of the present invention, so that the present invention is not limited by the specific embodiments disclosed below.
Embodiments of the present invention provide a method for measuring the backward attenuation coefficient of a water body by flash imaging lidar including the steps of:
The attenuation coefficient of a location in the water body will be obtained by combining the model of the relationship between the attenuation coefficient of the water body and the backscattering intensity along the front and back edges of the range-gate measured using the adjacent frame difference method.
A method for measuring the attenuation coefficient of a water body by flash imaging lidar according to the method documented in Example 1, in a preferred embodiment, as shown in
A method for measuring the backward attenuation coefficient of a water body by flash imaging lidar as documented in Example 2. Where α is the half divergence angle of the laser, γ is the half field of view of the receiver, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, d0 is the distance of the transmitting center and receiving center, A(r) is the irradiated area of the laser at distance of r, C is the attenuation coefficient of water bodies. The total backscattering luminous flux from the wale body can be yield as:
Where U is the backward scattering coefficient, AR is the effective receiving area.
The geometry sizes of the laser spot and received field-of-view interface varied with the distances between the laser source and target, which results in the changes on their overlap area. In order to analyze the backscattering with different overlap area between the laser spot and the received field of view, an overlap coefficient F is introduced, which is defined as the ratio of the laser spot area captured in the received field of view to the total laser irradiation area. Can be expressed as:
the angles corresponding to the overlapping area sections on the cross sections of the laser spot and received field of view are defined as ω1 and 1072, respectively. And r1 and r2 are the radius of the laser spot and the receiving field-of-view, respectively. d is the distance between the centers of laser spot and receiving field-of-view. The radii of the transmitting optical system and the receiving optical system are rL and rR, respectively. Then the detailed expression for the overlapping area at different positions can be calculated from the trigonometric relationship:
The introduction of F reduces the dimension of the integral of Eq.11 and turns the double integral into an ordinary integral to simplify calculations. According to Eq.11 and Eq.12, the total backscattering flux of the water body can be further expressed as:
Where AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient, e≈2.7182.
In S102, the horizontal distance of the target from the LIDAR system is solved to obtain the differential equation as:
Where dpb is the backscattering intensity, AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient.
In S103, the relationship between the attenuation coefficient of the water body and the backscattering intensity is modeled as:
Where dpb is the backscattering intensity, AR is the effective receiving area, P0 is the output laser optical power, half of the angle between the transmitting center axis and the receiving center axis is the θ, r is the distance of the target from the Lidar system, C is the attenuation coefficient of water bodies, U is the backscattering coefficient, F is the overlap coefficient. only r and
are unknown value. Therefore, the attenuation coefficient of the water body can be calculated as long as the
and corresponding distance r been measured.
As known to those skilled in the art, the method of measuring the attenuation coefficient of a water body by the flash imaging lidar provided by embodiments of the present invention first establishes the relationship between the attenuation coefficient of the water body and the backscattering intensity value measured in real time by the adjacent frame difference (AFD) method based on the backscatter model in the water.
Essentially, both the water body backscattering light and the lidar target return light are reflected signals from different objects in the water to the incident light. Therefore, based on the principle of underwater light transmission and reflection provided by the lidar equation, a computational model of backscattering light from the water body can be developed.
The AFD method is a new three-dimensional (3D) reconstruction algorithm based on the range-gated method. In this method, the distance information of the target at the front and back edges of the range-gated are obtained by thresholding the intensity difference between two adjacent frames achieved by the intensified-charge-coupled device (ICCD), respectively. According to the characteristics of the AFD method, it can be known that the AFD method can obtain the backscattering intensity corresponding to the front and back edges of the range-gate in real time.
The Δτgate is the width of the range-gate, Δτstep is the delay step size of the flash lidar system, the imaging distance is large than R1 which indicates F=1, the intensity image received by the ICCD is the backscattering intensity image of the full laser beam in the water body.
By assuming that the adjacent frames of the backscattering intensity images received by the lidar are kth frame and k+1th frame, the corresponded intensity difference image ΔIk can be acquired. The pixel intensity values of both front and back edge parts in ΔIk as the backscattering intensity at the corresponding distance. Assuming the backscattering intensities at the back and front edges of the range-gate are dPr
dP
r
=ΣΔI
(i,j)
k
, dP
r
=ΣΔI
(m,n)
k (17)
Obtaining the attenuation coefficient at a location in the water body includes: based on the Eq. (16) and Eq. (17) we can get:
where r1=v(Δτ+kΔτstep), r2=v[Δτ+(k+1)Δτstep], v is the light speed in the water, Δτ is the initial delay of the lidar system. C is the attenuation coefficient at a certain location in the water body.
Embodiments of the present invention also provide a computer device including: at least one processor, a memory, and a computer program stored in said memory and runnable on the at least one processor, the processor executing said computer program to implement the steps in any of the method embodiments described above.
Embodiments of the present invention also provide a computer-readable storage medium, the computer-readable storage medium storing a computer program, the computer program being executed by a processor to implement the steps in each of the method embodiments described above.
The embodiments of the invention also provide an information data processing terminal. The information data processing terminal is used to implement the steps in the embodiments of the methods mentioned above, and the information data processing terminal is not limited to mobile phone, computer and switch.
Embodiments of the invention also provide a server which, when implemented on an electronic device, provides a user input interface to implement the steps in embodiments of the methods described above.
Embodiments of the present invention provide a computer program product that, when the computer program product is run on an electronic device, causes the steps in each of the method embodiments described above to be implemented when the electronic device is executed.
The integrated units can be stored on a computer-readable storage medium if they are implemented as software functional units and sold or used as stand-alone products. Based on such an understanding, the present invention implements all or part of the processes in the method of the above embodiments, which may be accomplished by means of a computer program to instruct the relevant hardware, said computer program may be stored in a computer readable storage medium which, when executed by a processor, implements the steps of each of the above method embodiments. Where the computer program includes computer program code, the computer program code may be in the form of source code, in the form of object code, in the form of an executable file or in some intermediate form, etc. The computer readable medium may include at least: any entity or device capable of carrying computer program code to a photographic device/terminal device, a recording medium, computer memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media.
The embodiment of the present invention was compared with the average attenuation coefficient at 3-8 m in water obtained by AFD method at the delay step size of 1 ns and 2 ns and by Beer-Lambert Law by the above experiments, and the relative errors were 0.98% and 3.91%, respectively. The relative errors of the attenuation coefficients in water obtained by the AFD method at a delay step size of 1 ns are slightly smaller, which is in line with the analytical results of the present invention.
In addition, when the delay steps are 1 ns and 2 ns, respectively, the attenuation coefficients of the water bodies obtained from the fitting method are 0.1287 and 0.1232, respectively. their relative errors are 3.23% and 7.37%, respectively. The relative errors are larger compared to the AFD method. This is due to the fact that the fitting method needs to measure a sufficient number of points to ensure its accuracy, the AFD method can overcome this drawback. When the delay step is 2 ns, the relative error can be controlled to 3.91%, even though the number of measured points is only half of the number when the step is 1 ns.
The relative error between the attenuation coefficients measured at different propagating distances by the AFD method with delay step time of 1 ns and Beer-Lambert Law was further compared. the relative errors for the Beer-Lambert Law are 6.82%, 0.82%, 5.26%, 3.78%, 9.38% and the relative errors for the AFD method are 4.76%, 2.89%, 2.59%, 4.82% and 0.11% in the ranges of 3-4 m, 3-5 m, 3-6 m, 3-7 m and 3-8 m. Their average relative errors are 5.21% and 3.04%. The standard deviation of the attenuation coefficient of water bodies in different distance ranges measured by the AFD method is 0.0052, which is smaller than the one from the Beer-Lambert law as 0.0083. Compared with the one for the Beer-Lambert law, it can be clearly found that the absolute errors for the AFD method does not change abruptly as the range changes. Therefore, as the measurement range increases, the AFD method is more stable than the Beer-Lambert Law. These results fully demonstrated that the attenuation coefficient of water bodies measured by AFD method in real time is feasible and credible.
The experimental results show that the present invention deduces the relationship between the backscattering intensity and the attenuation coefficient of the water body by establishing the backscattering model of the water body. And a method for remote real-time measurement of the attenuation coefficient of water bodies is proposed in combination with the adjacent frame difference method. Based on this relationship and AFD method, the attenuation coefficient of the water body can be further calculated in real time. In order to verify the feasibility of this method, the self-developed flash lidar system was used to conduct a measurement study of the fresh water body in the distance range of 3-8 m. The attenuation coefficient of the water body is measured as 0.1343±0.02 and 0.1382±0.03 under 1 ns and 2 ns delay step time of the AFD method. Theoretical analysis shows that the measurement results are more accurate for a delay step size of 1 ns. Compared these values to the one achieved following the conventional Beer-Lambert Law (0.1330±0.02), which consists well with each other. The AFD method with 1 ns delay step size has an error of about 1% compared to Beer-Lambert Law. By comparing the Beer-Lambert law and the AFD method in different distance ranges, the experimental results and theoretical analysis fully verify the feasibility and reliability of the AFD method to measure the attenuation coefficient of water bodies. A new scheme is proposed for the real-time remote measurement of attenuation coefficients in water bodies.
The above mentioned is only a specific and superior implementation of the present invention, but the scope of protection of the present invention is not limited to it. Any modification, equivalent replacement and improvement, etc. made by any person skilled in the art within the technical scope disclosed by the present invention and within the spirit and principles of the present invention shall be covered by the scope of protection of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
202211204456.0 | Sep 2022 | CN | national |