WO2012165991A1 - Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample - Google Patents

Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample Download PDF

Info

Publication number
WO2012165991A1
WO2012165991A1 PCT/RU2011/000378 RU2011000378W WO2012165991A1 WO 2012165991 A1 WO2012165991 A1 WO 2012165991A1 RU 2011000378 W RU2011000378 W RU 2011000378W WO 2012165991 A1 WO2012165991 A1 WO 2012165991A1
Authority
WO
WIPO (PCT)
Prior art keywords
section
cross
distribution
grayness
tomographic image
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.)
Ceased
Application number
PCT/RU2011/000378
Other languages
French (fr)
Inventor
Valery Vasilievich Shako
Nikita llyich RYZHIKOV
Dmitry Nikolaevich MIKHAILOV
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.)
Schlumberger Canada Ltd
Services Petroliers Schlumberger SA
Schlumberger Technology BV
Schlumberger Holdings Ltd
Prad Research and Development Ltd
Original Assignee
Schlumberger Canada Ltd
Services Petroliers Schlumberger SA
Schlumberger Technology BV
Schlumberger Holdings Ltd
Prad Research and Development Ltd
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 Schlumberger Canada Ltd, Services Petroliers Schlumberger SA, Schlumberger Technology BV, Schlumberger Holdings Ltd, Prad Research and Development Ltd filed Critical Schlumberger Canada Ltd
Priority to US14/119,661 priority Critical patent/US9008372B2/en
Priority to PCT/RU2011/000378 priority patent/WO2012165991A1/en
Publication of WO2012165991A1 publication Critical patent/WO2012165991A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • 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/02Investigating 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 transmitting the radiation through the material
    • G01N23/04Investigating 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 transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating 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 transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • 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/419Imaging computed tomograph
    • 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
    • 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/646Specific applications or type of materials flaws, defects
    • 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/649Specific applications or type of materials porosity

Definitions

  • the present invention relates generally to the methods for determination of spatial distribution and concentration of contrast components (mineral grains, components inside pore space, inclusions) in materials by analysing X-ray tomography data and can be used, for example, for calculation of fraction of contrast minerals in rock sample, diagnostic of tumor (cancerous or innocent) sizes in medicine (after using of special contrast agents), crack detection of composite and geo- materials.
  • contrast components mineral grains, components inside pore space, inclusions
  • X-ray micro Computed Tomography This technique provides 3D object structure composed of cross- section images of the internal structure with acquisition of grayscale which represents the X-ray absorption distribution (attenuation coefficients) within the object.
  • grayscale represents the X-ray absorption distribution (attenuation coefficients) within the object.
  • Computer tomography technique is based on the interaction of X-rays with material. Passing through an object X-rays will be attenuated depending on the physical density and atomic number of the studied object and on the used X-ray energies. This attenuation information is collected on 2D XmCT image.
  • each point of the black- white slice of object is characterized by different grayness. Because attenuation coefficient depends on material, through X-ray has passed, different materials are characterized by different grayness that allow us to separate individual materials and estimate it's fraction on each X-ray section.
  • Thresholding is method to separate object of interest from "background” (i.e. other objects) which is based on choosing of optimal threshold level of grayness. All points ("pixels") of the black-white X-ray image of object which grayness is lower then the threshold value are supposed to belong to object (or background, depending what is more bright).
  • This invention describes new method which is based on solving of inverse problem while analysis of histogram and provide deterministic solution.
  • number of the components is defined as the number of peaks on the histogram of grayness distribution inside the area of interest on the reference cross-section.
  • i index of a component
  • / is "intensity" (total number of pixels) of grayness z
  • A) , B) , C) are adjustable parameters
  • Adjustable parameters for all components are crudely estimated from the histogram of grayness distribution inside the area of interest. Accurate estimation of adjustable parameters for the components is made by minimization of module of difference between the histogram of grayness distribution inside the area of interest of the reference cross-section and sum of normal distributions corresponding to individual components
  • j index of grayscale
  • M total number of grayscales
  • N total number of different components.
  • the sub-areas each containing only one individual component are selected inside the area of interest on the reference cross-section and histograms of grayness distribution of the individual components are obtained. All histograms are normalized by their area, i.e. by number of pixels. The histograms of all individual components are transformed to general scale. Histogram of the area of interest of the reference cross-section is approximated by sum of histograms of individual components with some weight coefficients, corresponding to the areas that individual components occupy on this cross-section of the computer tomographic image.
  • Weight coefficients b c 1 etc correspond to areas that individual components occupied on the reference cross-section of the histograms of grayness distribution in the cross-sections of the computer tomographic image.
  • Fig. 1 shows an example of selecting an area of interest on the cross- section of the computer tomographic image and obtaining the histogram of grayness distribution ;
  • Fig. 2 shows an example of selecting the sub-areas, containing only one individual material, inside the total area of interest and obtaining the histogram of grayness distribution;
  • Fig. 3 shows an approximation of histogram using Gaussians as probability density functions
  • Fig. 4 shows an example of profile for mixture of three different materials.
  • a histogram of grayness distribution inside the area of interest of the reference cross-section is obtained using special program (for, example, ImageJ - http://rsbweb.nih.gov/ij/), see Fig. 1.
  • the number of the components is defined as the number of peaks on the histogram of grayness distribution inside the area of interest on the reference cross-section.
  • Adjustable parameters for all components are crudely estimated from the histogram of grayness distribution inside the area of interest. Accurate estimation of adjustable parameters for the components is made by minimization of module of difference between the histogram of grayness distribution inside the area of interest of the reference cross-section and sum of normal distributions, corresponding to individual components
  • Typical example of reconstructed profile for mixture of three different components is shown on Fig.4.
  • a porous sample (artificial sample, consisting of sand grains and liquid glass as cement material) is scanned with X-ray and a computer tomographic image of the sample is obtained.
  • the sub-areas each containing only one individual component are selected inside the area of interest on the reference cross-section and histograms of grayness distribution of the individual components are obtained, see Fig.2.
  • All histograms are normalized by their area, i.e. by number of pixels. This way amount of pixels in every range of histogram must be divided by amount of all selected pixels.
  • the histograms of all individual components are transformed to general scale, for example, to scale of histogram of total area of interest. Histogram of the area of interest of the reference cross-section is approximated by sum of histograms of individual components with some weight coefficients, corresponding to the areas that individual components occupy on thiscross-section of the computer tomographic image.
  • Weight coefficients b c 1 etc correspond to areas that individual components occupied on the reference cross-section of the X-ray image.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

A method for determination of a spatial distribution and concentration of contrast components in a porous sample comprises the steps of scanning a sample with X-ray and obtaining a computer tomographic image of the sample. Then an area of interest inside the obtained computer tomographic image is selected and a first cross-section of the computer tomographic image is defined. Spatial distribution and concentration of contrast components inside the area of interest are determined by analyzing histograms of grayness distribution in the cross-sections of the computer tomographic image starting with the reference cross-section.

Description

Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample
Field of the invention
The present invention relates generally to the methods for determination of spatial distribution and concentration of contrast components (mineral grains, components inside pore space, inclusions) in materials by analysing X-ray tomography data and can be used, for example, for calculation of fraction of contrast minerals in rock sample, diagnostic of tumor (cancerous or innocent) sizes in medicine (after using of special contrast agents), crack detection of composite and geo- materials.
Background of the invention
Powerful nondestructive method for analysis of contrast components in heterogeneous or porous medium is X-ray micro Computed Tomography (XmCT). This technique provides 3D object structure composed of cross- section images of the internal structure with acquisition of grayscale which represents the X-ray absorption distribution (attenuation coefficients) within the object. The idea of correlation between the gray values distribution in X- ray CT image and material densities distribution is described in US Pat. No. 2005/0010106).
Computer tomography technique is based on the interaction of X-rays with material. Passing through an object X-rays will be attenuated depending on the physical density and atomic number of the studied object and on the used X-ray energies. This attenuation information is collected on 2D XmCT image.
Depending on X-ray attenuation coefficient, each point of the black- white slice of object is characterized by different grayness. Because attenuation coefficient depends on material, through X-ray has passed, different materials are characterized by different grayness that allow us to separate individual materials and estimate it's fraction on each X-ray section.
Most widespread method of image recognizing is thresholding. Thresholding is method to separate object of interest from "background" (i.e. other objects) which is based on choosing of optimal threshold level of grayness. All points ("pixels") of the black-white X-ray image of object which grayness is lower then the threshold value are supposed to belong to object (or background, depending what is more bright).
Special type of this method is histogram based thresholding which is developed for case when an image have only two principal ("dark" and "light") grayscale regions (see, for example, Gupta L., Sortrakul T. A Gaussian-mixture based image segmentation algorithm. Pattern Recognition, Vol. 31 , No 3, p.315— 325,1998).
Principle disadvantage of thresholding technique is sensitivity of results to selected value of threshold and thus requires a priori information about analyzed parametres.
This invention describes new method which is based on solving of inverse problem while analysis of histogram and provide deterministic solution.
Summary of the invention
A method for determination of spatial distribution and concentration of contrast components in a porous sample comprises the steps of scanning a sample with X-ray and obtaining a computer tomographic image of the sample. Then an area of interest inside the obtained computer tomographic image is selected and a first cross-section of the computer tomographic image (let us call it as reference cross-section and assign a number k = 1) is defined. Histograms of grayness distribution in the cross-sections of the computer tomographic image are obtained. Spatial distribution and concentration of contrast components inside the area of interest are determined by analyzing histograms of grayness distribution in the cross-sections of the computer tomographic image, starting with the reference cross-section.
According to the first embodiment of the invention number of the components is defined as the number of peaks on the histogram of grayness distribution inside the area of interest on the reference cross-section. The histogram l) {z) of grayness distribution for each component inside the reference cross-section (k = 1) is approximated by normal distribution (Gaussian function :
Figure imgf000004_0001
wherein i is index of a component; / is "intensity" (total number of pixels) of grayness z ; A) , B) , C) are adjustable parameters; upper index "1" corresponds to number of cross-section (k = 1). Adjustable parameters for all components are crudely estimated from the histogram of grayness distribution inside the area of interest. Accurate estimation of adjustable parameters for the components is made by minimization of module of difference between the histogram of grayness distribution inside the area of interest of the reference cross-section and sum of normal distributions corresponding to individual components
Figure imgf000004_0002
wherein j is index of grayscale; M is total number of grayscales; N is total number of different components.
The fraction a) of the individual components inside the reference cross- section inside of the com uter tomographic image is calculated as:
Figure imgf000005_0001
Obtained parameters of Gaussians A) , B] , C) of the reference cross- section are used as initial parameters for minimization of module of difference between the real histogram of grayness distribution and it's approximation by sum of normal distributions (1) for next cross-section of X- ray image (k = 2) and so an.
To reconstruct the distribution and concentration of the contrast components along the sample, the mininization of module of difference between the real histogram of grayness distribution and it's approximation by sum of normal distributions (1) and the expression (2) are applied to all cross- sections inside area of interest of the computer tomographic image (k =
1 . ..K).
According to the second embodiment of the invention the sub-areas each containing only one individual component are selected inside the area of interest on the reference cross-section and histograms of grayness distribution of the individual components are obtained. All histograms are normalized by their area, i.e. by number of pixels. The histograms of all individual components are transformed to general scale. Histogram of the area of interest of the reference cross-section is approximated by sum of histograms of individual components with some weight coefficients, corresponding to the areas that individual components occupy on this cross-section of the computer tomographic image. Weight coefficients are obtained by minimization of module of difference between the real histogram of the area of interest of the reference cross-section (k = 1) and sum of histograms of individual components:
Figure imgf000006_0001
(3) wherein A) , B) , C) , ... are vectors of values of histograms, bx , c1 etc are weight coefficients; j is index of grayscale; M is total number of grayscales; upper index "1" corresponds to number of cross-section (k = 1).
Weight coefficients b c1 etc correspond to areas that individual components occupied on the reference cross-section of the histograms of grayness distribution in the cross-sections of the computer tomographic image.
To reconstruct distribution and concentration of the contrast components along the sample, the procedure described above and including the mininization (3), is applied to all cross-sections inside the area of interest of the computer tomographic image (k = 1...K).
In the case of bad convergence of this problem, it is possible to extend the sub-areas containing individual components.
Brief description of the drawings
Fig. 1 shows an example of selecting an area of interest on the cross- section of the computer tomographic image and obtaining the histogram of grayness distribution ;
Fig. 2 shows an example of selecting the sub-areas, containing only one individual material, inside the total area of interest and obtaining the histogram of grayness distribution;
Fig. 3 shows an approximation of histogram using Gaussians as probability density functions; Fig. 4 shows an example of profile for mixture of three different materials.
Detailed description of the invention
According to the first embodiment of the invention a porous sample (artificial sample, consisting of sand grains and liquid glass as cement material) is scanned with X-ray and a computer tomographic image of the sample is obtained. Then an area of interest inside of this computer tomographic image is selected and a first cross-section of the computer tomographic image (let us call it as reference cross-section and assign a number k = 1) is defined. Under an area of interest we understand a sub-area of 3D X-ray computer tomographic image which is selected for detailed analysis. This area can be selected because it includes some specific features (microfractures, defect or specific inclusion) or as typical volume of image to reduce simulation time (if initial 3D X-ray image is too large for analysis).
A histogram of grayness distribution inside the area of interest of the reference cross-section is obtained using special program (for, example, ImageJ - http://rsbweb.nih.gov/ij/), see Fig. 1.
The number of the components is defined as the number of peaks on the histogram of grayness distribution inside the area of interest on the reference cross-section. The histogram l) {z) of grayness distribution for each component inside reference cross-section (k = 1) is approximated by normal distribution Gaussian function):
Figure imgf000007_0001
wherein / is index of a component; / is "intensity" (total number of pixels) of grayness z ; A) , B) , C) are adjustable parameters; upper index "1" corresponds to number of cross-section (k = 1).
Example of hystogram of selected individual material is given in Fig. 2.
Adjustable parameters for all components are crudely estimated from the histogram of grayness distribution inside the area of interest. Accurate estimation of adjustable parameters for the components is made by minimization of module of difference between the histogram of grayness distribution inside the area of interest of the reference cross-section and sum of normal distributions, corresponding to individual components
Figure imgf000008_0001
Figure imgf000008_0004
wherein j is index of grayscale; M is total number of grayscales; N is total number of different components. Approximation result is shown on Fig.3.
The fraction a) of the individual components inside the reference cross- section of the com uter tomographic image is calculated as:
Figure imgf000008_0002
Obtained parameters of Gaussians
Figure imgf000008_0003
, B) , C) of the reference cross- section are used as initial parameters for minimization of module of difference between the real histogram of grayness distribution and it's approximation by sum of normal distributions (1) for next cross-section of the computer tomographic image {k = 2) and so an. To reconstruct the distribution and concentration of the contrast components along the sample the mininization of module of difference between the real histogram of grayness distribution and it's approximation by sum of normal distributions (1) and the expression (2) are applied to all cross- sections inside the area of interest of the computer tomographic image (k = 1 ...K).
Typical example of reconstructed profile for mixture of three different components is shown on Fig.4.
According to the second embodiment of the invention a porous sample (artificial sample, consisting of sand grains and liquid glass as cement material) is scanned with X-ray and a computer tomographic image of the sample is obtained.
Then an area of interest inside of this computer tomographic image is selected and a first cross-section of the computer tomographic image (let us call it as reference cross-section and assign a number k = 1) is defined. A histogram of grayness distribution inside the area of interest of the reference cross-section is obtained using special program (for, example, ImageJ - http://rsbweb.nih.gov/ij/), see Fig. 1.
Then, the sub-areas each containing only one individual component are selected inside the area of interest on the reference cross-section and histograms of grayness distribution of the individual components are obtained, see Fig.2. All histograms are normalized by their area, i.e. by number of pixels. This way amount of pixels in every range of histogram must be divided by amount of all selected pixels. The histograms of all individual components are transformed to general scale, for example, to scale of histogram of total area of interest. Histogram of the area of interest of the reference cross-section is approximated by sum of histograms of individual components with some weight coefficients, corresponding to the areas that individual components occupy on thiscross-section of the computer tomographic image.
Weight coefficients are obtained by minimization of module of difference between the real histogram of the area of interest of reference cross-section k = 1) and sum of histograms of individual components:
Figure imgf000010_0001
wherein A) , B , C) , ... are vectors of values of histograms, b1 , c1 etc are weight coefficients; j is index of grayscale; M is total number of grayscales; upper index "1" corresponds to number of cross-section (k = 1).
Weight coefficients b c1 etc correspond to areas that individual components occupied on the reference cross-section of the X-ray image.
To reconstruct distribution and concentration of the contrast components along the sample, the procedure described above and including the mininization (3), is applied to all cross-sections inside area of interest of the computer tomographic image (k = I ...K).
In the case of bad convergence of this problem, it is possible to extend the sub-areas containing individual components or check that all principle materials are taken into account.

Claims

1. A method for determination of spatial distribution and concentration of contrast components in a porous or / and heterogeneous sample comprising the following steps:
- scanning the sample with X-ray,
- obtaining a computer tomographic image of the sample,
- defining an area of interest inside the obtained computer tomographic image,
- defining a reference cross-section of the obtained computer tomographic image ,
- obtaining histograms of grayness distribution in the cross-sections of the computer tomographic image,
- determining spatial distribution and concentration of contrast components along the sample by analyzing histograms of grayness distribution in the cross-sections of the computer tomographic image starting with a histogram of grayness distribution in the reference cross-section.
2. A method of claim 1 wherein a step of analyzing histograms of grayness distribution in the cross-sections of the computer tomographic image comprises the following:
- defining the number of the components as the number of peaks on the histogram of grayness distribution inside the area of interest of the reference cross-section,
- approximating the histogram l) {∑) of grayness distribution for each component inside the reference cross-section by normal distribution (Gaussian function):
Figure imgf000012_0001
wherein
/' is an index of a component;
/ is "intensity" (total number of pixels) of grayness z ;
A) , B] , C) are adjustable parameters corresponding to the reference cross-section;
upper index "1" corresponds to the reference cross-section,
- crudely estimating adjustable parameters for all components from the histogram of grayness distribution inside the area of interest of the reference cross-section,
- accurately estimating adjustable parameters for the components by minimization of module of difference between the histogram of grayness distribution inside the area of interest of the reference cross-section and sum of normal distributions corresponding to individual components
Figure imgf000012_0002
wherein
j is index of grayscale;
M is total number of grayscales;
N is total number of different components,
- using obtained parameters of Gaussians A) , B] , C) as initial parameters for minimization of module of difference between the histogram of grayness distribution and it's approximation by sum of normal distributions for the next cross-section of the computer tomographic image, - applying previous step to all cross-sections inside the area of interest of the computer tomographic image,
- defining the fraction of the individual components inside each cross section k b integration of individual Gaussians:
Figure imgf000013_0001
wherein k = 1...K is the number of cross-section.
3. The method of claim 1 wherein a step of analyzing grayness distribution in the cross-section of the computer tomographic image comprises the following:
- selecting the sub-areas inside the area of interest on the reference cross-section, each containing only one individual component,
- obtaining histograms of grayness distribution of the individual components,
- normalizing all histograms by their area,
- transforming the histograms of all individual components to general scale,
- approximating the histogram of the area of interest on the reference cross-section by sum of histograms of individual components with the weight coefficients, corresponding to the areas occupied by the individual components on this cross-section of the computer tomographic image,
- defining weight coefficients by minimization of module of difference between the histogram of the area of interest on the reference cross-section and sum of histograms of individual components:
Figure imgf000013_0002
wherein
A) , B) , C) , . .. are vectors of values of histograms, b1 , c1 , ... are weight coefficients;
j is index of grayscale;
M is total number of grayscales;
- applying previous step to all cross-sections inside the area of interest of computer tomographic image.
4. Method of claim 3 wherein the sub-areas containing individual components are expanded.
PCT/RU2011/000378 2011-05-31 2011-05-31 Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample Ceased WO2012165991A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US14/119,661 US9008372B2 (en) 2011-05-31 2011-05-31 Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample
PCT/RU2011/000378 WO2012165991A1 (en) 2011-05-31 2011-05-31 Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/RU2011/000378 WO2012165991A1 (en) 2011-05-31 2011-05-31 Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample

Publications (1)

Publication Number Publication Date
WO2012165991A1 true WO2012165991A1 (en) 2012-12-06

Family

ID=47259590

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/RU2011/000378 Ceased WO2012165991A1 (en) 2011-05-31 2011-05-31 Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample

Country Status (2)

Country Link
US (1) US9008372B2 (en)
WO (1) WO2012165991A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2515009C1 (en) * 2012-12-21 2014-05-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Ярославский государственный технический университет" Method of determining coefficient of heterogeneity of mixture of hard-to-separate granular materials
RU2620387C1 (en) * 2016-03-21 2017-05-25 Федеральное государственное бюджетное образовательное учреждение высшего образования "Ярославский государственный технический университет" (ФГБОУ ВО "ЯГТУ") Method of determining coefficient of hard-to-separate bulk materials heterogeneous mixture
US9951613B2 (en) 2015-06-11 2018-04-24 Schlumberger Technology Corporation Method for quantitative analysis of a distribution of solid particles of a suspension invaded in a porous medium during filtration
FR3097328A1 (en) * 2019-06-17 2020-12-18 Safran Method, device and diagnostic program for densification faults of a CMC part

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106528997B (en) * 2016-10-28 2020-08-14 江苏天瑞仪器股份有限公司 Method for drawing regional particle hourly concentration distribution map

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6738144B1 (en) * 1999-12-17 2004-05-18 University Of Central Florida Non-invasive method and low-coherence apparatus system analysis and process control
US20050010106A1 (en) * 2003-03-25 2005-01-13 Imaging Therapeutics, Inc. Methods for the compensation of imaging technique in the processing of radiographic images
US7319739B2 (en) * 2003-03-14 2008-01-15 Siemens Aktiengesellschaft Imaging method based on two different x-ray spectra
RU84548U1 (en) * 2009-03-19 2009-07-10 Общество С Ограниченной Ответственностью "Лаборатория Оптико-Электронных Приборов" DEVICE FOR MEASURING THE DISTRIBUTION OF THE SIZES AND CONCENTRATIONS OF NON-SPHERICAL NANOPARTICLES IN LIQUIDS AND GASES

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4991092A (en) * 1988-08-12 1991-02-05 The Regents Of The University Of California Image processor for enhancing contrast between subregions of a region of interest
US5027379A (en) * 1990-02-22 1991-06-25 Bp America Inc. Method for identifying drilling mud filtrate invasion of a core sample from a subterranean formation
JP3478566B2 (en) 1993-09-16 2003-12-15 株式会社東芝 X-ray CT scanner
GB0115615D0 (en) * 2001-06-27 2001-08-15 Univ Coventry Image segmentation
US6985612B2 (en) * 2001-10-05 2006-01-10 Mevis - Centrum Fur Medizinische Diagnosesysteme Und Visualisierung Gmbh Computer system and a method for segmentation of a digital image

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6738144B1 (en) * 1999-12-17 2004-05-18 University Of Central Florida Non-invasive method and low-coherence apparatus system analysis and process control
US7319739B2 (en) * 2003-03-14 2008-01-15 Siemens Aktiengesellschaft Imaging method based on two different x-ray spectra
US20050010106A1 (en) * 2003-03-25 2005-01-13 Imaging Therapeutics, Inc. Methods for the compensation of imaging technique in the processing of radiographic images
RU84548U1 (en) * 2009-03-19 2009-07-10 Общество С Ограниченной Ответственностью "Лаборатория Оптико-Электронных Приборов" DEVICE FOR MEASURING THE DISTRIBUTION OF THE SIZES AND CONCENTRATIONS OF NON-SPHERICAL NANOPARTICLES IN LIQUIDS AND GASES

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2515009C1 (en) * 2012-12-21 2014-05-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Ярославский государственный технический университет" Method of determining coefficient of heterogeneity of mixture of hard-to-separate granular materials
US9951613B2 (en) 2015-06-11 2018-04-24 Schlumberger Technology Corporation Method for quantitative analysis of a distribution of solid particles of a suspension invaded in a porous medium during filtration
RU2620387C1 (en) * 2016-03-21 2017-05-25 Федеральное государственное бюджетное образовательное учреждение высшего образования "Ярославский государственный технический университет" (ФГБОУ ВО "ЯГТУ") Method of determining coefficient of hard-to-separate bulk materials heterogeneous mixture
FR3097328A1 (en) * 2019-06-17 2020-12-18 Safran Method, device and diagnostic program for densification faults of a CMC part

Also Published As

Publication number Publication date
US20140192953A1 (en) 2014-07-10
US9008372B2 (en) 2015-04-14

Similar Documents

Publication Publication Date Title
Wang et al. Comparison of image segmentation methods in simulated 2D and 3D microtomographic images of soil aggregates
US7646902B2 (en) Computerized detection of breast cancer on digital tomosynthesis mammograms
Rabbani et al. Comparing three image processing algorithms to estimate the grain-size distribution of porous rocks from binary 2D images and sensitivity analysis of the grain overlapping degree
JP6475691B2 (en) Method and x-ray system for computer-aided detection of structures in x-ray images
Karimpouli et al. Conditional reconstruction: An alternative strategy in digital rock physics
Abera et al. The effectiveness of global thresholding techniques in segmenting two-phase porous media
CN104700389B (en) Object recognition method in dual-energy CT scan images
Kazantsev et al. Statistical detection of defects in radiographic images in nondestructive testing
US9008372B2 (en) Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample
Wu et al. Extraction of fractures in shale CT images using improved U-Net
Gupta et al. A new workflow of X-ray CT image processing and data analysis of structural features in rock using open-source software
Nguyen et al. An effective deep learning architecture combination for tissue microarray spots classification of H&E stained colorectal images
Wang et al. Machine learning and UNet++ based microfracture evaluation from CT images
Mitra et al. X-ray tomography applied to tsunami deposits: Optimized image processing and quantitative analysis of particle size, particle shape, and sedimentary fabric in 3D
Eddaoudi et al. Microcalcifications detection in mammographic images using texture coding
Ushizima et al. Statistical segmentation and porosity quantification of 3D x-ray microtomography
Kanagaraj et al. Pulmonary tumor detection by virtue of GLCM
Cardoso et al. Porosity features extraction based on image segmentation technique applying k-means clustering algorithm
Sparkman et al. A supervised learning approach for prediction of x-ray computed tomography data from ultrasonic testing data
Grandin et al. Implementation of automated 3D defect detection for low signal-to noise features in NDE data
Baravalle et al. Three‐dimensional multifractal analysis of trabecular bone under clinical computed tomography
Leitner et al. 3D-modeling of deformed halite hopper crystals by Object Based Image Analysis
Crandall et al. A New Workflow of X-Ray CT Image Processing and Data Analysis of Structural Features in Rock Using Open-Source Software
CN119478558B (en) SAM model-based scanning electron microscope mineral identification method and equipment
CN119723212B (en) Image data processing methods, apparatus, computer-readable storage media, and processors

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11866833

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 14119661

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 11866833

Country of ref document: EP

Kind code of ref document: A1