Baseline Setting Method

Abstract

The figure is fixed in a given position. If peaks are seen in the positive Y direction, the minimum value of the difference in height between the spectrum and the figure in the range where the figure is present on the X-axis. The minimum value and the height of the figure at the reference point are added. The figure is moved within a range containing the reference point, and the minimum value of the difference in height between the spectrum and the figure is added to the height of the figure at the reference point, at each point on the figure. A maximum value L(xi) of the calculated values is obtained, and the maximum value L(xi) is obtained as a baseline value at the X coordinate of the reference point.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to baseline setting methods, and more specifically, to an improvement of a baseline estimation method.

2. Description of the Related Art

In spectrum measurement by an infrared spectrophotometer or the like, the wavelength or wavenumber is indicated on the reference axis (X-axis), and the absorbance, reflectance, or the like is indicated on the Y-axis, in usual cases. Then, measured values are plotted with respect to points on the reference axis.

While the absorbance, reflectance, or the like is plotted when the wavelength is changed, a narrow mountain, known as a peak, appears in the spectrum. This occurs because the atom or molecule has the property of emitting or absorbing light at a specific wavelength.

However, the baseline of the spectrum depends on the characteristics of the instrument used to take the spectrum or the environment in which the spectrum measurement is being carried out. Sometimes, the occurrence of fluorescence and the like may also influence the baseline.

Therefore, the baseline setting for the spectrum and spectrum data correction based on the set baseline are essential techniques in spectrum data processing.

In particular, in infrared transmission spectrum measurement of a specimen having a rough surface or a specimen containing a pulverized inorganic compound, dispersion occurs on the surface of the specimen or in the specimen. This has a great influence on light with a short wavelength and can lower the transmission spectrum (% T) in the high wavenumber region. The ATR spectrum (% T) in the low wavenumber region of a specimen containing carbon black has a falling tendency.

Since these spectra have a large overall height, the peaks become relatively small, making it difficult to conduct a library search.

It is difficult to measure changes in the baseline directly, so that those changes in the baseline must be estimated.

For example, in Japanese Unexamined Patent Application Publication No. Hei-05-60614, a circle having a diameter not smaller than twice the full width at half maximum of a peak, an ellipse, or the like is moved in contact with but not intersecting the actually measured spectrum, and a part of the track is used as the baseline.

Although conventional methods related to baseline correction such as that described above have advantages, there are many problems in actual use, such as high computational load and a wide range of parameters to be specified, such as the shape and size of the figure.

A simple, low-computational-load baseline setting method that can be used for any shape of spectrum, rising toward the right, rising toward the left, rising in the middle, falling in the middle, or being wavy, has been desired.

SUMMARY OF THE INVENTION

In view of the related art, it is an object of the present invention to provide a simple, high-accuracy baseline setting method with low computational load.

If peaks are seen in the positive Y direction, a baseline can be specified by using a semicircle or a semi-ellipse, as described below.

A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value y_i identified with respect to a point x_i on the reference axis X, the baseline setting method including the steps of:

specifying a semicircle or semi-ellipse C_xi; centered on the reference axis X and an

X radius R, expressed by

x ² +a ² y ² =r ²

y≧0

f _Xi(x)={(r²−x²)^1/2}/a

if peaks of the measured curve are seen in the positive Y direction;

comparing individual points (x_i−r, f_xi(x_i−r)) to (x_i+r, f_xi(x_i+r)) of the semicircle or semi-ellipse C_Xi, centered at (x_i, 0) with the corresponding points (x_i−r, y_Xi−r) to (x_i+r, y_Xi+r) on the measured curve to calculate

$\begin{array}{c}{l}_{\mathrm{Xi}-r}=y_{\mathrm{Xi}-r}-f_{\mathrm{Xi}}\left(x_i-r\right)\\ ⋮\\ l_{\mathrm{Xi}+r}=y_{\mathrm{Xi}+r}-f_{\mathrm{Xi}}\left(x_i+r\right)\end{array}$

and to specify the minimum value of to 1_Xi−r to 1_Xi+r as 1_(Xi)min;

specifying the maximum value of

$\begin{array}{c}{l\left(x_i-r\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)\\ ⋮\\ {l\left(x_i+r\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)\end{array}$

obtained with respect to each of semicircles or semi-ellipses C_Xi−r to C_Xi+r centered at (x_i−r, 0) to (x_i+r, 0), as L_(xi);

setting a baseline point (x_i, L_(xi)) corresponding to a specific point (x_i, y_i) on the measured curve; and

connecting the baseline points corresponding to the individual points on the measured curve, to form a baseline.

If peaks are seen in the negative Y direction, a baseline can be specified as described below.

A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value y_i identified with respect to each point x_i on the reference axis X, the baseline setting method including the steps of:

specifying a semicircle or semi-ellipse C_xi centered on the reference axis X and an X radius, expressed by

x ² +a ² y ² =r ²

y≦0

f _Xi(x)={(r²−x²)^1/2}/a

if peaks of the measured curve are seen in the negative Y direction;

comparing individual points (x_i−r, f_Xi(x_i−r)) to (x_i+r, f_Xi(x_i+r)) of the semicircle or semi-ellipse C_Xi centered at (x_i, 0) with the corresponding points (x_i−r, y_Xi−r) to (x_i+r, y_Xi+r) on the measured curve to calculate

$\begin{array}{c}{l}_{\mathrm{Xi}-r}=f_{\mathrm{Xi}}\left(x_i-r\right)-y_{\mathrm{Xi}-r}\\ ⋮\\ l_{\mathrm{Xi}+r}=f_{\mathrm{Xi}}\left(x_i+r\right)-y_{\mathrm{Xi}+r}\end{array}$

and to specify the minimum value of 1_Xi−r to 1_Xi+r as 1_(Xi)min;

specifying the minimum value of

$f_{\mathrm{Xi}-r}\left(x_i\right)-1{\left(x_i-r\right)}_{\min }$ $⋮$ $f_{\mathrm{Xi}+r}\left(x_i\right)-1{\left(x_i+r\right)}_{\min }$

obtained with respect to each of semicircles or semi-ellipses C_Xi−r to C_Xi+r centered at (x_i−r, 0) to (x_i+r, 0), as L_(xi);

setting a baseline point (x_i, L_(xi)) corresponding to a specific point (x_i, y_i) on the measured curve; and

connecting the baseline points corresponding to the individual points on the measured curve, to form a baseline.

If peaks are seen in the positive Y direction, a baseline can be specified by using a quadratic curve, as described below.

A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value y_i identified with respect to each point x_i on the reference axis X, the baseline setting method including the steps of:

specifying a quadratic curve D_Xi centered on the reference axis X, expressed by

y=a(x−b)²+c

under the conditions:

b−M≦x≦b+M

M≧W

a<0

f(x)=a(x−b)²+c

where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks seen in the measured curve, if the peaks of the measured curve are seen in the positive Y direction;

comparing individual points (x_i−M, f_Xi(x_i−M)) to (x_i+M, f_Xi(x_i+M)) of the quadratic curve D_Xi having its vertex at (x_i, c) (x_i=b at the beginning of the measurement) with the corresponding points (x_i−M, y_Xi−M) to (x_i+M, y_Xi+M) on the measured curve to calculate

$1_{\mathrm{Xi}-M}=y_{\mathrm{Xi}-M}-f_{\mathrm{Xi}}\left(x_i-M\right)$ $⋮$ $1_{\mathrm{Xi}+M}=y_{\mathrm{Xi}+M}-f_{\mathrm{Xi}}\left(x_i+M\right)$

and to specify the minimum value of 1_Xi−M to 1_Xi+M as 1_(Xi)min;

specifying the maximum value of

$1{\left(x_i-M\right)}_{\min }+f_{\mathrm{Xi}-M}\left(x_i\right)$ $⋮$ $1{\left(x_i+M\right)}_{\min }+f_{\mathrm{Xi}+M}\left(x_i\right)$

obtained with respect to each of quadratic curves D_Xi−m to D_Xi+M having their vertices at (x_i−M, c) to (x_i+M, c), as L_(xi);

setting a baseline point (x_i, L_(xi) corresponding to a specific point (x_i, y_i) on the measured curve;

setting baseline points corresponding to the individual points on the measured curve by moving the vertex in the X direction; and

connecting the baseline points to form a baseline.

If peaks are seen in the negative Y direction, a baseline can be specified as described below.

A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value y_i identified with respect to each point x_i on the reference axis X, the baseline setting method including the steps of:

specifying a quadratic curve D_Xi centered on the reference axis X, expressed by

y=a(x−b)²+c

under the conditions

b−M≦x≦b+M

M≧W

a>0

f(x)=a(x−b)²+c

where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks seen in the measured curve, if the peaks of the measured curve are seen in the negative Y direction;

comparing individual points (x_i−M, f_Xi(x_i−M)) to (x_i+M, f_Xi(x_i+M)) of the quadratic curve D_Xi having its vertex at (x_i, c) (x_i=b at the beginning of the measurement) with the corresponding points (x_i−M, y_Xi−M) to (x_i+M, y_Xi+M) on the measured curve to calculate

$1_{\mathrm{Xi}-M}=f_{\mathrm{Xi}}\left(x_i-M\right)-y_{\mathrm{Xi}-M}$ $⋮$ $1_{\mathrm{Xi}+M}=f_{\mathrm{Xi}}\left(x_i+M\right)-y_{\mathrm{Xi}+M}$

and to specify the minimum value of 1_Xi−M to 1_Xi+M as 1_(xi)min;

specifying the minimum value of

$f_{\mathrm{Xi}-M}\left(x_i\right)-1{\left(x_i-M\right)}_{\min }$ $⋮$ $f_{\mathrm{Xi}+M}\left(x_i\right)-1{\left(x_i+M\right)}_{\min }$

obtained with respect to each of quadratic curves D_Xi−M to D_Xi+M with their vertices at (x_i−M, c) to (x_i+M, c), as L_(xi);

setting a baseline point (x_i, L_(xi) corresponding to a specific point (x_i, y_i) on the measured curve;

setting baseline points corresponding to the individual points on the measured curve by moving the vertex in the X direction; and

connecting the baseline points to form a baseline.

The track of a baseline is calculated by using a circle or ellipse, as described below.

First, in a plane having a reference axis X and a measured value axis Y extending orthogonally to the reference axis X, a semicircle or semi-ellipse C centered on the reference axis X, expressed by

x ² +a ² y ² =r ²

(f(x)={(r²−x²)^1/2}/a)

is specified, where y≧0 if peaks are seen in the positive Y direction or y≦0 if peaks are seen in the negative Y direction.

The values of “a” and “r” of the semicircle or semi-ellipse C are specified empirically.

If the X radius R of the semicircle or semi-ellipse is too small, the baseline is set at a high position in a peak and becomes too close to the spectrum in the peak.

Therefore, it is preferable to set the X radius R of the semicircle or semi-ellipse C to twice the full width at half maximum, W, or greater, where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measurement curve.

Optimum values of curvature “a” and radius “r” should be specified with the shape and inclination of the base of the measured curve and the peak width taken into account.

If peaks are seen in the positive Y direction, points (x_i−r, f_Xi(x_i−r)) to (x_i+r, f_Xi(x_i+r)) on the semicircle or semi-ellipse C_xi centered at point (x_i, 0) are compared with points (x_i−r, y_Xi−r) to (x_i+r, y_Xi+r) on the spectrum. Differences 1_Xi−r to 1_Xi+r corresponding to the individual points on the reference axis X are obtained as follows:

$1_{\mathrm{Xi}-r}=y_{\mathrm{Xi}-r}-f_{\mathrm{Xi}}\left(x_i-r\right)$ $⋮$ $1_{\mathrm{Xi}+r}=y_{\mathrm{Xi}+r}-f_{\mathrm{Xi}}\left(x_i+r\right)$

(If peaks are seen in the negative Y direction,

$1_{\mathrm{Xi}-r}=f_{\mathrm{Xi}}\left(x_i-r\right)-y_{\mathrm{Xi}-r}$ $⋮$ $1_{\mathrm{Xi}+r}=f_{\mathrm{Xi}}\left(x_i+r\right)-y_{\mathrm{Xi}+r}$

are calculated.) The minimum value among to is specified as 1_Xi(x_i)_min.

Then, even when the measured curve has its peak on point (x_i, 0) on the X-axis, corresponding to the center of the semicircle or semi-ellipse, if the X radius of the semicircle or semi-ellipse is greater than the maximum peak width appearing in the spectrum, 1_Xi(x_i)_min is calculated at an off-peak position on the X-axis.

Semicircles or semi-ellipses C_Xi−r to C_Xi+r are moved in the range of (x_i−r, 0) to (x_i+r, 0), 1_Xi−r(x_i)_min to 1_Xi+r(x_i)_min are obtained, by following the procedure described above, and

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)$

are calculated. (If peaks are seen in the negative Y direction,

$f_{\mathrm{Xi}-r}\left(x_i\right)-1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }$ $⋮$ $f_{\mathrm{Xi}+r}\left(x_i\right)-1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }$

are calculated.)

The maximum value (or the minimum value, if peaks are seen in the negative Y direction) among those values is specified as L_(xi).

Then, a point (x_i, L_(xi)) is set as a baseline point.

As has been described above, by specifying appropriate values of curvature “a” and radius “r” for f(x), a natural baseline can be set for any type of spectrum having a measured curve rising toward the right, rising toward the left, rising in the middle, falling in the middle, or being wavy.

The point (x_i, 0) is moved in the X direction, L_(x(i+1)) is calculated by following the procedure described above, and a baseline point corresponding to a specific point (x_(i+1), y_x(i+1) on the measured curve is set as (x_(i+1), L_(x(i+1))).

The baseline points corresponding to individual points on the measured curve are calculated in the same way.

By connecting these baseline points, the baseline of the measured curve can be specified.

The procedure for the baseline setting method using a quadratic curve is almost the same as the procedure for the baseline setting method using a semicircle or semi-ellipse. If peaks are seen in the positive Y direction, a quadratic curve expressed by

f(x)=a(x−b)²+c

a<0

is specified in the range of

b−M≦x≦b+M.

If peaks are seen in the negative Y direction, the quadratic curve is expressed by

f(x)=a(x−b)²+c

a>0.

First, a quadratic curve D centered at a point on the reference axis X, represented by

y=a(x−b)²+c

(f(x)=a(x−b)²+c)

is specified, under the following conditions:

b−M≦x≦b+M

M≧W

• • a<0 (a>0, if peaks are seen in the negative Y direction) where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measured curve. The value of “a” is specified empirically. The values of “b” and “c” are specified by the person performing the measurement as the initial position (b, c) of the vertex of the quadratic curve D.

The value of M should not be smaller than the full width at half maximum, W, of the peak having the greatest peak width among the plurality of peaks appearing in the measured curve. If the value of M is smaller, the baseline is specified at a high in a peak and becomes too close to the spectrum in the peak.

An optimum value of “a” should be specified empirically with the shape and inclination of the base of the measured curve and the peak width taken into account.

If peaks are seen in the positive Y direction, individual points (x_i−r, f_Xi(x_i−r)) to (x_i+r, f_Xi(x_i+r)) on the quadratic curve D_Xi having its vertex at point (x_i, c) (x_i=b at the beginning of measurement) are compared with points (x_i−r, y_Xi−r) to (x_i+r, y_Xi+r) on the spectrum. Differences 1_Xi−r to 1_Xi+r corresponding to the individual points on the reference axis X are obtained as follows:

$1_{\mathrm{Xi}-r}=y_{\mathrm{Xi}-r}-f_{\mathrm{Xi}}\left(x_i-r\right)$ $⋮$ $1_{\mathrm{Xi}+r}=y_{\mathrm{Xi}+r}-f_{\mathrm{Xi}}\left(x_i+r\right)$

(If peaks are seen in the negative Y direction,

$1_{\mathrm{Xi}-r}=f_{\mathrm{Xi}}\left(x_i-r\right)-y_{\mathrm{Xi}-r}$ $⋮$ $1_{\mathrm{Xi}+r}=f_{\mathrm{Xi}}\left(x_i+r\right)-y_{\mathrm{Xi}+r}$

are calculated.) The minimum value among 1_Xi−r to 1_Xi+r is specified as 1_Xi(x_i)_min.

Then, even when the measured curve has a peak at point (x_i, 0) on the X-axis, corresponding to the vertex of the quadratic curve, if the value of “a” is specified to set the baseline to an appropriate height with respect to the peak and if the width (=2M) of the quadratic curve is greater than the maximum peak width appearing on the spectrum, 1_Xi(x_i)_min is calculated at an off-peak position on the X-axis.

Quadratic curves D_Xi−r to D_Xi+r are moved in the range of (x_i−r, 0) to (x_i+r, 0), 1_Xi−r(x_i)_min to 1_Xi+r(x_i)_min are obtained by following the procedure as described above, and

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)$

are calculated. (If peaks are seen in the negative Y direction,

$f_{\mathrm{Xi}-r}\left(x_i\right)-1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }$ $⋮$ $f_{\mathrm{Xi}+r}\left(x_i\right)-1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }$

are calculated.)

The maximum value (or the minimum value, if peaks are seen in the negative Y direction) among those values is specified as L_(xi).

Then, a point (x_i, L_(xi) is set as a baseline point.

As has been described above, by specifying an appropriate value of “a” for f(x), a natural baseline can be set for any type of spectrum having a measured curve rising toward the right, rising toward the left, rising in the middle, falling in the middle, or being wavy.

The vertex is moved in the X direction, L_(x(i+1)) is calculated by following the procedure described above, and a baseline point corresponding to a specific point (x_(x+1), y_x(i+1)) on the measured curve is set as (x_(i+1), L_(x(i+1))).

The baseline points corresponding to the individual points on the measured curve are calculated in the same way.

By connecting these baseline points, the baseline of the measured curve can be specified.

As has been described above, if peaks are seen in the positive Y direction, by scanning a semicircle or semi-ellipse at any position on the X-axis about the reference point (x_i, 0), the minimum values 1_Xi−r(x_i)_min to 1_Xi+r(x_i)_min of difference in height between the spectrum and the figure at individual positions of the figure (having its center in the range of (x_i−r) to (x_i+r)) are added to the height f_Xi+r(x_i) to f_Xi+r(x_i) at the reference point (x_i, 0) of the individual positions of the figure, and

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)$

are obtained. The maximum value L_(xi) of the sum is obtained as a baseline value. The reference point is moved in the X direction, and the other baseline values are obtained by following the procedure as described above. By this simple method, a highly accurate baseline can be created automatically. Here, the figure should be moved in the X direction alone. Since complicated movements are not required, the computational load is low.

By adjusting the value of “r” to make the X radius R of the semicircle twice the full width at half maximum of a peak having the greatest peak width or greater, even if there is a peak at point (x_i, 0) on the X-axis, where the baseline point is going to be specified, the difference in height between the spectrum and the figure is calculated as the minimum value 1(x_i)_min at an off-peak X-axis position. Then, the figure is moved as described above, 1_Xi−r(x_i)_min to 1_Xi+r(x_i)_min are calculated, and the values are incorporated in the following calculation:

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right).$

Since the values of 1_Xi−r(x_i)_min to 1_Xi+r(x_i)_min are the differences in height between the spectrum and the figure at off-peak positions, they do not become too large, and L_(xi) is set to an appropriate value.

If the spectrum rises toward the right, rises toward the left, rises in the middle, falls in the middle, or is wavy, an appropriate baseline can be obtained by adjusting the values of curvature “a” and radius “r” in accordance with the shape and inclination of the entire spectrum and the width of the peak.

When a quadratic curve is used, a quadratic curve D having its vertex at point (x_i, c) (x_i=b at the beginning of the measurement) is formed under the following conditions

f(x)=a(x−b)²+c

b−M≦x≦b+M

M≧W

• • a<0 (if peaks are seen in the positive Y direction)where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measured curve, and a baseline is formed by following the procedure described above. By this simple method, a highly accurate baseline can be created automatically.

The values of “b” and “c” are specified as the initial position (b, c) of the quadratic curve.

The figure should be moved in the X direction alone. Since complicated movements are not required, the computational load is low.

If the figure such as a semicircle is moved along the spectrum, the complicated movement of the figure would require manual intervention. According to the present invention, the figure should be moved in the X direction alone. The baseline can be set through a very simple operation. After the initial parameters are specified, the baseline is set almost automatically without any other manual intervention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates the calculation of Y_mint+f_Xi(X_i) with a semicircle centered at a reference point (X_i, 0), when peaks are seen in the positive Y direction.

FIG. 2 illustrates the calculation of Y_min(t+1)+f_Xi+d(X_i) with a semicircle centered at point (X_i+d, 0) while the reference point is (X_i, 0), when peaks are seen in the positive Y direction.

FIG. 3 illustrates the calculation when the center of the semicircle is moved in the range of (X_i−r, 0) to (X_i+r, 0) while the reference point is (X_i, 0), when peaks are seen in the positive Y direction.

FIG. 4 illustrates the setting of baseline points by moving the reference point and using a semicircle at the reference point in accordance with a baseline setting method according to the present invention, when peaks are seen in the positive Y direction.

FIG. 5 illustrates the position of Y_min when the baseline is set by using a semicircle (y≦0), if peaks are seen in the negative Y direction.

FIG. 6 illustrates the calculation when the center of the semicircle (y≦0) is moved in the range of (X_i−r, 0) to (X_i+r, 0) while the reference point is (X_i, 0), if peaks are seen in the negative Y direction.

FIG. 7 shows the correspondence between positions of the figures and points plotted in FIG. 6, where the value of Y_min is subtracted from the height of each figure at the reference point (X_i, 0), and the calculated difference is plotted as height on x=X_i.

FIG. 8 shows an example of a baseline setting method according to the present invention, by using a semicircle.

FIG. 9 shows an example of baseline setting by using a line segment, and by adjusting the position of the line segment according to the method of the present invention.

FIG. 10 shows an example of the baseline setting method according to the present invention, by using a spectrum differing from that used in FIG. 8 and a semicircle.

FIG. 11 shows an example of baseline setting by using a spectrum differing from that used in FIG. 9 and a line segment, and by adjusting the position of the line segment according to the method of the present invention.

FIG. 12 shows an example of a baseline setting method according to the present invention, by using a semicircle, when peaks are seen in the negative Y direction.

FIG. 13 illustrates the calculation of Y_min in the baseline setting method according to the present invention, by using a semicircle with its circumference intersecting the measured curve.

FIG. 14 illustrates the calculation of Y_min1+f_Xi(X_i) with a quadratic curve having its vertex at a reference point (X_i, 0), when peaks are seen in the positive Y direction.

FIG. 15 illustrates the calculation of Y_min2+f_Xi+d(X_i) with a quadratic curve having its vertex at point (X_i+d, 0) while the reference point is (X_i, 0), when peaks are seen in the positive Y direction.

FIG. 16 illustrates the setting of a baseline point by moving the reference point and using a quadratic curve at the moved reference point in accordance with the baseline setting method according to the present invention, when peaks are seen in the positive Y direction.

FIG. 17 illustrates the calculation of f_Xi(X_i)+Y_min1 in an initial position (b, c) of the vertex of a quadratic curve D, the initial position being aligned with the x coordinate of a peak of a spectrum and the height of the peak, while the reference point is (Xi, 0).

FIG. 18 illustrates the calculation of f_Xi+r(X_i)+Y_min2 when the quadratic curve D is moved from the initial position (b, c) of the vertex by distance “r” in the X direction, the initial position being aligned with the x coordinate of the peak of the spectrum and the height of the peak, while the reference point is (X_i, 0).

DESCRIPTION OF THE PREFERRED EMBODIMENTS

With reference to FIGS. 1 to 4, a procedure for setting a baseline by using a semicircle at point (X, 0) when peaks are seen in the positive Y direction will be described.

With reference to FIGS. 5 to 7, a baseline point setting process according to the present invention, when the peaks are seen in the negative Y direction, will be described.

FIGS. 8 to 11 show baselines specified by using a line or semicircle. The baselines specified by using a semicircle or line will be compared while the adjustment of “a” is described.

FIG. 12 shows an example in which the baseline is specified by using a semicircle and the spectrum is baseline-corrected, when peaks are seen in the negative Y direction.

FIG. 13 illustrates a spectrum intersecting a semicircle having a large diameter.

With reference to FIGS. 14 to 16, a procedure for setting a baseline by using a quadratic curve with its vertex at point (X, 0) will be described.

FIGS. 17 and 18 show an example in which the initial position (b, c) of the vertex of a quadratic curve D is aligned with the x coordinate of a peak of the measured line and the height of the peak.

As a first embodiment of the present invention, a procedure for setting a baseline by using a semicircle when peaks are seen in the positive Y direction will be described with reference to FIGS. 1 to 4.

The operator first specifies a circle or ellipse with its center at point (X_i, 0).

The circle or ellipse is expressed by

X ² +a ² Y ² =r ²

(f(x)={(r²−x²)^1/2}/a)

and set by specifying the values of curvature “a” and radius “r”.

The values of “a” and “r” of the circle or ellipse are specified empirically. It is preferable to set the X radius R of the semicircle or semi-ellipse C to twice the full width at half maximum, W, or greater, where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measured curve, so that baseline points are set in an appropriate height.

From the spectrum S at each point in the horizontal axis region containing the figure, a measured value S at (x, 0) in the range of (X_i−r, 0) to (X_i+r, 0) is compared with {(r²−x²)^1/2}/a (y≧0) to calculate the difference Y between them:

Y=(S−{(r²−x²)^1/2}/a).

Then, the minimum value Y_mint of Y is calculated.

The calculated Y_mint is added to the height f_xi(X_i) of the figure at point (X_i, 0) to obtain (see FIG. 1)

Y_mint+f_xi(X_i).

The figure is moved by “d” parallel to the X-axis within a range including point (X_i, 0).

From the spectrum at each point in the horizontal axis region containing the figure, a measured value “S” in the range of (X_i+d−r, 0) to (X_i+d+r, 0) is compared with {(r²−x²)^1/2}/a (y≧0) to calculate the difference Y between them:

Y=(S−{(r²−x²)^1/2}/a).

Then, the minimum value Y_min(t+1) of Y is calculated.

The calculated Y_min(t+1) is added to the height f_Xi+d(X_i) of the figure at point (X_i, 0) to obtain (see FIG. 2)

Y_min(t+1)+f_xi+d(X_i).

By repeating this step while moving the figure within the range of (X−r, 0) to (X+r, 0), the following are calculated:

$Y_{\min 1}+f_{\mathrm{Xi}-r}\left(X_i\right)$ $⋮$ $Y_{\mathrm{mint}}+f_{\mathrm{Xi}}\left(X_i\right)$ $Y_{\min \left(t+1\right)}+f_{\mathrm{Xi}+d}\left(X_i\right)$ $⋮$ $Y_{\mathrm{minn}}+f_{\mathrm{xi}+r}\left(X_i\right).$

From the results, the maximum value P₁ is selected. A baseline point (X_i, P₁) at (X_i, 0) is specified (FIG. 3).

In FIG. 3, a peak on the measured curve is on the reference point (X_i, 0), but Y_min is set in an off-peak position in each of the five shown semicircles (Y_min values are calculated separately for the individual circles). By calculating

$Y_{\min 1}+f_{\mathrm{Xi}-r}\left(X_i\right)⋮Y_{\mathrm{mint}}+f_{\mathrm{Xi}}\left(X_i\right)Y_{\min \left(t+1\right)}+f_{\mathrm{Xi}+d}\left(X_i\right)⋮Y_{\mathrm{minn}}+f_{\mathrm{Xi}+r}\left(X_i\right)\}$

and selecting the maximum value P₁, the baseline point (X_i, P₁) is set to an appropriate height.

By moving the point (X_i, 0) in the X direction and repeating the same steps for point (X_(i+1), 0), the maximum value P₂ is set as a baseline point at point (X_(i+1), 0).

The same procedure is repeated to calculate baseline points P at (Xi, 0), (X_(i+1), 0), to (X_(i+n−1), 0). By connecting those points, a baseline BL is obtained (FIG. 4).

The present invention can also be used when peaks are seen in the negative Y direction. The procedure differs from the procedure when peaks are seen in the positive Y direction, and a supplementary description will be added.

When peaks are shown in the negative Y direction, a semicircle must be specified in the range of y≦0, as shown in FIG. 5.

Moreover, Y_min should be taken from a position where the absolute value of the distance between the circumference of the circle and the spectrum is the greatest. So, Ymin is set to f(x)−S, where f(x) is the value of height of the figure, and S is a measured value.

The height of the figure at (X_i, 0) is expressed as f(x_i), as shown in FIG. 6. While the figure is moved in the X direction, the following is obtained:

$f_{\mathrm{Xi}-r}\left(x_i\right)-Y_{\min 1}$ $⋮$ $f_{\mathrm{Xi}-1/2r}\left(x_i\right)-Y_{\mathrm{mini}}$ $⋮$ $f_{\mathrm{Xi}}\left(x_i\right)-Y_{\mathrm{mint}}$ $⋮$ $f_{\mathrm{Xi}+1/2r}\left(x_i\right)-Y_{\mathrm{minj}}$ $⋮$ $f_{\mathrm{Xi}+r}\left(x_i\right)-Y_{\mathrm{minn}}.$

In FIG. 6, a plurality of points are plotted, and the minimum value f_Xi+1/2r(x_i)−Y_minj becomes a baseline point.

FIG. 7 shows arrows indicating the correspondence between the points plotted in FIG. 6 and circles.

If the specimen emits light in absorbance measurement with an FT-IR or ultraviolet-visible photometer, the amount of light reaching the detector increases, decreasing the absorbance. In Raman scattering intensity measurement using a Raman spectrometer, if the specimen emits light, the spectrum can be raised.

In those cases, baseline correction using a circle, as shown in FIG. 8 or 10, can scale down the vertical axis, making it possible to detect a small peak that would not be detected on the previous scale of the vertical axis. Accordingly, the hit ratio in a library search can be improved.

For the comparison, FIGS. 9 and 11 show examples of baseline setting where a figure is adjusted according to the present invention, and a line segment is used.

(1) FIG. 8 shows an example using a circle, and FIG. 9 shows an example using a line segment, given for comparison. Although the baseline setting by using a line segment is not included in the scope of the present invention, it can be regarded as an example using an ellipse having a very small Y diameter because F(x), or the ellipse having a very small Y diameter, is close to a line segment.

The spectra shown in FIGS. 8(a) and 9(a) are rising toward the left. The scale of the Y-axis after baseline correction in FIG. 9(b) is smaller than that in FIG. 8(b).

That is because the baseline in FIG. 9(a) is parallel to the X-axis around a peak of the spectrum appearing in the range of 2800 to 3000 (cm⁻¹).

The baseline-corrected spectrum seen in the range of 2800 to 3000 (cm⁻¹) in FIG. 9(b) has a narrower shape than that in FIG. 8(b), and the corrected spectrum has a missing part. Therefore, it is understood that the baseline in FIG. 9(b) is inappropriate.

(2) FIG. 10 shows an example using a circle, and FIG. 11 shows an example using a line segment, given for comparison.

The baseline shown in FIG. 11(a) has a part parallel to the X-axis in the range of 500 to 1700 [cm⁻¹], and the shape of the baseline-corrected spectrum differs from that in FIG. 10(b).

Like the spectrum shown in FIG. 9(b), the baseline-corrected spectrum in FIG. 11(b) has a missing part. An ideal baseline provides a spectrum without a missing part, as shown in FIG. 10(a).

FIG. 12 shows an example of baseline setting according to the present invention, where a semicircle is used when peaks are seen in the negative Y direction.

The upper graph in FIG. 12 shows peaks in the negative Y direction and a spectrum declining toward the right. After baseline correction according to the present invention, the vertical axis representing transmittance (% T) is scaled down, and the peaks of the spectrum are emphasized.

FIG. 13 shows a case where the circumference of a circle intersects a measured curve.

When peaks are seen in the positive Y direction, if f(x)(={(r²−x²)^1/2}/a)(y≧0) is greater than a specific point y_i on the measured curve, 1(x_i)_min [=y_i−f(x)] becomes a negative value. Even in that case, the maximum value of

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)$

is selected as L_(xi), and a baseline point (x_i, L_(xi)) can be calculated normally in the same procedure.

As a second embodiment of the present invention, a procedure for setting a baseline by using a quadratic curve will be described with reference to FIGS. 14 to 16.

If peaks are seen in the positive Y direction, the operator first specifies a quadratic curve having its vertex at point (X_i, 0), under the following conditions:

f(x)=a(x−b)²+c

b−M≦x≦b+M

M≧W

a<0

where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measured curve.

The quadratic curve expressed by

y=a(x−b)²+c

(f(x)=a(x−b)²+c)

is set by specifying the values of “a”, “b”, and “c”.

The parameter “a” of the quadratic curve is set empirically, with the width of the peak taken into account. If peaks are seen in the positive Y direction, however, the value of “a” should be negative. The value of M is specified to satisfy W≦M, where W is the full width at half maximum of a peak having the greatest peak width among a plurality of peaks appearing in the measured curve. The range of the quadratic curve must be specified by setting the following:

b−M≦x≦b+M

The initial position (b, c) of the vertex of the quadratic curve is also specified.

In the second embodiment, “c” is set to zero (hereafter f(x)=a(x−b)²), but “c” is not confined to zero.

If peaks of the measured curve are seen in the negative Y direction, the value of “a” should be positive.

From the spectrum at each point in the horizontal axis region containing the figure, the measured value S in the range of (X_i−M, 0) to (X_i+M, 0) and the height a(x−b)² of the figure in that position are obtained to calculate the difference Y between them:

Y={S−a(x−b)²}.

The minimum value Y_min1 of Y is calculated next.

The height f_Xi(X_i) of the figure at point (X_i, 0) is added to Y_min1 to obtain (see FIG. 14)

Y_min1+f_Xi(X_i).

If peaks are seen in the negative Y direction,

Y={a(x−b)²−S}

is obtained, and the minimum value of Y should be set as Y_min1.

The figure is moved by “d” parallel to the X-axis within the range including point (X_i, 0).

From the spectrum at each point in the horizontal axis region including the figure, a measured value S in the range of (X_i+d−M, 0) to (X_i+d+M, 0) is compared with the height a(x−b)² of the figure in that position to calculate the difference Y between them:

Y={S−a(x−b)²}.

The minimum value of Y is calculated as Y_min2.

The height f_Xi+d(X_i) of the figure at point (X_i, 0) and Y_min2 are added to obtain (see FIG. 15)

Y_min2+f_Xi+d(X_i).

By repeating this step,

$Y_{\min 1}+f_{\mathrm{Xi}}\left(X_i\right)$ $Y_{\min 2}+f_{\mathrm{Xi}+d}\left(X_i\right)$ $⋮$ $Y_{\mathrm{minn}}+f_{\mathrm{Xi}+d\left(n-1\right)}\left(X_i\right)$

are calculated. The maximum value of those values is selected as P₁, and a baseline point (X_i, P₁) at point (X_i, 0) is specified.

The maximum value P₁ corresponds to the vertex (X_i, 0) of the figure in FIG. 14. That is,

Y_min1+f_xi(X_i)

becomes the maximum value P₁.

If peaks of the measured curve are seen in the negative Y direction, the minimum value of

$f_{\mathrm{Xi}}\left(X_i\right)-Y_{\min 1}$ $f_{\mathrm{Xi}+d}\left(X_i\right)-Y_{\min 2}$ $⋮$ $f_{\mathrm{Xi}+d\left(n-1\right)}\left(X_i\right)-Y_{\mathrm{minn}}$

should be set as P₁.

Point (X_i, 0) is moved in the X direction, the same operation is performed at point (X_(i+1), 0), and the maximum value P₂ is specified as a baseline point at point (X_(i+1), 0).

By repeating this operation, baseline points P at points (X_i, 0) through (X_(i+1), 0) to (X_(i+n−1), 0), 0) are calculated. By connecting those points, a baseline BL is obtained (FIG. 16). FIG. 16 shows three baseline points (X₁ to X₃), but more points are plotted in practice.

FIGS. 17 and 18 show examples in which the initial position (b, c) of the vertex of a quadratic curve D is aligned with the height of the peak and the x coordinate of the peak of the spectrum. The quadratic curve shown in FIG. 18 has been moved by distance “r” in the X direction from the position shown in FIG. 17.

If f(x) (=a(x−b)²+c) is greater than a specific point y_i on the measured curve, the value of 1(x_i)_min [=y_Xi−f(x)] becomes negative. Even in that case, the maximum value of

$1_{\mathrm{Xi}-r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}-r}\left(x_i\right)$ $⋮$ $1_{\mathrm{Xi}+r}{\left(x_i\right)}_{\min }+f_{\mathrm{Xi}+r}\left(x_i\right)$

is specified as L_(xi), and a baseline point (x_i, L_(xi)) can be calculated normally in the same procedure.

Claims

1. A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value yi identified with respect to a point xi on the reference axis X, the baseline setting method comprising the steps of: specifying a semicircle or semi-ellipse Cxi centered on the reference axis X and an X radius R, expressed by x2+a2y2=r2 y≧0fXi(x)={(r2−x2)1/2}/a

2. A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value yi identified with respect to each point xi on the reference axis X, the baseline setting method comprising the steps of: specifying a semicircle or semi-ellipse Cxi centered on the reference axis X and an X radius, expressed by x2+a2y2=r2 y≦0fXi(x)={(r2−x2)1/2}/a

3. The baseline setting method of claim 1, wherein the X radius R of the semicircle or semi-ellipse CXi is specified to be at least twice the full width at half maximum, W, of a peak having the greatest peak width among a plurality of peaks seen in the measured curve.

4. A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value yi identified with respect to each point x, on the reference axis X, the baseline setting method comprising the steps of: specifying a quadratic curve DXi centered on the reference axis X, expressed by y=a(x−b)2+c

5. A baseline setting method for a measured curve plotted on a reference axis X and a measured value axis Y extending in a direction orthogonal to the reference axis X, the measured curve having a single measured value yi identified with respect to each point xi on the reference axis X, the baseline setting method comprising the steps of: specifying a quadratic curve DXi centered on the reference axis X, expressed by y=a(x−b)2+c

6. The baseline setting method of claim 2, wherein the X radius R of the semicircle or semi-ellipse CXi is specified to be at least twice the full width at half maximum, W, of a peak having the greatest peak width among a plurality of peaks seen in the measured curve.