JP2020146204A - 情報処理装置、情報処理方法、プログラムおよび情報処理システム - Google Patents
情報処理装置、情報処理方法、プログラムおよび情報処理システム Download PDFInfo
- Publication number
- JP2020146204A JP2020146204A JP2019045867A JP2019045867A JP2020146204A JP 2020146204 A JP2020146204 A JP 2020146204A JP 2019045867 A JP2019045867 A JP 2019045867A JP 2019045867 A JP2019045867 A JP 2019045867A JP 2020146204 A JP2020146204 A JP 2020146204A
- Authority
- JP
- Japan
- Prior art keywords
- pixel value
- image
- unit
- view
- brain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 230000010365 information processing Effects 0.000 title claims abstract description 76
- 238000003672 processing method Methods 0.000 title claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims abstract description 106
- 230000002093 peripheral effect Effects 0.000 claims 2
- 210000004556 brain Anatomy 0.000 description 222
- 238000010586 diagram Methods 0.000 description 212
- 238000005259 measurement Methods 0.000 description 81
- 230000008859 change Effects 0.000 description 64
- 238000000034 method Methods 0.000 description 28
- 230000004048 modification Effects 0.000 description 28
- 238000012986 modification Methods 0.000 description 28
- 238000003860 storage Methods 0.000 description 19
- 230000000694 effects Effects 0.000 description 16
- 238000004891 communication Methods 0.000 description 15
- 238000009826 distribution Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 15
- 238000004040 coloring Methods 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 8
- 210000000857 visual cortex Anatomy 0.000 description 8
- 230000007177 brain activity Effects 0.000 description 6
- 210000003926 auditory cortex Anatomy 0.000 description 5
- 210000000337 motor cortex Anatomy 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000001028 reflection method Methods 0.000 description 5
- 210000004092 somatosensory cortex Anatomy 0.000 description 5
- 238000001356 surgical procedure Methods 0.000 description 5
- 238000002582 magnetoencephalography Methods 0.000 description 4
- 239000003086 colorant Substances 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 3
- 230000001771 impaired effect Effects 0.000 description 3
- 210000004705 lumbosacral region Anatomy 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 206010015037 epilepsy Diseases 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000014644 Brain disease Diseases 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 210000001787 dendrite Anatomy 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000000537 electroencephalography Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 210000001652 frontal lobe Anatomy 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 206010041232 sneezing Diseases 0.000 description 1
- 230000003238 somatosensory effect Effects 0.000 description 1
- 210000000278 spinal cord Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000005062 synaptic transmission Effects 0.000 description 1
- 210000003478 temporal lobe Anatomy 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/18—Image warping, e.g. rearranging pixels individually
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
【課題】生体画像の任意の位置に、当該生体画像への影響を抑えつつ、視認性を向上させた情報を重畳して表示することができる情報処理装置、情報処理方法、プログラムおよび情報処理システムを提供することを目的とする。【解決手段】生体の解析をするための情報処理装置であって、生体を示す画像である生体画像の画素値に基づいて、解析の結果を示す結果画像に用いる画素値を特定する特定部と、特定部により特定された画素値とした結果画像を、生体画像に重畳させる重畳部と、を備える。【選択図】図3
Description
本発明は、情報処理装置、情報処理方法、プログラムおよび情報処理システムに関する。
脳外科手術等を行うためには、切除対象となる脳の疾患部位であるターゲット部位と、切除せずに温存すべき温存部位とを特定することが必要となる。温存部位としては、視覚野、聴覚野、体性感覚野、運動野および言語野等がある。これらの温存部位を誤って切除してしまうと、人間の感覚および運動等が害されることになるので、脳外科手術等に際してのターゲット部位および温存部位の特定は非常に重要である。このような脳外科手術等のために、事前に脳の活動を調べるためには脳磁計、脳波計、fMRI(functional Magnetic Resonance Imaging)、またはfNIRS(functional Near−Infrared Spectroscopy)等を用いて脳内の物理現象を計測する。このうち、fMRIおよびfNIRSは、脳内を流れる血流を測定して生体信号を得るものであるが、血流という性質上、脳の活動を計測するのに精度的に制限がある。一方、脳磁計は、脳内の電気的活動により生じる磁場を測定し、脳波計は、脳内の電気的活動を測定して波形として生体信号を測定できる。
このような脳磁計、脳波計等による解析で得られた情報等を脳の画像等の生体画像に重畳して表示する技術が開示されている。このような生体画像に情報を表示する技術として、生体画像上の空き領域に文字列の情報を表示させる技術が開示されている(例えば特許文献1)。
しかしながら、例えば脳活動の解析結果を含む情報を重畳するような場合には、脳活動が認められた部位にマーク等を表示する必要がある。特許文献1に記載された技術では、生体画像上の空き領域に情報を表示するものであり、空き領域ではなく、例えば脳画像上に情報を表示させた場合に、当該情報の表示色が不適切な場合には、当該情報の視認性が損なわれるという問題があった。
本発明は、上記に鑑みてなされたものであって、生体画像の任意の位置に、当該生体画像への影響を抑えつつ、視認性を向上させた情報を重畳して表示することができる情報処理装置、情報処理方法、プログラムおよび情報処理システムを提供することを目的とする。
上述した課題を解決し、目的を達成するために、本発明は、生体の解析をするための情報処理装置であって、前記生体を示す画像である生体画像の画素値に基づいて、前記解析の結果を示す結果画像に用いる画素値を特定する特定部と、前記特定部により特定された画素値とした前記結果画像を、前記生体画像に重畳させる重畳部と、を備えたことを特徴とする。
本発明によれば、生体画像の任意の位置に、当該生体画像への影響を抑えつつ、視認性を向上させた情報を重畳して表示することができる。
以下に、図面を参照しながら、本発明に係る情報処理装置、情報処理方法、プログラムおよび情報処理システムの実施形態を詳細に説明する。また、以下の実施形態によって本発明が限定されるものではなく、以下の実施形態における構成要素には、当業者が容易に想到できるもの、実質的に同一のもの、およびいわゆる均等の範囲のものが含まれる。さらに、以下の実施の形態の要旨を逸脱しない範囲で構成要素の種々の省略、置換、変更および組み合わせを行うことができる。
(生体信号計測システムの概略)
図1は、実施形態に係る生体信号計測システムの概略図である。図1を参照しながら、本実施形態に係る生体信号計測システム1の概略について説明する。
図1は、実施形態に係る生体信号計測システムの概略図である。図1を参照しながら、本実施形態に係る生体信号計測システム1の概略について説明する。
生体信号計測システム1(情報処理システムの一例)は、被検者の複数種類の生体信号(例えば、脳磁(MEG:Magneto−encephalography)信号、および脳波(EEG:Electro−encephalography)信号)を計測し、表示するシステムである。なお、測定対象となる生体信号は、脳磁信号および脳波信号に限られるものではなく、例えば、心臓の活動に応じて発生する電気信号(心電図として表現可能な電気信号)であってもよい。図1に示すように、生体信号計測システム1は、被検者の1以上の生体信号を測定する測定装置3と、測定装置3で測定された1種類以上の生体信号を記録するサーバ40と、サーバ40に記録された1種類以上の生体信号を解析する情報処理装置50と、を含む。なお、図1では、サーバ40と情報処理装置50とが別々に記載されているが、例えば、サーバ40が有する機能の少なくとも一部が情報処理装置50に組み込まれる形態であってもよい。
図1の例では、被検者(被測定者)は、頭に脳波測定用の電極(またはセンサ)を付けた状態で測定テーブル4に仰向けで横たわり、測定装置3のデュワ31の窪み32に頭部を入れる。デュワ31は、液体ヘリウムを用いた極低温環境の保持容器であり、デュワ31の窪み32の内側には脳磁測定用の多数の磁気センサが配置されている。測定装置3は、電極からの脳波信号と、磁気センサからの脳磁信号とを収集し、収集した脳波信号および脳磁信号を含むデータ(以下、「測定データ」と称する場合がある)をサーバ40へ出力する。サーバ40へ出力された測定データは、情報処理装置50に読み出されて表示され、解析される。一般的に、磁気センサを内蔵するデュワ31および測定テーブル4は、磁気シールドルーム内に配置されているが、図1では便宜上、磁気シールドルームの図示を省略している。
情報処理装置50は、複数の磁気センサからの脳磁信号の波形と、複数の電極からの脳波信号の波形を、同じ時間軸上に同期させて表示する装置である。脳波信号とは、神経細胞の電気的な活動(シナプス伝達の際のニューロンの樹状突起で起きるイオン電荷の流れ)を電極間の電圧値として表される信号である。脳磁信号とは、脳の電気活動により生じた微小な電場変動を表す信号である。脳磁場は、高感度の超伝導量子干渉計(SQUID:Superconducting Quantum Interference Device)センサで検知される。これらの脳波信号および脳磁信号は、「生体信号」の一例である。
(情報処理装置のハードウェア構成)
図2は、実施形態に係る情報処理装置のハードウェア構成の一例を示す図である。図2を参照しながら、本実施形態に係る情報処理装置50のハードウェア構成について説明する。
図2は、実施形態に係る情報処理装置のハードウェア構成の一例を示す図である。図2を参照しながら、本実施形態に係る情報処理装置50のハードウェア構成について説明する。
図2に示すように、情報処理装置50は、CPU(Central Processing Unit)101と、RAM(Random Access Memory)102と、ROM(Read Only Memory)103と、補助記憶装置104と、ネットワークI/F105と、入力装置106と、表示装置107と、を有し、これらがバス108で相互に接続されている。
CPU101は、情報処理装置50の全体の動作を制御し、各種の情報処理を行う演算装置である。CPU101は、ROM103または補助記憶装置104に記憶された情報表示プログラムを実行して、測定収集画面および解析画面(時間周波数解析画面等)の表示動作を制御する。
RAM102は、CPU101のワークエリアとして用いられ、主要な制御パラメータおよび情報を記憶する揮発性の記憶装置である。ROM103は、基本入出力プログラム等を記憶する不揮発性の記憶装置である。例えば、上述の情報表示プログラムがROM103に記憶されているものとしてもよい。
補助記憶装置104は、HDD(Hard Disk Drive)またはSSD(Solid State Drive)等の記憶装置である。補助記憶装置104は、例えば、情報処理装置50の動作を制御する制御プログラム、ならびに、情報処理装置50の動作に必要な各種のデータおよびファイル等を記憶する。
ネットワークI/F105は、サーバ40等のネットワーク上の機器と通信を行うための通信インターフェースである。ネットワークI/F105は、例えば、TCP(Transmission Control Protocol)/IP(Internet Protocol)に準拠したNIC(Network Interface Card)等によって実現される。
入力装置106は、タッチパネルの入力機能、キーボード、マウスおよび操作ボタン等のユーザインターフェース等である。表示装置107は、各種の情報を表示するディスプレイ装置である。表示装置107は、例えば、タッチパネルの表示機能、液晶ディスプレイ(LCD:Liquid Crystal Display)または有機EL(Electro−Luminescence)等によって実現される。表示装置107は、測定収集画面および解析画面を表示し、入力装置106を介した入出力操作に応じて画面が更新される。
なお、図2に示す情報処理装置50のハードウェア構成は一例であり、これ以外の装置が備えられるものとしてもよい。また、図2に示す情報処理装置50は、例えば、PC(Personal Computer)を想定したハードウェア構成であるが、これに限定されるものではなく、タブレット等のモバイル端末であってもよい。この場合、ネットワークI/F105は、無線通信機能を有する通信インターフェースであればよい。
(情報処理装置の機能ブロック構成)
図3は、実施形態に係る情報処理装置の機能ブロック構成の一例を示す図である。図3を参照しながら、本実施形態に係る情報処理装置50の機能ブロック構成について説明する。
図3は、実施形態に係る情報処理装置の機能ブロック構成の一例を示す図である。図3を参照しながら、本実施形態に係る情報処理装置50の機能ブロック構成について説明する。
図3に示すように、情報処理装置50は、収集表示制御部201と、解析表示制御部202と、ピークリスト制御部203(ピーク制御部)と、通信部204と、センサ情報取得部205と、解析部206(算出部)と、記憶部207と、入力部208と、生成部221と、特定部222と、重畳部223と、出力部224と、を有する。
収集表示制御部201は、センサ情報の収集動作時の画面表示を、図5〜図10を参照して説明する手法で制御する機能部である。
解析表示制御部202は、センサ情報取得部205により取得されたセンサ情報(脳波信号または脳磁信号)から解析部206により算出された生体信号の信号強度等の画面表示を、図11〜図60を参照して説明する手法で制御する機能部である。解析表示制御部202は、図3に示すように、ヒートマップ表示制御部211と、立体表示制御部212と、断面表示制御部213と、再生表示制御部214と、を有する。
ヒートマップ表示制御部211は、後述する図11等に示す時間周波数解析画面601のヒートマップ611の画面表示を制御する機能部である。立体表示制御部212は、時間周波数解析画面601の立体図612の画面表示を制御する機能部である。断面表示制御部213は、時間周波数解析画面601の頭部三面図613の画面表示を制御する機能部である。再生表示制御部214は、時間周波数解析画面601の再生制御パネル615に対する操作入力に従って、再生表示を制御する機能部である。
ピークリスト制御部203は、設定された条件を満たす信号強度のピークを抽出して、後述する図11等に示す時間周波数解析画面601のピークリスト614に登録する機能部である。
通信部204は、測定装置3またはサーバ40等とデータ通信を行う機能部である。通信部204は、図2に示すネットワークI/F105によって実現される。
センサ情報取得部205は、測定装置3またはサーバ40から、通信部204を介して、センサ情報(脳波信号または脳磁信号)を取得する機能部である。解析部206は、センサ情報取得部205により取得されたセンサ情報(測定された信号)を解析し、脳内の各部における信号強度を示す信号(以下、当該信号についても「生体信号」と称する場合がある)を算出する機能部である。
記憶部207は、解析部206により算出された信号強度を示す生体信号のデータ等を記憶する機能部である。記憶部207は、図2に示すRAM102または補助記憶装置104によって実現される。
入力部208は、センサ情報に付加されるアノテーション情報の入力操作、および、時間周波数解析画面601に対する各種入力操作を受け付ける機能部である。入力部208は、図2に示す入力装置106によって実現される。
生成部221は、脳の断面図等の生体画像において、横軸を生体画像の各画素の画素値(濃度)、縦軸をその画素値を有する画素の数とするヒストグラム(第1情報の一例)(後述の図64参照)を生成する機能部である。
特定部222は、生成部221により生成されたヒストグラムにおいて、最も画素の数が少ない画素値を特定する機能部である。
重畳部223は、ピークリスト制御部203により抽出されたピークの位置、またはダイポール推定の位置等を示すマークを、特定部222により特定された画素値のマークにして、脳の断面図等の生体画像に重畳させる機能部である。
出力部224は、重畳部223によりマークが重畳された生体画像をデータ出力する機能部である。出力部224は、例えば、重畳部223によりマークが重畳された生体画像を、DICOM(Digital Imaging and COmmunications in Medicine)形式等のフォーマットのファイルに変換して、補助記憶装置104に記憶させたり、あるいは、通信部204を介してサーバ40へ送信する。なお、出力部224により出力される画像のフォーマットはDICOM形式に限定されるものではない。
上述の収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、重畳部223および出力部224は、CPU101がROM103等に記憶されたプログラムをRAM102に展開して実行することにより実現される。なお、収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、重畳部223および出力部224の一部または全部は、ソフトウェアであるプログラムではなく、ASIC(Application Specific Integrated Circuit)またはFPGA(Field−Programmable Gate Array)等のハードウェア回路によって実現されてもよい。
なお、図3に示した各機能部は、機能を概念的に示したものであって、このような構成に限定されるものではない。例えば、図3で独立した機能部として図示した複数の機能部を、1つの機能部として構成してもよい。一方、図3の1つの機能部が有する機能を複数に分割し、複数の機能部として構成するものとしてもよい。
(開始画面での動作)
図4は、実施形態に係る情報処理装置に表示される開始画面の一例を示す図である。図4を参照しながら、開始画面501での動作について説明する。
図4は、実施形態に係る情報処理装置に表示される開始画面の一例を示す図である。図4を参照しながら、開始画面501での動作について説明する。
開始画面501には、「測定収集」および「解析」の選択ボタンが表示される。脳波測定および脳磁測定の場合、データの測定収集とデータの解析とでは、別々の主体によって行われる場合が多い。例えば、測定技師(測定者)によって「測定収集」のボタンが選択されると、測定装置3で測定されたデータは順次、サーバ40に保存され、情報処理装置50に読み出されて表示される。測定収集の終了後、医師によって「解析」のボタンが選択されると、収集された測定データが読み出されて解析される。
(測定収集時の動作)
図5は、測定収集画面の一例を示す図である。測定収集画面502は、図5に示すように、測定された生体信号(ここでは、脳磁信号および脳波信号)の信号波形を表示する領域511aと、信号波形以外のモニタ情報を表示する領域511bと、を有する。信号波形を表示する領域511aは、測定者からみて画面の左側に配置され、信号波形以外のモニタ情報を表示する領域511bは、測定者からみて画面の右側に配置されている。これによって、リアルタイムで検出され表示される波形の動き(画面の左側から右側に向かって表示される)に合わせた測定者の視線の動きと、画面左側の領域511aから右側の領域511bへマウスを移動させるときの動きに無駄が生じず、作業効率が向上する。
図5は、測定収集画面の一例を示す図である。測定収集画面502は、図5に示すように、測定された生体信号(ここでは、脳磁信号および脳波信号)の信号波形を表示する領域511aと、信号波形以外のモニタ情報を表示する領域511bと、を有する。信号波形を表示する領域511aは、測定者からみて画面の左側に配置され、信号波形以外のモニタ情報を表示する領域511bは、測定者からみて画面の右側に配置されている。これによって、リアルタイムで検出され表示される波形の動き(画面の左側から右側に向かって表示される)に合わせた測定者の視線の動きと、画面左側の領域511aから右側の領域511bへマウスを移動させるときの動きに無駄が生じず、作業効率が向上する。
表示画面の領域511bでは、測定中に被測定者の様子を確認するためのモニタウィンドウ512が表示される。測定中の被測定者のライブ映像を表示することで、後述するように、信号波形のチェックおよび判断の信頼性を高めることができる。図5では、1台のモニタディスプレイ(表示装置107)の表示画面に、測定収集画面502の全体が表示される場合を示しているが、左側の領域511aおよび右側の領域511bを、2台またはそれ以上のモニタディスプレイに分けて別々に表示してもよい。
図6は、測定収集画面の左側の領域の拡大図である。領域511aは、信号検出の時間情報を画面の水平方向に表示する第1表示領域530と、信号検出に基づく複数の信号波形を画面の垂直方向に並列に表示する第2表示領域521〜523と、を有する。
第1表示領域530で表示される時間情報は、図6の例では、時間軸531に沿って付された時間表示を含むタイムラインであるが、時間(数字)を表示せずに、帯状の軸だけでもよく、軸を設けずに時間(数字)の表示だけであってもよい。また、画面上側の第1表示領域530の他に、第2表示領域523の下側に時間軸531を表示して、タイムラインを表示してもよい。
領域511aでは、同種の複数のセンサから取得される複数の信号波形、または複数種類のセンサ群から取得される複数種類の信号の波形が、同じ時間軸531で同期して表示される。図6に示す例では、第2表示領域521には被測定者の頭部右側から得られる複数の脳磁信号の波形が、第2表示領域522には被測定者の頭部左側から得られる複数の脳磁信号の波形が、それぞれ並列に表示されている。第2表示領域523には、複数の脳波信号の波形が並列に表示されている。これらの複数の脳波信号の波形は、各電極間で測定された電圧信号である。これらの複数の信号波形それぞれは、その信号が取得されたセンサの識別番号またはチャネル番号と対応付けられて表示されている。
測定が開始され各センサからの測定情報が収集されると、時間の経過と共に領域511aの各第2表示領域521〜523の左端から右方向に向けて信号波形が表示される。ライン532は、計測の時刻(現在)を示しており、画面の左から右へ向けて移動する。領域511aの右端(時間軸531の右端)まで信号波形が表示されると、その後は画面の左端から右へ向けて徐々に信号波形が消え、消えた位置に新しい信号波形が順次左から右方向に表示され、かつ、ライン532も左端から右へ向けて移動していく。これと共に、水平方向の第1表示領域530でも測定の進行に対応して、時間の経過が時間軸531上に表示される。測定収集は、終了ボタン539が押下されるまで継続される。
実施形態の特徴として、測定者(収集者)がデータ収集中に信号波形上で波形の乱れ、振幅の特異点等に気付いたときに、問題となる箇所または範囲を信号波形上でマークすることができる。マーキングの箇所または範囲は、マウスによるポインタ操作またはクリック操作で指定することができる。指定された箇所または範囲は、第2表示領域521〜523の信号波形上に強調表示されると共に、指定結果が対応する時刻位置または時間範囲で、第1表示領域530の時間軸531に沿って表示される。時間軸531上への表示を含むマーキングの情報は、信号波形データと共に保存される。指定された箇所は、ある時刻に対応し、指定された範囲は、ある時刻を含む一定範囲に対応する。
図6の例では、時刻t1で、第2表示領域523で1以上のチャネルを含む範囲が指定され、マーク523a−1で時刻t1を含んだ時間がハイライトで表示されている。マーク523a−1の表示と関連して第1表示領域530の対応する時刻位置に、指定結果を示すアノテーション530a−1が表示されている。時刻t2で、第2表示領域523で別の波形位置またはその近傍がマークされ、その位置(時刻t2)または近傍の領域(少なくとも時間範囲か複数の波形のいずれか一つが指示される)にマーク523a−2がハイライト表示されている。同時に、第1表示領域530の対応する時刻位置(時間範囲)に、アノテーション530a−2が表示される。ここで、アノテーションとは、あるデータに対して関連する情報を注釈として付与することを示す。本実施形態では、少なくとも指定された時間情報に基づいて注釈として表示されるものであって、少なくとも時間情報に基づく波形が表示される位置と結び付けて、注釈として表示するものである。また、複数のチャネルを表示する場合、対応するチャネル情報と結びつけ、注釈として表示してもよい。
時刻t1で第1表示領域530に追加されたアノテーション530a−1は、一例として、アノテーション識別番号と、波形の属性を示す情報と、を含む。この例では、アノテーション番号「1」と共に、波形の属性を表わすアイコンと「strong spike」(ストロングスパイク)というテキスト情報が表示されている。
時刻t2で、測定者が別の波形箇所またはその近傍領域を指定すると、指定された箇所でマーク523a−2がハイライト表示され、これと共に、第1表示領域530の対応する時刻位置に、アノテーション番号「2」が表示される。さらに、ハイライト表示された箇所に、属性選択のためのポップアップウィンドウ535が表示される。ポップアップウィンドウ535は、種々の属性を選択する選択ボタン535aと、コメントおよび追加情報を入力する入力ボックス535bと、を有する。選択ボタン535aには、波形の属性として「速波(fast activity)」、「眼球運動(eye motion)」、「体動(body motion)」、「スパイク(spike)」等、波形乱れの要因が示されている。測定者は、画面の領域511bのモニタウィンドウ512で被測定者の様子を確認することができるので、波形の乱れの原因を示す属性を適切に選択することができる。例えば、波形にスパイクが生じたときに、てんかんの症状を示すスパイクなのか、被測定者の体動(くしゃみ等)に起因するスパイクなのかを判断することができる。
時刻t1でも同じ操作が行われており、図6では、ポップアップウィンドウ535で「スパイク」の選択ボタン535aが選択され、入力ボックス535bに「strong spike」と入力されたことにより第1表示領域530にアノテーション530a−1が表示されている。このような表示態様により、同じ時間軸531上に多数の信号波形を同期して表示する際に、信号波形の注目箇所または範囲を視認により容易に特定することができ、かつ注目箇所の基本情報を容易に把握することができる。
なお、アノテーション530a−1の一部または全部、例えば、属性アイコンおよびテキスト情報の少なくとも一方を、第2表示領域523の信号波形上のマーク523a−1の近傍にも表示してもよい。ここで、信号波形上へのアノテーションの追加は、波形形状のチェックの妨げになる場合もあり得るので、第2表示領域521〜523の信号波形上にアノテーションを表示させる場合は、表示・非表示を選択可能にしておくことが望ましい。
カウンタボックス538は、スパイクのアノテーションの累積数を表示する。「スパイク」が選択される都度、カウンタボックス538のカウンタ値がインクリメントされ、収録開始から現在(ライン532)までのトータルのスパイク数が一目でわかるようになっている。
図7は、測定収集画面の右側の領域の拡大図である。図7は、図6と同じ時刻(ライン532の時点)での状態を示している。領域511bのモニタウィンドウ512では、頭部を測定装置3に入れて測定テーブル4に横たわっている被測定者の状態のライブ映像が表示される。領域511bでは、第2表示領域521、522、523の信号波形のそれぞれに対応する脳磁分布図541、542、脳波分布図550と、アノテーションリスト560と、が表示される。アノテーションリスト560は、図6の信号波形上でマークされたアノテーションの一覧である。第2表示領域521〜523で信号波形上の位置または範囲が指定されアノテーションが付される都度、対応する情報がアノテーションリスト560に順次追加される。測定収集画面502におけるアノテーションリスト560への追加および表示は、例えば、降順(新しいデータを上に表示)で行われるが、この例に限定されない。アノテーションリスト560の表示を昇順にしても構わないが、第1表示領域530で時間軸531に沿って表示されるアノテーションとの対応関係がわかるように表示する。さらに、表示順序を変更させたり、項目ごとにソートさせることも可能である。
図7の例では、アノテーションリスト560には、アノテーション番号「1」に対応する時刻情報と、付加されたアノテーション情報と、がリストされている。アノテーション情報として、「スパイク」を表わす属性アイコンと、「strong spike」というテキストが記録されている。また、マーク523a−1がハイライト表示された時点で、アノテーション番号「2」に対応する時刻情報がリストされている。ここでは、「アノテーション」は、アノテーション番号と、時刻情報と、アノテーション情報との組であると考えてもよいし、アノテーション情報のみであると考えてもよいし、アノテーション情報と、アノテーション番号または時刻情報との組であると考えてもよい。
また、アノテーションリスト560の近傍に、表示/非表示の選択ボックス560aが配置されている。選択ボックス560aで非表示が選択されると、第2表示領域521〜523で、信号波形上のハイライトマーク以外のアノテーションが非表示にされるが、第1表示領域530の時間軸531に沿ったアノテーションの表示は維持される。これにより、信号波形の視認性を阻害せずにアノテーション情報を認識可能にする。
図8は、アノテーション情報が入力された直後の状態を示す図である。具体的には、図8は、時刻t2でポップアップウィンドウ535の「スパイク」が選択され、「normal spike」(ノーマルスパイク)というテキストが入力された直後の画面を示す。図6で例示したポップアップウィンドウ535で「OK」ボタンが選択されると、ポップアップウィンドウ535が閉じて、図8に示すように第1表示領域530の対応する時刻位置にアノテーション530a−2が表示される。アノテーション番号「2」と対応付けて、「スパイク」を表わす属性アイコンと、「normal spike」のテキスト情報が表示される。これと同時に、カウンタボックス538の値がインクリメントされる。また、ハイライト表示されたマーク523a−2の近傍に、属性アイコン526−2が表示される。この例では、マーク523a−1の近傍にも属性アイコン526−1が表示されているが、上述したように、属性アイコン526−1、526−2の表示、非表示は選択可能である。マーク523a−1と属性アイコン526−1とを含むアノテーションA1、およびマーク523a−2と属性アイコン526−2とを含むアノテーションA2も、アノテーション情報に含まれる。
図9は、更新されたアノテーションリストの図である。測定収集画面502の画面左側の領域511aでマーク523a−2に対応するアノテーションが付加されたことにより、アノテーションリスト560が更新される。アノテーション番号「2」に「normal spike」というメモが追加される。
以下同様に、測定中に領域511aで信号波形上の特定の箇所または範囲が指定される都度、指定箇所が強調表示されると共に、第1表示領域530の時間軸531に沿ってアノテーション情報が表示される。領域511bでは、アノテーションリスト560にアノテーション情報が順次追加される。
アノテーションリスト560および信号波形表示の領域511aにおいて、アノテーション番号の表示は必須ではなく、用いなくてもよい。付加されたアノテーションを識別できる情報であれば任意の情報を識別情報として用いることができる。例えば、属性アイコンと、属性文字列(「strong spike」等)と、時間軸531の近傍での時刻とを対応付けて表示してもよい。さらに、領域511aにファイル番号(図9の「File」の項目に表示される番号)を併記して表示してもよい。
終了ボタン539(図8参照)が選択(押下)され測定が終了すると、第2表示領域521〜523で指定されたハイライト箇所は信号波形に対応付けて保存される。第1表示領域530の対応する時刻位置に表示されたアノテーション情報も、アノテーション番号と時刻とに対応付けて保存される。カウンタボックス538のカウンタ値、および、アノテーションリスト560の内容等の関連情報も保存される。これらの表示情報を保存することで、測定者と解析者とが異なる場合でも、解析者は容易に問題となる箇所を認識し、解析することができる。
図10は、実施形態に係る情報処理装置における測定収集動作の一例を示すフローチャートである。図10を参照しながら、本実形態に係る情報処理装置50での測定収集動作について説明する。
図4に示す開始画面501で「測定収集」が選択されると(ステップS11)、測定が開始され、複数の信号の波形が同じ時間軸に沿って同期して表示される(ステップS12)。ここで「複数の信号波形」という場合は、同一種類の複数のセンサで検知された信号波形、および、異なる種類のセンサのそれぞれで検知された複数の信号波形の両方を含む。この例では、複数の生体信号の波形は、被測定者の頭部右側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の頭部左側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の脳波測定用の電極から得られる脳波信号の波形と、から構成されるが、これに限られるものではない。なお、センサの選択については、左/右のセンサ群にとらわれず、頭頂部、前頭葉、側頭葉等のパートから、任意に選択可能である。例えば、図7等に示す「MEG Window Control 1」で頭頂部のセンサが選択された場合は、「MEG Window Control 2」でそれ以外のセンサが選択された状態となる。
情報処理装置50は、表示されている信号波形上で注目箇所または範囲が指定されたか否かを判断する(ステップS13)。注目箇所または範囲の指定があると(ステップS13:Yes)、信号波形の表示領域(第2表示領域521〜523)に指定箇所を強調表示すると共に、時間軸領域(第1表示領域530)の対応する時刻位置に指定結果を表示する(ステップS14)。指定結果には、指定がなされたこと自体を示す情報、または指定の識別情報が含まれる。時間軸領域への指定結果の表示と同時に、または前後して、アノテーションの入力要求の有無を判断する(ステップS15)。アノテーション入力の要求がある場合は(ステップS15:Yes)、入力されたアノテーション情報を時間軸領域の対応する時刻位置に表示すると共に、アノテーションリストに追加する(ステップS16)。その後、測定終了コマンドが入力されたか否かを判断する(ステップS17)。注目位置(領域)の指定がない場合(ステップS13:No)、および、アノテーションの入力要求がない場合(ステップS15:No)は、ステップS17へ移行して測定終了の判断を行う。測定が終了するまで(ステップS17:Yes)、ステップS13〜S16を繰り返す。
この情報表示方法により、複数のセンサからの信号を収集する際に、信号情報の視認性の高い測定収集画面502が提供される。
(時間周波数解析画面における解析動作)
図11は、時間周波数解析画面の一例を示す図である。図11を参照しながら、情報処理装置50に表示される時間周波数解析画面601における解析動作について説明する。
図11は、時間周波数解析画面の一例を示す図である。図11を参照しながら、情報処理装置50に表示される時間周波数解析画面601における解析動作について説明する。
上述の図4に示す開始画面501において「解析」ボタンが押下されると、解析部206は、上述の測定収集画面502での測定収集動作により収集されたセンサ情報(脳波信号または脳磁信号)を解析し、脳(生体部位の一例、発信源の一例)内の各位置における信号強度を示す生体信号を算出する。信号強度を算出する方法としては、例えば、公知の方法である空間フィルタ法等が挙げられるが、その他の方法を用いるものとしてもよい。解析表示制御部202は、上述の図4に示す開始画面501において、「解析」ボタンが選択されると、表示装置107に図11に示す時間周波数解析画面601を表示させる。時間周波数解析画面601は、図11に示すように、解析画面切替リスト605と、ヒートマップ611と、立体図612と、頭部三面図613と、ピークリスト614と、再生制御パネル615と、を表示する。時間周波数解析画面601を用いた解析および測定の主目的は、視覚野、聴覚野、体性感覚野、運動野および言語野等の人間の生存に欠かせない部位を特定して表示することである。ピークリスト614の右側に表示されているピークリスト設定ボタン614aは、ピークリスト614に登録されるピークの条件を設定するためのウィンドウを表示させるためのボタンである。ピークリスト設定ボタン614aの押下による上述の設定については後述する。また、ヒートマップ611、立体図612、頭部三面図613、ピークリスト614および再生制御パネル615の表示および動作の詳細は、順次後述する。
解析画面切替リスト605は、各種解析画面を選択するためのリストである。解析画面切替リスト605から選択できる解析画面としては、本実施形態の生体信号に対する時間および周波数での解析をするための時間周波数解析画面601の他、生体信号からてんかん等の部位を推定および解析を行うためのダイポール推定を行うための解析画面等がある。本実施形態では時間周波数解析画面601での解析動作について説明する。
<ヒートマップについて>
図13は、ヒートマップにおいて特定の位置が指定された状態の一例を示す図である。図14は、ヒートマップにおいてピークリストの3つのピークの位置を表示した状態の一例を示す図である。図15は、ヒートマップにおいてピークの情報に応じて表示形態を変えた状態の一例を示す図である。図16は、ヒートマップにおいて特定の範囲が指定された状態の一例を示す図である。図17は、ヒートマップにおいて特定の範囲が複数指定された状態の一例を示す図である。図18は、時間周波数解析画面において立体図および頭部三面図が追加された状態の一例を示す図である。図13〜図18を参照しながら、時間周波数解析画面601のヒートマップ611についての動作について説明する。
図13は、ヒートマップにおいて特定の位置が指定された状態の一例を示す図である。図14は、ヒートマップにおいてピークリストの3つのピークの位置を表示した状態の一例を示す図である。図15は、ヒートマップにおいてピークの情報に応じて表示形態を変えた状態の一例を示す図である。図16は、ヒートマップにおいて特定の範囲が指定された状態の一例を示す図である。図17は、ヒートマップにおいて特定の範囲が複数指定された状態の一例を示す図である。図18は、時間周波数解析画面において立体図および頭部三面図が追加された状態の一例を示す図である。図13〜図18を参照しながら、時間周波数解析画面601のヒートマップ611についての動作について説明する。
図11に示すように、ヒートマップ611は、解析部206により算出された脳内の各位置の信号強度を示す生体信号について時間周波数分解を施し、横軸を時刻(トリガ時刻からの時間)、縦軸を周波数、そして、時刻および周波数で特定される生体信号の信号強度の分布を色で表した図である。図11に示す例では、信号強度は、例えば、所定の基準に対する増減で示されている。ここで所定の基準とは、例えば、被験者に何も刺激を与えていない場合の信号強度の平均値を0[%]とする。なお、本実施形態では、0±100[%]で図示しているが、これに限られず、100[%]以上の場合には、例えば200[%]等、レンジを変えて表示してもよい。また、単位[%]に代えて、例えば図12に示すヒートマップ611のように、[db](デシベル)で表記してもよい。また、ヒートマップ611は、例えば、時刻0[ms]に被験者に対して何らかの刺激(物理的に衝撃を与える、腕を動かす、言葉を発する、音を聴かせる等)を与えた場合、その後の時刻では刺激を与えた後の脳の活動状態を示し、時刻0[ms]以前の時刻では、刺激を与える前の脳の活動状態を示すことになる。このヒートマップ611の表示動作は、ヒートマップ表示制御部211によって制御される。
図13に示すように、解析者による入力部208に対する操作入力(クリック操作またはタップ操作)によって、ヒートマップ611上の特定の位置(点)を指定することができる。図13に示すように、ヒートマップ表示制御部211は、指定された位置を、例えば、指定部621のように表示する。なお、指定部621は、白抜きの矩形で示しているが、これに限定されるものではなく、その他の表示態様で示すものとしてもよい。
また、図13に示すヒートマップ611では、入力部208に対する操作入力によって指定された位置を表示しているが、ピークリスト614に登録されたピークのうち選択されたピークの時刻および周波数でのヒートマップおよびそのピーク位置が表示されるものとしてもよい。また、例えば、ピークリスト614で選択されたピークから数えて上位N個のピークの位置がヒートマップ611に表示されるものとしてもよい。図14では、上位3個のピークの位置が表示された例を示す。なお、設定によりピークの位置をどのように表示させるかが決定されてもよく、例えば、上述の他、ピークそのものを表示させない設定、信号強度がM以上のピークを表示させる設定等を切り替えられるものとしてもよい。
また、図15に示すように、ヒートマップ611上に表示された複数のピークについて、当該ピークの属性情報に応じて、表示の形態を変えるものとしてもよい。図15では、表示されたピークごとに番号を付加し、番号の表示部分の色をそれぞれ異なる色となるように変えて表示させた例を示している。
上述のように、ヒートマップ611において特定の位置が指定されると、指定された位置に対応する時刻および周波数に対応する生体信号の信号強度の分布が、例えば図11に示すように、立体図612の脳の図では部位712a−1〜712a−5、712b−1〜712b−5のように、頭部三面図613の脳の図では部位713a−1、713a−2、713b、713c、713dのようにヒートマップ(ヒートマップ611におけるヒートマップとは異なる)として表示される。具体的には、ヒートマップ611において指定された位置に対応する時刻および周波数の生体信号の信号強度が、赤〜青のヒートマップとして表示される。
また、図16に示すように、解析者による入力部208に対するドラッグ操作またはスワイプ操作によって、ヒートマップ611上の特定の範囲を指定することができる。図16に示すように、ヒートマップ表示制御部211は、指定された範囲を、例えば、ドラッグ操作等で定まる面積を有する矩形形状の指定領域622のように表示する。なお、指定領域622は、白抜きの矩形範囲で示しているが、これに限定されるものではなく、丸形状、自由形状その他の表示態様で示すものとしてもよい。
このように、ヒートマップ611において特定の範囲が指定されると、指定された範囲に含まれる時刻および周波数に対応する生体信号の信号強度の平均の分布が、例えば図11に示すように、立体図612の脳の図では部位712a−1〜712a−5、712b−1〜712b−5のように、頭部三面図613の脳の図では部位713a−1、713a−2、713b、713c、713dのようにヒートマップ(ヒートマップ611におけるヒートマップとは異なる)として表示される。
また、図17に示すように、解析者による入力部208に対する追加の操作(例えば、右クリックによるドラッグ操作、または新たなスワイプ操作等)によって、すでに指定された指定領域622の他、追加で指定領域623のように指定することができる。この場合、新たに指定された指定領域623に対応する立体図および頭部三面図として、図18に示すように、それぞれ立体図612aおよび頭部三面図613aが表示される。そして、指定領域623に含まれる時刻および周波数に対応する生体信号の信号強度の平均の分布が、立体図612aおよび頭部三面図613aの脳の図にヒートマップ(ヒートマップ611におけるヒートマップとは異なる)として表示される。ヒートマップ611で複数指定した情報を受け取ると、この情報を受け取った順にこれに対応する立体図612および頭部三面図613が上から下に表示される。図18の例では、指定領域622、623の順に選択された例を示す。このように表示することで、直感的にわかりやすくなる。ここで、ヒートマップ611で複数選択した情報を受け取ると、この情報を受け取った順にこれに対応する立体図612および頭部三面図613が下から上に表示される形態としてもよい。この場合、ヒートマップ611の直下に、最新で指定された選択領域に対応する立体図612および頭部三面図613が表示されるため、ヒートマップ611、立体図612、頭部三面図613への視線の移動を少なくできる。なお、ヒートマップ611上で複数の位置を指定する場合、上述の指定領域622、623のような範囲に限らず、指定部621のように点位置を複数指定することも可能である。このように、ヒートマップ611上で複数の位置(点または領域)が指定されることによって、その指定された時刻および周波数に対応する生体信号の信号強度の分布を比較することができる。
<立体図について>
図19は、時間周波数解析画面における立体図の一例を示す図である。図20は、ヒートマップで指定された位置に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図21は、ヒートマップで指定された範囲に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図22は、立体図に表示された脳がヒートマップ上のどの時刻・周波数に対応するのかを線分で示した状態の一例を示す図である。図23は、立体図に表示された脳がヒートマップ上のどの時刻・周波数の範囲に対応するのかを矩形領域で示した状態の一例を示す図である。図24は、立体図上をドラッグすることによって立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図25は、立体図上のいずれかの脳の表示をクリックすることによって、立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図19〜図25を参照しながら、時間周波数解析画面601の立体図612の基本的な表示動作について説明する。
図19は、時間周波数解析画面における立体図の一例を示す図である。図20は、ヒートマップで指定された位置に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図21は、ヒートマップで指定された範囲に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図22は、立体図に表示された脳がヒートマップ上のどの時刻・周波数に対応するのかを線分で示した状態の一例を示す図である。図23は、立体図に表示された脳がヒートマップ上のどの時刻・周波数の範囲に対応するのかを矩形領域で示した状態の一例を示す図である。図24は、立体図上をドラッグすることによって立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図25は、立体図上のいずれかの脳の表示をクリックすることによって、立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図19〜図25を参照しながら、時間周波数解析画面601の立体図612の基本的な表示動作について説明する。
図19に示すように、立体図612は、所定の視点からの脳の立体画像(3D画像)を表示する図であり、ヒートマップ611上で指定された位置(点または範囲)、または、ピークリスト614で選択されたピークの位置に対応する生体信号の信号強度をヒートマップとして重畳表示する。図19に示すように、立体図612では、同じ行では同じ視点からの脳の立体画像が表示される。図19に示す例では、立体図612の上側の行である表示領域612−1に表示された脳の立体画像は、脳の左側面からの視点での画像であり、下側の行である表示領域612−2に表示された脳の立体画像は、脳の右側面からの視点での画像である。この立体図612の表示動作は、立体表示制御部212によって制御される。
なお、図19に示す立体図612では2つの視点から見た場合の脳の立体図、すなわち、2行で構成された脳の立体図であるが、これに限定されるものではなく、その他の数の行で表示するものとしてもよく、行の数を設定により変更できるものとしてもよい。例えば、脳の言語野の測定であれば、左右の際が重要な情報となるので、脳の左側面および右側面からの2つの視点からの脳の立体図を表示(2行で表示)するものとすればよい。ここで、測定対象と視点の対応の一例を下記の(表1)に示す。測定対象は、測定時に被験者へ与える刺激(刺激装置による付与であり(表1)のNo.1〜4に対応する)や、被験者の動作((表1)のNo.5)を示し、収集時に測定収集画面上で選択する項目である。この測定対象を選択することで、これに対応する視点の脳立体図が表示される。また、視点は、被験者の正面を起点としたときの方向を示す。なお、行数についても別途設定することができるものとしてもよい。また、図19に示す立体図は、(表1)におけるNo.2に相当する。以降の説明では、便宜上、立体図612は、2つの行(2つの視点)で構成されているものとして説明する。
また、図20に示すように、ヒートマップ611上で指定部621が指定されると、立体表示制御部212は、指定部621に対応する時刻を立体図612の表示領域の中心として、前後の時刻に対応した脳における信号強度のヒートマップを立体図612に表示する。図20の例では、ヒートマップ611上で時刻560[ms]が指定されたため、立体図612において560[ms]を中心として、550、555、560、565、570[ms]のように、隣り合う脳の表示時刻の間隔を5[ms]おきで表示している。ただし、隣り合う脳の表示時刻の間隔は、その他10[ms]、25[ms]等、設定で変更できるようにしてもよい。
なお、図21に示すように、ヒートマップ611上で範囲(指定領域622)が指定された場合、当該範囲内の信号強度を平均した信号強度のヒートマップを立体図612に表示するものとすればよい。この場合、選択された時刻の範囲に応じて、立体図612で表示される隣り合う立体画像の時刻を調整設定するものとすればよい。例えば、図21に示すように、指定領域622の範囲が450〜600[ms]の場合、立体図612で表示させる隣り合う立体画像の表示時刻の間隔を150[ms]に設定したとき、立体図612の中央に表示される立体画像は450〜600[ms]に対応する。立体図612の中央の立体画像の左の立体画像は300〜450[ms]に対応する図とし、右の立体画像は、600〜750[ms]に対応する図とすればよい。また、各立体画像に表示されるヒートマップは、各時間幅における平均値である。
次に、図22および図23を参照しながら、ヒートマップ611上の位置または範囲と、立体図612の各立体画像との対応付けについて説明する。このような対応付けとしては、まず、図22に示すように、ヒートマップ611上で特定の点位置である指定部621−1が指定された場合、その指定部621−1の時刻および周波数に対応する脳の立体画像が立体図612に表示される。この場合、その脳の立体画像を中心として前後の時刻に対応した脳の立体画像が表示されるので(図22に示す例では5つの時刻に対応する脳)、ヒートマップ表示制御部211は、それぞれの脳の時刻に対応する点を、ヒートマップ611上で対応部621−2〜621−5としてそれぞれ表示させる。この場合、対応部621−2〜621−5の周波数位置は、指定部621−1の周波数位置に合わせる。さらに、ヒートマップ表示制御部211は、図22に示すように、ヒートマップ611上の指定部621−1および対応部621−2〜621−5と、立体図612上の対応する脳の立体画像とをそれぞれ結ぶ線分631−1〜631−5で表示する。これによって、立体図612上で表示されている脳の状態が、ヒートマップ611上のどの位置に対応する状態であるのかを一目で確認することができる。なお、図22に示す例では、線分を用いているが、その他の対応付けの方法として、指定部621−1および対応部621−2〜621−5のマークと、立体図612の脳の背景色とを合わせる等、対応付けが確認できる方法であればよい。また、この場合、解析者が指定した指定部621−1は、対応部621−2〜621−5と区別する態様で表示するものとすればよい。
また、ヒートマップ611上で特定の範囲である指定領域622−1が指定された場合、まず、図23に示すように、その指定領域622−1の時刻および周波数に対応する脳の立体画像が立体図612に表示される。この場合、その脳の立体画像を中心として前後の時刻範囲に対応した脳の立体画像が表示されるので(図23に示す例では5つの時刻範囲に対応する脳)、ヒートマップ表示制御部211は、それぞれの脳の時刻範囲に対応する範囲を、ヒートマップ611上で対応領域622−2〜622−5としてそれぞれ表示させる。また、この場合、解析者が指定した指定領域622−1は、対応領域622−2〜622−5と区別する態様で表示するものとすればよく、例えば、指定領域622−1の矩形の枠の色だけ異なる色となるように表示するものとすればよい。さらに、立体表示制御部212は、図23に示すように、ヒートマップ611上の指定領域622−1および対応領域622−2〜622−5と同様な矩形を、立体図612上で対応する脳の立体画像を囲うように表示させる。これによって、立体図612上で表示されている脳の状態が、ヒートマップ611上のどの範囲に対応する状態であるのかを一目で確認することができる。なお、ヒートマップ611上で特定の範囲である指定領域622−1が指定された場合、指定領域622−1および対応領域622−2〜622−5の枠を表示すると共に、立体図612上で、枠722−1〜722−5およびヒートマップを表示するものとしてもよい。
次に、図24および図25を参照しながら、立体図612上でドラッグ操作、スワイプ操作またはカーソルキー操作をした場合に、立体図612の表示を左右に移動させる動作について説明する。図24では、立体図612上でドラッグ操作、スワイプ操作またはカーソルキー操作により、立体図612の立体画像を右側に移動させた状態の例を示す。この場合、図24に示すように、移動した結果、表示されている脳に応じて時刻表示も更新され、立体図612で中央に表示された脳の立体図が選択された状態となるように矩形が表示される。さらに、立体表示制御部212は、立体図612における立体画像に移動に従って、ヒートマップ611上の指定領域622−1および対応領域622−2〜622−5の表示も移動させる。
また、図25に示すように、解析者による操作によって、立体図612上の中央の立体画像以外の立体画像がクリック操作またはタップ操作されると、操作された脳の立体画像が立体図612の中心に移動する。なお、実装上は脳の画像はそのままで、オーバーラップするヒートマップのみが移動するものとしてもよい。この場合、図25に示すように、移動した結果、表示されている脳に応じて時刻表示も更新され、立体図612で中央に表示された脳の立体図が選択された状態となるように矩形が表示される。さらに、立体表示制御部212は、立体図612における立体画像に移動に従って、ヒートマップ611上の指定領域622−1および対応領域622−2〜622−5の表示も移動させる。
以上のように、立体図612における表示を自由に移動することが可能となることにより、時間の前後で脳の状態の変化を素早く確認することができる。
図26は、立体図に表示されたいずれかの脳の視点を変更した場合、同じ行のすべての脳の視点を変更する状態の一例を示す図である。図27は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の一例を示す図である。図28は、図27における視点の変更動作を具体的に説明する図である。図29は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の別の例を示す図である。図30は、図29における視点の変更動作を具体的に説明する図である。図31は、立体図にコメントを付加した状態の一例を示す図である。図26〜図31を参照しながら、時間周波数解析画面601の立体図612の特定の立体画像について視点変更を行った場合の動作について説明する。
立体図612に立体画像として表示されている脳は、解析者による操作(例えば、ドラッグ操作またはスワイプ操作等)によって視点を変更することが可能である。この場合、立体図612の特定の脳の立体画像の視点を変更した場合において、他の立体画像への反映方法の類型について以下に説明する。
まず、特定の立体画像の視点が変更された場合に、同じ行の他の立体画像について同じ視点変更がなされる場合について説明する。図26(a)に示すように、解析者は、2行で表示された立体図612において、表示領域612−1に表示された脳の立体画像のうち右端の立体画像(以下、「対象立体画像」と称する場合がある)に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部212は、図26(b)に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の背面の視点の立体画像の表示となるように視点変更する。このとき、脳の画像に重畳しているヒートマップも同様に視点変更する。そして、立体表示制御部212は、図26(c)に示すように、対象立体画像と同じ行(表示領域612−1)の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。これによって、特定の立体画像(対象立体画像)に対する視点変更を行うのみで、同じ行の他の立体画像についてもその視点変更が反映されるので、操作性が向上し、同じ視点での時間的に前後の脳の活動変化を確認することが容易になる。なお、立体画像の視点変更の操作として、例えば、視点変更したい立体画像上にマウスを合わせてドラック操作またはクリック操作等でもよく、ポップアップ表示によるパラメータの指定等でもよい。
次に、特定の立体画像の視点が変更された場合に、他のすべての立体画像について同じ視点変更がなされる場合について説明する。図27(a)に示すように、解析者は、2行で表示された立体図612において、表示領域612−1に表示された脳の立体画像のうち右端の対象立体画像に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部212は、図27(b)に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の背面の視点の立体画像の表示となるように視点変更する。そして、立体表示制御部212は、図27(c)に示すように、対象立体画像と同じ行(表示領域612−1)の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図28(a)に示す脳の左側面の視点を、表示領域612−1の他の立体画像について、図28(b)に示すように、脳の背面の視点となるように変更する。さらに、立体表示制御部212は、図27(c)に示すように、対象立体画像と異なる行(表示領域612−2)の立体画像に対しても、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図28(a)に示す脳の右側面の視点を、表示領域612−2の立体画像について、図28(c)に示すように、脳の前面の視点となるように変更する。なお、処理環境に余裕があれば図28(a)〜図28(c)の処理を高速に行い、見た目上すべての脳の視点が同時に変わっているようにしてもよいし、処理環境が貧弱な場合には、ユーザが動かしている画像のみについて視点を変えた後、それが確定したタイミング(例えばドラッグ操作で脳画像を回転させることで視点を変更している場合、マウスのボタンを離したタイミング)で他の画像の視点を変えてもよい。このとき、脳の画像に重畳しているヒートマップもそれぞれ同様に視点を変更する。これによって、特定の立体画像(対象立体画像)に対する視点変更を行うのみで、同じ行の他の立体画像、および他の行の立体画像についてもその視点変更が反映されるので、操作性が向上し、時間的に前後の脳の活動変化を確認することが容易になる。
そして、特定の立体画像の視点が変更された場合に、同じ行の他の立体画像については同じ視点変更がなされ、他の行の立体画像については、対応する視点変更(具体的には、脳の中心面(対象面)に対して面対称となる視点変更)がなされる場合について説明する。図29(a)に示すように、解析者は、2行で表示された立体図612において、表示領域612−1に表示された脳の立体画像のうち右端の対象立体画像に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部212は、図29(b)に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の左前側の視点の立体画像の表示となるように視点変更する。そして、立体表示制御部212は、図29(c)に示すように、対象立体画像と同じ行(表示領域612−1)の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図30(a)に示す脳の左側面の視点を、表示領域612−1の他の立体画像について、図30(b)に示すように、脳の左前側の視点となるように変更する。さらに、立体表示制御部212は、図29(c)に示すように、対象立体画像と異なる行(表示領域612−2)の立体画像に対して、対象立体画像の視点変更に対応する視点変更をして表示する。すなわち、図30(a)に示す脳の右側面の視点を、表示領域612−2の立体画像について、図30(c)に示すように、脳の中心面(対象面)に対して面対称となるような視点変更、すなわち、脳の右前側の視点となるように変更する。このとき、脳の画像に重畳しているヒートマップもそれぞれ同様に視点を変更する。これによって、特定の立体画像(対象立体画像)に対する視点変更を行うのみで、同じ行の他の立体画像についてはその視点変更が反映され、他の行の立体画像については、対応する視点変更が反映されるので、操作性が向上し、さらに、複数の行を見比べることによって、対応する視点から時間的に前後の脳の活動変化を確認することができる。
以上で説明した3つの他の立体画像への反映方法のうち、いずれの方法を採用してもよく、または、これらの方法のうちどの方法で反映するかを設定により切り替えることができるものとしてもよい。
なお、図26〜図30で上述したように、解析者が最初に視点変更の操作の対象とする対象立体画像を、表示領域612−1の右端の立体画像としたが、これに限定されるものではなく、表示領域612−1または表示領域612−2のいずれの立体画像が操作の対象であってもよい。
また、立体図612の特定の脳の立体画像の視点を変更した場合に、他の立体画像へ視点変更を反映する動作を説明したが、立体画像に対して変更を加える表示態様について視点だけではなく、例えば、拡大縮小、輝度の変更、または透明度の変更等であってもよい。これらの変更についても、上述の視点変更の趣旨を逸脱しない範囲で、他の立体画像に変更を反映するものとすればよい。
また、上述のように、立体図612に表示される立体画像に対して変更を加えた結果、図31に示すように、解析者による入力部208への操作入力により、特定の立体画像に対してメモ(例えば、図31に示すコメント635)を付加できるようにしてもよい。これによって、解析者(医者等)が脳の気になる活動部位に対するコメントをその立体画像に関連付けて保存しておくことができ、脳に対する外手術の他、当該脳の疾患についてのカンファレンスにおいても役立てることができる。
<頭部三面図について>
図32は、時間周波数解析画面における頭部三面図の一例を示す図である。図33は、頭部三面図の立体画像としてカットモデルを表示した状態の一例を示す図である。図34は、頭部三面図においてピークリストで選択されたピークの位置を示した状態の一例を示す図である。図35は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を示した状態の一例を示す図である。図36は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を色を変えて示した状態の一例を示す図である。図37は、頭部三面図の立体画像にダイポール推定の結果を重畳表示した状態の一例を示す図である。図38は、頭部三面図の立体画像に複数の測定対象の結果(ヒートマップ)を重畳表示した状態の一例を示す図である。図32〜図38を参照しながら、時間周波数解析画面601の頭部三面図613の基本的な表示動作について説明する。
図32は、時間周波数解析画面における頭部三面図の一例を示す図である。図33は、頭部三面図の立体画像としてカットモデルを表示した状態の一例を示す図である。図34は、頭部三面図においてピークリストで選択されたピークの位置を示した状態の一例を示す図である。図35は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を示した状態の一例を示す図である。図36は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を色を変えて示した状態の一例を示す図である。図37は、頭部三面図の立体画像にダイポール推定の結果を重畳表示した状態の一例を示す図である。図38は、頭部三面図の立体画像に複数の測定対象の結果(ヒートマップ)を重畳表示した状態の一例を示す図である。図32〜図38を参照しながら、時間周波数解析画面601の頭部三面図613の基本的な表示動作について説明する。
図32に示すように、頭部三面図613は、脳の特定の位置(点)における3方向の断面図(以下、「三面図」と総称する場合がある)、および3D画像である立体画像644を含む図である。図32に示す例では、頭部三面図613は、脳の特定の位置における3方向の断面図として、脳の前後方向に対して垂直な断面を示す断面図641、脳の左右方向に対して垂直な断面を示す断面図642、および、脳の上下方向に対して垂直な断面を示す断面図643を含む。断面図641は、上述の特定の位置を通るように基準線645aおよび基準線645bが引かれている。また、断面図642は、上述の特定の位置を通るように、基準線645aおよび基準線645cが引かれている。また、断面図643は、上述の特定の位置を通るように、基準線645bおよび基準線645dが引かれている。断面図641〜643には、それぞれ、ヒートマップ611で指定された位置(点または範囲)に対応する時刻および周波数の生体信号の信号強度の分布を示すヒートマップ(ヒートマップ611とは異なる)が重畳して表示される。この頭部三面図613の表示動作は、断面表示制御部213によって制御される。
基準線645aは、上述の脳の特定の位置の上下方向の位置を規定する線であるため、断面図641および断面図642にわたって連続した線として描画されている。また、基準線645bは、上述の脳の特定の位置の左右方向の位置を規定する線であるため、断面図641および断面図643にわたって連続した線として描画されている。また、基準線645cは、断面図642において、上述の脳の特定の位置の前後方向の位置を規定する線である。また、基準線645dは、断面図643において、上述の脳の特定の位置の前後方向の位置を規定する線である。なお、頭部三面図613における断面図641〜643の並びは、上述のように、基準線645aおよび基準線645bが複数の断面図に連続して描画することが可能であるために、図32に示す並びとしているが、これに限定されるものではなく、任意でよい。その場合、各断面図において、脳の特定の位置を通るように基準線がそれぞれ描画されるものとすればよい。また、基準線は必ずしも描画される必要はなく、各断面図において、脳の特定の位置を示すマーク等が表示されるものとしてもよい。
立体画像644は、脳の3D画像であり、後述するように、この立体画像644に対する操作に応じて、立体図612に描画された脳の立体画像の視点が変更される。また、立体画像644には、ヒートマップ611で指定された位置(点または範囲)に対応する時刻および周波数の生体信号の信号強度の分布を示すヒートマップ(ヒートマップ611とは異なる)が重畳して表示される。なお、立体画像644は、脳の特定の視点の立体画像を表示することに限定されるものではなく、例えば、図33に示すように、三面図で特定される脳の位置を中心に3次元方向で一部カットしたカットモデル画像であってもよい。
また、図32に示す頭部三面図613では、ピークリスト614に登録されたピークのうち選択されたピークの位置が特定される三面図を表示し、図34に示すように、立体画像644上に当該選択されたピーク位置を示すピーク点646を表示するものとしてもよい。また、例えば、ピークリスト614で選択されたピークから数えて上位N個のピークの位置が立体画像644に表示されるものとしてもよい。図35では、上位3個のピークの位置(ピーク点646、646a、646b)が表示された例を示す。また、図35では、上位3個ではなく、ピークリスト614で選択されたピークの前後の時刻におけるピーク位置をピーク点646、646a、646bとして表示(ピークの軌跡を表示)するものとしてもよい。なお、設定によりピークの位置をどのように表示させるかが決定されてもよく、例えば、上述の他、ピークそのものを表示させない設定、信号強度がM以上のピークを表示させる設定等を切り替えられるものとしてもよい。
また、図36に示すように、立体画像644上に表示された複数のピークについて、当該ピークの属性情報に応じて、表示の形態を変えるものとしてもよい。図36では、表示されたピークごとにマークの色をそれぞれ異なる色となるように変えて表示させた例を示している。
また、図37に示すように、断面表示制御部213は、立体画像644に対して、他の解析画面等でダイポール推定の結果であるダイポール647を重畳表示させるものとしてもよい。これによって、立体画像644が示す温存部位を示すヒートマップと、疾患部位(ターゲット部位)を示すダイポールとの位置関係を把握することができ、外科手術等に役立てることができる。
また、三面図のいずれかの断面図上で、解析者による入力部208を用いてクリック操作またはタップ操作によって、脳の三次元空間の特定の位置を指定することができる。このように、三面図において特定の位置が指定されると、指定された位置(点)に対応する時刻および周波数に関する生体信号の信号強度の分布がヒートマップ611に反映される。
また、三面図のいずれかの断面図上で、解析者による入力部208に対するドラッグ操作またはスワイプ操作によって、脳の三次元空間の特定の範囲を指定することができる。このように、三面図において特定の範囲が指定されると、指定された範囲に対応する時刻および周波数に関する生体信号の信号強度(その範囲での平均)の分布がヒートマップ611に反映される。
なお、立体画像644(および三面図)に描画されるヒートマップ(信号強度の強弱を示す等高線図)は、各脳の部位の活動が活発となるようなそれぞれの刺激の結果を重畳して表示するものとしてもよい。例えば、言語刺激(測定時)を行った信号に対する結果と、視覚刺激(測定時)を行った信号に対する結果を求めた後に、断面表示制御部213は、図38(a)に示す言語野を活性化させた場合のヒートマップと、図38(b)に示す視覚野を活性化させた場合のヒートマップとを、図38(c)に示すように、重畳させて表示可能とするものとしてもよい。このように図38(c)に示すように重畳されたヒートマップで示される部位が温存部位であるということが確認できる。重畳させるための操作方法としては、現在表示されている測定結果が言語野であるものとした場合、メニューによって異なる測定結果(例えば、視覚野)を選択することができるようにすればよい。なお、重畳する際には、測定対象より刺激に対する反応時間が異なる可能性もある。そこで、測定対象を追加する場合に時間のずれを設定できるようにすれば、より的確に重畳することができる。さらに、ヒートマップを重畳した結果である図38(c)に示す立体画像を、図38(d)に示すように反転表示させることによって、逆に温存部位ではない切除可能な部位を示すことが可能になる。
また、頭部三面図613に示す断面図は、断面方向の異なる3つの断面図である三面図としているが、これに限定されるものではなく、特定の断面方向の1つの断面図、または断面方向の異なる2つもしくは4つ以上の断面図であってもよい。
図39は、頭部三面図の立体画像に対して視点が変更される前の状態の一例を示す図である。図40は、頭部三面図の立体画像に対して視点が変更された場合に表示されるダイアログボックスを示す図である。図41は、立体画像の視点変更により立体図の1行目に同じ表示を反映するための設定例を示す図である。図42は、頭部三面図の立体画像に対して視点が変更された場合、立体図の1行目に同じ表示を反映した状態を示す図である。図43は、立体画像の視点変更により立体図の1、2行目に同じ変換を反映するための設定例を示す図である。図44は、頭部三面図の立体画像に対して視点が変更された場合、立体図の1、2行目に同じ変換を反映した状態を示す図である。図45は、立体画像の視点変更により立体図の1、2行目に、対応した変換を反映するための設定例を示す図である。図46は、頭部三面図の立体画像に対して視点が変更された場合、立体図の1、2行目に、対応した変換を反映した状態を示す図である。図47は、立体画像の視点変更により立体図に同じ表示の新たな行を追加するための設定例を示す図である。図48は、頭部三面図の立体画像に対して視点が変更された場合、立体図に同じ表示の新たな行を追加した状態を示す図である。図39〜図48を参照しながら、時間周波数解析画面601の頭部三面図613の立体画像に対して視点変更した場合における立体図612への反映動作について説明する。
頭部三面図613の立体画像644に表示されている脳は、立体図612と同様に、解析者による操作(例えば、ドラッグ操作またはスワイプ操作等)によって視点を変更することが可能である。この場合、立体画像644の脳の視点を変更した場合において、立体図612に表示された脳の視点に反映するものとしてもよく、以下、その反映方法の類型について以下に説明する。
図39に示すような時間周波数解析画面601に表示された頭部三面図613の立体画像644に対して、解析者による操作(例えば、ドラッグ操作またはスワイプ操作等)が行われると、断面表示制御部213は、図40に示すように、ダイアログボックス650を表示させる。ダイアログボックス650は、立体画像644の脳について視点変更された場合に、どのように立体図612に反映するかを決定するための画面である。ここで、例えば、「立体図は変更しない」ボタンが押下されると、立体図612の立体画像に対する視点変更はなされない。ここでは、解析者は、図40に示すように、脳の左側面の視点の立体画像644を、脳の背面の視点の立体画像の表示となるように視点変更したものとする。
まず、図41に示すように、ダイアログボックス650の「立体図の行の表示に反映する」ボタンが押下されたものとする。すると、断面表示制御部213は、図41に示すように、立体図612にどのように反映するかの詳細を設定するためのダイアログボックス651を表示させる。解析者は、図41に示すように、ダイアログボックス651において、変更する行として立体図612の1行目に選択し、「立体図の3Dと同じ表示」を選択したものとする。この場合、立体表示制御部212は、図42に示すように、立体図612の1行目(上の行)の立体画像に対して、立体画像644の変更された視点と同じ視点となるように表示する。
次に、図40のように視点変更された状態から、解析者は、図43に示すように、ダイアログボックス650の「立体図の行の表示に反映する」ボタンを押下し、ダイアログボックス651において、変更する行として立体図612の1行目および2行目に選択し、「立体図の3Dと同じ変換」を選択したものとする。この場合、立体表示制御部212は、図44に示すように、もともと立体画像644と同じ視点であった立体図612の1行目の立体画像に対して、立体画像644の視点変更と同じ視点変更をして表示する。すなわち、図44に示すように、脳の背面の視点となるように変更する。さらに、立体表示制御部212は、図44に示すように、もともと脳の右側面の視点であった立体図612の2行目の立体画像に対して、立体画像644の視点変更と同じ視点変更をして表示する。すなわち、図44に示すように、脳の前面の視点となるように変更する。なお、「立体図の行の表示に反映する」ボタン押下によるダイアログボックス651での選択を初期設定済みとするか初期設定可能なようにしておき、その上で、例えば、[ビューのリンク]/[ビューのリンク解除]ボタンを設けることで、その選択結果を表示するようにしてもよい。このような構成とすることで毎回の選択動作を省略(簡略化)できる。
次に、解析者は、図46に示すように、脳の左側面の視点の立体画像644を、脳の左前側の視点の立体画像の表示となるように視点変更したものとする。そして、この状態から、解析者は、図45に示すように、ダイアログボックス650の「立体図の行の表示に反映する」ボタンを押下し、ダイアログボックス651において、変更する行として立体図612の1行目および2行目に選択し、「立体図の3Dに対応した変換」を選択したものとする。この場合、立体表示制御部212は、図46に示すように、もともと立体画像644と同じ視点であった立体図612の1行目の立体画像に対して、立体画像644の視点変更と同じ視点変更をして表示する。すなわち、図46に示すように、脳の左前側の視点となるように変更する。さらに、立体表示制御部212は、図46に示すように、もともと脳の右側面の視点であった立体図612の2行目の立体画像に対して、立体画像644の視点変更に対応する視点変更をして表示する。すなわち、図46に示すように、脳の中心面(対象面)に対して面対称となるような視点変更、すなわち、脳の右前側の視点となるように変更する。
次に、図40のように視点変更された状態から、解析者は、図47に示すように、ダイアログボックス650の「立体図の行を追加する」ボタンを押下することによって表示されるダイアログボックス652において、「立体図の3Dと同じ表示」を選択したものとする。この場合、立体表示制御部212は、図48に示すように、立体画像644の視点変更の結果、同一の視点となる脳の立体画像を、立体図612の表示領域612−3に新たな行として表示する。すなわち、図48に示すように、表示領域612−3における新たな行に、脳の背面の視点の立体画像が表示される。
以上のように、頭部三面図613における立体画像644に対して行った視点変更を、各種設定に応じて、立体図612の時系列に並んだ脳の立体画像の視点に反映することができる。これによって、立体画像644の視点変更と同様の視点変更を、立体図612に対して改めて行う必要がないので、操作性が向上し、さらに、立体画像644で変更した視点と同一の視点、または対応する視点で、立体図612において脳の状態の変化を時系列に確認することができる。
なお、図41、図43、図45および図47で示すダイアログボックス650〜652によって設定される立体図612に表示された脳の視点への反映方法は、一例であり、その他の反映方法が設定できるものとしてもよい。
また、立体画像644の視点を変更した場合に、立体図612の立体画像へ視点変更を反映する動作を説明したが、立体画像644に対して変更を加える表示態様は視点だけではなく、例えば、拡大縮小、輝度の変更、または透明度の変更等であってもよい。これらの変更についても、上述の視点変更の趣旨を逸脱しない範囲で、立体図612の立体画像に変更を反映するものとすればよい。
<ピークリストについて>
図49は、ピークリストの設定例を示す図である。図50は、空間的なピークを説明する図である。図51は、時間・周波数的なピークを説明する図である。図52は、プルダウンされたピークリストから特定のピークを選択する状態を示す図である。図53は、プルダウンされたピークリストから選択されたピークについてヒートマップ、立体図および頭部三面図に反映させた状態を示す図である。図49〜図53を参照しながら、時間周波数解析画面601のピークリスト614の基本的な動作について説明する。
図49は、ピークリストの設定例を示す図である。図50は、空間的なピークを説明する図である。図51は、時間・周波数的なピークを説明する図である。図52は、プルダウンされたピークリストから特定のピークを選択する状態を示す図である。図53は、プルダウンされたピークリストから選択されたピークについてヒートマップ、立体図および頭部三面図に反映させた状態を示す図である。図49〜図53を参照しながら、時間周波数解析画面601のピークリスト614の基本的な動作について説明する。
ピークリスト614は、ピークリスト制御部203により抽出された、設定された条件を満たす信号強度のピークが登録されたリストである。図49に示すように、ピークリスト制御部203は、ピークリスト614でプルダウンされることによって登録された信号強度の一覧であるプルダウンリスト656を表示する。
また、上述のピークリスト制御部203により抽出される信号強度のピークの条件は、ピークリスト設定ボタン614aを押下することにより設定することができる。ピークリスト制御部203は、ピークリスト設定ボタン614aが押下されると、抽出される信号強度のピークの条件を設定するためのダイアログボックス655を表示させる。
ダイアログボックス655では、まず、ピークリスト614に登録されたピーク情報をどのようにソートするかを設定することができる。ダイアログボックス655において、「ピークの値が大きい順番」が選択された場合、ピークリスト制御部203は、ピークリスト614で登録されているピーク情報を、ピークの信号強度が大きい順にソートする。一方、ダイアログボックス655において、「ピークの高さ(頂点と谷点の差)が大きい順番」が選択された場合、ピークリスト制御部203は、ピークリスト614で登録されているピーク情報を、ピーク点の信号強度と、当該ピークの谷部分の信号強度との差が大きい順にソートする。
さらに、ダイアログボックス655では、ピークリスト614にどのようなピーク情報を登録(リストアップ)させるかを設定することができる。ダイアログボックス655において、「すべての空間的なピーク」が選択された場合、ピークリスト制御部203は、時刻・周波数平面の各時刻・各周波数においての脳全体における空間的なピークを抽出して、ピークリスト614に登録する。ここで、空間的なピークとは、図50に示すピーク部801のように、着目する時刻・周波数の生体信号の、脳全体における信号強度のピークであり、このピーク部801の信号強度は、周辺より信号強度が強い。
また、ダイアログボックス655において、「すべての時刻/周波数的なピーク」が選択された場合、ピークリスト制御部203は、脳全体の各位置においての時刻・周波数平面における時刻/周波数的なピークを抽出して、ピークリスト614に登録する。ここで、時刻/周波数的なピークとは、図51に示すピーク部802のように、着目する脳の位置での生体信号の、時刻・周波数平面における信号強度のピークであり、このピーク部802の信号強度は、周辺よりも信号強度が強い。
また、ダイアログボックス655において、「指定されている時刻/周波数における空間的なピーク」が選択された場合、ピークリスト制御部203は、時刻・周波数平面で指定された時刻・周波数においての脳全体における空間的なピークを抽出して、ピークリスト614に登録する。なお、指定されている時刻/周波数は一点とは限らず、範囲で選択される場合もある。
また、ダイアログボックス655において、「指定されている位置における時刻/周波数的なピーク」が選択された場合、ピークリスト制御部203は、指定された脳の位置においての時刻・周波数平面における時刻/周波数的なピークを抽出して、ピークリスト614に登録する。なお、指定されている位置は一点とは限らず、範囲で選択される場合もある。例えば、視覚野についてのピークを抽出する場合、後頭部全体の範囲を指定することによって、所望のピークを抽出しやすくなる。
次に、各ピーク情報が登録されたピークリスト614から特定のピーク情報を選択した場合の動作について説明する。解析者によって、ピークリスト614から表示されたプルダウンリスト656から特定のピーク情報(例えば、図52に示す「95%/9ms/70Hz/voxel:1736」)が選択されると、ヒートマップ表示制御部211は、当該ピーク情報が示す脳の特定の位置に対応するヒートマップ611を表示する。この場合、ヒートマップ表示制御部211は、図14で上述したように、ヒートマップ611においてピーク情報が示すピークの位置を具体的に示すものとすればよい。
また、立体表示制御部212は、選択されたピーク情報が示す時刻・周波数の脳の立体画像を立体図612の各行の中央に表示させ、さらに、その時刻の前後の脳の立体画像を表示させる。この場合、立体図612の脳の各立体画像に重畳されるヒートマップは、ピーク情報が示す周波数の生体信号の信号強度に対応するものとすればよい。
また、断面表示制御部213は、選択されたピーク情報が示す脳の位置を通る三面図を頭部三面図613に表示させる。さらに、断面表示制御部213は、選択されたピーク情報が示す時刻・周波数の生体信号の信号強度に対応するヒートマップを立体画像644の脳に重畳表示させるものとすればよい。なお、断面表示制御部213は、図53に示すように、選択されたピーク情報が示す脳の位置を中心に3次元方向で一部カットしたカットモデルを立体画像644に表示させるものとしてもよい。
以上のように、ピークリスト614に登録されたピーク情報から特定のピーク情報を選択することによって、当該ピーク情報に対応したヒートマップ611、立体図612および頭部三面図613が表示される。これによって、選択したピークが脳のどの位置、およびどの時刻・周波数であるかを瞬時に認識することができ、さらに、ヒートマップ611においては当該ピークおよびその周りの時刻/周波数における信号強度の状態、および、当該ピークの位置およびその周りの脳の信号強度の状態も把握することが可能となる。
<再生制御パネルについて>
図54は、再生制御パネルの操作によりヒートマップおよび立体図が再生表示される状態を示す図である。図55は、再生制御パネルの操作によりヒートマップおよび立体図がコマ戻しされる状態を示す図である。図56は、再生制御パネルの操作によりヒートマップおよび立体図がコマ送りされる状態を示す図である。図54〜図56を参照しながら、時間周波数解析画面601の再生制御パネル615を操作した場合の動作について説明する。
図54は、再生制御パネルの操作によりヒートマップおよび立体図が再生表示される状態を示す図である。図55は、再生制御パネルの操作によりヒートマップおよび立体図がコマ戻しされる状態を示す図である。図56は、再生制御パネルの操作によりヒートマップおよび立体図がコマ送りされる状態を示す図である。図54〜図56を参照しながら、時間周波数解析画面601の再生制御パネル615を操作した場合の動作について説明する。
再生制御パネル615は、解析者の操作によって、時間経過と共に、ヒートマップ611、立体図612および頭部三面図613の状態を再生表示するためのユーザインターフェースである。
例えば、解析者が、再生制御パネル615の「再生」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211に対して、図54(a)および図54(b)に示すように、ヒートマップ611で指定されている指定領域622−1、およびその周りの対応領域622−2〜622−5を、時間経過と共に、右方向(時間が進む方向)に移動させるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1および対応領域622−2〜622−5の移動に伴い、立体表示制御部212に対して、図54(a)および図54(b)に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、断面表示制御部213に対して、移動している指定領域622−1に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像644に表示させるように指示する。
また、解析者が、再生制御パネル615の「コマ戻し」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211に対して、図55(a)および図55(b)に示すように、ヒートマップ611で指定されている指定領域622−1、およびその周りの対応領域622−2〜622−5を、所定時間分だけ左方向(時間が戻る方向)に移動させるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1および対応領域622−2〜622−5の移動に伴い、立体表示制御部212に対して、図55(a)および図55(b)に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、断面表示制御部213に対して、移動した指定領域622−1に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像644に表示させるように指示する。
また、解析者が、再生制御パネル615の「コマ送り」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211に対して、図56(a)および図56(b)に示すように、ヒートマップ611で指定されている指定領域622−1、およびその周りの対応領域622−2〜622−5を、所定時間分だけ右方向(時間が進む方向)に移動させるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1および対応領域622−2〜622−5の移動に伴い、立体表示制御部212に対して、図56(a)および図56(b)に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、断面表示制御部213に対して、移動した指定領域622−1に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像644に表示させるように指示する。
また、解析者が、再生制御パネル615の「停止」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211、立体表示制御部212および断面表示制御部213に対して、ヒートマップ611、立体図612、および頭部三面図613におけるそれぞれの再生表示動作を停止するように指示する。
また、解析者が、再生制御パネル615の「先頭へ移動」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211に対して、ヒートマップ611で指定されている指定領域622−1を、時刻の先頭に移動させるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、立体表示制御部212に対して、指定領域622−1に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、断面表示制御部213に対して、移動した指定領域622−1に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像644に表示させるように指示する。
また、解析者が、再生制御パネル615の「末尾へ移動」ボタンを押下すると、再生表示制御部214は、ヒートマップ表示制御部211に対して、ヒートマップ611で指定されている指定領域622−1を、時刻の末尾に移動させるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、立体表示制御部212に対して、指定領域622−1に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部214は、ヒートマップ611における指定領域622−1の移動に伴い、断面表示制御部213に対して、移動した指定領域622−1に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像644に表示させるように指示する。
以上のように、再生表示させることによって、立体図612および頭部三面図613に表示される信号強度の分布(ヒートマップ)の時間経過による変化を動画として確認することができ、時間経過によるピークの移動等を視覚的に確認することができる。
<初期表示について>
図57は、ピークに対してどの視点からの図を初期表示させるかを説明する図である。図58は、2つのピークに対してどの視点からの図を初期表示させるかを説明する図である。図59は、図58に説明する視点からの図を立体図に初期表示させた状態を示す図である。図57〜図59を参照しながら、時間周波数解析画面601を起動した(開いた)場合におけるヒートマップ611、立体図612および頭部三面図613の初期表示について説明する。
図57は、ピークに対してどの視点からの図を初期表示させるかを説明する図である。図58は、2つのピークに対してどの視点からの図を初期表示させるかを説明する図である。図59は、図58に説明する視点からの図を立体図に初期表示させた状態を示す図である。図57〜図59を参照しながら、時間周波数解析画面601を起動した(開いた)場合におけるヒートマップ611、立体図612および頭部三面図613の初期表示について説明する。
解析者が、時間周波数解析画面601を起動した(開いた)場合に、ヒートマップ611、立体図612および頭部三面図613の初期表示としてどのような画像を表示させるかについての類型を説明する。
例えば、解析表示制御部202は、時刻・周波数全体を通じて、脳全体において信号強度が最大となる時刻・周波数および脳内の位置を求める。この場合、ヒートマップ表示制御部211は、解析表示制御部202により求められた脳内の位置でのヒートマップ611を表示させる。また、立体表示制御部212は、解析表示制御部202により求められた信号強度が最大となる時刻・周波数に対応する脳の立体画像を立体図612に表示させる。また、断面表示制御部213は、解析表示制御部202により求められた脳内の位置を通るような三面図を頭部三面図613に表示させ、解析表示制御部202により求められた信号強度が最大となる時刻・周波数のヒートマップを三面図および立体画像644に重畳させる。
また、解析表示制御部202は、時刻・周波数全体を通じての信号強度の平均が最大となる脳内の位置を求めるものとしてもよい。この場合、ヒートマップ表示制御部211は、解析表示制御部202により求められた脳内の位置でのヒートマップ611を表示させる。また、立体表示制御部212は、表示されたヒートマップ611のうち信号強度が最大となる時刻・周波数に対応する脳の立体画像を立体図612に表示させる。また、断面表示制御部213は、解析表示制御部202により求められた脳内の位置を通るような三面図を頭部三面図613に表示させ、表示されたヒートマップ611のうち信号強度が最大となる時刻・周波数のヒートマップを三面図および立体画像644に重畳させる。
また、解析表示制御部202は、脳全体で信号強度の平均値が最大となる時刻・周波数を求めるものとしてもよい。この場合、立体表示制御部212は、解析表示制御部202により求められた時刻・周波数に対応する脳の立体画像を立体図612に表示させる。また、ヒートマップ表示制御部211は、立体図612の立体画像に表示された、解析表示制御部202により求められた時刻・周波数に対応するヒートマップにおいて信号強度が最大となる脳内の位置を求め、当該位置でのヒートマップ611を表示させる。また、断面表示制御部213は、ヒートマップ表示制御部211により求められた脳内の位置を通るような三面図を頭部三面図613に表示させ、解析表示制御部202により求められた時刻・周波数のヒートマップを三面図および立体画像644に重畳させる。
また、立体表示制御部212は、ピークリスト614に登録されているピーク情報のうち、先頭のピーク情報が示す脳内の位置でのヒートマップ611を表示させるものとしてもよい。また、立体表示制御部212は、ピークリスト614に登録されているピーク情報のうち、先頭のピーク情報が示す時刻・周波数に対応する脳の立体画像を立体図612に表示させる。また、断面表示制御部213は、ピークリスト614に登録されているピーク情報のうち、先頭のピーク情報が示す脳内の位置を通るような三面図を頭部三面図613に表示させ、当該ピーク情報が示す時刻・周波数のヒートマップを三面図および立体画像644に重畳させる。
また、立体表示制御部212は、測定の対象(視覚野、聴覚野、体性感覚野、運動野および言語野等)に応じてプリセットされた脳内の位置でのヒートマップ611を表示させるものとしてもよい。また、立体表示制御部212は、測定の対象(視覚野、聴覚野、体性感覚野、運動野および言語野等)に応じてプリセットされた時刻・周波数に対応する脳の立体画像を立体図612に表示させる。また、断面表示制御部213は、測定の対象(視覚野、聴覚野、体性感覚野、運動野および言語野等)に応じてプリセットされた脳内の位置を通るような三面図を頭部三面図613に表示させ、当該ピーク情報が示す時刻・周波数のヒートマップを三面図および立体画像644に重畳させる。
次に、解析者が、時間周波数解析画面601を起動した(開いた)場合に、立体図612の脳の立体画像、および、頭部三面図613の立体画像644の表示の初期視点について説明する。
例えば、測定の対象(視覚野、聴覚野、体性感覚野、運動野および言語野等)に応じてプリセットされた視点を初期視点として用いるものとしてもよい。この場合、立体図612については、行(視点)の数もプリセットしておく。立体図612が2行である場合、2つの視点をプリセットしておく必要がある。例えば、言語野が測定対象であれば、脳の左側面および右側面の視点をプリセットしておく。
また、ピークリスト614において先頭に登録されているピークが最もよく見えるような視点を初期視点として用いるものとしてもよい。具体的には、図57に示すように、脳の中心とピークとを結んだ直線811上に初期視点として視点P0を設定するものとすればよい。
また、ピークリスト614において所定のパラメータ(例えば、図50に示すピークの値(信号強度)またはピークの高さ)が、所定の閾値を超えたピークを用いて設定した視点を初期視点として用いるものとしてもよい。例えば、閾値を超えたピークが2つある場合、立体図612の表示を2行表示とし、図58に示すように、脳の中心とそれぞれのピークとを結んだ直線812、813上に、初期視点としてそれぞれ視点P1、P2を設定するものとすればよい。この場合において、視点P1からの脳の立体画像を立体図612の上の行に、視点P2からの脳の立体画像を立体図612の下の行に表示した例を、図59に示す。
以上のように、脳内の特定の位置または特定の範囲における生体信号の時刻・周波数に関するヒートマップ611を表示するものとし、そのヒートマップ611上で指定された点または範囲に対応する脳の活動を示すヒートマップが重畳された立体画像を中心に、その時刻の前後の脳の活動を示す立体画像を表示、すなわち、脳の活動を示す静止画(ここでは立体画像)をコマ送り的またはコマ戻し的に表示するものとしている。これによって、脳の活動を示す静止画を適切かつ迅速に抽出することができ、脳の活動の解析を容易にすることができる。また、カンファレンス等で議論のベースにすることが容易となる。
また、ピークリスト614に登録されたピーク情報から特定のピーク情報を選択することによって、当該ピーク情報に対応したヒートマップ611、立体図612および頭部三面図613が表示される。これによって、選択したピークが脳のどの位置、およびどの時刻・周波数であるかを瞬時に認識することができ、さらに、ヒートマップ611においては当該ピークおよびその周りの時刻/周波数における信号強度の状態、および、当該ピークの位置およびその周りの脳の信号強度の状態も把握することが可能となる。
また、立体図612上で自由に脳の視点を変更することができ、さらに、その視点変更に基づく変更を同じ行または異なる行の脳に対しても反映することができる。これによって、特定の立体画像(対象立体画像)に対する視点変更を行うのみで、他の立体画像についてもその視点変更に基づく変更が反映されるので、操作性が向上し、さらに、複数の行を見比べることによって、対応する視点から時間的に前後の脳の活動変化を容易に確認することができる。また、立体画像で描画される脳の視点を自在に変えることによって、ある視点では見えなかった発火位置を確認することができる。
また、頭部三面図613における立体画像644に対して行った視点変更を、各種設定に応じて、立体図612の時系列に並んだ脳の立体画像の視点に反映することができる。これによって、立体画像644の視点変更と同様の視点変更を、立体図612に対して改めて行う必要がないので、操作性が向上し、さらに、立体画像644で変更した視点と同一の視点、または対応する視点で、立体図612において脳の状態の変化を時系列に確認することができる。
なお、上述の実施形態は、生体部位としての脳についての生体信号を取り扱うものとしたが、これに限定されるものではなく、例えば、脊髄、筋肉等の生体部位の生体信号に対して適用することが可能である。例えば、脳の図で説明した立体図612を、腰椎の場合、図60のように表示することが可能である。図60では、図60(a)〜図60(d)の順で時系列に腰椎の信号が上に伝わっていく状態を示している。
(解析結果を示すマークの重畳処理)
図61は、脳の断面図の一例を示す図である。図62は、マーク画像のサンプルを示す図である。図63は、断面図にマーク画像のサンプルを重畳させた図である。図64は、断面図の画素値についてのヒストグラムの一例を示す図である。図61〜図64を参照しながら、生体画像に解析結果を示すマークを重畳させる処理について説明する。
図61は、脳の断面図の一例を示す図である。図62は、マーク画像のサンプルを示す図である。図63は、断面図にマーク画像のサンプルを重畳させた図である。図64は、断面図の画素値についてのヒストグラムの一例を示す図である。図61〜図64を参照しながら、生体画像に解析結果を示すマークを重畳させる処理について説明する。
図61に示す断面図661は、例えば、上述の図32等で示した頭部三面図613のうちの断面図641等に対応する生体画像の一例としての脳の断面図であり、例えば、DICOM形式のフォーマットに基づく画像である。図61に示す断面図661は、例えば、MRI検査により得られるMRI画像である。
図61に示すように、MRI画像等である断面図661は、複雑な形状が入り組んでおり、また通常は視認性を高めるためにすべての濃度を用いて表示される。このような断面図661を構成する各画素は、通常1バイトの情報により256階調の濃度である画素値が表され、V=0〜255の値をとる。ここでは、V=0が黒、V=255が白を表すものとして説明する。
例えば、ピークリスト制御部203により抽出されピークリスト614に登録された信号強度のピークの位置、または、ダイポール推定の位置等を示すマーク等の情報(以下、単にピークの情報と称する)(結果情報の一例)は、医者による診察処理において重要である。したがって、断面図661に示される脳画像中に、ピークの情報を何らかのマーク等として重畳させて示すことが好ましい。ここで、図62に、断面図661に重畳させるためのマークの一例であるマーク画像761を示す。マーク画像761には、画像全体において4か所にマーク群852a〜852dが含まれている。図62ではすべてのマークが視認可能なように、背景の画素値をV=128としてグレー表示にしているが、実際には、マーク以外の部分は透明として扱い、マーク部分のみが断面図661等の生体画像に重畳される。また、マーク画像761に含まれるマーク群852a〜852dそれぞれには、4つのマークが含まれており、4つマークのうち、左上のマークがV=0(黒)、右上のマークがV=255(白)、右下のマークがV=64(濃いグレー)、左下のマークがV=192(薄いグレー)としている。なお、図62に示すマークの形状は十字形状としているが、これは例示であり、矢印形状または×形状等の他の形状のマークであってもよい。
図62に示すマーク群852a〜852dに含まれる各マークを、断面図661に重畳させた画像として、図63に断面図662として示す。図62に示すように、マークがいかなる濃度を用いても、ピークの位置によっては、脳画像を構成する複雑な模様および複雑な画素値構成により視認性が悪い、すなわち、マークが見づらいケースがあることが把握される。
ここで、実際に生体画像に重畳させるマークに用いる画素値を求める方法について説明する。図64に示すグラフは、図61に示した断面図661において、横軸を断面図661の各画素の画素値(濃度)、縦軸をその画素値を有する画素の数(画素数)とするヒストグラムである。情報処理装置50の生成部221は、このようなヒストグラムを生成する。なお、生成部221によるヒストグラムの生成とは、必ずしも図64のようにグラフとして表示させるためのヒストグラムを生成することに限られるものではなく、対象となる生体画像(ここでは断面図661)を構成する各画素の画素値と、当該画素値を有する画素の数とを関連付ける情報(第1情報の一例)を生成することを含む。
そして、情報処理装置50の特定部222は、生成部221により生成されたヒストグラムにおいて、最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いるものとする。図64に示すヒストグラムの例では、V=240の画素値の画素が最も少ないことが特定される。そして、情報処理装置50の重畳部223は、特定部222により特定された画素値(濃度)をマークの画素値にして、断面図661に重畳させる。この場合、同じ濃度を有する画素は少ないので、特定部222により特定された画素値(濃度)のマークを断面図661に重畳させた場合、マークの視認性を向上させることができる。特に、特定部222により特定された画素値(濃度)の画素数が0である場合には、当該画素値は断面図661内で用いられていない画素値なので、マークの視認性がさらに高まることになる。
そして、情報処理装置50の出力部224は、重畳部223により断面図661にマークが重畳されて得られた画像を、例えば、DICOM形式のフォーマットの画像ファイルとして出力する。また、出力部224により出力された画像ファイルについて、特定の濃度の画素を着色する機能を有するビューワを用いて表示させる場合、重畳部223により重畳されたマークと同じ濃度を有する画素は少ないので、当該マークの濃度に対して着色することによって、さらにマークの視認性を向上させることができる。
以上のように、本実施形態に係る情報処理装置50では、対象となる生体画像を構成する各画素の画素値と、当該画素値を有する画素の数とを関連付ける情報(例えば、ヒストグラム)を生成し、当該ヒストグラムにおいて、最も画素数が少ない画素値(濃度)を特定し、そして、特定した画素値(濃度)をマークの画素値にして、対象となる生体画像に重畳させるものとしている。これによって、対象となる生体画像において、特定された画素値と同じ画素値を有する画素は少ないので、重畳されたマークの視認性を向上させることができる。
なお、ピークの情報を示すマークを重畳させる生体画像として、図61に示すような脳の断面図である断面図661を例としたが、これに限定されるものではなく、上述の図32に示す脳の立体画像644等のその他の生体画像であってもよい。
(変形例1)
上述の実施形態では、生成部221が対象となる生体画像全体についてのヒストグラムを生成し、特定部222が当該ヒストグラムにおいて最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いることを説明した。本変形例では、対象となる生体画像のうち、マークの位置を含む所定の領域においてヒストグラムを生成する動作を説明する。
上述の実施形態では、生成部221が対象となる生体画像全体についてのヒストグラムを生成し、特定部222が当該ヒストグラムにおいて最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いることを説明した。本変形例では、対象となる生体画像のうち、マークの位置を含む所定の領域においてヒストグラムを生成する動作を説明する。
対象となる生体画像全体でヒストグラムを生成する場合、上述の実施形態で説明した方法でマークの濃度を特定した場合でも、たまたまマークの周囲の同じ濃度の画素が多数存在している場合もあり得る。この場合、特定した濃度でマークを重畳させても、同じ濃度を有する画素に埋もれて、当該マークの視認性が意図した通りには向上できない可能性もある。そこで、本変形例では、対象となる生体画像において、マークの位置を中心とする所定の領域を定め、生成部221は、当該所定の領域についてヒストグラムを生成し、特定部222は、当該所定の領域について生成されたヒストグラムにおいて、最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いるものとする。
これによって、マークに含まれる所定の領域には当該マークと同じ濃度を有する画素が少ないので、マークの視認性をさらに向上させることができる。また、ヒストグラムにおいて谷点が複数ある、すなわち、画素数が極小となる濃度が複数ある場合、同様の方法により、マーク付近に少ない濃度を選択することもできる。
なお、上述の所定の領域は、マークの位置を中心とする領域に限定されるものではなく、少なくとも当該マークを含む領域であってもよい。
(変形例2)
上述の実施形態では、生成部221が対象となる生体画像全体についてのヒストグラムを生成し、特定部222が当該ヒストグラムにおいて最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いることを説明した。本変形例では、特定した画素値(濃度)を有するマークと同じ画素値(濃度)を有する画素の当該画素値を変更する動作について説明する。
上述の実施形態では、生成部221が対象となる生体画像全体についてのヒストグラムを生成し、特定部222が当該ヒストグラムにおいて最も画素数が少ない画素値(濃度)を特定し、その画素値をマークの画素値(濃度)として用いることを説明した。本変形例では、特定した画素値(濃度)を有するマークと同じ画素値(濃度)を有する画素の当該画素値を変更する動作について説明する。
マークの視認性をさらに向上させる方法として、特定部222によりマークの画素値として特定されたものと、同じ画素値を有する画素を、対象とする生体画像から削除する方法がある。本変形例では、当該方法について説明する。ただし、生体画像の性質上、当該生体画像の見た目の影響を極力排除することが望まれる。
図65は、変形例2に係る情報処理装置の機能ブロック構成の一例を示す図である。図65を参照しながら、本変形例に係る情報処理装置50aの機能ブロック構成について説明する。
図65に示すように、情報処理装置50aは、収集表示制御部201と、解析表示制御部202と、ピークリスト制御部203と、通信部204と、センサ情報取得部205と、解析部206と、記憶部207と、入力部208と、生成部221と、特定部222と、変更部225と、重畳部223と、出力部224と、を有する。なお、このうち収集表示制御部201、解析表示制御部202、ピークリスト制御部203、通信部204、センサ情報取得部205、解析部206、記憶部207、入力部208、生成部221、特定部222および出力部224の動作は、上述の実施形態で説明した通りである。
変更部225は、対象となる生体画像において、特定部222により特定された画素値と同じ画素値を有する画素の当該画素値を変更する機能部である。例えば、変更部225は、特定部222により特定された画素値と同じ画素値を有する画素の当該画素値を±1だけ変更する。
重畳部223は、特定部222により特定された画素値をマークの画素値にして、変更部225により画素値が変更された後の生体画像に重畳させる。
上述の収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、変更部225、重畳部223および出力部224は、CPU101がROM103等に記憶されたプログラムをRAM102に展開して実行することにより実現される。なお、収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、変更部225、重畳部223および出力部224の一部または全部は、ソフトウェアであるプログラムではなく、ASICまたはFPGA等のハードウェア回路によって実現されてもよい。
なお、図65に示した各機能部は、機能を概念的に示したものであって、このような構成に限定されるものではない。例えば、図65で独立した機能部として図示した複数の機能部を、1つの機能部として構成してもよい。一方、図65の1つの機能部が有する機能を複数に分割し、複数の機能部として構成するものとしてもよい。
図66は、断面図にV=64のマークを重畳した場合の一例を示す図である。図67は、断面図のうちV=64の画素値を有する画素を着色した図である。図68は、断面図のうちV=64のマークのみ着色した図である。図66〜図68を参照しながら、特定部222により特定された画素値と同じ画素値を有する画素の当該画素値を変更する動作の詳細を説明する。
図66に示す断面図663は、V=64のマーク853を重畳させた脳の断面図を示す。変形例1で上述したように、ヒストグラムを用いて特定した濃度でマークを重畳させても、同じ濃度を有する画素に埋もれて、当該マークの視認性が意図した通りには向上できない可能性もある。すなわち、図66に示す断面図663にように、特定された画素値(濃度)がV=64であり、当該画素値でマークを重畳させたとしても、視認性が意図した通りには向上できない場合もあり得る。また、上述の実施形態で説明したように、重畳部223により重畳されたマークと同じ濃度を有する画素に対して着色した場合、図67に示す断面図663aのようになり、着色されたマーク853a以外にもV=64の画素値(濃度)を有する画素が散見され、この場合もマークの視認性が十分であるとは言い難い。そこで、本変形例では、以下のように画素値を変更する処理を施す。
本変形例では、特定部222によりヒストグラムを用いてマークの画素値が特定された後、変更部225が、対象となる生体画像において、例えば、特定部222により特定された画素値と同じ画素値を有する画素の当該画素値を±1だけ変更する。具体的には、画素値がV=64のマークを重畳させる場合には、変更部225は、対象となる生体画像において、V=64の画素の画素値をV=63、またはV=65に変更する。そして、重畳部223は、特定部222により特定された画素値をマークの画素値にして、変更部225により画素値が変更された後の生体画像に重畳させる。
これによって、対象となる生体画像の見た目の影響を極力排除しつつ、マークが重畳された生体画像において、マークの画素値(ここではV=64)と同じ画素値を有する画素は、マークを構成する画素のみとすることができる。この場合、マークの濃度に対して着色した画像の例を、図68に示す。図68に示す断面図664のように、着色された画素はマークを構成する画素のみであるので、着色されたマークであるマーク854について非常に視認性を高く認識することができる。
また、変更部225で、特定部222により特定された画素値と同じ画素値を有する画素の当該画素値を−1すなわちV=63とするか、+1すなわちV=65とするかについては、画素値を変更する画素の周辺の画素値に基づいて、生体画像の見た目への影響を小さくする方法で選択すればよい。例えば、変更部225は、画素値の変更対象となる画素の周囲26画素、すなわち、画素値を変更する画素を取り囲む3×3×3の画素から当該画素を除いた26画素の画素値を参照し、V=63に近い画素の数と、V=65に近い画素の数とを数え、画素の数が多い方の画素値に変更するという方法が挙げられる。これによって、変更部225により画素値が変更された生体画像について見た目の影響を極力排除することができる。
なお、変更部225により変更される画素値の幅は±1に限定されるものではなく、生体画像の見た目の影響を極力排除できる範囲であれば、画素値を変更する幅を広げてもよい。
(変形例3)
上述の実施形態および変形例1、2では、対象となる生体画像にピークの情報としてのマークを1つのみ重畳させる例を示した。本変形例では、生体画像に複数のマークを重畳させる動作について説明する。
上述の実施形態および変形例1、2では、対象となる生体画像にピークの情報としてのマークを1つのみ重畳させる例を示した。本変形例では、生体画像に複数のマークを重畳させる動作について説明する。
図69は、断面図の複数のマークを重畳した場合の一例を示す図である。図70は、断面図の複数のマークを着色した図である。図69および図70を参照しながら、生体画像に複数のマークを重畳させる動作について説明する。
生体画像に重畳を所望するマークは1つとは限られず、属性の異なる複数のマーク、例えば、信号強度が1番目に強いピークおよび2番目に強いピークをそれぞれ示すマーク、または、信号強度の上向きのピークおよび下向きのピークをそれぞれ示すマーク等のように複数のピークを重畳させることが望まれる場合もある。この場合、例えば、特定部222は、生成部221により生成されたヒストグラムを用いて、画素数が少ない画素値(濃度)を所望されるマークの数だけ特定する。例えば、特定部222は、ヒストグラムにおいて、複数の極小値から複数の画素値を特定すればよく、最も画素数が少ない画素値から±1だけ変更した画素値をマークの画素値として特定するものとしてもよい。そして、重畳部223は、特定部222により特定された複数の画素値(濃度)をマークの画素値にして、複数のマークを対象となる生体画像に重畳させる。この場合、変更部225により、生体画像を構成する画素のうち、特定部222により特定された複数の画素値と同じ画素値を有する画素の当該画素値が変更されるようにしてもよい。
以上のような動作で、重畳部223により画素値の異なる複数のマークが重畳された生体画像の一例を、図69に示す。図69に示す断面図665では、信号強度が1番強いピークを示すマーク856にV=66の画素値を割り当て、2番目に強いピークを示すマーク855にV=64を割り当てた例を示す。そして、特定の濃度の画素を着色する機能を有するビューワを用いて、V=66のマーク856を赤で着色してマーク856aとし、V=64のマーク855を青に着色してマーク855aとした例を、図70の断面図665aに示す。
以上のように、複数のマークを視認性高く生体画像に重畳させることができ、当該複数のマークを異なる色で着色することによって、さらに視認性を向上させることができる。
(変形例4)
上述の実施形態および変形例1〜3では、対象となる生体画像に対して、ピークリスト制御部203により抽出されピークリスト614に登録された信号強度のピークの位置、または、ダイポール推定の位置等を示すマークのみを重畳させる動作について説明した。本変形例では、生体画像に重畳したマークを着色するために必要な、当該マークの画素値(濃度)を含むテキストを生体画像に付加する動作について説明する。
上述の実施形態および変形例1〜3では、対象となる生体画像に対して、ピークリスト制御部203により抽出されピークリスト614に登録された信号強度のピークの位置、または、ダイポール推定の位置等を示すマークのみを重畳させる動作について説明した。本変形例では、生体画像に重畳したマークを着色するために必要な、当該マークの画素値(濃度)を含むテキストを生体画像に付加する動作について説明する。
対象となる生体画像に対してマークのみを重畳させたとしても、特定の濃度の画素を着色する機能を有するビューワを用いて当該マークを着色する場合に、当該マークの画素値を知得していなければ着色操作を行うことができない。そこで、本変形例では、生体画像に重畳したマークを着色するために必要な、当該マークの画素値(濃度)を含むテキストを生体画像に付加することによって、ユーザに対してマークの画素値を知得させる動作について説明する。
図71は、変形例4に係る情報処理装置の機能ブロック構成の一例を示す図である。図72は、断面図にテキストを付加した場合の一例を示す図である。図71および図72を参照しながら、本変形例に係る情報処理装置50bの機能ブロックの構成および動作について説明する。
図71に示すように、情報処理装置50bは、収集表示制御部201と、解析表示制御部202と、ピークリスト制御部203と、通信部204と、センサ情報取得部205と、解析部206と、記憶部207と、入力部208と、生成部221と、特定部222と、変更部225と、重畳部223と、付加部226と、出力部224と、を有する。なお、このうち収集表示制御部201、解析表示制御部202、ピークリスト制御部203、通信部204、センサ情報取得部205、解析部206、記憶部207、入力部208、生成部221、特定部222、変更部225、重畳部223および出力部224の動作は、上述の実施形態および変形例2で説明した通りである。また、情報処理装置50bは変更部225を有するものとしているが、これに限定されえるものではなく、変更部225を有さないものとしてもよい。
付加部226は、重畳部223によりマークが重畳された生体画像に対して、当該マークの画素値(濃度)を含むテキスト(第2情報の一例)を付加する機能部である。
上述の収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、変更部225、重畳部223、付加部226および出力部224は、CPU101がROM103等に記憶されたプログラムをRAM102に展開して実行することにより実現される。なお、収集表示制御部201、解析表示制御部202、ピークリスト制御部203、センサ情報取得部205、解析部206、生成部221、特定部222、変更部225、重畳部223、付加部226および出力部224の一部または全部は、ソフトウェアであるプログラムではなく、ASICまたはFPGA等のハードウェア回路によって実現されてもよい。
なお、図71に示した各機能部は、機能を概念的に示したものであって、このような構成に限定されるものではない。例えば、図71で独立した機能部として図示した複数の機能部を、1つの機能部として構成してもよい。一方、図71の1つの機能部が有する機能を複数に分割し、複数の機能部として構成するものとしてもよい。
上述の図69に示すマーク855、856が重畳された断面図665に対して、付加部226により、マーク855、856が示す情報の内容、および画素値を示すテキスト表示部866が付加された断面図666を、図72に示す。図72に示すように、テキスト表示部866は、断面図666のうち脳画像部分以外の部分に付加されることが望ましい。これによって、生体画像のうち生体を示す部分にテキストが付加されることにより、当該生体の描画を損なうこと防止することができる。また、図72に示すように、付加部226により付加されるテキストとしては、マークの画素値だけでなく、当該マークが示す情報の内容を含むものとすることによって、ユーザは、生体画像に重畳されているマークが何の情報であるのかを把握することができる。
以上のように、本変形例では、付加部226によって、重畳部223によりマークが重畳された生体画像に対して、当該マークの画素値(濃度)を含むテキストが付加されるものとしている。これによって、ユーザは、生体画像に重畳されているマークの画素値(濃度)を知得することができ、当該マークについて着色操作を行うことができる。
なお、上述の実施形態または各変形例では、出力部224により、重畳部223によりマークが重畳された生体画像をDICOM形式のフォーマットに変換したファイルとして出力するものとしているが、これに加えて、当該マークが重畳された生体画像を、例えば、図11に示す時間周波数解析画面601の頭部三面図613に表示させるものとしてもよい。さらに、この場合、マークが重畳された生体画像に、信号強度の分布であるヒートマップを重畳させて表示するものとしてもよい。
また、上述の実施形態および各変形例において、生体信号計測システム1の各機能部の少なくともいずれかがプログラムの実行によって実現される場合、そのプログラムは、ROM等に予め組み込まれて提供される。また、上述の実施形態および各変形例に係る生体信号計測システム1で実行されるプログラムは、インストール可能な形式または実行可能な形式のファイルでCD−ROM、フレキシブルディスク(FD)、CD−R(Compact Disk Recordable)、DVD(Digital Versatile Disc)等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよい。また、上述の実施形態および各変形例の生体信号計測システム1で実行されるプログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、上述の実施形態および各変形例の生体信号計測システム1で実行されるプログラムを、インターネット等のネットワーク経由で提供または配布するように構成してもよい。また、上述の実施形態および各変形例の生体信号計測システム1で実行されるプログラムは、上述した各機能部のうち少なくともいずれかを含むモジュール構成となっており、実際のハードウェアとしてはCPUがROM等からプログラムを読み出して実行することにより、上述の各機能部が主記憶装置上にロードされて生成されるようになっている。
1 生体信号計測システム
3 測定装置
4 測定テーブル
31 デュワ
32 窪み
40 サーバ
50、50a、50b 情報処理装置
101 CPU
102 RAM
103 ROM
104 補助記憶装置
105 ネットワークI/F
106 入力装置
107 表示装置
108 バス
201 収集表示制御部
202 解析表示制御部
203 ピークリスト制御部
204 通信部
205 センサ情報取得部
206 解析部
207 記憶部
208 入力部
211 ヒートマップ表示制御部
212 立体表示制御部
213 断面表示制御部
214 再生表示制御部
221 生成部
222 特定部
223 重畳部
224 出力部
225 変更部
226 付加部
501 開始画面
502 測定収集画面
511a、511b 領域
512 モニタウィンドウ
521〜523 第2表示領域
523a−1、523a−2 マーク
526−1、526−2 属性アイコン
530 第1表示領域
530a−1、530a−2 アノテーション
531 時間軸
532 ライン
535 ポップアップウィンドウ
535a 選択ボタン
535b 入力ボックス
538 カウンタボックス
539 終了ボタン
541、542 脳磁分布図
550 脳波分布図
560 アノテーションリスト
560a 選択ボックス
601 時間周波数解析画面
605 解析画面切替リスト
611 ヒートマップ
612、612a 立体図
612−1〜612−3 表示領域
613、613a 頭部三面図
614 ピークリスト
614a ピークリスト設定ボタン
615 再生制御パネル
621、621−1 指定部
621−2〜621−5 対応部
622、622−1 指定領域
622−2〜622−5 対応領域
623 指定領域
631−1〜631−5 線分
635 コメント
641〜643 断面図
644 立体画像
645a〜645d 基準線
646、646a、646b ピーク点
647 ダイポール
650〜652 ダイアログボックス
655 ダイアログボックス
656 プルダウンリスト
661〜663、663a、664、665、665a、666 断面図
712a−1〜712a−5 部位
712b−1〜712b−5 部位
713a−1、713a−2、713b〜713d 部位
722−1〜722−5 枠
761 マーク画像
801、802 ピーク部
811〜813 直線
852a〜852d マーク群
853、853a、854、855、855a、856、856a マーク
866 テキスト表示部
A1、A2 アノテーション
P0〜P2 視点
3 測定装置
4 測定テーブル
31 デュワ
32 窪み
40 サーバ
50、50a、50b 情報処理装置
101 CPU
102 RAM
103 ROM
104 補助記憶装置
105 ネットワークI/F
106 入力装置
107 表示装置
108 バス
201 収集表示制御部
202 解析表示制御部
203 ピークリスト制御部
204 通信部
205 センサ情報取得部
206 解析部
207 記憶部
208 入力部
211 ヒートマップ表示制御部
212 立体表示制御部
213 断面表示制御部
214 再生表示制御部
221 生成部
222 特定部
223 重畳部
224 出力部
225 変更部
226 付加部
501 開始画面
502 測定収集画面
511a、511b 領域
512 モニタウィンドウ
521〜523 第2表示領域
523a−1、523a−2 マーク
526−1、526−2 属性アイコン
530 第1表示領域
530a−1、530a−2 アノテーション
531 時間軸
532 ライン
535 ポップアップウィンドウ
535a 選択ボタン
535b 入力ボックス
538 カウンタボックス
539 終了ボタン
541、542 脳磁分布図
550 脳波分布図
560 アノテーションリスト
560a 選択ボックス
601 時間周波数解析画面
605 解析画面切替リスト
611 ヒートマップ
612、612a 立体図
612−1〜612−3 表示領域
613、613a 頭部三面図
614 ピークリスト
614a ピークリスト設定ボタン
615 再生制御パネル
621、621−1 指定部
621−2〜621−5 対応部
622、622−1 指定領域
622−2〜622−5 対応領域
623 指定領域
631−1〜631−5 線分
635 コメント
641〜643 断面図
644 立体画像
645a〜645d 基準線
646、646a、646b ピーク点
647 ダイポール
650〜652 ダイアログボックス
655 ダイアログボックス
656 プルダウンリスト
661〜663、663a、664、665、665a、666 断面図
712a−1〜712a−5 部位
712b−1〜712b−5 部位
713a−1、713a−2、713b〜713d 部位
722−1〜722−5 枠
761 マーク画像
801、802 ピーク部
811〜813 直線
852a〜852d マーク群
853、853a、854、855、855a、856、856a マーク
866 テキスト表示部
A1、A2 アノテーション
P0〜P2 視点
Claims (13)
- 生体の解析をするための情報処理装置であって、
前記生体を示す画像である生体画像の画素値に基づいて、前記解析の結果を示す結果画像に用いる画素値を特定する特定部と、
前記特定部により特定された画素値とした前記結果画像を、前記生体画像に重畳させる重畳部と、
を備えた情報処理装置。 - 前記生体画像に含まれる画素の画素値と、該画素値を有する画素の画素数とを関連付ける第1情報を生成する生成部を、さらに備え、
前記特定部は、前記第1情報に基づいて、前記結果画像に用いる画素値を特定する請求項1に記載の情報処理装置。 - 前記特定部は、前記第1情報から、画素数が最も少ない画素値を特定する請求項2に記載の情報処理装置。
- 前記生成部は、前記生体画像における前記結果画像の位置を含む所定の領域を構成する画素の画素値と、該画素値を有する画素の画素数とを関連付ける情報を前記第1情報として生成する請求項2または3に記載の情報処理装置。
- 前記生体画像において、前記特定部により特定された画素値と同じ画素値を有する画素の該画素値を変更する変更部を、さらに備え、
前記重畳部は、前記特定部により特定された画素値とした前記結果画像を、前記変更部により画素値が変更された前記生体画像に重畳させる請求項1〜4のいずれか一項に記載の情報処理装置。 - 前記変更部は、前記生体画像において、前記特定部により特定された画素値と同じ画素値を有する画素の周辺画素の画素値に基づいて、前記同じ画素値を有する画素の該画素値を変更する請求項5に記載の情報処理装置。
- 前記変更部は、前記同じ画素値を有する画素の画素値を、前記周辺画素の画素値のうち画素数が最も多い画素値へ近づくように変更する請求項6に記載の情報処理装置。
- 前記特定部は、前記結果画像が複数である場合、前記各結果画像に対してそれぞれ異なる画素値を特定し、
前記重畳部は、前記各結果画像を、前記特定部により特定されたそれぞれ異なる画素値として前記生体画像に重畳させる請求項1〜7のいずれか一項に記載の情報処理装置。 - 前記重畳部により前記結果画像が重畳された前記生体画像に対して、少なくとも該結果画像の画素値を示す第2情報を付加する付加部を、さらに備えた請求項1〜8のいずれか一項に記載の情報処理装置。
- 前記付加部は、前記重畳部により前記結果画像が重畳された前記生体画像に対して、前記結果画像に対応する前記解析の結果内容を含めた前記第2情報を付加する請求項9に記載の情報処理装置。
- 生体を示す画像である生体画像の画素値に基づいて、前記生体の解析の結果を示す結果画像に用いる画素値を特定する特定ステップと、
特定した画素値とした前記結果画像を、前記生体画像に重畳させる重畳ステップと、
を有する情報処理方法。 - コンピュータに、
生体を示す画像である生体画像の画素値に基づいて、前記生体の解析の結果を示す結果画像に用いる画素値を特定する特定ステップと、
特定した画素値とした前記結果画像を、前記生体画像に重畳させる重畳ステップと、
を実行させるためのプログラム。 - 被験者から1種類以上の生体信号を測定する測定装置と、
請求項1〜10のいずれか一項に記載の情報処理装置と、
を備えた情報処理システム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019045867A JP2020146204A (ja) | 2019-03-13 | 2019-03-13 | 情報処理装置、情報処理方法、プログラムおよび情報処理システム |
US16/795,828 US20200294189A1 (en) | 2019-03-13 | 2020-02-20 | Information processing device, information processing method, recording medium storing program code, and information processing system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019045867A JP2020146204A (ja) | 2019-03-13 | 2019-03-13 | 情報処理装置、情報処理方法、プログラムおよび情報処理システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2020146204A true JP2020146204A (ja) | 2020-09-17 |
Family
ID=72423471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019045867A Withdrawn JP2020146204A (ja) | 2019-03-13 | 2019-03-13 | 情報処理装置、情報処理方法、プログラムおよび情報処理システム |
Country Status (2)
Country | Link |
---|---|
US (1) | US20200294189A1 (ja) |
JP (1) | JP2020146204A (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6962188B2 (ja) * | 2017-12-28 | 2021-11-05 | 株式会社リコー | 情報処理装置、生体信号計測システム、及びプログラム |
US10818386B2 (en) * | 2018-11-21 | 2020-10-27 | Enlitic, Inc. | Multi-label heat map generating system |
JP2020146286A (ja) * | 2019-03-14 | 2020-09-17 | 株式会社リコー | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム |
KR20220073209A (ko) * | 2020-11-26 | 2022-06-03 | 조선대학교산학협력단 | 심전도 신호로부터 설명 가능한 히트맵 영상에 근거한 질병을 판별하는 시스템 및 판별 방법 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0461843A (ja) * | 1990-06-30 | 1992-02-27 | Shimadzu Corp | 脳磁計測装置 |
JPH06254094A (ja) * | 1993-03-08 | 1994-09-13 | Toshiba Corp | 超音波診断装置 |
JP2015116435A (ja) * | 2013-12-20 | 2015-06-25 | キヤノン株式会社 | 放射線画像表示装置及び放射線画像表示方法 |
-
2019
- 2019-03-13 JP JP2019045867A patent/JP2020146204A/ja not_active Withdrawn
-
2020
- 2020-02-20 US US16/795,828 patent/US20200294189A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0461843A (ja) * | 1990-06-30 | 1992-02-27 | Shimadzu Corp | 脳磁計測装置 |
JPH06254094A (ja) * | 1993-03-08 | 1994-09-13 | Toshiba Corp | 超音波診断装置 |
JP2015116435A (ja) * | 2013-12-20 | 2015-06-25 | キヤノン株式会社 | 放射線画像表示装置及び放射線画像表示方法 |
Also Published As
Publication number | Publication date |
---|---|
US20200294189A1 (en) | 2020-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108471942B (zh) | 定量心脏测试 | |
US11037349B2 (en) | Information displaying system and non-transitory recording medium | |
JP2020146204A (ja) | 情報処理装置、情報処理方法、プログラムおよび情報処理システム | |
US10679394B2 (en) | Information processing device, information processing method, computer program product, and biosignal measurement system | |
US11666262B2 (en) | Information display device, biological signal measurement system and computer-readable recording medium | |
JP6996203B2 (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
JP7009906B2 (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
US20180325483A1 (en) | Information processing apparatus, information processing method, computer-readable medium, and biological signal measurement system | |
US20200297231A1 (en) | Information processing device, information processing method, recording medium storing program code, and biomedical-signal measuring system | |
US11237712B2 (en) | Information processing device, biomedical-signal measuring system, display method, and recording medium storing program code | |
US11109810B2 (en) | Information display device, biological signal measurement system, and computer program product | |
US11457856B2 (en) | Information processing device, information processing method, recording medium storing program code, and biomedical-signal measuring system | |
US11138779B2 (en) | Information processing apparatus, information processing method, computer-readable medium, and biological signal measurement system | |
US11484268B2 (en) | Biological signal analysis device, biological signal measurement system, and computer-readable medium | |
JP7135845B2 (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
JP7176197B2 (ja) | 情報処理装置、生体信号計測システム、表示方法、及びプログラム | |
US11229389B2 (en) | Information processing device, biomedical-signal measuring system, and recording medium storing program code | |
JP2021145969A (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
JP2021145875A (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
JP2019155074A (ja) | 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム | |
JP2021019954A (ja) | 情報処理装置、生体信号表示装置およびプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211223 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220810 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220823 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20221111 |