A method, device, equipment, medium and product for digital reconstruction of gravity and magnetic anomaly color images

Abstract

The present application discloses a method, device, equipment, medium, and product for digitally reconstructing a gravity and magnetic anomaly color image, relating to the field of geophysical exploration technology. The method comprises: obtaining a gravity and magnetic anomaly color image; determining an actual geographic coordinate system based on projection parameters; converting the gravity and magnetic anomaly color image into an actual geographic coordinate system using pixel points to obtain anomaly grid data; determining discrete sample points based on the mapping scale of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data; editing the sample grid data based on the anomaly grid data to obtain initial contour lines of the gravity and magnetic anomaly color image; and reconstructing the anomaly grid data based on the actual anomaly values ​​and initial contour lines of the gravity and magnetic anomaly color image to complete the reconstruction of the gravity and magnetic anomaly color image. The present application can reconstruct a gravity and magnetic anomaly color image, and the reconstructed gravity and magnetic anomaly color image can directly extract gravity and magnetic anomaly data.

Description

Technical Field

[0001] The present application relates to the field of geophysical exploration technology, and in particular to a method, device, equipment, medium and product for digitally reconstructing color images of gravity and magnetic anomalies. Background Art

[0002] In gravity and magnetic exploration, the primary form of exploration results is gravity and magnetic data. Comprehensive analysis and research of these data require their visualization. Maps drawn from gravity and magnetic data are called gravity and magnetic anomaly (gravity and magnetic anomaly) maps. In practical work, cross-sectional and planar gravity and magnetic anomaly maps are the most commonly used. Planar contour maps of gravity and magnetic anomalies are crucial foundational maps for further geological research and analysis. They assess the changing state of gravity and magnetic anomalies through the sparseness of contour lines, their strike characteristics, and their amplitude. Furthermore, because planar color images facilitate the comparison and analysis of anomaly features, combining contour maps with color images, which can be used to create different presentation effects by varying color depth, brightness, and hue, is increasingly popular among geophysical researchers.

[0003] However, in related technologies, gravity and magnetic anomaly data cannot be directly extracted from planar color images, which causes many inconveniences in the comprehensive analysis and research of gravity and magnetic anomaly data. Summary of the Invention

[0004] The purpose of this application is to provide a method, device, equipment, medium and product for digital reconstruction of gravity and magnetic anomaly color images, which can reconstruct gravity and magnetic anomaly color images, and the reconstructed gravity and magnetic anomaly color images can directly extract gravity and magnetic anomaly data.

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

[0006] In a first aspect, the present application provides a method for digitally reconstructing a color image of gravity and magnetic anomalies, comprising:

[0007] Obtain color images of gravity and magnetic anomalies;

[0008] Determine the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection datum plane, and a projection plane unit;

[0009] Converting the gravity and magnetic anomaly color image into the actual geographic coordinate system using pixel points to obtain anomaly grid data;

[0010] Determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data;

[0011] Editing the sample point grid data according to the abnormal grid data to obtain initial contour lines of the gravity and magnetic anomaly color image;

[0012] The abnormal grid data is reconstructed according to the actual abnormal value of the gravity and magnetic anomaly color image and the initial isovalue line to complete the reconstruction of the gravity and magnetic anomaly color image.

[0013] Furthermore, the abnormal grid data includes an initial abnormal value of the abnormal point and an actual coordinate of the abnormal point in the actual geographic coordinate system;

[0014] The method of converting the gravity and magnetic anomaly color image into the actual geographic coordinate system using pixel points to obtain anomaly grid data includes:

[0015] Selecting a plurality of pixel points of the gravity and magnetic anomaly color image as image control points;

[0016] Determining the actual coordinates of the image control point according to the pixel coordinates of the image control point, and determining the conversion relationship between the pixel coordinates of the image control point and the actual coordinates;

[0017] According to the transformation relationship, all pixel points in the gravity and magnetic anomaly color image are transformed into an actual geographic coordinate system to obtain the actual coordinates of the anomaly points; wherein the anomaly points correspond to the pixel points one to one;

[0018] Calculate the initial outlier value of the outlier point according to the RGB value of the pixel point;

[0019] The abnormal grid data is obtained based on the actual coordinates of the abnormal point and the initial abnormal value of the abnormal point.

[0020] Furthermore, the calculation formula of the initial outlier value of the outlier point is:

[0021] A N =0.299*R+0.587*G+0.114*B;

[0022] Among them, A N is the initial outlier value of the outlier point, R is the component value of the red channel, G is the component value of the green channel, and B is the component value of the blue channel.

[0023] Furthermore, the sample point grid data includes actual coordinates of the discrete sample points and initial outlier values ​​of the discrete sample points;

[0024] The method of determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image and obtaining sample grid data of the discrete sample points by using the discrete sample points and the anomaly grid data comprises:

[0025] Determining discrete sample points according to the mapping scale, preset point spacing, and preset line spacing of the gravity and magnetic anomaly color image;

[0026] Determine the actual coordinates of the discrete sample points according to the position of the gravity and magnetic anomaly color image in the actual geographic coordinate system;

[0027] Based on the abnormal grid data, the radial basis function method is used to interpolate the GRD data of the abnormal point;

[0028] Extracting initial outlier values ​​of the discrete sample points from the GRD data of the outlier points using the actual coordinates of the discrete sample points;

[0029] The sample point grid data of the discrete sample points are obtained based on the actual coordinates of the discrete sample points and the initial abnormal values ​​of the discrete sample points.

[0030] Furthermore, the sample point grid data is edited according to the abnormal grid data to obtain the initial contour lines of the gravity and magnetic anomaly color image, including:

[0031] Draw the contour lines of the outlier points according to the GRD data of the outlier points;

[0032] Filtering the initial outliers of the discrete sample points;

[0033] Based on the actual coordinates of the discrete sample points and the initial outlier values ​​of the filtered discrete sample points, the improved Shepard method is used to interpolate and obtain the GRD data of the discrete sample points;

[0034] Draw the contour lines of the discrete sample points according to the GRD data of the discrete sample points;

[0035] The contour lines of the discrete sample points are edited according to the contour lines of the abnormal points to obtain the initial contour lines of the gravity and magnetic anomaly color image.

[0036] Furthermore, the abnormal grid data is reconstructed according to the actual abnormal value of the gravity and magnetic abnormal color image and the initial contour line, including:

[0037] Re-assigning the initial contour lines according to the actual abnormal values ​​of the gravity and magnetic anomaly color image to obtain contour lines of the actual abnormal values;

[0038] The abnormal grid data is reconstructed based on the contour lines of the abnormal true value.

[0039] In a second aspect, the present application provides a device for digitally reconstructing color images of gravity and magnetic anomalies, comprising:

[0040] An acquisition module is used to obtain color images of gravity and magnetic anomalies;

[0041] A coordinate system determination module is used to determine the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection reference plane, and a projection plane unit;

[0042] A conversion module is used to convert the gravity and magnetic anomaly color image into an actual geographic coordinate system using pixel points to obtain anomaly grid data;

[0043] a discrete module for determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data;

[0044] An editing module, configured to edit the sample point grid data according to the abnormal grid data to obtain initial contour lines of the gravity and magnetic anomaly color image;

[0045] A reconstruction module is used to reconstruct the abnormal grid data according to the actual abnormal value of the gravity and magnetic abnormal color image and the initial contour line.

[0046] In a third aspect, the present application provides a computer device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-mentioned method for digitally reconstructing color images of gravity and magnetic anomalies.

[0047] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned method for digitally reconstructing color images of gravity and magnetic anomalies.

[0048] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements the above-mentioned method for digitally reconstructing color images of gravity and magnetic anomalies.

[0049] According to the specific embodiments provided in this application, this application discloses the following technical effects:

[0050] The present application provides a method, apparatus, device, medium and product for digital reconstruction of gravity and magnetic anomaly color images, which obtain gravity and magnetic anomaly color images and determine their actual geographic coordinate system based on projection parameters; convert the gravity and magnetic anomaly color images into the actual geographic coordinate system, that is, convert the gravity and magnetic anomaly color images into a specific actual geographic coordinate system according to the research needs of gravity and magnetic anomaly data, and then obtain abnormal grid data; again, determine discrete sample points according to the mapping scale, discretize the gravity and magnetic anomaly color images, and obtain sample grid data using the abnormal grid data and discrete sample points; finally, edit the sample grid data according to the abnormal grid data to obtain initial contour lines of the gravity and magnetic anomaly color images; on this basis, reconstruct the abnormal grid data using the actual abnormal values ​​and the initial contour lines, so that the reconstructed abnormal grid data has abnormal true values, and complete the digital reconstruction of the gravity and magnetic anomaly color images. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0052] Figure 1 This is a diagram of the application environment of a method for digitally reconstructing color images of gravity and magnetic anomalies in one embodiment of the present application;

[0053] Figure 2 A schematic flow chart of a method for digitally reconstructing a color image of gravity and magnetic anomalies provided in one embodiment of the present application;

[0054] Figure 3 are the original image of the gravity and magnetic anomaly color image and the contour image of the gravity and magnetic anomaly color image; among them, (a) is the original image of the gravity and magnetic anomaly color image, (b) is the contour image of the gravity and magnetic anomaly color image;

[0055] Figure 4 is the initial contour line combination map of the edited gravity and magnetic anomaly color image; where (a) is the edited gravity and magnetic anomaly color image combination map, and (b) is the edited contour line image combination map of the gravity and magnetic anomaly color image;

[0056] Figure 5 A schematic diagram of a reconstructed gravity and magnetic anomaly map provided in one embodiment of the present application;

[0057] Figure 6 A schematic diagram of the functional modules of a device for digitally reconstructing color images of gravity and magnetic anomalies provided in one embodiment of the present application;

[0058] Figure 7 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0059] 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.

[0060] In order to make the purpose, features and advantages of this application more obvious and easy to understand, this application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0061] The method for digitally reconstructing color images of gravity and magnetic anomalies provided in the embodiment of the present application can be applied to Figure 1 In the application environment shown, the terminal 102 communicates with the server 104 via a network. The data storage system can store data that the server 104 needs to process. The data storage system can be set up separately, integrated on the server 104, or placed on the cloud or other servers. Terminal 102 can send the gravity and magnetic anomaly color image to server 104. After receiving the gravity and magnetic anomaly color image, server 104 determines the actual geographic coordinate system of the gravity and magnetic anomaly color image based on projection parameters, wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection datum, and a projection plane unit. The gravity and magnetic anomaly color image is converted into the actual geographic coordinate system using pixels to obtain anomaly grid data. Discrete sample points are determined based on the mapping scale of the gravity and magnetic anomaly color image, and sample grid data of the discrete sample points is obtained using the discrete sample points and the anomaly grid data. The sample grid data is edited based on the anomaly grid data to obtain initial contour lines of the gravity and magnetic anomaly color image. The anomaly grid data is reconstructed based on the actual anomaly values ​​of the gravity and magnetic anomaly color image and the initial contour lines, thereby completing the reconstruction of the gravity and magnetic anomaly color image. Server 104 can provide feedback of the reconstructed gravity and magnetic anomaly color image to terminal 102. In addition, in some embodiments, the method for digitally reconstructing the gravity and magnetic anomaly color image can also be implemented independently by the server 104 or the terminal 102. For example, the terminal 102 can directly perform digital reconstruction processing on the gravity and magnetic anomaly color image, or the server 104 can obtain the gravity and magnetic anomaly color image from the data storage system and perform digital reconstruction processing on the gravity and magnetic anomaly color image.

[0062] The terminal 102 may be, but is not limited to, various desktop computers, laptop computers, smart phones, tablet computers, etc. The server 104 may be implemented as an independent server or a server cluster consisting of multiple servers, or may be a cloud server.

[0063] In an exemplary embodiment, Figure 2 As shown, a method for digital reconstruction of color images of gravity and magnetic anomalies is provided. The method is executed by a computer device, specifically a computer device such as a terminal or a server, and can also be executed by a terminal and a server together. In the embodiment of the present application, the method is applied to Figure 1 The server 104 in the example is used as an example to illustrate the process, including the following steps 201 to 206.

[0064] Step 201: Acquire a color image of gravity and magnetic anomalies.

[0065] Specifically, according to the actual research work situation, the gravity and magnetic anomaly color image whose gravity and magnetic anomaly data needs to be extracted is found from the literature or reports, and the gravity and magnetic anomaly color image is stored in JPG format.

[0066] Step 202: determining the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection datum plane, and a projection plane unit.

[0067] Step 203: Using pixel points, the gravity and magnetic anomaly color image is converted into the actual geographic coordinate system to obtain anomaly grid data.

[0068] Step 204 : determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data.

[0069] Step 205 : Editing the sample point grid data according to the abnormal grid data to obtain the initial contour lines of the gravity and magnetic anomaly color image.

[0070] Step 206 : reconstructing the abnormal grid data according to the actual abnormal value of the gravity and magnetic anomaly color image and the initial contour line to complete the reconstruction of the gravity and magnetic anomaly color image.

[0071] By executing steps 201 to 206 , the gravity and magnetic anomaly color image can be digitally reconstructed, and the gravity and magnetic anomaly data can be extracted using the reconstructed gravity and magnetic anomaly color image.

[0072] In an exemplary embodiment, the abnormal grid data includes the initial abnormal value of the abnormal point and the actual coordinates of the abnormal point in the actual geographic coordinate system, and step 203 specifically includes steps 20301 to 20305:

[0073] Step 20301: Select multiple pixel points of the gravity and magnetic anomaly color image as image control points.

[0074] Step 20302: Determine the actual coordinates of the image control point based on the pixel coordinates of the image control point, and determine the conversion relationship between the pixel coordinates and the actual coordinates of the image control point.

[0075] Step 20303: transform all pixel points in the gravity and magnetic anomaly color image into the actual geographic coordinate system according to the transformation relationship to obtain the actual coordinates of the anomaly points; wherein the anomaly points and the pixel points correspond one to one.

[0076] Step 20304: Calculate the initial outlier value of the outlier point based on the RGB value of the pixel point.

[0077] The calculation formula of the initial outlier value of the outlier point is:

[0078] A N =0.299*R+0.587*G+0.114*B;

[0079] Among them, A N is the initial outlier value of the outlier point, R is the component value of the red channel, G is the component value of the green channel, and B is the component value of the blue channel.

[0080] Step 20305: Obtain the abnormal grid data based on the actual coordinates of the abnormal point and the initial abnormal value of the abnormal point.

[0081] In the specific implementation process, the Global Mapper software is used to convert the gravity and magnetic anomaly color image into actual geographic coordinates. Specifically, the Global Mapper software is used to open the JPG image of the gravity and magnetic anomaly color image obtained in step 201, and the projection parameters are set in the Global Mapper software according to the research needs of the gravity and magnetic anomaly data, so as to determine the actual geographic coordinates of the gravity and magnetic anomaly color image. On this basis, multiple pixel points are selected as image control points. The number of image control points is generally greater than 4 to achieve the purpose of accurate conversion of image coordinates. The specific number is not limited in this embodiment.

[0082] Input the pixel coordinates of the selected pixel point and the actual coordinates of the corresponding image control point into the GlobalMapper software to determine the conversion relationship between the pixel point and the image control point. The pixel point converted to the actual geographic coordinates is called an outlier.

[0083] According to this transformation relationship, all pixel points in the gravity and magnetic anomaly color image are transformed into actual geographic coordinates. Each pixel point corresponds to an anomaly point in the actual geographic coordinates, thereby obtaining the actual coordinates of the anomaly points corresponding to all pixel points, and obtaining a JPG image file with actual coordinates. Combined with the initial anomaly value of the anomaly point, the anomaly grid data is obtained.

[0084] In an exemplary embodiment, the sample point grid data includes the actual coordinates of the discrete sample points and the initial outlier values ​​of the discrete sample points. Then, step 204 specifically includes steps 20401 to 20405:

[0085] Step 20401: Determine discrete sample points according to the mapping scale, preset point spacing, and preset line spacing of the gravity and magnetic anomaly color image.

[0086] Specifically, the mapping scale of the gravity and magnetic anomaly color image can be obtained from the literature or report on obtaining the gravity and magnetic anomaly color image. The mapping scale is used to characterize the data accuracy of the gravity and magnetic anomaly color image. For example, the mapping scale may be 1:50,000 or 1:1,000,000. Then the corresponding preset point spacing and preset line spacing can be set to 500m×500m preset point spacing×preset line spacing, or 10,000m×10,000m preset point spacing×preset line spacing; different preset point spacings and preset line spacings can be used to determine the spacing of discrete sample points for different mapping scales. Generally, when the mapping scale is 1:50,000, the spacing of discrete sample points is determined to be 0.5km×0.5km. Other spacings of discrete sample points can also be used, which is not limited in this embodiment.

[0087] Step 20402: Determine the actual coordinates of discrete sample points based on the positions of the gravity and magnetic anomaly color image in the actual geographic coordinate system.

[0088] A uniform grid map is formed according to the preset point spacing and the preset line spacing. Each point in the grid map is regarded as a discrete sample point, and a plurality of discrete sample points are finally determined.

[0089] The grid map where the discrete sample points are located is overlapped with the gravity and magnetic anomaly color image in the actual geographic coordinate system, and the corresponding actual coordinates of the discrete sample points in the gravity and magnetic anomaly color image area are used as the actual coordinates of the discrete sample points, thereby determining the actual coordinates of the discrete sample points in the gravity and magnetic anomaly color image.

[0090] Step 20403: Based on the abnormal grid data, radial basis function method is used to interpolate to obtain GRD data of the abnormal point.

[0091] Specifically, based on the initial outlier value of the outlier point, the GlobalMapper software was used to convert the JPG image file with the actual coordinates and output it as a data file. The radial basis function method in the Surfer software was then used to interpolate the data file into a GRD file to obtain the GRD data of the outlier point.

[0092] Step 20404: extract the initial outlier value of the discrete sample point from the GRD data of the outlier point using the actual coordinates of the discrete sample point.

[0093] The initial outliers in GRD data change continuously, and the initial outliers corresponding to the actual coordinates of the discrete sample points are used as the initial outliers of the discrete sample points.

[0094] Step 20405: Obtain sample point grid data of the discrete sample points based on the actual coordinates of the discrete sample points and the initial abnormal values ​​of the discrete sample points.

[0095] In an exemplary embodiment, step 205 specifically includes steps 20501 to 20505:

[0096] Step 20501: Draw the contour line of the outlier point based on the GRD data of the outlier point.

[0097] Specifically, the contour line drawing interface of the Surfer software is used to draw the contour lines of the abnormal points.

[0098] Step 20502: Filter the initial outliers of the discrete sample points.

[0099] Specifically, a symmetrical small subdomain filtering method (gradient-preserving filtering method) is used to filter the initial anomalies of discrete sample points, which can better preserve the characteristics of the original anomaly field while filtering and denoising.

[0100] Step 20503: Based on the actual coordinates of the discrete sample points and the initial outlier values ​​of the filtered discrete sample points, an improved Shepard method is used to perform interpolation to obtain the GRD data of the discrete sample points.

[0101] Step 20504: Draw the contour lines of the discrete sample points based on the GRD data of the discrete sample points.

[0102] Specifically, the same method as step 20501 is used to draw the contour lines of the discrete sample points.

[0103] Step 20505: Edit the contour lines of the discrete sample points according to the contour lines of the abnormal points to obtain the initial contour lines of the gravity and magnetic anomaly color image.

[0104] The contour lines of the discrete sample points do not completely coincide with the contour lines of the outliers. Therefore, the contour lines of the discrete sample points are edited according to the contour lines of the outliers so that the contour lines of the discrete sample points coincide with the contour lines of the outliers one by one, which can correspond well.

[0105] Specifically, there are two types of images provided in literature or reports: one is when the original literature contains not only color images of gravity and magnetic anomalies but also contour images of these images, i.e., contour images of the anomaly points; the other is when the original literature contains only color images of gravity and magnetic anomalies. The contour editing methods for both cases are largely the same; however, editing and adjusting the contours is faster and easier when the contour images of the anomaly points are included.

[0106] like Figure 3 As shown, Figure 3 (a) is the original color image of gravity and magnetic anomaly provided in the literature or reports. Figure 3 (b) in the figure is the contour image of the color image of gravity and magnetic anomalies provided in the literature or reports.

[0107] like Figure 4 As shown, Figure 4 The edited gravity and magnetic anomaly color image composite image in (a) includes a color background image and a red curve, where the color background image is the gravity and magnetic anomaly color image with actual coordinates, and the red curve is the edited initial contour line.

[0108] Figure 4 The contour image combination of the edited gravity and magnetic anomaly color image in (b) includes a grayscale background image, a black curve and a red curve, wherein the image composed of the grayscale background image and the black curve is the contour image of the gravity and magnetic anomaly color image with actual coordinates, and the red curve is the initial contour line after editing.

[0109] In an exemplary embodiment, step 206 specifically includes steps 20601 and 20602:

[0110] Step 20601: re-assign the initial contour lines according to the actual abnormal values ​​of the gravity and magnetic anomaly color image to obtain the contour lines of the actual abnormal values.

[0111] Specifically, in literature or reports, the distribution of the true abnormal values ​​of the gravity and magnetic anomaly color image is generally marked by the image legend color bar. Therefore, by reassigning the initial contour lines according to the actual abnormal values ​​marked by the legend color bar of the gravity and magnetic anomaly color image, the contour lines of the true abnormal values ​​of the gravity and magnetic anomaly color image can be obtained.

[0112] Step 20602: reconstruct the abnormal grid data based on the contour lines of the abnormal true value.

[0113] Specifically, using Surfer software, the contour lines of the true anomaly values ​​are output one by one, obtaining multiple point data. Each point data point includes the horizontal coordinate, vertical coordinate, and true anomaly value in actual physical coordinates. Finally, all the output point data are continuously integrated into a single anomaly result file. In actual operation, before outputting the point data, the contour line assignment results of the true anomaly values ​​are checked and verified based on the characteristics of the contour lines to ensure the correctness and accuracy of the assignment results of each contour line. Among them, the characteristics of contour lines include that the true anomaly values ​​of adjacent contour lines generally change gradually, the true anomaly values ​​of contour lines corresponding to the same color boundary are consistent, and the true anomaly value of a closed color area gradually increases or decreases from the boundary to the center.

[0114] Finally, the radial basis function method of Surfer software is used to interpolate the abnormal result file to form a new GRD file, and the digital reconstruction of the gravity and magnetic anomaly color image is completed, such as Figure 5 As shown in FIG, the gravity and magnetic anomaly map is obtained after completing the digital reconstruction of the gravity and magnetic anomaly color image. Based on the gravity and magnetic anomaly map, the required gravity and magnetic anomaly data can be extracted.

[0115] Based on the same inventive concept, embodiments of the present application also provide a device for digitally reconstructing color images of gravity and magnetic anomalies, for implementing the aforementioned method for digitally reconstructing color images of gravity and magnetic anomalies. The solution provided by this device is similar to the solution described in the aforementioned method. Therefore, the specific limitations of one or more embodiments of the device for digitally reconstructing color images of gravity and magnetic anomalies provided below can be found in the aforementioned limitations of the method for digitally reconstructing color images of gravity and magnetic anomalies, and will not be further elaborated here.

[0116] In an exemplary embodiment, Figure 6 As shown, a device for digitally reconstructing color images of gravity and magnetic anomalies is provided, comprising:

[0117] The acquisition module 61 is used to acquire a color image of gravity and magnetic anomalies.

[0118] The coordinate system determination module 62 is used to determine the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include projection ellipsoid, projection zone number, projection datum plane and projection plane unit.

[0119] The conversion module 63 is used to convert the gravity and magnetic anomaly color image into an actual geographic coordinate system using pixel points to obtain anomaly grid data.

[0120] The discrete module 64 is configured to determine discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtain sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data.

[0121] The editing module 65 is configured to edit the sample point grid data according to the abnormal grid data to obtain the initial contour lines of the gravity and magnetic anomaly color image.

[0122] The reconstruction module 66 is used to reconstruct the abnormal grid data according to the actual abnormal value of the gravity and magnetic abnormal color image and the initial contour line.

[0123] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure of the computer device is as follows: Figure 7 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O) and a communication interface. The processor, memory and input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store color images of gravity and magnetic anomalies. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a method for digitally reconstructing color images of gravity and magnetic anomalies is implemented.

[0124] Those skilled in the art will understand that Figure 7 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

[0125] In an exemplary embodiment, a computer device is further provided, including a memory and a processor. The memory stores a computer program, and the processor implements the steps in the above method embodiments when executing the computer program.

[0126] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented.

[0127] In an exemplary embodiment, a computer program product is provided, including a computer program. When the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0128] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0129] Those skilled in the art will understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, database or other media used in the embodiments provided in this application may include at least one of non-volatile and volatile memory. Non-volatile memory may include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory may include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0130] The databases involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processors involved in the various embodiments provided herein may include, but are not limited to, general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic units, data processing logic units based on quantum computing, and the like.

[0131] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0132] This specification uses specific examples to illustrate the principles and implementation methods of this application. The above examples are only intended to help understand the method and core concept of this application. At the same time, those skilled in the art will find that the specific implementation methods and application scopes may vary based on the concept of this application. In summary, the contents of this specification should not be construed as limiting this application.

Claims

1. A method for digitally reconstructing color images of gravity and magnetic anomalies, characterized in that: The method for digitally reconstructing gravity and magnetic anomaly color images comprises: Obtain color images of gravity and magnetic anomalies; Determine the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection datum plane, and a projection plane unit; Converting the gravity and magnetic anomaly color image into the actual geographic coordinate system using pixel points to obtain anomaly grid data; Determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data; Editing the sample point grid data according to the abnormal grid data to obtain initial contour lines of the gravity and magnetic anomaly color image; The abnormal grid data is reconstructed according to the actual abnormal value of the gravity and magnetic anomaly color image and the initial isovalue line to complete the reconstruction of the gravity and magnetic anomaly color image.

2. The method for digitally reconstructing color images of gravity and magnetic anomalies according to claim 1, characterized in that: The abnormal grid data includes the initial abnormal value of the abnormal point and the actual coordinates of the abnormal point in the actual geographic coordinate system; The method of converting the gravity and magnetic anomaly color image into the actual geographic coordinate system using pixel points to obtain anomaly grid data includes: Selecting a plurality of pixel points of the gravity and magnetic anomaly color image as image control points; Determining the actual coordinates of the image control point according to the pixel coordinates of the image control point, and determining the conversion relationship between the pixel coordinates of the image control point and the actual coordinates; According to the transformation relationship, all pixel points in the gravity and magnetic anomaly color image are transformed into an actual geographic coordinate system to obtain the actual coordinates of the anomaly points; wherein the anomaly points correspond to the pixel points one to one; Calculate the initial outlier value of the outlier point according to the RGB value of the pixel point; The abnormal grid data is obtained based on the actual coordinates of the abnormal point and the initial abnormal value of the abnormal point.

3. The method for digitally reconstructing color images of gravity and magnetic anomalies according to claim 2, characterized in that: The calculation formula of the initial outlier value of the outlier point is: A N =0.299*R+0.587*G+0.114*B; Among them, A N is the initial outlier value of the outlier point, R is the component value of the red channel, G is the component value of the green channel, and B is the component value of the blue channel.

4. The method for digitally reconstructing color images of gravity and magnetic anomalies according to claim 2, characterized in that: The sample point grid data includes the actual coordinates of the discrete sample points and the initial abnormal values ​​of the discrete sample points; The method of determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image and obtaining sample grid data of the discrete sample points by using the discrete sample points and the anomaly grid data comprises: Determining discrete sample points according to the mapping scale, preset point spacing, and preset line spacing of the gravity and magnetic anomaly color image; Determine the actual coordinates of the discrete sample points according to the position of the gravity and magnetic anomaly color image in the actual geographic coordinate system; Based on the abnormal grid data, the radial basis function method is used to interpolate the GRD data of the abnormal point; Extracting initial outlier values ​​of the discrete sample points from the GRD data of the outlier points using the actual coordinates of the discrete sample points; The sample point grid data of the discrete sample points are obtained based on the actual coordinates of the discrete sample points and the initial abnormal values ​​of the discrete sample points.

5. The method for digitally reconstructing color images of gravity and magnetic anomalies according to claim 4, characterized in that: Editing the sample point grid data according to the abnormal grid data to obtain the initial contour lines of the gravity and magnetic anomaly color image includes: Draw the contour lines of the outlier points according to the GRD data of the outlier points; Filtering the initial outliers of the discrete sample points; Based on the actual coordinates of the discrete sample points and the initial outlier values ​​of the filtered discrete sample points, the improved Shepard method is used to interpolate and obtain the GRD data of the discrete sample points; Draw the contour lines of the discrete sample points according to the GRD data of the discrete sample points; The contour lines of the discrete sample points are edited according to the contour lines of the abnormal points to obtain the initial contour lines of the gravity and magnetic anomaly color image.

6. The method for digitally reconstructing color images of gravity and magnetic anomalies according to claim 1, characterized in that: Reconstructing the abnormal grid data according to the actual abnormal value of the gravity and magnetic anomaly color image and the initial contour line includes: Re-assigning the initial contour lines according to the actual abnormal values ​​of the gravity and magnetic anomaly color image to obtain contour lines of the actual abnormal values; The abnormal grid data is reconstructed based on the contour lines of the abnormal true value.

7. A device for digital reconstruction of color images of gravity and magnetic anomalies, characterized in that: The device for digitally reconstructing color images of gravity and magnetic anomalies comprises: An acquisition module is used to obtain color images of gravity and magnetic anomalies; A coordinate system determination module is used to determine the actual geographic coordinate system of the gravity and magnetic anomaly color image according to projection parameters; wherein the projection parameters include a projection ellipsoid, a projection zone number, a projection reference plane, and a projection plane unit; A conversion module is used to convert the gravity and magnetic anomaly color image into an actual geographic coordinate system using pixel points to obtain anomaly grid data; a discrete module for determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample grid data of the discrete sample points using the discrete sample points and the anomaly grid data; An editing module, configured to edit the sample point grid data according to the abnormal grid data to obtain initial contour lines of the gravity and magnetic anomaly color image; A reconstruction module is used to reconstruct the abnormal grid data according to the actual abnormal value of the gravity and magnetic abnormal color image and the initial contour line.

8. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for digitally reconstructing color images of gravity and magnetic anomalies according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for digitally reconstructing a gravity and magnetic anomaly color image according to any one of claims 1 to 6 is implemented.

10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the method for digitally reconstructing a gravity and magnetic anomaly color image according to any one of claims 1 to 6 is implemented.

Images