JP2011206286A - Apparatus and method for measuring biological information - Google Patents

Apparatus and method for measuring biological information Download PDF

Info

Publication number
JP2011206286A
JP2011206286A JP2010077211A JP2010077211A JP2011206286A JP 2011206286 A JP2011206286 A JP 2011206286A JP 2010077211 A JP2010077211 A JP 2010077211A JP 2010077211 A JP2010077211 A JP 2010077211A JP 2011206286 A JP2011206286 A JP 2011206286A
Authority
JP
Japan
Prior art keywords
histogram
ratio
predetermined
biological information
measurement data
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.)
Granted
Application number
JP2010077211A
Other languages
Japanese (ja)
Other versions
JP5531715B2 (en
Inventor
Kenji Hamaguri
謙治 蛤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Konica Minolta Opto Inc
Original Assignee
Konica Minolta Opto Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Konica Minolta Opto Inc filed Critical Konica Minolta Opto Inc
Priority to JP2010077211A priority Critical patent/JP5531715B2/en
Publication of JP2011206286A publication Critical patent/JP2011206286A/en
Application granted granted Critical
Publication of JP5531715B2 publication Critical patent/JP5531715B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

PROBLEM TO BE SOLVED: To provide an apparatus and a method for measuring biological information which acquire biological information while reducing noises by simple processing and reducing power consumption.SOLUTION: The apparatus for measuring biological information Sa includes a sensor section 1 for measuring a predetermined physiological phenomenon of a living body which is a measurement object and outputting measurement data, a histogram calculating section 35 for calculating a first histogram on the basis of the measurement data measured by the sensor section 1, and a statistic parameter calculating section 37a for acquiring a predetermined statistic parameter in the first histogram acquired by the histogram calculating section 35. The apparatus for measuring biological information acquires a second histogram on the basis of the first histogram acquired by the histogram calculating section 35 and the predetermined statistic parameter calculated by the statistic parameter calculating section 37a, and acquires predetermined biological information on the basis of the second histogram.

Description

本発明は、生体の生理的現象を測定することによって生体に関する生体情報を測定する生体情報測定装置に関し、特に、例えば生体の体動等によってノイズが生じた場合でも例えば血中酸素飽和度や脈拍数等を好適に測定し得る生体情報測定装置および生体情報測定方法に関する。   The present invention relates to a biological information measuring apparatus that measures biological information related to a living body by measuring physiological phenomena of the living body, and in particular, even when noise occurs due to, for example, body movement of the living body, for example, blood oxygen saturation or pulse rate The present invention relates to a biological information measuring device and a biological information measuring method capable of suitably measuring numbers and the like.

生体組織の酸素濃度を監視する意義が、臨床において極めて大きいことは周知の通りである。酸素は、生命活動維持のために最も重要な物質であり、酸素の供給が絶たれると生体組織細胞は、重大な傷害を受けることから、酸素供給に関するパラメータが重要であると考えられている。このため、酸素供給が不安定になり得る場合、例えば、麻酔中、術後、呼吸不全および循環不全等の患者を治療する場合には、特に酸素が適切に供給されているか否かをモニタすることが重要である。この生体組織への酸素供給は、動脈血によって行われる。そのため、生体組織への酸素供給が適切に行われているか否かを把握するために、例えば脈拍数や血中酸素飽和度等の生体に関する生体情報がモニタされる。   As is well known, the significance of monitoring the oxygen concentration in living tissue is extremely great in clinical practice. Oxygen is the most important substance for maintaining life activity, and biological tissue cells are seriously injured when the supply of oxygen is cut off. Therefore, parameters relating to oxygen supply are considered to be important. For this reason, when oxygen supply can become unstable, for example, when treating patients with anesthesia, post-surgery, respiratory failure and circulatory failure, monitor whether or not oxygen is being supplied appropriately This is very important. This oxygen supply to the living tissue is performed by arterial blood. Therefore, in order to grasp whether or not oxygen supply to the living tissue is appropriately performed, living body information such as a pulse rate and blood oxygen saturation is monitored.

これら脈拍数や血中酸素飽和度等の生体情報を測定する装置として、従来、パルスオキシメータと称される生体情報測定装置が知られている。この生体情報測定装置は、動脈血の脈動により生ずる生体組織の透過または反射光量における変動成分を利用するもので、脈動分の減光度を異なる2波長の光で測定してその比から酸素飽和度を求める装置である。この生体情報測定装置では、例えば生体の体動等によってノイズが生じるため、このノイズを除去または低減する必要があり、ノイズを除去する装置が例えば、特許文献1ないし特許文献3に開示されている。   Conventionally, a biological information measuring device called a pulse oximeter is known as a device for measuring biological information such as the pulse rate and blood oxygen saturation. This biological information measuring device uses a fluctuation component in the amount of transmitted or reflected light of living tissue caused by pulsation of arterial blood, and measures the attenuation of pulsation with two different wavelengths of light, and calculates the oxygen saturation from the ratio. It is a device to seek. In this biological information measuring apparatus, noise is generated due to, for example, body movements of the living body. Therefore, it is necessary to remove or reduce the noise. For example, Patent Documents 1 to 3 disclose apparatuses for removing noise. .

この特許文献1に開示のパルスフォトメータは、異なる2つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光部により得られる前記2つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする2次元直交座標系に展開した、各波長の脈波データを、各脈波データの平均値を中心として所定角度に回転させる回転行列を用いて、前記離散的時系列脈波データに含まれるノイズを除去処理した波形を求める波形取得手段と、前記波形取得手段により得られた波形の周波数解析により脈波の基本周波数または脈拍数を求める波形解析手段と、前記波形解析手段の出力から血中の酸素の濃度を求める血中酸素濃度演算手段とを備えている。   The pulse photometer disclosed in Patent Document 1 is a light emitting means for irradiating a living tissue with light of two different wavelengths, and light of each wavelength generated from the light emitting means and transmitted or reflected by the living tissue as an electrical signal. The light-receiving means for converting, and the pulse wave data of each wavelength developed in a two-dimensional orthogonal coordinate system with the discrete time-series pulse wave data of the two wavelengths obtained by the light-receiving unit as the vertical axis or the horizontal axis, Using a rotation matrix that rotates at a predetermined angle around the average value of each pulse wave data, a waveform acquisition unit that obtains a waveform obtained by removing noise included in the discrete time-series pulse wave data, and the waveform acquisition unit Waveform analysis means for obtaining the fundamental frequency or pulse rate of the pulse wave by frequency analysis of the obtained waveform, and blood oxygen concentration calculation means for obtaining the oxygen concentration in the blood from the output of the waveform analysis means Eteiru.

また、前記特許文献2および特許文献3に開示の装置では、赤色光を生体組織に照射することによって得られた赤色脈波信号Rch、および、赤外光を前記生体組織に照射することによって得られた赤外脈波信号IRchが、周波数ドメインへ変換され、この周波数ドメインでの複数のピーク近傍における赤色脈波信号Rchと赤外脈波信号IRchとの比が求められ、この比から信号成分およびノイズ成分が求められ、このノイズ成分を除去した後の脈波信号におけるピークを与える周波数が求められ、この周波数から脈拍数が求められている。   Further, in the devices disclosed in Patent Document 2 and Patent Document 3, the red pulse wave signal Rch obtained by irradiating the biological tissue with red light and the biological tissue obtained by irradiating the biological tissue with infrared light. The obtained infrared pulse wave signal IRch is converted to the frequency domain, and a ratio between the red pulse wave signal Rch and the infrared pulse wave signal IRch in the vicinity of a plurality of peaks in the frequency domain is obtained. And a noise component, a frequency that gives a peak in the pulse wave signal after removing the noise component is obtained, and a pulse rate is obtained from this frequency.

特許第4352315号公報Japanese Patent No. 4352315 米国特許第6002952号明細書US Pat. No. 6,0029,52 米国特許第6067462号明細書US Pat. No. 6,067,462

ところで、前記特許文献1に開示の装置では、ノイズの除去に最適な回転行列を探索するために繰り返し演算が必要であり、このため、その演算処理量が比較的多くなってしまう。また、前記特許文献2および特許文献3に記載の各装置では、周波数ドメインへ変換する処理でフーリエ変換が用いられるので、その演算処理量が比較的多くなってしまう。この結果、これら特許文献1ないし特許文献3に記載の各装置では、その演算処理量が比較的多いために、消費電力が多くなってしまう。   By the way, in the apparatus disclosed in Patent Document 1, iterative computation is required to search for a rotation matrix that is optimal for noise removal, and thus the amount of computation processing is relatively large. Moreover, in each apparatus described in Patent Document 2 and Patent Document 3, since Fourier transform is used in the process of converting to the frequency domain, the amount of calculation processing becomes relatively large. As a result, in each of the devices described in Patent Document 1 to Patent Document 3, the amount of calculation processing is relatively large, and thus power consumption increases.

本発明は、上述の事情に鑑みて為された発明であり、その目的は、より簡易な処理によって生体情報を求めることができ、消費電力を低減することができる生体情報測定装置および生体情報測定方法を提供することである。   The present invention has been made in view of the above-described circumstances, and its object is to obtain biological information by simpler processing and to reduce biological power information and biological information measurement. Is to provide a method.

本発明者は、種々検討した結果、上記目的は、以下の本発明により達成されることを見出した。すなわち、本発明の一態様にかかる生体情報測定装置は、測定対象の生体における所定の生理的現象を測定して測定データを出力する測定部と、前記測定部によって測定された測定データに基づいて第1ヒストグラムを求めるヒストグラム算出部と、前記ヒストグラム算出部によって求められた第1ヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出部と、前記ヒストグラム算出部によって求められた第1ヒストグラムと前記統計パラメータ算出部で算出された前記所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求める生体情報算出部とを備えることを特徴とする。そして、本発明の一態様にかかる生体情報測定方法は、測定対象の生体における所定の生理的現象を測定して測定データを出力する測定工程と、前記測定工程によって測定された測定データに基づいて第1ヒストグラムを求めるヒストグラム算出工程と、前記ヒストグラム算出工程によって求められた第1ヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出工程と、前記ヒストグラム算出工程によって求められた第1ヒストグラムと前記統計パラメータ算出工程で算出された前記所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求める生体情報算出工程とを備えることを特徴とする。   As a result of various studies, the present inventor has found that the above object is achieved by the present invention described below. That is, the biological information measuring device according to one aspect of the present invention is based on a measurement unit that measures a predetermined physiological phenomenon in a living body to be measured and outputs measurement data, and the measurement data measured by the measurement unit. A histogram calculation unit for obtaining a first histogram, a statistical parameter calculation unit for obtaining a predetermined statistical parameter in the first histogram obtained by the histogram calculation unit, a first histogram obtained by the histogram calculation unit, and the statistical parameter calculation And a biological information calculation unit that obtains a second histogram based on the predetermined statistical parameter calculated by the unit and obtains the predetermined biological information based on the second histogram. The biological information measurement method according to one aspect of the present invention is based on a measurement process of measuring a predetermined physiological phenomenon in a measurement target living body and outputting measurement data, and the measurement data measured by the measurement process A histogram calculation step for obtaining a first histogram, a statistical parameter calculation step for obtaining a predetermined statistical parameter in the first histogram obtained by the histogram calculation step, a first histogram obtained by the histogram calculation step, and the statistical parameter calculation A biometric information calculating step of obtaining a second histogram based on the predetermined statistical parameter calculated in the step, and obtaining the predetermined biometric information based on the second histogram.

このような構成の生体情報測定装置および生体情報測定方法では、測定データの第1ヒストグラムが求められ、その第1ヒストグラムの統計パラメーラが求められ、これら第1ヒストグラムおよび統計パラメータに基づいて第2ヒストグラムが求められ、その第2ヒストグラムに基づいて所定の生体情報が求められる。このように本発明にかかる生体情報測定装置および生体情報測定方法は、より簡易な処理によって生体情報を求めることができ、その結果、消費電力を低減することができる。   In the biological information measuring device and the biological information measuring method configured as described above, the first histogram of the measurement data is obtained, the statistical parameter of the first histogram is obtained, and the second histogram is obtained based on the first histogram and the statistical parameters. And predetermined biological information is obtained based on the second histogram. As described above, the biological information measuring apparatus and the biological information measuring method according to the present invention can obtain biological information by simpler processing, and as a result, can reduce power consumption.

また、他の一態様では、上述の生体情報測定装置において、好ましくは、前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、前記統計パラメータ算出部は、前記所定の統計パラメータとして前記第1ヒストグラムf(p)における期待値evを求め、前記生体情報算出部は、g1(p)=f(2×ev−p)とする際に、前記第3比pが前記期待値evよりも大きい場合、または、前記g1(p)−f(p)が0より大きい場合における{g1(p)−f(p)}/∫{g1(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求める。   In another aspect, in the above-described biological information measurement device, preferably, the measurement unit irradiates the biological tissue of the living body with first light having a predetermined first wavelength and transmits the biological tissue of the biological body. Alternatively, the first sensor unit that receives the reflected first light and outputs first measurement data, and irradiates the living tissue of the living body with the second light having a predetermined second wavelength different from the first wavelength, and the living body A sensor unit including a second sensor unit that receives second light transmitted or reflected from the living tissue and outputs second measurement data, wherein the histogram calculation unit is an alternating current for a DC component in the first measurement data. A histogram f (p) of a third ratio p of a first ratio of components and a second ratio of AC components to DC components in the second measurement data is obtained as the first histogram, and the statistical parameter calculation unit The order of The expected value ev in the first histogram f (p) is obtained as a parameter, and when the biological information calculation unit sets g1 (p) = f (2 × ev−p), the third ratio p is the expected value. {G1 (p) -f (p)} / ∫ {g1 (p) -f (p)} dp when the value is larger than the value ev or when g1 (p) -f (p) is larger than 0 Is obtained as the second histogram, and the pulse rate of the living body is obtained based on the second histogram.

この構成によれば、第1ヒストグラムf(p)によるpの期待値evをおよび第2ヒストグラムg1(p)に基づき前記生体の脈拍数を求める生体情報測定装置が提供される。   According to this configuration, there is provided a biological information measuring device that obtains the expected value ev of p based on the first histogram f (p) and the pulse rate of the living body based on the second histogram g1 (p).

また、他の一態様では、上述の生体情報測定装置において、好ましくは、前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、前記統計パラメータ算出部は、前記所定の統計パラメータとして前記第1ヒストグラムf(p)における期待値evおよび移動平均g2(p)を求め、前記生体情報算出部は、前記第3比pが前記期待値evよりも大きい場合、または、前記g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求める。   In another aspect, in the above-described biological information measurement device, preferably, the measurement unit irradiates the biological tissue of the living body with first light having a predetermined first wavelength and transmits the biological tissue of the biological body. Alternatively, the first sensor unit that receives the reflected first light and outputs first measurement data, and irradiates the living tissue of the living body with the second light having a predetermined second wavelength different from the first wavelength, and the living body A sensor unit including a second sensor unit that receives second light transmitted or reflected from the living tissue and outputs second measurement data, wherein the histogram calculation unit is an alternating current for a DC component in the first measurement data. A histogram f (p) of a third ratio p of a first ratio of components and a second ratio of AC components to DC components in the second measurement data is obtained as the first histogram, and the statistical parameter calculation unit The order of An expected value ev and a moving average g2 (p) in the first histogram f (p) are obtained as parameters, and the biological information calculation unit determines that the third ratio p is greater than the expected value ev, or the g2 {G2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp when (p) −f (p) is greater than 0 is obtained as the second histogram, The pulse rate of the living body is obtained based on the histogram.

この構成によれば、第1ヒストグラムf(p)によるpの期待値evおよび移動平均g2(p)を所定の統計パラメータとして求め、この第1ヒストグラムf(p)の期待値evおよび移動平均g2(p)に基づき前記生体の脈拍数を求める生体情報測定装置が提供される。   According to this configuration, the expected value ev of p and the moving average g2 (p) from the first histogram f (p) are obtained as predetermined statistical parameters, and the expected value ev and moving average g2 of the first histogram f (p) are obtained. A biological information measuring device for obtaining the pulse rate of the living body based on (p) is provided.

また、他の一態様では、上述の生体情報測定装置において、好ましくは、前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、前記統計パラメータ算出部は、前記第1ヒストグラムf(p)における移動平均g2(p)を求め、前記移動平均g2(p)のピークを与える第3比p_pを前記所定の統計パラメータとして求め、前記生体情報算出部は、前記第3比pが前記移動平均g2(p)のピークを与える第3比p_pよりも大きい場合、または、前記g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求める。   In another aspect, in the above-described biological information measurement device, preferably, the measurement unit irradiates the biological tissue of the living body with first light having a predetermined first wavelength and transmits the biological tissue of the biological body. Alternatively, the first sensor unit that receives the reflected first light and outputs first measurement data, and irradiates the living tissue of the living body with the second light having a predetermined second wavelength different from the first wavelength, and the living body A sensor unit including a second sensor unit that receives second light transmitted or reflected from the living tissue and outputs second measurement data, wherein the histogram calculation unit is an alternating current for a DC component in the first measurement data. A histogram f (p) of a third ratio p of a first ratio of components and a second ratio of AC components to a DC component in the second measurement data is obtained as the first histogram, and the statistical parameter calculator is 1 hiss A moving average g2 (p) in the gram f (p) is obtained, a third ratio p_p that gives a peak of the moving average g2 (p) is obtained as the predetermined statistical parameter, and the biological information calculating unit is configured to obtain the third ratio. {g2 (p) −f (p) when p is larger than a third ratio p_p giving the peak of the moving average g2 (p) or when g2 (p) −f (p) is greater than 0 } / ∫ {g2 (p) −f (p)} dp is obtained as the second histogram, and the pulse rate of the living body is obtained based on the second histogram.

この構成によれば、第1ヒストグラムf(p)の移動平均g2(p)およびこの移動平均g2(p)のピークを与える第3比p_pを所定の統計パラメータとして求め、この第1ヒストグラムf(p)の移動平均g2(p)に基づき前記生体の脈拍数を求める生体情報測定装置が提供される。   According to this configuration, the moving average g2 (p) of the first histogram f (p) and the third ratio p_p that gives the peak of the moving average g2 (p) are obtained as predetermined statistical parameters, and the first histogram f (p There is provided a biological information measuring device for obtaining the pulse rate of the living body based on the moving average g2 (p) of p).

また、他の一態様では、これら上述の生体情報測定装置において、好ましくは、前記生体情報算出部は、前記第1ヒストグラムの標準偏差の値に応じて前記生体の脈拍数を求める算出式を変更する。   In another aspect, in the above-described biological information measuring device, preferably, the biological information calculation unit changes a calculation formula for calculating the pulse rate of the living body according to a standard deviation value of the first histogram. To do.

この構成によれば、測定データに含まれるノイズの度合いに応じて前記生体の脈拍数を求める算出式を変更する生体情報測定装置が提供される。   According to this configuration, there is provided a biological information measuring device that changes a calculation formula for obtaining the pulse rate of the living body according to the degree of noise included in the measurement data.

また、他の一態様では、これら上述の生体情報測定装置において、好ましくは、前記生体情報算出部は、前記所定の生体情報としてさらに血中酸素飽和度を求める。   In another aspect, in the above-described biological information measuring device, preferably, the biological information calculation unit further calculates blood oxygen saturation as the predetermined biological information.

この構成によれば、脈拍数に加えてさらに血中酸素飽和度も求められる。   According to this configuration, in addition to the pulse rate, the blood oxygen saturation is also obtained.

本発明にかかる生体情報測定装置および生体情報測定方法は、より簡易な処理によって生体情報を求めることができ、その結果、消費電力を低減することができる。   The biological information measuring device and the biological information measuring method according to the present invention can obtain biological information by simpler processing, and as a result, can reduce power consumption.

第1実施形態における生体情報測定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the biometric information measuring apparatus in 1st Embodiment. 第1実施形態における生体情報測定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the biometric information measuring apparatus in 1st Embodiment. 第1実施形態における生体情報測定装置の動作を説明するための図である。It is a figure for demonstrating operation | movement of the biometric information measuring apparatus in 1st Embodiment. 第2実施形態における生体情報測定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the biometric information measuring apparatus in 2nd Embodiment. 第2実施形態における生体情報測定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the biometric information measuring apparatus in 2nd Embodiment. 第2実施形態における生体情報測定装置の動作を説明するための図である。It is a figure for demonstrating operation | movement of the biometric information measuring apparatus in 2nd Embodiment.

以下、本発明にかかる実施の一形態を図面に基づいて説明する。なお、各図において同一の符号を付した構成は、同一の構成であることを示し、適宜、その説明を省略する。   Hereinafter, an embodiment according to the present invention will be described with reference to the drawings. In addition, the structure which attached | subjected the same code | symbol in each figure shows that it is the same structure, The description is abbreviate | omitted suitably.

本実施形態における生体情報測定装置は、測定対象の生体における所定の生理的現象を測定して測定データを出力する測定部と、前記測定部によって測定された測定データに基づいて第1ヒストグラムを求めるヒストグラム算出部と、前記ヒストグラム算出部によって求められた第1ヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出部と、前記ヒストグラム算出部によって求められた第1ヒストグラムと前記統計パラメータ算出部で算出された前記所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求める生体情報算出部とを備えるものである。あるいは、本実施形態における生体情報測定装置は、測定対象の生体における所定の生理的現象を測定して測定データを出力する測定部と、前記測定部によって測定された測定データに基づいてヒストグラムを求めるヒストグラム算出部と、前記ヒストグラム算出部によって求められたヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出部と、前記統計パラメータ算出部で算出された前記所定の統計パラメータに基づいて前記ヒストグラムのクラスを複数の所定範囲に分け、これら複数の所定範囲の中の1つに対応する前記測定データに基づいて前記所定の生体情報を求める生体情報算出部とを備えるものである。以下、より具体的に本実施形態における生体情報測定装置について説明する。   The biological information measuring apparatus according to the present embodiment obtains a first histogram based on a measurement unit that measures a predetermined physiological phenomenon in a living body to be measured and outputs measurement data, and the measurement data measured by the measurement unit. Calculated by a histogram calculation unit, a statistical parameter calculation unit for obtaining a predetermined statistical parameter in the first histogram obtained by the histogram calculation unit, a first histogram obtained by the histogram calculation unit, and the statistical parameter calculation unit A biometric information calculation unit that obtains a second histogram based on the predetermined statistical parameter and obtains the predetermined biometric information based on the second histogram; Alternatively, the biological information measurement device according to the present embodiment obtains a histogram based on a measurement unit that measures a predetermined physiological phenomenon in a measurement target living body and outputs measurement data, and the measurement data measured by the measurement unit. A histogram calculation unit, a statistical parameter calculation unit for obtaining a predetermined statistical parameter in the histogram obtained by the histogram calculation unit, and a plurality of classes of the histogram based on the predetermined statistical parameter calculated by the statistical parameter calculation unit And a biological information calculation unit for obtaining the predetermined biological information based on the measurement data corresponding to one of the plurality of predetermined ranges. Hereinafter, the biological information measuring apparatus in the present embodiment will be described more specifically.

(第1実施形態)
図1は、第1実施形態における生体情報測定装置の構成を示すブロック図である。第1実施形態における生体情報測定装置Saは、測定対象である生体の生理的現象を測定することによって生体に関する生体情報を測定する装置である。生体情報は、例えば、脈拍数や血中酸素飽和度(血中酸素濃度)等である。このような生体情報は、動脈血の脈動により生ずる生体組織の透過光量または反射光量における変動成分を利用することによって求めることができ、その基本原理は、公知の常套手段であって、例えば、特公昭53−026437号公報等に開示されている。
(First embodiment)
FIG. 1 is a block diagram showing the configuration of the biological information measuring apparatus in the first embodiment. The biological information measuring device Sa in the first embodiment is a device that measures biological information related to a living body by measuring a physiological phenomenon of the living body that is a measurement target. The biological information is, for example, a pulse rate or blood oxygen saturation (blood oxygen concentration). Such biological information can be obtained by using a fluctuation component in the transmitted light amount or reflected light amount of a living tissue caused by pulsation of arterial blood, and the basic principle thereof is a well-known conventional means, for example, 53-026437 and the like.

このような生体情報測定装置Saは、例えば、図1に示すように、生体の所定の生理的現象を測定して測定データを出力するセンサ部1と、センサ部1で測定された測定データに基づいて例えば脈拍数や血中酸素飽和度等の生体情報を求める演算制御部3aと、演算制御部3aで測定した生体情報を外部から認識可能に表示する表示部5とを備えて構成される。   For example, as shown in FIG. 1, such a biological information measuring device Sa measures a predetermined physiological phenomenon of a living body and outputs measurement data, and the measurement data measured by the sensor unit 1. Based on, for example, a calculation control unit 3a that obtains biological information such as a pulse rate and blood oxygen saturation, and a display unit 5 that displays the biological information measured by the calculation control unit 3a so as to be recognized from the outside. .

センサ部1は、演算制御部3aに接続され、本実施形態では、生体の生理的現象として生体組織中の血液に関する情報を測定する装置である。より具体的には、センサ部1は、生体の生理的現象として、心拍による動脈血の脈動により生ずる生体組織の透過光量または反射光量における変動成分を測定する装置である。   The sensor unit 1 is connected to the calculation control unit 3a, and in this embodiment, is a device that measures information related to blood in a living tissue as a physiological phenomenon of the living body. More specifically, the sensor unit 1 is a device that measures a fluctuation component in a transmitted light amount or a reflected light amount of a living tissue caused by arterial blood pulsation due to a heartbeat as a physiological phenomenon of the living body.

このような生体の生理的現象を測定する方法は、例えば、生体組織のヘモグロビンの吸光特性を利用する方法を挙げることができる。酸素は、ヘモグロビンによって生体の各細胞に運ばれるが、ヘモグロビンは、肺で酸素と結合して酸化ヘモグロビンとなり、生体の細胞で酸素が消費されるとヘモグロビン(還元ヘモグロビン)に戻る。血中酸素飽和度は、血中(血液中)の酸化ヘモグロビンの割合として定義される。これらヘモグロビンおよび酸化ヘモグロビンの各吸光度は、波長依存性を有しており、ヘモグロビンは、例えば赤色光(赤色波長領域の光)に対し酸化ヘモグロビンより光を多く吸収し、一方、赤外光(赤外線波長領域の光)に対して酸化ヘモグロビンより光の吸収が少ない。すなわち、ヘモグロビンは、酸化されて酸化ヘモグロビンになると赤色光の吸収が減少して赤外光の吸収が増加し、逆に還元されてヘモグロビンに戻ると赤色光の吸収が増加して赤外光の吸収が減少するという光学的特性を有している。本実施形態の生体情報測定装置Saは、このようなヘモグロビンと酸化ヘモグロビンとの赤色光と赤外光とに対する吸光特性の違いを利用することによって、例えば脈拍数や血中酸素飽和度等の生体情報を求めるものである。すなわち、赤色光の吸光度が大きいほど、還元ヘモグロビンが多く、酸素飽和度が小さくなり、逆に、赤外光の吸光度が大きいほど、酸化ヘモグロビンが多く、酸素飽和度が大きくなる。   As a method for measuring the physiological phenomenon of the living body, for example, a method using the light absorption characteristic of hemoglobin in the living tissue can be mentioned. Oxygen is transported to each cell of the living body by hemoglobin, but hemoglobin combines with oxygen in the lung to become oxygenated hemoglobin, and returns to hemoglobin (reduced hemoglobin) when oxygen is consumed in the cells of the living body. Blood oxygen saturation is defined as the percentage of oxyhemoglobin in the blood (in the blood). The absorbance of each of these hemoglobin and oxyhemoglobin has a wavelength dependency, and hemoglobin absorbs more light than oxyhemoglobin for red light (red wavelength region light), for example, while infrared light (infrared light) Less light is absorbed than oxyhemoglobin. In other words, when hemoglobin is oxidized to oxygenated hemoglobin, the absorption of red light decreases and the absorption of infrared light increases. Conversely, when it is reduced to return to hemoglobin, the absorption of red light increases and the absorption of infrared light increases. It has an optical property that absorption is reduced. The biological information measuring device Sa according to the present embodiment uses such a difference in absorption characteristics of hemoglobin and oxyhemoglobin with respect to red light and infrared light, for example, a living body such as a pulse rate and blood oxygen saturation. It seeks information. That is, the greater the absorbance of red light, the greater the amount of reduced hemoglobin and the lower the oxygen saturation. On the contrary, the greater the absorbance of infrared light, the greater the amount of oxyhemoglobin and the greater the oxygen saturation.

このような方法によるため、本実施形態のセンサ部1は、例えば、赤色光に対する生体組織における吸光特性を測定する第1センサ部11と、赤外光に対する前記生体組織における吸光特性を測定する第2センサ部12とを備えて構成され、演算制御部3aに接続される。この第1センサ部11は、例えば、波長λ1の赤色光を前記生体組織に照射する例えば発光ダイオード等のR発光素子と、前記R発光素子で照射され前記生体組織を透過または反射した赤色光を受光する例えばシリコンホトダイオード等のR受光素子とを備えて構成され、この第2センサ部12は、例えば、前記波長λ1と異なる波長λ2の赤外色光を前記生体組織に照射する例えば発光ダイオード等のIR発光素子と、前記IR発光素子で照射され前記生体組織を透過または反射した赤外光を受光する例えばシリコンホトダイオード等のIR受光素子とを備えて構成される。センサ部1は、このような透過型または反射型のセンサを用いることができる。   Due to such a method, the sensor unit 1 according to the present embodiment includes, for example, a first sensor unit 11 that measures the light absorption characteristics in the biological tissue with respect to red light, and a first sensor that measures the light absorption characteristics in the biological tissue with respect to infrared light. 2 sensor section 12 and is connected to calculation control section 3a. The first sensor unit 11 includes, for example, an R light emitting element such as a light emitting diode that irradiates the biological tissue with red light having a wavelength λ1, and red light that is irradiated by the R light emitting element and transmitted or reflected by the biological tissue. For example, the second sensor unit 12 is configured to irradiate the living tissue with infrared light having a wavelength λ2 different from the wavelength λ1, for example, a light-emitting diode. An IR light emitting element and an IR light receiving element such as a silicon photodiode that receives infrared light irradiated by the IR light emitting element and transmitted or reflected by the living tissue are configured. The sensor unit 1 can use such a transmission type or reflection type sensor.

センサ部1は、例えば、手指や耳朶等や、乳幼児の場合の手の甲、手首、足の甲等の、所定の生体組織にセット可能な構造を有して前記所定の生体組織にセットされ、第1および第2センサ部11、12で測定された測定データを演算制御部3aへ出力する。より具体的には、このような構成の第1センサ部11では、前記R発光素子は、前記生体組織に対し赤色光を照射し、前記R受光素子は、このR発光素子によって前記生体組織に照射された赤色光の前記生体組織を透過または反射した赤色光Rを受光し、この受光した赤色光を光電変換することによって、その受光量に応じた電気信号を前記測定データとして演算制御部3aへ出力する。同様に、第2センサ部12では、前記IR発光素子は、前記生体組織に対し赤外光を照射し、前記IR受光素子は、このIR発光素子によって前記生体組織に照射された赤外光の前記生体組織を透過または反射した赤外光を受光し、この受光した赤外光を光電変換することによって、その受光量に応じた電気信号を前記測定データとして演算制御部3aへ出力する。   The sensor unit 1 has a structure that can be set in a predetermined biological tissue, such as a finger, an earlobe, or the back of the hand, wrist, or foot of an infant, and is set in the predetermined biological tissue. The measurement data measured by the first and second sensor units 11 and 12 are output to the arithmetic control unit 3a. More specifically, in the first sensor unit 11 having such a configuration, the R light emitting element emits red light to the living tissue, and the R light receiving element is applied to the living tissue by the R light emitting element. The red light R that is transmitted or reflected through the living tissue of the irradiated red light is received, and the received red light is photoelectrically converted, whereby an electric signal corresponding to the amount of received light is used as the measurement data as the operation control unit 3a. Output to. Similarly, in the second sensor unit 12, the IR light emitting element irradiates the biological tissue with infrared light, and the IR light receiving element emits infrared light irradiated onto the biological tissue by the IR light emitting element. Infrared light transmitted or reflected through the living tissue is received, and the received infrared light is photoelectrically converted to output an electrical signal corresponding to the amount of received light to the arithmetic control unit 3a as the measurement data.

演算制御部3aは、表示部5に接続され、センサ部1で測定された測定データに基づいて生体情報を求めるとともに、生体情報測定装置Sa全体の制御を司る装置である。演算制御部3aは、例えば、センサ部1で測定された測定データを所定のサンプリング周期(例えば周波数37.5Hz等)でサンプリングすることによって測定データの時系列データをセンサ部1から取得するものである。また例えば、演算制御部3aは、所定の周期でセンサ部1を駆動、すなわち、発光および受光の各動作を行わせることによって、時系列データとして測定データをセンサ部1から取得するものである。また例えば、センサ部1が所定のサンプリング周期でサンプリングすることによって前記生体組織から時系列データとして測定データを測定し、この時系列データの測定データを演算制御部3aへ出力するものである。この測定データは、アナログデータであってもよいが、本実施形態では、ディジタルデータであり、アナログデータからディジタルデータへの変換(AD変換)は、センサ部1または演算制御部3aで行われ、また、必要に応じて、前記AD変換前の測定データを増幅する増幅部をセンサ部1または演算制御部3aにさらに備えてよい。   The arithmetic control unit 3a is connected to the display unit 5 and is a device that obtains biological information based on measurement data measured by the sensor unit 1 and controls the entire biological information measuring device Sa. The arithmetic control unit 3a acquires time-series data of measurement data from the sensor unit 1 by sampling the measurement data measured by the sensor unit 1 at a predetermined sampling period (for example, a frequency of 37.5 Hz, for example). is there. Further, for example, the calculation control unit 3a acquires measurement data from the sensor unit 1 as time-series data by driving the sensor unit 1 with a predetermined cycle, that is, by causing each operation of light emission and light reception to be performed. Further, for example, the sensor unit 1 measures the measurement data as time series data from the living tissue by sampling at a predetermined sampling period, and outputs the measurement data of the time series data to the arithmetic control unit 3a. The measurement data may be analog data, but in the present embodiment, the measurement data is digital data. Conversion from analog data to digital data (AD conversion) is performed by the sensor unit 1 or the arithmetic control unit 3a. Further, if necessary, the sensor unit 1 or the operation control unit 3a may further include an amplification unit that amplifies the measurement data before AD conversion.

より具体的には、この演算制御部3aは、センサ部1によって測定された測定データに基づいて第1ヒストグラムを求め、この求められた第1ヒストグラムにおける所定の統計パラメータを求め、これら求められた第1ヒストグラムと所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求めるものである。言い換えれば、この演算制御部3aは、センサ部1によって測定された測定データに基づいてヒストグラムを求め、この求められたヒストグラムにおける所定の統計パラメータを求め、この求められた所定の統計パラメータに基づいて前記ヒストグラムのクラスを複数の所定範囲に分け、これら複数の所定範囲の中の1つに対応する前記測定データ(前記第2ヒストグラムの測定データ)に基づいて前記所定の生体情報を求めるものである。このような演算制御部3aは、例えば、マイクロプロセッサ、メモリおよびその周辺回路を備えるマイクロコンピュータによって構成される。前記メモリは、センサ部1で測定された測定データに基づいて生体情報を求めるための生体情報演算プログラムや、生体情報測定装置Sa全体を制御するための制御プログラム等の各種のプログラムや、センサ部1で測定された前記測定データや前記プログラムの実行に必要なデータ等の各種のデータを記憶する例えば書き換え可能な不揮発性の記憶素子であるEEPROM(Electrically Erasable Programmable Read Only Memory)や不揮発性の記憶素子であるROM(Read Only Memory)、および、前記マイクロプロセッサのいわゆるワーキングメモリとなる例えば揮発性の記憶素子であるRAM(Random Access Memory)等を備えて構成され、前記マイクロプロセッサは、いわゆるCPU(Central Processing Unit)等であり、前記プログラムを実行することにより、機能的に、例えば、第1前処理部31と、第2前処理部32と、第1バンドパスフィルタ部(第1BPF部)33と、第2バンドパスフィルタ部(第2BPF部)34と、ヒストグラム算出部35と、比算出部36と、統計パラメータ算出部37aと、pv算出部38aと、血中酸素飽和度算出部39と、脈拍数算出部40と、信頼度算出部41とを備える。   More specifically, the arithmetic control unit 3a obtains a first histogram based on measurement data measured by the sensor unit 1, obtains predetermined statistical parameters in the obtained first histogram, and obtains these. A second histogram is obtained based on the first histogram and a predetermined statistical parameter, and the predetermined biological information is obtained based on the second histogram. In other words, the arithmetic control unit 3a obtains a histogram based on the measurement data measured by the sensor unit 1, obtains a predetermined statistical parameter in the obtained histogram, and based on the obtained predetermined statistical parameter. The histogram class is divided into a plurality of predetermined ranges, and the predetermined biological information is obtained based on the measurement data (measurement data of the second histogram) corresponding to one of the plurality of predetermined ranges. . Such an arithmetic control unit 3a is constituted by, for example, a microcomputer including a microprocessor, a memory, and its peripheral circuits. The memory includes various programs such as a biological information calculation program for obtaining biological information based on measurement data measured by the sensor unit 1, a control program for controlling the entire biological information measuring device Sa, and a sensor unit. For example, an EEPROM (Electrically Erasable Programmable Read Only Memory), which is a rewritable nonvolatile memory element, or a nonvolatile memory that stores various data such as the measurement data measured in Step 1 and data necessary for the execution of the program. A ROM (Read Only Memory) which is an element, and a RAM (Random Access Memory) which is a so-called working memory of the microprocessor, for example, which is a volatile storage element, is configured. Central Processing Unit), etc., and functions by executing the program In addition, for example, a first pre-processing unit 31, a second pre-processing unit 32, a first band-pass filter unit (first BPF unit) 33, a second band-pass filter unit (second BPF unit) 34, and a histogram calculation Unit 35, ratio calculation unit 36, statistical parameter calculation unit 37 a, pv calculation unit 38 a, blood oxygen saturation calculation unit 39, pulse rate calculation unit 40, and reliability calculation unit 41.

第1および第2前処理部31、32は、センサ部1から入力された測定データに対し、所定の前処理を行うものである。より具体的には、第1前処理部31は、第1センサ部11から入力された赤色光に係わる測定データに対し、前記R受光素子での暗電流を補正するためのいわゆるダーク処理を行い、そして、その直流成分RDCに対する交流成分RACの第1比(赤色光交直比)R(=RAC/RDC)を算出し、この第1比Rを第1BPF部33へ通知(出力)する。この第1比Rは、ランバート・ビア(Lambert-Beer)の法則による近似から、赤色光に対する前記生体組織の吸光度の変化分に等しいとみなされる。また、第2前処理部32は、第2センサ部12から入力された赤外光に係わる測定データに対し、前記IR受光素子での暗電流を補正するためのいわゆるダーク処理を行い、そして、その直流成分IRDCに対する交流成分IRACの第2比(赤外光交直比)IR(=IRAC/IRDC)を算出し、この第2比IRを第1BPF部33へ通知(出力)する。この第2比IRは、前記ランバート・ビアの法則による近似から、赤外光に対する前記生体組織の吸光度の変化分に等しいとみなされる。前記ダーク処理は、公知の常套手段が用いられ、例えば、前記赤色光に係わる測定データから遮光状態のR受光素子から出力される出力値(暗電流値)Rdarkを減算するとともに、前記赤外光に係わる測定データから遮光状態のIR受光素子から出力される出力値(暗電流値)IRdarkを減算することによって行われる。これら遮光状態のR受光素子およびIR受光素子から出力される各出力値Rdark、IRdarkは、予め測定される。例えば、これら遮光状態のR受光素子およびIR受光素子から出力される各出力値Rdark、IRdarkは、測定の都度、予め測定されてもよく、また例えば、予め測定して記憶しておいてもよい。 The first and second preprocessing units 31 and 32 perform predetermined preprocessing on the measurement data input from the sensor unit 1. More specifically, the first preprocessing unit 31 performs so-called dark processing for correcting dark current in the R light receiving element on the measurement data related to red light input from the first sensor unit 11. Then, the first ratio (red light AC / DC ratio) R (= R AC / R DC ) of the AC component R AC to the DC component R DC is calculated, and the first BPF unit 33 is notified (output). ) The first ratio R is considered to be equal to the change in absorbance of the living tissue with respect to red light, based on the approximation by Lambert-Beer law. The second preprocessing unit 32 performs a so-called dark process for correcting dark current in the IR light receiving element on the measurement data related to the infrared light input from the second sensor unit 12, and A second ratio (infrared light AC / DC ratio) IR (= IR AC / IR DC ) of the AC component IR AC with respect to the DC component IR DC is calculated, and this second ratio IR is notified (output) to the first BPF unit 33. . This second ratio IR is considered to be equal to the change in absorbance of the living tissue with respect to infrared light based on the approximation by the Lambert-Beer law. For the dark process, known conventional means are used. For example, an output value (dark current value) Rdark output from the light-receiving R light receiving element is subtracted from the measurement data related to the red light, and the infrared light is subtracted. This is performed by subtracting the output value (dark current value) IRdark output from the IR light receiving element in the light-shielded state from the measurement data related to. The output values Rdark and IRdark output from the light-receiving R light-receiving element and the IR light-receiving element are measured in advance. For example, the output values Rdark and IRdark output from the R light receiving element and the IR light receiving element in the light shielding state may be measured in advance every measurement, or may be measured and stored in advance, for example. .

第1および第2BPF部33、34は、センサ部1によって測定された測定データから所定のノイズ成分を除去するフィルタであり、動脈血の脈動により生ずる生体組織の透過または反射光量における変動成分として通常含まれる周波数成分以外の周波数成分を除去するものである。第1BPF部33は、例えば、赤色光に対する、動脈血の脈動により生ずる生体組織の透過光量または反射光量における変動成分として通常含まれる周波数成分を含む所定の周波数帯域を通過帯域とするフィルタであり、第1比Rをフィルタ処理(フィルタリング)し、このフィルタ処理した後の第1比Rを比算出部36、血中酸素飽和度算出部39、脈拍数算出部40および信頼度算出部41の各部へ通知する。第2BPF部34は、例えば、赤外光に対する、動脈血の脈動により生ずる生体組織の透過光量または反射光量における変動成分として通常含まれる周波数成分を含む所定の周波数帯域を通過帯域とするフィルタであり、第2比IRをフィルタ処理(フィルタリング)し、このフィルタ処理した後の第2比IRを比算出部36、血中酸素飽和度算出部39、脈拍数算出部40および信頼度算出部41の各部へ通知する。   The first and second BPF units 33 and 34 are filters that remove a predetermined noise component from the measurement data measured by the sensor unit 1, and are usually included as fluctuation components in the amount of transmitted or reflected light of living tissue caused by pulsation of arterial blood. The frequency components other than the frequency components to be removed are removed. The first BPF unit 33 is, for example, a filter that uses a predetermined frequency band including a frequency component that is normally included as a fluctuation component in the transmitted light amount or reflected light amount of biological tissue caused by pulsation of arterial blood with respect to red light. 1 ratio R is filtered (filtered), and the first ratio R after the filtering process is sent to each of the ratio calculating section 36, blood oxygen saturation calculating section 39, pulse rate calculating section 40, and reliability calculating section 41. Notice. The second BPF unit 34 is, for example, a filter having a predetermined frequency band including a frequency component that is normally included as a fluctuation component in a transmitted light amount or a reflected light amount of living tissue caused by pulsation of arterial blood with respect to infrared light as a pass band, The second ratio IR is filtered (filtered), and the second ratio IR after the filter processing is subjected to the ratio calculation unit 36, the blood oxygen saturation calculation unit 39, the pulse rate calculation unit 40, and the reliability calculation unit 41. To notify.

比算出部36は、第2比IRに対する第1比Rの第3比(2波長変動比)p(=R/IR)を算出し、この第3比pをヒストグラム算出部35へ通知するものである。   The ratio calculation unit 36 calculates a third ratio (two-wavelength variation ratio) p (= R / IR) of the first ratio R with respect to the second ratio IR, and notifies the histogram calculation unit 35 of the third ratio p. It is.

ヒストグラム算出部35は、所定の時間範囲内においてセンサ部1によって測定された測定データから求められた第3比pの第1ヒストグラム(度数分布)f(p)を算出し、この第3比pの第1ヒストグラムf(p)を統計パラメータ算出部37aおよびpv算出部38aの各部へ通知するものである。前記所定の時間範囲は、第1ヒストグラムf(p)が統計的に有意となるような個数の測定データを取得することができる時間範囲であり、この測定データの個数は、例えば、300や400や500等である。ヒストグラム算出部35は、第3比pが取り得る数値範囲として想定される所定の範囲を所定の間隔で複数の階級(クラス)に階級分け(クラス分け)を行い、この第1ヒストグラムf(p)を作成する時点から前記所定の時間だけ前までに比算出部36から入力された複数の第3比pのそれぞれをこれらいずれかのクラスに振り分けることによって、各クラスにおける第3比pの頻度を求める。これによってヒストグラム算出部35は、第3比pの第1ヒストグラムf(p)を求める。   The histogram calculator 35 calculates a first histogram (frequency distribution) f (p) of the third ratio p obtained from the measurement data measured by the sensor unit 1 within a predetermined time range, and this third ratio p The first histogram f (p) is notified to each of the statistical parameter calculation unit 37a and the pv calculation unit 38a. The predetermined time range is a time range in which the number of pieces of measurement data in which the first histogram f (p) is statistically significant can be acquired. The number of pieces of measurement data is, for example, 300 or 400 Or 500 mag. The histogram calculation unit 35 classifies a predetermined range assumed as a numerical range that can be taken by the third ratio p into a plurality of classes (classes) at predetermined intervals, and the first histogram f (p ), The frequency of the third ratio p in each class by allocating each of the plurality of third ratios p input from the ratio calculation unit 36 to the one of these classes by the predetermined time before Ask for. Thereby, the histogram calculation unit 35 obtains the first histogram f (p) of the third ratio p.

統計パラメータ算出部37aは、ヒストグラム算出部35で算出された第3比pの第1ヒストグラムf(p)に対し、所定の統計パラメータを算出し、この算出した統計パラメータをpv算出部38aおよび血中酸素飽和度算出部39の各部へ通知するものである。前記所定の統計パラメータは、本実施形態では、第1ヒストグラムf(p)の期待値evおよび標準偏差σである。前記期待値evは、次式A1によって算出される。
期待値ev=∫p×f(p)dp ・・・(A1)
ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。
The statistical parameter calculator 37a calculates a predetermined statistical parameter for the first histogram f (p) of the third ratio p calculated by the histogram calculator 35, and uses the calculated statistical parameter as the pv calculator 38a and the blood This is notified to each part of the medium oxygen saturation calculating unit 39. In the present embodiment, the predetermined statistical parameter is an expected value ev and a standard deviation σ of the first histogram f (p). The expected value ev is calculated by the following equation A1.
Expected value ev = ∫p × f (p) dp (A1)
However, the integral ∫ is executed from the minimum class value a to the maximum class value b in the plurality of classes.

そして、統計パラメータ算出部37aは、この算出した期待値evを用いて、関数g1(p)=f(2×ev−p)を求める。   Then, the statistical parameter calculation unit 37a obtains a function g1 (p) = f (2 × ev−p) using the calculated expected value ev.

pv算出部38aは、統計パラメータ算出部37によって算出された第3比pの第1ヒストグラムf(p)における統計パラメータに基づいて、第3比pを信号成分比paとノイズ成分比pvとに分けた場合におけるノイズ成分比pvを算出し、この算出したノイズ成分比pvを血中酸素飽和度算出部39および脈拍数算出部40の各部へ通知するものである。   The pv calculation unit 38a converts the third ratio p into the signal component ratio pa and the noise component ratio pv based on the statistical parameter in the first histogram f (p) of the third ratio p calculated by the statistical parameter calculation unit 37. In this case, the noise component ratio pv is calculated, and the calculated noise component ratio pv is notified to each of the blood oxygen saturation calculator 39 and the pulse rate calculator 40.

ここで、生体組織の吸光度を測定する場合、その吸光度の変化分の信号成分をsとし、この信号成分に重畳されるノイズ成分をnとする場合に、第1比Rおよび第2比IRは、次の式B1および式B2が成り立つ。
IR=s+n ・・・(B1)
R=s×pa+n×pv ・・・(B2)
Here, when measuring the absorbance of a living tissue, when the signal component corresponding to the change in absorbance is s and the noise component superimposed on this signal component is n, the first ratio R and the second ratio IR are The following formulas B1 and B2 hold.
IR = s + n (B1)
R = s × pa + n × pv (B2)

この信号成分比paは、赤外光における吸光度の変化分の信号成分sと赤色光における吸光度の変化分の信号成分との比であり、この信号成分比paは、一般に、血中酸素飽和度と一対一で対応することが知られている。また、ノイズ成分比pvは、赤外光に対する信号成分sに重畳されるノイズ成分nと赤色光に対する信号成分に重畳されるノイズ成分との比である。   The signal component ratio pa is a ratio between the signal component s corresponding to the change in absorbance in infrared light and the signal component corresponding to the change in absorbance in red light. This signal component ratio pa is generally the oxygen saturation level in blood. It is known to correspond one-to-one. The noise component ratio pv is a ratio between the noise component n superimposed on the signal component s for infrared light and the noise component superimposed on the signal component for red light.

pv算出部38aは、より具体的には、本実施形態では、ヒストグラム算出部35によって求められた第1ヒストグラムf(p)と統計パラメータ算出部37aで算出された所定の統計パラメータとに基づいて第2ヒストグラムh1(p)を求め、この第2ヒストグラムh1(p)を用いてノイズ成分比pvを算出する。本実施形態では、前記関数g1(p)=f(2×ev−p)と定義した場合に、p≦期待値evでは、h1(p)=0と定義され、p>期待値evでは、h1(p)={g1(p)−f(p)}/∫{g1(p)−f(p)}dpと定義され、g1(p)−f(p)≦0では、h1(p)=0と定義され、g1(p)−f(p)>0では、h1(p)={g1(p)−f(p)}/∫{g1(p)−f(p)}dpと定義される。すなわち、
h1(p)=0;p≦期待値evまたはg1(p)−f(p)≦0
h1(p)={g1(p)−f(p)}/∫{g1(p)−f(p)}dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。ただし、p≦evおよびg1(p)−f(p)≦0のpの領域を除く);p>期待値evまたはg1(p)−f(p)>0 ・・・(A2)
そして、ノイズ成分比pvは、次式A3によって算出される。
ノイズ成分比pv=∫p×h1(p)dp ・・・(A3)
ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。
More specifically, in this embodiment, the pv calculation unit 38a is based on the first histogram f (p) obtained by the histogram calculation unit 35 and the predetermined statistical parameter calculated by the statistical parameter calculation unit 37a. A second histogram h1 (p) is obtained, and the noise component ratio pv is calculated using the second histogram h1 (p). In this embodiment, when the function g1 (p) = f (2 × ev−p) is defined, when p ≦ expected value ev, h1 (p) = 0 is defined, and when p> expected value ev, h1 (p) = {g1 (p) −f (p)} / ∫ {g1 (p) −f (p)} dp, and when g1 (p) −f (p) ≦ 0, h1 (p ) = 0 and g1 (p) −f (p)> 0, h1 (p) = {g1 (p) −f (p)} /} {g1 (p) −f (p)} dp Is defined. That is,
h1 (p) = 0; p ≦ expected value ev or g1 (p) −f (p) ≦ 0
h1 (p) = {g1 (p) −f (p)} / ∫ {g1 (p) −f (p)} dp (where the integral ∫ is the maximum from the minimum class value a in the plurality of classes. Up to class value b, except for p ≦ ev and g1 (p) −f (p) ≦ 0, p); p> expected value ev or g1 (p) −f (p )> 0 (A2)
The noise component ratio pv is calculated by the following equation A3.
Noise component ratio pv = ∫p × h1 (p) dp (A3)
However, the integral ∫ is executed from the minimum class value a to the maximum class value b in the plurality of classes.

血中酸素飽和度算出部39は、pv算出部38aによって算出されたノイズ成分比pvに基づいて酸素飽和度、例えば、静脈血の酸素飽和度の推定値SvOや動脈血の酸素飽和度の推定値SpOを求め、この求めた酸素飽和度SvO、SpOを表示部5へ出力するものである。 The blood oxygen saturation calculation unit 39 estimates the oxygen saturation based on the noise component ratio pv calculated by the pv calculation unit 38a, for example, the estimated value SvO 2 of venous blood oxygen saturation or the oxygen saturation of arterial blood. The value SpO 2 is obtained, and the obtained oxygen saturations SvO 2 and SpO 2 are output to the display unit 5.

より具体的には、例えば、血中酸素飽和度算出部39は、予め求められ記憶されている、ノイズ成分比pvと酸素飽和度との関係を表すルックアップテーブル(校正テーブル)を用いることにより、静脈血の酸素飽和度の推定値SvOを求める。また、血中酸素飽和度算出部39は、ノイズ成分比pvに基づいて、例えば、次式B3を用いることにより、信号成分比paを算出し、予め求められ記憶されている、信号成分比paと酸素飽和度との関係を表すルックアップテーブル(校正テーブル)を用いることにより、動脈血の酸素飽和度の推定値SpOを求める。
pa={pv×ΣR(i)×IR(i)−ΣR(i)×R(i)}/{pv×ΣIR(i)×IR(i)−ΣR(i)×IR(i)} ・・・(B3)
More specifically, for example, the blood oxygen saturation calculation unit 39 uses a lookup table (calibration table) representing a relationship between the noise component ratio pv and the oxygen saturation, which is obtained and stored in advance. Then, an estimated value SvO 2 of oxygen saturation of venous blood is obtained. Further, the blood oxygen saturation calculating unit 39 calculates the signal component ratio pa based on the noise component ratio pv, for example, by using the following equation B3, and the signal component ratio pa obtained and stored in advance. The estimated value SpO 2 of the arterial blood oxygen saturation is obtained by using a look-up table (calibration table) representing the relationship between the oxygen saturation and the oxygen saturation.
pa = {pv × ΣR (i) × IR (i) −ΣR (i) × R (i)} / {pv × ΣIR (i) × IR (i) −ΣR (i) × IR (i)} .. (B3)

脈拍数算出部40は、赤色光に係わる前記第1比R、赤外光に係わる前記第2比IRおよびpv算出部38aによって算出されたノイズ成分比pvに基づいて所定の時間内における脈拍数(例えば1分間の脈拍数)を算出し、この脈拍数を表示部5へ出力するものである。   The pulse rate calculation unit 40 calculates the pulse rate within a predetermined time based on the first ratio R related to red light, the second ratio IR related to infrared light, and the noise component ratio pv calculated by the pv calculation unit 38a. (For example, the pulse rate for 1 minute) is calculated, and this pulse rate is output to the display unit 5.

より具体的には、脈拍数算出部40は、赤色光に係わる前記第1比R(i)、赤外光に係わる前記第2比IR(i)およびpv算出部38aによって算出されたノイズ成分比pvを用いて、脈波の信号S(i)=R(i)−pv×IR(i)を求め、この信号S(i)の周期T_periodを求め、この信号S(i)の周期T_periodで60を割り算することによって脈拍数(=60/周期T_period)を求める。   More specifically, the pulse rate calculation unit 40 includes the first ratio R (i) related to red light, the second ratio IR (i) related to infrared light, and the noise component calculated by the pv calculation unit 38a. Using the ratio pv, a pulse wave signal S (i) = R (i) −pv × IR (i) is obtained, a period T_period of the signal S (i) is obtained, and a period T_period of the signal S (i) is obtained. The pulse rate (= 60 / period T_period) is obtained by dividing 60 by.

信頼度算出部41は、所定の信頼度を算出し、この算出した信頼度を表示部5へ出力するものである。   The reliability calculation unit 41 calculates a predetermined reliability and outputs the calculated reliability to the display unit 5.

信頼度は、生体情報に係わる算出値がどの程度信頼することができる値であるかを表す指標(度合い)である。このような信頼度は、例えば、次の式C1ないし式C6のいずれかによって求めることができ、信頼度算出部41は、赤色光に係わる前記第1比Rおよび赤外光に係わる前記第2比IRに基づいて、血中酸素飽和度算出部39によって算出された酸素飽和度に対する所定の信頼度を算出する。このような式によって得られる信頼度zでは、その値zの絶対値が大きいほど血中酸素飽和度の信頼度が低くなる。なお、これら各式において、Σは、iについて和を求める。
z=[Σ{R(i)×IR(i)}/Σ{IR(i)}]−[Σ{R(i)}/ΣR{(i)×IR(i)}] ・・・(C1)
z=[Σ{R(i)×IR(i)}/Σ{IR(i)}−[Σ{R(i)}/Σ{IR(i)}] ・・・(C2)
z=[Σ{R(i)}/Σ{R(i)×IR(i)}]−[Σ{R(i)}/(Σ{IR(i)}] ・・・(C3)
z=[(1/N)×ΣR(i)/IR(i)]−(Σ{R(i)})/(Σ{IR(i)}) ・・・(C4)
z=(ΣR(i)×IR(i))/(Σ{IR(i)})−(1/N)×ΣR(i)/IR(i) ・・・(C5)
z=(Σ{R(i)})/(ΣR(i)×IR(i))−(1/N)×ΣR(i)/IR(i) ・・・(C6)
The reliability is an index (degree) indicating how reliable the calculated value related to biological information is. Such reliability can be obtained by, for example, any one of the following formulas C1 to C6, and the reliability calculation unit 41 includes the first ratio R related to red light and the second ratio related to infrared light. Based on the ratio IR, a predetermined reliability for the oxygen saturation calculated by the blood oxygen saturation calculating unit 39 is calculated. In the reliability z obtained by such an equation, the reliability of the blood oxygen saturation level decreases as the absolute value of the value z increases. In each of these equations, Σ calculates the sum for i.
z = [Σ {R (i) × IR (i)} / Σ {IR (i)} 2 ] − [Σ {R (i)} 2 / ΣR {(i) × IR (i)}]・ (C1)
z = [Σ {R (i) × IR (i)} / Σ {IR (i)} 2 ] 2 − [Σ {R (i)} 2 / Σ {IR (i)} 2 ] ( C2)
z = [Σ {R (i)} 2 / Σ {R (i) × IR (i)}] 2 − [Σ {R (i)} 2 / (Σ {IR (i)} 2 ]... (C3)
z = [(1 / N) × ΣR (i) / IR (i)] 2 − (Σ {R (i)} 2 ) / (Σ {IR (i)} 2 ) (C4)
z = (ΣR (i) × IR (i)) / (Σ {IR (i)} 2 ) − (1 / N) × ΣR (i) / IR (i) (C5)
z = (Σ {R (i)} 2 ) / (ΣR (i) × IR (i)) − (1 / N) × ΣR (i) / IR (i) (C6)

表示部5は、この生体情報測定装置Saの動作状態や演算制御部3aによって求められた生体情報等を表示する装置であり、例えば、液晶表示装置(LCD)や有機EL表示装置等である。表示部5は、例えば、本実施形態では、脈拍数算出部40で算出された脈拍数を表示する脈拍数表示部51と、血中酸素飽和度算出部39で算出された酸素飽和度を表示する血中酸素飽和度表示部52と、信頼度算出部41で算出された信頼度を表示する信頼度表示部53とを備えている。   The display unit 5 is a device that displays the operating state of the biological information measuring device Sa, biological information obtained by the arithmetic control unit 3a, and the like, and is, for example, a liquid crystal display device (LCD) or an organic EL display device. In the present embodiment, for example, the display unit 5 displays the pulse rate display unit 51 that displays the pulse rate calculated by the pulse rate calculation unit 40 and the oxygen saturation calculated by the blood oxygen saturation calculation unit 39. The blood oxygen saturation level display unit 52 and the reliability level display unit 53 that displays the reliability level calculated by the reliability level calculation unit 41 are provided.

次に、第1実施形態の生体情報測定装置Saの動作について説明する。図2は、第1実施形態における生体情報測定装置の動作を示すフローチャートである。図3は、第1実施形態における生体情報測定装置の動作を説明するための図である。図3の比較的太い実線は、第3比pのヒストグラムf(p)を示し、破線は、関数g1(p)(=f(2×ev−p))を示し、比較的細い実線は、第2ヒストグラムh1(p)を表す。   Next, the operation of the biological information measuring device Sa of the first embodiment will be described. FIG. 2 is a flowchart showing the operation of the biological information measuring apparatus in the first embodiment. FIG. 3 is a diagram for explaining the operation of the biological information measuring apparatus according to the first embodiment. The relatively thick solid line in FIG. 3 indicates the histogram f (p) of the third ratio p, the broken line indicates the function g1 (p) (= f (2 × ev−p)), and the relatively thin solid line indicates Represents the second histogram h1 (p).

本実施形態の生体情報測定装置Saでは、例えば、図略の電源スイッチの投入や電源スイッチの投入後に図略の測定開始スイッチの投入によって、測定対象である生体の生体情報の測定が開始され、測定データが取得され、第1ヒストグラムf(p)が作成される。   In the biological information measuring apparatus Sa of the present embodiment, for example, measurement of biological information of a biological object to be measured is started by turning on a power switch (not shown) or turning on a measurement start switch after turning on the power switch, Measurement data is acquired and a first histogram f (p) is created.

より具体的には、まず、センサ部1の第1センサ部11によって、赤色光に係わる測定データRs(暗電流を含む)およびその暗電流Rdarkが測定され、アナログ信号からディジタル信号へ変換されるとともに、センサ部1の第2センサ部12によって、赤外光に係わる測定データIRs(暗電流を含む)およびその暗電流IRdarkが測定され、アナログ信号からディジタル信号へ変換される。続いて、演算制御部3aの第1前処理部31によって、センサ部1から入力された赤色光に係わる測定データRsに対し、ダーク処理(Rs−Rdark)が実行され、第1比Rが算出されるとともに、演算制御部3aの第2前処理部32によって、センサ部1から入力された赤外光に係わる測定データIRsに対し、ダーク処理(IRs−IRdark)が実行され、第2比IRが算出される。続いて、演算制御部3aの第1BPF部33によって、第1前処理部31から通知された第1比Rがフィルタリングされるとともに、演算制御部3aの第2BPF部34によって、第2前処理部32から通知された第2比IRがフィルタリングされる。続いて、演算制御部3aの比算出部36によって、第2比IRに対する第1比Rの第3比pが算出され、演算制御部3aのヒストグラム算出部35によって、所定の時間範囲での第3比pのヒストグラムf(p)が算出される。例えば、現在時刻tiから過去の時刻ti−(N+1)までの時間範囲の測定データについて、図3に実線で示すような第3比pのヒストグラムf(p)が求められる。   More specifically, first, the measurement data Rs (including the dark current) relating to the red light and the dark current Rdark are measured by the first sensor unit 11 of the sensor unit 1 and converted from an analog signal to a digital signal. At the same time, the second sensor unit 12 of the sensor unit 1 measures the measurement data IRs (including the dark current) related to the infrared light and the dark current IRdark, and converts them from an analog signal to a digital signal. Subsequently, dark processing (Rs-Rdark) is performed on the measurement data Rs related to red light input from the sensor unit 1 by the first preprocessing unit 31 of the arithmetic control unit 3a, and the first ratio R is calculated. In addition, dark processing (IRs-IRdark) is performed on the measurement data IRs related to the infrared light input from the sensor unit 1 by the second preprocessing unit 32 of the arithmetic control unit 3a, and the second ratio IR Is calculated. Subsequently, the first ratio R notified from the first preprocessing unit 31 is filtered by the first BPF unit 33 of the calculation control unit 3a, and the second preprocessing unit is filtered by the second BPF unit 34 of the calculation control unit 3a. The second ratio IR notified from 32 is filtered. Subsequently, a third ratio p of the first ratio R with respect to the second ratio IR is calculated by the ratio calculation unit 36 of the calculation control unit 3a, and the histogram calculation unit 35 of the calculation control unit 3a calculates the third ratio p in a predetermined time range. A histogram f (p) of 3 ratio p is calculated. For example, the histogram f (p) of the third ratio p as shown by the solid line in FIG. 3 is obtained for the measurement data in the time range from the current time ti to the past time ti− (N + 1).

このように第3比pの第1ヒストグラムf(p)が作成されると、図2において、まず、ステップS11では、演算制御部3aの統計パラメータ算出部37aは、所定の範囲、例えば、第1ヒストグラムf(p)のクラスにおける最小のクラスの値aから最大のクラスの値bまでの範囲で、第3比pのヒストグラムf(p)に対する期待値evおよび標準偏差σをそれぞれ算出する。   When the first histogram f (p) of the third ratio p is created in this way, in FIG. 2, first, in step S11, the statistical parameter calculation unit 37a of the calculation control unit 3a has a predetermined range, for example, the first The expected value ev and standard deviation σ for the histogram f (p) of the third ratio p are calculated in the range from the minimum class value a to the maximum class value b in the class of one histogram f (p).

続いて、ステップS12では、演算制御部3aは、この求めた標準偏差σが予め設定された所定値σth以下であるか否かを判断する(σ≦σthか?)。このステップS12は、第3比pの第1ヒストグラムf(p)に含まれるノイズの程度を判断するための処理であり、前記所定値σthは、閾値として、後述のノイズ除去処理を行うことなく、測定データから例えば仕様に応じた所定の精度で脈拍数を求めることができるか否かの観点から適宜に設定される。   Subsequently, in step S12, the arithmetic control unit 3a determines whether or not the obtained standard deviation σ is equal to or smaller than a predetermined value σth (whether σ ≦ σth?). This step S12 is a process for determining the degree of noise included in the first histogram f (p) of the third ratio p, and the predetermined value σth is set as a threshold without performing a noise removal process described later. From the viewpoint of whether or not the pulse rate can be obtained from the measurement data with a predetermined accuracy according to the specification, for example, it is set appropriately.

このステップS12の判断の結果、標準偏差σが所定値σth以下である場合(Yes)には、ステップS22の処理が実行され、一方、このステップS12の判断の結果、標準偏差σが所定値σth以下ではない場合(No)には、ステップS13の処理が実行される。   If the standard deviation σ is equal to or smaller than the predetermined value σth as a result of the determination in step S12 (Yes), the process of step S22 is executed, while the standard deviation σ is determined to be the predetermined value σth as a result of the determination in step S12. If not (No), the process of step S13 is executed.

ステップS13では、統計パラメータ算出部37aは、g1(p)(=f(2×ev−p))を求める。例えば、上述の例では、図3に実線で示す第3比pの第1ヒストグラムf(p)に対し、図3に破線で示すg1(p)が求められる。   In step S13, the statistical parameter calculation unit 37a calculates g1 (p) (= f (2 × ev−p)). For example, in the above example, g1 (p) indicated by a broken line in FIG. 3 is obtained with respect to the first histogram f (p) having the third ratio p indicated by a solid line in FIG.

続いて、ステップS14では、演算制御部3aのpv算出部38aは、第3比pが期待値ev以上であるか否かを判断する(p≧evか?)。このステップS14の判断の結果、第3比pが期待値ev以上ではない場合(No、p<ev)には、pv算出部38aは、ステップS16を実行し、当該第3比pに対するh1(p)=0とし、ステップS19の処理を実行する。一方、このステップS14の判断の結果、第3比pが期待値ev以上である場合(Yes、p≧ev)には、pv算出部38aは、ステップS15を実行し、g1(p)−f(p)が0以上であるか否かを判断する(g1(p)−f(p)≧0か?)。   Subsequently, in step S14, the pv calculation unit 38a of the calculation control unit 3a determines whether or not the third ratio p is equal to or higher than the expected value ev (p ≧ ev?). As a result of the determination in step S14, when the third ratio p is not equal to or greater than the expected value ev (No, p <ev), the pv calculation unit 38a executes step S16, and h1 ( p) = 0, and the process of step S19 is executed. On the other hand, as a result of the determination in step S14, when the third ratio p is equal to or greater than the expected value ev (Yes, p ≧ ev), the pv calculation unit 38a executes step S15, and g1 (p) −f It is determined whether (p) is 0 or more (whether g1 (p) −f (p) ≧ 0).

このステップS15の判断の結果、g1(p)−f(p)が0以上ではない場合(No、g1(p)−f(p)<0)には、pv算出部38aは、ステップS17を実行し、当該第3比pに対するh1(p)=0とし、ステップS19の処理を実行する。一方、このステップS15の判断の結果、g1(p)−f(p)が0以上である場合(Yes、g1(p)−f(p)≧0)には、pv算出部38aは、ステップS18を実行し、当該第3比pに対するh1(p)={g1(p)−f(p)}/∫{g1(p)−f(p)}dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。ただし、p≦evおよびg1(p)−f(p)≦0のpの領域を除く)を求める。   If g1 (p) -f (p) is not equal to or greater than 0 as a result of the determination in step S15 (No, g1 (p) -f (p) <0), the pv calculation unit 38a performs step S17. And h1 (p) = 0 for the third ratio p, and the process of step S19 is executed. On the other hand, if g1 (p) −f (p) is 0 or more as a result of the determination in step S15 (Yes, g1 (p) −f (p) ≧ 0), the pv calculation unit 38a S18 is executed, and h1 (p) = {g1 (p) −f (p)} / ∫ {g1 (p) −f (p)} dp for the third ratio p (provided that the integral ∫ Is performed from the minimum class value a to the maximum class value b in the class (except for the p region where p ≦ ev and g1 (p) −f (p) ≦ 0).

このようなステップS14からステップS18までの各処理が、a≦p≦bの範囲における第3比pのそれぞれについて実行され、第2ヒストグラムh1(p)が求められる。例えば、上述の例では、図3に比較的太い実線で示す第3比pの第1ヒストグラムf(p)に対し、図3に比較的細い実線で示す第2ヒストグラムh1(p)が求められる。   Each process from step S14 to step S18 is executed for each of the third ratios p in the range of a ≦ p ≦ b, and the second histogram h1 (p) is obtained. For example, in the above-described example, the second histogram h1 (p) indicated by the relatively thin solid line in FIG. 3 is obtained with respect to the first histogram f (p) having the third ratio p indicated by the relatively thick solid line in FIG. .

そして、ステップS19では、pv算出部38aは、ノイズ成分比pv=∫p×h1(p)dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。)を求める。   In step S19, the pv calculation unit 38a determines that the noise component ratio pv = ∫p × h1 (p) dp (where the integral ∫ is the minimum class value a to the maximum class value b in the plurality of classes). Will be executed by.)

続いて、ステップS20では、脈拍数算出部40は、上述の各処理によって求めたノイズ成分比pvを用いて、脈波の信号S(i)=R(i)−pv×IR(i)を求め、この信号S(i)の周期T_periodを求める。   Subsequently, in step S20, the pulse rate calculator 40 uses the noise component ratio pv obtained by the above-described processes to obtain a pulse wave signal S (i) = R (i) −pv × IR (i). The period T_period of the signal S (i) is obtained.

続いて、ステップS21では、脈拍数算出部40は、このように求めた信号S(i)の周期T_periodで60を割り算することによって脈拍数(=60/周期T_period)を求める。   Subsequently, in step S21, the pulse rate calculation unit 40 calculates the pulse rate (= 60 / period T_period) by dividing 60 by the period T_period of the signal S (i) thus determined.

一方、ステップS22では、脈拍数算出部40は、第2比IR(i)の周期T_periodを求め、続いて、ステップS21を実行し、このように求めた第2比IR(i)の周期T_periodで60を割り算することによって脈拍数(=60/周期T_period)を求める。   On the other hand, in step S22, the pulse rate calculation unit 40 obtains the cycle T_period of the second ratio IR (i), and subsequently executes step S21, and thus obtains the cycle T_period of the second ratio IR (i) thus obtained. The pulse rate (= 60 / period T_period) is obtained by dividing 60 by.

なお、σが所定値σth以下の場合、pv=[Σ{R(i)×IR(i)}/Σ{IR(i)}]+δ としてステップS20にジャンプしてもよい。ここで、δは、正の値であって、一定値でもよいし、[Σ{R(i)×IR(i)}/Σ{IR(i)}]に応じた値としてもよい。 If σ is equal to or smaller than the predetermined value σth, the process may jump to step S20 as pv = [Σ {R (i) × IR (i)} / Σ {IR (i)} 2 ] + δ. Here, δ is a positive value, may be a constant value, or may be a value according to [Σ {R (i) × IR (i)} / Σ {IR (i)} 2 ].

そして、演算制御部3aは、このような各処理によって求めた脈拍数を表示部5の脈拍数表示部51に表示する。   Then, the arithmetic control unit 3 a displays the pulse rate obtained by each of these processes on the pulse rate display unit 51 of the display unit 5.

なお、上述の処理において、期待値evを求める際の積分範囲とノイズ成分比pvを求める際の積分範囲とは、同一であっても異なっていてもよい。   In the above-described processing, the integration range for obtaining the expected value ev and the integration range for obtaining the noise component ratio pv may be the same or different.

また、演算制御部3aの血中酸素飽和度算出部39は、ノイズ成分比pvに基づいて静脈血の酸素飽和度SvOを求めるとともに、ノイズ成分比pvに基づいて信号成分比paを求め、この求めた信号成分比paに基づいて動脈血の酸素飽和度SpOを求める。そして、演算制御部3aは、これら求めた酸素飽和度SvO、SpOを表示部5の血中酸素飽和度表示部52に表示する。 Further, the blood oxygen saturation calculation unit 39 of the calculation control unit 3a calculates the venous blood oxygen saturation SvO 2 based on the noise component ratio pv, and calculates the signal component ratio pa based on the noise component ratio pv. Based on the obtained signal component ratio pa, the oxygen saturation SpO 2 of arterial blood is obtained. Then, the arithmetic control unit 3 a displays the obtained oxygen saturation SvO 2 and SpO 2 on the blood oxygen saturation display unit 52 of the display unit 5.

また、演算制御部3aの信頼度算出部41は、所定の信頼度を求め、演算制御部3aは、この求めた信頼度を表示部5の信頼度表示部53に表示する。   Further, the reliability calculation unit 41 of the calculation control unit 3 a obtains a predetermined reliability, and the calculation control unit 3 a displays the obtained reliability on the reliability display unit 53 of the display unit 5.

このように本実施形態の生体情報測定装置Saでは、測定データR、IRによる第3比pの第1ヒストグラムf(p)が求められ、その第1ヒストグラムf(p)の統計パラメーラとして期待値evが求められ、その期待値evに基づく所定範囲における測定データに基づいて所定の生体情報が求められる。このように本実施形態の生体情報測定装置Saは、より簡易な処理によって生体情報を求めることができ、その結果、消費電力を低減することができる。   As described above, in the biological information measuring device Sa of the present embodiment, the first histogram f (p) of the third ratio p based on the measurement data R and IR is obtained, and an expected value is used as a statistical parameter of the first histogram f (p). ev is obtained, and predetermined biological information is obtained based on measurement data in a predetermined range based on the expected value ev. As described above, the biological information measuring device Sa of the present embodiment can obtain the biological information by simpler processing, and as a result, can reduce power consumption.

そして、本実施形態では、第3比pの第1ヒストグラムf(p)における期待値evを所定の統計パラメータとして求め、この第1ヒストグラムf(p)によるpの期待値evに基づき生体の脈拍数を求め、さらに加えて血中酸素飽和度を求める生体情報測定装置Saが提供される。   In this embodiment, the expected value ev in the first histogram f (p) of the third ratio p is obtained as a predetermined statistical parameter, and the pulse of the living body is based on the expected value ev of p based on the first histogram f (p). A biological information measuring device Sa for obtaining the number and further obtaining the blood oxygen saturation is provided.

次に、別の実施形態について説明する。   Next, another embodiment will be described.

(第2実施形態)
第1実施形態における生体情報測定装置Saは、第3比pの第1ヒストグラムf(p)における期待値evおよび関数g1(p)=f(2×ev−p)を用いて脈拍数および血中酸素飽和度を求めたが、第2実施形態における生体情報測定装置Sbは、第3比pの第1ヒストグラムf(p)における期待値evおよび第3比pに関するその移動平均g2(p)を用いて脈拍数および血中酸素飽和度を求めるものである。
(Second Embodiment)
The biological information measuring apparatus Sa according to the first embodiment uses the expected value ev and the function g1 (p) = f (2 × ev−p) in the first histogram f (p) with the third ratio p, and the pulse rate and blood The biological oxygen measurement device Sb in the second embodiment calculates the intermediate oxygen saturation, but the expected value ev in the first histogram f (p) of the third ratio p and the moving average g2 (p) of the third ratio p. Is used to determine the pulse rate and blood oxygen saturation.

図4は、第2実施形態における生体情報測定装置の構成を示すブロック図である。図4において、このような第2実施形態の生体情報想定装置Sbは、センサ部1と、演算制御部3bと、表示部5とを備えて構成される。第2実施形態の生体情報想定装置Sbにおけるセンサ部1および表示部5は、それぞれ、第1実施形態の生体情報想定装置Saにおけるセンサ部1および表示部5と同様であるので、その説明を省略する。   FIG. 4 is a block diagram showing the configuration of the biological information measuring device in the second embodiment. In FIG. 4, the biological information assumption device Sb according to the second embodiment includes a sensor unit 1, an arithmetic control unit 3 b, and a display unit 5. Since the sensor unit 1 and the display unit 5 in the biological information assumption device Sb of the second embodiment are the same as the sensor unit 1 and the display unit 5 in the biological information assumption device Sa of the first embodiment, respectively, description thereof is omitted. To do.

そして、第2実施形態の生体情報想定装置Sbにおける演算制御部3bは、表示部5に接続され、センサ部1で測定された測定データに基づいて生体情報を求めるとともに、生体情報測定装置Sa全体の制御を司る装置であり、機能的に、例えば、第1前処理部31と、第2前処理部32と、第1バンドパスフィルタ部(第1BPF部)33と、第2バンドパスフィルタ部(第2BPF部)34と、ヒストグラム算出部35と、比算出部36と、統計パラメータ算出部37bと、pv算出部38bと、血中酸素飽和度算出部39と、脈拍数算出部40と、信頼度算出部41とを備える。第2実施形態の生体情報想定装置Sbにおける第1前処理部31、第2前処理部32、第1BPF部33、第2BPF部34、ヒストグラム算出部35、比算出部36、血中酸素飽和度算出部39、脈拍数算出部40および信頼度算出部41は、それぞれ、第1実施形態の生体情報想定装置Saにおける第1前処理部31、第2前処理部32、第1BPF部33、第2BPF部34、ヒストグラム算出部35、比算出部36、血中酸素飽和度算出部39、脈拍数算出部40および信頼度算出部41と同様であるので、その説明を省略する。   And the calculation control part 3b in the biometric information assumption apparatus Sb of 2nd Embodiment is connected to the display part 5, and while calculating | requiring biometric information based on the measurement data measured by the sensor part 1, the biometric information measurement apparatus Sa whole Functionally, for example, a first pre-processing unit 31, a second pre-processing unit 32, a first band-pass filter unit (first BPF unit) 33, and a second band-pass filter unit (Second BPF unit) 34, histogram calculation unit 35, ratio calculation unit 36, statistical parameter calculation unit 37b, pv calculation unit 38b, blood oxygen saturation calculation unit 39, pulse rate calculation unit 40, A reliability calculation unit 41; The first preprocessing unit 31, the second preprocessing unit 32, the first BPF unit 33, the second BPF unit 34, the histogram calculation unit 35, the ratio calculation unit 36, and the blood oxygen saturation in the biological information assumption device Sb of the second embodiment. The calculation unit 39, the pulse rate calculation unit 40, and the reliability calculation unit 41 are respectively a first preprocessing unit 31, a second preprocessing unit 32, a first BPF unit 33, and a first BPF unit 33 in the biological information assumption device Sa of the first embodiment. Since it is similar to the 2BPF unit 34, the histogram calculation unit 35, the ratio calculation unit 36, the blood oxygen saturation calculation unit 39, the pulse rate calculation unit 40, and the reliability calculation unit 41, description thereof is omitted.

統計パラメータ算出部37bは、ヒストグラム算出部35で算出された第3比pの第1ヒストグラムf(p)に対し、所定の統計パラメータを算出し、この算出した統計パラメータをpv算出部38bおよび血中酸素飽和度算出部39の各部へ通知するものである。前記所定の統計パラメータは、本実施形態では、第1ヒストグラムf(p)における期待値ev、標準偏差σおよび第3比pに関する第1ヒストグラムf(p)の移動平均g2(p)である。   The statistical parameter calculator 37b calculates a predetermined statistical parameter for the first histogram f (p) of the third ratio p calculated by the histogram calculator 35, and uses the calculated statistical parameter as the pv calculator 38b and the blood This is notified to each part of the medium oxygen saturation calculating unit 39. In the present embodiment, the predetermined statistical parameter is the moving average g2 (p) of the first histogram f (p) related to the expected value ev, the standard deviation σ, and the third ratio p in the first histogram f (p).

pv算出部38bは、ヒストグラム算出部35によって求められた第3比pの第1ヒストグラムf(p)に基づいて、ノイズ成分比pvを算出し、この算出したノイズ成分比pvを血中酸素飽和度算出部39および脈拍数算出部40の各部へ通知するものである。   The pv calculation unit 38b calculates a noise component ratio pv based on the first histogram f (p) of the third ratio p obtained by the histogram calculation unit 35, and uses the calculated noise component ratio pv for blood oxygen saturation. It notifies to each part of the degree calculation part 39 and the pulse rate calculation part 40.

pv算出部38bは、より具体的には、本実施形態では、ヒストグラム算出部35によって求められた第1ヒストグラムf(p)と統計パラメータ算出部37aで算出された所定の統計パラメータとに基づいて第2ヒストグラムh2(p)を求め、この第2ヒストグラムh2(p)を用いてノイズ成分比pvを算出する。本実施形態では、p≦期待値evでは、h2(p)=0と定義され、p>期待値evでは、h2(p)={g2(p)−f(p)}/∫{g2(p)−f(p)}dpと定義され、g2(p)−f(p)≦0では、h2(p)=0と定義され、g2(p)−f(p)>0では、h2(p)={g2(p)−f(p)}/∫{g2(p)−f(p)}dpと定義される。すなわち、
h2(p)=0;p≦期待値evまたはg2(p)−f(p)≦0
h2(p)={g2(p)−f(p)}/∫{g2(p)−f(p)}dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。ただし、p≦evおよびg1(p)−f(p)≦0のpの領域を除く);p>期待値evまたはg2(p)−f(p)>0 ・・・(A4)
そして、ノイズ成分比pvは、次式A5によって算出される。
ノイズ成分比pv=∫p×h2(p)dp ・・・(A5)
ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。
More specifically, in this embodiment, the pv calculator 38b is based on the first histogram f (p) obtained by the histogram calculator 35 and the predetermined statistical parameter calculated by the statistical parameter calculator 37a. A second histogram h2 (p) is obtained, and the noise component ratio pv is calculated using the second histogram h2 (p). In this embodiment, when p ≦ expected value ev, h2 (p) = 0 is defined, and when p> expected value ev, h2 (p) = {g2 (p) −f (p)} / ∫ {g2 ( p) −f (p)} dp, h2 (p) = 0 is defined when g2 (p) −f (p) ≦ 0, and h2 when g2 (p) −f (p)> 0. (P) = {g2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp. That is,
h2 (p) = 0; p ≦ expected value ev or g2 (p) −f (p) ≦ 0
h2 (p) = {g2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp (where the integral ∫ is the maximum from the minimum class value a in the plurality of classes. Up to the value b of the class of ## EQU2 ## except for p ≦ ev and g1 (p) −f (p) ≦ 0, where p> expected value ev or g2 (p) −f (p )> 0 (A4)
The noise component ratio pv is calculated by the following equation A5.
Noise component ratio pv = ∫p × h2 (p) dp (A5)
However, the integral ∫ is executed from the minimum class value a to the maximum class value b in the plurality of classes.

次に、第2実施形態の生体情報測定装置Sbの動作について説明する。図5は、第2実施形態における生体情報測定装置の動作を示すフローチャートである。図6は、第2実施形態における生体情報測定装置の動作を説明するための図である。図6の比較的太い実線は、第3比pのヒストグラムf(p)を示し、破線は、第3比pに関するヒストグラムf(p)の移動平均g2(p)を示し、比較的細い実線は、関数h2(p)を表す。   Next, the operation of the biological information measuring device Sb of the second embodiment will be described. FIG. 5 is a flowchart showing the operation of the biological information measuring apparatus in the second embodiment. FIG. 6 is a diagram for explaining the operation of the biological information measuring apparatus according to the second embodiment. The relatively thick solid line in FIG. 6 indicates the histogram f (p) of the third ratio p, the broken line indicates the moving average g2 (p) of the histogram f (p) with respect to the third ratio p, and the relatively thin solid line indicates , Represents the function h2 (p).

本実施形態の生体情報測定装置Sbでは、第1実施形態の生体情報測定装置Saと同様に、例えば、図略の電源スイッチの投入や電源スイッチの投入後に図略の測定開始スイッチの投入によって、測定対象である生体の生体情報の測定が開始され、測定データが取得され、第1ヒストグラムf(p)が作成される。例えば、現在時刻tiから過去の時刻ti−(N+1)までの時間範囲の測定データについて、図6に実線で示すような第3比pの第1ヒストグラムf(p)が求められる。   In the biological information measuring device Sb of the present embodiment, similarly to the biological information measuring device Sa of the first embodiment, for example, by turning on a power switch (not shown) or turning on a measurement start switch (not shown) after turning on the power switch, Measurement of biological information of a living body to be measured is started, measurement data is acquired, and a first histogram f (p) is created. For example, for the measurement data in the time range from the current time ti to the past time ti− (N + 1), the first histogram f (p) having the third ratio p as shown by the solid line in FIG. 6 is obtained.

この第3比pの第1ヒストグラムf(p)が作成されると、図5において、まず、ステップS31では、演算制御部3bの統計パラメータ算出部37は、ステップS11と同様に、所定の範囲、例えば、第1ヒストグラムf(p)のクラスにおける最小のクラスの値aから最大のクラスの値bまでの範囲で、第3比pの第1ヒストグラムf(p)に対する期待値evおよび標準偏差σをそれぞれ算出する。   When the first histogram f (p) of the third ratio p is created, in FIG. 5, first, in step S31, the statistical parameter calculation unit 37 of the calculation control unit 3b performs a predetermined range as in step S11. For example, the expected value ev and the standard deviation for the first histogram f (p) of the third ratio p in the range from the minimum class value a to the maximum class value b in the class of the first histogram f (p). Each σ is calculated.

続いて、ステップS32では、演算制御部3bは、ステップS12と同様に、この求めた標準偏差σが予め設定された所定値σth以下であるか否かを判断する(σ≦σthか?)。   Subsequently, in step S32, the arithmetic control unit 3b determines whether or not the obtained standard deviation σ is equal to or less than a predetermined value σth set in advance (step σ ≦ σth?), As in step S12.

このステップS32の判断の結果、標準偏差σが所定値σth以下である場合(Yes)には、ステップS42の処理が実行され、一方、このステップS32の判断の結果、標準偏差σが所定値σth以下ではない場合(No)には、ステップS33の処理が実行される。   If the standard deviation σ is equal to or smaller than the predetermined value σth as a result of the determination in step S32 (Yes), the process of step S42 is executed, while the standard deviation σ is determined to be the predetermined value σth as a result of the determination in step S32. If not (No), the process of step S33 is executed.

ステップS33では、統計パラメータ算出部37bは、第3比pに関する第1ヒストグラムf(p)の移動平均g2(p)=∫f(p)dp(ただし、積分∫は、現在時点から2△Pの範囲(p−△pからp+△pまで)で実行される)を求める。例えば、上述の例では、図6に実線で示す第3比pの第1ヒストグラムf(p)に対し、図6に破線で示すg2(p)が求められる。   In step S33, the statistical parameter calculator 37b moves the moving average g2 (p) = ∫f (p) dp of the first histogram f (p) with respect to the third ratio p (where the integral ∫ is 2ΔP from the current time point). (Executed in the range from p−Δp to p + Δp). For example, in the above example, g2 (p) indicated by a broken line in FIG. 6 is obtained with respect to the first histogram f (p) having the third ratio p indicated by a solid line in FIG.

続いて、ステップS34では、演算制御部3bのpv算出部38bは、第3比pが期待値ev以上であるか否かを判断する(p≧evか?)。このステップS34の判断の結果、第3比pが期待値ev以上ではない場合(No、p<ev)には、pv算出部38bは、ステップS36を実行し、当該第3比pに対するh2(p)=0とし、ステップS39の処理を実行する。一方、このステップS34の判断の結果、第3比pが期待値ev以上である場合(Yes、p≧ev)には、pv算出部38bは、ステップS35を実行し、g2(p)−f(p)が0以上であるか否かを判断する(g2(p)−f(p)≧0か?)。   Subsequently, in step S34, the pv calculation unit 38b of the arithmetic control unit 3b determines whether or not the third ratio p is equal to or larger than the expected value ev (p ≧ ev?). As a result of the determination in step S34, when the third ratio p is not equal to or greater than the expected value ev (No, p <ev), the pv calculation unit 38b executes step S36, and h2 ( p) = 0, and the process of step S39 is executed. On the other hand, if the result of determination in step S34 is that the third ratio p is equal to or greater than the expected value ev (Yes, p ≧ ev), the pv calculation unit 38b executes step S35, and g2 (p) −f It is determined whether (p) is 0 or more (whether g2 (p) −f (p) ≧ 0).

このステップS35の判断の結果、g2(p)−f(p)が0以上ではない場合(No、g2(p)−f(p)<0)には、pv算出部38bは、ステップS37を実行し、当該第3比pに対するh2(p)=0とし、ステップS39の処理を実行する。一方、このステップS35の判断の結果、g2(p)−f(p)が0以上である場合(Yes、g2(p)−f(p)≧0)には、pv算出部38bは、ステップS38を実行し、当該第3比pに対するh2(p)={g2(p)−f(p)}/∫{g2(p)−f(p)}dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。ただし、p≦evおよびg1(p)−f(p)≦0のpの領域を除く)を求める。   If g2 (p) -f (p) is not equal to or greater than 0 as a result of the determination in step S35 (No, g2 (p) -f (p) <0), the pv calculation unit 38b performs step S37. And h2 (p) = 0 for the third ratio p, and the process of step S39 is executed. On the other hand, if g2 (p) −f (p) is equal to or greater than 0 as a result of the determination in step S35 (Yes, g2 (p) −f (p) ≧ 0), the pv calculation unit 38b S38 is executed, h2 (p) = {g2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp for the third ratio p (where integral ∫ Is performed from the minimum class value a to the maximum class value b in the class (except for the p region where p ≦ ev and g1 (p) −f (p) ≦ 0).

このようなステップS34からステップS38までの各処理が、a≦p≦bの範囲における第3比pのそれぞれについて実行され、第2ヒストグラムh2(p)が求められる。例えば、上述の例では、図6に比較的太い実線で示す第3比pの第1ヒストグラムf(p)に対し、図6に比較的細い実線で示す第2ヒストグラムh2(p)が求められる。   Each processing from step S34 to step S38 is executed for each of the third ratios p in the range of a ≦ p ≦ b, and the second histogram h2 (p) is obtained. For example, in the above-described example, the second histogram h2 (p) indicated by the relatively thin solid line in FIG. 6 is obtained with respect to the first histogram f (p) having the third ratio p indicated by the relatively thick solid line in FIG. .

そして、ステップS39では、pv算出部38bは、ノイズ成分比pv=∫p×h2(p)dp(ただし、積分∫は、前記複数のクラスにおける最小のクラスの値aから最大のクラスの値bまでで実行される。)を求める。   In step S39, the pv calculator 38b determines that the noise component ratio pv = ∫p × h2 (p) dp (where the integral ∫ is the minimum class value a to the maximum class value b in the plurality of classes). Will be executed by.)

続いて、ステップS40では、脈拍数算出部40は、ステップS20と同様に、上述の各処理によって求めたノイズ成分比pvを用いて、脈波の信号S(i)=R(i)−pv×IR(i)を求め、この信号S(i)の周期T_periodを求める。   Subsequently, in step S40, the pulse rate calculation unit 40 uses the noise component ratio pv obtained by each of the above-described processes, similarly to step S20, to detect a pulse wave signal S (i) = R (i) −pv. XIR (i) is obtained, and the period T_period of the signal S (i) is obtained.

続いて、ステップS41では、脈拍数算出部40は、ステップS21と同様に、このように求めた信号S(i)の周期T_periodで60を割り算することによって脈拍数(=60/周期T_period)を求める。   Subsequently, in step S41, the pulse rate calculation unit 40 divides 60 by the period T_period of the signal S (i) obtained in this way, similarly to step S21, thereby dividing the pulse rate (= 60 / period T_period). Ask.

一方、ステップS42では、脈拍数算出部40は、ステップS22と同様に、第2比IR(i)の周期T_periodを求め、続いて、ステップS41を実行し、このように求めた第2比IR(i)の周期T_periodで60を割り算することによって脈拍数(=60/周期T_period)を求める。   On the other hand, in step S42, the pulse rate calculation unit 40 obtains the cycle T_period of the second ratio IR (i) in the same manner as in step S22, and subsequently executes step S41 to obtain the second ratio IR thus obtained. The pulse rate (= 60 / period T_period) is obtained by dividing 60 by the period T_period of (i).

そして、演算制御部3bは、このような各処理によって求めた脈拍数を表示部5の脈拍数表示部51に表示する。   Then, the arithmetic control unit 3 b displays the pulse rate obtained by each of these processes on the pulse rate display unit 51 of the display unit 5.

なお、上述の処理において、期待値evを求める際の積分範囲とノイズ成分比pvを求める際の積分範囲とは、同一であっても異なっていてもよい。   In the above-described processing, the integration range for obtaining the expected value ev and the integration range for obtaining the noise component ratio pv may be the same or different.

また、演算制御部3bの血中酸素飽和度算出部39は、ノイズ成分比pvに基づいて静脈血の酸素飽和度SvOを求めるとともに、ノイズ成分比pvに基づいて信号成分比paを求め、この求めた信号成分比paに基づいて動脈血の酸素飽和度SpOを求める。そして、演算制御部3bは、これら求めた酸素飽和度SvO、SpOを表示部5の血中酸素飽和度表示部52に表示する。 The blood oxygen saturation calculation unit 39 of the calculation control unit 3b calculates the venous blood oxygen saturation SvO 2 based on the noise component ratio pv, and calculates the signal component ratio pa based on the noise component ratio pv. Based on the obtained signal component ratio pa, the oxygen saturation SpO 2 of arterial blood is obtained. Then, the arithmetic control unit 3 b displays the obtained oxygen saturation SvO 2 and SpO 2 on the blood oxygen saturation display unit 52 of the display unit 5.

また、演算制御部3bの信頼度算出部41は、所定の信頼度を求め、演算制御部3bは、この求めた信頼度を表示部5の信頼度表示部53に表示する。   Further, the reliability calculation unit 41 of the calculation control unit 3 b obtains a predetermined reliability, and the calculation control unit 3 b displays the obtained reliability on the reliability display unit 53 of the display unit 5.

このように本実施形態の生体情報測定装置Sbでは、測定データR、IRによる第3比pの第1ヒストグラムf(p)が求められ、その第1ヒストグラムf(p)の統計パラメーラとして期待値evが求められ、その期待値evに基づく所定範囲における測定データに基づいて所定の生体情報が求められる。このように本実施形態の生体情報測定装置Sbは、より簡易な処理によって生体情報を求めることができ、その結果、消費電力を低減することができる。   As described above, in the biological information measuring device Sb of the present embodiment, the first histogram f (p) of the third ratio p based on the measurement data R and IR is obtained, and an expected value as a statistical parameter of the first histogram f (p). ev is obtained, and predetermined biological information is obtained based on measurement data in a predetermined range based on the expected value ev. As described above, the biological information measuring device Sb according to the present embodiment can obtain biological information through simpler processing, and as a result, can reduce power consumption.

そして、本実施形態では、第3比pの第1ヒストグラムf(p)における期待値evおよび移動平均g2(p)を所定の統計パラメータとして求め、この第1ヒストグラムf(p)によるpの期待値evおよび移動平均g2(p)に基づき生体の脈拍数を求め、さらに加えて血中酸素飽和度を求める生体情報測定装置Sbが提供される。   In this embodiment, the expected value ev and the moving average g2 (p) in the first histogram f (p) of the third ratio p are obtained as predetermined statistical parameters, and the expectation of p based on the first histogram f (p) A biological information measuring device Sb is provided that calculates the pulse rate of the living body based on the value ev and the moving average g2 (p), and further calculates the blood oxygen saturation.

なお、上述の第2実施形態では、統計パラメータ算出部37bは、前記所定の統計パラメータとして第3比pの第1ヒストグラムf(p)における期待値evおよびその移動平均g2(p)を求め、pv算出部38bは、第3比pが期待値evよりも大きい場合、または、g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムh2(p)として求め、この第2ヒストグラムh2(p)に基づいてノイズ成分比pvを求め、脈拍数算出部40は、このノイズ成分比pvに基づいて生体の脈拍数を求めたが、統計パラメータ算出部37bは、第3比pの第1ヒストグラムf(p)における移動平均g2(p)を求め、前記移動平均g2(p)のピークを与える第3比p_pを前記所定の統計パラメータとして求め、pv算出部38bは、第3比pが移動平均g2(p)のピークを与える第3比p_pよりも大きい場合、または、g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムh2(p)として求め、この第2ヒストグラムh2(p)に基づいてノイズ成分比pvを求め、脈拍数算出部40は、このノイズ成分比pvに基づいて生体の脈拍数を求めてもよい。   In the second embodiment described above, the statistical parameter calculation unit 37b obtains the expected value ev in the first histogram f (p) of the third ratio p and its moving average g2 (p) as the predetermined statistical parameter, When the third ratio p is larger than the expected value ev or when g2 (p) −f (p) is larger than 0, the pv calculation unit 38b {g2 (p) −f (p)} / ∫ { g2 (p) −f (p)} dp is obtained as the second histogram h2 (p), the noise component ratio pv is obtained based on the second histogram h2 (p), and the pulse rate calculating unit 40 Although the pulse rate of the living body is obtained based on the component ratio pv, the statistical parameter calculation unit 37b obtains the moving average g2 (p) in the first histogram f (p) of the third ratio p, and the moving average g2 (p ) Give a peak The ratio p_p is obtained as the predetermined statistical parameter, and the pv calculation unit 38b determines that the third ratio p is larger than the third ratio p_p that gives the peak of the moving average g2 (p), or g2 (p) −f ( {g2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp when p) is larger than 0 is obtained as the second histogram h2 (p), and this second histogram h2 Based on (p), the noise component ratio pv may be obtained, and the pulse rate calculator 40 may obtain the pulse rate of the living body based on the noise component ratio pv.

このような構成によっても生体情報測定装置Sは、より簡易な処理によって生体情報を求めることができ、その結果、消費電力を低減することができる。そして、この構成では、第3比pの第1ヒストグラムf(p)における移動平均g2(p)およびこの移動平均g2(p)のピークを与える第3比p_pを所定の統計パラメータとして求め、この第1ヒストグラムf(p)の移動平均g2(p)に基づき生体の脈拍数を求め、さらに加えて血中酸素飽和度を求める生体情報測定装置Sbが提供される。   Even with such a configuration, the biological information measuring device S can obtain biological information through simpler processing, and as a result, power consumption can be reduced. In this configuration, the moving average g2 (p) in the first histogram f (p) of the third ratio p and the third ratio p_p that gives the peak of the moving average g2 (p) are obtained as predetermined statistical parameters. A biological information measuring device Sb is provided that calculates the pulse rate of the living body based on the moving average g2 (p) of the first histogram f (p) and further calculates the blood oxygen saturation.

本発明を表現するために、上述において図面を参照しながら実施形態を通して本発明を適切且つ十分に説明したが、当業者であれば上述の実施形態を変更および/または改良することは容易に為し得ることであると認識すべきである。したがって、当業者が実施する変更形態または改良形態が、請求の範囲に記載された請求項の権利範囲を離脱するレベルのものでない限り、当該変更形態または当該改良形態は、当該請求項の権利範囲に包括されると解釈される。   In order to express the present invention, the present invention has been properly and fully described through the embodiments with reference to the drawings. However, those skilled in the art can easily change and / or improve the above-described embodiments. It should be recognized that this is possible. Therefore, unless the modifications or improvements implemented by those skilled in the art are at a level that departs from the scope of the claims recited in the claims, the modifications or improvements are not covered by the claims. To be construed as inclusive.

Sa、Sb 生体情報測定装置
1 センサ部
3a、3b 演算制御部
5 表示部
11 第1センサ部
12 第2センサ部
35 ヒストグラム算出部
36 比算出部
37a、37b 統計パラメータ算出部
38a、38b pv算出部
39 血中酸素飽和度算出部
40 脈拍数算出部
Sa, Sb Biological information measuring device 1 Sensor unit 3a, 3b Operation control unit 5 Display unit 11 First sensor unit 12 Second sensor unit 35 Histogram calculation unit 36 Ratio calculation unit 37a, 37b Statistical parameter calculation unit 38a, 38b pv calculation unit 39 Blood oxygen saturation calculator 40 Pulse rate calculator

Claims (7)

測定対象の生体における所定の生理的現象を測定して測定データを出力する測定部と、
前記測定部によって測定された測定データに基づいて第1ヒストグラムを求めるヒストグラム算出部と、
前記ヒストグラム算出部によって求められた第1ヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出部と、
前記ヒストグラム算出部によって求められた第1ヒストグラムと前記統計パラメータ算出部で算出された前記所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求める生体情報算出部とを備えること
を特徴とする生体情報測定装置。
A measurement unit that measures a predetermined physiological phenomenon in the living body to be measured and outputs measurement data; and
A histogram calculation unit for obtaining a first histogram based on the measurement data measured by the measurement unit;
A statistical parameter calculator for obtaining a predetermined statistical parameter in the first histogram obtained by the histogram calculator;
A second histogram is obtained based on the first histogram obtained by the histogram calculation unit and the predetermined statistical parameter calculated by the statistical parameter calculation unit, and the predetermined biological information is obtained based on the second histogram. A biological information measuring device comprising: a biological information calculating unit.
前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、
前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、
前記統計パラメータ算出部は、前記所定の統計パラメータとして前記第1ヒストグラムf(p)における期待値evを求め、
前記生体情報算出部は、g1(p)=f(2×ev−p)とする際に、前記第3比pが前記期待値evよりも大きい場合、または、前記g1(p)−f(p)が0より大きい場合における{g1(p)−f(p)}/∫{g1(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求めること
を特徴とする請求項1に記載の生体情報測定装置。
The measurement unit irradiates the living tissue of the living body with first light having a predetermined first wavelength, receives the first light transmitted or reflected through the living tissue of the living body, and outputs first measurement data. Second measurement data is received by irradiating the biological tissue of the living body with second light having a predetermined second wavelength different from the first wavelength, and receiving or transmitting the second light transmitted or reflected through the biological tissue of the biological body. A sensor unit including a second sensor unit that outputs
The histogram calculation unit calculates a histogram f (p) of a third ratio p between a first ratio of the AC component to the DC component in the first measurement data and a second ratio of the AC component to the DC component in the second measurement data. Obtained as the first histogram,
The statistical parameter calculation unit obtains an expected value ev in the first histogram f (p) as the predetermined statistical parameter,
When the third ratio p is larger than the expected value ev when g1 (p) = f (2 × ev−p), or the g1 (p) −f ( {g1 (p) −f (p)} / ∫ {g1 (p) −f (p)} dp when p) is larger than 0 is obtained as the second histogram, and the living body is calculated based on the second histogram. The biological information measuring device according to claim 1, wherein the pulse rate is calculated.
前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、
前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、
前記統計パラメータ算出部は、前記所定の統計パラメータとして前記第1ヒストグラムf(p)における期待値evおよび移動平均g2(p)を求め、
前記生体情報算出部は、前記第3比pが前記期待値evよりも大きい場合、または、前記g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求めること
を特徴とする請求項1に記載の生体情報測定装置。
The measurement unit irradiates the living tissue of the living body with first light having a predetermined first wavelength, receives the first light transmitted or reflected through the living tissue of the living body, and outputs first measurement data. Second measurement data is received by irradiating the biological tissue of the living body with second light having a predetermined second wavelength different from the first wavelength, and receiving or transmitting the second light transmitted or reflected through the biological tissue of the biological body. A sensor unit including a second sensor unit that outputs
The histogram calculation unit calculates a histogram f (p) of a third ratio p between a first ratio of the AC component to the DC component in the first measurement data and a second ratio of the AC component to the DC component in the second measurement data. Obtained as the first histogram,
The statistical parameter calculation unit obtains an expected value ev and a moving average g2 (p) in the first histogram f (p) as the predetermined statistical parameter,
The biometric information calculation unit {g2 (p) −f (p) when the third ratio p is larger than the expected value ev or when the g2 (p) −f (p) is larger than 0. } / 情報 {g2 (p) −f (p)} dp is obtained as the second histogram, and the pulse rate of the living body is obtained based on the second histogram. measuring device.
前記測定部は、所定の第1波長の第1光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第1光を受光して第1測定データを出力する第1センサ部と、前記第1波長と異なる所定の第2波長の第2光を前記生体の生体組織に照射し、前記生体の生体組織を透過または反射した第2光を受光して第2測定データを出力する第2センサ部とを備えるセンサ部を備え、
前記ヒストグラム算出部は、前記第1測定データにおける直流成分に対する交流成分の第1比と前記第2測定データにおける直流成分に対する交流成分の第2比との第3比pのヒストグラムf(p)を前記第1ヒストグラムとして求め、
前記統計パラメータ算出部は、前記第1ヒストグラムf(p)における移動平均g2(p)を求め、前記移動平均g2(p)のピークを与える第3比p_pを前記所定の統計パラメータとして求め、
前記生体情報算出部は、前記第3比pが前記移動平均g2(p)のピークを与える第3比p_pよりも大きい場合、または、前記g2(p)−f(p)が0より大きい場合における{g2(p)−f(p)}/∫{g2(p)−f(p)}dpを前記第2ヒストグラムとして求め、この第2ヒストグラムに基づいて前記生体の脈拍数を求めること
を特徴とする請求項1に記載の生体情報測定装置。
The measurement unit irradiates the living tissue of the living body with first light having a predetermined first wavelength, receives the first light transmitted or reflected through the living tissue of the living body, and outputs first measurement data. Second measurement data is received by irradiating the biological tissue of the living body with second light having a predetermined second wavelength different from the first wavelength, and receiving or transmitting the second light transmitted or reflected through the biological tissue of the biological body. A sensor unit including a second sensor unit that outputs
The histogram calculation unit calculates a histogram f (p) of a third ratio p between a first ratio of the AC component to the DC component in the first measurement data and a second ratio of the AC component to the DC component in the second measurement data. Obtained as the first histogram,
The statistical parameter calculation unit obtains a moving average g2 (p) in the first histogram f (p), obtains a third ratio p_p that gives a peak of the moving average g2 (p) as the predetermined statistical parameter,
When the third ratio p is larger than a third ratio p_p that gives the peak of the moving average g2 (p), or when the g2 (p) −f (p) is larger than 0 {G2 (p) −f (p)} / ∫ {g2 (p) −f (p)} dp is obtained as the second histogram, and the pulse rate of the living body is obtained based on the second histogram. The living body information measuring device according to claim 1 characterized by things.
前記生体情報算出部は、前記第1ヒストグラムの標準偏差の値に応じて前記生体の脈拍数を求める算出式を変更すること
を特徴とする請求項2ないし請求項4のいずれか1項に記載の生体情報測定装置。
The said biological information calculation part changes the calculation formula which calculates | requires the pulse rate of the said biological body according to the value of the standard deviation of the said 1st histogram, The any one of Claim 2 thru | or 4 characterized by the above-mentioned. Biological information measuring device.
前記生体情報算出部は、前記所定の生体情報としてさらに血中酸素飽和度を求めること
を特徴とする請求項2ないし請求項5のいずれか1項に記載の生体情報測定装置。
The biological information measuring apparatus according to any one of claims 2 to 5, wherein the biological information calculation unit further calculates blood oxygen saturation as the predetermined biological information.
測定対象の生体における所定の生理的現象を測定して測定データを出力する測定工程と、
前記測定工程によって測定された測定データに基づいて第1ヒストグラムを求めるヒストグラム算出工程と、
前記ヒストグラム算出工程によって求められた第1ヒストグラムにおける所定の統計パラメータを求める統計パラメータ算出工程と、
前記ヒストグラム算出工程によって求められた第1ヒストグラムと前記統計パラメータ算出工程で算出された前記所定の統計パラメータとに基づいて第2ヒストグラムを求め、この第2ヒストグラムに基づいて前記所定の生体情報を求める生体情報算出工程とを備えること
を特徴とする生体情報測定方法。
A measurement process for measuring a predetermined physiological phenomenon in the measurement target and outputting measurement data;
A histogram calculation step for obtaining a first histogram based on the measurement data measured by the measurement step;
A statistical parameter calculating step for obtaining a predetermined statistical parameter in the first histogram obtained by the histogram calculating step;
A second histogram is obtained based on the first histogram obtained by the histogram calculating step and the predetermined statistical parameter calculated by the statistical parameter calculating step, and the predetermined biological information is obtained based on the second histogram. A biological information measuring method comprising: a biological information calculating step.
JP2010077211A 2010-03-30 2010-03-30 Biological information measuring apparatus and method Expired - Fee Related JP5531715B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010077211A JP5531715B2 (en) 2010-03-30 2010-03-30 Biological information measuring apparatus and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010077211A JP5531715B2 (en) 2010-03-30 2010-03-30 Biological information measuring apparatus and method

Publications (2)

Publication Number Publication Date
JP2011206286A true JP2011206286A (en) 2011-10-20
JP5531715B2 JP5531715B2 (en) 2014-06-25

Family

ID=44938071

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010077211A Expired - Fee Related JP5531715B2 (en) 2010-03-30 2010-03-30 Biological information measuring apparatus and method

Country Status (1)

Country Link
JP (1) JP5531715B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014024626A1 (en) * 2012-08-09 2014-02-13 コニカミノルタ株式会社 Biological information measuring device and pulse oximeter
JP2017516533A (en) * 2014-05-22 2017-06-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and apparatus for optical sensing of tissue variation with increased accuracy
JP2017176265A (en) * 2016-03-28 2017-10-05 富士ゼロックス株式会社 Apparatus and program for measuring biological information
JP2017213389A (en) * 2012-10-26 2017-12-07 ナイキ イノベイト シーブイ Athletic performance monitoring system using heart rate information
JP2018501883A (en) * 2014-12-31 2018-01-25 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. System for performing histogram analysis of time-based capnography signals and method of operation thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001078990A (en) * 1999-09-13 2001-03-27 Nippon Koden Corp Apparatus for measuring light absorber in blood
JP2002516130A (en) * 1998-05-26 2002-06-04 ネルコー・ピューリタン・ベネット・インコーポレイテッド Method and apparatus for estimating physiological parameters using transform
JP2007532207A (en) * 2004-04-16 2007-11-15 マシモ コーポレーション Non-invasive measurement method for the second heart sound component
JP2009034427A (en) * 2007-08-03 2009-02-19 Konica Minolta Sensing Inc Pulse oxymeter
WO2010001235A1 (en) * 2008-06-30 2010-01-07 Nellcor Puritan Bennett Ireland Signal processing systems and methods for determining slope using an origin point

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002516130A (en) * 1998-05-26 2002-06-04 ネルコー・ピューリタン・ベネット・インコーポレイテッド Method and apparatus for estimating physiological parameters using transform
JP2001078990A (en) * 1999-09-13 2001-03-27 Nippon Koden Corp Apparatus for measuring light absorber in blood
JP2007532207A (en) * 2004-04-16 2007-11-15 マシモ コーポレーション Non-invasive measurement method for the second heart sound component
JP2009034427A (en) * 2007-08-03 2009-02-19 Konica Minolta Sensing Inc Pulse oxymeter
WO2010001235A1 (en) * 2008-06-30 2010-01-07 Nellcor Puritan Bennett Ireland Signal processing systems and methods for determining slope using an origin point

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014024626A1 (en) * 2012-08-09 2014-02-13 コニカミノルタ株式会社 Biological information measuring device and pulse oximeter
JPWO2014024626A1 (en) * 2012-08-09 2016-07-25 コニカミノルタ株式会社 Biological information measuring device and pulse oximeter
US9642567B2 (en) 2012-08-09 2017-05-09 Konica Minolta, Inc. Biological information measurement device and pulse oximeter
JP2017213389A (en) * 2012-10-26 2017-12-07 ナイキ イノベイト シーブイ Athletic performance monitoring system using heart rate information
US10206589B2 (en) 2012-10-26 2019-02-19 Nike, Inc. Athletic performance monitoring system utilizing heart rate information
US11478156B2 (en) 2012-10-26 2022-10-25 Nike, Inc. Athletic performance monitoring system utilizing heart rate information
JP2017516533A (en) * 2014-05-22 2017-06-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and apparatus for optical sensing of tissue variation with increased accuracy
US10433738B2 (en) 2014-05-22 2019-10-08 Koninklijke Philips N.V. Method and apparatus for optical sensing of tissue variation at increased accuracy
JP2018501883A (en) * 2014-12-31 2018-01-25 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. System for performing histogram analysis of time-based capnography signals and method of operation thereof
US10869614B2 (en) 2014-12-31 2020-12-22 Koninklijke Philips N.V. System for performing histogram analysis of the time-based capnography signals and method of operation thereof
JP2017176265A (en) * 2016-03-28 2017-10-05 富士ゼロックス株式会社 Apparatus and program for measuring biological information

Also Published As

Publication number Publication date
JP5531715B2 (en) 2014-06-25

Similar Documents

Publication Publication Date Title
US10271748B2 (en) Patient monitor for determining microcirculation state
US8494786B2 (en) Exponential sampling of red and infrared signals
US8457706B2 (en) Estimation of a physiological parameter using a neural network
US20120165629A1 (en) Systems and methods of monitoring a patient through frequency-domain photo migration spectroscopy
US8649838B2 (en) Wavelength switching for pulse oximetry
US9247896B2 (en) Systems and methods for determining respiration information using phase locked loop
EP3307162B1 (en) Pulse oximetry
US9078609B2 (en) Extraction of physiological measurements from a photoplethysmograph (PPG) signal
US20100081899A1 (en) System and Method for Photon Density Wave Pulse Oximetry and Pulse Hemometry
US20130137946A1 (en) Medical device with conditional power consumption
JP2004202190A (en) Biological information measuring device
JP2008188216A (en) Biological information measurement apparatus
JP2008532682A (en) Method and apparatus for determining blood perfusion in a body part
WO2008112582A2 (en) Detection of oximetry sensor sites based on waveform characteristics
JP5531715B2 (en) Biological information measuring apparatus and method
US8870783B2 (en) Pulse rate determination using Gaussian kernel smoothing of multiple inter-fiducial pulse periods
US20140275825A1 (en) Methods and systems for light signal control in a physiological monitor
US20140221852A1 (en) Systems and methods for determining respiration information using frequency demodulation
US8840562B2 (en) Signal processing warping technique
US20140066782A1 (en) System and method for determining a resting heart rate of an individual
JP5353790B2 (en) Biological information measuring apparatus and method
US20140187884A1 (en) Systems and methods for ensemble averaging in pulse oximetry
US9888871B2 (en) Methods and systems for determining a venous signal using a physiological monitor
JP5527277B2 (en) Signal processing apparatus, signal processing method, and biological information measuring apparatus
US9560978B2 (en) Systems and methods for determining respiration information from a physiological signal using amplitude demodulation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20121115

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20130418

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130711

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130723

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130910

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140325

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140407

R150 Certificate of patent or registration of utility model

Ref document number: 5531715

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees