JP6996203B2 - 情報処理装置、情報処理方法、プログラムおよび生体信号計測システム - Google Patents

Description

本発明は、情報処理装置、情報処理方法、プログラムおよび生体信号計測システムに関する。

従来、測定した被検者の生体信号に基づき推定した生体内の信号源を断層画像上に重畳表示するとともに、信号源に対応する生体信号の波形を画面上に並列表示する技術が知られている（例えば特許文献１参照）。この技術においては、オペレータは、表示されている生体信号の波形の所望の時間を指定すると、断層画像上の対応する信号源位置を特定する表示が行われる。

例えば、脳の神経活動を測定する脳磁計や脳波計では、測定した波形からてんかん特有の波形箇所（以下、特異点という。）を判別し、特異点から信号源の推定を行い、信号源を断層画像上に重畳表示させる。そして、この断層画像上の信号源の位置などに基づき、手術にて切除する位置（てんかんの原因となっている箇所）を特定することになる。

しかしながら、従来の技術では、症例の原因となる対象箇所を特定する確度が十分ではないという問題がある。

本発明は、上記に鑑みてなされたものであって、症例の原因となる対象箇所を特定する確度を向上させることが可能な情報処理装置、情報処理方法、プログラムおよび生体信号計測システムを提供することを目的とする。

上述した課題を解決し、目的を達成するために、本発明は、所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示可能な表示制御部を備え、該表示制御部は、前記複数スライスされた生体断層画像のうち、所定の条件に合致する信号源が重畳された生体断層画像を含めて、前記複数スライスされた生体断層画像を層順に並べて表示領域に表示するよう制御し、前記信号源が重畳されない生体断層画像の表示または非表示の指示に応じて、信号源が重畳されない生体断層画像の表示・非表示を切り替えることを特徴する、情報処理装置である。

本発明によれば、症例の原因となる対象箇所を特定する確度を向上させることができる。

図１は、実施形態の生体信号計測システムの概略図である。 図２は、情報処理装置に表示される開始画面の一例を示す図である。 図３は、測定収録画面の一例を示す図である。 図４は、測定収録画面の左側の領域の拡大図である。 図５は、測定収録画面の右側の領域の拡大図である。 図６は、アノテーション情報が入力された直後の画面を示す図である。 図７は、更新されたアノテーションリストの図である。 図８は、測定収録時の情報処理装置の動作を示すフローチャートである。 図９は、解析画面の一例を示す図である。 図１０は、解析画面の左側の領域の拡大図である。 図１１は、解析画面の右側の領域の拡大図である。 図１２は、解析画面で特定のアノテーションラインが選択された直後の画面を示す図である。 図１３は、図１２の画面の左側の領域の拡大図である。 図１４は、図１２の画面の右側の領域の拡大図である。 図１５は、解析時の情報処理装置の動作を示すフローチャートである。 図１６は、情報処理装置のハードウェア構成図である。 図１７は、情報処理装置の機能ブロック図である。 図１８は、表示制御部が有する機能の一例を示す図である。 図１９は、第１の実施形態において、マージボタンが押下されたときに表示される画面を示す図である。 図２０は、３次元方向に対応する３つのスライスの画像の関係を説明するための図である。 図２１は、「only dipole」ボタンが押された時の画面の一例を示す図である。 図２２は、第１の実施形態の変形例において、マージボタンが押下されたときに表示される画面を示す図である。 図２３は、第１の実施形態の変形例において、マージボタンが押下されたときに表示される画面を示す図である。 図２４は、第１の実施形態の変形例において、マージボタンが押下されたときに表示される画面を示す図である。 図２５は、第１の実施形態の変形例において、マージボタンが押下されたときに表示される画面を示す図である。 図２６は、第１の実施形態の変形例において、マージボタンが押下されたときに表示される画面を示す図である。 図２７は、第２の実施形態の情報処理装置の動作例（３回にわたる測定収録時の動作）を示すフローチャートである。 図２８は、第２の実施形態の情報処理装置の動作例を示すフローチャートである。 図２９は、第２の実施形態における解析画面の左側の領域を示す図である。 図３０は、第２の実施形態の変形例の解析画面の一例を示す図である。 図３１は、異なる範囲情報ごとの信号波形を区別するための表示方法を説明するための図である。 図３２は、複数の測定ファイルごとに解析ファイルを紐付けて管理することを表す模式図である。 図３３は、第２の実施形態の変形例における選択画面の一例を示す図である。

以下、添付図面を参照しながら、本発明に係る情報処理装置、情報処理方法、プログラムおよび生体信号計測システムの実施形態を詳細に説明する。

（第１の実施形態）
図１は、本実施形態の生体信号計測システム１の概略図である。生体信号計測システム１は、被検者の複数種類の生体信号、たとえば脳磁（ＭＥＧ：Magneto-encephalography）信号と脳波（ＥＥＧ：Electro-encephalography）信号を計測し、表示する。測定対象となる生体信号は、脳磁信号および脳波信号に限られるものではなく、例えば心臓の活動に応じて発生する電気信号（心電図として表現可能な電気信号）であってもよい。図１に示すように、生体信号計測システム１は、被検者の１以上の生体信号を測定する測定装置３と、測定装置３で測定された１以上の生体信号を記録するサーバ４０と、サーバ４０に記録された１以上の生体信号を解析する情報処置装置５０と、を含む。ここでは、サーバ４０と情報処理装置５０とは別々に記載されているが、例えばサーバ４０が有する機能の少なくとも一部が情報処理装置５０に組み込まれる形態であってもよい。

図１の例では、被検者（被測定者）は、頭に脳波測定用の電極（またはセンサ）を付けた状態で測定テーブル４に仰向けで横たわり、測定装置３のデュワ３０の窪み３１に頭部を入れる。デュワ３０は、液体ヘリウムを用いた極低温環境の保持容器であり、デュワ３０の窪み３１の内側には脳磁測定用の多数の磁気センサが配置されている。測定装置３は、電極からの脳波信号と、磁気センサからの脳磁信号を収集し、収集した脳波信号および脳磁信号を含むデータ（以下の説明では「測定データ」と称する場合がある）をサーバ４０に出力する。サーバ４０に収録された測定データは、情報処理装置５０に読み出されて表示され、解析される。一般的に、磁気センサを内蔵するデュワ３０と測定テーブル４は磁気シールドルーム内に配置されているが、図示の便宜上、磁気シールドルームを省略している。

情報処理装置５０は、複数の磁気センサからの脳磁信号の波形と、複数の電極からの脳波信号の波形を、同じ時間軸上に同期させて表示する。脳波信号は、神経細胞の電気的な活動（シナプス伝達の際にニューロンの樹状突起で起きるイオン電荷の流れ）を電極間の電圧値として表すものである。脳磁信号は、脳の電気活動により生じた微小な磁場変動を表す。脳磁場は高感度の超伝導量子干渉計（ＳＱＵＩＤ）センサで検知される。これらの脳波信号および脳磁信号は「生体信号」の一例である。

図２は、情報処理装置５０に表示される開始画面２０４の一例を示す。開始画面２０４には、「測定収録」と「解析」の選択ボタンが表示される。脳波及び／または脳磁測定の場合、データの測定収録とデータの解析は、別々の主体によって行われる場合が多い。例えば、測定技師（測定者）によって「測定収録」のボタンが選択されると、測定装置３で測定されたデータは順次サーバ４０に保存され、情報処理装置５０に読み出されて表示される。測定収録の終了後、医師によって「解析」ボタンが選択されると、収録された測定データが読み出され解析される。

＜測定収録時の動作＞
図３は、測定収録画面の一例を示す。画面のタブ１１１に「測定収録」画面であることが表示されている。測定収録画面は、測定された信号波形を表示する領域２０１Ａと、信号波形以外のモニタ情報を表示する領域２０１Ｂを有する。信号波形を表示する領域２０１Ａは、測定者からみて画面の左側に配置され、信号波形以外のモニタ情報を表示する領域２０１Ｂは、測定者からみて画面の右側に配置されている。リアルタイムで検出され表示される波形の動き（画面の左側から右側に向かって表示される）に合わせた測定者の視線の動きと、画面左側の領域２０１Ａから右側の領域２０１Ｂへマウスを移動させるときの動きに無駄が生じず、作業効率が向上する。

表示画面の領域２０１Ｂでは、測定中に被測定者の様子を確認するためのモニタウィンドウ１７０が表示される。測定中の被測定者のライブ映像を表示することで、後述するように、信号波形のチェック、判断の信頼性を高めることができる。図３では、ひとつのモニタディスプレイ（後述の表示装置２８）の表示画面に、測定収録画面の全体が表示される場合を示しているが、左側の領域２０１Ａと右側の領域２０１Ｂを、２台もしくはそれ以上のモニタディスプレイに分けて別々に表示してもよい。

図４は、図３の左側の領域２０１Ａの拡大図である。領域２０１Ａは、信号検出の時間情報を画面の水平方向（第１方向）に表示する第１の表示領域１１０と、信号検出に基づく複数の信号波形を画面の垂直方向（第２方向）に並列に表示する第２の表示領域１０１～１０３を有する。

第１の表示領域１１０で表示される時間情報は、図４の例では、時間軸１１２に沿って付された時間表示を含むタイムラインであるが、時間（数字）を表示せずに、帯状の軸だけでもよいし、軸を設けずに時間（数字）の表示だけであってもよい。また、画面上側の表示領域１１０の他に、表示領域１０３の下側に時間軸１１２を表示して、タイムラインを表示してもよい。

領域２０１Ａでは、同種の複数のセンサから取得される複数の信号波形、あるいは複数種類のセンサ群から取得される複数種類の信号の波形が、同じ時間軸１１２で同期して表示される。たとえば、表示領域１０１には被測定者の頭部右側から得られる複数の脳磁信号の波形が、表示領域１０２には被測定者の頭部左側から得られる複数の脳磁信号の波形が、それぞれ並列に表示されている。表示領域１０３には、複数の脳波信号の波形が並列に表示されている。これらの複数の脳波信号波形は、各電極間で測定された電圧信号である。これらの複数の信号波形の各々は、その信号が取得されたセンサの識別番号あるいはチャネル番号と対応付けられて表示されている。

測定が開始され各センサからの測定情報が収集されると、時間の経過とともに領域２０１Ａの各表示領域１０１～１０３の左端から右方向に向けて信号波形が表示される。ライン１１３は計測の時刻（現在）を示しており、画面の左から右に向けて移動する。領域２０１Ａの右端（時間軸１１２の右端）まで信号波形が表示されると、その後は画面の左端から右に向けて徐々に信号波形が消え、消えた位置に新しい信号波形が順次左から右方向に表示され且つライン１１３も左端から右に向けて移動していく。これとともに、水平方向の表示領域１１０でも測定の進行に対応して、時間の経過が時間軸１１２上に表示される。測定収録は、終了ボタン１１９が押されるまで継続される。

実施形態の特徴として、測定者（収録者）がデータ収録中に信号波形上で波形の乱れ、振幅の特異点等に気付いたときに、問題となる箇所または範囲を信号波形上でマークすることができる。マーキングの箇所または範囲は、マウスによるポインタ操作あるいはクリック操作で指定することができる。指定された箇所（または範囲）は、表示領域１０１～１０３の信号波形上に強調表示されるとともに、指定結果が対応する時刻位置又は時間範囲で、表示領域１１０の時間軸１１２に沿って表示される。時間軸１１２上への表示を含むマーキングの情報は、信号波形データとともに保存される。指定された箇所は或る時刻に対応し、指定された範囲は或る時刻を含む一定範囲に対応する。

図４の例では、時刻ｔ１で、表示領域１０３で１以上のチャネルを含む範囲が指定され、マーク１０３ａ－１で時刻ｔ１を含んだ時間がハイライト表示されている。マーク１０３ａ－１の表示と関連して表示領域１１０の対応する時刻位置に、指定結果を示すアノテーション１１０ａ－１が表示されている。時刻ｔ２で、表示領域１０３で別の波形位置またはその近傍がマークされ、その位置（時刻ｔ２）または近傍の領域（少なくとも時間範囲か複数の波形のいずれか一つが指示される）にマーク１０３ａ－２がハイライト表示されている。同時に、表示領域１１０の対応する時刻位置（時間範囲）に、アノテーション１１０ａ－２が表示される。なお、アノテーションとは、あるデータに対して関連する情報を注釈として付与することを指す。本実施例では、少なくとも指定された時間情報に基づいて注釈として表示されるものであって、少なくとも時間情報に基づく波形が表示される位置と結び付けて、注釈として表示するものである。また、複数のチャネルを表示する場合、対応するチャネル情報と結びつけ、注釈として表示しても良い。

時刻ｔ１で表示領域１１０に追加されたアノテーション１１０ａ－１は、一例として、アノテーション識別番号と、波形の属性を示す情報を含む。この例では、アノテーション番号「１」とともに、波形の属性を表わすアイコンと「strong spike」（ストロングスパイク）というテキスト情報が表示されている。

時刻ｔ２で、測定者が別の波形箇所またはその近傍領域を指定すると、指定された箇所でマーク１０３ａ－２がハイライト表示され、これとともに、表示領域１１０の対応する時刻位置に、アノテーション番号「２」が表示される。さらに、ハイライト表示された箇所に、属性選択のためのポップアップウィンドウ１１５が表示される。ポップアップウィンドウ１１５は、種々の属性を選択する選択ボタン１１５ａと、コメントや追加情報を入力する入力ボックス１１５ｂを有する。選択ボタン１１５ａには、波形の属性として「速波（fast activity）」、「眼球運動（eye motion）」、「体動（body motion）」、「スパイク（spike）」など、波形乱れの要因が示されている。測定者は、画面の領域２０１Ｂのモニタウィンドウ１７０で被測定者の様子を確認することができるので、波形の乱れの原因を示す属性を適切に選択することができる。たとえば、波形にスパイクが生じたときに、てんかんの症状を示すスパイクなのか、被測定者の体動（くしゃみ等）に起因するスパイクなのかを判断することができる。

時刻ｔ１でも同じ操作が行われており、図４では、ポップアップウィンドウ１１５で「スパイク」の選択ボタン１１５ａが選択され、入力ボックス１１５ｂに「strong spike」と入力されたことにより表示領域１１０にアノテーション１１０ａ―１が表示されている。このような表示態様により、同じ時間軸１１２上に多数の信号波形を同期して表示する際に、信号波形の注目箇所または範囲を視認により容易に特定することができ、かつ注目箇所の基本情報を容易に把握することができる。

アノテーション１１０ａ－１の一部または全部、たとえば、属性アイコンとテキストアノテーションの少なくとも一方を、表示領域１０３の信号波形上のマーク１０３ａ－１の近傍にも表示してもよい。信号波形上へのアノテーションの追加は、波形形状のチェックの妨げになる場合もあり得るので、表示領域１０１～１０３の信号波形上にアノテーションを表示させる場合は、表示・非表示を選択可能にしておくことが望ましい。

カウンタボックス１１８は、スパイクアノテーションの累積数を表示する。「スパイク」が選択される都度、カウンタボックス１１８のカウンタ値がインクリメントされ、収録開始から現在（ライン１１３）までのトータルのスパイク数が一目でわかるようになっている。

図５は、図３の右側の領域２０１Ｂの拡大図であり、図４と同じ時刻（ライン１１３の時点）での状態を示す。領域２０１Ｂのモニタウィンドウ１７０では、頭部を測定装置３に入れて測定テーブル４に横たわっている被測定者の状態のライブ映像が表示される。領域２０１Ｂでは、表示領域１０１、１０２、１０３の信号波形の各々に対応する分布図１４１、１４２、１３０と、アノテーションリスト１８０が表示される。アノテーションリスト１８０は、図４の信号波形上でマークされたアノテーションの一覧である。表示領域１０１～１０３で信号波形上の位置または範囲が指定されアノテーションが付される都度、対応する情報がアノテーションリスト１８０に順次追加される。測定収録画面におけるアノテーションリスト１８０への追加と表示は、たとえば降順（新しいデータを上に表示）で行われるが、この例に限定されない。アノテーションリスト１８０の表示を昇順にしてもかまわないが、表示領域１１０で時間軸１１２に沿って表示されるアノテーションとの対応関係がわかるように表示する。さらに、表示順序を変更させたり、項目ごとにソートさせることも可能である。

図５の例では、アノテーション番号「１」に対応する時刻情報と、付加されたアノテーション情報がリストされている。アノテーション情報として、「スパイク」を表わす属性アイコンと、「strong spike」というテキストが記録されている。また、マーク１０３ａ－１がハイライト表示された時点で、アノテーション番号「２」に対応する時刻情報がリストされている。ここでは、「アノテーション」は、アノテーション番号と、時刻情報と、アノテーション情報との組であると考えてもよいし、アノテーション情報のみであると考えてもよいし、アノテーション情報と、アノテーション番号または時刻情報との組であると考えてもよい。

また、アノテーションリスト１８０の近傍に、表示／非表示の選択ボックス１８０ａが配置されている。選択ボックス１８０ａで非表示が選択されると、表示領域１０１～１０３で、信号波形上のハイライトマーク以外のアノテーションが非表示にされるが、表示領域１１０の時間軸１１２に沿ったアノテーションの表示は維持される。これにより、信号波形の視認性を阻害せずにアノテーション情報を認識可能にする。

図６は、時刻ｔ２でポップアップウィンドウ１１５の「スパイク」が選択され、「normal spike」（ノーマルスパイク）というテキストが入力された直後の画面を示す。図４で例示したポップアップウィンドウ１１５で「ＯＫ」ボタンが選択されると、ポップアップウィンドウ１１５が閉じて、図６に示すように表示領域１１０の対応する時刻位置にアノテーション１１０ａ－２が表示される。アノテーション番号「２」と対応付けて、「スパイク」を表わす属性アイコンと、「normal spike」のテキスト情報が表示される。これと同時に、カウンタボックス１１８の値がインクリメントされる。また、ハイライト表示されたマーク１０３ａ－２の近傍に、属性アイコン１０６－２が表示される。この例では、マーク１０３ａ－１の近傍にも属性アイコン１０６－１が表示されているが、上述したように、属性アイコン１０６－１、１０６－２の表示、非表示は選択可能である。マーク１０３ａ－１と属性アイコン１０６－１を含むアノテーションＡ１、及びマーク１０３ａ－２と属性アイコン１０６－２を含むアノテーションＡ２も、アノテーション情報に含まれる。

図７は、アノテーションリスト１８０を示す。画面左側の領域２０１Ａでマーク１０３ａ－２に対応するアノテーションが付加されたことにより、アノテーションリスト１８０が更新される。アノテーション番号「２」に「normal spike」というメモが追加される。

以下同様に、測定中に領域２０１Ａで信号波形上の特定の箇所または範囲が指定される都度、指定箇所が強調表示されるとともに、表示領域１１０の時間軸１１２に沿ってアノテーション情報が表示される。領域２０１Ｂでは、アノテーションリスト１８０にアノテーション情報が順次追加される。

アノテーションリスト１８０及び信号波形表示の領域２０１Ａにおいて、アノテーション番号の表示は必須ではなく、用いなくてもよい。付加されたアノテーションを識別できる情報であれば任意の情報を識別情報として用いることができる。たとえば、属性アイコンと属性文字列（「strong spike」等）と、時間軸１１２の近傍に時刻と対応付けて表示してもよい。さらに、領域２０１Ａにファイル番号（図７の「Ｆｉｌｅ」の項目に表示される番号）を併記して表示してもよい。

終了ボタン１１９（図４に示す）が選択（押下）され測定が終了すると、表示領域１０１～１０３で指定されたハイライト箇所は信号波形に対応付けて保存される。表示領域１１０の対応する時刻位置に表示されたアノテーション情報も、アノテーション番号と時刻に対応付けて保存される。カウンタボックス１１８のカウンタ値、アノテーションリスト１８０の内容等の関連情報も保存される。これらの表示情報を保存することで、測定者と解析者が異なる場合でも、解析者は容易に問題となる箇所を認識し、解析することができる。

図８は、情報処理装置５０で行われる測定収録段階での情報表示処理のフローチャートである。図２に示す開始画面２０４で「測定収録」が選択されると（ステップＳ１１）、測定が開始され、複数の信号の波形が同じ時間軸に沿って同期して表示される（ステップＳ１２）。ここで「複数の信号波形」という場合は、同一種類の複数のセンサで検知された信号波形と、異なる種類のセンサの各々で検知された複数の信号波形の両方を含む。この例では、複数の生体信号の波形は、被測定者の頭部右側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の頭部左側に対応する磁気センサ群から得られる脳磁信号の波形と、被測定者の脳波測定用の電極から得られる脳波信号の波形と、から構成されるが、これに限られるものではない。なお、センサの選択については、左／右のセンサ群にとらわれず、頭頂部、前頭葉、側頭葉などのパートから、任意に選択可能である。例えば図５等に示す「MEG Window Control 1」で頭頂部のセンサが選択された場合は、「MEG Window Control 2」でそれ以外のセンサが選択された状態となる。

情報処理装置５０は、表示されている信号波形上で注目箇所または範囲が指定されたか否かを判断する（ステップＳ１３）。注目箇所または範囲の指定があると（ステップＳ１３：Ｙｅｓ）、信号波形の表示領域（表示領域１０１～１０３）に指定箇所を強調表示するとともに、時間軸領域（表示領域１１０）の対応する時刻位置に指定結果を表示する（ステップＳ１４）。指定結果には、指定がなされたこと自体を示す情報、または指定の識別情報が含まれる。時間軸領域への指定結果の表示と同時に、あるいは前後して、アノテーションの入力要求の有無を判断する（ステップＳ１５）。アノテーション入力の要求がある場合は（ステップＳ１５：Ｙｅｓ）、入力されたアノテーション情報を時間軸領域の対応する時刻位置に表示するとともに、アノテーションリストに追加する（ステップＳ１６）。その後、測定終了コマンドが入力されたか否かを判断する（ステップＳ１７）。注目位置（領域）の指定がない場合（ステップＳ１３：Ｎｏ）、及びアノテーションの入力要求がない場合（ステップＳ１５：Ｎｏ）は、ステップＳ１７へ飛んで測定終了の判断を行う。測定が終了するまで（ステップＳ１７：Ｙｅｓ）、ステップＳ１３～Ｓ１６を繰り返す。

この情報表示方法により、複数のセンサからの信号を収集する際に、信号情報の視認性の高い測定収録画面が提供される。

＜解析時の動作＞
図９は、解析時の情報処理装置５０の画面の一例を示す。解析画面は、図２の開始画面２０４で「解析」ボタンを選択することで表示させる。画面のタブ１１１に「解析」画面であることが表示されている。解析画面は、測定により得られた被検者の１以上の生体信号（この例では、被測定者の頭部右側に対応する磁気センサ群から得られる脳磁信号、被測定者の頭部左側に対応する磁気センサ群から得られる脳磁信号、被測定者の脳波測定用の電極から得られる脳波信号）の経時的変化を示す生体データと、測定時に生体データに対して入力されたアノテーションとを対応付けた画面である。本実施形態の情報処理装置５０は、この解析画面を表示部（後述の表示装置２８）に表示する制御を行う機能を有している。図９の例では、解析画面は、収録された３つの生体信号の経時的変化を示す波形（生体データに相当）をアノテーションとともに表示する領域２０２Ａと、解析情報を表示する領域２０２Ｂを有する。収録された信号波形とアノテーション情報を表示する領域２０２Ａは、測定者からみて画面の左側に配置され、解析表示する領域２０２Ｂは、測定者からみて右側に配置されている。これは、解析の際に、領域２０２Ａで信号波形をチェックまたは選択しながら、マウス等を操作して領域２０２Ｂで解析結果を確認または確定させる作業の効率が良いからである。

本例では、領域２０２Ａの第２の表示領域１０３の脳波信号の波形の画面上方に、第２の表示領域１０１，１０２の脳磁信号の波形を表示している。また、領域２０２Ａの右側の領域２０２Ｂでは、領域２０２Ａに近い側の画面領域で且つ画面上方に脳磁分布図１４１，１４２を表示し、脳波分布図１３０をその下方に表示している。このため、解析者は、第２の表示領域１０３の「脳波信号の波形」→第２の表示領域１０１，１０２の「脳磁信号の波形」→脳磁分布図１４１、１４２→脳波分布図１３０の順に（この場合は、時計回りで）視線移動を行える。そのため、解析者（または測定者）の視線移動が効率的になり、その結果、解析作業効率を向上させることが可能となる。なお、上記では、時計回りとして説明したがこの例に限定されない。

図１０は、図９の左側の領域２０２Ａの拡大図である。領域２０２Ａは、測定時の時間情報を画面の水平方向（第１方向）に表示する表示領域１１０及び表示領域１２０と、収録された信号波形を種類ごとに画面の垂直方向（第２方向）に並べて表示する表示領域１０１～１０３を有する。

表示領域１１０には、収録時の時間の経過を示す時間軸１１２と、時間軸１１２に沿って付加されたアノテーション１１０ａ－７、１１０ａ－８が表示される。表示領域１２０には収録時間の全体を示す時間軸１２２が表示される。時間軸１２２に沿って、アノテーションが付加された時刻位置を示すポインタマーク１２０ａと、表示領域１０１～１０３に現在表示されている信号波形が収録された時間帯を示すタイムゾーン１２０ｂが表示される。この表示により、解析者は、現在解析中の信号波形が、測定収録時のどの段階で取得された信号波形なのかを直感的に把握することができる。

解析者は、解析画面を開いた後に、たとえば、時間軸１２２のバー上でタイムゾーン１２０ｂをドラッグすることで所望の時間帯の信号波形を表示領域１０１～１０３に表示させることができる。あるいは、後述するように、アノテーションリスト１８０の中から所望のアノテーションを選択することで、そのアノテーションを含む前後の信号波形を表示領域１０１～１０３に表示させることができる。

表示領域１０１～１０３には、収録時に信号波形に付加されたアノテーションＡ７、Ａ８が表示されている。マーク１０３ａ－７、１０３ａ－８がハイライト表示され、マーク１０３ａ－７、１０３ａ－８の近傍に対応する属性アイコン１０６－７、１０６－８が表示されている。また、マーク１０３ａ－７、１０３ａ－８の時刻位置を示す縦のライン１１７－７、１１７－８が表示されている。ライン１１７－７、１１７－８が表示されることで、たとえば、表示領域１０３の所定の箇所の指定と関連してアノテーションが付加されたときに、異なる種類の信号表示エリアである表示領域１０２、１０１においても、指定の結果が容易に視認できる。ライン１１７－７、１１７－８は、アノテーション情報の視認を容易にするという意味でアノテーション情報に含めることができ、「アノテーションライン」と称してもよい。ライン１１７－７、１１７－８を選択することで、その時刻の前後の一定時間を含めて、信号波形が拡大表示される。この処理については後述する。

図１１は、図１０と同じ時間での右側の領域２０２Ｂの拡大図である。表示領域１０１と１０２に表示されている信号波形に対応する脳磁分布図１４１、１４２、及び表示領域１０３に表示されている信号波形に対応する脳波分布図１３０が表示されている。また、脳磁図（ＭＥＧ：Magneto-encephalograph）の等磁場図１５０と、脳波図（ＥＥＧ：Electro-encephalograph）のマップエリア１６０と、ＭＲＩ（磁気共鳴画像：Magnetic Resonance Imaging）で取得された被測定者の脳の断層画像の表示ウィンドウ１９０が表示されている。等磁場図１５０では、磁場の湧き出し領域と沈み込み領域が色分けされて表示され、電流の流れる方向が視覚的に把握される。等磁場図１５０とマップエリア１６０は測定完了後に得られる情報であり、ＭＲＩの断層画像は、別途検査で得られる情報である。

モニタウィンドウ１７０には、表示領域１０１～１０３の信号波形が取得された時刻に同期して、測定時の被測定者の映像が表示される。解析者は、モニタウィンドウ１７０を見て被測定者の状態を確認しながら信号波形を解析することができる。

アノテーションリスト１８０は、測定収録で付加されたすべてのアノテーションがリストされている。アノテーションリスト１８０には、アノテーション番号１８１と対応付けて付加されたアノテーション情報（属性アイコン、テキスト入力情報など）が記載されている。解析画面のアノテーションリスト１８０は、たとえば、付加されたアノテーションが昇順で（古いデータが上になるように）表示されるが、これに限定されない。測定収録画面と同様にアノテーション番号の使用は必須ではなく、時刻、ファイル名、属性等の組み合わせでアノテーションを識別することも可能である。また、アノテーションリスト１８０に含まれるアノテーションの表示順序の変更や、項目ごとのソートも可能である。所望のアノテーション番号１８１または行をクリックすることで、図１０の表示領域１０１～１０３に、そのアノテーションが付加された時刻位置を含む所定の時間帯の信号波形を表示させることができる。

測定収録画面と異なり、解析者がアノテーション部分の信号波形を確認して最終的に信号源の推定がなされたアノテーションには、推定完了マーク１８２（図１１に示す）が表示されている。

アノテーションの表示／非表示を選択する選択ボックス１８０ａで非表示の指定がされると、図１０の表示領域１０３の属性アイコン１０６－７、１０６－８が消える。表示／非表示の選択ボックス１８０ａで、ハイライトされたマーク１０３ａ－７、１０３ａ－８の非表示を選択できるようにしてもよい。

図１２は、図１０の解析画面でライン１１７－７が選択（たとえばダブルクリック）された直後の画面の全体図である。解析者がアノテーションＡ７に着目して、この領域の波形を解析するためにライン１１７－７を選択（たとえばダブルクリック）すると、ハイライトされた信号波形の近傍の信号波形が拡大表示領域２００に拡大表示される。領域１１４で示される一定の時間範囲にわたって、信号波形が時刻位置を示すライン２１７－７とともに拡大表示される。

図１３は、図１２の左側の領域２０３Ａ（信号波形の表示領域）の拡大図である。解析者は、拡大表示領域２００に信号波形を拡大表示することで、収録時に付加されたマークの妥当性を再確認し、あるいは計測収録時にチェックされていない波形部分をチェックすることができる。たとえば、ライン２１７－７を左右にドラッグすることで、問題となる波形の正確な個所を特定または変更することができる。拡大表示領域２００に、表示領域１０３でハイライト表示されているマーク１０３ａ及びまたは属性アイコン１０６を反映させてもよい。ただし、振幅の特異点を正確に判断する際の視認の妨げになることも考えられるので、拡大表示領域２００にハイライトされたマーク１０３ａや属性アイコン１０６を表示する場合は、表示と非表示の選択を可能にしておくのが望ましい。

拡大表示領域２００に表示する信号波形の種類と、チャネル範囲を指定することも可能である。たとえば、解析者は、視線を表示領域１０３でハイライトされたマーク１０３ａ－７から画面の上方へ移し、脳磁波形の表示領域１０１または１０２の波形に振幅の特異点がないかを確認する。この場合、ボックス１２５に表示領域１０１または１０２のターゲットのチャネル領域を入力することで、拡大表示領域２００にマーク１０３ａ－７と関連する脳磁波形を拡大表示することができる。

拡大表示領域２００の画面の下側には、確認ウィンドウ２１０が表示されている。確認ウィンドウ２１０は、信号波形の属性ボタン２１１と、信号源の推定ボタン２１２を含む。属性ボタン２１１は、測定収録画面のポップアップウィンドウ１１５に含まれる属性情報と対応し、収録時に付加された属性が誤っているときは、属性ボタン２１１を選択して正しい属性を選択することができる。信号波形の正しい位置、及び／あるいは属性の選択が確認されたならば、推定ボタン２１２をクリックすることで信号源の推定をアノテーションに反映することができる。つまり、本実施形態の情報処理装置５０は、解析画面から選択されたアノテーションに対応する信号源を推定する機能を有している。後述するように、推定された信号源は、複数のＭＲＩの断層画像のうち該推定された信号源に対応する断層画像上に重畳して表示されることになる。

図１４は、図１２の右側の領域２０３Ｂの拡大図である。図１３で、所望のアノテーションについて、信号波形位置及び／または属性が確認され、信号源の推定ボタン２１２が選択されると、アノテーションリスト１８０の対応するアノテーション（この例では、アノテーション番号「７」）に、推定完了マーク１８２が付加される。さらに、表示ウィンドウ１９０のＭＲＩ断層画像に、ダイポールの推定結果１９０ａが表示される。

解析者によって、表示領域１０１～１０３にハイライト表示されるマーク位置、及び／またはアノテーションの内容が変更される場合に、アノテーションリスト１８０の更新方法には、２通りの方法がある。解析者による最新の更新情報だけをアノテーションリスト１８０に反映する方法と、測定収録時のアノテーション情報を維持したまま、新たなアノテーション情報として追加する方法である。後者の方法を採用する場合は、アノテーション識別情報として、たとえば収録時のアノテーション番号からの枝番号を付けることができる。この場合、表示領域１１０にも新たなアノテーション情報を追加して、追加のアノテーション情報を時間軸に沿って異なる色で表示してもよい。

図１５は、情報処理装置５０で行われる解析段階での情報表示処理のフローチャートである。開始画面２０４（図２参照）で「解析」が選択されると（ステップＳ２１）、解析が開始され、解析画面が表示される（ステップＳ２２）。初期の解析画面は、信号波形が表示されていないブランク画面でもよいし、収録の先頭または最後の一定の時間範囲の信号波形であってもよい。解析画面が表示されると、特定のアノテーションが選択されたか否かが判断される（ステップＳ２３）。アノテーションの選択は、アノテーションリスト１８０内での特定のアノテーション番号または行の選択であってもよいし、表示領域１２０の時間軸１２２上のタイムゾーン１２０ｂを操作することによる時刻位置の指定であってもよい。アノテーションの選択があると（ステップＳ２３：Ｙｅｓ）、選択されたアノテーションの時刻位置を含む所定時間分の信号波長が表示される（ステップＳ２４）。

表示された場面で、ハイライト表示されたマークの時間位置を示すライン１１７が選択されたか否かが判断される（ステップＳ２５）。ライン１１７が選択されると（ステップＳ２５：Ｙｅｓ）、選択されたラインを含む一定の時間範囲の信号波形が拡大表示される（ステップＳ２６）。拡大表示は、必ずしもハイライト表示されたマークの近傍の信号波形に限定されず、同じ時間位置の異なる種類の信号波形を拡大表示してもよい。たとえば、脳波信号波形にハイライト表示されたマークが付されている場合に、同じ時間位置の脳磁信号波形を拡大表示させてもよい。また、すべてのチャネルの信号波形を拡大表示する代わりに、マークされた信号波形が取得されたチャネルを含む一定範囲のチャネルで取得された信号波形を拡大表示させてもよい。この場合、拡大表示させたい信号波形の種類、及び／またはチャネル範囲の指定入力の有無を判断してもよい。

次に、信号源の推定ボタン２１２が押下されたか否かが判断される（ステップＳ２７）。信号源の推定ボタン２１２が押下されると（ステップＳ２７：Ｙｅｓ）、信号源推定の演算が行われる。推定結果がＭＲＩ断層画面に表示されるとともに、推定完了マーク１８２がアノテーションリスト１８０に追加される（ステップＳ２８）。そして、アノテーションリスト１８０の下に配置されたマージボタン３００の押下を受け付けた場合（ステップＳ２９：Ｙｅｓ）、情報処理装置５０は後述の画面４００を表示し、画面４００に関する処理を行う（ステップＳ３０）。ステップＳ２９およびステップＳ３０の具体的な内容については後述する。マージボタン３００の押下を受け付けなかった場合（ステップＳ２９：Ｎｏ）、または、ステップＳ３０の後は、解析終了コマンドが入力されたか否かが判断される（ステップＳ３１）。アノテーションの選択がない場合（ステップＳ２３：Ｎｏ）、拡大表示のためのアノテーションラインの選択がない場合（ステップＳ２５：Ｎｏ）、及び信号源の推定ボタン２１２押下の入力がない場合（ステップＳ２７：Ｎｏ）は、ステップＳ３１へ飛んで解析終了の判断を行う。解析終了コマンドが入力されるまで（ステップＳ３１：Ｙｅｓ）、ステップＳ２３～Ｓ３０が繰り返される。

ステップＳ２６とＳ２７の間に、アノテーションが変更されたか否かの判断を行ってもよい。アノテーションが変更された場合は、アノテーションリスト１８０へ変更を反映して、ステップＳ２７の判断に移行する。

上述した表示処理動作により、視認性と操作性に優れた情報表示が実現する。

図１６は、情報処理装置５０のハードウェア構成図である。情報処理装置５０は、ＣＰＵ（Central Processing Unit：プロセッサ）２１、ＲＡＭ（Random Access Memory）２２、ＲＯＭ（Read Only Memory）２３、補助記憶装置２４、入出力インタフェース２５、及び表示装置２８を有し、これらがバス２７で相互に接続されている。

ＣＰＵ２１は、情報処理装置５０の全体の動作を制御し、各種の情報処理を行う。ＣＰＵ２１はまた、ＲＯＭ２３または補助記憶装置２４に格納された情報表示プログラムを実行して、測定収録画面と解析画面の表示動作を制御する。ＲＡＭ２２は、ＣＰＵ２１のワークエリアとして用いられ、主要な制御パラメータや情報を記憶する不揮発ＲＡＭを含んでもよい。ＲＯＭ２３は、基本入出力プログラム等を記憶する。本発明の情報表示プログラムもＲＯＭ２３に保存されてもよい。補助記憶装置２４は、ＳＳＤ（Solid State Drive）、ＨＤＤ（Hard Disk Drive）などの記憶装置であり、たとえば、情報処理装置５０の動作を制御する制御プログラムや、情報処理装置５０の動作に必要な各種のデータ、ファイル等を格納する。入出力インタフェース２５は、タッチパネル、キーボード、表示画面、操作ボタン等のユーザインタフェースと、各種センサあるいはサーバ４０からの情報を取り込み、他の電子機器に解析情報を出力する通信インタフェースの双方を含む。表示装置２８は各種の情報を表示するモニタディスプレイである。表示装置２８では測定収録画面と解析画面が表示され、入出力インタフェース２５を介した入出力操作に応じて画面が更新される。

図１７は、情報処理装置５０の機能ブロック図である。情報処理装置５０は、制御部２５０、解析部２５２、センサ情報取得部２５３、収録／解析情報保存部２５４、及びアノテーション入力部２５５を有する。制御部２５０は、情報処理装置５０の画面表示を制御する表示制御部２５１を含む。

センサ情報取得部２５３は、測定装置３あるいはサーバ４０から、センサ情報を取得する。アノテーション入力部２５５は、センサ情報に付加されるアノテーション情報を入力する。解析部２５２は、収集されたセンサ情報を解析する。センサ情報の解析には、信号波形の解析、振幅の特異点の解析、電流ダイポールの向きを含む脳磁場の解析が含まれる。つまり、この例では、解析部２５２は、解析画面から選択されたアノテーションに対応する信号源を推定する機能（推定部の機能）を有している。表示制御部２５１は、センサ情報の測定収録時、及び解析時の画面表示を図２～図１７を参照して説明した手法で制御する。収録／解析情報保存部２５４は、測定データと解析結果を保存する。測定収録時に信号波形にアノテーションが付加された場合は、信号波形が取得された時間情報と対応付けて、アノテーションも保存される。表示制御部２５１を含む制御部２５０の機能は、図１６のＣＰＵ２１がＲＯＭ２３等に格納されたプログラムをＲＡＭ２２上に展開して実行することによって実現される。解析部２５２の機能も、ＣＰＵ２１がＲＯＭ２３等に格納されたプログラムをＲＡＭ２２上に展開して実行することにより実現される。なお、これに限らず、例えば制御部２５０および解析部２５２の機能のうちの少なくとも一部が専用のハードウェア回路（半導体集積回路等）で実現される形態であってもよい。センサ情報取得部２５３とアノテーション入力部２５５の機能は入出力インタフェース２５によって実現される。収録／解析情報保存部２５４の機能は、ＲＯＭ２３または補助記憶装置２４によって実現される。

図１８は、表示制御部２５１が有する機能の一例を示す図である。図１８に示すように、表示制御部２５１は、第１の表示制御部２６１と、第２の表示制御部２６２と、を有する。図１８では、本発明に関する機能のみを例示しているが、表示制御部２５１が有する機能はこれらに限られるものではなく、それ以外の上述した機能も当然に有している。

第１の表示制御部２６１は、上述の解析画面を表示装置２８に表示する制御を行う。

前述したように、解析部２５２の解析により信号源の推定が順次に行われていき、アノテーションリスト１８０の下に配置されたマージボタン３００の押下を受け付けた場合、第２の表示制御部２６２は、所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示するとともに、該複数スライスされた生体断層画像のうち、所定の信号源が重畳された生体断層画像を表示領域に初期表示するよう制御する。ここでは、所定の信号源とは、所定の条件に合致する信号源である。この例では所定の条件は信号源数であるが、これに限られるものではない。本実施形態では、所定の信号源とは、複数スライスされた生体断層画像毎の信号源の数の中で最も多い信号源である。具体的な内容は後述するが、第２の表示制御部２６２は、信号源の数が最も多い生体断層画像を、表示領域の中央部付近に初期表示する。そして、第２の表示制御部２６２は、中央部付近の生体断層画像から左右に生体断層画像が層順に並ぶように、他の生体断層画像を配置して表示する。また、第２の表示制御部２６２は、信号源が重畳されない生体断層画像も表示することができるし、信号源が重畳されない生体断層画像は表示しないこともできる。以下、具体的な内容を説明する。

本実施形態では、第２の表示制御部２６２は、図１４に示すアノテーションリスト１８０の下に配置されたマージボタン３００の押下を受け付けると、図１９に示すような画面４００を表示装置２８に表示する制御を行う。画面４００は、複数の生体断層画像が左右方向に並んで表示される領域４０１Ａと、領域４０１Ａから選択された生体断層画像の断層位置が表示される領域４０１Ｂと、を有する。

領域４０１Ａは、上面からの断層画像（以下の説明では「スライス画像Ａ」と称する場合がある）を表示する表示領域４１０Ａと、側面方向からの断層画像（以下の説明では「スライス画像Ｂ」と称する場合がある）を表示する表示領域４１０Ｂと、背面方向からの断層画像（以下の説明では「スライス画像Ｃ」と称する場合がある）を表示する表示領域４１０Ｃとで構成されている。以下の説明では、スライス画像Ａ、Ｂ、Ｃを互いに区別しない場合は、単に「スライス画像」と称する場合がある。なお、領域４０１Ａにおける断層画像の上下並び方向は、本実施形態の態様に限られない。

ここで、図２０を参照して、３次元方向に対応するスライス画像Ａ、スライス画像Ｂ、スライス画像Ｃの関係を説明する。図２０において、スライス画像Ａ～Ｃは、３次元方向に位置がリンクしている。基準線１９０Ｅは、各スライス画像をまたがるように表示され、各基準線１９０Ｅの交点Ｏは各スライス画像のスライス位置を表示している。この例では、スライス画像Ｃは、水平方向（左右方向）の基準線１９０Ｅでスライス画像Ａを切断した断面を図２０に示すＡ方向から見た断面図である。また、スライス画像Ｂは、上下方向の基準線１９０Ｅでスライス画像Ａを切断した断面を図２０に示すＢ方向から見た断面図である。以下の説明では、スライス画像Ａに対応する視点を「Axial View」、スライス画像Ｂに対応する視点を「Sagittal View」、スライス画像Ｃに対応する視点を「Coronal View」と称する。

つまり、この例では、生体断層画像は、第１方向の断面である第１断層画像（例えばスライス画像Ａ）と、第１方向と直交する第２方向の断面である第２断層画像（例えばスライス画像Ｂ）と、第１方向および第２方向と直交する第３方向の断面である第３断層画像（例えばスライス画像Ｃ）と、を含む。

図１９に戻って説明を続ける。第２の表示制御部２６２は、重畳される信号源の数を示す情報を、対応する生体断層画像とともに表示する制御を行う。各表示領域４１０Ａ～４１０Ｃには、スライス画像ごとに、どの位置でスライスした画像であるかを示すスライス番号を示す情報４４０Ａと、該スライス画像に重畳される信号源の数（ダイポールの推定結果の数）を示す情報４４０Ｂとが表示される。

領域４０１Ｂは、表示領域４１０Ａに対応する表示領域４２０Ａと、表示領域４１０Ｂに対応する表示領域４２０Ｂと、表示領域４１０Ｃに対応する表示領域４２０Ｃとを有する。

表示領域４２０Ａには、表示領域４１０Ａに表示されるスライス画像Ａが、側面（表示領域４２０Ａの左側の画像）および背面（表示領域４２０Ａの右側の画像）のそれぞれから見た断層画像においてどの位置でスライスした画像であるかが表示され、この断層位置を示す断層位置ライン４５０が重畳して表示される。側面から見たスライス位置Ａと背面から見たスライス位置Ｂの各隣り合う断層位置ライン４５０は、図における上下方向の位置が一致している。また、各断層位置ライン４５０に対応するスライス番号が表示領域４１０Ａに表示されるスライス画像Ａと紐付けられている。たとえば、表示領域４２０Ａの下から上に向けてスライス番号が０１～１５と付与される。

同様に、表示領域４２０Ｂには、表示領域４１０Ｂに表示されるスライス画像Ｂが、上面（表示領域４２０Ｂの左側の画像）および背面（表示領域４２０Ｂの右側の画像）から見た断層位置においてどの位置でスライスした画像であるかが表示され、この断層位置を示す断層位置ライン４５０が重畳して表示される。上面から見たスライス位置Ｃと背面から見たスライス位置Ｄの断層位置ライン４５０は、図における左右方向の位置が一致している。また、各断層位置ライン４５０に対応するスライス番号が表示領域４１０Ｂに表示されるスライス画像Ｂと紐付けられている。たとえば、表示領域４２０Ｂの左から右に向けてスライス番号が０１～１４と付与される。

同様に、表示領域４２０Ｃには、表示領域４１０Ｃに表示されるスライス画像Ｃが、上面（表示領域４２０Ｃの左側の画像）および側面（表示領域４２０Ｃの右側の画像）から見た断層位置においてどの位置でスライスした画像であるかが表示され、この断層位置を示す断層位置ライン４５０が重畳して表示される。上面から見たスライス位置Ｅの断層位置ライン４５０の上から下に向けて、側面から見たスライス位置Ｆの断層位置ライン４５０の左から右に向けて、位置が一致している。また、各断層位置ライン４５０に対応するスライス番号が表示領域４１０Ｃに表示されるスライス画像Ａと紐付けられている。たとえば、表示領域４２０Ｃの上から下（左の上面断層画像の場合）または左から右（右側の側面断層画像の場合）に向けてスライス番号が０１～１５と付与されている。

つまり、本実施形態の第２の表示制御部２６２は、領域４０１Ａ（表示領域）に表示される生体断層画像の断層位置を示す情報を表示する制御を行う。この例では、第２の表示制御部２６２は、複数の生体断層画像（スライス画像）の中から選択された生体断層画像の断層位置を示す情報を表示する制御を行う。なお、上述の各断層位置ライン４５０と、上述のスライス番号を示す情報４４０Ａとは連携して記憶装置（補助記憶装置２４等）に格納されている。

この例では、各表示領域４１０Ａ～４１０Ｃで表示されるスライス画像は、重畳表示されるダイポールの推定結果の数が最も多いスライス画像が中央に配置されている。そして、中央のスライス画像から左右にスライス番号順に（層順に）並ぶように、他のスライス画像が配置される。例えば、表示領域４１０Ａには、スライス番号１０のスライス画像が中央に配置され、その右側にスライス番号１１，１２，１３（一部のみ）の順で配置されている。さらに、スライス番号１０の左側にはスライス番号９，８，７（一部のみ）の順で配置されている。また、表示領域４１０Ｂには、スライス番号１０のスライス画像が中央に配置され、その右側にスライス番号１１，１２，１３（一部のみ）の順で配置されている。さらに、スライス番号１０の左側にはスライス番号９，８，７（一部のみ）の順で配置されている。そして、表示領域４１０Ｃには、スライス番号７のスライス画像が中央に配置され、その右側にスライス番号８，９，１０（一部のみ）の順で配置されている。さらに、スライス番号７の左側にはスライス番号６，５，４（一部のみ）の順で配置されている。ここで、中央とは、領域４０１Ａ（「表示領域」に対応）の幅方向における中央である。また、解析者が視覚的に見つけやすくするために、例えば図１９のように、ダイポールの推定結果の数が最も多いスライス画像の上にタイトル（Axial View、Sagittal View、Coronal View）を表示させてもよい。この場合、解析者がスライス画像の初期表示状態から左右にスクロールして他のスライス画像を表示させたとき、これらタイトルが連動して移動してもよい。連動して移動することで、スクロール後でも、ダイポールの推定結果の数が最も多いスライス画像を見つけやすくできる。一方、どの方向の断層画像であるかに注目するのであれば、スクロールしてもタイトルは固定とすることが好ましい。

なお、本実施形態においては、表示領域４１０Ａに表示されたスライス画像Ａと、その直下に表示されたスライス画像Ｂ、および、さらにその直下に表示されたスライス画像Ｃとは３次元方向に対応していない。つまり、表示領域４１０Ａ～４１０Ｃごとに、該表示領域４１０に表示される複数のスライス画像のうち、重畳表示されるダイポールの推定結果１９０ａの数が最も多いスライス画像が中央に配置され、中央のスライス画像から左右にスライス番号順に並ぶように、他のスライス画像が配置される。このように表示することで、中央から左右へのダイポールの推定結果１９０ａの広がりを視認できる。

また、各スライス画像にはダイポールの推定結果の数を示す情報４４０Ｂも表示されるので、どのスライス画像にどれだけダイポールの推定結果が重畳されているかを見たいときに、併せて確認することができる。また、領域４０１Ｂの断層位置ライン４５０と、選択した（注目したい）スライス画像のスライス番号を示す情報４４０Ａから、所定範囲（たとえば１ｍｍ以内）にダイポールの推定結果が収まっているかどうかを把握できる。また、領域４０１に全てのスライス画像を表示できない場合は、たとえばマウスでスクロールすることで、スライス画像を左右方向に移動させて新たなスライス画像を表示することもできる。つまり、第２の表示制御部２６２は、スライス画像を送る／戻すための操作（スクロールするための操作）に応じて、スライス画像を左右方向に移動して、新たなスライス画像を表示する制御を行うことができる。

また、この例では、表示領域４１０Ａの上には、ダイポールの推定結果が重畳表示されているスライス画像のみを表示するモードを選択するための「Only Dipole」ボタン４３０Ｂと、ダイポールの推定結果が重畳されないスライス画像も含めた全てのスライス画像を表示するモードを選択するための「ALL」ボタン４３０Ａと、が配置されている。図１９は、「ALL」ボタン４３０Ａが押されたときの画面４００の一例を示す図である。図２１は、「only dipole」ボタン４３０Ｂが押された時の画面４００の一例を示す図であり、図２１に示すように、ダイポールの推定結果が重畳表示されているスライス画像のみが領域４０１Ａに表示され、ダイポールの推定結果が重畳されないスライス画像に対応する断層位置ライン４５０が領域４０１Ｂにて非表示になる。言い換えると、ダイポールの推定結果が重畳されたスライス画像に対応する断層位置ライン４５０のみが表示される。このように、注目したいスライス画像と、対応する断層位置ライン４５０とを同一画面上で比較することで、各ダイポールの推定結果がどのくらい離れた位置にそれぞれあるかが容易に認識できる。

解析者は、ダイポールの推定結果が重畳表示されているスライス画像から、どの位置に一番多くダイポールの推定結果が存在するかを検証することができる。そして、出力ボタン４６０を押下すると、ダイポールの推定結果が重畳表示されたスライス画像が出力（そのときの画面４００が出力）され、プリントアウトされる。このように、立体的な信号源（ダイポールの推定結果）の位置を従来に比べてより詳細に特定することができる。

以上に説明したように、本実施形態では、領域４０１Ａに初期表示するスライス画像の条件として、信号源の数が最も多いスライス画像であることを条件とし、該信号源の数が最も多いスライス画像を少なくとも初期表示する。ここでは、上述したように、第２の表示制御部２６２は、信号源の数が最も多いスライス画像を領域４０１Ａの中央部付近に初期表示し、中央部付近のスライス画像から左右にスライス画像が層順に並ぶように、他のスライス画像を配置して表示する。このように表示することで、中央から左右への信号源の広がりを視認できるので、解析者が、症例の原因となる対象箇所を特定する確度を向上させることができる。また、隣り合うスライス画像同士で信号源の有無を確認することができる。

なお、例えば、まず全てのスライス画像に対して信号源（群）を重畳して信号源重畳済みスライス画像を生成し、信号源数が最も多い信号源重畳済みスライス画像を選定してもよい。このような、信号源数が最も多い信号源重畳済みスライス画像を選定する機能は、第２の表示制御部２６２が有していてもよいが、これに限らず、例えば上記機能が、第２の表示制御部２６２とは別に設けられる形態であってもよい。つまり、所定の条件に合致する生体断層画像（領域４０１Ａに初期表示する生体断層画像）を選定する機能（選定部）が、第２の表示制御部２６２とは別に設けられる形態であってもよい。なお、上記機能（選定部）は、ソフトウェア的に実現（例えばＣＰＵ２１がプログラムを実行することにより実現）されてもよいし、専用のハードウェア回路で実現されてもよい。

また、例えば、まずすべてのスライス画像ごとに、該スライス画像に存在する信号源（群）を特定し、その結果から、信号源数が最も多いスライス画像を選定し、その後、その選定したスライス画像に信号源を重畳させて表示してもよい。なお、例えば信号源を重畳させずに、以上のようにして選定したスライス画像（信号源数が最も多いスライス画像）を初期表示しておき、任意のタイミングで、信号源（群）や数を示す情報を重畳表示する形態であってもよい。また、重畳表示に加えて、信号源（群）や数を示す情報を表示領域にスクロール表示してもよい。この場合、信号源が重畳表示されていない状態のスライス画像であっても、潜在的に信号源が紐付けられた状態であると考えることができるため、このスライス画像を表示する場合も、「所定の信号源が重畳された生体断層画像を表示領域に初期表示する」態様の一例であると考えることができる。

（第１の実施形態の変形例１）
本変形例では、図２２に示すように、各表示領域４１０Ａ～４１０Ｃで表示されるスライス画像は、重畳表示されるダイポールの推定結果の数が最も多いスライス画像が左側に配置されている。そして、左端のスライス画像から左右にスライス番号順に（層順に）並ぶように、他のスライス画像が配置されることは図１９に示す例と同じである。本変形例において、ダイポールの推定結果の数が最も多いスライス画像のスライス番号を図１９と同様にしたとき、図２２に示すようになる。例えば、表示領域４１０Ａには、スライス番号１０のスライス画像が左端に配置され、その右側にスライス番号１１，１２，１３，１４の順で配置されている。

また、表示領域４１０Ｂには、スライス番号１０のスライス画像が左端に配置され、その右側にスライス番号１１，１２，１３，１４の順で配置されている。そして、表示領域４１０Ｃには、スライス番号７のスライス画像が左端に配置され、その右側にスライス番号８，９，１０，１１，１２の順で配置されている。

このように、ダイポールの推定結果の数が最も多いスライス画像を左端に寄せて配置することで、ダイポールの推定結果の数が最も多いスライス画像を見つけやすくなる。

（第１の実施形態の変形例２）
例えば領域４０１Ａに初期表示するスライス画像の条件として、スライス画像に重畳される信号源の数が多い順に所定方向に沿ってスライス画像（生体断層画像）を並べて表示することもできる。つまり、第２の表示制御部２６２は、信号源の数が最も多いスライス画像（所定の信号源が重畳された生体断層画像）を基準として、信号源の数が多い順に所定方向に沿って他の生体断層画像を並べて表示することもできる。例えば図２３に示すように、第２の表示制御部２６２は、重畳される信号源の数が最も多いスライス画像を左端とし、左から右へ向けて（所定方向の一例）、重畳される信号源の数が少なくなるようにスライス画像を並べて表示する形態であってもよい。スライス画像上の信号源の数（ダイポールの推定結果の数）を確認するだけであれば、領域４０１Ａの左から右に向けて一方向に視線を移動するだけでよいので、図１９等に比べて視認性がよいという利点もある。また、図１９と同様に、信号源の数が最も多いスライス画像を中央に配置し、左右には信号源の数が徐々に少なくなるように並べて配置してもよい。

（第１の実施形態の変形例３）
図１９等の例では、領域４１０Ａ～領域４１０Ｃに３方向のスライス画像（スライス画像Ａ、スライス画像Ｂ、スライス画像Ｃ）を表示していた。本変形例では、第２の表示制御部２６２は、１方向のみのスライス画像を表示し、他の２方向のスライス画像は非表示にする。このときの領域４０１Ａに初期表示するスライス画像の条件としては、上記第１の実施形態や変形例１および２が適用可能である。例えば、領域４０１Ａにスライス画像を１列で表示すれば、１方向におけるスライス画像の視認性が向上する。また、スライス画像が多数ある場合は、図２４に示すように複数行にわたって１方向の断層画像を表示すれば、スクロールの回数を減らすことができ、より全体を把握でき、視認性が向上する。図２４は、上述の第１の実施形態の変形例２を適用した例で、信号源の数が最も多いスライス画像を中央に配置し、左右には信号源の数が徐々に少なくなるように並べたものであり、初期表示後にスクロールした状態を示している。なお、図２４の例では、１方向のみのスライス画像が表示される表示領域４１０の上に、スライス画像Ａに対応する方向（「Axial View」の方向）を選択するためのボタン５００Ａと、スライス画像Ｂに対応する方向（「Sagittal View」の方向）を選択するためのボタン５００Ｂと、スライス画像Ｃに対応する方向（「Coronal View」の方向）を選択するためのボタン５００Ｃとが配置される。第２の表示制御部２６２は、ボタン５００Ａの押下を受け付けた場合はスライス画像Ａ群のみを表示し、ボタン５００Ｂの押下を受け付けた場合はスライス画像Ｂ群のみを表示し、ボタン５００Ｃの押下を受け付けた場合はスライス画像Ｃ群のみを表示する。

（第１の実施形態の変形例４）
前述したように、第１の実施形態においては、表示領域４１０Ａに表示されたスライス画像Ａと、その直下に表示されたスライス画像Ｂ、および、さらにその直下に表示されたスライス画像Ｃとは３次元方向に対応していないが、本変形例のように、表示領域４１０Ａに表示されたスライス画像Ａと、その直下に表示されたスライス画像Ｂ、および、さらにその直下に表示されたスライス画像Ｃとが３次元方向に対応する形態であってもよい。この場合、表示領域４１０Ａ～４１０Ｃの何れかを基準とし、基準となる表示領域４１０に表示された複数のスライス画像のうち、重畳される信号源の数が最も多いスライス画像を中央に配置して表示し、中央のスライス画像から左右にスライス画像が層順（スライス番号順）に並ぶように、他のスライス画像を配置して表示した上で、他の表示領域４１０のスライス画像を対応させて表示することができる。例えば表示領域４１０Ａを基準とした場合、表示領域４１０Ａに表示されるスライス画像のうち、重畳される信号源の数が最も多いスライス画像を中央に配置して表示し、中央のスライス画像から左右にスライス画像が層順に並ぶように、他のスライス画像を配置して表示する。そして、他の表示領域４１０Ｂおよび表示領域４１０Ｃの各々に表示される複数のスライス画像の各々を、表示領域４１０Ａに表示された各スライス画像に対応させて表示するという具合である。このように、３次元的に対応する３つのスライス画像を上下方向に揃えて表示することで、立体的にダイポールの推定結果（信号源）の位置を把握できる。また、基準となる表示領域４１０に表示された複数のスライス画像のうちダイポールの推定結果の数が最も多いスライス画像を、上述の第１の実施形態の変形例１のように左端に寄せて配置してもよい。

（第１の実施形態の変形例５）
本変形例は、上述の第１の実施形態の変形例４において、基準となる表示領域４１０に表示された複数のスライス画像のうち、重畳される信号源の数が最も多いスライス画像を中央に配置して表示し、中央のスライス画像から左右に向けて信号源の数が少なくなるように他のスライス画像を配置して表示する。その上で、他の表示領域４１０のスライス画像を対応させて表示することができる。例えば表示領域４１０Ａを基準とした場合、表示領域４１０Ａに表示されるスライス画像のうち、重畳される信号源の数が最も多いスライス画像を中央に配置して表示し、中央のスライス画像から左右に向けて信号源の数が少なくなるようにスライス画像が並ぶように、他のスライス画像を配置して表示する。そして、他の表示領域４１０Ｂおよび表示領域４１０Ｃの各々に表示される複数のスライス画像の各々を、表示領域４１０Ａに表示された各スライス画像に対応させて表示するという具合である。このように、３次元的に対応する３つのスライス画像を上下方向に揃えて表示することで、立体的にダイポールの推定結果（信号源）の位置を把握できる。また、基準となる表示領域４１０に表示された複数のスライス画像のうちダイポールの推定結果の数が最も多いスライス画像を変形例１のように左端に寄せて配置してもよい。

（第１の実施形態の変形例６）
本変形例では、３次元方向に対応するスライス画像Ａ、スライス画像Ｂおよびスライス画像Ｃを１群とし、各群におけるダイポールの推定結果の数が最も多いスライス画像の群を中央に配置して表示し、中央のスライス画像の群から左右に向けて信号源の数が少なくなるように他のスライス画像の郡を配置して表示する。なお、ダイポールの推定結果の数が最も多いスライス画像の群を、上述の第１の実施形態の変形例１のように左端に寄せて配置してもよい。このように、３次元的に対応する３つのスライス画像の群を揃えて表示することで、立体的にダイポールの推定結果（信号源）の位置を把握できる。

（第１の実施形態の変形例７）
上述の第１の実施形態では、「only dipole」ボタン４３０Ｂと、「ALL」ボタン４３０Ａとにより、スライス画像の表示を切り替えていたが、例えばこれらのボタンが設けられずに、第２の表示制御部２６２は、ダイポールの推定結果が重畳表示されているスライス画像のみを表示する形態であってもよい。

（第１の実施形態の変形例８）
領域４０１に表示されるスライス画像には、スライス番号を示す情報４４０Ａを表示していたが、領域４０１で選択されたスライス画像に対応する断層位置ライン４５０が強調して表示される(例えば「赤」などの視認性の高い色で表示される等)形態であってもよい。また、スライス番号を示す情報４４０Ａに加えて、上記の強調表示する形態を組み合わせてもよい。

（第１の実施形態の変形例９）
上述した実施形態および各変形例では、１つのスライス画像に対応する信号源を当該スライス画像に重畳して領域４０１Ａに表示していた。本変形例では、複数のスライス画像にそれぞれ重畳されている各信号源を、所定のスライス画像に集約して重畳して領域４０１に表示する。例えば、複数のスライス画像についてスライス番号順に１０枚ずつ１つの群を構成し、群ごとに、当該群に含まれる各スライス画像に対応する信号源の数の合計を演算する。次に、各群の信号源の数の合計を比較し、図２５に示すように、その合計の数が最も多い群に含まれるスライス画像（のいずれか）を領域４０１Ａの中央に表示する。そして、中央のスライス画像から左右にスライス番号順に（層順に）並ぶように、他の群に含まれるスライス画像を配置して表示する。

ここで、領域４０１Ａに表示する各スライス画像は、群に含まれる各スライス画像のいずれかである。例えば、群に含まれる連続したスライス番号のスライス画像のうち、真中のスライス番号のスライス画像でもよく、スライス番号の一番小さい、または、一番大きいスライス画像であってもよい。そして、表示する画像として特定されたスライス画像に、当該スライス画像を含む群に含まれるスライス画像の全ての信号源を重畳する。また、表示するスライス画像のスライス番号が決まったら、その左右に表示するスライス画像のスライス番号は、群を構成するスライス画像の枚数だけシフトしたスライス番号とすればよい。例えば、１０枚のスライス画像で１つの群を構成する場合、左右に表示させるスライス画像のスライス番号は、中央に表示されたスライス画像のスライス番号に１０だけ加算したスライス番号、または減算したスライス番号とすればよい。

また、図２５に示すように、スライス番号を示す情報４４０Ａには、群を構成するスライス画像のスライス番号の範囲が表示される。また、信号源の数を示す情報４４０Ｂには、表示しているスライス画像を含む群に含まれる各スライス画像に対応する信号源の数の合計値が表示されている。なお、図２５では、領域４０１Ａの中央、および、その左右に１ずつのスライス画像のみを表示する例を示しているが、左右方向に複数のスライス画像が表示されていてもよい。本変形例によれば、領域４０１Ａに表示されるスライス画像の総数が減るため、上述した実施形態および各変形例に比べて閲覧性が向上する。なお、信号源の数の合計値が最も多い群のスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。

（第１の実施形態の変形例１０）
上述した第１の実施形態および各変形例では、スライス画像に重畳される信号源の数を、領域４０１Ａに初期表示するスライス画像の条件として用いたが、これに限られない。たとえば、信号源のベクトル方向や強度など、解析の目的にあった条件を適用することが可能である。

一例として、推定した信号源の妥当性もしくは信憑性、または信号源の近似の妥当性もしくは信憑性などを示す数値を用い、当該数値が最も高い信号源が重畳されたスライス画像を領域４０１Ａの中央に表示することができる。妥当性または信憑性など（以下、単に信憑性と総称する）を示す数値は、一例として、ＧＯＦ（Goodness Of Fit）を用いて算定することができる。そして、算定した信憑性を示す数値が所定の閾値を越えた信号源が重畳されているスライス画像を領域４０１Ａに表示する。そして、信憑性を示す数値が最も大きい（信憑性が最も高い）信号源が重畳されたスライス画像を領域４０１Ａに表示し、その左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示する。なお、信憑性を示す数値が最も大きいスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。

ここで、図２６は、上述の変形例９での図２５に示す表示形態を本変形例に適用した図を示している。本変形例では、情報４４０Ｂに示す数が信号源の数ではなく信憑性を示す数値である点が図２５と異なる。また、図２６に示す信憑性を示す数値は、１つの群に含まれるスライス画像に重畳する各信号源の信憑性の平均値であってもよく、または、１つの群に含まれるスライス画像に重畳する各信号源のうち、最も信憑性を示す数値が高い信号源の当該数値であってもよい。

本変形例によれば、信憑性の高い信号源を特定することで、症例の原因となる対象箇所を特定する確度がより向上する。また、情報４４０Ｂに示す数を、図２５に示すように信号源の数とし、信号源の表現（色、形、大きさなど）をＧＯＦの値に応じた表現とすることにより区別してもよい。

（第１の実施形態の変形例１１）
変形例１０では、信憑性を示す数値が所定の閾値以上の信号源をスライス画像に重畳させて表示させる動作を示したが、これに限られることはない。例えば、本変形例では、信憑性が低い画像も含めて表示させる。この場合、信憑性を示す数値が所定の閾値以上の信号源が重畳されたスライス画像は、信憑性が低いスライス画像と区別して表示させる。例えば、スライス画像の背景色を閾値以上と閾値未満とで変更させたり、スライス画像の外周枠の色を変更させたり、または、スライス画像内に注意喚起するマークを表示させてもよい。また、閾値以上の信号源と閾値未満の信号源とがスライス画像に重畳される場合は、閾値以上と閾値未満とで信号源の表現（色、形、大きさなど）を区別して表示させるものとしてもよい。

（第１の実施形態の変形例１２）
変形例１０とは別の例として、本変形例では、すべての信号源の平均座標を算出し、この平均座標に対応するスライス画像を領域４０１Ａの中央に配置する。そして、中央に表示したスライス画像の左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示する。これにより、各信号源の位置の平均座標からの広がりを視認しやすくなる。なお、平均座標に対応するスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。

さらに別の例として、信号源が存在する各スライス画像のスライス番号のうち、最大番号と最小番号との中心の番号を算出し、当該中心の番号に対応するスライス画像を領域４０１Ａの中央に表示し、その左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示するものとしてもよい。これにより、対象箇所の中心からの広がりを視認しやすくなる。この例においても、中心座標に対応するスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。

（第２の実施形態）
次に、第２の実施形態を説明する。上述の各実施形態と共通する部分については適宜に説明を省略する。基本的な装置構成は上述の第１の実施形態と同じである。上述の各実施形態では、解析画面には、所定の時間長分の生体データ（「１つの生体データ」とみなせる）が表示されているが、本実施形態では、第１の表示制御部２６１は、所定の時間長ごとに区切られた複数の生体データを表示対象とし、タイムゾーン１２０ｂに対応する何れかの生体データの信号波形を表示する。

また、解析部２５２（推定部）は、所定の時間長ごとに区切られた複数の生体データごとに、該生体データに対して入力済みの複数のアノテーションの中から選択されたアノテーションに対応する信号源を推定する。そして、第２の表示制御部２６２は、複数の生体データの一部に対応する信号源の数に基づいて、１以上の生体断層画像の表示を可変に制御する。この例では、第２の表示制御部２６２は、複数スライスされた生体断層画像に対して、所定の時間長ごとに区切られた複数の生体データの一部に対応する信号源を重畳して表示するとともに、上述の第１の実施形態と同様に、複数スライスされた生体断層画像のうち所定の信号源が重畳された生体断層画像を表示領域に初期表示するよう制御する。以下、具体的な内容を説明する。

＜測定収録時の動作＞
例えば、上述の第１の実施形態で説明したような測定動作が３回にわたって断続的に実施される場合を想定する。各測定の間には所定のインターバル（インターバルの時間は任意）が設けられるとする。なお、この「３回」という回数は一例であり、これに限られるものではない。要するに測定回数は検査目的に応じて任意に変更可能である。図２７は、この場合の情報処理装置５０の動作例（３回にわたる測定収録時の動作）を示すフローチャートである。図２７に示すように、ステップＳ４１において、情報処理装置５０は１回目の測定を行う。このときの動作は、図８のステップＳ１２～ステップＳ１７の処理と同様である。そして、１回目の測定が終了すると、情報処理装置５０は、１回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを、被検者を識別する被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４２）。

次に、情報処理装置５０は、２回目の測定を行う（ステップＳ４３）。このときの動作は、図８のステップＳ１２～ステップＳ１７の処理と同様である。そして、２回目の測定が終了すると、情報処理装置５０は、２回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４４）。

次に、情報処理装置５０は、３回目の測定を行う（ステップＳ４５）。このときの動作は、図８のステップＳ１２～ステップＳ１７の処理と同様である。そして、３回目の測定が終了すると、情報処理装置５０は、３回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４６）。

以上のようにして、１回の測定（所定の時間長にわたる測定）が完了するたびに、その測定結果を示す測定データがファイル単位で収録／解析情報保存部２５４に保存されていく。以下の説明では、収録／解析情報保存部２５４に保存された１つの測定データのファイルを「測定ファイル」と称する場合がある。この例では、３回にわたる測定が終了した後、収録／解析情報保存部２５４には３つの測定ファイルが記憶されることになる。以下では、１回目の測定に対応する測定ファイルを第１の測定ファイル、２回目の測定に対応する測定ファイルを第２の測定ファイル、３回目の測定に対応する測定ファイルを第３の測定ファイルと称する場合がある。上述したように、各測定ファイルに対しては被検者ＩＤが紐付けられて収録／解析情報保存部２５４に保存される。

＜解析時の動作＞
次に、解析時の動作を説明する。ここでは、図２の開始画面２０４で「解析」ボタンの押下を受け付けた場合、情報処理装置５０（第１の表示制御部２６１）は、測定により得られた測定ファイルを選択するための選択画面を表示装置２８に表示することを前提とする。図２８は、このときの情報処理装置５０の動作例を示すフローチャートである。

まず、情報処理装置５０（第１の表示制御部２６１）は、選択画面から測定ファイルを選択する操作を受け付ける（ステップＳ５１）。次に、情報処理装置５０（第１の表示制御部２６１）は、ステップＳ５１で選択された測定ファイルと、該測定ファイルに紐付く被検者ＩＤと同一の被検者ＩＤが紐付けられた他の１以上の測定ファイルと、を含む一連の測定ファイル（この例では上述の３つの測定ファイル）を読み込み、その読み込んだ一連の測定ファイルを反映させた解析画面を表示装置２８に表示する制御を行う（ステップＳ５２）。

図２９は、このときの解析画面の左側の領域２０２Ａの一例を示す図である。時間軸１２２上には、何れか１つの測定ファイルに対応する収録時間全体ではなく、一連の測定ファイル（第１の測定ファイル、第２の測定ファイル、第３の測定ファイル）の全ての収録時間を含む全体の時間が表示される。そして、時間軸１２２上には、第１の測定ファイルの収録期間を表す範囲情報９００ａと、第２の測定ファイルの収録期間を表す範囲情報９００ｂと、第３の測定ファイルの収録期間を表す範囲情報９００ｃと、が表示される。以下の説明では、範囲情報９００ａ、９００ｂ、９００ｃを互いに区別しない場合は単に「範囲情報９００」と称する場合がある。各範囲情報９００に対しては、対応する測定ファイルの名称を示す情報が付加される形態であってもよい。この例では、各測定は間隔を空けて実施されているため、時間軸１２２上において、範囲情報９００の間には隙間（空き領域）が設けられる。解析者は、マウス等によりタイムゾーン１２０ｂを移動する操作を行うことで、領域２０２Ａに表示される信号波形を切り替えることができる。この例では、タイムゾーン１２０ｂに対応する信号波形（何れかの測定ファイルの生体データのうちの一部）が領域２０２Ａに表示される。つまり、解析者は、時間軸１２２上でタイムゾーン１２０ｂを移動させることにより、測定ファイルを跨いで所望の時間帯の信号波形を領域２０２Ａに表示させることができる。また、ここでは、この解析画面の右側の領域２０２Ｂのアノテーションリスト１８０には、３つの測定ファイルの各々に含まれるアノテーションが全て表示される。また、例えば測定ファイルごとに対象となる検査の名称が紐付けられ、タイムゾーン１２０ｂに対応する検査の名称も併せて解析画面に表示する形態であってもよい。

図２８に戻って説明を続ける。情報処理装置５０（第１の表示制御部２６１）は、ステップＳ５２で解析画面を表示した後に、タイムゾーン１２０ｂの位置を変更する操作を受け付けた場合（ステップＳ５３：Ｙｅｓ）、変更後のタイムゾーン１２０ｂの位置に対応する信号波形が現在の領域２０２Ａに表示されているか否かを確認する（ステップＳ５４）。

ステップＳ５４の結果が否定の場合（ステップＳ５４：Ｎｏ）、情報処理装置５０（第１の表示制御部２６１）は、領域２０２Ａに表示される信号波形を、変更後のタイムゾーン１２０ｂの位置に対応する信号波形に切り替える（ステップＳ５５）。ステップＳ５４の結果が肯定の場合（ステップＳ５４：Ｙｅｓ）、または、ステップＳ５５の後、情報処理装置５０は、解析者の操作に応じて解析処理を実行する（ステップＳ５６）。この解析処理の内容は、図１５に示すステップＳ２３～ステップＳ３１の処理となる。なお、上述したように、アノテーションリスト１８０には、３つの測定ファイルの各々に含まれるアノテーションが全て表示されている。そして、ステップＳ２９において、マージボタン３００の押下を受け付けた場合、第２の表示制御部２６２は、アノテーションリスト１８０に表示された複数のアノテーション（複数の測定ファイルに跨る全てのアノテーション）のうち、推定完了マーク１８２が付与された複数のアノテーションと１対１に対応する複数の信号源の各々が重畳される複数の生体断層画像を、それぞれに重畳される信号源の数に応じて並べて表示する制御を行う。要するに、図１９に示す各スライス画像に重畳される信号源の集合は、複数の測定ファイルに跨る全てのアノテーションのうち、信号源の推定が行われた複数のアノテーションと１対１に対応する複数の信号源となる。

（第２の実施形態の変形例１）
例えば解析画面の時間軸１２２上には、何れか１つの測定ファイルに対応する範囲情報９００のみが表示され、解析者の操作に応じて、時間軸１２２上の範囲情報９００が測定ファイル単位で切り替わる形態であってもよい。図３０は、解析画面の一例を示す図である。図３０の例では、時間軸１２２上に、第１の測定ファイルの収録期間を表す範囲情報９００ａのみが表示されている。情報処理装置５０（第１の表示制御部２６１）は、解析者による範囲情報９００を切り替えるための操作を受け付けた場合、その受け付けた操作に応じて、時間軸１２２上の範囲情報９００を測定ファイル単位で切り替え、切り替え後の測定ファイルに対応させて表示領域２０２Ａおよび表示領域２０２Ｂの表示を切り替える。

（第２の実施形態の変形例２）
上述の第２の実施形態では、タイムゾーン１２０ｂの位置は、異なる範囲情報９００間を跨がないように設定される。例えばタイムゾーン１２０ｂが図２９に示す範囲情報９００ａの終点に位置している状態で、タイムゾーン１２０ｂの位置を１段階送る操作を受け付けた場合、情報処理装置５０（第１の表示制御部２６１）は、タイムゾーン１２０ｂが範囲情報９００ａと範囲情報９００ｂとの間を跨がずに次の範囲情報９００ｂの始点に位置するように、タイムゾーン１２０ｂの表示を切り替えるといった具合である。

ただし、これに限らず、本変形例では、タイムゾーン１２０ｂが異なる範囲情報９００間を跨ぐように配置されることも許容される。この場合、図３１の（Ａ）に示すように、タイムゾーン１２０ｂに対応する信号波形は、生体信号が存在しない空白領域（測定間の隙間に相当）が発生する。例えば異なる測定ファイル間の領域であることを視認し易くするため、図３１の（Ｂ）に示すように、情報処理装置５０（第１の表示制御部２６１）は、空白領域の背景を変更して表示する形態であってもよい。

また、例えばタイムゾーン１２０ｂが異なる範囲情報９００間を跨ぐように配置され、かつ、測定間の時間間隔が短い場合は、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形との隙間が殆ど無いことが想定される。このような場合、例えば図３１の（Ｃ）に示すように、情報処理装置５０（第１の表示制御部２６１）は、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形との間に、つなぎ目を表すライン（アノテーションラインとは異なるライン）を表示する形態であってもよい。さらに、例えば図３１の（Ｄ）に示すように、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形とを異なる表示としてもよいし、例えば図３１の（Ｅ）に示すように、何れかの測定ファイルに対応する信号波形の背景の色を変更する形態であってもよい。

（第２の実施形態の変形例３）
図３２に示すように、第２の実施形態の情報処理装置５０は、複数の測定ファイルごとに、解析結果（解析情報）を示す解析ファイルを紐付けて管理（保存）する。ここで、解析者は１人とは限らず、複数の解析者の各々が解析を行う場合も想定される。この場合、複数の解析者ごとに異なる解析ファイルが生成され、測定ファイルに紐付けられる。また、この例では、解析が終了するたびに、そのときの解析結果を示す解析ファイルが新たに測定ファイルに紐付けられていく。つまり、情報処理装置５０は、１回の解析のたびに、解析終了コマンドの入力を受け付けると、そのときの解析結果を示す解析ファイルを新たに測定ファイルに紐付けて保存する。例えば解析者が、解析画面上の「保存」または「終了」を示すボタンを押下することを契機として、解析終了コマンドが入力される形態であってもよい。

ここで、全てのファイルのアノテーションを解析画面のアノテーションリスト１８０に表示することとすると、現在ログインして解析を行っている解析者には、自身が解析を行っていないアノテーションも提示されることになり、全てのアノテーションの中から、自身が解析したアノテーションを探し出す必要があるために負担が大きい。

そこで、本変形例においては、情報処理装置５０は、各解析ファイルに対して、解析者と更新日（解析ファイルの作成日）と被検者ＩＤとを紐付けて管理する。そして、情報処理装置５０（第１の表示制御部２６１）は、図２の開始画面２０４で「解析」ボタンの押下を受け付けた場合、図３２のようなファイルリストに含まれる解析ファイルの中から、ログイン中の解析者に対応する解析ファイルを抽出し、複数の測定ファイルごとに該抽出した解析ファイルが紐付けられた選択画面を表示装置２８に表示する。なお、抽出した解析ファイルが紐付けられていない測定ファイルに関しては、測定ファイルのみが表示されることになる。例えばログイン中の解析者が解析者「Ａ」である場合、図３２のファイルリストに含まれる解析ファイルの中から、解析者「Ａ」に対応する解析ファイルが抽出され、図３３に示すような選択画面が表示されるといった具合である。

情報処理装置５０は、この選択画面から、何れかの測定ファイルの選択を受け付けた場合、その選択された測定ファイルと、該測定ファイルに紐付く被検者ＩＤと同一の被検者ＩＤが紐付けられた他の１以上の測定ファイルと、を含む一連の測定ファイルを反映させた解析画面を表示する。ここでは、この解析画面におけるアノテーションリスト１８０の表示に際して、情報処理装置５０（第１の表示制御部２６１）は、一連の測定ファイルごとに、該測定ファイルに紐付く解析ファイルを特定し、特定した解析ファイルに対応するアノテーションをアノテーションリスト１８０に表示する。また、１つの測定ファイルに対して、同じ解析者の解析ファイルが複数紐付けられている場合は、最新の更新日の解析ファイルのみを特定し、その特定した解析ファイルに対応するアノテーションをアノテーションリスト１８０に表示する。また、１つの測定ファイルに対して解析ファイルが紐付けられていない場合は、その測定ファイルに含まれるアノテーションをアノテーションリスト１８０に表示する。

また、例えば選択画面において、表示対象としたい測定ファイルを、解析者が全て選択する形態であってもよい。例えば図３３の選択画面において、解析者「Ａ」が測定ファイル２を選択した場合、情報処理装置５０（第１の表示制御部２６１）は、測定ファイル２のみを反映させた解析画面を表示する形態であってもよい。ここでは、この解析画面におけるアノテーションリスト１８０の表示に際して、情報処理装置５０（第１の表示制御部２６１）は、測定ファイル２に紐付く２つの解析ファイル２，２’のうち、最新の更新日の解析ファイル２’に対応するアノテーションをアノテーションリスト１８０に表示するといった具合である。

以上より、解析画面のアノテーションリスト１８０には、ログイン中の解析者に対応するアノテーションのみが適切に表示されることになるので、解析者の利便性が向上する。

なお、ログインとは、情報処理装置５０を使用する権限を有することを意味し、情報処理装置５０は、ユーザ（解析者）のログインの可否を判断する機能を有している。例えば情報処理装置５０の起動時には、ログインするための情報（例えばＩＤとパスワードとの組み合わせなどからなるログイン情報）の入力を促すログイン画面が表示され、解析者はログイン画面に対して、自身が所持するログイン情報の入力を行う。情報処理装置５０は、使用権限を有するユーザごとに、予め設定したログイン情報を紐付けて登録しておき、ログイン画面から入力されたログイン情報が、登録済みのログイン情報と一致する場合、ログイン情報の入力を行ったユーザ（解析者）のログインを許可し、一致しない場合はログインを許可しない。

以上、本発明に係る実施形態について説明したが、本発明は、上述の実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上述の実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、上述の実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態および変形例にわたる構成要素を適宜組み合わせてもよい。

また、上述した各実施形態の生体信号計測システム１で実行されるプログラムは、インストール可能な形式または実行可能な形式のファイルでＣＤ－ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ－Ｒ、ＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｋ）、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ Ｓｅｒｉａｌ Ｂｕｓ）等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよいし、インターネット等のネットワーク経由で提供または配布するように構成してもよい。また、各種プログラムを、ＲＯＭ等に予め組み込んで提供するように構成してもよい。

１ 生体信号計測システム
３ 測定装置
２１ ＣＰＵ
２２ ＲＡＭ
２３ ＲＯＭ
２４ 補助記憶装置
２５ 入出力インタフェース
２７ バス
２８ 表示装置
４０ サーバ
５０ 情報処理装置
２５０ 制御部
２５１ 表示制御部
２５２ 解析部
２５３ センサ情報取得部
２５４ 収録／解析情報保存部
２５５ アノテーション入力部
２６１ 第１の表示制御部
２６２ 第２の表示制御部

特開２０００－５１３３号公報

Claims (15)

1. 所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示可能な表示制御部を備え、
   該表示制御部は、
   前記複数スライスされた生体断層画像のうち、所定の条件に合致する信号源が重畳された生体断層画像を含めて、前記複数スライスされた生体断層画像を層順に並べて表示領域に表示するよう制御し、
   前記信号源が重畳されない生体断層画像の表示または非表示の指示に応じて、信号源が重畳されない生体断層画像の表示・非表示を切り替えることを特徴する、
   情報処理装置。
2. 前記所定の条件は、前記複数スライスされた生体断層画像毎の信号源の数が最も多いことを特徴とする、
   請求項１に記載の情報処理装置。
3. 前記表示制御部は、前記信号源の数が最も多い生体断層画像を、前記表示領域の中央部付近に表示することを特徴とする、
   請求項２に記載の情報処理装置。
4. 前記表示制御部は、前記中央部付近の生体断層画像から左右に生体断層画像が層順に並ぶように、他の生体断層画像を配置して表示する、
   請求項３に記載の情報処理装置。
5. 前記所定の条件は、前記複数スライスされた生体断層画像を、前記所定方向に所定数ずつの群に分けた場合において、前記群のうち、該群に含まれる生体断層画像に対応する各信号源の数の合計が最も多い群であることを特徴とする、
   請求項１に記載の情報処理装置。
6. 前記所定の条件は、前記複数スライスされた生体断層画像毎の信号源に対するＧＯＦ（Goodness Of Fit）が最も大きいことを特徴とする、
   請求項１に記載の情報処理装置。
7. 前記表示制御部は、前記表示の際に、前記複数スライスされた生体断層画像毎に重畳された信号源の座標の平均座標に対応する生体断層画像を前記表示領域に表示することを特徴とする、
   請求項１に記載の情報処理装置。
8. 前記表示制御部は、前記信号源の数が最も多い生体断層画像を基準として、前記信号源の数が多い順に所定方向に沿って他の生体断層画像を並べて表示する、
   請求項２に記載の情報処理装置。
9. 前記表示制御部は、前記信号源の数を示す情報を、対応する生体断層画像とともに表示する制御を行う、
   請求項１乃至８のうちの何れか１項に記載の情報処理装置。
10. 前記表示制御部は、前記表示領域に表示される生体断層画像の断層位置を示す情報を表示する制御を行う、
    請求項１乃至９のうちの何れか１項に記載の情報処理装置。
11. 前記複数の生体断層画像の各々は、
    前記所定方向としての第１方向の断面である第１断層画像と、前記第１方向と直交する第２方向の断面である第２断層画像と、前記第１方向および前記第２方向と直交する第３方向の断面である第３断層画像と、を含む、
    請求項１乃至１０のうちの何れか１項に記載の情報処理装置。
12. 前記表示制御部は、
    前記複数スライスされた生体断層画像に対して、所定の時間長ごとに区切られた複数の前記生体データの一部に対応する前記信号源を重畳して表示する、
    請求項１乃至１１のうちの何れか１項に記載の情報処理装置。
13. 所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示するにあたり、
    前記複数スライスされた生体断層画像のうち、所定の条件に合致する信号源が重畳された生体断層画像を含めて、層順に並べて表示領域に表示するよう制御する表示制御ステップを含み、
    前記信号源が重畳されない生体断層画像の前記信号源が重畳されない生体断層画像の表示または非表示の指示に応じて、信号源が重畳されない生体断層画像の表示・非表示を切り替えるステップを有することを特徴する、
    情報処理方法。
14. コンピュータに、
    所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示するにあたり、
    前記複数スライスされた生体断層画像のうち、所定の条件に合致する信号源が重畳された生体断層画像を含めて、前記複数スライスされた生体断層画像を層順に並べて表示領域に表示するよう制御する表示制御ステップと、前記信号源が重畳されない生体断層画像の表示または非表示の指示に応じて、信号源が重畳されない生体断層画像の表示・非表示を切り替えるステップを実行させるためのプログラム。
15. 被検者の１以上の生体信号を測定する測定装置と、前記測定装置で測定された１以上の前記生体信号を記録するサーバと、前記サーバに記録された１以上の前記生体信号を解析する情報処置装置と、を備えた生体信号計測システムであって、
    所定方向に複数スライスされた生体断層画像に対して、生体信号の経時的変化を示す生体データの一部に対応する信号源を重畳して表示可能な表示制御部を備え、
    該表示制御部は、
    前記複数スライスされた生体断層画像のうち、所定の条件に合致する信号源が重畳された生体断層画像を含めて、前記複数スライスされた生体断層画像を層順に並べて表示領域に表示するよう制御し、前記信号源が重畳されない生体断層画像の表示または非表示の指示に応じて、信号源が重畳されない生体断層画像の表示・非表示を切り替えることを特徴とする表示制御部を備える、
    生体信号計測システム。

Images