Layered automatic drawing method and system for pollution range and remediation range of soil pollutants

Abstract

The invention discloses a layered automatic drawing method and system for a pollution range and a restoration range of a soil pollutant, and the drawing method is based on a written Arcpy script, and mainly comprises the steps: obtaining the basic information of the soil pollutant; reading layered soil pollution information and a land parcel red line range in the basic information; performing inverse distance interpolation on the layered soil pollution information to obtain pollutant concentration distribution; according to the pollutant concentration distribution and the land parcel red line range, pollutant concentration distribution in the land parcel red line is obtained; and according to the pollutant remediation target value and the screening value, the concentration distribution of the pollutants in the land parcel red line is subjected to reclassification, and the pollution range and remediation range of the layered soil pollutants are obtained. According to the method, batch processing and analysis tasks of data can be automatically and rapidly completed, the working efficiency is greatly improved, great flexibility and convenience are provided for users, and the time and error rate of manual operation are remarkably reduced.

Description

Technical Field

[0001] The present invention belongs to the technical field of soil pollutant range mapping, and in particular relates to a method and system for automatically mapping the pollution range and remediation range of soil pollutants in layers. Background Art

[0002] Survey units usually use ArcGIS software to depict the scope of soil pollution and remediation. The application of ArcGIS not only improves the efficiency of data processing, but also enhances the depth and breadth of analysis, providing a scientific basis for soil pollution prevention and control.

[0003] Although ArcGIS provides a powerful graphical user interface (GUI) for users to perform interactive operations, when processing a large amount of pollutant data and performing the depiction of complex pollution ranges and remediation ranges, manually setting parameters, adjusting layer coordinate systems, drawing precise boundaries and other operations may be very time-consuming. In addition, human factors can easily lead to error accumulation, affecting the accuracy of the final results, and it is also not easy for beginners to get started. Therefore, it is urgent to solve the problem of low efficiency and poor accuracy in depicting the pollution range and remediation range of a large number of stratified soil pollutants. Summary of the invention

[0004] In order to solve the above technical problems, the present invention proposes a method and system for automatically drawing the pollution range and remediation range of soil pollutants in layers, which can solve the problems of low efficiency and poor accuracy when depicting the pollution range and remediation range of a large number of layers of soil pollutants.

[0005] The present invention provides a method for automatically drawing the pollution range and remediation range of soil pollutants in layers, based on the written Arcpy script, which mainly includes:

[0006] Based on the written Arcpy script, it mainly includes:

[0007] Obtaining basic information on soil pollutants, wherein the basic information on soil pollutants includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is collected, pollutant type, pollutant screening value, pollutant remediation target value, and concentration of each pollutant in the sample;

[0008] Reading the layered soil pollution information and the red line range of the land parcel in the basic information;

[0009] Performing inverse distance interpolation on the layered soil pollution information to obtain pollutant concentration distribution;

[0010] According to the pollutant concentration distribution and the red line range of the plot, obtaining the pollutant concentration distribution within the red line range of the plot;

[0011] According to the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line range of the plot is reclassified to obtain the stratified soil pollutant pollution range and remediation range.

[0012] Optionally, before reading the layered soil pollution information and the land red line range in the basic information, the method includes: preparing a land red line Shapefile vector file and placing it in an ArcGIS work path folder.

[0013] Optionally, obtaining the pollutant concentration distribution includes: performing inverse distance interpolation using an Idw function;

[0014] Using the Idw function to perform inverse distance difference includes:

[0015] Obtain pollutant concentration values ​​at known points;

[0016] A weighted estimate is made based on the distance between the pollutant concentration value at the known point and the unknown point to obtain the pollutant concentration distribution.

[0017] Optionally, according to the pollutant concentration distribution and the red line range of the plot, obtaining the pollutant concentration distribution within the red line range of the plot includes:

[0018] According to the plot redline Shapefile vector file, obtain the plot redline range;

[0019] Based on the red line range of the plot, the pollutant concentration distribution is clipped to obtain the pollutant concentration distribution within the red line range of the plot.

[0020] Optionally, based on the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line of the plot is reclassified to obtain the stratified soil pollutant contamination range and remediation range, including:

[0021] Classify the pollutant concentration distribution within the red line of the plot according to the pollutant screening value and the pollutant remediation target value;

[0022] The categories with less than the pollutant screening value are not reflected, and the pollution range of the stratified soil pollutants is obtained;

[0023] Categories with values ​​smaller than the pollutant remediation target are not reflected, and the stratified soil pollutant remediation range is obtained.

[0024] The present invention also provides a layered automatic drawing system for soil pollutant contamination range and remediation range, based on the written Arcpy script, which mainly includes: an information acquisition module, an information reading module, an interpolation module, a clipping module and a reclassification module;

[0025] The information acquisition module is used to obtain basic information on soil pollutants, which includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is located, pollutant type, pollutant screening value, pollutant remediation target value and concentration of each pollutant in the sample;

[0026] The information reading module is used to read the layered soil pollution information and the red line range of the plot in the basic information;

[0027] The interpolation module is used to perform inverse distance interpolation on the layered soil pollution information;

[0028] The trimming module is used to trim the pollutant concentration distribution;

[0029] The reclassification module is used to reclassify the pollutant concentration distribution within the red line range of the land plot.

[0030] Optionally, before reading the layered soil pollution information and the land red line range in the basic information, the method includes: preparing a land red line Shapefile vector file and placing it in an ArcGIS work path folder.

[0031] Optionally, performing inverse distance interpolation on the layered soil pollution information includes: performing inverse distance interpolation using an Idw function.

[0032] Optionally, the cropping module includes: a mask unit and a cropping processing unit;

[0033] The mask unit is used to obtain the red line range of the plot according to the red line Shapefile vector file of the plot;

[0034] The clipping processing unit is used to clip the pollutant concentration distribution within the red line range of the plot.

[0035] Optionally, the reclassification module includes: a classification unit and a category processing unit;

[0036] The classification unit is used to classify the pollutant concentration distribution within the red line range of the plot according to the pollutant screening value and the pollutant remediation target value;

[0037] The category processing unit is used to perform non-reflection processing on categories that are smaller than the pollutant screening value and the pollutant remediation target value.

[0038] Compared with the prior art, the present invention has the following advantages and technical effects:

[0039] 1. Through the Arcpy script of the present invention, the complex operation process of layered soil pollution scope and remediation scope can be automatically processed, which is very useful for scenarios that need to process a large amount of sampling point data and generate continuous surface models. It can quickly complete batch processing and analysis tasks of data, greatly improve work efficiency, provide users with great flexibility and convenience, and significantly reduce the time and error rate of manual operations.

[0040] 2. Interactive operation: The Arcpy script of the present invention can be interactively operated in the Python window of ArcGIS, so that the user can quickly input the established pollution information database and can also flexibly adjust the execution parameters therein.

[0041] 3. The Arcpy script of the present invention can support more complex depth applications, such as joint clipping based on multiple geometric figures, clipping based on specific attributes, etc. These complex clipping conditions may be difficult to implement directly in the ArcGIS graphical interface, but can be easily completed through the Arcpy script. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] The drawings constituting a part of the present application are used to provide a further understanding of the present application. The illustrative embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0043] Figure 1 It is a flow chart of a method for automatically drawing the contaminated area and remediation area of ​​soil pollutants in layers according to an embodiment of the present invention;

[0044] Figure 2 Schematic diagram of the pollutant pollution range and remediation range outputted by an embodiment of the present invention, wherein (a) is a schematic diagram of the first-layer selenium pollution range and remediation range Shapefile, (b) is a schematic diagram of the first-layer fluoride pollution range and remediation range Shapefile, (c) is a schematic diagram of the first-layer phenol pollution range and remediation range Shapefile, (d) is a schematic diagram of the second-layer selenium pollution range and remediation range Shapefile, (e) is a schematic diagram of the second-layer fluoride pollution range and remediation range Shapefile, and (f) is a schematic diagram of the second-layer phenol pollution range and remediation range Shapefile. DETAILED DESCRIPTION

[0045] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0046] It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and that, although a logical order is shown in the flowcharts, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0047] The present invention proposes a method for automatically drawing the pollution range and remediation range of soil pollutants in layers, based on the written Arcpy script, such as Figure 1 As shown, the specific steps include:

[0048] Obtain basic information on soil pollutants;

[0049] Read the layered soil pollution information and land red line range in the basic information;

[0050] Perform inverse distance interpolation on layered soil pollution information to obtain pollutant concentration distribution;

[0051] According to the pollutant concentration distribution and the red line range of the plot, obtain the pollutant concentration distribution within the red line range of the plot;

[0052] According to the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line of the plot is reclassified to obtain the stratified soil pollutant pollution range and remediation range.

[0053] Furthermore, the basic information of soil pollutants includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is collected, pollutant type, pollutant screening value, pollutant remediation target value and concentration of each pollutant in the sample;

[0054] Furthermore, before reading the layered soil pollution information and the land red line range in the basic information, the following steps are included: preparing a land red line Shapefile vector file and placing it in the ArcGIS work path folder.

[0055] Specifically, according to the survey results of the soil pollution status of the plot and the remediation requirements, the basic information of soil pollutants obtained will be filled in MuBan.xlsx (Excel custom file name), as shown in Table 1; in addition, a plot red line Shapefile (vector file) will be prepared and placed in the ArcGIS work path folder.

[0056] Table 1

[0057]

[0058] Furthermore, Arcpy reading was performed: as shown in Table 2, combined with the actual project requirements, the three variables f_workspace, Poldata_f, and red_boundary in the script were modified and clarified, and then all the script contents were copied and pasted into the Python window of ArcGIS software and the Enter key was pressed to run. ArcGIS software automatically identified the data frame coordinate system, pollution layer, sampling point coordinates, sample pollutant concentration and other information through computing libraries such as pandas, numpy, and arcpy, and the Excel database was converted into a layered soil pollution information Shapefile.

[0059] Table 2

[0060]

[0061]

[0062]

[0063] Furthermore, performing inverse distance interpolation on the layered soil pollution information includes: performing inverse distance interpolation using an Idw function.

[0064] In this step, as an embodiment that can be added, the inverse distance weighted (IDW) interpolation uses a linear weight combination of a set of sampling points to determine the pixel value, the weight is an inverse distance function, and the surface to be interpolated should be a surface with local dependent variables. The IDW method assumes that the mapped variable is reduced due to the influence of the distance between its sampling position, that is, the concentration value of the unknown point can be estimated by the weighted average of the concentration values ​​of the known points around it, and the weight depends on the distance between the unknown point and the known point.

[0065] Specifically, inverse distance interpolation: IDW is an interpolation method based on the similarity of sample points. It estimates the value of a cell by calculating the weighted average of the sample points in the adjacent area, and then interpolates to obtain a surface. According to the layered soil pollution information Shapefile output in step 2, ArcGIS software automatically performs inverse distance interpolation through the Idw function to obtain the concentration distribution of each pollutant in each layer of soil in the polluted plot, and outputs the inverse distance interpolation Tiff file.

[0066] Further, according to the pollutant concentration distribution and the layered soil pollution information, obtaining the pollutant concentration distribution within the mask range includes:

[0067] Get the red line range of the plot according to the red line Shapefile vector file;

[0068] Based on the red line range of the plot, the pollutant concentration distribution is clipped to obtain the pollutant concentration distribution within the red line range of the plot.

[0069] Specifically, raster file clipping: Based on the obtained red_boundary variable information, ArcGIS software automatically identifies the plot redline file in the working path, and automatically clips the output inverse distance interpolation Tiff file through the ExtractByMask function, and outputs the clipped Tiff file of the concentration distribution of each pollutant in each layer of soil within the red line range of the plot.

[0070] Furthermore, according to the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line of the plot is reclassified to obtain the stratified soil pollutant pollution range and remediation range including:

[0071] Classify the distribution of pollutant concentrations within the red line of the plot according to the pollutant screening value and the pollutant remediation target value;

[0072] The categories with less than the pollutant screening value are not reflected, and the pollution range of stratified soil pollutants is obtained;

[0073] Categories with values ​​smaller than the pollutant remediation target are not reflected, and the stratified soil pollutant remediation range is obtained.

[0074] Specifically, raster file reclassification: According to the pollutant screening value and pollutant remediation target value obtained in step one, ArcGIS software automatically reclassifies and changes the pixel values ​​of the cropped Tiff file through the RemapValue and Reclassify functions. Pixel values ​​less than the screening value or the remediation target value are classified into one category, and pixel values ​​greater than the screening value or the remediation target value are classified into one category, and the reclassified Tiff file is output.

[0075] Convert raster files to Shapefiles: ArcGIS software automatically converts the reclassified Tiff file output in step 5 into a Shapefile through the arcpy.RasterToPolygon_conversion function, and uses the arcpy.da.UpdateCursor and cursor.deleteRow functions to eliminate the categories that are smaller than the screening value or the restoration target value, and outputs the stratified soil pollutant pollution range and restoration range Shapefile, which is convenient for users to make subsequent maps, perform statistical analysis, and other work.

[0076] The present invention also provides a layered automatic drawing system for soil pollutant contamination range and remediation range, based on a written Arcpy script, including: an information acquisition module, an information reading module, an interpolation module, a clipping module and a reclassification module;

[0077] Information acquisition module, used to obtain basic information on soil pollutants;

[0078] An information reading module is used to read the layered soil pollution information and the red line range of the land parcel in the basic information;

[0079] Interpolation module, used to perform inverse distance interpolation of layered soil pollution information;

[0080] A clipping module is used to clip the pollutant concentration distribution;

[0081] The reclassification module is used to reclassify the pollutant concentration distribution within the red line of the land.

[0082] Furthermore, the basic information of soil pollutants includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is located, pollutant type, pollutant screening value, pollutant remediation target value and concentration of each pollutant in the sample; in addition, prepare the plot red line Shapefile (vector file) and place it in the ArcGIS work path folder.

[0083] Furthermore, performing inverse distance interpolation on the layered soil pollution information includes: performing inverse distance interpolation using an Idw function.

[0084] Further, the cropping module includes: a mask unit and a cropping processing unit;

[0085] The mask unit is used to obtain the red line range of the plot according to the red line Shapefile (vector file);

[0086] The clipping processing unit is used to clip the pollutant concentration distribution within the red line of the land plot.

[0087] Further, the reclassification module includes: a classification unit and a category processing unit;

[0088] The classification unit is used to classify the pollutant concentration distribution within the red line of the plot according to the pollutant screening value and the pollutant remediation target value;

[0089] The category processing unit is used to perform non-reflection processing on categories that are smaller than the pollutant screening value and the pollutant remediation target value.

[0090] This example is described in detail for a certain construction land:

[0091] The strata in a certain construction land plot are plain fill (0-1.5m), silt (1.5-3.5m) and silt clay (water-proof layer) from top to bottom. There are 16 survey points in the plot, numbered S01-S16. One sample is taken from each point in the plain fill layer, silt layer and silt clay layer for inspection. The test results show that there are points where selenium, fluoride and phenol exceed the second-class screening value in the plain fill layer and silt layer, and the silt clay layer shows that the above three pollutants do not exceed the second-class screening value. After risk assessment, the risk control values ​​of the three pollutants are greater than the second-class screening value, so the risk control value is selected as the remediation target value of each pollutant. In order to compare the pollution range and remediation range of each pollutant, the pollution range and remediation range of selenium, fluoride and phenol are automatically characterized for the plain fill layer and silt layer of the plot.

[0092] Step 1: Obtain basic information:

[0093] According to the survey results of the plot, the plain fill layer was set as the first layer, the silt layer was set as the second layer, the coordinates of 16 survey points were obtained, the sample number, sample depth, pollutant type, pollutant screening value, pollutant remediation target value and the concentration of each pollutant in the sample were collected and filled in MuBan.xlsx, as shown in Tables 3 and 4; in addition, the plot red line Shapefile "plot red line.shp" was prepared and placed in the ArcGIS working path "D:\desktopfile\Pollution and remediation scope characterization".

[0094] Table 3

[0095]

[0096] Table 4

[0097]

[0098] Step 2. Execute Arcpy reading: According to the actual project requirements, modify the three variables f_workspace, Poldata_f, and red_boundary in the script to "f_workspace = r"D:\desktopfile\Pollution and Remediation Scope Characterization", Poldata_f = "MuBan.xlsx", red_boundary = r"Plot Red Line.shp"", copy all the script contents and paste them into the ArcGIS software Python window and press Enter to run. ArcGIS software automatically identifies the data frame coordinate system, pollution layer, sampling point coordinates, sample pollutant concentration and other information through pandas, numpy, arcpy and other computing libraries, and converts the Excel database into a layered soil pollution information Shapefile.

[0099] Step 3: Inverse distance interpolation: Based on the layered soil pollution information Shapefile output in step 2, ArcGIS software automatically performs inverse distance interpolation through the Idw function to obtain the concentration distribution of each pollutant in each layer of soil in the polluted plot, and outputs the inverse distance interpolation Tiff file.

[0100] Step 4: Raster file clipping: Based on the red_boundary variable information obtained in step 2, ArcGIS software automatically identifies the land redline file in the working path, and automatically clips the inverse distance interpolation Tiff file output in step 3 through the ExtractByMask function, and outputs the clipped Tiff file of the concentration distribution of each pollutant in each layer of soil within the red line range of the land.

[0101] Step 5. Reclassify raster files: Based on the pollutant screening value and pollutant remediation target value obtained in step 1, ArcGIS software automatically reclassifies and changes the pixel values ​​of the clipped Tiff file output in step 4 through the RemapValue and Reclassify functions. Pixel values ​​less than the screening value or the remediation target value are classified into one category, and pixel values ​​greater than the screening value or the remediation target value are classified into one category, and the reclassified Tiff file is output.

[0102] Step 6: Convert raster files to Shapefiles: ArcGIS software automatically converts the reclassified Tiff files output in step 5 into Shapefiles through the arcpy.RasterToPolygon_conversion function, and uses arcpy.da.UpdateCursor and cursor.deleteRow functions to exclude the categories that are less than the screening value or the restoration target value. Finally, the Shapefiles of selenium pollution range and restoration range of the plain fill layer and silt layer, the Shapefiles of fluoride pollution range and restoration range, and the Shapefiles of phenol pollution range and restoration range are output. For details, see Figure 2 (a)-(f).

[0103] The above are only preferred specific implementations of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed in the present application should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for automatically drawing the pollution range and remediation range of soil pollutants in layers, characterized in that: Based on the written Arcpy script, it mainly includes: Obtaining basic information on soil pollutants, wherein the basic information on soil pollutants includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is collected, pollutant type, pollutant screening value, pollutant remediation target value, and concentration of each pollutant in the sample; Reading the layered soil pollution information and the red line range of the land parcel in the basic information; Performing inverse distance interpolation on the layered soil pollution information to obtain pollutant concentration distribution; According to the pollutant concentration distribution and the red line range of the plot, obtaining the pollutant concentration distribution within the red line range of the plot; According to the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line range of the plot is reclassified to obtain the stratified soil pollutant pollution range and remediation range.

2. The method for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 1 is characterized in that: Before reading the layered soil pollution information and the land red line range in the basic information, the following steps are included: preparing a land red line Shapefile vector file and placing it in an ArcGIS work path folder.

3. The method for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 1 is characterized in that: Obtaining the pollutant concentration distribution includes: using the Idw function to perform inverse distance interpolation; Using the Idw function to perform inverse distance difference includes: Obtain pollutant concentration values ​​at known points; A weighted estimate is made based on the distance between the pollutant concentration value at the known point and the unknown point to obtain the pollutant concentration distribution.

4. The method for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 1 is characterized in that: According to the pollutant concentration distribution and the red line range of the plot, obtaining the pollutant concentration distribution within the red line range of the plot includes: According to the plot redline Shapefile vector file, obtain the plot redline range; Based on the red line range of the plot, the pollutant concentration distribution is clipped to obtain the pollutant concentration distribution within the red line range of the plot.

5. The method for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 1 is characterized in that: According to the pollutant remediation target value and screening value, the pollutant concentration distribution within the red line of the plot is reclassified to obtain the stratified soil pollutant pollution range and remediation range including: Classify the distribution of pollutant concentrations within the red line of the plot according to the pollutant screening value and the pollutant remediation target value; The categories with less than the pollutant screening value are not reflected, and the pollution range of the stratified soil pollutants is obtained; Categories with values ​​smaller than the pollutant remediation target are not reflected, and the stratified soil pollutant remediation range is obtained.

6. A layered automatic drawing system for soil pollutant contamination range and remediation range, characterized by: Based on the written Arcpy script, it mainly includes: information acquisition module, information reading module, interpolation module, clipping module and reclassification module; The information acquisition module is used to obtain basic information on soil pollutants, which includes: survey point number, point coordinates, sample number, sample depth, soil layer number where the sample is located, pollutant type, pollutant screening value, pollutant remediation target value and concentration of each pollutant in the sample; The information reading module is used to read the layered soil pollution information and the red line range of the plot in the basic information; The interpolation module is used to perform inverse distance interpolation on the layered soil pollution information; The cutting module is used to cut the pollutant concentration distribution; The reclassification module is used to reclassify the pollutant concentration distribution within the red line range of the land plot.

7. The system for automatically drawing the soil pollutant pollution range and remediation range according to claim 6 is characterized in that: Before reading the layered soil pollution information and the land red line range in the basic information, the following steps are included: preparing a land red line Shapefile vector file and placing it in an ArcGIS work path folder.

8. The system for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 6 is characterized in that: Performing inverse distance interpolation on the layered soil pollution information includes: performing inverse distance interpolation using an Idw function.

9. The system for automatically drawing the contaminated area and remediation area of ​​soil pollutants according to claim 6 is characterized in that: The cropping module includes: a mask unit and a cropping processing unit; The mask unit is used to obtain the red line range of the plot according to the red line Shapefile vector file of the plot; The clipping processing unit is used to clip the pollutant concentration distribution within the red line range of the plot.

10. The system for automatically drawing the soil pollutant pollution range and remediation range according to claim 7 is characterized in that: The reclassification module includes: a classification unit and a category processing unit; The classification unit is used to classify the pollutant concentration distribution within the red line range of the plot according to the pollutant screening value and the pollutant remediation target value; The category processing unit is used to perform non-reflection processing on categories that are smaller than the pollutant screening value and the pollutant remediation target value.

Images