Patent Grant
9677934
This invention relates to the field of optical emission spectroscopy, and in particular to spectral data processing to provide background-corrected portions of spectral emission data. It may be applied, for example, in inductively coupled plasma or microwave optical emission spectrometry.
Inductively coupled plasma and microwave induced plasma optical emission spectroscopy (ICP-OES and MIP-OES respectively) are analytical techniques for determining the concentration of elements in a sample. A sample solution is injected into a plasma source, the temperature of which vaporises the sample, destroys any chemical bonds, ionises atoms and causes electronic excitation of the atoms and ions. A plasma spectrum consists of a continuum background and discrete spectral lines at wavelengths characteristic of any element contained in the sample. The intensity of any given spectral line is proportional to the concentration of the element in the sample and quantitative assessment of element concentration is then based upon a measurement of the intensity of light at a wavelength which is characteristic of the element.
A problem which exists in this and other analytical techniques involving emission spectra is that a continuum background can distort the shape of the spectral lines making an accurate assessment of element concentration difficult. Typical conventional techniques for providing some form of background correction are discussed in “ICP in Analytical Atomic Spectrometry”, Montaser & Golightly, VCH, 1987, section 6.6, and in “Median Filtering for Removal of Low-Frequency Background Drift”, Alvin W. Moore, James W. Jorgenson, Analytical Chemistry 1993, 65, 188-191. The simplest and a relatively unreliable technique is for the user to guess the wavelength value of one or two background points and use these wavelength values to apply a background correction in the determination of a concentration estimate. This initial choice of wavelength is typically used for all subsequent analyses of samples. Background correction is achieved at any particular wavelength of interest by linear interpolation between the background points. However this technique is time consuming and requires a relatively skilled technician to select suitable background points.
Similar problems exist in infra-red, gas chromatography, liquid chromatography and ultraviolet analytical instruments where changes in conditions affect the measurements made by the instrument.
Against this background the present invention has been made.
The present invention provides in a first aspect a method for deriving a background-corrected portion of a measured emission spectrum comprising the steps of identifying two or more background correction points from the portion of the measured emission spectrum, deriving a background correction function fitted to the identified background correction points, and applying the background correction function to the portion of the measured emission spectrum so as to produce a background-corrected portion of the emission spectrum; wherein the portion of measured emission spectrum comprises measured data points which consist of a series of n intensity values In, at discrete wavelengths or values corresponding to wavelength, λn, the series extending from a first measured data point to a last measured data point, and the background correction points are identified from the measured data points by the following steps:
identifying a first background correction point as the first measured data point;
calculating the gradients Gm of straight lines between the background correction point just identified and each subsequent measured data point in the series;
identifying a next background correction point as the measured data point which lies on the straight line which has the minimum gradient of all the calculated gradients Gm;
if the background correction point just identified does not lie at the last measured data point, repeat the procedure from step (2) until it does.
The method of the invention is applied to data points within a portion of the measured emission spectrum and the following description considers only those data points. Each portion of the measured emission spectrum comprises measured data points which consist of a series of n intensity values In, at discrete wavelengths or values corresponding to wavelength, λn, the series extending from a first measured data point to a last measured data point. The number of points, n, within each portion may differ from one portion to another. The measured emission spectrum may be in a form in which the values corresponding to wavelength are stored and manipulated as wavelengths, or as frequency, distance, pixel number, memory bin number, or some other numerical value.
The portion of the measured emission spectrum may be found manually by the user, or automatically. If only a single spectral line needs to be background corrected, selecting the first and last data points manually can give improved results for that spectral line. In this case, first and last data points for the portion of the emission spectrum are selected by the user and the method of the invention determines a background correction function from the selected portion of the emission spectrum. Generally however, preferably the portion of the measured emission spectrum is found using the following automated procedure.
The measured emission spectrum is analysed to find all peaks in intensity which lie above a predetermined intensity level and which comprise a minimum number of measured data points thereby being a predetermined spectral width. Some or all of the peaks thus found may then be chosen. Preferably peaks are chosen on the basis of where in the spectrum of measured wavelengths or values corresponding to wavelength they lie, so that peaks corresponding to chosen elements may be selected. Each such peak selected comprises a number of measured data points, and these together with a number of adjacent measured data points either side of the peak is selected to be a portion of the measured emission spectrum. The number of adjacent measured data points either side of the peak may be chosen on the basis of a fixed predetermined number, or by analysis of the intensities of the adjacent points, so as to select additional points which include signal minima, at least one minima at each side of the peak. Still more preferably the number of adjacent measured data points either side of the peak is chosen so as to include at least one signal minima at each side of the peak, and one or more data points beyond the signal minima. Hence the portion of the emission spectrum is identified from a larger set of measured spectral data by a method of peak detection.
Having selected a portion of the measured emission spectrum, the method starts by considering a first point in the portion of the measured emission spectrum, and then subsequent points, until reaching a last data point in the portion of the measured emission spectrum. The first data point may be the point corresponding to the lowest wavelength or value corresponding to wavelength in the portion of measured emission spectrum, in which case the last data point is the point corresponding to the highest wavelength or value corresponding to wavelength in the portion of measured emission spectrum. Conversely, the method of the invention may be worked by taking the first data point to be the point corresponding to the highest wavelength or value corresponding to wavelength in the portion of measured emission spectrum, in which case the last data point is the point corresponding to the lowest wavelength or value corresponding to wavelength in the portion of measured emission spectrum.
The background correction points are identified as follows. The first background correction point, BC1 is identified as the first measured data point. The next background correction point BC2 is determined to be a data point in the selected portion of measured emission spectrum at which a minimum gradient occurs, the gradient being the gradient of a straight line between the first background correction point and any other subsequent data point in the selected portion of measured emission spectrum.
If the background correction point BC2 is not the last data point in the selected portion of measured emission spectrum, one or more additional background correction points are then determined. A next background correction point BC3 is determined as a data point in the selected portion of measured emission spectrum at which a minimum gradient occurs, the gradient being the gradient of the line between the second background correction point BC2 and any other subsequent data point in the selected portion of measured emission spectrum. Subsequent data points are data points which lie between the second background correction point and the last point; data points preceding the last background correction point BC2 are not considered in this step. If the background correction point BC3 thus found is not the last point in the selected region, one or more additional background correction points are then determined in the same way as described above for finding BC3, where the gradient is taken of straight lines between the background correction point just identified and each subsequent measured data point in the series.
Hence the method of deriving background correction points may be described as follows, acting upon n points Pn, which make up the selected portion of measured emission spectrum, the points consisting of intensity values In at discrete wavelengths or values corresponding to wavelength λn:
Step (a): BC1=: P1
Step (b): i=: 1
Step (c): For each subsequent data point Pm, where m=i+1 . . . n, calculate the gradient Gm of the line between the previous background correction point BCi and Pm, where
and (Im,λm) are the intensity and value corresponding to wavelength for the point Pm, and (Iiλi) are the intensity and value corresponding to wavelength for the point BCi
Step (d) i=: i+1
Step (e): BCi=Pk, where Gk is the minimum gradient
Step (f): If Bi is not the last data point, i.e. BCi≠Pn, repeat from step (c); otherwise terminate.
The background correction points thus identified are then used in a background correction function which is applied to the portion of the measured emission spectrum to produce a background-corrected portion of the emission spectrum. The background correction function is preferably a mathematical function which is fitted to the background correction points.
The mathematical function may be a curve or it may be a straight line. In a preferred embodiment the background correction function is a combination of straight line functions taken between the each successive pair of background correction points. In this case, the background correction function is applied to any point Pk lying between a pair of background correction points by subtracting a background Bk from Ik, the intensity at point Pk. Bk is derived by linear interpolation of the straight line function taken between the pair of background correction points which lie either side of point Pk. This linear interpolation technique amounts to fitting a linear background correction function between adjacent background points which bound a range of wavelengths or values corresponding to wavelengths of interest. In this case it is assumed that an emission spectrum consists only of additive features and that over the width of a spectral line, the continuum background is modeled adequately by a straight line.
Where the background correction function is a curve, preferably the curve is fitted to three or more background correction points. The curve may be a spline fit, for example. Where the background correction function is a curve, the background correction function is applied to any point Pk lying within the wavelength or values corresponding to wavelength bounds of the fitted curve by subtracting a background Bk from Ik, Ik being the intensity at point Pk, and Bk being derived from the intensity ordinate of the fitted curve at the wavelength or value corresponding to wavelength λk.
Preferably all the points Pn which lie within the portion of the measured emission spectrum are thus background-corrected producing a background-corrected portion of the emission spectrum, which may then be further processed.
Preferably all selected portions of the emission spectrum are subjected to the above background correction method. The method may be applied to just the first detected spectrum in a series of analyses or to each subsequent spectrum as it is produced.
The method can be refined by allowing “roll-up” at the first and last data points in a portion of the emission spectrum to ensure that at least one intermediate background correction point is identified. The reason this may be required is that some spectral lines of interest lie beside molecular emission bands, and a graph of wavelength or values corresponding to wavelength against intensity shows a slight dip at either end. In this case only two background correction points are found by the method of the invention. Then, use of a linear interpolation function as the background correction function would effectively create a single straight line drawn from one edge of the portion of the measured emission spectrum to the other, giving an inaccurate estimate of the background, due to the presence of the molecular emission bands. The roll-up feature guarantees that at least one intermediate data point of the emission spectrum is used as a background correction point, eliminating this problem.
The roll-up method is used if the method of finding the background correction points only finds the first and last data points in the portion of measured emission spectrum. In this case, the method of finding the background correction points is followed again using a subset of the portion of measured spectral data, the subset being all points from and including the second data point to the next-to-last data point. Additional background correction points are then found and added to the original two background correction points found.
An alternative to the roll-up method is to apply the original method twice, once from the top of the peak to the first data point and once from the top of the peak to the last data point, to produce two sets of background correction points. The first background correction point in each set (corresponding to the data point on the peak top) is discarded, and the two sets of background correction points are then combined to provide a single set of background correction points. The background correction function is then derived from the background correction points in the single set.
Hence the invention further provides a method for deriving a background-corrected portion of a measured emission spectrum wherein if the background correction points identified consist only of the first and last data points in the portion of measured emission spectrum, the method further comprises:
The method of the present invention provides background-corrected portions of a measured emission spectrum which allow a good estimate of the true intensity value at a particular wavelength or value corresponding to wavelength of interest to be obtained. The method is rapid, robust, may be fully automated, and is less susceptible to noise. The method provides a means of automatically identifying background points to allow a background correction to be applied to a portion of an emission spectrum. The method may be applied to data measured in real time, or it may be applied to data gathered previously. It may therefore be applied to archived data to obtain improved estimates of spectral intensity from data recorded before the invention was made. It may also be applied to data acquired by instrumentation at locations remote from the location at which the method of the invention is performed. Accordingly independent aspects of the invention do not include the step of measuring the spectral data.
Further details of the means by which the invention provides these, and other, advantages are given in the description below.
The invention is preferably at least partially implemented on a computer, especially the data processing steps. Accordingly, the invention may utilise a computer readable medium carrying a computer program, the program having modules of program code for carrying out the method of the invention. A computer may be used to read the computer readable medium and perform some or all the steps of the method of the invention by using the modules of program code.
The method preferably further comprises outputting the background-corrected portion of the spectrum. Correspondingly, the system preferably further comprises an outputting device for outputting the background-corrected portion of the spectrum. The outputting device may comprise an electronic display device (e.g. VDU screen) and/or a printer, the outputting device preferably being under the control of the computer that is used to implement the invention. In this way, a tangible (user readable) form of the background-corrected portion of the spectrum is produced, e.g. on a display or on paper.
In a further aspect of the present invention, there is provided a system for determining a background-corrected portion of a measured emission spectrum, comprising a computer having an input for receiving a portion of measured emission spectrum data, and an output for outputting the background-corrected portion of a measured emission spectrum derived by the computer from the portion of measured emission spectrum data; wherein the computer is programmed with a program comprising one or more modules of program code for identifying two or more background correction points from the portion of the measured emission spectrum; one or more modules of program code for deriving a background correction function fitted to the identified background correction points, and one or more modules of program code for applying the background correction function to the portion of the measured emission spectrum so as to produce a background-corrected portion of the emission spectrum; wherein the portion of measured emission spectrum comprises measured data points which consist of a series of n intensity values In, at discrete wavelengths or values corresponding to wavelength λn, the series extending from a first measured data point to a last measured data point, and the background correction points are identified from the measured data points by the following steps:
While shown and described herein as a method and system for determining a background-corrected portion of a measured emission spectrum, it is understood that aspects of the invention further provide various alternative embodiments. For example, in one embodiment, the invention provides a computer program fixed in at least one computer-readable medium, which when executed, enables a computer system to determine a background-corrected portion of a measured emission spectrum. To this extent, the computer-readable medium includes program code which implements some or all of a process described herein. It is understood that the term “computer-readable medium” comprises one or more of any type of tangible medium of expression, now known or later developed, from which a copy of the program code can be perceived, reproduced, or otherwise communicated by a computing device. For example, the computer-readable medium can comprise: one or more portable storage articles of manufacture; one or more memory/storage components of a computing device; paper; and/or the like.
In order to further understand the invention, embodiments will now be described in more detail with reference to the accompanying drawings. The embodiments described are examples and for illustration only and are not intended to, and do not, limit the scope of the invention.
Spectrometer 200 comprises a lens system and polychromator 210, a detector 220, a data processing unit 230 and a spectrometer controller 240. The lens system captures the portion of optical emission intercepted by the spectrometer 200 and directs it onto the polychromator. The polychromator disperses the optical emission as a function of the wavelength of the photons. Dispersed photons are intercepted by detector 220 which comprises an array of light sensitive detection elements. Information on the intensity at different wavelengths is received from the detector 220 by data processing unit 230. Spectrometer 200 is controlled by spectrometer controller 240.
Optical emission source 100 and spectrometer 200 are controlled by a computer 300. Computer 300 is in communication with source control system 120 and spectrometer controller 240. Computer 300 also receives information from data processing unit 230.
System 1 is configured for determining a background-corrected portion of a measured emission spectrum in accordance with the present invention. Computer 300 comprises a computer-readable medium which includes program code which when performed by computer 300 implements the method of the invention.
A background correction function is derived, fitted to the identified background correction points, the background correction function being in this example a combination of straight line functions, A and B, taken between the each successive pair of background correction points, as shown in
An alternative background correction function could be used, in which the background correction function is a curve fitted to the three background correction points. The curve may be a spline fit, for example. In this case, the background correction function is applied to any point Pk lying within the wavelength or values corresponding to wavelength bounds of the fitted curve by subtracting a background Bk from Ik, Ik being the intensity at point Pk, and Bk being derived from the intensity ordinate of the fitted curve at the wavelength or value corresponding to wavelength λk.
As used herein, including in the claims, unless the context indicates otherwise, singular forms of the terms herein are to be construed as including the plural form and vice versa.
Throughout the description and claims of this specification, the words “comprise”, “including”, “having” and “contain” and variations of the words, for example “comprising” and “comprises” etc, mean “including but not limited to”, and are not intended to (and do not) exclude other components.
It will be appreciated that variations to the foregoing embodiments of the invention can be made while still falling within the scope of the invention. Each feature disclosed in this specification, unless stated otherwise, may be replaced by alternative features serving the same, equivalent or similar purpose. Thus, unless stated otherwise, each feature disclosed is one example only of a generic series of equivalent or similar features.
The use of any and all examples, or exemplary language (“for instance”, “such as”, “for example” and like language) provided herein, is intended merely to better illustrate the invention and does not indicate a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention.
Any steps described in this specification may be performed in any order or simultaneously unless stated or the context requires otherwise.
All of the features disclosed in this specification may be combined in any combination, except combinations where at least some of such features and/or steps are mutually exclusive. In particular, the preferred features of the invention are applicable to all aspects of the invention and may be used in any combination. Likewise, features described in non-essential combinations may be used separately (not in combination).
| Number | Date | Country | Kind |
|---|---|---|---|
| 1315284.8 | Aug 2013 | GB | national |
| Number | Name | Date | Kind |
|---|---|---|---|
| 5121337 | Brown | Jun 1992 | A |
| 5789754 | Cathey et al. | Aug 1998 | A |
| 6564114 | Toprac et al. | May 2003 | B1 |
| 7468790 | Sogan et al. | Dec 2008 | B2 |
| 7768639 | Chan et al. | Aug 2010 | B1 |
| 20040114139 | Florek | Jun 2004 | A1 |
| 20050283058 | Choo-Smith | Dec 2005 | A1 |
| 20090014635 | Neitsch | Jan 2009 | A1 |
| 20090086205 | Grun et al. | Apr 2009 | A1 |
| Number | Date | Country |
|---|---|---|
| 101692045 | Apr 2010 | CN |
| 59221626 | Dec 1984 | JP |
| Number | Date | Country | |
|---|---|---|---|
| 20150062577 A1 | Mar 2015 | US |
