JP4473779B2 - 超音波診断装置及びその画像処理方法 - Google Patents

超音波診断装置及びその画像処理方法 Download PDF

Info

Publication number
JP4473779B2
JP4473779B2 JP2005150233A JP2005150233A JP4473779B2 JP 4473779 B2 JP4473779 B2 JP 4473779B2 JP 2005150233 A JP2005150233 A JP 2005150233A JP 2005150233 A JP2005150233 A JP 2005150233A JP 4473779 B2 JP4473779 B2 JP 4473779B2
Authority
JP
Japan
Prior art keywords
displacement
image
series
time
subject
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2005150233A
Other languages
English (en)
Other versions
JP2006325686A (ja
Inventor
正英 西浦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2005150233A priority Critical patent/JP4473779B2/ja
Priority to US11/438,584 priority patent/US20070016036A1/en
Publication of JP2006325686A publication Critical patent/JP2006325686A/ja
Application granted granted Critical
Publication of JP4473779B2 publication Critical patent/JP4473779B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995Combining images from different aspect angles, e.g. spatial compounding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は超音波診断装置に係わり、超音波診断装置で生成する画像のノイズを低減させる機能を有する超音波診断装置装置及びその画像処理方法に関する。
超音波診断装置では、超音波パルス波を体内に送信し、生体組織間で音響インピーダンスの差があることにより反射されたエコー信号を受信し、そのエコー信号の強度を画像化し、体内の組織形状等をみることができる。
エコー信号には種々の原因によるノイズが含まれ、これらのノイズ量が多いと診断上の障害となる。そのため、ノイズを低減させることは重要な課題である。すなわち、信号対雑音比S/N比を向上させることは重要な課題である。
種々のノイズの中でも自己相関性の低いランダムなノイズについてみると、複数の信号を加算することでS/N比を向上させることができる。これは、n個の信号を加算すると、自己相関性の高い本来の信号成分(生体からの反射信号)はn倍されるが、ランダムなノイズは期待値として√n倍しかされないという原理に基づいている。
したがって、超音波診断装置で得られる時系列の信号から複数の時相の信号を加算すれば、S/N比を向上させることができる。
しかし、この方法は複数の時相で対象物が同じ位置になければ効果がない。つまり対象物が動く場合にはこの方法を用いることはできない。
これに対して特許文献1、特許文献2では、周期的に運動する被検体について、動きが同位相となる複数の画像フレームを加算する方法を開示している。
特開平10−118061号公報 特開2001−170047公報
超音波診断装置で得られる画像のノイズが多いと診断上の障害となる。信号対雑音比(S/N比)を向上させるために、複数の信号を加算する方法が考えられるが、対象物に動きがある場合には効果がない。
また、周期的に運動する被検体について、動きが同位相となる複数の画像フレームを加算する方法では、対象が周期性をもった動きをしている必要があり、周期性のない動きをした対象物に用いることはできない。さらに、一周期分以上の画像データを保持しておくメモリが必要な他に、同位相となる画像フレームの対応づけが必要である。
そこで、本発明は、上記問題点を解決するためになされたものであり、対象物の運動の有無、運動の周期性の有無にかかわらずS/N比の向上した画像を得ることのできる超音波診断装置及びその方法を提供することを目的とする。
本発明は、超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、前記時系列の受信信号を時系列の変位検出用画像に変換する第1画像生成手段と、前記変換された時系列の変位検出用画像間の変位に関する情報を検出する変位検出手段と、前記変位に関する情報を逆変換して、その逆変換した変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正信号を複数生成する変位補正手段と、前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算手段と、前記加算受信信号を表示用画像に変換する第2画像生成手段と、を有すること特徴とする超音波診断装置である。
本発明は、超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、前記時系列の受信信号を、表示用画像と異なる解像度の時系列の変位検出用画像に変換する第1画像生成手段と、前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出手段と、前記時系列の受信信号を、時系列の前記表示用画像に変換する第2画像生成手段と、前記表示用画像と前記変位検出用画像との解像度に基づいて、前記変位に関する情報を変換し、前記表示用画像における前記被検体の動きによる位置ずれを補正した変位補正画像を生成する変位補正手段と、前記複数の変位補正画像を重み付け加算して加算画像を生成する加算手段と、を有することを特徴とする超音波診断装置である。
本発明は、超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、前記時系列の受信信号間の変位に関する情報を検出する変位検出手段と、前記変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正受信信号を複数生成する変位補正手段と、前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算手段と、前記加算受信信号を表示用画像に変換する画像生成手段と、を有することを特徴とする超音波診断装置である。
本発明によれば、超音波診断装置において、対象物の運動の有無、運動の周期性の有無にかかわらずS/N比の向上した画像を得ることができる。
以下、本発明の実施形態について図面に基づいて説明する。
[第1の実施形態]
本発明の第1の実施形態に係わる超音波診断装置10について図1から図5に基づいて説明する。
(1)超音波診断装置10の構成
図1は、本実施形態に係わる超音波診断装置10を示すブロック図である。
超音波診断装置10は、超音波を送受信する超音波プローブ12と、超音波プローブ12を駆動し超音波を送信すると共に、反射した超音波を受信しA/D変換処理を行い受信信号に変換する受信部14と、変換された受信信号を保持しておく受信バッファ16と、受信信号に対してフィルタ操作やスキャン変換操作を行い、画像データを生成する第1画像生成部18と、画像データに対して被検体の局所毎のフレーム間変位を検出する変位検出部20と、検出された変位情報を用いて複数フレームの受信信号を重み付け加算する信号加算部22と、加算された受信信号に対してフィルタ操作やスキャン変換操作を行い、表示用画像データを生成する第2画像生成部24と、表示用画像データを表示する表示部26とで構成される。
(2)受信部14と受信バッファ16
受信部14により駆動された超音波プローブ12から超音波ビームが射出される。運動する被検体で反射した超音波は超音波プローブ12で受波されてセクタスキャンされ(図2左参照)、受信部14で検波やA/D変換処理により時系列の受信信号に変換される(図2中参照)。時系列の受信信号は、超音波ビーム方向毎のデータの羅列であり、一連のビーム方向のデータにより一フレーム分のデータとなる。受信バッファ16は複数フレームにわたる時系列の受信信号を保持する。
(3)第1画像生成部18
受信バッファ16に保持された時系列の受信信号は、第1画像生成部18によりフィルタ操作やスキャン変換操作が施され時系列の変位検出用画像データとなる(図2参照)。第1画像生成部18では、生成する画像データの画像特性が後段の変位検出に適したものとなるように、フィルタ操作やスキャン変換操作のパラメータを設定すると良い。生成された画像データは第1画像生成部18に設けられた画像バッファに記憶される。この画像バッファは複数フレームの変位検出用画像データを保持する容量を持つ。
(4)変位検出部20
次に、変位検出部20では、画像バッファに保持されている時系列の変位検出用画像データを用いて時間的に隣接するフレーム間での運動する被検体の局所毎の変位を検出する。変位検出は、例えば局所毎のブロックマッチングによって検出することができる。
全ての位置で変位を検出しても良いが、目的の解像度に応じて離散位置で変位を検出し補間により全ての位置での変位を算出できるように構成しても良い。
変位検出処理により、変位検出用画像フレーム間での被検体の局所毎の動きが明らかになり、被検体が移動や変形をしても局所毎の対応する位置がわかることになる。
変位検出処理は、ブロックマッチングによる方法に限らず他の手法を用いることも可能である。
図3は変位検出用画像フレーム間の変位の例を示している。時間的に連続するi−1番目のフレーム画像とi番目のフレーム画像の間、及び、i番目のフレーム画像とi+1番目のフレーム画像の間での被検体の変位検出の例である。被検体の形状やテクスチャの有無等の画像の性質によって変位検出の困難さが異なるため、変位を検出しやすい位置での変位を求めると精度面や処理速度の点から見て効率が良い。
図3は、被検体(この図の例では心臓)の局所毎の変位が特徴点上で検出される様子を模式的に図示しており、図中の四角形が変位を検出する位置を示し、四角形から伸びる線分が変位の方向と大きさを示している。この特徴点は、画像上でのエッジ処理が行われ、画像上でエッジがある位置を示している。なお、この特徴点は画像上でコーナー検出あるいはピーク検出によって求められた位置を示しても構わない。
但し、変位検出を行う位置は必ずしも特徴点上で行われる必要はなく、全ての画素位置で行うことも可能である。
(5)信号加算部22
次に、信号加算部22での処理を行う。信号加算部22の処理は、変位逆変換処理、変位補正受信信号の生成処理、信号加算処理との3つに大別される。
(5−1)第1の処理
第1の処理としては、変位検出部20で検出された変位検出用画像フレーム間での変位情報をスキャン変換される前の受信信号の座標系での変位に逆変換する(図4左上参照)。
図3で例示した処理方法では、変位検出部20で検出される変位情報は特徴点位置のみのものである。そこで、特徴点上で得られた変位情報を空間的に補間することで任意の位置での変位を算出できるようにする(図4右上参照)。例えば、断面像の場合は2次元補間、3次元データの場合は変位情報も3次元ベクトルとなり補間も3次元的に行う。
そして、受信信号から画像データへのスキャン変換の逆変換によって受信信号の座標系での変位を生成する(図4左下参照)。
なお、変位情報の補間を画像データ上で行わずに、特徴点位置での変位情報を受信信号座標系に変換した後、受信信号上で空間的に補間を行う構成とすることも可能である。
(5−2)第2の処理
第2の処理としては、受信信号の座標系での変位を用いて変位を補正した変位補正受信信号を生成する(図4右下参照)。
i番目フレームのある位置(位置A)がi−1番目フレームのどの位置(位置B)であったかを受信信号の座標系での変位から読み出し、i−1番目フレームの受信信号の位置Bでのデータをi−1番目フレームに対応する変位補正受信信号の位置Aでのデータとする。
同様にi番目フレームのある位置(位置A)がi+1番目フレームのどの位置(位置B)になるかを受信信号の座標系での変位から読み出し、i+1番目フレームの受信信号の位置Bでのデータをi+1番目フレームに対応する変位補正受信信号の位置Aでのデータとする。
データ位置が整数でない場合は周囲のデータを補間したものを用いると良い。
このように変位補正受信信号を生成することで、i番目のフレームにおける被検体の位置にあわせた(動きによるずれが補正された)i−1番目及びi+1番目フレームの受信信号データから構成される変位補正受信信号を生成することができる。
(5−3)第3の処理
第3の処理としては、複数フレームの変位補正受信信号を重み付け加算し加算受信信号を生成する(図5参照)。
nフレームの信号を加算することで自己相関性の高い本来の信号成分(生体からの反射信号)はn倍されるが、ランダムなノイズは期待値として√n倍しかされないという原理によって、加算受信信号はS/N比が向上した信号となる。また、加算時に被検体の動きを補正しているので、被検体が移動・変形する場合でも効果を有する。
なお、図3、図5ではn=3の場合について例示している。より多くのフレームについて加算する構成とすることも可能である。その場合、連続したフレーム間の変位を求めてその変位を累積することで2フレーム以上離れたフレーム間の変位を算出することも可能であるが、注目するフレームと対象フレームとの間で直接変位を検出しても良い。
加算時の重み係数wについては、例えば1/nとすると各フレームの平均値をとることとなる。あるいは、変位検出結果の精度によって重み係数を変化させることも可能である。例えば、変位の検出精度が良い場合は各フレームの重み係数を均等にし、精度が悪い場合には注目フレーム(上記の例ではi番目フレーム)の重みを他のフレームに比べて大きくする。このようにすることで、変位の検出精度が悪い場合に生じる加算画像のぼけなどの悪影響を低減することができる。
同様に、加算に用いるフレーム数nを変位検出の精度あるいは変位の大きさによって動的に変化させることも可能である。例えば、検出精度の良いときは加算フレーム数nを大きく、検出精度の悪いときには加算フレーム数nを小さくすると、加算画像のぼけ等の誤差による悪影響を抑制しながら、S/N向上の利点を効果的に得ることができる。あるいは、変位の大/小によって加算フレーム数をそれぞれ小/大とすることで、より大きな効果を得ることができる。
なお、複数フレームを加算するときに、被検体の動きによっては対応位置がデータの範囲外となる場合が生じる。この場合には、その位置に関してデータの無いフレームを加算の対象外とし重み係数を調整すると良い。
(6)第2画像生成部24と表示部26
第2画像生成部24で、加算受信信号にフィルタ操作やスキャン変換操作を行い表示用画像データを生成する。生成された表示用画像データは表示部26によりユーザに提示される。
以上、本実施形態によれば、被検体が移動あるいは変形している場合においても、S/Nの向上した画像を得ることができる。また、受信信号を加算するため、表示用画像を生成する際のフィルタ操作や輝度変換等の影響を受けずに済むという利点があるさらに、画像データが多くの場合8ビットの256階調で表現されるのに対して、受信信号は一般的により多くのビット数(例えば、16ビット)で表現されているので、変位補正受信信号を重み付け加算し加算受信信号を生成することにより、S/Nの向上した画像を得ることができる。。
(7)変更例
第1画像生成部18で生成される画像データは変位検出に用いる変位検出用画像データであるので、第2画像生成部24で表示用画像データを生成する時と同じパラメータで画像生成する必要はない。変位検出に適した解像度や、フィルタ操作をパラメータとするように設定することができる。
[第2の実施形態]
本発明の第2の実施形態に係わる超音波診断装置10について図6から図8に基づいて説明する。本実施形態は、信号加算を画像データについて行う点で第1の実施形態と異なるものである。
図6は、本実施形態に係わる超音波診断装置10を示すブロック図である。
超音波診断装置10は、超音波を送受信する超音波プローブ12と、超音波プローブ12を駆動し超音波を送信すると共に、反射した超音波を受信しA/D変換処理を行い受信信号に変換する受信部14と、変換された受信信号を保持しておく受信バッファ16と、受信信号に対してフィルタ操作やスキャン変換操作を行い変位検出用画像データを生成する第1画像生成部18と、画像データに対して被検体の局所毎のフレーム間変位を検出する変位検出部20と、受信信号に対してフィルタ操作やスキャン変換操作を行い表示用画像データを生成する第2画像生成部24と、検出された変位情報を用いて複数フレームの表示用画像データを重み付け加算する信号加算部22と、加算された表示用画像データを表示する表示部26とで構成される。
超音波の送受信、受信信号から変位検出用画像データの生成、変位の検出部分については第1の実施形態と同様の処理の流れである。
本実施形態では、信号加算部22で加算されるデータは変位検出用画像データである。変位検出処理が行われる一方で、第2画像生成部24において表示用画像を生成する。第2画像生成部24では、受信バッファ16に保持されている受信信号からフィルタ操作やスキャン変換操作によって表示用画像データを生成する。
信号加算部22では、まず変位検出部20で検出された変位検出用画像フレーム間での変位情報を表示用画像フレームの座標系での変位に変換し、その変換された変位を用いて変位補正画像データを得る。図7は、変位検出用画像データと表示用画像データの解像度が異なる場合について変位の変換の例を示している。本実施形態においても、変位検出部20で検出される変位情報は特徴点位置のみのものであるので、特徴点上で得られた変位情報を空間的に補間することで任意の位置での変位を算出できるようにする(図7右上参照)。すなわち、断面像の場合は2次元補間、3次元データの場合は変位情報も3次元ベクトルとなり補間も3次元的に行う。その後変位情報を表示用画像データ上での変位情報に変換し(図7左下参照)、変位補正画像を生成する(図7右下参照)。
さらに、信号加算部22では、複数フレームの変位補正画像を重み付け加算し加算画像を生成する(図8参照)。画像データに対して加算を行った場合でも、nフレームの信号を加算することで自己相関性の高い本来の信号成分(生体からの反射信号)はn倍されるが、ランダムなノイズは期待値として√n倍しかされないという同様の原理によって、加算画像はS/N比が向上した画像となる。
最後に、加算画像を表示部26によってユーザに提示する。
本実施形態のような画像データを加算する構成とすると、変位の座標系変換がスケールの変換だけで済むことになるという利点を持つ。
なお、本実施形態では第1画像生成部18と第2画像生成部24を分けて構成したが、これら2つの画像生成部を共通化し、変位検出用画像データと表示用画像データを同一にする構成としても良い。このようにすると装置の構成を簡略化できる。
[第3の実施形態]
本発明の第3の実施形態に係わる超音波診断装置10について図9に基づいて説明する。本実施形態は、変位検出を受信信号上で行う点で第1の実施形態と異なるものである。
図9に本実施形態に係わる超音波診断装置10を示すブロック図を示す。
超音波診断装置10は、超音波を送受信する超音波プローブ12と、超音波プローブ12を駆動し超音波を送信すると共に、反射した超音波を受信しA/D変換処理を行い受信信号に変換する受信部14と、変換された受信信号を保持しておく受信バッファ16と、受信信号において被検体の局所毎のフレーム間変位を検出する変位検出部20と、検出された変位情報を用いて複数フレームの受信信号を重み付け加算する信号加算部22と、加算された加算受信信号に対してフィルタ操作やスキャン変換操作を行い表示用画像データを生成する画像生成部28と、表示用画像データを表示する表示部26とで構成される。
本実施形態では、超音波を送受信し受信信号を受信バッファ16に格納するまでは第1の実施形態と同様である。本実施形態では、被検体の変位検出を受信信号上で行う。画像データが多くの場合8ビットの256階調で表現されるのに対して、受信信号は一般的により多くのビット数で表現されている。また構成によっては、超音波信号を位相検波した後の複素信号として表現するように構成することも可能である。そのためより多くの情報を用いて変位を検出することが可能である。変位の検出は、受信信号を2次元データの配列とみなし、画像の場合と同様にブロックマッチング等の方法を用いて行うことができる。
信号加算部22では、検出された変位を用いて変位補正受信信号を生成し、複数フレームの変位補正受信信号を重み付け加算し加算受信信号を生成する。本実施形態では、変位検出する座標系と変位補正受信信号の座標系が同一であるので変換の必要はなく、任意の位置での変位を算出するための補間処理が必要なだけである。
最後に、加算受信信号にフィルタ操作やスキャン変換操作を行い表示用画像を生成する。生成された表示用画像は表示部26によりユーザに提示される。
本実施形態は、変位を変換する処理が不要であることに加え、画像生成部28も1つで済むためより簡単な構成とすることができる。
[第4の実施形態]
本発明の第4の実施形態に係わる超音波診断装置10について図10に基づいて説明する。本実施形態は、変位検出を受信信号上で行い、信号加算を画像データに対して行う点で第1の実施形態と異なるものである。
図10に本実施形態に係わる超音波診断装置10を示すブロック図を示す。
超音波診断装置10は、超音波を送受信する超音波プローブ12と、超音波プローブ12を駆動し超音波を送信すると共に、反射した超音波を受信しA/D変換処理を行い受信信号に変換する受信部14と、変換された受信信号を保持しておく受信バッファ16と、受信信号において被検体の局所毎のフレーム間変位を検出する変位検出部20と、受信信号に対してフィルタ操作やスキャン変換操作を行い表示用画像データを生成する画像生成部30と、検出された変位情報を用いて複数フレームの画像データを重み付け加算する信号加算部22と、加算された表示用画像データを表示する表示部26とで構成される。
本実施形態における変位検出処理までは第3の実施形態と同様の手順であり、被検体の変位検出を受信信号上で行う。
本実施形態では、信号加算部22で加算されるデータは画像データである。変位検出処理が行われる一方で、画像生成部30において表示用画像を生成する。画像生成部30では、受信バッファ16に保持されている受信信号からフィルタ操作やスキャン変換操作によって表示用画像データを生成する。
本実施形態における信号加算部22では、まず変位検出部20で検出された受信信号上での変位情報を表示用画像フレームの座標系での変位に変換し、その変換された変位を用いて変位補正画像データを得る。受信信号上では特徴点位置でのみ変位情報が検出されているので、任意の位置での変位を算出できるように空間的な補間処理を行う。その後、補間された変位情報を表示用画像データの座標系に変換し、これを用いて変位補正画像を生成する。空間的な補間処理は表示用画像データの座標系で行ってもよい。
さらに信号加算部22では、複数フレームの変位補正画像を重み付け加算し加算画像を生成する。
最後に、加算画像を表示部26によってユーザに提示する。
このような構成にすることで、変位情報の変換は受信信号の座標系から表示用画像データの座標系への変換となり、この変換は表示用の画像生成で通常必要となる変換であるので、この変換機能部分を画像生成部30と共用することができ、構成要素を減らすことが可能である。
[変更例]
本発明は上記各実施形態に限らず、その主旨を逸脱しない限り種々に変更することができる。
本発明の第1の実施形態に係わる超音波診断装置の構成を示すブロック図である。 超音波診断装置での超音波信号から画像データを生成する手順を説明した図である。 第1の実施形態に係わる超音波診断装置の変位検出部における変位検出処理を例示した図である。 第1の実施形態に係わる超音波診断装置の信号加算部における変位補正信号の生成手順を示した図である。 第1の実施形態に係わる超音波診断装置の信号加算部におけるフレームデータ加算演算を示した図である。 第2の実施形態に係わる超音波診断装置の構成を示すブロック図である。 第2の実施形態に係わる超音波診断装置の信号加算部における変位補正画像の生成手順を示した図である。 第2の実施形態に係わる超音波診断装置の信号加算部におけるフレーム画像加算演算を示した図である。 第3の実施形態に係わる超音波診断装置の構成を示すブロック図である。 第4の実施形態に係わる超音波診断装置の構成を示すブロック図である。
符号の説明
10 超音波診断装置
12 超音波プローブ
14 受信部
16 受信バッファ
18 第1画像生成部
20 変位検出部
22 信号加算部
24 第2画像生成部
26 表示部

Claims (13)

  1. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、
    前記時系列の受信信号を時系列の変位検出用画像に変換する第1画像生成手段と、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出手段と、
    前記変位に関する情報を逆変換して、その逆変換した変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正受信信号を複数生成する変位補正手段と、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算手段と、
    前記加算受信信号を表示用画像に変換する第2画像生成手段と、
    を有すること特徴とする超音波診断装置。
  2. 前記変位検出用画像と前記表示用画像をそれぞれ表示する表示手段を有する
    こと特徴とする請求項1記載の超音波診断装置。
  3. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、
    前記時系列の受信信号を、表示用画像と異なる解像度の時系列の変位検出用画像に変換する第1画像生成手段と、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出手段と、
    前記時系列の受信信号を、時系列の前記表示用画像に変換する第2画像生成手段と、
    前記表示用画像と前記変位検出用画像との解像度に基づいて、前記変位に関する情報を変換し、前記表示用画像における前記被検体の動きによる位置ずれを補正した変位補正画像を生成する変位補正手段と、
    前記複数の変位補正画像を重み付け加算して加算画像を生成する加算手段と、
    を有することを特徴とする超音波診断装置。
  4. 前記時系列の画像間の変位に関する情報は、前記画像上の複数の特徴点の変位である
    こと特徴とする請求項1または3記載の超音波診断装置。
  5. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置において、
    前記時系列の受信信号間の変位に関する情報を検出する変位検出手段と、
    前記変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正受信信号を複数生成する変位補正手段と、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算手段と、
    前記加算受信信号を表示用画像に変換する画像生成手段と、
    を有することを特徴とする超音波診断装置。
  6. 前記重み付けの重み係数は、加算する画像または受信信号のフレーム数nの逆数である
    こと特徴とする請求項1から5のいずれか一項に記載の超音波診断装置。
  7. 前記重み付けの重み係数は、前記変位検出結果の精度によって決定する
    こと特徴とする請求項1から5のいずれか一項に記載の超音波診断装置。
  8. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理方法において、
    前記時系列の受信信号を時系列の変位検出用画像に変換する第1画像生成ステップと、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出ステップと、
    前記変位に関する情報を逆変換して、その逆変換した変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正信号を複数生成する変位補正ステップと、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算ステップと、
    前記加算受信信号を表示用画像に変換する第2画像生成ステップと、
    を有すること特徴とする超音波診断装置における画像処理方法。
  9. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理方法において、
    前記時系列の受信信号を、表示用画像と異なる解像度の時系列の変位検出用画像に変換する第1画像生成ステップと、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出ステップと、
    前記時系列の受信信号を、時系列の前記表示用画像に変換する第2画像生成ステップと、
    前記表示用画像と前記変位検出用画像との解像度に基づいて、前記変位に関する情報を変換し、前記表示用画像における前記被検体の動きによる位置ずれを補正した変位補正画像を生成する変位補正ステップと、
    前記複数の変位補正画像を重み付け加算して加算画像を生成する加算ステップと、
    を有することを特徴とする超音波診断装置における画像処理方法。
  10. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理方法において、
    前記時系列の受信信号間の変位に関する情報を検出する変位検出ステップと、
    前記変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正受信信号を複数生成する変位補正ステップと、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算ステップと、
    前記加算受信信号を表示用画像に変換する画像生成ステップと、
    を有することを特徴とする超音波診断装置における画像処理方法。
  11. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理方法をコンピュータによって実現するプログラムにおいて、
    コンピュータに、
    前記時系列の受信信号を時系列の変位検出用画像に変換する第1画像生成機能と、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出機能と、
    前記変位に関する情報を逆変換して、その逆変換した変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正信号を複数生成する変位補正機能と、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算機能と、
    前記加算受信信号を表示用画像に変換する第2画像生成機能と、
    を実現させるための超音波診断装置における画像処理ログラム。
  12. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理ログラムにおいて、
    コンピュータに、
    前記時系列の受信信号を、表示用画像と異なる解像度の時系列の変位検出用画像に変換する第1画像生成機能と、
    前記時系列の変位検出用画像間の変位に関する情報を検出する変位検出機能と、
    前記時系列の受信信号を、時系列の前記表示用画像に変換する第2画像生成機能と、
    前記表示用画像と前記変位検出用画像との解像度に基づいて、前記変位に関する情報を変換し、前記表示用画像における前記被検体の動きによる位置ずれを補正した変位補正画像を生成する変位補正機能と、
    前記複数の変位補正画像を重み付け加算して加算画像を生成する加算機能と、
    を実現させるための超音波診断装置における画像処理ログラム。
  13. 超音波を送受信する超音波プローブを駆動し超音波を被検体に送信すると共に、前記被検体から反射した超音波を時系列に受信し、その時系列の受信信号を時系列の画像に変換して表示する超音波診断装置における画像処理ログラムにおいて、
    コンピュータに、
    前記時系列の受信信号間の変位に関する情報を検出する変位検出機能と、
    前記変位に関する情報に基づき、前記被検体の動きによる位置ずれを補正した変位補正受信信号を複数生成する変位補正機能と、
    前記複数の変位補正受信信号を重み付け加算して加算受信信号を生成する加算機能と、
    前記加算受信信号を表示用画像に変換する画像生成機能と、
    を実現させるための超音波診断装置における画像処理ログラム。
JP2005150233A 2005-05-23 2005-05-23 超音波診断装置及びその画像処理方法 Expired - Fee Related JP4473779B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2005150233A JP4473779B2 (ja) 2005-05-23 2005-05-23 超音波診断装置及びその画像処理方法
US11/438,584 US20070016036A1 (en) 2005-05-23 2006-05-22 Ultrasonic diagnostic apparatus and image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005150233A JP4473779B2 (ja) 2005-05-23 2005-05-23 超音波診断装置及びその画像処理方法

Publications (2)

Publication Number Publication Date
JP2006325686A JP2006325686A (ja) 2006-12-07
JP4473779B2 true JP4473779B2 (ja) 2010-06-02

Family

ID=37548035

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005150233A Expired - Fee Related JP4473779B2 (ja) 2005-05-23 2005-05-23 超音波診断装置及びその画像処理方法

Country Status (2)

Country Link
US (1) US20070016036A1 (ja)
JP (1) JP4473779B2 (ja)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090306505A1 (en) * 2006-02-22 2009-12-10 Hideki Yoshikawa Ultrasonic diagnostic apparatus
KR100983770B1 (ko) 2007-03-23 2010-09-28 주식회사 메디슨 초음파 탄성영상을 형성하기 위한 시스템 및 방법
JP2008256462A (ja) * 2007-04-03 2008-10-23 Fanuc Ltd 形状データの画像表示方法
WO2009104525A1 (ja) * 2008-02-18 2009-08-27 株式会社 日立メディコ 超音波診断装置、超音波弾性情報処理方法及び超音波弾性情報処理プログラム
JP5388492B2 (ja) * 2008-07-04 2014-01-15 株式会社東芝 画像表示装置および医用画像診断装置
WO2012008217A1 (ja) * 2010-07-14 2012-01-19 株式会社日立メディコ 超音波画像の画像復元方法およびその装置並びに超音波診断装置
JP2012196308A (ja) * 2011-03-22 2012-10-18 Fujifilm Corp 光音響画像生成装置及び方法
KR101982149B1 (ko) 2011-09-05 2019-05-27 삼성전자주식회사 의료 영상의 일부 정보를 활용한 장기 영상 생성 방법 및 장치
CN102727255B (zh) * 2012-07-13 2014-04-30 深圳市理邦精密仪器股份有限公司 一种超声图像空间复合成像的方法及装置
US10143381B2 (en) * 2013-04-19 2018-12-04 Canon Kabushiki Kaisha Object information acquiring apparatus and control method therefor
US20170315219A1 (en) * 2014-10-28 2017-11-02 Super Sonic Imagine Method and system for adjusting image gain
EP3235437A4 (en) * 2014-12-19 2018-09-12 Olympus Corporation Ultrasonic observation device
CN107847141B (zh) * 2015-07-09 2020-08-28 佳能株式会社 用于获取与多个图像数据集的位置位移有关的信息的设备、方法和程序
JP7336760B2 (ja) * 2019-02-04 2023-09-01 国立大学法人富山大学 超音波断層像生成方法、超音波断層像生成装置、およびプログラム

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5091782A (en) * 1990-04-09 1992-02-25 General Instrument Corporation Apparatus and method for adaptively compressing successive blocks of digital video
EP0477571B1 (en) * 1990-08-29 2001-03-21 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus capable of acquiring high quality image by correcting phase distortion contained in ultrasonic pulses
US5799111A (en) * 1991-06-14 1998-08-25 D.V.P. Technologies, Ltd. Apparatus and methods for smoothing images
US5447158A (en) * 1992-06-16 1995-09-05 Kabushiki Kaisha Toshiba Ultrasonic imaging method and system capable of displaying B-mode image and color flow mapping image over wide field
US5538004A (en) * 1995-02-28 1996-07-23 Hewlett-Packard Company Method and apparatus for tissue-centered scan conversion in an ultrasound imaging system
SE515535C2 (sv) * 1996-10-25 2001-08-27 Ericsson Telefon Ab L M En transkoder
US5873830A (en) * 1997-08-22 1999-02-23 Acuson Corporation Ultrasound imaging system and method for improving resolution and operation
CN1306649A (zh) * 1998-06-19 2001-08-01 赤道技术公司 把具有第一分辨率的编码图像直接解码为具有第二分辨率的解码图像
US7339617B1 (en) * 1999-03-12 2008-03-04 Sony Corporation Image providing device and its providing method, image processing device and processing method, and storage medium
US6775400B1 (en) * 1999-10-29 2004-08-10 Acuson Corporation Medical diagnostic ultrasonic imaging method and apparatus for suppressing electronic noise
US7085318B2 (en) * 2000-06-15 2006-08-01 Sony Corporation Image processing system, image processing method, program, and recording medium
US6901157B2 (en) * 2001-01-15 2005-05-31 Fuji Photo Film Co., Ltd. Ultrasonic diagnostic apparatus
JP3750797B2 (ja) * 2001-06-20 2006-03-01 ソニー株式会社 画像処理方法および装置
JP4958348B2 (ja) * 2001-09-06 2012-06-20 株式会社日立メディコ 超音波撮像装置
US7092584B2 (en) * 2002-01-04 2006-08-15 Time Warner Entertainment Company Lp Registration of separations
US7379496B2 (en) * 2002-09-04 2008-05-27 Microsoft Corporation Multi-resolution video coding and decoding
JP4214459B2 (ja) * 2003-02-13 2009-01-28 ソニー株式会社 信号処理装置および方法、記録媒体、並びにプログラム
KR100624751B1 (ko) * 2003-04-25 2006-09-19 (주)마크텍 영상에 워터마크를 삽입하는 방법 및 상기 방법을 이용한디지털 비디오 저장장치
US7542622B1 (en) * 2003-06-02 2009-06-02 The Trustees Of Columbia University In The City Of New York Spatio-temporal treatment of noisy images using brushlets
WO2005020821A1 (ja) * 2003-09-01 2005-03-10 Matsushita Electric Industrial Co., Ltd. 生体信号モニタ装置
US7004906B1 (en) * 2004-07-26 2006-02-28 Siemens Medical Solutions Usa, Inc. Contrast agent imaging with agent specific ultrasound detection
ITVR20060113A1 (it) * 2006-06-07 2008-01-07 Giglio Antonio Del Dispositivo per il trattamento del tessuto adiposo sottocutaneo mediante shockwaves non foicalizzate e contrapposte

Also Published As

Publication number Publication date
JP2006325686A (ja) 2006-12-07
US20070016036A1 (en) 2007-01-18

Similar Documents

Publication Publication Date Title
JP4473779B2 (ja) 超音波診断装置及びその画像処理方法
US8157737B2 (en) Method and apparatus for real time spatial compound imaging
JP5063515B2 (ja) 超音波診断装置
KR100971433B1 (ko) 사이드 로브의 영향을 제거하는 방법
CN101491450B (zh) 用于造影剂增强的医学诊断超声成像的同步组合
CN101396281B (zh) 用于彩色超声成像的扫描变换方法及其装置
US20130258805A1 (en) Methods and systems for producing compounded ultrasound images
US20200341142A1 (en) Motion-compensated processing
JP5065629B2 (ja) 超音波診断装置および超音波診断装置の制御プログラム
JP2004507298A (ja) 空間的に合成されたパノラマ画像を持つ超音波診断システム
US20130172747A1 (en) Estimating motion of particle based on vector doppler in ultrasound system
US9384568B2 (en) Method and system for enhanced frame rate upconversion in ultrasound imaging
JP2012105950A (ja) 超音波診断装置及びプログラム
KR20110104451A (ko) 초음파 진단 장치
KR100971425B1 (ko) 영상 움직임 추정 및 보상 장치와 그 방법
JP4426472B2 (ja) 超音波診断装置
JP5159480B2 (ja) 超音波診断装置および超音波診断装置の制御プログラム
US11151697B2 (en) Ultrasonic diagnosis device and program
WO2013150917A1 (ja) 超音波診断装置及び超音波画像の超解像生成方法
JP5665304B2 (ja) 超音波システム及び周期的に動く対象体の体積情報を提供する方法
JP5829008B2 (ja) 超音波診断装置
JP5786857B2 (ja) 超音波診断装置
JP5190248B2 (ja) 超音波診断装置
WO2016194432A1 (ja) 超音波診断装置
JP6132665B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061020

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090901

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090908

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091102

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20100209

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100305

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130312

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130312

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130312

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140312

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees