A method for extracting aviation release characteristics in a uranium resource sample set

Abstract

The present application relates to a kind of extraction methods of uranium resource sample set in flight release feature, specifically to a kind of extraction methods of uranium resource sample set in flight release feature.The present application includes the following steps: step 1, known uranium deposit is converted into polygonal uranium mine area information;Step 2, sample label data in uranium mine area is generated according to fixed distance;Step 3, generate the flight release data of grid format;Step 4, the sample label corresponding to the feature value of flight release data generated by polygonal uranium mine area is calculated;Step 5, the calculation result is saved as uranium resource flight release sample data set.The present application can more accurately reflect the situation of flight release feature in uranium mineralization range.

Description

TECHNICAL FIELD

[0001] The application relates to a method for extracting aeromagnetic features in a uranium resource sample set. BACKGROUND

[0002] In recent years, with the explosive growth of geological data and the rapid development of computer science and technology, new methods and technologies represented by machine learning are gradually applied to the field of metallogenic prediction. At present, some scholars have used machine learning methods to carry out metallogenic modeling analysis. The key to the application of uranium resource exploration big data is to combine massive data in the field of uranium resource exploration with machine learning and data mining technologies. Among them, the establishment of a large number of uranium resource sample sets is the basis for the application of uranium resource big data in the later period.

[0003] Machine learning needs a large number of sample data with labels to train the model, and the data used should contain two parts of information, that is, "features" and "labels". Among them, the features are used to represent things and map a series of information of the labels, and the labels are the results of prediction or classification. In order to realize the application of uranium resource machine learning such as data mining method and machine learning prediction, the data related to mineralization such as geology, geophysical exploration, geochemical exploration and remote sensing need to be converted into sample features that can be fused with uranium mineralization label information. The aeromagnetic information is an important prospecting feature for the uranium resource sample set, and it is closely related to the uranium mineralization with a shallow burial depth.

[0004] The conventional method for extracting aeromagnetic features generally converts the aeromagnetic data format into a raster form, calculates the pixel value of the aeromagnetic data at the position of the known uranium mineralization through spatial analysis method, and directly uses the calculated pixel value as the aeromagnetic information feature corresponding to the uranium mineralization label. The advantage of this method is that the calculation is simple and convenient, and the aeromagnetic features of different levels of uranium mineralization in the working area can be quickly obtained. However, the conventional method also has obvious disadvantages. Because the uranium mineralization data is often the projection of the mine or ore body on the plane, the pixel value of the aeromagnetic data at the position of the projection point of the uranium mineralization may not be able to obtain the true aeromagnetic information corresponding to the uranium mineralization. How to relatively truly reflect and calculate the aeromagnetic uranium, thorium, potassium and uranium-thorium ratio information corresponding to the uranium mineralization is the problem to be solved by the application. SUMMARY

[0005] The purpose of the application is to provide a method for extracting aeromagnetic features in a uranium resource sample set, which enlarges the range of the projected uranium mineral point, replaces the projected point form of the uranium mineralization with the polygon uranium mine area range, generates a plurality of sample label points with equal spacing within the uranium mine area range, calculates the aeromagnetic pixel value at the position of the sample label points as the feature value corresponding to the sample label, and can more accurately reflect the situation of the aeromagnetic features in the uranium mineralization range.

[0006] The technical scheme adopted by the present application is as follows:

[0007] A method for extracting aerial and ground features in a uranium resource sample set, comprising the following steps: step 1, converting known uranium ore spots into polygonal uranium mine area information; step 2, generating sample label data within the uranium mine area according to a fixed distance; step 3, generating aerial and ground data in a grid format; step 4, calculating feature values of the aerial and ground data corresponding to the sample labels generated by the polygonal uranium mine area; and step 5, saving the calculation results as a uranium resource aerial and ground sample data set.

[0008] The step 1 comprises the following steps:

[0009] Step 1.1, reading in uranium ore spot data using ArcGIS software, with the data name being U.shp;

[0010] Step 1.2, performing buffer analysis processing on the uranium ore spot data U.shp using the buffer analysis function in the spatial analysis module of the ArcGIS software.

[0011] In the step 1.2, a uranium mineralization square region is generated according to a length of 900 meters as the side length and the uranium mineralization point as the center.

[0012] The step 2 comprises the following steps:

[0013] Step 2.1, obtaining sample label points bq; step 2.2, saving the sample label points bq obtained in step 2.1 as a point file, denoted as bq.shp; step 2.3, opening the attribute table of bq.shp and adding a mineralization level attribute field; and adding a belonging mineralization name attribute field to the attribute table of bq.shp.

[0014] In the step 2.1, an interval distance of 300 meters is used to generate equidistant points within the square region obtained in step 1.2 using the equidistant point generation function in the spatial analysis of arcgis, i.e., each uranium mineralization point is expanded into 9 label points, and these equidistant label points are denoted as bq, which are the sample label points.

[0015] The step 3 comprises the following steps:

[0016] Step 3.1, reading in the aerial and ground uranium contour line file, denoted as U.shp; reading in the aerial and ground thorium contour line file, denoted as Th.shp; and reading in the aerial and ground K contour line file, denoted as K.shp;

[0017] Step 3.2, using the contour data conversion to raster data function in arcgis, respectively, U.shp, Th.shp, K.shp is converted to raster format data, the converted flight U grid data is RasterU.tif, the converted flight thorium grid data is recorded as RasterTh.if, the converted flight potassium grid data is recorded as RasterK.tif;

[0018] Step 3.3, calculate the uranium thorium ratio data.

[0019] In step 3.3, the ratio of the flight U grid data RasterU.tif and the flight thorium grid data RasterTh.tif is calculated using the raster calculator in the GIS software, and the calculation result is the uranium thorium ratio data in the grid format, recorded as RasterUth.tif.

[0020] The step 4 includes the following steps:

[0021] Step 4.1, read the sample label point file bq.shp obtained in step 2.3, read the flight U grid data RasterU.tif, the flight thorium grid data RasterTh.tif, the flight potassium grid data RasterK.tif obtained in step 3.2, and the uranium thorium ratio data RasterUth.tif obtained in step 3.3;

[0022] Step 4.2, using the raster value extraction to point function in arcgis spatial analysis, the pixel value of the flight U grid data RasterU corresponding to each point position in the label point file bq.shp is calculated, an attribute field is generated in bq.shp, recorded as the attribute field name recorded as flight U, the calculated pixel value is written into the attribute field, and the bq.shp file after calculation is recorded as bq_u.shp

[0023] Step 4.3, using the raster value extraction to point function in GIS spatial analysis, the pixel value of the flight thorium grid data RasterTh corresponding to each point position in the label point file bq_u.shp is calculated, an attribute field is generated in bq_u.shp, recorded as the attribute field name recorded as flight Th, the calculated pixel value is written into the attribute field, and the bq_u.shp file after calculation is recorded as bq_u_th.shp;

[0024] Step 4.4, the pixel value of the raster data RasterK corresponding to the position of each point in the label point file bq_u_th.shp is calculated by using the grid value extraction to point function in the GIS spatial analysis, an attribute field is generated in bq_u_th.shp, and the attribute field name is recorded as the flight K, the calculated pixel value is written into the attribute field, and the bq_u_th.shp file after the calculation is completed is recorded as bq_u_th_k.shp;

[0025] Step 4.5, the pixel value of the raster data RasterUth corresponding to the position of each point in the label point file bq_u_th_k.shp is calculated by using the grid value extraction to point function in the GIS spatial analysis, an attribute field is generated in bq_u_th_k.shp, and the attribute field name is recorded as the flight Uth, the calculated pixel value is written into the attribute field, and the bq_u_th_k.shp file after the calculation is completed is recorded as bq_u_th_k_uth.shp.

[0026] The attribute table of the file bq_u_th_k_uth.shp is saved as a table document format, and is recorded as bq_u_th_k_uth.xlsx.

[0027] The table contains six fields, namely the mineralization name, the mineralization level, the flight U, the flight Th, the flight K and the flight Uth.

[0028] Compared with the prior art, the beneficial effects of the present application are that:

[0029] (1) The method for extracting the flight characteristics in the uranium resource sample set provided by the present application can expand the sample label from the uranium mineralization projection point to the polygon range of the uranium mineralization area, that is, the polygon uranium mineral area range is used to replace the uranium mineralization in the form of the projection point, so that the flight characteristics in the known uranium mineralization range can be more accurately reflected, and the real flight information characteristics in the uranium mineralization range are met.

[0030] (2) The method for extracting the flight characteristics in the uranium resource sample set provided by the present application can effectively increase the sample quantity of the sample set, and provide sample data for more accurate machine learning modeling. DETAILED DESCRIPTION

[0031] Figure 1 It is a distribution diagram of known uranium points;

[0032] Figure 2 It is a schematic diagram of the buffer analysis result of the uranium points;

[0033] Figure 3 It is a schematic diagram of the expansion of the uranium points into label points;

[0034] Figure 4 Overlaying the tag points with the airborne uranium data. DETAILED DESCRIPTION

[0035] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative effort fall within the protection scope of the present application.

[0036] In the description of the present application, it should be noted that the orientations or positional relationships indicated by the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", and the like are based on the orientations or positional relationships shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and therefore cannot be understood as indicating or implying that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation on the present application. In addition, the terms "first", "second", "third" are only for the purpose of description, and cannot be understood as indicating or implying relative importance.

[0037] In the description of the present application, it should be noted that unless otherwise explicitly specified and limited, the terms "mounting", "connection", "connection" should be understood in a broad sense, for example, it can be fixed connection, or detachable connection, or integrally connected; it can be mechanical connection, or electrical connection; it can be directly connected, or indirectly connected through an intermediate medium, or the internal communication of two elements. For a person of ordinary skill in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0038] The uranium resource sample set feature extraction method provided by the present application includes the following steps: step 1, converting the known uranium ore points into polygonal uranium mine area information. This step includes the following sub-steps:

[0039] Step 1.1, read in the uranium ore point data using ArcGIS software, the data name is U.shp, the file contains 4 uranium mineralization, and the 4 uranium mineralization names are A, B, C, and D, respectively, as shown in Figure 1 The mineralization levels of points A, B, C, and D are ore deposit, ore point, mineralization point, and anomaly point, respectively.

[0040] Step 1.2, using the buffer analysis function in the spatial analysis module of the ArcGIS software to perform buffer analysis processing on the uranium ore point data U.shp, generating a uranium mineralization square area with a length of 900 meters as the side length and the uranium mineralization point as the center. That is, the 4 uranium mineralization points in the U.shp file are converted into 4 square regions with a side length of 900 meters, as shown inFigure 2 .

[0041] Step 2, generate sample label data in the uranium mine area according to fixed distance.

[0042] Step 2.1, use the "generate equidistant points" function in arcgis spatial analysis to generate equidistant points in the 4 square areas obtained in step 1.2 at an interval distance of 300 meters, i.e. each uranium mineralization point is expanded to 9 label points, and a total of 4 uranium mineralization points are expanded to 36 label points. These equidistant label points are denoted as bq, which are sample label points, such as Figure 3 .

[0043] Step 2.2, save the sample label points bq obtained in step 2.1 as a point file, denoted as bq.shp.

[0044] Step 2.3 open the attribute table of bq.shp, add a "mineralization level" attribute field, and add a label value "mineralization level" to the sample label points bq according to the mineralization levels of uranium mineralization A, B, C and D, i.e. the mineralization level of uranium mineralization point A is "ore deposit", and the mineralization level of the 9 label points generated from A is also "ore deposit", and the mineralization levels of the label points expanded from B, C and D are also assigned according to the corresponding mineralization levels. Add an "associated mineralization name" attribute field to the attribute table of bq.shp, i.e. the name of uranium mineralization point A is "A", and the "associated mineralization name" attribute of the 9 label points generated from A is assigned "A", and the names of uranium mineralization B, C and D are added to the corresponding expanded sample label points. The attribute table after completing step 2.3 is shown in Table 1.

[0045] Table 1 attribute table of expanded sample label points

[0046]

[0047]

[0048] Step 3, generate raster format aerial data. This step includes the following substeps:

[0049] Step 3.1, read in the aerial uranium contour file, denoted as U.shp; read in the aerial thorium contour file, denoted as Th.shp; read in the aerial K contour file, denoted as K.shp.

[0050] Step 3.2, using the contour data conversion to raster data function in arcgis, respectively, U.shp, Th.shp, K.shp is converted to raster format data, the converted flight and put uranium grid data is RasterU.tif, the converted flight and put thorium grid data is recorded as RasterTh.if, the converted flight and put potassium grid data is recorded as RasterK.tif.

[0051] Step 3.3, calculate the ratio of uranium thorium data. Using the raster calculator in GIS software, the ratio of flight and put uranium grid data RasterU.tif and flight and put thorium grid data RasterTh.tif is calculated, and the calculation result is the uranium thorium ratio data in raster format, recorded as RasterUth.tif.

[0052] Step 4, calculate the sample label corresponding to the flight and put data feature value generated by the polygon uranium mine area. This step includes the following sub steps:

[0053] Step 4.1, read in the sample label point file bq.shp obtained in step 2.3, read in the flight and put uranium grid data RasterU.tif, flight and put thorium grid data RasterTh.tif, flight and put potassium grid data RasterK.tif, and uranium thorium ratio data RasterUth.tif obtained in step 3.3. The schematic diagram of the sample label point file bq.shp and the flight and put uranium grid data RasterU.tif is as follows Figure 4 .

[0054] Step 4.2, using the "raster value extraction to point" function in arcgis spatial analysis, calculate the pixel value of each point position in the label point file bq.shp corresponding to the flight and put uranium grid data RasterU, generate an attribute field in bq.shp, recorded as the attribute field name is recorded as "flight and put U", write the calculated pixel value into the attribute field, and record the bq.shp file after calculation as bq_u.shp.

[0055] Step 4.3, using the "raster value extraction to point" function in GIS spatial analysis, calculate the pixel value of each point position in the label point file bq_u.shp corresponding to the flight and put thorium grid data RasterTh, generate an attribute field in bq_u.shp, recorded as the attribute field name is recorded as "flight and put Th", write the calculated pixel value into the attribute field, and record the bq_u.shp file after calculation as bq_u_th.shp.

[0056] Step 4.4, using the "raster value to point" function in the GIS spatial analysis, the pixel value of the each point position in the label point file bq_u_th.shp corresponding to the aerial potassium grid data RasterK is calculated, an attribute field is generated in bq_u_th.shp, and the attribute field name is recorded as "aerial K", the calculated pixel value is written into the attribute field, and the bq_u_th.shp file after the calculation is recorded as bq_u_th_k.shp.

[0057] Step 4.5, using the "raster value to point" function in the GIS spatial analysis, the pixel value of the each point position in the label point file bq_u_th_k.shp corresponding to the aerial potassium grid data RasterUth is calculated, an attribute field is generated in bq_u_th_k.shp, and the attribute field name is recorded as "aerial Uth", the calculated pixel value is written into the attribute field, and the bq_u_th_k.shp file after the calculation is recorded as bq_u_th_k_uth.shp.

[0058] Step 5, the calculation result is saved as the uranium resource aerial sample data set. That is, the attribute table of the file bq_u_th_k_uth.shp is saved as a table document format, recorded as bq_u_th_k_uth.xlsx, and the table contains six fields respectively as the mineralization name, the mineralization level, the aerial U, the aerial Th, the aerial K and the aerial Uth, and the data set is as shown in Table 2.

[0059] Table 2 uranium resource aerial sample data table

[0060]

[0061]

[0062] It is apparent to those skilled in the art that the present application is not limited to the details of the foregoing exemplary embodiments, but can be implemented in other specific forms without departing from the spirit or essential characteristics of the present application. Therefore, the embodiments should be considered in all respects as illustrative and not restrictive, the scope of the present application being defined by the appended claims rather than the foregoing description, and it is intended that all changes and modifications which come within the meaning and range of equivalency of the claims are resolutely intended to be embraced therein. Any reference signs in the claims should not be construed as limiting the claims to the figures in which the reference signs are used.

[0063] Furthermore, it should be understood that although the specification is described in terms of embodiments, not every embodiment includes every feature described. The specification can include implicit combinations of explicitly mentioned features and / or explicit combinations of implicitely mentioned features. Each embodiment depends on the explicit combinations of features and / or the implicit combinations of features made specifically within that embodiment, and each such embodiment can be combined with every other such embodiment to create further embodiments.

Claims

1. A method for extracting airborne characteristics in a set of uranium resource samples, characterized in that: It comprises the following steps: step 1, converting known uranium ore spots into polygonal uranium mine information; step 2, generating sample label data in the uranium mine according to a fixed distance; step 3, generating flight and release data in a raster format; step 4, calculating characteristic values of flight and release data corresponding to the sample label generated by the polygonal uranium mine; and step 5, saving the calculation results as a flight and release sample data set of uranium resources. The step 2 comprises the following steps: Step 2.1, obtaining sample label points bq; step 2.2, saving the sample label points bq obtained in step 2.1 as a point file, denoted as bq.shp; and step 2.3, opening the attribute table of bq.shp and adding a mineralization level attribute field; and adding a belonging mineralization name attribute field to the attribute table of bq.shp. In the step 2.1, an equal-interval point function in spatial analysis of arcgis is used to generate equal-interval points in the square region obtained in step 1.2 at an interval distance of 300 meters, i.e., each uranium mineralization point is expanded into 9 label points, and these equal-interval label points are denoted as bq, which are the sample label points.

2. The method of claim 1, wherein the method comprises: The step 1 comprises the following steps: Step 1.1, reading in uranium ore point data using ArcGIS software, and the data is named as U.shp; Step 1.2, performing buffer analysis processing on the uranium ore point data U.shp using a buffer analysis function in a spatial analysis module of the ArcGIS software.

3. The method of claim 2, wherein the method further comprises: determining the number of the uranium resource sample sets; and determining the number of the uranium resource sample sets based on the number of the uranium resource sample sets. In the step 1.2, a square region of uranium mineralization is generated with a length of 900 meters as the side length and the uranium mineralization point as the center.

4. The method of claim 1, wherein the method further comprises: The step 3 comprises the following steps: Step 3.1, reading in flight and release uranium contour line files, denoted as U.shp; reading in flight and release thorium contour line files, denoted as Th.shp; and reading in flight and release K contour line files, denoted as K.shp; Step 3.2, converting U.shp, Th.shp and K.shp into raster format data using an equal contour line data conversion function in arcgis, and the converted flight and release uranium raster data is denoted as RasterU.GIF, the converted flight and release thorium raster data is denoted as RasterTh.GIF, and the converted flight and release potassium raster data is denoted as RasterK.GIF; Step 3.3, calculating uranium thorium ratio data.

5. The method of claim 4, wherein the method further comprises: In the step 3.3, a raster calculator in the GIS software is used to calculate the ratio of the flight and release uranium raster data RasterU.GIF and the flight and release thorium raster data RasterTh.GIF, and the calculation result is uranium thorium ratio data in a raster format, denoted as RasterUth.GIF.

6. The method of claim 5, wherein the method further comprises: The step 4 comprises the following steps: Step 4.1, reading in the sample label point file bq.shp obtained in step 2.3, reading in the flight and release uranium raster data RasterU.GIF, the flight and release thorium raster data RasterTh.GIF, the flight and release potassium raster data RasterK.GIF and the uranium thorium ratio data RasterUth.GIF obtained in step 3.3; Step 4.2, using the raster value extraction to point function in the arcgis spatial analysis, calculate the pixel value of the RasterU corresponding to the position of each point in the label point file bq.shp, generate an attribute field in bq.shp, and record the attribute field name as the flight U. Write the calculated pixel value to the attribute field. The bq.shp file after calculation is recorded as bq_u.shp; Step 4.3, using the raster value extraction to point function in the arcgis spatial analysis, calculate the pixel value of the RasterTh corresponding to the position of each point in the label point file bq_u.shp, generate an attribute field in bq_u.shp, and record the attribute field name as the flight Th. Write the calculated pixel value to the attribute field. The bq_u.shp file after calculation is recorded as bq_u_th.shp; Step 4.4, using the raster value extraction to point function in the arcgis spatial analysis, calculate the pixel value of the RasterK corresponding to the position of each point in the label point file bq_u_th.shp, generate an attribute field in bq_u_th.shp, and record the attribute field name as the flight K. Write the calculated pixel value to the attribute field. The bq_u_th.shp file after calculation is recorded as bq_u_th_k.shp; Step 4.5, using the raster value extraction to point function in the arcgis spatial analysis, calculate the pixel value of the RasterUth corresponding to the position of each point in the label point file bq_u_th_k.shp, generate an attribute field in bq_u_th_k.shp, and record the attribute field name as the flight Uth. Write the calculated pixel value to the attribute field. The bq_u_th_k.shp file after calculation is recorded as bq_u_th_k_uth.shp.

7. The method of claim 6, wherein the method further comprises: determining the number of the uranium resource sample sets; and determining the number of the uranium resource sample sets based on the number of the uranium resource sample sets. The step 5 is to save the attribute table of the file bq_u_th_k_uth.shp as a table document format, recorded as bq_u_th_k_uth.xlsx.

8. The method of claim 7, wherein the method further comprises: determining the number of the uranium resource sample sets; and determining the number of the uranium resource sample sets based on the number of the uranium resource sample sets. The table contains 6 fields respectively as the name of the mineralization, the level of mineralization, the flight U, the flight Th, the flight K and the flight Uth.

Images