JP5159893B2 - 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 - Google Patents
生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 Download PDFInfo
- Publication number
- JP5159893B2 JP5159893B2 JP2010531878A JP2010531878A JP5159893B2 JP 5159893 B2 JP5159893 B2 JP 5159893B2 JP 2010531878 A JP2010531878 A JP 2010531878A JP 2010531878 A JP2010531878 A JP 2010531878A JP 5159893 B2 JP5159893 B2 JP 5159893B2
- Authority
- JP
- Japan
- Prior art keywords
- dimensional
- probe
- light
- coordinate
- image
- 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 139
- 238000001514 detection method Methods 0.000 title claims description 110
- 238000000034 method Methods 0.000 title claims description 106
- 239000000523 sample Substances 0.000 claims description 128
- 230000008569 process Effects 0.000 claims description 66
- 238000012545 processing Methods 0.000 claims description 62
- 238000004364 calculation method Methods 0.000 claims description 60
- 238000006243 chemical reaction Methods 0.000 claims description 49
- 230000009466 transformation Effects 0.000 claims description 47
- 230000003287 optical effect Effects 0.000 claims description 11
- 210000004556 brain Anatomy 0.000 claims description 9
- 230000001678 irradiating effect Effects 0.000 claims 1
- 210000003128 head Anatomy 0.000 description 156
- 238000010586 diagram Methods 0.000 description 43
- 230000006870 function Effects 0.000 description 23
- 230000014509 gene expression Effects 0.000 description 18
- 238000012876 topography Methods 0.000 description 15
- 102000001554 Hemoglobins Human genes 0.000 description 12
- 108010054147 Hemoglobins Proteins 0.000 description 12
- 230000008859 change Effects 0.000 description 10
- 239000013307 optical fiber Substances 0.000 description 10
- 238000004070 electrodeposition Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 3
- 210000001061 forehead Anatomy 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 230000001131 transforming effect Effects 0.000 description 3
- 0 C*C[C@](I)[I-]=[N+]*C Chemical compound C*C[C@](I)[I-]=[N+]*C 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000001936 parietal effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002123 temporal effect 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
-
- 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/14553—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 specially adapted for cerebral tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6814—Head
-
- 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Neurology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は生体光計測装置、及び光照射位置及び光検出位置、又は測定チャンネルの位置表示方法に係り、特に、近赤外光を生体に照射し、生体内部を通過或いは生体内部で反射した光を計測し生体内部の血液循環、血行動態及びヘモグロビン量変化を計測する生体光計測装置、及び光照射位置及び光検出位置、又は測定チャンネルの位置表示方法に関するものである。
生体光計測装置は、生体内部の血液循環、血行動態及びヘモグロビン量変化を、簡便に被検体に対し低拘束で且つ害を与えずに計測できる装置であり、近年、多チャンネル装置による測定データの画像化が実現され、臨床への応用が期待されている。
生体光計測装置の原理として、可視から近赤外の波長の光を生体に照射し、照射位置から10-50mm程度離れた位置での反射光から生体内部を測定する装置が知られている。また、多チャンネルを用いた生体光計測装置において、測定対象における計測部位と光照射位置及び光検出位置の位置関係を簡単に2次元上において、表示設定する機能を備えたものも知られている。
特許文献1には、測定対象における計測部位と光照射位置及び光検出位置並びに生体通過光強度画像との3次元的な位置関係を明確にするために、光照射位置及び光検出位置を計測し、計測された光照射位置及び光検出位置に基づいて補正された生体通過光強度画像と、計測部位の目安となる測定対象の形状画像とを重ね合わせて表示する機能を有した生体光計測装置が記載されている。測定対象の形状画像としては、測定対象のワイヤーフレーム画像、断層画像、CT画像、MRI画像が挙げられている。
上記公知文献には、光照射位置及び光検出位置を計測しなければいけない、という問題がある。
光照射位置及び光検出位置を計測しなければいけないという事は、計測者にとって、その計測手段を有する装置を使用するというコストがかかる。
また、生体光計測時に光照射位置及び光検出位置を計測する事自体は手間と時間を要し、計測者及び被検体に対し負担がかかる。特に、通常の生体光計測は多くのチャンネル数を用いるため、計測しなければいけない光照射位置及び光検出位置の個数も多くなり、大きな負担となる。
一方、測定対象の形状画像に対し、光照射位置及び光検出位置を大まかに把握し、重ね合わせ表示を行いたいというニーズが存在する。例えば、過去の生体光計測の経験から、計測部位が前頭部、後頭部といった特定の位置に概ね定まっており、詳細な位置情報を必要としない場合である。
このニーズに対しては、光照射位置及び光検出位置又は測定チャンネルの正確な3次元位置の計測結果は過剰な情報であり、上記のコスト及び負担は目的に比して過剰な負担となっている。
本発明の目的は、光照射位置及び光検出位置、又は測定チャンネルの3次元位置を実際に計測する事なく、3次元頭部画像上での前記光照射位置及び光検出位置、又は前記測定チャンネルの算出位置を表示可能にすることができる生体光計測装置を提供することにある。
上記課題を解決するために、本発明は以下の様に構成される。本発明の生体光計測装置は、近赤外光を照射する光源部と、被検体の2次元計測点における前記近赤外光の通過光強度を計測し測定点毎の通過光強度に対応する信号を測定チャンネル毎の測定データとして出力する2次元プローブと、該2次元プローブの測定データを処理して画像化する信号処理部と、画像化した測定データを表示する表示部とを備えた生体光計測装置において、表示用の頭部形状に関するデータを記憶する記憶部と、前記頭部形状に関するデータから選択された2次元頭部画像上において設定された前記2次元プローブの位置情報を座標変換し3次元頭部画像上での光照射位置及び光検出位置、又は前記測定チャンネルの位置を算出する座標変換部を有した制御部とを備えたことを特徴とする。
以上、本発明によれば、計測者が光照射位置及び光検出位置の3次元位置を計測する事なく、3次元頭部画像上での光照射位置及び光検出位置、又は測定チャンネルの位置を表示する事が可能となる。
本発明の代表的な実施例によれば、生体光計測装置がプローブ位置簡易入力部を備えており、計測者が2次元上に表示した頭部画像上で光照射位置及び光検出位置を設定し、制御部が、設定された情報に従って3次元上での光照射位置及び光検出位置を計算する。これにより、光照射位置及び光検出位置の3次元位置を計測する事無く、生体通過光強度画像と計測部位の目安となる測定対象の形状画像とを3次元画像上で重ね合わせて表示する事が可能となる。
以下、本発明の生体光計測装置について、全体的な装置構成、具体的な構成例など、図面を参照しながら説明する。
(装置構成)
本発明の生体光計測装置は、近赤外光を生体内に照射し、生体の表面近傍から反射或いは生体内を通過した光(以下、単に通過光という)を検出し、光の強度に対応する電気信号を発生する装置である。この生体光計測装置100は、図1に示すように、近赤外光を照射する光源部101と、通過光を計測し、電気信号に変換する光計測部102と、光源部101及び光計測部102の駆動を制御する制御部103とを備えている。
本発明の生体光計測装置は、近赤外光を生体内に照射し、生体の表面近傍から反射或いは生体内を通過した光(以下、単に通過光という)を検出し、光の強度に対応する電気信号を発生する装置である。この生体光計測装置100は、図1に示すように、近赤外光を照射する光源部101と、通過光を計測し、電気信号に変換する光計測部102と、光源部101及び光計測部102の駆動を制御する制御部103とを備えている。
光源部101は、所定の波長の光を放射する半導体レーザ104と、半導体レーザ104が発生する光を複数の異なる周波数で変調するための変調器を備えた複数の光モジュール105とを備え、各光モジュール105の出力光はそれぞれ光ファイバ106を介して被検体107の所定の計測領域、例えば頭部の複数箇所から照射される。なお、プローブホルダ108は被検体107に取り付けられており、光ファイバ106や検出用光ファイバ109がプローブホルダ108に固定されている。本発明では、プローブホルダ108内の光ファイバ106の光照射位置及び検出用光ファイバ109の光検出位置の略中間点を、測定チャンネルと定義する。
プローブホルダ108は、被検体107の頭部形状に合わせて光照射位置及び光検出位置が所定の間隔で正方格子状に配列されており、それぞれ隣接する光照射位置と光検出位置との中間位置である計測位置(測定チャンネル)において、被検体107の脳に対する刺激が加えられていない時と脳に刺激が加えられている時との酸素化ヘモグロビン濃度変化及び脱酸素化ヘモグロビン濃度変化並びにヘモグロビン濃度総量の変化を求める。
光計測部102は、計測領域の複数の計測箇所から検出用光ファイバ109を介して誘導された通過光をそれぞれ光量に対応する電気量に変換するフォトダイオード等の光電変換素子110と、光電変換素子110からの電気信号を入力し、光照射位置に対応した変調信号を選択的に検出するロックインアンプ111と、ロックインアンプ111の出力信号をデジタル信号に変換するA/D変換器112とからなる。
光源部101は、n個(nは自然数)の光モジュールを備える。光の波長は生体内の注目物質の分光特性によるが、HbとHbO2の濃度から酸素飽和度や血液量を計測する場合には600nm〜1400nmの波長範囲の光の中から1あるいは複数波長選択して用いる。
光源部101は酸素化ヘモグロビンと脱酸素化ヘモグロビンの2種類の測定対象に対応して2種類の波長、例えば780nm及び830nmの光を発生するように構成され、これら二波長の光は合成され一つの照射位置から照射される。ロックインアンプ111は光照射位置とこれら二波長に対応した変調信号を選択的に検出する。光照射位置と検出位置との間の点(計測点)の数の2倍のチャンネル数のヘモグロビン量変化信号が得られる。制御部103には、さらに、信号処理部113、記憶部115、表示部114やマウス117を含む入出力部116が接続されている。さらに、被検体107の近傍には刺激呈示部118が設置され、制御部103等で生成された所定の刺激が被検体107に提示(表示)される。
記憶部115は、信号処理部113の処理に必要なデータや処理結果を記憶している。表示用の頭部形状に関するデータとして、例えば、被検体107の年齢や性別に応じた複数の表示用の「頭部形状の標準的なパターン」が、2次元の平面図形と3次元の形状(ワイヤーフレーム像)のデータが各々対をなす形でテーブル化され、記憶部115に記憶されている。この頭部形状の標準的なパターンやプローブの形状は、過去の生体光計測のデータ例を基に、被検体107に対する計測目的に応じて計測者(使用者)が被検体の特定の計測部位、例えば、前頭部、後頭部といった計測位置を把握できるに足る位置情報として予め準備されたものである。表示用の頭部形状に関するデータとして、さらに、生体光計測装置に備え付けられた複数のプローブホルダの各形状や、各プローブホルダ内の光ファイバや検出用光ファイバの位置関係のデータがセットになって、所定のパターンを有する「プローブ」のデータとしてテーブル化され、記憶部115に記憶されている。
制御部103や信号処理部113は、CPU、メモリ、ソフトウェア等により実現にされる各種の演算処理機能や画像処理表示機能を有している。演算処理機能により、例えば、デジタル信号に変換されたヘモグロビン量の変化信号を処理する。この処理された情報に基づいて、信号処理部113で酸素化ヘモグロビン濃度変化、脱酸素化ヘモグロビン濃度変化、全ヘモグロビン濃度変化などをチャンネル毎に示すグラフやそれを被検体の2次元画像上にプロットした画像を作成する。この信号処理部113の処理結果は、入出力部116の表示部114に表示される。
入出力部116は、計測者が装置の動作に必要な種々の指令を入力したり、測定結果や処理結果を出力するため機能を備えている。表示部114の表示画面はグラフィカルユーザーインターフェース(GUI)機能を有し、使用者が表示画面上のアイコンや「プローブ」位置等をマウス117などで操作することで各種の測定や処理に必要な指令や座標位置等の情報を入力するための入力部として機能する。例えば、計測者が、被検体107のプローブホルダ108の位置、すなわち光照射位置及び光検出位置、又は測定チャンネルの3次元位置を実際に計測する事なく、これに相当する情報として、マウス操作で2次元の「頭部形状」を選択し、さらに、この「頭部形状」に対する2次元の「プローブ」の位置を設定することで、「頭部形状」上におけるプローブの位置座標の情報を入力するのに使用される。なお、2次元の「プローブ」の位置座標入力のための表示部114と、信号処理部113の処理結果を頭部形状のプローブの位置に重ねて表示する表示部とは、必要に応じて、別個に設けても良いことは言うまでもない。
図2に、本発明の生体光計測装置100における制御部103や表示部114等が備える主な機能を示す。生体光計測装置100の表示部114は、記憶部115の「標準的な頭部形状」等のデータから被検体に最適な「頭部形状の標準的なパターン」や「プローブ」を使用者に選択させるための頭部形状データ選択部201、選択された頭部形状の標準的なパターン等に基づき生成、表示された「頭部形状」の画面像上で被検体のプローブ位置を使用者に入力、設定させるためのプローブ位置簡易入力部202、座標変換された画像や、頭部形状上の光照射位置及び光検出位置、又は測定チャンネルの算出位置にトポグラフィ像を重ね合わせた画像を表示する表示部209を備えている。
生体光計測装置100の制御部103は、入力されたプローブ位置の座標変換処理を行い3次元頭部画像上での光照射位置及び光検出位置、又は測定チャンネルの位置を算出する座標変換部203、被検体に刺激を呈示する刺激提示部206、光照射位置に基づきトポグラフィ像を生成するトポグラフィ像生成部207、及び入力されたプローブ位置に基づき頭部形状上の光照射位置及び光検出位置、又は測定チャンネルの算出位置にトポグラフィ像を重ね合わせて生成する重ね合わせ処理部208を備えている。座標変換部203は、プローブ位置簡易入力部202で設定された2次元の位置情報に従って3次元頭部画像上での光照射位置及び光検出位置、又は測定チャンネルの位置を計算する機能を備えている。
なお、座標変換部203は、2次元座標から3次元座標に変換する2次元座標−3次元座標変換部204と、3次元座標から2次元座標に変換する3次元座標−2次元座標変換部205とがあり、例えば、2次元頭部画像から前記3次元頭部画像へ、途中で近似した形状間での座標変換処理を順次繰り返すことにより、光照射位置及び光検出位置、又は測定チャンネルの算出位置を求める。後で詳細に述べるように、座標変換部203は、使用者の情報入力のため、あるいは演算結果の出力、表示のために、プローブ位置は2次元表示または3次元表示のいずれも変換可能である。なお、本発明の生体光計測装置100が備える上記主な機能は、ソフトウェアを用いて制御部103や表示部114等の全体で実現されるものであり、図2に示したブロック構成に限定されるものでないことは言うまでもない。
換言すると、本発明の生体光計測装置100は、GUI機能を有する表示部を利用した光照射位置及び光検出位置、又は測定チャンネルの3次元位置の簡易表示、設定機能を備えている。すなわち、計測者がプローブホルダ108を被検体107の頭部に装着した後、被検体に装着された実際の光照射位置及び光検出位置を計測することなく、表示画面の2次元頭部画像上で被検体の光照射位置及び光検出位置を簡易設定し、このようにして設定された情報に基づき座標変換などの演算処理により3次元頭部画像上での光照射位置及び光検出位置を推定し、実際に計測されたトポグラフィ像を演算処理で得られた頭部画像上の算出位置に重ねて表示する機能を実現することができる。これにより、生体光計測時に被検体に装着された実際の光照射位置及び光検出位置を計測する手間と時間を省略し、計測者及び被検体に対する負担を大幅に低減することができる。
(1．0 全体の具体的な構成例)
次に、本発明の生体光計測装置100の各機能を実現するより具体的な構成について、図3のフローチャート及び図4〜図10の表示画面に表示される操作画面等の例を用いて説明する。
次に、本発明の生体光計測装置100の各機能を実現するより具体的な構成について、図3のフローチャート及び図4〜図10の表示画面に表示される操作画面等の例を用いて説明する。
まず、使用者が頭部形状データ選択部201により、表示画面を利用して被検体に最適な「頭部形状の標準的なパターン」や「プローブ」を選択する。図4(図4A、図4B、図4C)は本発明における表示画面の例を示している。
本実施例においては、操作の簡便性を優先した操作画面2011を想定する。すなわち、表示部114の表示画面には、図4Aに示すような操作画面2011が表示される。この操作画面2011には、光照射位置及び光検出位置、又は測定チャンネルすなわち2次元頭部画像2012、複数の光照射位置及び光検出位置を全体として一つの塊とした一つのプローブ(2次元プローブ)2013、2次元プローブの中心座標位置2014、2次元頭部画像の中心座標位置2015が表示される。2019は各種の指令やデータの入力欄、プルダウン形式の操作ボタンなどのアイコンが配置された操作入力部である。操作画面2011における2次元頭部画像2012、2次元プローブ2013の形状、位置や大きさは、マウス操作等により使用者が、被検体107の年齢や、測定目的、使用するプローブホルダ108の種類などに応じて、任意に選択し、あるいは変更して表示可能である。
図4Aに示された2次元頭部画像2012は、頭部を上から見た平面画像である。使用されるプローブホルダ108の形状に対応して、一つの正方形のプローブに複数の光照射位置及び光検出位置が規則的に配列され、図形の中心にプローブ中心点が存在している。なお、頭部上に装着された状態でのプローブホルダ108の上面図に対応するように、プローブホルダ108の設置位置によって、光照射位置や光検出位置の間隔を不等間隔とし、図形の中心位置もプローブ内でシフトして2次元頭部画像上に表示するようにしても良い。
また、2次元頭部画像として、測定目的に応じて頭部を種々の方向から見た画像、例えば、左右の側頭部、正面、あるいは後頭部の2次元頭部画像も位置設定用の2次元画像として使用者が選択できるように、データベースを構成しても良い。
次に、選択された2次元頭部画像2012上でのプローブ2013の位置を使用者がソフトウェア上で設定する(S2010)。
次に、計測者が簡易操作で設定した2次元頭部画像上でのプローブ位置の設定位置情報を用いて、3次元頭部画像上での光照射位置及び光検出位置の計算を行なう(S2110)。
まず、使用者がソフトウェア上で設定を行なった2次元頭部画像上でのプローブ2013の位置を基に、座標変換処理を行う。この座標変換処理は、2次元頭部画像から3次元頭部画像へ、途中で近似した形状間での座標変換処理を順次繰り返しながら実行され、光照射位置及び光検出位置、又は測定チャンネルの算出位置を求める。すなわち、プローブ中心点の位置に関して、座標変換処理A(S2020)、座標変換処理B(S2040)、座標変換処理C(S2060)を行い、3次元半楕円体面上でのプローブ中心点の位置を求める(S2060)。
さらに、楕円弧の長さの計算処理、3次元半楕円体面上でのプローブ中心点の座標位置から光照射位置、光検出位置を計算する計算処理(S2080)が実行される。これらの計算処理の結果に基づき、さらに、3次元頭部画像上の光照射位置及び光検出位置、又は測定チャンネルすなわち光照射位置及び光検出位置の計算(S2110)を行なう。なお、途中の近似した形状の種類や座標変換処理の回数は上記例に限定されるものではない。例えば、途中の上記座標変換処理C(S2060)を省略しても良い。あるいは、2次元頭部画像から座標変換処理B、Cまでを行い3次元頭部画像に代えて3次元半球画像を最終的な画像とし、この3次元半球画上でのプローブ中心点の座標位置から光照射位置、光検出位置、又は測定チャンネルを計算するようにしても良い。あるいは3次元頭部画像に代えて3次元半球画像を最終的な画像とし、光照射位置及び光検出位置、又は測定チャンネルを計算するようにしても良い。
そして、3次元頭部画像上の光照射位置及び光検出位置、又は測定チャンネル(光照射位置及び光検出位置)が計算(S2110)された後、生体通過光強度画像と計測部位の目安となる測定対象の形状画像とを3次元上で重ね合わせて表示する(S2120)。
なお、図4Bに示すように、プローブ位置を使用者がソフトウェア上で設定するステップ(S2010)において、操作画面2011には、2次元頭部画像2012に重ねて3次元頭部画像(ワイヤーフレーム像)3100も同時に表示可能である。これは、2次元と3次元の形状データが各々対をなしてテーブル化され記憶装置に保持されていることや、後で詳細に述べる2次元と3次元の座標変換の双方向処理機能に基き、演算処理の結果が直ちに操作画面上に反映される結果によるものである。重ね合わせる測定対象の3次元頭部形状画像としては、本実施例においては、被検体の頭部の形状フレーム画像を用いる。図4Bにおいて、2次元プローブ2013内の円及び矩形角内の数値は、被検体の左頭部における照射位置及び光検出位置を示している。また、3次元プローブ2017には、3次元頭部画像3100が表示されている。3次元頭部画像3100における3次元プローブ2018の位置は、2次元頭部画像2012における2次元プローブ2013の位置に対応した左頭部にある。
図4Cは、2次元プローブ2013の形状として横長の矩形がマウス操作で選択され、位置も前頭部にある状態を示しており、このとき、3次元頭部画像3100内における3次元プローブ2018の位置も、前頭部にある。座標変換処理の双方向処理機能により、マウス操作等により使用者が3次元プローブ2018の形状、位置を変えることで、2次元頭部画像2012の形状、位置を変えることもできる。また、3次元頭部画像3100をマウス操作などで回転させ、3次元プローブ2018の形状、位置を見やすくすることも可能である。このとき、2次元頭部画像内における2次元プローブの位置も連動して、変形、移動表示される。このように、操作画面2011に、2次元頭部画像と3次元頭部画像を同時に表示することで、計測者が簡易操作で設定した表示画面の頭部画像上における被検体の光照射位置及び光検出位置をより正確に確認し、必要に応じて修正することを可能にしている。状況に応じて、計測者がマウス等の簡易操作により3次元頭部画像上で3次元プローブ2018を操作し、プローブホルダ108の3次元位置座標を入力するようにしても良い。
ここで、2次元頭部画像上でのプローブ位置2010を使用者がソフトウェア上で設定する処理S2010について、説明を加える。
まず、図4Aの操作画面2011に表示された2次元頭部画像2012上でのプローブ位置2010を使用者が、マウス操作などにより設定する。なお、本実施例においては、2次元頭部画像2012が光照射位置及び光検出位置のおおよその位置を示すために用いられ、2次元頭部画像2012上での光照射位置及び光検出位置の位置情報の正確さがそれほど無い場合を想定する。
以上の想定条件において、2次元頭部画像上のプローブの中心座標位置2014を2次元頭部画像2012上での使用者の設定位置情報として利用する事とする。
次に、図3のフローにおける［順方向の処理］、すなわち2次元頭部画像上→3次元頭部画像上の座標変換処理A〜D、及び計算処理A(S2080)について説明する。
順方向の処理では、まず、座標変換処理A(S2020)を行い、2次元頭部画像2012上のプローブの中心座標位置2014から、図5に示す、2次元円画像2032上でのプローブ中心点の座標位置2031を求める(S2030)。なお、2次元円画像2032の中心は、2次元頭部画像の中心座標位置2015と一致する。
そして、座標変換処理B(S2040)を行い、2次元円画像2012上でのプローブ中心点の座標位置2031から、図6に示す、3次元半球面2052上でのプローブ中心点の座標位置2051を求める(S2050)。
そして、座標変換処理C(S2060)を行い、3次元半球面上でのプローブ中心点の座標位置2051から、図7に示す、3次元半楕円体面2072上でのプローブ中心点の座標位置2071を求める(S2070)。
さらに、計算処理A(S2080)により、3次元半楕円体面上でのプローブ中心点の座標位置2071から、図8に示す、3次元半楕円体面2072上での光照射位置及び光検出位置2091を求める(S2090)。計算処理A(S2080)の詳細に関しては後に説明する。
そして、座標変換処理D(S2100)を用いて、3次元半楕円体面上での光照射位置及び光検出位置2090から、図9に示す、3次元頭部画像3100上での光照射位置及び光検出位置2111を求める(S2110)。座標変換の途中で3次元半楕円体面2072上での位置の算出を行うのは、3次元半楕円体面2072の3次元形状が、3次元半球面2052よりも、より3次元頭部画像3100の3次元形状に近いものであるため、座標変換の処理が容易になる利点があるためである。
なお、2次元頭部画像上での使用者による設定位置の情報としては、2次元頭部画像上での光照射位置及び光検出位置の位置情報を用いて、本実施例と同様の座標変換処理を行ない、3次元頭部画像上での光照射位置及び光検出位置を算出する(S2110)事も可能である。
3次元頭部画像上の光照射位置及び光検出位置2111が計算された後、図10に示すように、生体通過光強度画像2121と計測部位の目安となる測定対象の形状画像とを3次元上で重ね合わせて表示する(S2120)。
次に、形状画像との重ね合わせ表示の処理(S2120)について、説明を補足する。光照射位置に基づくトポグラフィ像生成部207は、読み込んだ光照射位置及び光検出位置2111からそれぞれ隣接する光照射位置及び光検出位置の三次元空間における距離ｄを各計測位置毎に計算する。このとき本実施の形態では、計測位置数は16個となるので、計算によって算出された光照射位置及び光検出位置の三次元空間における距離ｄは24個となる。次に、トポグラフィ像生成部207は、補正前のトポグラフィ像である二次元の光照射位置及び光検出位置に基づいたトポグラフィ像Hb data(original)を下記の式(1)に従って補正し、補正後のトポグラフィ像である三次元の光照射位置及び光検出位置に基づいたトポグラフィ像Hb data(new)を生成する。具体的には、3次スプライン補間によって得られたトポグラフィ像の各画素毎の画素値を、その画素に最も近い計測位置の距離ｄで補正する。なお、このときの補正は、下記の式(1)に従う。
Hb data(new)＝Hb data(original)×Distance・・・(1)
ただし、Distance＝d/30
ただし、Distance＝d/30
次に、トポグラフィ像生成部207は、補正後のトポグラフィ像を先に生成したワイヤーフレーム像に沿ったトポグラフィ像に変換する。具体的には、トポグラフィ像を構成する各画素の画素値を、頭部ワイヤフレーム座標に合わせて3次元マッピング補間を行うことにより、三次元トポグラフィ像画像へ変換する。
この後に、トポグラフィ像生成部207は、ワイヤーフレーム像にトポグラフィ像を重ねた三次元像を生成し、その三次元像を入力手段から指示された視点位置から見た二次元の画像(「三次元的トポグラフィ像」)に変換し、表示手段の表示面上に表示させる。
次に、図3のフローチャートにおける［逆方向の処理］、すなわち、3次元頭部画像上→2次元頭部画像上の処理について、説明する。
先にも述べたとおり、座標変換処理A(S2020)、座標変換処理B(S2040)、座標変換処理C(S2060)、座標変換処理D(S2100)は可逆な座標変換処理であるため、3次元頭部画像上で3次元頭部画像上でのプローブ位置を使用者が設定する事で、2次元頭部画像上でのプローブ位置を計算し、表示させる事も可能である。
まず、例えば図4Bに示した状態において、3次元頭部画像3100上で、3次元頭部画像上でのプローブ位置2017を使用者がソフトウェア上で設定する。ここでは、簡単に、図9に示す、3次元頭部画像3100上でのプローブ2111の中心点に相当する位置が設定されるものとする。
次に、3次元頭部画像上でのプローブ中心点から座標変換処理D(S2100)を用いて、図7に示す、3次元半楕円体面2072上でのプローブ中心点の座標位置2071が求まる。
次に、3次元半楕円体面上でのプローブ中心点の座標位置2070から座標変換処理C(S2060)を用いて、図6に示す、3次元半球面2052上でのプローブ中心点の座標位置2051が求まる。さらに、3次元半球面上でのプローブ中心点の座標位置2051から座標変換処理B(S2040)を用いて、図5に示す、2次元円画像2032上でのプローブ中心点の座標位置2031が求まる。この2次元円画像上でのプローブ中心点の座標位置2031から、座標変換処理B(S2020)を用いて、2次元頭部画像上でのプローブ2013の位置が求まる。
3次元半楕円体面2072上での光照射位置及び光検出位置2091、3次元頭部画像3100上での光照射位置及び光検出位置2111は、先に説明した通り、3次元半楕円体面2072上でのプローブ中心点の座標位置2071から求められる。
以上、本発明によれば、計測者が2次元上に表示した頭部画像上で光照射位置及び光検出位置を設定し、設定された情報に従って3次元上での光照射位置及び光検出位置を計算する機能を実現する事で、光照射位置及び光検出位置の3次元座標を計測する事無く、生体通過光強度画像と計測部位の目安となる測定対象の形状画像とを3次元上で重ね合わせて表示する事が可能となる。すなわち、被検体の頭部における光照射位置及び光検出位置の3次元座標を計測者が計測する必要なしに、被検体の生体通過光強度画像と計測部位の目安となる測定対象の形状画像とを3次元画像上で重ね合わせて表示する事が可能となる。光照射位置及び光検出位置の正確な計測結果を必要としないため、位置計測の省略による生体光計測全体の時間を短縮し、コストの節減及び計測者の負担軽減を図ることができる。
次に、本発明の実施例における、座標変換処理について、より詳細に説明する。
(2．1．2次元頭部画像⇔2次元円画像)
まず、座標変換処理A(S2020)に関し、図11を用いて、説明する。
まず、座標変換処理A(S2020)に関し、図11を用いて、説明する。
(2．1．1．座標変換処理A 2次元頭部画像→2次元円画像)
座標変換処理A(S2020)により、2次元頭部画像上の点を2次元円画像上の点に一意に対応させる計算手順を説明する。2次元円画像は半径Rc、中心点(0 0)の円とする。
座標変換処理A(S2020)により、2次元頭部画像上の点を2次元円画像上の点に一意に対応させる計算手順を説明する。2次元円画像は半径Rc、中心点(0 0)の円とする。
各座標点、偏角を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
(2．1．2． 座標変換処理A 2次元円画像→2次元頭部画像)
座標変換処理A(S2020)により、2次元円画像上の点を2次元頭部画像上の点に一意に対応させる計算手順を説明する。2次元円画像は半径Rc、中心点(0 0)の円とする。
各座標点、偏角を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
座標変換処理A(S2020)により、2次元円画像上の点を2次元頭部画像上の点に一意に対応させる計算手順を説明する。2次元円画像は半径Rc、中心点(0 0)の円とする。
各座標点、偏角を下記のように表記する事とする。
(2．2． 2次元円画像⇔3次元半球面)
座標変換処理B(S2040)に関し、図12を用いて、説明する。
座標変換処理B(S2040)に関し、図12を用いて、説明する。
(2．2．1． 座標変換処理B 2次元円画像→3次元半球面)
座標変換処理B(S2040)により、2次元円画像上の点を3次元半球面上の点に一意に対応させる計算手順を説明する。
以下のように座標変換の計算を行なう。
座標変換処理B(S2040)により、2次元円画像上の点を3次元半球面上の点に一意に対応させる計算手順を説明する。
以下のように座標変換の計算を行なう。
(2．2．2． 座標変換処理B 3次元半球面→2次元円画像)
座標変換処理B(S2040)により、3次元半球面上の点を2次元円画像上の点に一意に対応させる計算手順を説明する。
座標変換処理B(S2040)により、3次元半球面上の点を2次元円画像上の点に一意に対応させる計算手順を説明する。
2次元円画像は半径Rc、中心点(0 0)の円とし、3次元半球面は半径Rs、中心点(0 0 0)の半球面とする。
各座標点を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
以下のように座標変換の計算を行なう。
(2．3． 3次元半球面⇔3次元半楕円体面)
座標変換処理C(S2060)に関し、図13を用いて、説明する。
座標変換処理C(S2060)に関し、図13を用いて、説明する。
(2．3．1． 座標変換処理C 3次元半球面→3次元半楕円体面)
座標変換処理C(S2060)により、3次元半球面上の点を3次元半楕円体面上の点に一意に対応させる計算手順を説明する。
座標変換処理C(S2060)により、3次元半球面上の点を3次元半楕円体面上の点に一意に対応させる計算手順を説明する。
3次元半球面は半径Rs、中心点(0 0 0)の半球面とする。3次元半楕円体面は半径Rvx, Rvy, Rvz、中心点(0 0 0)の半楕円体面とする。
各座標点を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
以下のように座標変換の計算を行なう。
(2．3．2． 座標変換処理C 3次元半楕円体面→3次元半球面)
座標変換処理C(S2060)により、3次元半楕円体面上の点を3次元半球面上の点に一意に対応させる計算手順を説明する。
3次元半楕円体面は半径Rvx, Rvy, Rvz、中心点(0 0 0)の半楕円体面とする。3次元半球面は半径Rs、中心点(0 0 0)の半球面とする。
各座標点を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
座標変換処理C(S2060)により、3次元半楕円体面上の点を3次元半球面上の点に一意に対応させる計算手順を説明する。
3次元半楕円体面は半径Rvx, Rvy, Rvz、中心点(0 0 0)の半楕円体面とする。3次元半球面は半径Rs、中心点(0 0 0)の半球面とする。
各座標点を下記のように表記する事とする。
以下のように座標変換の計算を行なう。
(2．4． 3次元半楕円体面⇔3次元頭部画像)
座標変換処理D(S2100)に関し、図14を用いて、説明する。
座標変換処理D(S2100)に関し、図14を用いて、説明する。
(2．4．1． 3次元半楕円体面→3次元頭部画像)
座標変換処理D(S2100)により、3次元半楕円体面上の点を3次元頭部画像上の点に一意に対応させる。
本計算処理はポリゴンの形状に依存しないが、本実施例では三角形ポリゴンで構成されているものとする。
座標変換処理D(S2100)により、3次元半楕円体面上の点を3次元頭部画像上の点に一意に対応させる。
本計算処理はポリゴンの形状に依存しないが、本実施例では三角形ポリゴンで構成されているものとする。
各々の三角形ポリゴンを構成する3点を
と表記する。各座標点を下記のように表記する事とする。
直線OvPvと3次元頭部画像との交点を点Phとして定義し、以下のように計算を行なう。
(Ｉ)各々の三角形ポリゴンに対し、下記の計算を行う。
なお、図15に、三角形ポリゴンと直線の関係を示す。
(Ｉ)各々の三角形ポリゴンに対し、下記の計算を行う。
なお、図15に、三角形ポリゴンと直線の関係を示す。
(2．4．2． 3次元頭部画像→3次元半楕円体面)
座標変換処理D(S2100)により、3次元頭部画像上の点を3次元半楕円体面上の点に一意に対応させる。3次元半楕円体面は半径Rvx, Rvy, Rvz、中心点(0 0 0)の半楕円体面とする。3次元頭部画像の中心点は(0 0 0)であるものとする。
各座標点を下記のように表記する事とする。
直線OvPhと3次元半楕円体面との交点を点Pvとして定義し、以下のように求める。
座標変換処理D(S2100)により、3次元頭部画像上の点を3次元半楕円体面上の点に一意に対応させる。3次元半楕円体面は半径Rvx, Rvy, Rvz、中心点(0 0 0)の半楕円体面とする。3次元頭部画像の中心点は(0 0 0)であるものとする。
各座標点を下記のように表記する事とする。
直線OvPhと3次元半楕円体面との交点を点Pvとして定義し、以下のように求める。
(3． 3次元頭部画像の3次元半楕円体への近似計算)
3次元頭部画像を3次元半楕円体を用いて近似する。近似の方法は幾つか考えられるが、本実施例では単純に3次元頭部画像から3次元半楕円体の3つの径Rvx、Rvy、Rvzを求める事とする。
3次元頭部画像中の右耳Ph,R、左耳Ph,L、NasionPh,N、頭頂Ph,Tの4点を用いて、3次元半楕円体の3つの径を求める手順を、図16を用いて説明する。図16は、3次元頭部画像の3次元半楕円体への近似計算処理を示す図である。
3次元頭部画像を3次元半楕円体を用いて近似する。近似の方法は幾つか考えられるが、本実施例では単純に3次元頭部画像から3次元半楕円体の3つの径Rvx、Rvy、Rvzを求める事とする。
3次元頭部画像中の右耳Ph,R、左耳Ph,L、NasionPh,N、頭頂Ph,Tの4点を用いて、3次元半楕円体の3つの径を求める手順を、図16を用いて説明する。図16は、3次元頭部画像の3次元半楕円体への近似計算処理を示す図である。
(4． 楕円弧の長さの計算)
楕円弧の長さの計算方法について、図17に基づき説明する。図17は楕円弧の長さの計算例を示す図である。
楕円弧の長さを求める問題は、楕円積分の計算を行なう問題として知られており、初等関数によって求める事が出来なく、数値積分により、近似的に計算する事が可能である。数値積分は幾つかの手法が提案されているが、本実施例においては、拡張中点則を用いた区分求積法により計算を行なう。
楕円の弧の長さＬは次の定積分にて表現される。
楕円弧の長さの計算方法について、図17に基づき説明する。図17は楕円弧の長さの計算例を示す図である。
楕円弧の長さを求める問題は、楕円積分の計算を行なう問題として知られており、初等関数によって求める事が出来なく、数値積分により、近似的に計算する事が可能である。数値積分は幾つかの手法が提案されているが、本実施例においては、拡張中点則を用いた区分求積法により計算を行なう。
楕円の弧の長さＬは次の定積分にて表現される。
(II)他の場合
下表のように、定積分を分解し、分解された定積分を(Ｉ)の近似計算に従って計算する。
下表のように、定積分を分解し、分解された定積分を(Ｉ)の近似計算に従って計算する。
(5． 楕円弧の長さが指定距離となる点座標の計算)
楕円上のある点から指定距離となる点座標を計算する方法について、図17に基づき説明する。
楕円上のある点から指定距離となる点座標を計算する方法について、図17に基づき説明する。
楕円弧の長さは先に述べた通り、演繹的に計算する事は出来なく、近似計算により求める事になるため、可能な範囲における偏角βの候補値に対し楕円弧の長さの計算を行ない、指定距離に最も近くなる偏角を採用する事とする。具体的には以下のように計算を行なう。
(6． 実際の計算処理)
使用者がソフトウェア上で設定を行なった2次元頭部画像上でのプローブ位置(S2010)を元に3次元頭部画像上の光照射位置及び光検出位置(S2110)の計算を行なう手順を以降、説明する。
使用者がソフトウェア上で設定を行なった2次元頭部画像上でのプローブ位置(S2010)を元に3次元頭部画像上の光照射位置及び光検出位置(S2110)の計算を行なう手順を以降、説明する。
(プローブ中心点の座標位置から光照射位置、光検出位置を計算する計算処理A(S2080))
計算処理A(S2080)に関して図18、図19、図20を用いて以下に説明する。図18は、2次元頭部画像上でのプローブ位置の設定例を説明する図、図19は、2次元平面上におけるプローブの形状情報(回転前)を説明する図、図20は、2次元平面上におけるプローブの形状情報(回転後)を説明する図である。
計算処理A(S2080)に関して図18、図19、図20を用いて以下に説明する。図18は、2次元頭部画像上でのプローブ位置の設定例を説明する図、図19は、2次元平面上におけるプローブの形状情報(回転前)を説明する図、図20は、2次元平面上におけるプローブの形状情報(回転後)を説明する図である。
3次元半楕円体面上におけるプローブ中心点の座標位置から3次元半楕円体面上における光照射位置、光検出位置を計算する。
本実施例では、例えば、2次元頭部画像上でプローブを、図18のように配置したとする。
プローブの中心点と2次元頭部画像の中心点を結んだ基準線(901)の傾きをαuとし、基準線(901)に対するプローブの中心線(902)の傾きをθuとする。
3次元半楕円体面上におけるプローブ中心点の座標位置から3次元半楕円体面上における光照射位置、光検出位置を計算するためには、そのプローブにおける光照射位置、光検出位置のプローブ中心点との関係情報を既知のプローブの形状情報として用いる事とする。
プローブの形状情報は、通常ある一定の個数の光照射位置、光検出位置を一組にしたシート状のプローブを想定し、その形状情報を用いればよい。
完全に頭部形状に合致する形状を持ったプローブを作成する事は、被検体の頭部形状の個体差や、取り付け位置に依存した形状の違いにより、非常に困難であり、通常は、多くの被検体の頭部形状に可能な範囲で合致するような形状を有したプローブが使用される。
本実施例では、プローブの形状情報として、2次元平面上の形状と、3次元球面上に配置した形状の2つを想定した場合における計算処理A(S2080)を説明する。
(6．1． 2次元平面上の形状を利用した場合)
まず、プローブの形状情報として、2次元平面上のシート状の上に各光照射位置、光検出位置を配置した場合を想定する。例えば、プローブの形状情報として、図19に示す形状情報を用いたとする。光照射位置、光検出位置間の距離を全てWとした。
まず、プローブの形状情報として、2次元平面上のシート状の上に各光照射位置、光検出位置を配置した場合を想定する。例えば、プローブの形状情報として、図19に示す形状情報を用いたとする。光照射位置、光検出位置間の距離を全てWとした。
3次元半楕円体上での各光照射位置、光検出位置を求める際に、このプローブの形状情報から、プローブの中心から基準点までの距離(式2)と基準点から各光照射位置、光検出位置までの距離(式3)を用いる事とする。
Miの3次元半楕円体上での位置を求めるために、次の2つの計算を行なう。
(Ｉ)3次元半楕円体上での基準線に沿って、3次元半楕円体上のプローブの中心点から3次元半楕円体上で距離が(式4)となる3次元半楕円体上での基準点Kiの座標を求める。
(II)3次元半楕円体上での基準点Kiから、3次元半楕円体上で距離が(式5)となる座標点を求め、3次元半楕円体上でのMiとする。
以降、(Ｉ)、(II)の計算処理の詳細を説明する。
(Ｉ)3次元半楕円体上での基準線に沿って、3次元半楕円体上のプローブの中心点から3次元半楕円体上で距離が(式4)となる3次元半楕円体上での基準点Kiの座標を求める。
(II)3次元半楕円体上での基準点Kiから、3次元半楕円体上で距離が(式5)となる座標点を求め、3次元半楕円体上でのMiとする。
以降、(Ｉ)、(II)の計算処理の詳細を説明する。
(Ｉ)3次元半楕円体上での基準点Kiの座標を求める。
3次元半楕円体上でのプローブの中心点Pvと半楕円体の頂点を通る基準線上の点として、基準点Kiを求める手順を図21、図22を用いて説明する。図21は、基準点の座標位置の計算(3次元半楕円体)を説明する図、図22は、基準点の座標位置の計算(切断楕円面 S1)を説明する図である。
3次元半楕円体上でのプローブの中心点Pvと半楕円体の頂点を通る基準線上の点として、基準点Kiを求める手順を図21、図22を用いて説明する。図21は、基準点の座標位置の計算(3次元半楕円体)を説明する図、図22は、基準点の座標位置の計算(切断楕円面 S1)を説明する図である。
図21において、プローブの中心点Pvを次のように表記する。
基準点Kiは、基準線上の点として、
基準点Kiは、基準線上の点として、
3次元半楕円体上でのプローブの中心点Pvと半楕円体の頂点を通る基準線を含む切断面において形成される楕円上で、プローブの中心点Pvから距離が(式4)となる点を求める。(図22参照)
切断楕円面S1上では点Pv、Kiは以下のように表記される。
これに対し、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、偏角ψKiを求める事で、プローブの中心点Pvから距離が(式4)となる基準点Kiが求まる。
切断楕円面S1上では点Pv、Kiは以下のように表記される。
これに対し、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、偏角ψKiを求める事で、プローブの中心点Pvから距離が(式4)となる基準点Kiが求まる。
また、ここで、基準線をXvYv座標平面へ射影した線分がXv軸と成す角度αを以下の通りに求めておく。
(II)3次元半楕円体上におけるMiの座標を求める。
3次元半楕円体上におけるMiの座標を求める方法を以下に説明する。
点Miは、「楕円体上の点Kiにおける接平面」に垂直であり、かつ、「点Kiを通り楕円体の底面に垂直な断面(断面S1)」に垂直な断面(断面S2)と、楕円体の交線として構成される曲線(曲線L2)上で、基準点Kiからの距離が(式5)となる点として求める。その計算手順は以下の通りとなる。
3次元半楕円体上におけるMiの座標を求める方法を以下に説明する。
点Miは、「楕円体上の点Kiにおける接平面」に垂直であり、かつ、「点Kiを通り楕円体の底面に垂直な断面(断面S1)」に垂直な断面(断面S2)と、楕円体の交線として構成される曲線(曲線L2)上で、基準点Kiからの距離が(式5)となる点として求める。その計算手順は以下の通りとなる。
(i)断面S1に含まれる点Kiにおける楕円体への接線L1の式を求める。
(ii)断面S2の式を求める。
(iii)曲線L2の式を求める。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
基準点Kiの座標値(式6)に依存して、各計算手順の詳細が異なる。
基準点Kiの座標値が以下の6通りの条件の各々に関し、計算手順(i)〜(iv)の詳細に関し、以下に説明する。
(ii)断面S2の式を求める。
(iii)曲線L2の式を求める。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
基準点Kiの座標値(式6)に依存して、各計算手順の詳細が異なる。
(1)条件1の場合
直線L1は断面S1上のXsYs座標系において
と表記される。なお、図23は、条件1における3次元半楕円体上の断面S 1 、S 2 と曲線L 1 、L 2 、L 3 の関係を示す図、図24は、条件1における断面S 1 上の直線L 1 を示す図である。
直線L1は断面S1上のXsYs座標系において
と表記される。なお、図23は、条件1における3次元半楕円体上の断面S 1 、S 2 と曲線L 1 、L 2 、L 3 の関係を示す図、図24は、条件1における断面S 1 上の直線L 1 を示す図である。
図24中に示される断面S1上の座標系と図23中に示される3次元座標系の間には以下の関係式が成り立つ。
(ii)断面S2の式を求める。
図23に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
図23に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
(iii)曲線L2の式を求める。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。図25は、条件1における3次元半楕円体上の曲線L 2 とXY平面上曲線L 3 の関係を示す図である。
以上を用いて、曲線L2上の任意の点Qは以下のように表記される。
ただし、媒介変数θQは図25に示す通り、曲線L2をXY平面上に射影した曲線L3において定義した。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。図25は、条件1における3次元半楕円体上の曲線L 2 とXY平面上曲線L 3 の関係を示す図である。
以上を用いて、曲線L2上の任意の点Qは以下のように表記される。
ただし、媒介変数θQは図25に示す通り、曲線L2をXY平面上に射影した曲線L3において定義した。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
この定積分は演繹的に解を得る事が出来ないので、拡張中点則を用いた区分求積法により、近似的に解を得る事とする。拡張中点則においては、適当な分割数Nを用い、下記のように計算を行なう。
以上により、曲線L2上で基準点Kiからの距離が(式5)となる点として、点Miの座標が求まる。
(2)条件2の場合(式7の場合)
(i)断面S1に含まれる点Kiにおける楕円体への接線L1の式を求める
(i)断面S1に含まれる点Kiにおける楕円体への接線L1の式を求める
(ii)断面S2の式を求める。
図26に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
図26に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
(iii)曲線L2の式を求める。
楕円体の式は下記の通りとなる。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。
以上を用いて、曲線L2上の任意の点Qは媒介変数θQを用いて以下のように表記される。
楕円体の式は下記の通りとなる。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
この定積分は演繹的に解を得る事が出来ないので、拡張中点則を用いた区分求積法により、近似的に解を得る事とする。拡張中点則においては、適当な分割数Nを用い、下記のように計算を行なう。
以上により、曲線L2上で基準点Kiからの距離が(式5)となる点として、点Miの座標が求まる。
この定積分は演繹的に解を得る事が出来ないので、拡張中点則を用いた区分求積法により、近似的に解を得る事とする。拡張中点則においては、適当な分割数Nを用い、下記のように計算を行なう。
以上により、曲線L2上で基準点Kiからの距離が(式5)となる点として、点Miの座標が求まる。
(3)条件3の場合(ZKi=0)
(i)断面S1に含まれる点Kiにおける楕円体への接線L1の式を求める
接線L1 X=X1, Y=Y1
(ii)断面S2の式を求める。
半楕円体の底面が断面S2となるので断面S2の式は Z=0 となる。
(iii)曲線L2の式を求める。
半楕円体の底面の楕円が曲線L2なるので曲線L2の式は(式8)となる。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
曲線L2は楕円になるので、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いる事で、曲線L2上で基準点Kiからの距離が(式5)となる点の座標値を求める。
(i)断面S1に含まれる点Kiにおける楕円体への接線L1の式を求める
接線L1 X=X1, Y=Y1
(ii)断面S2の式を求める。
半楕円体の底面が断面S2となるので断面S2の式は Z=0 となる。
(iii)曲線L2の式を求める。
半楕円体の底面の楕円が曲線L2なるので曲線L2の式は(式8)となる。
曲線L2は楕円になるので、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いる事で、曲線L2上で基準点Kiからの距離が(式5)となる点の座標値を求める。
(4)条件4の場合
直線L1は
と求められる。
直線L1は
と求められる。
(ii)断面S2の式を求める。
図29に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
図29に示すように、直線L1に垂直な平面として、断面S2が下式の通り、求まる。
(iii)曲線L2の式を求める。
楕円体の式は下記の通りとなる。
断面S2の式を楕円体の式に代入し、断面S2と楕円体の交線で構成される断面線を求める。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。
以上を用いて、曲線L2上の任意の点Qは媒介変数θQを用いて以下のように表記される。
楕円体の式は下記の通りとなる。
断面S2の式を楕円体の式に代入し、断面S2と楕円体の交線で構成される断面線を求める。
上式を変形後、下式を得る。下式は曲線L2をXY平面上に射影して定義される曲線L3を示し、図23、図25に示す通りである。
以上を用いて、曲線L2上の任意の点Qは媒介変数θQを用いて以下のように表記される。
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
この定積分は演繹的に解を得る事が出来ないので、拡張中点則を用いた区分求積法により、近似的に解を得る事とする。拡張中点則においては、適当な分割数Nを用い、下記のように計算を行なう。
以上により、曲線L2上で基準点Kiからの距離が(式5)となる点として、点Miの座標が求まる。
(5)条件5の場合
(ii)曲線L2の式を求める。
(iii)楕円体における曲線L2上の点Qの座標は
(iv)曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、曲線L2上で基準点Kiからの距離が(式5)となる点を求め、その点の偏角θcを求める事で、その点の3次元半楕円体上の座標値が求まる。
先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、曲線L2上で基準点Kiからの距離が(式5)となる点を求め、その点の偏角θcを求める事で、その点の3次元半楕円体上の座標値が求まる。
(6)条件6の場合
これに対し、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
これに対し、先に説明した「楕円弧の長さが指定距離となる点座標の計算」方法を用いて、曲線L2上で基準点Kiからの距離が(式5)となる点を求める。
(6．2． 3次元球面上に配置した形状を利用した場合)
次に、プローブの形状情報として、3次元球面上に配置した形状を利用した場合を想定する。
次に、プローブの形状情報として、3次元球面上に配置した形状を利用した場合を想定する。
3次元球の半径をRsとする。プローブの中心点をP(Rs 0 0)とする。例えば、プローブの形状情報として、図35に示す形状情報を用いたとする。3次元上のx軸を中心として、3次元球面上に光照射位置、光検出位置を配置する。N点の光照射位置、光検出位置の3次元上の座標点を
とする。
とする。
次に2次元頭部画像上で設定されたプローブの回転を3次元球面上のプローブの形状情報において反映する。
3次元上のx軸を中心として、回転を行なう事とする。(図36参照)
N点の光照射位置、光検出位置の3次元上の座標点は
となり、図36のようになる。
N点の光照射位置、光検出位置の3次元上の座標点は
光照射位置、光検出位置の各点において、プローブの基準線上への足となる点を求め、光照射位置、光検出位置の各点に対する基準点K1,・・・Ki,・・・,KNとする。例えば、Miに対し、基準点Kiは
と計算される。
と計算される。
そして、3次元球面上における、プローブの中心点と各基準点の間の弧長(式2)、各基準点と光照射位置、光検出位置の各点の間の弧長(式3)を求める。例えば、Mi及びその基準点Kiに対し、(式4)、(式5)を求めると以下のようになる。
Miを通り、yz平面に水平な平面で切断した断面が図37のようになるので
と求まる。Miを通り、xy平面に水平な平面で切断した断面が図38のようになるので
と求まる。
と求まる。Miを通り、xy平面に水平な平面で切断した断面が図38のようになるので
と求まる。
プローブの中心点から基準点までの弧長(式2)と基準点から光照射位置、光検出位置の各点までの弧長(式3)をプローブの形状情報として利用して、2次元平面上の形状を利用した場合と同様に、光照射位置、光検出位置の3次元楕円半球上での座標位置が求まる。
(7． 脳波形電極位置の表示機能)
本実施例において、座標変換処理A、B、C、D(S2020、S2040、S2060、S2100)を用いて脳波形電極位置を2次元頭部画像上、2次元円画像上、3次元半球面上、3次元半楕円体面上及び3次元頭部画像上に表示する事が可能である。脳波形電極位置の表示機能に関し、図39及び図40(40A，40B，40C)を用いて以下に説明する。
図39は脳波形電極位置の表示機能のフローチャートを示したものである。
本実施例において、座標変換処理A、B、C、D(S2020、S2040、S2060、S2100)を用いて脳波形電極位置を2次元頭部画像上、2次元円画像上、3次元半球面上、3次元半楕円体面上及び3次元頭部画像上に表示する事が可能である。脳波形電極位置の表示機能に関し、図39及び図40(40A，40B，40C)を用いて以下に説明する。
図39は脳波形電極位置の表示機能のフローチャートを示したものである。
脳波形電極位置として、例えば10-20電極配置法が一般的に知られている。10-20電極配置法の詳細は、例えば非特許文献1(新臨床検査技師講座 7 臨床生理学(第3版)、医学書院、pp.174-176)などに記載されている。本実施例ではこの10-20電極配置法を用いる事とする。定められた配置法に従って、3次元半球上での各電極の座標位置(301)を求める。求めた結果は図40Aの302のようになる。
3次元半球上での各電極の座標位置(301)に対し、座標変換処理B(S2040)を施し、2次元円画像上での各電極の座標位置(303)を得る事が出来る。
さらに、2次元円画像上での各電極の座標位置(303)に対し、座標変換処理A(S2020)を施し、2次元頭部画像上での各電極の座標位置(305)を得る事が出来る。
2次元頭部画像上での各電極の座標位置(305)を用いて2次元頭部画像2012との重ね合わせ画像(306)を、図40Bのように得、呈示する事ができる。
また、3次元半球上での各電極の座標位置(301)に対し、座標変換処理C(S2060)を施し、3次元半楕円体上での各電極の座標位置(307)を得る事が出来る。さらに、3次元半楕円体上での各電極の座標位置(307)に対し、座標変換処理D(S2080)を施し、3次元頭部画像3100上での各電極の座標位置(308)を得る事が出来る。
得られた3次元頭部画像上での各電極の座標位置(308)を用いて形状画像との重ね合わせ画像(310)を、図40Cのように得、呈示する事ができる。
また、以上の各画像(3次元頭部画像、3次元半楕円体画像、3次元半球画像、2次元円画像、2次元頭部画像)上における計算結果を用いて、脳波形電極位置と生体光計測結果を重ね合わせ表示する事も可能である。
100 生体光計測装置、101 光源部、102 光計測部、103 制御部、104 半導体レーザ、105 光モジュール、106 光ファイバ、107 被検体、108 プローブホルダ、109 検出用光ファイバ、110 光電変換素子、111 ロックインアンプ、112 A/D変換器、113 信号処理部、114 表示部、115 記憶部、116 入出力部、201 頭部形状データ選択部、202 プローブ位置簡易入力部、203 座標変換部、204 2次元座標−3次元座標変換部、205 3次元座標−2次元座標変換部、206 刺激提示部、207 光照射位置に基づくトポグラフィ像生成部、208 トポグラフィ像と頭部形状の重ね合わせ処理部、209 トポグラフィ像と頭部形状の重ね合わせ表示部、2011 操作画面、2012 2次元頭部画像、2013 2次元プローブ、3100 3次元頭部画像
Claims (11)
- 近赤外光を照射する光源部と、被検体の2次元計測点における前記近赤外光の通過光強度を計測し測定点毎の通過光強度に対応する信号を測定チャンネル毎の測定データとして出力する2次元プローブと、該2次元プローブの測定データを処理して画像化する信号処理部と、画像化した測定データを表示する表示部とを備えた生体光計測装置において、
表示用の頭部形状に関するデータを記憶する記憶部と、
前記頭部形状に関するデータから選択された2次元頭部画像上において設定された2次元プローブの位置情報を座標変換し、3次元頭部画像上での光照射位置及び光検出位置、又は前記測定チャンネルの位置を算出する座標変換部を有した制御部とを備えることを特徴とする生体光計測装置。 - 請求項1において、
前記表示部は、前記2次元プローブにおける前記光照射位置及び光検出位置、又は前記測定チャンネルを構成する3次元プローブに対応する形状を表示することを特徴とする生体光計測装置。 - 請求項1において、
前記座標変換部は、
前記2次元頭部画像から前記3次元頭部画像へ、途中で近似した形状間での座標変換処理を順次繰り返し前記光照射位置及び光検出位置、又は前記測定チャンネルの算出位置を求めることを特徴とする生体光計測装置。 - 請求項3において、
前記座標変換部は、
前記2次元頭部画像上の前記プローブの中心座標位置から、2次元円画像上での2次元プローブ中心点の座標位置を求め、
該2次元円画像上でのプローブ中心点の座標位置から3次元半球面上でのプローブ中心点の座標位置を求め、
該3次元半球面上でのプローブ中心点の座標位置から3次元半楕円体面上での3次元プローブ中心点の座標位置を求め、
該3次元半楕円体面上での3次元プローブ中心点の座標位置から、3次元半楕円体面上での前記光照射位置及び光検出位置、又は前記測定チャンネルの位置を求め、前記3次元頭部画像上での前記光照射位置及び光検出位置、又は前記測定チャンネルの算出位置とすることを特徴とする生体光計測装置。 - 請求項1において、
前記表示部は、前記2次元頭部画像に重ねて前記2次元プローブと、対応する3次元頭部画像及び3次元プローブを同時に表示することを特徴とする生体光計測装置。 - 請求項1において、
前記2次元プローブからの測定データの処理結果を、前記3次元頭部画像上の前記光照射位置及び光検出位置、又は前記測定チャンネルにおける2次元プローブの算出位置に重ね合わせ、前記表示部に表示する重ね表示部を備えたことを特徴とする生体光計測装置。 - 請求項4において、
前記座標変換部は、
前記3次元頭部画像上でのプローブの中心点に相当する位置が設定されたとき、該3次元頭部画像上での3次元プローブ中心点から座標変換処理を用いて3次元半楕円体面上でのプローブ中心点の座標位置を求め、
前記3次元半楕円体面上での前記3次元プローブ中心点の座標位置から座標変換処理を用いて、前記3次元半球面上でのプローブ中心点の座標位置を求め、
前記3次元半球面上でのプローブ中心点の座標位置から座標変換処理を用いて、前記2次元円画像上での2次元プローブ中心点の座標位置を求め、
該2次元円画像上でのプローブ中心点の座標位置から、座標変換処理を用いて、前記2次元頭部画像上での前記2次元プローブの位置を求めることを特徴とする生体光計測装置。 - 請求項4において、
前記測定チャンネルに複数の脳波形電極を備え、
前記制御部は、 前記3次元半球上での各電極の座標位置に対し、座標変換処理を施し、前記2次元円画像上での前記各電極の座標位置を得て、
前記2次元円画像上での各電極の座標位置に対し、座標変換処理を施し、前記2次元頭部画像上での各電極の座標位置を得て、
前記2次元頭部画像上での前記各電極の座標位置を用いて前記2次元頭部画像との重ね合わせ画像を生成して、表示し、前記3次元半球上での前記各電極の座標位置に対し、座標変換処理を施し、前記3次元半楕円体上での前記各電極の座標位置を得て、
前記3次元半楕円体上での前記各電極の座標位置に対し、座標変換処理を施し、前記3次元頭部画像上での前記各電極の座標位置を得て、
得られた前記3次元頭部画像上での前記各電極の座標位置を用いて形状画像との重ね合わせ画像を生成して、前記表示部に表示することを特徴とする生体光計測装置。 - 請求項8において、
前記制御部は、
前記3次元頭部画像、前記3次元半楕円体画像、前記3次元半球画像、前記2次元円画像、前記2次元頭部画像の各画像上における計算結果を用いて、前記脳波形電極の位置と生体光計測結果を重ね合わせ前記表示部に表示することを特徴とする生体光計測装置。 - 生体光計測装置における光照射位置及び光検出位置、又は測定チャンネルの位置表示方法であって、
前記生体光計測装置は、近赤外光を照射する光源部と、被検体の2次元計測点における通過光強度を計測し測定点毎の通過光強度に対応する信号を測定チャンネル毎の測定データとして出力する2次元プローブと、該2次元プローブの測定データを処理して画像化する信号処理部と、表示用の頭部形状に関するデータを記憶する記憶部と、画像化した測定データを表示する表示部と、座標変換部を有した制御部とを備えており、
前記表示部に、2次元頭部画像と2次元プローブを表示し、
表示された前記2次元頭部画像上における前記2次元プローブの位置情報を座標変換し3次元頭部画像上での3次元プローブ及び前記光照射位置及び光検出位置、又は測定チャンネルの位置を算出し、
前記表示部に、前記2次元頭部画像に重ねて前記2次元プローブと、対応する前記3次元頭部画像及び3次元プローブを表示することを特徴とする生体光計測装置における光照射位置及び光検出位置、又は測定チャンネルの位置表示方法。 - 請求項10において、
前記2次元頭部画像から前記3次元頭部画像へ、途中で近似した形状間での座標変換処理を順次繰り返し前記3次元プローブの前記光照射位置及び光検出位置、又は前記測定チャンネルの算出位置を求め、 前記2次元プローブからの測定データの処理結果を、前記3次元頭部画像上の前記3次元プローブの前記光照射位置及び光検出位置、又は前記測定チャンネルの算出位置に重ね合わせ、前記表示部に表示することを特徴とする生体光計測装置における光照射位置及び光検出位置、又は測定チャンネルの位置表示方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010531878A JP5159893B2 (ja) | 2008-10-01 | 2009-09-30 | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 |
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008256510 | 2008-10-01 | ||
| JP2008256510 | 2008-10-01 | ||
| JP2010531878A JP5159893B2 (ja) | 2008-10-01 | 2009-09-30 | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 |
| PCT/JP2009/067023 WO2010038774A1 (ja) | 2008-10-01 | 2009-09-30 | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2010038774A1 JPWO2010038774A1 (ja) | 2012-03-01 |
| JP5159893B2 true JP5159893B2 (ja) | 2013-03-13 |
Family
ID=42073534
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010531878A Active JP5159893B2 (ja) | 2008-10-01 | 2009-09-30 | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US8565501B2 (ja) |
| EP (1) | EP2335580A4 (ja) |
| JP (1) | JP5159893B2 (ja) |
| CN (1) | CN102170826B (ja) |
| WO (1) | WO2010038774A1 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101661934B1 (ko) * | 2010-07-29 | 2016-10-04 | 삼성전자주식회사 | 영상 처리 장치 및 방법 |
| JP2014030621A (ja) * | 2012-08-03 | 2014-02-20 | Sony Corp | 情報処理装置、プログラム及び生体信号測定セット |
| US20150342461A1 (en) * | 2012-11-14 | 2015-12-03 | Shimadzu Corporation | Optical biometric device and position measuring device used therein |
| US9091628B2 (en) | 2012-12-21 | 2015-07-28 | L-3 Communications Security And Detection Systems, Inc. | 3D mapping with two orthogonal imaging views |
| JP7415856B2 (ja) | 2020-09-02 | 2024-01-17 | 株式会社デンソーウェーブ | 物体認識装置、物体認識方法及びプログラム |
| CN114246557A (zh) * | 2022-03-01 | 2022-03-29 | 慧创科仪(北京)科技有限公司 | 用于近红外脑功能成像装置的定位方法、设备及存储介质 |
| CN114246556B (zh) * | 2022-03-01 | 2022-05-24 | 慧创科仪(北京)科技有限公司 | 用于近红外脑功能成像装置的定位方法、设备和存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0919408A (ja) * | 1995-07-05 | 1997-01-21 | Hitachi Ltd | 生体光計測装置及びその装置における画像作成方法 |
| JP2001079008A (ja) * | 1999-09-14 | 2001-03-27 | Hitachi Medical Corp | 生体光計測装置 |
Family Cites Families (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5803909A (en) * | 1994-10-06 | 1998-09-08 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
| US6240309B1 (en) * | 1995-10-06 | 2001-05-29 | Hitachi, Ltd. | Optical measurement instrument for living body |
| US20060184047A1 (en) * | 1995-11-17 | 2006-08-17 | Yuichi Yamashita | Optical measurement instrument for living body |
| JP4006826B2 (ja) * | 1998-04-28 | 2007-11-14 | 株式会社日立製作所 | 生体光計測装置 |
| JP4076003B2 (ja) * | 1999-02-19 | 2008-04-16 | 株式会社日立製作所 | 生体光計測装置 |
| US7228166B1 (en) * | 1999-09-14 | 2007-06-05 | Hitachi Medical Corporation | Biological light measuring instrument |
| JP3839202B2 (ja) * | 1999-10-28 | 2006-11-01 | 株式会社日立製作所 | 生体光計測装置及びこの装置を機能させるためのプログラム |
| JP2003088528A (ja) * | 2001-09-18 | 2003-03-25 | Hitachi Medical Corp | 生体光計測装置 |
| JP4071475B2 (ja) * | 2001-11-12 | 2008-04-02 | 株式会社日立メディコ | 生体光計測装置 |
| JP3635332B2 (ja) * | 2003-03-20 | 2005-04-06 | 独立行政法人情報通信研究機構 | 頭部装着具 |
| US7974671B2 (en) * | 2003-09-19 | 2011-07-05 | Hitachi Medical Corporation | Living body information signal processing system combining living body optical measurement apparatus and brain wave measurement apparatus and probe device used for the same |
| JP4403453B2 (ja) * | 2003-11-13 | 2010-01-27 | 株式会社島津製作所 | 頭表座標を脳表座標に変換する方法と、その変換データを利用する経頭蓋的脳機能測定装置 |
| EP1749476A1 (en) * | 2004-05-26 | 2007-02-07 | Matsushita Electric Industrial Co., Ltd. | Optical element for measuring bioinformation and bioinformation measuring device |
| EP1818016B1 (en) * | 2004-07-20 | 2019-02-20 | Toshinori Kato | Biofunction diagnosis device, biofunction diagnosis method, bioprobe, bioprobe wearing tool, bioprobe support tool, and bioprobe wearing assisting tool |
| JP4590555B2 (ja) * | 2004-09-02 | 2010-12-01 | 国立大学法人長岡技術科学大学 | 感性状態判別方法及び装置 |
| JP2006122086A (ja) * | 2004-10-26 | 2006-05-18 | Hitachi Ltd | 生体光計測装置 |
| US8038625B2 (en) * | 2005-09-15 | 2011-10-18 | St. Jude Medical, Atrial Fibrillation Division, Inc. | System and method for three-dimensional mapping of electrophysiology information |
| JP4702107B2 (ja) * | 2006-03-03 | 2011-06-15 | 株式会社日立製作所 | 生体光計測装置 |
-
2009
- 2009-09-30 WO PCT/JP2009/067023 patent/WO2010038774A1/ja active Application Filing
- 2009-09-30 CN CN2009801391048A patent/CN102170826B/zh active Active
- 2009-09-30 US US13/121,309 patent/US8565501B2/en active Active
- 2009-09-30 EP EP09817806.4A patent/EP2335580A4/en not_active Withdrawn
- 2009-09-30 JP JP2010531878A patent/JP5159893B2/ja active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0919408A (ja) * | 1995-07-05 | 1997-01-21 | Hitachi Ltd | 生体光計測装置及びその装置における画像作成方法 |
| JP2001079008A (ja) * | 1999-09-14 | 2001-03-27 | Hitachi Medical Corp | 生体光計測装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2010038774A1 (ja) | 2012-03-01 |
| EP2335580A1 (en) | 2011-06-22 |
| CN102170826B (zh) | 2013-07-24 |
| CN102170826A (zh) | 2011-08-31 |
| WO2010038774A1 (ja) | 2010-04-08 |
| US8565501B2 (en) | 2013-10-22 |
| US20110176713A1 (en) | 2011-07-21 |
| EP2335580A4 (en) | 2015-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5159893B2 (ja) | 生体光計測装置及び、光照射位置及び光検出位置又は測定チャンネルの位置表示方法 | |
| JP4968167B2 (ja) | 光生体測定装置及びそれに用いられるホルダ配置支援システム | |
| JP7196606B2 (ja) | 眼科情報処理装置および眼科情報処理プログラム | |
| JP6273207B2 (ja) | 光生体計測装置及びそれに用いられる位置計測装置 | |
| JP5347448B2 (ja) | 生体測定装置 | |
| JP4661688B2 (ja) | 光生体計測装置、光生体計測装置用プログラム及び光生体計測方法 | |
| US7228166B1 (en) | Biological light measuring instrument | |
| US9138146B2 (en) | Biological optical measurement instrument and reference image display method | |
| JP4071475B2 (ja) | 生体光計測装置 | |
| JP4876808B2 (ja) | 脳機能データ制御装置 | |
| JP2006218196A (ja) | 生体光計測装置 | |
| JP2013244343A (ja) | 生体情報提示装置および生体情報提示方法 | |
| KR20220131207A (ko) | 유방암 진단 시스템 | |
| JP4266453B2 (ja) | 生体光計測装置 | |
| JP6011636B2 (ja) | 光生体計測装置 | |
| JP2017038878A (ja) | 生体光計測装置 | |
| JP5729490B2 (ja) | 光生体計測装置 | |
| JP2005328933A (ja) | 経頭蓋的脳機能測定装置 | |
| JP5853777B2 (ja) | 光生体計測装置 | |
| JP5484786B2 (ja) | 生体光計測装置 | |
| JP2014110844A (ja) | 生体計測装置及びそれに用いられる位置計測装置 | |
| JP2013059437A (ja) | 光生体測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120323 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121105 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121106 |
|
| 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: 20121203 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121211 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5159893 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151221 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 |
|
| 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 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |