JP2014042725A - 被検体情報取得装置、表示方法、及びプログラム - Google Patents

Abstract

【課題】 適応型信号処理を用いて生成した画像のみを表示した場合、違和感が大きくなる可能性がある。
【解決手段】 本発明の被検体情報取得装置は、被検体内の各位置で反射した反射波を受信して複数の受信信号に変換する複数の変換素子と、前記複数の受信信号を用いて固定の重みで加算処理し、前記被検体内の各位置からの反射波に対応する信号を複数取得し、第１の分布情報を取得する固定型信号処理手段と、前記複数の受信信号を用いて、適応型信号処理を行い、第２の分布情報を取得する適応型信号処理手段と、前記第１の分布情報と、前記第２の分布情報と、が入力され、表示手段に画像情報を出力する表示制御手段と、を有し、
前記表示制御手段は、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記第１と第２の分布情報が合成された合成画像と、を同一画面内に並列表示させるための画像情報を出力する。
【選択図】 図１

Description

本発明は、被検体情報取得装置、表示方法、及びプログラムに関する。特に、被検体に弾性波を送信し、被検体からの反射波を受信することにより得られる分布情報を表示するための技術に関する。

被検体情報取得装置である超音波診断装置において、パルスエコー法によって画像データを形成する場合の深さ方向の空間分解能は、超音波の波長をλ、送信波数をｎとすると、（ｎλ）／２で一般的に表すことが可能である。例えば、１２ＭＨｚの中心周波数の超音波を２波長分送信した場合は約０．１３ｍｍ程度となる。

パルスエコー法について説明する。まず超音波パルス（弾性波）を被検体に送信すると、被検体内での音響インピーダンス差に応じて超音波が反射されて戻ってくる。次に、この反射波を受信し、反射波の受信信号を用いて画像データを生成する。代表的には、受信信号を用いて整相加算した後、包絡線を取得し、この包絡線を輝度値に変換して画像データを生成する。被検体内の複数の方向もしくは位置に対して超音波の送受信を繰り返すことで、超音波を送受信した方向の複数の走査線上の輝度情報を取得できる。この複数の走査線上の輝度情報を並べることで被検体内の画像化が可能となる。

なお超音波診断装置においては、超音波を電気信号に変換する複数の変換素子を用い、それぞれの素子間の受信信号波形に時間的なずれを加えることで、送信受信ともに被検体内でフォーカスさせるのが一般的である。

一方、レーダーの分野で発展してきた適応型信号処理という手法を超音波に組み合わせることで、空間解像度を向上させることができる。非特許文献１には、深さ方向に垂直な方向（走査線方向に垂直な方向）の空間解像度を上げるために、適応型信号処理であるＣａｐｏｎ法を用いる手法が記載されている。

また、深さ方向（走査線方向）の空間分解能を向上する技術として、非特許文献２では周波数領域干渉計法（ＦＤＩ法：Ｆｒｅｑｕｅｎｃｙ Ｄｏｍａｉｎ Ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ）と、適応型信号処理であるＣａｐｏｎ法と、を適用し、血管壁の層構造を画像化した結果が示されている。受信信号にＦＤＩ法及びＣａｐｏｎ法を適用することで、深さ方向の空間分解能を向上させることができる。ただし、ＦＤＩの処理を行うために切り出した深さ方向の信号の範囲（処理レンジ内）には、複数の反射層が存在することが想定される。また、近接した反射層からの複数の反射波は、互いに高い相関性を有している可能性が高い。このような高い相関性を有する複数の反射波の受信信号に対してＣａｐｏｎ法などの適応型信号処理をそのまま適用すると、所望の信号を打ち消すなどの予期しない動作を行うことが知られている。このような相関性を有する信号（相関性干渉波）による影響を低減（抑圧）するため、周波数平均法（ｆｒｅｑｕｅｎｃｙ ａｖｅｒａｇｉｎｇ ｔｅｃｈｎｉｑｕｅ）を合わせて用いることで、反射波の受信信号に対してＦＤＩ法とＣａｐｏｎ法が適用可能となる。

Ｐｒｏｃ．Ａｃｏｕｓｔｉｃｓ，Ｓｐｅｅｃｈ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓ．ｐｐ．４８９−４９２（Ｍａｒ．２００５） Ｈｉｒｏｆｕｍｉ Ｔａｋｉ，Ｋｏｕｓｕｋｅ Ｔａｋｉ，Ｔａｋｕｙａ Ｓａｋａｍｏｔｏ，Ｍａｋｏｔｏ Ｙａｍａｋａｗａ，Ｔｓｕｙｏｓｈｉ Ｓｈｉｉｎａ ａｎｄ Ｔｏｒｕ Ｓａｔｏ：Ｃｏｎｆ Ｐｒｏｃ ＩＥＥＥ Ｅｎｇ Ｍｅｄ Ｂｉｏｌ Ｓｏｃ．２０１０；１：５２９８−５３０１．

Ｃａｐｏｎ法や、ＦＤＩ法及びＣａｐｏｎ法を組み合わせる手法等の、適応型信号処理を適用することで、画像の空間分解能を向上させることができる。しかしながら、このような新規な技術で生成された画像を表示すると、ユーザー（特に医師）は、従来のＢモード画像（複数の受信信号を整相加算した後、包絡線を取得し、この包絡線を輝度値に変換した画像）を見慣れているため、違和感を生じる可能性がある。特に、適応型信号処理を用いて生成した画像のみを表示した場合、その違和感が大きくなる可能性がある。

そこで本発明は、上記課題に鑑み、適応型信号処理を用いて生成した画像を表示する場合において、使い勝手の良い表示方法及び被検体情報取得装置を提供することを目的とする。

本発明の被検体情報取得装置は、被検体に弾性波を送信し、被検体内の各位置で反射した反射波を受信して複数の受信信号に変換する複数の変換素子と、
前記複数の受信信号を用いて固定の重みで加算処理し、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得し、前記複数の走査線信号を用いて第１の分布情報を取得する固定型信号処理手段と、
前記複数の受信信号を用いて、前記受信信号に応じて適応的に変化する重みで適応型信号処理を行い、第２の分布情報を取得する適応型信号処理手段と、
前記第１の分布情報と、前記第２の分布情報と、が入力され、表示手段に画像情報を出力する表示制御手段と、
を有し、
前記表示制御手段は、
前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記第１と第２の分布情報が合成された合成画像と、を同一画面内に並列表示させるための画像情報を出力することを特徴とする。

また、本発明の表示方法は、被検体情報取得装置で取得された分布情報を用いて画像を表示する表示方法であって、
前記取得された分布情報は、
被検体に弾性波を送信し被検体内で反射した反射波を受信することにより得られる受信信号を用いて固定の重みで加算処理し、前記被検体内の各位置からの反射波に対応する信号である走査線信号を複数取得し、複数の前記走査線信号を用いることで得られる第１の分布情報と、
前記複数の受信信号を用いて前記受信信号に応じて適応的に変化する重みで適応型信号処理を行うことにより得られる第２の分布情報と、
を含み、
前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記第１と第２の分布情報が合成された合成画像と、を同一画面内に並列表示するステップを有することを特徴とする。

本発明によれば、適応型信号処理を用いて生成した画像を表示する場合において、使い勝手の良い表示方法及び被検体情報取得装置を提供することができる。

本発明の適用可能な被検体情報取得装置のシステム概要を示す模式図である。 固定型信号処理ブロックの構成を示す模式図である。 適応型信号処理ブロックの構成を示す模式図である。 第１の実施形態の表示方法のフローを示すフローチャートである。 第１の実施形態により表示部に表示される画面の一例を示す図である。 第２の実施形態により表示部に表示される画面の一例を示す図である。 第３の実施形態の表示方法のフローを示すフローチャートである。 第３の実施形態により表示部に表示される画面の一例を示す図である。 第４の実施形態において拡大率と合成率の関係を示す図である。

以下、本発明の実施の形態について、図面を用いて説明する。同一の構成要素には原則として同一の符号を付して、説明を省略する。

なお、本発明において、弾性波とは、典型的には超音波であり、音波、超音波、音響波、と呼ばれる弾性波を含む。本発明の被検体情報取得装置とは、被検体に弾性波を送信し、被検体内部で反射した反射波（反射した弾性波）を受信して、被検体内の分布情報を画像データとして取得する装置を含む。取得される被検体内の分布情報とは、被検体内部の組織の音響インピーダンスの違いを反映した情報である。また、本発明において走査線とは、探触子から送信される弾性波の進行方向に形成される仮想的な線を示す。

＜第１の実施形態＞
本実施形態では、本発明の基本的な装置構成及び処理フローについて説明する。

（被検体情報取得装置の基本的な構成）
本発明の被検体情報取得装置の構成について図１を用いて説明する。図１は本発明が適用できる被検体情報取得装置のシステム概要を示す模式図である。本実施形態の被検体情報取得装置は、複数の変換素子００２を有する探触子００１、受信回路系００５、送信回路系００３、固定型信号処理ブロック００６、適応型信号処理ブロック００７、表示制御手段００８を備える。さらに本発明の被検体情報取得装置は、表示手段００９、入力手段０１０、システム制御部００４を備える。

探触子００１は弾性波を被検体内の複数位置に送信して、反射波を受信する送受信器であり、弾性波を電気信号に変換する変換素子００２を複数有する。

送信回路系００３は、システム制御部００４からの制御信号に従って、注目位置や注目方向に応じた遅延時間や振幅を有する複数の送信信号を生成する送信信号生成手段である。この送信信号は複数の変換素子００２によって弾性波に変換され、弾性波ビームとして探触子００１から被検体へと送信される。被検体内部の対象物（反射界面や反射体）で反射された弾性波（反射波）は、複数の変換素子００２によって受信され複数の受信信号に変換される。受信信号は受信回路系００５に入力される。

受信回路系００５は、複数の受信信号を増幅し、複数のデジタル信号（デジタル化された受信信号）に変換する受信信号処理手段である。ここで、本発明では、変換素子００２が出力したアナログの受信信号だけでなく、増幅やデジタル変換等の処理を行った信号も受信信号と表現する。受信回路系００５から出力された複数のデジタル信号は、固定型信号処理ブロック００６及び適応型信号処理ブロック００７に入力される。

固定型信号処理ブロック００６は、本発明における固定型信号処理手段に相当する。図２は、固定型信号処理ブロック００６内の構成を示している。固定型信号処理ブロック００６内では、まず、整相加算手段である遅延加算回路０１１が、弾性波を送信した方向や位置に応じて、複数のデジタル信号に対する遅延処理を行い、さらに遅延処理後の複数のデジタル信号に対して加算処理を行う。つまり整相加算処理を実行する。この整相加算処理により、走査線信号が複数取得される。なお、この加算処理の際に複数のデジタル信号にそれぞれ重みを乗じた上で加算処理することも可能である。この重みは観察位置や送受信の条件によって変化するが、あらかじめ決められた（固定された）重みを用いることが多い。整相加算処理は、被検体内の各位置で反射した反射波の音圧に対応する信号を走査線信号として生成する処理である。次に、包絡線検波手段である包絡線検波回路０１２が、複数の走査線信号を用いて包絡線検波を行い、第１の分布情報を取得する。固定型信号処理ブロック００６で得られた第１の分布情報は、表示制御手段００８に入力される。

適応型信号処理ブロック００７は、本発明のおける適応型信号処理手段に相当する。適応型信号処理は、受信信号に応じて、その処理パラメータを適応的に変化させる処理である。特に、適応型信号処理の一つであるＣａｐｏｎ法（拘束付電力最小化規範（ＣＭＰ：Ｃｏｎｓｔｒａｉｎｅｄ Ｍｉｎｉｍｉｚａｔｉｏｎ ｏｆ Ｐｏｗｅｒ）ともいう）は、複数の入力信号に対して、注目方向や注目位置に関する感度を固定した状態で電力を最小化するように処理する方法である。このような適応型信号処理は空間解像度を向上させる効果がある。適応型信号処理ブロック００７は、深さ方向と、深さ方向に垂直な方向と、のうち少なくとも一つの方向の解像度を向上させた電力強度分布を第２の分布情報として出力する。ここで、深さ方向とは、探触子から送信される弾性波（超音波ビーム）の進行方向であり、走査線方向である。適応型信号処理の詳細は、図３を用いて後述する。

表示制御手段００８は、固定型信号処理ブロック００６から第１の分布情報が入力され、適応型信号処理ブロック００７から第２の分布情報が入力される。そして、表示制御手段００８は表示手段００９に表示用の画像情報を出力する。この表示制御手段００８から出力される画像情報を基に被検体内の分布情報を示す画像が表示手段００９に表示される。表示制御手段００８が行う詳細な処理内容は図４を用いて後述する。なお、表示制御手段００８は、第１の分布情報の画像情報、第２の分布情報の画像情報、第１と第２の分布情報の合成画像情報に対して、エッジ強調やコントラスト調整などの各種画像処理を行い、輝度データの画像情報を出力する。

表示手段００９は、表示制御手段００８から入力された画像情報をもとに画像を表示する。ＬＣＤ（Ｌｉｑｕｉｄ Ｃｒｙｓｔａｌ Ｄｉｓｐｌａｙ）やＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）、有機ＥＬディスプレイ等で構成される。

入力手段０１０は、ユーザーが指定領域の入力を行うための手段である。ユーザーは表示手段００９に表示された第１の分布情報の画像を参照しながら、入力手段０１０を用いて所定の領域を指定する。入力手段０１０としては、マウスやキーボード等のポインティングデバイスや、ペンタブレットタイプ、表示手段００９表面に取り付けたタッチパッド等で構成される。なお、表示手段００９、入力手段０１０は、本発明の被検体情報処理装置が有する構成とはせずに、別に用意してから夫々を本発明の被検体情報取得装置に接続しても良い。

（適応型信号処理の詳細）
ここで、本発明の適応型信号処理ブロック００７が行う処理について説明する。図３は、適応型信号処理ブロック００７内の構成を示している。図３（ａ），（ｂ），（ｃ）は、それぞれ適応型信号処理ブロック００７内の構成が異なる。以下、図３（ａ），（ｂ），（ｃ）を用いて、それぞれ本発明の適用可能な適応型信号処理ブロック００７内の構成の例を説明する。

まず、図３（ａ）は、深さ方向（探触子から送信される弾性波（超音波ビーム）の進行方向）に垂直な方向の解像度を上げるための適応型信号処理ブロック００７の構成を示す。Ｐｒｏｃ．Ａｃｏｕｓｔｉｃｓ，Ｓｐｅｅｃｈ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓ．ｐｐ．４８９−４９２（Ｍａｒ．２００５）には、このような深さ方向に垂直な方向の解像度を上げるための適応型信号処理の手法が記載されている。

ここで、複数の受信信号に対して適応型信号処理を適用した場合の処理についてＣａｐｏｎ法を例にして述べる。

まず複数の受信信号から相関行列を算出するところまでを説明する。最初に、遅延処理回路２０１が、複数の変換素子から出力された複数の受信信号に対してヒルベルト変換、注目位置に応じた遅延処理（整相処理）を実施する。このようにして算出されるのが、複素表現された受信信号である。ここでｋ番目の素子からの受信信号を処理して得られた信号のｓサンプル目をｘｋ［ｓ］として、ｓサンプル目の入力ベクトルＸ［ｓ］を以下のように定義する。なお、ここでＭは素子数である。

次にこの入力ベクトルＸ［ｓ］を用いて適応型信号処理部であるＣａｐｏｎ回路２０２が、相関行列Ｒｘｘを算出する。

式中の右肩のＨは複素共役転置を表し、右肩の＊は複素共役を表す。Ｅ［・］は時間平均を算出する処理であり、サンプルの番号（ここではｓ）を変化させ、その平均を算出することを意味する。

次に、Ｃａｐｏｎ回路２０２は、注目方向以外から探触子に到達する相関性干渉波による影響を抑圧するために、相関行列Ｒｘｘに対して空間平均法を適用し、平均相関行列Ｒ’ｘｘを求める。

ここでＲ^ｎｘｘは相関行列Ｒｘｘの中の部分行列を表しており、Ｒｘｘの対角成分上を移動し、Ｒｘｘの（ｎ、ｎ）成分をその１番目の対角成分とする位置にあるＫ×Ｋのサイズの行列である。Ｚｎはそれぞれの部分行列を加算する際の係数であり、Ｚｎの総和が１になるように調整される。

Ｃａｐｏｎ法では、ある拘束条件下で出力電力を最小化するための複素ウェイトを求める。複素ウェイトは、複素ベクトルで表現されるウェイト（重み）である。注目方向からの弾性波の受信信号に対する感度を１に拘束した状態で、出力電力を最小化するための最適な複素ウェイトＷｏｐｔはＣａｐｏｎ法においては以下の式で求められる。

Ｃは拘束ベクトルであり、素子の位置と注目方向に応じて変化するものである。ただし、受信信号に対して整相遅延処理を実施している場合は、平均相関行列のサイズ（この場合はＫ）において、すべての値が１であるベクトルとして構わない。

また、複素ウェイトＷｏｐｔを用いて、算出された電力Ｐｍｉｎは以下のように求まる。この算出された電力Ｐｍｉｎが、本例における、被検体内部の組織の音響インピーダンスの違いを反映した分布（音響特性に関する分布）情報を示している。

このように、Ｃａｐｏｎ回路２０２は受信信号から相関行列、さらには平均相関行列を求め、その逆行列を用いて複素ウェイトや複素ウェイトを用いた場合の電力分布を取得できる。この複素ウェイトや複素ウェイトを用いた場合の電力は、注目方向からの弾性波の信号に対して感度を１にし、それ以外の方向から到達する弾性波の信号を抑圧した場合のウェイトや電力である。つまり、Ｃａｐｏｎ法では注目方向からの弾性波の信号を選択的に抽出することが可能で、その結果として深さ方向に直交する方向の空間分解能を向上することが出来る。なお、逆行列を直接求めずに、平均相関行列に対するＱＲ分解と後退代入処理によっても、電力は算出可能である。

このように、適応型信号処理ブロック００７は、複数の受信信号を用いて、受信信号に応じて適応的に変化する重みで適応型信号処理（ここではＣａｐｏｎ法を用いた）を行う。その結果、深さ方向に垂直な方向の解像度が向上した電力強度分布（第２の分布情報に相当する）を出力する。

次に、図３（ｂ）を用いて、適応型信号処理ブロック００７の構成の２つ目の例を示す。
図３（ｂ）は、深さ方向（探触子から送信される弾性波（超音波ビーム）の進行方向）の解像度を上げるための適応型信号処理ブロック００７の構成を示す。このような深さ方向の空間分解能を向上するため、適応型信号処理に、周波数領域干渉計法（ＦＤＩ法：Ｆｒｅｑｕｅｎｃｙ Ｄｏｍａｉｎ Ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ）を組み合わせる手法がある。Ｈｉｒｏｆｕｍｉ Ｔａｋｉ，Ｋｏｕｓｕｋｅ Ｔａｋｉ，Ｔａｋｕｙａ Ｓａｋａｍｏｔｏ，Ｍａｋｏｔｏ Ｙａｍａｋａｗａ，Ｔｓｕｙｏｓｈｉ Ｓｈｉｉｎａ ａｎｄ Ｔｏｒｕ Ｓａｔｏ：Ｃｏｎｆ Ｐｒｏｃ ＩＥＥＥ Ｅｎｇ Ｍｅｄ Ｂｉｏｌ Ｓｏｃ．２０１０；１：５２９８−５３０１．に、周波数領域干渉計法（ＦＤＩ法）と、適応型信号処理であるＣａｐｏｎ法と、を適用した手法について記載されている。

ＦＤＩ法は、受信信号を周波数ごとに分解し、分解された信号の位相を注目位置に応じて変化させることで、注目位置における受信電力を推定する方法である。なお、位相の変化量はある基準位置から注目位置までの距離と周波数に対応した波数の積とからあらかじめ決定できる。

つまり、ＦＤＩ法と適応型信号処理とを組み合わせることは、各周波数成分に分解された受信信号に対して、あらかじめ決定された固定の位相変化量・重みではなく、適応型信号処理によって信号に応じて算出された位相変化量・重みを用いて、注目位置における受信電力を推定することになる。

さらに、パルス波のように広い周波数帯域を有する弾性波の受信信号に対して周波数平均法を用いる際には、参照信号によって受信信号のホワイトニングを行うとよい。

図３（ｂ）では、まず、整相加算手段である遅延加算回路３０１が、弾性波を送信した方向や位置に応じて、複数のデジタル信号に対する遅延処理を行い、さらに遅延処理後の複数のデジタル信号の加算処理を行う。つまり整相加算処理を実行する。この整相加算処理は、固定型信号処理ブロック００６内の整相加算と同じであり、被検体内の各位置で反射した反射波の音圧に対応する信号を走査線信号として生成する処理である。

次に、ＦＤＩ適応処理部であるＦＤＩ・Ｃａｐｏｎ回路３０２が、遅延加算回路３０１から出力された複数の走査線信号を入力信号として受け取る。そして、この複数の走査線信号から、１回で処理する時間分、つまり処理レンジ分の信号を抽出する。

次に、ＦＤＩ・Ｃａｐｏｎ回路３０２は、抽出した信号をフーリエ変換し周波数ごとの成分（Ｘｓ１、Ｘｓ２、Ｘｓ３、・・・、ＸｓＮ）に分割する。一方で、ＦＤＩ・Ｃａｐｏｎ回路３０２には、不図示の参照信号記憶部から少なくとも１つの参照信号が入力される。そして、ＦＤＩ・Ｃａｐｏｎ回路３０２は、参照信号のフーリエ変換を行い、参照信号を周波数ごとの成分（Ｘｒ１、Ｘｒ２、Ｘｒ３、・・・、ＸｒＮ）に分割する。

次に、ＦＤＩ・Ｃａｐｏｎ回路３０２は、下記式に示すホワイトニング処理を行う。

ここでＸｗｋ（ｋ＝１，２，・・・，Ｎ）はホワイトニング処理後の周波数ごとの成分、ηは安定化のための微小量、＊は複素共役を意味する。ホワイトニング処理された各周波数成分からなるベクトルＸｆを用いて、相関行列Ｒを算出する。
Ｘｆ＝［Ｘ_Ｗ１，Ｘ_Ｗ２，・・・，Ｘ_ＷＮ］^Ｔ
Ｒ＝ＸｆＸｆ^Ｔ＊
なおＴは転置を意味する。ここで相関行列ＲはＮ×Ｎのサイズを有する行列となる。次に相関行列Ｒから部分行列を抽出し、それらを平均化する周波数平均法を適用する。

Ｒ’は周波数平均相関行列、ＲｍはＲｍｉｊを要素に持つ相関行列Ｒの部分行列である。このようにして周波数平均相関行列Ｒ’が算出される。

次に、拘束ベクトルＣが、ＦＤＩ・Ｃａｐｏｎ回路３０２に入力される。拘束ベクトルＣは、処理レンジ内での位置ｒに応じて変化するベクトルであり、以下の式で定義される。
Ｃ＝［ｅｘｐ（ｊｋ_１ｒ），ｅｘｐ（ｊｋ_２ｒ），・・・，ｅｘｐ（ｊｋ_{（Ｎ−Ｍ＋１）}ｒ）］
これらの周波数平均相関行列Ｒ’ならびに拘束ベクトルＣを用いて、処理レンジ内の電力強度分布Ｐ（ｒ）を算出する。この算出された電力Ｐ（ｒ）が、本例における、被検体内部の組織の音響インピーダンスの違いを反映した分布（音響特性に関する分布）情報を示している。

η’Ｅは逆行列算出を安定させるために加算した対角行列である。

このように、本例では、適応型信号処理ブロック００７は、複数の受信信号を用いてＦＤＩ法及び適応型信号処理（ここではＣａｐｏｎ法を用いた）を行う。その結果、深さ方向の解像度が向上した電力強度分布（第２の分布情報に相当する）を出力する。

次に、図３（ｃ）を用いて、適応型信号処理ブロック００７の構成の３つ目の例を示す。遅延処理回路４０１は、複数の変換素子から出力された複数の受信信号に対してヒルベルト変換、注目位置に応じた遅延処理を実施しデジタル信号を出力する。Ｃａｐｏｎ回路４０２は、遅延処理されたデジタル信号を入力としＣａｐｏｎ処理を行う。先ほどと同様の処理は省略するが、最終的にＣａｐｏｎ回路４０２は以下の式で算出される信号Ｙ［ｓ］を出力する。ここでＸ‘［ｓ］はｓサンプル目の入力ベクトルＸ［ｓ］からＷｏｐｔのサイズに合わせたベクトルを抽出したものである。

この出力Ｙ［ｓ］は注目位置に応じた反射波形の位相情報まで保持しているため、この後のＦＤＩ・ＣＡＰＯＮ処理を行うことが可能となる。ＦＤＩ・ＣＡＰＯＮ回路３０２では入力された信号Ｙ［ｓ］に対してＦＤＩ・ＣＡＰＯＮ処理を行い電力強度分布を出力する。

このような処理を行うことで、深さ方向に直交する方向ならびに深さ方向の解像度が向上した電力強度分布を得ることが出来る。

なおここでは適応型信号処理の例としてＣａｐｏｎ法の処理について述べたが、他の適応型信号処理である、ＭＵＳＩＣ法やＥＳＰＲＩＴ法などにおいても本発明の効果が同様に得られる。

（表示方法）
以下、本実施形態の表示方法のフローについて図４を用いて説明する。図４は、本実施形態の表示方法を説明するためのフローチャートである。

まず、表示制御手段００８は、Ｓ１０１でユーザーからの、被検体内の分布情報を示す画像の並列表示の指示があるかどうか（並列表示の指示情報が入力されているかどうか）を確認する。例えば、表示手段００９内の画面に表示される並列表示ボタン（図５（ａ）（ｂ）参照）がユーザーによりクリックされることで、並列表示の指示が入力されるようにするとよい。

Ｓ１０１でＹＥＳの場合、Ｓ１０２のステップにおいて表示制御手段００８は、第１の分布情報の画像と、第２の分布情報の画像又は第１と第２の分布情報が合成された合成画像と、を同一画面内に並列に表示するための画像情報を出力する。この画像情報により表示手段００９には画像が並列表示される。本実施形態では、第１の分布情報の画像と、第２の分布情報の画像と、を並列表示する場合について説明する。第１の分布情報と第２の分布情報との合成画像を表示する例は後述の第４の実施形態で説明する。

図５に、本実施形態により表示手段００９に表示される画面の例を示す。図５は、血管壁の層構造を示している。図５の例では、第２の分布情報の画像を得るため、適応型信号処理として、ＦＤＩ法とＣａｐｏｎ法を組み合わせた処理（図３（ｂ）の例）を行った。

図５の例では、画面上に並列表示ボタンが表示されている。並列表示ボタンがクリックされ、並列表示モードがオンになると、第１の分布情報の画像と第２の分布情報の画像とを並列表示するための指示が、表示制御手段００８に入力される。並列表示ボタンは、表示制御手段００８に第１の分布情報及び第２の分布情報が入力されている状態で、ユーザーがクリックすることにより並列表示モードをオンにすることが可能である。また、第１の分布情報及び第２の分布情報が入力される前に並列表示ボタンをクリックする構成としてもよい。

図５に示されているように、右図の第２の分布情報の画像（ＦＤＩ＋Ｃａｐｏｎ）のほうが、左図の第１の分布情報の画像（通常Ｂモード）よりも、深さ方向において高解像度となっていることが分かる。具体的には、血管壁がシャープになっていることが分かる。

本実施形態のように、第２の分布情報の画像だけでなく、第１の分布情報の画像も並列表示することで、ユーザーの違和感を低減することができる。

また、本実施形態では、第１の分布情報の画像内の位置と第２の分布情報の画像内の位置とを対応付けるためのガイドを表示している。つまり、第１の分布情報の画像内と第２の分布情報の画像内とで、被検体内の同じ位置を示すガイドを表示している。具体的に、図５（ａ）では、第１の分布情報の画像内の所定位置と、第２の分布情報の画像内の前記所定位置に対応する位置と、を対応づけるためのガイドとして十字のマークを表示し、被検体内の同じ位置を示している。また、図５（ｂ）では、点線のガイドにより各画像内における被検体の同じ位置を示している。

ただし、本実施形態では、位置を対応づけるためのガイドがない方が、画像の視認性を高めることができる場合があるため、位置を対応づけるためのガイドを表示するモードと表示しないモードとをユーザーが選択できるようにするとよい。例えば、表示画面上に、ガイド表示をオン、オフするためのボタンを設けるとよい。

また、表示制御手段００８は、上述の画像を並列表示する並列表示モードに加え、単独表示モードと重畳表示モードのうちいずれか一方のモードをさらに有するとよい。つまり、並列表示モードを、単独表示モードや重畳表示モードに切り替えることが可能であり、最初に表示する表示モードをいずれかから選択的に実行可能である。単独表示モードは、第１の分布情報の画像、第２の分布情報の画像、第１及び第２の分布情報の合成画像、のいずれかの画像を単独で表示するモードである。重畳表示モードは、第１の分布情報の画像と、第２の分布情報の画像（又は第１及び第２の分布情報の合成画像）と、を同じ表示領域に重ねて表示するモードである。このように、並列表示モード以外の異なる表示モードも有することで、ユーザーの使い勝手がよくなる。モード切替え等の場合は、例えば、表示手段００９の画面上にモード切替えボタンを表示し、このモード切替えボタンをクリックすることにより、ユーザーは、どの表示モードで表示するかを選択できる。表示制御手段００８は、ユーザーからのモード選択の入力を受け、選択された表示モードの画像情報を出力する。

また、上述の例では、第１の分布情報と並列表示する画像として、第２の分布情報の画像を表示したが、第１の分布情報と第２の分布情報との合成画像を、第１の分布情報の画像と並列表示しても本発明の効果は得られる。

＜第２の実施形態＞
本実施形態は、表示手段００９に表示される画像の並べ方が、適応型信号処理に応じて決められていることが第１の実施形態と異なる。本実施形態の被検体情報取得装置は、図１で示した装置と同様の構成の装置を用いる。また、表示方法の概略は図４で説明したフローと基本的に同じであるため、以下では、図６を用いて第１の実施形態とは異なる表示部分に絞って説明する。

本実施形態では、第２の分布情報を生成する適応型信号処理が、深さ方向の解像度を上げるための処理か、深さ方向に直交する方向の解像度を上げるための処理か、によって、並列表示する画像の並べ方を変える。つまり、用いた適応型信号処理が、複数の受信信号を用いて、注目方向に関する感度を固定した状態で電力を最小化するように処理する方法（Ｃａｐｏｎ法）か、深さ方向の注目位置に関する感度を固定した状態で電力を最小化するように処理する方法（Ｃａｐｏｎ法とＦＤＩ法とを組み合わせた方法）か、により並列表示の仕方を変える。

図６は、本実施形態により表示手段００９に表示される画面の一例である。図６では、深さ方向に直交する方向の解像度を向上するために、適応型信号処理としてＣａｐｏｎ法を用いて得られた第２の分布情報の画像を、第１の分布情報の画像の下側に表示している。つまり、第１の分布情報の画像と第２の分布情報の画像とを、ユーザーから見て画面上の縦方向に並ぶように並列表示している。ここで、本実施形態における縦方向とは、第２の分布情報の画像内（もしくは合成画像内）における深さ方向と同じ方向である。これは、用いている適応型信号処理が、深さ方向に直交する方向の解像度を向上するＣａｐｏｎ法であるため、縦方向に並ぶよう表示したほうが、ユーザーの視認性があがるためである。例えば、第１の分布情報と第２の分布情報との同じ位置の画像を比較したい場合、深さ方向に直交する方向の解像度が高い場合は、２つの分布情報を横に並べるよりも、縦に並べた方が深さ方向に直交する方向の位置の比較が行いやすい場合が多い。

また、適応型信号処理としてＣａｐｏｎ法とＦＤＩ法とを組み合わせた方法を用いた場合は、得られる第２の分布情報の画像は、図５（ａ）（ｂ）のように、第１の分布情報の画像の横側に表示するとよい。本実施形態における横方向とは、第２の分布情報の画像内（もしくは合成画像内）における深さ方向に交差（典型的には直交）する方向である。つまり、第１の分布情報の画像と第２の分布情報の画像とを、ユーザーから見て画面上の横方向に並ぶように並列表示するとよい。これは、用いている適応型信号処理が、深さ方向の解像度を向上するＣａｐｏｎ法とＦＤＩ法とを組み合わせた方法であるため、横方向に並ぶよう表示したほうが、ユーザーの視認性があがるためである。例えば、第１の分布情報と第２の分布情報との同じ位置の画像を比較したい場合、深さ方向の解像度が高い場合は、２つの分布情報を縦に並べるよりも、横に並べた方が深さ方向の位置の比較が行いやすい場合が多い。

このように、適応型信号処理の処理内容に応じて並列表示の並べ方を変えることで、ユーザーの視認性がより向上する。

また、本実施形態においても、図６に示すように、第１の分布情報の画像内の位置と第２の分布情報の画像内の位置とを対応付けるためのガイドとして点線を表示している。このガイドは第１の実施形態と同様に、位置を対応づけるためのガイドを表示するモードと表示しないモードとをユーザーが選択できるようにするとよい。

さらに、本実施形態においても、表示制御手段００８は、上述の画像を並列表示する並列表示モードに加え、単独表示モードと重畳表示モードのうちいずれか一方のモードをさらに有するとよい。

また、第１の分布情報と並列表示する画像として、第１の分布情報と第２の分布情報との合成画像を、第１の分布情報の画像と並列表示しても本発明の効果は得られる。

＜第３の実施形態＞
本実施形態は、表示制御手段００８が、まずは第１の分布情報の画像を表示する。そして、第１の分布情報の画像内における指定領域の情報を受けた場合に、指定領域の位置の第１の分布情報の画像と、指定領域に対応する位置の第２の分布情報の画像又は第１及び第２の分布情報の合成画像と、を拡大して並列表示することを特徴とする。それ以外は第１の実施形態及び第２の実施形態と同じである。本実施形態の被検体情報取得装置は、図１で示した装置と同様の構成の装置を用いる。

本実施形態の表示方法のフローについて図７を用いて説明する。図７は、本実施形態の表示方法を説明するためのフローチャートである。本実施形態では、第１の分布情報の画像と、第２の分布情報の画像と、を拡大して並列表示する場合について説明する。

まず、表示制御手段００８は、Ｓ３０１でユーザーからの、並列表示の指示があるかどうか（並列表示の指示情報が入力されているかどうか）を確認する。Ｓ３０１でＹＥＳの場合、実施形態１と同様に、Ｓ３０４のステップにおいて、第１の分布情報の画像と、第２の分布情報の画像と、を同一画面内に並列に表示するための画像情報を出力する。

Ｓ３０１でＮＯの場合、Ｓ３０２のステップにおいて、第１の分布情報の画像を単独で表示するための画像情報を表示手段００９に出力する。この画像情報により表示手段００９には第１の分布情報の画像が表示される。

次に、Ｓ３０３でユーザーからの指定領域の指示があるかどうか（指定領域の情報が入力されたかどうか）を確認する。ユーザーは、表示手段００９に表示された第１の分布情報の画像を見ながら、マウス等の入力手段０１０を用いて拡大したい領域として指定領域を入力する。システム制御部００４には入力手段０１０から指定領域の情報が入力され、システム制御部００４は表示制御手段００８に、指定領域の情報をユーザーからの拡大指示に関する情報として出力する。このように、拡大領域を指定する場合は、ユーザーが第１の分布情報の画像内における所定の領域を指定することによって拡大指示となる。そして、表示制御手段００８は、ユーザーにより指定された領域の大きさと表示手段００９の表示領域の大きさとの関係で拡大画像の拡大倍率を決定する。また、拡大領域が指定された後、表示手段００９内の画面に表示される拡大ボタン（図８参照）がユーザーによりクリックされることで、拡大開始の指示が入力されるようにしてもよい。

Ｓ３０３でＹＥＳの場合、Ｓ３０４のステップで、指定領域の位置の第１の分布情報の拡大画像と、指定領域に対応する位置の第２の分布情報の拡大画像とを並列表示する。図８は、本実施形態により表示手段００９に表示される画面の一例である。

図８は、血管壁の層構造を示している。画面内において、上図は拡大前の画像（第１の分布情報の画像）であり、左下図は指定領域が拡大された第１の分布情報の画像（通常Ｂモード）、右下図は指定領域が拡大された第２の分布情報の画像（ＦＤＩ＋Ｃａｐｏｎ）である。図８の例では、第２の分布情報の画像を得るため、適応型信号処理として、ＦＤＩ法とＣａｐｏｎ法を組み合わせた処理（図３（ｂ）の例）を行った。このように、第１の分布情報の画像と、第２の分布情報の画像と、を拡大して並列表示することによりさらに使い勝手の良い表示とすることができる。これは、拡大率によっては、単純に第１の分布情報の画像を拡大しても解像度が悪く視認性が高まらない場合がある。この場合、適応型信号処理を行って取得した解像度の良い第２の分布情報の拡大画像を表示することにより視認性があがるだけでなく、第１の分布情報の拡大画像と並列表示することにより両者を比較しやすくなる。

また、図８では、画面上には上述した拡大ボタンが表示されている。ユーザーが拡大領域を指定完了している状態（表示制御手段００８に拡大指示に関する情報が入力されている状態）で、拡大ボタンをクリックすると、拡大開始の指示が表示制御手段００８に入力される。表示制御手段００８は、入力された拡大開始指示により、第１の分布情報の拡大画像と第２の分布情報の拡大画像とを並列表示させるための画像情報を、表示手段００９に出力する。この画像情報により、表示手段００９の画面が切り替わるようにするとよい。

さらに、本実施形態では、拡大画像が表示された画面内の別の表示領域（図８では上側）に、拡大前の第１の分布情報の画像もサムネイル表示している。そして、拡大画像が示している、第１の分布情報の画像上の位置（第１の分布情報の画像上における拡大領域の位置）を点線で囲まれた矩形で示している。この矩形のように、拡大している領域の位置を示すガイドを表示することで、ユーザーはどの位置の拡大画像を見ているのかが把握しやすくなる。

拡大領域の位置を示すガイドをユーザーが移動させることで、拡大領域の位置を変更可能としてもよい。入力手段によりユーザーからのガイドの移動の指示がシステム制御部００４に入力されると、システム制御部００４は、表示制御手段００８にガイドの移動情報を出力する。表示制御手段００８は、このガイドの移動情報を受けて、画面上でガイドを移動させるとともに、変更後の拡大領域の拡大画像を表示手段００９に表示させる。

さらに、ガイドのサイズをユーザーが変更することで、拡大領域のサイズ（つまり拡大倍率）を変更可能としてもよい。入力手段によりユーザーからのガイドのサイズ変更の指示がシステム制御部００４に入力されると、システム制御部００４は、表示制御手段００８にガイドのサイズ変更情報を出力する。表示制御手段００８は、このガイドのサイズ変更情報を受けて、画面上でガイドのサイズを変更させるとともに、変更後の拡大領域の拡大画像を表示手段００９に表示させる。このように、ガイドの位置やサイズの変更を可能とすることにより、ユーザーが拡大して観察したい被検体内の領域の位置やサイズを容易に変更でき、操作性が向上する。

また、本実施形態においても、第１の分布情報の画像内の位置と第２の分布情報の画像内の位置とを対応付けるためのガイドを表示してもよい。もちろん、この位置を対応づけるためのガイドは第１の実施形態と同様に、ガイドを表示するモードと表示しないモードとをユーザーが選択できるようにするとよい。

さらに、本実施形態においても、表示制御手段００８は、上述の画像を並列表示する並列表示モードに加え、単独表示モードと重畳表示モードのうちいずれか一方のモードをさらに有するとよい。

また、第１の分布情報と並列表示する画像として、第１の分布情報と第２の分布情報との合成画像を、第１の分布情報の画像と並列表示しても本発明の効果は得られる。

＜第４の実施形態＞
本実施形態は、表示制御手段００８が、並列表示の指示情報が入力された場合に、第１の分布情報の画像と、第１及び第２の分布情報との合成画像と、を並列表示することを特徴とする。それ以外は第１乃至３の実施形態と同じである。本実施形態の被検体情報取得装置は、図１で示した装置と同様の構成の装置を用いる。また、表示方法の概略は図４で説明したフローと基本的に同じである。

本実施形態では、ユーザーからの並列表示の指示情報を受け、図４のＳ１０３のステップで、第１の分布情報と、第１及び第２の分布情報との合成画像と、を表示する。第１の分布情報の画像と第２の分布情報の画像との合成率は５０：５０のように予め決められていても良いし、ユーザーが任意に設定しても良い。

また、第３の実施形態のように、指定領域の画像を拡大して表示する場合は、その拡大率に応じて合成率を変化させても良い。図９に、拡大率と合成率との関係の一例を示す。図９では、拡大率が第１の所定値以下の場合は、第１と第２の分布情報の合成率は一定である。この場合、拡大率が低いので、第２の分布情報の画像の合成率（合成画像における第１の分布情報の画像の割合が高く、第２の分布情報の画像の割合）は低い。そして、拡大率が第１の所定値より高く第２の所定値より低い場合は、拡大率が高くなるに従って第２の分布情報の画像の合成率（合成画像中の第１の分布情報の画像に対する第２の分布情報の画像の割合）を高くしている。さらに、拡大率が第２の所定値以上の場合は、第１と第２の分布情報の合成率を一定としている。この場合、拡大率が高いため、第２の分布情報の画像の合成率を高くしている。

＜第５の実施形態＞
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した各実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。

００１ 探触子
００２ 変換素子
００３ 送信回路系
００４ システム制御部
００５ 受信回路系
００６ 固定型信号処理ブロック
００７ 適応型信号処理ブロック
００８ 表示制御手段
００９ 表示手段
０１０ 入力手段

Claims (17)

1. 被検体に弾性波を送信し、被検体内の各位置で反射した反射波を受信して複数の受信信号に変換する複数の変換素子と、
   前記複数の受信信号を用いて固定の重みで加算処理し、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得し、前記複数の走査線信号を用いて第１の分布情報を取得する固定型信号処理手段と、
   前記複数の受信信号を用いて、前記受信信号に応じて適応的に変化する重みで適応型信号処理を行い、第２の分布情報を取得する適応型信号処理手段と、
   前記第１の分布情報と、前記第２の分布情報と、が入力され、表示手段に画像情報を出力する表示制御手段と、
   を有し、
   前記表示制御手段は、
   前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記第１と第２の分布情報が合成された合成画像と、を同一画面内に並列表示させるための画像情報を出力することを特徴とする被検体情報取得装置。
2. 前記表示制御手段は、前記第１の分布情報の画像内における所定の位置と、前記第２の分布情報の画像内又は前記合成画像内における前記所定の位置に対応する位置と、を対応づけるためのガイドを表示させることを特徴とする請求項１に記載の被検体情報取得装置。
3. 前記表示制御手段は、
   前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像と、を並列表示するモードと、
   前記第１の分布情報の画像と、前記第２の分布情報の画像と、前記合成画像と、のうちいずれか一つの画像を単独で表示するモードと、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像と、を重畳表示するモードと、のうち少なくともいずれか一方のモードと、
   を選択的に実行可能に構成されていることを特徴とする請求項１又は２に記載の被検体情報取得装置。
4. 前記表示制御手段は、前記適応型信号処理の方法に応じて、前記第１の分布情報の画像と、前記第２の分布情報の画像と、を並列表示する場合の並べ方を変えることを特徴とする請求項１乃至３のいずれか１項に記載の被検体情報取得装置。
5. 前記適応型信号処理手段は、前記複数の受信信号を用いて、注目方向に関する感度を固定した状態で電力を最小化するように処理し、
   前記表示制御手段は、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像とを、前記第２の分布情報の画像内又は前記合成画像内における深さ方向と同じ方向に並ぶように並列表示させることを特徴とする請求項１乃至４のいずれか１項に記載の被検体情報取得装置。
6. 前記適応型信号処理手段は、前記複数の受信信号を用いて、深さ方向の注目位置に関する感度を固定した状態で電力を最小化するように処理し、
   前記表示制御手段は、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像とを、前記第２の分布情報の画像内又は前記合成画像内における深さ方向と交差する方向に並ぶように並列表示させることを特徴とする請求項１乃至４のいずれか１項に記載の被検体情報取得装置。
7. 前記表示制御手段は、
   前記第１の分布情報の画像と前記第２の分布情報の画像とのうち前記第１の分布情報のみが表示されている状態で、ユーザーにより入力された、前記第１の分布情報の画像内の指定領域の情報を受け、
   前記指定領域の位置の前記第１の分布情報の画像と、前記指定領域に対応する位置の前記第２の分布情報の画像又は前記合成画像と、を夫々拡大した拡大画像を並列表示させることを特徴とする請求項１乃至６のいずれか１項に記載の被検体情報取得装置。
8. 前記表示制御手段は、
   前記ユーザーにより入力された、前記第１の分布情報の画像内の指定領域の情報を受け、
   前記第１の分布情報の画像内の指定領域を示すガイドを表示させることを特徴とする請求項７に記載の被検体情報取得装置。
9. 被検体情報取得装置で取得された分布情報を用いて画像を表示手段に表示する表示方法であって、
   前記取得された分布情報は、
   被検体に弾性波を送信し被検体内で反射した反射波を受信することにより得られる受信信号を用いて固定の重みで加算処理し、前記被検体内の各位置からの反射波に対応する信号である走査線信号を複数取得し、複数の前記走査線信号を用いることで得られる第１の分布情報と、
   前記複数の受信信号を用いて前記受信信号に応じて適応的に変化する重みで適応型信号処理を行うことにより得られる第２の分布情報と、
   を含み、
   前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記第１と第２の分布情報が合成された合成画像と、を同一画面内に並列表示するステップを有することを特徴とする表示方法。
10. 前記第１の分布情報の画像内における所定の位置と、前記第２の分布情報の画像内又は前記合成画像内における前記所定の位置に対応する位置と、を対応づけるためのガイドを表示するステップをさらに有することを特徴とする請求項９に記載の表示方法。
11. 前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像と、を並列表示するモードと、
    前記第１の分布情報の画像と、前記第２の分布情報の画像と、前記合成画像と、のうちいずれか一つの画像を単独で表示するモードと、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像と、を重畳表示するモードと、のうち少なくともいずれか一方のモードと、
    を選択的に実行可能であることを特徴とする請求項９又は１０に記載の表示方法。
12. 前記並列表示するステップでは、前記適応型信号処理の方法に応じて、前記第１の分布情報の画像と、前記第２の分布情報の画像と、を並列表示する場合の並べ方を変えることを特徴とする請求項９乃至１１のいずれか１項に記載の表示方法。
13. 前記第２の分布情報は、前記複数の受信信号を用いて、注目方向に関する感度を固定した状態で電力を最小化するように適応型信号処理が行われることにより得られる分布情報であり、
    前記並列表示するステップでは、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像とを、前記第２の分布情報の画像内又は前記合成画像内における深さ方向と同じ方向に並ぶように並列表示することを特徴とする請求項９乃至１２のいずれか１項に記載の表示方法。
14. 前記第２の分布情報は、前記複数の受信信号を用いて、深さ方向の注目位置に関する感度を固定した状態で電力を最小化するように適応型信号処理が行われることにより得られる分布情報であり、
    前記並列表示するステップでは、前記第１の分布情報の画像と、前記第２の分布情報の画像又は前記合成画像とを、前記第２の分布情報の画像内又は前記合成画像内における深さ方向と交差する方向に並ぶように並列表示することを特徴とする請求項９乃至１２のいずれか１項に記載の表示方法。
15. 前記並列表示するステップの前に、前記第１の分布情報の画像と前記第２の分布情報の画像とのうち前記第１の分布情報の画像のみを表示するステップを有し、
    前記並列表示するステップでは、
    ユーザーにより入力された、前記第１の分布情報の画像内の指定領域の位置の、前記第１の分布情報の画像の拡大画像と、
    前記指定領域に対応する位置の、前記第２の分布情報の画像又は前記合成画像の拡大画像と、
    を並列表示することを特徴とする請求項９乃至１４のいずれか１項に記載の表示方法。
16. 前記ユーザーにより入力された、前記第１の分布情報の画像内の指定領域を示すガイドを表示するステップをさらに有することを特徴とする請求項１５に記載の表示方法。
17. 請求項９乃至１６のいずれか１項に記載の表示方法の各ステップをコンピュータに実行させるためのプログラム。