1. A method for acquiring sea pollution space distribution is characterized by comprising the following steps:

calculating the path distances between a monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, wherein the data of the sea area to be researched comprises the coordinates of the monitoring station, the flow velocity of seawater and the flow field data of seawater;

determining an interpolation reference station meeting set conditions according to the path distances of the position points;

calculating distance weights of points to be measured of the sea area to be researched according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station, wherein the points to be measured are position points in the plurality of position points of the sea area to be researched;

and calculating according to the distance weight to obtain the sea area pollution space distribution of the sea area to be researched.

2. The method according to claim 1, wherein the calculating the path distances between the monitoring station of the sea area to be researched and the plurality of position points of the sea area to be researched according to the data of the sea area to be researched comprises:

calculating a surface impedance factor of the pollution diffusion of the sea area to be researched according to the flow velocity of the seawater;

calculating a horizontal impedance factor of pollution diffusion of the sea area to be researched according to the flow direction of the seawater;

and calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor.

3. The method according to claim 2, wherein the calculation formula of the surface impedance factor of the pollution diffusion of the sea area to be researched according to the flow velocity of the seawater is as follows:

wherein, F (x)_i) Is a point x to be measured_iThe surface impedance factor at the point of the beam,is a point x to be measured_iFlow velocity of seawater at a point, S_minIs the minimum value of sea water flow velocity, S_maxIs the maximum value of the sea water flow velocity in the sea area.

4. The method according to claim 2, wherein the calculation formula of the horizontal impedance factor of the pollution spread of the sea area to be researched according to the seawater flow direction data is as follows:

wherein, HF (x)_i) Is a point x to be measured_iThe horizontal impedance factor at the point of the point,is the included angle between the diffusion direction of the pollutants and the flow direction of the seawater in the seawater flow field data.

5. The method according to claim 2, wherein the calculation formula for calculating the path distances between the monitoring station of the sea area to be researched and the plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor is as follows:

d (A, B) is the path distance from the monitoring station A to the monitoring station B, N is the number of grid pixels passing between the monitoring station A and the monitoring station B, and D (x)_i,x_i+1) Is an adjacent point x to be measured_iAnd x_i+1Grid cell path distance between, F (x)_i) Is a point x to be measured_iSurface impedance factor of the spot, HF (x)_i) Is a point x to be measured_iHorizontal impedance factor of (2), F (x)_i) Is a point x to be measured_iPoint adjacent site x_i+1Surface impedance factor of, HF (x)_i) Is a point x to be measured_iAdjacent site x_i+1Horizontal impedance factor of (x)₁A) indicates that the starting point of the path distance from monitored station a to monitored station B is monitored station a.

6. The method according to claim 1, wherein the determining an interpolation reference station satisfying a set condition according to the path distances of the plurality of position points comprises:

judging whether the current reference monitoring station meets a set condition or not based on a convex hull algorithm, wherein the current reference monitoring station is a specified number of monitoring stations closest to the station to be detected when the current reference monitoring station is judged for the first time;

if the preset conditions are met, taking the current reference monitoring station as an interpolation reference station;

and if the set condition is not met, updating the current reference monitoring station to obtain an updated current reference position point.

7. The method according to claim 1, wherein the calculating the sea area pollution spatial distribution of the sea area to be researched according to the distance weight comprises:

calculating according to the distance weight to obtain the concentration of the pollutants at the point to be measured;

and obtaining the spatial distribution of the concentration of the pollutants in the sea area according to the concentration of the pollutants at the point to be detected.

8. A sea area polluted space interpolation device is characterized by comprising:

a first calculation module: the system comprises a data acquisition unit, a data processing unit and a data processing unit, wherein the data acquisition unit is used for acquiring data of a sea area to be researched, calculating path distances between a monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, and the data of the sea area to be researched comprises coordinates of the monitoring station, seawater flow velocity and seawater flow field data;

a judging module: the interpolation reference station is used for determining an interpolation reference station meeting a set condition according to the path distances of the position points;

a second calculation module: the distance weight of a point to be measured of the sea area to be researched is calculated according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station, wherein the point to be measured is a position point in the plurality of position points of the sea area to be researched;

a third calculation module: and the sea area pollution space distribution of the sea area to be researched is obtained through calculation according to the distance weight.

9. An electronic device, comprising: a processor, a memory storing machine-readable instructions executable by the processor, the machine-readable instructions when executed by the processor performing the steps of the method of any of claims 1 to 7 when the electronic device is run.

10. A computer-readable storage medium, having stored thereon a computer program which, when being executed by a processor, is adapted to carry out the steps of the method according to any one of claims 1 to 7.

Background

For the spatial interpolation of marine pollution, the traditional spatial interpolation method adopts Euclidean distance in weight calculation, and a reference point selection method is a fixed range or fixed point number, so that two problems are easily ignored when the reference point selection method is used for the marine pollution interpolation: (1) the sea has a soft natural boundary, the diffusion of pollutants in the sea is influenced by factors such as ocean currents, the diffusion path is not a straight line, and the diffusion path is not only influenced by the blocking of lands such as peninsula and islands, but also influenced by the speed and direction of the ocean currents to a great extent. (2) The diffusion distribution of the ocean pollution sources generally conforms to a certain trend, namely if the monitoring result of a certain station near the shore is heavy pollution, and the monitoring results of the station near the shore and far sea are better in pollution condition, the water quality of the farther sea area is certain to be lighter in pollution degree. The influence of surrounding reference points on interpolation points is exaggerated by a traditional space interpolation Euclidean distance weighting method, the reference points are fixedly selected, interpolation is seriously distorted in an area with large water quality span, a pollution diffusion rule is violated, and the final interpolation result cannot reflect the real distribution condition and the diffusion trend of sea area pollution.

Disclosure of Invention

In view of the above, an object of the embodiments of the present application is to provide a method and an apparatus for obtaining a spatial distribution of marine pollution, an electronic device, and a readable storage medium. The method can effectively express the acquisition method of the sea area pollution space distribution under the action of the horizontal motion of ocean currents.

In order to achieve the above purpose, the technical solutions provided in the embodiments of the present application are as follows:

in a first aspect, an embodiment of the present application provides a method for obtaining a spatial distribution of sea pollution, including: calculating the path distances between a monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, wherein the data of the sea area to be researched comprises the coordinates of the monitoring station, the flow velocity of seawater and the flow field data of seawater; determining an interpolation reference station meeting set conditions according to the path distances of the position points; calculating distance weights of points to be measured of the sea area to be researched according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station, wherein the points to be measured are position points in the plurality of position points of the sea area to be researched; and calculating according to the distance weight to obtain the sea area pollution space distribution of the sea area to be researched.

According to the embodiment of the application, the interpolation reference point is selected according to the path distance between the monitoring station and the point to be measured, and the distribution of the sea area pollution space with the sea area to be researched is calculated according to the distance weight of the path distance of the plurality of position points of the interpolation reference point meeting the standard, so that the distribution trend of the sea area pollution can be directly and accurately reflected.

With reference to the first aspect, an embodiment of the present application provides a first possible implementation manner of the first aspect, where: calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, wherein the path distances comprise the following steps: calculating a surface impedance factor of the pollution diffusion of the sea area to be researched according to the flow velocity of the seawater; calculating a horizontal impedance factor of pollution diffusion of the sea area to be researched according to the seawater flow direction data; and calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor.

According to the method and the device, the surface impedance factor of pollution diffusion is calculated through the seawater flow velocity data, the horizontal impedance factor of pollution diffusion is calculated through the seawater flow velocity data, the path distance can be calculated according to the obtained surface impedance factor and the horizontal impedance factor, the path distances between the monitoring station and the multiple position points of the sea area to be researched can be reflected according to the actual seawater flow velocity and the flow direction, the obtained path distance data is consistent with the real-time seawater data, and the accuracy of the data is guaranteed.

With reference to the first possible implementation manner of the first aspect, an embodiment of the present application provides a second possible implementation manner of the first aspect, where: and calculating a calculation formula of the surface impedance factor of the pollution diffusion of the sea area to be researched according to the flow velocity of the seawater, wherein the calculation formula comprises the following steps:

wherein, F (x)_i) Is a point x to be measured_iThe surface impedance factor at the point of the beam,is a point x to be measured_iFlow velocity of seawater at a point, S_minIs the minimum value of sea water flow velocity, S_maxIs the maximum value of the sea water flow velocity in the sea area.

With reference to the second possible implementation manner of the first aspect, the present application provides a third possible implementation manner of the first aspect, where a calculation formula for calculating a horizontal impedance factor of pollution diffusion of the sea area to be studied according to the seawater flow direction data is as follows:

wherein, HF (x)_i) Is a point x to be measured_iThe horizontal impedance factor at the point of the point,is the included angle between the diffusion direction of the pollutants and the flow direction of the seawater in the seawater flow field data.

With reference to the third possible implementation manner of the first aspect, an embodiment of the present application provides a fourth possible implementation manner of the first aspect, where a calculation formula of path distances between a monitoring station of the sea area to be researched and a plurality of location points of the sea area to be researched, which is obtained by calculation according to the surface impedance factor and the horizontal impedance factor, is as follows:

wherein D (A, B) is the path distance from the monitoring station A to the monitoring station B, N is the number of grid pixels passing between the monitoring station A and the monitoring station B, and D (x)_i,x_i+1) Is an adjacent point x to be measured_iAnd x_i+1Grid in betweenUnit path distance, F (x)_i) Is a point x to be measured_iSurface impedance factor of the spot, HF (x)_i) Is a point x to be measured_iHorizontal impedance factor of (2), F (x)_i) Is a point x to be measured_iPoint adjacent site x_i+1Surface impedance factor of, HF (x)_i) Is a point x to be measured_iAdjacent site x_i+1Horizontal impedance factor of (x)₁A) indicates that the starting point of the path distance from monitored station a to monitored station B is monitored station a.

With reference to the fourth possible implementation manner of the first aspect, an embodiment of the present application provides a fifth possible implementation manner of the first aspect, where determining an interpolation reference station that meets a set condition according to path distances of the multiple position points includes: judging whether the current reference monitoring station meets a set condition or not based on a convex hull algorithm, wherein the current reference monitoring station is a specified number of monitoring stations closest to the station to be detected when the current reference monitoring station is judged for the first time; if the preset conditions are met, taking the current reference monitoring station as an interpolation reference station; and if the set condition is not met, updating the current reference monitoring station to obtain an updated current reference position point.

Whether the monitoring station meets the conditions or not is judged through a convex hull algorithm, if yes, the monitoring station is determined to be an interpolation reference station, if not, the monitoring station is updated to be judged again, until the monitoring station meeting the conditions is determined to be used as the interpolation reference station, the interpolation reference station meeting the set conditions is selected in a circulating mode, the selection of the interpolation reference station is enabled to be more in line with the actual requirements, and the actual seawater pollution data are more fit.

With reference to the fifth possible implementation manner of the first aspect, an embodiment of the present application provides a sixth possible implementation manner of the first aspect, where the calculating, according to the distance weight, to obtain a sea area pollution spatial distribution of the sea area to be researched includes: calculating according to the distance weight to obtain the concentration of the pollutants at the point to be measured; and obtaining the spatial distribution of the concentration of the pollutants in the sea area according to the concentration of the pollutants at the point to be detected.

According to the method and the device, the pollutant concentration is calculated according to the distance weight obtained by calculating the actual seawater data, and then interpolation calculation is carried out on the concentration, so that the obtained sea area pollutant concentration space distribution situation can effectively reflect the distribution trend of sea area pollution under the action of ocean currents.

In a second aspect, an embodiment of the present application further provides a sea area polluted space interpolation apparatus, including: a first calculation module: the system comprises a data acquisition unit, a data processing unit and a data processing unit, wherein the data acquisition unit is used for acquiring data of a sea area to be researched, calculating path distances between a monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, and the data of the sea area to be researched comprises coordinates of the monitoring station, seawater flow velocity and seawater flow direction data; a judging module: the interpolation reference station is used for determining an interpolation reference station meeting a set condition according to the path distances of the position points; a second calculation module: the distance weight of a point to be measured of the sea area to be researched is calculated according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station, wherein the point to be measured is a position point in the plurality of position points of the sea area to be researched; a third calculation module: and the sea area pollution space distribution of the sea area to be researched is obtained through calculation according to the distance weight.

In a third aspect, an embodiment of the present application further provides an electronic device, including: a processor, a memory storing machine-readable instructions executable by the processor, the machine-readable instructions, when executed by the processor, performing the steps of the method of the first aspect described above, or any possible implementation of the first aspect, when the electronic device is run.

In a fourth aspect, this embodiment of the present application further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of the method in the first aspect or any one of the possible implementation manners of the first aspect.

In order to make the aforementioned objects, features and advantages of the present application more comprehensible, embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.

Fig. 1 is a block diagram of an electronic device according to an embodiment of the present disclosure.

Fig. 2 is a flowchart of a method for obtaining sea area pollution space distribution according to an embodiment of the present application.

Fig. 3 is a detailed flowchart of step 201 of a method for obtaining sea area pollution space distribution according to an embodiment of the present application.

Fig. 4 is a flowchart illustrating a step 202 of a method for obtaining a spatial distribution of marine pollution according to an embodiment of the present application.

Fig. 5 is a functional module schematic diagram of a sea area pollution space distribution obtaining apparatus according to an embodiment of the present application.

Detailed Description

The technical solution in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.

It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only for distinguishing the description, and are not to be construed as indicating or implying relative importance.

The method comprises a diffusion simulation model based on a sewage drain and a spatial interpolation based on monitoring station data, wherein the former needs a large amount of previous work to obtain parameters such as hydrology, meteorology and runoff to establish a pollution diffusion numerical simulation model, the calculation is complex, the construction period is long, sea area environment parameters are usually obtained by monitoring stations, and the distribution of the sea area monitoring stations is often sparse and uneven, so that the method is not suitable for the research on the spatial distribution of the pollution concentration in a large-scale sea area. A space interpolation method based on Geographic statistics and GIS (Geographic Information System) technology completely depends on environment monitoring data, fully excavates useful Information of the environment monitoring data, and avoids complex calculation. The inverse distance weight interpolation method has the advantages of simple principle, convenience in calculation, conformity with the first law of geography and the like, and is widely applied to environmental pollution simulation.

Based on the research, the sea area pollution space distribution obtaining method, the sea area pollution space distribution obtaining device, the electronic equipment and the readable storage medium improve the traditional space interpolation method in weight calculation and reference point selection, and can effectively simulate the space distribution and diffusion trend of pollutants under the action of ocean currents.

Example one

For the convenience of understanding the present embodiment, first, an electronic device for performing a method for acquiring a sea pollution spatial distribution disclosed in the embodiments of the present application will be described in detail.

As shown in fig. 1, is a block schematic diagram of an electronic device. The electronic device 100 may include a memory 111, a memory controller 112, a processor 113, a peripheral interface 114, an input-output unit 115, and a display unit 116. It will be understood by those of ordinary skill in the art that the structure shown in fig. 1 is merely exemplary and is not intended to limit the structure of the electronic device 100. For example, electronic device 100 may also include more or fewer components than shown in FIG. 1, or have a different configuration than shown in FIG. 1.

The above-mentioned elements of the memory 111, the memory controller 112, the processor 113, the peripheral interface 114, the input/output unit 115 and the display unit 116 are electrically connected to each other directly or indirectly, so as to implement data transmission or interaction. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The processor 113 is used to execute the executable modules stored in the memory.

The Memory 111 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like. The memory 111 is configured to store a program, and the processor 113 executes the program after receiving an execution instruction, and the method executed by the electronic device 100 defined by the process disclosed in any embodiment of the present application may be applied to the processor 113, or implemented by the processor 113.

The processor 113 may be an integrated circuit chip having signal processing capability. The Processor 113 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the Integrated Circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, or a discrete hardware component. The various methods, steps, and logic blocks disclosed in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.

The peripheral interface 114 couples various input/output devices to the processor 113 and memory 111. In some embodiments, the peripheral interface 114, the processor 113, and the memory controller 112 may be implemented in a single chip. In other examples, they may be implemented separately from the individual chips.

The input/output unit 115 is used to provide input data to the user. The input/output unit 115 may be, but is not limited to, a mouse, a keyboard, and the like.

The display unit 116 provides an interactive interface (e.g., a user operation interface) between the electronic device 100 and the user or is used for displaying image data to the user for reference. In this embodiment, the display unit may be a liquid crystal display or a touch display. In the case of a touch display, the display can be a capacitive touch screen or a resistive touch screen, which supports single-point and multi-point touch operations. The support of single-point and multi-point touch operations means that the touch display can sense touch operations simultaneously generated from one or more positions on the touch display, and the sensed touch operations are sent to the processor for calculation and processing.

The electronic device 100 in this embodiment may be configured to perform each step in each method provided in this embodiment. The following describes in detail an implementation procedure of the sea pollution spatial distribution acquisition method by several embodiments.

Example two

Please refer to fig. 2, which is a flowchart illustrating a method for obtaining a spatial distribution of marine pollution according to an embodiment of the present application. The specific process shown in fig. 2 will be described in detail below.

Step 201, calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched.

Optionally, the data of the sea area to be studied includes coordinates of monitoring stations, flow velocity of seawater and flow direction data of seawater.

Optionally, before step 201, the method may further include: and acquiring original data of the sea area to be researched and preprocessing the data to obtain the data of the sea area to be researched.

Optionally, the data of the sea area to be studied may further include: recording attribute data of various water quality monitoring index concentration information and a boundary range of a sea area to be researched.

Optionally, the various water quality monitoring index concentration information includes: information on the concentration of various contaminants.

Optionally, the determining of the boundary range of the sea area to be studied includes: starting from two end points of a land area shoreline, extending rays to a leading sea baseline in the offshore direction, and taking a sea area range defined by the shoreline, the two rays and the leading sea baseline as a boundary of a research sea area.

Optionally, the preprocessing the raw data of the sea area to be researched includes: and selecting station data in the offshore sea area of the research area according to the original data of the sea area to be researched, and rejecting other station data.

Optionally, the preprocessing the raw data of the sea area to be researched further includes: and screening the data of each station in the original data of the sea area to be researched to remove the data of the station without the record of the pollutant concentration information.

Optionally, the preprocessing the raw data of the sea area to be researched further includes: and converting the flow velocity and flow direction information recorded by the sampling points into raster data based on the sampling point data in the path distance of the plurality of position points.

Optionally, the preprocessing the raw data of the sea area to be researched further includes: and carrying out projection and cutting processing on the obtained boundary range of the sea area to be researched and the grid to form projection coordinates.

In one embodiment, as shown in fig. 3, step 201 may include the following steps 2011 to 2013.

And 2011, calculating a surface impedance factor of the pollution diffusion of the sea area to be researched according to the flow velocity of the seawater.

In one embodiment, the calculation formula for calculating the surface impedance factor of the pollution spread of the sea area to be studied from the flow velocity of the sea water is:

wherein, F (x)_i) Is a point x to be measured_iThe surface impedance factor at the point of the beam,is a point x to be measured_iFlow velocity of seawater at a point, S_minIs the minimum value of sea water flow velocity, S_maxIs the maximum value of the sea water flow velocity in the sea area.

Alternatively, the flow velocity of the seawater at the point x of the point to be measured can be obtained from the flow field data.

Optionally, the sea area sea water flow velocity minimum value and the sea area sea water flow velocity maximum value may be obtained by sorting flow velocities in the flow field data.

Alternatively, the surface impedance may be calculated using a calculation formula of the surface impedance factor using a grid calculation tool based on the grid data.

Step 2012, calculating the horizontal impedance factor of the pollution diffusion of the sea area to be studied according to the seawater flow field data.

Optionally, the seawater flow field data includes: the flow velocity of the seawater and the flow direction of the seawater.

In one embodiment, the calculation formula for calculating the horizontal impedance factor of the pollution spread of the sea area to be studied according to the seawater flow direction data is as follows:

wherein, HF (x)_i) Is a point x to be measured_iThe horizontal impedance factor at the point of the point,is the included angle between the diffusion direction of the pollutants and the flow direction of the seawater in the seawater flow field data.

Optionally, the direction of the spread of the contaminant is the direction of the current seawater flow.

Optionally, the included angle between the seawater flow direction at the current position and the seawater flow directions at a plurality of positions around the current position is the included angle between the pollutant diffusion direction and the water flow direction.

Alternatively, the surface impedance may be calculated using a calculation formula of the surface impedance factor based on the water flow direction data.

Step 2013, calculating and obtaining the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor.

In one embodiment, the calculation formula for calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor is as follows:

wherein D (A, B) is the path distance from the monitoring station A to the monitoring station B, N is the number of grid pixels passing between the monitoring station A and the monitoring station B, and D (x)_i,x_i+1) Is an adjacent point x to be measured_iAnd x_i+1Grid cell path distance between, F (x)_i) Is a point x to be measured_iSurface impedance factor of the spot, HF (x)_i) Is a point x to be measured_iHorizontal impedance factor of (2), F (x)_i) Is a point x to be measured_iPoint adjacent site x_i+1Surface impedance factor of, HF (x)_i) Is a point x to be measured_iAdjacent site x_i+1Horizontal impedance factor of (x)₁A) indicates that the starting point of the path distance from monitored station a to monitored station B is monitored station a.

Optionally, the path distance may be calculated from the monitoring point to the monitoring point by using a path distance calculation module according to the surface impedance factor grid, the horizontal impedance factor grid, and the water quality monitoring point data.

Step 202, determining an interpolation reference station meeting the set conditions according to the path distances of the plurality of position points.

In one embodiment, as shown in FIG. 4, step 202 may include the following steps.

And judging whether the current reference monitoring station meets the set conditions or not based on a convex hull algorithm.

And if the preset conditions are met, taking the current reference monitoring station as an interpolation reference station.

And if the set condition is not met, updating the current reference monitoring station to obtain an updated current reference position point.

Optionally, when the first judgment is made, the current reference monitoring station is a specified number of monitoring stations closest to the station to be detected.

Alternatively, the reference site needs to satisfy the following condition: (1) the polygon enclosed by the reference sites is a convex polygon. (2) The point to be measured is located in the convex polygon. For example, the convex polygon may be a convex quadrangle, a convex pentagon, a convex hexagon, etc., and the convex polygon is not a convex quadrangle as an example for description.

Optionally, before the convex hull algorithm is performed, the longitude and latitude position information of the monitored station can be converted into an array.

Optionally, whether the monitored site meets the condition is determined based on the following function method:

a＝(m2-m1)*(n-n1)-(n2-n1)*(m-m1)

b＝(m3-m2)*(n-n2)-(n3-n2)*(m-m2)

c＝(m4-m3)*(n-n3)-(n4-n3)*(m-m3)

d＝(m1-m4)*(n-n4)-(n1-n4)*(m-m4)

wherein m1 and n1 are longitude and latitude coordinates of a first point of the convex quadrilateral respectively; m2 and n2 are longitude and latitude coordinates of a second point of the convex quadrilateral respectively; m3 and n3 are longitude and latitude coordinates of a third point of the convex quadrilateral respectively; m4 and n4 are longitude and latitude coordinates of a fourth point of the convex quadrilateral respectively; and m and n are respectively longitude and latitude coordinates of the point to be measured.

Optionally, if the monitored site satisfies the condition, the function returns 1, otherwise returns 0.

Optionally, the above function method is used for each point to be measured, and whether the monitored stations around the point to be measured meet the condition for interpolation is judged.

Optionally, the number of the current reference monitoring stations may be 4, 5, or 6, and the specific number may be adjusted according to an actual situation, which is not limited in this embodiment of the application.

Illustratively, if the number of the current reference monitoring stations is 4, selecting 4 monitoring stations closest to the point to be measured as the previous reference monitoring stations, if the judgment result does not meet the set condition, selecting a 5 th monitoring point close to the point to be measured to replace the point farthest from the point to be measured in the original 4 monitoring stations, and judging again until the monitoring station meeting the condition is determined to be the interpolation reference station.

Illustratively, if the number of the current reference monitoring stations is 5, selecting 5 monitoring stations closest to the point to be measured as the previous reference monitoring stations, if the judgment does not meet the set condition, selecting the 6 th monitoring point close to the point to be measured to replace the point farthest from the point to be measured in the original 5 monitoring stations, and judging again until the monitoring station meeting the condition is determined to be the interpolation reference station.

Illustratively, if the number of the current reference monitoring stations is 6, 6 monitoring stations closest to the point to be measured are selected as the previous reference monitoring stations, if the judgment result does not meet the set condition, the 7 th monitoring point close to the point to be measured is selected to replace the point farthest from the point to be measured in the original 6 monitoring stations, and the judgment is carried out again until the monitoring station meeting the condition is determined to be the interpolation reference station.

And 203, calculating the distance weight of the point to be measured in the sea area to be researched according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station.

Optionally, the point to be measured is a location point of a plurality of location points of the sea area to be studied.

And 204, calculating according to the distance weight to obtain the sea area pollution space distribution of the sea area to be researched.

Optionally, step 204 comprises: and calculating the pollutant concentration of the point to be measured according to the distance weight from the interpolation reference point to the point to be measured and the pollutant concentration of the interpolation reference point.

Optionally, step 204 further comprises: and carrying out spatial interpolation based on the pollutant concentration of the interpolation reference point to obtain the pollutant concentration of the point to be measured.

Optionally, the interpolation formula of the measured point pollutant concentration is as follows:

wherein S is the concentration of the pollutants at the point to be measured; lambda [ alpha ]_iAn inverse distance weight for each interpolated reference point; p_iThe contaminant concentration at the site is monitored and N is the number of interpolation points.

Optionally, the calculation formula of the inverse distance weight of each interpolation reference point is:

wherein λ is_iInverse distance weights, D, for respective interpolated reference points_iThe distance of a path from a monitoring station to a point to be measured is obtained; p is an index of distance; n is the number of interpolation points.

Optionally, the pollutant concentration of the measuring point is mapped in a coordinate system to obtain the spatial distribution of the pollutant concentration of the sea area.

In the embodiment, when the monitoring station data is used for obtaining the distribution condition of the sea area pollution space, the space interpolation method based on the traditional inverse distance weight is improved. The vast majority of marine pollutants originate from land, land-source pollutants are gathered at a river mouth and are diffused in the sea, the sea is a flowing space and has a soft natural boundary, the land-source pollutants are not only simply accumulated along the river mouth line or a coastline but also are unevenly diffused under the action of ocean currents in the space distribution of the sea, the diffusion path is not straight, and the diffusion of the marine pollution has a specific tendency.

Based on the principle, the invention uses an improved sea area pollution spatial interpolation method to calculate the pollutant diffusion path distance under the action of ocean currents, and the pollutant diffusion path distance is used as the weight calculation basis of spatial interpolation, so that the interpolation result can effectively reflect the distribution trend of the ocean current pollution in the sea area.

EXAMPLE III

Based on the same application concept, the embodiment of the present application further provides a sea area pollution space distribution device corresponding to the sea area pollution space distribution method, and as the principle of solving the problem of the device in the embodiment of the present application is similar to that of the embodiment of the sea area pollution space distribution method, the implementation of the device in the embodiment of the present application can refer to the description in the embodiment of the above method, and repeated parts are not described again.

Please refer to fig. 5, which is a schematic diagram of functional modules of a marine pollution space distribution apparatus according to an embodiment of the present application. Each module in the sea area polluted space distribution device in the embodiment is used for executing each step in the above method embodiment. A sea area pollution space distribution device comprises a first calculation module 301, a judgment module 302, a second calculation module 303 and a third calculation module 304; wherein the content of the first and second substances,

the first calculation module 301: the method is used for calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the data of the sea area to be researched, wherein the data of the sea area to be researched comprises the coordinates of the monitoring station, the flow speed of the seawater and the flow direction data of the seawater.

The judging module 302: and the interpolation reference station is used for determining the interpolation reference station meeting the set conditions according to the path distances of the plurality of position points.

The second calculation module 303: and the distance weight of the point to be measured of the sea area to be researched is calculated according to the interpolation reference station and the path distances of the plurality of position points corresponding to the interpolation reference station, wherein the point to be measured is a position point in the plurality of position points of the sea area to be researched.

The third calculation module 304: and the method is used for calculating and obtaining the sea area pollution space distribution of the sea area to be researched according to the distance weight.

In a possible implementation, the first computing module 301 is further configured to: calculating a surface impedance factor of pollution diffusion of a sea area to be researched according to the flow velocity of the seawater; calculating a horizontal impedance factor of pollution diffusion of a sea area to be researched according to the flow direction data of the seawater; (ii) a And calculating the path distances between the monitoring station of the sea area to be researched and a plurality of position points of the sea area to be researched according to the surface impedance factor and the horizontal impedance factor.

In a possible implementation manner, the determining module 302 is specifically configured to: whether the current reference monitoring station meets the set conditions or not is judged based on a convex hull algorithm, and when the current reference monitoring station is judged for the first time, the current reference monitoring station is a specified number of monitoring stations closest to the point to be measured. And if the preset conditions are met, taking the current reference monitoring station as an interpolation reference station. And if the set condition is not met, updating the current reference monitoring station to obtain an updated current reference position point.

In a possible implementation, the third calculating module 304 is further configured to: calculating according to the distance weight to obtain the concentration of the pollutants at the point to be measured; and obtaining the spatial distribution of the concentration of the pollutants in the sea area according to the concentration of the pollutants at the point to be detected.

In addition, an embodiment of the present application further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program performs the steps of the method for acquiring the distribution of marine pollution spaces described in the above method embodiments.

The computer program product of the sea area pollution space distribution obtaining method provided in the embodiment of the present application includes a computer readable storage medium storing a program code, where instructions included in the program code may be used to execute steps of the sea area pollution space distribution obtaining method described in the above method embodiment, which may be specifically referred to in the above method embodiment, and details are not described here.

In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

In addition, functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes. It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.

The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.

The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.