(57)【要約】【課題】注目する生体信号についての強度分布の領域を容易に抽出することができる情報処理装置、情報処理方法、プログラムおよび生体信号計測システムを提供する。【解決手段】表示部に表示される、特定の発信源からの生体信号の第１強度分布の特定情報を取得する取得部と、前記取得部により取得された前記特定情報に基づいて、前記発信源を示す画像に重畳させる前記生体信号の第２強度分布を求めるための前記第１強度分布の指定領域を抽出する抽出部と、を備える。【選択図】図３
2020146206.tif 000002
【特許請求の範囲】
【請求項１】
表示部に表示される、特定の発信源からの生体信号の第１強度分布の特定情報を取得する取得部と、
前記取得部により取得された前記特定情報に基づいて、前記発信源を示す画像に重畳させる前記生体信号の第２強度分布を求めるための前記第１強度分布の指定領域を抽出する抽出部と、
を備えた情報処理装置。
【請求項２】
前記第１強度分布は、前記生体信号の時刻および周波数についての二次元分布であり、
前記取得部は、前記特定情報として少なくとも前記二次元分布の特定の時刻および周波数を取得する請求項１に記載の情報処理装置。
【請求項３】
前記第１強度分布は、前記生体信号の信号強度の時刻および周波数についての二次元分布であり、
前記取得部は、前記特定情報として、さらに前記特定の時刻および周波数における前記信号強度の値を取得する請求項２に記載の情報処理装置。
【請求項４】
前記第１強度分布は、前記生体信号の信号強度の確度、分散または標準偏差の時刻および周波数についての二次元分布であり、
前記取得部は、前記特定情報として、さらに前記特定の時刻および周波数における前記信号強度の確度、分散または標準偏差の値を取得する請求項２に記載の情報処理装置。
【請求項５】
前記取得部により取得された前記特定情報としての前記値に基づいて、前記第１強度分布を二値化した二値化画像を得る二値化部を、さらに備え、
前記抽出部は、前記二値化画像における連結した領域であって、前記特定の時刻および周波数に対応する位置を含む領域を、前記指定領域として抽出する請求項３または４に記載の情報処理装置。
【請求項６】
前記二値化部は、前記特定情報としての前記値を閾値として、前記第１強度分布の各画素値に対する閾値判定により前記二値化画像を得る請求項５に記載の情報処理装置。
【請求項７】
前記二値化部は、前記特定情報としての前記値の所定の割合を閾値として、前記第１強度分布の各画素値に対する閾値判定により前記二値化画像を得る請求項５に記載の情報処理装置。
【請求項８】
前記取得部は、入力部からの操作入力により指定された前記第１強度分布での位置に対応する前記特定情報を取得する請求項１〜７のいずれか一項に記載の情報処理装置。
【請求項９】
前記取得部は、前記第１強度分布における予め定められた所定領域において前記特定情報を取得する請求項１〜７のいずれか一項に記載の情報処理装置。
【請求項１０】
前記取得部は、前記第１強度分布における前記所定領域の代表点に対応する前記特定情報を取得する請求項９に記載の情報処理装置。
【請求項１１】
前記抽出部は、前記取得部により取得された前記特定情報に基づいて、前記第２強度分布を求めるために抽出した前記第１強度分布の領域に外接する矩形領域を前記指定領域とする請求項１〜１０のいずれか一項に記載の情報処理装置。
【請求項１２】
前記取得部は、前記第１強度分布において複数の前記特定情報を取得し、
前記抽出部は、前記各特定情報に基づいて、前記第１強度分布の前記指定領域をそれぞれ抽出し、
前記抽出部により抽出された前記各指定領域に基づく前記第２強度分布それぞれを、別々の前記発信源を示す画像に重畳させて前記表示部に表示させる表示部制御部を、さらに備えた請求項１〜１１のいずれか一項に記載の情報処理装置。
【請求項１３】
前記取得部は、前記第１強度分布において複数の前記特定情報を取得し、
前記抽出部は、前記各特定情報に基づいて、前記第１強度分布の前記指定領域をそれぞれ抽出し、
前記抽出部により抽出された前記各指定領域のうち、入力部からの操作入力により選択された前記指定領域に基づく前記第２強度分布を、前記発信源を示す画像に重畳させて前記表示部に表示させる表示部制御部を、さらに備えた請求項１〜１１のいずれか一項に記載の情報処理装置。
【請求項１４】
表示部に表示される、特定の発信源からの生体信号の第１強度分布の特定情報を取得する取得ステップと、
取得した前記特定情報に基づいて、前記発信源を示す画像に重畳させる前記生体信号の第２強度分布を求めるための前記第１強度分布の指定領域を抽出する抽出ステップと、
を有する情報処理方法。
【請求項１５】
コンピュータに
表示部に表示される、特定の発信源からの生体信号の第１強度分布の特定情報を取得する取得ステップと、
取得した前記特定情報に基づいて、前記発信源を示す画像に重畳させる前記生体信号の第２強度分布を求めるための前記第１強度分布の指定領域を抽出する抽出ステップと、
を実行させるためのプログラム。
【請求項１６】
被験者から１種類以上の前記生体信号を測定する測定装置と、
請求項１〜１３のいずれか一項に記載の情報処理装置と、
を備えた生体信号計測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、情報処理方法、プログラムおよび生体信号計測システムに関する。
【背景技術】
【０００２】
脳外科手術等を行うためには、切除対象となる脳の疾患部位であるターゲット部位と、切除せずに温存すべき温存部位とを特定することが必要となる。温存部位としては、視覚野、聴覚野、体性感覚野、運動野および言語野等がある。これらの温存部位を誤って切除してしまうと、人間の感覚および運動等が害されることになるので、脳外科手術等に際してのターゲット部位および温存部位の特定は非常に重要である。このような脳外科手術等のために、事前に脳の活動を調べるためには脳磁計、脳波計、ｆＭＲＩ（ｆｕｎｃｔｉｏｎａｌ  Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）、またはｆＮＩＲＳ（ｆｕｎｃｔｉｏｎａｌ  Ｎｅａｒ−Ｉｎｆｒａｒｅｄ  Ｓｐｅｃｔｒｏｓｃｏｐｙ）等を用いて脳内の物理現象を計測する。このうち、ｆＭＲＩおよびｆＮＩＲＳは、脳内を流れる血流を測定して生体信号を得るものであるが、血流という性質上、脳の活動を計測するのに精度的に制限がある。一方、脳磁計は、脳内の電気的活動により生じる磁場を測定し、脳波計は、脳内の電気的活動を測定して波形として生体信号を測定できる。このような生体信号は、例えば時間周波数解析が行われて表示される。
【０００３】
このような、時間周波数解析を用いて脳活動の変化を測定する技術として、加算平均法では解析不能な脳活動の変化を解析する技術が開示されている（特許文献１参照）。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１に記載された技術では、脳等の生体の活動が行われている時間（時刻）および周波数等についての生体信号の強度分布の領域を容易に抽出することができないという問題がある。
【０００５】
本発明は、上記に鑑みてなされたものであって、注目する生体信号についての強度分布の領域を容易に抽出することができる情報処理装置、情報処理方法、プログラムおよび生体信号計測システムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
上述した課題を解決し、目的を達成するために、本発明は、表示部に表示される、特定の発信源からの生体信号の第１強度分布の特定情報を取得する取得部と、前記取得部により取得された前記特定情報に基づいて、前記発信源を示す画像に重畳させる前記生体信号の第２強度分布を求めるための前記第１強度分布の指定領域を抽出する抽出部と、を備えたことを特徴とする。
【発明の効果】
【０００７】
本発明によれば、注目する生体信号についての強度分布の領域を容易に指定することができる。
【図面の簡単な説明】
【０００８】
【図１】図１は、実施形態に係る生体信号計測システムの概略図である。
【図２】図２は、実施形態に係る情報処理装置のハードウェア構成の一例を示す図である。
【図３】図３は、実施形態に係る情報処理装置の機能ブロック構成の一例を示す図である。
【図４】図４は、実施形態に係る情報処理装置に表示される開始画面の一例を示す図である。
【図５】図５は、測定収集画面の一例を示す図である。
【図６】図６は、測定収集画面の左側の領域の拡大図である。
【図７】図７は、測定収集画面の右側の領域の拡大図である。
【図８】図８は、アノテーション情報が入力された直後の状態を示す図である。
【図９】図９は、更新されたアノテーションリストの図である。
【図１０】図１０は、実施形態に係る情報処理装置における測定収集動作の一例を示すフローチャートである。
【図１１】図１１は、時間周波数解析画面の一例を示す図である。
【図１２】図１２は、レンジをデシベルで表記したヒートマップの一例を示す図である。
【図１３】図１３は、ヒートマップにおいて特定の位置が指定された状態の一例を示す図である。
【図１４】図１４は、ヒートマップにおいてピークリストの３つのピークの位置を表示した状態の一例を示す図である。
【図１５】図１５は、ヒートマップにおいてピークの情報に応じて表示形態を変えた状態の一例を示す図である。
【図１６】図１６は、ヒートマップにおいて特定の範囲が指定された状態の一例を示す図である。
【図１７】図１７は、領域指定処理の一例を示すフローチャートである。
【図１８】図１８は、ヒートマップの一点を選択した状態を示す図である。
【図１９】図１９は、二値化された画像の一例を示す図である。
【図２０】図２０は、連続した領域を指定領域とすることを説明する図である。
【図２１】図２１は、クロージング処理を説明する図である。
【図２２】図２２は、連続した領域に外接する矩形領域を指定領域とすることを説明する図である。
【図２３】図２３は、ヒートマップにおいて特定の範囲が複数指定された状態の一例を示す図である。
【図２４】図２４は、時間周波数解析画面において立体図および頭部三面図が追加された状態の一例を示す図である。
【図２５】図２５は、予め定まった領域で指定領域を抽出することを説明する図である。
【図２６】図２６は、時間周波数解析画面における立体図の一例を示す図である。
【図２７】図２７は、ヒートマップで指定された位置に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。
【図２８】図２８は、ヒートマップで指定された範囲に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。
【図２９】図２９は、立体図に表示された脳がヒートマップ上のどの時刻・周波数に対応するのかを線分で示した状態の一例を示す図である。
【図３０】図３０は、立体図に表示された脳がヒートマップ上のどの時刻・周波数の範囲に対応するのかを矩形領域で示した状態の一例を示す図である。
【図３１】図３１は、立体図上をドラッグすることによって立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。
【図３２】図３２は、立体図上のいずれかの脳の表示をクリックすることによって、立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。
【図３３】図３３は、立体図に表示されたいずれかの脳の視点を変更した場合、同じ行のすべての脳の視点を変更する状態の一例を示す図である。
【図３４】図３４は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の一例を示す図である。
【図３５】図３５は、図３４における視点の変更動作を具体的に説明する図である。
【図３６】図３６は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の別の例を示す図である。
【図３７】図３７は、図３６における視点の変更動作を具体的に説明する図である。
【図３８】図３８は、立体図にコメントを付加した状態の一例を示す図である。
【図３９】図３９は、時間周波数解析画面における頭部三面図の一例を示す図である。
【図４０】図４０は、頭部三面図の立体画像としてカットモデルを表示した状態の一例を示す図である。
【図４１】図４１は、頭部三面図においてピークリストで選択されたピークの位置を示した状態の一例を示す図である。
【図４２】図４２は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を示した状態の一例を示す図である。
【図４３】図４３は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を色を変えて示した状態の一例を示す図である。
【図４４】図４４は、頭部三面図の立体画像にダイポール推定の結果を重畳表示した状態の一例を示す図である。
【図４５】図４５は、頭部三面図の立体画像に複数の測定対象の結果（ヒートマップ）を重畳表示した状態の一例を示す図である。
【図４６】図４６は、頭部三面図の立体画像に対して視点が変更される前の状態の一例を示す図である。
【図４７】図４７は、頭部三面図の立体画像に対して視点が変更された場合に表示されるダイアログボックスを示す図である。
【図４８】図４８は、立体画像の視点変更により立体図の１行目に同じ表示を反映するための設定例を示す図である。
【図４９】図４９は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１行目に同じ表示を反映した状態を示す図である。
【図５０】図５０は、立体画像の視点変更により立体図の１、２行目に同じ変換を反映するための設定例を示す図である。
【図５１】図５１は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１、２行目に同じ変換を反映した状態を示す図である。
【図５２】図５２は、立体画像の視点変更により立体図の１、２行目に、対応した変換を反映するための設定例を示す図である。
【図５３】図５３は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１、２行目に、対応した変換を反映した状態を示す図である。
【図５４】図５４は、立体画像の視点変更により立体図に同じ表示の新たな行を追加するための設定例を示す図である。
【図５５】図５５は、頭部三面図の立体画像に対して視点が変更された場合、立体図に同じ表示の新たな行を追加した状態を示す図である。
【図５６】図５６は、ピークリストの設定例を示す図である。
【図５７】図５７は、空間的なピークを説明する図である。
【図５８】図５８は、時間・周波数的なピークを説明する図である。
【図５９】図５９は、プルダウンされたピークリストから特定のピークを選択する状態を示す図である。
【図６０】図６０は、プルダウンされたピークリストから選択されたピークについてヒートマップ、立体図および頭部三面図に反映させた状態を示す図である。
【図６１】図６１は、再生制御パネルの操作によりヒートマップおよび立体図が再生表示される状態を示す図である。
【図６２】図６２は、再生制御パネルの操作によりヒートマップおよび立体図がコマ戻しされる状態を示す図である。
【図６３】図６３は、再生制御パネルの操作によりヒートマップおよび立体図がコマ送りされる状態を示す図である。
【図６４】図６４は、ピークに対してどの視点からの図を初期表示させるかを説明する図である。
【図６５】図６５は、２つのピークに対してどの視点からの図を初期表示させるかを説明する図である。
【図６６】図６６は、図６５に説明する視点からの図を立体図に初期表示させた状態を示す図である。
【図６７】図６７は、時系列で、腰椎の信号が上に伝わっていく状態の一例を示す図である。
【発明を実施するための形態】
【０００９】
以下に、図面を参照しながら、本発明に係る情報処理装置、情報処理方法、プログラムおよび生体信号計測システムの実施形態を詳細に説明する。また、以下の実施形態によって本発明が限定されるものではなく、以下の実施形態における構成要素には、当業者が容易に想到できるもの、実質的に同一のもの、およびいわゆる均等の範囲のものが含まれる。さらに、以下の実施の形態の要旨を逸脱しない範囲で構成要素の種々の省略、置換、変更および組み合わせを行うことができる。
【００１０】
（生体信号計測システムの概略）
図１は、実施形態に係る生体信号計測システムの概略図である。図１を参照しながら、本実施形態に係る生体信号計測システム１の概略について説明する。
【００１１】
生体信号計測システム１は、被検者の複数種類の生体信号（例えば、脳磁（ＭＥＧ：Ｍａｇｎｅｔｏ−ｅｎｃｅｐｈａｌｏｇｒａｐｈｙ）信号、および脳波（ＥＥＧ：Ｅｌｅｃｔｒｏ−ｅｎｃｅｐｈａｌｏｇｒａｐｈｙ）信号）を計測し、表示するシステムである。なお、測定対象となる生体信号は、脳磁信号および脳波信号に限られるものではなく、例えば、心臓の活動に応じて発生する電気信号（心電図として表現可能な電気信号）であってもよい。図１に示すように、生体信号計測システム１は、被検者の１以上の生体信号を測定する測定装置３と、測定装置３で測定された１種類以上の生体信号を記録するサーバ４０と、サーバ４０に記録された１種類以上の生体信号を解析する情報処理装置５０と、を含む。なお、図１では、サーバ４０と情報処理装置５０とが別々に記載されているが、例えば、サーバ４０が有する機能の少なくとも一部が情報処理装置５０に組み込まれる形態であってもよい。
【００１２】
図１の例では、被検者（被測定者）は、頭に脳波測定用の電極（またはセンサ）を付けた状態で測定テーブル４に仰向けで横たわり、測定装置３のデュワ３１の窪み３２に頭部を入れる。デュワ３１は、液体ヘリウムを用いた極低温環境の保持容器であり、デュワ３１の窪み３２の内側には脳磁測定用の多数の磁気センサが配置されている。測定装置３は、電極からの脳波信号と、磁気センサからの脳磁信号とを収集し、収集した脳波信号および脳磁信号を含むデータ（以下、「測定データ」と称する場合がある）をサーバ４０へ出力する。サーバ４０へ出力された測定データは、情報処理装置５０に読み出されて表示され、解析される。一般的に、磁気センサを内蔵するデュワ３１および測定テーブル４は、磁気シールドルーム内に配置されているが、図１では便宜上、磁気シールドルームの図示を省略している。
【００１３】
情報処理装置５０は、複数の磁気センサからの脳磁信号の波形と、複数の電極からの脳波信号の波形を、同じ時間軸上に同期させて表示する装置である。脳波信号とは、神経細胞の電気的な活動（シナプス伝達の際のニューロンの樹状突起で起きるイオン電荷の流れ）を電極間の電圧値として表される信号である。脳磁信号とは、脳の電気活動により生じた微小な電場変動を表す信号である。脳磁場は、高感度の超伝導量子干渉計（ＳＱＵＩＤ：Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ  Ｑｕａｎｔｕｍ  Ｉｎｔｅｒｆｅｒｅｎｃｅ  Ｄｅｖｉｃｅ）センサで検知される。これらの脳波信号および脳磁信号は、「生体信号」の一例である。
【００１４】
（情報処理装置のハードウェア構成）
図２は、実施形態に係る情報処理装置のハードウェア構成の一例を示す図である。図２を参照しながら、本実施形態に係る情報処理装置５０のハードウェア構成について説明する。
【００１５】
図２に示すように、情報処理装置５０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１０１と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）１０２と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）１０３と、補助記憶装置１０４と、ネットワークＩ／Ｆ１０５と、入力装置１０６と、表示装置１０７と、を有し、これらがバス１０８で相互に接続されている。
【００１６】
ＣＰＵ１０１は、情報処理装置５０の全体の動作を制御し、各種の情報処理を行う演算装置である。ＣＰＵ１０１は、ＲＯＭ１０３または補助記憶装置１０４に記憶された情報表示プログラムを実行して、測定収集画面および解析画面（時間周波数解析画面等）の表示動作を制御する。
【００１７】
ＲＡＭ１０２は、ＣＰＵ１０１のワークエリアとして用いられ、主要な制御パラメータおよび情報を記憶する揮発性の記憶装置である。ＲＯＭ１０３は、基本入出力プログラム等を記憶する不揮発性の記憶装置である。例えば、上述の情報表示プログラムがＲＯＭ１０３に記憶されているものとしてもよい。
【００１８】
補助記憶装置１０４は、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）またはＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）等の記憶装置である。補助記憶装置１０４は、例えば、情報処理装置５０の動作を制御する制御プログラム、ならびに、情報処理装置５０の動作に必要な各種のデータおよびファイル等を記憶する。
【００１９】
ネットワークＩ／Ｆ１０５は、サーバ４０等のネットワーク上の機器と通信を行うための通信インターフェースである。ネットワークＩ／Ｆ１０５は、例えば、ＴＣＰ（Ｔｒａｎｓｍｉｓｓｉｏｎ  Ｃｏｎｔｒｏｌ  Ｐｒｏｔｏｃｏｌ）／ＩＰ（Ｉｎｔｅｒｎｅｔ  Ｐｒｏｔｏｃｏｌ）に準拠したＮＩＣ（Ｎｅｔｗｏｒｋ  Ｉｎｔｅｒｆａｃｅ  Ｃａｒｄ）等によって実現される。
【００２０】
入力装置１０６は、タッチパネルの入力機能、キーボード、マウスおよび操作ボタン等のユーザインターフェース等である。表示装置１０７は、各種の情報を表示するディスプレイ装置である。表示装置１０７は、例えば、タッチパネルの表示機能、液晶ディスプレイ（ＬＣＤ：Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）または有機ＥＬ（Ｅｌｅｃｔｒｏ−Ｌｕｍｉｎｅｓｃｅｎｃｅ）等によって実現される。表示装置１０７は、測定収集画面および解析画面を表示し、入力装置１０６を介した入出力操作に応じて画面が更新される。
【００２１】
なお、図２に示す情報処理装置５０のハードウェア構成は一例であり、これ以外の装置が備えられるものとしてもよい。また、図２に示す情報処理装置５０は、例えば、ＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）を想定したハードウェア構成であるが、これに限定されるものではなく、タブレット等のモバイル端末であってもよい。この場合、ネットワークＩ／Ｆ１０５は、無線通信機能を有する通信インターフェースであればよい。
【００２２】
（情報処理装置の機能ブロック構成）
図３は、実施形態に係る情報処理装置の機能ブロック構成の一例を示す図である。図３を参照しながら、本実施形態に係る情報処理装置５０の機能ブロック構成について説明する。
【００２３】
図３に示すように、情報処理装置５０は、収集表示制御部２０１と、解析表示制御部２０２と、ピークリスト制御部２０３と、通信部２０４と、センサ情報取得部２０５と、解析部２０６と、記憶部２０７と、入力部２０８と、を有する。
【００２４】
収集表示制御部２０１は、センサ情報の収集動作時の画面表示を、図５〜図１０を参照して説明する手法で制御する機能部である。
【００２５】
解析表示制御部２０２は、センサ情報取得部２０５により取得されたセンサ情報（脳波信号または脳磁信号）から解析部２０６により算出された生体信号の信号強度等の画面表示を、図１１〜図６７を参照して説明する手法で制御する機能部である。解析表示制御部２０２は、図３に示すように、ヒートマップ表示制御部２１１と、立体表示制御部２１２と（表示制御部の一例）、断面表示制御部２１３（表示制御部の一例）と、再生表示制御部２１４と、を有する。
【００２６】
ヒートマップ表示制御部２１１は、後述する図１１等に示す時間周波数解析画面６０１のヒートマップ６１１の画面表示を制御する機能部である。ヒートマップ表示制御部２１１は、図３に示すように、取得部２２１と、二値化部２２２と、抽出部２２３と、を有する。
【００２７】
取得部２２１は、後述する図１１等に示すヒートマップ６１１において、ユーザが入力部２０８を介して指定した位置の時刻および周波数（特定情報の一例）と、当該時刻および周波数での生体信号の信号強度の値（特定情報の一例）と、を取得する機能部である。
【００２８】
二値化部２２２は、取得部２２１により取得された信号強度の値を用いて閾値判定することにより、後述するヒートマップ６１１の画像を二値化する機能部である。すなわち、二値化部２２２は、閾値判定の結果に基づいて、例えばヒートマップ６１１の各画素の画素値を０または１に設定することによって二値化する。
【００２９】
抽出部２２３は、二値化部２２２により得られた二値化画像から、取得部２２１により取得された時刻および周波数が示す位置の画素を含み、値（例えば１）を有する画素が連結した領域を抽出する機能部である。
【００３０】
立体表示制御部２１２は、時間周波数解析画面６０１の立体図６１２の画面表示を制御する機能部である。断面表示制御部２１３は、時間周波数解析画面６０１の頭部三面図６１３の画面表示を制御する機能部である。再生表示制御部２１４は、時間周波数解析画面６０１の再生制御パネル６１５に対する操作入力に従って、再生表示を制御する機能部である。
【００３１】
ピークリスト制御部２０３は、設定された条件を満たす信号強度のピークを抽出して、後述する図１１等に示す時間周波数解析画面６０１のピークリスト６１４に登録する機能部である。
【００３２】
通信部２０４は、測定装置３またはサーバ４０等とデータ通信を行う機能部である。通信部２０４は、図２に示すネットワークＩ／Ｆ１０５によって実現される。
【００３３】
センサ情報取得部２０５は、測定装置３またはサーバ４０から、通信部２０４を介して、センサ情報（脳波信号または脳磁信号）を取得する機能部である。解析部２０６は、センサ情報取得部２０５により取得されたセンサ情報（測定された信号）を解析し、脳内の各部における信号強度を示す信号（以下、当該信号についても「生体信号」と称する場合がある）を算出する機能部である。
【００３４】
記憶部２０７は、解析部２０６により算出された信号強度を示す生体信号のデータ等を記憶する機能部である。記憶部２０７は、図２に示すＲＡＭ１０２または補助記憶装置１０４によって実現される。
【００３５】
入力部２０８は、センサ情報に付加されるアノテーション情報の入力操作、および、時間周波数解析画面６０１に対する各種入力操作を受け付ける機能部である。入力部２０８は、図２に示す入力装置１０６によって実現される。
【００３６】
上述の収集表示制御部２０１、解析表示制御部２０２、ピークリスト制御部２０３、センサ情報取得部２０５および解析部２０６は、ＣＰＵ１０１がＲＯＭ１０３等に記憶されたプログラムをＲＡＭ１０２に展開して実行することにより実現される。なお、収集表示制御部２０１、解析表示制御部２０２、ピークリスト制御部２０３、センサ情報取得部２０５および解析部２０６の一部または全部は、ソフトウェアであるプログラムではなく、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）またはＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）等のハードウェア回路によって実現されてもよい。
【００３７】
また、図３に示した各機能部は、機能を概念的に示したものであって、このような構成に限定されるものではない。例えば、図３で独立した機能部として図示した複数の機能部を、１つの機能部として構成してもよい。一方、図３の１つの機能部が有する機能を複数に分割し、複数の機能部として構成するものとしてもよい。
【００３８】
（開始画面での動作）
図４は、実施形態に係る情報処理装置に表示される開始画面の一例を示す図である。図４を参照しながら、開始画面５０１での動作について説明する。
【００３９】
開始画面５０１には、「測定収集」および「解析」の選択ボタンが表示される。脳波測定および脳磁測定の場合、データの測定収集とデータの解析とでは、別々の主体によって行われる場合が多い。例えば、測定技師（測定者）によって「測定収集」のボタンが選択されると、測定装置３で測定されたデータは順次、サーバ４０に保存され、情報処理装置５０に読み出されて表示される。測定収集の終了後、医師によって「解析」のボタンが選択されると、収集された測定データが読み出されて解析される。
【００４０】
（測定収集時の動作）
図５は、測定収集画面の一例を示す図である。測定収集画面５０２は、図５に示すように、測定された生体信号（ここでは、脳磁信号および脳波信号）の信号波形を表示する領域５１１ａと、信号波形以外のモニタ情報を表示する領域５１１ｂと、を有する。信号波形を表示する領域５１１ａは、測定者からみて画面の左側に配置され、信号波形以外のモニタ情報を表示する領域５１１ｂは、測定者からみて画面の右側に配置されている。これによって、リアルタイムで検出され表示される波形の動き（画面の左側から右側に向かって表示される）に合わせた測定者の視線の動きと、画面左側の領域５１１ａから右側の領域５１１ｂへマウスを移動させるときの動きに無駄が生じず、作業効率が向上する。
【００４１】
表示画面の領域５１１ｂでは、測定中に被測定者の様子を確認するためのモニタウィンドウ５１２が表示される。測定中の被測定者のライブ映像を表示することで、後述するように、信号波形のチェックおよび判断の信頼性を高めることができる。図５では、１台のモニタディスプレイ（表示装置１０７）の表示画面に、測定収集画面５０２の全体が表示される場合を示しているが、左側の領域５１１ａおよび右側の領域５１１ｂを、２台またはそれ以上のモニタディスプレイに分けて別々に表示してもよい。
【００４２】
図６は、測定収集画面の左側の領域の拡大図である。領域５１１ａは、信号検出の時間情報を画面の水平方向に表示する第１表示領域５３０と、信号検出に基づく複数の信号波形を画面の垂直方向に並列に表示する第２表示領域５２１〜５２３と、を有する。
【００４３】
第１表示領域５３０で表示される時間情報は、図６の例では、時間軸５３１に沿って付された時間表示を含むタイムラインであるが、時間（数字）を表示せずに、帯状の軸だけでもよく、軸を設けずに時間（数字）の表示だけであってもよい。また、画面上側の第１表示領域５３０の他に、第２表示領域５２３の下側に時間軸５３１を表示して、タイムラインを表示してもよい。
【００４４】
領域５１１ａでは、同種の複数のセンサから取得される複数の信号波形、または複数種類のセンサ群から取得される複数種類の信号の波形が、同じ時間軸５３１で同期して表示される。図６に示す例では、第２表示領域５２１には被測定者の頭部右側から得られる複数の脳磁信号の波形が、第２表示領域５２２には被測定者の頭部左側から得られる複数の脳磁信号の波形が、それぞれ並列に表示されている。第２表示領域５２３には、複数の脳波信号の波形が並列に表示されている。これらの複数の脳波信号の波形は、各電極間で測定された電圧信号である。これらの複数の信号波形それぞれは、その信号が取得されたセンサの識別番号またはチャネル番号と対応付けられて表示されている。
【００４５】
測定が開始され各センサからの測定情報が収集されると、時間の経過と共に領域５１１ａの各第２表示領域５２１〜５２３の左端から右方向に向けて信号波形が表示される。ライン５３２は、計測の時刻（現在）を示しており、画面の左から右へ向けて移動する。領域５１１ａの右端（時間軸５３１の右端）まで信号波形が表示されると、その後は画面の左端から右へ向けて徐々に信号波形が消え、消えた位置に新しい信号波形が順次左から右方向に表示され、かつ、ライン５３２も左端から右へ向けて移動していく。これと共に、水平方向の第１表示領域５３０でも測定の進行に対応して、時間の経過が時間軸５３１上に表示される。測定収集は、終了ボタン５３９が押下されるまで継続される。
【００４６】
実施形態の特徴として、測定者（収集者）がデータ収集中に信号波形上で波形の乱れ、振幅の特異点等に気付いたときに、問題となる箇所または範囲を信号波形上でマークすることができる。マーキングの箇所または範囲は、マウスによるポインタ操作またはクリック操作で指定することができる。指定された箇所または範囲は、第２表示領域５２１〜５２３の信号波形上に強調表示されると共に、指定結果が対応する時刻位置または時間範囲で、第１表示領域５３０の時間軸５３１に沿って表示される。時間軸５３１上への表示を含むマーキングの情報は、信号波形データと共に保存される。指定された箇所は、ある時刻に対応し、指定された範囲は、ある時刻を含む一定範囲に対応する。
【００４７】
図６の例では、時刻ｔ１で、第２表示領域５２３で１以上のチャネルを含む範囲が指定され、マーク５２３ａ−１で時刻ｔ１を含んだ時間がハイライトで表示されている。マーク５２３ａ−１の表示と関連して第１表示領域５３０の対応する時刻位置に、指定結果を示すアノテーション５３０ａ−１が表示されている。時刻ｔ２で、第２表示領域５２３で別の波形位置またはその近傍がマークされ、その位置（時刻ｔ２）または近傍の領域（少なくとも時間範囲か複数の波形のいずれか一つが指示される）にマーク５２３ａ−２がハイライト表示されている。同時に、第１表示領域５３０の対応する時刻位置（時間範囲）に、アノテーション５３０ａ−２が表示される。ここで、アノテーションとは、あるデータに対して関連する情報を注釈として付与することを示す。本実施形態では、少なくとも指定された時間情報に基づいて注釈として表示されるものであって、少なくとも時間情報に基づく波形が表示される位置と結び付けて、注釈として表示するものである。また、複数のチャネルを表示する場合、対応するチャネル情報と結びつけ、注釈として表示してもよい。
【００４８】
時刻ｔ１で第１表示領域５３０に追加されたアノテーション５３０ａ−１は、一例として、アノテーション識別番号と、波形の属性を示す情報と、を含む。この例では、アノテーション番号「１」と共に、波形の属性を表わすアイコンと「ｓｔｒｏｎｇ  ｓｐｉｋｅ」（ストロングスパイク）というテキスト情報が表示されている。
【００４９】
時刻ｔ２で、測定者が別の波形箇所またはその近傍領域を指定すると、指定された箇所でマーク５２３ａ−２がハイライト表示され、これと共に、第１表示領域５３０の対応する時刻位置に、アノテーション番号「２」が表示される。さらに、ハイライト表示された箇所に、属性選択のためのポップアップウィンドウ５３５が表示される。ポップアップウィンドウ５３５は、種々の属性を選択する選択ボタン５３５ａと、コメントおよび追加情報を入力する入力ボックス５３５ｂと、を有する。選択ボタン５３５ａには、波形の属性として「速波（ｆａｓｔ  ａｃｔｉｖｉｔｙ）」、「眼球運動（ｅｙｅ  ｍｏｔｉｏｎ）」、「体動（ｂｏｄｙ  ｍｏｔｉｏｎ）」、「スパイク（ｓｐｉｋｅ）」等、波形乱れの要因が示されている。測定者は、画面の領域５１１ｂのモニタウィンドウ５１２で被測定者の様子を確認することができるので、波形の乱れの原因を示す属性を適切に選択することができる。例えば、波形にスパイクが生じたときに、てんかんの症状を示すスパイクなのか、被測定者の体動（くしゃみ等）に起因するスパイクなのかを判断することができる。
【００５０】
時刻ｔ１でも同じ操作が行われており、図６では、ポップアップウィンドウ５３５で「スパイク」の選択ボタン５３５ａが選択され、入力ボックス５３５ｂに「ｓｔｒｏｎｇ  ｓｐｉｋｅ」と入力されたことにより第１表示領域５３０にアノテーション５３０ａ−１が表示されている。このような表示態様により、同じ時間軸５３１上に多数の信号波形を同期して表示する際に、信号波形の注目箇所または範囲を視認により容易に特定することができ、かつ注目箇所の基本情報を容易に把握することができる。
【００５１】
なお、アノテーション５３０ａ−１の一部または全部、例えば、属性アイコンおよびテキスト情報の少なくとも一方を、第２表示領域５２３の信号波形上のマーク５２３ａ−１の近傍にも表示してもよい。ここで、信号波形上へのアノテーションの追加は、波形形状のチェックの妨げになる場合もあり得るので、第２表示領域５２１〜５２３の信号波形上にアノテーションを表示させる場合は、表示・非表示を選択可能にしておくことが望ましい。
【００５２】
カウンタボックス５３８は、スパイクのアノテーションの累積数を表示する。「スパイク」が選択される都度、カウンタボックス５３８のカウンタ値がインクリメントされ、収録開始から現在（ライン５３２）までのトータルのスパイク数が一目でわかるようになっている。
【００５３】
図７は、測定収集画面の右側の領域の拡大図である。図７は、図６と同じ時刻（ライン５３２の時点）での状態を示している。領域５１１ｂのモニタウィンドウ５１２では、頭部を測定装置３に入れて測定テーブル４に横たわっている被測定者の状態のライブ映像が表示される。領域５１１ｂでは、第２表示領域５２１、５２２、５２３の信号波形のそれぞれに対応する脳磁分布図５４１、５４２、脳波分布図５５０と、アノテーションリスト５６０と、が表示される。アノテーションリスト５６０は、図６の信号波形上でマークされたアノテーションの一覧である。第２表示領域５２１〜５２３で信号波形上の位置または範囲が指定されアノテーションが付される都度、対応する情報がアノテーションリスト５６０に順次追加される。測定収集画面５０２におけるアノテーションリスト５６０への追加および表示は、例えば、降順（新しいデータを上に表示）で行われるが、この例に限定されない。アノテーションリスト５６０の表示を昇順にしても構わないが、第１表示領域５３０で時間軸５３１に沿って表示されるアノテーションとの対応関係がわかるように表示する。さらに、表示順序を変更させたり、項目ごとにソートさせることも可能である。
【００５４】
図７の例では、アノテーションリスト５６０には、アノテーション番号「１」に対応する時刻情報と、付加されたアノテーション情報と、がリストされている。アノテーション情報として、「スパイク」を表わす属性アイコンと、「ｓｔｒｏｎｇ  ｓｐｉｋｅ」というテキストが記録されている。また、マーク５２３ａ−１がハイライト表示された時点で、アノテーション番号「２」に対応する時刻情報がリストされている。ここでは、「アノテーション」は、アノテーション番号と、時刻情報と、アノテーション情報との組であると考えてもよいし、アノテーション情報のみであると考えてもよいし、アノテーション情報と、アノテーション番号または時刻情報との組であると考えてもよい。
【００５５】
また、アノテーションリスト５６０の近傍に、表示／非表示の選択ボックス５６０ａが配置されている。選択ボックス５６０ａで非表示が選択されると、第２表示領域５２１〜５２３で、信号波形上のハイライトマーク以外のアノテーションが非表示にされるが、第１表示領域５３０の時間軸５３１に沿ったアノテーションの表示は維持される。これにより、信号波形の視認性を阻害せずにアノテーション情報を認識可能にする。
【００５６】
図８は、アノテーション情報が入力された直後の状態を示す図である。具体的には、図８は、時刻ｔ２でポップアップウィンドウ５３５の「スパイク」が選択され、「ｎｏｒｍａｌ  ｓｐｉｋｅ」（ノーマルスパイク）というテキストが入力された直後の画面を示す。図６で例示したポップアップウィンドウ５３５で「ＯＫ」ボタンが選択されると、ポップアップウィンドウ５３５が閉じて、図８に示すように第１表示領域５３０の対応する時刻位置にアノテーション５３０ａ−２が表示される。アノテーション番号「２」と対応付けて、「スパイク」を表わす属性アイコンと、「ｎｏｒｍａｌ  ｓｐｉｋｅ」のテキスト情報が表示される。これと同時に、カウンタボックス５３８の値がインクリメントされる。また、ハイライト表示されたマーク５２３ａ−２の近傍に、属性アイコン５２６−２が表示される。この例では、マーク５２３ａ−１の近傍にも属性アイコン５２６−１が表示されているが、上述したように、属性アイコン５２６−１、５２６−２の表示、非表示は選択可能である。マーク５２３ａ−１と属性アイコン５２６−１とを含むアノテーションＡ１、およびマーク５２３ａ−２と属性アイコン５２６−２とを含むアノテーションＡ２も、アノテーション情報に含まれる。
【００５７】
図９は、更新されたアノテーションリストの図である。測定収集画面５０２の画面左側の領域５１１ａでマーク５２３ａ−２に対応するアノテーションが付加されたことにより、アノテーションリスト５６０が更新される。アノテーション番号「２」に「ｎｏｒｍａｌ  ｓｐｉｋｅ」というメモが追加される。
【００５８】
以下同様に、測定中に領域５１１ａで信号波形上の特定の箇所または範囲が指定される都度、指定箇所が強調表示されると共に、第１表示領域５３０の時間軸５３１に沿ってアノテーション情報が表示される。領域５１１ｂでは、アノテーションリスト５６０にアノテーション情報が順次追加される。
【００５９】
アノテーションリスト５６０および信号波形表示の領域５１１ａにおいて、アノテーション番号の表示は必須ではなく、用いなくてもよい。付加されたアノテーションを識別できる情報であれば任意の情報を識別情報として用いることができる。例えば、属性アイコンと、属性文字列（「ｓｔｒｏｎｇ  ｓｐｉｋｅ」等）と、時間軸５３１の近傍での時刻とを対応付けて表示してもよい。さらに、領域５１１ａにファイル番号（図９の「Ｆｉｌｅ」の項目に表示される番号）を併記して表示してもよい。
【００６０】
終了ボタン５３９（図８参照）が選択（押下）され測定が終了すると、第２表示領域５２１〜５２３で指定されたハイライト箇所は信号波形に対応付けて保存される。第１表示領域５３０の対応する時刻位置に表示されたアノテーション情報も、アノテーション番号と時刻とに対応付けて保存される。カウンタボックス５３８のカウンタ値、および、アノテーションリスト５６０の内容等の関連情報も保存される。これらの表示情報を保存することで、測定者と解析者とが異なる場合でも、解析者は容易に問題となる箇所を認識し、解析することができる。
【００６１】
図１０は、実施形態に係る情報処理装置における測定収集動作の一例を示すフローチャートである。図１０を参照しながら、本実形態に係る情報処理装置５０での測定収集動作について説明する。
【００６２】
図４に示す開始画面５０１で「測定収集」が選択されると（ステップＳ１１）、測定が開始され、複数の信号の波形が同じ時間軸に沿って同期して表示される（ステップＳ１２）。ここで「複数の信号波形」という場合は、同一種類の複数のセンサで検知された信号波形、および、異なる種類のセンサのそれぞれで検知された複数の信号波形の両方を含む。この例では、複数の生体信号の波形は、被測定者の頭部右側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の頭部左側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の脳波測定用の電極から得られる脳波信号の波形と、から構成されるが、これに限られるものではない。なお、センサの選択については、左／右のセンサ群にとらわれず、頭頂部、前頭葉、側頭葉等のパートから、任意に選択可能である。例えば、図７等に示す「ＭＥＧ  Ｗｉｎｄｏｗ  Ｃｏｎｔｒｏｌ  １」で頭頂部のセンサが選択された場合は、「ＭＥＧ  Ｗｉｎｄｏｗ  Ｃｏｎｔｒｏｌ  ２」でそれ以外のセンサが選択された状態となる。
【００６３】
情報処理装置５０は、表示されている信号波形上で注目箇所または範囲が指定されたか否かを判断する（ステップＳ１３）。注目箇所または範囲の指定があると（ステップＳ１３：Ｙｅｓ）、信号波形の表示領域（第２表示領域５２１〜５２３）に指定箇所を強調表示すると共に、時間軸領域（第１表示領域５３０）の対応する時刻位置に指定結果を表示する（ステップＳ１４）。指定結果には、指定がなされたこと自体を示す情報、または指定の識別情報が含まれる。時間軸領域への指定結果の表示と同時に、または前後して、アノテーションの入力要求の有無を判断する（ステップＳ１５）。アノテーション入力の要求がある場合は（ステップＳ１５：Ｙｅｓ）、入力されたアノテーション情報を時間軸領域の対応する時刻位置に表示すると共に、アノテーションリストに追加する（ステップＳ１６）。その後、測定終了コマンドが入力されたか否かを判断する（ステップＳ１７）。注目位置（領域）の指定がない場合（ステップＳ１３：Ｎｏ）、および、アノテーションの入力要求がない場合（ステップＳ１５：Ｎｏ）は、ステップＳ１７へ移行して測定終了の判断を行う。測定が終了するまで（ステップＳ１７：Ｙｅｓ）、ステップＳ１３〜Ｓ１６を繰り返す。
【００６４】
この情報表示方法により、複数のセンサからの信号を収集する際に、信号情報の視認性の高い測定収集画面５０２が提供される。
【００６５】
（時間周波数解析画面における解析動作）
図１１は、時間周波数解析画面の一例を示す図である。図１１を参照しながら、情報処理装置５０に表示される時間周波数解析画面６０１における解析動作について説明する。
【００６６】
上述の図４に示す開始画面５０１において「解析」ボタンが押下されると、解析部２０６は、上述の測定収集画面５０２での測定収集動作により収集されたセンサ情報（脳波信号または脳磁信号）を解析し、脳（発信源の一例）内の各位置における信号強度を示す生体信号を算出する。信号強度を算出する方法としては、例えば、公知の方法である空間フィルタ法等が挙げられるが、その他の方法を用いるものとしてもよい。解析表示制御部２０２は、上述の図４に示す開始画面５０１において、「解析」ボタンが選択されると、表示装置１０７に図１１に示す時間周波数解析画面６０１を表示させる。時間周波数解析画面６０１は、図１１に示すように、解析画面切替リスト６０５と、ヒートマップ６１１と、立体図６１２と、頭部三面図６１３と、ピークリスト６１４と、再生制御パネル６１５と、を表示する。時間周波数解析画面６０１を用いた解析および測定の主目的は、視覚野、聴覚野、体性感覚野、運動野および言語野等の人間の生存に欠かせない部位を特定して表示することである。ピークリスト６１４の右側に表示されているピークリスト設定ボタン６１４ａは、ピークリスト６１４に登録されるピークの条件を設定するためのウィンドウを表示させるためのボタンである。ピークリスト設定ボタン６１４ａの押下による上述の設定については後述する。また、ヒートマップ６１１、立体図６１２、頭部三面図６１３、ピークリスト６１４および再生制御パネル６１５の表示および動作の詳細は、順次後述する。
【００６７】
解析画面切替リスト６０５は、各種解析画面を選択するためのリストである。解析画面切替リスト６０５から選択できる解析画面としては、本実施形態の生体信号に対する時間および周波数での解析をするための時間周波数解析画面６０１の他、生体信号からてんかん等の部位を推定および解析を行うためのダイポール推定を行うための解析画面等がある。本実施形態では時間周波数解析画面６０１での解析動作について説明する。
【００６８】
＜ヒートマップについて＞
図１２は、レンジをデシベルで表記したヒートマップの一例を示す図である。図１３は、ヒートマップにおいて特定の位置が指定された状態の一例を示す図である。図１４は、ヒートマップにおいてピークリストの３つのピークの位置を表示した状態の一例を示す図である。図１５は、ヒートマップにおいてピークの情報に応じて表示形態を変えた状態の一例を示す図である。図１２〜図１５を参照しながら、時間周波数解析画面６０１のヒートマップ６１１についての動作について説明する。
【００６９】
図１１に示すように、ヒートマップ６１１は、解析部２０６により算出された脳内の各位置の信号強度を示す生体信号について時間周波数分解を施し、横軸を時刻（トリガ時刻からの時間）、縦軸を周波数、そして、時刻および周波数で特定される生体信号の信号強度の二次元分布（第１強度分布）を色で表した図である。すなわち、ヒートマップ６１１は、各画素について信号強度の強弱を示すための色（例えばＲＧＢ）を示す画素値により色を表示しているが、実質的には、各画素には生体信号の信号強度が対応付けられている。図１１に示す例では、信号強度は、例えば、所定の基準に対する増減で示されている。ここで所定の基準とは、例えば、被験者に何も刺激を与えていない場合の信号強度の平均値を０［％］とする。
【００７０】
なお、本実施形態では、０±１００［％］で図示しているが、これに限られず、１００［％］以上の場合には、例えば２００［％］等、レンジを変えて表示してもよい。また、単位［％］に代えて、例えば図１２に示すヒートマップ６１１のように、［ｄｂ］（デシベル）で表記してもよい。また、ヒートマップ６１１は、例えば、時刻０［ｍｓ］に被験者に対して何らかの刺激（物理的に衝撃を与える、腕を動かす、言葉を発する、音を聴かせる等）を与えた場合、その後の時刻では刺激を与えた後の脳の活動状態を示し、時刻０［ｍｓ］以前の時刻では、刺激を与える前の脳の活動状態を示すことになる。
【００７１】
このヒートマップ６１１の表示動作は、ヒートマップ表示制御部２１１によって制御される。
【００７２】
図１３に示すように、解析者による入力部２０８に対する操作入力（クリック操作またはタップ操作）によって、ヒートマップ６１１上の特定の位置（点）を指定することができる。図１３に示すように、ヒートマップ表示制御部２１１は、指定された位置を、例えば、指定部６２１のように表示する。なお、指定部６２１は、白抜きの矩形で示しているが、これに限定されるものではなく、その他の表示態様で示すものとしてもよい。
【００７３】
また、図１３に示すヒートマップ６１１では、入力部２０８に対する操作入力によって指定された位置を表示しているが、ピークリスト６１４に登録されたピークのうち選択されたピークの時刻および周波数でのヒートマップおよびそのピーク位置が表示されるものとしてもよい。また、例えば、ピークリスト６１４で選択されたピークから数えて上位Ｎ個のピークの位置がヒートマップ６１１に表示されるものとしてもよい。図１４では、上位３個のピークの位置が表示された例を示す。なお、設定によりピークの位置をどのように表示させるかが決定されてもよく、例えば、上述の他、ピークそのものを表示させない設定、信号強度がＭ以上のピークを表示させる設定等を切り替えられるものとしてもよい。
【００７４】
また、図１５に示すように、ヒートマップ６１１上に表示された複数のピークについて、当該ピークの属性情報に応じて、表示の形態を変えるものとしてもよい。図１５では、表示されたピークごとに番号を付加し、番号の表示部分の色をそれぞれ異なる色となるように変えて表示させた例を示している。
【００７５】
上述のように、ヒートマップ６１１において特定の位置が指定されると、指定された位置に対応する時刻および周波数に対応する生体信号の信号強度の分布が、例えば図１１に示すように、立体図６１２（発信源を示す画像の一例）の脳の図では部位７１２ａ−１〜７１２ａ−５、７１２ｂ−１〜７１２ｂ−５のように、頭部三面図６１３（発信源を示す画像の一例）の脳の図では部位７１３ａ−１、７１３ａ−２、７１３ｂ、７１３ｃ、７１３ｄのようにヒートマップ（ヒートマップ６１１におけるヒートマップとは異なる）（第２強度分布）として表示される。具体的には、ヒートマップ６１１において指定された位置に対応する時刻および周波数の生体信号の信号強度が、赤〜青のヒートマップとして表示される。
【００７６】
＜ヒートマップ上での領域指定＞
図１６は、ヒートマップにおいて特定の範囲が指定された状態の一例を示す図である。図１７は、領域指定処理の一例を示すフローチャートである。図１８は、ヒートマップの一点を選択した状態を示す図である。図１９は、二値化された画像の一例を示す図である。図２０は、連続した領域を指定領域とすることを説明する図である。図２１は、クロージング処理を説明する図である。図２２は、連続した領域に外接する矩形領域を指定領域とすることを説明する図である。図２３は、ヒートマップにおいて特定の範囲が複数指定された状態の一例を示す図である。図２４は、時間周波数解析画面において立体図および頭部三面図が追加された状態の一例を示す図である。図２５は、予め定まった領域で指定領域を抽出することを説明する図である。図１６〜図２５を参照しながら、ヒートマップ上での範囲を指定する処理について説明する。
【００７７】
上述の図１３〜図１５では、ヒートマップ６１１の特定の位置（点）を指定することによって、当該位置に対応する時刻および周波数に対応する生体信号の信号強度の分布が、立体図６１２および頭部三面図６１３にヒートマップとして表示されることを説明した。さらに、ヒートマップ６１１上で時刻および周波数についての特定の範囲を指定し、当該範囲に対応する生体信号の信号強度の分布を、ヒートマップとして表示させたいニーズがある。しかし、例えば、解析者による入力部２０８に対するドラッグ操作またはスワイプ操作等によって矩形領域等の範囲を指定することができるものとしても、信号強度の分布を示すヒートマップとして観察したい領域以外の領域も含まれてしまう。すなわち、例えば、図１６に示すように、ヒートマップ６１１において生体信号の信号強度が強い領域としての指定領域６２２のような領域を指定したいというニーズがある。指定領域６２２は、同じ脳活動により強い生体信号が発生している時刻および周波数の範囲であると期待される。しかし、図１６に示す指定領域６２２のような領域を、解析者が手作業で指定することは非常に煩わしく困難である。そこで、本実施形態では、図１６に示す指定領域６２２のような領域を自動で抽出する処理について説明する。具体的には、図１６に示すように、解析者による入力部２０８に対する操作入力（クリック操作またはタップ操作）によって、ヒートマップ６１１上の特定の位置（点）が指定部６２１として指定されると、当該指定部６２１を含み、かつ所定の条件を満たす領域としての指定領域６２２を自動的に抽出処理について、図１７に示すフローチャートの流れに沿って説明する。図１６のような指定領域６２２が自動で抽出されると、指定領域６２２に含まれる時刻および周波数に対応する生体信号の信号強度の平均の分布が、例えば図１１に示すように、立体図６１２（発信源を示す画像の一例）の脳の図では部位７１２ａ−１〜７１２ａ−５、７１２ｂ−１〜７１２ｂ−５のように、頭部三面図６１３（発信源を示す画像の一例）の脳の図では部位７１３ａ−１、７１３ａ−２、７１３ｂ、７１３ｃ、７１３ｄのようにヒートマップ（ヒートマップ６１１におけるヒートマップとは異なる）（第２強度分布）として表示されることになる。
【００７８】
＜＜ステップＳ２１＞＞
解析者は、例えば、図１８に示すヒートマップ６１１が拡大されたヒートマップにおいて、入力部２０８に対する操作入力（クリック操作またはタップ操作）によって、特定の位置（画素）である指定部６２１を指定したものとする。図１８に示すヒートマップにおいては、各画素の画素値である生体信号の信号強度が強いほど、密度の濃いハッチングを施した画素で示している。すると、取得部２２１は、指定部６２１の位置の時刻および周波数と、当該時刻および周波数での生体信号の信号強度の値と、を取得する。そして、ステップＳ２２へ移行する。
【００７９】
＜＜ステップＳ２２＞＞
二値化部２２２は、取得部２２１により取得された信号強度の値を用いて閾値判定することにより、ヒートマップ６１１の画像を二値化する。例えば、二値化部２２２は、取得部２２１により取得された信号強度の値以上の信号強度を有する画素に対応する画素値を「１」（黒）、取得部２２１により取得された信号強度の値未満の信号強度を有する画素に対応する画素値を「０」（白）とした図１９に示すような二値化画像を生成する。
【００８０】
なお、上述では、二値化部２２２によって、取得部２２１により取得された信号強度の値を閾値として閾値判定するものとしたが、これに限定されるものではない。例えば、取得部２２１により取得された信号強度の値の所定割合（例えば８０％）以上であるか否かにより閾値判定することにより、「０」、「１」の値を振り分けて二値化画像を得るものとしてもよい。この場合、所定割合については、スライドバー等で容易に変更できることが望ましい。また、取得部２２１により取得された時刻および周波数で特定される位置の周辺で探索された最も信号強度が大きい画素の当該信号強度についての所定割合（例えば７０％）を用いて閾値判定することにより、「０」、「１」の値を振り分けて二値化画像を得るものとしてもよい。また、信号強度に対する固定の閾値を用いて閾値判定をすることにより、「０」、「１」の値を振り分けて二値化画像を得るものとしてもよい。ただし、この場合、取得部２２１により取得された信号強度の値が０の場合、ダイアログボックスを表示して、他の時刻および周波数の指定を促す等の例外処理を行う必要がある。また、信号強度の値ではなく、別途、統計的に算出することが可能な信号強度についての確度（尤もらしさ）、分散または標準偏差について閾値判定をすることにより、「０」、「１」の値を振り分けて二値化画像を得るものとしてもよい。この場合、ヒートマップ６１１で表示するのは信号強度のヒートマップではなく、確度（尤もらしさ）、分散または標準偏差のヒートマップとして表示するものとすればよい。
【００８１】
そして、ステップＳ２３へ移行する。
【００８２】
＜＜ステップＳ２３＞＞
抽出部２２３は、二値化部２２２により得られた二値化画像から、取得部２２１により取得された時刻および周波数が示す位置の画素を含み、かつ、画素値が「１」である画素が連結した領域を、図２０に示す指定領域６２２として抽出する。
【００８３】
以上のステップＳ２１〜Ｓ２３の領域指定処理によって、図１６に示すように、解析者によって指定された指定部６２１から自動的に指定領域６２２を抽出することができる。すなわち、解析者によりヒートマップ６１１上の１点をクリック操作等により指定するのみで、自動的に複雑な形状を呈した指定領域６２２を指定することが可能となる。
【００８４】
なお、抽出部２２３は、二値化部２２２により得られた二値化画像において、図２１に示すクロージング処理を施した上で、指定領域を抽出するものとしてもよい。例えば、図２１（ａ）に示す画像が二値化部２２２により得られた二値化画像であるものとする。クロージング処理では、まず、図２１（ａ）に示す二値化画像に対して膨張処理を施すことによって図２１（ｂ）に示すような画像を得る。この図２１（ｂ）に示した画像に対して収縮処理を施することによって、図２１（ｃ）に示すような画像を得る。これによって、元の二値化画像における小さなギャップを介した連結領域同士を繋げることができ、二値化部２２２により二値化を行う際の閾値の設定の精度が不十分な場合でも、適切に指定領域を抽出することが可能となる。なお、抽出部２２３は、クロージング処理の代わりに、またはクロージング処理に加えて、収縮処理を行った後に膨張処理を行うオープニング処理を施するものとしてもよい。これによって、例えば、指定領域に含まれることを望まない細かい連結領域を排除することができる。
【００８５】
また、図２０に示すように、抽出部２２３により複雑な形状をした指定領域６２２を抽出することに限定されるものではなく、例えば、図２２に示すように、抽出部２２３は、上述の処理により求めた指定領域６２２に相当する領域に外接する矩形領域を、指定領域６２２ａとして抽出するものとしてもよい。より簡便には、抽出部２２３は、ヒートマップ６１１上でクリック操作等により指定された位置から上下左右に信号強度の値を確認し、指定された位置の信号強度の値の所定割合（例えば８０％）以下の値になるところまでの矩形範囲を指定領域として抽出するものとしてもよい。
【００８６】
また、図２３に示すように、解析者による入力部２０８に対する追加の操作（例えば、新たなクリック操作またはタップ操作）によって、既に指定された指定領域６２２の他、指定部６２１ａを指定することにより、当該指定部６２１ａに対する上述の図１７に示す領域指定処理により、追加で指定領域６２３を指定することができる。この場合、新たに指定された指定領域６２３に対応する立体図および頭部三面図として、図２４に示すように、それぞれ立体図６１２ａおよび頭部三面図６１３ａが表示される。そして、指定領域６２３に含まれる時刻および周波数に対応する生体信号の信号強度の平均の分布が、立体図６１２ａおよび頭部三面図６１３ａの脳の図にヒートマップ（ヒートマップ６１１におけるヒートマップとは異なる）として表示される。ヒートマップ６１１で複数指定した情報を受け取ると、この情報を受け取った順にこれに対応する立体図６１２および頭部三面図６１３が上から下に表示される。図２４の例では、指定領域６２２、６２３の順に選択された例を示す。このように表示することで、直感的にわかりやすくなる。ここで、ヒートマップ６１１で複数選択した情報を受け取ると、この情報を受け取った順にこれに対応する立体図６１２および頭部三面図６１３が下から上に表示される形態としてもよい。この場合、ヒートマップ６１１の直下に、最新で指定された選択領域に対応する立体図６１２および頭部三面図６１３が表示されるため、ヒートマップ６１１、立体図６１２、頭部三面図６１３への視線の移動を少なくできる。なお、ヒートマップ６１１上で複数の位置を指定する場合、上述の指定領域６２２、６２３のような範囲に限らず、指定部６２１のように点位置を複数指定することも可能である。このように、ヒートマップ６１１上で複数の位置（点または領域）が指定されることによって、その指定された時刻および周波数に対応する生体信号の信号強度の分布を比較することができる。
【００８７】
また、例えば、ヒートマップ６１１の任意の位置をクリック操作またはタップ操作等がされたことをトリガとして、予め決められた閾値を用いて、二値化部２２２により二値化された画像から、抽出部２２３によりすべての連結領域がそれぞれ指定領域として抽出されるものとしてもよい。この場合、抽出された各指定領域に対応する立体図６１２、６１２ａ、・・・、および頭部三面図６１３、６１３ａ、・・・にそれぞれヒートマップが重畳された状態で表示されるものとすればよい。または、抽出部２２３により抽出された複数の指定領域から解析者により選択操作された指定領域に対応するヒートマップを、立体図６１２および頭部三面図６１３に表示させるものとしてもよい。この場合、抽出部２２３により複数の指定領域が抽出された直後の状態としては、例えば、指定領域内の信号強度の平均が最も大きいものに対応するヒートマップを、立体図６１２および頭部三面図６１３に表示させるものとしてもよい。
【００８８】
また、人間の特定の行動（例えば、手を握る等）に応じてヒートマップ６１１において着目すべき時刻および周波数の領域が予め定まっている場合がある。例えば、対象となる人間の行動が運動である場合には、図２５に示すヒートマップ６１１における予め定められた領域６２５〜６２７に注目して解析を行う場合が多い。この場合、二値化部２２２は、領域６２５〜６２７について二値化を行い、各領域にそれぞれ対応する二値化画像から連結した領域を抽出して指定領域を求めるものとしてもよい。この場合、領域６２５〜６２７の具体的な位置は、例えば、ヒートマップ６１１における時刻０の位置から定められる。ただし、脳活動においては、ある程度、領域６２５〜６２７の範囲が分かっていても、個人差が存在する。そこで、予め大雑把な範囲は設定した上で、その範囲内で時刻および周波数について精査することが好ましい。そこで、例えば、領域６２５〜６２７の中心位置（代表点の一例）を予め定めておき、当該各中心位置がそれぞれクリック操作またはタップ操作されたものとして、上述の図１７に示す領域指定処理により、各中心位置に対応する指定領域をそれぞれ抽出することができる。これによって、個人差を吸収してより適切な指定領域を抽出することができる。
【００８９】
＜立体図について＞
図２６は、時間周波数解析画面における立体図の一例を示す図である。図２７は、ヒートマップで指定された位置に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図２８は、ヒートマップで指定された範囲に対応する脳の状態を立体図において中央に表示させた状態の一例を示す図である。図２９は、立体図に表示された脳がヒートマップ上のどの時刻・周波数に対応するのかを線分で示した状態の一例を示す図である。図３０は、立体図に表示された脳がヒートマップ上のどの時刻・周波数の範囲に対応するのかを矩形領域で示した状態の一例を示す図である。図３１は、立体図上をドラッグすることによって立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図３２は、立体図上のいずれかの脳の表示をクリックすることによって、立体図の表示、およびヒートマップ上の矩形領域の表示が移動する状態の一例を示す図である。図２６〜図３２を参照しながら、時間周波数解析画面６０１の立体図６１２の基本的な表示動作について説明する。
【００９０】
図２６に示すように、立体図６１２は、所定の視点からの脳の立体画像（３Ｄ画像）を表示する図であり、ヒートマップ６１１上で指定された位置（点または範囲）、または、ピークリスト６１４で選択されたピークの位置に対応する生体信号の信号強度をヒートマップ（第２強度分布）として重畳表示する。図２６に示すように、立体図６１２では、同じ行では同じ視点からの脳の立体画像が表示される。図２６に示す例では、立体図６１２の上側の行である表示領域６１２−１に表示された脳の立体画像は、脳の左側面からの視点での画像であり、下側の行である表示領域６１２−２に表示された脳の立体画像は、脳の右側面からの視点での画像である。この立体図６１２の表示動作は、立体表示制御部２１２によって制御される。
【００９１】
なお、図２６に示す立体図６１２では２つの視点から見た場合の脳の立体図、すなわち、２行で構成された脳の立体図であるが、これに限定されるものではなく、その他の数の行で表示するものとしてもよく、行の数を設定により変更できるものとしてもよい。例えば、脳の言語野の測定であれば、左右の際が重要な情報となるので、脳の左側面および右側面からの２つの視点からの脳の立体図を表示（２行で表示）するものとすればよい。ここで、測定対象と視点の対応の一例を下記の（表１）に示す。測定対象は、測定時に被験者へ与える刺激（刺激装置による付与であり（表１）のＮｏ．１〜４に対応する）や、被験者の動作（（表１）のＮｏ．５）を示し、収集時に測定収集画面上で選択する項目である。この測定対象を選択することで、これに対応する視点の脳立体図が表示される。また、視点は、被験者の正面を起点としたときの方向を示す。なお、行数についても別途設定することができるものとしてもよい。また、図２６に示す立体図は、（表１）におけるＮｏ．２に相当する。以降の説明では、便宜上、立体図６１２は、２つの行（２つの視点）で構成されているものとして説明する。
【表１】
2020146206.tif 000003
【００９２】
また、図２７に示すように、ヒートマップ６１１上で指定部６２１が指定されると、立体表示制御部２１２は、指定部６２１に対応する時刻を立体図６１２の表示領域の中心として、前後の時刻に対応した脳における信号強度のヒートマップを立体図６１２に表示する。図２７の例では、ヒートマップ６１１上で時刻５６０［ｍｓ］が指定されたため、立体図６１２において５６０［ｍｓ］を中心として、５５０、５５５、５６０、５６５、５７０［ｍｓ］のように、隣り合う脳の表示時刻の間隔を５［ｍｓ］おきで表示している。ただし、隣り合う脳の表示時刻の間隔は、その他１０［ｍｓ］、２５［ｍｓ］等、設定で変更できるようにしてもよい。
【００９３】
なお、図２８に示すように、ヒートマップ６１１上で範囲（指定領域６２２）が指定された場合、当該範囲内の信号強度を平均した信号強度のヒートマップを立体図６１２に表示するものとすればよい。この場合、選択された時刻の範囲に応じて、立体図６１２で表示される隣り合う立体画像の時刻を調整設定するものとすればよい。例えば、図２８に示すように、指定領域６２２の時刻の範囲が４５０〜６００［ｍｓ］の場合、立体図６１２で表示させる隣り合う立体画像の表示時刻の間隔を１５０［ｍｓ］に設定したとき、立体図６１２の中央に表示される立体画像は４５０〜６００［ｍｓ］に対応する。立体図６１２の中央の立体画像の左の立体画像は３００〜４５０［ｍｓ］に対応する図とし、右の立体画像は、６００〜７５０［ｍｓ］に対応する図とすればよい。また、各立体画像に表示されるヒートマップは、例えば、各時間幅における平均値とすればよい。
【００９４】
次に、図２９および図３０を参照しながら、ヒートマップ６１１上の位置または範囲と、立体図６１２の各立体画像との対応付けについて説明する。このような対応付けとしては、まず、図２９に示すように、ヒートマップ６１１上で特定の点位置である指定部６２１−１が指定された場合、その指定部６２１−１の時刻および周波数に対応する脳の立体画像が立体図６１２に表示される。この場合、その脳の立体画像を中心として前後の時刻に対応した脳の立体画像が表示されるので（図２９に示す例では５つの時刻に対応する脳）、ヒートマップ表示制御部２１１は、それぞれの脳の時刻に対応する点を、ヒートマップ６１１上で対応部６２１−２〜６２１−５としてそれぞれ表示させる。この場合、対応部６２１−２〜６２１−５の周波数位置は、指定部６２１−１の周波数位置に合わせる。さらに、ヒートマップ表示制御部２１１は、図２９に示すように、ヒートマップ６１１上の指定部６２１−１および対応部６２１−２〜６２１−５と、立体図６１２上の対応する脳の立体画像とをそれぞれ結ぶ線分６３１−１〜６３１−５で表示する。これによって、立体図６１２上で表示されている脳の状態が、ヒートマップ６１１上のどの位置に対応する状態であるのかを一目で確認することができる。なお、図２９に示す例では、線分を用いているが、その他の対応付けの方法として、指定部６２１−１および対応部６２１−２〜６２１−５のマークと、立体図６１２の脳の背景色とを合わせる等、対応付けが確認できる方法であればよい。また、この場合、解析者が指定した指定部６２１−１は、対応部６２１−２〜６２１−５と区別する態様で表示するものとすればよい。
【００９５】
また、ヒートマップ６１１上で特定の範囲である指定領域６２２−１が指定された場合、まず、図３０に示すように、その指定領域６２２−１の時刻および周波数に対応する脳の立体画像が立体図６１２に表示される。このとき、指定領域６２２−１は、例えば、抽出部２２３により抽出された連結領域に外接する矩形領域であるものとする。この場合、その脳の立体画像を中心として前後の時刻範囲に対応した脳の立体画像が表示されるので（図３０に示す例では５つの時刻範囲に対応する脳）、ヒートマップ表示制御部２１１は、それぞれの脳の時刻範囲に対応する範囲を、ヒートマップ６１１上で対応領域６２２−２〜６２２−５としてそれぞれ表示させる。また、この場合、解析者が指定した指定領域６２２−１は、対応領域６２２−２〜６２２−５と区別する態様で表示するものとすればよく、例えば、指定領域６２２−１の矩形の枠の色だけ異なる色となるように表示するものとすればよい。さらに、立体表示制御部２１２は、図３０に示すように、ヒートマップ６１１上の指定領域６２２−１および対応領域６２２−２〜６２２−５と同様な矩形を、立体図６１２上で対応する脳の立体画像を囲うように表示させる。これによって、立体図６１２上で表示されている脳の状態が、ヒートマップ６１１上のどの範囲に対応する状態であるのかを一目で確認することができる。なお、ヒートマップ６１１上で特定の範囲である指定領域６２２−１が指定された場合、指定領域６２２−１および対応領域６２２−２〜６２２−５の枠を表示すると共に、立体図６１２上で、枠７２２−１〜７２２−５およびヒートマップを表示するものとしてもよい。
【００９６】
次に、図３１および図３２を参照しながら、立体図６１２上でドラッグ操作、スワイプ操作またはカーソルキー操作をした場合に、立体図６１２の表示を左右に移動させる動作について説明する。図３１では、立体図６１２上でドラッグ操作、スワイプ操作またはカーソルキー操作により、立体図６１２の立体画像を右側に移動させた状態の例を示す。この場合、図３１に示すように、移動した結果、表示されている脳に応じて時刻表示も更新され、立体図６１２で中央に表示された脳の立体図が選択された状態となるように矩形が表示される。さらに、立体表示制御部２１２は、立体図６１２における立体画像に移動に従って、ヒートマップ６１１上の指定領域６２２−１および対応領域６２２−２〜６２２−５の表示も移動させる。
【００９７】
また、図３２に示すように、解析者による操作によって、立体図６１２上の中央の立体画像以外の立体画像がクリック操作またはタップ操作されると、操作された脳の立体画像が立体図６１２の中心に移動する。なお、実装上は脳の画像はそのままで、オーバーラップするヒートマップのみが移動するものとしてもよい。この場合、図３２に示すように、移動した結果、表示されている脳に応じて時刻表示も更新され、立体図６１２で中央に表示された脳の立体図が選択された状態となるように矩形が表示される。さらに、立体表示制御部２１２は、立体図６１２における立体画像に移動に従って、ヒートマップ６１１上の指定領域６２２−１および対応領域６２２−２〜６２２−５の表示も移動させる。
【００９８】
以上のように、立体図６１２における表示を自由に移動することが可能となることにより、時間の前後で脳の状態の変化を素早く確認することができる。
【００９９】
図３３は、立体図に表示されたいずれかの脳の視点を変更した場合、同じ行のすべての脳の視点を変更する状態の一例を示す図である。図３４は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の一例を示す図である。図３５は、図３４における視点の変更動作を具体的に説明する図である。図３６は、立体図に表示されたいずれかの脳の視点を変更した場合、すべての行のすべての脳の視点を変更する状態の別の例を示す図である。図３７は、図３６における視点の変更動作を具体的に説明する図である。図３８は、立体図にコメントを付加した状態の一例を示す図である。図３３〜図３８を参照しながら、時間周波数解析画面６０１の立体図６１２の特定の立体画像について視点変更を行った場合の動作について説明する。
【０１００】
立体図６１２に立体画像として表示されている脳は、解析者による操作（例えば、ドラッグ操作またはスワイプ操作等）によって視点を変更することが可能である。この場合、立体図６１２の特定の脳の立体画像の視点を変更した場合において、他の立体画像への反映方法の類型について以下に説明する。
【０１０１】
まず、特定の立体画像の視点が変更された場合に、同じ行の他の立体画像について同じ視点変更がなされる場合について説明する。図３３（ａ）に示すように、解析者は、２行で表示された立体図６１２において、表示領域６１２−１に表示された脳の立体画像のうち右端の立体画像（以下、「対象立体画像」と称する場合がある）に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部２１２は、図３３（ｂ）に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の背面の視点の立体画像の表示となるように視点変更する。このとき、脳の画像に重畳しているヒートマップも同様に視点変更する。そして、立体表示制御部２１２は、図３３（ｃ）に示すように、対象立体画像と同じ行（表示領域６１２−１）の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。これによって、特定の立体画像（対象立体画像）に対する視点変更を行うのみで、同じ行の他の立体画像についてもその視点変更が反映されるので、操作性が向上し、同じ視点での時間的に前後の脳の活動変化を確認することが容易になる。なお、立体画像の視点変更の操作として、例えば、視点変更したい立体画像上にマウスを合わせてドラック操作またはクリック操作等でもよく、ポップアップ表示によるパラメータの指定等でもよい。
【０１０２】
次に、特定の立体画像の視点が変更された場合に、他のすべての立体画像について同じ視点変更がなされる場合について説明する。図３４（ａ）に示すように、解析者は、２行で表示された立体図６１２において、表示領域６１２−１に表示された脳の立体画像のうち右端の対象立体画像に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部２１２は、図３４（ｂ）に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の背面の視点の立体画像の表示となるように視点変更する。そして、立体表示制御部２１２は、図３４（ｃ）に示すように、対象立体画像と同じ行（表示領域６１２−１）の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図３５（ａ）に示す脳の左側面の視点を、表示領域６１２−１の他の立体画像について、図３５（ｂ）に示すように、脳の背面の視点となるように変更する。さらに、立体表示制御部２１２は、図３４（ｃ）に示すように、対象立体画像と異なる行（表示領域６１２−２）の立体画像に対しても、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図３５（ａ）に示す脳の右側面の視点を、表示領域６１２−２の立体画像について、図３５（ｃ）に示すように、脳の前面の視点となるように変更する。なお、処理環境に余裕があれば図３５（ａ）〜図３５（ｃ）の処理を高速に行い、見た目上すべての脳の視点が同時に変わっているようにしてもよいし、処理環境が貧弱な場合には、ユーザが動かしている画像のみについて視点を変えた後、それが確定したタイミング（例えばドラッグ操作で脳画像を回転させることで視点を変更している場合、マウスのボタンを離したタイミング）で他の画像の視点を変えてもよい。このとき、脳の画像に重畳しているヒートマップもそれぞれ同様に視点を変更する。これによって、特定の立体画像（対象立体画像）に対する視点変更を行うのみで、同じ行の他の立体画像、および他の行の立体画像についてもその視点変更が反映されるので、操作性が向上し、時間的に前後の脳の活動変化を確認することが容易になる。
【０１０３】
そして、特定の立体画像の視点が変更された場合に、同じ行の他の立体画像については同じ視点変更がなされ、他の行の立体画像については、対応する視点変更（具体的には、脳の中心面（対象面）に対して面対称となる視点変更）がなされる場合について説明する。図３６（ａ）に示すように、解析者は、２行で表示された立体図６１２において、表示領域６１２−１に表示された脳の立体画像のうち右端の対象立体画像に対して、視点を変更する操作を行ったものとする。この場合、立体表示制御部２１２は、図３６（ｂ）に示すように、脳の左側面の視点の対象立体画像を、解析者の操作に従って、脳の左前側の視点の立体画像の表示となるように視点変更する。そして、立体表示制御部２１２は、図３６（ｃ）に示すように、対象立体画像と同じ行（表示領域６１２−１）の他の立体画像に対して、対象立体画像の視点変更と同じ視点変更をして表示する。すなわち、図３７（ａ）に示す脳の左側面の視点を、表示領域６１２−１の他の立体画像について、図３７（ｂ）に示すように、脳の左前側の視点となるように変更する。さらに、立体表示制御部２１２は、図３６（ｃ）に示すように、対象立体画像と異なる行（表示領域６１２−２）の立体画像に対して、対象立体画像の視点変更に対応する視点変更をして表示する。すなわち、図３７（ａ）に示す脳の右側面の視点を、表示領域６１２−２の立体画像について、図３７（ｃ）に示すように、脳の中心面（対象面）に対して面対称となるような視点変更、すなわち、脳の右前側の視点となるように変更する。このとき、脳の画像に重畳しているヒートマップもそれぞれ同様に視点を変更する。これによって、特定の立体画像（対象立体画像）に対する視点変更を行うのみで、同じ行の他の立体画像についてはその視点変更が反映され、他の行の立体画像については、対応する視点変更が反映されるので、操作性が向上し、さらに、複数の行を見比べることによって、対応する視点から時間的に前後の脳の活動変化を確認することができる。
【０１０４】
以上で説明した３つの他の立体画像への反映方法のうち、いずれの方法を採用してもよく、または、これらの方法のうちどの方法で反映するかを設定により切り替えることができるものとしてもよい。
【０１０５】
なお、図３３〜図３７で上述したように、解析者が最初に視点変更の操作の対象とする対象立体画像を、表示領域６１２−１の右端の立体画像としたが、これに限定されるものではなく、表示領域６１２−１または表示領域６１２−２のいずれの立体画像が操作の対象であってもよい。また、表示領域６１２−１に含まれる一群の立体画像、および表示領域６１２−２に含まれる一群の立体画像が、それぞれ本発明の「形状画像」または「第３画像」に相当する。
【０１０６】
また、立体図６１２の特定の脳の立体画像の視点を変更した場合に、他の立体画像へ視点変更を反映する動作を説明したが、立体画像に対して変更を加える表示態様について視点だけではなく、例えば、拡大縮小、輝度の変更、または透明度の変更等であってもよい。これらの変更についても、上述の視点変更の趣旨を逸脱しない範囲で、他の立体画像に変更を反映するものとすればよい。
【０１０７】
また、上述のように、立体図６１２に表示される立体画像に対して変更を加えた結果、図３８に示すように、解析者による入力部２０８への操作入力により、特定の立体画像に対してメモ（例えば、図３８に示すコメント６３５）を付加できるようにしてもよい。これによって、解析者（医者等）が脳の気になる活動部位に対するコメントをその立体画像に関連付けて保存しておくことができ、脳に対する外手術の他、当該脳の疾患についてのカンファレンスにおいても役立てることができる。
【０１０８】
＜頭部三面図について＞
図３９は、時間周波数解析画面における頭部三面図の一例を示す図である。図４０は、頭部三面図の立体画像としてカットモデルを表示した状態の一例を示す図である。図４１は、頭部三面図においてピークリストで選択されたピークの位置を示した状態の一例を示す図である。図４２は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を示した状態の一例を示す図である。図４３は、頭部三面図においてピークリストで選択されたピークの位置、および時間的な前後のピークの位置を色を変えて示した状態の一例を示す図である。図４４は、頭部三面図の立体画像にダイポール推定の結果を重畳表示した状態の一例を示す図である。図４５は、頭部三面図の立体画像に複数の測定対象の結果（ヒートマップ）を重畳表示した状態の一例を示す図である。図３９〜図４５を参照しながら、時間周波数解析画面６０１の頭部三面図６１３の基本的な表示動作について説明する。
【０１０９】
図３９に示すように、頭部三面図６１３は、脳の特定の位置（点）における３方向の断面図（以下、「三面図」と総称する場合がある）、および３Ｄ画像である立体画像６４４を含む図である。図３９に示す例では、頭部三面図６１３は、脳の特定の位置における３方向の断面図として、脳の前後方向に対して垂直な断面を示す断面図６４１、脳の左右方向に対して垂直な断面を示す断面図６４２、および、脳の上下方向に対して垂直な断面を示す断面図６４３を含む。断面図６４１は、上述の特定の位置を通るように基準線６４５ａおよび基準線６４５ｂが引かれている。また、断面図６４２は、上述の特定の位置を通るように、基準線６４５ａおよび基準線６４５ｃが引かれている。また、断面図６４３は、上述の特定の位置を通るように、基準線６４５ｂおよび基準線６４５ｄが引かれている。断面図６４１〜６４３には、それぞれ、ヒートマップ６１１で指定された位置（点または範囲）に対応する時刻および周波数の生体信号の信号強度の分布を示すヒートマップ（ヒートマップ６１１とは異なる）が重畳して表示される。この頭部三面図６１３の表示動作は、断面表示制御部２１３によって制御される。
【０１１０】
基準線６４５ａは、上述の脳の特定の位置の上下方向の位置を規定する線であるため、断面図６４１および断面図６４２にわたって連続した線として描画されている。また、基準線６４５ｂは、上述の脳の特定の位置の左右方向の位置を規定する線であるため、断面図６４１および断面図６４３にわたって連続した線として描画されている。また、基準線６４５ｃは、断面図６４２において、上述の脳の特定の位置の前後方向の位置を規定する線である。また、基準線６４５ｄは、断面図６４３において、上述の脳の特定の位置の前後方向の位置を規定する線である。なお、頭部三面図６１３における断面図６４１〜６４３の並びは、上述のように、基準線６４５ａおよび基準線６４５ｂが複数の断面図に連続して描画することが可能であるために、図３９に示す並びとしているが、これに限定されるものではなく、任意でよい。その場合、各断面図において、脳の特定の位置を通るように基準線がそれぞれ描画されるものとすればよい。また、基準線は必ずしも描画される必要はなく、各断面図において、脳の特定の位置を示すマーク等が表示されるものとしてもよい。
【０１１１】
立体画像６４４は、脳の３Ｄ画像であり、後述するように、この立体画像６４４に対する操作に応じて、立体図６１２に描画された脳の立体画像の視点が変更される。また、立体画像６４４には、ヒートマップ６１１で指定された位置（点または範囲）に対応する時刻および周波数の生体信号の信号強度の分布を示すヒートマップ（ヒートマップ６１１とは異なる）が重畳して表示される。なお、立体画像６４４は、脳の特定の視点の立体画像を表示することに限定されるものではなく、例えば、図４０に示すように、三面図で特定される脳の位置を中心に３次元方向で一部カットしたカットモデル画像であってもよい。
【０１１２】
また、図３９に示す頭部三面図６１３では、ピークリスト６１４に登録されたピークのうち選択されたピークの位置が特定される三面図を表示し、図４１に示すように、立体画像６４４上に当該選択されたピーク位置を示すピーク点６４６を表示するものとしてもよい。また、例えば、ピークリスト６１４で選択されたピークから数えて上位Ｎ個のピークの位置が立体画像６４４に表示されるものとしてもよい。図４２では、上位３個のピークの位置（ピーク点６４６、６４６ａ、６４６ｂ）が表示された例を示す。また、図４２では、上位３個ではなく、ピークリスト６１４で選択されたピークの前後の時刻におけるピーク位置をピーク点６４６、６４６ａ、６４６ｂとして表示（ピークの軌跡を表示）するものとしてもよい。なお、設定によりピークの位置をどのように表示させるかが決定されてもよく、例えば、上述の他、ピークそのものを表示させない設定、信号強度がＭ以上のピークを表示させる設定等を切り替えられるものとしてもよい。
【０１１３】
また、図４３に示すように、立体画像６４４上に表示された複数のピークについて、当該ピークの属性情報に応じて、表示の形態を変えるものとしてもよい。図４３では、表示されたピークごとにマークの色をそれぞれ異なる色となるように変えて表示させた例を示している。
【０１１４】
また、図４４に示すように、断面表示制御部２１３は、立体画像６４４に対して、他の解析画面等でダイポール推定の結果であるダイポール６４７を重畳表示させるものとしてもよい。これによって、立体画像６４４が示す温存部位を示すヒートマップと、疾患部位（ターゲット部位）を示すダイポールとの位置関係を把握することができ、外科手術等に役立てることができる。
【０１１５】
また、三面図のいずれかの断面図上で、解析者による入力部２０８を用いてクリック操作またはタップ操作によって、脳の三次元空間の特定の位置を指定することができる。このように、三面図において特定の位置が指定されると、指定された位置（点）に対応する時刻および周波数に関する生体信号の信号強度の分布がヒートマップ６１１に反映される。
【０１１６】
また、三面図のいずれかの断面図上で、解析者による入力部２０８に対するドラッグ操作またはスワイプ操作によって、脳の三次元空間の特定の範囲を指定することができる。このように、三面図において特定の範囲が指定されると、指定された範囲に対応する時刻および周波数に関する生体信号の信号強度（その範囲での平均）の分布がヒートマップ６１１に反映される。
【０１１７】
なお、立体画像６４４（および三面図）に描画されるヒートマップ（信号強度の強弱を示す等高線図）は、各脳の部位の活動が活発となるようなそれぞれの刺激の結果を重畳して表示するものとしてもよい。例えば、言語刺激（測定時）を行った信号に対する結果と、視覚刺激（測定時）を行った信号に対する結果を求めた後に、断面表示制御部２１３は、図４５（ａ）に示す言語野を活性化させた場合のヒートマップと、図４５（ｂ）に示す視覚野を活性化させた場合のヒートマップとを、図４５（ｃ）に示すように、重畳させて表示可能とするものとしてもよい。このように図４５（ｃ）に示すように重畳されたヒートマップで示される部位が温存部位であるということが確認できる。重畳させるための操作方法としては、現在表示されている測定結果が言語野であるものとした場合、メニューによって異なる測定結果（例えば、視覚野）を選択することができるようにすればよい。なお、重畳する際には、測定対象より刺激に対する反応時間が異なる可能性もある。そこで、測定対象を追加する場合に時間のずれを設定できるようにすれば、より的確に重畳することができる。さらに、ヒートマップを重畳した結果である図４５（ｃ）に示す立体画像を、図４５（ｄ）に示すように反転表示させることによって、逆に温存部位ではない切除可能な部位を示すことが可能になる。
【０１１８】
また、頭部三面図６１３に示す断面図は、断面方向の異なる３つの断面図である三面図としているが、これに限定されるものではなく、特定の断面方向の１つの断面図、または断面方向の異なる２つもしくは４つ以上の断面図であってもよい。
【０１１９】
図３９は、頭部三面図の立体画像に対して視点が変更される前の状態の一例を示す図である。図４７は、頭部三面図の立体画像に対して視点が変更された場合に表示されるダイアログボックスを示す図である。図４８は、立体画像の視点変更により立体図の１行目に同じ表示を反映するための設定例を示す図である。図４９は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１行目に同じ表示を反映した状態を示す図である。図５０は、立体画像の視点変更により立体図の１、２行目に同じ変換を反映するための設定例を示す図である。図５１は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１、２行目に同じ変換を反映した状態を示す図である。図５２は、立体画像の視点変更により立体図の１、２行目に、対応した変換を反映するための設定例を示す図である。図５３は、頭部三面図の立体画像に対して視点が変更された場合、立体図の１、２行目に、対応した変換を反映した状態を示す図である。図５４は、立体画像の視点変更により立体図に同じ表示の新たな行を追加するための設定例を示す図である。図５５は、頭部三面図の立体画像に対して視点が変更された場合、立体図に同じ表示の新たな行を追加した状態を示す図である。図４６〜図５５を参照しながら、時間周波数解析画面６０１の頭部三面図６１３の立体画像に対して視点変更した場合における立体図６１２への反映動作について説明する。
【０１２０】
頭部三面図６１３の立体画像６４４に表示されている脳は、立体図６１２と同様に、解析者による操作（例えば、ドラッグ操作またはスワイプ操作等）によって視点を変更することが可能である。この場合、立体画像６４４の脳の視点を変更した場合において、立体図６１２に表示された脳の視点に反映するものとしてもよく、以下、その反映方法の類型について以下に説明する。
【０１２１】
図４６に示すような時間周波数解析画面６０１に表示された頭部三面図６１３の立体画像６４４に対して、解析者による操作（例えば、ドラッグ操作またはスワイプ操作等）が行われると、断面表示制御部２１３は、図４７に示すように、ダイアログボックス６５０を表示させる。ダイアログボックス６５０は、立体画像６４４の脳について視点変更された場合に、どのように立体図６１２に反映するかを決定するための画面である。ここで、例えば、「立体図は変更しない」ボタンが押下されると、立体図６１２の立体画像に対する視点変更はなされない。ここでは、解析者は、図４７に示すように、脳の左側面の視点の立体画像６４４を、脳の背面の視点の立体画像の表示となるように視点変更したものとする。
【０１２２】
まず、図４８に示すように、ダイアログボックス６５０の「立体図の行の表示に反映する」ボタンが押下されたものとする。すると、断面表示制御部２１３は、図４８に示すように、立体図６１２にどのように反映するかの詳細を設定するためのダイアログボックス６５１を表示させる。解析者は、図４８に示すように、ダイアログボックス６５１において、変更する行として立体図６１２の１行目に選択し、「立体図の３Ｄと同じ表示」を選択したものとする。この場合、立体表示制御部２１２は、図４９に示すように、立体図６１２の１行目（上の行）の立体画像に対して、立体画像６４４の変更された視点と同じ視点となるように表示する。
【０１２３】
次に、図４７のように視点変更された状態から、解析者は、図５０に示すように、ダイアログボックス６５０の「立体図の行の表示に反映する」ボタンを押下し、ダイアログボックス６５１において、変更する行として立体図６１２の１行目および２行目に選択し、「立体図の３Ｄと同じ変換」を選択したものとする。この場合、立体表示制御部２１２は、図５１に示すように、もともと立体画像６４４と同じ視点であった立体図６１２の１行目の立体画像に対して、立体画像６４４の視点変更と同じ視点変更をして表示する。すなわち、図５１に示すように、脳の背面の視点となるように変更する。さらに、立体表示制御部２１２は、図５１に示すように、もともと脳の右側面の視点であった立体図６１２の２行目の立体画像に対して、立体画像６４４の視点変更と同じ視点変更をして表示する。すなわち、図５１に示すように、脳の前面の視点となるように変更する。なお、「立体図の行の表示に反映する」ボタン押下によるダイアログボックス６５１での選択を初期設定済みとするか初期設定可能なようにしておき、その上で、例えば、［ビューのリンク］／［ビューのリンク解除］ボタンを設けることで、その選択結果を表示するようにしてもよい。このような構成とすることで毎回の選択動作を省略（簡略化）できる。
【０１２４】
次に、解析者は、図５３に示すように、脳の左側面の視点の立体画像６４４を、脳の左前側の視点の立体画像の表示となるように視点変更したものとする。そして、この状態から、解析者は、図５２に示すように、ダイアログボックス６５０の「立体図の行の表示に反映する」ボタンを押下し、ダイアログボックス６５１において、変更する行として立体図６１２の１行目および２行目に選択し、「立体図の３Ｄに対応した変換」を選択したものとする。この場合、立体表示制御部２１２は、図５３に示すように、もともと立体画像６４４と同じ視点であった立体図６１２の１行目の立体画像に対して、立体画像６４４の視点変更と同じ視点変更をして表示する。すなわち、図５３に示すように、脳の左前側の視点となるように変更する。さらに、立体表示制御部２１２は、図５３に示すように、もともと脳の右側面の視点であった立体図６１２の２行目の立体画像に対して、立体画像６４４の視点変更に対応する視点変更をして表示する。すなわち、図５３に示すように、脳の中心面（対象面）に対して面対称となるような視点変更、すなわち、脳の右前側の視点となるように変更する。
【０１２５】
次に、図４７のように視点変更された状態から、解析者は、図５４に示すように、ダイアログボックス６５０の「立体図の行を追加する」ボタンを押下することによって表示されるダイアログボックス６５２において、「立体図の３Ｄと同じ表示」を選択したものとする。この場合、立体表示制御部２１２は、図５５に示すように、立体画像６４４の視点変更の結果、同一の視点となる脳の立体画像を、立体図６１２の表示領域６１２−３に新たな行として表示する。すなわち、図５５に示すように、表示領域６１２−３における新たな行に、脳の背面の視点の立体画像が表示される。
【０１２６】
以上のように、頭部三面図６１３における立体画像６４４に対して行った視点変更を、各種設定に応じて、立体図６１２の時系列に並んだ脳の立体画像の視点に反映することができる。これによって、立体画像６４４の視点変更と同様の視点変更を、立体図６１２に対して改めて行う必要がないので、操作性が向上し、さらに、立体画像６４４で変更した視点と同一の視点、または対応する視点で、立体図６１２において脳の状態の変化を時系列に確認することができる。
【０１２７】
なお、図４８、図５０、図５２および図５４で示すダイアログボックス６５０〜６５２によって設定される立体図６１２に表示された脳の視点への反映方法は、一例であり、その他の反映方法が設定できるものとしてもよい。
【０１２８】
また、立体画像６４４の視点を変更した場合に、立体図６１２の立体画像へ視点変更を反映する動作を説明したが、立体画像６４４に対して変更を加える表示態様は視点だけではなく、例えば、拡大縮小、輝度の変更、または透明度の変更等であってもよい。これらの変更についても、上述の視点変更の趣旨を逸脱しない範囲で、立体図６１２の立体画像に変更を反映するものとすればよい。
【０１２９】
＜ピークリストについて＞
図５６は、ピークリストの設定例を示す図である。図５７は、空間的なピークを説明する図である。図５８は、時間・周波数的なピークを説明する図である。図５９は、プルダウンされたピークリストから特定のピークを選択する状態を示す図である。図６０は、プルダウンされたピークリストから選択されたピークについてヒートマップ、立体図および頭部三面図に反映させた状態を示す図である。図５６〜図６０を参照しながら、時間周波数解析画面６０１のピークリスト６１４の基本的な動作について説明する。
【０１３０】
ピークリスト６１４は、ピークリスト制御部２０３により抽出された、設定された条件を満たす信号強度のピークが登録されたリストである。図５６に示すように、ピークリスト制御部２０３は、ピークリスト６１４でプルダウンされることによって登録された信号強度の一覧であるプルダウンリスト６５６を表示する。
【０１３１】
また、上述のピークリスト制御部２０３により抽出される信号強度のピークの条件は、ピークリスト設定ボタン６１４ａを押下することにより設定することができる。ピークリスト制御部２０３は、ピークリスト設定ボタン６１４ａが押下されると、抽出される信号強度のピークの条件を設定するためのダイアログボックス６５５を表示させる。
【０１３２】
ダイアログボックス６５５では、まず、ピークリスト６１４に登録されたピーク情報をどのようにソートするかを設定することができる。ダイアログボックス６５５において、「ピークの値が大きい順番」が選択された場合、ピークリスト制御部２０３は、ピークリスト６１４で登録されているピーク情報を、ピークの信号強度が大きい順にソートする。一方、ダイアログボックス６５５において、「ピークの高さ（頂点と谷点の差）が大きい順番」が選択された場合、ピークリスト制御部２０３は、ピークリスト６１４で登録されているピーク情報を、ピーク点の信号強度と、当該ピークの谷部分の信号強度との差が大きい順にソートする。
【０１３３】
さらに、ダイアログボックス６５５では、ピークリスト６１４にどのようなピーク情報を登録（リストアップ）させるかを設定することができる。ダイアログボックス６５５において、「すべての空間的なピーク」が選択された場合、ピークリスト制御部２０３は、時刻・周波数平面の各時刻・各周波数においての脳全体における空間的なピークを抽出して、ピークリスト６１４に登録する。ここで、空間的なピークとは、図５７に示すピーク部８０１のように、着目する時刻・周波数の生体信号の、脳全体における信号強度のピークであり、このピーク部８０１の信号強度は、周辺より信号強度が強い。
【０１３４】
また、ダイアログボックス６５５において、「すべての時刻／周波数的なピーク」が選択された場合、ピークリスト制御部２０３は、脳全体の各位置においての時刻・周波数平面における時刻／周波数的なピークを抽出して、ピークリスト６１４に登録する。ここで、時刻／周波数的なピークとは、図５８に示すピーク部８０２のように、着目する脳の位置での生体信号の、時刻・周波数平面における信号強度のピークであり、このピーク部８０２の信号強度は、周辺よりも信号強度が強い。
【０１３５】
また、ダイアログボックス６５５において、「指定されている時刻／周波数における空間的なピーク」が選択された場合、ピークリスト制御部２０３は、時刻・周波数平面で指定された時刻・周波数においての脳全体における空間的なピークを抽出して、ピークリスト６１４に登録する。なお、指定されている時刻／周波数は一点とは限らず、範囲で選択される場合もある。
【０１３６】
また、ダイアログボックス６５５において、「指定されている位置における時刻／周波数的なピーク」が選択された場合、ピークリスト制御部２０３は、指定された脳の位置においての時刻・周波数平面における時刻／周波数的なピークを抽出して、ピークリスト６１４に登録する。なお、指定されている位置は一点とは限らず、範囲で選択される場合もある。例えば、視覚野についてのピークを抽出する場合、後頭部全体の範囲を指定することによって、所望のピークを抽出しやすくなる。
【０１３７】
次に、各ピーク情報が登録されたピークリスト６１４から特定のピーク情報を選択した場合の動作について説明する。解析者によって、ピークリスト６１４から表示されたプルダウンリスト６５６から特定のピーク情報（例えば、図５９に示す「９５％／９ｍｓ／７０Ｈｚ／ｖｏｘｅｌ：１７３６」）が選択されると、ヒートマップ表示制御部２１１は、当該ピーク情報が示す脳の特定の位置に対応するヒートマップ６１１を表示する。この場合、ヒートマップ表示制御部２１１は、図１４で上述したように、ヒートマップ６１１においてピーク情報が示すピークの位置を具体的に示すものとすればよい。
【０１３８】
また、立体表示制御部２１２は、選択されたピーク情報が示す時刻・周波数の脳の立体画像を立体図６１２の各行の中央に表示させ、さらに、その時刻の前後の脳の立体画像を表示させる。この場合、立体図６１２の脳の各立体画像に重畳されるヒートマップは、ピーク情報が示す周波数の生体信号の信号強度に対応するものとすればよい。
【０１３９】
また、断面表示制御部２１３は、選択されたピーク情報が示す脳の位置を通る三面図を頭部三面図６１３に表示させる。さらに、断面表示制御部２１３は、選択されたピーク情報が示す時刻・周波数の生体信号の信号強度に対応するヒートマップを立体画像６４４の脳に重畳表示させるものとすればよい。なお、断面表示制御部２１３は、図６０に示すように、選択されたピーク情報が示す脳の位置を中心に３次元方向で一部カットしたカットモデルを立体画像６４４に表示させるものとしてもよい。
【０１４０】
以上のように、ピークリスト６１４に登録されたピーク情報から特定のピーク情報を選択することによって、当該ピーク情報に対応したヒートマップ６１１、立体図６１２および頭部三面図６１３が表示される。これによって、選択したピークが脳のどの位置、およびどの時刻・周波数であるかを瞬時に認識することができ、さらに、ヒートマップ６１１においては当該ピークおよびその周りの時刻／周波数における信号強度の状態、および、当該ピークの位置およびその周りの脳の信号強度の状態も把握することが可能となる。
【０１４１】
＜再生制御パネルについて＞
図６１は、再生制御パネルの操作によりヒートマップおよび立体図が再生表示される状態を示す図である。図６２は、再生制御パネルの操作によりヒートマップおよび立体図がコマ戻しされる状態を示す図である。図６３は、再生制御パネルの操作によりヒートマップおよび立体図がコマ送りされる状態を示す図である。図６１〜図６３を参照しながら、時間周波数解析画面６０１の再生制御パネル６１５を操作した場合の動作について説明する。
【０１４２】
再生制御パネル６１５は、解析者の操作によって、時間経過と共に、ヒートマップ６１１、立体図６１２および頭部三面図６１３の状態を再生表示するためのユーザインターフェースである。
【０１４３】
例えば、解析者が、再生制御パネル６１５の「再生」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１に対して、図６１（ａ）および図６１（ｂ）に示すように、ヒートマップ６１１で指定されている指定領域６２２−１、およびその周りの対応領域６２２−２〜６２２−５を、時間経過と共に、右方向（時間が進む方向）に移動させるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１および対応領域６２２−２〜６２２−５の移動に伴い、立体表示制御部２１２に対して、図６１（ａ）および図６１（ｂ）に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、断面表示制御部２１３に対して、移動している指定領域６２２−１に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像６４４に表示させるように指示する。
【０１４４】
また、解析者が、再生制御パネル６１５の「コマ戻し」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１に対して、図６２（ａ）および図６２（ｂ）に示すように、ヒートマップ６１１で指定されている指定領域６２２−１、およびその周りの対応領域６２２−２〜６２２−５を、所定時間分だけ左方向（時間が戻る方向）に移動させるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１および対応領域６２２−２〜６２２−５の移動に伴い、立体表示制御部２１２に対して、図６２（ａ）および図６２（ｂ）に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、断面表示制御部２１３に対して、移動した指定領域６２２−１に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像６４４に表示させるように指示する。
【０１４５】
また、解析者が、再生制御パネル６１５の「コマ送り」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１に対して、図６３（ａ）および図６３（ｂ）に示すように、ヒートマップ６１１で指定されている指定領域６２２−１、およびその周りの対応領域６２２−２〜６２２−５を、所定時間分だけ右方向（時間が進む方向）に移動させるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１および対応領域６２２−２〜６２２−５の移動に伴い、立体表示制御部２１２に対して、図６３（ａ）および図６３（ｂ）に示すように、各領域に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、断面表示制御部２１３に対して、移動した指定領域６２２−１に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像６４４に表示させるように指示する。
【０１４６】
また、解析者が、再生制御パネル６１５の「停止」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１、立体表示制御部２１２および断面表示制御部２１３に対して、ヒートマップ６１１、立体図６１２、および頭部三面図６１３におけるそれぞれの再生表示動作を停止するように指示する。
【０１４７】
また、解析者が、再生制御パネル６１５の「先頭へ移動」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１に対して、ヒートマップ６１１で指定されている指定領域６２２−１を、時刻の先頭に移動させるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、立体表示制御部２１２に対して、指定領域６２２−１に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、断面表示制御部２１３に対して、移動した指定領域６２２−１に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像６４４に表示させるように指示する。
【０１４８】
また、解析者が、再生制御パネル６１５の「末尾へ移動」ボタンを押下すると、再生表示制御部２１４は、ヒートマップ表示制御部２１１に対して、ヒートマップ６１１で指定されている指定領域６２２−１を、時刻の末尾に移動させるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、立体表示制御部２１２に対して、指定領域６２２−１に対応した脳の立体画像の表示に切り替えるように指示する。また、再生表示制御部２１４は、ヒートマップ６１１における指定領域６２２−１の移動に伴い、断面表示制御部２１３に対して、移動した指定領域６２２−１に対応する時刻・周波数の範囲に対応する信号強度のヒートマップを、三面図および立体画像６４４に表示させるように指示する。
【０１４９】
以上のように、再生表示させることによって、立体図６１２および頭部三面図６１３に表示される信号強度の分布（ヒートマップ）の時間経過による変化を動画として確認することができ、時間経過によるピークの移動等を視覚的に確認することができる。
【０１５０】
＜初期表示について＞
図６４は、ピークに対してどの視点からの図を初期表示させるかを説明する図である。図６５は、２つのピークに対してどの視点からの図を初期表示させるかを説明する図である。図６６は、図６５に説明する視点からの図を立体図に初期表示させた状態を示す図である。図６４〜図６６を参照しながら、時間周波数解析画面６０１を起動した（開いた）場合におけるヒートマップ６１１、立体図６１２および頭部三面図６１３の初期表示について説明する。
【０１５１】
解析者が、時間周波数解析画面６０１を起動した（開いた）場合に、ヒートマップ６１１、立体図６１２および頭部三面図６１３の初期表示としてどのような画像を表示させるかについての類型を説明する。
【０１５２】
例えば、解析表示制御部２０２は、時刻・周波数全体を通じて、脳全体において信号強度が最大となる時刻・周波数および脳内の位置を求める。この場合、ヒートマップ表示制御部２１１は、解析表示制御部２０２により求められた脳内の位置でのヒートマップ６１１を表示させる。また、立体表示制御部２１２は、解析表示制御部２０２により求められた信号強度が最大となる時刻・周波数に対応する脳の立体画像を立体図６１２に表示させる。また、断面表示制御部２１３は、解析表示制御部２０２により求められた脳内の位置を通るような三面図を頭部三面図６１３に表示させ、解析表示制御部２０２により求められた信号強度が最大となる時刻・周波数のヒートマップを三面図および立体画像６４４に重畳させる。
【０１５３】
また、解析表示制御部２０２は、時刻・周波数全体を通じての信号強度の平均が最大となる脳内の位置を求めるものとしてもよい。この場合、ヒートマップ表示制御部２１１は、解析表示制御部２０２により求められた脳内の位置でのヒートマップ６１１を表示させる。また、立体表示制御部２１２は、表示されたヒートマップ６１１のうち信号強度が最大となる時刻・周波数に対応する脳の立体画像を立体図６１２に表示させる。また、断面表示制御部２１３は、解析表示制御部２０２により求められた脳内の位置を通るような三面図を頭部三面図６１３に表示させ、表示されたヒートマップ６１１のうち信号強度が最大となる時刻・周波数のヒートマップを三面図および立体画像６４４に重畳させる。
【０１５４】
また、解析表示制御部２０２は、脳全体で信号強度の平均値が最大となる時刻・周波数を求めるものとしてもよい。この場合、立体表示制御部２１２は、解析表示制御部２０２により求められた時刻・周波数に対応する脳の立体画像を立体図６１２に表示させる。また、ヒートマップ表示制御部２１１は、立体図６１２の立体画像に表示された、解析表示制御部２０２により求められた時刻・周波数に対応するヒートマップにおいて信号強度が最大となる脳内の位置を求め、当該位置でのヒートマップ６１１を表示させる。また、断面表示制御部２１３は、ヒートマップ表示制御部２１１により求められた脳内の位置を通るような三面図を頭部三面図６１３に表示させ、解析表示制御部２０２により求められた時刻・周波数のヒートマップを三面図および立体画像６４４に重畳させる。
【０１５５】
また、立体表示制御部２１２は、ピークリスト６１４に登録されているピーク情報のうち、先頭のピーク情報が示す脳内の位置でのヒートマップ６１１を表示させるものとしてもよい。また、立体表示制御部２１２は、ピークリスト６１４に登録されているピーク情報のうち、先頭のピーク情報が示す時刻・周波数に対応する脳の立体画像を立体図６１２に表示させる。また、断面表示制御部２１３は、ピークリスト６１４に登録されているピーク情報のうち、先頭のピーク情報が示す脳内の位置を通るような三面図を頭部三面図６１３に表示させ、当該ピーク情報が示す時刻・周波数のヒートマップを三面図および立体画像６４４に重畳させる。
【０１５６】
また、立体表示制御部２１２は、測定の対象（視覚野、聴覚野、体性感覚野、運動野および言語野等）に応じてプリセットされた脳内の位置でのヒートマップ６１１を表示させるものとしてもよい。また、立体表示制御部２１２は、測定の対象（視覚野、聴覚野、体性感覚野、運動野および言語野等）に応じてプリセットされた時刻・周波数に対応する脳の立体画像を立体図６１２に表示させる。また、断面表示制御部２１３は、測定の対象（視覚野、聴覚野、体性感覚野、運動野および言語野等）に応じてプリセットされた脳内の位置を通るような三面図を頭部三面図６１３に表示させ、当該ピーク情報が示す時刻・周波数のヒートマップを三面図および立体画像６４４に重畳させる。
【０１５７】
次に、解析者が、時間周波数解析画面６０１を起動した（開いた）場合に、立体図６１２の脳の立体画像、および、頭部三面図６１３の立体画像６４４の表示の初期視点について説明する。
【０１５８】
例えば、測定の対象（視覚野、聴覚野、体性感覚野、運動野および言語野等）に応じてプリセットされた視点を初期視点として用いるものとしてもよい。この場合、立体図６１２については、行（視点）の数もプリセットしておく。立体図６１２が２行である場合、２つの視点をプリセットしておく必要がある。例えば、言語野が測定対象であれば、脳の左側面および右側面の視点をプリセットしておく。
【０１５９】
また、ピークリスト６１４において先頭に登録されているピークが最もよく見えるような視点を初期視点として用いるものとしてもよい。具体的には、図６４に示すように、脳の中心とピークとを結んだ直線８１１上に初期視点として視点Ｐ０を設定するものとすればよい。
【０１６０】
また、ピークリスト６１４において所定のパラメータ（例えば、図５７に示すピークの値（信号強度）またはピークの高さ）が、所定の閾値を超えたピークを用いて設定した視点を初期視点として用いるものとしてもよい。例えば、閾値を超えたピークが２つある場合、立体図６１２の表示を２行表示とし、図６５に示すように、脳の中心とそれぞれのピークとを結んだ直線８１２、８１３上に、初期視点としてそれぞれ視点Ｐ１、Ｐ２を設定するものとすればよい。この場合において、視点Ｐ１からの脳の立体画像を立体図６１２の上の行に、視点Ｐ２からの脳の立体画像を立体図６１２の下の行に表示した例を、図６６に示す。
【０１６１】
以上のように、本実施形態では、生体信号の信号強度の分布を示すヒートマップ６１１上で特定の時刻および周波数が指定されることにより、指定された位置の信号強度に基づく閾値判定を行うことによって、所定の条件を満たす信号強度を有する時刻および周波数の領域を自動的に抽出するものとしている。これによって、注目する生体信号についての時刻および周波数の詳細な領域を容易に抽出することができる。また、このような詳細な領域を自動的に指定できることによって、より適切なヒートマップを立体図６１２および頭部三面図６１３に重畳することができる。
【０１６２】
また、脳内の特定の位置または特定の範囲における生体信号の時刻・周波数に関するヒートマップ６１１を表示するものとし、そのヒートマップ６１１上で指定された点または範囲に対応する脳の活動を示すヒートマップが重畳された立体画像を中心に、その時刻の前後の脳の活動を示す立体画像を表示、すなわち、脳の活動を示す静止画（ここでは立体画像）をコマ送り的またはコマ戻し的に表示するものとしている。これによって、脳の活動を示す静止画を適切かつ迅速に抽出することができ、脳の活動の解析を容易にすることができる。また、カンファレンス等で議論のベースにすることが容易となる。
【０１６３】
また、ピークリスト６１４に登録されたピーク情報から特定のピーク情報を選択することによって、当該ピーク情報に対応したヒートマップ６１１、立体図６１２および頭部三面図６１３が表示される。これによって、選択したピークが脳のどの位置、およびどの時刻・周波数であるかを瞬時に認識することができ、さらに、ヒートマップ６１１においては当該ピークおよびその周りの時刻／周波数における信号強度の状態、および、当該ピークの位置およびその周りの脳の信号強度の状態も把握することが可能となる。
【０１６４】
また、立体図６１２上で自由に脳の視点を変更することができ、さらに、その視点変更に基づく変更を同じ行または異なる行の脳に対しても反映することができる。これによって、特定の立体画像（対象立体画像）に対する視点変更を行うのみで、他の立体画像についてもその視点変更に基づく変更が反映されるので、操作性が向上し、さらに、複数の行を見比べることによって、対応する視点から時間的に前後の脳の活動変化を容易に確認することができる。また、立体画像で描画される脳の視点を自在に変えることによって、ある視点では見えなかった発火位置を確認することができる。
【０１６５】
また、頭部三面図６１３における立体画像６４４に対して行った視点変更を、各種設定に応じて、立体図６１２の時系列に並んだ脳の立体画像の視点に反映することができる。これによって、立体画像６４４の視点変更と同様の視点変更を、立体図６１２に対して改めて行う必要がないので、操作性が向上し、さらに、立体画像６４４で変更した視点と同一の視点、または対応する視点で、立体図６１２において脳の状態の変化を時系列に確認することができる。
【０１６６】
なお、上述の実施形態は、生体部位としての脳についての生体信号を取り扱うものとしたが、これに限定されるものではなく、例えば、脊髄、筋肉等の生体部位の生体信号に対して適用することが可能である。例えば、脳の図で説明した立体図６１２を、腰椎の場合、図６７のように表示することが可能である。図６７では、図６７（ａ）〜図６７（ｄ）の順で時系列に腰椎の信号が上に伝わっていく状態を示している。
【０１６７】
また、上述の実施形態において、生体信号計測システム１の各機能部の少なくともいずれかがプログラムの実行によって実現される場合、そのプログラムは、ＲＯＭ等に予め組み込まれて提供される。また、上述の実施形態に係る生体信号計測システム１で実行されるプログラムは、インストール可能な形式または実行可能な形式のファイルでＣＤ−ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ−Ｒ（Ｃｏｍｐａｃｔ  Ｄｉｓｋ  Ｒｅｃｏｒｄａｂｌｅ）、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｃ）等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよい。また、上述の実施形態の生体信号計測システム１で実行されるプログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、上述の実施形態の生体信号計測システム１で実行されるプログラムを、インターネット等のネットワーク経由で提供または配布するように構成してもよい。また、上述の実施形態の生体信号計測システム１で実行されるプログラムは、上述した各機能部のうち少なくともいずれかを含むモジュール構成となっており、実際のハードウェアとしてはＣＰＵがＲＯＭ等からプログラムを読み出して実行することにより、上述の各機能部が主記憶装置上にロードされて生成されるようになっている。
【符号の説明】
【０１６８】
１  生体信号計測システム
３  測定装置
４  測定テーブル
３１  デュワ
３２  窪み
４０  サーバ
５０  情報処理装置
１０１  ＣＰＵ
１０２  ＲＡＭ
１０３  ＲＯＭ
１０４  補助記憶装置
１０５  ネットワークＩ／Ｆ
１０６  入力装置
１０７  表示装置
１０８  バス
２０１  収集表示制御部
２０２  解析表示制御部
２０３  ピークリスト制御部
２０４  通信部
２０５  センサ情報取得部
２０６  解析部
２０７  記憶部
２０８  入力部
２１１  ヒートマップ表示制御部
２１２  立体表示制御部
２１３  断面表示制御部
２１４  再生表示制御部
２２１  取得部
２２２  二値化部
２２３  抽出部
５０１  開始画面
５０２  測定収集画面
５１１ａ、５１１ｂ  領域
５１２  モニタウィンドウ
５２１〜５２３  第２表示領域
５２３ａ−１、５２３ａ−２  マーク
５２６−１、５２６−２  属性アイコン
５３０  第１表示領域
５３０ａ−１、５３０ａ−２  アノテーション
５３１  時間軸
５３２  ライン
５３５  ポップアップウィンドウ
５３５ａ  選択ボタン
５３５ｂ  入力ボックス
５３８  カウンタボックス
５３９  終了ボタン
５４１、５４２  脳磁分布図
５５０  脳波分布図
５６０  アノテーションリスト
５６０ａ  選択ボックス
６０１  時間周波数解析画面
６０５  解析画面切替リスト
６１１  ヒートマップ
６１２、６１２ａ  立体図
６１２−１〜６１２−３  表示領域
６１３、６１３ａ  頭部三面図
６１４  ピークリスト
６１４ａ  ピークリスト設定ボタン
６１５  再生制御パネル
６２１、６２１−１  指定部
６２１−２〜６２１−５  対応部
６２１ａ  指定部
６２２、６２２−１  指定領域
６２２−２〜６２２−５  対応領域
６２２ａ  指定領域
６２３  指定領域
６２５〜６２７  領域
６３１−１〜６３１−５  線分
６３５  コメント
６４１〜６４３  断面図
６４４  立体画像
６４５ａ〜６４５ｄ  基準線
６４６、６４６ａ、６４６ｂ  ピーク点
６４７  ダイポール
６５０〜６５２  ダイアログボックス
６５５  ダイアログボックス
６５６  プルダウンリスト
７１２ａ−１〜７１２ａ−５  部位
７１２ｂ−１〜７１２ｂ−５  部位
７１３ａ−１、７１３ａ−２、７１３ｂ〜７１３ｄ  部位
７２２−１〜７２２−５  枠
８０１、８０２  ピーク部
８１１〜８１３  直線
Ａ１、Ａ２  アノテーション
Ｐ０〜Ｐ２  視点
【先行技術文献】
【特許文献】
【０１６９】
【特許文献１】特開２００６−１１０２３４号公報
【図１】
2020146206.tif 000004
【図２】
2020146206.tif 000005
【図３】
2020146206.tif 000006
【図４】
2020146206.tif 000007
【図５】
2020146206.tif 000008
【図６】
2020146206.tif 000009
【図７】
2020146206.tif 000010
【図８】
2020146206.tif 000011
【図９】
2020146206.tif 000012
【図１０】
2020146206.tif 000013
【図１１】
2020146206.tif 000014
【図１２】
2020146206.tif 000015
【図１３】
2020146206.tif 000016
【図１４】
2020146206.tif 000017
【図１５】
2020146206.tif 000018
【図１６】
2020146206.tif 000019
【図１７】
2020146206.tif 000020
【図１８】
2020146206.tif 000021
【図１９】
2020146206.tif 000022
【図２０】
2020146206.tif 000023
【図２１】
2020146206.tif 000024
【図２２】
2020146206.tif 000025
【図２３】
2020146206.tif 000026
【図２４】
2020146206.tif 000027
【図２５】
2020146206.tif 000028
【図２６】
2020146206.tif 000029
【図２７】
2020146206.tif 000030
【図２８】
2020146206.tif 000031
【図２９】
2020146206.tif 000032
【図３０】
2020146206.tif 000033
【図３１】
2020146206.tif 000034
【図３２】
2020146206.tif 000035
【図３３】
2020146206.tif 000036
【図３４】
2020146206.tif 000037
【図３５】
2020146206.tif 000038
【図３６】
2020146206.tif 000039
【図３７】
2020146206.tif 000040
【図３８】
2020146206.tif 000041
【図３９】
2020146206.tif 000042
【図４０】
2020146206.tif 000043
【図４１】
2020146206.tif 000044
【図４２】
2020146206.tif 000045
【図４３】
2020146206.tif 000046
【図４４】
2020146206.tif 000047
【図４５】
2020146206.tif 000048
【図４６】
2020146206.tif 000049
【図４７】
2020146206.tif 000050
【図４８】
2020146206.tif 000051
【図４９】
2020146206.tif 000052
【図５０】
2020146206.tif 000053
【図５１】
2020146206.tif 000054
【図５２】
2020146206.tif 000055
【図５３】
2020146206.tif 000056
【図５４】
2020146206.tif 000057
【図５５】
2020146206.tif 000058
【図５６】
2020146206.tif 000059
【図５７】
2020146206.tif 000060
【図５８】
2020146206.tif 000061
【図５９】
2020146206.tif 000062
【図６０】
2020146206.tif 000063
【図６１】
2020146206.tif 000064
【図６２】
2020146206.tif 000065
【図６３】
2020146206.tif 000066
【図６４】
2020146206.tif 000067
【図６５】
2020146206.tif 000068
【図６６】
2020146206.tif 000069
【図６７】
2020146206.tif 000070