GIS-based rapid exploration method, electronic equipment and system for rare earth minerals

Abstract

This application discloses a GIS-based rapid exploration method, electronic equipment, and system for rare earth minerals, relating to the field of data processing technology. The GIS-based rapid exploration method for rare earth minerals includes: obtaining a grid based on the target mining area; obtaining multiple known points based on the grid; obtaining a concentration change network for each known point based on each known point; obtaining a deviation coefficient of the rare earth element concentration change in each direction of the corresponding known point based on the concentration change network; obtaining an interpolation weight for each known point based on the corresponding deviation coefficient in each direction of each known point and the point to be evaluated; and obtaining the rare earth element concentration of the point to be evaluated based on the interpolation weight for each known point. Compared to the prior art that only uses the physical distance to the point to be evaluated to obtain the interpolation weight, the interpolation weight obtained by the deviation coefficient in this application is more accurate, which means that the rare earth element concentration of the point to be evaluated obtained subsequently is also more accurate.

Description

Technical Field

[0001] The present application relates to the field of data processing technology, specifically to a GIS-based rapid exploration method, electronic equipment, and system for rare earth minerals. Background Art

[0002] Rare earth minerals refer to ores containing rare earth elements (REEs). Geochemical exploration is a common method for rare earth mineral exploration. Surface soil and rock samples are collected, analyzed for REE content, and REE distribution maps are created to identify ore bodies and locate areas of rare earth mineral concentration. Rare earth elements are typically present in the Earth's crust at low concentrations and are relatively dispersed, often presenting as small ore bodies or mineral points. Therefore, accurate identification of potential mining areas is crucial to provide a scientific basis for resource development and utilization. Geographic Information Systems (GIS) can currently be used to build models based on geochemical data and, combined with known mineral data, assess the resource potential of mining areas. Mineral resource exploration typically involves sampling observations, and interpolation methods are often required to analyze the mineral distribution across an entire region. GIS-based spatial distribution analysis of geochemical exploration data models and, using inverse distance weighted interpolation, can be used to infer the REE concentration distribution in unsampled areas based on the known geochemical data distribution. This enriches the REE concentration distribution dataset and allows for more accurate identification of potential mining areas.

[0003] Based on the rare earth element concentrations in known areas, GIS is used to estimate the rare earth element concentrations in unknown areas through inverse distance weighted interpolation. The basic idea is that the closer the distance between the known area and the unknown area, the greater the influence of the rare earth element concentration in the known area on the rare earth element concentration in the unknown area. When the rare earth elements are evenly distributed in various areas, inverse distance weighted interpolation can obtain more accurate estimation results. However, in reality, the distribution of rare earth elements in the land is not uniform. Therefore, when the inverse distance weighted interpolation method is used only based on the physical distance between the known point and the point to be evaluated, the result obtained is not accurate enough. Summary of the Invention

[0004] The purpose of this application is to provide a GIS-based rapid exploration method, electronic equipment and system for rare earth minerals to solve the technical problem in the existing technology that the rare earth element concentration at the point to be evaluated cannot be accurately estimated.

[0005] To achieve the above objectives, this application provides the following technical solutions:

[0006] In the first aspect, the present application proposes a GIS-based rapid exploration method for rare earth minerals, which includes:

[0007] Based on the target mining area, a division grid is obtained; the target mining area is obtained in advance;

[0008] Based on the divided grid, a plurality of known points are obtained; the known points are sampling points of any one of the divided grids, and the rare earth element concentration of each known point is known;

[0009] Based on each known point, obtain the concentration change network of each known point;

[0010] Based on the concentration change network, a deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point is obtained; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value;

[0011] Based on the deviation coefficients corresponding to the respective directions of the respective known points and the points to be evaluated, the interpolation weights of the respective known points are obtained in advance;

[0012] Based on the interpolation weights of the known points, the rare earth element concentration of the point to be evaluated is obtained.

[0013] As a specific solution in the technical solution of this application, the method of obtaining the concentration change network of each known point based on each known point includes:

[0014] Based on the known points, a first known point is obtained; the first known point is any one of the known points;

[0015] Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed with the first known point as the center;

[0016] Based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point;

[0017] Based on the rare earth element concentration change rates in various directions, a concentration change network of the first known point is obtained.

[0018] As a specific solution in the technical solution of the present application, the method of 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:

[0019] Based on each known point, a second known point is obtained; the second known point is any one of the known points;

[0020] Based on the second known point, a first reference point is obtained; the first reference point is any reference point among the reference points corresponding to the second known point;

[0021] Based on the second known point and the first reference point, a plurality of verification points are obtained; the verification point is any known point located between the second known point and the first reference point among the known points;

[0022] Based on each verification point, a deviation coefficient of the rare earth element concentration change of the second known point along a first direction is obtained; the first direction points from the second known point to the first reference point.

[0023] As a specific solution in the technical solution of the present application, the acquiring of multiple verification points based on the second known point and the first reference point includes:

[0024] Based on the second known point and the first reference point, a first line segment is obtained; 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;

[0025] Based on the first line segment, a first straight line and a second straight line are obtained; 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;

[0026] Each known point between the first straight line and the second straight line is used as a verification point.

[0027] As a specific solution in the technical solution of the present application, obtaining the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on each verification point includes:

[0028] Based on each verification point, a first verification point is obtained; the first verification point is any one of the verification points;

[0029] Based on the second known point and the first verification point, obtaining concentration change components of the second known point and the first verification point along a first direction;

[0030] Based on the concentration change component, a deviation coefficient of the rare earth element concentration change at the second known point along the first direction is obtained.

[0031] As a specific solution in the technical solution of this application, the calculation formula for obtaining the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point is as follows:

[0032] ;

[0033] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

[0034] As a specific solution in the technical solution of the present application, the calculation formula for obtaining the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component is as follows:

[0035] ;

[0036] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

[0037] As a specific solution in the technical solution of the present application, the interpolation weight of each known point is obtained based on the corresponding deviation coefficient in each direction of each known point and the point to be evaluated, including:

[0038] Based on each known point, a third known point is obtained; the third known point is any one of the known points;

[0039] Based on the deviation coefficients of the third known point in various directions, an effective deviation coefficient is obtained; the effective deviation coefficient is a deviation coefficient that is less than a preset value among the deviation coefficients;

[0040] Based on each effective deviation coefficient, obtaining a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient is used to at least 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;

[0041] Based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0042] In a second aspect, the present 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 comprising:

[0043] A processing module is used to obtain a partitioning grid based on a target mining area; the target mining area is obtained in advance;

[0044] A reading module is used 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;

[0045] The processing module is further configured to obtain a concentration variation network of each known point based on each known point;

[0046] And, based on the concentration change network, obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value;

[0047] and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point;

[0048] And, based on the interpolation weight of each known point, the rare earth element concentration of each unknown point is obtained.

[0049] As a specific solution in the technical solution of the present application, the processing module is further used to obtain a first known point based on each known point; the first known point is any one of the known points;

[0050] Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed with the first known point as the center;

[0051] And, based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point;

[0052] And, based on the rare earth element concentration change rate in each direction, a concentration change network of the first known point is obtained.

[0053] As a specific solution in the technical solution of the present application, the processing module is further used to obtain a second known point based on each known point; the second known point is any one of the known points;

[0054] And, based on the second known point, obtaining a first reference point; the first reference point is any reference point among the reference points corresponding to the second known point;

[0055] And, based on the second known point and the first reference point, a plurality of verification points are obtained; the verification point is any known point located between the second known point and the first reference point among the known points;

[0056] And, based on each verification point, a deviation coefficient of the rare earth element concentration change of the second known point along a first direction is obtained; the first direction points from the second known point to the first reference point.

[0057] As a specific solution in the technical solution of the present application, the processing module is further used 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; and the first reference point is the other end of the first line segment;

[0058] And, based on the first line segment, obtaining 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;

[0059] Furthermore, each known point between the first straight line and the second straight line is used as a verification point.

[0060] 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;

[0061] and, based on the second known point and the first verification point, obtaining concentration variation components of the second known point and the first verification point along a first direction;

[0062] And, based on the concentration change component, a deviation coefficient of the rare earth element concentration change at the second known point along the first direction is obtained.

[0063] As a specific solution in the technical solution of this application, the processing module obtains the calculation formula of the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point as follows:

[0064] ;

[0065] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

[0066] As a specific solution in the technical solution of the present application, the processing module obtains the calculation formula of the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component as follows:

[0067] ;

[0068] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

[0069] As a specific solution in the technical solution of the present application, the processing module is further used to obtain a third known point based on each known point; the third known point is any one of the known points;

[0070] And, based on the deviation coefficients of the third known point in each direction, obtaining an effective deviation coefficient; the effective deviation coefficient is a deviation coefficient less than a preset value among the deviation coefficients;

[0071] and obtaining, based on each effective deviation coefficient, a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient being used to at least characterize the degree of uniform change in the rare earth element concentration of the third known point along a direction from the third known point to the point to be evaluated;

[0072] And, based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0073] In a third aspect, the present application proposes a GIS-based rapid exploration system for rare earth minerals, which includes:

[0074] A processor is configured to obtain a partitioning grid based on a target mining area; the target mining area is obtained in advance;

[0075] A reader, configured to obtain a plurality of known points based on the divided grids; the known points are sampling points of any one of the divided grids, and the rare earth element concentration of each known point is known;

[0076] The processor is further configured to obtain a concentration variation network of each known point based on each known point;

[0077] And, based on the concentration change network, obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value;

[0078] and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point;

[0079] And, based on the interpolation weight of each known point, the rare earth element concentration of each unknown point is obtained.

[0080] 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;

[0081] Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed with the first known point as the center;

[0082] And, based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point;

[0083] And, based on the rare earth element concentration change rate in each direction, a concentration change network of the first known point is obtained.

[0084] 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;

[0085] And, based on the second known point, obtaining a first reference point; the first reference point is any reference point among the reference points corresponding to the second known point;

[0086] And, based on the second known point and the first reference point, a plurality of verification points are obtained; the verification point is any known point located between the second known point and the first reference point among the known points;

[0087] And, based on each verification point, a deviation coefficient of the rare earth element concentration change of the second known point along a first direction is obtained; the first direction points from the second known point to the first reference point.

[0088] 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; and the first reference point is the other end of the first line segment;

[0089] And, based on the first line segment, obtaining 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;

[0090] Furthermore, each known point between the first straight line and the second straight line is used as a verification point.

[0091] 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;

[0092] and, based on the second known point and the first verification point, obtaining concentration variation components of the second known point and the first verification point along a first direction;

[0093] And, based on the concentration change component, a deviation coefficient of the rare earth element concentration change at the second known point along the first direction is obtained.

[0094] As a specific solution in the technical solution of the present application, the processor obtains the calculation formula of the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point as follows:

[0095] ;

[0096] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

[0097] As a specific solution in the technical solution of the present application, the processor obtains the calculation formula of the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component as follows:

[0098] ;

[0099] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

[0100] 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;

[0101] And, based on the deviation coefficients of the third known point in each direction, obtaining an effective deviation coefficient; the effective deviation coefficient is a deviation coefficient less than a preset value among the deviation coefficients;

[0102] and obtaining, based on each effective deviation coefficient, a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient being used to at least characterize the degree of uniform change in the rare earth element concentration of the third known point along a direction from the third known point to the point to be evaluated;

[0103] And, based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0104] Compared with the prior art, the present invention has the following advantages:

[0105] This application calculates the rate of change of rare earth element concentration at each known point in each direction by taking several reference directions. The actual concentration change in each direction is then calculated using the concentrations of several known points in the corresponding directions to determine whether the actual concentration change is uniform. Interpolation weights are then obtained for each known point using the directions where the concentration change is relatively uniform. Compared to existing methods that use only the physical distance to the point to be evaluated to obtain interpolation weights, the interpolation weights obtained in this application are more accurate, which means that the rare earth element concentrations of the points to be evaluated obtained subsequently are also more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0106] Figure 1 A schematic diagram of a process for rapid exploration of rare earth minerals based on GIS proposed in an embodiment of the present application;

[0107] Figure 2 This is a schematic diagram of the structure of a GIS-based rare earth mineral rapid exploration electronic device proposed in an embodiment of the present application;

[0108] Figure 3 This is a schematic diagram of the structure of a GIS-based rare earth mineral rapid exploration system proposed in an embodiment of the present application;

[0109] Figure 4 This is a schematic diagram of setting sampling points using a serpentine point distribution method in dividing grid units proposed in an embodiment of the present application;

[0110] Figure 5 A schematic diagram of selecting a reference point according to an embodiment of the present application;

[0111] Figure 6 A schematic diagram of selecting verification points proposed in an embodiment of the present application;

[0112] Figure 7 This is a schematic diagram of the third known point and the point to be evaluated proposed in an embodiment of the present application. DETAILED DESCRIPTION

[0113] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0114] The terms "first", "second", etc. in the description of the embodiments of the present application and the above-mentioned drawings are used to distinguish similar objects and are not necessarily 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 the known points used in this way can be interchanged where appropriate so that the embodiments described herein can be implemented in an order other than that shown or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or modules is not necessarily limited to those steps or modules clearly listed, but may include other steps or modules that are not clearly listed or inherent to these processes, methods, products or devices. The division of modules that appears in the embodiments of the present application is only a logical division. In actual applications, there may be other division methods. For example, multiple modules can be combined or integrated into another system, or some features can be ignored or not performed. In addition, the coupling or direct coupling or communication connection between each other shown or discussed 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. Moreover, 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 into multiple circuit modules, and some or all of the modules may be selected according to actual needs to achieve the purpose of the embodiment of the present application.

[0115] In order to solve the technical problem that the prior art cannot accurately estimate the rare earth element concentration of the point to be evaluated, this application proposes a GIS-based rapid exploration method for rare earth minerals. Specifically, Figure 1 As shown, the GIS-based rapid exploration method for rare earth minerals includes steps S100 to S600.

[0116] Step S100: Based on the target mining area, a division grid is obtained.

[0117] In this embodiment, the target mining area is acquired in advance. That is, in this embodiment, the target mining area can be any area where rare earth element concentration estimation is required.

[0118] It should be noted that if a model needs to be built based on geochemical data, it is necessary to collect soil and rock samples from the surface. When sampling for mineral exploration, in order to make the sampling points distributed as evenly as possible in the mining area and avoid excessive concentration or omissions in some areas, the target mining area is usually divided into several equal-sized grids, and sampling points are selected within the grids. For example, the target mining area can be divided into several 50m×50m grids. The shape of the target mining area may not be completely divided into regular, equal-sized grids. When dividing, the number of regular grids should be as large as possible. For irregular areas, the sampling points 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 can be [ 】, 【】 denotes rounding of the values in brackets. Since the target mining area is usually a large area with uneven terrain and uneven soil, the serpentine method can be used as the grid point distribution method. Figure 4 As shown, the sampling points are set using the snake-shaped point distribution method in a single regular grid (i.e. Figure 4 Schematic diagram of the diamond hollow dots, diamond black dots and circular hollow dots).

[0119] Step S200: Based on the grid division, a plurality of known points are obtained.

[0120] In this embodiment, the known point is a sampling point of any grid in the divided grid (that is, Figure 4 The REE concentration at each known point is known. Samples are collected from all sampling points within the target mining area's grid. Each sampling point is marked as a known point. The REE content of all samples is measured to obtain the REE concentration at all known points. Testing for REE concentration in samples is a mature technology and will not be described in detail here.

[0121] It should be noted that because the distribution of rare earth elements in soil varies unevenly, interpolating directly based on the physical distances between several adjacent known points and the point to be evaluated approximates that the rare earth element concentration at each known point varies uniformly in all directions. Therefore, the resulting rare earth element concentration at the point to be evaluated is inaccurate. To perform interpolation more accurately, it is necessary to analyze the trend of rare earth element concentration changes at each known point along all directions (also referred to as the concentration change network below). This trend of rare earth element concentration changes along the line connecting each known point and the point to be evaluated can then be determined. This, combined with the distance between the known points and the point to be evaluated, allows for more accurate determination of interpolation weights and accurate determination of the rare earth element concentration at the point to be evaluated. The following examples of this application illustrate the present invention using only a single-unit grid. This does not imply that the embodiments proposed in this application are only applicable to obtaining a concentration change network for each known point based on a single-unit grid. It should be understood that a concentration change network for each known point can be obtained based on a single-unit grid or a multi-unit grid.

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

[0123] As previously mentioned, a concentration variation network refers to the trend of rare earth element concentration changes along various directions at corresponding known points. In embodiments of the present application, any suitable method can be used to obtain the concentration variation network for each known point. For example, in embodiments of the present application, step S300, based on each known point, obtaining the concentration variation network for each known point, may include steps S310 to S340.

[0124] Step S310: Based on various known points, obtain a first known point.

[0125] In this embodiment, the first known point is any one of the known points. That is, in this embodiment, each known point can use the same method as the first known point to obtain its corresponding concentration change network.

[0126] Step S320: Based on the first known point, multiple reference points are obtained from each known point.

[0127] It should be clear that in order 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 embodiment of the present application, the reference points are distributed with the first known point as the center. In this embodiment, a random selection method can be adopted 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 point distribution method when sampling, in order to make the direction of the concentration change analysis more suitable for the sampling situation, all the extreme points of the serpentine point distribution method (that is, as shown in FIG. Figure 4 The diamond hollow point shown in the figure is used as a reference point. The extreme point is the point closest to the boundary in the serpentine point distribution trajectory. Figure 4 As shown, it is assumed that the circular point is the first known point.

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

[0129] In this embodiment, each direction is from the first known point to the corresponding reference point. Figure 4 It can be seen that the number of reference points is 4, and the first known point pointing to any reference point can form a direction, that is, in this embodiment, the number of reference points is equal to the number of directions in the concentration change network. Figure 5 As shown, the four reference points can form four directions with the first known point. In other words, if there are four directions, the rare earth element concentration change rate of the first known point in the four directions needs to be obtained in step S340.

[0130] In an embodiment of the present application, any suitable method can be used to obtain the rare earth element concentration change rate in each direction based on the first known point and each reference point. For example, after obtaining the reference point, the concentration difference between the known point and each reference point can be compared to obtain the concentration change of the known point in multiple directions. Taking the i-th known point (that is, the first known point mentioned above) as an example, the concentration difference between the known point and each reference point is compared. The rare earth element concentration change rate from the i-th known point to the t-th reference point The calculation formula can be:

[0131] ;

[0132] in, Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration at the t-th reference point; represents the Euclidean distance between the i-th known point and the t-th reference point; It 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 a negative value.

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

[0134] It should be noted that, in this embodiment, the rare earth element concentration change rate of the first known point along all directions is the concentration change network of the first known point. Figure 5 As shown, the line connecting the i-th known point (also known as the first known point mentioned above) and each reference point can be recorded as each edge of the i-th known point. The value of each edge is the rate of change of rare earth element concentration from the i-th known point to the corresponding reference point. All edges of the i-th known point constitute the concentration change network of the i-th known point.

[0135] Step S400: 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.

[0136] It should be clear that when the existing technology only obtains interpolation weights based on distance for interpolation, the changes in rare earth element concentration in all directions of all known points are regarded as uniform changes. This application constructs a concentration change network, which splits the rare earth element concentration changes of each known point into multiple directions, so that the concentration changes in multiple directions can be discussed separately in subsequent analysis, and the rare earth element concentration changes of each known point are divided more finely, so that the subsequent measurement of interpolation weights is more accurate.

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

[0138] It should be noted that in this embodiment, any reasonable method can be used to obtain the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point based on the concentration change network. For example, step S400, 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, may include steps S410 to S440.

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

[0140] 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, in this embodiment, the deviation coefficient of the rare earth element concentration change in each direction can be obtained for each known point using the same method as the second known point.

[0141] Step S420: Acquire a first reference point based on the second known point.

[0142] In this embodiment, the first reference point is any one of the reference points corresponding to the second known point. As previously mentioned, the second known point and any one of the reference points can form a direction. That is, in this embodiment, the deviation coefficients of the rare earth element concentration change at the second known point along various directions can be obtained using the same method as the deviation coefficients of the rare earth element concentration change in the direction from the second known point to the first reference point.

[0143] Step S430: Acquire multiple verification points based on the second known point and the first reference point.

[0144] As previously mentioned, when calculating the rate of change of rare earth element concentration 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 assumed to be uniform. However, in reality, the concentration change may not be uniform. Therefore, it is necessary to analyze the reliability of the calculated rare earth element concentration rate of change result based on several other known points (i.e., verification points) between the two known points. In this embodiment, the verification points can be any known point located between the second known point and the first reference point.

[0145] In this embodiment, multiple verification points may be obtained from each known point in a random manner. In another embodiment of the present application, step S430, obtaining multiple verification points based on the second known point and the first reference point, may include steps S431 to S433.

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

[0147] 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 follows: Figure 6 As shown in the line segment 1, the circular hollow point in the line segment 1 is the second known point, and the diamond hollow point in the line segment 1 is the first reference point.

[0148] Step S432: Based on the first line segment, obtain a first straight line and a second straight line.

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

[0150] Step S433: using each known point between the first straight line and the second straight line as a verification point.

[0151] Specifically, such as Figure 6 As shown, a gray area is formed between the first straight line and the second straight line. That is, as long as the known point falls in the gray area (that is, Figure 6 The diamond-shaped black dots in the gray area shown can all be used as verification points between the second known point and the first reference point.

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

[0153] In this embodiment, the first direction is from the second known point to the first reference point. As mentioned above, the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration of the corresponding known point in the first direction and the uniform change value. In this embodiment, the deviation coefficient of the change in the rare earth element concentration of the second known point along the first direction can be obtained in any suitable manner based on each verification point. For example, step S440, based on each verification point, obtains the deviation coefficient of the change in the rare earth element concentration of the second known point along the first direction, which can include steps S441 to S443.

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

[0155] In this embodiment, the first verification point is any one of the verification points. As can be seen from step S442, each verification point needs to calculate the concentration change component along the first direction. In other words, in this embodiment, the concentration change component along the first direction of each verification point can refer to the first verification point, which will not be further described.

[0156] Step S442: 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 along the first direction.

[0157] In the embodiment of the present application, any reasonable method can be used to obtain the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point. For example, in step S442, the calculation formula for obtaining the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point can be as follows:

[0158] ;

[0159] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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. Figure 6 The angle formed by line segments 1 and 4 shown ); Represents the cosine function.

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

[0161] It's important to note that the closer the concentration change component corresponding to each verification point is to the rate of change of the rare earth element concentration at the second known point along the first direction, the more uniform the rare earth element concentration at the second known point along the first direction varies. In other words, when subsequently obtaining interpolation weights, the weight value in that direction should be increased. Conversely, when subsequently obtaining interpolation weights, the weight value in that direction should be decreased.

[0162] In this embodiment, any reasonable method can be used to obtain the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component. For example, in one embodiment of the present application, the variance or standard deviation of each concentration change component corresponding to the second known point and the rare earth element concentration change rate along the first direction can be used as the deviation coefficient. In another embodiment of the present application, step S443, based on the concentration change component, the calculation formula for obtaining the deviation coefficient of the rare earth element concentration change at the second known point along the first direction is as follows:

[0163] ;

[0164] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith 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; when the denominator is 0, it is replaced by 0.1; represents the sine function; represents the linear normalization function.

[0165] Step S500: Based on the deviation coefficient corresponding to each known point in each direction and the point to be evaluated, the interpolation weight of each known point is obtained.

[0166] In this embodiment, the point to be evaluated is pre-acquired. The point to be evaluated can be any point in the divided grid with an unknown rare earth element concentration. In the embodiment of the present application, any reasonable method can be used to obtain the interpolation weight of each known point based on the deviation coefficient corresponding to each known point in each direction and the point to be evaluated. For example, step S500, which obtains the interpolation weight of each known point based on the deviation coefficient corresponding to each known point in each direction and the point to be evaluated, can include steps S510 to S540.

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

[0168] In this embodiment, the third known point is any one of the known points. As previously mentioned, when calculating the rare earth element concentration at the point to be evaluated, it is necessary to assign an interpolation weight to each known point. In other words, in this embodiment, the interpolation weights for each known point can be obtained using the same method as for obtaining the interpolation weights for the third known point.

[0169] Step S520: Obtain an effective deviation coefficient based on the deviation coefficients of the third known point in various directions.

[0170] It should be noted that the rare earth element concentration variation along all directions of the third known point 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 rare earth element concentration variation of the third known point in that direction. If the inverse distance interpolation weight result obtained using the third known point in that direction has a large error, the interpolation result obtained for directions with smaller deviation estimates will be more accurate. Points with excessively large deviation estimates are screened using a preset threshold to avoid large deviations during interpolation. In an embodiment of the present application, a deviation coefficient with a smaller value among the deviation coefficients (i.e., an effective deviation coefficient) can be selected to participate in the subsequent interpolation weight acquisition. That is, in this embodiment, the effective deviation coefficient is a deviation coefficient among the deviation coefficients that is less than a preset value. In this embodiment, there is no restriction on the size of the preset value and it can be set according to needs. For example, the preset value can be 0.6 or 0.7, etc.

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

[0172] It should be clear that when using the inverse distance interpolation weight to obtain interpolation, the interpolation search radius must first be set. This application selects the nine-square grid range centered on the grid where the point to be evaluated is located as the search radius when obtaining the interpolation. According to the previous analysis, the smaller the deviation coefficient of a certain direction of the known point, the more uniform the rare earth element concentration of the known point along that direction, and the greater the interpolation weight of the known point along that direction should be. Using the deviation coefficient of each known point in each direction, combined with the Euclidean distance of each known point to the point to be evaluated, the interpolation weight of each known point to the point to be evaluated can be obtained.

[0173] In this embodiment, the uniformity adjustment coefficient is used to at least 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 embodiment of the present application, any reasonable method can be used to obtain the uniformity adjustment coefficients of 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 uniformity adjustment coefficients of the third known point and the point to be evaluated based on each effective deviation coefficient is as follows:

[0174] ;

[0175] in, The uniform adjustment coefficient of the i-th known point (also the third known point); Indicates the number of effective deviation coefficients of the i-th known point; Represents the angle between the point to be evaluated and the i-th known point in the m-th direction (i.e. Figure 7 Angle shown ); represents the deviation evaluation of the mth direction of the i-th known point, represents the inverse normalization function; As weights, calculate the weighted mean of the deviation evaluation, and use the inverse of the weighted mean as the uniform adjustment coefficient; The larger the value, the smaller the angle between the point to be evaluated and the mth effective side of the i-th known point. Here, +1 is to ensure The object is a positive number so that the final inverse normalization result is correctly mapped to the range of [0,1]. Figure 7 The points corresponding to the hollow triangles in are the points to be evaluated.

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

[0177] It should be noted that the uniform adjustment coefficient reflects the uniformity of the change of the i-th known point along the direction of the line connecting the points to be evaluated. The larger the value, the greater the reference degree of the i-th known point for interpolation, that is, the greater the interpolation weight. The uniform adjustment coefficient is used to adjust the interpolation weight method commonly used for distance calculation to obtain the interpolation weight of each known point. The interpolation weight of the i-th known point is The calculation method is as follows:

[0178] ;

[0179] Where: Represents the interpolation weight of the i-th known point for the point to be evaluated; represents the uniform adjustment coefficient of the i-th known point; Represents the Euclidean distance between the i-th known point and the point to be evaluated; Represents the exponential parameter, which affects the smoothness of the interpolation, here k=2; represents the linear normalization function.

[0180] Step S600: obtaining the rare earth element concentration of the point to be evaluated based on the interpolation weight of each known point.

[0181] In this embodiment, several points to be evaluated may be set, interpolation weights of all known points and each point to be evaluated may be calculated, and each point to be evaluated may be interpolated to obtain rare earth element concentrations of all points to be evaluated.

[0182] In this application, an embodiment of a GIS-based rapid rare earth mineral exploration method proposed by the present invention calculates the rate of change of rare earth element concentration in each direction at each known point by taking several reference directions. The concentrations of several known points in the corresponding directions are then used to calculate whether the actual concentration change in each direction is uniform. Interpolation weights are then obtained for each known point using the directions in which the concentration change is relatively uniform. Compared to existing methods that use only the physical distance from the point to be evaluated to obtain interpolation weights, the interpolation weights obtained in this application are more accurate, which in turn results in more accurate rare earth element concentrations at the point to be evaluated.

[0183] After introducing the GIS-based rapid exploration method for rare earth minerals proposed in the embodiment of this application, the following introduces an embodiment of the GIS-based rapid exploration electronic device for rare earth minerals proposed in this application, such as Figure 2 As shown, the GIS-based rare earth mineral rapid exploration electronic device 10 includes:

[0184] The processing module 11 is used to obtain a partitioning grid based on a target mining area; the target mining area is obtained in advance;

[0185] A reading module 12 is configured to obtain a plurality of known points based on the divided grids; the known points are sampling points of any one of the divided grids, and the rare earth element concentration of each known point is known;

[0186] The processing module 11 is further used to obtain a concentration change network of each known point based on each known point;

[0187] And, based on the concentration change network, obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value;

[0188] and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point;

[0189] And, based on the interpolation weight of each known point, the rare earth element concentration of each unknown point is obtained.

[0190] As a specific embodiment of the present 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;

[0191] Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed with the first known point as the center;

[0192] And, based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point;

[0193] And, based on the rare earth element concentration change rate in each direction, a concentration change network of the first known point is obtained.

[0194] As a specific embodiment of 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;

[0195] And, based on the second known point, obtaining a first reference point; the first reference point is any reference point among the reference points corresponding to the second known point;

[0196] And, based on the second known point and the first reference point, a plurality of verification points are obtained; the verification point is any known point located between the second known point and the first reference point among the known points;

[0197] And, based on each verification point, a deviation coefficient of the rare earth element concentration change of the second known point along a first direction is obtained; the first direction points from the second known point to the first reference point.

[0198] As a specific embodiment of 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; and the first reference point is the other end of the first line segment.

[0199] And, based on the first line segment, obtaining 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;

[0200] Furthermore, each known point between the first straight line and the second straight line is used as a verification point.

[0201] As a specific embodiment of 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;

[0202] and, based on the second known point and the first verification point, obtaining concentration variation components of the second known point and the first verification point along a first direction;

[0203] And, based on the concentration change component, a deviation coefficient of the rare earth element concentration change at the second known point along the first direction is obtained.

[0204] As a specific embodiment of the present application, the processing module 11 obtains the calculation formula of the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point as follows:

[0205] ;

[0206] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

[0207] As a specific embodiment of the present application, the processing module 11 obtains the calculation formula of the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component as follows:

[0208] ;

[0209] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

[0210] As a specific embodiment of 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;

[0211] And, based on the deviation coefficients of the third known point in each direction, obtaining an effective deviation coefficient; the effective deviation coefficient is a deviation coefficient less than a preset value among the deviation coefficients;

[0212] and obtaining, based on each effective deviation coefficient, a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient being used to at least characterize the degree of uniform change in the rare earth element concentration of the third known point along a direction from the third known point to the point to be evaluated;

[0213] And, based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0214] The embodiment of the GIS-based electronic equipment for rapid rare earth mineral exploration proposed in this application calculates the rate of change of rare earth element concentration in each direction at each known point by taking several reference directions. The concentrations of several known points in the corresponding directions are then used to calculate whether the actual concentration change in each direction is uniform. Interpolation weights are then obtained for each known point using the directions where the concentration change is relatively uniform. Compared to the prior art method of obtaining interpolation weights based solely on the physical distance to the point to be evaluated, the interpolation weights obtained in this application are more accurate, which means that the rare earth element concentrations of the points to be evaluated obtained subsequently are also more accurate.

[0215] After introducing the GIS-based rare earth mineral rapid exploration electronic equipment proposed in the embodiment of this application, the following introduces the embodiment of the GIS-based rare earth mineral rapid exploration system proposed in this application, such as Figure 3 As shown, the GIS-based rare earth mineral rapid exploration electronic system 20 includes:

[0216] The processor 21 is configured to obtain a partitioning grid based on a target mining area; the target mining area is obtained in advance;

[0217] The reader 22 is configured to obtain a plurality of known points based on the divided grids; the known points are sampling points of any one of the divided grids, and the rare earth element concentration of each known point is known;

[0218] The processor 21 is further configured to obtain a concentration variation network of each known point based on each known point;

[0219] And, based on the concentration change network, obtaining the deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value;

[0220] and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point;

[0221] And, based on the interpolation weight of each known point, the rare earth element concentration of each unknown point is obtained.

[0222] As a specific embodiment of the present 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;

[0223] Based on the first known point, a plurality of reference points are obtained from each known point; each reference point is distributed with the first known point as the center;

[0224] And, based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point;

[0225] And, based on the rare earth element concentration change rate in each direction, a concentration change network of the first known point is obtained.

[0226] As a specific embodiment of the present 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;

[0227] And, based on the second known point, obtaining a first reference point; the first reference point is any reference point among the reference points corresponding to the second known point;

[0228] And, based on the second known point and the first reference point, a plurality of verification points are obtained; the verification point is any known point located between the second known point and the first reference point among the known points;

[0229] And, based on each verification point, a deviation coefficient of the rare earth element concentration change of the second known point along a first direction is obtained; the first direction points from the second known point to the first reference point.

[0230] As a specific embodiment of the present 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; and the first reference point is the other end of the first line segment.

[0231] And, based on the first line segment, obtaining 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;

[0232] Furthermore, each known point between the first straight line and the second straight line is used as a verification point.

[0233] As a specific embodiment of the present 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;

[0234] and, based on the second known point and the first verification point, obtaining concentration variation components of the second known point and the first verification point along a first direction;

[0235] And, based on the concentration change component, a deviation coefficient of the rare earth element concentration change at the second known point along the first direction is obtained.

[0236] As a specific embodiment of the present application, the processor 21 obtains the calculation formula of the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point as follows:

[0237] ;

[0238] in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

[0239] As a specific embodiment of the present application, the processor 21 obtains the calculation formula of the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component as follows:

[0240] ;

[0241] in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

[0242] As a specific embodiment of 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;

[0243] And, based on the deviation coefficients of the third known point in each direction, obtaining an effective deviation coefficient; the effective deviation coefficient is a deviation coefficient less than a preset value among the deviation coefficients;

[0244] and obtaining, based on each effective deviation coefficient, a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient being used to at least characterize the degree of uniform change in the rare earth element concentration of the third known point along a direction from the third known point to the point to be evaluated;

[0245] And, based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

[0246] The embodiment of the GIS-based rapid rare earth mineral exploration system proposed in this application calculates the rate of change of rare earth element concentration in each direction at each known point by taking several reference directions. The concentrations of several known points in the corresponding directions are then used to calculate whether the actual concentration change in each direction is uniform. Interpolation weights are then obtained for each known point using the directions in which the concentration change is relatively uniform. Compared to the prior art, which uses only the physical distance to the point to be evaluated to obtain interpolation weights, the interpolation weights obtained in this application are more accurate, which means that the rare earth element concentrations of the points to be evaluated obtained later are also more accurate.

[0247] It should be understood that computer-readable storage media in this application include permanent and non-permanent, removable and non-removable media that can be used to store information 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 technology, read-only compact disc read-only memory, digital versatile disc or other optical storage, magnetic cassettes, disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include temporary computer-readable media such as modulated data signals and carrier waves.

[0248] In the above embodiments, the description of each embodiment has its own focus. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0249] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the methods, devices and equipment described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0250] In the several embodiments provided in 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 device embodiments described above are merely schematic. For example, the division of the modules is merely a logical function division. In actual implementation, there may be other division methods, such as 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 mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or modules, which can be electrical, mechanical or other forms.

[0251] The modules described as separate components may or may not be physically separate, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed across multiple network modules. Some or all of the modules may be selected to achieve the purpose of the present embodiment according to actual needs.

[0252] In addition, the functional modules in each embodiment of the present application can be integrated into a processing module, or each module can exist physically separately, or two or more modules can be integrated into a module. The above-mentioned integrated modules can be implemented in the form of hardware or in the form of software functional modules. If the integrated modules are implemented in the form of software functional modules and sold or used as independent products, they can be stored in a computer-readable storage medium.

[0253] In the above embodiments, all or part of the embodiments may be implemented by software, hardware, firmware, or any combination thereof. When implemented by software, all or part of the embodiments may be implemented in the form of a computer program product.

[0254] 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 accordance with the embodiments of the present application are fully or partially generated. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. 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 via a wired (e.g., coaxial cable, optical fiber, digital subscriber line) or wireless (e.g., infrared, wireless, microwave, etc.) method. The computer-readable storage medium may be any available medium that can be stored on a computer or a data storage device such as a server or data center that integrates one or more available media. The available medium may be a magnetic medium (e.g., a floppy disk, hard disk, or magnetic tape), an optical medium (e.g., a digital versatile disk), or a semiconductor medium (e.g., a solid-state drive (SSD)).

[0255] Although the embodiments of the present application have been shown and described, it will be appreciated by those skilled in the art that many changes, modifications, substitutions, and alterations may 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 by: include: Based on the target mining area, obtain the divided grid; The target mining area is obtained in advance; Based on the grid division, a plurality of known points are obtained; The known points are sampling points of any one of the divided grids, 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, a deviation coefficient of the rare earth element concentration change in each direction corresponding to the known point is obtained; the deviation coefficient is used to at least characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction corresponding to the known point and the uniform change value; Based on the deviation coefficients corresponding to the respective directions of the respective known points and the points to be evaluated, the interpolation weights of the respective known points are obtained in advance; Obtaining the rare earth element concentration of the point to be evaluated based on the interpolation weight of each known point; The step of 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 reference point among 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 point is any known point located between the second known point and the first reference point among the known points; Based on each verification point, obtaining a deviation coefficient of a change in rare earth element concentration at the second known point along a first direction, wherein the first direction is from the second known point to the first reference point; The acquiring 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, a first line segment is obtained; 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, a first straight line and a second straight line are obtained; 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; using each known point between the first straight line and the second straight line as each verification point; The step of obtaining the interpolation weight of each known point based on the corresponding deviation coefficient in each direction of each known point and the point to be evaluated includes: Based on each known point, a third known point is obtained; the third known point is any one of the known points; Based on the deviation coefficients of the third known point in various directions, an effective deviation coefficient is obtained; the effective deviation coefficient is a deviation coefficient that is less than a preset value among the deviation coefficients; Based on each effective deviation coefficient, obtaining a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient is used to at least 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; Based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

2. The GIS-based rapid exploration method for rare earth minerals according to claim 1, characterized in that: The step of obtaining a concentration change network of each known point based on each known point includes: Based on the known points, 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 with the first known point as the center; Based on the first known point and each reference point, obtaining the rare earth element concentration change rate in each direction; each direction is pointed from the first known point to the corresponding reference point; Based on the rare earth element concentration change rates in various directions, a concentration change network of the first known point is obtained.

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

4. The GIS-based rapid exploration method for rare earth minerals according to claim 3, characterized in that: The calculation formula for obtaining the concentration change component of the second known point and the first verification point along the first direction based on the second known point and the first verification point is as follows: ； in, The concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; represents the rare earth element concentration at the i-th known point; represents the rare earth element concentration of the yth verification point of the tth reference point of the i-th known point; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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 cosine function.

5. The GIS-based rapid exploration method for rare earth minerals according to claim 4, characterized in that: The calculation formula for obtaining the deviation coefficient of the rare earth element concentration change at the second known point along the first direction based on the concentration change component is as follows: ； in, It represents the deviation coefficient of the rare earth element concentration change from the i-th known point to the t-th reference point; represents the variance of the concentration change component 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 concentration change rate component of the yth verification point of the tth reference point of the ith known point along the direction from the ith known point to the tth reference point; Indicates the concentration change rate of the ith known point in the concentration change network pointing to the tth reference point; Indicates taking the absolute value; represents the Euclidean distance between the ith known point and the yth verification point of the tth reference point in the ith known point; 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.

6. GIS-based rare earth mineral rapid exploration electronic equipment, characterized by: include: A processing module is used to obtain a divided grid based on the target mining area; The target mining area is obtained in advance; A reading module is used 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 variation network of each known point based on each known point; and, based on the concentration variation network, obtaining deviation coefficients of rare earth element concentration variations in various directions corresponding to known points; The deviation coefficient is used to characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction of the corresponding known point and the uniform change value; and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point; and, obtaining the rare earth element concentration of each unknown point based on the interpolation weight of each known point; The step of 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 reference point among 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 point is any known point located between the second known point and the first reference point among the known points; Based on each verification point, obtaining a deviation coefficient of a change in rare earth element concentration at the second known point along a first direction, wherein the first direction is from the second known point to the first reference point; The acquiring 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, a first line segment is obtained; 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, a first straight line and a second straight line are obtained; 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; using each known point between the first straight line and the second straight line as each verification point; The step of obtaining the interpolation weight of each known point based on the corresponding deviation coefficient in each direction of each known point and the point to be evaluated includes: Based on each known point, a third known point is obtained; the third known point is any one of the known points; Based on the deviation coefficients of the third known point in various directions, an effective deviation coefficient is obtained; the effective deviation coefficient is a deviation coefficient that is less than a preset value among the deviation coefficients; Based on each effective deviation coefficient, obtaining a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient is used to at least 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; Based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

7. The GIS-based rare earth mineral rapid exploration system is characterized by: include: A processor, configured to obtain a partitioning 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 grids; the known points are sampling points of any one of the divided grids, and the rare earth element concentration of each known point is known; The processor is further configured to obtain a concentration variation network of each known point based on each known point; and, based on the concentration variation network, obtaining deviation coefficients of rare earth element concentration variations in various directions corresponding to known points; The deviation coefficient is used to characterize the degree of deviation between the actual change value of the rare earth element concentration in each direction of the corresponding known point and the uniform change value; and, obtaining an interpolation weight of the known point based on the corresponding deviation coefficient in each direction of the known point; and, obtaining the rare earth element concentration of each unknown point based on the interpolation weight of each known point; The step of 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 reference point among 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 point is any known point located between the second known point and the first reference point among the known points; Based on each verification point, obtaining a deviation coefficient of a change in rare earth element concentration at the second known point along a first direction, wherein the first direction is from the second known point to the first reference point; The acquiring 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, a first line segment is obtained; 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, a first straight line and a second straight line are obtained; 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; using each known point between the first straight line and the second straight line as each verification point; The step of obtaining the interpolation weight of each known point based on the corresponding deviation coefficient in each direction of each known point and the point to be evaluated includes: Based on each known point, a third known point is obtained; the third known point is any one of the known points; Based on the deviation coefficients of the third known point in various directions, an effective deviation coefficient is obtained; the effective deviation coefficient is a deviation coefficient that is less than a preset value among the deviation coefficients; Based on each effective deviation coefficient, obtaining a uniform adjustment coefficient for the third known point and the point to be evaluated; the uniform adjustment coefficient is used to at least 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; Based on the uniform adjustment coefficient, an interpolation weight of the third known point is obtained.

Images