JPH0994248A - 超音波診断装置におけるスペックルノイズ判定方法及びスペックルノイズ判定除去回路を備えた超音波診断装置 - Google Patents
超音波診断装置におけるスペックルノイズ判定方法及びスペックルノイズ判定除去回路を備えた超音波診断装置Info
- Publication number
- JPH0994248A JPH0994248A JP7275107A JP27510795A JPH0994248A JP H0994248 A JPH0994248 A JP H0994248A JP 7275107 A JP7275107 A JP 7275107A JP 27510795 A JP27510795 A JP 27510795A JP H0994248 A JPH0994248 A JP H0994248A
- Authority
- JP
- Japan
- Prior art keywords
- image data
- speckle noise
- ultrasonic
- image
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000015654 memory Effects 0.000 claims abstract description 79
- 238000004364 calculation method Methods 0.000 claims description 37
- 239000000523 sample Substances 0.000 claims description 16
- 239000006185 dispersion Substances 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000009434 installation Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000005242 cardiac chamber Anatomy 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 206010027339 Menstruation irregular Diseases 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52077—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】
【課題】 超音波診断装置におけるスペックルノイズ判
定方法において、超音波画像に混入しているスペックル
ノイズを有効に判定する。 【解決手段】 超音波診断装置の超音波本体部のメモリ
部から連続した複数フレームの画像データを並列に読み
出し、これら複数の画像データを統計処理することによ
り、対象画像データがスペックルノイズであるか否かを
判定することができる。特に、上記複数の画像データの
統計処理として、複数フレームの画像データについて同
一番地のデータの平均輝度と平均分散とを求め、この平
均輝度と平均分散との比率を演算することにより、超音
波画像に混入しているスペックルノイズを有効に判定す
ることができる。
定方法において、超音波画像に混入しているスペックル
ノイズを有効に判定する。 【解決手段】 超音波診断装置の超音波本体部のメモリ
部から連続した複数フレームの画像データを並列に読み
出し、これら複数の画像データを統計処理することによ
り、対象画像データがスペックルノイズであるか否かを
判定することができる。特に、上記複数の画像データの
統計処理として、複数フレームの画像データについて同
一番地のデータの平均輝度と平均分散とを求め、この平
均輝度と平均分散との比率を演算することにより、超音
波画像に混入しているスペックルノイズを有効に判定す
ることができる。
Description
【０００１】
【発明の属する技術分野】本発明は、超音波を利用して
被検体内の診断部位について断層像を計測し表示する超
音波診断装置において、超音波の波長に比べて十分小さ
い生体組織の反射体群によりさまざまな位相で発生する
散乱波が干渉してランダムに小輝点群が見えるスペック
ルノイズを判定する方法及びそのスペックルノイズの判
定除去回路を備えた超音波診断装置に関する。
被検体内の診断部位について断層像を計測し表示する超
音波診断装置において、超音波の波長に比べて十分小さ
い生体組織の反射体群によりさまざまな位相で発生する
散乱波が干渉してランダムに小輝点群が見えるスペック
ルノイズを判定する方法及びそのスペックルノイズの判
定除去回路を備えた超音波診断装置に関する。
【０００２】
【従来の技術】一般に、超音波診断装置で得られる超音
波画像には、スペックルノイズと呼ばれる周期が不規則
な雑音が混入している。このスペックルノイズは、超音
波の波長に比べて十分小さい生体組織の反射体群により
さまざまな位相で発生する散乱波が干渉してランダムに
小輝点群が出現すると考えられている。そして、このス
ペックルノイズは、Ｂモード像の全域に本来の反射エコ
ー信号に重畳して出現することが知られている。すなわ
ち、例えば心臓の超音波画像において、心壁等の構造物
の反射エコー信号にスペックルノイズが重なって本来の
信号を乱し画像のチラツキ等の画質劣化を引き起こすと
共に、心腔等の反射体が無い部位等にもスペックルノイ
ズが出現して画質を劣化させるものであった。また、一
般に、スペックルノイズは、生体内に分布する微小散乱
体の配置により形状が大きく変化するとされており、生
体組織が動くと上記スペックルノイズのパターンも変わ
ることが知られている。しかし、そのパターンの変化の
度合、方向、輝度などの変わり方は一様ではない。この
ことから、生体が動いている場合は、スペックルノイズ
は不規則な変化をすることが確認されている。
波画像には、スペックルノイズと呼ばれる周期が不規則
な雑音が混入している。このスペックルノイズは、超音
波の波長に比べて十分小さい生体組織の反射体群により
さまざまな位相で発生する散乱波が干渉してランダムに
小輝点群が出現すると考えられている。そして、このス
ペックルノイズは、Ｂモード像の全域に本来の反射エコ
ー信号に重畳して出現することが知られている。すなわ
ち、例えば心臓の超音波画像において、心壁等の構造物
の反射エコー信号にスペックルノイズが重なって本来の
信号を乱し画像のチラツキ等の画質劣化を引き起こすと
共に、心腔等の反射体が無い部位等にもスペックルノイ
ズが出現して画質を劣化させるものであった。また、一
般に、スペックルノイズは、生体内に分布する微小散乱
体の配置により形状が大きく変化するとされており、生
体組織が動くと上記スペックルノイズのパターンも変わ
ることが知られている。しかし、そのパターンの変化の
度合、方向、輝度などの変わり方は一様ではない。この
ことから、生体が動いている場合は、スペックルノイズ
は不規則な変化をすることが確認されている。
【０００３】このような事情に対処して、従来の超音波
診断装置においては、図１０に示すようなフレーム相関
処理回路により経時的に重みを付けた加算平均化処理を
行い、動画像のスペックルノイズを低減していた。この
フレーム相関処理回路は、従来一般に「フレーム相関」
又は「ｓｃｃ（スキャンコリレーション）走査相関」と
呼ばれる手法でスペックルノイズを低減処理するもの
で、図１０に示すように、例えば重み付け加算器から成
る演算処理回路１と、Ｂモード像のデータを記録する機
能を有するフレームメモリ２と、このフレームメモリ２
の入力端子Ｄに入力される超音波画像のデータに同期し
て上記フレームメモリ２の番地を順次発生させる制御回
路３とから構成されている。
診断装置においては、図１０に示すようなフレーム相関
処理回路により経時的に重みを付けた加算平均化処理を
行い、動画像のスペックルノイズを低減していた。この
フレーム相関処理回路は、従来一般に「フレーム相関」
又は「ｓｃｃ（スキャンコリレーション）走査相関」と
呼ばれる手法でスペックルノイズを低減処理するもの
で、図１０に示すように、例えば重み付け加算器から成
る演算処理回路１と、Ｂモード像のデータを記録する機
能を有するフレームメモリ２と、このフレームメモリ２
の入力端子Ｄに入力される超音波画像のデータに同期し
て上記フレームメモリ２の番地を順次発生させる制御回
路３とから構成されている。
【０００４】上記フレーム相関処理回路において、フレ
ームメモリ２は、制御回路３により選択された番地に記
録されていたデータ、つまり１フレーム前のデータを出
力端子Ｅから出力する。この出力端子Ｅからの出力デー
タは、前記演算処理回路１のＢ入力に入力する。次に、
演算処理回路１は、図示外の超音波本体部から送出され
Ａ入力に入力してくる現時点でのデータと、上記Ｂ入力
に入力する１フレーム前のデータとを相関処理し、その
結果を出力端子Ｃより出力する。すると、この出力端子
Ｃからの出力データは、新しい時相のデータとしてフレ
ームメモリ２の入力端子Ｄへ送られる。次に、制御回路
３は、上記フレームメモリ２に新しい時相のデータが入
力された後に、これを前回と同一番地に記録するよう指
示を出す。
ームメモリ２は、制御回路３により選択された番地に記
録されていたデータ、つまり１フレーム前のデータを出
力端子Ｅから出力する。この出力端子Ｅからの出力デー
タは、前記演算処理回路１のＢ入力に入力する。次に、
演算処理回路１は、図示外の超音波本体部から送出され
Ａ入力に入力してくる現時点でのデータと、上記Ｂ入力
に入力する１フレーム前のデータとを相関処理し、その
結果を出力端子Ｃより出力する。すると、この出力端子
Ｃからの出力データは、新しい時相のデータとしてフレ
ームメモリ２の入力端子Ｄへ送られる。次に、制御回路
３は、上記フレームメモリ２に新しい時相のデータが入
力された後に、これを前回と同一番地に記録するよう指
示を出す。
【０００５】ここで、上記演算処理回路１で行われてい
る相関処理について説明する。いま、演算処理回路１の
Ａ入力に入力される数値をａで表し、Ｂ入力に入力され
る数値をｂで表し、出力端子Ｃに出現する演算結果をｃ
とすると、 ａ・Ｘ＋ｂ・Ｙ＝ｃ …（１） の関係があるＸ，Ｙの相関係数により、重み付け加算平
均の処理が行われる場合が多い。ここでは、Ｘ＜１，Ｙ
＜１であり、 Ｘ＋Ｙ＝１ …（２） の関係がある。一般には、被検体の診断部位や雑音の大
小により上記相関係数Ｘ，Ｙを選定して使用している。
る相関処理について説明する。いま、演算処理回路１の
Ａ入力に入力される数値をａで表し、Ｂ入力に入力され
る数値をｂで表し、出力端子Ｃに出現する演算結果をｃ
とすると、 ａ・Ｘ＋ｂ・Ｙ＝ｃ …（１） の関係があるＸ，Ｙの相関係数により、重み付け加算平
均の処理が行われる場合が多い。ここでは、Ｘ＜１，Ｙ
＜１であり、 Ｘ＋Ｙ＝１ …（２） の関係がある。一般には、被検体の診断部位や雑音の大
小により上記相関係数Ｘ，Ｙを選定して使用している。
【０００６】
【発明が解決しようとする課題】しかし、このような従
来の超音波診断装置におけるフレーム相関処理回路にお
いては、図１０に示す演算処理回路１のＢ入力に入力さ
れる画像データは１フレーム前の時相で既に演算処理さ
れたものであるので、上記演算処理回路１の出力端子Ｃ
からの演算出力は、経時的に旧フレームの影響が長時間
残留してしまうものであった。また、相関をとるフレー
ムの画像データ、例えば輝度レベル等の高低により相関
処理後の依存度が大幅に変化してしまうことがあった。
すなわち、輝度の高いフレームの画像データは長く影響
が残るが、輝度の低いフレームの画像データは非常に短
時間で無関係となってしまうものであった。これでは、
フレーム間で相関処理する時間間隔が不明となり、フレ
ームの指定が不可能となる。従って、上記相関処理の効
果は完全には推定できず、不確定要素となり、スペック
ルノイズ除去による画質の安定性は不十分であった。さ
らに、図１０に示すフレーム相関処理回路においては、
元の超音波画像の輝度レベル等のデータが演算処理回路
１の相関処理により順次変化して行き、輝度に応じたフ
ィルタ処理等が実行できないことがあった。そのため、
輝度情報により処理のアルゴリズムを変更することがで
きず、前述のように単純な重み付け加算程度しか行えな
かった。
来の超音波診断装置におけるフレーム相関処理回路にお
いては、図１０に示す演算処理回路１のＢ入力に入力さ
れる画像データは１フレーム前の時相で既に演算処理さ
れたものであるので、上記演算処理回路１の出力端子Ｃ
からの演算出力は、経時的に旧フレームの影響が長時間
残留してしまうものであった。また、相関をとるフレー
ムの画像データ、例えば輝度レベル等の高低により相関
処理後の依存度が大幅に変化してしまうことがあった。
すなわち、輝度の高いフレームの画像データは長く影響
が残るが、輝度の低いフレームの画像データは非常に短
時間で無関係となってしまうものであった。これでは、
フレーム間で相関処理する時間間隔が不明となり、フレ
ームの指定が不可能となる。従って、上記相関処理の効
果は完全には推定できず、不確定要素となり、スペック
ルノイズ除去による画質の安定性は不十分であった。さ
らに、図１０に示すフレーム相関処理回路においては、
元の超音波画像の輝度レベル等のデータが演算処理回路
１の相関処理により順次変化して行き、輝度に応じたフ
ィルタ処理等が実行できないことがあった。そのため、
輝度情報により処理のアルゴリズムを変更することがで
きず、前述のように単純な重み付け加算程度しか行えな
かった。
【０００７】そこで、本発明は、このような問題点に対
処し、超音波画像に混入しているスペックルノイズを有
効に判定すると共にこのスペックルノイズのみを除去し
て画質を向上することができる超音波診断装置における
スペックルノイズ判定方法及びスペックルノイズ判定除
去回路を提供することを目的とする。
処し、超音波画像に混入しているスペックルノイズを有
効に判定すると共にこのスペックルノイズのみを除去し
て画質を向上することができる超音波診断装置における
スペックルノイズ判定方法及びスペックルノイズ判定除
去回路を提供することを目的とする。
【０００８】
【課題を解決するための手段】上記目的を達成するため
に、第一の発明による超音波診断装置におけるスペック
ルノイズ判定方法は、被検体内に超音波を送受信する探
触子と、この探触子を駆動して超音波を発生させると共
に受信した反射エコー信号を処理する超音波送受信部及
びこの超音波送受信部からの反射エコー信号をディジタ
ル化し運動組織を含む被検体内の超音波画像のデータを
時系列に複数フレームに記録するメモリ部を備えた超音
波本体部と、上記メモリ部からのディジタル信号を超音
波ビームの走査線毎に書き込んで画像データを形成する
ディジタルスキャンコンバータ部と、このディジタルス
キャンコンバータ部からの画像データをアナログ変換し
画像として表示する画像表示部とを有する超音波診断装
置において、上記超音波本体部のメモリ部から連続した
複数フレームの画像データを並列に読み出し、これら複
数の画像データを統計処理することにより、対象画像デ
ータがスペックルノイズであるか否かを判定するもので
ある。
に、第一の発明による超音波診断装置におけるスペック
ルノイズ判定方法は、被検体内に超音波を送受信する探
触子と、この探触子を駆動して超音波を発生させると共
に受信した反射エコー信号を処理する超音波送受信部及
びこの超音波送受信部からの反射エコー信号をディジタ
ル化し運動組織を含む被検体内の超音波画像のデータを
時系列に複数フレームに記録するメモリ部を備えた超音
波本体部と、上記メモリ部からのディジタル信号を超音
波ビームの走査線毎に書き込んで画像データを形成する
ディジタルスキャンコンバータ部と、このディジタルス
キャンコンバータ部からの画像データをアナログ変換し
画像として表示する画像表示部とを有する超音波診断装
置において、上記超音波本体部のメモリ部から連続した
複数フレームの画像データを並列に読み出し、これら複
数の画像データを統計処理することにより、対象画像デ
ータがスペックルノイズであるか否かを判定するもので
ある。
【０００９】そして、上記複数の画像データの統計処理
は、複数フレームの画像データについて同一番地のデー
タの平均輝度と平均分散とを求め、この平均輝度と平均
分散との比率を演算するものとしてもよい。
は、複数フレームの画像データについて同一番地のデー
タの平均輝度と平均分散とを求め、この平均輝度と平均
分散との比率を演算するものとしてもよい。
【００１０】また、第二の発明による超音波診断装置
は、被検体内に超音波を送受信する探触子と、この探触
子を駆動して超音波を発生させると共に受信した反射エ
コー信号を処理する超音波送受信部及びこの超音波送受
信部からの反射エコー信号をディジタル化し運動組織を
含む被検体内の超音波画像のデータを時系列に複数フレ
ームに記録するメモリ部を備えた超音波本体部と、上記
メモリ部からのディジタル信号を超音波ビームの走査線
毎に書き込んで画像データを形成するディジタルスキャ
ンコンバータ部と、このディジタルスキャンコンバータ
部からの画像データをアナログ変換し画像として表示す
る画像表示部とを有する超音波診断装置において、上記
超音波本体部又はディジタルスキャンコンバータ部から
出力される連続した複数フレームの画像データを順次書
き込む複数個並列の画像メモリと、この画像メモリから
読み出した複数の画像データを取り込みそれらの画像デ
ータを統計処理することにより対象画像データがスペッ
クルノイズであるか否か判定すると共にそのスペックル
ノイズを除去する演算処理回路と、上記画像メモリ及び
演算処理回路の動作を制御する制御回路とを備え、上記
演算処理回路から前記ディジタルスキャンコンバータ部
又は画像表示部へ処理後の画像データを送出するスペッ
クルノイズ判定除去回路を備えたものである。
は、被検体内に超音波を送受信する探触子と、この探触
子を駆動して超音波を発生させると共に受信した反射エ
コー信号を処理する超音波送受信部及びこの超音波送受
信部からの反射エコー信号をディジタル化し運動組織を
含む被検体内の超音波画像のデータを時系列に複数フレ
ームに記録するメモリ部を備えた超音波本体部と、上記
メモリ部からのディジタル信号を超音波ビームの走査線
毎に書き込んで画像データを形成するディジタルスキャ
ンコンバータ部と、このディジタルスキャンコンバータ
部からの画像データをアナログ変換し画像として表示す
る画像表示部とを有する超音波診断装置において、上記
超音波本体部又はディジタルスキャンコンバータ部から
出力される連続した複数フレームの画像データを順次書
き込む複数個並列の画像メモリと、この画像メモリから
読み出した複数の画像データを取り込みそれらの画像デ
ータを統計処理することにより対象画像データがスペッ
クルノイズであるか否か判定すると共にそのスペックル
ノイズを除去する演算処理回路と、上記画像メモリ及び
演算処理回路の動作を制御する制御回路とを備え、上記
演算処理回路から前記ディジタルスキャンコンバータ部
又は画像表示部へ処理後の画像データを送出するスペッ
クルノイズ判定除去回路を備えたものである。
【００１１】そして、上記スペックルノイズ判定除去回
路の演算処理回路は、上記画像メモリから読み出した複
数フレームの画像データについて同一番地のデータの平
均輝度を求める平均輝度算出回路と、上記複数フレーム
の画像データについて同一番地のデータの平均分散を求
める平均分散算出回路と、上記それぞれの回路で求めた
算出結果を入力して平均輝度と平均分散との比率を演算
し対象画像データがスペックルノイズであるか否かを判
定するスペックルノイズ判定演算回路と、この判定結果
を入力して前記超音波本体部又はディジタルスキャンコ
ンバータ部からの画像データについてスペックルノイズ
を除去するスペックルノイズ除去処理回路とを有して成
るものとしてもよい。
路の演算処理回路は、上記画像メモリから読み出した複
数フレームの画像データについて同一番地のデータの平
均輝度を求める平均輝度算出回路と、上記複数フレーム
の画像データについて同一番地のデータの平均分散を求
める平均分散算出回路と、上記それぞれの回路で求めた
算出結果を入力して平均輝度と平均分散との比率を演算
し対象画像データがスペックルノイズであるか否かを判
定するスペックルノイズ判定演算回路と、この判定結果
を入力して前記超音波本体部又はディジタルスキャンコ
ンバータ部からの画像データについてスペックルノイズ
を除去するスペックルノイズ除去処理回路とを有して成
るものとしてもよい。
【００１２】
【発明の実施の形態】以下、本発明の実施の形態を添付
図面に基づいて詳細に説明する。第一の発明による超音
波診断装置におけるスペックルノイズ判定方法は、被検
体内に超音波を送受信する探触子と、この探触子を駆動
して超音波を発生させると共に受信した反射エコー信号
を処理する超音波送受信部及びこの超音波送受信部から
の反射エコー信号をディジタル化し運動組織を含む被検
体内の超音波画像のデータを時系列に複数フレームに記
録するメモリ部を備えた超音波本体部と、上記メモリ部
からのディジタル信号を超音波ビームの走査線毎に書き
込んで画像データを形成するディジタルスキャンコンバ
ータ部と、このディジタルスキャンコンバータ部からの
画像データをアナログ変換し画像として表示する画像表
示部とを有する超音波診断装置において、上記超音波本
体部のメモリ部から図１に示すように連続した複数フレ
ームｆ1，ｆ2，…，ｆnの画像データを同時に並列に読
み出し、これらの画像データをｆ1，ｆ2，…，ｆnのよ
うに正しく経時的に配列した後、上記複数の画像データ
（ｆ1〜ｆn）を統計処理することにより、対象画像デー
タがスペックルノイズであるか否かを判定するものであ
る。
図面に基づいて詳細に説明する。第一の発明による超音
波診断装置におけるスペックルノイズ判定方法は、被検
体内に超音波を送受信する探触子と、この探触子を駆動
して超音波を発生させると共に受信した反射エコー信号
を処理する超音波送受信部及びこの超音波送受信部から
の反射エコー信号をディジタル化し運動組織を含む被検
体内の超音波画像のデータを時系列に複数フレームに記
録するメモリ部を備えた超音波本体部と、上記メモリ部
からのディジタル信号を超音波ビームの走査線毎に書き
込んで画像データを形成するディジタルスキャンコンバ
ータ部と、このディジタルスキャンコンバータ部からの
画像データをアナログ変換し画像として表示する画像表
示部とを有する超音波診断装置において、上記超音波本
体部のメモリ部から図１に示すように連続した複数フレ
ームｆ1，ｆ2，…，ｆnの画像データを同時に並列に読
み出し、これらの画像データをｆ1，ｆ2，…，ｆnのよ
うに正しく経時的に配列した後、上記複数の画像データ
（ｆ1〜ｆn）を統計処理することにより、対象画像デー
タがスペックルノイズであるか否かを判定するものであ
る。
【００１３】この場合、図２に示す例えば心臓の断層像
５において、心壁６の内部の心腔内に出現するスペック
ルノイズ７は、上記心壁６などの組織よりも輝度が低
い。また、上記スペックルノイズ７の輝度レベルの変化
は不規則である。これを横軸を時間（フレーム番号）と
し縦軸を輝度として時間変化に対する輝度変化の状況を
グラフで示すと、図３のようになる。すなわち、心壁６
などの組織の運動等においてはカーブＣ1に示すように
輝度が高くそろっているが、スペックルノイズ７におい
てはカーブＣ2に示すように輝度が低くバラツイてい
る。なお、図２において、符号８は生体組織上に出現す
るスペックルノイズを示している。図３のグラフからわ
かるように、時系列的に並んだ複数枚（複数フレーム）
の画像データをもとに統計的手法を用いてスペックルノ
イズを判定することができる。
５において、心壁６の内部の心腔内に出現するスペック
ルノイズ７は、上記心壁６などの組織よりも輝度が低
い。また、上記スペックルノイズ７の輝度レベルの変化
は不規則である。これを横軸を時間（フレーム番号）と
し縦軸を輝度として時間変化に対する輝度変化の状況を
グラフで示すと、図３のようになる。すなわち、心壁６
などの組織の運動等においてはカーブＣ1に示すように
輝度が高くそろっているが、スペックルノイズ７におい
てはカーブＣ2に示すように輝度が低くバラツイてい
る。なお、図２において、符号８は生体組織上に出現す
るスペックルノイズを示している。図３のグラフからわ
かるように、時系列的に並んだ複数枚（複数フレーム）
の画像データをもとに統計的手法を用いてスペックルノ
イズを判定することができる。
【００１４】ここで、前記複数の画像データの統計処理
は、図１において正しく経時的に配列された複数フレー
ムｆ1，ｆ2，…，ｆnの画像データについて同一番地
（ｘ，ｙ）のデータの平均輝度と平均分散とを求め、こ
の平均輝度と平均分散との比率を演算するものである。
すなわち、図１に示すように、ある決められた枚数のフ
レームメモリの同一番地（ｘ，ｙ）、つまり断層像５の
視野の同じ位置の反射エコー信号を輝度変調した情報量
からその平均輝度及び平均分散を算出し、これらを利用
して演算する。いま、ｎフレームの同一番地（ｘ，ｙ）
のデータについての平均輝度をＰとし、平均分散をＱと
すると、スペックルノイズの判定関数Ｓは、 Ｓ＝Ｐ／Ｑ …（３） で与えられる。そして、この判定関数Ｓが所定値より大
きいか否かでスペックルノイズであるかどうかが判定可
能となる。すなわち、Ｓが大きい場合は固定部組織など
のデータであり、Ｓが小さい場合はスペックルノイズで
あると判定する。換言すると、平均輝度Ｐが小さくて平
均分散Ｑが大きい場合は、Ｓの値が小さくスペックルノ
イズと判定できる。この場合、平均分散Ｑが大きいとい
うことは、図３のカーブＣ2に示すようにデータが不規
則な変化をしていることを意味する。なお、式（３）の
判定関数Ｓに対する所定値は、対象物により自由に選択
して用いればよい。
は、図１において正しく経時的に配列された複数フレー
ムｆ1，ｆ2，…，ｆnの画像データについて同一番地
（ｘ，ｙ）のデータの平均輝度と平均分散とを求め、こ
の平均輝度と平均分散との比率を演算するものである。
すなわち、図１に示すように、ある決められた枚数のフ
レームメモリの同一番地（ｘ，ｙ）、つまり断層像５の
視野の同じ位置の反射エコー信号を輝度変調した情報量
からその平均輝度及び平均分散を算出し、これらを利用
して演算する。いま、ｎフレームの同一番地（ｘ，ｙ）
のデータについての平均輝度をＰとし、平均分散をＱと
すると、スペックルノイズの判定関数Ｓは、 Ｓ＝Ｐ／Ｑ …（３） で与えられる。そして、この判定関数Ｓが所定値より大
きいか否かでスペックルノイズであるかどうかが判定可
能となる。すなわち、Ｓが大きい場合は固定部組織など
のデータであり、Ｓが小さい場合はスペックルノイズで
あると判定する。換言すると、平均輝度Ｐが小さくて平
均分散Ｑが大きい場合は、Ｓの値が小さくスペックルノ
イズと判定できる。この場合、平均分散Ｑが大きいとい
うことは、図３のカーブＣ2に示すようにデータが不規
則な変化をしていることを意味する。なお、式（３）の
判定関数Ｓに対する所定値は、対象物により自由に選択
して用いればよい。
【００１５】また、図３に示す時間変化に対する輝度変
化の状況を示すグラフにおいて、カーブＣ1に示す高輝
度を保っていた生体内の組織からの輝度データが特定フ
レームのみ低い場合は、データの欠落であると判定す
る。さらに、図３において、カーブＣ2に示すスペック
ルノイズの輝度レベルを保っていた輝度データが突然に
特定フレームのみ上昇した場合は、その他のノイズであ
ると判定する。以上のことを基準として、複数フレーム
のデータを比較してスペックルノイズであるか否かの判
定を行うことができる。
化の状況を示すグラフにおいて、カーブＣ1に示す高輝
度を保っていた生体内の組織からの輝度データが特定フ
レームのみ低い場合は、データの欠落であると判定す
る。さらに、図３において、カーブＣ2に示すスペック
ルノイズの輝度レベルを保っていた輝度データが突然に
特定フレームのみ上昇した場合は、その他のノイズであ
ると判定する。以上のことを基準として、複数フレーム
のデータを比較してスペックルノイズであるか否かの判
定を行うことができる。
【００１６】次に、これらの判定基準により（ａ）スペ
ックルノイズ、（ｂ）データの欠落、（ｃ）その他のノ
イズの判定を、複数フレームのデータを図４に示すよう
にグラフの形状で分類して行う場合について説明する。
この図４は、例えば５個のフレームｆ1〜ｆ5の輝度デー
タを用い、（ａ）,（ｂ）,（ｃ）図はスペックルノイズ
の判定を目的とした各フレームｆ1〜ｆ5の同一番地
（ｘ，ｙ）のデータを経時的に示したものである。初め
に、図４（ａ）は、５個のデータｆ1〜ｆ5が総て低い輝
度レベルの範囲に分布している場合を示している。ま
た、それぞれの輝度値はランダムで不規則である。この
（ａ）に示すグラフのパターンは、図２に示すように、
心壁６の内部のように反射対象物の存在しない心腔内な
どに発生しているスペックルノイズ７を特徴付けると判
定できる。次に、図４（ｂ）は、特定のフレーム、例え
ば第３フレームｆ3のみがスペックルノイズのデータで
ある場合を示す。この（ｂ）の場合は、ｆ3の前後のフ
レームのデータの大きさからして、電源回路等のゆらぎ
などに起因するいわゆる黒抜けノイズと判定される。次
に、図４（ｃ）は、特定フレーム、例えば第３フレーム
ｆ3のみが高輝度で他のフレームはスペックルノイズの
輝度レベルである場合を示す。この（ｃ）の場合は、上
記第３フレームｆ3のデータの大きさからして、電気回
路や外乱ノイズによるペッパーノイズと判定される。な
お、図４に示す（ａ）,（ｂ）,（ｃ）の経時的なデータ
によるグラフ形状は、実験的に求められるものである。
また、図４は超音波画像の視野のある１点に注目したグ
ラフであるが、実際には得られた超音波画像の視野全体
について処理を行う。
ックルノイズ、（ｂ）データの欠落、（ｃ）その他のノ
イズの判定を、複数フレームのデータを図４に示すよう
にグラフの形状で分類して行う場合について説明する。
この図４は、例えば５個のフレームｆ1〜ｆ5の輝度デー
タを用い、（ａ）,（ｂ）,（ｃ）図はスペックルノイズ
の判定を目的とした各フレームｆ1〜ｆ5の同一番地
（ｘ，ｙ）のデータを経時的に示したものである。初め
に、図４（ａ）は、５個のデータｆ1〜ｆ5が総て低い輝
度レベルの範囲に分布している場合を示している。ま
た、それぞれの輝度値はランダムで不規則である。この
（ａ）に示すグラフのパターンは、図２に示すように、
心壁６の内部のように反射対象物の存在しない心腔内な
どに発生しているスペックルノイズ７を特徴付けると判
定できる。次に、図４（ｂ）は、特定のフレーム、例え
ば第３フレームｆ3のみがスペックルノイズのデータで
ある場合を示す。この（ｂ）の場合は、ｆ3の前後のフ
レームのデータの大きさからして、電源回路等のゆらぎ
などに起因するいわゆる黒抜けノイズと判定される。次
に、図４（ｃ）は、特定フレーム、例えば第３フレーム
ｆ3のみが高輝度で他のフレームはスペックルノイズの
輝度レベルである場合を示す。この（ｃ）の場合は、上
記第３フレームｆ3のデータの大きさからして、電気回
路や外乱ノイズによるペッパーノイズと判定される。な
お、図４に示す（ａ）,（ｂ）,（ｃ）の経時的なデータ
によるグラフ形状は、実験的に求められるものである。
また、図４は超音波画像の視野のある１点に注目したグ
ラフであるが、実際には得られた超音波画像の視野全体
について処理を行う。
【００１７】次に、第二の発明による超音波診断装置に
ついて、図５〜図８を参照して説明する。図５は、本発
明のスペックルノイズ判定除去回路を備えた超音波診断
装置の全体構成を示すブロック図である。この超音波診
断装置は、超音波を利用して被検体内の診断部位につい
て断層像を計測し表示するもので、図５に示すように、
探触子１０と、超音波本体部１１と、ディジタルスキャ
ンコンバータ部（以下「ＤＳＣ部」と略称する）１２
と、画像表示部１３とを有し、さらにスペックルノイズ
判定除去回路１４を備えて成る。
ついて、図５〜図８を参照して説明する。図５は、本発
明のスペックルノイズ判定除去回路を備えた超音波診断
装置の全体構成を示すブロック図である。この超音波診
断装置は、超音波を利用して被検体内の診断部位につい
て断層像を計測し表示するもので、図５に示すように、
探触子１０と、超音波本体部１１と、ディジタルスキャ
ンコンバータ部（以下「ＤＳＣ部」と略称する）１２
と、画像表示部１３とを有し、さらにスペックルノイズ
判定除去回路１４を備えて成る。
【００１８】上記探触子１０は、被検体内に超音波を送
受信するもので、図示省略したがその内部には超音波の
発生源であると共に反射エコーを受信して電気信号に変
換する超音波振動子を有している。超音波本体部１１
は、上記探触子１０を駆動して超音波を送受信させると
共に得られた超音波画像のデータを記録するもので、図
示省略したがその内部には、上記探触子１０を駆動して
超音波を発生させると共に受信した反射エコー信号を処
理する超音波送受信部と、この超音波送受信部からの反
射エコー信号をディジタル化し運動組織を含む被検体内
の超音波画像のデータを時系列に複数フレームに記録す
るメモリ部とを有している。また、ＤＳＣ部１２は、上
記メモリ部からのディジタル信号を超音波ビームの走査
線毎に書き込んで画像データを形成するもので、その内
部にはディジタルの超音波走査線メモリ群と、そのメモ
リ制御器とを有している。さらに、画像表示部１３は、
上記ＤＳＣ部１２から出力される画像データを入力して
アナログ変換し画像として表示するもので、その内部に
はＤ／Ａ変換器と、テレビモニタとを有している。
受信するもので、図示省略したがその内部には超音波の
発生源であると共に反射エコーを受信して電気信号に変
換する超音波振動子を有している。超音波本体部１１
は、上記探触子１０を駆動して超音波を送受信させると
共に得られた超音波画像のデータを記録するもので、図
示省略したがその内部には、上記探触子１０を駆動して
超音波を発生させると共に受信した反射エコー信号を処
理する超音波送受信部と、この超音波送受信部からの反
射エコー信号をディジタル化し運動組織を含む被検体内
の超音波画像のデータを時系列に複数フレームに記録す
るメモリ部とを有している。また、ＤＳＣ部１２は、上
記メモリ部からのディジタル信号を超音波ビームの走査
線毎に書き込んで画像データを形成するもので、その内
部にはディジタルの超音波走査線メモリ群と、そのメモ
リ制御器とを有している。さらに、画像表示部１３は、
上記ＤＳＣ部１２から出力される画像データを入力して
アナログ変換し画像として表示するもので、その内部に
はＤ／Ａ変換器と、テレビモニタとを有している。
【００１９】ここで、本発明においては、上記超音波本
体部１１とＤＳＣ部１２との間にスペックルノイズ判定
除去回路１４が設けられている。このスペックルノイズ
判定除去回路１４は、対象画像データについてスペック
ルノイズであるか否かを判定すると共にスペックルノイ
ズと判定したらそのデータを除去するもので、その内部
構成は図６に示すように、データ分配器１５と、複数個
の画像メモリ１６ａ〜１６ｅと、データ選別器１７と、
演算処理回路１８と、制御回路１９とから成る。
体部１１とＤＳＣ部１２との間にスペックルノイズ判定
除去回路１４が設けられている。このスペックルノイズ
判定除去回路１４は、対象画像データについてスペック
ルノイズであるか否かを判定すると共にスペックルノイ
ズと判定したらそのデータを除去するもので、その内部
構成は図６に示すように、データ分配器１５と、複数個
の画像メモリ１６ａ〜１６ｅと、データ選別器１７と、
演算処理回路１８と、制御回路１９とから成る。
【００２０】上記データ分配器１５は、前記超音波本体
部１１から出力された画像データを入力し、制御回路１
９の制御により複数個の画像メモリ１６ａ〜１６ｅに対
して１フレーム毎に切り換えて上記画像データを送出す
るものである。複数個の画像メモリ１６ａ〜１６ｅは、
上記データ分配器１５で分配された連続する複数フレー
ムの画像データを１フレーム毎に順次入力して書き込む
もので、図６では例えば５個並列に設けられている。な
お、この画像メモリは、５個に限られず、処理の目的及
び精度等により他の任意の個数としてもよい。データ選
別器１７は、上記複数個の画像メモリ１６ａ〜１６ｅか
ら読み出された画像データを入力し、制御回路１９の制
御により経時的に古いものから新しいものへと正しく並
び換えて画像データを送出するものである。演算処理回
路１８は、上記データ選別器１７で正しく経時的に配列
された複数フレームの画像データを取り込み、それらの
画像データを統計処理することにより対象画像データが
スペックルノイズであるか否か判定すると共にそのスペ
ックルノイズを除去するもので、その演算処理後の画像
データを前記ＤＳＣ部１２へ送出するようになってい
る。そして、制御回路１９は、上記の各構成要素の動作
を制御するものである。
部１１から出力された画像データを入力し、制御回路１
９の制御により複数個の画像メモリ１６ａ〜１６ｅに対
して１フレーム毎に切り換えて上記画像データを送出す
るものである。複数個の画像メモリ１６ａ〜１６ｅは、
上記データ分配器１５で分配された連続する複数フレー
ムの画像データを１フレーム毎に順次入力して書き込む
もので、図６では例えば５個並列に設けられている。な
お、この画像メモリは、５個に限られず、処理の目的及
び精度等により他の任意の個数としてもよい。データ選
別器１７は、上記複数個の画像メモリ１６ａ〜１６ｅか
ら読み出された画像データを入力し、制御回路１９の制
御により経時的に古いものから新しいものへと正しく並
び換えて画像データを送出するものである。演算処理回
路１８は、上記データ選別器１７で正しく経時的に配列
された複数フレームの画像データを取り込み、それらの
画像データを統計処理することにより対象画像データが
スペックルノイズであるか否か判定すると共にそのスペ
ックルノイズを除去するもので、その演算処理後の画像
データを前記ＤＳＣ部１２へ送出するようになってい
る。そして、制御回路１９は、上記の各構成要素の動作
を制御するものである。
【００２１】このような構成のスペックルノイズ判定除
去回路１４において、複数個の画像メモリ１６ａ〜１６
ｅに記録される画像データは周期的に更新され、制御回
路１９は指定した画像メモリに常に最新のデータが記録
されるように制御している。このとき、上記複数個の画
像メモリ１６ａ〜１６ｅはそのアドレスが総て共通とさ
れ、記録時には同一番地のデータが出力に現れる方式と
されている。すなわち、画像メモリ１６ａ〜１６ｅへの
データ記録時には新旧フレームの同一番地のデータが同
時に出力され、データ選別器１７へ入力する。しかし、
このデータは上記複数個の画像メモリ１６ａ〜１６ｅの
うちどれが最新であるかは、制御回路１９により指示を
受けなければわからない。そこで、上記データ選別器１
７は、制御回路１９から指示を受けて古いデータから新
しいデータへと経時的に再配列し、この再配列後のデー
タを演算処理回路１８へ送る。なお、上記データ選別器
１７は、データ長を例えば６ビットとすると、６ビット
×５個で30ビットのデータを選別する能力を要する。従
って、図６の実施例では、５×５のマトリックススイッ
チが６個あればデータ選別を実現することができる。そ
して、上記演算処理回路１８は、データ選別器１７から
送られた再配列後のデータを入力し、正しく経時的に配
列された複数フレームの画像データを比較して統計処理
する。この場合は、従来手法と異なり、設定したフレー
ムの枚数、すなわち心時相等の範囲のデータが正しく相
関等の演算処理に用いられ、対象画像データがスペック
ルノイズであるか否か判定すると共にそのスペックルノ
イズを除去することができる。なお、図６において、デ
ータ選別器１７は必ずしも設けなくてもよい。
去回路１４において、複数個の画像メモリ１６ａ〜１６
ｅに記録される画像データは周期的に更新され、制御回
路１９は指定した画像メモリに常に最新のデータが記録
されるように制御している。このとき、上記複数個の画
像メモリ１６ａ〜１６ｅはそのアドレスが総て共通とさ
れ、記録時には同一番地のデータが出力に現れる方式と
されている。すなわち、画像メモリ１６ａ〜１６ｅへの
データ記録時には新旧フレームの同一番地のデータが同
時に出力され、データ選別器１７へ入力する。しかし、
このデータは上記複数個の画像メモリ１６ａ〜１６ｅの
うちどれが最新であるかは、制御回路１９により指示を
受けなければわからない。そこで、上記データ選別器１
７は、制御回路１９から指示を受けて古いデータから新
しいデータへと経時的に再配列し、この再配列後のデー
タを演算処理回路１８へ送る。なお、上記データ選別器
１７は、データ長を例えば６ビットとすると、６ビット
×５個で30ビットのデータを選別する能力を要する。従
って、図６の実施例では、５×５のマトリックススイッ
チが６個あればデータ選別を実現することができる。そ
して、上記演算処理回路１８は、データ選別器１７から
送られた再配列後のデータを入力し、正しく経時的に配
列された複数フレームの画像データを比較して統計処理
する。この場合は、従来手法と異なり、設定したフレー
ムの枚数、すなわち心時相等の範囲のデータが正しく相
関等の演算処理に用いられ、対象画像データがスペック
ルノイズであるか否か判定すると共にそのスペックルノ
イズを除去することができる。なお、図６において、デ
ータ選別器１７は必ずしも設けなくてもよい。
【００２２】図７は上記演算処理回路１８の具体的な内
部構成例を図示したスペックルノイズ判定除去回路１
４′のブロック図である。この実施例では、図６に示す
データ分配器１５の機能を制御回路１９で実行させるも
のとし、該データ分配器１５は省略してある。そして、
上記演算処理回路１８は、平均輝度算出回路２０と、複
数個の個別分散回路２１ａ〜２１ｅと、平均分散算出回
路２２と、スペックルノイズ判定演算回路２３と、スペ
ックルノイズ除去処理回路２４とから成る。
部構成例を図示したスペックルノイズ判定除去回路１
４′のブロック図である。この実施例では、図６に示す
データ分配器１５の機能を制御回路１９で実行させるも
のとし、該データ分配器１５は省略してある。そして、
上記演算処理回路１８は、平均輝度算出回路２０と、複
数個の個別分散回路２１ａ〜２１ｅと、平均分散算出回
路２２と、スペックルノイズ判定演算回路２３と、スペ
ックルノイズ除去処理回路２４とから成る。
【００２３】上記平均輝度算出回路２０は、複数個たと
えば５個の画像メモリ１６ａ〜１６ｅから読み出し経時
的に配列された複数フレームの画像データについて同一
番地のデータの平均輝度Ｐ（前述の式(３)参照）を求め
るものである。個別分散回路２１ａ〜２１ｅは、上記平
均輝度算出回路２０で求めた平均輝度Ｐと、上記画像メ
モリ１６ａ〜１６ｅの各々から読み出した各画像データ
とを入力して、それぞれ両者の差の絶対値を算出して個
別分散値を求めるもので、画像メモリ１６ａ〜１６ｅの
個数に合わせて例えば５個設けられている。平均分散算
出回路２２は、上記個別分散回路２１ａ〜２１ｅからそ
れぞれ出力された個別分散値を入力して全体の平均値を
算出することにより経時的に配列された複数フレームの
画像データについて同一番地のデータの平均分散Ｑ（前
述の式(３)参照）を求めるものである。また、スペック
ルノイズ判定演算回路２３は、上記平均輝度算出回路２
０と平均分散算出回路２２とで、求めた算出結果を入力
して平均輝度Ｐと平均分散Ｑとの比率Ｓ（前述の式(３)
参照）を演算し対象画像データがスペックルノイズであ
るか否かを判定するものである。さらに、スペックルノ
イズ除去処理回路２４は、上記スペックルノイズ判定演
算回路２３の判定結果Ｒを入力して前記超音波本体部１
１から入力する現時点の画像データについてスペックル
ノイズを除去するものである。そして、このスペックル
ノイズ除去処理回路２４からの出力データは、前記ＤＳ
Ｃ部１２へ送られるようになっている。
えば５個の画像メモリ１６ａ〜１６ｅから読み出し経時
的に配列された複数フレームの画像データについて同一
番地のデータの平均輝度Ｐ（前述の式(３)参照）を求め
るものである。個別分散回路２１ａ〜２１ｅは、上記平
均輝度算出回路２０で求めた平均輝度Ｐと、上記画像メ
モリ１６ａ〜１６ｅの各々から読み出した各画像データ
とを入力して、それぞれ両者の差の絶対値を算出して個
別分散値を求めるもので、画像メモリ１６ａ〜１６ｅの
個数に合わせて例えば５個設けられている。平均分散算
出回路２２は、上記個別分散回路２１ａ〜２１ｅからそ
れぞれ出力された個別分散値を入力して全体の平均値を
算出することにより経時的に配列された複数フレームの
画像データについて同一番地のデータの平均分散Ｑ（前
述の式(３)参照）を求めるものである。また、スペック
ルノイズ判定演算回路２３は、上記平均輝度算出回路２
０と平均分散算出回路２２とで、求めた算出結果を入力
して平均輝度Ｐと平均分散Ｑとの比率Ｓ（前述の式(３)
参照）を演算し対象画像データがスペックルノイズであ
るか否かを判定するものである。さらに、スペックルノ
イズ除去処理回路２４は、上記スペックルノイズ判定演
算回路２３の判定結果Ｒを入力して前記超音波本体部１
１から入力する現時点の画像データについてスペックル
ノイズを除去するものである。そして、このスペックル
ノイズ除去処理回路２４からの出力データは、前記ＤＳ
Ｃ部１２へ送られるようになっている。
【００２４】次に、このように構成された演算処理回路
１８の動作について、図８に示すタイミング線図を参照
して説明する。図８（ａ）は、図７に示す例えば５個の
画像メモリ１６ａ〜１６ｅに対するフレーム更新信号の
タイミングを示しており、フレームの書換開始時に出力
される。図８（ｂ）は、上記各画像メモリ１６ａ〜１６
ｅのフレームデータであり、例えば８ビットのデータが
５フレーム使用されており、図中でｆ1，ｆ2，ｆ3，…
はフレーム番号を示している。図８（ｃ）は、第一の画
像メモリ１６ａに入力する書き込み許可信号Ｗ1を示し
ている。図８（ｄ）は第二の画像メモリ１６ｂに入力す
る書き込み許可信号Ｗ2を示し、図８（ｅ）は第三の画
像メモリ１６ｃに入力する書き込み許可信号Ｗ3を示
し、…、図８（ｇ）は第五の画像メモリ１６ｅに入力す
る書き込み許可信号Ｗ5を示している。これらの書き込
み許可信号Ｗ1〜Ｗ5は、図７に示す制御回路１９から各
画像メモリ１６ａ〜１６ｅに送出されるメモリ制御信号
Ｓ1に含まれる。この図８（ａ）〜（ｇ）からわかるよ
うに、各画像メモリ１６ａ〜１６ｅは、１フレーム毎に
経時的にデータが書き込まれて行く。この場合、画像メ
モリ１６ａ〜１６ｅは例えば５フレームとされており、
５フレーム毎に各画像メモリ１６ａ〜１６ｅが書き込み
許可される。
１８の動作について、図８に示すタイミング線図を参照
して説明する。図８（ａ）は、図７に示す例えば５個の
画像メモリ１６ａ〜１６ｅに対するフレーム更新信号の
タイミングを示しており、フレームの書換開始時に出力
される。図８（ｂ）は、上記各画像メモリ１６ａ〜１６
ｅのフレームデータであり、例えば８ビットのデータが
５フレーム使用されており、図中でｆ1，ｆ2，ｆ3，…
はフレーム番号を示している。図８（ｃ）は、第一の画
像メモリ１６ａに入力する書き込み許可信号Ｗ1を示し
ている。図８（ｄ）は第二の画像メモリ１６ｂに入力す
る書き込み許可信号Ｗ2を示し、図８（ｅ）は第三の画
像メモリ１６ｃに入力する書き込み許可信号Ｗ3を示
し、…、図８（ｇ）は第五の画像メモリ１６ｅに入力す
る書き込み許可信号Ｗ5を示している。これらの書き込
み許可信号Ｗ1〜Ｗ5は、図７に示す制御回路１９から各
画像メモリ１６ａ〜１６ｅに送出されるメモリ制御信号
Ｓ1に含まれる。この図８（ａ）〜（ｇ）からわかるよ
うに、各画像メモリ１６ａ〜１６ｅは、１フレーム毎に
経時的にデータが書き込まれて行く。この場合、画像メ
モリ１６ａ〜１６ｅは例えば５フレームとされており、
５フレーム毎に各画像メモリ１６ａ〜１６ｅが書き込み
許可される。
【００２５】次に、上記各画像メモリ１６ａ〜１６ｅか
ら読み出されるデータは以下のようになる。まず、図８
（ｈ）に示すように、第一の画像メモリ１６ａにフレー
ム番号ｆ1のデータが書き込まれてから５フレームの間
は、そのデータｆ1が出力されつづける。同様に、図８
（ｉ）に示すように、第二の画像メモリ１６ｂにはフレ
ーム番号ｆ2のデータが書き込まれてから５フレームの
間はそのデータｆ2が出力されつづけ、図８（ｊ）に示
すように、第三の画像メモリ１６ｃにはフレーム番号ｆ
3のデータが書き込まれてから５フレームの間はそのデ
ータｆ3が出力されつづけ、…、図８（ｌ）に示すよう
に、第五の画像メモリ１６ｅにはフレーム番号ｆ5のデ
ータが書き込まれてから５フレームの間はそのデータｆ
5が出力されつづける。なお、図８（ｍ）は、各画像メ
モリ１６ａ〜１６ｅのアドレス信号を示しており、各フ
レームとも共通であり、符号Ａdで示すデータが毎フレ
ーム繰り返されている。このアドレス信号Ａdも、図７
に示すメモリ制御信号Ｓ1に含まれている。
ら読み出されるデータは以下のようになる。まず、図８
（ｈ）に示すように、第一の画像メモリ１６ａにフレー
ム番号ｆ1のデータが書き込まれてから５フレームの間
は、そのデータｆ1が出力されつづける。同様に、図８
（ｉ）に示すように、第二の画像メモリ１６ｂにはフレ
ーム番号ｆ2のデータが書き込まれてから５フレームの
間はそのデータｆ2が出力されつづけ、図８（ｊ）に示
すように、第三の画像メモリ１６ｃにはフレーム番号ｆ
3のデータが書き込まれてから５フレームの間はそのデ
ータｆ3が出力されつづけ、…、図８（ｌ）に示すよう
に、第五の画像メモリ１６ｅにはフレーム番号ｆ5のデ
ータが書き込まれてから５フレームの間はそのデータｆ
5が出力されつづける。なお、図８（ｍ）は、各画像メ
モリ１６ａ〜１６ｅのアドレス信号を示しており、各フ
レームとも共通であり、符号Ａdで示すデータが毎フレ
ーム繰り返されている。このアドレス信号Ａdも、図７
に示すメモリ制御信号Ｓ1に含まれている。
【００２６】次に、図８（ｂ）のフレームデータがｆ5
になった時間に注目すると、図８（ｈ）〜（ｌ）に示す
各画像メモリ１６ａ〜１６ｅの出力としてフレームデー
タｆ1〜ｆ5が揃っていることがわかる。そして、これら
のデータｆ1〜ｆ5が図７に示す画像メモリ１６ａ〜１６
ｅから同時に平均輝度算出回路２０へ送出される。この
とき、制御回路１９からは平均輝度算出回路制御信号Ｓ
2が入力して、内部の加算平均回路のタイミングを制御
する。これにより、上記平均輝度算出回路２０は、画像
メモリ１６ａ〜１６ｅからのデータが入力される毎に平
均輝度Ｐを算出する。次に、個別分散回路２１ａ〜２１
ｅは、同様に上記画像メモリ１６ａ〜１６ｅからデータ
が出力される毎に、この画像メモリ１６ａ〜１６ｅから
のデータと上記平均輝度算出回路２０からの平均輝度Ｐ
とを入力して両者の差を計算し、その絶対値を求める。
その後、平均分散算出回路２２は、上記個別分散回路２
１ａ〜２１ｅからの個別分散値が入力する毎にそれらを
加算平均して、平均分散Ｑを算出する。このとき、制御
回路１９からは個別分散回路制御信号Ｓ3及び平均分散
算出回路制御信号Ｓ4がそれぞれ該当する回路に出力さ
れ、それぞれの内部回路のタイミングを制御する。
になった時間に注目すると、図８（ｈ）〜（ｌ）に示す
各画像メモリ１６ａ〜１６ｅの出力としてフレームデー
タｆ1〜ｆ5が揃っていることがわかる。そして、これら
のデータｆ1〜ｆ5が図７に示す画像メモリ１６ａ〜１６
ｅから同時に平均輝度算出回路２０へ送出される。この
とき、制御回路１９からは平均輝度算出回路制御信号Ｓ
2が入力して、内部の加算平均回路のタイミングを制御
する。これにより、上記平均輝度算出回路２０は、画像
メモリ１６ａ〜１６ｅからのデータが入力される毎に平
均輝度Ｐを算出する。次に、個別分散回路２１ａ〜２１
ｅは、同様に上記画像メモリ１６ａ〜１６ｅからデータ
が出力される毎に、この画像メモリ１６ａ〜１６ｅから
のデータと上記平均輝度算出回路２０からの平均輝度Ｐ
とを入力して両者の差を計算し、その絶対値を求める。
その後、平均分散算出回路２２は、上記個別分散回路２
１ａ〜２１ｅからの個別分散値が入力する毎にそれらを
加算平均して、平均分散Ｑを算出する。このとき、制御
回路１９からは個別分散回路制御信号Ｓ3及び平均分散
算出回路制御信号Ｓ4がそれぞれ該当する回路に出力さ
れ、それぞれの内部回路のタイミングを制御する。
【００２７】次に、図７において、スペックルノイズ判
定演算回路２３は、上記平均輝度算出回路２０からの平
均輝度Ｐと平均分散算出回路２２からの平均分散Ｑとを
入力する毎に、これらの比率を前述の式（３）によって
演算し（Ｓ＝Ｐ／Ｑ）、この判定関数Ｓが所定値より大
きいか否かでスペックルノイズであるかどうか判定す
る。その後、スペックルノイズ除去処理回路２４は、上
記スペックルノイズ判定演算回路２３からの判定結果Ｒ
を入力して超音波本体部１１から入力する現時点の画像
データについてスペックルノイズ除去の画像処理を行
う。このとき、制御回路１９からはスペックルノイズ判
定演算回路制御信号Ｓ5及びスペックルノイズ除去処理
回路制御信号Ｓ6がそれぞれ該当する回路に出力され、
それぞれの内部回路のタイミングを制御する。なお、図
７及び図８においては、画像メモリを５個とし５フレー
ムのデータ処理としたが、これに限らず、例えば８個の
画像メモリとし８フレームのデータ処理としてもよい。
定演算回路２３は、上記平均輝度算出回路２０からの平
均輝度Ｐと平均分散算出回路２２からの平均分散Ｑとを
入力する毎に、これらの比率を前述の式（３）によって
演算し（Ｓ＝Ｐ／Ｑ）、この判定関数Ｓが所定値より大
きいか否かでスペックルノイズであるかどうか判定す
る。その後、スペックルノイズ除去処理回路２４は、上
記スペックルノイズ判定演算回路２３からの判定結果Ｒ
を入力して超音波本体部１１から入力する現時点の画像
データについてスペックルノイズ除去の画像処理を行
う。このとき、制御回路１９からはスペックルノイズ判
定演算回路制御信号Ｓ5及びスペックルノイズ除去処理
回路制御信号Ｓ6がそれぞれ該当する回路に出力され、
それぞれの内部回路のタイミングを制御する。なお、図
７及び図８においては、画像メモリを５個とし５フレー
ムのデータ処理としたが、これに限らず、例えば８個の
画像メモリとし８フレームのデータ処理としてもよい。
【００２８】図９は第二の発明の超音波診断装置におけ
るスペックルノイズ判定除去回路１４の他の設置例を示
す全体構成のブロック図である。この設置例は、ＤＳＣ
部１２と画像表示部１３との間にスペックルノイズ判定
除去回路１４を設けたものである。図５の設置例では、
超音波本体部１１とＤＳＣ部１２との間にスペックルノ
イズ判定除去回路１４を設け、上記超音波本体部１１よ
り得たデータを超音波信号の送受信タイミングに同期し
て演算処理していたのに対し、図９の設置例では、ＤＳ
Ｃ部１２から出力されるフレーム画として完成されたデ
ータ、すなわちテレビモニタなどに表示するため走査変
換した画像データを画像表示部１３の表示のタイミング
に同期して演算処理するものである。この場合は、図６
に示すスペックルノイズ判定除去回路１４の内部構成に
おいて、複数個の画像メモリ１６ａ〜１６ｅへの書き込
み、読み出しのタイミングが走査変換後のタイミング、
すなわち画像表示部１３のテレビモニタなどへの表示タ
イミングに従っている。また、上記画像表示部１３の表
示用のフレーム画として上記複数個の画像メモリ１６ａ
〜１６ｅの画像データを使用している。
るスペックルノイズ判定除去回路１４の他の設置例を示
す全体構成のブロック図である。この設置例は、ＤＳＣ
部１２と画像表示部１３との間にスペックルノイズ判定
除去回路１４を設けたものである。図５の設置例では、
超音波本体部１１とＤＳＣ部１２との間にスペックルノ
イズ判定除去回路１４を設け、上記超音波本体部１１よ
り得たデータを超音波信号の送受信タイミングに同期し
て演算処理していたのに対し、図９の設置例では、ＤＳ
Ｃ部１２から出力されるフレーム画として完成されたデ
ータ、すなわちテレビモニタなどに表示するため走査変
換した画像データを画像表示部１３の表示のタイミング
に同期して演算処理するものである。この場合は、図６
に示すスペックルノイズ判定除去回路１４の内部構成に
おいて、複数個の画像メモリ１６ａ〜１６ｅへの書き込
み、読み出しのタイミングが走査変換後のタイミング、
すなわち画像表示部１３のテレビモニタなどへの表示タ
イミングに従っている。また、上記画像表示部１３の表
示用のフレーム画として上記複数個の画像メモリ１６ａ
〜１６ｅの画像データを使用している。
【００２９】
【発明の効果】第一の発明によるスペックルノイズ判定
方法は以上のように構成されたので、超音波診断装置の
超音波本体部のメモリ部から連続した複数フレームの画
像データを並列に読み出し、これら複数の画像データを
統計処理することにより、対象画像データがスペックル
ノイズであるか否かを判定することができる。特に、上
記複数の画像データの統計処理として、複数フレームの
画像データについて同一番地のデータの平均輝度と平均
分散とを求め、この平均輝度と平均分散との比率を演算
することにより、超音波画像に混入しているスペックル
ノイズを有効に判定することができる。
方法は以上のように構成されたので、超音波診断装置の
超音波本体部のメモリ部から連続した複数フレームの画
像データを並列に読み出し、これら複数の画像データを
統計処理することにより、対象画像データがスペックル
ノイズであるか否かを判定することができる。特に、上
記複数の画像データの統計処理として、複数フレームの
画像データについて同一番地のデータの平均輝度と平均
分散とを求め、この平均輝度と平均分散との比率を演算
することにより、超音波画像に混入しているスペックル
ノイズを有効に判定することができる。
【００３０】また、第二の発明による超音波診断装置
は、超音波本体部又はディジタルスキャンコンバータ部
から出力される連続した複数フレームの画像データを順
次書き込む複数個並列の画像メモリと、この画像メモリ
から読み出した複数の画像データを取り込みそれらの画
像データを統計処理することにより対象画像データがス
ペックルノイズであるか否か判定すると共にそのスペッ
クルノイズを除去する演算処理回路と、上記画像メモリ
及び演算処理回路の動作を制御する制御回路とを備え、
上記演算処理回路から前記ディジタルスキャンコンバー
タ部又は画像表示部へ処理後の画像データを送出するス
ペックルノイズ判定除去回路を備えたことにより、超音
波画像に混入しているスペックルノイズのみを除去する
ことができる。このスペックルノイズの除去により、画
質を向上することができると共に画質の安定性を増すこ
とができる。
は、超音波本体部又はディジタルスキャンコンバータ部
から出力される連続した複数フレームの画像データを順
次書き込む複数個並列の画像メモリと、この画像メモリ
から読み出した複数の画像データを取り込みそれらの画
像データを統計処理することにより対象画像データがス
ペックルノイズであるか否か判定すると共にそのスペッ
クルノイズを除去する演算処理回路と、上記画像メモリ
及び演算処理回路の動作を制御する制御回路とを備え、
上記演算処理回路から前記ディジタルスキャンコンバー
タ部又は画像表示部へ処理後の画像データを送出するス
ペックルノイズ判定除去回路を備えたことにより、超音
波画像に混入しているスペックルノイズのみを除去する
ことができる。このスペックルノイズの除去により、画
質を向上することができると共に画質の安定性を増すこ
とができる。
【００３１】また、上記スペックルノイズ判定除去回路
の演算処理回路として、上記画像メモリから読み出した
複数フレームの画像データについて同一番地のデータの
平均輝度を求める平均輝度算出回路と、上記複数フレー
ムの画像データについて同一番地のデータの平均分散を
求める平均分散算出回路と、上記それぞれの回路で求め
た算出結果を入力して平均輝度と平均分散との比率を演
算し対象画像データがスペックルノイズであるか否かを
判定するスペックルノイズ判定演算回路と、この判定結
果を入力して前記超音波本体部又はディジタルスキャン
コンバータ部からの画像データについてスペックルノイ
ズを除去するスペックルノイズ除去処理回路とを有して
成るものとした場合は、例えば画像メモリが８個で８ビ
ットのデータを８フレーム使用して統計処理するときで
も、スペックルノイズ判定除去回路を小形かつ低コスト
で実現することができる。
の演算処理回路として、上記画像メモリから読み出した
複数フレームの画像データについて同一番地のデータの
平均輝度を求める平均輝度算出回路と、上記複数フレー
ムの画像データについて同一番地のデータの平均分散を
求める平均分散算出回路と、上記それぞれの回路で求め
た算出結果を入力して平均輝度と平均分散との比率を演
算し対象画像データがスペックルノイズであるか否かを
判定するスペックルノイズ判定演算回路と、この判定結
果を入力して前記超音波本体部又はディジタルスキャン
コンバータ部からの画像データについてスペックルノイ
ズを除去するスペックルノイズ除去処理回路とを有して
成るものとした場合は、例えば画像メモリが８個で８ビ
ットのデータを８フレーム使用して統計処理するときで
も、スペックルノイズ判定除去回路を小形かつ低コスト
で実現することができる。
【図１】第一の発明によるスペックルノイズ判定方法の
原理を示す説明図である。
原理を示す説明図である。
【図２】被検体内の診断部位について収集した超音波断
層像におけるスペックルノイズの分布状況を示す説明図
である。
層像におけるスペックルノイズの分布状況を示す説明図
である。
【図３】収集した超音波画像の時間変化に対する輝度変
化の状況を示すグラフである。
化の状況を示すグラフである。
【図４】スペックルノイズであるか否かの判定を複数フ
レームのデータをグラフの形状で分類して行う場合を説
明するグラフである。
レームのデータをグラフの形状で分類して行う場合を説
明するグラフである。
【図５】第二の発明によるスペックルノイズ判定除去回
路を備えた超音波診断装置の全体構成を示すブロック図
である。
路を備えた超音波診断装置の全体構成を示すブロック図
である。
【図６】スペックルノイズ判定除去回路の内部構成を示
すブロック図である。
すブロック図である。
【図７】上記スペックルノイズ判定除去回路における演
算処理回路の具体的な内部構成例を示すブロック図であ
る。
算処理回路の具体的な内部構成例を示すブロック図であ
る。
【図８】上記演算処理回路の動作を説明するためのタイ
ミング線図である。
ミング線図である。
【図９】第二の発明の超音波診断装置におけるスペック
ルノイズ判定除去回路の他の設置例を示す全体構成のブ
ロック図である。
ルノイズ判定除去回路の他の設置例を示す全体構成のブ
ロック図である。
【図１０】従来の超音波診断装置においてスペックルノ
イズを除去するためのフレーム相関処理回路を示すブロ
ック図である。
イズを除去するためのフレーム相関処理回路を示すブロ
ック図である。
１０…探触子 １１…超音波本体部 １２…ＤＳＣ部 １３…画像表示部 １４，１４′…スペックルノイズ判定除去回路 １５…データ分配器 １６ａ〜１６ｅ…画像メモリ １７…データ選別器 １８…演算処理回路 １９…制御回路 ２０…平均輝度算出回路 ２１ａ〜２１ｅ…個別分散回路 ２２…平均分散算出回路 ２３…スペックルノイズ判定演算回路 ２４…スペックルノイズ除去処理回路
Claims (4)
- 【請求項１】 被検体内に超音波を送受信する探触子
と、この探触子を駆動して超音波を発生させると共に受
信した反射エコー信号を処理する超音波送受信部及びこ
の超音波送受信部からの反射エコー信号をディジタル化
し運動組織を含む被検体内の超音波画像のデータを時系
列に複数フレームに記録するメモリ部を備えた超音波本
体部と、上記メモリ部からのディジタル信号を超音波ビ
ームの走査線毎に書き込んで画像データを形成するディ
ジタルスキャンコンバータ部と、このディジタルスキャ
ンコンバータ部からの画像データをアナログ変換し画像
として表示する画像表示部とを有する超音波診断装置に
おいて、上記超音波本体部のメモリ部から連続した複数
フレームの画像データを並列に読み出し、これら複数の
画像データを統計処理することにより、対象画像データ
がスペックルノイズであるか否かを判定することを特徴
とするスペックルノイズ判定方法。 - 【請求項２】 上記複数の画像データの統計処理は、複
数フレームの画像データについて同一番地のデータの平
均輝度と平均分散とを求め、この平均輝度と平均分散と
の比率を演算することを特徴とする請求項１記載の超音
波診断装置におけるスペックルノイズ判定方法。 - 【請求項３】 被検体内に超音波を送受信する探触子
と、この探触子を駆動して超音波を発生させると共に受
信した反射エコー信号を処理する超音波送受信部及びこ
の超音波送受信部からの反射エコー信号をディジタル化
し運動組織を含む被検体内の超音波画像のデータを時系
列に複数フレームに記録するメモリ部を備えた超音波本
体部と、上記メモリ部からのディジタル信号を超音波ビ
ームの走査線毎に書き込んで画像データを形成するディ
ジタルスキャンコンバータ部と、このディジタルスキャ
ンコンバータ部からの画像データをアナログ変換し画像
として表示する画像表示部とを有する超音波診断装置に
おいて、上記超音波本体部又はディジタルスキャンコン
バータ部から出力される連続した複数フレームの画像デ
ータを順次書き込む複数個並列の画像メモリと、この画
像メモリから読み出した複数の画像データを取り込みそ
れらの画像データを統計処理することにより対象画像デ
ータがスペックルノイズであるか否か判定すると共にそ
のスペックルノイズを除去する演算処理回路と、上記画
像メモリ及び演算処理回路の動作を制御する制御回路と
を備え、上記演算処理回路から前記ディジタルスキャン
コンバータ部又は画像表示部へ処理後の画像データを送
出するスペックルノイズ判定除去回路を備えたことを特
徴とする超音波診断装置。 - 【請求項４】 上記スペックルノイズ判定除去回路の演
算処理回路は、上記画像メモリから読み出した複数フレ
ームの画像データについて同一番地のデータの平均輝度
を求める平均輝度算出回路と、上記複数フレームの画像
データについて同一番地のデータの平均分散を求める平
均分散算出回路と、上記それぞれの回路で求めた算出結
果を入力して平均輝度と平均分散との比率を演算し対象
画像データがスペックルノイズであるか否かを判定する
スペックルノイズ判定演算回路と、この判定結果を入力
して前記超音波本体部又はディジタルスキャンコンバー
タ部からの画像データについてスペックルノイズを除去
するスペックルノイズ除去処理回路とを有して成るもの
であることを特徴とする請求項３記載の超音波診断装
置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27510795A JP3542858B2 (ja) | 1995-09-29 | 1995-09-29 | 超音波診断装置 |
| US08/721,715 US5664572A (en) | 1995-09-29 | 1996-09-27 | Method for discriminating speckle noises from signals in ultrasonic tomography apparatus and ultrasonic tomography apparatus including speckle noise removing circuit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27510795A JP3542858B2 (ja) | 1995-09-29 | 1995-09-29 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0994248A true JPH0994248A (ja) | 1997-04-08 |
| JP3542858B2 JP3542858B2 (ja) | 2004-07-14 |
Family
ID=17550847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP27510795A Expired - Fee Related JP3542858B2 (ja) | 1995-09-29 | 1995-09-29 | 超音波診断装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5664572A (ja) |
| JP (1) | JP3542858B2 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007135884A1 (ja) | 2006-05-19 | 2007-11-29 | Hitachi Medical Corporation | 超音波診断装置及び超音波診断方法 |
| JP2009090104A (ja) * | 2007-09-18 | 2009-04-30 | Fujifilm Corp | 超音波診断方法及び装置 |
| JP2009095665A (ja) * | 2007-09-28 | 2009-05-07 | Fujifilm Corp | 超音波診断方法及び装置 |
| EP2169425A2 (en) | 2008-09-25 | 2010-03-31 | Fujifilm Corporation | Ultrasound diagnostic apparatus and ultrasound diagnostic method |
| JP2014028283A (ja) * | 2007-09-18 | 2014-02-13 | Fujifilm Corp | 超音波診断方法及び装置 |
| JP2014210121A (ja) * | 2013-04-22 | 2014-11-13 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置及びその制御プログラム |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5876343A (en) * | 1997-09-23 | 1999-03-02 | Scimed Life Systems, Inc. | Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system |
| JP3242887B2 (ja) * | 1998-09-21 | 2001-12-25 | 松下電器産業株式会社 | 超音波診断装置 |
| JP3330092B2 (ja) * | 1998-09-30 | 2002-09-30 | 松下電器産業株式会社 | 超音波診断装置 |
| US6484191B1 (en) | 1999-07-02 | 2002-11-19 | Aloka Co., Ltd. | Apparatus and method for the real-time calculation of local variance in images |
| US6454715B2 (en) | 2000-04-11 | 2002-09-24 | Scimed Life Systems, Inc. | Methods and apparatus for blood speckle detection in an intravascular ultrasound imaging system |
| ITSV20020040A1 (it) * | 2002-08-13 | 2004-02-14 | Esaote Spa | Metodo e dispositivo per il rilevamento di immagini ecografiche |
| US20040077946A1 (en) * | 2002-10-15 | 2004-04-22 | Jun Ohmiya | Image processing apparatus, method and program |
| US8280482B2 (en) * | 2004-04-19 | 2012-10-02 | New York University | Method and apparatus for evaluating regional changes in three-dimensional tomographic images |
| JP4625663B2 (ja) * | 2004-08-18 | 2011-02-02 | 株式会社東芝 | 超音波診断装置 |
| WO2010125789A1 (ja) * | 2009-04-28 | 2010-11-04 | 株式会社日立メディコ | 超音波画像の画質改善方法、超音波診断装置、及び画質改善プログラム |
| JP2012182678A (ja) * | 2011-03-01 | 2012-09-20 | Jvc Kenwood Corp | 立体映像処理装置および立体映像処理方法 |
| GB2502997B (en) * | 2012-06-13 | 2014-09-03 | Crystalview Medical Imaging Ltd | Suppression of reverberations and/or clutter in ultrasonic imaging systems |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5477858A (en) * | 1986-07-30 | 1995-12-26 | Siemens Medical Systems, Inc. | Ultrasound blood flow/tissue imaging system |
| US5479926A (en) * | 1995-03-10 | 1996-01-02 | Acuson Corporation | Imaging system display processor |
-
1995
- 1995-09-29 JP JP27510795A patent/JP3542858B2/ja not_active Expired - Fee Related
-
1996
- 1996-09-27 US US08/721,715 patent/US5664572A/en not_active Expired - Lifetime
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007135884A1 (ja) | 2006-05-19 | 2007-11-29 | Hitachi Medical Corporation | 超音波診断装置及び超音波診断方法 |
| JP4879263B2 (ja) * | 2006-05-19 | 2012-02-22 | 株式会社日立メディコ | 超音波診断装置及び超音波診断方法 |
| US8300909B2 (en) | 2006-05-19 | 2012-10-30 | Hitachi Medical Corporation | Ultrasonographic device and ultrasonographic method |
| JP2009090104A (ja) * | 2007-09-18 | 2009-04-30 | Fujifilm Corp | 超音波診断方法及び装置 |
| JP2014028283A (ja) * | 2007-09-18 | 2014-02-13 | Fujifilm Corp | 超音波診断方法及び装置 |
| JP2009095665A (ja) * | 2007-09-28 | 2009-05-07 | Fujifilm Corp | 超音波診断方法及び装置 |
| EP2169425A2 (en) | 2008-09-25 | 2010-03-31 | Fujifilm Corporation | Ultrasound diagnostic apparatus and ultrasound diagnostic method |
| US8398548B2 (en) | 2008-09-25 | 2013-03-19 | Fujifilm Corporation | Ultrasound diagnostic apparatus and ultrasound diagnostic method |
| JP2014210121A (ja) * | 2013-04-22 | 2014-11-13 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置及びその制御プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| US5664572A (en) | 1997-09-09 |
| JP3542858B2 (ja) | 2004-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0994248A (ja) | 超音波診断装置におけるスペックルノイズ判定方法及びスペックルノイズ判定除去回路を備えた超音波診断装置 | |
| JP4817728B2 (ja) | 超音波診断装置 | |
| JP3732613B2 (ja) | 超音波診断装置 | |
| EP0406909A2 (en) | Ultrasonic diagnosing apparatus | |
| JP3878273B2 (ja) | 超音波診断装置 | |
| JP2000287972A (ja) | 超音波診断装置および超音波診断装置の画像生成方法 | |
| JPH04200457A (ja) | 超音波診断装置 | |
| JPH10243945A5 (ja) | ||
| JPS6214294B2 (ja) | ||
| JP2006130009A (ja) | 超音波診断装置 | |
| JP3889121B2 (ja) | 超音波診断装置 | |
| JPH11226013A (ja) | 超音波診断装置 | |
| JPH0975350A (ja) | 超音波診断装置 | |
| JPH06209942A (ja) | 超音波診断装置 | |
| JP3047478B2 (ja) | 超音波診断装置のディジタルスキャンコンバータ | |
| JP2000316846A (ja) | 超音波診断装置 | |
| JPH08266538A (ja) | 超音波診断装置 | |
| JP2939640B2 (ja) | 被検体内の血流速分布像を二次元表示する超音波診断装置 | |
| JP2011217826A (ja) | 超音波診断装置及び多重検出プログラム | |
| JPH0221262B2 (ja) | ||
| JPH08308840A (ja) | 超音波診断装置 | |
| JP2791578B2 (ja) | 超音波ドプラ装置 | |
| JPH0614924A (ja) | 超音波診断装置 | |
| JPH03297454A (ja) | 超音波診断装置 | |
| JPH06154215A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20031216 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040210 |
|
| 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: 20040309 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040402 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090409 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100409 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100409 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110409 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120409 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120409 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130409 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140409 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |