CN112907579A - Mineralogy parameter fitting analysis method based on multiple Mapping images - Google Patents

Mineralogy parameter fitting analysis method based on multiple Mapping images Download PDF

Info

Publication number
CN112907579A
CN112907579A CN202110324309.6A CN202110324309A CN112907579A CN 112907579 A CN112907579 A CN 112907579A CN 202110324309 A CN202110324309 A CN 202110324309A CN 112907579 A CN112907579 A CN 112907579A
Authority
CN
China
Prior art keywords
mapping
image
fitting
imaging
mirror
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.)
Granted
Application number
CN202110324309.6A
Other languages
Chinese (zh)
Other versions
CN112907579B (en
Inventor
宋昊
姚畅
池国祥
张刚阳
李圻
徐争启
王泽鑫
李娜
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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202110324309.6A priority Critical patent/CN112907579B/en
Publication of CN112907579A publication Critical patent/CN112907579A/en
Application granted granted Critical
Publication of CN112907579B publication Critical patent/CN112907579B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G06T7/337Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

The invention discloses a mineralogy parameter fitting analysis method based on multiple Mapping images, which comprises the following steps of: s1, scanning and imaging the sample by multiple Mapping, and photographing under a mirror to form a picture; s2, adjusting the Mapping scanned image into a gray mosaic image; s3, correction processing is carried out; s4, processing image torque array data; s5, converting the matrix data into an XYZ three-column data worksheet; s6, merging the worksheets; s7, carrying out error elimination processing to obtain a multi-imaging fitting analysis summary table; s8, fitting a regular linear or curve relation according to the imaging data in the multi-imaging fitting analysis summary table to obtain a pixel value equation; and S9, converting the pixel point equation into a relation equation between Mapping parameters according to the pixel values and respective conversion formulas to obtain a rule equation of real parameters. The parameter fitting analysis method is simple to operate, and compared with the traditional dotting analysis, the data volume is rich, and the position fitting is accurate.

Description

Mineralogy parameter fitting analysis method based on multiple Mapping images
Technical Field
The invention relates to a mineralogical parameter fitting analysis method based on multiple Mapping images.
Background
Mapping analysis can be more extensive compared with dotting analysis, comprehensively reflects the physical and chemical information of each region of the whole target sample, and is widely applied to the fields of geology, biology, materials and the like at present.
However, in the prior imaging analysis of Mapping, rules and phenomena are often obtained through naked eye distinguishing and contrast observation, and when a more specific fitting equation relation is involved, the analysis is still performed by means of instrument dotting, and for the data in Mapping, the dotting is still adopted as a main analysis fitting method because real values cannot be derived or other reasons cannot be well applied. However, compared with Mapping imaging analysis, dotting analysis has the disadvantages of complex operation, difficulty in aligning point locations, and too few point locations, so that a novel technical mode with fast and accurate Mapping data and reliable result is required.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide the mineralogy parameter fitting analysis method based on the various Mapping images, which is simple to operate, abundant in data quantity, accurate in position fitting and capable of comprehensively analyzing data relations among the various Mapping images and obtaining a quantitative equation.
The purpose of the invention is realized by the following technical scheme: the method for fitting and analyzing the mineralogical parameters based on various Mapping images comprises the following steps:
s1, scanning and imaging the sample by multiple Mapping, and photographing under a mirror to form a picture;
s2, adjusting the multiple Mapping scanning images into gray mosaic images;
s3, correcting the gray mosaic image and the under-mirror photographed image by using Photoshop software to keep the position, size and shape of minerals in each image consistent;
s4, carrying out image torque matrix data processing on the corrected image in Origin software;
s5, converting the matrix data into an XYZ three-column data worksheet, wherein X, Y is data position, and Z is data value;
s6, combining the working tables of various Mapping and under-mirror images to obtain a multi-imaging fitting summary initial table;
s7, carrying out error elimination processing on the multi-imaging fitting initial table in Origin to obtain a multi-imaging fitting analysis summary table;
s8, fitting a regular linear or curve relation according to the imaging data parameters in the multi-imaging fitting analysis summary table to obtain a pixel value equation;
and S9, converting the pixel point equation into a relation equation among Mapping parameters according to the pixel value and the respective conversion formulas of multiple Mapping, and obtaining a rule equation of real parameters.
Further, in the step S1, the photographing under the mirror needs to be performed under reflected light, so that the boundary and the hole or other mineral inclusion can be clearly distinguished when the photographing under the mirror is formed. The number of pixel points contained in each Mapping scanning imaging graph is not less than 200.
Further, in step S3, before the images are overlapped, the Photoshop software needs to take each Mapping image and the under-mirror photograph to the target mineral in the image and perform the edge smoothing process, and convert the background of the image into white.
Further, the specific implementation method of step S4 is as follows: and (3) photographing each Mapping image and the image under the mirror to convert the Mapping image into a gray matrix: and converting according to the color depth of each pixel point, wherein black is 0, white is 255, and intermediate transition gray is located in a (0, 255) interval.
Further, the specific implementation method of the error elimination processing in step S7 is as follows: deleting pixel points with the pixel value <1 of the picture under the mirror or the pixel value >244 of the picture under the mirror; and simultaneously deleting the pixel points with the burrs at the edges, namely eliminating the edge sawtooth effect.
Further, in step S8, the coefficient of variation of the fitting equation is less than 0.6.
The invention has the beneficial effects that: the invention provides a mineralogy-geochemistry parameter fitting analysis method based on multiple Mapping image pixel values, which is simple to operate, has abundant data volume and accurate position fitting compared with the traditional dotting analysis, and can comprehensively analyze the data relation in multiple Mapping images and obtain a quantitative equation.
Drawings
FIG. 1 is a flow chart of a mineralogical parameter fitting analysis method of the present invention;
FIG. 2 is a summary chart of the effect of the mineralogy-geochemistry parameter fitting analysis method based on various Mapping image pixel values in this embodiment;
fig. 3 is a fitting graph showing the experimental group and the fitting comparison with the control group.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
As shown in FIG. 1, the method for fitting and analyzing mineralogical parameters based on multiple Mapping images comprises the following steps:
s1, scanning and imaging the sample by multiple Mapping, and photographing under a mirror to form a picture; shooting under the mirror needs to be carried out under reflected light, so that the boundary, holes or other mineral bag bodies can be clearly distinguished when the picture is shot under the mirror; the areas tested by various Mapping need to be photographed under the mirror to form a picture and are matched.
The scanning imaging of various Mapping mainly comprises Raman imaging, scanning of electronic probe elements, scanning of EBSD imaging and the like, and the number of pixel points contained in each Mapping scanning imaging graph is not less than 200.
S2, adjusting the multiple Mapping scanning images into gray mosaic images;
s3, correcting the gray mosaic image and the under-mirror photographed image by using Photoshop software to keep the position, size and shape of minerals in each image consistent; before image registration, Photoshop software needs to take each Mapping image and the under-mirror picture of the target mineral in the image to carry out smooth edge cutting processing, and convert the image background into white.
S4, carrying out image torque matrix data processing on the corrected image in Origin software; the specific implementation method comprises the following steps: and (3) photographing each Mapping image and the image under the mirror to convert the Mapping image into a gray matrix: and converting according to the color depth of each pixel point, wherein black is 0, white is 255, and intermediate transition gray is located in a (0, 255) interval.
S5, converting the matrix data into an XYZ three-column data worksheet, wherein X, Y is data position, and Z is data value; the conversion modes and formats of the Mapping images and the under-mirror images are the same, an XYZ column mode is adopted, and the X is firstly arranged, so that missing values are not eliminated.
S6, combining the working tables of various Mapping and under-mirror images to obtain a multi-imaging fitting summary initial table; the multiple imaging fit summary table should contain the XY position as well as the parameter information for multiple imaging, i.e. the Z values in each table.
S7, carrying out error elimination processing on the multi-imaging fitting initial table in Origin to obtain a multi-imaging fitting analysis summary table; the specific implementation method of the error elimination processing is as follows: deleting pixel points with the pixel value <1 of the picture under the mirror or the pixel value >244 of the picture under the mirror; and simultaneously deleting the pixel points with the burrs at the edges, namely eliminating the edge sawtooth effect.
S8, fitting a regular linear or curve relation according to the imaging data parameters in the multi-imaging fitting analysis summary table to obtain a pixel value equation; the coefficient of variation of the fitting relation equation is required to be less than 0.6 so as to ensure the authenticity of data.
And S9, converting the pixel point equation into a relation equation among Mapping parameters according to the pixel value and the respective conversion formulas of multiple Mapping, and obtaining a rule equation of real parameters.
The following experiments further verify the recognition effect of the present invention.
Control group: and (3) carrying out element dotting analysis test (13 points) on Co in the pyrite sample by adopting an electronic probe, and simultaneously carrying out Raman dotting at the same dotting position, wherein the sample sources of the control group and the experimental group are the same region of the same sample.
Experimental groups: a 556nm laser light source is selected, a 10-fold focal length lens is adopted as an objective lens, after a target visual field is found, Raman surface scanning is carried out on the same area of the same sample, and the Raman displacement scanning range is 300-400 cm-1. Placing the dotted pyrite sample slice under the field of view of an objective lens of an objective table, shooting in the same area under reflected light, obviously distinguishing boundaries and holes or other mineral inclusion, and finally carrying out imaging scanning on the same area of an electronic probe Co, wherein the step length is 0.8 mu m, and the area range is the same as that of the former two.
After the test was completed, the data was processed by the peak clipping method in labspec6 of HORIBA JOBIN YVON corporation to obtain a mosaic gray image with respect to raman Ag intensity, and a Co gray scan image was obtained by a JEOL-JXA-8230 micro-area X-ray spectrometer attached instrument manufactured by japan electronics corporation.
And (3) carrying out coincidence and correction processing on the Raman intensity imaging image, the Co imaging image and the under-mirror image on various Mapping and under-mirror images in Photoshop software, so that pixel points of the three images coincide and are equal in number.
And (4) carrying out respective image conversion matrix data processing on all corrected images in Origin software, wherein the processing process method is kept consistent.
The matrices of the various Mapping and under-mirror maps were each converted to XYZ three-column data worksheets (XY is data position, Z is data value). And combining the working tables of various Mapping and under-mirror images to obtain a multi-imaging fitting summary initial table.
A series of error elimination treatments are carried out on the assembly image fitting primary table in Origin to obtain a multi-imaging fitting analysis summary table, wherein the error elimination treatments comprise that a picture pixel value under a mirror is deleted to be less than 1or 'a picture pixel value under the mirror' >244, namely holes, cracks, backgrounds and mineral inclusion in a sample to be detected are eliminated, for example, a picture under the mirror is { x I x >244} or a picture under the mirror is { x I x <1}, meanwhile, pixel points with burrs at the edges are deleted, namely, an XY equation about positions is deleted, and edge sawtooth effects are eliminated. For example { (x, y) | y > -0.0857x +523.2} - { (x, y) | x <222}, { (x, y) | y >0.2016x +458.84} - { (x, y) | x >222}, and the like, finally, an effect diagram (fig. 2) with respect to Co and raman intensity is obtained.
And (3) fitting a regular linear or curve relation by adopting a median fitting method for Raman intensity according to imaging data parameters in the multi-imaging fitting analysis summary table to obtain a pixel value equation chart 3.
And converting the pixel point equation into a relation equation among Mapping parameters according to the pixel value and respective conversion formulas of multiple Mapping to obtain a rule equation of real parameters.
Fig. 3 shows that the raman Intensity (Intensity) and the Co content in the area of the pyrite sample exhibit a relatively obvious inverse relationship through an experimental composition method, so that visual information obtained in the picture is well decoded, and further, fitting information is specifically given. Meanwhile, the results obtained by the control group by using dotting analysis and the experimental group by using the mineralogy-geochemistry parameter fitting analysis method based on the pixel values of the various Mapping images are approximately consistent, and the results of the mineralogy-geochemistry parameter fitting analysis method based on the pixel values of the various Mapping images are reliable.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Those skilled in the art can make various other specific changes and combinations based on the teachings of the present invention without departing from the spirit of the invention, and these changes and combinations are within the scope of the invention.

Claims (7)

1. The method for fitting and analyzing the mineralogical parameters based on various Mapping images is characterized by comprising the following steps of:
s1, scanning and imaging the sample by multiple Mapping, and photographing under a mirror to form a picture;
s2, adjusting the multiple Mapping scanning images into gray mosaic images;
s3, correcting the gray mosaic image and the under-mirror photographed image by using Photoshop software to keep the position, size and shape of minerals in each image consistent;
s4, carrying out image torque matrix data processing on the corrected image in Origin software;
s5, converting the matrix data into an XYZ three-column data worksheet, wherein X, Y is data position, and Z is data value;
s6, combining the working tables of various Mapping and under-mirror images to obtain a multi-imaging fitting summary initial table;
s7, carrying out error elimination processing on the multi-imaging fitting initial table in Origin to obtain a multi-imaging fitting analysis summary table;
s8, fitting a regular linear or curve relation according to the imaging data parameters in the multi-imaging fitting analysis summary table to obtain a pixel value equation;
and S9, converting the pixel point equation into a relation equation among Mapping parameters according to the pixel value and the respective conversion formulas of multiple Mapping, and obtaining a rule equation of real parameters.
2. The method for fitting analysis of mineralogical parameters based on multiple Mapping images of claim 1, wherein in the step S1, the photographing under the mirror needs to be performed under reflected light, so as to ensure that the boundary and the hole or other mineral inclusion can be clearly distinguished when the photographing under the mirror is taken.
3. The method for fitting analysis of mineralogical parameters based on multiple Mapping images of claim 1, wherein in the step S1, the number of pixel points included in each Mapping scan image is not less than 200.
4. The method for mineralogical parameter fitting analysis based on multiple Mapping images as claimed in claim 1, wherein in step S3, before image registration, Photoshop software is required to perform smooth trimming on each Mapping image and the target mineral in the image photographed under the mirror, and convert the background of the image into white.
5. The method for fitting and analyzing mineralogical parameters based on multiple Mapping images as claimed in claim 1, wherein the step S4 is implemented by: and (3) photographing each Mapping image and the image under the mirror to convert the Mapping image into a gray matrix: and converting according to the color depth of each pixel point, wherein black is 0, white is 255, and intermediate transition gray is located in a (0, 255) interval.
6. The method for fitting and analyzing mineralogical parameters based on Mapping images of claim 1, wherein the error elimination in the step S7 is realized by the following steps: deleting pixel points with the pixel value <1 of the picture under the mirror or the pixel value >244 of the picture under the mirror; and simultaneously deleting the pixel points with the burrs at the edges, namely eliminating the edge sawtooth effect.
7. The method for fitting analysis of mineralogical parameters based on Mapping images of claim 1, wherein in the step S8, the coefficient of variation of the fitting relation equation is less than 0.6.
CN202110324309.6A 2021-03-26 2021-03-26 Mineralogy parameter fitting analysis method based on multiple Mapping images Active CN112907579B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110324309.6A CN112907579B (en) 2021-03-26 2021-03-26 Mineralogy parameter fitting analysis method based on multiple Mapping images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110324309.6A CN112907579B (en) 2021-03-26 2021-03-26 Mineralogy parameter fitting analysis method based on multiple Mapping images

Publications (2)

Publication Number Publication Date
CN112907579A true CN112907579A (en) 2021-06-04
CN112907579B CN112907579B (en) 2022-06-21

Family

ID=76108729

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110324309.6A Active CN112907579B (en) 2021-03-26 2021-03-26 Mineralogy parameter fitting analysis method based on multiple Mapping images

Country Status (1)

Country Link
CN (1) CN112907579B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6675064B1 (en) * 1999-09-27 2004-01-06 University Of Kentucky Research Foundation Process for the physical segregation of minerals
CN106815805A (en) * 2017-01-17 2017-06-09 湖南优象科技有限公司 Rapid distortion bearing calibration based on Bayer images
CN109839369A (en) * 2019-03-11 2019-06-04 成都理工大学 A method of the graphite degree of order is measured based on LR laser raman Mapping

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6675064B1 (en) * 1999-09-27 2004-01-06 University Of Kentucky Research Foundation Process for the physical segregation of minerals
CN106815805A (en) * 2017-01-17 2017-06-09 湖南优象科技有限公司 Rapid distortion bearing calibration based on Bayer images
CN109839369A (en) * 2019-03-11 2019-06-04 成都理工大学 A method of the graphite degree of order is measured based on LR laser raman Mapping

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JOHANNES TRAN-GIA 等: "Model-based Acceleration of Parameter mapping (MAP) for saturation prepared radially acquired data", 《MAGNETIC ROSONANCE IN MEDICINE》 *
李娜: "褪色文物模型色彩重建技术研究", 《中国博士学位论文全文数据库 (信息科技辑)》 *

Also Published As

Publication number Publication date
CN112907579B (en) 2022-06-21

Similar Documents

Publication Publication Date Title
US7065236B2 (en) Method for quantitative video-microscopy and associated system and computer software program product
JP5940383B2 (en) Microscope system
Franke et al. Measurement of microdisplacements by machine vision photogrammetry (DISMAP)
Levanič ATRICS–A new system for image acquisition in dendrochronology
CN111626936B (en) Quick panoramic stitching method and system for microscopic images
ES2542106T3 (en) Method for preparing quantitative videomicroscopy and associated system
AU2002334590A1 (en) Method quantitative video-microscopy and associated system and computer software program product
US10746980B2 (en) Calibration of microscopy systems
MX2007016046A (en) System and method for re-locating an object in a sample on a slide with a microscope imaging device.
JP7146635B2 (en) Autofocus method and system for digital imaging using multispectral trajectories
RU2009140389A (en) DEVICE AND METHODS OF EXAMINATION, EVALUATION AND CLASSIFICATION OF PRECIOUS STONES
CN108109157B (en) Rock mass evaluation analysis method based on digital panoramic borehole image
CN112907579B (en) Mineralogy parameter fitting analysis method based on multiple Mapping images
EP1492152A2 (en) Electron microscope
US20220189048A1 (en) Generating composite image from multiple images captured for subject
WO2001016884A1 (en) Gel electrophoresis image combining for improved dynamic range
CN115439477B (en) 24 color card positioning method and device, electronic equipment and storage medium
Nederbragt et al. Image calibration, filtering, and processing
US11644359B2 (en) Method of reading the result of an electrophoretic assay comprising a digital image indicating the intensity of light emitted by chemiluminescence from the output medium of the electrophoretic assay
CN108520499A (en) A kind of image shift antidote based on leukorrhea micro-imaging
WO2007054804A2 (en) Digital inspection of the physical quality of plain surfaces
CN112884821A (en) Method for making super-long train template image
CN111553871B (en) Method for judging concentration of metal element in binary alloy by using gray value of synchronous radiation in-situ imaging image
CN113295666A (en) Quantitative analysis method for As element in pyrite by using mineral Raman parameters
US20030043285A1 (en) Image data forming method and apparatus

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
GR01 Patent grant
GR01 Patent grant