JPH07113743A - 生体信号処理方法 - Google Patents
生体信号処理方法Info
- Publication number
- JPH07113743A JPH07113743A JP5280293A JP28029393A JPH07113743A JP H07113743 A JPH07113743 A JP H07113743A JP 5280293 A JP5280293 A JP 5280293A JP 28029393 A JP28029393 A JP 28029393A JP H07113743 A JPH07113743 A JP H07113743A
- Authority
- JP
- Japan
- Prior art keywords
- power spectrum
- blood flow
- light
- tissue
- red blood
- 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
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
する。 【構成】 生体に対しレーザー光を照射して、反射散乱
光を収集し、前記反射散乱光からパワースペクトルを検
出し、多重回帰処理を施す。
Description
小循環血流計測等の為の生体情報を処理する方法に関す
る。
ックダイアグラムを示す。図1(2)は、図1(1)の
(01)部分を拡大したものである。図1(1)に於て、(1
1)はプローブであり、送光用ファイバー(12)と受光用フ
ァイバー(13)を束ねたものである。プローブ(11)は図1
(2)で示す様に生体組織(MA)の表面に接触して使用
される。送光用ファイバー(12)はレーザー光源(14)と接
続し、受光用ファイバー(13)はフォトダイオード(15)と
接続している。フォトダイオード(15)は更にアンプ(16)
に接続される。アンプ(16)の出力は、入力信号に対しA
(ω)演算処理を施すバンドパスフィルタ(17)、入力信号
に対し〈I(t)〉処理を施すローパスフィルタ(18)と接
続し、バンドパスフィルタ(17)は更に入力信号A(ω)に
対し、
し、P(ω)処理を施す2乗演算回路(21)に接続してい
る。ラティスフィルタ(20)の出力は入力信号
2乗回路(22)の出力は入力信号(ωP(ω))に対し、
力は、入力信号(P(ω))に対し、
(18)の出力は入力信号〈I(t)〉に対し、〈I(t)〉2
処理を施す2乗回路(19)と接続する。積分回路(23)と2
乗回路(19)との出力は、これら2入力に対し
の出力は、出力端(MA)並びに除算回路(27)の一端と
接続する。2乗回路(19)の出力と積分回路(24)の出力
は、これら2入力に対し、
の出力は出力端(FL)並びに除算回路(27)の他端と接
続する。除算回路(27)は、これら2入力に対し、
接続する。レーザー光源(14)から出力されたレーザー光
は、送光用ファイバー(12)を通して組織(MA)に照射
される。図1(2)で示す組織中の赤血球(MB)等で
散乱された光の一部が受光用の光ファイバー(13)を通じ
てフォトダイオード(15)に導かれる。フォトダイオード
(15)で光電変換され、アンプ(16)で増幅された信号は、
上述した信号処理回路で処理され組織血流量に相当する
量を出力端(FL)、組織血液量に相当する量を出力端
(MA)、血流速度に相当する量を出力端(VE)に出
力する。
で測定される量は、被測定体積内の平均値である。つま
り、血流速度は、その組織内を流れる血流の平均速度で
ある。従って、組織中には毛細血管と細動静脈が存在す
るが、これらを分けて測定できない。
定生体体積内の血流速度を分けて測定するために、アン
プ出力の信号をFFT処理し、これを多重回帰分析する
ことで血流速度分布を求めることを実現した。
光手段を供えたプローブと前記受光手段で得られた散乱
光を光電変換して電気信号に変換する手段及び前記電気
信号を増幅する増幅手段、パワースペクトル検出手段、
多重回帰分析手段、より構成される。レーザー光出力手
段は例えば、半導体レーザーを用いたもの、受光手段は
例えば、受光ダイオード等の光電変換素子を用いたもの
であり、レーザー光出力手段と受光手段は図1の様に一
体化されたプローブを形成されるものであるが、その
他、別体である場合もある。又、プローブと生体皮膚表
面とは、接触、非接触のどちらの形態でもかまわない。
パワースペクトル検出手段、多重回帰分析手段は、1つ
乃至複数のマイクロコンピュータ、アナログ演算処理手
段あるいはDSP等で構成された信号処理手段で実施さ
れるものであるが、具体的処理過程について次に説明す
る。アンプ出力の信号をFFT処理すると、得られる信
号は図2のようになる。これは人の前腕の皮膚血流を測
定した場合である。スペクトル波形がピークを持たない
理由は、光子が赤血球に衝突する前に静止組織によって
何回も散乱されるため、入射光ベクトルと赤血球の速度
ベクトルの方向の関係がランダムになるためである。こ
のパワースペクトル信号を信号処理手段で多重回帰して
分解することで、それぞれの速度の要素に分けることが
できる。次に多重回帰について詳細に説明する。 モデル理論:500から900nm波長域では生体組織
は強い散乱体であるので、組織に入射した光は拡散され
る。一般に組織中の赤血球数密度は低いので、送光点と
受光点間の距離が短ければ、ほとんどの光子は静止組織
のみによって散乱され、わずかな光子がその散乱過程で
1回、または数回散乱される。静止組織によって散乱さ
れた光は周波数変調されないが、流れている赤血球によ
って散乱された光は周波数変調される。静止組織からの
散乱光の割合が多いため、散乱光の一部を受光して信号
処理するときにはヘテロダイン法として扱うことができ
る。光子は赤血球に衝突する前に静止組織によって何回
か散乱されるため、入射光ベクトルの方向はランダムと
して考えられる。赤血球からの散乱光の自己相関関数、
赤血球の散乱光強度で、θは散乱角度である。F(q,
t)は散乱関数である。赤血球の大きさが、波長よりも
かなり大きい場合には、|f(θ)|2は次式で近似でき
る。 |f(θ)|2=(α/π)σsexp(−αθ2)・・・・・(2) ここでσsは赤血球の散乱断面積、αは散乱光分布から
得られる係数である。赤血球の速度分布がガウス分布だ
と仮定すると、F(q,t)は、 F(q,t)=exp(−q2〈V2〉t2/6)・・・・・(3) として表される。ここでqは散乱ベクトルで、q=(4
π/λ)sin(θ/2)の大きさを持つ。λは波長で、
Vは赤血球の速度である。ここで速度分布が2つあり、
それぞれが独立であるならば、F(q,t)はそれぞれ
の合成として考えられるので、 Ft(q,t)=A exp(−q2〈V1 2〉t2/6)+B exp(−q2〈V2 2〉t2/6) ・・・・・(4) として表すことができる。AとBはそれぞれV1とV2k
速度を持つ赤血球の数密度に依存する係数で、A+B=
1である。したがって、この場合(1)式は、
〈V1 2〉t2/λ2,τ2=(2/3)π2〈V2 2〉t2/
λ2である。この相関関数のパワースペクトルS(ω)は
次式で表すことができる。
(5),(6)式より、 S(ω)=(〈N〉/π)[A exp(−aω)+B exp(−bω)]・・(7) となる。ここで、a=[(3αλ2)/(2π
2〈V1 2〉)]1/2,b=[(3αλ2)/(2π
2〈V1 2〉)]1/2である。ここで最小自乗法を用いて回
帰すれば、A,B,a,bが求まる。1/aと1/bは
平均周波数であり、速度に比例する。また、それぞれの
速度を持つ赤血球からの散乱光強度(IIPS)は、
(πb)である。IIPSは赤血球数密度に比例するの
で、A/aとB/bの比から、それぞれの速度成分を持
つ赤血球量の割合が分かる。
かめた。実験は図3に示す組織血流モデルを用いて行っ
た。図3(2)は、図3(1)の、プローブ(11)周辺の
断面を示している。このモデルは、血管の代わりである
2本のマイクロチューブ(211)(212)と、静止組織の代わ
りである2枚のポリアセタール板(213)(214)から構成さ
れている。ポリアセタール板(213)の厚さは、約0.3m
mである。赤血球の代わりとして、マイクロチューブ(21
1)(212)の中にポリスチレン粒子(215)をマイクロフィー
ダー(216)(217)を用いて流した。粒子の直径は約2μm
である。送光用ファイバー(12)と受光用ファイバー(13)
を束ねたプローブ(11)をポリアセタール板(213)上に接
触させる。送光用ファイバー(12)はレーザー光源(14)、
受光用ファイバー(13)は、フォトダイオード(15)と接続
する。フォトダイオード(15)は、増幅器(16)と接続す
る。2本のチューブ内(211)(212)の粒子速度を0.3から
3.4mm/sの速度で階段的に変えて、パワースペクトル
を測定した。パワースペクトルの測定にはFFTアナラ
イザーを用いた。FFTアナライザーは、増幅器(16)の
出力端(30)と接続される。解析は自作のプログラムで行
った。このとき、パワースペクトルをそのまま回帰する
と、高周波成分での誤差が低く見積もられてしまうた
め、パワースペクトルに周波数の重みをかけたデータと
して解析した。(f×S(f)) 図4(a)は、実際にヒトの人差し指から図3(1)で
示すプローブ(11)を用いて血流を測定し、そのパワ
ースペクトルを算出して描いたグラフであり、図4
(b)は、上述した実験用モデルシステムに於てポリス
チレン粒子の流れをプローブ(11)を用いて測定し、その
パワースペクトルを算出して描いたグラフである。尚、
本モデルの妥当性は既に発表されている。(JAPANESE J
ournal ofAPPLIED PHYSICS Vol.32, pp.2177-2182 (199
3) ジャパニーズジャーナル オブアプライド フィジ
クス) 図5に、設定した速度と、回帰した後のパワースペクト
ルの平均周波数の関係を示す。(a)は粒子濃度が、0.0
2%、(b)は粒子濃度が、0.04%である。パワースペ
クトルの平均周波数は、粒子速度に比例する。自然対数
の多重回帰は何本でも可能であるが、本実験では2本で
回帰したときの残差2乗和の値が、1本で回帰したとき
の残差2乗和とあまり変わらなければ(本実験では0.7
倍以上小さくならなければ)、1本の回帰結果を用い
た。したがって、2本のチューブ(211)(212)内の速度が
0.3と0.6mm/sなど、あまり差がない場合には、2種類
の速度として回帰できず、1種類の速度として回帰され
た。この実験結果からは、高い相関が得られた。本実験
より、2つの速度にある程度差があれば、本方法で分離
して測定できると考えられた。実際に前腕と指の血流を
測定し、本方法で回帰した結果を表1に示す。表1は、
前腕と指の血流を回帰したときのパワースペクトルの平
均周波数とパワースペクトルの積分強度を示している。
血流速度は、2つに分けた場合、前腕、指ともに10倍程
度速度が異なる成分に分けられた。
密度に比例する。この値を見ると、前腕では早い速度の
赤血球の数は全体の14%程度であるが、指では早い速
度の成分が38%程度ある。この結果から、指のほうが
前腕に比べて太い血管の割合が高いと推測される。
生体組織に照射し、その反射散乱先に対しパワースペク
トル処理を施すことによって微小循環血流情報が得ら
れ、しかもこのパワースペクトル信号に多重回帰分析処
理を施すことによって血流速度分布を得ることができ
る。又、血流速度は血管の種類と太さに関係しているた
め、速度分布からその組織中の血管網構築が推定でき
る。組織中の毛細血管の割合は、病気などと関連があ
り、例えば糖尿病の患者では指先の毛細血管が病気の進
行とともに消失してくる。また、やけどの治り具合など
も推定できると考えられる。したがって、血流量の多少
だけでは判断できないことも、血流速度分布が分かれば
組織の状態を知ることができる。
Claims (1)
- 【請求項1】 生体に対しレーザー光を照射して、反射
散乱光を収集し、前記反射散乱光からパワースペクトル
を検出し、多重回帰処理を施すことを特徴とする生体信
号処理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP28029393A JP3478346B2 (ja) | 1993-10-14 | 1993-10-14 | 生体信号処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP28029393A JP3478346B2 (ja) | 1993-10-14 | 1993-10-14 | 生体信号処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH07113743A true JPH07113743A (ja) | 1995-05-02 |
JP3478346B2 JP3478346B2 (ja) | 2003-12-15 |
Family
ID=17622972
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP28029393A Expired - Fee Related JP3478346B2 (ja) | 1993-10-14 | 1993-10-14 | 生体信号処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3478346B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286870B2 (en) | 1994-10-06 | 2007-10-23 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
JP2012210321A (ja) * | 2011-03-31 | 2012-11-01 | Pioneer Electronic Corp | 光検出装置及び流体計測装置 |
WO2013183401A1 (ja) * | 2012-06-04 | 2013-12-12 | 株式会社アドバンテスト | 光音響診断装置、方法、プログラム、記録媒体 |
-
1993
- 1993-10-14 JP JP28029393A patent/JP3478346B2/ja not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286870B2 (en) | 1994-10-06 | 2007-10-23 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US7440794B2 (en) | 1994-10-06 | 2008-10-21 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US7715904B2 (en) | 1994-10-06 | 2010-05-11 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US8050744B2 (en) | 1994-10-06 | 2011-11-01 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
JP2012210321A (ja) * | 2011-03-31 | 2012-11-01 | Pioneer Electronic Corp | 光検出装置及び流体計測装置 |
WO2013183401A1 (ja) * | 2012-06-04 | 2013-12-12 | 株式会社アドバンテスト | 光音響診断装置、方法、プログラム、記録媒体 |
CN104168836A (zh) * | 2012-06-04 | 2014-11-26 | 株式会社爱德万测试 | 光声诊断装置、方法、程序、记录介质 |
Also Published As
Publication number | Publication date |
---|---|
JP3478346B2 (ja) | 2003-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nilsson et al. | A new instrument for continuous measurement of tissue blood flow by light beating spectroscopy | |
US4446239A (en) | Light scattering immunoassay involving particles with selective frequency band apparatus | |
Bonner et al. | Principles of laser-Doppler flowmetry | |
JP3577335B2 (ja) | 散乱吸収体計測方法及び装置 | |
JP3107914B2 (ja) | 散乱吸収体内部の吸収情報計測装置及び方法 | |
EP0102726B1 (en) | Method and apparatus for characterizing microparticles or measuring their response to their environment | |
CN100512752C (zh) | 光程可选择的无创人体成分检测方法及检测装置 | |
US5825488A (en) | Method and apparatus for determining analytical data concerning the inside of a scattering matrix | |
JP2505873B2 (ja) | 散乱光による粒子分析装置 | |
US5348002A (en) | Method and apparatus for material analysis | |
EP0703445B1 (en) | Method and apparatus for measuring concentration of absorptive constituent in scattering medium | |
EP0888531A1 (en) | Multi-parameter hematology apparatus and method | |
JPH11506202A (ja) | 非侵入的検査及びイメージングにおける散乱を最小化し且つ組織サンプリングを改善する方法 | |
Pedersen et al. | Transport calculations for light scattering in blood | |
US20020165456A1 (en) | Estimation of the average size of white light scatterers in normal and cancerous tissue using light scattering spectrum | |
Newson et al. | Laser Doppler velocimetry: the problem of fibre movement artefact | |
JP4018799B2 (ja) | 散乱吸収体の吸収成分の濃度計測方法及び装置 | |
EP1796535A1 (en) | A method of non-invasive measurement of sugar in blood and construction for its realisation | |
JPH07113743A (ja) | 生体信号処理方法 | |
US6660995B1 (en) | Particle size analysis in a turbid media with a single-fiber, optical probe while using a visible spectrometer | |
JPH1189799A (ja) | 特定成分の濃度測定装置および濃度測定方法 | |
CN109520998A (zh) | 一种用于医疗检测的光纤拉曼系统 | |
JP2641575B2 (ja) | グルコース無侵襲計測装置 | |
JP3694291B2 (ja) | 血糖値の無侵襲測定装置 | |
JP3114312B2 (ja) | 組織酸素流量計 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081003 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081003 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091003 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091003 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101003 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111003 Year of fee payment: 8 |
|
LAPS | Cancellation because of no payment of annual fees |