(57)【要約】【課題】複数時刻のダイポールの推定結果について視認性を高める。【解決手段】生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御部を備え、前記表示制御部は、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する。【選択図】図３５
2020146286.tif 000002
【特許請求の範囲】
【請求項１】
生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御部を備え、
前記表示制御部は、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する、
ことを特徴とする情報処理装置。
【請求項２】
前記表示制御部は、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果はグループごとに識別可能な色または形態で表示する、
ことを特徴とする請求項１に記載の情報処理装置。
【請求項３】
前記表示制御部は、グループ化されたダイポールの推定結果を表示する場合に、グループを代表するダイポールの推定結果を表示する、
ことを特徴とする請求項１に記載の情報処理装置。
【請求項４】
前記表示制御部は、グループを代表するダイポールの推定結果とともに、当該グループに関する情報を表示する、
ことを特徴とする請求項３に記載の情報処理装置。
【請求項５】
前記表示制御部は、グループに関する情報の表示を、ダイポールの推定結果の色または形態を用いて表現する、
ことを特徴とする請求項４に記載の情報処理装置。
【請求項６】
前記表示制御部は、ダイポールの推定結果をグループ化する際に、ダイポールの推定結果の信頼体積を用いてダイポールの推定結果間の位置の相違度を算出する、
ことを特徴とする請求項１ないし５の何れか一項に記載の情報処理装置。
【請求項７】
生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御ステップを含み、
前記表示制御ステップは、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する、
ことを特徴とする情報処理方法。
【請求項８】
コンピュータに、
生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御ステップを実行させ、
前記表示制御ステップは、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する、
ためのプログラム。
【請求項９】
被検者の１以上の生体信号を測定する測定装置と、前記測定装置で測定された１以上の前記生体信号を記録するサーバと、前記サーバに記録された１以上の前記生体信号を解析する情報処置装置と、を備えた生体信号計測システムであって、
生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御部を備え、
前記表示制御部は、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する、
ことを特徴とする生体信号計測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、情報処理方法、プログラムおよび生体信号計測システムに関する。
【背景技術】
【０００２】
従来、測定した被検者の生体信号に基づき推定した生体内の信号源を断層画像上に重畳表示するとともに、信号源に対応する生体信号の波形を画面上に並列表示する技術が知られている（例えば特許文献１参照）。この技術においては、オペレータは、表示されている生体信号の波形の所望の時間を指定すると、断層画像上の対応する信号源位置を特定する表示が行われる。
【０００３】
例えば、脳の神経活動を測定する脳磁計や脳波計では、測定した波形からてんかん特有の波形箇所（以下、特異点という。）を判別し、特異点から信号源の推定を行い、信号源のダイポールの推定結果を断層画像上に重畳表示させる。そして、この断層画像上の信号源の位置などに基づき、手術にて切除する位置（てんかんの原因となっている箇所）を特定することになる。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の技術では、推定した複数時刻のダイポールの推定結果を表示する際に、すべてのダイポールの推定結果を重畳表示しており、解析が困難な場合があるという問題があった。
【０００５】
本発明は、上記に鑑みてなされたものであって、複数時刻のダイポールの推定結果について視認性を高めることを目的とする。
【課題を解決するための手段】
【０００６】
上述した課題を解決し、目的を達成するために、本発明は、生体信号の経時的変化を示す生体データの一部に対応する信号源のダイポールの推定結果について同じ方向を向いたダイポールの推定結果をグループ化し、所定方向に複数スライスされた生体断層画像に対して、グループ化されたダイポールの推定結果を重畳して表示する表示制御部を備え、前記表示制御部は、グループ化されないダイポールの推定結果を表示する場合に、グループ化されたダイポールの推定結果とはダイポールの推定結果の方向に応じて色または形態を変えて表示する、ことを特徴とする。
【発明の効果】
【０００７】
本発明によれば、複数時刻のダイポールの推定結果について、グループ化を行って結果の視認性を高めることができる、という効果を奏する。
【図面の簡単な説明】
【０００８】
【図１】図１は、実施形態の生体信号計測システムの概略図である。
【図２】図２は、情報処理装置に表示される開始画面の一例を示す図である。
【図３】図３は、測定収録画面の一例を示す図である。
【図４】図４は、測定収録画面の左側の領域の拡大図である。
【図５】図５は、測定収録画面の右側の領域の拡大図である。
【図６】図６は、アノテーション情報が入力された直後の画面を示す図である。
【図７】図７は、更新されたアノテーションリストの図である。
【図８】図８は、測定収録時の情報処理装置の動作を示すフローチャートである。
【図９】図９は、解析画面の一例を示す図である。
【図１０】図１０は、解析画面の左側の領域の拡大図である。
【図１１】図１１は、解析画面の右側の領域の拡大図である。
【図１２】図１２は、解析画面で特定のアノテーションラインが選択された直後の画面を示す図である。
【図１３】図１３は、図１２の画面の左側の領域の拡大図である。
【図１４】図１４は、図１２の画面の右側の領域の拡大図である。
【図１５】図１５は、解析時の情報処理装置の動作を示すフローチャートである。
【図１６】図１６は、情報処理装置のハードウェア構成図である。
【図１７】図１７は、情報処理装置の機能ブロック図である。
【図１８】図１８は、表示制御部が有する機能の一例を示す図である。
【図１９】図１９は、第１の実施形態において、マージボタンが押下されたときに表示される画面を示す図である。
【図２０】図２０は、３次元方向に対応する３つのスライスの画像の関係を説明するための図である。
【図２１】図２１は、「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」ボタン４３０Ａとにより、スライス画像の表示を切り替えていたが、例えばこれらのボタンが設けられずに、第２の表示制御部２６２は、ダイポールの推定結果が重畳表示されているスライス画像のみを表示する形態であってもよい。
【００９４】
（第１の実施形態の変形例８）
領域４０１に表示されるスライス画像には、スライス番号を示す情報４４０Ａを表示していたが、領域４０１で選択されたスライス画像に対応する断層位置ライン４５０が強調して表示される（例えば「赤」などの視認性の高い色で表示される等）形態であってもよい。また、スライス番号を示す情報４４０Ａに加えて、上記の強調表示する形態を組み合わせてもよい。
【００９５】
（第１の実施形態の変形例９）
上述した実施形態および各変形例では、１つのスライス画像に対応する信号源を当該スライス画像に重畳して領域４０１Ａに表示していた。本変形例では、複数のスライス画像にそれぞれ重畳されている各信号源を、所定のスライス画像に集約して重畳して領域４０１に表示する。例えば、複数のスライス画像についてスライス番号順に１０枚ずつ１つの群を構成し、群ごとに、当該群に含まれる各スライス画像に対応する信号源の数の合計を演算する。次に、各群の信号源の数の合計を比較し、図２５に示すように、その合計の数が最も多い群に含まれるスライス画像（のいずれか）を領域４０１Ａの中央に表示する。そして、中央のスライス画像から左右にスライス番号順に（層順に）並ぶように、他の群に含まれるスライス画像を配置して表示する。
【００９６】
ここで、領域４０１Ａに表示する各スライス画像は、群に含まれる各スライス画像のいずれかである。例えば、群に含まれる連続したスライス番号のスライス画像のうち、真中のスライス番号のスライス画像でもよく、スライス番号の一番小さい、または、一番大きいスライス画像であってもよい。そして、表示する画像として特定されたスライス画像に、当該スライス画像を含む群に含まれるスライス画像の全ての信号源を重畳する。また、表示するスライス画像のスライス番号が決まったら、その左右に表示するスライス画像のスライス番号は、群を構成するスライス画像の枚数だけシフトしたスライス番号とすればよい。例えば、１０枚のスライス画像で１つの群を構成する場合、左右に表示させるスライス画像のスライス番号は、中央に表示されたスライス画像のスライス番号に１０だけ加算したスライス番号、または減算したスライス番号とすればよい。
【００９７】
また、図２５に示すように、スライス番号を示す情報４４０Ａには、群を構成するスライス画像のスライス番号の範囲が表示される。また、信号源の数を示す情報４４０Ｂには、表示しているスライス画像を含む群に含まれる各スライス画像に対応する信号源の数の合計値が表示されている。なお、図２５では、領域４０１Ａの中央、および、その左右に１ずつのスライス画像のみを表示する例を示しているが、左右方向に複数のスライス画像が表示されていてもよい。本変形例によれば、領域４０１Ａに表示されるスライス画像の総数が減るため、上述した実施形態および各変形例に比べて閲覧性が向上する。なお、信号源の数の合計値が最も多い群のスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。
【００９８】
（第１の実施形態の変形例１０）
上述した第１の実施形態および各変形例では、スライス画像に重畳される信号源の数を、領域４０１Ａに初期表示するスライス画像の条件として用いたが、これに限られない。たとえば、信号源のベクトル方向や強度など、解析の目的にあった条件を適用することが可能である。
【００９９】
一例として、推定した信号源の妥当性もしくは信憑性、または信号源の近似の妥当性もしくは信憑性などを示す数値を用い、当該数値が最も高い信号源が重畳されたスライス画像を領域４０１Ａの中央に表示することができる。妥当性または信憑性など（以下、単に信憑性と総称する）を示す数値は、一例として、ＧＯＦ（Good  Of  Fitness）を用いて算定することができる。そして、算定した信憑性を示す数値が所定の閾値を越えた信号源が重畳されているスライス画像を領域４０１Ａに表示する。そして、信憑性を示す数値が最も大きい（信憑性が最も高い）信号源が重畳されたスライス画像を領域４０１Ａに表示し、その左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示する。なお、信憑性を示す数値が最も大きいスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。
【０１００】
ここで、図２６は、上述の変形例９での図２５に示す表示形態を本変形例に適用した図を示している。本変形例では、情報４４０Ｂに示す数が信号源の数ではなく信憑性を示す数値である点が図２５と異なる。また、図２６に示す信憑性を示す数値は、１つの群に含まれるスライス画像に重畳する各信号源の信憑性の平均値であってもよく、または、１つの群に含まれるスライス画像に重畳する各信号源のうち、最も信憑性を示す数値が高い信号源の当該数値であってもよい。
【０１０１】
本変形例によれば、信憑性の高い信号源を特定することで、症例の原因となる対象箇所を特定する確度がより向上する。また、情報４４０Ｂに示す数を、図２５に示すように信号源の数とし、信号源の表現（色、形、大きさなど）をＧＯＦの値に応じた表現とすることにより区別してもよい。
【０１０２】
（第１の実施形態の変形例１１）
変形例１０では、信憑性を示す数値が所定の閾値以上の信号源をスライス画像に重畳させて表示させる動作を示したが、これに限られることはない。例えば、本変形例では、信憑性が低い画像も含めて表示させる。この場合、信憑性を示す数値が所定の閾値以上の信号源が重畳されたスライス画像は、信憑性が低いスライス画像と区別して表示させる。例えば、スライス画像の背景色を閾値以上と閾値未満とで変更させたり、スライス画像の外周枠の色を変更させたり、または、スライス画像内に注意喚起するマークを表示させてもよい。また、閾値以上の信号源と閾値未満の信号源とがスライス画像に重畳される場合は、閾値以上と閾値未満とで信号源の表現（色、形、大きさなど）を区別して表示させるものとしてもよい。
【０１０３】
（第１の実施形態の変形例１２）
変形例１０とは別の例として、本変形例では、すべての信号源の平均座標を算出し、この平均座標に対応するスライス画像を領域４０１Ａの中央に配置する。そして、中央に表示したスライス画像の左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示する。これにより、各信号源の位置の平均座標からの広がりを視認しやすくなる。なお、平均座標に対応するスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。
【０１０４】
さらに別の例として、信号源が存在する各スライス画像のスライス番号のうち、最大番号と最小番号との中心の番号を算出し、当該中心の番号に対応するスライス画像を領域４０１Ａの中央に表示し、その左右にはスライス番号順に（層順に）並ぶように、他のスライス画像を配置して表示するものとしてもよい。これにより、対象箇所の中心からの広がりを視認しやすくなる。この例においても、中心座標に対応するスライス画像を領域４０１Ａの中央に表示する形態に限られず、上述した各変形例に示す表示形態を本変形例に適用することもできる。
【０１０５】
（第２の実施形態）
次に、第２の実施形態を説明する。上述の各実施形態と共通する部分については適宜に説明を省略する。基本的な装置構成は上述の第１の実施形態と同じである。上述の各実施形態では、解析画面には、所定の時間長分の生体データ（「１つの生体データ」とみなせる）が表示されているが、本実施形態では、第１の表示制御部２６１は、所定の時間長ごとに区切られた複数の生体データを表示対象とし、タイムゾーン１２０ｂに対応する何れかの生体データの信号波形を表示する。
【０１０６】
また、解析部２５２（推定部）は、所定の時間長ごとに区切られた複数の生体データごとに、該生体データに対して入力済みの複数のアノテーションの中から選択されたアノテーションに対応する信号源を推定する。そして、第２の表示制御部２６２は、複数の生体データの一部に対応する信号源の数に基づいて、１以上の生体断層画像の表示を可変に制御する。この例では、第２の表示制御部２６２は、複数スライスされた生体断層画像に対して、所定の時間長ごとに区切られた複数の生体データの一部に対応する信号源を重畳して表示するとともに、上述の第１の実施形態と同様に、複数スライスされた生体断層画像のうち所定の信号源が重畳された生体断層画像を表示領域に初期表示するよう制御する。以下、具体的な内容を説明する。
【０１０７】
＜測定収録時の動作＞
例えば、上述の第１の実施形態で説明したような測定動作が３回にわたって断続的に実施される場合を想定する。各測定の間には所定のインターバル（インターバルの時間は任意）が設けられるとする。なお、この「３回」という回数は一例であり、これに限られるものではない。要するに測定回数は検査目的に応じて任意に変更可能である。図２７は、この場合の情報処理装置５０の動作例（３回にわたる測定収録時の動作）を示すフローチャートである。図２７に示すように、ステップＳ４１において、情報処理装置５０は１回目の測定を行う。このときの動作は、図８のステップＳ１２〜ステップＳ１７の処理と同様である。そして、１回目の測定が終了すると、情報処理装置５０は、１回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを、被検者を識別する被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４２）。
【０１０８】
次に、情報処理装置５０は、２回目の測定を行う（ステップＳ４３）。このときの動作は、図８のステップＳ１２〜ステップＳ１７の処理と同様である。そして、２回目の測定が終了すると、情報処理装置５０は、２回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４４）。
【０１０９】
次に、情報処理装置５０は、３回目の測定を行う（ステップＳ４５）。このときの動作は、図８のステップＳ１２〜ステップＳ１７の処理と同様である。そして、３回目の測定が終了すると、情報処理装置５０は、３回目の測定により得られた生体データと、入力されたアノテーションとを含む測定データを被検者ＩＤと紐付けて収録／解析情報保存部２５４に保存する（ステップＳ４６）。
【０１１０】
以上のようにして、１回の測定（所定の時間長にわたる測定）が完了するたびに、その測定結果を示す測定データがファイル単位で収録／解析情報保存部２５４に保存されていく。以下の説明では、収録／解析情報保存部２５４に保存された１つの測定データのファイルを「測定ファイル」と称する場合がある。この例では、３回にわたる測定が終了した後、収録／解析情報保存部２５４には３つの測定ファイルが記憶されることになる。以下では、１回目の測定に対応する測定ファイルを第１の測定ファイル、２回目の測定に対応する測定ファイルを第２の測定ファイル、３回目の測定に対応する測定ファイルを第３の測定ファイルと称する場合がある。上述したように、各測定ファイルに対しては被検者ＩＤが紐付けられて収録／解析情報保存部２５４に保存される。
【０１１１】
＜解析時の動作＞
次に、解析時の動作を説明する。ここでは、図２の開始画面２０４で「解析」ボタンの押下を受け付けた場合、情報処理装置５０（第１の表示制御部２６１）は、測定により得られた測定ファイルを選択するための選択画面を表示装置２８に表示することを前提とする。図２８は、このときの情報処理装置５０の動作例を示すフローチャートである。
【０１１２】
まず、情報処理装置５０（第１の表示制御部２６１）は、選択画面から測定ファイルを選択する操作を受け付ける（ステップＳ５１）。次に、情報処理装置５０（第１の表示制御部２６１）は、ステップＳ５１で選択された測定ファイルと、該測定ファイルに紐付く被検者ＩＤと同一の被検者ＩＤが紐付けられた他の１以上の測定ファイルと、を含む一連の測定ファイル（この例では上述の３つの測定ファイル）を読み込み、その読み込んだ一連の測定ファイルを反映させた解析画面を表示装置２８に表示する制御を行う（ステップＳ５２）。
【０１１３】
図２９は、このときの解析画面の左側の領域２０２Ａの一例を示す図である。時間軸１２２上には、何れか１つの測定ファイルに対応する収録時間全体ではなく、一連の測定ファイル（第１の測定ファイル、第２の測定ファイル、第３の測定ファイル）の全ての収録時間を含む全体の時間が表示される。そして、時間軸１２２上には、第１の測定ファイルの収録期間を表す範囲情報９００ａと、第２の測定ファイルの収録期間を表す範囲情報９００ｂと、第３の測定ファイルの収録期間を表す範囲情報９００ｃと、が表示される。以下の説明では、範囲情報９００ａ、９００ｂ、９００ｃを互いに区別しない場合は単に「範囲情報９００」と称する場合がある。各範囲情報９００に対しては、対応する測定ファイルの名称を示す情報が付加される形態であってもよい。この例では、各測定は間隔を空けて実施されているため、時間軸１２２上において、範囲情報９００の間には隙間（空き領域）が設けられる。解析者は、マウス等によりタイムゾーン１２０ｂを移動する操作を行うことで、領域２０２Ａに表示される信号波形を切り替えることができる。この例では、タイムゾーン１２０ｂに対応する信号波形（何れかの測定ファイルの生体データのうちの一部）が領域２０２Ａに表示される。つまり、解析者は、時間軸１２２上でタイムゾーン１２０ｂを移動させることにより、測定ファイルを跨いで所望の時間帯の信号波形を領域２０２Ａに表示させることができる。また、ここでは、この解析画面の右側の領域２０２Ｂのアノテーションリスト１８０には、３つの測定ファイルの各々に含まれるアノテーションが全て表示される。また、例えば測定ファイルごとに対象となる検査の名称が紐付けられ、タイムゾーン１２０ｂに対応する検査の名称も併せて解析画面に表示する形態であってもよい。
【０１１４】
図２８に戻って説明を続ける。情報処理装置５０（第１の表示制御部２６１）は、ステップＳ５２で解析画面を表示した後に、タイムゾーン１２０ｂの位置を変更する操作を受け付けた場合（ステップＳ５３：Ｙｅｓ）、変更後のタイムゾーン１２０ｂの位置に対応する信号波形が現在の領域２０２Ａに表示されているか否かを確認する（ステップＳ５４）。
【０１１５】
ステップＳ５４の結果が否定の場合（ステップＳ５４：Ｎｏ）、情報処理装置５０（第１の表示制御部２６１）は、領域２０２Ａに表示される信号波形を、変更後のタイムゾーン１２０ｂの位置に対応する信号波形に切り替える（ステップＳ５５）。ステップＳ５４の結果が肯定の場合（ステップＳ５４：Ｙｅｓ）、または、ステップＳ５５の後、情報処理装置５０は、解析者の操作に応じて解析処理を実行する（ステップＳ５６）。この解析処理の内容は、図１５に示すステップＳ２３〜ステップＳ３１の処理となる。なお、上述したように、アノテーションリスト１８０には、３つの測定ファイルの各々に含まれるアノテーションが全て表示されている。そして、ステップＳ２９において、マージボタン３００の押下を受け付けた場合、第２の表示制御部２６２は、アノテーションリスト１８０に表示された複数のアノテーション（複数の測定ファイルに跨る全てのアノテーション）のうち、推定完了マーク１８２が付与された複数のアノテーションと１対１に対応する複数の信号源の各々が重畳される複数の生体断層画像を、それぞれに重畳される信号源の数に応じて並べて表示する制御を行う。要するに、図１９に示す各スライス画像に重畳される信号源の集合は、複数の測定ファイルに跨る全てのアノテーションのうち、信号源の推定が行われた複数のアノテーションと１対１に対応する複数の信号源となる。
【０１１６】
（第２の実施形態の変形例１）
例えば解析画面の時間軸１２２上には、何れか１つの測定ファイルに対応する範囲情報９００のみが表示され、解析者の操作に応じて、時間軸１２２上の範囲情報９００が測定ファイル単位で切り替わる形態であってもよい。図３０は、解析画面の一例を示す図である。図３０の例では、時間軸１２２上に、第１の測定ファイルの収録期間を表す範囲情報９００ａのみが表示されている。情報処理装置５０（第１の表示制御部２６１）は、解析者による範囲情報９００を切り替えるための操作を受け付けた場合、その受け付けた操作に応じて、時間軸１２２上の範囲情報９００を測定ファイル単位で切り替え、切り替え後の測定ファイルに対応させて表示領域２０２Ａおよび表示領域２０２Ｂの表示を切り替える。
【０１１７】
（第２の実施形態の変形例２）
上述の第２の実施形態では、タイムゾーン１２０ｂの位置は、異なる範囲情報９００間を跨がないように設定される。例えばタイムゾーン１２０ｂが図２９に示す範囲情報９００ａの終点に位置している状態で、タイムゾーン１２０ｂの位置を１段階送る操作を受け付けた場合、情報処理装置５０（第１の表示制御部２６１）は、タイムゾーン１２０ｂが範囲情報９００ａと範囲情報９００ｂとの間を跨がずに次の範囲情報９００ｂの始点に位置するように、タイムゾーン１２０ｂの表示を切り替えるといった具合である。
【０１１８】
ただし、これに限らず、本変形例では、タイムゾーン１２０ｂが異なる範囲情報９００間を跨ぐように配置されることも許容される。この場合、図３１の（Ａ）に示すように、タイムゾーン１２０ｂに対応する信号波形は、生体信号が存在しない空白領域（測定間の隙間に相当）が発生する。例えば異なる測定ファイル間の領域であることを視認し易くするため、図３１の（Ｂ）に示すように、情報処理装置５０（第１の表示制御部２６１）は、空白領域の背景を変更して表示する形態であってもよい。
【０１１９】
また、例えばタイムゾーン１２０ｂが異なる範囲情報９００間を跨ぐように配置され、かつ、測定間の時間間隔が短い場合は、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形との隙間が殆ど無いことが想定される。このような場合、例えば図３１の（Ｃ）に示すように、情報処理装置５０（第１の表示制御部２６１）は、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形との間に、つなぎ目を表すライン（アノテーションラインとは異なるライン）を表示する形態であってもよい。さらに、例えば図３１の（Ｄ）に示すように、一方の測定ファイルに対応する信号波形と、他方の測定ファイルに対応する信号波形とを異なる表示としてもよいし、例えば図３１の（Ｅ）に示すように、何れかの測定ファイルに対応する信号波形の背景の色を変更する形態であってもよい。
【０１２０】
（第２の実施形態の変形例３）
図３２に示すように、第２の実施形態の情報処理装置５０は、複数の測定ファイルごとに、解析結果（解析情報）を示す解析ファイルを紐付けて管理（保存）する。ここで、解析者は１人とは限らず、複数の解析者の各々が解析を行う場合も想定される。この場合、複数の解析者ごとに異なる解析ファイルが生成され、測定ファイルに紐付けられる。また、この例では、解析が終了するたびに、そのときの解析結果を示す解析ファイルが新たに測定ファイルに紐付けられていく。つまり、情報処理装置５０は、１回の解析のたびに、解析終了コマンドの入力を受け付けると、そのときの解析結果を示す解析ファイルを新たに測定ファイルに紐付けて保存する。例えば解析者が、解析画面上の「保存」または「終了」を示すボタンを押下することを契機として、解析終了コマンドが入力される形態であってもよい。
【０１２１】
ここで、全てのファイルのアノテーションを解析画面のアノテーションリスト１８０に表示することとすると、現在ログインして解析を行っている解析者には、自身が解析を行っていないアノテーションも提示されることになり、全てのアノテーションの中から、自身が解析したアノテーションを探し出す必要があるために負担が大きい。
【０１２２】
そこで、本変形例においては、情報処理装置５０は、各解析ファイルに対して、解析者と更新日（解析ファイルの作成日）と被検者ＩＤとを紐付けて管理する。そして、情報処理装置５０（第１の表示制御部２６１）は、図２の開始画面２０４で「解析」ボタンの押下を受け付けた場合、図３２のようなファイルリストに含まれる解析ファイルの中から、ログイン中の解析者に対応する解析ファイルを抽出し、複数の測定ファイルごとに該抽出した解析ファイルが紐付けられた選択画面を表示装置２８に表示する。なお、抽出した解析ファイルが紐付けられていない測定ファイルに関しては、測定ファイルのみが表示されることになる。例えばログイン中の解析者が解析者「Ａ」である場合、図３２のファイルリストに含まれる解析ファイルの中から、解析者「Ａ」に対応する解析ファイルが抽出され、図３３に示すような選択画面が表示されるといった具合である。
【０１２３】
情報処理装置５０は、この選択画面から、何れかの測定ファイルの選択を受け付けた場合、その選択された測定ファイルと、該測定ファイルに紐付く被検者ＩＤと同一の被検者ＩＤが紐付けられた他の１以上の測定ファイルと、を含む一連の測定ファイルを反映させた解析画面を表示する。ここでは、この解析画面におけるアノテーションリスト１８０の表示に際して、情報処理装置５０（第１の表示制御部２６１）は、一連の測定ファイルごとに、該測定ファイルに紐付く解析ファイルを特定し、特定した解析ファイルに対応するアノテーションをアノテーションリスト１８０に表示する。また、１つの測定ファイルに対して、同じ解析者の解析ファイルが複数紐付けられている場合は、最新の更新日の解析ファイルのみを特定し、その特定した解析ファイルに対応するアノテーションをアノテーションリスト１８０に表示する。また、１つの測定ファイルに対して解析ファイルが紐付けられていない場合は、その測定ファイルに含まれるアノテーションをアノテーションリスト１８０に表示する。
【０１２４】
また、例えば選択画面において、表示対象としたい測定ファイルを、解析者が全て選択する形態であってもよい。例えば図３３の選択画面において、解析者「Ａ」が測定ファイル２を選択した場合、情報処理装置５０（第１の表示制御部２６１）は、測定ファイル２のみを反映させた解析画面を表示する形態であってもよい。ここでは、この解析画面におけるアノテーションリスト１８０の表示に際して、情報処理装置５０（第１の表示制御部２６１）は、測定ファイル２に紐付く２つの解析ファイル２，２’のうち、最新の更新日の解析ファイル２’に対応するアノテーションをアノテーションリスト１８０に表示するといった具合である。
【０１２５】
以上より、解析画面のアノテーションリスト１８０には、ログイン中の解析者に対応するアノテーションのみが適切に表示されることになるので、解析者の利便性が向上する。
【０１２６】
なお、ログインとは、情報処理装置５０を使用する権限を有することを意味し、情報処理装置５０は、ユーザ（解析者）のログインの可否を判断する機能を有している。例えば情報処理装置５０の起動時には、ログインするための情報（例えばＩＤとパスワードとの組み合わせなどからなるログイン情報）の入力を促すログイン画面が表示され、解析者はログイン画面に対して、自身が所持するログイン情報の入力を行う。情報処理装置５０は、使用権限を有するユーザごとに、予め設定したログイン情報を紐付けて登録しておき、ログイン画面から入力されたログイン情報が、登録済みのログイン情報と一致する場合、ログイン情報の入力を行ったユーザ（解析者）のログインを許可し、一致しない場合はログインを許可しない。
【０１２７】
（第３の実施の形態）
次に、第３の実施の形態について説明する。
【０１２８】
第３の実施の形態は、数時刻のダイポール推定の結果について、グループ化を行って結果の視認性を高めるとともに、グループ化の手段を複数用意して表示を切り替え可能とする点が、第１の実施の形態ないし第２の実施の形態と異なる。以下、第３の実施の形態の説明では、第１の実施の形態ないし第２の実施の形態と同一部分の説明については省略し、第１の実施の形態ないし第２の実施の形態と異なる箇所について説明する。
【０１２９】
例えば図２０においては、複数の時刻において算出されたダイポールの推定結果１９０ａをすべて表示した例を示した。図２０のスライス画像Ａ〜Ｃについては、全て同じ方向のダイポールの推定結果１９０ａが表示されている。しかしながら、例えば図２０のスライス画像Ａ〜Ｃに対して、図３４にａとして示すようなダイポールの推定結果１９０ａが追加された場合、図３４に示すように非常に視認性が悪く、見落とす可能性がある。
【０１３０】
ここで、図３５は第３の実施の形態にかかるダイポールの推定結果１９０ａを表示したスライスの画像の一例を示す図である。上述のような事態を防止すべく、本実施の形態において、第２の表示制御部２６２は、図３５に示すように、ダイポールの推定結果１９０ａの方向により色を変えて表示する。なお、図３５は、分かりやすさのために平面での例を示している。第２の表示制御部２６２は、ダイポールの推定結果１９０ａの向きが図３５の右下に示す色環のどこを示すかを求め、示した先の色環の色でダイポールの推定結果１９０ａを着色する。なお、第２の表示制御部２６２は、３次元の場合でも同様に色立体を用いてダイポールの推定結果１９０ａの色を定めれば良い。
【０１３１】
ただし、図３５のように表示した場合でも、ダイポールの推定結果１９０ａを重畳する上下関係によっては、図３６に示すように、一部のダイポールの推定結果１９０ａの視認性が悪くなる場合がある。
【０１３２】
そこで、第２の表示制御部２６２は、多くのダイポールの推定結果１９０ａが重畳されることで生じる視認性の低下を防ぐために、ダイポールの推定結果１９０ａのグループ化処理を行う。図３６に示す例の場合、多くのダイポールの推定結果１９０ａが同じ方向を向いており、これをまとめて表示することで視認性の向上を図ることができる。
【０１３３】
ここで、図３７はグループ化処理を行ったダイポールの推定結果１９０ａを表示したスライスの画像の一例を示す図である。図３７に示した例は、図３６に示した同じ方向を向いたダイポールの推定結果１９０ａをグループ化したものである。
【０１３４】
なお、グループ化（クラスタリング）については様々な手法が提案されており、如何なる手法を用いてグループ化を行っても構わない。診察のためには以下が肝要となる。
・グループ化の手法が選択できること。
・グループ化に用いるパラメータが調整できること。すなわち、更に纏めたり、分散させたりすることができること。
【０１３５】
本実施の形態においては、一例としてダイポールの推定結果１９０ａの信頼体積を用いたグループ化を採用する（米国特許７，８４０，０３９号など参照）（。ここで、図３８はダイポールの推定結果１９０ａの信頼体積を用いたグループ化を説明する図である。個々のダイポールの推定結果１９０ａの位置については、統計的にそのもっともらしさを算出することが出来る。例えば、ダイポールの推定結果１９０ａの位置の９５％存在範囲（９５％の確率でその範囲にダイポールの推定結果１９０ａの正しい位置が存在する）を計算することが可能で、その範囲を図３８に示すように楕円で表示することが可能である。この信頼体積を用いてダイポールの推定結果１９０ａの位置の相違を算出することで、単純に物理的な距離を用いて位置の相違を算出するのと比べて、より適切に（統計的に）位置の相違を評価できる。
【０１３６】
図３９は、信頼体積を用いて位置の相違（距離）を算出する一例を示す図である。図３９では、ダイポールの推定結果１９０ａの方向を示す棒線を割愛しており、中心位置を示す円（塗りつぶし）と信頼体積を示す楕円（線分のみ）を記載している。通常、ダイポールの推定結果１９０ａ間の位置の相違は、Ａ−Ｂの距離を持って評価している。しかしながら、ダイポールの推定結果１９０ａの位置の信頼性はダイポールの推定結果１９０ａごとに異なる。この信頼性を示すものが、図３９において楕円で示されている信頼体積である。本実施の形態では、信頼体積を用いて相違度を算出することでより確かな位置の相違を評価する。
【０１３７】
図３９に示すように、線分Ａ−ＢまたはＡ−Ｂの延長線と信頼体積の楕円が交わる点をＡｂ，Ｂａとする。ここでは相違度Ｄを以下の式で定義する。
Ｄ＝（Ａ−Ｂ）／（Ａ−Ａｂ）＋（Ａ−Ｂ）／（Ｂ−Ｂａ）
ただし、Ｘ−Ｙは点Ｘと点Ｙの間の距離
【０１３８】
そして、相違度Ｄを位置の相違度としてグループ化を行う。グループ化の方法についてはＫ−ｍｅａｎ法、Ｗａｒｄ法等さまざまなものがあり、そのいずれを用いても構わない。Ｋ−ｍｅａｎ法のような非階層的手法の場合には、初期シード点を設定する必要がある。この初期シード点の設定には、ダイポールの推定結果１９０ａのもっともらしさ（ＧＯＦ）を参照して、ＧＯＦの大きなダイポールの推定結果１９０ａを初期シードとする方法などが考えられる。
【０１３９】
また、グループ化にあたってはダイポールの推定結果１９０ａの位置、方向、強さの情報を用いことができる。脳活動解析の観点からは、同じ位置のダイポールの推定結果１９０ａであっても、方向が異なるものは別のグループと判断されることが好ましい。これは、脳の構造的には同じ位置であれば同じ方向のダイポールの推定結果１９０ａが検出されるはずであるが、位置推定の精度が悪いために異なる方向のダイポールの推定結果１９０ａが同じ位置に検出されていると解釈されるためである。
【０１４０】
そこで、本実施の形態においては、ダイポールの推定結果１９０ａのグループ化にあたり、同じ位置であっても方向が異なる場合には別のグループとする。例えば、以下の手順で段階的にグループ分けを行えば、同じ位置のダイポールの推定結果１９０ａでも別のグループに分類できる。
（１）．ダイポールの推定結果１９０ａの方向に着目してグループ化を行う。
（２）．（１）各グループに対して、ダイポールの推定結果１９０ａの位置に着目してグループ化を行う。（１）で別のグループに入ったダイポールの推定結果１９０ａは最終的に同じグループとはならない。
【０１４１】
なお、グループ分けの処理に用いるパラメータ（例えば、各種クラスタリングの終息条件）は、ユーザが適宜調整できることが好ましい。
【０１４２】
次に、第２の表示制御部２６２によるグループ分けを行ったダイポールの推定結果１９０ａの表示方法について説明する。
【０１４３】
第２の表示制御部２６２は、一例として、図３５に示した色分けをグループに基づいて行う。ただし、図４０に示すように、発生源の異なるダイポールの推定結果１９０ａが同じ（あるいは見分けが困難な）色で表示される場合がある。
【０１４４】
そこで、第２の表示制御部２６２は、グループごとに識別しやすい色を割り当てるようにする。ここで、図４１はグループごとに識別しやすい色を割り当てた例を示す図である。図４１に示すように、グループ分けを行い、そのグループごとに識別しやすい色を割り当てたならば、視認性よく異なる発生源のダイポールの推定結果１９０ａを見分けることが出来る。
【０１４５】
また、第２の表示制御部２６２は、全てのダイポールの推定結果１９０ａを表示するのではなく、グループを代表するダイポールの推定結果１９０ａのみを表示するようにしてもよい。代表するダイポールの推定結果１９０ａは、グループ内で位置的に最もグループの中心にあるダイポールの推定結果１９０ａや、グループ内で最も信頼性の高いダイポールの推定結果１９０ａ等とすれば良い。または、グループ内のダイポールの推定結果１９０ａについて位置、方向それぞれ平均をとって、平均位置に平均方向を持つダイポールの推定結果１９０ａにすれば良い。いずれにしても、第２の表示制御部２６２は、グループを代表して表示されるダイポールの推定結果１９０ａとして適切なダイポールの推定結果１９０ａを選択／算出する。
【０１４６】
ここで、図４２はグループを代表するダイポールの推定結果１９０ａのみを表示した例を示す図、図４３はグループを代表するダイポールの推定結果１９０ａのみを表示した別の例を示す図である。図４２に示すように、第２の表示制御部２６２は、ダイポールの推定結果１９０ａの位置を示す円の中に、グループに属するダイポールの推定結果１９０ａの数を表示する。さらに、図４３に示すように、第２の表示制御部２６２は、ダイポールの推定結果１９０ａの位置を示す円の大きさを、グループ内のダイポールの推定結果１９０ａの個数により変更して表示する。
【０１４７】
上述したように、グループごとに代表のダイポールの推定結果１９０ａのみを表示する場合には、グループに関する情報（ダイポールの推定結果１９０ａの数、位置のばらつき、方向のばらつき、信頼性（ＧＯＦ）の平均、ばらつきなど）を適宜ダイポールの推定結果１９０ａと紐づけて表示することが好ましい。特に、ダイポールの推定結果１９０ａの方向の情報は診断の際に重要であり、方向のばらつき等の統計的な情報は分かりやすく表示されていることが必須である。
【０１４８】
例えば、論文「Sensorimotor  seizures  of  pediatric  onset  with  unusual  posteriorly  oriented  rolandic  spikes.  Kakisaka  Y1,  Nakasato  N,  Haginoya  K,  Kanno  A,  Tsuchiya  S.  Epilepsy  Res.  2009  Apr;84（2-3）:153-8.  doi:  10.1016/j.eplepsyres.2009.01.012.  Epub  2009  Feb  28.」では、ダイポールの推定結果１９０ａの方向によりてんかんの良性／悪性を判断している。
【０１４９】
図４４は、グループ内の方向のばらつきの表示例を示す図である。図４４に示すように、第２の表示制御部２６２は、ダイポールの推定結果１９０ａの方向を示す線分（台形）の先端の広がり具合でグループ内の方向のばらつきを示している。
【０１５０】
以上、説明を簡単にするために２次元の３面図を用いて説明したが、３Ｄ画像上にダイポールの推定結果１９０ａを重畳する場合でも同様に処理を行うことが出来る。３Ｄ画像上にダイポールの推定結果１９０ａを重畳した場合の例を図４５に示す。
【０１５１】
また、第２の表示制御部２６２は、グループ内のアノテーションが存在する全てのスライスに代表ダイポールの推定結果１９０ａを表示するようにしてもよい。この場合には、スライス上にどの程度のアノテーションが乗っていたかが分かることが好ましい。ここで、図４６はスライス上にどの程度のダイポールの推定結果１９０ａが乗っているかを示す図である。図４６に示す例では、表示されている代表ダイポールの推定結果１９０ａについて、グループ内の６つのダイポールの推定結果１９０ａのうち２つのダイポールの推定結果１９０ａがこのスライス上に乗っていたことを、ダイポールの推定結果１９０ａの位置を示す円内の「２／６」の表示で示している。
【０１５２】
また、図４７は図４６に示す表示を図１９に示した画面に反映した例を示す図である。なお、図２１〜図２６に示した画面についても、同様に代表ダイポールの推定結果１９０ａの表示に変更できる。
【０１５３】
なお、グループの代表ダイポールの推定結果１９０ａを表示する場合、図１９〜２６の表示については以下のように対応すれば良い。第２の表示制御部２６２は、スライス上に乗るダイポールの推定結果１９０ａの数等を用いて中央に配置するスライスやスライスのソート順序を決定している。代表のダイポールの推定結果１９０ａのみを表示する場合には、以下のように取り扱えば同様の処理を行うことができる。
（１）代表のダイポールの推定結果１９０ａの位置に、グループ内に存在するダイポールの推定結果１９０ａの数が存在するものとする。
（２）応用として、グループ内のばらつきが小さい代表ダイポールの推定結果１９０ａが乗っているスライスを中央に配置するなど、グループに関する情報を利用して中央に表示するスライスや、スライスのソート順を決めることも可能となる。第１の実施の形態では、１つのアノテーションの情報のみを利用していたが、グループ内の複数のアノテーションに関する統計的な情報が同様に利用可能となる。
【０１５４】
このように本実施の形態によれば、複数時刻のダイポールの推定結果について、グループ化を行って結果の視認性を高めることができる。また、グループ化の手段を複数用意して表示を切り替え可能とすることにより、解析者が臨機応変に表示を切り替えることで、複数時刻のダイポールの推定結果に対する解析を容易かつ適切に行うことができる。
【０１５５】
以上、本発明に係る実施形態について説明したが、本発明は、上述の実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上述の実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、上述の実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態および変形例にわたる構成要素を適宜組み合わせてもよい。
【０１５６】
また、上述した各実施形態の生体信号計測システム１で実行されるプログラムは、インストール可能な形式または実行可能な形式のファイルでＣＤ−ＲＯＭ、フレキシブルディスク（ＦＤ）、ＣＤ−Ｒ、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｋ）、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ  Ｓｅｒｉａｌ  Ｂｕｓ）等のコンピュータで読み取り可能な記録媒体に記録して提供するように構成してもよいし、インターネット等のネットワーク経由で提供または配布するように構成してもよい。また、各種プログラムを、ＲＯＭ等に予め組み込んで提供するように構成してもよい。
【符号の説明】
【０１５７】
１      生体信号計測システム
３      測定装置
２１    ＣＰＵ
２２    ＲＡＭ
２３    ＲＯＭ
２４    補助記憶装置
２５    入出力インタフェース
２７    バス
２８    表示装置
４０    サーバ
５０    情報処理装置
２５０  制御部
２５１  表示制御部
２５２  解析部
２５３  センサ情報取得部
２５４  収録／解析情報保存部
２５５  アノテーション入力部
２６１  第１の表示制御部
２６２  第２の表示制御部
【先行技術文献】
【特許文献】
【０１５８】
【特許文献１】特開２０００−５１３３号公報
【図１】
2020146286.tif 000003
【図２】
2020146286.tif 000004
【図３】
2020146286.tif 000005
【図４】
2020146286.tif 000006
【図５】
2020146286.tif 000007
【図６】
2020146286.tif 000008
【図７】
2020146286.tif 000009
【図８】
2020146286.tif 000010
【図９】
2020146286.tif 000011
【図１０】
2020146286.tif 000012
【図１１】
2020146286.tif 000013
【図１２】
2020146286.tif 000014
【図１３】
2020146286.tif 000015
【図１４】
2020146286.tif 000016
【図１５】
2020146286.tif 000017
【図１６】
2020146286.tif 000018
【図１７】
2020146286.tif 000019
【図１８】
2020146286.tif 000020
【図１９】
2020146286.tif 000021
【図２０】
2020146286.tif 000022
【図２１】
2020146286.tif 000023
【図２２】
2020146286.tif 000024
【図２３】
2020146286.tif 000025
【図２４】
2020146286.tif 000026
【図２５】
2020146286.tif 000027
【図２６】
2020146286.tif 000028
【図２７】
2020146286.tif 000029
【図２８】
2020146286.tif 000030
【図２９】
2020146286.tif 000031
【図３０】
2020146286.tif 000032
【図３１】
2020146286.tif 000033
【図３２】
2020146286.tif 000034
【図３３】
2020146286.tif 000035
【図３４】
2020146286.tif 000036
【図３５】
2020146286.tif 000037
【図３６】
2020146286.tif 000038
【図３７】
2020146286.tif 000039
【図３８】
2020146286.tif 000040
【図３９】
2020146286.tif 000041
【図４０】
2020146286.tif 000042
【図４１】
2020146286.tif 000043
【図４２】
2020146286.tif 000044
【図４３】
2020146286.tif 000045
【図４４】
2020146286.tif 000046
【図４５】
2020146286.tif 000047
【図４６】
2020146286.tif 000048
【図４７】
2020146286.tif 000049