JPH1024024A - 脳硬膜下の記録脳波による脳内双極子の追跡装置 - Google Patents

脳硬膜下の記録脳波による脳内双極子の追跡装置

Info

Publication number
JPH1024024A
JPH1024024A JP8203298A JP20329896A JPH1024024A JP H1024024 A JPH1024024 A JP H1024024A JP 8203298 A JP8203298 A JP 8203298A JP 20329896 A JP20329896 A JP 20329896A JP H1024024 A JPH1024024 A JP H1024024A
Authority
JP
Japan
Prior art keywords
brain
dipole
dimensional shape
information
dimensional
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.)
Granted
Application number
JP8203298A
Other languages
English (en)
Other versions
JP3191858B2 (ja
Inventor
Yoshio Okamoto
良夫 岡本
Toshio Saijo
寿夫 西条
Keiichi Miyamoto
啓一 宮本
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.)
Chuo Electronics Co Ltd
CHUO DENSHI KK
Original Assignee
Chuo Electronics Co Ltd
CHUO DENSHI KK
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 Chuo Electronics Co Ltd, CHUO DENSHI KK filed Critical Chuo Electronics Co Ltd
Priority to JP20329896A priority Critical patent/JP3191858B2/ja
Publication of JPH1024024A publication Critical patent/JPH1024024A/ja
Application granted granted Critical
Publication of JP3191858B2 publication Critical patent/JP3191858B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

(57)【要約】 【課題】 患者の負担の軽減と、脳外科手術箇所の決定
精度を向上させる。 【解決手段】 脳硬膜下に挿入した記録電極によって脳
波情報を収集する脳波収集装置2と、X線CT写真とM
RI写真からの情報を読み取る透過型イメージ・スキャ
ナ1と、頭蓋骨の内側立体形状と脳実質の立体形状を作
成し、記録電極位置を座標数値化して伝達行列を計算
し、さらに、推定双極子位置、電位の大きさと方向を脳
実質の立体形状上にスーパー・インポーズするコンピュ
ータ3と、スーパー・インポーズされた双極子位置等を
表示する表示装置4とによって構成した。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、脳外科手術箇所を
決定するために、脳硬膜下に挿入した記録電極を用いた
脳内双極子の追跡装置に関する。
【0002】
【従来の技術】従来、癲癇等の脳外科手術においては、
異常脳波の脳内発生箇所を決定するために脳硬膜下に脳
波記録電極を挿入しておき、長期間に亙って収集した記
録電極による脳波情報に基づいて、脳外科手術箇所を決
定していた。
【0003】
【発明が解決しようとする課題】脳硬膜下に記録電極を
挿入するためには頭蓋骨に穴を開ける手術が必要であ
り、また、記録電極による脳波情報に基づいて手術箇所
を決定するためには長期間に亙る脳波情報を収集する必
要があり、このため患者の身体を長期間拘束しなくては
ならず、患者の負担は大きかった。
【0004】本発明は、上述した従来手術による患者の
負担を軽減させるように、記録電極による脳波記録収集
期間を短縮させ、かつ、脳外科手術箇所の決定精度を向
上させることを目的とするものである。
【0005】
【課題を解決するための手段】本発明による脳硬膜下の
記録脳波による脳内双極子の追跡装置は、脳硬膜下に挿
入した記録電極によって脳電位情報を収集する脳波収集
装置と、X線CT写真およびMRI写真による脳断層写
真情報を読み取る透過型イメージ・スキャナと、この透
過型イメージ・スキャナが読み取ったX線CT写真情報
とMRI写真情報からそれぞれ構築した頭蓋骨の内側立
体形状と脳実質の立体形状を作成すると共に、前記立体
形状情報および3次元座標化した記録電極の3次元座標
値に基づいて脳2層構造での伝達行列計算を行い、さら
に、前記記録電極によって記録した脳波の1個又は2個
の等価双極子計算を行ってその計算結果を3次元立体形
状上にスーパー・インポーズするコンピュータと、これ
を表示する表示装置とによって構成した。
【0006】
【発明の実施の形態】本発明の実施例を、図面を参照し
ながら説明する。図2は、本発明の脳硬膜下の記録脳波
による脳内双極子の追跡装置の構成を示すブロック図で
ある。図2において、透過型イメージ・スキャナ1は、
X線CT写真情報やMRI写真情報を読み取ってコンピ
ュータ3の記録装置に格納させるものであり、脳波収集
装置2は、脳硬膜下に挿入された記録電極によって収集
された脳波情報をコンピュータ3の記憶装置に格納させ
るものである。コンピュータ3は、上述した各種情報に
基づいて頭蓋骨の内側立体形状と脳実質の立体形状を作
成すると共に、記録電極位置の座標数値化、伝達行列計
算、双極子位置の推定計算を行い、推定双極子の位置、
電位の大きさとその方向および測定電位と計算電位の近
似度を3次元立体形状上にスーパー・インポーズして表
示装置4に表示させる。
【0007】以下、本発明の脳硬膜下の記録脳波による
脳内双極子の追跡装置を用いて、双極子位置を推定する
手順を、図1を参照しながら説明する。
【0008】ステップ1(S1 )において、頭蓋骨の内
側立体形状を作成する。X線CT写真は頭蓋底から頭頂
まで、OMラインに平行に撮影した2〜3mm幅の複数枚
の脳断層写真を用いる。夫々の断層写真を、透過型イメ
ージ・スキャナ1によりコンピュータ3に読み込む。読
み込まれた夫々の断層写真データを8ビット/画素のグ
レー・スケールで表わし、検出しようとする頭蓋骨内側
の境界の画素の平均値を閾値として2値化し、2値化し
た画像データの周辺をトレースして境界を座標数値化
し、頭蓋骨の内側立体形状を作成する。頭蓋骨の内側立
体形状の内容としては、三角形ワイヤ・メッシュで形状
を現すため、各三角形の頂点の3次元座標値とその座標
値の結線情報で構成され、コンピュータ3の記録装置に
格納する。座標系原点は、OMラインを含む面上に置
く。この作成作業は、1被験者に対して1度必要であ
り、その実施例を図3に示す。
【0009】ステップ2(S2 )において、脳実質の立
体形状を作成する。MRI断層写真は頭蓋底から頭頂ま
で、OMラインに平行に撮影した2〜3mm幅の複数枚の
脳断層写真を用いる。夫々の断層写真を、透過型イメー
ジ・スキャナ1によりコンピュータ3に読み込む。読み
込まれた夫々の断層写真データを8ビット/画素のグレ
ー・スケールで表わし、検出しようとする脳実質内側の
境界の画素の平均値を閾値として2値化し、2値化した
画像データの周辺をトレースして境界を座標数値化し、
脳実質の立体形状を作成する。脳実質の立体形状の内容
としては、三角形ワイヤ・メッシュで形状を現すため、
各三角形の頂点の3次元座標値とその座標値の結線情報
で構成され、コンピュータ3の記録装置に格納する。座
標系元点はOMラインを含む面上に置く。この作成作業
は1被験者に1度行う必要があり、その実施例を図4に
示す。
【0010】ステップ(S3 )において、記録電極位置
を座標数値化する。ステップ1(S1 )においてコンピ
ュータ3に読み込んだX線CT写真上の脳硬膜下の記録
電極位置を座標数値化し、コンピュータ3の記録装置に
格納する。図5に脳硬膜下に挿入した記録電極位置を示
す。
【0011】ステップ4(S4 )において、2層構造
(脳脊髄液と脳実質)での伝達行列の計算を行う。境界
要素法を用いて双極子位置の推定計算を行うため、ステ
ップ1(S1 )、ステップ2(S2 )、ステップ3(S
3 )で求めた数値データを使って伝達行列計算を行う
が、その詳細は後述する。なお、この計算は1被験者に
1度行う必要がある。
【0012】ステップ5(S5 )において、記録した電
位に基づいて双極子位置の推定計算を行う。脳内の複数
箇所(基本的には6箇所)に双極子の初期値を置き、そ
の初期値ごとにシンプレックス法を用いて双極子の位置
を移動し、記録電極位置における記録電位値と双極子が
発生する計算電位値の二乗誤差が最小になる双極子の位
置とベクトル成分を求める。推定計算については後述す
る。双極子の初期値の配置を、図6(1双極子)と図7
(2双極子)に示す。因みに、双極子の初期値は位置パ
ラメータのみである。
【0013】ステップ6(S6 )において、双極子の推
定位置を表示する。ステップ5(S5 )において求めた
双極子の位置、電位の大きさとその方向をステップ2
(S2 )によって作成した脳実質の立体形状上にスーパ
ー・インポーズし、表示装置4に表示する。図8に実施
例を示すが、測定電位と計算電位の近似度は別途表示す
る。
【0014】ステップ4(S4 )における2層構造(脳
脊髄液と脳実質)での伝達行列計算を説明すると次の通
りである。 脳実質領域をΩ1 、その導電率をδ1 脳脊髄液領域をΩ2 、その導電率をδ2 脳内起電力が発生する電位をu、脳内起電力が導電率δ
1 の無限一様媒質内で発生する電位を
【0015】
【数1】
【0016】とすると、境界要素法の関係式(1)が成
立する。
【0017】
【数2】
【0018】電位及び電流密度の法線成分の連続性から
(2)式が成立するので、(3)式を定義して(1)式
を(3)式に代入すると、(4)式が得られる。
【0019】
【数3】
【0020】
【数4】
【0021】
【数5】
【0022】(4)式から(5)式が得られるので、
(5)式を(1)式に代入して整理すると(6)式が得
られる。
【0023】
【数6】
【0024】
【数7】
【0025】(6)式における行列HとGを(7)式に
よって表わし、Moore-Penrose の一般行列H+ を使用す
ると(8)式が得られる。
【0026】
【数8】
【0027】
【数9】
【0028】(6)式を(5)式に代入して整理すると
(9)式が得られ、脳髄液表面の電位をucsf とすると
(10)式が得られる。
【0029】
【数10】
【0030】
【数11】
【0031】脳実質領域の導電率δ1 は脳脊髄液領域の
導電率δ2 に比べて極めて小さいとすると、q2 ≒0で
ある。従って、(6)式,(7)式,(8)式,(9)
式より(11)式が得られ、(11)式における伝達行
列は(12)式によって求めることができる。
【0032】
【数12】
【0033】
【数13】
【0034】次に、ステップ5(S5 )における双極子
推定計算について説明する。実測電位をVmes 、双極子
の発生する電位をVcal とし、二乗誤差をSとすると
(13)式によって表わせる。
【0035】
【数14】
【0036】任意の双極子位置に対し、二乗誤差Sを最
小にする双極子ベクトルは(14)式によって表わせ
る。(14)式より、二乗誤差Sは双極子の位置のみの
関数として表わせることが判る。
【0037】
【数15】
【0038】(15)式を最小にする双極子位置は(1
6)式によって求められる。
【0039】
【数16】
【0040】
【数17】
【0041】2双極子の場合におけるVcal は(17)
式によって表わせる。(18)式を最小にする双極子位
置を求めると、(19)式が得られ、(19)式によっ
て2双極子の場合の双極子ベクトルを決定する。
【0042】
【数18】
【0043】
【数19】
【0044】
【数20】
【0045】
【発明の効果】以上説明したように、本発明による脳内
双極子の追跡装置によると、脳硬膜下に挿入して脳波収
集を行う期間を短縮化できるので患者の負担を軽減させ
ることができ、また、脳外科手術箇所を正確に決定でき
る効果がある。
【図面の簡単な説明】
【図1】本発明による双極子位置推定方法を示すフロー
チャート。
【図2】本発明による脳内双極子の追跡装置の構成を示
すブロック図。
【図3】頭蓋骨の内側立体形状。
【図4】脳実質の立体形状。
【図5】脳硬膜下の記録電極位置。
【図6】初期値の配置(1双極子)。
【図7】初期値の配置(2双極子)。
【図8】双極子の表示例。
【符号の説明】
1 透過型イメージ・スキャナ 2 脳波収集装置 3 コンピュータ 4 表示装置
───────────────────────────────────────────────────── フロントページの続き (72)発明者 西条 寿夫 富山県射水郡小杉町南太閤山2−2 富山 医科薬科大学職員宿舎 5−303 (72)発明者 宮本 啓一 東京都八王子市元本郷町1丁目9番9号 中央電子株式会社内

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】 脳硬膜下に挿入した記録電極によって脳
    電位情報を収集する脳波収集装置と、 X線CTによる複数の脳断層写真およびMRIによる複
    数の脳断層写真を読み取り、読み取った前記写真情報を
    コンピュータに格納する透過型イメージ・スキャナと、 前記透過型イメージ・スキャナが読み取った複数のX線
    CT写真情報から構築した頭蓋骨の内側立体形状と、複
    数のMRI写真情報から構築した脳実質の立体形状とを
    作成すると共に、前記X線CT写真情報から3次元座標
    化した記録電極位置の3次元座標値と前記立体形状情報
    とに基づいて脳2層構造での伝達行列計算を行い、さら
    に、前記記録電極によって記録した脳波の1個又は2個
    の等価双極子計算を行い、その計算結果を脳実質の3次
    元立体形状上にスーパー・インポーズするコンピュータ
    と、 前記3次元立体形状上にスーパー・インポーズされた計
    算結果である双極子の位置、電位の大きさと方向、およ
    び測定電位と計算電位との近似度を表示する表示装置
    と、 によって構成したことを特徴とする脳硬膜下の記録脳波
    による脳内双極子の追跡装置。
JP20329896A 1996-07-12 1996-07-12 脳硬膜下の記録脳波による脳内双極子の追跡装置 Expired - Fee Related JP3191858B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP20329896A JP3191858B2 (ja) 1996-07-12 1996-07-12 脳硬膜下の記録脳波による脳内双極子の追跡装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP20329896A JP3191858B2 (ja) 1996-07-12 1996-07-12 脳硬膜下の記録脳波による脳内双極子の追跡装置

Publications (2)

Publication Number Publication Date
JPH1024024A true JPH1024024A (ja) 1998-01-27
JP3191858B2 JP3191858B2 (ja) 2001-07-23

Family

ID=16471722

Family Applications (1)

Application Number Title Priority Date Filing Date
JP20329896A Expired - Fee Related JP3191858B2 (ja) 1996-07-12 1996-07-12 脳硬膜下の記録脳波による脳内双極子の追跡装置

Country Status (1)

Country Link
JP (1) JP3191858B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103505203A (zh) * 2013-09-30 2014-01-15 西安交通大学 基于脑电源定位方法的人体精神状态检测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103505203A (zh) * 2013-09-30 2014-01-15 西安交通大学 基于脑电源定位方法的人体精神状态检测方法

Also Published As

Publication number Publication date
JP3191858B2 (ja) 2001-07-23

Similar Documents

Publication Publication Date Title
Liu et al. Detecting large‐scale networks in the human brain using high‐density electroencephalography
US7778686B2 (en) Method and apparatus for medical intervention procedure planning and location and navigation of an intervention tool
CN106056537A (zh) 一种医学图像拼接方法及装置
US20040087850A1 (en) Method and apparatus for medical intervention procedure planning
JP6175071B2 (ja) 胸部画像の処理及び表示
US8538106B2 (en) Three-dimensional esophageal reconstruction
CN103228214B (zh) 图像处理装置及图像处理方法
US9600895B2 (en) System and method for three-dimensional nerve segmentation using magnetic resonance imaging
JP5658686B2 (ja) 貫壁性のかん流の勾配の画像分析
TWI586326B (zh) 具有顱內電極的個人腦結構之顯示裝置及其顯示方法
US9730631B1 (en) Diagnostic imaging method
US20200411180A1 (en) Biological data processing apparatus, biological data measurement system, and recording medium
WO2019167731A1 (ja) 画像処理装置、画像処理方法およびプログラム
JP4299781B2 (ja) ソース再構成のための信頼区間を表示する方法及びシステム
US8644608B2 (en) Bone imagery segmentation method and apparatus
US20070092124A1 (en) System for and method of displaying subtraction image and computer program for the system
US20200305837A1 (en) System and method for guided ultrasound imaging
JP4051142B2 (ja) 関心領域設定方法、画像処理装置および医用画像処理装置
US8331638B2 (en) Creation of motion compensated MRI M-mode images of the myocardial wall
JPS62327A (ja) 頭蓋内手術部位の位置検出装置
CN112641466A (zh) 一种超声人工智能辅助诊断方法及装置
JP3191858B2 (ja) 脳硬膜下の記録脳波による脳内双極子の追跡装置
JP2021083479A (ja) 生体情報表示装置、生体情報表示方法及び表示プログラム
Ren et al. A knowledge-based automatic cephalometric analysis method
JP3191862B2 (ja) バックグラウンド脳電位の双極子推定位置固定による脳内双極子追跡装置

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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

Free format text: PAYMENT UNTIL: 20090525

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20090525

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20100525

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees