JPH02164352A - 超音波診断画像処理装置 - Google Patents

超音波診断画像処理装置

Info

Publication number
JPH02164352A
JPH02164352A JP63317918A JP31791888A JPH02164352A JP H02164352 A JPH02164352 A JP H02164352A JP 63317918 A JP63317918 A JP 63317918A JP 31791888 A JP31791888 A JP 31791888A JP H02164352 A JPH02164352 A JP H02164352A
Authority
JP
Japan
Prior art keywords
value
pixel
calculated
dispersion value
picture element
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
Application number
JP63317918A
Other languages
English (en)
Other versions
JPH0751134B2 (ja
Inventor
Masaru Murashita
賢 村下
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.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP63317918A priority Critical patent/JPH0751134B2/ja
Publication of JPH02164352A publication Critical patent/JPH02164352A/ja
Publication of JPH0751134B2 publication Critical patent/JPH0751134B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 [産業上の利用分野] 本発明は超音波診断画像処理装置、特に超音波を被検体
内に送受波し反射エコーから得られた情報を画像表示す
るためにデジタル処理する装置の回路構成に関する。
[従来の技術] 超音波を生体などの被検体内に送受波して被検体の画像
診断を行う超音波診断装置が周知であり、この装置は、
例えば一定繰り返し周期で出力されるパルス超音波を被
検体内に送波し、被検体内から反射するエコー信号を検
波し、これを所定領域について行うことにより一画面分
の情報を得るものである。そして、この情報はCRT表
示器に表示され、この表示器により例えば胆嚢などの腹
部臓器や心臓などの断層像を観察することができる。
このような超音波診断装置で画像処理を行う装置では、
エコー信号をアナログ的に輪郭強調することが行われて
いる。すなわち、検波後のエコー信号に対して微分処理
を施し、この微分信号を元の信号に加算することにより
、エコー信号中の輪郭情報を強調しており、これによっ
て被検体内の形状や構造を明瞭に画像表示することがで
きる。
[発明が解決しようとする課題] しかしながら、従来のような微分処理による輪郭強調操
作では、処理が一律に固定されており、自由度の高い処
理ができず、また微分処理であるため、一連の情報とな
る信号の立上り部は強調されるが、立下り部は強調され
ないという問題がある。
すなわち、第6図に示されるように、(a)のエコー信
号を微分処理する場合には、(b)のように立上り部は
輪郭信号100として検出されるが、立下り部は検出さ
れないことになる。従って、輪郭を正確に強調する場合
には、(b)の信号をもう一度微分して信号101の立
上り部を輪郭信号として検出しなければならない等の不
便がある。
また、このような微分処理は抽出したり情報以外の不必
要な信号、例えば雑音やスペックルなどの信号について
も行われることになり、不必要な信号まで強調してしま
う不都合がある。
更に、微分処理は画素単位に一次元的に行っているので
、全体的に見た場合には他の画素間でのつながりがなく
、二次元的な処理を行うことができないという問題があ
る。
発明の目的 本発明は前記問題点を解決することを課題としてなされ
たもので、その目的は、輪郭強調を精度良く行い、かつ
画素間での情報のつながりのある二次元的な処理を行う
ことができる超音波診断画像処理装置を提供することに
ある。
また、他の目的は、前記処理を実時間で行うことのでき
る超音波診断画像処理装置を提供することにある。
[課題を解決するための手段] 前記目的を達成するために、本発明は、注目画素を中心
とした所定領域の画素群を選択しこの画素群の画素平均
値を演算する平均値演算回路と、前記画素群の分散値を
演算する分散値演算回路と、この分散値演算回路から出
力された分散値が基準分散値よりも大きいか否かを判定
する比較器と、この比較器により分散値が大きいと判定
された場合には輪郭強調処理を行い注目画素値をその値
が前記画素平均値より大きい場合はより大きく設定し画
素平均値より小さい場合はより小さく設定する画素値強
調回路と、を含むことを特徴とする。
また、他の発明は、前記平均値演算回路及び分散値演算
回路はシフトレジスタを含み、バイブライン方式により
複数の処理を並列して順次行うことを特徴とする。
[作用] 以上の構成によれば、まず注目画素を中心とした、例え
ば縦横NXNの画素領域が選択され、この所定領域内の
画素平均値が演算される。そして、これと同時に分散値
が演算されるが、この分散値は各画素の信号の2乗平均
値z2を求め、この2乗平均値から前記画素平均値を2
乗(Σ)2を差し引(ことにより求められる。
そして、この分散値は比較器にて基準分散値と比較され
、分散値が基準分散値より大きい場合にのみ強調処理が
行われる。すなわち、この輪郭強調処理は、注目画素が
前記画素平均値より大きい場合には注目画素値をより大
きく設定し、画素平均値より小さい場合は注目画素値を
より小さく設定する。従って、診断画像において濃い部
分はより濃く、淡い部分はより淡くなり、輪郭が極めて
はっきりと表現されることになる。
また、この輪郭強調処理は分散値が基準分散値より小さ
い場合には行われず、そのままの画素値が出力されるの
で、ノイズやスペックルの信号は強調されないことにな
る。
[実施例] 以下、図面に基づいて本発明の好適な実施例を説明する
第1図には、実施例に係る超音波診断画像処理装置が示
され、第2図には、第1図装置を超音波診断装置に適用
した場合の概略構成が示されており、まず全体構成を説
明する。
第2図において、先端部に振動子が配列された探触子1
Gには送受信部12が接続されており、この送受信器1
2により超音波を被検体内に送信するとともに、その反
射エコーを受信する。例えば、断層像を画像表示する場
合には、受信した信号は検波され、信号振幅の減衰情報
を示す信号が形成される。
この送受信器12にはメモリ14が接続され、このメモ
リ14にはCRT表示器を含む表示部16が接続される
。前記メモリ14は、表示部16内のCRTの走査速度
に合せるために画像情報を一旦保持するもので、画素毎
の情報が格納されている。
本発明は、このメモリ14の画素情報に輪郭強調処理を
施すものであり、この輪郭強調処理を精度よくかつ二次
元的に行うための画像処理装置18に関するものである
第1図において、画像処理装置18は平均値演算回路、
分散値演算回路、比較器及び画素値強調回路から構成さ
れるが、本発明では、実時間処理を行うために複数の処
理を並列して順次行うバイブライン方式を採用する。
すなわち、減算処理する所定の領域の画素がNXN個で
あるとすると、メモリ14内の画素メモリをN列で出力
させることになり、まずN個のシフトレジスタ20と加
減算器22が設けられ、この加減算器22には加算器2
4及び除算器(十N)26が接続され、前記シフトレジ
スタ20から除算器26までの回路にて平均値演算回路
が構成される。
一方、前記シフトレジスタ20の他方の出力側にはN個
の2乗算器30が接続され、この2乗算器30の後段に
はN個のシフトレジスタ32及び加減算器34が接続さ
れ、更に加算器36及び除算器(+N)3gが接続され
ており、前記2乗算器30から除算器38までの回路に
て2乗平均値z2が演算される。また、この2乗平均値
演算と同時に前記平均値演算回路から出力された平均値
の2乗も演算されており、これは2乗算器28にて行わ
れる。なお、前記除算器24.38は、実施例では入力
情報のビットをシフトさせるだけで割り算を行うことが
できる。
そして、前記除算器38の後段には減算器40が設けら
れており、前記除算器38から出力された画素数の2乗
平均値から前記2乗算器28から出力された画素平均値
の2乗値を減算し、次式に示す分散値σを求める。
σ−z2− (H) 2        ・・・(1)
以上のようにして、所定領域の画素群の画素平均値と分
散値が求められることになるが、演算処理の対象となる
所定領域は、第3図に示されるように、実施例では画面
内においてNXN個で形成される正方形領域であり、そ
の中心位置にある画素が注目画素となり、この注目画素
の画素値決定処理を行うことになる。この場合の所定領
域の大きさは、3×3.5×5.7X7の画素群で形成
される大きさが好適である。
本発明は、処理演算をバイブライン方式にて行うことも
特徴としており、このために、シフトレジスタ20及び
32を設け、メモリ内の複数の情報を並列して順次処理
している。
すなわち、第4図には、N−5の場合のシフトレジスタ
20(32)と加減算器22(34)の動作が示されて
おり、メモリ14のメモリ1ラインから’11.112
・・・ilBの情報が、メモリ2ラインから’21”2
2・・・’2B・・・の情報が順次出力され、同じよう
にメモリ3、メモリ4、メモリ5のラインから画素情報
が出力される。そして、図示の状態では加減算器22(
34)はシフトレジスタ20(32)内の位置1〜5に
格納されている’11〜i  あるいは’ 21= ’
 25の画素値をそれぞれ15ゝ 加算することになり、これら加算値は更に第1図の加算
器24(36)にて加算され、所定位置の注目画素につ
いての処理が行われる。そして、次の注目画素について
処理を行う場合には、前記画素値を第4図の右側にシフ
トすることになり、画素値i とi は外され、’1B
や’2Bが新たに加わって加算演算が行われることにな
る。従って、シフトレジスタ20(32)及び加減算器
22(34)においては、メモリ14から順次画素情報
を読み出し画素値を1個ずつずらしながら単純に加算す
るだけで、短時間に多数の情報を加算処理することが可
能となる。
なお、バイブライン方式にて複数の処理を行うため既に
演算された値を所定時間だけ保持するラッチ回路42が
各ステージに設けられており、ステージ1においてメモ
リ14から得られた注目画素の画素値2はラッチ回路4
2a〜42eによりラッチされる。
次に、前記減算器40には比較器44が接続されており
、この比較器44は演算された分散値を基準分散値σ 
と比較し、この基準分散値よりも大きい場合にのみ後の
処理を行うようにしている。
すなわち、本発明では輪郭強調すべき情報を選択してお
り、ノイズやスペックルにより生じる信号成分について
は輪郭強調処理を除外し、輪郭強調の精度を高めている
。従って、分散値が基準分散値よりも小さい場合には注
目画素値をそのまま出力する。
そして、分散値が基準分散値よりも大きくなり輪郭強調
を行う場合には、まず所定の係数値Kを次式にて求める
K−kEX  (1−(σ /σ)) ・・・ (σ〉
σ )Cc 鴎0                ・・・ (σ≦
σ )但し、kE:輪郭強調定数     ・・・(2
)すなわち、比較器44には除算器48、減算器50及
び定数掛算器52が設けられ、除算器48にて前記(2
)式の(σ /σ)が演算され、減算器50にて(1−
(σ /σ))が演算され、最後に乗算器52により定
数kEが乗算されて係数Kが求められる。
そして、次段の乗算器54により、KX (z −2)
を演算する。すなわち、除算器26の後段には減算器4
6が接続されており、この減算器46はラッチ42aに
てラッチされている注目画素値から除算器26にて出力
される画素平均値を差し引いて、z−Eを演算し、ラッ
チ回路42f、42gを介して乗算器54に出力する。
従って、乗算器54ではKX (Z−E)が演算される
。この場合、(z−E)の値は正値と負値の両者があり
、正値の場合は画素値を増加させる方向に、負値の場合
は画素値を減少させる方向に作用することになり、これ
により輪郭強調が良好に行われる。
この乗算器54の後段には加算器56が接続され、この
加算器56ではラッチ42eの出力を加算して次式が演
算される。
Z−z+KX (z−z)       −(3)なお
、この加算器56の出力はラッチ回路42hを介して注
目画素の値としてメモリ14に供給する。
実施例は以上の構成からなり、以下にその作用を説明す
る。
まず、第2図の送受信部12の制御により探触子10か
ら超音波が被検体内に送信され、被検体内から反射され
るエコー信号が受信される。そして、受信信号は送信部
12内で検波され、被検体内の走査ライン順の情報がメ
モリ14に格納される。
このようにしてメモリ14内に格納された情報は、画像
処理装置18内のシフトレジスタ2Gに順次供給される
ことになり、第1図に示される加減算器22から除算器
26までの回路にて画素平均値iが演算される。この画
素平均値Σは2乗算器28にて2乗される。
一方、シフトレジスタ20の出力は2乗算器30に供給
されており、シフトレジスタ32から除算器38までの
回路を介して2乗平均値z2が演算され、減算器40で
は除算器38の出力から2乗算器28の出力を引いて分
散値σ−22−(Σ)2を演算する。
次に、比較器44において、前記分散値σが基準分散値
σ より大きいか否かが判断され、σ〉σ である場合
には前記(2)式に従って係数Kが演算され、σ≦σ 
である場合には係数にはOとなる。
従って、分散値が基準値より大きい場合にのみ輪郭強調
処理が行われることになり、除算器48から加算器56
までの回路及び減算器46等にて前記(3)式が演算さ
れることになる。
この場合、除算器46の出力z−Eが正値なら前記(3
)式の最終的な画素値は大きな値になり、z−zが負値
なら最終的な画素値はより小さな値となる。すなわち、
第5図(b)に示されるように、注目画素値が画素平均
値よりも大きい場合にはより大きく強調され、注目画素
値が画素平均値より小さい場合にはより小さく強調され
ることになり、これにより輪郭強調が精度良く7行われ
る。
このようにして、本発明ではシフトレジスタ20(32
)を用いパイプライン方式にて処理するので、輪郭強調
処理を実時間で行うことができる。
すなわち、通常、メモリ14内に存在する二次元情報を
表示部に表示させるにはCRTの走査線1本分83.5
μSecの間に600個前後のデータをメモリ14から
読み出すことになる。このことは、一画素の情報を読み
出してCRTに出力するのに、約1000μsecかか
ることを意味しており、この1000μsecの間にメ
モリ14からNxN個の情報を読出してさらに処理を施
すということは不可能である。従って、本発明のように
メモリ14内から画素情報を複数個並列に出力し、同時
に処理すれば、実時間で精度の良い輪郭強調処理を行う
ことが可能であることが理解される。
前記実施例で説明した輪郭強調の際の係数K(−kE)
は、注目画素値が画素平均値よりも大きいか否かで異な
る値に設定することも良い。
例えば、z>Hの場合はに−4、z<zの場合はに−2
というように、濃度の高い方をより強調することにより
、より鮮明な画像を得ることが可能である。
[発明の効果] 以上説明したように、本発明によれば、所定領域の画素
群の平均値及び分散値を求め、分散値が基準分散値より
も大きい場合に輪郭強調処理を行い、注目画素値をその
値が画素平均値より大きい場合はより大きくし、画素平
均値よりも小さい場合はより小さく設定するようにした
ので、各部の濃淡差をはっきり表現することができる。
また、基準分散値よりも大きい場合にのみ輪郭強調処理
を行うようにしたので、ノイズやスペックルの信号を強
調することがなく、濃淡差を精度良く求めることができ
る。
更に、所定領域の画素平均値を考慮して輪郭強調を行う
ので、−次元的処理ではなく、二次元的な処理を行うこ
とができ、精度の良い輪郭強調を行うことが可能である
また、他の発明によれば、本発明の輪郭強調処理をパイ
プライン方式にて行うようにしたので、多くの処理を短
時間に行い、実時間に処理することができる。
【図面の簡単な説明】
第1図は本実施例に係る超音波診断画像処理装置を示す
回路ブロック図、 第2図は本発明装置を超音波診断装置に適用した場合の
概略を示すブロック図、 第3図は1画面内において演算処理する所定領域を示す
説明図、 第4図はシフトレジスタ及び加減算器の作用を示す説明
図、 第5図は実施例装置において輪郭強調される状態を示す
波形図、 第6図は従来の装置において行われる微分処理を示す波
形図である。 10 ・・・ 探触子 14 ・・・ メモリ 16 ・・・ 表示部 18 ・・・ 画像処理装置 20、 2 ・・・ シフトレジスタ 22、 4 ・・・ 加減算器 24、  6.56  ・・・ 加算器26、  8.
48  ・・・ 除算器28、 0 ・・・ 2乗算器 40、  6.50  ・・・ 減算器44 ・・・ 
比較器 52.54 乗算器。

Claims (2)

    【特許請求の範囲】
  1. (1)注目画素を中心とした所定領域の画素群を選択し
    この画素群の画素平均値を演算する平均値演算回路と、
    前記画素群の分散値を演算する分散値演算回路と、この
    分散値演算回路から出力された分散値が基準分散値より
    も大きいか否かを判定する比較器と、この比較器により
    分散値が大きいと判定された時に輪郭強調処理を行い注
    目画素値をその値が前記画素平均値より大きい場合はよ
    り大きく設定し画素平均値より小さい場合はより小さく
    設定する画素値強調回路と、を含む超音波診断画像処理
    装置。
  2. (2)請求項(1)記載の装置において、前記平均値演
    算回路及び分散値演算回路はシフトレジスタを含み、パ
    イプライン方式により複数の処理を並列して順次行うこ
    とを特徴とする超音波診断画像処理装置。
JP63317918A 1988-12-16 1988-12-16 超音波診断画像処理装置 Expired - Fee Related JPH0751134B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP63317918A JPH0751134B2 (ja) 1988-12-16 1988-12-16 超音波診断画像処理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP63317918A JPH0751134B2 (ja) 1988-12-16 1988-12-16 超音波診断画像処理装置

Publications (2)

Publication Number Publication Date
JPH02164352A true JPH02164352A (ja) 1990-06-25
JPH0751134B2 JPH0751134B2 (ja) 1995-06-05

Family

ID=18093500

Family Applications (1)

Application Number Title Priority Date Filing Date
JP63317918A Expired - Fee Related JPH0751134B2 (ja) 1988-12-16 1988-12-16 超音波診断画像処理装置

Country Status (1)

Country Link
JP (1) JPH0751134B2 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04122359A (ja) * 1990-09-14 1992-04-22 Aloka Co Ltd 画像輪郭強調処理装置
JPH04141158A (ja) * 1990-10-02 1992-05-14 Aloka Co Ltd 超音波画像処理装置
JPH057579A (ja) * 1991-06-18 1993-01-19 Konica Corp 放射線画像の照射野領域抽出装置
JPH08206111A (ja) * 1995-02-06 1996-08-13 Ge Yokogawa Medical Syst Ltd 超音波画像表示方法および超音波診断装置
JPH09179962A (ja) * 1995-12-22 1997-07-11 Matsushita Electric Ind Co Ltd 画像処理装置
JP2000316851A (ja) * 1999-05-06 2000-11-21 Aloka Co Ltd 超音波診断装置
JP2002183709A (ja) * 2000-12-05 2002-06-28 Ge Medical Systems Global Technology Co Llc 画像処理方法および装置、記録媒体並びに画像撮影装置
US9817119B2 (en) 2011-05-25 2017-11-14 Hitachi, Ltd. Ultrasound image processing apparatus and program
JP2018072266A (ja) * 2016-11-02 2018-05-10 株式会社プロアシスト 超音波アレイセンサシステム

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04122359A (ja) * 1990-09-14 1992-04-22 Aloka Co Ltd 画像輪郭強調処理装置
JPH04141158A (ja) * 1990-10-02 1992-05-14 Aloka Co Ltd 超音波画像処理装置
JPH057579A (ja) * 1991-06-18 1993-01-19 Konica Corp 放射線画像の照射野領域抽出装置
JPH08206111A (ja) * 1995-02-06 1996-08-13 Ge Yokogawa Medical Syst Ltd 超音波画像表示方法および超音波診断装置
JPH09179962A (ja) * 1995-12-22 1997-07-11 Matsushita Electric Ind Co Ltd 画像処理装置
JP2000316851A (ja) * 1999-05-06 2000-11-21 Aloka Co Ltd 超音波診断装置
JP2002183709A (ja) * 2000-12-05 2002-06-28 Ge Medical Systems Global Technology Co Llc 画像処理方法および装置、記録媒体並びに画像撮影装置
US9817119B2 (en) 2011-05-25 2017-11-14 Hitachi, Ltd. Ultrasound image processing apparatus and program
JP2018072266A (ja) * 2016-11-02 2018-05-10 株式会社プロアシスト 超音波アレイセンサシステム

Also Published As

Publication number Publication date
JPH0751134B2 (ja) 1995-06-05

Similar Documents

Publication Publication Date Title
JP4722283B2 (ja) 連続データ獲得を用いた超音波フロー・イメージングにおける運動の可視化のための方法および装置
US5170792A (en) Adaptive tissue velocity compensation for ultrasonic Doppler imaging
JP3182479B2 (ja) 弾性計測装置
KR100742467B1 (ko) 촬상 시스템, 촬상 방법 및 혈액 움직임 촬상 시스템
JPH0824678B2 (ja) 超音波診断装置
US6402694B1 (en) Color imaging display method and apparatus for color doppler imaging system
JPH05220138A (ja) 超音波診断装置
JP3410821B2 (ja) 超音波診断装置
JPH02164352A (ja) 超音波診断画像処理装置
US5865752A (en) Method and apparatus for ultrasound imaging using normalized difference between successive frames
US6500125B1 (en) Ultrasound b/color priority threshold calculation
JP4562656B2 (ja) 超音波ドプラ血流測定装置
JP2001299765A (ja) 超音波画像表示方法および超音波診断装置
US6045504A (en) Method and apparatus for polynomial approximation of nonlinear operations in medical ultrasound imaging
KR0178095B1 (ko) 초음파 칼라 도플러 영상시스템을 위한 필터링방법
JP3000823B2 (ja) 超音波ドプラ血流計
JP2723458B2 (ja) 超音波ドプラ診断装置
JP2563656B2 (ja) 超音波ドプラ映像装置
JPH04336056A (ja) 超音波ドプラ映像装置
JP2565917B2 (ja) 画像処理表示装置
JPH0523334A (ja) 超音波ドプラ映像装置
JPS6274344A (ja) 超音波診断装置
JPH08173427A (ja) 超音波診断装置
JPH0556972A (ja) 超音波血流イメージング装置
JPH0966055A (ja) 超音波診断装置

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees