Three-dimensional visualization method and system for real-time mapping of multiple types of geochemical anomaly data

Abstract

The present invention discloses a three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data, comprising the following steps: S1, interpolating and gridding a certain original geochemical element data to be mapped to generate an image IMG, and binding the image IMG as a texture to a preset mapping model; S2, establishing a color mapping palette image with a height of at least 1 pixel, wherein the leftmost side of the image represents the color of the minimum chemical element intensity and the rightmost side represents the color of the maximum chemical element intensity; S3, creating a fragment shader in a GPU, wherein the fragment shader uses a model UV to sample a value W of the image IMG, the sample value being the chemical element intensity value, and the sample value W being the input value IN of the fragment shader; S4, rendering the mapping model using the fragment shader to obtain a color-mapped three-dimensional visualization image. The present invention can adjust the color mapping method in real time, and the adjustment result can be seen in a short time, which is convenient and fast.

Description

Technical Field

[0001] The present invention relates to the field of geographic data image display, and in particular to a three-dimensional visualization method and system for real-time mapping of multiple types of geochemical anomaly data. Background Art

[0002] During geological data processing, most data ultimately results in maps. For example, geochemical data processing involves processing, analyzing, and interpreting geochemical data to extract information relevant to geological prospecting and other useful information during geochemical exploration and geochemical research. This processing typically involves using statistical tools to screen collected sample data, identify sample points with "anomalies," and delineate the extent of these anomalies to create an anomaly map. Areas containing anomalies represent areas with potential for prospecting. Similar processes include remote sensing and geophysical anomaly processing and mapping.

[0003] Geological data maps generated using traditional methods are fixed and unchanging. Users typically judge the mineralization potential of an area based on its color and area. However, these areas are affected by pre-selected parameters. If inappropriate image parameters are selected, users may misjudge the scope and shape of the mineralization potential area, missing out on important prospecting information. Summary of the Invention

[0004] The main purpose of the present invention is to provide a 3D visualization method and system for real-time mapping of multiple types of geochemical anomaly data, in which users can adjust the color of 3D visualization images in real time and see the adjustment results in a very short time.

[0005] The technical solution adopted in the present invention is:

[0006] A three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data is provided, comprising the following steps:

[0007] S1. Interpolate and grid the original geochemical element data to be mapped to generate an image IMG, and bind the image IMG as a texture to a preset mapping model;

[0008] S2. Create a color mapping palette image with a height of at least 1 pixel, wherein the leftmost side of the color represents the minimum intensity of the chemical element, and the rightmost side of the color represents the maximum intensity of the chemical element;

[0009] S3. Create a fragment shader in the GPU. The fragment shader uses the model UV to sample the image IMG as W, which is the intensity value of the chemical element. The sample value W is used as the input value IN of the fragment shader. The execution logic of the fragment shader is as follows:

[0010] 1) Set Calculate U = f(IN; MIN, MAX), where MIN and MAX are determined based on the values entered by the user;

[0011] 2) Using (U, 0) as the UV coordinate value, sample the color map palette image and output the resulting color;

[0012] S4. Use a fragment shader to render the image model to obtain a three-dimensional visualization image after color mapping.

[0013] Following the above technical solution, users can dynamically adjust MIN and MAX according to their needs and obtain the re-color-mapped 3D visualization image in real time.

[0014] Following the above technical solution, the number of channels of the image IMG is at least 1, and the data type of the channel should be a floating point number.

[0015] The present invention also provides a three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data, comprising the following steps:

[0016] S1. Interpolate and grid the original geochemical element data to be mapped to generate an image IMG, and bind the image IMG as a texture to a preset mapping model;

[0017] S2. Sort the original geochemical element data in ascending order according to the element intensity value to obtain array A, and count the number of elements as C, the minimum value VMin, and the maximum value Vmax;

[0018] S3. Create an image LUT_1 with a height of at least 1 pixel, and use the image LUT_1 to store percentiles.

[0019] S4. Traverse each element in array A and perform the following operations on each element:

[0020] a) Let the element value be X, and calculate the normalized element value

[0021] Xn=(X-VMin) / (VMax-VMin);

[0022] b) Let the sequence number of the element in the array be N, and calculate the normalized sequence number Nn = N / C, where the starting sequence number of the elements in the array is 1;

[0023] c) Set the value of each pixel with U coordinate Xn in image LUT_1 to Nn;

[0024] S5. Traverse each pixel in LUT_1 and perform the following operations on each pixel:

[0025] a) Let the pixel's U coordinate be X and its value be P. If P is greater than or equal to 0, skip; if it is less than 0, proceed to the next step.

[0026] b) Find the first pixel to the left of the pixel that is not less than 0, whose value is PMin and whose U coordinate is

[0027] UMin.

[0028] c) Find the first pixel to the right of the pixel that is not less than 0, whose value is PMax and whose U coordinate is

[0029] UMax.

[0030] d) Assign P = PMin + (X - UMin) / (UMax - UMin) * (PMax - PMin);

[0031] S6. The user provides a color mapping palette image LUT_2, which has a height of at least 1 pixel. The leftmost side of the image LUT_2 represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value.

[0032] S7, the user inputs the minimum percentile MIN and maximum percentile MAX of the desired map;

[0033] S8. Create a 3D model based on the user's actual situation, and bind IMG as a texture to the 3D model.

[0034] S9. Create a fragment shader in the GPU, pass VMin and VMax into the shader as uniform variables, and set the shader to use the model UV to sample the IMG value as W. The main logic of the shader is as follows:

[0035] a) Calculate the normalized value Wn = (W - VMin) / (VMax - VMin);

[0036] b) Using (Wn, 0) as UV, sample LUT_1 to obtain the corresponding normalized percentile Pt;

[0037] c) Set Calculate U = f(Pt;MIN,MAX);

[0038] d) Use (U, 0) as the UV coordinate, sample the image LUT_2, and output the resulting color;

[0039] S10. Rendering the three-dimensional model using a fragment shader to obtain a three-dimensional visualization image after color mapping.

[0040] Following the above technical solution, users can dynamically adjust MIN and MAX according to their needs and obtain the re-color-mapped 3D visualization image in real time.

[0041] Following the above technical solution, the pixel format in the image LUT_1 is Float16 or Float32, and the pixel value is initialized to -1.

[0042] The present invention also provides a three-dimensional visualization system for real-time mapping of multiple types of geochemical anomaly data, comprising:

[0043] The interpolation module is used to interpolate and grid the original geochemical element data to be mapped, generate an image IMG, and bind the image IMG as a texture to the preset mapping model;

[0044] The palette construction module is used to create a color mapping palette image with a height of at least 1 pixel. The leftmost side of the palette image represents the color corresponding to the minimum intensity value MIN of the chemical element, and the rightmost side represents the color corresponding to the maximum intensity MAX of the chemical element. The values of MIN and MAX are determined according to the values input by the user.

[0045] The shader construction module is used to create a fragment shader. The fragment shader uses the model UV to sample the image IMG as W, which is the intensity value of the chemical element. The sample value W is used as the input value IN of the fragment shader. The execution logic of the fragment shader is:

[0046] 1) Set Calculate U = f(IN; MIN, MAX);

[0047] 2) Using (U, 0) as the UV coordinate value, sample the color map palette image and output the resulting color;

[0048] The rendering module is used to render the image model using a fragment shader to obtain a three-dimensional visual image after color mapping.

[0049] Following the above technical solution, the number of channels of the image IMG is at least 1, and the data type of the channel should be a floating point number.

[0050] Following the above technical solution, the system also includes a setting module for obtaining the MIN and MAX values input by the user, and dynamically adjusting the MIN and MAX according to user needs to obtain a three-dimensional visualization image after re-color mapping in real time.

[0051] The present invention also provides a computer storage medium storing a computer program executable by a processor, wherein the computer program executes the three-dimensional visualization method for multi-type real-time mapping of geochemical anomaly data described in the above technical solution.

[0052] The beneficial effects of the present invention are as follows: after interpolating and gridding traditional geological data, especially geochemical data, the present invention uses a GPU to display it in three-dimensional space according to a user-defined color mapping method, and the user can adjust the color mapping method in real time and see the adjustment results in a short time, which is convenient and fast.

[0053] Furthermore, the method of this patent can generate both single-element geochemical maps and element combination geochemical maps, etc., to achieve multi-type three-dimensional display of geochemical anomaly data. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following is a brief introduction to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0055] Figure 1 This is a flow chart of a three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data according to Example 1 of the present invention;

[0056] Figure 2 This is a diagrammatic flow chart of Example 2 of the present invention;

[0057] Figure 3 is the TIFF image outputted in step 1 of embodiment 3 of the present invention;

[0058] Figure 4 is the palette image generated in step 2 of embodiment 3 of the present invention;

[0059] Figure 5 This is the effect diagram after real-time rendering of Example 3 of the present invention;

[0060] Figure 6 This is a rendering after adjusting parameters in real time according to embodiment 3 of the present invention;

[0061] Figure 7 It is a diagrammatic flow chart of a three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data according to Example 4 of the present invention. DETAILED DESCRIPTION

[0062] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0063] Example 1

[0064] like Figure 1As shown, the three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data in this embodiment includes the following steps:

[0065] S1. Interpolate and grid the original geochemical element data to be mapped to generate an image IMG, and bind the image IMG as a texture to a preset mapping model;

[0066] S2. Create a color mapping palette image with a height of at least 1 pixel, wherein the leftmost side of the color represents the minimum intensity of the chemical element, and the rightmost side of the color represents the maximum intensity of the chemical element;

[0067] S3. Create a fragment shader in the GPU. The fragment shader uses the model UV to sample the image IMG as W, which is the intensity value of the chemical element. The sample value W is used as the input value IN of the fragment shader. The execution logic of the fragment shader is as follows:

[0068] 1) Set Calculate U = f(IN; MIN, MAX), where MIN and MAX are determined based on the values entered by the user;

[0069] 2) Using (U, 0) as the UV coordinate value, sample the color map palette image and output the resulting color;

[0070] S4. Use a fragment shader to render the image model to obtain a three-dimensional visualization image after color mapping.

[0071] This embodiment interpolates and grids traditional geological data, especially geochemical data, and uses a GPU to display it in three-dimensional space according to a user-defined color mapping method. The user can adjust the color mapping method in real time by inputting the minimum value (MIN) and maximum value (MAX) of the data for the desired map, that is, the maximum and minimum values of the chemical element intensity value when the fragment shader is normalized. The adjustment results can be seen in a very short time (within 1 second), which is convenient and fast.

[0072] Example 2

[0073] This embodiment is based on embodiment 1. Figure 2 As shown, this embodiment includes the following steps:

[0074] 1. Interpolate and grid the data to be mapped to generate the image IMG. The image must have at least 1 channel, and the channel data type must be floating point.

[0075] 2. The user provides a color map palette image with a height of at least 1 pixel. The leftmost side of the image represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value.

[0076] 3. The user enters the minimum value MIN and the maximum value MAX of the data to be mapped.

[0077] 4. Create the model that the final image depends on based on the user's actual situation, and bind IMG as a texture to the model.

[0078] 5. Create a fragment shader (also called pixel shader), set the shader input value to IN, and the main logic of the shader is as follows:

[0079] a. Calculate U = f(IN; MIN, MAX)

[0080] b. Using (U, 0) as UV, sample the palette image and output the resulting color.

[0081] 6. Use the above fragment shader to render the model in step 4 to obtain the color mapped data.

[0082] As you can see, users can dynamically adjust the MIN and MAX values to achieve real-time color mapping of the map. By adjusting the maximum and minimum values, users can quickly identify high-value areas (anomalies) in the map during the interactive process. Based on the morphological changes of the anomaly areas, users can understand the distribution of data within the area, thereby assisting in the final selection of the appropriate prospecting target area.

[0083] Example 3

[0084] This embodiment is based on Example 2, taking the copper element geochemical intensity map of the Juneau area as an example, and uses the method of this patent to generate a dynamic image:

[0085] 1. Use Pandas to read the Excel spreadsheet of geochemical exploration data in the Juneau area, and then use GDAL's gdal_grid to perform grid interpolation on the copper element intensity values. The output is a TIFF image, single channel,

[0086] Float32 type, resolution 512×391. Figure 3 As shown, the data generated in step 1 is directly observed using the image viewing tool DJV. Since the original data will be overexposed, the exposure value is set to EV = -10 during observation.

[0087] 2. Using OpenCV, generate a palette TIFF image based on COLORMAP_JET with dark blue on the left and dark red on the right (such as Figure 4 As shown), RGB three channels, Uint32 type, resolution 512×16.

[0088] 3. Use GLSL to write a fragment shader as follows, where texture[0] is bound to the element intensity image generated in step 1, and texture[1] is bound to the palette image generated in step 2:

[0089]

[0090] 4. The steps of creating polygonal planes and binding fragment shader variables are consistent with the traditional GLSL usage.

[0091] I will omit this here. The image after real-time rendering is as follows Figure 5 As shown, Min=5, Max=8457.

[0092] 5. When the user modifies min and max, the color range of the output image will change in real time, such as Figure 6 As shown, where Min=213, Max=1220.

[0093] Example 4

[0094] The inventive concept of this embodiment is similar to that of embodiment 1, except that the chemical element intensity values are converted into quantiles in advance.

[0095] The three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data mainly includes the following steps:

[0096] 1. Interpolate and grid the raw geochemical element data to generate an image (IMG). The image must have at least one channel, and the channel data type must be floating point. This step converts the geochemical element data into an image file through interpolation, allowing the GPU to process the data as a texture.

[0097] 2. Sort the raw geochemical element data in ascending order by element intensity to create array A. Count the number of elements as C, with the minimum value VMin and the maximum value VMax. This step prepares the data for generating the next image.

[0098] 3. Create an image (LUT_1) with a height of at least 1 pixel. The image width W can be arbitrary, but a larger width will produce smoother color transitions. This image is used to store percentiles. The pixel format is Float16 or Float32, and the pixel values are initialized to -1. This image is used to establish the correspondence between element values and percentiles. The image's U coordinate is the element value normalized between VMin and VMax, and the pixel value is the corresponding percentile. The image is initialized to -1 because the normalized element values in array A may not correspond to all pixels in LUT_1, resulting in missing pixels, and we need to be able to identify missing pixels. For example, in the array 0, 10, 10, 10, 10, where 0 corresponds to U coordinate 0 and 10 corresponds to U coordinate 1, an image with a width of 3 will only have the leftmost and rightmost pixels filled in, skipping the pixels in the middle. Because we initialized it with -1, we can identify these missing pixels in the following steps and fill them using interpolation.

[0099] 4. Traverse each element in array A and perform the following operations on each element:

[0100] a. Let the element value be X, and calculate the normalized element value

[0101] Xn=(X-VMin) / (VMax-VMin). This step determines the image U coordinate corresponding to the value.

[0102] b. Let the sequence number of the element in array A be N (the starting sequence number of elements in array A is 1), and calculate the normalized sequence number Nn = N / C. This step determines the percentile corresponding to the value.

[0103] c. Set the value of each pixel with U coordinate Xn in LUT_1 to Nn, thereby filling the pixels.

[0104] 5. Traverse each pixel in LUT_1 and perform the following operations on each pixel (mainly used to find vacant pixels and interpolate to fill them):

[0105] a. Let the pixel's U coordinate be X and its value be P. If P is greater than or equal to 0, skip it.

[0106] If it is less than 0, proceed to the next step.

[0107] b. Find the first pixel to the left of the pixel that is not less than 0. Its value is PMin and its U coordinate is UMin.

[0108] c. Find the first pixel to the right of the pixel that is not less than 0. Its value is PMax and its U coordinate is UMax.

[0109] d. Assign value P = PMin + (X - UMin) / (UMax - UMin) * (PMax - PMin). This step mainly performs linear interpolation based on the effective values around the missing pixel.

[0110] 6. The user provides a color mapping palette image (LUT_2) with a height of at least 1 pixel. The leftmost side of the image represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value. This image (LUT_2) is mainly used to color the percentiles.

[0111] 7. The user enters the minimum (MIN) and maximum (MAX) percentile values of the desired graph.

[0112] 8. Create the model that the final image depends on based on the user's actual situation, and bind IMG as a texture to the model.

[0113] 9. Create a fragment shader (also called a pixel shader), pass VMin and VMax into the shader as uniform variables, and set the shader to use the model UV to sample the IMG value as W. The main logic of the shader is as follows:

[0114] a. Calculate the normalized value Wn = (W-VMin) / (VMax-VMin).

[0115] b. Using (Wn, 0) as the UV, sample LUT_1 to obtain the corresponding normalized percentile Pt. (LUT_1 is used to convert element values in geochemical images into percentiles in real time.)

[0116] c. Set Calculate U = f(Pt;MIN,MAX)

[0117] (The percentiles are clipped and normalized based on user input so they can be mapped to the palette image LUT_2.)

[0118] d. Use (U, 0) as UV, sample LUT_2, and output the resulting color.

[0119] 10. Use the above fragment shader to render the model in step 4 to obtain the color mapped data.

[0120] Example 5

[0121] This system embodiment is mainly used to implement method embodiments 1, 2, and 3. The three-dimensional visualization system for multi-type real-time mapping of geochemical anomaly data in this embodiment includes:

[0122] The interpolation module is used to interpolate and grid the original geochemical element data to be mapped, generate an image IMG, and bind the image IMG as a texture to the preset mapping model;

[0123] The palette construction module is used to create a color mapping palette image with a height of at least 1 pixel. The leftmost side of the palette image represents the color corresponding to the minimum intensity value MIN of the chemical element, and the rightmost side represents the color corresponding to the maximum intensity MAX of the chemical element. The values of MIN and MAX are determined according to the values input by the user.

[0124] The shader construction module is used to create a fragment shader. The fragment shader uses the model UV to sample the image IMG as W, which is the intensity value of the chemical element. The sample value W is used as the input value IN of the fragment shader. The execution logic of the fragment shader is:

[0125] 1) Set Calculate U = f(IN; MIN, MAX);

[0126] 2) Using (U, 0) as the UV coordinate value, sample the color map palette image and output the resulting color;

[0127] The rendering module is used to render the image model using a fragment shader to obtain a three-dimensional visual image after color mapping.

[0128] Example 6

[0129] This system embodiment is mainly used to implement method embodiment 4. The three-dimensional visualization system for multi-type real-time mapping of geochemical anomaly data in this embodiment includes:

[0130] The interpolation module is used to interpolate and grid the original geochemical element data to be mapped, generate an image IMG, and bind the image IMG as a texture to the preset mapping model;

[0131] The percentile image creation module is used to create an image LUT_1 with a height of at least 1 pixel. The image LUT_1 is used to store percentiles. It is also used to traverse each pixel in LUT_1 and perform the following operations on each pixel:

[0132] a) Let the pixel's U coordinate be X and its value be P. If P is greater than or equal to 0, skip; if it is less than 0, proceed to the next step.

[0133] b) Find the first pixel to the left of the pixel that is not less than 0, whose value is PMin and whose U coordinate is

[0134] UMin.

[0135] c) Find the first pixel to the right of the pixel that is not less than 0, whose value is PMax and whose U coordinate is

[0136] UMax.

[0137] d) Assign P = PMin + (X - UMin) / (UMax - UMin) * (PMax - PMin);

[0138] The array setting module is used to sort the original geochemical element data in ascending order according to the element intensity value to obtain array A, and count the number of elements as C, the minimum value VMin, and the maximum value Vmax; and is used to traverse each element in array A and perform the following operations on each element:

[0139] a) Let the element value be X, and calculate the normalized element value Xn = (X-VMin) / (VMax-VMin);

[0140] b) Let the sequence number of the element in the array be N, and calculate the normalized sequence number Nn = N / C, where

[0141] The starting number of the elements in the array is 1;

[0142] c) Set the value of each pixel with U coordinate Xn in image LUT_1 to Nn;

[0143] The palette creation module is used to provide a color mapping palette image LUT_2 with a height of at least 1 pixel; the leftmost side of the image LUT_2 represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value;

[0144] The setting module is used to input the minimum value MIN and the maximum value MAX of the percentile of the expected map by the user;

[0145] The 3D model creation module is used to create the 3D model that the final drawing depends on according to the user's actual situation, and bind the IMG as a texture to the 3D model;

[0146] The fragment shader creation module is used to create a fragment shader in the GPU, pass VMin and VMax as uniform variables to the shader, and set the shader to use the model UV to sample the IMG value as W. The main logic of the shader is as follows:

[0147] a) Calculate the normalized value Wn = (W - VMin) / (VMax - VMin);

[0148] b) Using (Wn, 0) as UV, sample LUT_1 to obtain the corresponding normalized percentile Pt;

[0149] c) Set Calculate U = f(Pt;MIN,MAX);

[0150] d) Use (U, 0) as the UV coordinate, sample the image LUT_2, and output the resulting color;

[0151] The rendering module is used to render the three-dimensional model using a fragment shader to obtain a three-dimensional visual image after color mapping.

[0152] Example 7

[0153] This application also provides a computer-readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card-type memory (e.g., SD or DX memory), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic storage device, a disk, an optical disk, a server, an App store, etc., storing a computer program that, when executed by a processor, implements a corresponding function. The computer-readable storage medium of this embodiment, when executed by a processor, implements a three-dimensional visualization method for real-time mapping of multiple types of geochemical exploration anomaly data in the aforementioned multiple method embodiments.

[0154] It should be pointed out that, according to the needs of implementation, the various steps / components described in this application can be split into more steps / components, or two or more steps / components or partial operations of steps / components can be combined into new steps / components to achieve the purpose of the present invention.

[0155] The size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0156] It should be understood that those skilled in the art can make improvements or changes based on the above description, and all such improvements and changes should fall within the scope of protection of the appended claims of the present invention.

Claims

1. A three-dimensional visualization method for real-time mapping of multiple types of geochemical anomaly data, characterized by: The following steps are involved: S1. Interpolate and grid the original geochemical element data to be mapped to generate an image IMG, and bind the image IMG as a texture to a preset mapping model; S2. Sort the original geochemical element data in ascending order according to the element intensity value to obtain array A, and count the number of elements as C, the minimum value VMin, and the maximum value VMax; S3. Create an image LUT_1 with a height of at least 1 pixel, and use the image LUT_1 to store percentiles. S4. Traverse each element in array A and perform the following operations on each element: a. Let the element value be X, and calculate the normalized element value Xn = (X-VMin) / (VMax-VMin); b. Let the element's sequence number in the array be N, and calculate the normalized sequence number Nn = N / C, where the starting sequence number of the element in the array is 1; c. Set the value of each pixel with U coordinate Xn in image LUT_1 to Nn; S5. Traverse each pixel in LUT_1 and perform the following operations on each pixel: a) Let the pixel's U coordinate be X and its value be P. If P is greater than or equal to 0, skip; if it is less than 0, proceed to the next step. b) Find the first pixel to the left of the pixel that is not less than 0, whose value is PMin and whose U coordinate is UMin; c) Find the first pixel to the right of the pixel that is not less than 0, whose value is PMax and whose U coordinate is UMax; d) Assign P = PMin + (X - UMin) / (UMax - UMin) · (PMax - PMin); S6. The user provides a color mapping palette image LUT_2, which has a height of at least 1 pixel. The leftmost side of the image LUT_2 represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value. S7, the user inputs the minimum percentile MIN and maximum percentile MAX of the desired map; S8. Create a 3D model based on the user's actual situation, and bind IMG as a texture to the 3D model. S9. Create a fragment shader in the GPU, pass VMin and VMax into the shader as uniform variables, and set the shader to use the model UV to sample the IMG value as W. The main logic of the shader is as follows: a) Calculate the normalized value Wn = (W-VMin) / (VMax-VMin); b) Use As UV, sample LUT_1 and obtain the corresponding normalized percentile Pt; c) Set ,calculate ; d) Use As UV coordinates, sample the image LUT_2 and output the resulting color; S10. Rendering the three-dimensional model using a fragment shader to obtain a three-dimensional visualization image after color mapping.

2. The three-dimensional visualization method for multi-type real-time mapping of geochemical anomaly data according to claim 1 is characterized in that: Users can dynamically adjust MIN and MAX according to their needs and obtain the re-color-mapped 3D visualization map in real time.

3. The three-dimensional visualization method for multi-type real-time mapping of geochemical anomaly data according to claim 1 is characterized in that: The pixel format in image LUT_1 is Float16 or Float32, and the pixel value is initialized to -1.

4. A three-dimensional visualization system for real-time mapping of multiple types of geochemical anomaly data, characterized by: include: The interpolation module is used to interpolate and grid the original geochemical element data to be mapped, generate an image IMG, and bind the image IMG as a texture to the preset mapping model; The percentile image creation module is used to create an image LUT_1 with a height of at least 1 pixel. The image LUT_1 is used to store percentiles. It is also used to traverse each pixel in LUT_1 and perform the following operations on each pixel: a) Let the pixel's U coordinate be X and its value be P. If P is greater than or equal to 0, skip; if it is less than 0, proceed to the next step. b) Find the first pixel to the left of the pixel that is not less than 0, whose value is PMin and whose U coordinate is UMin; c) Find the first pixel to the right of the pixel that is not less than 0, whose value is PMax and whose U coordinate is UMax; d) Assign P = PMin + (X - UMin) / (UMax - UMin) * (PMax - PMin); The array setting module is used to sort the original geochemical element data in ascending order according to the element intensity value to obtain array A, and count the number of elements as C, the minimum value VMin, and the maximum value VMax; and is used to traverse each element in array A and perform the following operations on each element: a) Let the element value be X, and calculate the normalized element value Xn = (X-VMin) / (VMax-VMin); b) Let the sequence number of the element in the array be N, and calculate the normalized sequence number Nn = N / C, where the starting sequence number of the elements in the array is 1; c) Set the value of each pixel with U coordinate Xn in image LUT_1 to Nn; The palette creation module is used to provide a color mapping palette image LUT_2 with a height of at least 1 pixel; the leftmost side of the image LUT_2 represents the color of the minimum data value, and the rightmost side represents the color of the maximum data value; The setting module is used to input the minimum value MIN and the maximum value MAX of the percentile of the expected map by the user; The 3D model creation module is used to create the 3D model that the final drawing depends on according to the user's actual situation, and bind the IMG as a texture to the 3D model; The fragment shader creation module is used to create a fragment shader in the GPU, pass VMin and VMax as uniform variables to the shader, and set the shader to use the model UV to sample the IMG value as W. The main logic of the shader is as follows: a) Calculate the normalized value Wn = (W-VMin) / (VMax-VMin); b) Use As UV, sample LUT_1 and obtain the corresponding normalized percentile Pt; c) Set ,calculate ; d) Use As UV coordinates, sample the image LUT_2 and output the resulting color; The rendering module is used to render the three-dimensional model using a fragment shader to obtain a three-dimensional visual image after color mapping.

5. The three-dimensional visualization system for real-time mapping of multiple types of geochemical anomaly data according to claim 4 is characterized in that: Users can dynamically adjust MIN and MAX according to their needs and obtain the re-color-mapped 3D visualization map in real time.

6. The three-dimensional visualization system for real-time mapping of multiple types of geochemical anomaly data according to claim 4 is characterized in that: The pixel format in image LUT_1 is Float16 or Float32, and the pixel value is initialized to -1.

7. A computer storage medium, characterized in that A computer program executable by a processor is stored therein, and the computer program executes the three-dimensional visualization method for multi-type real-time mapping of geochemical anomaly data according to any one of claims 1 to 3.

Images