Automation of ingredient-specific particle sizing employing raman chemical imaging

• US 8,374,801 B2
• Filed: 01/08/2010
• Issued: 02/12/2013
• Est. Priority Date: 01/09/2009
• Status: Active Grant

• Current Assignee: Chemimage Corporation
• Original Assignee: Chemimage Corporation
• Inventors: Fuhrman, Michael, Priore, Ryan, Olkhovyk, Oksana, Klueva, Oksana
• Primary Examiner(s): BHAT, ADITYA S
• Application Number:
  US12/684,495
• Publication Number:
  US 20100179770A1
• Time in Patent Office: 1,131 Days
• Field of Search: 702/28
• US Class Current: 702/28
• CPC Class Codes: G01J 3/02 Details
• CPC Class Codes: G01J 3/027 Control of working procedures of a spectrometer; Failure detection; Bandwidth calculation
• CPC Class Codes: G01J 3/44 Raman spectrometry; Scattering spectrometry ; Fluorescence spectrometry
• CPC Class Codes: G01N 15/0227 using imaging; using holography
• CPC Class Codes: G01N 21/65 Raman scattering

Abstract

A system and method for determining at least one geometric property of a particle in a sample. A sample is irradiated to thereby generate Raman scattered photons. These photons are collected to generate a Raman chemical image. A first threshold is applied wherein the first threshold is such that all particles in the sample are detected. A particle in the sample is selected and a second threshold is applied so that at least one geometric property of the selected particle can be determined. At least one spectrum representative of the selected particle is analyzed to determine whether or not it is a particle of interest. The step of determining a second threshold may be iterative and automated via software so that candidate second thresholds are applied until a satisfactory result is achieved.

Claims

1. A method comprising:

configuring a spectroscopic imaging device to;

(a) irradiate a sample comprising at least one unknown particle to produce Raman scattered photons, and (b) collect the Raman scattered photons to generate a Raman chemical image representative of the sample;

integrating the intensity within the spectral peak at each pixel of the Raman chemical image to generate a working image; and

configuring a processor to;

(c) apply a first threshold to the working image wherein the first threshold is such that all particles in the sample are detected;

(d) select one of the detected particles;

(e) apply a second threshold to the working image to determine at least one geometric property of the selected particle, wherein the second threshold is unique to the selected particle such that the geometric property can be determined;

(f) analyze at least one spectrum, representative of the selected particle, wherein analyzing further comprises comparing at least one spectrum representative of the selected particle to at least one reference spectrum representative of a particle of interest; and

if comparing results in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as a particle of interest, and

if comparing does not result in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as not a particle of interest.

2. The method of claim 1 further comprising repeating steps (d)-(f) for at least one other unknown particle present in the sample.

3. The method of claim 1 wherein the at least one geometric property comprises the size of the selected particle.

4. The method of claim 1 wherein the geometric property of the selected particle comprises at least one of:

an area, a perimeter, a feret diameter, a maximum chord length, a shape factor, an aspect ratio, and combinations thereof.

5. The method of claim 1 wherein the geometric property of the selected particle comprises at least one of:

particle size, morphology, spatial distribution, and combinations thereof.

6. The method of claim 1 further comprising fusing the working image with a brightfield image representative of the sample to generate a fused image representative of the sample.

7. The method of claim 6 further comprising analyzing the fused image to determine at least one geometric property of at least one unknown particle in the sample.

8. The method of claim 7 wherein the geometric property further comprises at least one of:

particle size, morphology, spatial distribution, and combinations thereof.

9. The method of claim 7 wherein the geometric property further comprises at least one of:

an area, a perimeter, a feret diameter, a maximum chord length, a shape factor, an aspect ratio, and combinations thereof.

10. The method of claim 1 wherein the second threshold comprises a fraction of the integrated peak intensity of the Raman spectrum corresponding to the selected particle.

11. The method of claim 10 wherein the fraction comprises one half.

12. The method of claim 1 wherein the second threshold comprises the integrated peak intensity of at least one Raman spectrum corresponding to at least one edge of the selected particle.

13. The method of claim 1 wherein the second threshold is determined by averaging the peak intensities of two or more Raman spectra corresponding to at least one edge or the selected particle.

14. The method of claim 1 wherein the second threshold is determining by a method comprising:

(a) applying a candidate second threshold to the working image;

assessing the effectiveness of the candidate second threshold; and

(i) if, based on the assessment, the candidate second threshold is effective, identifying the candidate second threshold as a second threshold unique to the selected particle, such that the geometric property can be determined, and

(ii) if based on the assessment, the candidate second threshold is not effective, repeating steps (a)-(b) for at least one other candidate second threshold.

15. The method of claim 14 wherein a candidate second threshold is determined to be effective when the integrated peak intensity of the selected particle is five standard deviations above the average background.

16. The method of claim 14 wherein the second threshold determination method is automated via software.

17. The method of claim 1 wherein comparing further comprises applying a chemometric technique to the working image.

18. The method of claim 17 wherein the chemometric technique further comprises at least one of:

cluster analysis, principal component analysis (PCA), Cosine Correlation Analysis (CCA), Euclidian distance analysis (EDA), multivariate curve resolution (MCR), band target entropy method (BTEM), Mahalanobis distance (MD), adaptive subspace detector (ASO), multivariate curve resolution (MCR), and combinations thereof.

19. The method of claim 1 wherein the method is automated via software.

20. The method of claim 1 wherein the spectrum representative of the selected particle further comprises an average spectrum.

21. A method comprising:

configuring a spectroscopic imaging device to;

(a) irradiate a sample comprising at least one unknown particle of interest to thereby produce interacted photons wherein the interacted photons are selected from the group consisting of;

Raman scattered by the sample, reflected by the sample, emitted by the sample, absorbed by the sample, and combinations thereof; and

configuring a processor to;

(b) collect the interacted photons to thereby generate a chemical image representative of the sample;

(c) apply a first threshold to the chemical image wherein the first threshold is such that all particles in the sample are detected;

(d) select one of the particles detected as a result of applying the first threshold;

(e) apply a second threshold to the chemical image to thereby determine at least one geometric property of the selected particle, wherein the second threshold is unique to the selected particle such that the geometric property can be determined; and

(f) analyze at least one spectrum representative of the selected particle, wherein analyzing further comprises comparing at least one spectrum representative of the selected particle to at least one reference spectrum representative of a particle of interest, and

if comparing results in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as a particle of interest; and

if comparing does not result in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as not a particle of interest.

22. The method of claim 21 wherein comparing further comprises applying a chemometric technique to the working image.

23. The method of claim 22 wherein the chemometric technique further comprises at least one of cluster analysis, principal component analysis (PCA), Cosine Correlation Analysis (CCA), Euclidian distance analysis (EDA), multivariate curve resolution (MCR), band target entropy method (BTEM), Mahalanobis distance (MD), adaptive subspace detector (ASD), multivariate curve resolution (MCR), and combinations thereof.

24. A method comprising:

configuring a spectroscopic imaging device to;

(a) irradiate a sample comprising at least one unknown particle to produce Raman scattered photons, and (b) collect the Raman scattered photons to generate a Raman chemical image representative of the sample;

integrating the intensity within the spectral peak at each pixel of the Raman chemical image corresponding to at least two edges of the particle to generate a working image; and

configuring a processor to;

(c) apply a first threshold to the working image wherein the first threshold is such that all particles in the sample are detected;

(d) select one of the detected particles;

(e) apply a second threshold to the working image, wherein the second threshold is unique to the selected particle such that the geometric property can be determined;

(f) analyze at least one spectrum representative of the selected particle, wherein analyzing further comprises comparing at least one spectrum representative of the selected particle to at least one reference spectrum representative of a particle of interest; and

if comparing results in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as a particle of interest, and

if comparing does not result in a match between the spectrum representative of the selected particle and the reference spectrum representative of the particle of interest, identifying the selected particle as not a particle of interest.

25. The method of claim 24 wherein comparing further comprises applying a chemometric technique to the working image.

26. The method of claim 25 wherein the chemometric technique further comprises at least one of:

cluster analysis, principal component analysis (PCA), Cosine Correlation Analysis (CCA), Euclidian distance analysis (EDA), multivariate curve resolution (MCR), band target entropy method (BTEM), Mahalanobis distance (MD), adaptive subspace detector (ASD), multivariate curve resolution (MCR), and combinations thereof.