Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4559417B2 - 生体光計測装置およびプログラム - Google Patents
[go: Go Back, main page]

JP4559417B2 - 生体光計測装置およびプログラム - Google Patents

生体光計測装置およびプログラム Download PDF

Info

Publication number
JP4559417B2
JP4559417B2 JP2006514547A JP2006514547A JP4559417B2 JP 4559417 B2 JP4559417 B2 JP 4559417B2 JP 2006514547 A JP2006514547 A JP 2006514547A JP 2006514547 A JP2006514547 A JP 2006514547A JP 4559417 B2 JP4559417 B2 JP 4559417B2
Authority
JP
Japan
Prior art keywords
measurement data
period
measurement
subject
concentration
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
Application number
JP2006514547A
Other languages
English (en)
Other versions
JPWO2005120349A1 (ja
Inventor
大樹 佐藤
敦 牧
剛 山本
雅史 木口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Publication of JPWO2005120349A1 publication Critical patent/JPWO2005120349A1/ja
Application granted granted Critical
Publication of JP4559417B2 publication Critical patent/JP4559417B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient; User input means
    • A61B5/742Details of notification to user or communication with user or patient; User input means using visual displays
    • A61B5/7425Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient; User input means
    • A61B5/742Details of notification to user or communication with user or patient; User input means using visual displays
    • A61B5/743Displaying an image simultaneously with additional graphical information, e.g. symbols, charts, function plots
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Artificial Intelligence (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • General Business, Economics & Management (AREA)
  • Epidemiology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

本発明は生体内部の情報、特に光吸収物質の濃度変化を、光によって計測する生体光計測装置に関し、特に生体光計測で計測したデータを用いて健康状態の診断を支援する情報を提供する生体光計測装置に関する。
生体に対する透過性が高い、可視から近赤外領域に光強度ピーク波長を持つ光を用いると、生体内部の情報を無侵襲に計測することが可能である。これは、計測される光信号の対数値が光路長と濃度の積に比例することを示したLambert-Beer則に基づく。この法則を発展させ、生体中の「酸素化ヘモグロビン(Hb)」と「脱酸素化Hb」、および「総Hb(酸素化Hbと脱酸素化Hbの総和)」の相対的濃度変化(以下、濃度変化と呼ぶ)を計測する技術が開発されてきた。
特に、この技術を用いて人間の大脳皮質におけるHb濃度変化を無侵襲に多点同時計測する技術が提案され(非特許文献1)、研究および臨床面において広まりつつある。上記文献では、大脳皮質のHb濃度変化を計測することにより、人の脳機能を計測する方法が開示されている。具体的には、人の知覚機能や運動機能の賦活に伴い、その機能を司る大脳皮質領野の血液量が局所的に増加し、該当部位の酸素化Hbや脱酸素化Hbの濃度が変化するため、脳の活動状況が評価できる。
この技術は、被験者に対し無侵襲且つ低拘束で、簡便に脳機能を計測できる特長を有し、また、血行動態や血液循環の状態などを評価できるため、従来にない医療診断への応用(臨床応用)が期待されている。これまでの臨床応用の例としては、てんかん発作の診断(非特許文献2)や言語有意半球の同定(非特許文献3)における有用性が報告されている。
最近では、精神医学への応用が進められており、うつ病や統合失調症などの精神疾患患者において、生体光計測により得られたHb濃度変化が健常者とは異なるパターンを示すことが報告された(非特許文献4)。ここでは、酸素化Hbを指標として用い、各疾病を持つ被験者に特有な波形が示されている。しかし、各個人内における再現性は明らかではなく、疾病の有無あるいは病状を各被験者において明確に診断することは容易ではない。また、各Hb濃度変化の個人差が、課題へ取り組む姿勢(心理的要因)によって生じるものか、また血行動態の違い(生理的要因)によるものか、についても解明されていない。
Maki, A et al., (1995) "Spatial and temporal analysis of human motor activity using noninvasive NIR topography・Medical Physics 22, 1997-2005. Watanabe, E et al., (1998) "Non-invasive cerebral blood volume measurement during seizures using multi-channel near-infrared spectroscopic topography・J. epilepsy 11, 335-340. Watanabe, E et al., (1998) "Non-invasive assessment of language dominance with near-infrared spectroscopic mapping・Neurosci. Lett. 256, 49-52. Suto, T et al., (2004) "Multichannel near-infrared spectroscopy in depression and schizophrenia: cognitive brain activation study・Biol. Psychiatry 55(5),501-511.
従来、脳活動信号を用いた診断支援では、複数被験者の平均値や分散値と比較することが主流であった。しかし、その場合、個人差を吸収することが難しく、診断支援技術としての精度が低かった。この診断支援技術の精度を高めることが、脳活動信号を用いた診断支援技術において重要な課題である。
本発明者は、31名の健常者を対象に、指の運動課題を用いて計測を行った。単純な運動課題による運動野の活動は、より原初的な活動例であるため、個人差の主たる原因となる心理的要因の影響は少ない。4名の被験者の結果(酸素化Hb濃度変化)を図3に示す。課題期間に伴い増加する傾向は共通するものの、その波形は被験者により大きく異なることが分かる。また、同様に脱酸素化Hb濃度変化や総Hb濃度変化の波形にも被験者間のばらつきが見られる。このように、比較的個人差が小さいと考えられる単純な運動野活動の例をとっても、健常者群においてさえ、Hb濃度変化のパターンは被験者間で異なる。
被験者群を全体として評価した場合、ある程度の特徴を見出せる可能性はあるが、診断とは各個人においてなされるものであり、集団における統計的な類型化とは異なる。現状の光脳機能計測データに見られるHb濃度変化の個人差を考慮すると、集団データから抽出された典型波形(テンプレート)とのマッチングを用いた個人診断は、精度が充分でないと考えられる。
一般に、脳活動に伴い、酸素化Hbの増加、脱酸素化Hbの減少、および総Hbの増加が計測される場合が多いが、我々の計測結果では、図4に示すような各Hb濃度変化の多様性と、変化傾向(増減)の個人差が明らかになった。表1に、各Hb濃度変化の増減によるパターン分類結果を示す。酸素化Hb、脱酸素化Hb、総Hbの各Hbに対し、増加した場合を「↑」、減少した場合を「↓」、有意な変化がなかった場合を「―」として表し、その組み合わせに従い、変化パターンとして分類した。
Figure 0004559417
ここでは、3種類のHb濃度(酸素化Hb,脱酸素化Hb、総酸素化Hb)を3つの変化(増加↑、減少↓、変化なし―)に分類しているため、27通りのパターンが可能であるが、便宜的に出現頻度の高い8パターンに1〜8の数字を割り当てて示した。右端の列が、被験者31名中、各変化パターンを示した割合(%)を示す。これまでの知見と一致する変化(運動野活動に伴う、酸素化Hb濃度の増加、脱酸素化Hb濃度の減少、および総Hb濃度の増加)が見られるケースが多いが、各Hb濃度は独立に増加や減少、あるいは無変化を示すため、被験者によって各Hb濃度変化の組み合わせパターン(以下、Hb変化パターンと略す)は異なる。
各Hb濃度変化信号の再現性を検討するため、7人の被験者を対象に同じ計測を約6ヶ月の期間を空けて繰り返した。その結果、上記Hb変化パターンはほぼ再現する。また、各Hb濃度変化量(信号強度)の再現性は低いものの(図5)、変化領域(活動領域)の中心部位(活性中心)は約2cm以内に再現する(表2)。生体光計測の空間分解能は2〜3cmであることを考慮すると、空間情報の再現性は許容範囲内である。更に、活動領域における各Hb濃度の時間変化(活動波形)を比較すると、非常に高い相関を示す(図6、表3)。これは、活動波形が生理状態の個人特性を表す指標として有効であることを示すものである。
Figure 0004559417
Figure 0004559417
以上の結果を考慮すると、各被験者から得られたデータを、集団から得られた典型的な変化パターン(例:変化量、波形)と比較する方法より、同一被験者に対して同じ計測方法で繰り返し計測を行い、個人内の時間経過に伴う変化を検討する方法の方が、診断に有用な情報を提供し得る。特に、個人内における脳活動に伴う血行動態変化の再現性を調べるためには、(1)複数のHb濃度の増減を組み合わせたHb変化パターンの分類、(2)各Hb濃度変化の活動波形(時間情報)の比較、が効果的である。
また、(3)各Hb濃度変化の活性中心(空間情報)の比較も有用である。なお、ここでは、酸素化ヘモグロビン、脱酸素化ヘモグロビン、総ヘモグロビンの3つの計測データを用いているが、酸素化ヘモグロビンおよび脱酸素化ヘモグロビンの2つの増減に関するデータを用いても診断に有用な情報を提供し得る。
本発明の生体光計測装置は、生体光計測装置で得られたデータにおける時間経過に伴う個人内の変化に注目し、また、その特性を利用し、各個人内における状態変化を示す有効な情報を作成および表示することによって、各個人に対する効果的な診断支援を行う。本生体光計測装置では、計測した光信号からHb濃度変化を算出・記録する手段と、異なる時間に計測した複数の記録データを同時に読み出し、表示する手段を有する。また、各データ(各Hb濃度変化)から、2つ以上のHb濃度変化(例:酸素化Hb濃度変化と脱酸素化Hb濃度変化)を組み合わせて、1つのHb変化パターン(例:パターン1=酸素化Hb濃度の増加、脱酸素化Hb濃度の減少、および総Hb濃度の増加)として分類し、記号、画像、あるいはグラフなどで表示する手段と、1つのHb濃度変化、あるいは2つ以上のHb濃度変化から算出した指標(例:総Hb=酸素化Hb+脱酸素化Hb)に対し、計測時刻の異なる2つ以上のデータを、画像、グラフもしくは数値として、同時に表示する手段を備え、圧縮・簡素化した情報を検査者に与える。
本発明によれば、複数の異なるデータの情報を統合し新たな情報を与えるため、検査者の処理が軽くなり診断が容易または的確になる。
また、同一被験者に対して異なる時刻に計測した複数のデータセットを統合し、その情報を分かりやすく検査者に表示できるため、各被験者における状態変化の容易な診断支援を可能とする。
また個人内の経時変化を効果的に把握できるため、リハビリテーションやイメージトレーニング効果の評価・診断などへも応用できる。
以下、本発明の実施形態について図面を参照して説明する。図1は、本発明による生体光計測装置の概要を示すブロック図である。生体光計測部は、パーソナルコンピュータやワークステーションに代表される電子計算機から構成される制御装置103と、異なる波長にピーク波長を持つ2つのレーザダイオード109と110と、前記2つのレーザダイオードを異なった周波数で変調するための信号を生成する発振器107と108と、前記ピーク波長の異なる2つの光を混合する光混合器111と、前記光混合器111からの光を光ファイバ経由で被検体上の光照射位置に照射する光照射手段と、前記光照射手段から適度に離れた光検出位置(本実施例では約3cm離れた点)から混合光を検出する光検出器112と、前記発振器からの変調周波数が参照信号として入力されたロックインアンプ105および106と、ロックインアンプの出力である各波長帯の光の透過光信号をアナログ信号からデジタル信号へ変換するアナログ−デジタル変換器104を備える。
ここで、光源と検出器の略中点の位置がもっとも多くの代謝物質情報を含むため、便宜的にこの略中点位置を計測点とするが、実質的には光源と検出器の対により計測することのできる感度分布にもとづき計測点(計測領域)は決められる。図1では、計測位置を1つのみ記したが、発振器を用いて複数の光信号を分離することが可能であるため、実際には図2のように複数の光照射位置と光検出位置を配列した生体光計測装置を用いて、複数の位置を計測することが出来る。
本実施例では発振器を用いて複数の光信号を分離しているが、発振器は使わずにパルス光を用いて点灯タイミングで光信号を分離することも可能である。各波長帯の光の透過光信号はアナログ−デジタル変換器104でアナログ−デジタル変換された後、制御装置103に入力・記憶される。制御装置103では透過光信号を元に各計測部位におけるHb濃度変化が算出され、元信号(透過光信号)と共にデータ蓄積サーバ102に転送され記憶される。上記計測データは、データ蓄積サーバ以外にも、ICカードなどの記憶媒体に記憶させ、各個人が携帯することが望ましい。
診断支援部101の記憶部では、103で計測された計測データと、データ蓄積サーバ102もしくはICカードなどの記憶媒体に記憶された同一被験者の過去の計測データの両方を読み込み、記憶する。そして診断支援部101の解析部では、これら複数のデータから、それぞれの活動(各Hb濃度変化)の空間情報(活動マップ)、時間情報(活動波形)を算出し、記号、数値、画像、グラフなどを用いてディスプレイ113上に表示する。
図7および図8には、脳活動に伴うHb濃度変化を多点同時計測した場合の表示内容の例を示す。また、図9には表示までのデータ処理フローを示す。なお、このデータ処理フローに示されたデータ処理は、従来の生体光計測装置にこれらの処理を実行させるためのプログラムをインストール等することによってもできる。
最初の処理では、まずデータの比較形式を選択する。直前に計測したデータと過去の任意のデータを比較する場合の表示画面を図7に示し(図7方式)、過去のデータから任意に選択し比較する場合の表示画面を図8に示した(図8方式)。ここでは、比較するデータセットを2つに設定しているが、同じ方式で2つ以上のデータ比較も可能である。ステップ901のデータの選択・入力処理では、図7プルダウンメニュー701のようなインターフェースにより、今回計測したデータと比較するデータを、データ蓄積サーバあるいはICカードのような記憶媒体に保存されている過去のデータ群から選択し読み込むことが出来る。
同様に、ステップ902では、図8のプルダウンメニュー801、802に示したように、2つのデータをデータ蓄積サーバもしくはICカードのような記憶媒体からそれぞれ選択し読み込むことが可能である。前記プルダウンメニューで「平均」を選択すると、同一被験者の複数データの平均が算出され、その被験者の一般傾向(過去のデータ全体)と比較することが出来るため、有用である。過去のデータ全体から複数のデータを任意に選択し平均することも可能である。選択され診断装置101に入力するデータ形式は任意だが、元信号である透過光強度データの場合、ステップ903でHb濃度変化を算出する。
その後、ステップ904では安静期間と活動期間を設定する。設定方法は、ユーザが任意に設定する以外に、自動的に信号波形から最適な活動期間を設定することも可能である。例えば、活動期間になり得る全ての範囲と活動期間の長さだけを設定し、平均Hb濃度変化量(絶対値)が最大になる期間を算出し設定する。安静期間および活動期間のパラメータを用いて、ステップ905では活動期間における酸素化Hb、脱酸素化Hb、総Hbの濃度変化が有意であるかどうか検定する。
例えば、課題が繰り返し複数回行われた場合、各回における安静期間の平均Hb濃度変化量(ベースライン)と、課題(活動)期間の平均Hb濃度変化量(信号強度)を算出し、各計測点において複数のベースライン期間と活動期間における平均値の差を、t検定や分散分析などの統計手法によって評価する。また、活動期間に予測されるHb濃度の時間変化との相関係数や、ノンパラメトリックな統計手法によっても、有意なHb濃度変化(活動)の有無を評価することが可能であり、ここでは統計手法を限定する必要はない。また、統計手法によっては安静期間と活動期間の設定を必要としない。活動(有意性)を検定した後、目的とする表示内容に従い処理を進めていく。
本実施形態の表示例では、比較する両データのHb変化パターン(図7、702;図8、803)、活動波形(図7、704;図8、805)、および活動マップ(図7、703;図8、804)を示した。これらは各期間における計測データを同期して表示したり、平均値を表示したりすることができる。それぞれの詳細を以下に記す。
表示内容がHb変化パターンあるいは活動波形である場合、はじめに、各計測において代表的な変化傾向(各Hb濃度の増減)を決定する必要がある。なぜならば、本実施例では多点同時計測を行っており、全計測点が同じ変化傾向を示すことは稀であるからである。代表的活動の選択方法(評価に用いる変化傾向)の決定方法には数通りあり、例えば、ステップ906のように各Hbで最も強い信号強度(統計値、平均変化量、最大変化量など)を示した計測点の変化を採用する方法がある。
あるいは、ステップ907のように全計測点の中で最も多くの計測点が示した変化(増加か減少)を採用し、どの計測点でも有意な変化がなかった場合に「変化なし」という分類を採用する、という方法も可能である。また、予め検討目的とする部位が決まっている場合、その部位に対応した計測点をそのまま選択し、その変化を採用する方法もある(ステップ908)。
その後、Hb変化パターンを表示する場合は、ステップ911で各Hb濃度変化の組み合わせによるHb変化パターンの分類を行う。例えば、課題期間における各Hb濃度変化の平均値が安静期間より有意に小さい場合を減少(↓)、有意に大きい場合を増加(↑)、有意差がない場合を変化なし(−)として、各Hb濃度変化の増減の組み合わせを記号で表示する(図7の702、図8の803)。本データと比較データ、あるいはデータ1とデータ2のHb変化パターンを並べて表示することにより、基本的なHb変化パターンの被験者内経時変化の有無が判定できる。
診断例として、従来、酸素化Hbの増加と脱酸素化Hbの減少を示していた被験者が、酸素化Hbと脱酸素化Hb両方の増加を示した場合を考える。この場合、今までは血液中の酸素消費の速さ以上に血流速が速くなり脱酸素化Hbが減少していたが、今回の結果では血流速の上昇が少なく、酸素消費による脱酸素化の増加が観察されたと考えられる。血流速が遅くなったという観点から、血液の流れが悪くなったり血管の弾力性が低下したりという症状が考えられ、被験者を診断する材料となり得る。
また、活動波形を表示する場合も、上記ステップ906、907、908で選択した変化を示す。本実施例の表示画面(図7:704、図8:805)では、太い線で今回の計測結果あるいはデータ1を、細い線で過去に計測した比較データあるいはデータ2を示した。縦軸が、Hb濃度変化量(mM・mm)、横軸が時間(秒)を表す。表示には、選択した1つの計測点における活動波形を用いてもよいし、複数計測点のデータを平均した活動波形を用いてもよい(ステップ913)。
複数データの比較を容易にするためには、例えば、図10(a)のように過去の複数データおよび今回のデータを、重ねてプロットし表示する方法がある。この方法により、今回のデータがこれまでの結果と一致するものかどうか容易に判断することが出来る。また、複数データから平均波形と標準偏差を算出して表示し(図10(b))、検討する今回のデータを重ねてプロットすることにより、これまでの傾向からどの程度離れているか、診断に有用な情報を与える。更にステップ914では、活動波形の比較結果を具体的に数値などで表示する。
本実施形態の表示例では波形間の類似性として相関係数をグラフ右に示した。相関係数以外にも、波形の特徴量(潜時、立ち上がり時間、角度、振幅など)を数値化して示すことが可能である。また、これまでに得られているデータの平均波形(および分散)と今回のデータ波形を統計的に比較し、波形の再現度を数値化して表示することも有効である。これらの活動波形の再現性評価に関して、健常状態における範囲を自動あるいは手動で設定しておき、その範囲から逸脱した場合、何らかの警告を発する機能が有効である。
これらの波形類似性の比較は、前記のHb活動パターンを更に詳細に比較したものであり、同じHb活動パターンを示した場合でも、波形パターンの詳細な違いにより、血行動態の変化を見出すことが出来る。極端な例としては、図11に示したようにこれまでの平均範囲に酸素化Hbの活動波形が示されていた状態(図11(a))から、平均範囲から大きく外れた波形を描く状態(図11(b))へ変化する場合がある。この場合、脳活動に伴う血液量の増加が弱まっているため、血管の柔軟性の低下や、脳活動に伴う生理変化の変容(例:血液量増加が主→酸素消費が主)の可能性がある。
更に、この同一被験者内の活動波形の類似性を利用して、信号処理方法を改善することも可能である。図12に示したように各被験者のHb濃度変化の波形には個人差が生じる。しかし、その個人差には、脳活動パターンの個人差だけでなく装置系雑音や他の生理雑音も含まれているため、適切な処理方法で脳活動由来の変化を抽出することが望ましい。具体的には、計測される元信号(透過光強度)から各Hb濃度変化を算出する時に設定する幾つかのパラメータを、活動波形の類似性を利用して最適化できる。例えば、Hb濃度変化の算出に必要なベースラインを、刺激前数秒間の値と刺激後数秒間の値からフィッティングカーブを求めて用いる場合、「刺激前の期間」「刺激後の期間」「フィッティングカーブの次数」を、各刺激期間のHb濃度変化の相関係数が最大となるよう設定する。また、任意の数値をベースラインとして各Hb濃度変化を算出後、周波数フィルタリングにより補正する場合には、各刺激期間のHb濃度変化の相関係数を基準にフィルターを作成することが出来る。
次に、活動マップについて説明する。本実施例で表示した活動マップ(図7、703;図8、804)では、左右の列に異なる日に計測した同一被験者の結果を示した(上から酸素化Hb、脱酸素化Hb、総Hb)。活動の強さを各計測点の信号強度(課題期間の平均Hb濃度変化量)とし、右側のカラーバーに従い図示した(ステップ909)。これは、活動期間における最大Hb濃度変化量などを示してもよい。また、前記Hb変化パターンの分類のように、課題に伴う有意な各Hb濃度変化を統計的に評価し、その統計値を示す方法や、統計的に有意な変化が認められた計測点のみを選択して、その信号強度を示す方法も有効である。ステップ910では、活動マップの比較結果を具体的に表示する。
例えば、活動マップ上に活動の中心部位(活動中心または活性中心)を表示すると、時間経過に伴う活性中心の変化が分かる。活性中心の決定方法としては、最も強い活動を示した計測点を選択する方法や、有意な変化を示した活動領域から信号強度や統計値を基準に重心を求める方法などがある。本実施例の表示画面では、信号強度を基準に求めた重心(活性中心)を「+」で表示し、各データの活性中心間の距離(cm)をマップの下に表示した。活動波形と同様に、活性中心の位置再現性に関して、健常状態における再現範囲を自動あるいは手動で設定しておき、その範囲から逸脱した場合、何らかの警告を発する機能が有効である。健常者の場合、この活性中心はほぼ同じ位置に再現するため(表2)、正しい部位を計測しているかどうかの確認、および位置補正にも使える。
また、脳梗塞患者など、ある機能が一時的に障害を受けた被験者の場合、リハビリテーションによって障害を受けた機能が他の脳部位で代償される場合があるので、その代償部位を観察することができる。また、イメージトレーニングを含む何らかの機能訓練中には、その機能の上達に伴う、活動部位の広がりが観察でき、訓練の効果をフィードバックできる。
個人内の経時変化を容易且つ効果的に把握できるため、リハビリテーションやイメージトレーニング効果の評価・診断などへ応用できる。
本発明の生体光計測装置の一実施形態を示すブロック図である。 本発明の生体光計測装置における実施形態の例(多点同時計測)を示す模式図である。 酸素化Hb濃度変化の波形における個人差の例を示す図である。 酸素化Hb、脱酸素化Hb、および総Hbの濃度変化パターンにおける個人差の例を示す図である。 酸素化Hb、脱酸素化Hb、および総Hb濃度における変化量(信号強度)の再現性を示す図である。 活動領域における酸素化Hb、脱酸素化Hb、および総Hb濃度の時間変化(活動波形)を示す図である。 本発明の生体光計測装置における表示画面の一例を示す図である。 本発明の生体光計測装置における表示画面の一例を示す図である。 診断支援部における主な処理を示したフローチャートを示す図である。 波形比較に関する表示画面の例を示す図であり、(a)は複数の比較データを重ねてプロットした図、(b)は複数データから平均波形と標準偏差を算出して表示した図である。 同一被験者内において酸素化Hb濃度の変化パターンが大きく異なる例((a)から(b)への変化)を示す図である。 装置系雑音や他の生理雑音を含んだ、Hb濃度変化波形における個人差を示す図である。
符号の説明
101…診断支援部、102…データ蓄積サーバ、103…制御装置、104…アナログデジタル変換器、105…ロックインアンプ、106…ロックインアンプ、107…発振器、108…発振器、109…λ1光源、110…λ2光源、111…光混合器、112…光検出器、113…ディスプレイ、701…比較データを選択するプルダウンメニュー、702…Hb変化パターンの表示部、703…活動マップの表示部、704…活動波形の表示部、801…データ1を選択するプルダウンメニュー、802…データ2を選択するプルダウンメニュー、803…Hb変化パターンの表示部、804…活動マップの表示部、805…活動波形の表示部。

Claims (14)

  1. 被検体に光を照射する光源と、
    前記光源から照射され前記被検体内を伝播した光を検出する検出器と、
    前記検出器で検出された光から前記被検体内部の酸素化ヘモグロビンおよび脱酸素化ヘモグロビンの濃度変化についての計測データを算出する演算部と、
    前記被検体の第一の期間における前記計測データと、前記第一の期間とは異なる第二の期間における前記計測データを記憶する記憶部を有し、
    前記演算部は前記計測データを、前記酸素化ヘモグロビンの濃度の増減、前記脱酸素化ヘモグロビンの濃度の増減ごとに分類し、
    前記第一の期間における前記分類と、前記第二の期間における前記分類を表示する表示部を有する、
    ことを特徴とする生体光計測装置。
  2. 前記第一の期間における前記計測データは前記記憶部に記憶された前記被検体の過去の計測データであり、前記第二の期間における前記計測データは前記被検体の直前に計測された計測データであることを特徴とする請求項1記載の生体光計測装置。
  3. 前記表示部は、前記過去の計測データを過去の複数の計測データのうちから選択するための選択部を有することを特徴とする請求項2記載の生体光計測装置。
  4. 前記記憶部は複数の期間における前記計測データを記憶し、
    前記表示部は前記複数の期間のうちどの期間における前記計測データを表示するかを選択する選択部を有し、選択部で選択された期間における計測データの前記分類を表示することを特徴とする請求項1記載の生体光計測装置。
  5. 前記演算部は前記複数の期間における前記計測データの平均値を算出し、
    前記表示部は前記平均値を表示することを特徴とする請求項4記載の生体光計測装置。
  6. 前記表示部は前記第一の期間における活動中心位置と、前記第二の期間における活動中心位置の距離を表示することを特徴とする請求項1記載の生体光計測装置。
  7. 前記演算部は前記第一の期間における前記計測データと、前記第二の期間における前記計測データの相関係数を算出し、
    前記表示部は前記相関係数を表示することを特徴とする請求項1記載の生体光計測装置。
  8. 前記第一および第二の期間は、前記被検体が安静にする安静期間と、前記被検体が活動する活動期間を有することを特徴とする請求項1記載の生体光計測装置。
  9. 前記表示部は、前記安静期間と前記活動期間の長さを設定する設定部有することを特徴とする請求項8記載の生体光計測装置。
  10. 前記演算部は、最も強い信号強度を示した計測点における前記計測データの増減、または最も多くの計測点が示した前記計測データの増減、または所定の計測点における前記計測データの増減のいずれか一の基準にもとづき、前記酸素化ヘモグロビン、前記脱酸素化ヘモグロビンの増減を分類することを特徴とする請求項1記載の生体光計測装置。
  11. 前記表示部は前記第一の期間における前記計測データと、前記第二の期間における前記計測データを同期して表示することを特徴とする請求項1記載の生体光計測装置。
  12. 前記計測データは、前記被検体内部の酸素化ヘモグロビン、脱酸素化ヘモグロビンおよび、前記酸素化ヘモグロビンおよび脱酸素化ヘモグロビンの総濃度の濃度変化についての計測データであり、
    前記演算部は前記計測データを、前記酸素化ヘモグロビンの濃度の増減、前記脱酸素化ヘモグロビンの濃度の増減、前記酸素化ヘモグロビンおよび脱酸素化ヘモグロビンの総濃度の増減ごとに分類することを特徴とする請求項1記載の生体光計測装置。
  13. 被検体に光を照射し、前記被検体内を伝播した光を検出することで前記被検体内部の情報を計測する生体光計測装置によって計測される計測データに適用されるプログラムであって、
    コンピュータに
    第一の期間における前記被検体の計測データを記憶する第一のステップと、
    第二の期間における前記被検体の計測データを記憶する第二のステップと、
    前記計測データを、前記酸素化ヘモグロビンの濃度の増減、前記脱酸素化ヘモグロビンの濃度の増減ごとに分類する第三のステップと、
    前記第一の期間における前記分類と、前記第二の期間における前記分類を表示する第四のステップ、
    を実行させることを特徴とするプログラム。
  14. コンピュータに
    前記第一のステップにおける前記計測データと前記第二のステップにおける前記計測データを同期して表示するステップを実行させることを特徴とする請求項13記載のプログラム。
JP2006514547A 2004-06-14 2005-06-08 生体光計測装置およびプログラム Expired - Fee Related JP4559417B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2004174925 2004-06-14
JP2004174925 2004-06-14
PCT/JP2005/010501 WO2005120349A1 (ja) 2004-06-14 2005-06-08 生体光計測装置、生体光計測方法およびプログラム

Publications (2)

Publication Number Publication Date
JPWO2005120349A1 JPWO2005120349A1 (ja) 2008-04-03
JP4559417B2 true JP4559417B2 (ja) 2010-10-06

Family

ID=35502774

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006514547A Expired - Fee Related JP4559417B2 (ja) 2004-06-14 2005-06-08 生体光計測装置およびプログラム

Country Status (4)

Country Link
US (2) US7516021B2 (ja)
EP (1) EP1762175B1 (ja)
JP (1) JP4559417B2 (ja)
WO (1) WO2005120349A1 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4702107B2 (ja) * 2006-03-03 2011-06-15 株式会社日立製作所 生体光計測装置
JP5244027B2 (ja) * 2008-05-12 2013-07-24 俊徳 加藤 脳機能解析支援装置及びプログラム
JP5386594B2 (ja) * 2009-11-27 2014-01-15 株式会社日立製作所 生体光計測装置
US20140206962A1 (en) * 2011-04-15 2014-07-24 Hitachi Medical Corporation Biophotonic measurement device, biophotonic measurement device operating method, and biophotonic measurement data analysis and display method
US9131889B2 (en) * 2011-05-31 2015-09-15 Nagoya Institute Of Technology Cognitive impairment determination apparatus, cognitive impairment determination system and program
JP5809864B2 (ja) * 2011-07-07 2015-11-11 株式会社日立メディコ 生体光計測装置、波形解析方法
JP5280494B2 (ja) * 2011-07-13 2013-09-04 株式会社日立製作所 複数脳賦活観測システム
AU2017355315B2 (en) 2016-11-07 2023-12-14 Climate Llc Work layer imaging and analysis for implement monitoring, control and operator feedback
JP6787069B2 (ja) * 2016-11-16 2020-11-18 株式会社デンソー 生体情報解析装置、及び生体情報解析方法
KR102444393B1 (ko) * 2020-09-25 2022-09-20 주식회사 올리브헬스케어 유방암 진단 시스템
WO2023188698A1 (ja) * 2022-03-31 2023-10-05 株式会社島津製作所 評価方法、評価用装置、およびプログラム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0998972A (ja) * 1995-10-06 1997-04-15 Hitachi Ltd 生体光計測装置及び画像作成方法
JP2002119511A (ja) * 2000-10-13 2002-04-23 Hitachi Ltd 意思表示装置及び意思表示方法。
WO2002032317A1 (en) * 2000-10-16 2002-04-25 Hitachi Medical Corporation Organism optical measurement instrument

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
USRE45607E1 (en) * 1998-10-13 2015-07-14 Covidien Lp Multi-channel non-invasive tissue oximeter
JP4167860B2 (ja) * 2002-07-08 2008-10-22 株式会社日立製作所 生体光計測装置
JP3772220B2 (ja) * 2002-11-01 2006-05-10 独立行政法人情報通信研究機構 習得度測定装置及び習得度測定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0998972A (ja) * 1995-10-06 1997-04-15 Hitachi Ltd 生体光計測装置及び画像作成方法
JP2002119511A (ja) * 2000-10-13 2002-04-23 Hitachi Ltd 意思表示装置及び意思表示方法。
WO2002032317A1 (en) * 2000-10-16 2002-04-25 Hitachi Medical Corporation Organism optical measurement instrument

Also Published As

Publication number Publication date
JPWO2005120349A1 (ja) 2008-04-03
WO2005120349A1 (ja) 2005-12-22
US20090163786A1 (en) 2009-06-25
EP1762175A4 (en) 2010-12-01
EP1762175B1 (en) 2012-03-21
EP1762175A1 (en) 2007-03-14
US20080040049A1 (en) 2008-02-14
US7890270B2 (en) 2011-02-15
US7516021B2 (en) 2009-04-07

Similar Documents

Publication Publication Date Title
Yücel et al. Best practices for fNIRS publications
US7890270B2 (en) Optical system for measuring metabolism in a body, method and program
Novi et al. Functional near-infrared spectroscopy for speech protocols: characterization of motion artifacts and guidelines for improving data analysis
Wang et al. HemaApp: noninvasive blood screening of hemoglobin using smartphone cameras
Lai et al. Functional near-infrared spectroscopy in psychiatry
US20190159675A1 (en) Point-of-care tele monitoring device for neurological disorders and neurovascular diseases and system and method thereof
CN100444793C (zh) 生物体光测量装置
JP5982483B2 (ja) 生体状態評価装置およびそのためのプログラム
Fantini et al. Perspective: Prospects of non-invasive sensing of the human brain with diffuse optical imaging
JP5319960B2 (ja) 生体光計測装置
US9167970B2 (en) Non-invasive optical imaging for measuring pulse and arterial elasticity in the brain
Dadgostar et al. Functional connectivity of the PFC via partial correlation
Kober et al. Hemodynamic signal changes during saliva and water swallowing: a near-infrared spectroscopy study
Ewerlöf et al. Estimation of skin microcirculatory hemoglobin oxygen saturation and red blood cell tissue fraction using a multispectral snapshot imaging system: a validation study
US10835169B2 (en) Brain function index computing device and brain function index computing method
JP4702107B2 (ja) 生体光計測装置
Karen et al. Cerebral hemodynamic responses in preterm-born neonates to visual stimulation: classification according to subgroups and analysis of frontotemporal–occipital functional connectivity
Roy et al. Exploring the impact and influence of melanin on frequency-domain near-infrared spectroscopy measurements
Al-Shargie Quantification of mental stress using fNIRS signals
Nambu et al. Transient increase in systemic interferences in the superficial layer and its influence on event-related motor tasks: a functional near-infrared spectroscopy study
YAMADA Continuous wave functional near-infrared spectroscopy: various signal components and appropriate management
Amyot et al. Cerebrovascular reactivity measurement with functional near infrared spectroscopy
Li et al. Correlation and asynchronization of electroencephalogram and cerebral blood flow in active and passive stimulations
Mieles Freire et al. A non-invasive portable solution to estimate hemoglobin levels in the blood
Kim et al. Frontal lobe hemodynamics detected by functional near-infrared spectroscopy during head-up tilt table tests in patients with electrical burns

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100420

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100616

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: 20100706

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100722

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130730

Year of fee payment: 3

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

LAPS Cancellation because of no payment of annual fees