AU2015229467B9 - Blend modes for mineralogy images - Google Patents

Blend modes for mineralogy images Download PDF

Info

Publication number
AU2015229467B9
AU2015229467B9 AU2015229467A AU2015229467A AU2015229467B9 AU 2015229467 B9 AU2015229467 B9 AU 2015229467B9 AU 2015229467 A AU2015229467 A AU 2015229467A AU 2015229467 A AU2015229467 A AU 2015229467A AU 2015229467 B9 AU2015229467 B9 AU 2015229467B9
Authority
AU
Australia
Prior art keywords
image
value
sample
determining
top layer
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.)
Active
Application number
AU2015229467A
Other versions
AU2015229467B2 (en
AU2015229467A1 (en
Inventor
Eric GOERGEN
David HABERLAH
Garth Howell
Michael James Owen
Leonardo SALAZAR
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
FEI Co
Original Assignee
FEI Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Priority claimed from US14/204,485 external-priority patent/US9183656B2/en
Priority claimed from EP14173457.4A external-priority patent/EP2960865A1/en
Application filed by FEI Co filed Critical FEI Co
Publication of AU2015229467A1 publication Critical patent/AU2015229467A1/en
Application granted granted Critical
Publication of AU2015229467B2 publication Critical patent/AU2015229467B2/en
Publication of AU2015229467B9 publication Critical patent/AU2015229467B9/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/22Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
    • G01N23/225Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion
    • G01N23/2251Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
    • G01N23/2252Measuring emitted X-rays, e.g. electron probe microanalysis [EPMA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/50Lighting effects
    • G06T15/503Blending, e.g. for anti-aliasing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/401Imaging image processing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/405Imaging mapping of a material property
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/418Imaging electron microscope
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/606Specific applications or type of materials texture
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/616Specific applications or type of materials earth materials

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

Optimized blending mode for mineralogy images. A luminosity value is determined for a pixel in a base layer or top layer mineralogy image. An image weighting value is determined from the luminosity value and an optional mixing parameter. A multiply value is determined by multiplying the base and top layer pixel values. An overlay value is determined from twice the multiply value if the value of one of the base layer or top layer pixel values is over a threshold, otherwise it is determined by inverting twice the product of the inverted top layer pixel value with the inverted base layer pixel value. A blended image pixel value is determined by adding the multiply value weighted with the image weighting value and the overlay value weighted with the inverted image weighting value.

Description

Blend Modes for Mineralogy Images
Technical Field of the Invention [001] The present invention relates to the generation and display of mineralogy images, and in particular to a method for automatically blending mineralogy images of a sample, including textural, elemental and spatial mineral identification images.
Background of the Invention [002] To properly understand and interpret the mineralogy of a sample, knowledge of both the textural information of the sample and the spatial distribution of minerals in the sample is required.
[003] The textural information of a geological sample refers to the size, shape and arrangement of grains or crystals within the sample, the homogeneity of the sample, its degree of isotropy, and the distribution of elements within the minerals The textural information can be measured, for example, by scanning the sample with a high energy beam of charged particles, and measuring at each scanned position the intensity of back-scattered electrons (“BSE’s”) that are received at an electron detector. The measured BSE intensities are indicators of the average atomic number of elements in the sample as a function of scan position. A plot of the intensity versus scan position can be presented as a grey scale image, where each pixel in the image is a shade of grey that corresponds to the measured BSE intensity at the scanned position. The image so obtained will be darkest where the measured BSE intensity is least, and lightest where the measured BSE intensity is greatest. Of course, this color scheme can be inverted so that the image will be lightest where the measured BSE intensity is least, and darkest where the measured BSE intensity is greatest.
[004] In addition to BSE intensity, other forms of textural information can include representation of the rate of x-ray production at each pixel. This is measured by scanning the sample with a high energy beam of charged particles, and measuring the intensity of x-rays produced by the sample and received by an x-ray detector. The x-rays are produced when the charged particles dislodge electrons in the sample, and electrons from higher energy shells fall into the energy shells vacated by the dislodged electrons. The rate at which the x-rays are produced is roughly proportional to the volume and atomic number of the material in the sample. A plot of the x-ray generation rate versus scan position can be presented as a grey scale image, where each pixel in the image is a shade of grey corresponding to the measured rate of x-ray production at each scan position. The image so obtained will be darkest where the measured x-ray production rate is lowest, and brightest where the measured x-ray production rate is highest.
[005] The spatial distribution of minerals in the sample is a plot of the minerals identified in the sample as a function of sample position. It can be generated by scanning the sample with a high energy beam, and measuring the energy distribution of x-rays emitted from the sample as a function of scan position. On a per pixel basis, these energy distributions can be fitted and/or compared to a catalog of energy distributions obtained from pure elements or pure minerals in order to identify the elements and/or minerals in the sample at each scanned position. Techniques for identifying minerals based on a catalog of elemental x-ray spectra are disclosed, for example, in U.S. Patent Application No. 12/866,697, filed on Feb. 6, 2009, which is herein incorporated by reference in its entirety. Techniques for identifying minerals based on a catalog of mineral x-ray spectra are disclosed, for example, in U.S. Patent Application No. 14/073,523, filed on Nov. 6, 2013, which is herein incorporated by reference in its entirety. Different colors can be assigned to different minerals in the catalog, and an image of the spatial mineral distribution in the sample can be generated by plotting the colors of identified minerals as a function of scanned position.
[006] Since proper understanding of the mineralogy of a sample requires knowledge of both the textural information of the sample and the spatial distribution of minerals in the sample, methods are needed to combine the textural and spatial mineral distribution images into a single composite image that reveals both the spatial mineral distribution and textural information about the sample.
Summary of the Invention [007] An object of the invention is to provide an automated mechanism to allow data from multiple detectors to be simultaneously visualized without losing information from either data source.
[008] In one aspect, the invention features a computer implemented method for blending mineralogy images, including a base layer mineralogy image and a top layer mineralogy image. For each pixel in the base or top layer mineralogy images a base layer value is determined from the pixel in the base layer mineralogy image, and a top layer value is determined from the pixel in the top layer mineralogy image. A luminosity value is determined from one of the base layer or top layer values. A multiply value is determined by multiplying the base and top layer values. An overlay value is determined by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value. An image blending weight is determined from the luminosity value. A blended image value is determined as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight.
[009] Implementations of the invention may include one or more of the following features. The base layer mineralogy image can be a textural mineralogy image and the top layer mineralogy image can be a spatial mineral distribution image. The textural mineralogy image can be an image of the intensity of back-scattered electrons or an image of the intensity of generated x-rays. The image blending weight can be determined by multiplying the luminosity value by a mixing parameter. A blended image, where each pixel of the blended image has a blended image value determined from the base layer value and top layer value, can be displayed on a display device.
[0010] In another aspect, the invention features a computer implemented method for blending mineralogy images of a sample, including an elemental distribution image of the sample, and a mineral distribution image of the sample. A blended image of the sample can be generated by blending the elemental and mineral distribution images of the sample.
[0011] Implementations of the invention may include one or more of the following features. Each pixel in the blended image of the sample can be generated from each pixel in the elemental and mineral distribution images. A base layer value can be determined from the pixel in the elemental distribution image. A top layer value can be determined from the pixel in the mineral distribution image. A luminosity value can be determined from one of the base layer or top layer values. A multiply value can be determined by multiplying the base and top layer values. An overlay value can be determined by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value. An image blending weight can be determined from the luminosity value. A blended image value can be determined as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight. The image blending weight can be determined by multiplying the luminosity value by a mixing parameter.
[0012] The foregoing has outlined rather broadly the features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of the invention will be described hereinafter. It should be appreciated by those skilled in the art that the conception and specific embodiments disclosed may be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims.
[0013] Implementations of the invention may also include one or more of the following features. Each pixel in the blended image of the sample can be generated by determining a top layer value from the pixel in the mineral distribution image, determining an opacity value from the pixel in the elemental distribution image; and alpha blending the top layer value with a black pixel value using the opacity value. The pixel in the blended image can be rendered black when the corresponding pixel in the elemental distribution image has a value that is less than an elemental threshold.
Brief Description of the Drawings [0014] The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by the Office upon request and payment of the necessary fee. For a more thorough understanding of the present invention, and advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which: [0015] FIG. 1 is a flow chart depicting a method for blending textural and spatial mineral distribution images.
[0016] FIG. 2A is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the multiply blend mode.
[0017] FIG. 2B is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the multiply blend mode.
[0018] FIG. 2C is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the overlay blend mode.
[0019] FIG. 2B is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the special blend mode.
[0020] FIG. 3 is a flow chart depicting a method for blending an elemental distribution image and a spatial mineral distribution image.
[0021] FIG. 4 is an illustration of a blended image made by blending an elemental distribution image with a spatial mineral distribution image.
[0022] FIG. 5 is an illustration of a mineral identification and analysis system. Detailed Description of Preferred Embodiments [0023] Proper understanding of the mineralogy of a sample requires knowledge of both the texture of the sample and the spatial distribution of minerals in the sample. Methods are disclosed for forming a blended image by combining textural and spatial mineral distribution images of a sample to facilitate an understanding of its mineralogy. In one embodiment, the two images can be blended together using known image blending techniques, such as by compositing the two images as layers in a blended image. The degree of blending can be controlled by assigning a global opacity value to the alpha channel of the spatial mineral distribution image. The textural and spatial mineral distribution images can then be blended using a conventional alpha blending algorithm. This approach, however, tends to wash out the colors in the spatial mineral distribution layer.
[0024] Other image blending techniques from the image processing arts can be used to combine the textural and spatial mineral distribution images. For example, the color values of pixels in the textural and spatial mineral distribution images can be added together, subtracted, multiplied or divided by one another to obtain the color values of corresponding pixels in the blended image. Since each of these image blending techniques is globally applied (i.e., identically applied to each pixel in the blended image), each technique produces a global effect on the blended image. For example, multiplying images creates a blended image that is everywhere darker than either of the combined images, while dividing images creates a blended image that is everywhere lighter than either of the combined images.
[0025] Some image blending modes address these global compositing issues by thresholding the image compositing function. That is, different image compositing functions are used depending on whether the color value in one of the layers to be blended exceeds a threshold. Examples of such blending modes include the overlay, hard light, and soft light blending. However, depending on the threshold selected, portions of the image far from the threshold may appear too dark or too light.
[0026] To optimize the information presented in a blended image consisting of textural and spatial mineral distribution images from a sample, a custom blending method has been developed. This method combines elements from conventional multiply, screen and overlay blending modes, and mixes these modes on a per pixel basis as a function of the luminosity of the pixel color in the textural or spatial mineral distribution images.
[0027] FIG. 1 is a flow chart depicting a method for blending textural and spatial mineral distribution images. The method begins by receiving a textural image and a spatial mineral distribution image for a sample (100). In one implementation, an optional mixing parameter k that ranges in value between 0 and 1 is also received (110). In some implementations, the mixing parameter k can be received from a computer memory, where it is set to a default value such as 0.6. In other implementations, the mixing parameter can be received as user input, e.g., from a text box or as a radio button selection in a graphical user interface presented to the user. Setting the mixing parameter k to a default value of lis equivalent to having no mixing parameter at all, as further explained below.
[0028] The method proceeds by looping through all of the pixels in the textural and spatial mineral distribution images (120). If the images have different resolutions, one or both images can be renormalized so they have the same resolution or number of pixels. For each pixel, the luminance of the color in either the textural or the spatial mineral distribution image is determined (130). For example, in one implementation, the luminance of the textural image is determined. If this distribution is a grey scale image (e.g., representing the intensity of BSE’s), the luminance is simply the grey value for that pixel. In another implementation, the luminance of the spatial mineral distribution image is determined. If the pixel color in the spatial mineral distribution image is stored in the RGB or red, green, blue color space, the luminance Y can be determined as:
If the pixel color is stored in other color spaces, well known formulas can be used to determine the luminance from the pixel’s color values in those color spaces. By normalization, the luminance lies within a range of 0 to 1.
[0029] Next, pixel values for the conventional multiply, screen and overlay blending modes are obtained from the pixel values of the textural and spatial mineral distribution images. The multiply blend value is simply the product of the pixel values from these two images: Μ^(1) = Τ·^(ΐ) * SM(Eq. 2) where My(l), 7)//), and SMDy(l) are respectively the multiply blend value, the textural image value, and the spatial mineral distribution image value for the Ith color component of pixel (ij). The pixel value in the multiply blend will be darker than the pixel values in either the textural or spatial mineral distribution images. This mode is therefore useful for blending high intensity areas of the textural and spatial mineral distribution images; however, low intensity areas often come out too dark to provide useful information about the mineralogy of the sample. The screen blend pixel value SCRy(l) is obtained by inverting the product of the inverted pixel values from the textural and spatial mineral distribution images. It can be computed as:
The pixel value in the screen blend will be lighter than the pixel values in either the textural or spatial mineral distribution images. This mode is therefore useful for low intensity areas of the textural and spatial mineral distribution images; however, high intensity areas often come out too light to provide useful information about the mineralogy of the sample. Finally, the overlay blend value Oy(l) can be computed as:
The overlay blend mode attempts to address the shortcomings in the multiply and screen blend modes when dealing with very bright or dark compositing layers by selecting the multiply mode if the pixel value is less than a threshold and the screen mode if the pixel value is greater than the threshold.
[0030] Next, a special blend value is calculated for the pixel in the blended image from the multiply and overlay pixel values, the mixing parameter k, and the luminosity of the pixel color in either the textural or the spatial mineral distribution image. The special blend value is computed as a weighted average of the multiply and overlay pixel values, where the weights are determined by the product r = k * Y, where k is the optional mixing parameter and Y is the luminosity of the pixel color in the textural or spatial mineral distribution image. Thus, the special blend SBy(l) value of the pixel is given by:
Since both the mixing parameter k and the pixel intensity Y are normalized to lie in the range of 0 to 1, the weight r also lies in the range of 0 to 1. Note that when the special mixing parameter k = 1, the weight r reduces to the luminosity Y as if there were no mixing parameter k.
[0031] As a weighted average, the special blend provides for continuous blending of the textural and spatial mineral distribution images, with no thresholding. It produces a blended image that is optimized for both low and high luminosity areas in the spatial mineral distribution image, when blended with both low and high luminosity areas in the textural image. High intensity colors in the spatial mineral distribution image (e.g., bright reds and yellows) retain their chrominance in the blended image and are easily discernible, even when blended with high intensity areas of the textural image. Moreover, high intensity colors in the spatial mineral distribution image retain their luminance in the blended image, even when blended with low intensity areas of the textural image.
[0032] FIG. 2A is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the multiply blend mode. The blended textural and spatial mineral distribution images were made for a sample containing the minerals pyrite (bright yellow) and rutile (red). As seen in the illustration, while the multiply blend image preservers the colors of the pyrite and rutile minerals, it is everywhere else too dark, and fails to preserve useful information from the textural and the spatial mineral distribution images.
[0033] FIG. 2B is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the screen blend mode. As seen in the illustration, the screen blend image is lighter than the multiply blend image shown in FIG. 2A. As a result, more information from both the textural and spatial mineral distribution images is revealed. Nonetheless, the high luminosity colors for pyrite and rutile are significantly faded, resulting in a loss of information from the spatial mineral distribution image.
[0034] FIG. 2C is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the overlay blend mode. As seen in the illustration, the overlay blend image is lighter than the multiply blend image shown in FIG. 2A, but darker than the screen blend image shown in FIG. 2B. As a result, more information from both the textural and spatial mineral distribution images is revealed in the overlay blended image. Nonetheless, the high luminosity colors for pyrite and rutile remain faded in the blended image, resulting in a loss of information from the spatial mineral distribution image.
[0035] FIG. 2D is an illustration of a blended image made by blending a textural image with a spatial mineral distribution image using the special blend mode. As seen in the illustration, the special blend image is lighter than the multiply blend image shown in FIG. 2A, but darker than the overlay blend image shown in FIG. 2C. As a result, more information from both the textural and spatial mineral distribution images is revealed. Moreover, the high luminosity colors for pyrite and rutile continue to stand out in the special blend image, preserving information about these minerals in the sample recorded in the spatial mineral distribution image [0036] The methods disclosed herein can be used to blend any type of textural and spatial mineral distribution information. In some embodiments, the textural image can be a per pixel record of the total number of x-rays emitted by the sample at the pixel scan position. In this embodiment, the blended image can provide a user with information regarding confidence in the mineral identification for the sample. For example, if a sample’s spatial mineral distribution image contained a large monotone red area (indicating the identification of rutile in some area of the sample), and the sample’s textural image (composed of per-pixel x-ray counts) showed large counts in the center of the red area but smaller counts toward its edges, the corresponding area in the blended image would be different shades of red corresponding to the different x-ray pixel counts in the textural image. Where the x-ray count is larger (e.g., in the center), the red in the blended image would be a brighter shade than in those areas where the x-ray count is smaller (e.g., along the edges). Since the likelihood of correctly identifying a mineral is proportional to the x-ray count, areas of the blended image having a brighter shade of a pure mineral color (e.g., red, blue) will be areas of higher mineral identification confidence, while areas of the blended image having a darker shade of a pure mineral color will be areas of lower mineral identification confidence.
[0037] FIG. 3 is a flow chart depicting, in yet another embodiment of the invention, a method for blending elemental and mineral distribution images. Both the elemental distribution image and the mineral distribution image can be obtained by analyzing, on a per pixel basis, the energy distribution of x-rays emitted from a sample when bombarded with charged particles. As shown in FIG. 3, the process begins by raster scanning the sample (300). Next, for each pixel in the image (310), an x-ray emission spectrum is obtained (320). This sample x-ray emission spectrum is then analyzed to determine both the elemental (330) and the mineral (340) composition of the sample at the pixel scan position. This process is repeated (350), for all pixels or scan positions for the sample.
[0038] As disclosed in U.S. Patent Application No. 12/866,697, filed on Feb. 6, 2009, which is herein incorporated by reference, the elemental distribution image for the sample can be obtained by fitting, on a per pixel basis, the sample’s x-ray emission spectrum with elemental x-ray emission spectra to determine the number, kind, and percentage of elements present in the sample. The elemental composition so determined can then be used to determine the mineral composition and the mineral distribution of the sample. Moreover, for each element present in the sample, an elemental distribution image can be generated by assigning a grayscale value to each pixel in the image corresponding to the percentage of the element found in the sample at that pixel location or scan position.
[0039] As disclosed in U.S. Patent Application No. 14/073,523, filed on Nov. 6, 2013, which is herein incorporated by reference in its entirety, the spatial mineral distribution image for the sample can also be obtained by fitting, on a per pixel basis, the sample’s x-ray emission spectra with mineral x-ray emission spectra to determine the number, kind, and percentage of minerals present in the sample. Each mineral determined to be present in the sample can be assigned a color, and when more than one mineral is present the colors can be blended together, with the amount of blend determined by the percentage of mineral present in the sample.
[0040] Referring again to FIG. 3, once the elemental and mineral distribution images are obtained, an elemental selection can be received. For example, the process can present a user with a list of elements and minerals identified in the sample, and with the percentages of each. The user can then select a particular element for forming a blended image. For example, the user can select the element by entering its name in a text box, or by selecting a radio button or check box associated with the element in a graphical user interface. Once selected, the elemental distribution image and the mineral distribution image can be combined.
[0041] In one embodiment, the elemental and mineral distribution images can be combined using the special blend mode discussed above in reference to FIG. 1. In other embodiments, the elemental and mineral distribution images can be combined using any conventional image blending function. For example, the grayscale value recorded in the elemental distribution image can be used to determine the opacity of the mineral distribution image. In this way, the blended image will only contain areas of the mineral distribution image that overlap with areas in the elemental distribution image. All other areas of the blended image will appear black. Moreover, the shades of colors assigned to minerals identified in the mineral distribution image will vary in the blended image depending on the percentage of the selected element found in the elemental distribution image. Identified mineral colors in the blended image will be brightest in those areas corresponding to elemental distribution image areas having a high percentage of the selected element, and darkest in those areas corresponding to elemental distribution image areas having a low percentage of the selected element.
[0042] In some embodiments, when an element selection is received in step 360, an elemental threshold can be optionally received. In these embodiments, pixels in the elemental distribution image will only be blended with pixels in the mineral distribution image if they are above the elemental threshold. For example, if the received elemental threshold is 12%, the blended image will only contain pixels in the mineral distribution image blended with pixels in the elemental distribution image that have a grayscale value greater than 0.12. Pixels in the elemental distribution image having a grayscale value less than 0.12 will appear black in the blended image. In this way, regions of interest, defined as regions having an elemental composition greater than a specified threshold value, can be easily defined and displayed.
[0043] FIG. 4 is an illustration of a blended image made by blending an elemental distribution image with a mineral distribution image. In FIG. 4, the selected element is magnesium (Mg), and the blended image shows areas of a sample where Mg is present in either dolomite (blue) or clay (green). As seen from the image, different shades of green indicate areas containing clays (identified in the mineral distribution image) having different percentages of Mg (determined from elemental distribution image).
[0044] FIG. 5 is an illustration of a mineral identification and analysis system 200. The mineral identification system 200 includes a scanning electron beam system 241, an x-ray detector 240, and a secondary or backscatter electron detector 243. An electron beam 232 emitted from a cathode 253 is accelerated toward an anode 254. Electron beam 232 is subsequently focused to a fine spot by means of a condensing lens 256 and an objective lens 258, and can be deflected across a sample 202 by means of a deflection coil 260 to perform a two-dimensional raster scan of the sample. The condensing lens 256, objective lens 258, and deflection coil 260 are supplied current by a power supply 245 operated under the control of a system controller 233. The sample 202 is preferably mounted on a movable X-Y stage 204 within a lower vacuum chamber 210. The vacuum chamber 210 is evacuated to high vacuum by a mechanical pumping system 269 and an ion pump 268 operated under the control of vacuum controller 232.
[0045] When the electron beam 232 strikes the sample 202, several forms of radiation are emitted, including backscattered electrons from the electron beam 232, secondary electrons produced by interactions between the electron beam 232 and the sample, and x-rays produced by interactions between the electron beam 232 and the sample that are characteristic of the elements in the sample. The back scattered electrons are detected by an electron detector 242, which outputs a signal indicative of the flux or intensity of back scattered electrons. This signal is received by processor 220. The x-rays emitted from the sample are detected by x-ray detector 240, which preferably outputs a signal indicative of the energy of the detected x-rays. To that end, x-ray detector 240 is preferably an energy dispersive detector such as a silicon drift detector. The output signal of x-ray detector 240 can be amplified and received by processor 220. The back-scattered or secondary electrons are detected by an electron detector, such as a [0046] Processor 220 can be programmed to store, for each scanned pixel, a count of the number of back scattered electrons detected, a count of the number of x-rays detected, and a histogram counting the number of x-rays detected in each of a plurality of energy bins over some energy range. Typically, the energy range is on the order of 0 - lOkeV, and is subdivided into energy bins of 10-20 eV, for a total of 500 to 1000 energy bins or channels per pixel.
[0047] System 200 also includes a display 244 for displaying the results of the textural and mineral analysis and other information by way of video circuit 242; a program memory 222 for storing executable computer program code to program the processor 220, and a data memory 223 /224 for storing data, such as per-pixel BSE counts, x-ray counts, and x-ray emission spectra recorded from sample 202, and a library of standardized elemental or mineral x-ray emission spectra. Program memory 222 can include computer storage media in the form of removable and/or non-removable, volatile and/or nonvolatile memory and can provide storage of computer-readable instructions, data structures, program modules and other data. Generally, the processor 220 is programmed by means of instructions stored in the various computer-readable storage media. Programs and operating systems are typically distributed, for example, on floppy disks or CD-ROMs. From there, they are installed or loaded into the secondary memory of a computer. At execution, they are loaded at least partially into the computer's primary electronic memory. The invention described herein includes these and other various types of computer-readable storage media when such media contain instructions or programs for implementing the steps described above in conjunction with a microprocessor or other data processor. The invention also includes the computer itself when programmed according to the methods and techniques described herein.
[0048] While the embodiment shown uses a scanning electron microscope to generate x-rays from sample 202, other embodiments could employ a transmission electron microscope or a scanning transmission electron microscope. An x-ray fluorescence system could also be used to generate x-rays from sample 202. In other embodiments, different forms of characteristic radiation emitted from the sample, such as gamma rays, may be detected.
[0049] Some embodiments of the invention provide a method of analyzing a sample, comprising: [0050] directing an electron beam toward multiple points on the surface of a sample; [0051] detecting x-rays emitted from the multiple points on the sample surface in response to the electron beam; [0052] forming a measured x-ray spectrum from the x-rays emitted from each of the multiple points on the surface; [0053] comparing the measured x-ray spectrum from each of the multiple points with reference x-ray spectra to determine the composition at each of the points on the surface; [0054] determining textural information at the multiple points on the surface; [0055] blending information from the multiple points on the surface to produce a blended sample image composed of pixels, the pixels corresponding to the multiple points on the sample, the color and brightness of each pixel in the blended sample image being determined by a combination of the composition and the textural information at the multiple points.
[0056] In some embodiments, the method further comprises displaying the blended sample image.
[0057] In some embodiments, determining textural information includes determining textural information from back scattered electrons produced by the impact of the electrons on the multiple points on the surface.
[0058] In some embodiments, determining textural information includes determining textural information from the intensity of x-ray produced by the impact of the electrons on the multiple points on the surface.
[0059] In some embodiments, comparing the measured x-ray spectrum with reference x-ray spectra to determine the composition of the points on the surface includes producing a top layer image composed of a top layer value for each of the multiple points and determining textural information includes producing a base layer image composed of base layer values, and these embodiments may further comprise: [0060] determining a luminosity value from one of the base layer or top layer values; [0061] determining a multiply value by multiplying the base and top layer values; [0062] determining an overlay value from the multiply value or from the top layer value and the bottom layer value; [0063] determining an image blending weight from the luminosity value; [0064] determining a blended image value from the multiply value and the overlay value; and [0065] displaying an image in which the pixels of the image correspond to the blended image values.
[0066] In some embodiments, determining an overlay value comprises determining an overlay value by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value; and [0067] determining a blended image value comprises determining a blended image value as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight.
[0068] Some embodiments of the invention provide a method of analyzing a sample, comprising: [0069] directing an electron beam toward multiple points on the surface of a sample; [0070] detecting x-rays emitted from the points on the sample surface in response to the electron beam; [0071] forming a measured x-ray spectrum from x-rays emitted from each of the points on the surface; [0072] comparing the measured x-ray spectra with reference x-ray spectra to determine the elemental composition of the points on the surface and to determine the mineral composition of the points on the surface; and [0073] generating a blended image corresponding to the multiple points on the sample, the blended image including pixels corresponding to the multiple points on the sample, the color and the brightness of the pixels determined in part by the elemental composition at the corresponding point on the sample and determined in part by the mineral composition at the corresponding points on the sample.
[0074] In some embodiments, comparing the measured x-ray spectra with reference x-ray spectra includes determining an elemental distribution map of the sample and determining a mineral distribution map of the sample, each of the maps composed of pixels; and [0075] generating a blended image corresponding to the multiple points on the sample includes blending the pixels of the elemental distribution map and the pixels of the mineral distribution map [0076] In some embodiments, the method further comprises receiving an elemental threshold and generating a blended image comprises rendering a pixel of the blended image black when a corresponding pixel of the elemental distribution image has a value that is less than the elemental threshold.
[0077] Some embodiments of the invention provide a computer implemented method for blending mineralogy images, comprising: [0078] receiving a base layer mineralogy image and a top layer mineralogy image; and [0079] sequentially, via a computer processor and for each pixel in the base or top layer mineralogy images: [0080] determining a base layer value from the pixel in the base layer mineralogy image; [0081] determining a top layer value from the pixel in the top layer mineralogy image; [0082] determining a luminosity value from one of the base layer or top layer values; [0083] determining a multiply value by multiplying the base and top layer values; [0084] determining an overlay value by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value; [0085] determining an image blending weight from the luminosity value; and [0086] determining a blended image value as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight.
[0087] In some embodiments, the base layer mineralogy image is a textural mineralogy image and the top layer mineralogy image is a spatial mineral distribution image.
[0088] In some embodiments the textural mineralogy image is an image of the intensity of back-scattered electrons.
[0089] In some embodiments the base layer mineralogy image is an image of the intensity of generated x-rays and the top layer mineralogy image is a spatial mineral distribution image.
[0090] In some embodiments the luminosity value is determined from the top layer value.
[0091] In some embodiments the overlay value is determined by doubling the multiply value if the top layer value is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value.
[0092] Some embodiments further comprise displaying the blended image on a display device.
[0093] In some embodiments the base layer mineralogy image is a spatial elemental distribution mineralogy image and the top layer mineralogy image is a spatial mineral distribution image.
[0094] Some embodiments further comprise: receiving a mixing parameter; and determining the image blending weight by multiplying the luminosity value by the mixing parameter.
[0095] Some embodiments of the invention provide a computer implemented method for blending mineralogy images of a sample, comprising: [0096] receiving an elemental distribution image of the sample; [0097] receiving a mineral distribution image of the sample; [0098] generating, via a computer processor, a blended image of the sample by blending the elemental and mineral distribution images of the sample; and [0099] displaying the blended image.
[00100] Some embodiments further comprise: receiving an elemental threshold; and wherein generating a blended image of the sample by blending the elemental and mineral distribution images of the sample comprises rendering a pixel of the blended image black when a corresponding pixel of the elemental distribution image has a value that is less than the elemental threshold.
[00101] Some embodiments provide a computer implemented method for blending mineralogy images of a sample, comprising: [00102] receiving a compositional image of the sample; [00103] receiving a textural image of the sample; [00104] generating, via a computer processor, a blended image of the sample by blending the compositional and textural images of the sample; and [00105] displaying the blended image.
[00106] Some embodiments of the invention provide a charged particle beam apparatus, comprising: [00107] a charged particle beam source for generating a beam of charged particles and directing the beam toward a portion of the sample; [00108] at least one detector for detecting radiation emitted from the illuminated portion of the sample; [00109] a program memory for storing computer instruction to carry out any of methods described above; and [00110] a controller configured to control the charged particle beam system in accordance with the computer instructions stored in the computer memory.
[00111] Some embodiments provide a computer program product, stored on a nontransitory computer readable medium, comprising instructions operable to cause a programmable processor to perform the any of the methods described above.
[00112] Although much of the previous description is directed at mineral samples from drill cores, the invention could be used to prepare samples of any suitable material. The terms “work piece,” “sample,” “substrate,” and “specimen” are used interchangeably in this application unless otherwise indicated. Further, whenever the terms “automatic,” “automated,” or similar terms are used herein, those terms will be understood to include manual initiation of the automatic or automated process or step.
[00113] In the following discussion and in the claims, the terms "including" and "comprising" are used in an open-ended fashion, and thus should be interpreted to mean "including, but not limited to . . . ." To the extent that any term is not specially defined in this specification, the intent is that the term be given its plain and ordinary meaning. The accompanying drawings are intended to aid in understanding the present invention and, unless otherwise indicated, are not drawn to scale. Particle beam systems suitable for carrying out the present invention are commercially available, for example, from FEI Company, the assignee of the present application.
[00114] Although the present invention and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made to the embodiments described herein without departing from the spirit and scope of the invention as defined by the appended claims. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.
[00115] We claim as follows:

Claims (23)

  1. CLAIMS:
    1. A method of analyzing a sample, comprising: directing an electron beam toward multiple points on the surface of a sample; detecting x-rays emitted from the multiple points on the sample surface in response to the electron beam; forming a measured x-ray spectrum from the x-rays emitted from each of the multiple points on the surface; comparing the measured x-ray spectrum from each of the multiple points with reference x-ray spectra to determine the composition at each of the points on the surface; determining textural information at the multiple points on the surface; blending information from the multiple points on the surface to produce a blended sample image composed of pixels, the pixels corresponding to the multiple points on the sample, the color and brightness of each pixel in the blended sample image being determined by a combination of the composition and the textural information at the multiple points.
  2. 2. The method of claim 1 further comprising displaying the blended sample image.
  3. 3. The method of claim 1 or claim 2 in which determining textural information includes determining textural information from back scattered electrons produced by the impact of the electrons on the multiple points on the surface.
  4. 4. The method of claim 1 or claim 2 in which determining textural information includes determining textural information from the intensity of x-ray produced by the impact of the electrons on the multiple points on the surface.
  5. 5. The method of any of the above claims in which comparing the measured x-ray spectrum with reference x-ray spectra to determine the composition of the points on the surface includes producing a top layer image composed of a top layer value for each of the multiple points and in which determining textural information includes producing a base layer image composed of base layer values, and further comprising; determining a luminosity value from one of the base layer or top layer values; determining a multiply value by multiplying the base and top layer values; determining an overlay value from the multiply value or from the top layer value and the bottom layer value; determining an image blending weight from the luminosity value; determining a blended image value from the multiply value and the overlay value; and displaying an image in which the pixels of the image correspond to the blended image values.
  6. 6. The method of claim 5 in which: determining an overlay value comprises determining an overlay value by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value; and determining a blended image value comprises determining a blended image value as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight.
  7. 7. A method of analyzing a sample, comprising: directing an electron beam toward multiple points on the surface of a sample; detecting x-rays emitted from the points on the sample surface in response to the electron beam; forming a measured x-ray spectrum from x-rays emitted from each of the points on the surface; comparing the measured x-ray spectra with reference x-ray spectra to determine the elemental composition of the points on the surface and to determine the mineral composition of the points on the surface; and generating a blended image corresponding to the multiple points on the sample, the blended image including pixels corresponding to the multiple points on the sample, the color and the brightness of the pixels determined in part by the elemental composition at the corresponding point on the sample and determined in part by the mineral composition at the corresponding points on the sample.
  8. 8. The method of claim 7 in which: comparing the measured x-ray spectra with reference x-ray spectra includes determining an elemental distribution map of the sample and determining a mineral distribution map of the sample, each of the maps composed of pixels; and generating a blended image corresponding to the multiple points on the sample includes blending the pixels of the elemental distribution map and the pixels of the mineral distribution map.
  9. 9. The method of claim 7 or claim 8 further comprising receiving an elemental threshold and wherein generating a blended image comprises rendering a pixel of the blended image black when a corresponding pixel of the elemental distribution image has a value that is less than the elemental threshold.
  10. 10. A computer implemented method for blending mineralogy images, comprising: receiving a base layer mineralogy image and a top layer mineralogy image; and sequentially, via a computer processor and for each pixel in the base or top layer mineralogy images: determining a base layer value from the pixel in the base layer mineralogy image; determining a top layer value from the pixel in the top layer mineralogy image; determining a luminosity value from one of the base layer or top layer values; determining a multiply value by multiplying the base and top layer values; determining an overlay value by doubling the multiply value if one of the base layer or top layer values is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value; determining an image blending weight from the luminosity value; and determining a blended image value as a weighted average of the multiply value and the overlay value, wherein the multiply value is weighted by the image blending weight and the overlay value is weighted by the inverted image blending weight.
  11. 11. The computer implemented method of claim 10, wherein the base layer mineralogy image is a textural mineralogy image and the top layer mineralogy image is a spatial mineral distribution image.
  12. 12. The computer implemented method of claim 11, wherein the textural mineralogy image is an image of the intensity of back-scattered electrons.
  13. 13. The computer implemented method of claim 10, wherein the base layer mineralogy image is an image of the intensity of generated x-rays and the top layer mineralogy image is a spatial mineral distribution image.
  14. 14. The computer implemented method of any of claims 10 to 13, wherein the luminosity value is determined from the top layer value.
  15. 15. The computer implemented method of any of claims 10 to 14, wherein the overlay value is determined by doubling the multiply value if the top layer value is less than a threshold, and by otherwise inverting twice the product of the inverted top layer value with the inverted bottom layer value.
  16. 16. The computer implemented method of any of claims 10 to 15, further comprising displaying the blended image on a display device.
  17. 17. The computer implemented method of claim 10, wherein the base layer mineralogy image is a spatial elemental distribution mineralogy image and the top layer mineralogy image is a spatial mineral distribution image.
  18. 18. The computer implemented method of any of claims 10 to 17, further comprising: receiving a mixing parameter; and determining the image blending weight by multiplying the luminosity value by the mixing parameter.
  19. 19. A computer implemented method for blending mineralogy images of a sample, comprising: receiving an elemental distribution image of the sample; receiving a mineral distribution image of the sample; generating, via a computer processor, a blended image of the sample by blending the elemental and mineral distribution images of the sample; and displaying the blended image.
  20. 20. The computer implemented method of claim 19, further comprising: receiving an elemental threshold; and wherein generating a blended image of the sample by blending the elemental and mineral distribution images of the sample comprises rendering a pixel of the blended image black when a corresponding pixel of the elemental distribution image has a value that is less than the elemental threshold.
  21. 21. A computer implemented method for blending mineralogy images of a sample, comprising: receiving a compositional image of the sample; receiving a textural image of the sample; generating, via a computer processor, a blended image of the sample by blending the compositional and textural images of the sample; and displaying the blended image.
  22. 22. A charged particle beam apparatus, comprising: a charged particle beam source for generating a beam of charged particles and directing the beam toward a portion of the sample; at least one detector for detecting radiation emitted from the illuminated portion of the sample; a program memory for storing computer instruction to carry out any of claims 1 through 21; and a controller configured to control the charged particle beam system in accordance with the computer instructions stored in the computer memory.
  23. 23. A computer program product, stored on a non-transitory computer readable medium, comprising instructions operable to cause a programmable processor to perform the method of any of claims 1 to 21.
AU2015229467A 2014-03-11 2015-03-11 Blend modes for mineralogy images Active AU2015229467B9 (en)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US14/204,485 US9183656B2 (en) 2014-03-11 2014-03-11 Blend modes for mineralogy images
US14/204,485 2014-03-11
EP14173457.4 2014-06-23
EP14173457.4A EP2960865A1 (en) 2014-06-23 2014-06-23 Blend modes for mineralogy images
PCT/US2015/019989 WO2015138619A1 (en) 2014-03-11 2015-03-11 Blend modes for mineralogy images

Publications (3)

Publication Number Publication Date
AU2015229467A1 AU2015229467A1 (en) 2016-09-08
AU2015229467B2 AU2015229467B2 (en) 2017-08-03
AU2015229467B9 true AU2015229467B9 (en) 2017-08-10

Family

ID=52727458

Family Applications (1)

Application Number Title Priority Date Filing Date
AU2015229467A Active AU2015229467B9 (en) 2014-03-11 2015-03-11 Blend modes for mineralogy images

Country Status (2)

Country Link
AU (1) AU2015229467B9 (en)
WO (1) WO2015138619A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3772085A1 (en) * 2019-08-01 2021-02-03 FEI Company A method of examining a sample using a charged particle microscope
CN115656053B (en) * 2022-10-19 2024-05-31 山东大学 Rock mineral content testing method and system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100198638A1 (en) * 2007-11-27 2010-08-05 Max Deffenbaugh Method for determining the properties of hydrocarbon reservoirs from geophysical data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2223842B (en) * 1988-09-06 1993-02-03 Shell Int Research Automated mineral identification and rock characterization process
EP2605005A1 (en) * 2011-12-14 2013-06-19 FEI Company Clustering of multi-modal data

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100198638A1 (en) * 2007-11-27 2010-08-05 Max Deffenbaugh Method for determining the properties of hydrocarbon reservoirs from geophysical data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Pointelectronic, "ColorMAP - Realtime processing of 8 elemental mapping images and one analogue image", Halle, Germany, (20050310), pages 1 - 2, URL: http://www.pointelectronic.de/pdf/ColorMAP_Productsheet.pdf, (20150610) *

Also Published As

Publication number Publication date
WO2015138619A1 (en) 2015-09-17
AU2015229467B2 (en) 2017-08-03
AU2015229467A1 (en) 2016-09-08

Similar Documents

Publication Publication Date Title
US9495786B2 (en) Blend modes for mineralogy images
AU2014346694B2 (en) Sub-pixel analysis and display of fine grained mineral samples
Friel et al. Tutorial review: X-ray mapping in electron-beam instruments
EP2546638B1 (en) Clustering of multi-modal data
AU2013335211B2 (en) Automated mineral classification
US11002692B2 (en) Method of examining a sample using a charged particle microscope
US9719950B2 (en) Sample-specific reference spectra library
US11327032B2 (en) Method of examining a sample using a charged particle microscope
US11598733B2 (en) Method of examining a sample using a charged particle microscope
AU2015229467B9 (en) Blend modes for mineralogy images
US20190035597A1 (en) Scanning electron microscope and image processing apparatus
EP2960865A1 (en) Blend modes for mineralogy images
EP3073253A1 (en) Sample-specific reference spectra library
Schwarzer et al. Backscattered electron imaging with an EBSD detector
JP2003240739A (en) X-ray analyzer having region division support function in composition distribution
US20210042943A1 (en) Scanning electron microscope
US20240027377A1 (en) Method of examining a sample using a charged particle microscope
JP2018179862A (en) X-ray image display device, x-ray image display method, and x-ray image display program
Schwarzer et al. Microscopy Education

Legal Events

Date Code Title Description
SREP Specification republished
FGA Letters patent sealed or granted (standard patent)