JP5944511B2 - 生体光計測装置およびそれを用いる生体光計測方法 - Google Patents
生体光計測装置およびそれを用いる生体光計測方法 Download PDFInfo
- Publication number
- JP5944511B2 JP5944511B2 JP2014532866A JP2014532866A JP5944511B2 JP 5944511 B2 JP5944511 B2 JP 5944511B2 JP 2014532866 A JP2014532866 A JP 2014532866A JP 2014532866 A JP2014532866 A JP 2014532866A JP 5944511 B2 JP5944511 B2 JP 5944511B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- distance
- measurement
- subject
- biological
- 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 208
- 238000000691 measurement method Methods 0.000 title claims description 4
- 230000003287 optical effect Effects 0.000 claims description 50
- 238000001514 detection method Methods 0.000 claims description 36
- 238000004458 analytical method Methods 0.000 claims description 32
- 210000004884 grey matter Anatomy 0.000 claims description 23
- 210000004556 brain Anatomy 0.000 claims description 16
- 238000000926 separation method Methods 0.000 claims description 16
- 239000010410 layer Substances 0.000 claims description 11
- 239000002344 surface layer Substances 0.000 claims description 5
- 230000001902 propagating effect Effects 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 3
- 230000031700 light absorption Effects 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 238000010586 diagram Methods 0.000 claims 1
- 239000000523 sample Substances 0.000 description 73
- 238000000034 method Methods 0.000 description 29
- 239000013307 optical fiber Substances 0.000 description 27
- 238000004497 NIR spectroscopy Methods 0.000 description 15
- 230000008859 change Effects 0.000 description 14
- 238000003780 insertion Methods 0.000 description 14
- 230000037431 insertion Effects 0.000 description 14
- 230000008326 skin blood flow Effects 0.000 description 11
- 230000035945 sensitivity Effects 0.000 description 10
- 238000002834 transmittance Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 210000004761 scalp Anatomy 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 238000012880 independent component analysis Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 230000003925 brain function Effects 0.000 description 4
- 239000000835 fiber Substances 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 102000001554 Hemoglobins Human genes 0.000 description 3
- 108010054147 Hemoglobins Proteins 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 210000003128 head Anatomy 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000003727 cerebral blood flow Effects 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000036992 cognitive tasks Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 210000001214 frontal sinus Anatomy 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000008557 oxygen metabolism Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- OGIDPMRJRNCKJF-UHFFFAOYSA-N titanium oxide Inorganic materials [Ti]=O OGIDPMRJRNCKJF-UHFFFAOYSA-N 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/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
- A61B5/1455—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 using optical sensors, e.g. spectral photometrical oximeters
-
- 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
- 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/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
- A61B5/4064—Evaluating the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/743—Displaying an image simultaneously with additional graphical information, e.g. symbols, charts, function plots
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/02—Operational features
-
- 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
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
-
- 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
- A61B5/1455—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 using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—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 using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
- A61B5/14556—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 using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases by fluorescence
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Neurology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Psychology (AREA)
- Neurosurgery (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は、可視光または近赤外光を用いた生体光計測装置に関する。
光トポグラフィ法を含む、近赤外分光法（Near−infrared spectroscopy: NIRS）を用いた無侵襲光脳機能イメージングで得られる光検出信号および生体信号（以下、NIRS信号）においては、光を頭皮上から照射するために、頭皮における皮膚血流変動の影響を受ける可能性があるとの報告がある。このような皮膚血流の影響を考慮して、それらの成分を抽出・除去する方法が研究されている。これらの多くは複数の照射−検出器（送光器−受光器）間距離（以下、SD（Source−Detector）距離）方式により深さの異なる部位からの信号成分を取得しており、これらを用いて浅い層部分の計測データに影響すると考えられる皮膚血流信号を除去しようとするものである。複数のSD距離において計測する方式を、以下ではマルチSD方式と呼ぶ。
マルチSD方式を利用した皮膚血流除去の方法としては、従来では、特許文献１において、皮膚血流等の信号を除去するため、被検者および計測部位に依存して、長いSD距離での計測信号から、短いSD距離での計測信号に適切な係数（scaling factor）を乗じたものを減算する方法が開示されている。しかし、この方法は、マルチSD計測データが使用可能であることが前提であり、各SD距離で適切な信号が取得できていない場合、例えば光量が大きいことによる検出器の飽和時には適用困難となる。
また、特許文献２において、光照射部から、計測対象によって吸収される波長を有する光を、その強度を時間的に変化させながら被検体内に向けて照射し、被検体を透過した光は、走査方向に沿って配列される複数の検出素子により検出する方法が開示されている。この方法は、１次元方向の計測に適しているが、脳機能計測等で用いられる２次元イメージングには適用困難である。
さらに、特許文献３において、NIRS信号に含まれる皮膚血流の影響を分離、除去し、脳または大脳皮質由来の信号を抽出するために、複数のSD距離による計測を実現するよう、かつ、各受光器が受光する光が灰白質を伝播するよう送光器および受光器を配置し、各計測点におけるデータを使用して独立成分分析（Independent component analysis、ICA）を実施し、各分離成分の重み値のSD距離依存性を用いることで、各々の独立成分が脳由来か皮膚由来かを判定する方法が開示されている。この方法においても、マルチSD計測データが使用可能であることが前提であり、安定した信号取得が課題となる。
非特許文献３記載のように、マルチSD方式で、ある送光器あるいは受光器を、異なる複数のSD距離の計測に共通に使用するときには、時分割で受光感度を調整する方法、および時分割で光源のパワーを変化させる方法もあるが、マルチSDではプローブが高密度に配置されることになるので、時分割で計測していては、時間分解能が低下するという課題がある。さらに、時分割で計測する場合でも、例えばSD距離5 mmと30 mmの計測を同時に行う場合には、図３に示すように４桁以上の検出光量の違いがあり、電磁気的作用による制御では、高精度で多段の増幅回路等が必要となり、生産コストの増大を引き起こす。
さらに、マルチSD計測時には、SD距離の種類が複数あるとともに、各SD距離の計測点が空間的に均等に配置されることで、各計測点における感度一様化の設計が効率化する。計測対象を２次元あるいは３次元に画像化する解析においても、計測点が一様に分布することは重要であり、そのような計測を実現するための送光器および受光器の適切な配置が必要である。
光の透過率は、SD距離、被検者（年齢、性別）、計測部位等の各種条件により異なる。特に、ヒトの前額部においては、副鼻腔（もしくは前頭洞）の個人差の影響により、平均光路長および透過率が大きく変化する（特許文献４参照）。マルチSD方式による光脳機能計測においては、計測時の各条件に応じて、光量を調整したり、光学フィルタまたは光減衰器挿入等により対処したりする必要があった。光学フィルタを用いない場合、検出器のダイナミックレンジ以上の複数光源からの検出光量差に対し、検出器の適切なゲイン設定ができない可能性があった。光学フィルタを用いる場合には、光ファイバをプローブ側または装置本体側で一度外す必要があり、手間がかかるという問題があった。さらに、一つの検出器が、同時に複数のSD距離の位置にある光源からの光を受光する場合、光量の異なる複数の信号を安定して検出するために、ダイナミックレンジの広い検出器を用いるか、光学フィルタをSD距離もしくは検出光量に応じて光源側に挿入する必要があった。
A. Maki et al.、"Spatial and temporal analysis of human mot or activity using noninvasive NIR topography"、Medical Physics、 Vol.22、 No.12、 pp.1997−2005（1995）
T. Funane et al.、マルチディスタンスプローブ及び独立成分分析を用いた光脳機能計測における皮膚血流影響の分離手法（Discrimination of Skin Blood Flow Using Multi−Distance Probe and Independent Component Analysis in Optical Brain Function Monitoring）、電気学会研究会資料（The papers of Technical Meeting on Optical and Quantum Devices、IEE Japan）、OQD−11−033、pp. 17−22 (2011)
I. Oda et al.、"Near Infrared Imager with a Flexible Source−Detector Arrangement and a New Detection Gain Control"、Optical Review 10(5)、pp. 422−426 (2003).
しかし、光学フィルタを送光器もしくは検出器に用いると、その送光器または検出器は、光の減衰のため長いSD距離では信号対雑音比(S/N)が小さく計測が困難となり、短いSD距離の計測にしか用いることができないという問題が生じる。また、光学フィルタを装着したプローブ（送光器、または受光器）を用いて長いSD距離の計測を行う場合には、経時的に光検出器の感度を変化させる（非特許文献３参照）、または経時的に光量を変化させる必要がある。
本発明は、近赤外分光技術（Near−infrared spectroscopy: NIRS）に基づき、生体組織の表層および深層における光吸収変化を分離するために複数の照射−検出器間距離を有するマルチディスタンスプローブを用いた生体光計測装置において、送光または受光する光の減衰量を調整することで、受光感度を各被検者・各部位で安定化することを目的とする。
上記課題を解決するため、本発明の生体光計測装置は、光減衰量調整手段を有し、送光または受光する光の減衰量を調整することで、受光感度を安定化させる。
本発明の生体光計測装置の一例を挙げるならば、被検体に光を照射するための１つまたは複数の光照射手段と、前記光照射手段から前記被検体上の照射点に照射され、前記被検体内を伝播してきた光を前記被検体上の検出点において検出するための１つまたは複数の光検出手段と、前記光照射手段に含まれる光源素子と前記被検体との間、もしくは前記光検出手段に含まれる受光素子と前記被検体との間の光の伝播路上に設置される光減衰量調整手段と、前記光検出手段で得られる信号を解析するための解析部と、前記解析部での解析結果を表示するための表示部とを有し、前記光照射手段と前記光検出手段の各々は、前記被検体上における、前記照射点と前記検出点間の距離として定義されるSD距離が２種以上となるように前記被検体上に配置され、前記解析部は、少なくとも１つの前記光検出手段が検出する２種以上のSD距離の位置にある前記光照射手段からの光による信号を各々解析し、各々の受光量を所定の範囲内とするための光減衰調整量を算出し、前記光減衰量調整手段は、前記光検出手段が検出する２種以上のSD距離の位置にある前記光照射手段からの光の減衰量を、前記光検出手段に入射する光量を変化させることにより、各々調整可能としたことを特徴とするものである。
本発明の生体光計測装置を用いた生体計測方法は、ＳＤ距離が１０ｍｍ未満、好ましくは８ｍｍ以下の少なくとも１つの計測点、およびＳＤ距離が１０ｍｍ以上、好ましくは１２ｍｍ以上の２つ以上の計測点において被検体内を伝播してきた光を計測するステップと、ＳＤ距離と信号分離手法により得られた各成分の重み値とをプロットする図において、ＳＤ距離が１０ｍｍ以上の２つ以上の計測点の重み値を用いて１次直線を求め、重みゼロに対応するＳＤ距離を求めるステップと、ＳＤ距離が１０ｍｍ未満の計測点の重み値を用いて軸に平行な直線を求めるステップと、前記軸に平行な直線と前記１次直線との交点のＳＤ距離を、灰白質到達最小ＳＤ距離として求めるステップと、前記重みゼロに対応するＳＤ距離と前記灰白質到達最小ＳＤ距離から、脳寄与率を算出するステップとを含むものである。
本発明によれば、送光または受光する光の減衰量を調整することで、被検者個人間や計測部位間の透過率の違いを考慮した上で、受光信号の安定化を実現できる。また、各送光器および受光器の導波路における光減衰量を、自動で、もしくは装置のオペレータが容易に手動で調整可能となり、時間的、人的コストを低減した効率的な信号取得を実現できる。
以下、図面を用いて本発明の実施の形態を説明する。
図１に、本発明の生体光計測装置の装置構成の一例を示す。光を生体に入射し、生体内で散乱・吸収され伝播して出てきた光を検出する生体光計測装置において、装置本体２０に含まれる１つまたは複数の光源１０１から照射される光３０を、導波路４０を介して、被検者１０に入射させる。光３０は、照射点１２から被検者１０内に入射し、被検者１０内を透過、伝播した後は、照射点１２とは離れた位置にある検出点１３から導波路４０を介して、１つまたは複数の光検出器１０２で検出される。SD距離は、前述のように照射点１２と検出点１３間の距離で定義される。
ここで、１つまたは複数の光源１０１は半導体レーザ（LD）や発光ダイオード（LED）等であり、１つまたは複数の光検出器はアバランシェフォトダイオード（APD）やフォトダイオード（PD）、光電子増倍管（PMT）等である。また、導波路４０は光ファイバ、ガラス、ライトガイド等である。
光源１０１は、光源駆動部１０３により駆動され、１つまたは複数の光検出器１０２のゲインは制御・解析部１０６により制御される。制御・解析部１０６は、光源駆動部１０３の制御も行い、入力部１０７からの条件等の入力を受ける。
光検出器１０２で光電変換した電気信号は、増幅器１０４で増幅され、アナログ−デジタル変換器１０５でアナログ−デジタル変換され、制御・解析部１０６へ送られ、処理される。
ここで、１つまたは複数の光源１０１は半導体レーザ（LD）や発光ダイオード（LED）等であり、１つまたは複数の光検出器はアバランシェフォトダイオード（APD）やフォトダイオード（PD）、光電子増倍管（PMT）等である。また、導波路４０は光ファイバ、ガラス、ライトガイド等である。
光源１０１は、光源駆動部１０３により駆動され、１つまたは複数の光検出器１０２のゲインは制御・解析部１０６により制御される。制御・解析部１０６は、光源駆動部１０３の制御も行い、入力部１０７からの条件等の入力を受ける。
光検出器１０２で光電変換した電気信号は、増幅器１０４で増幅され、アナログ−デジタル変換器１０５でアナログ−デジタル変換され、制御・解析部１０６へ送られ、処理される。
制御・解析部１０６では、光検出器１０２で検出された信号に基づき解析を実行する。具体的には、アナログ−デジタル変換器１０５で変換して得られたデジタル信号を入力し、当該デジタル信号をもとに、例えば非特許文献１に記載されている方法に基づいて、検出光量変化もしくは吸光度変化から、酸素化、脱酸素化ヘモグロビン濃度長変化（oxy−Hb、deoxy−Hb）を算出する。ここで、濃度長変化とは濃度と光路長の積の変化量である。
ここでは、制御・解析部１０６は光源１０１の駆動、光検出器１０２のゲイン制御、アナログ−デジタル変換器１０５からの信号処理を全て行うことを想定して記述したが、それぞれ別個の制御部を有し、さらにそれらを統合する手段を有することでも同機能を実現できる。
また、計測データおよびヘモグロビン濃度長変化算出結果は、記憶部１０８に保存され、解析結果および／または保存データに基づいて表示部１０９で計測結果を表示することが可能である。
ここでは、制御・解析部１０６は光源１０１の駆動、光検出器１０２のゲイン制御、アナログ−デジタル変換器１０５からの信号処理を全て行うことを想定して記述したが、それぞれ別個の制御部を有し、さらにそれらを統合する手段を有することでも同機能を実現できる。
また、計測データおよびヘモグロビン濃度長変化算出結果は、記憶部１０８に保存され、解析結果および／または保存データに基づいて表示部１０９で計測結果を表示することが可能である。
図２に示す送光器５０、受光器６０は、図１には記載していないが、送光器５０は、例えば光源１０１側の導波路４０を含み、被検者１０に接触あるいは接触に近い状態で設置され、受光器６０は、例えば光検出器１０２側の導波路４０を含み、被検者１０に接触あるいは接触に近い状態で設置される。
図２に、本発明の生体光計測装置により照射された光が生体組織内を伝播する状況について説明する。送光器５０が１つ、受光器６０が２つ、頭皮上に配置された例である。SD距離は１５ｍｍ、３０ｍｍの２種類である。SD距離の違いにより、照射光の到達深度（penetration depth）が異なる。つまり、平均実効光路は、SD距離１５ｍｍで浅く、SD距離３０ｍｍで深くなる。このような配置のプローブ（以下、送光器５０、受光器６０を合わせてプローブと呼び、さらに、複数のSD距離を有するプローブを特に、マルチディスタンスプローブまたはマルチSDプローブと呼ぶ。）を用いることにより、深層信号から表層信号の影響を除去することが可能となる。
図３に、典型的な頭部モデルにおけるSD距離と光子透過率との関係を示す。図３に示すようにSD距離により光子透過率は指数関数的に変化する。生体計測時には、光検出器１０２が飽和せずに、かつ十分な信号対雑音比で検出可能な光量とする必要がある。さらに、計測される被検者、計測部位により、組織構造が異なるため、光の透過率は変化する。多くの被検者および部位に対して使用可能な装置とするために、本発明の生体光計測装置は、図１に示す光減衰量調整手段１４を有し、複数の前記光検出器１０２で検出される光量を所定の範囲内になるように光減衰量を調整する。
図４に、被検者および光減衰量調整量表示を含めた本発明の実施の形態を示す。外部筺体に光源（半導体レーザ）、光検出器（フォトダイオード）を含む場合である。光源および検出器内臓回路２０２内の光源１０１が発した光３０は光ファイバ４０３を通りプローブホルダ５０３に固定された送光器５０より被検者１０に対して照射される。被検者１０内を伝播した光３０は、プローブホルダ５０３に固定された受光器６０から光ファイバ４０３を通り光源および検出器内臓回路２０２内の光検出器１０２で受光される。受光信号（検出光量）は電気信号２０１として装置本体２０に送られ、図に示していない制御・解析部１０６において検出光量が所定の範囲内かどうかを判定し、光減衰量の調整量２０３（例えば、＋１、−１、＋２、０、等の相対値）を表示部１０９により表示する。＋１等は、例えば１０の１乗、＋２は１０の２乗を表し、光減衰量調整手段１４において、これだけの調整量を自動または手動で操作し、光減衰量を調整可能とする。自動で光減衰量を調整する場合には、装置本体２０から電気信号２０１を介して光減衰量調整手段１４に制御信号を送る。手動で光減衰量を調整する場合には、オペレータは表示部１０９を見ながら、光減衰量調整手段１４を手動で操作する。この構成により、プローブホルダ５０３を再度着脱することなく、高速に光減衰量を調整することができる。さらにオペレータおよび被検者への負担軽減の効果もある。ここでは相対値の表示について説明したが、絶対値の表示でも良い。さらに、部位に応じて、典型的な光減衰量設定値、もしくは当該被検者の過去の設定値等の所定の設定値を自動で設定、または手動で設定できるよう、表示部１０９で表示しても良い。 このように、光減衰量調整手段を用いることで、プローブの再着脱のみでは実現困難な大きな光減衰を容易に調整可能となり、計測時間、コストの低減につながるという効果がある。
図５に、被検者および光減衰量調整量表示を含めた本発明の別の実施の形態を示す。外部筺体に光ファイバの中間コネクタを有する場合である。図４では、光源および検出器内臓回路２０２内で光の生成、検出光の光電変換を行ったが、図５の構成では、光減衰量調整手段１４が装置本体側およびプローブホルダ側の光ファイバ４０３に挟まれるように、一方または両方の光ファイバ４０３の接合部に配置される。光減衰量調整手段１４は送光ファイバ、受光ファイバのどちらに設置されても良い。光減衰量調整手段１４は、光散乱、吸収、反射等の光学的作用により、光ファイバ４０３を伝播する光量を調整する。他の構成は図４の場合と同様である。この構成により、既存装置の光ファイバに対し、独立したモジュールとして光減衰量調整手段１４を付加することができる。
次に、図６を用いて図４に対応する光減衰量調整手段１４を、説明する。図６に、光減衰量調整部と周囲のコンポーネントとの関連を示す。筺体２０４内において、光減衰量調整手段１４は、同じ筺体２０４内に配置される半導体レーザ２０６またはフォトダイオード２０７等の光学素子と光ファイバ４０３との間に設置され、各々の素子に接続する光ファイバ４０３との間の接続損失を個別に調整する。ここでは、電磁的手段として、自動または手動で調整可能なピエゾアクチュエータもしくは圧電素子を備え、光ファイバ４０３間の端面位置関係を微小に変化させる。このように装置本体２０とは異なる外部の筺体２０４に光学素子、および光減衰量調整手段１４がまとめられることで、より被検者１０に近い位置で、プローブ装着と同時に受光感度調整を行うことができ、計測を効率化できるという効果がある。なお、ここでは電磁的手段について説明したが、モータ、歯車、もしくはラック・ピニオン機構、油圧、空気圧等を用いた機械的手段でも良い。また、自動、手動のどちらの手段でもよい。
図７に、光減衰量調整手段（手動調整つまみ）と表示部（発光ダイオードによる色表示）の具体例を示す。図６の構成をより具体化したもので、光減衰量調整手段１４は、手動調整つまみ２０５を有しており、さらに、表示部１０９として、調整量を色による表示２０８で示したものである。なお、色による表示２０８では、白黒もしくは所定の色の濃淡、つまみの方向を表す矢印の長さ等を表示してもよい。オペレータが操作する手動調整つまみ２０５が含まれる同じ筺体２０４に表示部１０９を有することで、本構成の生体光計測装置は、操作方法が単純化されオペレータの操作時のミスを防ぐ効果がある。さらに、操作結果を早急にオペレータに知らせることができるため、計測を効率化する効果がある。
図８に、光減衰量調整手段の一例として、複数の挿入口を持つ光ファイバ接続用の中間コネクタの構造を示す。中間コネクタ２０９は、２本の光ファイバ４０３を光学的に結合するために用いられ、光ファイバ挿入口（減衰量が小さい）２１０、光ファイバ挿入口（減衰量が中程度）２１１、光ファイバ挿入口（減衰量が大きい）２１２をそれぞれ２つずつ、対面させて備えており、両側に光ファイバ４０３を１本ずつ挿入することで、光ファイバ端面間の距離および位置関係を変化させ、当該光ファイバ４０３における光減衰量を調整することが可能となる。この構成により、手動で容易に光減衰量を調整可能となる。中間コネクタの各挿入口の種類を増やすことで、より広い範囲で光減衰量を調整することも可能となる。
図９に、光減衰量調整手段の他の一例として、段階的に挿入できる中間コネクタの構造を示す。中間コネクタ２０９は、挿入深さ調整部２１３を有し、光ファイバ４０３間の端面間距離を数段階で調整可能である。ここでは、段階的に光ファイバ４０３を挿入する代わりに、ネジ状構造を設けて回転させながら連続的に端面間距離を調整する方法でもよい。この構成により、より小型の装置で光減衰量の調整が可能となる。
図１０に、光減衰量調整手段の他の一例として、穴サイズにより光減衰量を調整する機構を示す。図１０の中間コネクタ２０９は、挿入口２１４を有する。導波孔（減衰量が大きい）２１６、導波孔（減衰量が中程度）２１７、導波孔（減衰量が小さい）２１８を有する媒体もしくは金属部材を挿入口２１４より挿入し、光ファイバ４０３を伝播する光の導波路上に穴の径の異なる導波孔を配置することで、光減衰量を制御する。この操作は自動でも手動でもよい。
図１１に、光減衰量調整手段の他の一例として、媒質により光減衰量を調整する機構を示す。図１１の中間コネクタ２０９には、図１０における孔を有する媒体の代わりに、吸収または散乱媒体２１５を挿入する。吸収または散乱媒体２１５は複数準備しておき、実現したい光減衰量に相当するものを使用する。この操作は自動でも手動でもよい。ここでの媒体としては、エポキシ樹脂に、染料等の吸収物質および酸化チタン等の散乱物質を混合したもので良く、同一部材内の所定の面上における位置により光学特性が滑らかに変化する様、各物質濃度のグラデーションのある媒体としても良い。吸収物質および散乱物質の濃度を調整することで、光減衰量を調整可能である。さらに、ここでの媒体の代わりに、電磁的に調整可能で、光の透過率を変えられる液晶画面、またはエレクトロクロミック材料からなる調光ガラス及び調光フィルムを光路上に設置し、各画素のON/OFFを切り替えて適当な大きさおよび形状の模様を作ることで、光透過度を調整してもよい。さらに、熱もしくは温度で透過率を制御できるサーモクロミック材料を、温度制御手段と共に用いてもよい。
これら図１０、１１の構成により、広い範囲で任意の光減衰量を実現することができ、被検者、計測部位によらずに安定した受光感度が得られるという効果がある。
図１２に、検出光量または調整量を表示するフローチャートを示す。図１２のフローチャートを参照しながら、実施例１に従った調整量表示手順を説明する。まず、プローブを被検者１０に装着する（ステップS１０１）。次に、レーザを照射し検出光量を計測する（ステップS１０２）。検出光量が所定の光量（例えば、A/D変換後に0.1-2 V）の範囲内かを判定する（ステップS１０３）。検出光量が0.1-2 Vの範囲内と判定された場合（YES）には、検出光量が適正であるため、表示部において、「正常」と表示する（ステップS１０４）。ステップS１０３において、検出光量が0.1 V未満か2 Vより大きい場合（NO）、検出光量が適正でないため、さらに閾値（例えば、2 V）より大きいかを判定する（ステップS１０５）。検出光量が2 Vより大きいと判定された場合（YES）には、検出光量を小さくする必要があるため、それを示すために「大きい」と表示する（ステップS１０６）。検出光量が2 V以下であると判定された場合（NO）、ステップS１０３での判定と合わせて検出光量は0.1 V未満であるから、検出光量を大きくする必要があり、それを示すために「小さい」と表示する（ステップS１０７）。これらの表示では、ステップS１０３およびS１０５の判定結果を、色の違いで表示してもよい。さらには、どれだけ大きいか、または小さいかを、段階的に数値（例えば、＋１、＋３、−２等）および色の違いで表示してもよい。さらに、これらは光減衰量調整手段による調整量として表示してもよい。
図１２に、検出光量または調整量を表示するフローチャートを示す。図１２のフローチャートを参照しながら、実施例１に従った調整量表示手順を説明する。まず、プローブを被検者１０に装着する（ステップS１０１）。次に、レーザを照射し検出光量を計測する（ステップS１０２）。検出光量が所定の光量（例えば、A/D変換後に0.1-2 V）の範囲内かを判定する（ステップS１０３）。検出光量が0.1-2 Vの範囲内と判定された場合（YES）には、検出光量が適正であるため、表示部において、「正常」と表示する（ステップS１０４）。ステップS１０３において、検出光量が0.1 V未満か2 Vより大きい場合（NO）、検出光量が適正でないため、さらに閾値（例えば、2 V）より大きいかを判定する（ステップS１０５）。検出光量が2 Vより大きいと判定された場合（YES）には、検出光量を小さくする必要があるため、それを示すために「大きい」と表示する（ステップS１０６）。検出光量が2 V以下であると判定された場合（NO）、ステップS１０３での判定と合わせて検出光量は0.1 V未満であるから、検出光量を大きくする必要があり、それを示すために「小さい」と表示する（ステップS１０７）。これらの表示では、ステップS１０３およびS１０５の判定結果を、色の違いで表示してもよい。さらには、どれだけ大きいか、または小さいかを、段階的に数値（例えば、＋１、＋３、−２等）および色の違いで表示してもよい。さらに、これらは光減衰量調整手段による調整量として表示してもよい。
本実施例によれば、近赤外分光技術（Near−infrared spectroscopy: NIRS）に基づき、生体組織の表層および深層における光吸収変化を分離するために複数の照射−検出器間距離を有するマルチディスタンスプローブを用いた生体光計測装置において、送光または受光する光の減衰量を調整することで、受光感度を各被検者・各部位で安定化することができる。特に、光源１個に対して光検出器を２個用いる計測や、光検出器１個に対して光源を２個用いる計測や、両者が混在した計測の場合のような、複数のSD距離の計測点が混在する際に、各計測点間での干渉を防いで受光感度を安定化することができる。このように共通の光源または光検出器を用いて、SD距離の異なる複数の計測点を用いて計測するときに、共通して用いられる光源のパワー調整や光検出器のゲイン調整では、マルチSD計測で生じる大きい信号振幅比に対応することは困難であるが、本実施例では、各光源および各光検出器の導波路における光減衰量を、自動で、もしくはオペレータが手動で調整可能となり、効率的な信号取得が実現できる。
次に、制御・解析部１０６で行う、マルチSD計測データを用いた脳寄与率算出方法について説明する。特許文献３においては、信号振幅もしくは信号分離手法により得られた成分の重みのSD距離依存性を利用した、脳および皮膚由来信号の分離手法について記載されている。図１３に、 モンテカルロシミュレーションにより求めた、(a) SD距離と、頭皮および灰白質の部分平均光路長、(b) SD距離と、頭骨、脳脊髄液層（CSF）の部分平均光路長および総光路長との関係をそれぞれ示す。横軸はSD距離(mm)、縦軸は光路長(mm)である。NIRS計測においては、SD距離は照射光の到達深度に影響するため、SD距離を変えることで、対象とする生体組織深度を変えることができる。マルチSD計測では、各SD距離で、浅部および深部信号が、異なる重みで得られることから、複数点でデータ取得することで、信号分離手法により、各層由来の信号の寄与を算出可能となる。具体的には、頭皮の部分光路長には強いSD距離依存性は見られないが、灰白質には線形的なSD距離依存性が見られる。この方法は、NIRS信号強度は血流変化が生じる部位の部分光路長に比例する（例えば、非特許文献１参照）という原理に基づき、SD距離が大きくなると、NIRS計測信号における脳由来成分は大きくなるが、皮膚由来成分は変化しないことが予想される。
図１４に、３種類のSD距離における計測データから灰白質到達最小SD距離を推定する方法を幾何的に示す。この分離手法の適用には、マルチSD計測データが必要であり、このような計測においても、実施例１で説明した光減衰量調整手段により、安定した信号を得ることが可能となる。計測データを用いた解析時において、SD距離８ｍｍ以下の計測点１点、SD距離１２ｍｍ以上の計測点２点において計測されたNIRSデータに独立成分分析を適用し、得られた独立成分の重み６０１をSD距離に対してプロットする。このとき、信号分離手法を用いずに、所定の認知課題遂行中のNIRS信号の振幅等、代表的な振幅をプロットしても良い。図１３（ａ）に示す平均実効光路長のモデルを仮定し、近似的にSD距離約8 mm以降で皮膚の平均実効光路長は一定値、SD距離約10 mm以降で灰白質の平均実効光路長は1次関数的に変化すると仮定する。さらに、各分離成分が灰白質（脳）と皮膚（頭皮）における血流変化の両方の影響を受けている（混合成分）と仮定すれば、得られる信号の振幅は、両層の光路長に、各層の濃度変化で重みづけされた重ね合わせで表されるため、SD距離約8 mmから約10 mmまで一定値で、約10 mm（灰白質到達最小SD距離６０３）以降は線形的に増加することが期待される。このような、混合成分のSD距離依存性に対する定性的な特徴を利用し、計測データから灰白質到達最小SD距離６０３を推定する手法を以下に説明する。
まず、図１４において、SD距離８ｍｍ以下の計測点については傾き０の直線（ｘ軸に平行な直線、Ｌ２）６０７を引き、SD距離１２ｍｍ以上の計測点に１次直線（Ｌ１）６０６をあてはめ、それらの２直線の交点のx座標がＸ sgray（灰白質到達最小SD距離６０３）となる。この灰白質到達最小SD距離６０３をx軸上で通り、傾きが前記1次直線６０６と同じ直線は、灰白質平均実効光路長に比例する直線６０４であり、前記1次直線６０６と直線６０４の振幅比から、各SD距離での脳寄与率を算出できる。このとき、使用するSD距離は１点が８ｍｍ以下、２点が１２ｍｍ以上であるが、１点が１０ｍｍ未満、２点を１０ｍｍ以上から選択してもほぼ同様の解析が可能である。ただし、灰白質到達最小SD距離６０３が標準的なヒトの場合に１０ｍｍ前後のため、より一般的に述べれば、１点を、光が灰白質に到達しないSD距離から、２点を、光が灰白質に到達するSD距離から選択する必要がある。ここで設定される不使用SD距離帯６０５は、図１３（ａ）において灰白質の平均実効光路長が0 mmに近く、皮膚の平均実効光路長の変化が小さい、例えば８〜１２ｍｍが設定されるが、脳寄与率算出時の精度に関連する重要なパラメータとなるため、より正確な幅を決定するには、各被検者のMRI構造データ、もしくはＸ 線CT（コンピュータ断層撮影、Computed tomography）データ等を利用することが重要となる。
ここで述べた構成により、従来法では灰白質到達最小SD距離６０３を仮定する必要があったところが、実測より被検者毎、部位毎に推定できるようになり、より高精度な脳寄与率算出が期待できる。なお、ここでは信号分離手法として、独立成分分析を用いる方法を説明したが、主成分分析、因子分析等の方法でも良い。また、信号分離手法も用いなくても、検出光量またはヘモグロビン変化の代表振幅を用いる方法でも良い。
図１５に、３種類のSD距離における計測データから推定される灰白質到達最小SD距離Ｘsgrayを用いて、NIRS信号の脳寄与率を算出するフローチャートを示す。SD12 mm程度以上の重み値を用いて直線（L1）を計算し、x切片をＸ sとする（ステップS２０１）。次に、SD8 mm程度以下の重み値を用いてx軸に平行な直線（L2）を計算する（ステップS２０２）。次に、L1、L2の交点のx座標をＸ sgrayとする（ステップS２０３）。最後に、Ｘ sとＸ sgrayから、当該成分の脳寄与率を算出する（ステップS２０４）。脳寄与率(r)算出時には、例えば、非特許文献２に記載されているように、
r ＝ (x − Ｘ sgray)/(x − Ｘ s)
の式を用いる。これは、前記１次直線６０６と直線６０４の振幅比を式で表したものである。
r ＝ (x − Ｘ sgray)/(x − Ｘ s)
の式を用いる。これは、前記１次直線６０６と直線６０４の振幅比を式で表したものである。
本発明における光減衰量調整手段１４を効果的に適用可能である、マルチSD計測を実現するための具体的なプローブ配置の例を以下に図を用いて説明する。なお、プローブ配置はここに示したものに限らず、複数のSD距離で計測可能であれば良く、ユーザの目的に応じてあらゆるプローブ配置が考えられる。
図１６に、２×１１配置をベースとし、SD距離5、15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。
図１７に、２×１１配置をベースとし、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。
図１８、図１９、図２０に、２×８配置をベースとし、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。できるだけ、各SD距離の計測点が空間的に均一に分布するように考慮した配置である。
図２１に、SD距離30 mmの光源、検出器のペアを、菱形の頂点上（もしくは水平線に対して４５度傾いた正方形の頂点上）に配置し、さらに横一列に同形の菱形が接するようにプローブを並べて配置してSD距離15 mm、30 mmの計測点をそれぞれ24、12点計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。図２１の配置は、図１６〜図２０の配置に比べて、計測領域は狭くなるが、SD距離30 mmの各計測点に2個ずつSD距離15 mmの計測点があり、さらにその計測点間距離が全て7.5 mmと近くなっている。よって狭範囲の計測に適している。
図２２に、SD距離30 mmの光源、検出器のペアを、菱形の頂点上（もしくは水平線に対して４５度傾いた正方形の頂点上）に配置し、さらに横一列に同形の菱形が接するようにプローブを並べて配置してSD距離15 mm、30 mmの計測点をそれぞれ20、16点計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。
図２３に、３×５配置をベースとし、SD距離5、15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。
図２４に、４×４配置をベースとし、SD距離5、15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）計測点配置を示す。
図１６〜図２４において、各々の送光器―受光器で計測される領域において、SD距離30 mmの計測点５０７と、SD距離15 mmの計測点５０８と、SD距離5 mmの計測点５０９とができるだけ均等に配置されるよう、送光用プローブ５０１、受光用プローブ５０２、送光用プローブと受光用プローブの一体型プローブ５０４が配置されている。各々の計測部位において、所定の範囲を有する領域を計測単位と定義して、計測単位内の各SD距離のデータを同時に用いて解析するため、計測単位内では、できるだけ近い位置に各SD距離の計測点を配置するのが良い。この場合、あるSD距離（例えば、３０ｍｍ）を優先SD距離と定め、当該SD距離に対応する送光器５０、受光器６０で計測される領域または代表される位置（例えば、送光器５０と受光器６０の、被検体１０表面における中点位置、これを優先計測点と呼ぶ）からの距離が、所定範囲内（例えば、２０ｍｍ以内）に他のSD距離の計測点（従属計測点と呼ぶ）が配置されるように送光用プローブ５０１、受光用プローブ５０２の配置を決定する。ここでの所定範囲は、同一範囲を計測している計測点で解析するのが望ましいので、できるだけ小さい（つまり、優先計測点と従属計測点の距離が小さい）方が望ましい。しかし、装置のプローブ数の制約等のため、離れて配置せざるを得ない場合も考えられ、その場合には、４０ｍｍ以内等に設定する。図１６〜図２４では、異なるSD距離の計測点間距離の最大値は、それぞれ33.5 mm、30.1 mm、17 mm、16.8 mm、15 mm、7.5 mm、19 mm、15.81 mm、15.81 mmである。また、図１６でSD距離5 mm、15 mm、30 mmでの計測点数はそれぞれ５、１２、３１点、図２３でそれぞれ８、５４、２２点、図２４で８、６０、２４点である。図１７でSD距離15 mm、30 mmでの計測点数はそれぞれ１６、３１点、図１８でそれぞれ３１、２２点、図１９でそれぞれ２４、２２点、図２１でそれぞれ２４、１２点、図２２でそれぞれ２０、１６点、図２３でそれぞれ８、５４、２２点、図２４でそれぞれ８、６０、２４点である。図２０では15 mm、30 mmでの計測点数とSD距離16.8 mmの計測点５１０の数がそれぞれ１２、２２、８点である。
なお、SD距離15 mm、30 mmのみのプローブ配置の場合でも、例えば検出器を1個追加して、5 mmの計測点を1点追加する構成にしても良い。このSD距離5 mmの計測点は、脳、皮膚由来信号分離後に、分離性能の評価に用いることができる。また、皮膚血流が対象となる計測範囲で一様と仮定したときに、SD距離5 mmの計測点を参照信号として利用して、皮膚血流成分を分離する解析に用いることができる。
入力部１０７で、光照射手段と光検出手段の位置関係を入力し、解析部１０６は、入力された光照射手段と光検出手段の位置関係から、優先計測点と従属計測点の組み合わせを算出する。そして、表示部１０９は、同一計測単位内の優先計測点と従属計測点の組み合わせを表等にして表示し、さらに、１つまたは複数の優先計測点と１つまたは複数の従属計測点からなる計測単位をデータセットとして、各データセットでの解析結果を表示する。これにより、計測単位の範囲の妥当性判断や、優先計測点と従属計測点との間の距離の妥当性の判断に役立つとともに、深部信号の寄与率算出から、被検者１０の構造的特徴を取得することができるという効果がある。
以下に、SD距離15、30 mmの計測点間の距離を小さくすることを目的とした各種プローブ配置の例を示す。
図２５に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置A）。
図２６に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置B）。
図２７に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置C）。
図２８に、SD距離7.5、15、30 mmを計測するときの、（ａ）プローブ配置、（ｂ）各計測単位における各SD距離の計測点個数の表を示す。
なお、図２５〜図２８において、符号５０１は送光用プローブを、符号５０２は受光用プローブを、符号５０５は送光用プローブ(減衰：大）を、符号５０６は受光用プローブ(減衰：大）を示す。また、符号５０７はＳＤ距離３０ｍｍの計測点を、符号５０８はＳＤ距離１５ｍｍの計測点を示す。
図２５に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置A）。
図２６に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置B）。
図２７に、SD距離15、30 mmを計測するときの（ａ）プローブ配置、（ｂ）SD距離30 mmの計測点配置、（ｃ）SD距離15 mmの計測点配置を示す（配置C）。
図２８に、SD距離7.5、15、30 mmを計測するときの、（ａ）プローブ配置、（ｂ）各計測単位における各SD距離の計測点個数の表を示す。
なお、図２５〜図２８において、符号５０１は送光用プローブを、符号５０２は受光用プローブを、符号５０５は送光用プローブ(減衰：大）を、符号５０６は受光用プローブ(減衰：大）を示す。また、符号５０７はＳＤ距離３０ｍｍの計測点を、符号５０８はＳＤ距離１５ｍｍの計測点を示す。
図２５〜図２８において、異なるSD距離の計測点間距離の最大値は、それぞれ7.5 mm、5 mm、15 mm、11 mmとなっている。図２８では特に、各計測単位（ｃｈ）において、7. 5、15、30 mmの各SDの計測点が最低1個含まれる配置となっており、脳寄与率算出等の目的には有用である。図１６から図２８に示したプローブ配置は、システムが有するプローブ数、対象となる生体の領域、適用する解析手法によって適切に選択する必要がある。例えば、光脳機能計測の場合には、脳血流および皮膚血流の影響が広域性を有するものであれば、異なるSD距離の計測点間距離の最大値は大きくても良く、脳血流および皮膚血流の変化が局在性を有するものであれば、異なるSD距離の計測点間距離の最大値は、活動領域の期待される大きさより小さくすることが望ましい。
図２９に、表示の例として、各計測点およびSD距離における検出光量結果を示す。上側にSD距離30 mmの計測点における自動ゲイン設定結果１３９、下側にSD距離15 mmの計測点における自動ゲイン設定結果１４０が示されている。表示法は凡例１３５に示す通りである。計測位置のセルを黒く塗りつぶした表示とした、検出光量が強いことを示す表示１３６、計測位置のセルを灰色に塗りつぶし、中に丸（○）表示をした、検出光量が適正であることを示す表示１３７、計測位置のセルを白く塗りつぶした表示とした、検出光量が弱いことを示す表示１３８を用いた。検出光量の結果は、プローブの装着状態に大きく依存するため、一部の計測点で検出光量が小さくなったときには光減衰量調整手段１４を用いることで、改善が可能となる。その場合には、光減衰量を変更した後に、ゲイン調整の再試行ボタン１３４を押下することで、再度検出器ゲインを調整することが可能となる。
本発明により、医療用および研究用機器として、あるいは教育効果・リハビリ効果の確認、家庭における健康管理、商品モニタ等の市場調査、さらに、同様の手法により組織酸素飽和度計測や筋肉の酸素代謝計測に用いることができる生体計測装置において、受光感度を各被検者・各部位で安定化することに利用できる。特に表層および深層成分を分離する等のための、複数の照射―検出器間距離（ＳＤ距離）を有するプローブを用いた装置において、計測精度の向上の手段として利用できる。
１０：被検者
１２：照射点
１３：検出点
１４：光減衰量調整手段
２０：装置本体
３０：光
４０：導波路
５０：送光器
６０：受光器
１０１：光源
１０２：光検出器
１０３：光源駆動部
１０４：増幅器
１０５：アナログ−デジタル変換器
１０６：制御・解析部
１０７：入力部
１０８：記憶部
１０９：表示部
１１３：OKボタン
１１４：キャンセルボタン
１３４：ゲイン調整の再試行ボタン
１３５：凡例
１３６：検出光量が強いことを示す表示
１３７：検出光量が適正であることを示す表示
１３８：検出光量が弱いことを示す表示
１３９：SD距離30 mmの計測点における自動ゲイン設定結果
１４０：SD距離15 mmの計測点における自動ゲイン設定結果
２０１：電気信号
２０２：光源および検出器内臓回路
２０３：調整量
２０４：筺体
２０５：手動調整つまみ
２０６：半導体レーザ（Laser Diode、 LD）
２０７：フォトダイオード（Photodiode、 PD）
２０８：色による表示
２０９：中間コネクタ
２１０：光ファイバ挿入口（減衰量が小さい）
２１１：光ファイバ挿入口（減衰量が中程度）
２１２：光ファイバ挿入口（減衰量が大きい）
２１３：挿入深さ調整部
２１４：挿入口
２１５：吸収または散乱媒体
２１６：導波孔（減衰量が大きい）
２１７：導波孔（減衰量が中程度）
２１８：導波孔（減衰量が小さい）
４０３：光ファイバ
５０１：送光用プローブ
５０２：受光用プローブ
５０３：プローブホルダ
５０４：送光用プローブと受光用プローブの一体型プローブ
５０５：送光用プローブ（減衰：大）
５０６：受光用プローブ（減衰：大）
５０７：SD距離30 mmの計測点
５０８：SD距離15 mmの計測点
５０９：SD距離5 mmの計測点
５１０：SD距離16.8 mmの計測点
６０１：実測の成分重み
６０２：実測データから算出される重み直線のＸ 切片
６０３：灰白質到達最小SD距離
６０４：灰白質平均実効光路長
６０５：不使用SD距離帯
６０６：SD距離１２ｍｍ以上の計測点における重みにフィッティングした１次直線
６０７：傾き０の直線
１２：照射点
１３：検出点
１４：光減衰量調整手段
２０：装置本体
３０：光
４０：導波路
５０：送光器
６０：受光器
１０１：光源
１０２：光検出器
１０３：光源駆動部
１０４：増幅器
１０５：アナログ−デジタル変換器
１０６：制御・解析部
１０７：入力部
１０８：記憶部
１０９：表示部
１１３：OKボタン
１１４：キャンセルボタン
１３４：ゲイン調整の再試行ボタン
１３５：凡例
１３６：検出光量が強いことを示す表示
１３７：検出光量が適正であることを示す表示
１３８：検出光量が弱いことを示す表示
１３９：SD距離30 mmの計測点における自動ゲイン設定結果
１４０：SD距離15 mmの計測点における自動ゲイン設定結果
２０１：電気信号
２０２：光源および検出器内臓回路
２０３：調整量
２０４：筺体
２０５：手動調整つまみ
２０６：半導体レーザ（Laser Diode、 LD）
２０７：フォトダイオード（Photodiode、 PD）
２０８：色による表示
２０９：中間コネクタ
２１０：光ファイバ挿入口（減衰量が小さい）
２１１：光ファイバ挿入口（減衰量が中程度）
２１２：光ファイバ挿入口（減衰量が大きい）
２１３：挿入深さ調整部
２１４：挿入口
２１５：吸収または散乱媒体
２１６：導波孔（減衰量が大きい）
２１７：導波孔（減衰量が中程度）
２１８：導波孔（減衰量が小さい）
４０３：光ファイバ
５０１：送光用プローブ
５０２：受光用プローブ
５０３：プローブホルダ
５０４：送光用プローブと受光用プローブの一体型プローブ
５０５：送光用プローブ（減衰：大）
５０６：受光用プローブ（減衰：大）
５０７：SD距離30 mmの計測点
５０８：SD距離15 mmの計測点
５０９：SD距離5 mmの計測点
５１０：SD距離16.8 mmの計測点
６０１：実測の成分重み
６０２：実測データから算出される重み直線のＸ 切片
６０３：灰白質到達最小SD距離
６０４：灰白質平均実効光路長
６０５：不使用SD距離帯
６０６：SD距離１２ｍｍ以上の計測点における重みにフィッティングした１次直線
６０７：傾き０の直線
Claims (11)
- 被検体に光を照射するための１つまたは複数の光照射手段と、
前記光照射手段から前記被検体上の照射点に照射され、前記被検体内を伝播してきた光を前記被検体上の検出点において検出するための１つまたは複数の光検出手段と、
前記光照射手段に含まれる光源素子と前記被検体との間、もしくは前記光検出手段に含まれる受光素子と前記被検体との間の光の伝播路上に設置される光減衰量調整手段と、
前記光検出手段で得られる信号を解析するための解析部と、
表層及び深層由来信号の分離性能評価手段と、
前記解析部での解析結果を表示するための表示部とを有し、
前記光照射手段と前記光検出手段の各々は、前記被検体上における、前記照射点と前記検出点間の距離として定義されるＳＤ距離が２種以上となるように前記被検体上に配置され、
前記解析部は、少なくとも１つの前記光検出手段が検出する２種以上のＳＤ距離の位置にある前記光照射手段からの光による信号を各々解析し、各々の受光量を所定の範囲内とするための光減衰調整量を算出し、
前記解析部は、２種以上のＳＤ距離における計測データから前記被検体の表層および深層における光吸収変化もしくは生体信号を分離する機能を備え、
前記光減衰量調整手段は、前記光検出手段が検出する２種以上のSD距離の位置にある前記光照射手段からの光の減衰量を、前記光検出手段に入射する光量を変化させることにより、各々調整可能としたことを特徴とする生体光計測装置。 - 請求項１に記載の生体光計測装置において、
前記表示部は、前記解析部で算出した前記光減衰調整量を表示することを特徴とする生体光計測装置。 - 請求項１に記載の生体光計測装置において、
前記分離性能評価手段は、ＳＤ距離8 mm以下の計測点を構成する光照射手段及び光検出手段、もしくは当該計測点で得られたデータであることを特徴とする生体光計測装置。 - 請求項２に記載の生体光計測装置において、
前記表示部は、前記光減衰調整量を、色の違いとして表示することを特徴とする生体光計測装置。 - 請求項２または３に記載の生体光計測装置において、
前記表示部は、前記光減衰量調整手段と同一の筺体に配置されることを特徴とする生体光計測装置。 - 請求項１〜請求項３の何れか１つに記載の生体光計測装置において、
前記光減衰量調整手段は、前記ＳＤ距離もしくは前記解析部で算出した前記光減衰調整量をもとに、前記１つまたは複数の光検出手段における検出光量を自動調整することを特徴とする生体光計測装置。 - 請求項１〜請求項６の何れか１つに記載の生体光計測装置において、
前記解析部は、前記光照射手段と前記光検出手段との組み合わせにより計測された複数の計測データから信号分離手法を用いて１つまたは複数の分離成分を抽出することを特徴とする生体光計測装置。 - 請求項１〜請求項５の何れか１つに記載の生体光計測装置において、
所定のＳＤ距離を優先ＳＤ距離として定義し、１つの光照射手段と１つの光検出手段により計測する領域もしくは領域を代表する位置を当該光照射手段と光検出手段に対応する計測点と定義したとき、
前記優先SD距離に対応する計測点（優先計測点）から４０ｍｍ以内の距離、より好ましくは２０ｍｍ以内の距離に、他のＳＤ距離に対応する計測点（従属計測点）が位置するよう前記光照射手段と前記光検出手段が配置され、
前記解析部は、前記優先計測点と、その優先計測点に対応する前記従属計測点をまとめて１つのデータセットとして解析し、
前記表示部は、１つまたは複数のデータセットの解析結果を、データセット毎に表示することを特徴とする生体光計測装置。 - 請求項８に記載の生体光計測装置において、
前記光照射手段と前記光検出手段の位置関係を入力するための入力手段を有し、
前記解析部は、前記光照射手段と前記光検出手段の位置関係から、前記優先計測点と前記従属計測点の組み合わせを算出し、
前記表示部は、前記組み合わせを図表により表示することを特徴とする生体計測装置。 - 請求項１〜請求項６の何れか１つに記載の生体光計測装置において、
前記解析部は、３種以上のＳＤ距離の条件において検出された光量およびその時間変化、もしくは信号分離手法により得られた各成分の、各SD距離における重みおよび振幅から、前記被検体の頭部構造を推定することを特徴とする生体光計測装置。 - 請求項１に記載の生体光計測装置を用いる生体光計測方法であって、
ＳＤ距離が１０ｍｍ未満、好ましくは８ｍｍ以下の少なくとも１つの計測点、およびＳＤ距離が１０ｍｍ以上、好ましくは１２ｍｍ以上の２つ以上の計測点において被検体内を伝播してきた光を計測するステップと、
ＳＤ距離と信号分離手法により得られた各成分の重み値とをプロットする図において、ＳＤ距離が１０ｍｍ以上の２つ以上の計測点の重み値を用いて１次直線を求め、重みゼロに対応するＳＤ距離を求めるステップと、
ＳＤ距離が１０ｍｍ未満の計測点の重み値を用いて軸に平行な直線を求めるステップと、前記軸に平行な直線と前記１次直線との交点のＳＤ距離を、灰白質到達最小ＳＤ距離として求めるステップと、
前記重みゼロに対応するＳＤ距離と前記灰白質到達最小ＳＤ距離から、脳寄与率を算出するステップと
を含む生体光計測装置を用いる生体光計測方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012191039 | 2012-08-31 | ||
| JP2012191039 | 2012-08-31 | ||
| PCT/JP2013/068788 WO2014034285A1 (ja) | 2012-08-31 | 2013-07-09 | 生体光計測装置およびそれを用いる生体光計測方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP5944511B2 true JP5944511B2 (ja) | 2016-07-05 |
| JPWO2014034285A1 JPWO2014034285A1 (ja) | 2016-08-08 |
Family
ID=50183106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2014532866A Active JP5944511B2 (ja) | 2012-08-31 | 2013-07-09 | 生体光計測装置およびそれを用いる生体光計測方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20150223694A1 (ja) |
| EP (1) | EP2891457A4 (ja) |
| JP (1) | JP5944511B2 (ja) |
| CN (1) | CN104582588B (ja) |
| WO (1) | WO2014034285A1 (ja) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015092872A1 (ja) * | 2013-12-18 | 2015-06-25 | 株式会社日立製作所 | 生体光計測装置およびそれを用いた生体光計測方法 |
| US20170340252A1 (en) * | 2014-12-22 | 2017-11-30 | Hitachi, Ltd. | Device for optical measurement of living body, analysis device, and analysis method |
| EP3243444A4 (en) * | 2015-01-07 | 2018-09-19 | National Institute of Advanced Industrial Science and Technology | Brain function measurement device and brain function measurement method |
| JP6609737B2 (ja) * | 2015-09-25 | 2019-11-27 | 株式会社ＮｅＵ | 生体光計測装置、情報処理プログラム、情報処理方法 |
| EP3184043A1 (en) * | 2015-12-22 | 2017-06-28 | IMEC vzw | Sensor, system and holder arrangement for biosignal activity measurement |
| CN105615895A (zh) * | 2015-12-25 | 2016-06-01 | 广州视源电子科技股份有限公司 | 血氧饱和度测量仪的调光方法和系统、测量仪及测量方法 |
| JP2019506205A (ja) * | 2015-12-31 | 2019-03-07 | ウェア２ビー リミテッド | 生理学的測定値の非侵襲的監視のための装置、システム、および方法 |
| JP6951348B2 (ja) * | 2016-02-12 | 2021-10-20 | ブライトシード・エルエルシーＢｒｉｔｅｓｅｅｄ，Ｌｌｃ | 組織特性を決定するためのシステム及び方法 |
| US10791981B2 (en) * | 2016-06-06 | 2020-10-06 | S Square Detect Medical Devices | Neuro attack prevention system, method, and apparatus |
| CN108209867B (zh) * | 2016-12-15 | 2022-03-18 | 松下知识产权经营株式会社 | 摄像装置 |
| EP3360464A1 (en) * | 2017-02-10 | 2018-08-15 | Carag AG | Apparatus and method for measuring the blood oxygen saturation in a subject's tissue |
| EP3610796B1 (en) * | 2017-04-14 | 2024-09-25 | National Institute of Advanced Industrial Science and Technology | Brain function measurement device and method |
| US11911130B2 (en) * | 2019-02-27 | 2024-02-27 | National Institute Of Advanced Industrial Science And Technology | Brain function measurement device, brain function measurement method, and probe |
| US20230366669A1 (en) | 2020-09-25 | 2023-11-16 | National University Corporation Shizuoka University | Measurement Sensitivity Calculation Method, Measurement Sensitivity Calculation Device, Measurement Sensitivity Calculation Program, and Optical Measurement Device |
| EP4122379A1 (en) * | 2021-07-23 | 2023-01-25 | Newmanbrain, S.L. | A method to obtain a near-infrared spectroscopy cerebral signal |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003207443A (ja) * | 2002-01-11 | 2003-07-25 | Shimadzu Corp | 光計測装置 |
| JP2004102172A (ja) * | 2002-09-12 | 2004-04-02 | Sumitomo Electric Ind Ltd | 光デバイス及び多チャンネル光伝送装置 |
| JP2004309791A (ja) * | 2003-04-07 | 2004-11-04 | Sumitomo Electric Ind Ltd | 光可変減衰器 |
| JP2004333568A (ja) * | 2003-04-30 | 2004-11-25 | Fujikura Ltd | 可変光アッテネータ |
| WO2012005303A1 (ja) * | 2010-07-06 | 2012-01-12 | 株式会社日立メディコ | 生体光計測装置およびそれを用いた生体光計測方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09273961A (ja) * | 1996-04-05 | 1997-10-21 | Tohoku Electric Power Co Inc | 物理量多点計測システム |
| JP3796086B2 (ja) * | 1999-12-27 | 2006-07-12 | 株式会社日立製作所 | 生体光計測装置 |
| US20030045798A1 (en) * | 2001-09-04 | 2003-03-06 | Richard Hular | Multisensor probe for tissue identification |
| US6628882B2 (en) * | 2001-09-13 | 2003-09-30 | Megasense, Inc. | VOA device and attenuation method with improved linearity |
| JP4327738B2 (ja) | 2005-01-18 | 2009-09-09 | 株式会社東芝 | 生体光計測装置及び生体光計測方法 |
| WO2007060583A2 (en) * | 2005-11-28 | 2007-05-31 | Koninklijke Philips Electronics N.V. | Method and apparatus for determining concentrations of analytes in a turbid medium |
| JP5324999B2 (ja) | 2009-04-09 | 2013-10-23 | 株式会社日立製作所 | 生体光計測装置及び生体光計測方法 |
| WO2010150751A1 (ja) | 2009-06-24 | 2010-12-29 | 株式会社日立製作所 | 生体計測装置 |
| JP5528739B2 (ja) * | 2009-08-12 | 2014-06-25 | 株式会社ジャパンディスプレイ | 検出装置、表示装置、および物体の近接距離測定方法 |
| JP5662700B2 (ja) * | 2010-05-17 | 2015-02-04 | 株式会社日立製作所 | 生体光計測装置及び生体光計測方法 |
-
2013
- 2013-07-09 US US14/424,572 patent/US20150223694A1/en not_active Abandoned
- 2013-07-09 WO PCT/JP2013/068788 patent/WO2014034285A1/ja active Application Filing
- 2013-07-09 JP JP2014532866A patent/JP5944511B2/ja active Active
- 2013-07-09 CN CN201380045180.9A patent/CN104582588B/zh active Active
- 2013-07-09 EP EP13833379.4A patent/EP2891457A4/en not_active Ceased
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003207443A (ja) * | 2002-01-11 | 2003-07-25 | Shimadzu Corp | 光計測装置 |
| JP2004102172A (ja) * | 2002-09-12 | 2004-04-02 | Sumitomo Electric Ind Ltd | 光デバイス及び多チャンネル光伝送装置 |
| JP2004309791A (ja) * | 2003-04-07 | 2004-11-04 | Sumitomo Electric Ind Ltd | 光可変減衰器 |
| JP2004333568A (ja) * | 2003-04-30 | 2004-11-25 | Fujikura Ltd | 可変光アッテネータ |
| WO2012005303A1 (ja) * | 2010-07-06 | 2012-01-12 | 株式会社日立メディコ | 生体光計測装置およびそれを用いた生体光計測方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2891457A4 (en) | 2016-04-20 |
| WO2014034285A1 (ja) | 2014-03-06 |
| JPWO2014034285A1 (ja) | 2016-08-08 |
| EP2891457A1 (en) | 2015-07-08 |
| CN104582588A (zh) | 2015-04-29 |
| US20150223694A1 (en) | 2015-08-13 |
| CN104582588B (zh) | 2017-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5944511B2 (ja) | 生体光計測装置およびそれを用いる生体光計測方法 | |
| JP5895025B2 (ja) | 生体光計測装置 | |
| JP5658993B2 (ja) | 生体計測装置 | |
| US10912504B2 (en) | Near-infrared spectroscopy and diffuse correlation spectroscopy device and methods | |
| EP2389862B1 (en) | Optical probe for optical imaging system | |
| CN106132299B (zh) | 测量装置、测量方法、程序和记录媒介 | |
| JP6285576B2 (ja) | 生体光計測装置 | |
| US11872022B2 (en) | System and method for an optical blood flow measurement | |
| JP4136704B2 (ja) | 光計測装置用のプローブとそれを用いたマルチチャンネル光計測装置 | |
| JP6211628B2 (ja) | 生体光計測装置およびそれを用いた生体光計測方法 | |
| JP6412956B2 (ja) | 生体光計測装置、解析装置、及び方法 | |
| Chi-Cheng | Studies of optical fiber biosensors | |
| Chen et al. | Simultaneously imaging of brain functions with Phased Array and PMS systems |
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: 20160426 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160525 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5944511 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |