GIS-based rare earth mineral rapid exploration method, electronic device and system

Abstract

The invention discloses a GIS-based rare earth mineral rapid exploration method, electronic equipment and system, and relates to the technical field of data processing. The GIS-based rare earth mineral rapid exploration method comprises the following steps: based on a target mining area, obtaining a division grid; obtaining a plurality of known points based on the divided grids; based on each known point, acquiring a concentration change network of each known point; based on the concentration change net, deviation coefficients of rare earth element concentration changes in all directions of the corresponding known point are obtained; obtaining the interpolation weight of each known point based on the corresponding deviation coefficient of each known point in each direction and the to-be-evaluated point; and obtaining the rare earth element concentration of the to-be-evaluated point based on the interpolation weight of each known point. Compared with the prior art that the interpolation weight is obtained only by adopting the physical distance with the to-be-evaluated point, the interpolation weight obtained through the deviation coefficient is more accurate, and the subsequently obtained rare earth element concentration of the to-be-evaluated point is more accurate.

Description

Technical Field

[0001] This application relates to the technical field of data processing, and specifically to a rapid exploration method, electronic device, and system for rare earth minerals based on GIS. Background Art

[0002] Rare earth minerals refer to ores containing rare earth elements. Geochemical exploration is one of the common exploration methods for rare earth minerals. By collecting surface soil and rock samples, analyzing the content of rare earth elements in them, and drawing a rare earth element distribution map, the location of the ore body can be identified, and the enrichment area of rare earth minerals can be found. Rare earth elements usually exist in the earth's crust at a relatively low concentration and are relatively dispersed, often showing small ore bodies or ore spots. Therefore, it is necessary to accurately identify potential mining areas to provide a scientific basis for resource development and utilization. In the prior art, the Geographic Information System (GIS) can establish a model based on geochemical data and combine known mineral data to evaluate the resource potential of mining areas. When collecting information during mineral resource exploration, sampling observations are usually used. In order to analyze the mineral distribution of the entire area, interpolation methods are usually required. Based on GIS, spatial distribution analysis of the geochemical exploration data model can be carried out. Through the inverse distance weighted interpolation method, according to the distribution of known geochemical data, the concentration distribution of rare earth elements in the unsampled area can be inferred, enriching the rare earth element concentration distribution data set to more accurately identify potential mining areas.

[0003] Estimating the concentration of rare earth elements in the unknown area through inverse distance weighted interpolation based on the concentration of rare earth elements in the known area by GIS, the basic idea is that if the distance between the known area and the unknown area is closer, the influence of the rare earth element concentration in the known area on the rare earth element concentration in the unknown area is greater. In the case of uniform change distribution of rare earth elements in each area, inverse distance weighted interpolation can obtain a relatively accurate estimation result. However, in fact, the distribution of rare earth elements in the land is not uniform. Therefore, when estimating using the inverse distance weighted interpolation method only based on the physical distance between the known points and the points to be evaluated, the obtained result is not accurate enough. Summary of the Invention

[0004] The purpose of this application is to provide a rapid exploration method, electronic device, and system for rare earth minerals based on GIS to solve the technical problem in the prior art that the concentration of rare earth elements at the points to be evaluated cannot be accurately estimated.

[0005] To achieve the above purpose, this application provides the following technical solutions: In the first aspect, this application proposes a rapid exploration method for rare earth minerals based on GIS. The rapid exploration method for rare earth minerals based on GIS includes: Based on the target mining area, obtain a divided grid; the target mining area is obtained in advance; Based on the divided grid, a plurality of known points are obtained; the known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known; Based on each known point, a concentration change network of each known point is obtained; Based on the concentration change network, the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point is obtained; the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration at the corresponding known point in each direction from the uniform change value; Based on the deviation coefficients corresponding to each direction of each known point and the point to be evaluated, the interpolation weight of each known point is obtained; the point to be evaluated is obtained in advance; Based on the interpolation weights of each known point, the rare earth element concentration of the point to be evaluated is obtained.

[0006] As a specific solution in the technical solution of the present application, the obtaining a concentration change network of each known point based on each known point includes: Based on each known point, a first known point is obtained; the first known point is any one of the known points; Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed centered on the first known point; Based on the first known point and each reference point, the rare earth element concentration change rate in each direction is obtained; each direction points from the first known point to the corresponding reference point; Based on the rare earth element concentration change rates in each direction, the concentration change network of the first known point is obtained.

[0007] As a specific solution in the technical solution of the present application, the obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point based on the concentration change network includes: Based on each known point, a second known point is obtained; the second known point is any one of the known points; Based on the second known point, a first reference point is obtained; the first reference point is any one of the reference points corresponding to the second known point; Based on the second known point and the first reference point, a plurality of verification points are obtained; the verification points are any known points among the known points located between the second known point and the first reference point; Based on each verification point, the deviation coefficient of the rare earth element concentration change in the first direction of the second known point is obtained; the first direction points from the second known point to the first reference point.

[0008] As a specific solution in the technical solution of this application, obtaining a plurality of verification points based on the second known point and the first reference point includes: Based on the second known point and the first reference point, obtain a first line segment; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; Based on the first line segment, obtain a first straight line and a second straight line; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; Take each known point located between the first straight line and the second straight line as each verification point.

[0009] As a specific solution in the technical solution of this application, obtaining the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on each verification point includes: Based on each verification point, obtain a first verification point; the first verification point is any one of the verification points; Based on the second known point and the first verification point, obtain the concentration change components of the second known point and the first verification point in the first direction; Based on the concentration change components, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction.

[0010] As a specific solution in the technical solution of this application, the calculation formula for obtaining the concentration change components of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ; Wherein, represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the rare earth element concentration of the i-th known point; represents the rare earth element concentration of the y-th verification point of the t-th reference point of the i-th known point; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the cosine function.

[0011] As a specific solution in the technical solution of this application, the calculation formula for obtaining the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component is as follows: ; Wherein, represents the deviation coefficient of the rare earth element concentration change of the i-th known point along the direction from the i-th known point to the t-th reference point; represents the variance of the concentration change components of all verification points between the i-th known point and the t-th reference point; represents the total number of all verification points between the i-th known point and the t-th reference point; represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point along the direction from the i-th known point to the t-th reference point; represents the concentration change rate in the direction from the i-th known point to the t-th reference point in the concentration change network of the i-th known point; represents taking the absolute value; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the sine function, represents the linear normalization function.

[0012] As a specific solution in the technical solution of this application, obtaining the interpolation weights of each known point based on the deviation coefficients corresponding to each direction of each known point and the point to be evaluated includes: Obtaining a third known point based on each known point; the third known point is any one of the known points; Obtaining an effective deviation coefficient based on the deviation coefficients in each direction of the third known point; the effective deviation coefficient is the deviation coefficient less than a preset value among the deviation coefficients; Obtaining a uniform adjustment coefficient of the third known point and the point to be evaluated based on each effective deviation coefficient; the uniform adjustment coefficient is at least used to characterize the degree of uniform change of the rare earth element concentration of the third known point along the direction from the third known point to the point to be evaluated; Obtaining the interpolation weight of the third known point based on the uniform adjustment coefficient.

[0013] In a second aspect, this application proposes a GIS-based electronic device for rapid exploration of rare earth minerals. The GIS-based electronic device for rapid exploration of rare earth minerals includes: A processing module, configured to obtain a divided grid based on a target mining area; the target mining area is obtained in advance. A reading module, configured to obtain a plurality of known points based on the divided grid; the known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known. The processing module is further configured to obtain a concentration change network of each known point based on each known point. And, based on the concentration change network, obtain the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration at the corresponding known point in each direction from the uniform change value. And, based on the deviation coefficients corresponding to each direction of the known points, obtain the interpolation weight of the known point. And, based on the interpolation weights of each known point, obtain the rare earth element concentration of each unknown point.

[0014] As a specific solution in the technical solution of the present application, the processing module is further configured to obtain a first known point based on each known point; the first known point is any one of the known points. Based on the first known point, obtain a plurality of reference points from each known point; each reference point is distributed centered on the first known point. And, based on the first known point and each reference point, obtain the rare earth element concentration change rate in each direction; each direction points from the first known point to the corresponding reference point. And, based on the rare earth element concentration change rates in each direction, obtain the concentration change network of the first known point.

[0015] As a specific solution in the technical solution of the present application, the processing module is further configured to obtain a second known point based on each known point; the second known point is any one of the known points. And, based on the second known point, obtain a first reference point; the first reference point is any one of the reference points corresponding to the second known point. And, based on the second known point and the first reference point, obtain a plurality of verification points; the verification points are any known points among the known points located between the second known point and the first reference point. And, based on each verification point, obtain the deviation coefficient of the rare earth element concentration change along the first direction of the second known point; the first direction points from the second known point to the first reference point.

[0016] As a specific solution in the technical solution of the present application, the processing module is further configured to obtain a first line segment based on the second known point and the first reference point; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; and, obtain a first straight line and a second straight line based on the first line segment; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; and, take each known point located between the first straight line and the second straight line as each verification point.

[0017] As a specific solution in the technical solution of the present application, the processing module is further configured to obtain a first verification point based on each verification point; the first verification point is any one of the verification points; and, obtain the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point; and, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component.

[0018] As a specific solution in the technical solution of the present application, the calculation formula for the processing module to obtain the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ; wherein, represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the rare earth element concentration of the i-th known point; represents the rare earth element concentration of the y-th verification point of the t-th reference point of the i-th known point; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the cosine function.

[0019] As a specific solution in the technical solution of the present application, the calculation formula for the processing module to obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component is as follows: ; wherein, The deviation coefficient representing the change in rare earth element concentration of the \(i\)th known point along the direction from the \(i\)th known point to the \(t\)th reference point; The variance of the concentration change components of all verification points between the \(i\)th known point and the \(t\)th reference point; The total number of all verification points between the \(i\)th known point and the \(t\)th reference point; The component of the concentration change rate in the direction from the \(i\)th known point to the \(t\)th reference point of the \(y\)th verification point of the \(t\)th reference point of the \(i\)th known point; The concentration change rate in the direction from the \(i\)th known point to the \(t\)th reference point in the concentration change network of the \(i\)th known point; Denote taking the absolute value; The Euclidean distance between the \(i\)th known point and the \(y\)th verification point of the \(t\)th reference point in the \(i\)th known point; The included angle between the line segment connecting the \(i\)th known point and the \(t\)th reference point and the line segment connecting the \(i\)th known point and the \(y\)th verification point; Denote the sine function, Denote the linear normalization function.

[0020] As a specific solution in the technical solution of this application, the processing module is further configured to obtain a third known point based on each known point; the third known point is any one of the known points; And, obtain an effective deviation coefficient based on the deviation coefficients in each direction of the third known point; the effective deviation coefficient is the deviation coefficient less than a preset value among the deviation coefficients; And, obtain a uniform adjustment coefficient of the third known point and the point to be evaluated based on each effective deviation coefficient; the uniform adjustment coefficient is at least used to characterize the degree of uniform change in rare earth element concentration of the third known point along the direction from the third known point to the point to be evaluated; And, obtain the interpolation weight of the third known point based on the uniform adjustment coefficient.

[0021] In a third aspect, this application proposes a rapid rare earth mineral exploration system based on GIS. The rapid rare earth mineral exploration system based on GIS includes: A processor, configured to obtain a divided grid based on a target mining area; the target mining area is obtained in advance; A reader, configured to obtain a plurality of known points based on the divided grid; the known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known; The processor is further configured to obtain a concentration change network of each known point based on each known point; Further, based on the concentration change network, obtain the deviation coefficients of the rare earth element concentration changes in various directions corresponding to each known point; the deviation coefficients are at least used to characterize the deviation degree of the actual change value of the rare earth element concentration at each known point in each direction from the uniform change value; Further, based on the deviation coefficients corresponding to various directions of each known point, obtain the interpolation weights of the known points; Further, based on the interpolation weights of each known point, obtain the rare earth element concentrations of each unknown point.

[0022] As a specific solution in the technical solution of the present application, the processor is further configured to obtain a first known point based on each known point; the first known point is any one of the known points; Based on the first known point, obtain multiple reference points from each known point; each reference point is distributed centered on the first known point; Further, based on the first known point and each reference point, obtain the rare earth element concentration change rates in various directions; each direction points from the first known point to the corresponding reference point; Further, based on the rare earth element concentration change rates in various directions, obtain the concentration change network of the first known point.

[0023] As a specific solution in the technical solution of the present application, the processor is further configured to obtain a second known point based on each known point; the second known point is any one of the known points; Further, based on the second known point, obtain a first reference point; the first reference point is any one of the reference points corresponding to the second known point; Further, based on the second known point and the first reference point, obtain multiple verification points; the verification points are any known points among the known points located between the second known point and the first reference point; Further, based on each verification point, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction; the first direction points from the second known point to the first reference point.

[0024] As a specific solution in the technical solution of the present application, the processor is further configured to obtain a first line segment based on the second known point and the first reference point; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; Further, based on the first line segment, obtain a first straight line and a second straight line; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; Further, each known point located between the first straight line and the second straight line is used as each verification point.

[0025] As a specific solution in the technical solution of the present application, the processor is further configured to obtain a first verification point based on each verification point; the first verification point is any one of the verification points; and, based on the second known point and the first verification point, obtain a concentration change component of the second known point and the first verification point in the first direction; and, based on the concentration change component, obtain a deviation coefficient of the rare earth element concentration change of the second known point in the first direction.

[0026] As a specific solution in the technical solution of the present application, the calculation formula for the processor to obtain the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ; wherein, represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the rare earth element concentration of the i-th known point; represents the rare earth element concentration of the y-th verification point of the t-th reference point of the i-th known point; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the cosine function.

[0027] As a specific solution in the technical solution of the present application, the calculation formula for the processor to obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component is as follows: ; wherein, represents the deviation coefficient of the rare earth element concentration change of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the variance of the concentration change components of all verification points between the i-th known point and the t-th reference point; represents the total number of all verification points between the i-th known point and the t-th reference point; represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the concentration change rate of the \(i\)-th known point in the direction pointing from the \(i\)-th known point to the \(t\)-th reference point in the concentration change network of the \(i\)-th known point; represents taking the absolute value; represents the Euclidean distance between the \(i\)-th known point and the \(y\)-th verification point of the \(t\)-th reference point among the \(i\)-th known points; represents the angle between the line segment connecting the \(i\)-th known point and the \(t\)-th reference point and the line segment connecting the \(i\)-th known point and the \(y\)-th verification point; represents the sine function, represents the linear normalization function.

[0028] As a specific solution in the technical solution of the present application, the processor is further configured to obtain a third known point based on each known point; the third known point is any one of the known points; and, based on the deviation coefficients in each direction of the third known point, obtain an effective deviation coefficient; the effective deviation coefficient is the deviation coefficient less than a preset value among the deviation coefficients; and, based on each effective deviation coefficient, obtain a uniform adjustment coefficient between the third known point and the point to be evaluated; the uniform adjustment coefficient is at least used to characterize the degree of uniform change of the rare earth element concentration in the direction from the third known point to the point to be evaluated; and, based on the uniform adjustment coefficient, obtain the interpolation weight of the third known point.

[0029] Compared with the prior art, the beneficial effects of the present application are: In the present application, by taking several reference directions at each known point, the concentration change rate of the rare earth element at the known point in each direction is calculated respectively, and whether the actual concentration change in each direction is uniform is calculated through the concentrations of several known points in the corresponding direction. The interpolation weight of each known point is obtained by using the directions with relatively uniform concentration changes. Compared with the prior art that only uses the physical distance from the known point to the point to be evaluated to obtain the interpolation weight, the interpolation weight obtained in the present application is more accurate, that is, the rare earth element concentration of the point to be evaluated obtained subsequently is also more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 is a schematic flow chart of a method for rapid exploration of rare earth minerals based on GIS proposed in an embodiment of the present application; Figure 2 is a schematic structural diagram of an electronic device for rapid exploration of rare earth minerals based on GIS proposed in an embodiment of the present application; Figure 3 is a schematic structural diagram of a system for rapid exploration of rare earth minerals based on GIS proposed in an embodiment of the present application; Figure 4Schematic diagram of setting sampling points by the snake-shaped point distribution method in dividing grid cells proposed in the embodiments of the present application; Figure 5 Schematic diagram of selecting a reference point proposed in the embodiments of the present application; Figure 6 Schematic diagram of selecting a verification point proposed in the embodiments of the present application; Figure 7 Schematic diagram of the third known point and the point to be evaluated proposed in the embodiments of the present application. Detailed implementation manners

[0031] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

[0032] Terms such as "first" and "second" in the specification of the embodiments of the present application and the above accompanying drawings are used to distinguish similar objects and do not necessarily need to be used to describe a specific order or sequence. For example, the first known point and the second known point proposed below belong to different known points. It should be understood that such known points can be interchanged under appropriate circumstances so that the embodiments described here can be implemented in an order other than that shown or described here. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or modules does not necessarily need to be limited to those steps or modules clearly listed, but may include other steps or modules not clearly listed or inherent to these processes, methods, products, or devices. The division of modules in the embodiments of the present application is only a logical division, and there may be other division methods in actual implementation. For example, multiple modules can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the shown or discussed coupling or direct coupling or communication connection between each other can be through some interfaces, and the indirect coupling or communication connection between modules can be electrical or other similar forms, which are not limited in the embodiments of the present application. And the modules or sub-modules described as separate components may or may not be physically separated, may or may not be physical modules, or may be distributed to multiple circuit modules. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of the embodiments of the present application.

[0033] To solve the technical problem in the prior art that the rare earth element concentration of the point to be evaluated cannot be accurately estimated in the prior art, the present application proposes a rapid exploration method for rare earth minerals based on GIS. Specifically, as Figure 1 shown, the rapid exploration method for rare earth minerals based on GIS includes steps S100 to S600.

[0034] Step S100: Obtain a divided grid based on the target mining area.

[0035] In this embodiment, the target mining area is obtained in advance. That is to say, in this embodiment, the target mining area can be any area where the rare earth element concentration needs to be estimated.

[0036] It should be noted that if a model needs to be established based on geochemical data, soil and rock samples on the surface need to be collected. When sampling for mineral exploration, in order to make the sampling points distributed in the mining area as uniform as possible and avoid over-concentration or omission of local sampling, the target mining area is usually divided into several grids of equal size, and points are selected for sampling within the grids. For example, the target mining area can be evenly divided into several grids of 50m×50m. The shape of the target mining area may not be completely divided into regular and equal-sized grids. When dividing, the number of regular grids should be as large as possible. For irregular areas, the sampling points therein can be adjusted according to the area. For example: if the number of sampling points in a regular grid of 50m×50m is 20, and the area of some irregular areas is S, then the number of sampling points therein can be [ ], [ ] indicates rounding the value in the brackets. Since the target mining area is usually an area with a large area, uneven terrain, and uneven soil, the point layout method within the grid can choose the serpentine layout method. As Figure 4 shown, it is a schematic diagram of setting sampling points (i.e., Figure 4 the diamond hollow points, diamond black points, and circular hollow points in it) using the serpentine layout method in a single regular grid.

[0037] Step S200: Obtain a plurality of known points based on the divided grid.

[0038] In this embodiment, the known points are the sampling points of any one grid in the divided grid (i.e., Figure 4 the square hollow points, diamond black points, or circular hollow points in it), and the rare earth element concentration of each known point is known. Samples of all sampling points are obtained within all grids of the target mining area, each sampling point is recorded as a known point, the rare earth element content of all samples is measured, and the rare earth element concentration of all known points is obtained. Measuring the concentration of rare earth elements in the test samples is a mature technology and will not be elaborated here.

[0039] It should be noted that since the distribution of rare earth elements in the land is not uniform, when directly interpolating based on the physical distances between several adjacent known points and the point to be evaluated in the sampling, it is approximately considered that the rare earth element concentration of each known point changes uniformly in all directions. Therefore, the result of the rare earth element concentration of the point to be evaluated obtained is not accurate. To perform interpolation more accurately, it is necessary to analyze the change trend of the rare earth element concentration of each known point in all directions (i.e., the concentration change network in the following text), and then obtain the change trend of the rare earth element concentration in the direction connecting each known point and the point to be evaluated. Combining with the distance between the known point and the point to be evaluated, the interpolation weight can be determined more accurately to accurately obtain the rare earth element concentration of the point to be evaluated. In the embodiments described below in this application, only a division grid of one unit is used to illustrate the embodiments of this application, which does not mean that the embodiments proposed in this application are only applicable to obtaining the concentration change network of each known point based on a division grid of one unit. It should be understood that the concentration change network of each known point can be obtained based on a division grid of one unit, or can be obtained based on a division grid of multiple units.

[0040] Step S300: Based on each known point, obtain the concentration change network of each known point.

[0041] As can be seen from the foregoing, the concentration change network refers to the change trend of the rare earth element concentration of the corresponding known point in all directions. In the embodiments of this application, any suitable method can be used to obtain the concentration change network of each known point. For example, in the embodiments of this application, step S300, based on each known point, obtaining the concentration change network of each known point, may include steps S310 to S340.

[0042] Step S310: Based on each known point, obtain a first known point.

[0043] In this embodiment, the first known point is any one of the known points. That is to say, in this embodiment, the concentration change network corresponding to each known point can be obtained by the same method as that of the first known point.

[0044] Step S320: Based on the first known point, obtain multiple reference points from each known point.

[0045] It should be clear that to obtain the concentration change network corresponding to the first known point, it is first necessary to select reference points in several directions, and then compare the concentration changes between the reference points and the known points to obtain the concentration changes of each known point in multiple directions. In the embodiments of the present application, the respective reference points are distributed centered on the first known point. In this embodiment, a random selection method can be used to obtain multiple reference points from each known point. As can be seen from the foregoing, since the method selected in this embodiment is the serpentine sampling method during sampling, in order to make the direction of analyzing the concentration change more conform to the sampling situation, all the extreme points of the serpentine sampling method (i.e., the rhombus hollow points as shown in Figure 4 can be selected as reference points, and the extreme points are the points closest to the boundary in the trajectory of the serpentine sampling. In this embodiment, as shown in Figure 4 , it is assumed that the circular point is the first known point.

[0046] Step S330: Based on the first known point and each reference point, obtain the concentration change rate of rare earth elements in each direction.

[0047] In this embodiment, each direction points from the first known point to the corresponding reference point. As can be seen from Figure 4 , the number of reference points is 4, and a direction can be formed by pointing from the first known point to any one of the reference points. That is to say, in this embodiment, the number of reference points is equal to the number of directions in the concentration change network. That is, as shown in Figure 5 , 4 reference points can form 4 directions with the first known point. In other words, if there are 4 directions, then in step S340, it is necessary to obtain the concentration change rate of rare earth elements of the first known point in 4 directions.

[0048] In the embodiments of the present application, any suitable method can be used to obtain the concentration change rate of rare earth elements in each direction based on the first known point and each reference point. For example, after obtaining the reference points, the concentration differences between the known points and each reference point can be compared, and then the concentration change situations of the known points in multiple directions can be obtained. Taking the i-th known point (i.e., the first known point in the foregoing) as an example, the concentration differences between this known point and each reference point are compared. The calculation formula for the concentration change rate of the rare earth elements from the i-th known point to the t-th reference point can be: ; where represents the concentration change rate in the direction from the i-th known point to the t-th reference point in the concentration change network of the i-th known point; represents the rare earth element concentration of the i-th known point; represents the rare earth element concentration of the t-th reference point; represents the Euclidean distance between the i-th known point and the t-th reference point; Indicates the degree to which the rare earth element concentration of the \(i\)-th known point is higher than that of the \(t\)-th reference point. When the rare earth element concentration of the \(t\)-th reference point is greater, the concentration change rate is negative.

[0049] Step S340: Based on the rare earth element concentration change rates in each direction, obtain the concentration change network of the first known point.

[0050] It should be clear that in this embodiment, the rare earth element concentration change rates of the first known point in each direction are the concentration change network of the first known point. Specifically, as Figure 5 shown, the connection lines between the \(i\)-th known point (i.e., the first known point in the above text) and each reference point can be denoted as each side of the \(i\)-th known point. The value of each side is the rare earth element concentration change rate from the \(i\)-th known point to the corresponding reference point, and all the sides of the \(i\)-th known point form the concentration change network of the \(i\)-th known point.

[0051] Step S400: Based on the concentration change network, obtain the deviation coefficients of the rare earth element concentration changes in each direction of the corresponding known point.

[0052] It should be clear that when the prior art only obtains the interpolation weights based on the distance for interpolation, it regards the rare earth element concentration changes in each direction of all known points as uniform changes. In this application, a concentration change network is constructed, and the rare earth element concentration change of each known point is split into multiple directions to discuss the concentration change situations in multiple directions separately in subsequent analysis, which more finely divides the rare earth element concentration change of each known point. Therefore, the subsequent measurement of the interpolation weights is more accurate.

[0053] In this embodiment, the deviation coefficient is at least used to characterize the deviation degree between the actual change value and the uniform change value of the rare earth element concentration of the corresponding known point in each direction.

[0054] It should be noted that in this embodiment, any reasonable method can be adopted to obtain the deviation coefficients of the rare earth element concentration changes in each direction of the corresponding known point based on the concentration change network. For example: Step S400, to obtain the deviation coefficients of the rare earth element concentration changes in each direction of the corresponding known point based on the concentration change network, it may include steps S410 to S440.

[0055] Step S410: Based on each known point, obtain a second known point.

[0056] In this embodiment, the second known point can be any one of the known points. For example: the second known point can be the first known point. That is to say, in this embodiment, each known point can adopt the same method as the second known point to obtain the deviation coefficients of the rare earth element concentration changes in each direction.

[0057] Step S420: Obtain a first reference point based on the second known point.

[0058] In this embodiment, the first reference point is any one of the reference points corresponding to the second known point. As can be seen from the foregoing, a direction can be formed between the second known point and any one reference point. That is to say, in this embodiment, the deviation coefficients of the rare earth element concentration changes in each direction of the second known point can all be obtained by the same method as the deviation coefficient of the rare earth element concentration change in the direction from the second known point to the first reference point.

[0059] Step S430: Obtain a plurality of verification points based on the second known point and the first reference point.

[0060] As can be seen from the foregoing, since when calculating the rare earth element concentration change rate from the i-th known point to the t-th reference point, the concentration change between the i-th known point and the t-th reference point is regarded as a uniform change to obtain the rare earth element concentration change rate, but the actual situation may not be uniform. Therefore, it is necessary to analyze whether the calculation result of the rare earth element concentration change rate is reliable according to several other known points (i.e., verification points) between these two known points. In this embodiment, the verification points can be any known points among all the known points that are located between the second known point and the first reference point.

[0061] In this embodiment, a plurality of verification points can be randomly obtained from all the known points. In another embodiment of the present application, step S430, obtaining a plurality of verification points based on the second known point and the first reference point, may include steps S431 to S433.

[0062] Step S431: Obtain a first line segment based on the second known point and the first reference point.

[0063] In this embodiment, the second known point is one end of the first line segment, and the first reference point is the other end of the first line segment. In a specific embodiment of the present application, the first line segment can be as Figure 6 shown by line segment 1 in, where the circular hollow point in line segment 1 is the second known point, and the diamond hollow point in line segment 1 is the first reference point.

[0064] Step S432: Obtain a first straight line and a second straight line based on the first line segment.

[0065] In this embodiment, the first straight line passes through the second known point and is perpendicular to the first line segment, that is, the straight line 2 as Figure 6 shown. The second straight line passes through the first reference point and is perpendicular to the first line segment, that is, the straight line 3 as Figure 6 shown.

[0066] Step S433: Use each known point located between the first straight line and the second straight line as each verification point.

[0067] Specifically, as Figure 6 shown, a gray area is formed between the first straight line and the second straight line. That is to say, as long as the known points falling within the gray area (i.e., the diamond black dots in the gray area as Figure 6 shown) can be used as the verification points between the second known point and the first reference point.

[0068] Step S440: Based on each verification point, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction.

[0069] In this embodiment, the first direction points from the second known point to the first reference point. As can be seen from the foregoing, the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration of the corresponding known point in the first direction from the uniform change value. In this embodiment, based on each verification point, any suitable method can be used to obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction. For example, step S440, based on each verification point, obtaining the deviation coefficient of the rare earth element concentration change of the second known point in the first direction, may include steps S441 to S443.

[0070] Step S441: Based on each verification point, obtain a first verification point.

[0071] In this embodiment, the first verification point is any one of the verification points. As can be seen from step S442, the concentration change component in the first direction needs to be calculated for each verification point. That is to say, in this embodiment, the concentration change components of each verification point in the first direction can all refer to the first verification point, and will not be elaborated later.

[0072] Step S442: Based on the second known point and the first verification point, obtain the concentration change components of the second known point and the first verification point in the first direction.

[0073] In the embodiments of the present application, any reasonable method can be used to obtain the concentration change components of the second known point and the first verification point in the first direction based on the second known point and the first verification point. For example: step S442, the calculation formula for obtaining the concentration change components of the second known point and the first verification point in the first direction based on the second known point and the first verification point can be as follows: ; where The component of the rate of change of concentration in the direction from the $i$-th known point to the $t$-th reference point for the $y$-th verification point of the $t$-th reference point of the $i$-th known point; Represents the rare earth element concentration of the $i$-th known point; Represents the rare earth element concentration of the $y$-th verification point of the $t$-th reference point of the $i$-th known point; Represents the Euclidean distance between the $i$-th known point and the $y$-th verification point of the $t$-th reference point in the $i$-th known point; Represents the angle between the line segment connecting the $i$-th known point and the $t$-th reference point and the line segment connecting the $i$-th known point and the $y$-th verification point (i.e., the angle formed by line segment 1 and line segment 4 as shown in Figure 6 ); ); Represents the cosine function.

[0074] Step S443: Based on the concentration change component, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction.

[0075] It should be clear that if the concentration change components corresponding to each verification point are closer to the rate of change of the rare earth element concentration of the second known point in the first direction, it means that the rare earth element concentration of the second known point changes more uniformly in the first direction. That is, when obtaining the interpolation weight subsequently, it is more necessary to increase the weight value in this direction. On the contrary, when obtaining the interpolation weight subsequently, it is more necessary to decrease the weight value in this direction.

[0076] In this embodiment, any reasonable method can be used to obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component. For example, in an embodiment of the present application, the variance or standard deviation of the concentration change components corresponding to the second known point and the rate of change of the rare earth element concentration in the first direction can be used as the deviation coefficient. In another embodiment of the present application, in step S443, the calculation formula for obtaining the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component is as follows: ; Where, Represents the deviation coefficient of the change in rare earth element concentration from the $i$-th known point to the $t$-th reference point; Represents the variance of the concentration change components of all verification points between the $i$-th known point and the $t$-th reference point; Represents the total number of all verification points between the $i$-th known point and the $t$-th reference point; The component of the rate of change of concentration in the direction from the $i$-th known point to the $t$-th reference point for the $y$-th verification point of the $t$-th reference point of the $i$-th known point; represents the concentration change rate in the direction from the \(i\)-th known point to the \(t\)-th reference point in the concentration change network of the \(i\)-th known point; represents taking the absolute value; represents the Euclidean distance between the \(i\)-th known point and the \(y\)-th verification point of the \(t\)-th reference point among the \(i\)-th known points; represents the angle between the line segment connecting the \(i\)-th known point and the \(t\)-th reference point and the line segment connecting the \(i\)-th known point and the \(y\)-th verification point; when the denominator gets 0, it is replaced by 0.1; represents the sine function; represents the linear normalization function.

[0077] Step S500: Based on the deviation coefficients corresponding to each known point in each direction and the point to be evaluated, obtain the interpolation weights of each known point.

[0078] In this embodiment, the point to be evaluated is obtained in advance. The point to be evaluated can be any point with unknown rare earth element concentration in the divided grid. In the embodiments of the present application, any reasonable method can be used to obtain the interpolation weights of each known point based on the deviation coefficients corresponding to each known point in each direction and the point to be evaluated. For example: Step S500, based on the deviation coefficients corresponding to each known point in each direction and the point to be evaluated, obtaining the interpolation weights of each known point may include Step S510 to Step S540.

[0079] Step S510: Based on each known point, obtain a third known point.

[0080] In this embodiment, the third known point is any one of the known points. As can be seen from the foregoing, when calculating the rare earth element concentration of the point to be evaluated, interpolation weights need to be assigned to each known point. That is to say, in this embodiment, the method for obtaining the interpolation weights of each known point can adopt the method for obtaining the interpolation weights of the third known point.

[0081] Step S520: Based on the deviation coefficients in each direction of the third known point, obtain the effective deviation coefficients.

[0082] It should be noted that the change in the rare earth element concentration of the third known point in each direction is not uniform. It is easy to understand that the larger the deviation coefficient of the third known point in one direction, the more uneven the change in the rare earth element concentration of the third known point in that direction. If there is a large error in the inverse distance interpolation weight result obtained through the third known point in this direction, and for the direction with a smaller deviation evaluation, the obtained interpolation result is more accurate. By presetting a threshold to screen out points with too large deviation evaluations, the situation of large deviations during interpolation can be avoided. In the embodiments of the present application, the deviation coefficient with a smaller value (i.e., the effective deviation coefficient) can be selected to participate in the subsequent acquisition of the interpolation weight. That is to say, in this embodiment, the effective deviation coefficient is the deviation coefficient less than the preset value among all the deviation coefficients. In this embodiment, the size of the preset value is not limited and can be set according to requirements. For example, the preset value can be 0.6 or 0.7, etc.

[0083] Step S530: Based on each effective deviation coefficient, obtain the uniform adjustment coefficient between the third known point and the point to be evaluated.

[0084] It should be clear that when using the inverse distance interpolation weight to obtain the interpolation, the search radius of the interpolation needs to be set first. In the present application, when obtaining the interpolation, the nine-grid range centered on the grid where the point to be evaluated is located is selected as the search radius. According to the previous analysis, the smaller the deviation coefficient of a known point in a certain direction, the more uniform the change in the rare earth element concentration of the known point along that direction, and then the larger the interpolation weight of the known point along that direction. By using the deviation coefficients of each known point in each direction and combining the Euclidean distance between each known point and the point to be evaluated, the interpolation weight of each known point for the point to be evaluated can be obtained.

[0085] In this embodiment, the uniform adjustment coefficient is at least used to characterize the degree of uniform change in the rare earth element concentration of the third known point along the direction from the third known point to the point to be evaluated. In the embodiments of the present application, any reasonable method can be adopted to obtain the uniform adjustment coefficient between the third known point and the point to be evaluated based on each effective deviation coefficient. For example: in step S530, the calculation formula for obtaining the uniform adjustment coefficient between the third known point and the point to be evaluated based on each effective deviation coefficient is as follows: ; where The uniform adjustment coefficient of the i-th known point (i.e., the third known point); represents the number of effective deviation coefficients of the i-th known point; represents the angle of the m-th direction of the line connecting the point to be evaluated and the i-th known point (i.e., the angle Figure 7 shown as ); Denote the deviation evaluation in the m-th direction of the i-th known point. Denote the inverse proportional normalization function. As the weight, calculate the weighted mean of the deviation evaluations, and take the inverse of the weighted mean as the uniform adjustment coefficient. The larger it is, the smaller the angle between the point to be evaluated and the m-th effective edge of the i-th known point. Here, +1 is to ensure the object is a positive number, so that the final result of inverse proportional normalization is correctly mapped to the range of [0, 1]. Figure 7 The point corresponding to the hollow triangle in

[0086] Step S540: Based on the uniform adjustment coefficient, obtain the interpolation weight of the third known point.

[0087] It should be clear that the uniform adjustment coefficient reflects the degree of uniformity of the change of the i-th known point along the connection direction of the point to be evaluated. The larger this value is, the greater the reference degree of the i-th known point for interpolation, that is, the larger the interpolation weight. Use the uniform adjustment coefficient to adjust the method of calculating the interpolation weight usually using distance, and obtain the interpolation weight of each known point. The interpolation weight of the i-th known point is calculated as follows: ; In the formula: Denote the interpolation weight of the i-th known point for the point to be evaluated; Denote the uniform adjustment coefficient of the i-th known point; Denote the Euclidean distance between the i-th known point and the point to be evaluated; Denote the exponential parameter, which affects the smoothness of interpolation. Here, k = 2; Denote the linear normalization function.

[0088] Step S600: Based on the interpolation weights of each known point, obtain the rare earth element concentration of the point to be evaluated.

[0089] In this embodiment, several points to be evaluated can be set, calculate the interpolation weights of all known points and each point to be evaluated, perform interpolation on each point to be evaluated, and obtain the rare earth element concentrations of all points to be evaluated.

[0090] An embodiment of the GIS-based rapid exploration method for rare earth minerals proposed in this application calculates the change rate of rare earth element concentration at each known point in several reference directions respectively at each known point, and calculates whether the actual concentration change in each direction is uniform through the concentrations of several known points in the corresponding direction. The interpolation weights of each known point are obtained by using the directions with relatively uniform concentration changes. Compared with the prior art that only uses the physical distance from the point to be evaluated to obtain the interpolation weight, the interpolation weights obtained in this application are more accurate, that is, the rare earth element concentration of the point to be evaluated obtained subsequently is also more accurate.

[0091] After introducing the GIS-based rapid exploration method for rare earth minerals proposed in the embodiment of this application, the following introduces the embodiment of the GIS-based rapid exploration electronic device proposed in this application. As Figure 2 shown, the GIS-based rapid exploration electronic device 10 includes: A processing module 11, configured to obtain a divided grid based on a target mining area; the target mining area is obtained in advance; A reading module 12, configured to obtain a plurality of known points based on the divided grid; the known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known; The processing module 11 is further configured to obtain a concentration change network of each known point based on each known point; And, based on the concentration change network, obtain the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration at the corresponding known point in each direction from the uniform change value; And, based on the deviation coefficients corresponding to each direction of the known point, obtain the interpolation weight of the known point; And, based on the interpolation weights of each known point, obtain the rare earth element concentration of each unknown point.

[0092] As a specific embodiment in this application, the processing module 11 is further configured to obtain a first known point based on each known point; the first known point is any one of the known points; Based on the first known point, obtain a plurality of reference points from each known point; each reference point is distributed centered on the first known point; And, based on the first known point and each reference point, obtain the change rate of the rare earth element concentration in each direction; each direction points from the first known point to the corresponding reference point; And, based on the change rate of the rare earth element concentration in each direction, obtain the concentration change network of the first known point.

[0093] As a specific embodiment in the present application, the processing module 11 is further configured to obtain a second known point based on each known point; the second known point is any one of the known points; and obtain a first reference point based on the second known point; the first reference point is any one of the reference points corresponding to the second known point; and obtain a plurality of verification points based on the second known point and the first reference point; the verification point is any known point among the known points located between the second known point and the first reference point; and obtain a deviation coefficient of the change in the rare earth element concentration of the second known point in the first direction based on each verification point; the first direction points from the second known point to the first reference point.

[0094] As a specific embodiment in the present application, the processing module 11 is further configured to obtain a first line segment based on the second known point and the first reference point; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; and obtain a first straight line and a second straight line based on the first line segment; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; and use each known point located between the first straight line and the second straight line as each verification point.

[0095] As a specific embodiment in the present application, the processing module 11 is further configured to obtain a first verification point based on each verification point; the first verification point is any one of the verification points; and obtain a concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point; and obtain a deviation coefficient of the change in the rare earth element concentration of the second known point in the first direction based on the concentration change component.

[0096] As a specific embodiment in the present application, the formula for the processing module 11 to obtain the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ; where represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the rare earth element concentration of the i-th known point; The rare earth element concentration of the y-th verification point of the t-th reference point of the i-th known point; The Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; The included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; Represents the cosine function.

[0097] As a specific embodiment in the present application, the processing module 11 obtains the calculation formula for the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component as follows: ; Wherein, Represents the deviation coefficient of the rare earth element concentration change of the i-th known point along the direction from the i-th known point to the t-th reference point; Represents the variance of the concentration change components of all verification points between the i-th known point and the t-th reference point; Represents the total number of all verification points between the i-th known point and the t-th reference point; Represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; Represents the concentration change rate in the direction from the i-th known point to the t-th reference point in the concentration change network of the i-th known point; Represents taking the absolute value; Represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; Represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; Represents the sine function, Represents the linear normalization function.

[0098] As a specific embodiment in the present application, the processing module 11 is further configured to obtain a third known point based on each known point; the third known point is any one of the known points; And, obtain an effective deviation coefficient based on the deviation coefficients in each direction of the third known point; the effective deviation coefficient is the deviation coefficient less than a preset value among the deviation coefficients; And, obtain a uniform adjustment coefficient of the third known point and the point to be evaluated based on each effective deviation coefficient; the uniform adjustment coefficient is at least used to characterize the degree of uniform change of the rare earth element concentration of the third known point along the direction from the third known point to the point to be evaluated; Further, based on the uniform adjustment coefficient, obtain the interpolation weight of the third known point.

[0099] In the embodiment of the GIS-based rapid exploration electronic device for rare earth minerals proposed in this application, by taking several reference directions at each known point, calculate the change rate of rare earth element concentration of the known point in each direction, and calculate whether the actual concentration change in each direction is uniform through the concentrations of several known points in the corresponding direction. Use the directions with relatively uniform concentration changes to obtain the interpolation weight of each known point. Compared with the prior art that only uses the physical distance from the point to be evaluated to obtain the interpolation weight, the interpolation weight obtained in this application is more accurate, that is, the rare earth element concentration of the point to be evaluated obtained subsequently is also more accurate.

[0100] After introducing the GIS-based rapid exploration electronic device for rare earth minerals proposed in the embodiment of this application, the embodiment of the GIS-based rapid exploration system for rare earth minerals proposed in this application will be introduced below. As Figure 3 shown, the GIS-based rapid exploration electronic system 20 for rare earth minerals includes: A processor 21, configured to obtain a divided grid based on a target mining area; the target mining area is obtained in advance; A reader 22, configured to obtain a plurality of known points based on the divided grid; the known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known; The processor 21 is further configured to obtain a concentration change network of each known point based on each known point; Further, based on the concentration change network, obtain the deviation coefficient of the rare earth element concentration change in each direction of the corresponding known point; the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration of the corresponding known point in each direction from the uniform change value; Further, based on the deviation coefficients corresponding to each direction of the known points, obtain the interpolation weight of the known point; Further, based on the interpolation weights of each known point, obtain the rare earth element concentration of each unknown point.

[0101] As a specific embodiment of this application, the processor 21 is further configured to obtain a first known point based on each known point; the first known point is any one of the known points; Based on the first known point, obtain a plurality of reference points from each known point; each reference point is distributed centered on the first known point; Further, based on the first known point and each reference point, obtain the change rate of the rare earth element concentration in each direction; each direction points from the first known point to the corresponding reference point; And, based on the rate of change of the rare earth element concentration in each direction, obtain the concentration change network of the first known point.

[0102] As a specific embodiment in this application, the processor 21 is further configured to obtain a second known point based on each known point; the second known point is any one of the known points; And, based on the second known point, obtain a first reference point; the first reference point is any one of the reference points corresponding to the second known point; And, based on the second known point and the first reference point, obtain a plurality of verification points; the verification points are any known points among the known points located between the second known point and the first reference point; And, based on each verification point, obtain the deviation coefficient of the change in the rare earth element concentration of the second known point in the first direction; the first direction points from the second known point to the first reference point.

[0103] As a specific embodiment in this application, the processor 21 is further configured to obtain a first line segment based on the second known point and the first reference point; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; And, based on the first line segment, obtain a first straight line and a second straight line; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; And, use each known point located between the first straight line and the second straight line as each verification point.

[0104] As a specific embodiment in this application, the processor 21 is further configured to obtain a first verification point based on each verification point; the first verification point is any one of the verification points; And, based on the second known point and the first verification point, obtain the concentration change component of the second known point and the first verification point in the first direction; And, based on the concentration change component, obtain the deviation coefficient of the change in the rare earth element concentration of the second known point in the first direction.

[0105] As a specific embodiment in this application, the formula for the processor 21 to obtain the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ; Wherein, Denote the component of the concentration change rate in the direction from the $i$-th known point to the $t$-th reference point for the $y$-th verification point of the $t$-th reference point of the $i$-th known point; Denote the rare earth element concentration of the $i$-th known point; Denote the rare earth element concentration of the $y$-th verification point of the $t$-th reference point of the $i$-th known point; Denote the Euclidean distance between the $i$-th known point and the $y$-th verification point of the $t$-th reference point in the $i$-th known point; Denote the angle between the line segment connecting the $i$-th known point and the $t$-th reference point and the line segment connecting the $i$-th known point and the $y$-th verification point; Denote the cosine function.

[0106] As a specific embodiment in the present application, the processor 21 obtains the calculation formula of the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component as follows: ; wherein, Denote the deviation coefficient of the rare earth element concentration change in the direction from the $i$-th known point to the $t$-th reference point of the $i$-th known point; Denote the variance of the concentration change components of all verification points between the $i$-th known point and the $t$-th reference point; Denote the total number of all verification points between the $i$-th known point and the $t$-th reference point; Denote the component of the concentration change rate in the direction from the $i$-th known point to the $t$-th reference point for the $y$-th verification point of the $t$-th reference point of the $i$-th known point; Denote the concentration change rate in the direction from the $i$-th known point to the $t$-th reference point in the concentration change network of the $i$-th known point; Denote taking the absolute value; Denote the Euclidean distance between the $i$-th known point and the $y$-th verification point of the $t$-th reference point in the $i$-th known point; Denote the angle between the line segment connecting the $i$-th known point and the $t$-th reference point and the line segment connecting the $i$-th known point and the $y$-th verification point; Denote the sine function, Denote the linear normalization function.

[0107] As a specific embodiment in the present application, the processor 21 is further configured to obtain a third known point based on each known point; the third known point is any one of the known points; and, obtain an effective deviation coefficient based on the deviation coefficients in each direction of the third known point; the effective deviation coefficient is the deviation coefficient less than a preset value among the deviation coefficients; Further, based on each effective deviation coefficient, a uniform adjustment coefficient of the third known point and the point to be evaluated is obtained; the uniform adjustment coefficient is at least used to characterize the degree of uniform change in the rare earth element concentration of the third known point along the direction from the third known point to the point to be evaluated. Further, based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0108] In the embodiment of the GIS-based rapid exploration system for rare earth minerals proposed in this application, by taking a number of reference directions at each known point, the change rate of the rare earth element concentration of the known point in each direction is calculated respectively, and whether the actual concentration change in each direction is uniform is calculated through the concentrations of a number of known points in the corresponding direction. The interpolation weight of each known point is obtained by using the directions with relatively uniform concentration changes. Compared with the prior art that only uses the physical distance from the point to be evaluated to obtain the interpolation weight, the interpolation weight obtained in this application is more accurate, which also makes the rare earth element concentration of the point to be evaluated obtained subsequently more accurate.

[0109] It should be clear that the computer-readable storage medium in this application includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory, static random access memory, dynamic random access memory, other types of random access memory, read-only memory, electrically erasable programmable read-only memory, flash memory or other memory technologies, compact disc read-only memory, digital versatile disc or other optical storage, magnetic cassette tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transient computer-readable media such as modulated data signals and carrier waves.

[0110] In the above embodiments, the descriptions of the respective embodiments have their own emphases. For the parts not detailed in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0111] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the methods, apparatuses, and devices described above can refer to the corresponding processes in the foregoing method embodiments, and will not be elaborated herein.

[0112] In several embodiments provided by the embodiments of the present application, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For example, the division of the modules is only a logical function division. In actual implementation, there may be other division methods. For example, multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of the devices or modules can be in electrical, mechanical, or other forms.

[0113] The modules described as separate components may or may not be physically separated. The components displayed as modules may or may not be physical modules, that is, they can be located in one place or distributed to multiple network modules. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0114] In addition, in each embodiment of the embodiments of the present application, the various functional modules can be integrated in a processing module, or each module can exist physically alone, or two or more modules can be integrated in one module. The above integrated modules can be implemented in the form of hardware or in the form of software functional modules. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it can be stored in a computer-readable storage medium.

[0115] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product.

[0116] The computer program product includes one or more computer instructions. When the computer program is loaded and executed on a computer, the processes or functions described in the embodiments of the present application are generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by wire (such as coaxial cable, optical fiber, digital subscriber line) or wirelessly (such as infrared, wireless, microwave, etc.). The computer-readable storage medium may be any available medium that a computer can store or a data storage device such as a server or data center that includes one or more integrated available media. The available medium may be a magnetic medium (such as a floppy disk, hard disk, magnetic tape), an optical medium (such as a digital video disc), or a semiconductor medium (such as a solid-state disk (SSD)).

[0117] Although the embodiments of the present application have been shown and described, those of ordinary skill in the art can understand that various changes, modifications, substitutions, and variations can be made to these embodiments without departing from the principles of the present application.

Claims

1. A rapid exploration method for rare earth minerals based on GIS, characterized in that Including: Based on the target mining area, obtain a divided grid; The target mining area is obtained in advance; Based on the divided grid, obtain multiple known points; The known points are sampling points of any grid in the divided grid, and the rare earth element concentration of each known point is known; Based on each known point, obtain the concentration change network of each known point; Based on the concentration change network, obtain the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is at least used to characterize the deviation degree of the actual change value of the rare earth element concentration at the corresponding known point in each direction from the uniform change value; Based on the deviation coefficients corresponding to each direction of each known point and the point to be evaluated, obtain the interpolation weights of each known point; the point to be evaluated is obtained in advance; Based on the interpolation weights of each known point, obtain the rare earth element concentration of the point to be evaluated.

2. The rapid exploration method for rare earth minerals based on GIS according to claim 1, wherein The obtaining the concentration change network of each known point based on each known point includes: Based on each known point, obtain a first known point; the first known point is any one of the known points; Based on the first known point, obtain multiple reference points from each known point; each reference point is distributed centered on the first known point; Based on the first known point and each reference point, obtain the rare earth element concentration change rate in each direction; each direction points from the first known point to the corresponding reference point; Based on the rare earth element concentration change rate in each direction, obtain the concentration change network of the first known point.

3. The rapid exploration method of rare earth minerals based on GIS according to claim 1, characterized in that, The obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point based on the concentration change network includes: Based on each known point, obtain a second known point; the second known point is any one of the known points; Based on the second known point, obtain a first reference point; the first reference point is any one of the reference points corresponding to the second known point; Based on the second known point and the first reference point, obtain multiple verification points; the verification points are any known points among the known points located between the second known point and the first reference point; Based on each verification point, obtain the deviation coefficient of the rare earth element concentration change in the first direction of the second known point; the first direction points from the second known point to the first reference point.

4. The rapid exploration method of rare earth minerals based on GIS according to claim 3, characterized in that The obtaining multiple verification points based on the second known point and the first reference point includes: Based on the second known point and the first reference point, obtain a first line segment; the second known point is one end of the first line segment; the first reference point is the other end of the first line segment; Based on the first line segment, obtain a first straight line and a second straight line; the first straight line passes through the second known point and is perpendicular to the first line segment; the second straight line passes through the first reference point and is perpendicular to the first line segment; Take each known point located between the first straight line and the second straight line as each verification point.

5. The rapid exploration method of rare earth minerals based on GIS according to claim 3, characterized in that The obtaining the deviation coefficient of the rare earth element concentration change in the first direction of the second known point based on each verification point includes: Based on each verification point, obtain the first verification point; the first verification point is any one of the verification points; Based on the second known point and the first verification point, obtain the concentration change component of the second known point and the first verification point in the first direction; Based on the concentration change component, obtain the deviation coefficient of the rare earth element concentration change of the second known point in the first direction.

6. The rapid exploration method of rare earth minerals based on GIS according to claim 5, characterized in that, The calculation formula for obtaining the concentration change component of the second known point and the first verification point in the first direction based on the second known point and the first verification point is as follows: ； Among them, represents the concentration change rate component of the y-th verification point of the t-th reference point of the i-th known point in the direction from the i-th known point to the t-th reference point; represents the rare earth element concentration of the i-th known point; represents the rare earth element concentration of the y-th verification point of the t-th reference point of the i-th known point; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the cosine function.

7. The rapid exploration method of rare earth minerals based on GIS according to claim 6, characterized in that The calculation formula for obtaining the deviation coefficient of the rare earth element concentration change of the second known point in the first direction based on the concentration change component is as follows: ； Among them, represents the deviation coefficient of the rare earth element concentration change of the i-th known point along the direction from the i-th known point to the t-th reference point; represents the variance of the concentration change components of all verification points between the i-th known point and the t-th reference point; represents the total number of all verification points between the i-th known point and the t-th reference point; represents the component of the concentration change rate of the y-th verification point of the t-th reference point of the i-th known point along the direction from the i-th known point to the t-th reference point; represents the concentration change rate of the i-th known point in the concentration change network of the i-th known point along the direction to the t-th reference point; represents taking the absolute value; represents the Euclidean distance between the i-th known point and the y-th verification point of the t-th reference point in the i-th known point; represents the included angle between the line segment connecting the i-th known point and the t-th reference point and the line segment connecting the i-th known point and the y-th verification point; represents the sine function, represents the linear normalization function.

8. The rapid exploration method for rare earth minerals based on GIS according to any one of claims 1 to 7, characterized in that Obtaining the interpolation weights of each known point based on the deviation coefficients corresponding to each direction of each known point and the point to be evaluated includes: Based on each known point, obtain the third known point; the third known point is any one of the known points; Based on the deviation coefficients in each direction of the third known point, obtain the effective deviation coefficient; the effective deviation coefficient is the deviation coefficient less than the preset value among the deviation coefficients; Based on each effective deviation coefficient, obtain the uniform adjustment coefficient of the third known point and the point to be evaluated; the uniform adjustment coefficient is at least used to characterize the degree of uniform change of the rare earth element concentration of the third known point in the direction from the third known point to the point to be evaluated; Based on the uniform adjustment coefficient, obtain the interpolation weight of the third known point.

9. The rapid exploration electronic equipment for rare earth minerals based on GIS is characterized in that, Including: A processing module, configured to obtain a divided grid based on the target mining area; The target mining area is obtained in advance; A reading module, configured to obtain a plurality of known points based on the divided grid; the known points are the sampling points of any one grid in the divided grid, and the rare earth element concentration of each known point is known; The processing module is further configured to obtain the concentration change network of each known point based on each known point; And, based on the concentration change network, obtain the deviation coefficients of the rare earth element concentration changes in each direction of the corresponding known point; The deviation coefficient is at least used to characterize the deviation degree of the actual change value and the uniform change value of the rare earth element concentration of the corresponding known point in each direction; And, based on the deviation coefficients corresponding to each direction of the known point, obtain the interpolation weight of the known point; And, based on the interpolation weights of each known point, obtain the rare earth element concentrations of each unknown point.

10. The rapid exploration system for rare earth minerals based on GIS is characterized in that, Including: A processor, configured to obtain a divided grid based on the target mining area; The target mining area is obtained in advance; A reader, configured to obtain a plurality of known points based on the divided grid; the known points are the sampling points of any one grid in the divided grid, and the rare earth element concentration of each known point is known; The processor is further configured to obtain the concentration change network of each known point based on each known point; And, based on the concentration change network, obtain the deviation coefficients of the rare earth element concentration changes in each direction of the corresponding known point; The deviation coefficient is at least used to characterize the deviation degree of the actual change value and the uniform change value of the rare earth element concentration of the corresponding known point in each direction; And, based on the deviation coefficients corresponding to each direction of the known points, obtain the interpolation weights of the known points; And, based on the interpolation weights of each known point, obtain the rare earth element concentrations of each unknown point.

Images