JP2003235819A - Signal processing method and pulse wave signal processing method - Google Patents
Signal processing method and pulse wave signal processing methodInfo
- Publication number
- JP2003235819A JP2003235819A JP2002361800A JP2002361800A JP2003235819A JP 2003235819 A JP2003235819 A JP 2003235819A JP 2002361800 A JP2002361800 A JP 2002361800A JP 2002361800 A JP2002361800 A JP 2002361800A JP 2003235819 A JP2003235819 A JP 2003235819A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- spectrum
- calculated
- ratio
- processing method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、信号処理方法に関
し、特にパルスフォトメトリによる酸素飽和度測定等に
おいて用いることができる、脈波のノイズ除去に関す
る。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a signal processing method and, more particularly, to pulse wave noise removal that can be used in oxygen saturation measurement by pulse photometry.
【0002】動脈血の酸素飽和度を非観血的に連続測定
するには、従来からパルスオキシメータが用いられてい
る。このパルスオキシメータでは、プローブを被験者の
指先や耳朶に装着し、プローブから赤と赤外の異なる波
長の光を生体に時分割に照射して、異なる2波長の透過
光または反射光から得られる吸光度の脈動成分の比Φか
ら酸素飽和度Sを測定するものである。赤色光には、例
えば660nmの基準波長が用いられると共に、赤外光
には、例えば940nmの波長が用いられ、プローブ内
にはこれらの波長を発生する2つの発光ダイオードと受
光用の1つのフォトダイオードが内蔵されている。い
ま、赤外光の波長の吸光度の脈動成分をS1、赤色光の
波長の吸光度の脈動成分をS2とすると、異なる2波長
の吸光度の比Φは、次式で与えられる。
Φ=S2/S1 (1)
S=f(Φ) (2)A pulse oximeter has been conventionally used for non-invasive continuous measurement of oxygen saturation of arterial blood. In this pulse oximeter, a probe is attached to a fingertip or an earlobe of a subject, and light of different wavelengths of red and infrared is radiated from a probe to a living body in a time-sharing manner to obtain transmitted or reflected light of two different wavelengths The oxygen saturation S is measured from the ratio Φ of the pulsating component of the absorbance. A reference wavelength of, for example, 660 nm is used for red light, and a wavelength of, for example, 940 nm is used for infrared light. Two light emitting diodes that generate these wavelengths and one photo detector for receiving light are used in the probe. It has a built-in diode. Now, assuming that the pulsation component of the absorbance of the infrared light wavelength is S1 and the pulsation component of the absorbance of the red light wavelength is S2, the ratio Φ of the two different absorbances is given by the following equation. Φ = S2 / S1 (1) S = f (Φ) (2)
【0003】このようなパルスオキシメータでは、測定
中に患者に体動等が起きると、プローブで検出される脈
波にノイズ成分が混入してしまい、正確に酸素飽和度S
を測定できなくなる。そこで、このようなノイズの影響
を除去する試みが従来から種々なされている。In such a pulse oximeter, when a patient moves during measurement, a noise component is mixed in the pulse wave detected by the probe, and the oxygen saturation S is accurately measured.
Cannot be measured. Therefore, various attempts have conventionally been made to remove the influence of such noise.
【0004】例えば、酸素飽和度を測定するに当り、測
定された脈波信号S1、S2には、信号成分s1、s2
とノイズ成分n1、n2が含まれると考え、信号成分s
1とノイズ成分n1の相関が最小になるようなs1とs
2の比を求めるように構成した信号処理方法が提案され
ている(特許文献1参照)。For example, in measuring the oxygen saturation, the measured pulse wave signals S1 and S2 include signal components s1 and s2.
And noise components n1 and n2 are included, the signal component s
1 and s1 and s such that the correlation between the noise component n1 is minimized.
There has been proposed a signal processing method configured to obtain a ratio of 2 (see Patent Document 1).
【0005】しかしながら、前記特許文献1に記載され
るような信号処理方法には、次のような問題があった。
すなわち、信号成分s1とノイズ成分n1の相関が最小
になるようなs1とs2の比を求めるために、順次その
比の値を変えながら演算する必要があり、演算量が膨大
となって、測定装置に相当の演算負荷が掛かる他、演算
処理時間が掛かり、即時演算性に問題があった。However, the signal processing method as described in Patent Document 1 has the following problems.
That is, in order to obtain the ratio of s1 and s2 that minimizes the correlation between the signal component s1 and the noise component n1, it is necessary to sequentially change the value of the ratio, and the amount of calculation becomes enormous. In addition to a considerable calculation load on the device, there was a problem with the immediate calculation due to the long calculation processing time.
【0006】[0006]
【特許文献1】米国特許第5632272号明細書(特
許請求の範囲、明細書全文)[Patent Document 1] US Pat. No. 5,632,272 (Claims, full description)
【0007】[0007]
【発明が解決しようとする課題】本発明は、上記問題点
を解決するためになされたものであり、その目的は、演
算量も少なく、即時にノイズ成分を測定信号から除去す
ることにある。特に、体動等によるノイズが重畳された
脈波であっても、ノイズを除去し本来の脈波を再現する
ことにある他、パルスフォトメトリ法による酸素飽和度
を精度良く測定する方法を提供することである。SUMMARY OF THE INVENTION The present invention has been made to solve the above problems, and an object thereof is to immediately remove a noise component from a measurement signal with a small amount of calculation. In particular, even if it is a pulse wave on which noise due to body motion etc. is superimposed, it is to remove the noise and reproduce the original pulse wave, as well as provide a method for accurately measuring oxygen saturation by the pulse photometry method It is to be.
【0008】[0008]
【課題を解決するための手段】本発明に係る同じ基本周
波数を有し連続する第1の信号IRおよび第2の信号R
を処理する信号処理方法は、前記第1の信号IRおよび
前記第2の信号Rのそれぞれについて、所定期間におい
て、周波数スペクトルまたは周波数パワースペクトルを
演算し、前記第1の信号IRの第1スペクトルと、前記
第2の信号Rの第2スペクトルを演算するスペクトル演
算ステップと、前記スペクトル演算ステップにより演算
された前記第1スペクトルおよび前記第2スペクトルを
用い、周波数軸上で、相互の差および和の比の演算、ま
たはそれと等価の演算をするスペクトル正規化演算ステ
ップと、前記スペクトル正規化演算ステップにおいて演
算された結果における、前記基本周波数における演算値
を用いて、前記第1の信号IRの信号成分S1の振幅お
よび前記第2の信号Rの信号成分S2の振幅の比ΦAを
求める信号成分振幅比演算ステップと、を有することを
特徴とする(請求項1)。これをもって、前記第1の信
号IRおよび前記第2の信号Rにノイズが含まれる場合
であっても、その信号成分S1の振幅および信号成分S
2の振幅の比ΦAを求めることを目的とする。SUMMARY OF THE INVENTION A first signal IR and a second signal R which have the same fundamental frequency and are continuous according to the invention.
In the signal processing method for processing the first signal IR and the second signal R, a frequency spectrum or a frequency power spectrum is calculated in a predetermined period for each of the first signal IR and the second signal R, and , Using a spectrum calculation step of calculating a second spectrum of the second signal R, and the first spectrum and the second spectrum calculated in the spectrum calculation step, on the frequency axis, A signal component of the first signal IR is calculated by using a spectrum normalization calculation step for calculating a ratio or an equivalent calculation, and a calculation value at the fundamental frequency in the result calculated in the spectrum normalization calculation step. A signal component shake for obtaining a ratio ΦA of the amplitude of S1 and the amplitude of the signal component S2 of the second signal R. And having a ratio calculating step (claim 1). With this, even when the first signal IR and the second signal R include noise, the amplitude of the signal component S1 and the signal component S
The purpose is to find the ratio ΦA of the two amplitudes.
【0009】また、本発明に係る信号処理方法は、さら
に、前記第1の信号IRおよび前記第2の信号Rは、パ
ルスフォトメトリ法により測定された脈波に関するデー
タであることを特徴とする(請求項2)。これをもっ
て、パルスフォトメトリ法により測定された2つの脈波
の吸光度比ΦAを求めることを目的とする。Further, the signal processing method according to the present invention is further characterized in that the first signal IR and the second signal R are data relating to a pulse wave measured by a pulse photometry method. (Claim 2). With this, the purpose is to obtain the absorbance ratio ΦA of the two pulse waves measured by the pulse photometry method.
【0010】また、本発明に係る信号処理方法は、さら
に、前記第1の信号IRは、パルスフォトメトリ法によ
り測定された赤外光脈波に関するデータであり、前記第
2の信号Rは、パルスフォトメトリ法により測定された
赤色光脈波に関するデータであり、さらに、前記信号成
分振幅比演算ステップにより演算された比ΦAを用い
て、酸素飽和度を測定することを特徴とする(請求項
3)。これをもって、パルスフォトメトリ法により測定
された赤外光および赤色光の脈波の吸光度比ΦAを求
め、精度のよい酸素飽和度を測定することを目的とす
る。Further, in the signal processing method according to the present invention, the first signal IR is data relating to an infrared light pulse wave measured by a pulse photometry method, and the second signal R is The oxygen saturation is measured by using the ratio ΦA calculated in the signal component amplitude ratio calculation step, which is data regarding the red light pulse wave measured by the pulse photometry method. 3). With this, the purpose is to obtain the absorbance ratio ΦA of the pulse wave of infrared light and the pulse wave of red light measured by the pulse photometry method, and to measure the oxygen saturation with high accuracy.
【0011】また、本発明に係る同じ基本周波数を有し
連続する第1の信号IRおよび第2の信号Rを処理する
信号処理方法は、前記第1の信号IRおよび前記第2の
信号Rのそれぞれについて、所定期間において、周波数
スペクトルまたは周波数パワースペクトルを演算し、前
記第1の信号IRの第1スペクトルと、前記第2の信号
Rの第2スペクトルを演算するスペクトル演算ステップ
と、前記スペクトル演算ステップにより演算された前記
第1スペクトルおよび前記第2スペクトルを用い、周波
数軸上で、相互の差および和の比の演算、またはそれと
等価の演算をするスペクトル正規化演算ステップと、前
記スペクトル正規化演算ステップにおいて演算された結
果における、前記基本周波数における演算値を用いて、
前記第1の信号IRの信号成分S1の振幅および前記第
2の信号Rの信号成分S2の振幅の比ΦAを求める信号
成分振幅比演算ステップと、前記スペクトル正規化演算
ステップにおいて演算された結果における、ノイズ基本
周波数における演算値を用いて、前記第1の信号IRの
ノイズ成分N1の振幅および前記第2の信号Rのノイズ
成分N2の振幅の比ΦNを求めるノイズ成分振幅比演算
ステップと、前記信号成分振幅比演算ステップにおいて
求めた前記比ΦAおよび前記ノイズ成分振幅比演算ステ
ップにおいて求めた前記比ΦNを用いて、前記第1の信
号IRまたは前記第2の信号Rからノイズを除去する処
理を行うノイズ除去処理ステップと、を有することを特
徴とする(請求項4)。これをもって、前記第1の信号
IRおよび前記第2の信号Rにノイズが含まれる場合で
あっても、ノイズ成分を除去し、その信号成分S1の振
幅および信号成分S2の振幅の比ΦAを求めることを目
的とする。A signal processing method for processing a continuous first signal IR and a second signal R having the same fundamental frequency according to the present invention is a method for processing the first signal IR and the second signal R. For each of them, a spectrum calculation step of calculating a frequency spectrum or a frequency power spectrum and calculating a first spectrum of the first signal IR and a second spectrum of the second signal R in a predetermined period, and the spectrum calculation. Using the first spectrum and the second spectrum calculated in the step, a spectrum normalization calculation step of performing a calculation of a mutual difference and a sum ratio on the frequency axis, or a calculation equivalent thereto, and the spectrum normalization Using the calculated value at the fundamental frequency in the result calculated in the calculation step,
In the signal component amplitude ratio calculation step for obtaining the ratio ΦA of the amplitude of the signal component S1 of the first signal IR and the amplitude of the signal component S2 of the second signal R, and the result calculated in the spectrum normalization calculation step. A noise component amplitude ratio calculation step for obtaining a ratio ΦN between the amplitude of the noise component N1 of the first signal IR and the amplitude of the noise component N2 of the second signal R using the calculated value at the noise fundamental frequency; A process for removing noise from the first signal IR or the second signal R using the ratio ΦA obtained in the signal component amplitude ratio calculation step and the ratio ΦN obtained in the noise component amplitude ratio calculation step. And a noise removal processing step to be performed (claim 4). With this, even when the first signal IR and the second signal R include noise, the noise component is removed and the ratio ΦA of the amplitude of the signal component S1 and the amplitude of the signal component S2 is obtained. The purpose is to
【0012】また、本発明に係る信号処理方法は、さら
に、前記第1の信号IRおよび前記第2の信号Rは、パ
ルスフォトメトリ法により測定された脈波に関するデー
タであることを特徴とする(請求項5)。これをもっ
て、パルスフォトメトリ法により測定された脈波からノ
イズ成分を除去し、信号成分を求めることを目的とす
る。Further, the signal processing method according to the present invention is further characterized in that the first signal IR and the second signal R are data relating to a pulse wave measured by a pulse photometry method. (Claim 5). The purpose is to remove the noise component from the pulse wave measured by the pulse photometry method and obtain the signal component.
【0013】また、本発明に係る信号処理方法は、さら
に、前記ノイズ除去処理ステップによりノイズが除去さ
れた前記第1の信号IRまたは前記第2の信号Rを、表
示部に表示させることを特徴とする(請求項6)。これ
をもって、パルスフォトメトリ法により測定された脈波
からノイズ成分を除去された脈波を、表示部に表示させ
ることを目的とする。Further, the signal processing method according to the present invention is characterized in that the first signal IR or the second signal R from which noise has been removed by the noise removal processing step is displayed on a display section. (Claim 6). With this, it is an object to display a pulse wave in which a noise component is removed from the pulse wave measured by the pulse photometry method on the display unit.
【0014】[0014]
【発明の実施の形態】以下、本発明に係る信号処理方法
としての脈波信号処理方法およびそれを利用したパルス
オキシメータの実施の形態を、図面を参照しながら詳細
に説明する。図1は、本発明の実施の形態のパルスオキ
シメータである。プローブ1は、発光部2と受光部3を
備え、これらにより指先(生体組織)4を挾持する構成
となっている。発光部2は、第1波長光である赤外光
(波長λ1:940nm)と第2波長光である赤色光
(波長λ2:660nm)とをそれぞれ発光する2つの
発光ダイオードを備えている。発光部2は、発光部駆動
回路5により駆動されるものであり、赤色光と赤外光は
交互に発光される。受光部3は、フォトダイオードを備
え、指先の透過光を受光し、その透過光強度に応じた電
気信号を出力するものである。受光部3の出力信号は、
受光信号増幅回路6で増幅され、復調回路7で復調され
るようになっている。復調回路7は、赤色光と赤外光に
応じたそれぞれの信号を分けて出力する。これらの信号
は、増幅器9a、9bで増幅され、A/D変換器10
a、10bでデジタル信号に変換されて、CPU(セン
トラルプロセッシングユニット)8に入力されるように
なっている。CPU8は、復調回路7と発光部駆動回路
5を制御すると共に、A/D変換器10a、10bから
与えられる信号を処理して、その結果を表示部11に出
力する。本発明では、この信号処理における脈波のノイ
ズ除去処理に特徴がある。表示部11に出力される結果
は、ノイズが除去された脈波波形、脈拍数、SpO2 値
(酸素飽和度S)である。BEST MODE FOR CARRYING OUT THE INVENTION Hereinafter, embodiments of a pulse wave signal processing method as a signal processing method according to the present invention and a pulse oximeter using the same will be described in detail with reference to the drawings. FIG. 1 is a pulse oximeter according to an embodiment of the present invention. The probe 1 includes a light emitting unit 2 and a light receiving unit 3, and is configured to hold a fingertip (living tissue) 4 by these. The light emitting unit 2 includes two light emitting diodes that emit infrared light (wavelength λ1: 940 nm) that is the first wavelength light and red light (wavelength λ2: 660 nm) that is the second wavelength light, respectively. The light emitting section 2 is driven by the light emitting section drive circuit 5, and red light and infrared light are emitted alternately. The light receiving section 3 includes a photodiode, receives the transmitted light from the fingertip, and outputs an electric signal according to the transmitted light intensity. The output signal of the light receiving unit 3 is
The received light signal amplification circuit 6 amplifies the signal and the demodulation circuit 7 demodulates it. The demodulation circuit 7 separately outputs signals corresponding to red light and infrared light. These signals are amplified by the amplifiers 9a and 9b, and the A / D converter 10
It is adapted to be converted into a digital signal by a and 10b and input to a CPU (Central Processing Unit) 8. The CPU 8 controls the demodulation circuit 7 and the light emitting unit drive circuit 5, processes the signals given from the A / D converters 10 a and 10 b, and outputs the result to the display unit 11. The present invention is characterized by pulse wave noise removal processing in this signal processing. The results output to the display unit 11 are the pulse wave waveform from which noise is removed, the pulse rate, and the SpO2 value (oxygen saturation S).
【0015】以下、図2に示したフローチャートに沿っ
て説明する。
(ステップ1)の「スタート」において、測定装置を起
動させる。
(ステップ2)「生体組織を透過または反射した赤外光
(IR)・赤色光(R)の脈波測定」については、図1
で説明したように、赤外光(IR)・赤色光(R)の脈
波が測定され、CPU8に入力される。図3に示すよう
な脈波波形が測定される。ここで、酸素飽和度Sは、前
述したように、吸光度の比Φの関数fとして産出するこ
とができる。これを説明すると、まず、
赤外光の透過光または反射光から得られる脈波の直流成
分に対する脈動成分の比:IR
赤色光の透過光または反射光から得られる脈波の直流成
分に対する脈動成分の比:R
とする。これらの比により、それぞれの波長における吸
光度の脈動成分が近似される。また、それぞれの脈動成
分には、本来の脈動による信号成分SiとノイズNiが
含まれる。つまり、
IR=S1+N1 (3)
R=S2+N2 (4)
ここで、
S1:信号成分、N1:ノイズ成分
S2:信号成分、N2:ノイズ成分
である。また、ここで次のように定義する。
ΦA=S2/S1 (5)
ΦN=N2/N1 (6)
このとき、酸素飽和度Sとの間で、
S=f(ΦA) (7)
の関係が成り立つのである。なお、ノイズ成分が除去さ
れない状態で、式(5)の代わりに、
Φ=R/IR (8)
を用いると、酸素飽和度Sを演算した結果も測定誤差を
含むことになる。The operation will be described below with reference to the flow chart shown in FIG. In "Start" of (Step 1), the measuring device is started. (Step 2) For “measurement of pulse wave of infrared light (IR) / red light (R) transmitted or reflected by living tissue”, see FIG.
As described above, the pulse waves of infrared light (IR) and red light (R) are measured and input to the CPU 8. A pulse wave waveform as shown in FIG. 3 is measured. Here, the oxygen saturation S can be produced as a function f of the absorbance ratio Φ, as described above. To explain this, first, the ratio of the pulsating component to the DC component of the pulse wave obtained from the transmitted or reflected light of infrared light: IR The pulsating component to the DC component of the pulse wave obtained from the transmitted or reflected light of red light Ratio: R These ratios approximate the pulsating component of the absorbance at each wavelength. Further, each pulsation component includes a signal component Si and noise Ni due to the original pulsation. That is, IR = S1 + N1 (3) R = S2 + N2 (4) Here, S1: signal component, N1: noise component, S2: signal component, and N2: noise component. In addition, here, the definition is as follows. ΦA = S2 / S1 (5) ΦN = N2 / N1 (6) At this time, the relationship of S = f (ΦA) (7) holds with the oxygen saturation S. If Φ = R / IR (8) is used instead of the equation (5) in the state where the noise component is not removed, the result of calculating the oxygen saturation S also includes a measurement error.
【0016】ここで、上記各式を用いて式変形を試み
る。まず、測定された信号から信号成分を取り除くた
め、次のような定義をし、式変形をする。
ここで、式(5)よりS2=ΦA・S1であるから、さ
らに式変形をして、
N′=(S2−S2)+(N2−ΦA・N1)
=N2−ΦA・N1
=N1・(N/N2−ΦA)
また、式(6)より、
N′=N1・(ΦN−ΦA)
となる。つまり、
N′=R−ΦA・IR
=N1・(ΦN−ΦA)
であるから、
N1=(R−ΦA・IR)/(ΦN−ΦA) (10)
=N′/(ΦN−ΦA) (11)
N2=N1・ΦN
と式変形できる。Here, an equation transformation is attempted by using the above equations. First, in order to remove the signal component from the measured signal, the following definition is made and the formula is modified. Here, since S2 = ΦA · S1 from the formula (5), the formula is further modified so that N ′ = (S2-S2) + (N2-ΦA · N1) = N2-ΦA · N1 = N1 · ( N / N2-ΦA) Further, from the equation (6), N '= N1 · (ΦN-ΦA). That is, since N ′ = R−ΦA · IR = N1 · (ΦN−ΦA), N1 = (R−ΦA · IR) / (ΦN−ΦA) (10) = N ′ / (ΦN−ΦA) ( 11) The equation can be transformed into N2 = N1.ΦN.
【0017】一方、測定されたIR、Rをフーリエ変換
し、スペクトルの絶対値ないしパワースペクトルを演算
することで、脈波の特徴的な情報が得られることを説明
する。測定された第1波長の赤外光(IR)および第2
波長の赤色光(R)についての脈波が、図3に示すよう
に得られる。そして、所定期間毎(期間[i]、期間
[i+1]、期間[i+2]、・・・)に、それぞれの
データサンプルを用いてフーリエ変換を行い、スペクト
ルの絶対値ないしパワースペクトルを求める。図4
(a)、(b)に、赤外光(IR)のスペクトル(Spc.
IR)、赤色光(R)のスペクトル(Spc.R)を示す。
さらに、周波数に合わせて、下記の演算をすることによ
り、図5に示すようになる。
(Spc.IR−Spc.R)/(Spc.IR+Spc.R) (12)
なお、この等価の演算式として次の式(13)を用いて
もよい。
{1−(Spc.R/Spc.IR)}/{1+(Spc.R/Spc.IR)} (13)
ここで、Spc.R/Spc.IRは、Φに等しいと置けるの
で、
(1−Φ)/(1+Φ) (14)
とも表すことができる。従って、脈波の基本周波数fs
における値は、
(1−ΦA)/(1+ΦA) (15)
である。また、ノイズfnにおける値は、
(1−ΦN)/(1+ΦN) (16)
である。On the other hand, it will be explained that characteristic information of the pulse wave can be obtained by Fourier-transforming the measured IR and R and calculating the absolute value of the spectrum or the power spectrum. Measured first wavelength infrared light (IR) and second
A pulse wave for red light (R) of wavelength is obtained as shown in FIG. Then, for each predetermined period (period [i], period [i + 1], period [i + 2], ...), Fourier transform is performed using each data sample to obtain the absolute value of the spectrum or the power spectrum. Figure 4
In (a) and (b), the spectrum of infrared light (IR) (Spc.
IR) and the spectrum (Spc.R) of red light (R) are shown.
Further, by performing the following calculation according to the frequency, it becomes as shown in FIG. (Spc.IR-Spc.R) / (Spc.IR + Spc.R) (12) The following equation (13) may be used as the equivalent arithmetic expression. {1- (Spc.R / Spc.IR)} / {1+ (Spc.R / Spc.IR)} (13) Here, since Spc.R / Spc.IR can be said to be equal to Φ, (1 It can also be expressed as −Φ) / (1 + Φ) (14). Therefore, the fundamental frequency fs of the pulse wave
The value at is (1-ΦA) / (1 + ΦA) (15). The value of the noise fn is (1-ΦN) / (1 + ΦN) (16).
【0018】図5に示されるように、演算式(12)〜
(14)の結果には、脈波の基本周波数(心拍周期に対
応)fsにピークPs1が現れる。このピークPs1の
値をξsとする。また、fsの高調波の周波数2fs、
3fs、・・・にピークPs2、Ps3、・・・が現れ
る。また、Ps1付近にノイズによるピークPnが現れ
る。ピークPnの値をξnとし、その周波数をfnとす
る。As shown in FIG. 5, the arithmetic expressions (12)-
As a result of (14), a peak Ps1 appears at the fundamental frequency (corresponding to the heartbeat period) fs of the pulse wave. The value of this peak Ps1 is ξs. In addition, the frequency 2fs of the harmonic of fs,
The peaks Ps2, Ps3, ... Appear at 3fs. Further, a peak Pn due to noise appears near Ps1. The value of the peak Pn is ξn, and its frequency is fn.
【0019】そうとすると、式(15)により、
(1−ΦA)/(1+ΦA)=ξs (17)
ΦA=(1−ξs)/(1+ξs) (18)
が成り立ち、また式(16)により、
(1−ΦN)/(1+ΦN)=ξn (19)
ΦN=(1−ξn)/(1+ξn) (20)
が成り立つ。これにより、ΦA、ΦNがそれぞれ演算で
きる。Then, according to equation (15), (1-ΦA) / (1 + ΦA) = ξs (17) ΦA = (1-ξs) / (1 + ξs) (18) holds, and according to formula (16) , (1-ΦN) / (1 + ΦN) = ξn (19) ΦN = (1-ξn) / (1 + ξn) (20) holds. Thereby, ΦA and ΦN can be calculated respectively.
【0020】酸素飽和度Sは、式(18)により求めら
れたΦAを用いて、直接式(7)より求めることができ
る。また、式(3)、(4)、(10)、(6)より、
S1=IR−N1
=IR−(R−ΦA・IR)/(ΦN−ΦA) (21)
S2=R−N2
=R−N1・ΦN (22)
が成り立つ。よって、測定されたIR、Rおよび演算さ
れたΦA、ΦNを、式(21)、(22)に代入して演
算することにより、信号成分S1、S2を求めることが
できる。これを表示部(ディスプレイ)に表示させるこ
とができる。The oxygen saturation S can be directly obtained from the equation (7) by using ΦA obtained from the equation (18). Further, from the formulas (3), (4), (10) and (6), S1 = IR-N1 = IR- (R-ΦA · IR) / (ΦN-ΦA) (21) S2 = R-N2 = R−N1 · ΦN (22) holds. Therefore, the signal components S1 and S2 can be obtained by substituting the measured IR and R and the calculated .PHI.A and .PHI.N into the equations (21) and (22) for calculation. This can be displayed on the display unit (display).
【0021】そこで、(ステップ3)において、所定期
間における脈波データをフーリエ変換し、スペクトル
(スペクトルの絶対値またはパワースペクトル)を演算
する。そして、(ステップ4)では、(ステップ3)に
おいて得られたスペクトルを用いて、周波数毎に、式
(15)、(16)に基づき、基本周波数、ノイズ周波
数における演算値ξs、ξnを演算する。次に、(ステ
ップ5)では、(ステップ4)において得られた演算値
ξs、ξnから、式(18)、(20)を用いて、Φ
A、ΦNを演算する。さらに、(ステップ6)では、
(ステップ5)において得られたΦA、ΦNから、式
(7)、(21)または(22)を用いて、酸素飽和度
S、脈波の信号成分S1またはS2を演算して、表示部
(ディスプレイ)に、酸素飽和度S、脈波の信号成分S
1またはS2を表示する。演算する期間は、フーリエ変
換に用いたデータ期間分(単位期間:期間[i])のデ
ータでよい。そして、(ステップ2)に戻り、次の所定
期間(フーリエ変換する単位期間:期間[i+1])の
データを用いて、(ステップ2)〜(ステップ6)の処
理を繰り返す。「終了」させる場合(ステップ7)は、
測定装置をオフ状態とする。なお、単位期間は、脈波に
急激な変化が現れない期間、例えば6秒程度がよい。ま
た、本発明の適用は、酸素飽和度の測定に限られず、血
中にある吸光物質の濃度測定にも適用可能である。Therefore, in (step 3), the pulse wave data in a predetermined period is Fourier transformed to calculate the spectrum (absolute value of spectrum or power spectrum). Then, in (step 4), the calculated values ξs and ξn at the fundamental frequency and the noise frequency are calculated for each frequency based on the equations (15) and (16) using the spectrum obtained in (step 3). . Next, in (step 5), using equations (18) and (20) from the calculated values ξs and ξn obtained in (step 4), Φ
Calculate A and ΦN. Furthermore, in (step 6),
The oxygen saturation S and the pulse wave signal component S1 or S2 are calculated from ΦA and ΦN obtained in (Step 5) using the formula (7), (21) or (22), and the display unit ( Display), oxygen saturation S, pulse wave signal component S
1 or S2 is displayed. The calculation period may be data for the data period (unit period: period [i]) used in the Fourier transform. Then, returning to (step 2), the processing of (step 2) to (step 6) is repeated using the data of the next predetermined period (unit period for Fourier transform: period [i + 1]). If you want to "end" (step 7),
Turn off the measuring device. It should be noted that the unit period is preferably a period in which a rapid change does not appear in the pulse wave, for example, about 6 seconds. Further, the application of the present invention is not limited to the measurement of oxygen saturation, but can be applied to the measurement of the concentration of a light-absorbing substance in blood.
【0022】[0022]
【発明の効果】以上、詳記したように、請求項1に係る
信号処理方法によれば、前記第1の信号IRおよび前記
第2の信号Rにノイズが含まれる場合であっても、その
信号成分S1の振幅および信号成分S2の振幅の比ΦA
を求めることができる。As described above in detail, according to the signal processing method of the first aspect, even when the first signal IR and the second signal R contain noise, Ratio ΦA of amplitude of signal component S1 and amplitude of signal component S2
Can be asked.
【0023】また、請求項2に係る信号処理方法によれ
ば、パルスフォトメモリ法により測定された2つの脈波
の信号成分についての吸光度比ΦAを求めることができ
る。Further, according to the signal processing method of the second aspect, it is possible to obtain the absorbance ratio ΦA for the two pulse wave signal components measured by the pulse photomemory method.
【0024】また、請求項3に係る信号処理方法によれ
ば、パルスフォトメモリ法により測定された赤外光およ
び赤色光の脈波の信号成分についての吸光度比ΦAを求
め、精度の良い酸素飽和度を測定することができる。Further, according to the signal processing method of the third aspect, the absorbance ratio ΦA for the signal components of the pulse waves of infrared light and red light measured by the pulse photomemory method is obtained, and oxygen saturation with high accuracy is obtained. The degree can be measured.
【0025】また、請求項4に係る信号処理方法によれ
ば、前記第1の信号IRおよび前記第2の信号Rにノイ
ズが含まれる場合であっても、ノイズ成分を除去し、そ
の信号成分S1の振幅および信号成分S2の振幅の比Φ
Aを求めることができる。According to the signal processing method of the fourth aspect, even when the first signal IR and the second signal R contain noise, the noise component is removed and the signal component is removed. Ratio Φ of amplitude of S1 and amplitude of signal component S2
You can ask for A.
【0026】また、請求項5に係る信号処理方法によれ
ば、パルスフォトメモリ法により測定された脈波からノ
イズ成分を除去し、信号成分を求めることができる。According to the signal processing method of the fifth aspect, it is possible to remove the noise component from the pulse wave measured by the pulse photomemory method and obtain the signal component.
【0027】また、請求項6に係る信号処理方法によれ
ば、パルスフォトメモリ法により測定された脈波からノ
イズ成分を除去された脈波を、表示部に表示させること
ができる。According to the signal processing method of the sixth aspect, the pulse wave obtained by removing the noise component from the pulse wave measured by the pulse photomemory method can be displayed on the display unit.
【図1】本実施の形態に係る脈波信号処理を行う酸素飽
和度測定装置の構成を示すブロック回路図である。FIG. 1 is a block circuit diagram showing a configuration of an oxygen saturation measuring apparatus that performs pulse wave signal processing according to the present embodiment.
【図2】本実施の形態に係る脈波信号処理の操作プログ
ラムを示すフローチャート図である。FIG. 2 is a flowchart showing an operation program for pulse wave signal processing according to the present embodiment.
【図3】本実施の形態に係る脈波信号処理の処理波形を
示すタイムチャート図である。FIG. 3 is a time chart diagram showing a processing waveform of pulse wave signal processing according to the present embodiment.
【図4】(a)は赤外光脈波のスペクトルを示す波形図
である。(b)は赤色光脈波のスペクトルを示す波形図
である。FIG. 4A is a waveform diagram showing a spectrum of an infrared light pulse wave. (B) is a waveform diagram showing a spectrum of a red light pulse wave.
【図5】スペクトルの正規化の演算結果を示す波形図で
ある。FIG. 5 is a waveform diagram showing a calculation result of spectrum normalization.
【符号の説明】 1 プローブ 2 発光部 3 受光部 4 生体組織(指) 5 発光部駆動回路 6 受光信号増幅回路 7 復調回路 8 CPU 9a 増幅器 9b 増幅器 10a A/D変換器 10b A/D変換器 11 表示部[Explanation of symbols] 1 probe 2 Light emitting part 3 Light receiving part 4 Living tissue (fingers) 5 Light emission part drive circuit 6 Light receiving signal amplification circuit 7 Demodulation circuit 8 CPU 9a amplifier 9b amplifier 10a A / D converter 10b A / D converter 11 Display
Claims (6)
号IRおよび第2の信号Rを処理する信号処理方法にお
いて、 前記第1の信号IRおよび前記第2の信号Rのそれぞれ
について、所定期間において、周波数スペクトルまたは
周波数パワースペクトルを演算し、前記第1の信号IR
の第1スペクトルと、前記第2の信号Rの第1スペクト
ルを演算するスペクトル演算ステップと、 前記スペクトル演算ステップにより演算された前記第1
スペクトルおよび前記第2スペクトルを用い、周波数軸
上で、相互の差および和の比の演算、またはそれと等価
の演算をするスペクトル正規化演算ステップと、 前記スペクトル正規化演算ステップにおいて演算された
結果における、前記基本周波数における演算値を用い
て、前記第1の信号IRの信号成分S1の振幅および前
記第2の信号Rの信号成分S2の振幅の比ΦAを求める
信号成分振幅比演算ステップと、を有することを特徴と
する信号処理方法。1. A signal processing method for processing a continuous first signal IR and second signal R having the same fundamental frequency, wherein a predetermined value is set for each of the first signal IR and the second signal R. A frequency spectrum or a frequency power spectrum is calculated in the period, and the first signal IR is calculated.
Of the first spectrum of the second signal R, and a spectrum calculation step of calculating the first spectrum of the second signal R, and the first spectrum calculated by the spectrum calculation step.
A spectrum normalization operation step of performing a calculation of a mutual difference and a sum ratio on the frequency axis using the spectrum and the second spectrum, or an operation equivalent thereto; and a result calculated in the spectrum normalization operation step. , A signal component amplitude ratio calculating step for obtaining a ratio ΦA of the amplitude of the signal component S1 of the first signal IR and the amplitude of the signal component S2 of the second signal R using the calculated value at the fundamental frequency. A signal processing method comprising:
フォトメトリ法により測定された脈波に関するデータで
あることを特徴とする信号処理方法。2. The signal processing method according to claim 1, wherein the first signal IR and the second signal R are data relating to a pulse wave measured by a pulse photometry method. Processing method.
定された赤外光脈波に関するデータであり、 前記第2の信号Rは、パルスフォトメトリ法により測定
された赤色光脈波に関するデータであり、 さらに、前記信号成分振幅比演算ステップにより演算さ
れた比ΦAを用いて、酸素飽和度を測定することを特徴
とする信号処理方法。3. The signal processing method according to claim 2, wherein the first signal IR is data relating to an infrared light pulse wave measured by a pulse photometry method, and the second signal R is a pulse. A signal processing method, which is data relating to a red light pulse wave measured by a photometry method, and further uses the ratio ΦA calculated in the signal component amplitude ratio calculation step to measure oxygen saturation.
号IRおよび第2の信号Rを処理する信号処理方法にお
いて、 前記第1の信号IRおよび前記第2の信号Rのそれぞれ
について、所定期間において、周波数スペクトルまたは
周波数パワースペクトルを演算し、前記第1の信号IR
の第1スペクトルと、前記第2の信号Rの第2スペクト
ルを演算するスペクトル演算ステップと、 前記スペクトル演算ステップにより演算された前記第1
スペクトルおよび前記第2スペクトルを用い、周波数軸
上で、相互の差および和の比の演算、またはそれと等価
の演算をするスペクトル正規化演算ステップと、 前記スペクトル正規化演算ステップにおいて演算された
結果における、前記基本周波数における演算値を用い
て、前記第1の信号IRの信号成分S1の振幅および前
記第2の信号Rの信号成分S2の振幅の比ΦAを求める
信号成分振幅比演算ステップと、 前記スペクトル正規化演算ステップにおいて演算された
結果における、ノイズ基本周波数における演算値を用い
て、前記第1の信号IRのノイズ成分N1の振幅および
前記第2の信号Rのノイズ成分N2の振幅の比ΦNを求
めるノイズ成分振幅比演算ステップと、 前記信号成分振幅比演算ステップにおいて求めた前記比
ΦAおよび前記ノイズ成分振幅比演算ステップにおいて
求めた前記比ΦNを用いて、前記第1の信号IRまたは
前記第2の信号Rからノイズを除去する処理を行うノイ
ズ除去処理ステップと、を有することを特徴とする信号
処理方法。4. A signal processing method for processing a continuous first signal IR and second signal R having the same fundamental frequency, wherein a predetermined value is set for each of the first signal IR and the second signal R. A frequency spectrum or a frequency power spectrum is calculated in the period, and the first signal IR is calculated.
A first spectrum of the second signal R and a spectrum calculation step of calculating a second spectrum of the second signal R, and the first spectrum calculated by the spectrum calculation step.
A spectrum normalization operation step of performing a calculation of a mutual difference and a sum ratio on the frequency axis using the spectrum and the second spectrum, or an operation equivalent thereto; and a result calculated in the spectrum normalization operation step. A signal component amplitude ratio calculation step of obtaining a ratio ΦA of the amplitude of the signal component S1 of the first signal IR and the amplitude of the signal component S2 of the second signal R using the calculated value at the fundamental frequency, A ratio ΦN between the amplitude of the noise component N1 of the first signal IR and the amplitude of the noise component N2 of the second signal R is calculated by using the calculated value at the noise fundamental frequency in the result calculated in the spectrum normalization calculation step. And a ratio ΦA obtained in the signal component amplitude ratio calculation step. A noise removal process step of performing a process of removing noise from the first signal IR or the second signal R using the ratio ΦN obtained in the noise component amplitude ratio calculation step. Signal processing method.
フォトメトリ法により測定された脈波に関するデータで
あることを特徴とする信号処理方法。5. The signal processing method according to claim 4, wherein the first signal IR and the second signal R are data regarding a pulse wave measured by a pulse photometry method. Processing method.
方法において、 さらに、前記ノイズ除去処理ステップによりノイズが除
去された前記第1の信号IRまたは前記第2の信号R
を、表示部に表示させることを特徴とする信号処理方
法。6. The signal processing method according to claim 4 or 5, further comprising the first signal IR or the second signal R from which noise has been removed by the noise removal processing step.
Is displayed on the display unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002361800A JP2003235819A (en) | 2001-12-14 | 2002-12-13 | Signal processing method and pulse wave signal processing method |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001381374 | 2001-12-14 | ||
JP2001-381374 | 2001-12-14 | ||
JP2002361800A JP2003235819A (en) | 2001-12-14 | 2002-12-13 | Signal processing method and pulse wave signal processing method |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2003235819A true JP2003235819A (en) | 2003-08-26 |
Family
ID=27790859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002361800A Pending JP2003235819A (en) | 2001-12-14 | 2002-12-13 | Signal processing method and pulse wave signal processing method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2003235819A (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004202218A (en) * | 2002-10-31 | 2004-07-22 | Nippon Koden Corp | Signal processing method/apparatus and pulse photometer using it |
WO2008043864A1 (en) * | 2006-10-10 | 2008-04-17 | Universidad De Cádiz | System for determining and monitoring desaturation indices and instantaneous respiratory rate |
JP2009022484A (en) * | 2007-07-19 | 2009-02-05 | Nippon Koden Corp | Method for measuring pulse rate and oxygen saturation during exercise and measuring apparatus |
JP2009072639A (en) * | 2002-10-31 | 2009-04-09 | Nippon Koden Corp | Pulse photometer |
WO2010073908A1 (en) * | 2008-12-26 | 2010-07-01 | コニカミノルタセンシング株式会社 | Biological information signal processing apparatus, biological information signal processing method and biological information measuring apparatus |
JP2011050745A (en) * | 2009-09-03 | 2011-03-17 | Swatch Group Research & Development Ltd | Method and device for measuring pulsation by light waves of two wavelengths |
JP2012235890A (en) * | 2011-05-11 | 2012-12-06 | Konica Minolta Optics Inc | Signal processing unit, signal processing method, and biological information measuring apparatus |
JP2013059439A (en) * | 2011-09-13 | 2013-04-04 | Nippon Koden Corp | Biological signal measurement device |
JP2014509231A (en) * | 2011-02-09 | 2014-04-17 | マサチューセッツ インスティテュート オブ テクノロジー | Ear-mounted multiple vital signs monitor |
JP2014180417A (en) * | 2013-03-19 | 2014-09-29 | Seiko Epson Corp | Signal processing device, pulse wave measuring device and signal processing method |
JP2014180418A (en) * | 2013-03-19 | 2014-09-29 | Seiko Epson Corp | Signal processing device, pulse wave measuring device and signal processing method |
WO2015045109A1 (en) * | 2013-09-27 | 2015-04-02 | パイオニア株式会社 | Measuring instrument |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05269116A (en) * | 1992-02-07 | 1993-10-19 | Boc Health Care Inc | Improved artery blood monitor device |
JPH06178767A (en) * | 1992-07-06 | 1994-06-28 | Mark Ries Robinson | Highly reliable and non-invasion type method for measuring blood gas |
JPH11128184A (en) * | 1997-10-30 | 1999-05-18 | Nippon Colin Co Ltd | Reflection type photoelectric pulse wave detection device |
JP2000296117A (en) * | 1999-04-15 | 2000-10-24 | Nippon Colin Co Ltd | Fingertip photoelectric pulse wave sensor |
JP2001078990A (en) * | 1999-09-13 | 2001-03-27 | Nippon Koden Corp | Apparatus for measuring light absorber in blood |
-
2002
- 2002-12-13 JP JP2002361800A patent/JP2003235819A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05269116A (en) * | 1992-02-07 | 1993-10-19 | Boc Health Care Inc | Improved artery blood monitor device |
JPH06178767A (en) * | 1992-07-06 | 1994-06-28 | Mark Ries Robinson | Highly reliable and non-invasion type method for measuring blood gas |
JPH11128184A (en) * | 1997-10-30 | 1999-05-18 | Nippon Colin Co Ltd | Reflection type photoelectric pulse wave detection device |
JP2000296117A (en) * | 1999-04-15 | 2000-10-24 | Nippon Colin Co Ltd | Fingertip photoelectric pulse wave sensor |
JP2001078990A (en) * | 1999-09-13 | 2001-03-27 | Nippon Koden Corp | Apparatus for measuring light absorber in blood |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004202218A (en) * | 2002-10-31 | 2004-07-22 | Nippon Koden Corp | Signal processing method/apparatus and pulse photometer using it |
JP2009072639A (en) * | 2002-10-31 | 2009-04-09 | Nippon Koden Corp | Pulse photometer |
WO2008043864A1 (en) * | 2006-10-10 | 2008-04-17 | Universidad De Cádiz | System for determining and monitoring desaturation indices and instantaneous respiratory rate |
ES2302446A1 (en) * | 2006-10-10 | 2008-07-01 | Universidad De Cadiz | System for determining and monitoring desaturation indices and instantaneous respiratory rate |
US20100094108A1 (en) * | 2006-10-10 | 2010-04-15 | Universidad De Cadiz | System for determining and monitoring desaturation indices and instantaneous respiratory rate |
JP2009022484A (en) * | 2007-07-19 | 2009-02-05 | Nippon Koden Corp | Method for measuring pulse rate and oxygen saturation during exercise and measuring apparatus |
JPWO2010073908A1 (en) * | 2008-12-26 | 2012-06-14 | コニカミノルタセンシング株式会社 | Biological information signal processing device, biological information signal processing method, and biological information measuring device |
WO2010073908A1 (en) * | 2008-12-26 | 2010-07-01 | コニカミノルタセンシング株式会社 | Biological information signal processing apparatus, biological information signal processing method and biological information measuring apparatus |
JP2011050745A (en) * | 2009-09-03 | 2011-03-17 | Swatch Group Research & Development Ltd | Method and device for measuring pulsation by light waves of two wavelengths |
JP2014509231A (en) * | 2011-02-09 | 2014-04-17 | マサチューセッツ インスティテュート オブ テクノロジー | Ear-mounted multiple vital signs monitor |
JP2012235890A (en) * | 2011-05-11 | 2012-12-06 | Konica Minolta Optics Inc | Signal processing unit, signal processing method, and biological information measuring apparatus |
JP2013059439A (en) * | 2011-09-13 | 2013-04-04 | Nippon Koden Corp | Biological signal measurement device |
US9125605B2 (en) | 2011-09-13 | 2015-09-08 | Nihon Kohden Corporation | Biological signal measuring apparatus |
JP2014180417A (en) * | 2013-03-19 | 2014-09-29 | Seiko Epson Corp | Signal processing device, pulse wave measuring device and signal processing method |
JP2014180418A (en) * | 2013-03-19 | 2014-09-29 | Seiko Epson Corp | Signal processing device, pulse wave measuring device and signal processing method |
WO2015045109A1 (en) * | 2013-09-27 | 2015-04-02 | パイオニア株式会社 | Measuring instrument |
JPWO2015045109A1 (en) * | 2013-09-27 | 2017-03-02 | パイオニア株式会社 | Measuring instrument |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6780158B2 (en) | Signal processing method and pulse wave signal processing method | |
JP4830693B2 (en) | Oxygen saturation measuring apparatus and measuring method | |
US9538927B2 (en) | Optical measurement device and a method for an optical measurement | |
JP4352315B2 (en) | Signal processing method / apparatus and pulse photometer using the same | |
US7361148B2 (en) | Cuff volumetric pulse wave obtaining apparatus, cuff volumetric pulse wave analyzing apparatus, pressure pulse wave obtaining apparatus, and pressure pulse wave analyzing apparatus | |
JP3116252B2 (en) | Pulse oximeter | |
Nabeel et al. | Arterial blood pressure estimation from local pulse wave velocity using dual-element photoplethysmograph probe | |
US6650918B2 (en) | Cepstral domain pulse oximetry | |
JP2001526073A (en) | Artifact reduction in optical volume fluctuation recording | |
US20130137938A1 (en) | Systems and methods for determining differential pulse transit time from the phase difference of two analog plethysmographs | |
JP5568037B2 (en) | Cardiopulmonary resuscitation monitoring device | |
JP2001078990A (en) | Apparatus for measuring light absorber in blood | |
JP2003235819A (en) | Signal processing method and pulse wave signal processing method | |
JP3924636B2 (en) | Pulse wave signal processing method | |
Vinciguerra et al. | Progresses towards a processing pipeline in photoplethysmogram (PPG) based on SiPMs | |
JP3107630B2 (en) | Pulse oximeter | |
US8005523B2 (en) | Signal processing for pulse oximetry | |
JP2001112728A (en) | Pulsimeter | |
US10772513B2 (en) | Blood pressure ratio calculation device, blood pressure ratio calculation method, blood pressure ratio calculation program, and recording medium recording said program | |
JP6636784B2 (en) | Blood pressure information calculation device, blood pressure information calculation method, blood pressure information calculation program, and recording medium for recording the program | |
TW201440725A (en) | Denoising method and apparatus of pulse wave signal and pulse oximetry | |
JP4470056B2 (en) | Signal processing method used for pulse photometer and pulse photometer using the same | |
JP2019037713A (en) | Pulse oximeter | |
US20210085197A1 (en) | Physiological information measurement method, physiological information measurement program, and physiological information measurement apparatus | |
JP4807598B2 (en) | Pulse photometer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041025 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060417 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060711 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20061212 |