JP6011636B2 - 光生体計測装置 - Google Patents
光生体計測装置 Download PDFInfo
- Publication number
- JP6011636B2 JP6011636B2 JP2014546785A JP2014546785A JP6011636B2 JP 6011636 B2 JP6011636 B2 JP 6011636B2 JP 2014546785 A JP2014546785 A JP 2014546785A JP 2014546785 A JP2014546785 A JP 2014546785A JP 6011636 B2 JP6011636 B2 JP 6011636B2
- Authority
- JP
- Japan
- Prior art keywords
- measurement
- light
- correlation coefficient
- displayed
- brain
- 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.)
- Active
Links
- 238000005259 measurement Methods 0.000 title claims description 153
- 230000003287 optical effect Effects 0.000 title claims description 58
- 210000004556 brain Anatomy 0.000 claims description 77
- 239000000523 sample Substances 0.000 claims description 63
- 210000004761 scalp Anatomy 0.000 claims description 61
- 230000005540 biological transmission Effects 0.000 claims description 19
- 108010064719 Oxyhemoglobins Proteins 0.000 description 13
- 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 11
- 108010002255 deoxyhemoglobin Proteins 0.000 description 11
- 108010054147 Hemoglobins Proteins 0.000 description 6
- 102000001554 Hemoglobins Human genes 0.000 description 6
- 238000000034 method Methods 0.000 description 6
- 230000007177 brain activity Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 210000003484 anatomy Anatomy 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000002835 absorbance Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000258 photobiological effect Effects 0.000 description 2
- 240000005561 Musa balbisiana Species 0.000 description 1
- 235000018290 Musa x paradisiaca Nutrition 0.000 description 1
- 230000003925 brain function Effects 0.000 description 1
- 230000003727 cerebral blood flow Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0233—Special features of optical sensors or probes classified in A61B5/00
- A61B2562/0238—Optical sensor arrangements for performing transmission measurements on body tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Neurology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、光生体計測装置に関し、さらに詳細には非侵襲で脳活動を測定する光生体計測装置に関する。
近年、脳の活動状況を観察するために、光を用いて簡便に非侵襲で測定する光脳機能イメージング装置が開発されている。このような光脳機能イメージング装置では、被検者の頭皮表面上に配置した送光プローブにより、異なる3種類の波長λ1、λ2、λ3(例えば、780nmと805nmと830nm)の近赤外光を脳に照射するとともに、頭皮表面上に配置した受光プローブにより、脳から放出された各波長λ1、λ2、λ3の近赤外光の強度変化(受光量情報)ΔA(λ1)、ΔA(λ2)、ΔA(λ3)をそれぞれ検出する。
そして、このようにして得られた受光量情報ΔA(λ1)、ΔA(λ2)、ΔA(λ3)から、脳血流中のオキシヘモグロビンの濃度変化・光路長積[oxyHb]と、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]とを求めるために、例えば、ModifiedBeerLambert則を用いて関係式(1)(2)(3)に示す連立方程式を作成して、この連立方程式を解いている。さらには、オキシヘモグロビンの濃度変化・光路長積[oxyHb]と、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]とから総ヘモグロビンの濃度変化・光路長積([oxyHb]+[deoxyHb])を算出している。
ΔA(λ1)=EO(λ1)×[oxyHb]+Ed(λ1)×[deoxyHb]・・・(1)
ΔA(λ2)=EO(λ2)×[oxyHb]+Ed(λ2)×[deoxyHb]・・・(2)
ΔA(λ3)=EO(λ3)×[oxyHb]+Ed(λ3)×[deoxyHb]・・・(3)
なお、EO(λm)は、波長λmの光におけるオキシヘモグロビンの吸光度係数であり、Ed(λm)は、波長λmの光におけるデオキシヘモグロビンの吸光度係数である。
そして、このようにして得られた受光量情報ΔA(λ1)、ΔA(λ2)、ΔA(λ3)から、脳血流中のオキシヘモグロビンの濃度変化・光路長積[oxyHb]と、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]とを求めるために、例えば、ModifiedBeerLambert則を用いて関係式(1)(2)(3)に示す連立方程式を作成して、この連立方程式を解いている。さらには、オキシヘモグロビンの濃度変化・光路長積[oxyHb]と、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]とから総ヘモグロビンの濃度変化・光路長積([oxyHb]+[deoxyHb])を算出している。
ΔA(λ1)=EO(λ1)×[oxyHb]+Ed(λ1)×[deoxyHb]・・・(1)
ΔA(λ2)=EO(λ2)×[oxyHb]+Ed(λ2)×[deoxyHb]・・・(2)
ΔA(λ3)=EO(λ3)×[oxyHb]+Ed(λ3)×[deoxyHb]・・・(3)
なお、EO(λm)は、波長λmの光におけるオキシヘモグロビンの吸光度係数であり、Ed(λm)は、波長λmの光におけるデオキシヘモグロビンの吸光度係数である。
ここで、送光プローブと受光プローブとの間の距離と、測定部位との関係について説明する。図4は、一対の送光プローブ及び受光プローブと、測定部位との関係を示す図である。送光プローブ12が被検者の頭皮表面の送光点Tに押し当てられるとともに、受光プローブ13が被検者の頭皮表面の受光点Rに押し当てられる。そして、送光プローブ12から光を照射させるとともに、受光プローブ13に頭皮表面から放出される光を入射させる。このとき、頭皮表面の送光点Tから照射された光のうちで、バナナ形状(測定領域)を通過した光が頭皮表面の受光点Rに到達する。これにより、測定領域の中でも、特に送光点Tと受光点Rとを被検者の頭皮表面に沿って最短距離で結んだ線の中点(頭皮表面上の測定関連点)Mから、送光点Tと受光点Rとを被検者の頭皮表面に沿って最短距離で結んだ線の距離の半分の深さである被検者の測定部位(脳表面上の測定関連点)Sに関する受光量情報A(λ1)、A(λ2)、A(λ3)が得られるとしている。
また、光脳機能イメージング装置では、脳の複数箇所の測定部位に関するオキシヘモグロビンの濃度変化・光路長積[oxyHb]、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]及び総ヘモグロビンの濃度変化・光路長積([oxyHb]+[deoxyHb])をそれぞれ測定することが行われている。
このような光脳機能イメージング装置においては、15個の送光プローブと15個の受光プローブとを所定の配列で被検者の頭皮表面に接触させるために、ホルダが使用される。図5は、15個の送光プローブと15個の受光プローブとが挿入されるホルダの一例を示す平面図である。
送光プローブ12T1〜12T15と受光プローブ13R1〜13R15とは、縦方向に5個と横方向に6個とに交互となるように配置されている。このとき、送光プローブ12T1〜12T15と受光プローブ13R1〜13R15との間の間隔は、30mmとなっている。これにより、脳の49箇所の測定部位に関する受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)(n=1、2、・・・、49)を得ている。
このような光脳機能イメージング装置においては、15個の送光プローブと15個の受光プローブとを所定の配列で被検者の頭皮表面に接触させるために、ホルダが使用される。図5は、15個の送光プローブと15個の受光プローブとが挿入されるホルダの一例を示す平面図である。
送光プローブ12T1〜12T15と受光プローブ13R1〜13R15とは、縦方向に5個と横方向に6個とに交互となるように配置されている。このとき、送光プローブ12T1〜12T15と受光プローブ13R1〜13R15との間の間隔は、30mmとなっている。これにより、脳の49箇所の測定部位に関する受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)(n=1、2、・・・、49)を得ている。
そして、49個の受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)を所定時間間隔Δtで得ていくことで、関係式(1)(2)(3)を用いて、オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Xn(t)、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]の経時変化(測定データ)Yn(t)及び総ヘモグロビンの濃度変化・光路長積([oxyHb]+[deoxyHb])の経時変化(測定データ)Zn(t)(n=1、2、・・・、49)を求めている。
このとき、医師や検査技師等が、被検者へ刺激(以下、「負荷」や「タスク」という)を与えて、被検者の脳を活動させたときの測定データXn(t)、Yn(t)、Zn(t)を得ることも行われている(例えば、特許文献1参照)。
そして、医師や検査技師等は、被検者へタスクを与えたときの脳の活動領域を調べて、どの領域が何の機能を担っているのかを解明しようとする「脳の機能局在研究」も行われている。
そして、医師や検査技師等は、被検者へタスクを与えたときの脳の活動領域を調べて、どの領域が何の機能を担っているのかを解明しようとする「脳の機能局在研究」も行われている。
さらに、脳の各領域の機能だけでなく、複数の領域間の関係を調べて、脳内ネットワークの機能を解明しようとする「脳領域間の接続性研究」も行われている。
例えば、下記式(4)を用いて、第一の測定部位S1で得られた測定データX1(t)と、第二の測定部位S2で得られた測定データX2(t)との相関の強さを示す相関係数値α1−2を求めている。
例えば、下記式(4)を用いて、第一の測定部位S1で得られた測定データX1(t)と、第二の測定部位S2で得られた測定データX2(t)との相関の強さを示す相関係数値α1−2を求めている。
ここで、図6は、第一の測定部位S1で得られた測定データX1(t)と、他の測定部位Snで得られた測定データXn(t)(n=2、・・・、49)との間の相関係数値α1−nを示した表示画面の一例である。図6に示す表示画面では、図5に示した平面図において、送光プローブ12T1〜12T15と受光プローブ13R1〜13R15とを最短距離で結んだ線の各中点が、その送光プローブ12から照射させた光を、その受光プローブ13で検出させたときに得られた測定点Snとなるようにしている。そして、表示画面に測定点Snを表示するとともに、測定データX1(t)が得られた測定点S1と、測定データXn(t)(n=2、・・・、49)が得られた各測定点Snとを、相関係数値α1−nに応じた太さの直線L1−nで結んでいる。
しかしながら、図6に示す表示画面では脳を示す表示がないため、脳のどの領域とどの領域との間の相関係数値α1−nが大きいかを把握することが困難であり、さらに、一の直線L1−nと他の直線L1−n’とのオーバーラップにより、医師や検査技師等の観察において誤解を生じることがあった。また、脳の解剖学的構造には個人差があり、脳の形状が各人で違っていることから、脳の領域を正確に判別することができないという問題も生じていた。
そこで、本発明は、脳の領域間の相関関係値を正確に把握することができる光生体計測装置を提供することを目的とする。
そこで、本発明は、脳の領域間の相関関係値を正確に把握することができる光生体計測装置を提供することを目的とする。
上記課題を解決するためになされた本発明の光生体計測装置は、被検者の頭皮表面上に配置される複数個の送光プローブと、当該頭皮表面上に配置される複数個の受光プローブと、前記送光プローブが頭皮表面に光を照射するとともに、前記受光プローブが頭皮表面から放出される光を検出するように制御することで、複数個の測定部位に関する複数個の受光量情報を取得する送受光用制御部と、複数個の受光量情報に基づいて、複数個の測定データを取得する演算部と、各測定データ間の相関の強さを示す相関係数値を算出する相関係数算出部とを備える光生体計測装置であって、前記被検者を撮影することで取得された3次元脳表面画像を表示する3次元画像表示制御部と、前記3次元脳表面画像上に、前記測定データが得られる測定関連点を表示するとともに、前記相関係数値の大きさを示す線分で測定関連点と測定関連点とを結ぶように表示する相関係数表示制御部とを備え、前記3次元脳表面画像は、所望の方向から見た画像となるように、当該方向を変更可能として表示されることを特徴としている。
ここで、「測定データ」とは、受光プローブで検出された受光量情報の経時変化自体であってもよく、受光量情報から算出されたオキシヘモグロビン濃度の経時変化やデオキシヘモグロビン濃度の経時変化や全ヘモグロビン濃度の経時変化であってもよい。
また、「3次元脳表面画像」とは、例えば、MRIやCT画像等により作成された被検者の映像データから、脳表面を示す映像データを抽出することにより作成された3次元画像等のことをいう。
また、「3次元脳表面画像」とは、例えば、MRIやCT画像等により作成された被検者の映像データから、脳表面を示す映像データを抽出することにより作成された3次元画像等のことをいう。
本発明の光生体計測装置によれば、3次元脳表面画像上に、相関係数値の大きさを示す線分が重畳されて表示されるので、医師や検査技師等は、脳の領域間の相関係数値を正確に把握することができる。
(その他の課題を解決するための手段及び効果)
また、本発明の光生体計測装置においては、前記線分は、前記相関係数値の大きさに応じて太さが変化する直線であるようにしてもよい。
また、本発明の光生体計測装置においては、前記線分は、前記相関係数値の大きさに応じて太さが変化する直線であるようにしてもよい。
また、本発明の光生体計測装置においては、前記測定関連点は、脳表面上では、前記送光プローブの位置と前記受光プローブの位置とを結んだ線の垂直二等分線と交差する位置とする脳表面上の測定関連点とするようにしてもよい。
そして、本発明の光生体計測装置においては、前記3次元脳表面画像は、所望の方向から見た画像となるように、当該方向を変更可能として表示されるようにしている。
本発明の光生体計測装置によれば、3次元脳表面画像を所望の方向から見た画像となるように表示することができるので、特に注目したい領域が中心にくるようにして表示することができる。
本発明の光生体計測装置によれば、3次元脳表面画像を所望の方向から見た画像となるように表示することができるので、特に注目したい領域が中心にくるようにして表示することができる。
さらに、本発明の光生体計測装置においては、前記3次元脳表面画像は、前記被検者を撮影することで取得されたものであるようにしている。
本発明の光生体計測装置によれば、被検者の3次元脳表面画像を表示するので、脳の解剖学的構造の個人差にかかわらず、脳活動を詳細に考察することができる。
本発明の光生体計測装置によれば、被検者の3次元脳表面画像を表示するので、脳の解剖学的構造の個人差にかかわらず、脳活動を詳細に考察することができる。
以下、本発明の実施形態について図面を用いて説明する。なお、本発明は、以下に説明するような実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の態様が含まれる。
図1は、本発明の一実施形態である光生体計測装置の構成を示すブロック図である。光生体計測装置1は、光を出射する光源2と、光源2を駆動する光源駆動機構4と、光を検出する光検出器3と、A/D(A/Dコンバータ)5と、送受光用制御部21と、演算部22と、相関係数算出部23と、相関係数表示制御部24と、3次元画像表示制御部32と、ポインタ表示制御部33と、測定関連点算出部34と、メモリ25とを備えるとともに、15個の送光プローブ12T1〜12T15と、15個の受光プローブ13R1〜13R15と、表示装置26と、入力装置27とを備える。
ここで、本発明に係る光生体計測装置1により表示される表示画面について説明する。
図2は、本発明に係る光生体計測装置1により表示される表示画面の一例を示す図である。図2に示す表示画面には、3次元頭皮表面画像41と、3次元脳表面画像42と、49個の測定関連点Mn(n=1,2,・・・,49)と、直線L1−nとが表示されている。49個の測定関連点Mnは、3次元頭皮表面画像41上の所定位置(頭皮表面上の測定関連点に対応する位置)にそれぞれ表示されている。そして、測定関連点M1と測定関連点M2とを結ぶ直線L1−2が太さ5mmで表示され、測定関連点M1と測定関連点M3とを結ぶ直線L1−3が太さ3mmで表示されるように、測定関連点M1と他の各測定関連点Mn(n=2,・・・,49)とを結ぶ直線L1−nがそれぞれの太さで表示されている。なお、頭皮表面画像41は、半透明で表示されている。
図2は、本発明に係る光生体計測装置1により表示される表示画面の一例を示す図である。図2に示す表示画面には、3次元頭皮表面画像41と、3次元脳表面画像42と、49個の測定関連点Mn(n=1,2,・・・,49)と、直線L1−nとが表示されている。49個の測定関連点Mnは、3次元頭皮表面画像41上の所定位置(頭皮表面上の測定関連点に対応する位置)にそれぞれ表示されている。そして、測定関連点M1と測定関連点M2とを結ぶ直線L1−2が太さ5mmで表示され、測定関連点M1と測定関連点M3とを結ぶ直線L1−3が太さ3mmで表示されるように、測定関連点M1と他の各測定関連点Mn(n=2,・・・,49)とを結ぶ直線L1−nがそれぞれの太さで表示されている。なお、頭皮表面画像41は、半透明で表示されている。
また、図3は、本発明に係る光生体計測装置1により表示される表示画面の他の一例を示す図である。図3に示す表示画面には、3次元脳表面画像42と、49個の測定関連点Sn(n=1,2,・・・,49)と、直線L1−nとが表示されている。49個の測定関連点Snは、3次元脳表面画像42上の所定位置(脳表面上の測定関連点に対応する位置)にそれぞれ表示されている。そして、測定関連点S1と測定関連点S2とを結ぶ直線L1−2が太さ5mmで表示され、測定関連点S1と測定関連点S3とを結ぶ直線L1−3が太さ3mmで表示されるように、測定関連点S1と他の各測定関連点Sn(n=1,2,4,・・・,49)とを結ぶ直線L1−nがそれぞれの太さで表示されている。
次に、図2や図3の表示画面を表示するための光生体計測装置1の構成について説明する。
光源駆動機構4は、送受光用制御部21から入力された駆動信号により光源2を駆動する。光源2は、例えば異なる3種類の波長λ1、λ2、λ3の近赤外光を出射することができる半導体レーザLD1、LD2、LD3等である。
光検出器3は、例えば光電子増倍管等であり、15個の受光プローブ13R1〜13R15で受光した近赤外光を個別に検出することにより、15個の受光量情報ΔA(λ1)、ΔA(λ2)、ΔA(λ3)をA/D5を介して送受光用制御部21に出力する。
光源駆動機構4は、送受光用制御部21から入力された駆動信号により光源2を駆動する。光源2は、例えば異なる3種類の波長λ1、λ2、λ3の近赤外光を出射することができる半導体レーザLD1、LD2、LD3等である。
光検出器3は、例えば光電子増倍管等であり、15個の受光プローブ13R1〜13R15で受光した近赤外光を個別に検出することにより、15個の受光量情報ΔA(λ1)、ΔA(λ2)、ΔA(λ3)をA/D5を介して送受光用制御部21に出力する。
送受光用制御部21は、所定の時間に1個の送光プローブ12に光を送光する駆動信号を光源駆動機構4に出力するとともに、受光プローブ13で受光された受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)(n=1,2,・・・,49)を光検出器3で検出する制御を行う。具体的には、まず5ミリ秒間、送光プローブ12T1に波長780nmの光を送光させ、次の5ミリ秒間、送光プローブ12T1に波長805nmの光を送光させ、次の5ミリ秒間、送光プローブ12T1に波長830nmの光を送光させ、次の5ミリ秒間、送光プローブ12T2に波長780nmの光を送光させるように、所定のタイミングで1個の送光プローブ12T1〜12T15に光を順番に送光させていく。このとき、いずれか1個の送光プローブ12T1〜12T15に光を送光させるごとに、15個の受光プローブ13R1〜13R15で受光量情報を検出することになるが、所定のタイミングで検出した所定(光を照射した送光プローブと隣接する)の受光プローブ13R1〜13R15の受光量情報をメモリ25に記憶させる。これにより、合計49個の受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)(n=1,2,・・・,49)の収集が行われる。
演算部22は、メモリ25に記憶された49個の受光量情報ΔAn(λ1)、ΔAn(λ2)、ΔAn(λ3)に基づいて、関係式(1)(2)(3)を用いて、オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Xn(t)、デオキシヘモグロビンの濃度変化・光路長積[deoxyHb]の経時変化(測定データ)Yn(t)及び総ヘモグロビンの濃度変化・光路長積([oxyHb]+[deoxyHb])の経時変化(測定データ)Zn(t)(n=1,2,・・・,49)を求める制御を行う。
相関係数算出部23は、各測定データXn(t)、Yn(t)、Zn(t)間の相関の強さを示す相関係数値αn−n’を算出する制御を行う。例えば、前記式(4)を用いて、第一の測定部位S1で得られたオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)X1(t)と、第二の測定部位S2で得られたオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)X2(t)との間の相関係数値α1−2を算出し、第一の測定部位S1で得られたオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)X1(t)と、第三の測定部位S3で得られたオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)X3(t)との間の相関係数値α1−3を算出するように、1176個(49C48個)の相関係数値αn−n’を算出する。また、式(4)を用いて、第一の測定部位S1で得られたデオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Y1(t)と、第二の測定部位S2で得られたデオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Y2(t)との間の相関係数値β1−2を算出し、第一の測定部位S1で得られたデオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Y1(t)と、第三の測定部位S3で得られたデオキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Y3(t)との間の相関係数値β1−3を算出するように、1176個(49C48個)の相関係数値βn−n’を算出する。さらに、式(4)を用いて、第一の測定部位S1で得られた総オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Z1(t)と、第二の測定部位S2で得られた総オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Z2(t)との間の相関係数値γ1−2を算出し、第一の測定部位S1で得られた総オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Z1(t)と、第三の測定部位S3で得られた総オキシヘモグロビンの濃度変化・光路長積[oxyHb]の経時変化(測定データ)Z3(t)との間の相関係数値γ1−3を算出するように、1176個(49C48個)の相関係数値γn−n’を算出する。
3次元画像表示制御部32は、3次元画像取得部32dと、頭皮表面画像表示制御部32aと、脳表面画像表示制御部32bと、画像切替部32cとを有する。
3次元画像取得部32dは、計測前にMRI100により作成された映像データを取得することにより、頭皮表面を示す映像データを抽出することで3次元頭皮表面画像データを取得するとともに、脳表面を示す映像データを抽出することで3次元脳表面画像データを取得して、3次元頭皮表面画像データと3次元脳表面画像データとをメモリ25に記憶させる制御を行う。ここで、MRI100は、3方向の2次元画像を示す映像データを作成するものである。なお、映像データは、頭皮表面と脳表面とを含む被険者を示すものである。また、映像データは、MRI信号の強度情報や位相情報等の数値を有する複数のピクセルから構成される。そして、上述した抽出する方法としては、例えば、MRI信号の強度情報や位相情報等の数値を有する複数のピクセルを用いることにより、領域拡張法、領域併合法、ヒューリスティック法等の画像領域分割方法、境界要素を連結して領域を抽出する方法、閉曲線を変形させて領域を抽出する方法等を利用することが挙げられる。
3次元画像取得部32dは、計測前にMRI100により作成された映像データを取得することにより、頭皮表面を示す映像データを抽出することで3次元頭皮表面画像データを取得するとともに、脳表面を示す映像データを抽出することで3次元脳表面画像データを取得して、3次元頭皮表面画像データと3次元脳表面画像データとをメモリ25に記憶させる制御を行う。ここで、MRI100は、3方向の2次元画像を示す映像データを作成するものである。なお、映像データは、頭皮表面と脳表面とを含む被険者を示すものである。また、映像データは、MRI信号の強度情報や位相情報等の数値を有する複数のピクセルから構成される。そして、上述した抽出する方法としては、例えば、MRI信号の強度情報や位相情報等の数値を有する複数のピクセルを用いることにより、領域拡張法、領域併合法、ヒューリスティック法等の画像領域分割方法、境界要素を連結して領域を抽出する方法、閉曲線を変形させて領域を抽出する方法等を利用することが挙げられる。
頭皮表面画像表示制御部32aは、メモリ25に記憶された3次元頭皮表面画像データに基づいて、図2の3次元頭皮表面画像41を表示装置26に表示する制御を行う。なお、操作者が入力装置27を用いて、所望の方向から見た3次元頭皮表面画像41となるように、方向を変更して表示することができるようになっている。また、3次元頭皮表面画像41は、半透明で表示されたり、着色して表示されたりすることができるようになっている。
脳表面画像表示制御部32bは、メモリ25に記憶された3次元脳表面画像データに基づいて、3次元脳表面画像42を表示装置26に表示する制御を行う。なお、操作者が入力装置27を用いて、所望の方向から見た3次元脳表面画像42となるように、方向を変更して表示することができるようになっている。
脳表面画像表示制御部32bは、メモリ25に記憶された3次元脳表面画像データに基づいて、3次元脳表面画像42を表示装置26に表示する制御を行う。なお、操作者が入力装置27を用いて、所望の方向から見た3次元脳表面画像42となるように、方向を変更して表示することができるようになっている。
画像切替部32cは、操作者が入力装置27を用いて操作信号を入力することにより、3次元頭皮表面画像41を頭皮表面画像表示制御部32aに表示させるように決定したり、3次元脳表面画像42を脳表面画像表示制御部32bに表示させるように決定したり、3次元頭皮表面画像41を頭皮表面画像表示制御部32aに表示させると同時に3次元脳表面画像42を脳表面画像表示制御部32bに表示させるように決定したりする制御を行う。なお、3次元頭皮表面画像41を表示させると同時に3次元脳表面画像42を表示させる場合には、3次元頭皮表面画像41と3次元脳表面画像42とは位置を合わせた状態で重ねて表示されるようにしている。
ポインタ表示制御部33は、表示装置26にポインタ(図示せず)を表示するとともに、入力装置27から出力された操作信号に基づいて、表示装置26に表示されたポインタを移動したり、ポインタで画像上の位置を指定したりする制御を行う。
測定関連点算出部34は、表示装置26に表示された3次元頭皮表面画像41や3次元脳表面画像42の所定の位置がポインタで指定されることにより、メモリ25に記憶された3次元頭皮表面画像データと3次元脳表面画像データとを用いて、被検者の脳表面上や頭皮表面上との測定関連点Sn、Mnを決定する制御を行う。例えば、3次元頭皮表面画像41の所定の第一の位置がポインタで指定されると、指定された第一の位置が脳表面上の第一の測定関連点S1となり、第一の測定関連点S1から最短距離にある頭皮表面上の位置が頭皮表面上の第一の測定関連点M1となる。そして、3次元頭皮表面画像41上の第一の測定関連点M1が中点となる送光プローブ12の位置と受光プローブ13の位置とを算出して表示する。また、3次元頭皮表面画像41の所定の第二の位置がポインタで指定されると、指定された第二の位置が脳表面上の第二の測定関連点S2となり、第二の測定関連点S2から最短距離にある頭皮表面上の位置が頭皮表面上の第一の測定関連点M2となる。そして、3次元頭皮表面画像41上の第二の測定関連点M2が中点となる送光プローブ12の位置と受光プローブ13の位置とを算出して表示する。
このように、頭皮表面上のN個の測定関連点Mn(n=1,2,・・・,49)と、脳表面上のN個の測定関連点Sn(n=1,2,・・・,49)とが決定され、送光プローブ12の位置と受光プローブ13の位置とが表示される。これにより、医師や検査技師等は、被検者の頭皮表面上における送光プローブ12の配置位置と受光プローブ13の配置位置とを把握することができる。
このように、頭皮表面上のN個の測定関連点Mn(n=1,2,・・・,49)と、脳表面上のN個の測定関連点Sn(n=1,2,・・・,49)とが決定され、送光プローブ12の位置と受光プローブ13の位置とが表示される。これにより、医師や検査技師等は、被検者の頭皮表面上における送光プローブ12の配置位置と受光プローブ13の配置位置とを把握することができる。
相関係数表示制御部24は、3次元頭皮表面画像41や3次元脳表面画像42上に、測定関連点算出部34で決定された測定関連点Sn、Mnを表示するとともに、相関係数値αn−n’の大きさを示す直線Ln−n’で測定関連点Sn、Mnと測定関連点Sn、Mnとを結ぶように表示する制御を行う。このとき、相関係数値αn−n’の大きさと直線Ln−n’の太さとの対応関係を示すテーブルに基づいて、相関係数値αn−n’の大きさを直線Ln−n’の太さで表現する。例えば、相関係数値αn−n’が0.7以上であるときには、直線Ln−n’の太さは5mmとなり、相関係数値αn−n’が0.4以上0.7未満であるときには、直線Ln−n’の太さは3mmとなり、相関係数値αn−n’が0.1以上0.4未満であるときには、直線Ln−n’の太さは1mmとなり、相関係数αn−n’が0.1未満であるときには、直線Ln−n’の太さは0mm(直線は非表示)となるようにする。
そして、例えば、画像切替部32cで3次元頭皮表面画像41と3次元脳表面画像42とを表示させ、測定データX1(t)と他の各測定データXn(t)(n=2,・・・,49)との間の相関係数値α1−nを表示させるように操作者が入力装置27を用いて指示した場合には、相関係数表示制御部24は、3次元頭皮表面画像41上に49個の測定関連点Mnを表示して、測定データX1(t)と測定データX2(t)との間の相関係数値α1−2が0.7であるので、測定関連点M1と測定関連点M2とを結ぶ直線L1−2を太さ5mmで表示させ、測定データX1(t)と測定データX3(t)との間の相関係数値α1−3が0.3であるので、測定関連点M1と測定関連点M3とを結ぶ直線L1−3を太さ1mmで表示させるように、測定関連点#M1と他の各測定関連点Mn(n=2,・・・,49)とを結ぶ直線L1−nをそれぞれの太さで表示させる。(図2参照)。
これにより、3次元頭皮表面画像41と3次元脳表面画像42上に、相関係数値α1−nの大きさを示す直線L1−nが重畳して表示されるので、医師や検査技師等は、脳の領域間の相関係数値α1−nを正確に把握することができる。
これにより、3次元頭皮表面画像41と3次元脳表面画像42上に、相関係数値α1−nの大きさを示す直線L1−nが重畳して表示されるので、医師や検査技師等は、脳の領域間の相関係数値α1−nを正確に把握することができる。
また、画像切替部32cで3次元脳表面画像42を表示させ、測定データY3(t)と他の各測定データYn(t)(n=1,2,4・・・,49)との間の相関係数値β3−nを表示させるように操作者が入力装置27を用いて指示した場合には、相関係数表示制御部24は、3次元脳表面画像42上に49個の測定関連点Snを表示して、測定関連点S1と測定関連点S2とを結ぶ直線L1−2を太さ5mmで表示させ、測定関連点S1と測定関連点S3とを結ぶ直線L1−3を太さ3mmで表示させるように、測定関連点S1と他の各測定関連点Mn(n=1,2,4,・・・,49)とを結ぶ直線L1−nをそれぞれの太さで表示させる(図3参照)。
これにより、3次元脳表面画像42上に、相関係数値β1−nの大きさを示す直線L1−nが重畳されて表示されるので、医師や検査技師等は、脳の領域間の相関係数値β1−nを正確に把握することができる。
さらに、光生体計測装置1では、3次元頭皮表面画像41や3次元脳表面画像42を所望の方向から見た画像となるように、方向を変更可能として表示することができるので、相関係数表示制御部24は、49個の測定関連点Sn、Mnと直線Ln−n’とを表示する位置も移動させることになる。
これにより、3次元脳表面画像42上に、相関係数値β1−nの大きさを示す直線L1−nが重畳されて表示されるので、医師や検査技師等は、脳の領域間の相関係数値β1−nを正確に把握することができる。
さらに、光生体計測装置1では、3次元頭皮表面画像41や3次元脳表面画像42を所望の方向から見た画像となるように、方向を変更可能として表示することができるので、相関係数表示制御部24は、49個の測定関連点Sn、Mnと直線Ln−n’とを表示する位置も移動させることになる。
以上のように、本発明の光生体計測装置1によれば、3次元頭皮表面画像41や3次元脳表面画像42上に、相関係数値αn−n’、βn−n’、γn−n’の大きさを示す直線Ln−n’が重畳して表示されるので、医師や検査技師等は、脳の領域間の相関係数値αn−n’、βn−n’、γn−n’を正確に把握することができる。また、3次元頭皮表面画像41や3次元脳表面画像42を所望の方向から見た画像となるように表示することができるので、特に注目したい領域が表示画面の中心にくるように表示することができる。さらに、被検者自身の3次元頭皮表面画像41や3次元脳表面画像42を表示することができるので、脳の解剖学的構造の個人差にかかわらず、脳活動を詳細に考察することができる。
<他の実施形態>
(1)上述した光生体計測装置1では、測定データX1(t)と他の各測定データXn(t)(n=2,・・・,49)との間の相関係数値α1−nを表示させる構成を示したが、全ての測定データXn(t)間の相関係数値αn−n’を表示させる構成としてもよい。つまり、1176個(49C48個)の直線Ln−n’を表示してもよい。
(1)上述した光生体計測装置1では、測定データX1(t)と他の各測定データXn(t)(n=2,・・・,49)との間の相関係数値α1−nを表示させる構成を示したが、全ての測定データXn(t)間の相関係数値αn−n’を表示させる構成としてもよい。つまり、1176個(49C48個)の直線Ln−n’を表示してもよい。
(2)上述した光生体計測装置1では、3次元頭皮表面画像41や3次元脳表面画像42上に、測定関連点Sn、Mnと直線Ln−n’とを表示する構成を示したが、3次元テンプレート上に、測定関連点Sn、Mnと直線Ln−n’とを表示する構成としてもよい。
本発明は、非侵襲で脳活動を測定する光生体計測装置等に利用することができる。
1: 光生体計測装置
12: 送光プローブ
13: 受光プローブ
21: 送受光用制御部
22: 演算部
23: 相関係数算出部
12: 送光プローブ
13: 受光プローブ
21: 送受光用制御部
22: 演算部
23: 相関係数算出部
Claims (3)
- 被検者の頭皮表面上に配置される複数個の送光プローブと、
当該頭皮表面上に配置される複数個の受光プローブと、
前記送光プローブが頭皮表面に光を照射するとともに、前記受光プローブが頭皮表面から放出される光を検出するように制御することで、複数個の測定部位に関する複数個の受光量情報を取得する送受光用制御部と、
複数個の受光量情報に基づいて、複数個の測定データを取得する演算部と、
各測定データ間の相関の強さを示す相関係数値を算出する相関係数算出部とを備える光生体計測装置であって、
前記被検者を撮影することで取得された3次元脳表面画像を表示する3次元画像表示制御部と、
前記3次元脳表面画像上に、前記測定データが得られる測定関連点を表示するとともに、前記相関係数値の大きさを示す線分で測定関連点と測定関連点とを結ぶように表示する相関係数表示制御部とを備え、
前記3次元脳表面画像は、所望の方向から見た画像となるように、当該方向を変更可能として表示されることを特徴とする光生体計測装置。 - 前記線分は、前記相関係数値の大きさに応じて太さが変化する直線であることを特徴とする請求項1に記載の光生体計測装置。
- 前記測定関連点は、脳表面上では、前記送光プローブの位置と前記受光プローブの位置とを結んだ線の垂直二等分線と交差する位置とする脳表面上の測定関連点とすることを特徴とする請求項1又は請求項2に記載の光生体計測装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2012/079654 WO2014076795A1 (ja) | 2012-11-15 | 2012-11-15 | 光生体計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2014076795A1 JPWO2014076795A1 (ja) | 2016-09-08 |
JP6011636B2 true JP6011636B2 (ja) | 2016-10-19 |
Family
ID=50730735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014546785A Active JP6011636B2 (ja) | 2012-11-15 | 2012-11-15 | 光生体計測装置 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6011636B2 (ja) |
WO (1) | WO2014076795A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6661511B2 (ja) * | 2016-10-12 | 2020-03-11 | 株式会社日立製作所 | 脳結合性分析システム、及び脳結合性分析方法 |
WO2020162063A1 (ja) * | 2019-02-08 | 2020-08-13 | 株式会社島津製作所 | 脳機能計測装置 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5295491A (en) * | 1991-09-26 | 1994-03-22 | Sam Technology, Inc. | Non-invasive human neurocognitive performance capability testing method and system |
EP1272105A2 (en) * | 2000-03-30 | 2003-01-08 | The General Hospital Corporation | Method and apparatus for objectively measuring pain, pain treatment and other related techniques |
JP4230729B2 (ja) * | 2002-07-19 | 2009-02-25 | 株式会社日立メディコ | 生体光計測装置 |
JP2004184402A (ja) * | 2002-11-19 | 2004-07-02 | Hitachi Ltd | 生体光計測装置 |
JP2004187743A (ja) * | 2002-12-09 | 2004-07-08 | Hitachi Medical Corp | 医用3次元画像表示装置 |
JP4403453B2 (ja) * | 2003-11-13 | 2010-01-27 | 株式会社島津製作所 | 頭表座標を脳表座標に変換する方法と、その変換データを利用する経頭蓋的脳機能測定装置 |
JP4590555B2 (ja) * | 2004-09-02 | 2010-12-01 | 国立大学法人長岡技術科学大学 | 感性状態判別方法及び装置 |
US9101276B2 (en) * | 2006-07-06 | 2015-08-11 | Regents Of The University Of Minnesota | Analysis of brain patterns using temporal measures |
WO2008102825A1 (ja) * | 2007-02-20 | 2008-08-28 | Articell Systems Corporation | 遺伝子発現パターンから遺伝子を分類する方法 |
JP2009095380A (ja) * | 2007-10-12 | 2009-05-07 | Shimadzu Corp | 光生体計測装置 |
US10299710B2 (en) * | 2009-07-21 | 2019-05-28 | Shimadzu Corporation | Organism optical measurement device |
US9138146B2 (en) * | 2009-10-14 | 2015-09-22 | Hitachi Medical Corporation | Biological optical measurement instrument and reference image display method |
JP5693984B2 (ja) * | 2011-02-02 | 2015-04-01 | 株式会社トプコン | 測定装置、測定方法、およびプログラム |
TW201236646A (en) * | 2011-03-15 | 2012-09-16 | Ming-Da Wu | Method for allocating positive answer zone and negative answer zone of cerebral cortex |
-
2012
- 2012-11-15 JP JP2014546785A patent/JP6011636B2/ja active Active
- 2012-11-15 WO PCT/JP2012/079654 patent/WO2014076795A1/ja active Application Filing
Also Published As
Publication number | Publication date |
---|---|
JPWO2014076795A1 (ja) | 2016-09-08 |
WO2014076795A1 (ja) | 2014-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10299710B2 (en) | Organism optical measurement device | |
JP5896025B2 (ja) | 光生体計測装置及びそれを用いた解析方法 | |
JP5347448B2 (ja) | 生体測定装置 | |
US10052040B2 (en) | Photobiological measurement apparatus | |
JP6273207B2 (ja) | 光生体計測装置及びそれに用いられる位置計測装置 | |
JP4661688B2 (ja) | 光生体計測装置、光生体計測装置用プログラム及び光生体計測方法 | |
JP5248758B2 (ja) | 光計測装置 | |
JP5159893B2 (ja) | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 | |
JP6011636B2 (ja) | 光生体計測装置 | |
JP5790877B2 (ja) | 光生体測定システム及びその使用方法 | |
JP5262102B2 (ja) | 光計測装置 | |
JP4978485B2 (ja) | 光生体測定装置 | |
JP5946028B2 (ja) | 光生体計測装置 | |
JP7536186B2 (ja) | 乳癌診断システム及び乳癌診断システムの作動方法 | |
JP6154613B2 (ja) | 断面画像計測装置及び計測方法 | |
JP5729490B2 (ja) | 光生体計測装置 | |
JP5853777B2 (ja) | 光生体計測装置 | |
JP2013059437A (ja) | 光生体測定装置 | |
JP2014110844A (ja) | 生体計測装置及びそれに用いられる位置計測装置 | |
JP6281628B2 (ja) | 光計測システム | |
JP6069885B2 (ja) | 光計測システム及びその使用方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20160823 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160905 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6011636 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |