Gravity and magnetic anomaly color image digital reconstruction method, device, equipment, medium and product

Abstract

The invention discloses a gravity and magnetic anomaly color image digital reconstruction method and device, equipment, a medium and a product, and relates to the technical field of geophysical exploration, and the method comprises the steps: obtaining a gravity and magnetic anomaly color image; determining an actual geographic coordinate system according to the projection parameters; converting the gravity and magnetic abnormal color image into an actual geographic coordinate system by using the pixel points to obtain abnormal grid data; determining discrete sample points according to the charting proportion of the gravity and magnetic abnormal color image, and obtaining sample point grid data of the discrete sample points by using the discrete sample points and the abnormal grid data; editing the sampling point grid data according to the abnormal grid data to obtain an initial contour line of the gravity and magnetic abnormal color image; and reconstructing the abnormal grid data according to the actual abnormal value and the initial contour line of the gravity and magnetic abnormal color image to complete reconstruction of the gravity and magnetic abnormal color image. According to the method, the gravity and magnetic anomaly color image can be reconstructed, and gravity and magnetic anomaly data can be directly extracted from the reconstructed gravity and magnetic anomaly color image.

Description

Technical Field

[0001] The present application relates to the technical field of geophysical exploration, and particularly to a method, device, equipment, medium and product for digital reconstruction of gravity and magnetic anomaly color images. Background Art

[0002] In gravity exploration and magnetic exploration, the main forms of exploration results are gravity and magnetic data. To comprehensively analyze and study gravity and magnetic data, it is necessary to draw them in the form of images. The maps drawn with gravity and magnetic data are called gravity and magnetic anomaly (gravity anomaly, magnetic anomaly) maps. In actual work, the cross-section anomaly map and the plane anomaly map of gravity and magnetic anomalies are the most commonly used. Among them, the plane contour map of gravity and magnetic anomalies is an important basic map for further geological research and analysis. The change situation of gravity and magnetic anomalies is judged through aspects such as the sparse change, trend characteristics, and high and low amplitudes of the contour lines. In addition, since the expression form of the plane color image is more convenient for feature comparison and analysis of anomalies, the contour map and the color image are used in combination, and different display effects are presented by changing settings such as color depth, brightness, and color, which are increasingly favored and used by geophysical workers.

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

[0004] The purpose of the present application is to provide a method, device, equipment, medium and product for digital reconstruction of gravity and magnetic anomaly color images, which can reconstruct the 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 purpose, the present application provides the following solutions:

[0006] In the first aspect, the present application provides a method for digital reconstruction of gravity and magnetic anomaly color images, including:

[0007] Obtain a gravity and magnetic anomaly color image;

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

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

[0010] Determine discrete sample points according to the mapping scale of the gravity and magnetic anomaly color image, and obtain sample point grid data of the discrete sample points by using the discrete sample points and the anomaly grid data;

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

[0012] Reconstruct the abnormal grid data according to the actual anomaly values of the gravity and magnetic anomaly color image and the initial contour lines to complete the reconstruction of the gravity and magnetic anomaly color image.

[0013] Furthermore, the abnormal grid data includes the initial anomaly values of the abnormal points and the actual coordinates of the abnormal points in the actual geographic coordinate system;

[0014] Among them, converting the gravity and magnetic anomaly color image into the actual geographic coordinate system by using pixel points to obtain abnormal grid data, including:

[0015] Select multiple pixel points of the gravity and magnetic anomaly color image as image control points;

[0016] Determine the actual coordinates of the image control points according to the pixel coordinates of the image control points, and determine the conversion relationship between the pixel coordinates and the actual coordinates of the image control points;

[0017] Convert all pixel points in the gravity and magnetic anomaly color image into the actual geographic coordinate system according to the conversion relationship to obtain the actual coordinates of the abnormal points; among them, the abnormal points and the pixel points correspond one by one;

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

[0019] Obtain the abnormal grid data based on the actual coordinates of the abnormal points and the initial anomaly values of the abnormal points.

[0020] Furthermore, the calculation formula for the initial anomaly value of the abnormal point is:

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

[0022] Among them, A N is the initial anomaly value of the abnormal 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 the actual coordinates of the discrete sample points and the initial anomaly values of the discrete sample points;

[0024] Among them, determine the discrete sample points according to the mapping scale of the gravity and magnetic anomaly color image, and use the discrete sample points and the abnormal grid data to obtain the sample point grid data of the discrete sample points, including:

[0025] Determine discrete sample points according to the mapping scale, preset point distance, and preset line distance 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 anomaly grid data, use the radial basis function method to interpolate to obtain the GRD data of the anomaly points;

[0028] Extract the initial anomaly values of the discrete sample points from the GRD data of the anomaly points using the actual coordinates of the discrete sample points;

[0029] Obtain the sample point grid data of the discrete sample points based on the actual coordinates and initial anomaly values of the discrete sample points.

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

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

[0032] Filter the initial anomaly values of the discrete sample points;

[0033] Based on the actual coordinates of the discrete sample points and the filtered initial anomaly values of the discrete sample points, use the improved Shepard method to interpolate to 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] Edit the contour lines of the discrete sample points according to the contour lines of the anomaly points to obtain the initial contour lines of the gravity and magnetic anomaly color image.

[0036] Further, reconstruct the anomaly grid data according to the actual anomaly values of the gravity and magnetic anomaly color image and the initial contour lines, including:

[0037] Reassign the initial contour lines according to the actual anomaly values of the gravity and magnetic anomaly color image to obtain the contour lines of the true anomaly values;

[0038] Reconstruct the anomaly grid data based on the contour lines of the true anomaly values.

[0039] In a second aspect, the present application provides a device for digitizing and reconstructing a gravity and magnetic anomaly color image, including:

[0040] An acquisition module for acquiring a gravity and magnetic anomaly color image;

[0041] A coordinate system determination module, configured 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, configured to convert the gravity and magnetic anomaly color image into the actual geographic coordinate system by using pixel points to obtain anomaly grid data;

[0043] A discretization module, configured to determine discrete sample points according to the mapping scale of the gravity and magnetic anomaly color image, and obtain sample point grid data of the discrete sample points by 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 anomaly grid data to obtain an initial isoline of the gravity and magnetic anomaly color image;

[0045] A reconstruction module, configured to reconstruct the anomaly grid data according to the actual anomaly value of the gravity and magnetic anomaly color image and the initial isoline.

[0046] In a third aspect, the present application provides a computer device, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor executes the computer program to implement the above-mentioned gravity and magnetic anomaly color image digital reconstruction method.

[0047] In a fourth aspect, the present application provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the above-mentioned gravity and magnetic anomaly color image digital reconstruction method is implemented.

[0048] In a fifth aspect, the present application provides a computer program product, including a computer program, and when the computer program is executed by a processor, the above-mentioned gravity and magnetic anomaly color image digital reconstruction method is implemented.

[0049] According to the specific embodiments provided by the present application, the following technical effects are disclosed in the present application:

[0050] The present application provides a method, apparatus, device, medium and product for digital reconstruction of gravity and magnetic anomaly color images. The gravity and magnetic anomaly color images are obtained, and their actual geographic coordinate systems are determined according to the projection parameters. The gravity and magnetic anomaly color images are transformed into the actual geographic coordinate systems, that is, the gravity and magnetic anomaly color images are transformed into specific actual geographic coordinate systems according to the research needs of the gravity and magnetic anomaly data, so as to obtain anomaly grid data. Again, discrete sample points are determined according to the mapping scale, the gravity and magnetic anomaly color images are discretized, and sample point grid data are obtained by using the anomaly grid data and the discrete sample points. Finally, the sample point grid data are edited according to the anomaly grid data to obtain the initial contour lines of the gravity and magnetic anomaly color images. On this basis, the anomaly grid data are reconstructed by using the actual anomaly values and the initial contour lines, so that the reconstructed anomaly grid data have the true anomaly values, and the digital reconstruction of the gravity and magnetic anomaly color images is completed. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required to be used in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0052] Figure 1 It is an application environment diagram of a method for digital reconstruction of gravity and magnetic anomaly color images in an embodiment of the present application;

[0053] Figure 2 It is a schematic flowchart of a method for digital reconstruction of gravity and magnetic anomaly color images provided in an embodiment of the present application;

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

[0055] Figure 4 It is a combined diagram of the initial contour lines of the edited gravity and magnetic anomaly color image; among them, (a) is the combined diagram of the edited gravity and magnetic anomaly color image, and (b) is the combined diagram of the contour line images of the edited gravity and magnetic anomaly color image;

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

[0057] Figure 6 It is a schematic diagram of the functional modules of a device for digital reconstruction of gravity and magnetic anomaly color images provided in an embodiment of the present application;

[0058] Figure 7 A structural schematic diagram of a computer device provided by an embodiment of the present application. Specific implementation manners

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

[0060] To make the purpose, features, and advantages of the present application more obvious and understandable, the present application will be further described in detail below with reference to the accompanying drawings and specific implementation manners.

[0061] The method for digitizing and reconstructing the gravity and magnetic anomaly color image provided by the embodiment of the present application can be applied to an application environment as Figure 1 shown. Among them, the terminal 102 communicates with the server 104 through a network. The data storage system can store the data that the server 104 needs to process. The data storage system can be set separately, integrated on the server 104, placed on the cloud or other servers. The terminal 102 can send the gravity and magnetic anomaly color image to the server 104. After receiving the gravity and magnetic anomaly color image, for the gravity and magnetic anomaly color image, the server 104 determines the actual geographic coordinate system of the gravity and magnetic anomaly color image according to the projection parameters; wherein, the projection parameters include the projection ellipsoid, the projection zone number, the projection reference plane, and the projection plane unit; the gravity and magnetic anomaly color image is transformed into the actual geographic coordinate system by using pixel points to obtain anomaly grid data; discrete sample points are determined according to the mapping scale of the gravity and magnetic anomaly color image, and the sample point grid data of the discrete sample points is obtained by using the discrete sample points and the anomaly grid data; the initial isogram of the gravity and magnetic anomaly color image is obtained by editing the sample point grid data according to the anomaly grid data; the anomaly grid data is reconstructed according to the actual anomaly value and the initial isogram of the gravity and magnetic anomaly color image to complete the reconstruction of the gravity and magnetic anomaly color image. The server 104 can feedback the reconstructed gravity and magnetic anomaly color image to the terminal 102. In addition, in some embodiments, the method for digitizing and reconstructing the gravity and magnetic anomaly color image can also be implemented by the server 104 or the terminal 102 alone. For example, the terminal 102 can directly perform digitizing and reconstructing 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 digitizing and reconstructing processing on the gravity and magnetic anomaly color image.

[0062] Among them, the terminal 102 can be, but is not limited to, various desktop computers, laptop computers, smart phones, tablet computers, etc. The server 104 can be implemented by an independent server or a server cluster composed of multiple servers, and can also be a cloud server.

[0063] In an exemplary embodiment, as Figure 2 shown, a method for digital reconstruction of gravity and magnetic anomaly color images is provided. This method is executed by a computer device, and specifically can be executed alone by a computer device such as a terminal or a server, or can be jointly executed by a terminal and a server. In the embodiments of the present application, taking this method applied to Figure 1 the server 104 in it as an example for illustration, it includes the following steps 201 to step 206. Among them:

[0064] Step 201, obtain the gravity and magnetic anomaly color image.

[0065] Specifically, according to the actual research work situation, find the gravity and magnetic anomaly color image that needs to extract gravity and magnetic anomaly data from literature or reports, and store the gravity and magnetic anomaly color image in JPG format.

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

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

[0068] Step 204, determine discrete sample points according to the mapping scale of the gravity and magnetic anomaly color image, and use the discrete sample points and the anomaly grid data to obtain sample point grid data of the discrete sample points.

[0069] Step 205, edit the sample point grid data according to the anomaly grid data to obtain the initial isogram of the gravity and magnetic anomaly color image.

[0070] Step 206, reconstruct the anomaly grid data according to the actual anomaly value and the initial isogram of the gravity and magnetic anomaly color image to complete the reconstruction of the gravity and magnetic anomaly color image.

[0071] Executing steps 201 to 206 can digitally reconstruct the gravity and magnetic anomaly color image, 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. Then, 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 points according to the pixel coordinates of the image control points, and determine the conversion relationship between the pixel coordinates and the actual coordinates of the image control points.

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

[0076] Step 20304: Calculate the initial abnormal value of the abnormal point according to the RGB values of the pixel points.

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

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

[0079] wherein, A N is the initial abnormal value of the abnormal 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 and the initial abnormal value of the abnormal point.

[0081] In the specific implementation process, use Global Mapper software to convert the gravity and magnetic anomaly color image into the actual geographic coordinates. Specifically: Open the JPG image of the gravity and magnetic anomaly color image obtained in step 201 using GlobalMapper software, and set the projection parameters in Global Mapper software according to the research needs of the gravity and magnetic anomaly data, so as to determine the actual geographic coordinates for the conversion of the gravity and magnetic anomaly color image. On this basis, select multiple pixel points as image control points. The number of image control points is generally more than 4, aiming to accurately convert the image coordinates. The specific number is not limited in this embodiment.

[0082] Input the pixel coordinates of the selected pixel points and the actual coordinates of the corresponding image control points in the GlobalMapper software, so as to determine the conversion relationship between the pixel points and the image control points. The pixel points are converted into the actual geographical coordinates and are called abnormal points.

[0083] According to this conversion relationship, all the pixel points in the gravity and magnetic anomaly color image are converted into the actual geographical coordinates. Each pixel point corresponds to an abnormal point in the actual geographical coordinates, so as to obtain the actual coordinates of the abnormal points corresponding to all the pixel points, and obtain a JPG image file with actual coordinates. Then, combined with the initial anomaly values of the abnormal points, abnormal grid data is obtained.

[0084] In an exemplary embodiment, the sample point grid data includes the actual coordinates of discrete sample points and the initial anomaly values of 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 distance and preset line distance 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 for acquiring 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 distance and preset line distance can be set as the preset point distance × preset line distance of 500m × 500m, or as the preset point distance × preset line distance of 10,000m × 10,000m. Different preset point distances and preset line distances 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 are not limited in this embodiment.

[0087] Step 20402, determine the actual coordinates of discrete sample points according to the position of the gravity and magnetic anomaly color image in the actual geographical coordinate system.

[0088] Form a uniform grid map according to the preset point distance and preset line distance. Each point in the grid map is used as a discrete sample point, and finally multiple discrete sample points are determined.

[0089] The grid map where the discrete sample points are located overlaps with the gravity and magnetic anomaly color image in the actual geographical coordinate system. The corresponding actual coordinates of the discrete sample points within the area of the gravity and magnetic anomaly color image are used as the actual coordinates of the discrete sample points, and further the actual coordinates of the discrete sample points of the gravity and magnetic anomaly color image are determined.

[0090] Step 20403: Based on the abnormal grid data, use the radial basis function method to interpolate the GRD data of the abnormal points.

[0091] Specifically, based on the initial abnormal values of the abnormal points, use GlobalMapper software to convert and output the JPG image file with actual coordinates into a data file, and then use the radial basis function method in Surfer software to interpolate the data file into a GRD file to obtain the GRD data of the abnormal points.

[0092] Step 20404: Extract the initial abnormal values of the discrete sample points from the GRD data of the abnormal points using the actual coordinates of the discrete sample points.

[0093] The initial abnormal values in the GRD data change continuously. Take the initial abnormal value corresponding to the actual coordinates of the discrete sample points as the initial abnormal value of the discrete sample points.

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

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

[0096] Step 20501: Draw the isograms of the abnormal points according to the GRD data of the abnormal points.

[0097] Specifically, use the isogram drawing interface of Surfer software to draw the isograms of the abnormal points.

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

[0099] Specifically, use the symmetric small sub-domain filtering method (gradient-preserving filtering method) to filter the initial abnormality of the discrete sample points, which can better maintain the characteristics of the original abnormal field while filtering out noise.

[0100] Step 20503: Based on the actual coordinates of the discrete sample points and the filtered initial abnormal values of the discrete sample points, use the improved Shepard method to interpolate the GRD data of the discrete sample points.

[0101] Step 20504: Draw the isograms of the discrete sample points according to the GRD data of the discrete sample points.

[0102] Specifically, use the same method as in Step 20501 to draw the isograms of the discrete sample points.

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

[0104] The isolines of discrete sample points do not completely coincide with the isolines of abnormal points. Therefore, editing the isolines of discrete sample points according to the isolines of abnormal points so that the isolines of discrete sample points coincide with the isolines of abnormal points one by one can achieve a better correspondence.

[0105] Specifically, the images provided in the literature or report are divided into two cases. One is that in the original literature, there are not only the gravity and magnetic anomaly color images, but also the isoline images of the gravity and magnetic anomaly color images, that is, the isoline images of abnormal points; the other is that in the original literature, there is only the gravity and magnetic anomaly color image. For these two cases, the methods of isoline editing are roughly the same; however, when there are isoline images of abnormal points, the editing and adjustment will be faster and relatively easier.

[0106] As Figure 3 shown, Figure 3 (a) in is the original image of the gravity and magnetic anomaly color image provided in the literature or report, Figure 3 (b) in is the isoline image of the gravity and magnetic anomaly color image provided in the literature or report.

[0107] As Figure 4 shown, Figure 4 The combined image of the edited gravity and magnetic anomaly color image in (a) includes a color background image and a red curve. Among them, the color background image is the gravity and magnetic anomaly color image with actual coordinates, and the red curve is the initial isoline after editing.

[0108] Figure 4 The combined isoline image of the edited gravity and magnetic anomaly color image in (b) includes a gray background image, a black curve and a red curve. Among them, the image composed of the gray background image and the black curve is the isoline image of the gravity and magnetic anomaly color image with actual coordinates, and the red curve is the initial isoline after editing.

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

[0110] Step 20601, reassign the initial isoline according to the actual anomaly value of the gravity and magnetic anomaly color image to obtain the isoline of the true anomaly value.

[0111] Specifically, generally in the literature or report, the distribution of the true anomaly value of the gravity and magnetic anomaly color image is marked through the legend color bar of the image. Therefore, reassigning the initial isoline according to the actual anomaly value marked by the legend color bar of the gravity and magnetic anomaly color image can obtain the isoline of the true anomaly value of the gravity and magnetic anomaly color image.

[0112] Step 20602: Reconstruct the abnormal grid data based on the isolines of the abnormal true values.

[0113] Specifically, use Surfer software to output the isolines of the abnormal true values one by one to obtain a plurality of point data. Each point data includes the abscissa, ordinate in the actual physical coordinates, and the abnormal true value. Finally, continuously integrate all the output point data into an abnormal result file. During the actual operation process, before outputting the point data, check and approve the assignment results of the isolines of the abnormal true values according to the characteristics of the isolines to ensure the correctness and accuracy of the assignment results of each isoline. Among them, the characteristics of the isolines include that the abnormal true values of adjacent isolines are generally gradually changing, the abnormal true values of the isolines corresponding to the same color boundary are the same, and the abnormal true value of a closed color area gradually increases or decreases from the boundary to the center.

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

[0115] Based on the same inventive concept, the embodiment of the present application also provides a gravity and magnetic anomaly color image digital reconstruction device for implementing the above-mentioned gravity and magnetic anomaly color image digital reconstruction method. The solution provided by this device to solve the problem is similar to the solution described in the above method. Therefore, the specific limitations in one or more embodiments of the gravity and magnetic anomaly color image digital reconstruction device provided below can refer to the limitations on the gravity and magnetic anomaly color image digital reconstruction method in the above text, and will not be repeated here.

[0116] In an exemplary embodiment, as Figure 6 shown, a gravity and magnetic anomaly color image digital reconstruction device is provided, including:

[0117] An acquisition module 61, configured to acquire a gravity and magnetic anomaly color image.

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

[0119] A conversion module 63, configured to convert the gravity and magnetic anomaly color image into the actual geographic coordinate system by using pixel points to obtain abnormal grid data.

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

[0121] An editing module 65 is configured to edit the sample grid data according to the anomaly grid data to obtain an initial contour line of the gravity and magnetic anomaly color image.

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

[0123] In an exemplary embodiment, a computer device is provided. The computer device can be a server or a terminal, and its internal structure is as Figure 7 shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is configured 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 configured to store the gravity and magnetic anomaly color image. The input / output interface of the computer device is configured to exchange information between the processor and external devices. The communication interface of the computer device is configured to communicate with an external terminal through a network connection. When the computer program is executed by the processor, it implements a method for digital reconstruction of a gravity and magnetic anomaly color image.

[0124] Those skilled in the art can understand that Figure 7 the structure shown in

[0125] 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 those shown in the figure, or combine certain components, or have a different component layout.

[0126] In an exemplary embodiment, a computer device is further provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

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

[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 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 need to comply with relevant regulations.

[0129] Those of ordinary skill in the art can understand that all or part of the processes of implementing the above method embodiments can be completed by instructing relevant hardware through a computer program. 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 above method embodiments. Among them, any reference to a memory, database, or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memories. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0130] The databases involved in the embodiments provided in this application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., and are not limited thereto. The processors involved in the embodiments provided in this application can be general-purpose processors, central processors, graphics processors, digital signal processors, programmable logic devices, data processing logics based on quantum computing, etc., and are not limited thereto.

[0131] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, 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, it should be considered as the scope described in this specification.

[0132] Specific examples are used in this specification to elaborate on the principles and implementation manners of the present application. The descriptions of the above embodiments are only used to help understand the method and its core idea of the present application; at the same time, for those of ordinary skill in the art, according to the idea of the present application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A method for digital reconstruction of gravity and magnetic anomaly color images, characterized in that: The method for digitally reconstructing a gravity and magnetic anomaly color image 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 reference plane, and a projection plane unit; The gravity and magnetic anomaly color image is converted into the actual geographic coordinate system using pixel points to obtain anomaly grid data; Determine discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtain sample point grid data of the discrete sample points by using the discrete sample points and the anomaly grid data; 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; 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 digital reconstruction of gravity and magnetic anomaly color images according to claim 1 is characterized in that: The abnormal grid data includes an initial abnormal value of the abnormal point and the actual coordinates of the abnormal point in the actual geographic coordinate system; The pixel points are used to transform the gravity and magnetic anomaly color image into the actual geographic coordinate system to obtain anomaly grid data, including: Selecting a plurality of pixel points of the gravity and magnetic anomaly color image as image control points; Determine the actual coordinates of the image control point according to the pixel coordinates of the image control point, and determine the transformation 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 the actual geographic coordinate system to obtain the actual coordinates of the anomaly points; wherein the anomaly points correspond to the pixel points one by 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 digital reconstruction of gravity and magnetic anomaly color images according to claim 2 is 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 digital reconstruction of gravity and magnetic anomaly color images according to claim 2 is 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 comprises: determining discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtaining sample point grid data of the discrete sample points by using the discrete sample points and the anomaly grid data. Determine discrete sample points according to the mapping ratio, preset point spacing and preset line spacing of the gravity and magnetic anomaly color image; Determine the actual coordinates of 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, a 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 digital reconstruction of gravity and magnetic anomaly color images according to claim 4 is characterized in that: 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: Draw contour lines of the abnormal point according to the GRD data of the abnormal point; Filtering the initial outliers of the discrete sample points; Based on the actual coordinates of the discrete sample points and the initial outliers of the filtered discrete sample points, the GRD data of the discrete sample points are obtained by interpolation using the improved Shepard method; 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 digital reconstruction of gravity and magnetic anomaly color images according to claim 1 is 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 line according to the actual abnormal value of the gravity and magnetic anomaly color image to obtain the contour line of the actual abnormal value; 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 acquire color images of gravity and magnetic anomalies; A coordinate system determination module, 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 belt 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, used to determine discrete sample points according to the mapping ratio of the gravity and magnetic anomaly color image, and obtain sample point grid data of the discrete sample points by using the discrete sample points and the anomaly grid data; An editing module, used for 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; 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