METHOD OF CORRECTING FOR AN AMPLITUDE CHANGE IN A SPECTROMETER

Abstract

A method of correcting for an amplitude change in a spectrometric instrument output includes: exposing a sample in a sample holder to electromagnetic radiation at a plurality of wavenumbers; detecting electromagnetic absorption intensities in the sample at the plurality of wavenumbers; providing to a computer device the detected electromagnetic absorption intensities indexed against wavenumber as spectral data; applying in the computer device a mathematical transform (Icorr) to the spectral data to correct for an amplitude change in the spectrometric instrument's output and calculated by determining a difference (Δ(SBZ)/) between first derivatives of a logarithmic transformation of spectral data (SB.sub.Z) from the zero material sample at two different wavenumber ranges (log.sub.10(SB.sub.Z(x.sub.1)).sup./ and log.sub.10(SB.sub.Z(x.sub.2)).sup./); and calculating the mathematical transform (Icorr) as a function inversely dependent on the determined difference (Δ(SBZ)/).

Claims

1. A method of correcting for an amplitude change in an output of a spectrometric instrument due to variations in an optical path length through a sample holder, the method comprising: exposing an unknown sample in the sample holder electromagnetic radiation at a plurality of wavenumbers; detecting, by the spectrometric instrument, electromagnetic absorption intensities in the unknown sample at the plurality of wavenumbers; providing, to a computer device, associated with the spectrometric instrument the detected electromagnetic absorption intensities indexed against wavenumber as spectral data; generating, by the computer device, absorbance values relative to a zero material sample from the spectral data; applying, by the computer device, a mathematical transform (I.sub.corr) to the generated absorbance values to correct for the amplitude change in the output of the spectrometric instrument; wherein the method further includes calculating, in the computer device, the mathematical transform (I.sub.corr) by determining a difference (Δ(SB.sub.Z).sup./) between a first derivative of a logarithmic transformation of spectral data (SB.sub.Z) from the zero material sample at a first wavenumber range (x.sub.1) (log.sub.10(SB.sub.Z)(x.sub.1)).sup./), and a first derivative of a logarithmic transformation of the spectral data (SB.sub.Z) from the zero material sample at a second wavenumber range (x.sub.2) (log.sub.10(SB.sub.Z)(x.sub.2)).sup./), and calculating, in the computer device, the mathematical transform (I.sub.corr) as a function inversely dependent on the determined difference (Δ(SB.sub.Z).sup./).

2. The method according to claim 1, wherein calculating the mathematical transform (I.sub.corr) includes correcting for humidity present in air in a light path between a radiation device and a detector of the spectrometric instrument by determining a humidity correction factor (corr.sub.hum) as a slope/intercept on selected humidity regions of the spectral data (SB.sub.Z) from the zero material sample; and calculating the mathematical transform (I.sub.corr) as a function inversely dependent on a sum of the determined difference (Δ(SB.sub.Z).sup./) and the humidity correction factor (corr.sub.hum).

3. The method according to claim 1, wherein the zero material sample comprises a nominally pure water sample.

4. The method according to claim 1, wherein the detected electromagnetic absorption intensities are detected by Fourier transform spectroscopy.

5. The method according to claim 1, wherein the spectral data (SB.sub.Z) from the zero material sample wavenumber standardized spectral data.

6. An apparatus, comprising: a spectrometric instrument configured to perform Fourier transform spectroscopy on a sample in a sample holder; and an associated computing device configured to receive spectral data from the spectrometric instrument and to apply thereto a mathematical transform stored in a memory accessible to the associated computing device, which transform corrects the spectral data for amplitude change in electromagnetic absorption intensities detected by the spectrometric instrument; wherein the associated computing device is programmed to operate to cause the apparatus to perform the method according to claim 1.

7. The method according to claim 2, wherein the mathematical transform (I.sub.corr) is calculated as a reciprocal of an inverse intensity correction (1/I.sub.corr) that is calculated according to equation (1),
(1/I.sub.corr)=c.sub.1.Math.Δ(SB.sub.Z).sup./+c.sub.2+(corr.sub.hum) (1) wherein, in equation (1), c.sub.1 is a first empirically-determined constant that is based solely on an absorption of water relative to air at a particular nominal optical path length L.sub.0, c.sub.2 is a second empirically-determined constant that is based on a contribution to a background due to dry air, and (corr.sub.hum) is calculated according to equation (2), $\begin{matrix} (c o r r_{h u m}) = - \log_{10} [\frac{P e a k_{1}}{P e a k_{2}}] c_{3} + c_{4} & (2) \end{matrix}$ wherein, in equation (2), Peak.sub.1 and Peak.sub.2 are respective intensities of an adjacent valley and peak in a fringe pattern in an infrared region of the spectral data (SB.sub.Z) from the zero material sample that is associated with a presence of humidity in the light path, c.sub.3 is a third empirically-determined constant that is based on an influence of humidity on the mathematical transform (Lon), and c.sub.4 is a fourth empirically-determined constant that indicates an offset between the intensities at Peak.sub.1 and Peak.sub.2 when no humidity is present in the light path.

8. A computing device, comprising: a memory storing a program of instructions; and a processor configured to execute the program of instructions to perform a method of correcting for an amplitude change in an output of a spectrometric instrument due to variations in an optical path length through a sample holder, the method including generating, by the computer device, absorbance values relative to a zero material sample from spectral data received from the spectrometric instrument, the spectral data indicating electromagnetic absorption intensities that are detected, by the spectrometric instrument, in an unknown sample in that is exposed to electromagnetic radiation in the sample holder at a plurality of wavenumbers, the spectral data indicating the electromagnetic absorption intensities indexed against wavenumber, and applying, by the computer device, a mathematical transform (I.sub.corr) to the generated absorbance values to correct for the amplitude change in the output of the spectrometric instrument, wherein the method further includes calculating, in the computer device, the mathematical transform (I.sub.corr) by determining a difference (Δ(SB.sub.Z).sup./) between a first derivative of a logarithmic transformation of spectral data (SB.sub.Z) from the zero material sample at a first wavenumber range (x.sub.1) (log.sub.10(SB.sub.Z(x.sub.1)).sup./), and a first derivative of a logarithmic transformation of the spectral data (SB.sub.Z) from the zero material sample at a second wavenumber range (x.sub.2) (log.sub.10(SB.sub.Z(x.sub.2)).sup./), and calculating, in the computer device, the mathematical transform (I.sub.corr) as a function inversely dependent on the determined difference (Δ(SB.sub.Z).sup./).

9. The computing device according to claim 8, wherein calculating the mathematical transform (I.sub.corr) includes Correcting for humidity present in air in a light path between a radiation device and a detector of the spectrometric instrument by determining a humidity correction factor (corr.sub.hum) as a slope/intercept on selected humidity regions of the spectral data (SB.sub.Z) from the zero material sample; and calculating the mathematical transform (I.sub.corr) as a function inversely dependent on a sum of the determined difference (Δ(SB.sub.Z).sup./) and the humidity correction factor (corr.sub.hum).

10. The computing device according to claim 8, wherein the zero material sample comprises a nominally pure water sample.

11. The computing device according to claim 8, wherein the electromagnetic absorption intensities are detected by Fourier transform spectroscopy.

12. The computing device according to claim 8, wherein the spectral data (SB.sub.Z) from the zero material sample comprises wavenumber standardized spectral data.

13. The computing device according to claim 9, wherein the mathematical transform (I.sub.corr) is calculated as a reciprocal of an inverse intensity correction (1/I.sub.corr) that is calculated according to equation (1),
(1/I.sub.corr)=c.sub.1.Math.Δ(SB.sub.Z).sup./+c.sub.2+(corr.sub.hum) (1) wherein, in equation (1), c.sub.1 is a first empirically-determined constant that is based solely on an absorption of water relative to air at a particular nominal optical path length L.sub.0, c.sub.2 is a second empirically-determined constant that is based on a contribution to a background due to dry air, and (corr.sub.hum) is calculated according to equation (2), $\begin{matrix} (c o r r_{h u m}) = - \log_{10} [\frac{P e a k_{1}}{P e a k_{2}}] c_{3} + c_{4} & (2) \end{matrix}$ wherein, in equation (2), Peak.sub.1 and Peak.sub.2 are respective intensities of an adjacent valley and peak in a fringe pattern in an infrared region of the spectral data (SB.sub.Z) from the zero material sample that is associated with a presence of humidity in the light path, c.sub.3 is a third empirically-determined constant that is based on an influence of humidity on the mathematical transform (I.sub.corr), and c.sub.4 is a fourth empirically-determined constant that indicates an offset between the intensities at Peak.sub.1 and Peak.sub.2 when no humidity is present in the light path.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0016] These and other advantages associated with the present invention will become apparent from a consideration of the following description of aspects of non-limiting exemplary embodiments of the present invention which is made with reference to the accompanying figures, of which:

[0017] FIG. 1 Illustrates schematically an embodiment of an inventive apparatus according to the present invention;

[0018] FIG. 2 Illustrates a schematic cross-sectional top view of the sample holder illustrated in FIG. 1;

[0019] FIG. 3 is a block diagram illustrating a method of determining a correction factor according to an embodiment of the present invention;

[0020] FIG. 4 is a graphical representation of first derivative single beam spectra of water;

[0021] FIG. 5 is a graphical representation of a region of the single beam spectra of water illustrated in FIG. 7 showing wavenumber positions for humidity correction;

[0022] FIG. 6 is a graphical representation of uncorrected absorbance spectra of milk;

[0023] FIG. 7 is a graphical representation of single beam spectra of water; and

[0024] FIG. 8 is a graphical representation of absorbance spectra of milk corrected according to the method of the present invention.

DETAILED DESCRIPTION

[0025] In the following, an embodiment of the inventive apparatus 100 will be described with reference to FIGS. 1 and 2 in the context of absorption spectroscopy. The apparatus 100 comprises a radiation device 200, an interferometric arrangement 300, a detector 400 and a measuring device 500. Also, a sample holder 600 for holding a sample to be analysed is arranged to be placed in the apparatus 100.

[0026] The radiation device 200 comprises a radiation source 210 which is arranged to emit polychromatic infrared radiation in the direction as indicated by the letter R in FIGS. 1 and 2.

[0027] The interferometric arrangement 300 comprises necessary equipment for implementing Fourier transform spectroscopy, as is well-known to a person skilled in the art. For example, the interferometric arrangement 300 comprises a collimator which collimates the infrared radiation and additional equipment comprised in an interferometer, for example optical components such as mirrors and lenses.

[0028] The detector 400 is arranged to detect incoming infrared radiation which is transmitted through the sample holder 600, see further below.

[0029] The measuring device 500 comprises a computer 510 which is connected to the detector 400 for collecting unprocessed data about the detected infrared radiation. By means of this connection, the measuring device 500 is arranged to determine a transmittance in a discrete number of channels positioned equidistantly along a wavenumber axis. The computer 510 comprises a processor for processing the collected data, suitable computing software, as well as additional equipment well-known to a person skilled in the art. Moreover, the computer 510 is arranged to store the collected data and the processed data in a memory. According to the present embodiment, a routine using Fourier transform algorithms is used in order to transform the unprocessed data from the detector 400 into data about the intensity as a function of the wavenumber. Moreover, the computer 510 is arranged to present the data graphically in terms of two-dimensional plots, see FIGS. 4-8 referred to below.

[0030] The radiation device 200, the interferometric arrangement 300, the detector 400 and the measuring device 500 will in the following be referred to as an FTIR spectrometer, or simply a spectrometer. Further below, a method for correcting intensity deviations (or amplitude changes) of this FTIR spectrometer will be described.

[0031] The sample holder 600 is placed between the interferometric arrangement 300 and the detector 400. Furthermore, the sample holder 600 is arranged to hold a liquid sample which is to be spectrally analysed by letting infrared radiation be transmitted through it. For instance, the liquid sample may be milk or wine. In the present embodiment, the liquid sample predominantly comprises water 610 which serves as a reference or so-called “zero” fluid and is used in order to perform corrections of cuvette pathlength deviations, see further below. The water sample 610 is placed in a cuvette 620 which is in part made out of calcium fluoride. The outer surface of the cuvette 620 is shaped as a rectangular parallelepiped. The cuvette 620 comprises inner walls 630, window elements 640, spacers 650, cavities 660 and a sample space 622 for holding the sample 610, see the cross-sectional top view in FIG. 2. It is clear that the inner walls 630 and the window elements 640 are transparent to the infrared radiation which is sent through the sample 610. It is noted that the spacers 650 do not need to be transparent. For example, the spacers 650 may be comprised out of a plastic. The volume of the sample space 622 may be varied by varying the extension of the spacers 650. Indeed, the spacers 650 create a pathlength of the cuvette 620. Furthermore, there is an inlet 670 for introducing the sample 610 into the sample space 622 and an outlet 680 for removing the sample 610 from the sample space 622. According some embodiments, the sample 610 is kept in motion during the measurement, flowing from the inlet 670 to the outlet 680 via the sample space 622, as indicated by the arrows in FIG. 2. In other embodiments, however, the sample 610 is kept stationary in the sample space 622 during the measurement, in which embodiments the inlet 670 and the outlet 680 may be omitted.

[0032] The distance covered by the infrared radiation in the sample space 622 is referred to as a pathlength. Since the radiation is transmitted through the sample 610 at right angles with respect to a side edge of the cuvette 620, in the direction R in FIG. 1 and FIG. 2, the pathlength L coincides with an inner length extension of the cuvette 620, between the window elements 640. If the cuvette 620 wears down, the pathlength L will change (increase).

[0033] In fact, since the window elements 640 making contact with the water sample 610 are made from calcium fluoride, they will be dissolved over time. During its lifetime, the cuvette 620 may also have been deteriorated by other chemicals. For example, the thickness T (see FIG. 2) of the window elements 640 will become smaller over time. Consequently, the pathlength L will increase over time, giving rise to pathlength deviations. In addition, it is noted that cuvettes placed in different apparatuses of the same type 100 by default have different pathlengths. For instance, differing pathlengths may have resulted from having dissolved the cuvettes to various degrees, even if the cuvettes have been substantially similar at some point in time. Moreover, the extension of the spacers 650 may vary between different cuvettes 620, thereby giving rise to varying pathlengths. Therefore, in order to make the characteristics of different apparatuses of the same type 100 more similar and the characteristics of a same apparatus 100 more stable over time, the variation in/of pathlengths need to be compensated for.

[0034] An exemplary embodiment of the method of correcting for intensity deviations in the apparatus 100 (here for example an FTIR spectrometer) according to the present invention is described below with reference to the block diagram of FIG. 3. As will be understood further below, a correction of a pathlength deviation also means a correction of an intensity deviation. According to the present exemplary embodiment the method utilizes the single beam spectrum SB.sub.Z of the zero liquid sample, which is here a nominally pure water sample (the water sample may contain small amounts of other components such as around 0.01% by volume detergent), for detecting deviations in cuvette pathlengths. After the spectrometer 100 has been corrected using measurements on the water sample, it may be used for measurements on other liquid samples, such as milk or wine in order to make quantitative determinations of components of interest in these samples in a manner well known in the art.

[0035] According to the present method an intensity correction, I.sub.corr, may be determined which corrects the absorbance values (A.sub.samp,L.sup.w) of a sample relative to the zero material (here water) measured in a sample holder of a pathlength L to those values (A.sub.samp,L0.sup.w) measured in a sample holder of standard pathlength L.sub.0. This may be described by:

A.sub.samp,L0.sup.w=I.sub.corr.Math.A.sub.samp,L.sup.w (1)

From this it follows that the intensity correction is a ratio between the two pathlengths as:

I.sub.corr=L.sub.0/L (2)

[0036] An aim of the present invention is to provide a method by which I.sub.corr may be determined without the use of a reference material having a known absorbance.

[0037] A logarithmic transformation of the single beam spectrum of the zero material SB.sub.Z, here water, converts the intensity values (y axis values) into absorbance units and the single beam spectrum SB.sub.Z may then be decomposed into its different components, as:

log.sub.10(SB.sub.Z)=log.sub.10(SB.sub.air)−A.sub.w,L.sup.air (3)

where SB.sub.air is the single beam spectrum of air and A.sub.w,L.sup.air is the absorbance of water relative to air at the pathlength L. From equation (3) it follows that:

log.sub.10(SB.sub.Z)=log.sub.10(SB.sub.air)−(L/L.sub.0).Math.A.sub.w,L0.sup.air (4)

where A.sub.w,L0.sup.air is the absorbance of water relative to air at the pathlength L.sub.0 From equation (2) the equation (4) may be re-written as:

log.sub.10(SB.sub.Z)=log.sub.10(SB.sub.air)−(1/I.sub.corr).Math.A.sub.w,L0.sup.air (5)

[0038] In order to reduce the effects of the background, typically to a negligible level, according to the present method a first derivative is calculated. The derivative of the background SB.sub.air varies slowly whilst the derivative of the water absorption SB.sub.Z will vary much more quickly, particularly at the rising and the trailing edges of the water absorption bands in the in the single beam spectrum of water SB.sub.Z, thus from equation (5):

log.sub.10(SB.sub.Z).sup./=log.sub.10(SB.sub.air).sup./−((1/I.sub.corr).Math.A.sub.w,L0.sup.air).sup./ (6)

By obtaining the derivatives at two different wavenumbers (or ranges), x.sub.1 and x.sub.2 lying close to one another on the x axis (wavenumber axis) it is possible to also neglect any slope in the background. It is preferable that the two different wavenumbers (or ranges), x.sub.1 and x.sub.2 are selected such that the slopes in the water spectrum relative to air at these points (or ranges) are very different. This may be expressed mathematically as:

log.sub.10(SB.sub.Z(x.sub.1)).sup./−log.sub.10(SB.sub.Z(x.sub.2)).sup./=log.sub.10(SB.sub.air(x.sub.1)).sup./−log.sub.10(SB.sub.air(x.sub.2)).sup./−(1/I.sub.corr).Math.(A.sub.w,L0.sup.air/(x.sub.1)−A.sub.w,L0.sup.air/(x.sub.2)) (7)

Which may be more simply re-written as:

Δ(SB.sub.Z).sup./=Δ(SB.sub.air).sup./−(1/I.sub.corr).Math.ΔA.sub.L0.sup./ (8)

where:

Δ(SB.sub.Z).sup./=log.sub.10(SB.sub.Z(x.sub.1)).sup./−log.sub.10(SB.sub.Z(x.sub.2)).sup./

is the difference in slope in the measured single beam zero spectrum of water

Δ(SB.sub.air)/=log.sub.10(SB.sub.air(x.sub.1)).sup./−log.sub.10(SB.sub.air(x.sub.2)).sup./

is the difference in slope in the background of the measured single beam zero spectrum of water

ΔA.sub.L0.sup./=(A.sub.w,L0.sup.air/(x.sub.1)−A.sub.w,L0.sup.air/(x.sub.2))

Is the difference in slope of the absorption spectrum of water relative to air at the nominal pathlength

[0039] Next, the background may be decomposed into a component due to dry air in the light path, SB.sub.air,dry, and optionally a component due to humidity in the air in the light path SB.sub.hum, thus:

Δ(SB.sub.air).sup./=Δ(SB.sub.air,dry).sup./+Δ(SB.sub.hum).sup./ (9)

From equation (8) this gives:

Δ(SB.sub.Z).sup./=Δ(SB.sub.air,dry).sup./+Δ(SB.sub.hum).sup./−(1/I.sub.corr).Math.ΔA.sub.L0.sup./ (10)

The intensity correct I.sub.corr may then be expressed as:

(1/I.sub.corr)=(Δ(SB.sub.Z).sup.//ΔA.sub.L0.sup./)+(Δ(SB.sub.air,dry).sup.//ΔA.sub.L0.sup./)+(Δ(SB.sub.hum).sup.//ΔA.sub.L0.sup./) (11)

or:

(1/I.sub.corr)=c.sub.1.Math.Δ(SB.sub.Z).sup./+c.sub.2+corr.sub.hum (12)

The values c.sub.1, c.sub.2 and a formula for the humidity correction (corr.sub.hum), when used, needs to be determined experimentally: [0040] the constant c.sub.1 is dependent solely on the absorption of water relative to air at the nominal pathlength L.sub.0; [0041] the constant c.sub.2 is the contribution to the background due to dry air, particularly but not essentially this may be a mean contribution over a population of apparatus 100; and [0042] the correction, corr.sub.hum, is a correction based on the amplitude of the humidity signature in the single beam spectrum, SB.sub.Z, of the zero material (here water) and in some cases where its effects are negligible, may be ignored.
According to the present invention corr.sub.hum is calculated as a slope/intercept on selected humidity regions, as:

[00001] $\begin{matrix} c o r r_{h u m} = - \log_{10} [\frac{P e a k_{1}}{P e a k_{2}}] c_{3} + c_{4} & (13) \end{matrix}$ [0043] the constant c.sub.3 describes the influence of humidity on the intensity correction I.sub.corr; and [0044] the constant c.sub.4 is the offset between the intensity at the two peak positions when no humidity is present. This approaches zero as Peak.sub.1 and Peak.sub.2 approach one another.

[0045] The humidity (from water in the gas phase, i.e. water vapour) produces an infrared spectrum with a fringe pattern. The ratio between a valley (Peak.sub.2) and an adjacent peak (Peak.sub.1) in this fringe pattern can be used as a measure for the amount of water vapour in the light path. The constants c.sub.1, c.sub.2, c.sub.3 and c.sub.4 are determined empirically by measuring the single beam zero spectrum SB.sub.z of water on one or more apparatus 100 of the same type using cuvettes with pathlengths covering an expected variation, such as 50 μm to 60 μm or from 37 μm to 44 μm. The effect on the spectrum of a known sample (e.g. milk or wine or glycerol or other chemical solution) is collected and the constants are adjusted until all spectra are identical. This task need only be performed once for a given apparatus type and used in the future for all apparatus of this type.)

[0046] In some embodiments the x axis (or wavenumber scale) of the single beam spectrum is standardized before the y axis (amplitude) correction is performed. This may be achieved in a manner that is well known in the art by applying a mathematical transform to the spectrum by which transform measured data is standardized along the x axis. In the present exemplary embodiment, the x axis standardization is based on the CO.sub.2 peak in the infrared range.

[0047] As is known, this x axis standardisation comprises standardising the wavenumber scale of an optical spectrum recorded by the apparatus 100 by providing an optical spectrum recorded by the apparatus 100 and comprising spectral patterns originating from constituents of atmospheric air in the light path in the apparatus 100; selecting a spectral pattern originating from constituents of atmospheric air in the apparatus 100, here CO.sub.2 in air; determining one or more wavenumber dependent position values associated with the selected spectral pattern; constructing a mathematical transform based on a difference between the determined value or values and a corresponding reference value or values of the selected spectral pattern and applying the mathematical transform to optical spectra subsequently recorded by the apparatus 100 to standardise the wavenumber scale.

[0048] Additionally the humidity correction corr.sub.hum may be employed in other methods of determining an intensity correction, such as for example the method disclosed in U.S. Pat. No. 9,874,515, to compensate for background effects.

[0049] The exemplary embodiment of the method according to the present invention will now be further described with reference to the block diagram illustrated in FIG. 3. The method (Box 700) comprises an exposure of the zero liquid having a known amount of at least one component (here water sample) 610 to polychromatic infrared radiation (Box 710) from the radiation device 200. The radiation is illustrated by wavy lines in FIG. 1 and FIG. 2. The detector 400 detects the incoming infrared radiation which has been transmitted through the interferometric arrangement 300, the water sample 610 as well as the cuvette 620, thereby determining (Box 720) the intensity levels for wavenumbers in the range between around 900 cm.sup.−1 and around 3500 cm.sup.−1, utilizing the measuring device 500. More specifically, the intensity levels for a discrete set of, typically equidistantly distributed, wavenumbers in this range are determined. The intensity data indexed against wavenumber data are stored in the memory of the computer 510 as the single beam spectrum of the zero material, SB.sub.Z. The computer 510 processes the stored data using a suitable mathematical transform to obtain (Box 730) log.sub.10-transformed intensity (or Absorbance) levels standardised along the x axis (the wavenumber scale).

[0050] The first derivative, log.sub.10(SB.sub.Z).sup./, of this absorbance spectrum is calculated (Box 740), at least in the region of the spectrum characterised by water absorption at around 1650 cm.sup.−1. This may be achieved in the computer 510 using the known Savitzky-Golay algorithm. A difference between two ranges in the first derivative, Δ(SB.sub.Z).sup./, is calculated (Box 740) where the two ranges are characterised by being close to the water band absorption at around 1650 cm.sup.−1, preferably on the higher wavenumber shoulder of the water band. In the present embodiment the ranges 1740-1746 cm.sup.−1 and 1844-1850 cm.sup.−1 are used. This is illustrated by the broken line constructions in FIG. 4 by which is illustrated a representative first derivative single beam spectra of water obtained by a plurality of instruments 100 of the same type (here four: labelled as Instrument 1, Instrument 2, Instrument 3 and Instrument 4 in the FIG. 4).

[0051] The humidity correction, corr.sub.hum, when used, is calculated (Box 750) according to equation (13) to compensate for the effect of environmental humidity on the intensity correction I.sub.corr. In the present embodiment the ranges 1832-1840 cm.sup.−1 (Peak.sub.1) and 1814-1822 cm.sup.−1 (Peak.sub.2) are used. This is illustrated by the broken line constructions in FIG. 5 by which is illustrated representative single beam spectra of water obtained by the same plurality (here four) of instruments 100 referenced above with respect to FIG. 4.

[0052] The intensity correction, I.sub.corr, is then calculated (Box 760) as the reciprocal of the inverse intensity correction, 1/I.sub.corr, calculated according to equation (12). This intensity correction, I.sub.corr, may then be applied (Box 770) to optical spectra subsequently recorded by the apparatus 100 to standardise the absorbance intensity scale (y axis).

[0053] An example of the application of this method will now be described in relation to the standardization of the output of a plurality of instruments 100 of the same type, here the four instruments referenced above with respect to FIG. 4 and FIG. 5, measuring on a same milk sample. It will be appreciated that in the following figures which are stated as illustrating measurements made by each of these four instruments the variations between instruments may be so small as to be visually represented as overlapping. Differences between these four instruments may then be better discerned from the numerical values presented in the accompanying tables, Table 1 to Table 3 inclusively.

[0054] The resulting absorbance of milk measured on the four different instruments 100 of the same type are plotted in FIG. 6 versus the wavenumbers. Each of the spectra in FIG. 6, is manifested as an interpolating curve in a two-dimensional plot, with magnitude of absorbance on the vertical axis and a corresponding wavenumber on the horizontal spectral axis. According to an alternatively graphical presentation, the plot may be a scatter plot. As can be seen the absorbance values at the same wavenumber differ slightly between instruments. This will lead to differences in the amount of components of the same milk sample being determined chemometrically, from these absorbance spectra in a known manner, in each of the instruments. When applying prediction models for fat, protein, lactose, total solids (TS) and solids non-fat (SNF) to uncorrected spectra of a set of, here fifteen, standard milk samples the following Root Mean Square Errors (RMSEPs) between the instruments and their common mean, are obtained—results in g/100 mL:

TABLE-US-00001 TABLE 1 Predicted milk components using corrected spectra for each instrument Fat Protein Lactose TS SNF Instrument 1 0.147 0.132 0.188 0.479 0.349 Instrument 2 0.134 0.141 0.195 0.490 0.352 Instrument 3 0.040 0.041 0.068 0.146 0.117 Instrument 4 0.320 0.313 0.450 1.115 0.818 Mean 0.190 0.185 0.265 0.659 0.482 Relative error* 5.30% 5.08% 5.28% 5.11% 5.15% *The mean error relative to the average concentration of the parameter.
This level of performance is not acceptable to a user requiring virtually identical results from different instruments.

[0055] Following the method according to the present invention a single beam spectrum of water (SB.sub.Z) is obtained on each of the four instruments (see FIG. 7) and a first derivative of the absorbance spectrum, for each of these single beam spectra, standardised in the wavenumber axis, is determined in the computer 510. These first derivatives are illustrated in FIG. 4, in which the two regions at which the difference Δ(SB.sub.Z).sup./ is to be calculated are also depicted (broken lines). It can be seen from FIG. 7 and FIG. 4 that the y axis values at each wavenumber are slightly different for each of the four instruments 100. The recorded single beam spectrum is processed by the computer 510 of the instrument 100 which generates the respective spectrum in order to determine the humidity correction, corr.sub.hum, for each instrument from the intensity values at the two regions illustrated by the broken lines in FIG. 5. In each instrument the associated computer 510 then calculates the intensity correction factor I.sub.corr, specific to each instrument of the four.

TABLE-US-00002 TABLE 2 Intensity correction factors for each instrument Instrument Intensity correction (/.sub.corr) 1 0.9979 2 0.9953 3 0.9740 4 0.8823

[0056] The correction factor I.sub.corr for each instrument can then be applied to the absorbance data collected by each instrument (for example the absorbance data for the same milk sample illustrated in FIG. 6) and the corrected (or y axis standardised) absorbance data for each instrument is brought to closer conformity. This is illustrated in FIG. 8 which shows the absorbance spectra of milk that are illustrated in FIG. 6 but each corrected with the appropriate intensity correction factor, I.sub.corr, for the instrument which records the associated milk spectrum (see Table 2).

[0057] When applying this correction to a set of, here fifteen, standard milk samples the following Root Mean Square Errors (RMSEPs) between the instruments and their common mean, are found when predicting fat, protein, lactose, total solids (TS) and solids non-fat (SNF)—results in g/100 mL:

TABLE-US-00003 TABLE 3 Predicted milk components using corrected spectra for each instrument Fat Protein Lactose TS SNF Instrument 1 0.007 0.011 0.006 0.013 0.010 Instrument 2 0.007 0.010 0.017 0.034 0.022 Instrument 3 0.016 0.014 0.008 0.041 0.021 Instrument 4 0.012 0.017 0.010 0.022 0.012 Mean 0.011 0.014 0.011 0.030 0.017 Relative error* 0.32% 0.39% 0.23% 0.24% 0.19% *The mean error relative to the average concentration of the parameter.
These errors are very low compared to the typical prediction error of 0.8% relative against the values of these components as determined by known wet chemistry methods.

[0058] It will be appreciated by the skilled artisan that the application of the method according to the present invention permits greater standardisation of the outputs between different instruments 100 of the same type as well as the standardisation of outputs of the same instrument 100 as the pathlength through the cuvette 620 changes.

[0059] As will also be appreciated by the skilled artisan, instead of expressing the spectral information about the electromagnetic radiation in terms of a wavenumber, one may instead use a wavelength or a frequency without departing from the invention as claimed.

METHOD OF CORRECTING FOR AN AMPLITUDE CHANGE IN A SPECTROMETER

Assignee

Inventors

Cpc classification

Classification Explorer

G01N2201/129

PHYSICS

Classification Explorer

G01N2201/1214

PHYSICS

Classification Explorer

G01N21/3577

PHYSICS

Classification Explorer

G01N33/04

PHYSICS

Classification Explorer

G01N2021/3595

PHYSICS

Classification Explorer

G01N21/274

PHYSICS

International classification

Classification Explorer

G01N21/27

PHYSICS

Classification Explorer

G01N21/3577

PHYSICS

Abstract

Claims

Description