JP5658979B2 - 生体光計測装置及び生体における吸収係数の変動を推定する方法 - Google Patents
生体光計測装置及び生体における吸収係数の変動を推定する方法 Download PDFInfo
- Publication number
- JP5658979B2 JP5658979B2 JP2010253474A JP2010253474A JP5658979B2 JP 5658979 B2 JP5658979 B2 JP 5658979B2 JP 2010253474 A JP2010253474 A JP 2010253474A JP 2010253474 A JP2010253474 A JP 2010253474A JP 5658979 B2 JP5658979 B2 JP 5658979B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- living body
- distance
- absorption coefficient
- sensitivity
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 175
- 238000010521 absorption reaction Methods 0.000 title claims description 67
- 239000000523 sample Substances 0.000 claims description 99
- 238000005259 measurement Methods 0.000 claims description 77
- 230000008859 change Effects 0.000 claims description 71
- 230000035945 sensitivity Effects 0.000 claims description 70
- 238000001514 detection method Methods 0.000 claims description 53
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 13
- 230000003044 adaptive effect Effects 0.000 claims description 12
- 230000031700 light absorption Effects 0.000 claims description 10
- 239000004065 semiconductor Substances 0.000 claims description 8
- 239000013307 optical fiber Substances 0.000 claims description 7
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 230000017531 blood circulation Effects 0.000 description 56
- 102100023006 Basic leucine zipper transcriptional factor ATF-like 2 Human genes 0.000 description 53
- 101000903615 Homo sapiens Basic leucine zipper transcriptional factor ATF-like 2 Proteins 0.000 description 53
- 239000002344 surface layer Substances 0.000 description 39
- 230000003287 optical effect Effects 0.000 description 23
- 238000004088 simulation Methods 0.000 description 22
- 210000001519 tissue Anatomy 0.000 description 20
- 230000008326 skin blood flow Effects 0.000 description 18
- 239000006096 absorbing agent Substances 0.000 description 14
- 230000005540 biological transmission Effects 0.000 description 13
- 230000003925 brain function Effects 0.000 description 13
- 238000012876 topography Methods 0.000 description 13
- 238000009792 diffusion process Methods 0.000 description 10
- 210000004369 blood Anatomy 0.000 description 8
- 239000008280 blood Substances 0.000 description 8
- 210000004556 brain Anatomy 0.000 description 8
- 108010054147 Hemoglobins Proteins 0.000 description 7
- 102000001554 Hemoglobins Human genes 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 210000003710 cerebral cortex Anatomy 0.000 description 6
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 5
- 108010064719 Oxyhemoglobins Proteins 0.000 description 5
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 5
- 108010002255 deoxyhemoglobin Proteins 0.000 description 5
- 210000003625 skull Anatomy 0.000 description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 230000003727 cerebral blood flow Effects 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 4
- 210000003743 erythrocyte Anatomy 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 230000004913 activation Effects 0.000 description 3
- 230000007177 brain activity Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000006213 oxygenation reaction Methods 0.000 description 3
- 241000234295 Musa Species 0.000 description 2
- 235000018290 Musa x paradisiaca Nutrition 0.000 description 2
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 210000003403 autonomic nervous system Anatomy 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000004793 Polystyrene Substances 0.000 description 1
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000000427 antigen Substances 0.000 description 1
- 102000036639 antigens Human genes 0.000 description 1
- 108091007433 antigens Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- YPHMISFOHDHNIV-FSZOTQKASA-N cycloheximide Chemical compound C1[C@@H](C)C[C@H](C)C(=O)[C@@H]1[C@H](O)CC1CC(=O)NC(=O)C1 YPHMISFOHDHNIV-FSZOTQKASA-N 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006392 deoxygenation reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000008338 local blood flow Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229920002223 polystyrene Polymers 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
深部(深さ10mmの面上)に、独立した3個の吸収係数変化(ターゲット信号1〜3)を与え、これについて上述の3つの方法でシミュレーションを行いトポグラフィ画像を作成した。なお、以下の説明において、深部に与えた吸収係数変化をターゲット信号と呼び、表層部に与えた吸収係数変化を表層妨害信号と呼ぶ。本シミュレーションにおいては、以下のように第1〜第3ターゲット信号を定義する。第1〜第3ターゲット信号は、それぞれ、座標(x,y,z)=(15,−15,−10),(22.5,15,−10),(−15,0,−10)に位置し、それぞれの位置の吸収散乱媒質の体積1mm3の領域に、1%の吸収係数の変化があったとする。第1ターゲット信号は、y方向に沿った30mm間隔のオプトードの丁度中間点に位置しており、最も感度の高い部位に相当する。第2ターゲット信号は、2組の30mm間隔のオプトード間の感度の谷間になっている部位に位置している。また、第3ターゲット信号はx軸に沿った感度の高い部位に位置している。
次に、オプトード直下の座標(0,0,−2)の位置に表層妨害信号(以下、第1表装妨害信号と呼ぶ)がある場合について、シミュレーションを行った。この位置はオプトードの直下の浅い位置であるため、感度が極めて高いといえる。なお、表層妨害信号の強度(吸収係数の変化の大きさ)は上述の第1〜第3ターゲット信号の場合の10分の1である。図7(a)〜(c)は、それぞれ、SBP法、SRI法、SARI法により作成されたトポグラフィ画像を示す。なお、図7(a)〜(c)においては表示スケールを同一に設定している。
上述の第1ターゲット信号と第1表層妨害信号とを同時に与えた場合のシミュレーション結果を図8(a)〜(c)に示す。また、図9(a)〜(c)は、図8(a)〜(c)の表示スケールを同一に設定したものである。図8(a)及び図9(a)に示されるように、SBP法においてはターゲット信号と表層妨害信号とを区別できないことが分かった。これに対して、SRI法及びSARI法においては、若干の振動的な悪影響が現れているものの、表層妨害信号を効果的に抑制することができるとともに、画像コントラストがSBP法よりも優れていることが分かった。
表層妨害信号が空間の1点に分布している場合ではなく、ある面上に広がって分布している場合を考える。ここでは、上述の第1ターゲット信号に加えて、4点(−20,−30,−2)、(−20,20,−2)、(20,20,−2)、(20,−30,−2)で囲まれる領域に、1mm3あたりターゲット信号の100分の1の強度の表層妨害信号が分布しているとした。このような条件の元でのシミュレーション結果を図10(a)〜(c)に示す。また、図11(a)〜(c)は、図10(a)〜(c)の表示スケールを同一に設定したものである。
一般に、逆問題では雑音(ノイズ)の影響による解の不安定性が問題になることが知られている。そこで、ここでは前述の第1ターゲット信号に加えて、ショットノイズを想定した正規分布のランダムノイズを加えた場合についてのシミュレーションを行い、その結果を図12(a)〜12(c)に示す。また、図13(a)〜(c)は、図12(a)〜(c)の表示スケールを同一に設定したものである。ここで、ランダムノイズの標準偏差は、最大受光パワーの1%に相当するように設定されている。
前述の第1ターゲット信号及び第1表層妨害信号に加えて、さらに上述のランダムノイズを与えた場合についてシミュレーションを行い、その結果を図14(a)〜14(c)に示す。なお、また、図15(a)〜(c)は、図14(a)〜(c)の表示スケールを同一に設定したものである。図15(a)〜(c)を比較すると、SBP法においては、第1ターゲット信号と第1表層妨害信号とを区別することができず、SRI法及びSARI法においては第1表層信号の影響を抑制して第1ターゲット信号を再現できることが分かった。なお、SRI法とSARI法による結果を比べると、SARI法による結果の方が画像コントラストが高いことが分かる。
なお、上述のように、SRI法、SARI法においては、深部に対応する4636個のボクセルに推定値が返されるだけでなく、表層部に配置された18個のボクセルに表層妨害信号の推定値が返される。つまり、数14において、深部領域のターゲット信号sig1だけでなく、表層部の表層妨害信号sig2も推定することができるため、この表層妨害信号sig2をマッピングすることにより、表層妨害信号を現すトポグラフィ画像を作成することができる。上述の、「1個のターゲット信号と1個の表層妨害信号とランダムノイズ」を与えた場合における、SRI法、SARI法により得られた表層妨害信号のトポグラフィ画像をそれぞれ図16(a),(b)に示す。同様に、「1個のターゲット信号と広がった表層妨害信号」を与えた場合における、SRI法、SARI法により得られた表層妨害信号のトポグラフィ画像をそれぞれ図16(c),(d)に示し、「ターゲット信号とランダムノイズ」を与えた場合における、SRI法、SARI法により得られた表層妨害信号のトポグラフィ画像をそれぞれ図16(e),(f)に示す。
次に、ファントム(擬似体)を用いた予備実験について説明する。予備実験に使用した生体光測定装置は、前述の脳機能計測装置1と同様であるため、その構成についての説明はここでは省略する。図17に示されるように、ファントム40は、透明容器41と、透明容器41に満たされた分散液42と、分散液42中に沈められた、1辺の長さが4mmの立方体形状の黒色アクリル製吸収体Lと、同じく分散液42中に沈められた、1.4mm×1.4mm×2mmの大きさの、黒色アクリル製微小吸収体Sとを備える。ここで、分散液42は、染料(オリエント化学:WB−3)と散乱粒子(綜研化学:ポリスチレン球直径1.6μm)を、換算散乱係数μS’=1[mm−1]、吸収係数μa=0.02[mm−1]となるように調整されている。なお、吸収体L及び微小吸収体Sは、いずれも、直径0.5mmのステンレス製の棒材43により、背面側(図17におけるz軸方向のマイナス側)から支持されて、それぞれ独立に移動できるように構成されている。また、透明容器41の正面41a(図17においてz軸方向に垂直な面であって、z軸方向プラス側に位置する面)には、脳機能計測装置1の送光プローブ12及び受光プローブ13とが、前述のように配置されている(図2参照)。
また、図21(a),(b)は、同様の場合における、SRI法、SARI法により作成された表層領域のトポグラフィ画像である。図20(a)によれば、表層部に配置された微小吸収体Sの影響を受けて、(x,y)=(30,0)の位置にあるオプトードの周囲に大きなノイズが出現している。これに対して、図20(b)によれば、表層妨害信号の源となる微小吸収体Sの影響が取り除かれていることが分かる。さらに、(x,y)=(−10,−15)の近くに、吸収体Lに起因する吸収係数変化を示す像が現れているが、この像が2つに分離していることが分かった。これに対して、図20(c)によれば、同じく微小吸収体Sの影響が取り除かれているとともに、(x,y)=(−10,−15)の近くに、吸収体Lに起因する吸収係数変化を示す像が現れていることが分かった。なお、この場合においては、像は2つに分離せず単一の像として現れている。
11 半導体レーザ
12 送光プローブ
13 受光プローブ
14 アバランシェフォトダイオード
15 A/Dコンバータ
Claims (6)
- 生体内を通過した光を計測する生体光計測装置であって、
前記光を発するとともに、前記生体の表面の、複数の照射点に光学的に接続された光源と、
前記生体の表面の、前記照射点の1つからの距離が第1距離である第1の検出点及び第1距離よりも短い第2距離である第2の検出点を含む複数の検出点に、それぞれ光学的に接続されるとともに、前記光源から発せられた光のうち、前記生体の内部を通過してきた通過光の光量を検出する複数の検出部と、
前記検出部により検出された前記通過光の光量の変動に基づいて、前記生体内部の、第1の深さの領域に対応する複数の第1ボクセルにおける前記光の吸収係数の変動と、第2の深さの領域に対応する複数の第2ボクセルにおける前記光の吸収係数の変動とをそれぞれ推定するデータ処理部とを備え、
前記データ処理部は、Aを感度行列とし、λをティホノフの正則化パラメターとし、γを深部及び浅部の相対感度を調整するために設定される深さ感度パラメターとし、Iを単位行列としたときに、(ATA+λI(ATA)γ)−1ATの形で表される深さ感度可変感度適応型正則化を用いて逆問題を解くことにより前記第1、第2ボクセルにおける前記吸収係数の変動を推定する、生体光計測装置。 - さらに、前記光源と前記生体の表面の前記複数の照射点とをそれぞれ光学的に接続する、複数の光ファイバ製の送光プローブと、
前記検出部と前記生体の表面の前記検出点とをそれぞれ光学的に接続する、複数の光ファイバ製の受光プローブとを備え、
前記複数の送光プローブと前記複数の受光プローブとは、前記生体の表面において、互いに間隔の異なる2種類以上の送光プローブ及び受光プローブの組を形成するように、格子状に並べられている請求項1に記載の生体光計測装置。 - 前記複数の照射点及び前記複数の検出点は、第1距離だけ隔てて配置された複数の照射点及び検出点の組と、第2距離だけ隔てて配置された複数の照射点及び検出点の組とを含むように配置され、
第2距離だけ隔てて配置された複数の照射点及び検出点の組の数は、第1距離だけ隔てて配置された複数の照射点及び検出点の組の数よりも少ない、請求項1又は2に記載の生体光計測装置。 - 前記光源は、近赤外光を照射する半導体レーザ又はLEDである請求項1〜3のいずれか一項に記載の生体光計測装置。
- 前記検出部は、フォトダイオード又は光電子増倍管を有する請求項1〜4のいずれか一項に記載の生体光計測装置。
- 生体内を通過した光を計測して前記生体における吸収係数の変動を推定する方法であって、
光を発する光源を、前記生体の表面の、複数の照射点に光学的に接続することと、
複数の検出部を、前記生体の表面の、前記照射点の1つからの距離が第1距離である第1の検出点及び第1距離よりも短い第2距離である第2の検出点を含む複数の検出点に、それぞれ光学的に接続することと、
前記光源から前記光を照射し、前記生体の内部を通過してきた通過光の光量を前記検出部により検出することと、
前記生体内部の、第1の深さの領域に対応する複数の第1ボクセルと、第2の深さの領域に対応する複数の第2ボクセルとを設定することと、
前記検出部により検出された前記通過光の光量の変動に基づいて、前記第1、第2ボクセルにおける前記光の吸収係数の変動をそれぞれ推定することとを備え、
前記光の吸収係数の変動を推定する際に、Aを感度行列とし、λをティホノフの正則化パラメターとし、γを深部及び浅部の相対感度を調整するために設定される深さ感度パラメターとし、Iを単位行列としたときに、(ATA+λI(ATA)γ)−1ATの形で表される深さ感度可変感度適応型正則化を用いて逆問題を解くことにより前記2つの領域において前記光の前記吸収係数の変動をそれぞれ推定する、生体における吸収係数の変動を推定する方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010253474A JP5658979B2 (ja) | 2010-11-12 | 2010-11-12 | 生体光計測装置及び生体における吸収係数の変動を推定する方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010253474A JP5658979B2 (ja) | 2010-11-12 | 2010-11-12 | 生体光計測装置及び生体における吸収係数の変動を推定する方法 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2012100961A JP2012100961A (ja) | 2012-05-31 |
JP2012100961A5 JP2012100961A5 (ja) | 2013-12-19 |
JP5658979B2 true JP5658979B2 (ja) | 2015-01-28 |
Family
ID=46392117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010253474A Expired - Fee Related JP5658979B2 (ja) | 2010-11-12 | 2010-11-12 | 生体光計測装置及び生体における吸収係数の変動を推定する方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5658979B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6003762B2 (ja) * | 2013-03-27 | 2016-10-05 | セイコーエプソン株式会社 | 分光計測方法、分光計測器、および変換行列の生成方法 |
US20160025620A1 (en) * | 2014-07-22 | 2016-01-28 | Olympus Corporation | Optical Compositional Analysis of Mixtures |
CN111460966B (zh) * | 2020-03-27 | 2024-02-02 | 中国地质大学(武汉) | 基于度量学习和近邻增强的高光谱遥感图像分类方法 |
-
2010
- 2010-11-12 JP JP2010253474A patent/JP5658979B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2012100961A (ja) | 2012-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nguyen et al. | Bundled-optode implementation for 3D imaging in functional near-infrared spectroscopy | |
Strangman et al. | Depth sensitivity and source-detector separations for near infrared spectroscopy based on the Colin27 brain template | |
Patil et al. | Experimental investigation of NIRS spatial sensitivity | |
US9538926B2 (en) | Speckle contrast optical tomography | |
Zhou et al. | Diffuse optical correlation tomography of cerebral blood flow during cortical spreading depression in rat brain | |
US9498134B1 (en) | Diffuse optical tomography methods and system for determining optical properties | |
Zhan et al. | Image quality analysis of high-density diffuse optical tomography incorporating a subject-specific head model | |
JP5782314B2 (ja) | 生体計測装置および画像作成方法 | |
Kamran et al. | Differential path-length factor's effect on the characterization of brain's hemodynamic response function: a functional near-infrared study | |
Kannan et al. | A computational model to detect and quantify a primary blast lung injury using near‐infrared optical tomography | |
JP5658979B2 (ja) | 生体光計測装置及び生体における吸収係数の変動を推定する方法 | |
Heiskala et al. | Significance of background optical properties, time-resolved information and optode arrangement in diffuse optical imaging of term neonates | |
Causin et al. | Elastic net regularization in diffuse optical tomography applications | |
JP6472145B2 (ja) | 脳機能計測装置およびそれを用いた光照射・光検出方法 | |
Shoaib et al. | Approach to optimize 3-dimensional brain functional activation image with high resolution: a study on functional near-infrared spectroscopy | |
CN115516849A (zh) | 来自光相干的光学成像 | |
Brigadoi et al. | Evaluating real-time image reconstruction in diffuse optical tomography using physiologically realistic test data | |
JP4543774B2 (ja) | 生体光計測装置 | |
JP2005237618A (ja) | 生体光計測装置および信号処理方法 | |
JP2005160783A (ja) | 非侵襲脳活動計測方法 | |
Liu et al. | Suppressing physiological interferences and physical noises in functional diffuse optical tomography via tandem inversion filtering and LSTM classification | |
Meng et al. | Noninvasive optical monitoring of pulmonary embolism: a Monte Carlo study on visible Chinese human thoracic tissues | |
Chance et al. | Photon dynamics in tissue imaging | |
Lin et al. | Anatomical Modeling and Optimization of Speckle Contrast Optical Tomography | |
Zhan | Model-based high-density functional diffuse optical tomography of human brain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20131030 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131101 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140331 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140415 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140609 |
|
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: 20141111 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141201 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5658979 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |