CN116912438A - Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration - Google Patents
Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration Download PDFInfo
- Publication number
- CN116912438A CN116912438A CN202310802525.6A CN202310802525A CN116912438A CN 116912438 A CN116912438 A CN 116912438A CN 202310802525 A CN202310802525 A CN 202310802525A CN 116912438 A CN116912438 A CN 116912438A
- Authority
- CN
- China
- Prior art keywords
- value
- image
- color
- pixel
- mapping
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000013507 mapping Methods 0.000 title claims abstract description 76
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 26
- 238000007794 visualization technique Methods 0.000 title claims abstract description 17
- 239000000126 substance Substances 0.000 title claims abstract description 14
- 239000012634 fragment Substances 0.000 claims abstract description 47
- 229910052729 chemical element Inorganic materials 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims abstract description 24
- 238000009877 rendering Methods 0.000 claims abstract description 18
- 238000003384 imaging method Methods 0.000 claims abstract description 17
- 230000027455 binding Effects 0.000 claims abstract description 16
- 238000009739 binding Methods 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 16
- 238000012800 visualization Methods 0.000 claims description 13
- 230000000007 visual effect Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 5
- 230000001174 ascending effect Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000000523 sample Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 229910052802 copper Inorganic materials 0.000 description 2
- 239000010949 copper Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- MCSXGCZMEPXKIW-UHFFFAOYSA-N 3-hydroxy-4-[(4-methyl-2-nitrophenyl)diazenyl]-N-(3-nitrophenyl)naphthalene-2-carboxamide Chemical compound Cc1ccc(N=Nc2c(O)c(cc3ccccc23)C(=O)Nc2cccc(c2)[N+]([O-])=O)c(c1)[N+]([O-])=O MCSXGCZMEPXKIW-UHFFFAOYSA-N 0.000 description 1
- 208000025174 PANDAS Diseases 0.000 description 1
- 208000021155 Paediatric autoimmune neuropsychiatric disorders associated with streptococcal infection Diseases 0.000 description 1
- 240000000220 Panda oleosa Species 0.000 description 1
- 235000016496 Panda oleosa Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The application discloses a three-dimensional visualization method for multi-type real-time mapping of abnormal data of chemical detection, which comprises the following steps: s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model; s2, establishing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color of the minimum intensity value of the chemical element, and the rightmost side of the palette image represents the color of the maximum intensity value of the chemical element; s3, creating a fragment shader IN the GPU, wherein the fragment shader uses a sampling value of a model UV to an image IMG as W, the sampling value is a chemical element intensity value, and the sampling value W is used as an input value IN of the fragment shader; and S4, rendering the imaging model by using a fragment shader to obtain a three-dimensional visualized drawing piece after color mapping. The application can adjust the mapping mode of the color in real time, and can see the adjustment result in a short time, thereby being convenient and quick.
Description
Technical Field
The application relates to the field of geographic data image display, in particular to a three-dimensional visualization method and a system for multi-type real-time mapping of abnormal data by chemical detection.
Background
During the processing of geological data, most of the data eventually forms a map. For example, the chemical exploration data processing work is the processing, analysis and interpretation work of geochemical data in geochemical exploration and geochemical research in order to extract information having geological prospecting significance and other useful information. The process typically uses statistical tools to screen the collected sample data, determine which sample points are "abnormal", and define the scope of the abnormality to form an abnormality map. The areas with anomalies are areas with potential for prospecting. Also similar to this are remote sensing, geophysical prospecting exception handling and mapping.
The geological data map output based on the traditional method is fixed, and a user usually judges the ore-forming perspective of the region in the map through the color and the area of the region. However, the color and area of the area are affected by the parameters selected in advance, if the image parameters are not properly selected, the user may misjudge the range and form of the mine-forming remote scenic spot, and miss important mine-searching information.
Disclosure of Invention
The application mainly aims to provide a three-dimensional visualization method and a system for multi-type real-time mapping of abnormal data of visualization, wherein the mapping mode of the colors of the three-dimensional visualization map can be adjusted in real time by a user, and the adjustment result can be seen in a short time.
The technical scheme adopted by the application is as follows:
the three-dimensional visualization method for multi-type real-time mapping of the abnormal data of the chemical detection comprises the following steps:
s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
s2, establishing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color of the minimum intensity value of the chemical element, and the rightmost side of the palette image represents the color of the maximum intensity value of the chemical element;
s3, creating a fragment shader IN the GPU, wherein the fragment shader uses a sampling value of a model UV to an image IMG as W, the sampling value is a chemical element intensity value, the sampling value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is as follows:
1) Is provided withCalculating u=f (IN; MIN, MAX), the MIN, MAX being determined from the value entered by the user;
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and S4, rendering the imaging model by using a fragment shader to obtain a three-dimensional visualized drawing piece after color mapping.
By adopting the technical scheme, the user dynamically adjusts MIN and MAX according to the requirements, and the three-dimensional visual map after the re-color mapping is obtained in real time.
With the above technical solution, the number of channels of the image IMG is at least 1, and the data type of the channels should be floating point number.
The application also provides a three-dimensional visualization method for multi-type real-time mapping of the abnormal data, which comprises the following steps:
s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
s2, carrying out ascending order on original geochemical element data according to element intensity values to obtain an array A, and counting the number of elements as C, a minimum value VMin and a maximum value Vmax;
s3, establishing an image LUT_1 with the height of at least 1 pixel, wherein the image LUT_1 is used for storing the percentile;
s4, traversing each element in the array A, and carrying out the following operation on each element:
a) Setting the element value as X, calculating the normalized element value
Xn＝(X-VMin)/(VMax-VMin)；
b) Setting the serial number of the element in the array as N, and calculating the normalized serial number nn=N/C, wherein the initial serial number of the element in the array is 1;
c) Setting the pixel with U coordinate of Xn in the image LUT_1 as an Nn;
s5, traversing each pixel in the LUT_1, and carrying out the following operation on each pixel:
a) Setting the U coordinate of the pixel as X and the value as P, and skipping if the P is more than or equal to 0; if the ratio is less than 0, the next step is carried out;
b) Finding the first pixel at the left side of the pixel, wherein the value of the first pixel is PMin, and the U coordinate is
UMin。
c) Finding the first pixel not less than 0 on the right side of the pixel, wherein the value is PMax, and the U coordinate is
UMax。
d) Assignment p=pmin+ (X-omin)/(UMax-omin) × (PMax-PMin);
s6, providing an image LUT_2 of a color mapping palette by a user, wherein the height of the image LUT_2 is at least 1 pixel; the leftmost side of the image LUT_2 represents the data minimum color, and the rightmost side represents the data maximum color;
s7, inputting a minimum value MIN and a maximum value MAX of a percentile expected to be mapped by a user;
s8, creating a three-dimensional model on which the final mapping depends according to the actual situation of the user, and binding the IMG as a mapping to the three-dimensional model;
s9, creating a fragment shader in the GPU, taking VMin and VMax as uniform variables to be transmitted into the shader, setting the sampled value of the shader on the IMG as W by using a model UV, and the main logic of the shader is as follows:
a) Calculating a normalized value wn= (W-VMin)/(VMax-VMin);
b) Using (Wn, 0) as UV, sampling lut_1 to obtain a corresponding normalized percentile Pt;
c) Is provided withCalculating u=f (Pt; MIN, MAX);
d) Using (U, 0) as UV coordinates, sampling the image LUT_2, and outputting the obtained color;
and S10, rendering the three-dimensional model by using a fragment shader to obtain a three-dimensional visualized drawing after color mapping.
By adopting the technical scheme, the user dynamically adjusts MIN and MAX according to the requirements, and the three-dimensional visual map after the re-color mapping is obtained in real time.
With the above technical solution, the pixel format in the image lut_1 is either Float16 or Float32, and the pixel value is initialized to-1.
The application also provides a three-dimensional visualization system for multi-type real-time mapping of the abnormal data, which comprises the following steps:
the interpolation module is used for carrying out interpolation gridding on certain original geochemical element data needing to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
the palette construction module is used for constructing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color corresponding to the minimum value MIN of the intensity of the chemical element, and the rightmost side of the palette image represents the color corresponding to the maximum value MAX of the intensity of the chemical element; the values of MIN and MAX are determined according to the values input by a user;
the shader constructing module is configured to create a fragment shader, where the fragment shader uses a model UV to sample a value of an image IMG as W, the sample value is a chemical element intensity value, the sample value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is:
1) Is provided withCalculating u=f (IN; MIN, MAX);
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and the rendering module is used for rendering the imaging model by using the fragment shader to obtain a three-dimensional visual drawing piece after color mapping.
With the above technical solution, the number of channels of the image IMG is at least 1, and the data type of the channels should be floating point number.
By adopting the technical scheme, the system further comprises a setting module, wherein the setting module is used for acquiring MIN and MAX values input by a user, and dynamically adjusting the MIN and MAX according to the user requirements so as to obtain the three-dimensional visual map after re-color mapping in real time.
The application also provides a computer storage medium, in which a computer program executable by a processor is stored, and the computer program executes the three-dimensional visualization method for multi-type real-time mapping of the abnormal data by the technical scheme.
The application has the beneficial effects that: according to the application, after interpolation gridding is carried out on traditional geological data, particularly chemical detection data, the traditional geological data, particularly chemical detection data, is displayed in a three-dimensional space by using the GPU according to a user-defined color mapping mode, and a user can adjust the color mapping mode in real time, so that an adjustment result can be seen in a short time, and the method is convenient and quick.
Furthermore, by the method, not only can a single element geochemical diagram be generated, but also an element combination geochemical diagram and the like can be generated, and multi-type three-dimensional display of the chemical detection abnormal data is realized.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the present application, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a three-dimensional visualization method for multi-type real-time mapping of abnormal data in the embodiment 1 of the present application;
FIG. 2 is a schematic flow chart of embodiment 2 of the present application;
FIG. 3 is a TIFF image output by step 1 of example 3 of the present application;
FIG. 4 is a palette image generated in step 2 of embodiment 3 of the present application;
FIG. 5 is a graph of effects of embodiment 3 of the present application after real-time rendering;
FIG. 6 is a rendering chart after adjusting parameters in real time according to embodiment 3 of the present application;
fig. 7 is a flow chart illustrating a three-dimensional visualization method for multi-type real-time mapping of abnormal data in the embodiment 4 of the present application.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application.
Example 1
As shown in fig. 1, the three-dimensional visualization method for multi-type real-time mapping of abnormal data in the embodiment comprises the following steps:
s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
s2, establishing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color of the minimum intensity value of the chemical element, and the rightmost side of the palette image represents the color of the maximum intensity value of the chemical element;
s3, creating a fragment shader IN the GPU, wherein the fragment shader uses a sampling value of a model UV to an image IMG as W, the sampling value is a chemical element intensity value, the sampling value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is as follows:
1) Is provided withCalculating u=f (IN; MIN, MAX), the MIN, MAX being determined from the value entered by the user;
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and S4, rendering the imaging model by using a fragment shader to obtain a three-dimensional visualized drawing piece after color mapping.
According to the embodiment, after interpolation meshing is carried out on traditional geological data, particularly chemical detection data, the traditional geological data, particularly chemical detection data, is displayed in a three-dimensional space by using a GPU according to a user-defined color mapping mode, and a user can adjust the mapping mode of colors in real time by inputting a minimum value (MIN) and a maximum value (MAX) of data expected to be mapped, namely the maximum value and the minimum value of the intensity value of a chemical element when the chemical element is normalized by the pixel shader, so that an adjustment result can be seen in a short time (within 1 second), and convenience and rapidness are realized.
Example 2
This embodiment is based on embodiment 1, as shown in fig. 2, which comprises the following steps:
1. and carrying out interpolation gridding on the data needing to be imaged to generate an image IMG. The number of channels of the image is at least 1, and the data type of the channels should be floating point numbers.
2. The user provides an image of the color mapping palette having a height of at least 1 pixel. The leftmost side of the image represents the data minimum color and the rightmost side represents the data maximum color.
3. The user inputs the data minimum value MIN and maximum value MAX desired to be mapped.
4. And creating a model on which the final mapping depends according to the actual situation of the user, and binding the IMG as a mapping to the model.
5. A fragment shader (or pixel shader) is created, with shader input values IN, whose main logic is as follows:
a. is provided withCalculate u=f (IN; MIN, MAX)
b. Using (U, 0) as UV, the palette image is sampled and the resulting color is output.
6. Rendering the model in the step 4 by using the fragment shader, so as to obtain the data after color mapping.
Therefore, the user can dynamically adjust MIN and MAX according to the requirements, and the mapped drawing can be obtained in real time. By adjusting the maximum value and the minimum value, a user can quickly find out a high-value area (abnormal area) in the graph in the interactive process, and know the distribution condition of data in the area according to the morphological change of the abnormal area so as to assist in finally selecting a proper range of the prospecting target area.
Example 3
This example is based on example 2, taking the copper element geochemical intensity map in the cinnological region as an example, and dynamic images are generated by using the method of this patent:
1. excel tables of cinnological localization probe data were read using Pandas, and then grid interpolation was performed on copper element intensity values using GDAL's gdal_grid. The output is a TIFF image, a single channel,
the flow 32 type, resolution 512 x 391. As shown in fig. 3, in order to observe the data generated in step 1 directly using the image viewing tool DJV, since the original data is overexposed, the exposure value ev= -10 is set at the time of observation.
2. Using OpenCV, a palette TIFF image (as shown in fig. 4) with deep blue on the far left and deep red on the far right is generated based on color_jet, three channels of RGB, uint32 type, and resolution 512×16.
3. A fragment shader is written according to the following codes by using GLSL language, wherein texture [0] is bound to the element intensity image generated in step 1, and texture [1] is bound to the palette image in step 2:
4. the steps of creating polygon planes, fragment shader variable bindings and the like are consistent with the traditional GLSL usage,
and will be omitted herein. The real-time rendered image is shown in fig. 5, where min=5, max= 8457.
5. When the user modifies Min and max, the color range of the output image changes in real time, as shown in fig. 6, where min=213 and max=1220.
Example 4
The inventive concept of this embodiment is similar to embodiment 1, except that the chemical element intensity values are converted into quantiles in advance.
The three-dimensional visualization method for multi-type real-time mapping of the abnormal data of the visual detection mainly comprises the following steps:
1. the raw geochemical element data is interpolated and gridded to produce an Image (IMG). The number of channels of the image is at least 1, and the data type of the channels should be floating point numbers. This step is to convert the geochemical element data into an image file by interpolation to enable the GPU to process the data as a map. )
2. And (3) carrying out ascending order sequencing on the original geochemical element data according to the element intensity values to obtain an array A, and counting the number of elements as C, a minimum value VMin and a maximum value VMax. This step prepares the data for the generation of the next image.
3. An image (LUT_1) with a height of at least 1 pixel is created, the image width W being arbitrary, but a higher map width will have a smoother color change. The image is used to store the percentile. The pixel format is either Float16 or Float32, and the pixel value is initialized to-1. The image is used for establishing a corresponding relation between the element value and the percentile, the U coordinate of the image is the element value normalized between VMin and VMax, and the pixel value is the corresponding percentile. The reason for initializing the image to-1 here is that the element values in array a do not necessarily correspond to all pixels in the upper lut_1 after normalization, resulting in a missing pixel, which we need to be able to identify. For example, in arrays 0, 10, 10, 0 corresponds to U coordinate 0, 10 corresponds to U coordinate 1, then a width 3 image will have only the leftmost pixel and the rightmost pixel filled, and the middle pixel will be skipped. Since we have initialized with-1, these empty pixels will be found out in the later step, which is then filled in using interpolation methods.
4. Traversing each element in array A, and performing the following operation for each element:
a. setting the element value as X, calculating normalized element value
Xn= (X-VMin)/(VMax-VMin). This step determines the image U coordinate to which the value corresponds.
b. Let the element's sequence number in array a be N (the element's starting sequence number in array a be 1), calculate normalized sequence number nn=n/C. This step determines the percentile to which the value corresponds.
c. Pixels having U coordinates Xn in lut_1 are each set to value Nn, thereby filling the pixels.
5. Traversing each pixel in lut_1, for each pixel, performing the following operations (mainly for finding the empty pixel and interpolating the fill):
a. setting the U coordinate of the pixel as X and the value as P, and skipping if the P is more than or equal to 0;
if it is less than 0, the next step is performed.
b. Finding the first pixel at the left side of the pixel, wherein the value of the first pixel is PMin, and the U coordinate is UMin.
c. Finding the first pixel not less than 0 on the right side of the pixel, wherein the value is PMax, and the U coordinate is UMax.
d. Assignment p=pmin+ (X-PMin)/(UMax-PMin) ×x-PMin. The step mainly carries out linear interpolation according to the effective values around the empty pixel.
6. The user provides an image (lut_2) of the color mapping palette having a height of at least 1 pixel. The leftmost side of the image represents the data minimum color and the rightmost side represents the data maximum color. The image (lut_2) is mainly used to color the percentile.
7. The user inputs the Minimum (MIN) and Maximum (MAX) percentile desired to be mapped.
8. And creating a model on which the final mapping depends according to the actual situation of the user, and binding the IMG as a mapping to the model.
9. Creating a fragment shader (or called a pixel shader), taking VMin and VMax as uniform variables to enter the shader, and setting the sampled value of IMG by the shader to be W by using a model UV, wherein the main logic of the shader is as follows:
a. the normalized value wn= (W-VMin)/(VMax-VMin) is calculated.
b. Lut_1 is sampled using (Wn, 0) as UV, obtaining the corresponding normalized percentile Pt. (conversion of element values of geochemical images into percentiles in real time by LUT_1.)
c. Is provided withCalculate u=f (Pt; MIN, MAX)
(clipping and normalizing the percentile based on user input so that it can be mapped to the palette image LUT_2.)
d. Using (U, 0) as UV, lut_2 is sampled and the resulting color is output.
10. Rendering the model in the step 4 by using the fragment shader, so as to obtain the data after color mapping.
Example 5
The system embodiment is mainly used for realizing method embodiments 1, 2 and 3, and the three-dimensional visualization system for multi-type real-time mapping of the abnormal data by the chemical detection of the embodiment comprises:
the interpolation module is used for carrying out interpolation gridding on certain original geochemical element data needing to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
the palette construction module is used for constructing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color corresponding to the minimum value MIN of the intensity of the chemical element, and the rightmost side of the palette image represents the color corresponding to the maximum value MAX of the intensity of the chemical element; the values of MIN and MAX are determined according to the values input by a user;
the shader constructing module is configured to create a fragment shader, where the fragment shader uses a model UV to sample a value of an image IMG as W, the sample value is a chemical element intensity value, the sample value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is:
1) Is provided withCalculating u=f (IN; MIN, MAX);
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and the rendering module is used for rendering the imaging model by using the fragment shader to obtain a three-dimensional visual drawing piece after color mapping.
Example 6
The system embodiment is mainly used for realizing the method embodiment 4, and the three-dimensional visualization system for multi-type real-time mapping of the abnormal data by the embodiment comprises the following components:
the interpolation module is used for carrying out interpolation gridding on certain original geochemical element data needing to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
the percentile image creation module is used for creating an image LUT_1 with the height of at least 1 pixel, wherein the image LUT_1 is used for storing the percentile; and is used to traverse each pixel in lut_1, for each pixel:
a) Setting the U coordinate of the pixel as X and the value as P, and skipping if the P is more than or equal to 0; if the ratio is less than 0, the next step is carried out;
b) Finding the first pixel at the left side of the pixel, wherein the value of the first pixel is PMin, and the U coordinate is
UMin。
c) Finding the first pixel not less than 0 on the right side of the pixel, wherein the value is PMax, and the U coordinate is
UMax。
d) Assignment p=pmin+ (X-omin)/(UMax-omin) × (PMax-PMin);
the array setting module is used for carrying out ascending order on the original geochemical element data according to the element intensity values to obtain an array A, and counting the number of elements as C, a minimum value VMin and a maximum value Vmax; and is used to traverse each element in array a, for each element:
a) Assuming the element value as X, calculating normalized element value Xn= (X-VMin)/(VMax-VMin);
b) Setting the serial number of the element in the array as N, calculating the normalized serial number nn=N/C, wherein
The element initial sequence number in the array is 1;
c) Setting the pixel with U coordinate of Xn in the image LUT_1 as an Nn;
a palette creation module for providing an image lut_2 of a color mapping palette having a height of at least 1 pixel; the leftmost side of the image LUT_2 represents the data minimum color, and the rightmost side represents the data maximum color;
the setting module is used for inputting a minimum MIN and a maximum MAX of the percentile expected to be mapped through a user;
the three-dimensional model creation module is used for creating a three-dimensional model on which the final mapping depends according to the actual situation of a user, and binding the IMG as a map to the three-dimensional model;
the fragment shader creation module is used for creating a fragment shader in the GPU, taking VMin and VMax as uniform variables to be transmitted into the shader, setting the sampled value of the shader using a model UV to IMG as W, and the main logic of the shader is as follows:
a) Calculating a normalized value wn= (W-VMin)/(VMax-VMin);
b) Using (Wn, 0) as UV, sampling lut_1 to obtain a corresponding normalized percentile Pt;
c) Is provided withCalculating u=f (Pt; MIN, MAX);
d) Using (U, 0) as UV coordinates, sampling the image LUT_2, and outputting the obtained color;
and the rendering module is used for rendering the three-dimensional model by using the fragment shader to obtain a three-dimensional visualized drawing piece after color mapping.
Example 7
The present application also provides a computer readable storage medium such as a flash memory, a hard disk, a multimedia card, a card memory (e.g., SD or DX memory, etc.), 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 memory, a magnetic disk, an optical disk, a server, an App application store, etc., on which a computer program is stored that when executed by a processor performs a corresponding function. The computer readable storage medium of the present embodiment, when executed by a processor, implements the three-dimensional visualization method for multi-type real-time mapping of the abnormal data by the above-described method embodiments.
It should be noted that each step/component described in the present application may be split into more steps/components, or two or more steps/components or part of operations of the steps/components may be combined into new steps/components, according to the implementation needs, to achieve the object of the present application.
The sequence numbers of the steps in the above embodiments do not mean the execution sequence, and the execution sequence of the processes should be determined according to the functions and internal logic, and should not limit the implementation process of the embodiments of the present application.
It will be understood that modifications and variations will be apparent to those skilled in the art from the foregoing description, and it is intended that all such modifications and variations be included within the scope of the following claims.
Claims (10)
1. A three-dimensional visualization method for multi-type real-time mapping of abnormal data is characterized by comprising the following steps:
s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
s2, establishing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color of the minimum intensity value of the chemical element, and the rightmost side of the palette image represents the color of the maximum intensity value of the chemical element;
s3, creating a fragment shader IN the GPU, wherein the fragment shader uses a sampling value of a model UV to an image IMG as W, the sampling value is a chemical element intensity value, the sampling value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is as follows:
for determining from the value entered by the user;
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and S4, rendering the imaging model by using a fragment shader to obtain a three-dimensional visualized drawing piece after color mapping.
2. The three-dimensional visualization method for multi-type real-time mapping of abnormal data in visualization according to claim 1, wherein the user dynamically adjusts MIN and MAX according to the requirements to obtain the three-dimensional visualization map after re-color mapping in real time.
3. The method for three-dimensional visualization of multi-type real-time mapping of abnormal data of claim 1, wherein the number of channels of the image IMG is at least 1, and the data type of the channels is a floating point number.
4. A three-dimensional visualization method for multi-type real-time mapping of abnormal data is characterized by comprising the following steps:
s1, carrying out interpolation gridding on certain original geochemical element data to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
s2, carrying out ascending order on original geochemical element data according to element intensity values to obtain an array A, and counting the number of elements as C, a minimum value VMin and a maximum value Vmax;
s3, establishing an image LUT_1 with the height of at least 1 pixel, wherein the image LUT_1 is used for storing the percentile;
s4, traversing each element in the array A, and carrying out the following operation on each element:
a. assuming the element value as X, calculating normalized element value Xn= (X-VMin)/(VMax-VMin);
b. setting the serial number of the element in the array as N, calculating the normalized serial number nn=N/C, wherein
The element initial sequence number in the array is 1;
c. setting the pixel with U coordinate of Xn in the image LUT_1 as an Nn;
s5, traversing each pixel in the LUT_1, and carrying out the following operation on each pixel:
a) Setting the U coordinate of the pixel as X and the value as P, and skipping if the P is more than or equal to 0; if the ratio is less than 0, the next step is carried out;
b) Finding the first pixel at the left side of the pixel, wherein the value of the first pixel is PMin, and the U coordinate is UMin.
c) Finding the first pixel not less than 0 on the right side of the pixel, wherein the value is PMax, and the U coordinate is UMax.
d) Assignment p=pmin+ (X-omin)/(UMax-omin) × (PMax-PMin);
s6, providing an image LUT_2 of a color mapping palette by a user, wherein the height of the image LUT_2 is at least 1 pixel; the leftmost side of the image LUT_2 represents the data minimum color, and the rightmost side represents the data maximum color;
s7, inputting a minimum value MIN and a maximum value MAX of a percentile expected to be mapped by a user;
s8, creating a three-dimensional model on which the final mapping depends according to the actual situation of the user, and binding the IMG as a mapping to the three-dimensional model;
s9, creating a fragment shader in the GPU, taking VMin and VMax as uniform variables to be transmitted into the shader, setting the sampled value of the shader on the IMG as W by using a model UV, and the main logic of the shader is as follows:
a) Calculating a normalized value wn= (W-VMin)/(VMax-VMin);
b) Using (Wn, 0) as UV, sampling lut_1 to obtain a corresponding normalized percentile Pt;
c) Is provided withCalculating u=f (Pt; MIN, MAX);
d) Using (U, 0) as UV coordinates, sampling the image LUT_2, and outputting the obtained color;
and S10, rendering the three-dimensional model by using a fragment shader to obtain a three-dimensional visualized drawing after color mapping.
5. The three-dimensional visualization method for multi-type real-time mapping of abnormal data in visualization according to claim 4, wherein the user dynamically adjusts MIN and MAX according to the requirements to obtain the three-dimensional visualization map after re-color mapping in real time.
6. The method for three-dimensional visualization of multi-type real-time mapping of anomaly data of claim 4, wherein the pixel format in the image LUT_1 is Float16 or Float32, and the pixel value is initialized to-1.
7. A three-dimensional visualization system for multi-type real-time mapping of abnormal data of a chemical detection, comprising:
the interpolation module is used for carrying out interpolation gridding on certain original geochemical element data needing to be imaged to generate an image IMG, and binding the image IMG as a map to a preset imaging model;
the palette construction module is used for constructing a color mapping palette image with the height of at least 1 pixel, wherein the leftmost side of the palette image represents the color corresponding to the minimum value MIN of the intensity of the chemical element, and the rightmost side of the palette image represents the color corresponding to the maximum value MAX of the intensity of the chemical element; the values of MIN and MAX are determined according to the values input by a user;
the shader constructing module is configured to create a fragment shader, where the fragment shader uses a model UV to sample a value of an image IMG as W, the sample value is a chemical element intensity value, the sample value W is used as an input value IN of the fragment shader, and execution logic of the fragment shader is:
1) Is provided withCalculating u=f (IN; MIN, MAX);
2) Sampling the color mapping palette image using (U, 0) as a UV coordinate value, and outputting the obtained color;
and the rendering module is used for rendering the imaging model by using the fragment shader to obtain a three-dimensional visual drawing piece after color mapping.
8. The system of claim 7, wherein the number of channels of the image IMG is at least 1, and the data type of the channels is a floating point number.
9. The system of claim 7, further comprising a setting module for obtaining MIN and MAX values input by a user, and dynamically adjusting the MIN and MAX according to the user's needs to obtain the three-dimensional visual map after the re-color mapping in real time.
10. A computer storage medium, in which a computer program executable by a processor is stored, the computer program performing the method for three-dimensional visualization of multi-type real-time mapping of chemical detection anomaly data according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310802525.6A CN116912438A (en) | 2023-06-30 | 2023-06-30 | Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310802525.6A CN116912438A (en) | 2023-06-30 | 2023-06-30 | Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116912438A true CN116912438A (en) | 2023-10-20 |
Family
ID=88350346
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310802525.6A Pending CN116912438A (en) | 2023-06-30 | 2023-06-30 | Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116912438A (en) |
-
2023
- 2023-06-30 CN CN202310802525.6A patent/CN116912438A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8605951B2 (en) | Method for converting a digital image into a multi-dimensional geo-referenced data structure | |
| US7327363B2 (en) | Image processing apparatus, and computer program | |
| US7245300B2 (en) | Architecture for real-time texture look-up's for volume rendering | |
| US20100207936A1 (en) | Fusion of a 2d electro-optical image and 3d point cloud data for scene interpretation and registration performance assessment | |
| US8477135B2 (en) | Method and apparatus for volume rendering using depth weighted colorization | |
| JP4820004B2 (en) | Method and system for filtering image data to obtain samples mapped to pixel subcomponents of a display device | |
| CN109887073B (en) | Method and device for building three-dimensional digital model of rock core | |
| Marschallinger | Three-dimensional reconstruction and visualization of geological materials with IDL—examples and source code | |
| CN105096385A (en) | Method for three-dimensional display of two-dimensional earthquake profile | |
| WO2014072810A1 (en) | Method for processing an image | |
| US20090303236A1 (en) | Method and system for explicit control of lighting type in direct volume rendering | |
| CN110428504B (en) | Text image synthesis method, apparatus, computer device and storage medium | |
| US6400370B1 (en) | Stochastic sampling with constant density in object space for anisotropic texture mapping | |
| US20070188492A1 (en) | Architecture for real-time texture look-up's for volume rendering | |
| Görtler et al. | Stippling of 2D scalar fields | |
| US6967653B2 (en) | Apparatus and method for semi-automatic classification of volume data | |
| CN116912438A (en) | Three-dimensional visualization method and system for multi-type real-time mapping of abnormal data of chemical exploration | |
| CN117409157A (en) | Online sectioning and graph cutting method for three-dimensional geological model | |
| CN116758206A (en) | Vector data fusion rendering method and device, computer equipment and storage medium | |
| Gavrilescu et al. | Custom interface elements for improved parameter control in volume rendering | |
| US5821942A (en) | Ray tracing through an ordered array | |
| CN113126944A (en) | Depth map display method, display device, electronic device, and storage medium | |
| CN113593028A (en) | Three-dimensional digital earth construction method for avionic display control | |
| WO2019236064A1 (en) | Automatically updating a georeferencing graphical user interface for navigation line adjustments | |
| Kleefeld et al. | Processing Multispectral Images via Mathematical Morphology. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |