(57)【特許請求の範囲】
【請求項１】
被検体から収集した一連の時系列データを設定された解析期間にわたって所定サンプリング時間毎に解析して複数の解析データを求める手段と、
前記解析に関するパラメータとしての電流源の強度、深さ、および方向のうち少なくとも１つの度合いと、表示形態としての色相、輝度、およびコントラストのうち少なくとも１つとを関連付ける手段と、
前記度合いの序列を、段階的な前記表示形態の変化として示すバー又はサークルを表示すると共に、前記複数の解析データのそれぞれを前記パラメータと関連付けられた表示形態で表示する手段と、を設け、
前記時系列データを、磁気センサにより前記被検体としての生体から計測した磁場波形のデータとし、前記解析データを、前記磁場波形のデータを解析して得た前記生体内の電流源のデータとすることを特徴とする医用画像表示装置。
【請求項２】
請求項１に記載の発明において、
前記表示処理手段は、磁気共鳴イメージング装置、Ｘ線ＣＴスキャナを含む医用モダリティの中の特定のモダリティにより収集された前記被検体の形態画像上に前記複数の解析データを、前記所定サンプリング時間に対応して位置合わせして重畳して表示することを特徴とする医用画像表示装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、医用画像を表示する表示装置に係り、とくに、生体磁気計測装置やファンクショナルＭＲＩ（ｆＭＲＩ）装置で計測した医用情報の表示に好適な医用画像の表示装置に関する。
【０００２】
【従来の技術】
医用モダリティで収集した医用画像をいかに表示して見せるかは、読影能率や診断能の向上、さらには、医師の操作上、診断上の労力軽減の観点からも重要なファクタになっている。
【０００３】
このことは、生体磁気計測装置についても言えることである。生体磁気計測装置は、生体の脳や心臓などの部位内の電流源から発生する微弱な磁気（磁場）信号を捕らえて、それらの部位の機能診断に利用するもので、近年、その開発や研究が盛んに行われている。従来良く知られているように、例えば人の心臓は、その動きに伴って電流を発生する。この電気信号は非常に微弱であるが、種々の方法により非侵襲的に測定することができる。その１つの方法が生体磁気計測であり、これは生体中を流れる電流に因って生じる磁気（磁場）を測定するものである。とくに、脳を対象とした生体磁気計測を脳磁界計測と、心臓を対象としたそれを心磁界計測と別けて呼ぶことがある。
【０００４】
この生体磁気計測装置により得られた磁場データから電流源の強さ（電流の強さ）、電流源の向き、電流源の場所（位置）、および／または、それらの時間変化が推定（解析）され、表示される。これにより、生体内のどの部分が活性化しているかを知ることができ、診断に有益な情報を得ることができる。
【０００５】
生体磁気計測装置に係る従来の表示法の内、第１の表示法としては、生体磁場の時間変化の様子を表す表示法が採用されている。時間軸上で数ｍｓｅｃ毎にサンプリングしたときの磁場源解析データを用いて、この結果を形態画像上に矢印や線で連結しながらプロットしていく表示法である。一例として、図２４に示すように、丸（○）印に棒（−）線を付した模擬的な電流源モデルを採用している。棒（−）線の部分が電流の方向を示し、丸（○）印の部分が電流源の位置を示している。
【０００６】
また、第２の表示法としては、生体内を仮想的に格子状に区切った格子点上に電流モーメントを有する微小な電流源が在ると仮定し、それらの電流源の大きさや向きなどを推定し、表示する方法が知られている。
【０００７】
これらの推定結果は、ＭＲＩ（磁気共鳴イメージング）像やＸ線ＣＴ像に重畳して表示される。
【０００８】
【発明が解決しようとする課題】
しかしながら、上述した表示法の内、第１の表示法によれば、矢印や線などが邪魔になりがちで、脳や心臓の解析に係る一連の活動がどこから始まってどこで終わっているのか認識し難いという問題があった。つまり、時間的に変化する磁場波形の時間的前後関係を一目で把握することができなかったり、電流源の時間的な動きを観察しづらいという問題があった。
【０００９】
本発明は、上述した第１の表示法に関わる問題に鑑みてなされたもので、磁場の変化の時間的前後関係、詳細には、電流源の位置、深さ、強度、および／または向きの時間的前後関係を容易に且つ迅速に把握できる表示法を提供することを、その目的とする。
【００１０】
【課題を解決するための手段】
上記目的を達成するため、本発明の医用画像表示装置は、被検体から収集した一連の時系列データを設定された解析期間にわたって所定サンプリング時間毎に解析して複数の解析データを求める手段と、前記解析に関するパラメータとしての電流源の強度、深さ、および方向のうち少なくとも１つの度合いと、表示形態としての色相、輝度、およびコントラストのうち少なくとも１つとを関連付ける手段と、前記度合いの序列を、段階的な前記表示形態の変化として示すバー又はサークルを表示すると共に、前記複数の解析データのそれぞれを前記パラメータと関連付けられた表示形態で表示する手段と、を設け、前記時系列データを、磁気センサにより前記被検体としての生体から計測した磁場波形のデータとし、前記解析データを、前記磁場波形のデータを解析して得た前記生体内の電流源のデータとすることを特徴とする。
【００１１】
好適には、前記収集データは、磁気センサにより前記被検体としての生体から計測した磁場波形のデータであり、前記解析データは前記磁場波形のデータを解析して得た前記生体内の電流源のデータである。
【００１３】
本発明のそのほかの特徴および効果は、以下に説明する発明の実施の形態および添付の図面により明らかになる。
【００１４】
【発明の実施の形態】
以下、本発明の実施の形態を、添付図面を参照しながら説明する。
【００１５】
以下の各実施形態は、本発明に係る医用画像の表示装置を生体磁気計測装置に適用した例に関する。最初に、この生体磁気計測装置の全体構成の概略を図１に基づき説明する。
【００１６】
図１に示す生体磁気計測装置は、単チャンネルのｄｃ−ＳＱＵＩＤ磁束計をセンサの基本構成とし、これを複数チャンネル化した構成を採用している。なお、図１には、単チャンネルに関する構成を中心に示す。
【００１７】
各チャンネルのｄｃ−ＳＱＵＩＤ磁束計は、一例として、変調型ＦＬＬ（Flux Locked Loop）回路方式で構成されており、検出部としてのＳＱＵＩＤ磁束センサ１１と、このＳＱＵＩＤ磁束センサ１１を駆動して測定出力を得る駆動回路１２とを備える。駆動回路１２の出力はプリプロセッサ１３、Ａ／Ｄ変換器１４、および解析装置１５を介して表示装置１６に送られる。駆動回路１２、プリプロセッサ１３、Ａ／Ｄ変換器１４、解析装置１５などのユニットの動作に必要な情報は制御装置１７から与えられる。制御装置１７には、オペレータが必要な指令を与えるための入力装置１８が接続されている。
【００１８】
ＳＱＵＩＤ磁束センサ１１は、液体ヘリウムを入れたデュア内に浸され、超伝導状態におかれる。このＳＱＵＩＤ磁束センサ１１は、超伝導リングを有するセンサ本体部２１と、生体が発生する微弱な磁界を検知するピックアップコイル２２（検出コイル）とを備える。センサ本体部２１は、超伝導リングとしてのＳＱＵＩＤリング２１ａ、磁束を入力させるインプットコイル２１ｂ、およびフィードバックコイル２１ｃを備え、これらが一体的に形成されている。インプットコイル２１ｂにはピックアップコイル２２が超伝導接続される。ピックアップコイル２２の配線材としてはＮｂＴｉ線などが用いられ、この線材が種々の形式に基づき巻装・設置される。
【００１９】
駆動回路１２は、直流バイアス電源ＥB 、増幅器２３、同期検波器２４、励振信号発信器２５を備える。直流バイアス電源ＥB はＳＱＵＩＤリング２１ａに直流バイアス電流を供給する。励振信号発信器２５の出力信号ｆ0 はフィードバックコイル２１ｃに供給される。ＳＱＵＩＤリング２１ａの出力電圧は、増幅器２３を介して同期検波器２４に送られる。これにより、ＳＱＵＩＤリング２１ａのジョセフソン接合部分の電圧が変化しないようにフィードバックが掛けられ、このフィードバック信号が生体磁気の計測信号（つまり、その計測チャンネルにより計測された磁場の時間変化波形の信号：アナログ量）として取り出される。
【００２０】
このように非線形なＳＱＵＩＤ磁束センサ１１のセンサ信号が駆動回路１２を介して線形化される。駆動回路１２の計測データはプリプロセッサ１３に送られ、このプロセッサにより所望のゲインが与えられ、且つ、所望帯域に設定される。このプリプロセッサ１３からの出力信号はＡ／Ｄ変換器１４でデジタル信号に変換され、解析装置１５に送られる。
【００２１】
解析装置１５は例えばＣＰＵなどの演算素子を有して構成される。この解析装置１５において、計測データは、各サンプリング時刻における電流源（磁場源）の情報の推定演算（つまり、電流源解析処理）、および、本特許に係る画像表示の処理・演算（つまり、画像表示処理）に付されるとともに、必要に応じて、平滑化処理や空間フィルタ処理などの後処理に付される。制御装置１７は、上述した画像表示処理を行うために、プリプロセッサ１３およびＡ／Ｄ変換器１４にサンプリング時間を指定したり、プリプロセッサ１３のゲインや帯域を制御するなど、各ユニットに必要な情報を提供する。
【００２２】
第１の実施形態
第１の実施形態に係る生体磁気計測装置を上述した図１、及び、図２〜図６に基づき説明する。
【００２３】
上述した解析装置１５は、画像表示処理として、図２，３，６に示す処理を行う。具体的にはまず、図２において、入力装置１８を介してオペレータから与えられる情報に基づき、表示方法を決めるとともに、上述したサンプリング時間、ゲイン、帯域などの動作上の必要項目を設定する（図２、ステップＳ１）。
【００２４】
また、計測データの時系列上での解析範囲を選択する（ステップＳ２）。これは、画面上で解析範囲を設定するか、または、解析範囲を示すデータを入力器１８から与えることで行われる。
【００２５】
次いで解析装置１５は、設定した解析範囲の計測データを読み出し、その範囲内の指定サンプリング時間毎に、種々のアルゴリズムにより電流源に関わる諸量（電流強度、向き、位置、深さ）を解析（推定演算）する（ステップＳ３）。次いで、各サンプリング時間毎の解析位置の位置座標を決める（ステップＳ４）。つまり、電流源の推定位置はサンプリング時間毎に変化するから、その変化する位置をベクトル合成などにより決める。
【００２６】
さらに、解析装置１５は、選択された表示法に基づき、画像表示処理を実行して表示データを作成する（ステップＳ５）。この処理の一例を図３に示す。画像表示処理が済むと、表示データは表示装置１６に供給され、目視可能に画面表示される（ステップＳ６）。
【００２７】
本実施形態における画像表示処理を図３に例示する。解析装置１５は、最初に、計測した磁場波形における経過時間と色相とを一義的に対応させたテーブルのデータを読み出す（ステップＳ５−１）。このデータは予め設定され、メモリに記憶している。次いで、このテーブルを参照して、その磁場波形のデータを、これに対応した色相データに変換する（ステップＳ５−２）。この色相データは、電流源モデルとしてのある一定半径の円を埋めるとともに、各サンプリング時間に解析された電流源の強度に応じて、その輝度を変調させたデータに変換される（ステップＳ５−３）。これにより、例えば図４に示すように、ある磁場波形の時間変化（同図（ａ）参照）が、テーブルデータ（同図（ｂ）参照）を介して、時間ｔの変化に応じて変わる電流源モデルのデータ（同図（ｃ）参照）に変換される。なお、電流源の強度に応じて輝度を変える代わりに、コントラストを変えてもよい。
【００２８】
このため、表示装置１６の画面には、図５に示す如く、例えばＭＲＩ装置やＸ線ＣＴスキャナなどの他のモダリティで得られた頭部ＨＤの２次元断面像上に、ステップＳ４で決定された各サンプリング時間毎の位置に、ステップＳ５で加工された一連の電流源モデルＭＤが重畳表示される。同時に、画面上の適宜な位置に時間経過と色相とを対応させたカラーバーＣＢが分割表示される。
【００２９】
この画面上で、例えば赤系の最高の色相階調度で表示される解析開始位置（または活動の開始位置）から、例えば青系の最高の色相階調度で表示される解析終了位置（または活動の終了位置）まで経過時間と伴に色相が徐々に変わって表示される。また、生体内電流源（磁場源）の位置変化に応じて表示位置も徐々に移り変わる。この一連の時間変化が色相で色分けして表示されるので、その開始位置から終了位置までの視認または目視による追跡も容易になり、その強度の変化も輝度（またはコントラスト）から容易に把握できる。
【００３０】
さらに、この表示状態において、制御装置１７および解析装置１５は図６に示す処理を実行する。制御装置１７は、入力装置１８からオペレータが操作したカーソルやボタンの操作信号を入力し、監視している（ステップＳ１１，Ｓ１２）。操作信号からオペレータが表示画面上のカラーバーＣＢの特定の色相領域がクリックされたと判断されると、そのクリック情報を解析装置１５に送る。これに応答した解析装置１５は、クリックされた色相（すなわち、経過時間）に対応する電流源モデルが特定される（ステップＳ１３）。次いで、解析装置１５は、その電流源モデルの輝度（またはコントラスト）を変化させるデータを作成し、これを表示装置１６に送る（ステップＳ１４）。この輝度（またはコントラスト）の変化は、制御装置１７がその表示解除の指示をオペレータから受けるまで続けられる（ステップＳ１５〜Ｓ１７）。
【００３１】
この結果、例えば、カラーバーＣＢの橙色が指定され、この磁場波形において時刻ｔｎに相当していたとする（図４（ａ），（ｂ）参照）。この場合、時刻ｔｎに相当する橙色の電流源モデルが特定される（図４（ｃ）参照）。このため、画面上において、この色相を持つ円状の電流源モデル（例えば図５における円ＯＲ）内の輝度が例えば明滅する。これにより、オペレータがカラーバーＣＢ上で指定した任意の色相（すなわち時刻）が、活動部位の変化の軌跡線上においてどの位置の電流源モデルに対応するのか、容易に認識することができる。したがって、診断的に興味ある時間に対応した部位を容易に探すことができる。
【００３２】
変形例（その１）
第１の実施形態に係る生体磁気計測装置による表示法の１つの変形例を図７，８に基づき説明する。この変形例は、電流源モデルの分布表示に関する。
【００３３】
上述の実施形態では、電流源の推定位置を全て同一の大きさの円で表わすように構成したが、これを更に変形することもできる。一例として、電流源モデルを、電流源がその場所に存在する確率を反映させた大きさ又は領域で表すこともできる。
【００３４】
例えば、ダイポラリティ（Ｄ）という概念を用いることができる。計測した磁界分布をＢmea 、電流源推定を行った結果として得た電流源から計算した磁界分布をＢcal とすると、ダイポラリティＤは、
【数１】
0004229489.tif 000002
で表される。ここで、Ｂmea ＝Ｂcal のとき、Ｄ＝１となり、推定した電流源分布が測定した電流源分布に完全に一致する。このダイポラリティという概念を用いて電流源の分布表示を行うことができる。これに代えて、電流源の存在確率密度という概念を用いて電流源の分布表示を行うようにしてもよい。
【００３５】
これを行うには、解析装置１５は図７の処理を行う。同図の処理において、ステップＳ５−１１，１２，１４の処理は前述した図３のステップＳ５−１〜３の処理と同一である。解析装置１５は、ステップＳ５−１３にて、上述した電流源分布表示のための例えばダイポラリティＤを各サンプリング時間毎に演算し、演算結果に応じて電流源モデルの大きさを変更する。
【００３６】
この結果、例えば図８に示すように、電流源モデルＭＤの大きさがダイポラリティＤに応じて変更して画面表示される。このとき、電流源モデルＭＤと時間経過との関係および電流源の磁場強度と電流源モデルの輝度との関係は前述と同様である。したがって、この変形例によって、電流源の存在確率までをも取り込んだ表示が可能になり、表示情報の豊富化を図ることができる。
【００３７】
変形例（その２）
第１の実施形態に係る生体磁気計測装置による表示法の別の変形例を図９に基づき説明する。
【００３８】
この変形例は、磁場波形に対する各サンプリング時間毎の電流源の分布を等高線表示したもので、分布、すなわち等高線の変化に応じて色相、コントラストの濃淡、または輝度を一定刻みで連続的に変化させたグラデーション（階調）を持たせたものである。図９に表示例を示す。このための処理は、解析装置１５によって実施される。
【００３９】
なお、等高線表示のみならず、ドットの分布によって電流源分布を表示するようにしてもよい（図１０参照）。
【００４０】
第２の実施形態
本発明の第２の実施形態に係る生体磁気計測装置を図１１〜図１３に基づき説明する。なお、以下の実施形態およびその変形例において、第１の実施形態の装置と同一または同等の構成要素には同一符号を用いて、その説明を省略または簡略化する。
【００４１】
この実施形態の装置は、電流源の強度と色相の関係を表示に取り入れている。これを行うため、解析装置１５は、前述した図２の処理を行い、そのステップＳ５における画像表示処理の詳細な一例として、図１１に記載の処理を行う。また解析装置１５は前述した図６に記載の処理も合わせて実行する。
【００４２】
図１１の画像表示処理によれば、解析装置１５は、電流源の電流強度と色相とを一義的に対応させた対応関係を表すテーブルのデータを読み出す（ステップＳ５−２１）。このデータは予め設定され、メモリに記憶している。次いで、このテーブルを参照して、ステップＳ３で解析された電流源強度の経時的なデータを、これに対応した色相データに変換する（ステップＳ５−２２）。この色相データは表示用データとして表示装置１６に送られる。
【００４３】
これにより、例えば図１２に示すように、推定された電流強度の時間変化（同図（ａ）参照）が、テーブルデータ（同図（ｂ）参照）を介して、強度の変化に応じて変わる色相データ（同図（ｃ）参照）に変換される。この色相データは、表示時には、電流源モデルとしての一定半径の円内を埋めるデータとして用いられる。
【００４４】
このため、表示装置１６の画面には、図１３に示す如く、例えばＭＲＩ装置やＸ線ＣＴスキャナなどの他のモダリティで得られた頭部ＨＤの２次元断面像上に、ステップＳ４で決定された各サンプリング時間毎の推定位置に、ステップＳ５（ステップＳ５−２１，２２）で加工された一連の電流源モデルＭＤが重畳表示される。同時に、画面上の適宜な位置に電源源の強度と色相とを対応させたカラーバーＣＢが参照用として分割表示される。
【００４５】
これにより、画面上では、解析開始位置（または活動の開始位置）から解析終了位置（または活動の終了位置）まで強度に応じて色相が変わって表示される。また、生体内電流源（磁場源）の位置変化に応じて電流源モデルの表示位置も徐々に移り変わる。通常、電流源の位置は変化し且つその強度も変化するので、画面上の電流源モデルもその位置を変えながらかつ表示色相も経時的に変えられる。したがって、その開始位置から終了位置までの視認または目視による追跡も容易になり、その強度の変化も色相から容易に把握できる。
【００４６】
さらに、この表示状態において、制御装置１７および解析装置１５は図６に示す処理を実行するので、第１の実施形態のときと同様に、カラーバーＣＢ内の指定色相（すなわち、強度）に呼応して、その色相に相当する電流源モデル円の種々の情報を表す輝度、コントラストなどを変化させ、その位置を容易に把握できる。
【００４７】
第３の実施形態
第３の実施形態に係る生体磁気計測装置を図１４〜１５に基づき説明する。
【００４８】
この装置は、電流源の深さ（深さ方向の位置）と色相の関係を表示に取り入れている。これを行うため、解析装置１５は、前述した図２の処理を行い、そのステップＳ５における画像表示処理の詳細な一例として、図１４に記載の処理を行う。また解析装置１５は前述した図６に記載の処理も合わせて実行する。
【００４９】
図１４の画像表示処理によれば、解析装置１５は、電流源の深さと色相とを一義的に対応させた対応関係を表すテーブルのデータを読み出す（ステップＳ５−３１）。このデータは予め設定され、メモリに記憶している。次いで、このテーブルを参照して、ステップＳ３で解析された電流源の深さの経時的なデータを色相データに変換する（ステップＳ５−３２）。この色相データは表示用データとして表示装置１６に送られる。
【００５０】
これにより、推定された電流源の深さのデータが、テーブルデータを介して、深さの変化に応じて変わる色相データに変換される。この色相データは、表示時に、電流源モデルとしての一定半径の円内を埋めるデータとして用いられる。
【００５１】
このため、表示装置１６の画面には、図１５に示す如く、例えばＭＲＩ装置やＸ線ＣＴスキャナなどの他のモダリティで得られた頭部ＨＤの２次元断面像上に、ステップＳ４で決定された各サンプリング時間毎の位置に、ステップＳ５（ステップＳ５−３１，３２）で加工された一連の電流源モデルＭＤが重畳表示される。同時に、画面上の適宜な位置に電流源の方向と色相とを対応させたカラーバーＣＢが参照用として分割表示される。
【００５２】
これにより、画面上では、解析開始位置（または活動の開始位置）から解析終了位置（または活動の終了位置）まで電流源深さに応じて色相が変わって表示される。また、生体内電流源（磁場源）の位置変化に応じて電流源モデルの表示位置も徐々に移り変わる。通常、電流源の位置（深さ）は変化し且つその強度も変化するので、画面上の電流源モデルもその位置を変えながらかつ表示色相も経時的に変えられる。したがって、その開始位置から終了位置までの視認または目視による追跡も容易になり、その深さ位置の変化も色相から容易に把握できる。
【００５３】
さらに、この表示状態において、制御装置１７および解析装置１５は図６に示す処理を実行するので、第１の実施形態のときと同様に、カラーバーＣＢ内の指定色相（すなわち深さ位置）に呼応して、その色相に相当する電流原モデル円の種々の情報を表す輝度、コントラストなどを変化させ、その位置を容易に把握できる。
【００５４】
第４の実施形態
第４の実施形態に係る生体磁気計測装置を図１６〜１８に基づき説明する。
【００５５】
この装置は、電流源の方向と色相の関係を表示に取り入れている。これを行うため、解析装置１５は、前述した図２の処理を行い、そのステップＳ５における画像表示処理の詳細な一例として、図１６に記載の処理を行う。また解析装置１５と制御装置１７との共働による図１７に記載の処理も合わせて実行する。
【００５６】
図１６の画像表示処理によれば、解析装置１５は、電流源の方向と色相とを一義的に対応させた対応関係を表わすカラーサークルのデータを読み出す（ステップＳ５−４１）。このカラーサークルと計測された磁場波形とは一義的に対応している。このカラーサークルのデータは予め設定され、メモリに記憶している。次いで、このデータを参照して、ステップＳ３で解析された磁場源の方向の経時的データを色相データに変換する（ステップＳ５−４２）。この色相データは表示用データとして表示装置１６に送られる。
【００５７】
これにより、電流源の方向の時間変化データが、カラーサークルのデータを介して色相データに変換される。なお、方向を表すのに、その色相の輝度やコントラストを変えるように表示用データを作成してもよい。
【００５８】
このため、表示装置１６の画面には、図１８に示す如く、例えばＭＲＩ装置やＸ線ＣＴスキャナなどの他のモダリティで得られた頭部ＨＤの２次元断面像上にステップＳ５（ステップＳ５−４１，４２）で加工されたカラーサークルＣＣが重畳表示される。同時に、画面上の適宜な位置に、計測された磁場波形ＷＢが参照用として分割表示される。
【００５９】
この状態で、オペレータが入力装置１８を介して、磁場波形ＷＢ中の所望時間をカーソルで指定すると（図１７、ステップＳ５１，５２）、その指定を制御装置１７が検出して解析装置１５に伝える。このため、解析装置１５は、指定された時間を特定し（ステップＳ５３）、この特定した時間に対応した電流源の方向を決め、その方向に応じた色相を明滅させるデータを生成して、これを表示装置１６に送る（ステップＳ５４）。この一連の処理Ｓ５１〜Ｓ５４は、オペレータが入力装置１８に終了指令を与えるまで繰り返される（ステップＳ５５）。
【００６０】
このため、画面上では、オペレータが磁場波形ＷＢ中の所望時間を指定する度に、その時間に対応した方向が探索され、その方向の色相部分がカラーサークル上で例えば明滅（輝度又はコントラストの変更でも可）する。これにより、オペレータは自分が指定した時間に相当する電流源の方向を画面上で視覚的に容易に認識することができる。
【００６１】
つまり、波形の始まり位置および波形の終了位置を指定することで、解析開始位置（または活動の開始位置）および解析終了位置（または活動の終了位置）を電流原の方向の観点から容易に捕らえることができる。もちろん、時間波形中の診断的に興味ある任意の時間位置に対しても同様の把握ができる。
【００６２】
第５の実施形態
第５の実施形態を図１９に基づき説明する。この実施形態は、電流源の位置が不連続に変化するときの表示法に関する。
【００６３】
上述した実施形態およびその変形例では、解析結果である、電流源の強度の時間変化や、深さ、方向などの変化に応じて色相、コントラスト、輝度を変えて表示してきた。しかし、これにはさらに改善の余地がある。つまり、一連の解析結果の中で、電流源の位置が不連続に変化する場合である。この場合、上述してきた表示法によれば、電流源モデルの繋がりが複数個できるので、いずれの位置から解析が開始されてどの位置で解析が終了したかの判断にまごつくことにもなりかねない。
【００６４】
そこで、解析装置１５による画像表示処理において、ＳＱＵＩＤセンサからの出力信号を解析した結果、一連の解析結果において電流源の位置が空間的に不連続に変化していることが分かったとする。例えば、図１９に示す磁場波形について、電流源の位置がその前半と後半とでそれぞれ別個に局在し、不連続になっていたとする。
【００６５】
この場合、解析装置１５は、前半の位置の一群をグループ１として例えば赤を主原色とする色相で表し、例えば強度に応じてコントラストまたは輝度を変化させる一方で、後半の位置の一群をグループ２として別の色相、例えば緑を主原色とする色相で表し、例えば強度に応じてコントラストまたは輝度を変化させる電流源モデルの表示データを作成する。同時に、各グループにおいて、主原色の系統内で色相を経過時間に応じて変える。
【００６６】
この表示例を図１９に示す。同図にはカラーバーを図示していないが、グループ別に経過時間に対応したカラーバーを併せて表示する。これにより、電流源の位置が不連続に変化する場合でも、その変化自体を容易に認識でき、かつ、不連続な複数グループの電流源モデルにおける解析開始位置および解析終了位置を容易に観察できる。
【００６７】
この実施形態においても、前述した如く、磁場波形の時間変化グラフ上の所望時間位置をカーソルなどで指定すると、その指定時間に対応した電流源モデルＭＤの輝度、コントラストなどを変化させるようにしてもよい。
【００６８】
第６の実施形態
第６の実施形態を図２０に基づき説明する。この実施形態は、電流源の挙動の時間的変化の履歴を線、好適には矢印線で表示する方法に関する。
【００６９】
図２０には、この表示例を示す。解析装置１５により、ＳＱＵＩＤセンサからの磁場波形データが解析される。そこで、画像表示処理において、解析装置１５は、解析された磁場源の位置変化の履歴を一定幅の矢印線ＡＲ（図２０参照）で表す表示データを作成する。この表示データには、強度変化に応じて色相、コントラストおよび／または輝度を徐々に変化させる情報が組み込まれる。また、矢印線ＡＲの先頭部位が時間的に若い方を差すように作成される。
【００７０】
また、制御装置１７との共働によって、オペレータが画面を見ながら磁場波形の任意位置をカーソルなどで指定すると、その位置に対応した矢印線ＡＲ上の位置の色相、コントラストおよび／または輝度が変化するように構成される。
【００７１】
この実施形態によっても、電流源の一連の解析開始位置および解析終了位置を容易に且つ確実に把握できるとともに、変化の履歴の進む方向を矢印線ＡＲで表しているので、変化方向の読影作業の能率向上に寄与できる。
【００７２】
なお、この上述した変化履歴を表す矢印線の幅は一定限度内で任意であり、線そのものに近い矢印線から、かなり幅を持たせた面に近い矢印線まで適宜に設定可能である。
【００７３】
第７の実施形態
第７の実施形態を図２１に基づき説明する。この実施形態は電流源の表示モデルの形状に特徴を有し、その形状を表示の次元に応じて変えるものである。
【００７４】
解析装置１５は、画像表示処理を行って表示データを作成するが、その表示データに含まれる電流源モデルの形状を表示の次元に応じて変える処理を併せて行うようにする。例えば、指定されている表示法が３次元表示の場合には、図２１（ａ）〜（ｃ）に示すように、三角錐、多角錐、円錐の形状を成す電流源モデルを使用する。
【００７５】
一方、２次元表示の場合には、同図（ｄ）〜（ｉ）に例示する如く、３角形や、３角形と４角形を組み合わせた形、円に３角形を組み合わせた形などを使用する。いずれの形状もある一方向が突起状になっていることを特徴とする。
【００７６】
これらの形状データは、解析装置１５内のメモリに予め格納されており、どの形状を使用するかの入力装置１８からの選択信号に応じて読み出される。しかも、例えば、あるサンプリング時間における電流源の方向に、電流源モデルの突起の方向を合わせるように画像表示処理を行う。好適には、各サンプリング時間における電流源の方向に応じて色相を変えたり、その強度に応じて電流源モデル内の輝度を変えたりする態様が採られる。
【００７７】
これにより、表示画面には、各サンプリング時間における電流源の方向がそのモデルの突起の方向により逐一示されるから、電流源の方向の変化をより詳細に観察できるという利点がある。
【００７８】
変形例（その３）
上述の実施形態の変形例を図２２に基づき説明する。この変形例は、上述した第４の実施形態および第７の実施形態を組み合わせたものである。つまり、第４の実施形態で説明したカラーサークルＣＣと、第７の実施形態で説明した２次元表示用の３角形の電流源モデルＭＤとを分割表示するように構成したものである。この内、電流源モデルＭＤの方を頭部ＨＤの形態画像と重畳表示している。これにより、電流源の方向などの変化をより多角的に表示することができ、表示法の豊富化によって診断能の向上に寄与可能になる。
【００７９】
第８の実施形態
第８の実施形態を図２３に基づき説明する。この実施形態は、インデックス画像と個別画像と呼ばれる２種類の画像を同一画面上に表示することに特徴を有する。
【００８０】
以下に、電流源の強度を異なる色相で変調して表示する場合を説明するが、上述してきた種々の実施形態およびその変形例に基づき、その電流源強度を異なる輝度、またはコントラストで表す態様、電流源強度の時間変化を異なる色相、輝度、またはコントラストで表す態様、電流源の向き（方向）を異なる色相、輝度、またはコントラストで表す態様、さらには、電流源の深さを異なる色相、輝度、またはコントラストで表す態様などを実施してもよい。
【００８１】
解析装置１５は、その画像表示処理において、例えば電流源の強度と色相とを一義的に対応付けたテーブルを読み出す。次いで、一定の解析時間内における全サンプリング時間に存在する電流源による磁場分布に対し、全ての電流源をその強度に応じた色相の電流源モデルに変換し、それらの全モデルを形態画像上に重畳したインデックス画像データを作成する。これとともに、各サンプリング時間における個々の電流源をその強度に応じた色相で表した電流源モデルに変換し、それらを個別に形態画像上に重畳した個別画像データを用意する。
【００８２】
そして、インデックス画像および個別画像の表示データは、例えば図２３に示す如く、表示装置１６の同一画面上の分割態様で表示される。これにより、画面上部にはインデックス画像ＩＸが、その下側に複数個の個別画像ＩＶが同時に表示される。
【００８３】
オペレータが入力装置１８を介してインデックス画像ＩＸ上の興味ある１個又は複数個の電流源を例えばＲＯＩで指定すると、制御装置１７はこの指定内容を解析装置１５に伝える。このため、解析装置１５は、指定された電流源を個別に表した個別画像ＩＶを選択し、その画像の色相、輝度、コントラストなどを変化させて表示する。同時に、その個別画像を図２３に示すようにフレームＦＭで囲って明示してもよい。
【００８４】
反対に、個別画像ＩＶの中から任意の画像が指定されたときには、解析装置１５はインデックス画像ＩＸ内の対応する電流源の輝度、コントラスト、色相などを変化させて表示させる。
【００８５】
このようにインデックス画像ＩＸと個別画像ＩＶとの間で双方向若しくは一方向に指定し合うことができるので、電流源の解析時間内の全体的動きと個々のサンプリング時間における電流源の位置とを容易に確認できる。
【００８６】
なお、上述した実施形態およびその変形例にあっては、本発明の医用画像表示装置を生体磁気計測装置に適用した例を説明したが、この医用画像表示装置は生体から生じる磁気を用いて機能診断するファンクショナルＭＲＩに適用してもよい。
【００８７】
【発明の効果】
以上説明したように、本発明の医用画像表示装置は、収集データの経時的変化、または、解析期間内の各サンプリング時間に解析された解析データそれぞれを色相、輝度、またはコントラストの少なくとも１つを変化させて表示する表示処理手段を設けたので、とくに、生体磁気計測装置やファンクショナルＭＲＩ装置による収集データの表示処理に好適であり、磁場源（生体内電流源）の活動部位の動きの時間的変化、とくに、活動の始まる部位（解析を始めた部位）と活動の終わった部位（解析を終えた部位）を目視で良く観察できるとともに、時間的に変化する磁場波形の前後関係を簡単に理解できる。
【図面の簡単な説明】
【図１】本発明に係る生体磁気計測装置の全体構成の一例を示すブロック図。
【図２】この装置の解析装置で実施される解析処理および表示画像処理の概要を説明するための概略フローチャート。
【図３】第１の実施形態に係る画像表示処理を説明する概略フローチャート。
【図４】その画像表示処理を説明する説明図。
【図５】第１の実施形態に係る表示画面の一例を示す図。
【図６】制御装置および解析装置の共働により、色相指定に応じて対応する電流原モデルの表示状態を変化させる処理を示す概略フローチャート。
【図７】第１の実施形態の変形例に係る画像表示処理を説明するための概略フローチャート。
【図８】その変形例による表示画面を示す図。
【図９】第１の実施形態の別の変形例に係る表示画面のを示す図。
【図１０】第１の実施形態のさらに別の変形例に係る表示画面を示す図。
【図１１】第２の実施形態に係る画像表示処理を説明する概略フローチャート。
【図１２】その画像表示処理を説明する説明図。
【図１３】第２の実施形態に係る表示画面を示す図。
【図１４】第３の実施形態に係る画像表示処理を説明する概略フローチャート。
【図１５】第３の実施形態に係る表示画面を示す図。
【図１６】第４の実施形態に係る画像表示処理を説明する概略フローチャート。
【図１７】同様に第４の実施形態に係る画像表示処理を説明する概略フローチャート。
【図１８】第４の実施形態における表示画面の一例を示す図。
【図１９】第５の実施形態における表示画面の一例を示す図。
【図２０】第６の実施形態における表示画面の一例を示す図。
【図２１】第７の実施形態で用いられる電流源表示モデルの種々の形状を表す図。
【図２２】１つの変形例に係る表示画面を示す図。
【図２３】第８の実施形態で用いられる電流源表示モデルの種々の形状を表す図。
【図２４】従来技術の一例を説明する図。
【符号の説明】
１１  ＳＱＵＩＤ磁束センサ（磁気センサ）
１２  駆動回路
１３  プリプロセッサ
１４  Ａ／Ｄ変換器
１５  解析装置
１６  表示装置
１７  制御装置
１８  入力装置
【図１】
0004229489.tif 000003
【図２】
0004229489.tif 000004
【図３】
0004229489.tif 000005
【図４】
0004229489.tif 000006
【図５】
0004229489.tif 000007
【図６】
0004229489.tif 000008
【図７】
0004229489.tif 000009
【図８】
0004229489.tif 000010
【図９】
0004229489.tif 000011
【図１０】
0004229489.tif 000012
【図１１】
0004229489.tif 000013
【図１２】
0004229489.tif 000014
【図１３】
0004229489.tif 000015
【図１４】
0004229489.tif 000016
【図１５】
0004229489.tif 000017
【図１６】
0004229489.tif 000018
【図１７】
0004229489.tif 000019
【図１８】
0004229489.tif 000020
【図１９】
0004229489.tif 000021
【図２０】
0004229489.tif 000022
【図２１】
0004229489.tif 000023
【図２２】
0004229489.tif 000024
【図２３】
0004229489.tif 000025
【図２４】
0004229489.tif 000026