JP2723453B2 - 超音波血流観測装置 - Google Patents
超音波血流観測装置Info
- Publication number
- JP2723453B2 JP2723453B2 JP5255953A JP25595393A JP2723453B2 JP 2723453 B2 JP2723453 B2 JP 2723453B2 JP 5255953 A JP5255953 A JP 5255953A JP 25595393 A JP25595393 A JP 25595393A JP 2723453 B2 JP2723453 B2 JP 2723453B2
- Authority
- JP
- Japan
- Prior art keywords
- pixel
- window
- pixel value
- blood flow
- value
- 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
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、超音波パルスドプラ法
に基づいて生体内の血流情報を表示する超音波血流観測
装置に関し、特に前記血流情報の信号処理による血流情
報の高精度化に関する。
に基づいて生体内の血流情報を表示する超音波血流観測
装置に関し、特に前記血流情報の信号処理による血流情
報の高精度化に関する。
【０００２】
【従来の技術】医用分野で用いられる超音波血流観測装
置は、生体内の臓器を白黒のＢモード断層像で実時間表
示すると共に、生体内に流れている血流の情報、つまり
血流の平均速度又は血流の速度分散又はドプラ信号の強
度をＢモード断層像に重ねてカラーで実時間表示する。
このため、このような装置の信号処理回路は、Ｂモード
像を得るための断層データ処理系と血流情報を得るため
のドプラ信号処理系とを有する。このドプラ信号処理系
において、血流情報を得るためのドプラ信号処理におけ
る血流情報は、超音波パルスドプラ法に基づいて血球か
らのエコー信号のドプラシフト量を計測することによっ
て得られる。
置は、生体内の臓器を白黒のＢモード断層像で実時間表
示すると共に、生体内に流れている血流の情報、つまり
血流の平均速度又は血流の速度分散又はドプラ信号の強
度をＢモード断層像に重ねてカラーで実時間表示する。
このため、このような装置の信号処理回路は、Ｂモード
像を得るための断層データ処理系と血流情報を得るため
のドプラ信号処理系とを有する。このドプラ信号処理系
において、血流情報を得るためのドプラ信号処理におけ
る血流情報は、超音波パルスドプラ法に基づいて血球か
らのエコー信号のドプラシフト量を計測することによっ
て得られる。
【０００３】超音波パルスドプラ法は、生体に対してス
キャンしている超音波パルスビームの各ビーム方向につ
いて超音波パルスを所定の周期で複数回生体内に順次送
り、そして各々その生体内から返ってくるエコー信号を
受信し、パルス周期にわたるエコー信号の位相変化（ド
プラシフト周波数）又は相関関数などを算出し、血流情
報へと変換する方法である。そしてこの方法による血流
情報の抽出をスキャンしているすべての超音波ビーム方
向について行えば、二次元の血流情報が得られ、画像と
して表示できる。
キャンしている超音波パルスビームの各ビーム方向につ
いて超音波パルスを所定の周期で複数回生体内に順次送
り、そして各々その生体内から返ってくるエコー信号を
受信し、パルス周期にわたるエコー信号の位相変化（ド
プラシフト周波数）又は相関関数などを算出し、血流情
報へと変換する方法である。そしてこの方法による血流
情報の抽出をスキャンしているすべての超音波ビーム方
向について行えば、二次元の血流情報が得られ、画像と
して表示できる。
【０００４】パルスドプラ法により順次受信されたエコ
ー信号には、血流からのドプラシフトを受けたドプラ成
分と、血流以外の静止又は動きの少ない臓器からのクラ
ッタ成分と、ノイズとが含まれている。エコー信号にお
けるクラッタ成分を抑圧し、ノイズを除去し、ドプラ成
分のみから血流情報を精度良く抽出するために、従来の
超音波血流情報装置ではドプラ信号処理系は図４に示す
ように構成されている。以下その原理を説明する。
ー信号には、血流からのドプラシフトを受けたドプラ成
分と、血流以外の静止又は動きの少ない臓器からのクラ
ッタ成分と、ノイズとが含まれている。エコー信号にお
けるクラッタ成分を抑圧し、ノイズを除去し、ドプラ成
分のみから血流情報を精度良く抽出するために、従来の
超音波血流情報装置ではドプラ信号処理系は図４に示す
ように構成されている。以下その原理を説明する。
【０００５】まず、探触子４１及び送受信回路４２によ
り、超音波パルスが生体内に送信され、生体内から返っ
てくるエコー信号が受信される。受信回路で所定の処理
を受けたエコー信号は位相検波回路４３に入力される。
位相検波されたエコー信号がＡ／Ｄコンバータ４４でデ
ジタル信号に変換された後にＭＴＩ（ＭｏｖｉｎｇＴａ
ｒｇｅｔ Ｉｎｄｉｃａｔｉｏｎ）フィルタ４５に通さ
れ、ここでクラッタ成分が抑圧される。クラッタ成分が
抑圧されたエコー信号は、相関処理回路４６で解析され
る。この解析では、パルス周期にわたるエコー信号のド
プラシフト周波数又は相関係数などが算出され、これに
より血流の平均速度又は血流の速度分散又はドプラ信号
の強度の血流情報が得られる。そして、この血流情報
は、ノイズ判別処理回路４７に入力され、ここでノイズ
が判別されノイズが除去される。ノイズが除去された血
流情報は、ＤＳＣ（Ｄｉｇｉｔａｌ Ｓｃａｎ Ｃｏｎ
ｖｅｒｔｅｒ）を含む画像表示回路４８に入力され、こ
こでカラーコーディングされた後に、カラー画像として
表示器４９に表示される。
り、超音波パルスが生体内に送信され、生体内から返っ
てくるエコー信号が受信される。受信回路で所定の処理
を受けたエコー信号は位相検波回路４３に入力される。
位相検波されたエコー信号がＡ／Ｄコンバータ４４でデ
ジタル信号に変換された後にＭＴＩ（ＭｏｖｉｎｇＴａ
ｒｇｅｔ Ｉｎｄｉｃａｔｉｏｎ）フィルタ４５に通さ
れ、ここでクラッタ成分が抑圧される。クラッタ成分が
抑圧されたエコー信号は、相関処理回路４６で解析され
る。この解析では、パルス周期にわたるエコー信号のド
プラシフト周波数又は相関係数などが算出され、これに
より血流の平均速度又は血流の速度分散又はドプラ信号
の強度の血流情報が得られる。そして、この血流情報
は、ノイズ判別処理回路４７に入力され、ここでノイズ
が判別されノイズが除去される。ノイズが除去された血
流情報は、ＤＳＣ（Ｄｉｇｉｔａｌ Ｓｃａｎ Ｃｏｎ
ｖｅｒｔｅｒ）を含む画像表示回路４８に入力され、こ
こでカラーコーディングされた後に、カラー画像として
表示器４９に表示される。
【０００６】
【発明が解決しようとする課題】上記のように、探触子
４１及び送受信回路４２により受信されたエコー信号に
は、ドプラ成分、クラッタ成分、ノイズ成分が含まれて
おり、高域通過フィルタであるＭＴＩフィルタ４５によ
り、低周波成分であるクラッタ成分は抑圧される。しか
し、完全に抑圧しきれない場合、残存クラッタ成分及び
前記ノイズ成分は、ノイズ判別処理回路４７に入力され
る。
４１及び送受信回路４２により受信されたエコー信号に
は、ドプラ成分、クラッタ成分、ノイズ成分が含まれて
おり、高域通過フィルタであるＭＴＩフィルタ４５によ
り、低周波成分であるクラッタ成分は抑圧される。しか
し、完全に抑圧しきれない場合、残存クラッタ成分及び
前記ノイズ成分は、ノイズ判別処理回路４７に入力され
る。
【０００７】このノイズ判別処理回路４７により、入力
信号の強度が所定の閾値より大きい場合は、ドプラ成分
と判断し、入力信号をそのまま出力する。また、入力信
号の強度が所定の閾値より小さい場合は、ノイズ成分で
あると判断し、入力信号を除去する。
信号の強度が所定の閾値より大きい場合は、ドプラ成分
と判断し、入力信号をそのまま出力する。また、入力信
号の強度が所定の閾値より小さい場合は、ノイズ成分で
あると判断し、入力信号を除去する。
【０００８】このため、強度の小さいドプラ成分である
血流情報は失われてしまう。また、強度の大きい残存ク
ラッタ成分、又はノイズ成分は除去されることなく出力
される場合が生じる。このように従来の超音波血流観測
装置においては、血流情報の欠落、クラッタ成分又はノ
イズ成分の残存があるため、局所的に画質を劣化させる
という問題があった。
血流情報は失われてしまう。また、強度の大きい残存ク
ラッタ成分、又はノイズ成分は除去されることなく出力
される場合が生じる。このように従来の超音波血流観測
装置においては、血流情報の欠落、クラッタ成分又はノ
イズ成分の残存があるため、局所的に画質を劣化させる
という問題があった。
【０００９】本発明は、上記のような問題点を解決する
ためになされたもので、信号欠落の補填、残存ノイズ成
分の除去により前記血流画像の画質を改善し、画質の良
い超音波血流観測装置を提供することを目的とする。
ためになされたもので、信号欠落の補填、残存ノイズ成
分の除去により前記血流画像の画質を改善し、画質の良
い超音波血流観測装置を提供することを目的とする。
【００１０】
【課題を解決するための手段】上記目的を達成するため
に本発明は、超音波ビームをスキャンしながら、超音波
パルスを生体内に送信する送信装置と、前記送信された
超音波パルスによる生体からのエコーを受信する受信装
置と、前記受信装置からの受信信号により血流情報を得
るドプラ処理手段と、前記血流情報を超音波画像として
表示する手段と、を含む超音波血流観測装置において、
前記超音波画像を構成する各画素を注目画素とし、その
注目画素を中心とする複数の画素で構成される窓を設定
し、その窓内の各画素の画素値を抽出する窓抽出手段
と、前記窓抽出手段により抽出された前記窓内の各画素
の画素値を多値化する多値化手段と、前記多値化された
注目画素の画素値に応じて及び前記多値化された前記窓
内における各画素値の個数に応じて、前記注目画素が実
信号であるか、実信号の欠落であるか、又はノイズであ
るかを判別する判別手段と、前記実信号であると判別さ
れた場合には前記注目画素の更新画素値としてそのまま
の画素値を出力し、前記実信号の欠落であると判別され
た場合には前記窓内の他の画素の画素値を利用して前記
注目画素の更新画素値を決定し、前記ノイズであると判
別された場合には前記注目画素の更新画素値として所定
値を出力する画素値選択手段と、を含むことを特徴とす
る。
に本発明は、超音波ビームをスキャンしながら、超音波
パルスを生体内に送信する送信装置と、前記送信された
超音波パルスによる生体からのエコーを受信する受信装
置と、前記受信装置からの受信信号により血流情報を得
るドプラ処理手段と、前記血流情報を超音波画像として
表示する手段と、を含む超音波血流観測装置において、
前記超音波画像を構成する各画素を注目画素とし、その
注目画素を中心とする複数の画素で構成される窓を設定
し、その窓内の各画素の画素値を抽出する窓抽出手段
と、前記窓抽出手段により抽出された前記窓内の各画素
の画素値を多値化する多値化手段と、前記多値化された
注目画素の画素値に応じて及び前記多値化された前記窓
内における各画素値の個数に応じて、前記注目画素が実
信号であるか、実信号の欠落であるか、又はノイズであ
るかを判別する判別手段と、前記実信号であると判別さ
れた場合には前記注目画素の更新画素値としてそのまま
の画素値を出力し、前記実信号の欠落であると判別され
た場合には前記窓内の他の画素の画素値を利用して前記
注目画素の更新画素値を決定し、前記ノイズであると判
別された場合には前記注目画素の更新画素値として所定
値を出力する画素値選択手段と、を含むことを特徴とす
る。
【００１１】又は、上記に、この注目画素更新手段によ
り更新された注目画素の値を線形又は非線形フィルタに
よりフィルタリングするフィルタ手段とを備える。
り更新された注目画素の値を線形又は非線形フィルタに
よりフィルタリングするフィルタ手段とを備える。
【００１２】
【作用】本発明において、以上の手段により、画面の各
画素を順次注目画素とし、その注目画素毎に以下の処理
を行う。注目画素の値を血流情報と判断した時は、注目
画素の値をそのまま注目画素の値とする。また、注目画
素の値を欠落情報と判断した時は、窓内の注目画素以外
の全画素の重み付き平均値を注目画素の値とする。更
に、注目画素の値をノイズであると判断した時は、注目
画素の値を画像情報がない時の画素の値とする。
画素を順次注目画素とし、その注目画素毎に以下の処理
を行う。注目画素の値を血流情報と判断した時は、注目
画素の値をそのまま注目画素の値とする。また、注目画
素の値を欠落情報と判断した時は、窓内の注目画素以外
の全画素の重み付き平均値を注目画素の値とする。更
に、注目画素の値をノイズであると判断した時は、注目
画素の値を画像情報がない時の画素の値とする。
【００１３】又は、上記注目画素の最終結果を更にフィ
ルタリングすることにより各画素の値を平滑化する。
ルタリングすることにより各画素の値を平滑化する。
【００１４】
【実施例】図１に、本発明に係る超音波血流観測装置の
概略構成図を示す。図１の探触子４１及び送受信回路４
２により、超音波が生体内に送信され、生体内から返っ
てくるエコー信号が受信される。受信回路で所定の処理
を受けたエコー信号は、断層データ処理系１３及びドプ
ラ信号処理系１１に入力される。断層データ処理系１３
はＢモード断層像を得るための回路である。その出力は
画像表示回路４８に入力され、表示器４９に入力され
る。一方、ドプラ信号処理系１１は、図４の点線枠内に
示す回路構成であり、上述したように血流の平均速度又
は血流の速度分散又はドプラ信号の強度の血流情報を得
るための回路である。その出力は画像処理回路１２に入
力される。本実施例において特徴的な前記画像処理回路
１２は、画像処理により信号（血流情報）の欠落を補填
し、残存ノイズを除去する。この画像処理回路１２の内
部構成を図２に示す。以下にその動作を説明する。
概略構成図を示す。図１の探触子４１及び送受信回路４
２により、超音波が生体内に送信され、生体内から返っ
てくるエコー信号が受信される。受信回路で所定の処理
を受けたエコー信号は、断層データ処理系１３及びドプ
ラ信号処理系１１に入力される。断層データ処理系１３
はＢモード断層像を得るための回路である。その出力は
画像表示回路４８に入力され、表示器４９に入力され
る。一方、ドプラ信号処理系１１は、図４の点線枠内に
示す回路構成であり、上述したように血流の平均速度又
は血流の速度分散又はドプラ信号の強度の血流情報を得
るための回路である。その出力は画像処理回路１２に入
力される。本実施例において特徴的な前記画像処理回路
１２は、画像処理により信号（血流情報）の欠落を補填
し、残存ノイズを除去する。この画像処理回路１２の内
部構成を図２に示す。以下にその動作を説明する。
【００１５】初めに、ドプラ信号処理系１１による血流
情報を画像メモリ２０に格納し、この画像メモリ２０の
画像情報を原画像とする。次に、アドレス発生部２１か
ら発生するアドレスと各画素を一対一に対応させ、この
各画素をアドレス発生部２１により選択されるアドレス
毎に順次注目画素とし、この注目画素を中心とするｉ×
ｊの画素により構成されるマトリックスを窓とする。こ
の窓１個分の画素による情報が順次取り出され、絶対値
化回路２２と注目画素値抽出部２５と平均値算出部２６
の各々に入力される。ｉとｊの方向は各々超音波ビーム
のスキャン方向と超音波ビームのビーム方向に対応す
る。本実施例ではｉとｊの値を共に３とし、３×３の場
合を示す。
情報を画像メモリ２０に格納し、この画像メモリ２０の
画像情報を原画像とする。次に、アドレス発生部２１か
ら発生するアドレスと各画素を一対一に対応させ、この
各画素をアドレス発生部２１により選択されるアドレス
毎に順次注目画素とし、この注目画素を中心とするｉ×
ｊの画素により構成されるマトリックスを窓とする。こ
の窓１個分の画素による情報が順次取り出され、絶対値
化回路２２と注目画素値抽出部２５と平均値算出部２６
の各々に入力される。ｉとｊの方向は各々超音波ビーム
のスキャン方向と超音波ビームのビーム方向に対応す
る。本実施例ではｉとｊの値を共に３とし、３×３の場
合を示す。
【００１６】図２において、ドプラ信号処理系１１より
得られる血流の平均速度には、速度の方向を示す正又は
負の符号を有するが、絶対値化回路２２では、この平均
速度の絶対値化を行う。ただし、前記画像メモリ２０の
情報が、血流の平均速度画像である時、前記絶対値化回
路２２を必要とするが、画像メモリ２０の情報が、血流
の速度分散又はドプラ信号の強度である時、この強度は
常に正値であり、絶対値化回路２２を省略できる。
得られる血流の平均速度には、速度の方向を示す正又は
負の符号を有するが、絶対値化回路２２では、この平均
速度の絶対値化を行う。ただし、前記画像メモリ２０の
情報が、血流の平均速度画像である時、前記絶対値化回
路２２を必要とするが、画像メモリ２０の情報が、血流
の速度分散又はドプラ信号の強度である時、この強度は
常に正値であり、絶対値化回路２２を省略できる。
【００１７】この絶対値化回路２２の出力は二値化回路
２３に入力される。二値化回路２３では、絶対値化回路
２２より得られる窓内の各画素の信号の値が、閾値より
大きい値又はそれ以外の値であるかの判定により、窓内
の画素のすべてを“１”又は“０”の二値に分類する二
値化処理を行う。本実施例では、閾値を０とする。つま
り窓内の各画素の値は０であれば、“０”とし、それ以
外は“１”とする。この二値化処理より、窓の各画素が
二値で表現され、３×３画素の二値パターンが得られ
る。この窓の二値パターンには２9 ＝５１２通りの組合
せが存在する。
２３に入力される。二値化回路２３では、絶対値化回路
２２より得られる窓内の各画素の信号の値が、閾値より
大きい値又はそれ以外の値であるかの判定により、窓内
の画素のすべてを“１”又は“０”の二値に分類する二
値化処理を行う。本実施例では、閾値を０とする。つま
り窓内の各画素の値は０であれば、“０”とし、それ以
外は“１”とする。この二値化処理より、窓の各画素が
二値で表現され、３×３画素の二値パターンが得られ
る。この窓の二値パターンには２9 ＝５１２通りの組合
せが存在する。
【００１８】上記３×３画素の二値パターンのデータ、
つまり窓は、９ビットで表現でき、この９ビットによる
窓のデータより、注目画素の属性を決定する。この注目
画素の属性の決定方法は図２の記憶装置を備えるパター
ンテーブル２４との比較により行う。
つまり窓は、９ビットで表現でき、この９ビットによる
窓のデータより、注目画素の属性を決定する。この注目
画素の属性の決定方法は図２の記憶装置を備えるパター
ンテーブル２４との比較により行う。
【００１９】図３に、窓内のパターンと注目画素の属性
との関係を示す。注目画素値が“１”及び“０”の時
の、窓内全体の“１”の数と注目画素の属性との関係
を、図３にまとめる。また、図３の中で特に、注目画素
値が“１”、かつ窓内全体の“１”が２個の時、注目画
素の属性を“血流情報”と判定する場合の全パターンを
図３の（注）にまとめる。
との関係を示す。注目画素値が“１”及び“０”の時
の、窓内全体の“１”の数と注目画素の属性との関係
を、図３にまとめる。また、図３の中で特に、注目画素
値が“１”、かつ窓内全体の“１”が２個の時、注目画
素の属性を“血流情報”と判定する場合の全パターンを
図３の（注）にまとめる。
【００２０】前記パターンテーブル２４による属性の判
定の仕方を前記図３を用い、以下に説明する。図３に示
すように、窓の注目画素が“１”、かつ窓内の“１”の
数が２以上ならば、注目画素の属性を“血流情報”と判
定する。また、注目画素のみが“１”の時、注目画素の
属性は“ノイズ”であると判定する。更に、注目画素が
“１”、かつ窓内の注目画素以外に“１”が１つある場
合、図３の（注）に示す窓のパターン時のみ“血流情
報”と判定し、それ以外は“ノイズ”であると判定す
る。
定の仕方を前記図３を用い、以下に説明する。図３に示
すように、窓の注目画素が“１”、かつ窓内の“１”の
数が２以上ならば、注目画素の属性を“血流情報”と判
定する。また、注目画素のみが“１”の時、注目画素の
属性は“ノイズ”であると判定する。更に、注目画素が
“１”、かつ窓内の注目画素以外に“１”が１つある場
合、図３の（注）に示す窓のパターン時のみ“血流情
報”と判定し、それ以外は“ノイズ”であると判定す
る。
【００２１】更に図３に示すように、注目画素が
“０”、かつ窓内の“１”の数が４以上ならば、注目画
素の属性を“信号の欠落”であると判定する。また、注
目画素が“０”、かつ、窓内の“１”の数が４以下なら
ば、注目画素の属性は“ノイズ”であると判定する。
“０”、かつ窓内の“１”の数が４以上ならば、注目画
素の属性を“信号の欠落”であると判定する。また、注
目画素が“０”、かつ、窓内の“１”の数が４以下なら
ば、注目画素の属性は“ノイズ”であると判定する。
【００２２】上記のようにパターンテーブル２４と窓の
パターンとの比較によるパターン認識により、注目画素
の属性を判定する。この注目画素の属性を示す情報がセ
レクタ２８に入力される。セレクタ２８の入力側には、
前記属性情報以外に注目画素値抽出部２５と平均値算出
部２６と０値レジスタ２７の出力部が入力される。この
注目画素値抽出部２５では、３×３の画素で構成される
窓の中心に位置する注目画素の値のみが抽出される。ま
た、平均値算出部２６では、窓内の注目画素以外の各画
素による重み付き平均値が算出される。また、０値レジ
スタ２７は、０データを保持している。
パターンとの比較によるパターン認識により、注目画素
の属性を判定する。この注目画素の属性を示す情報がセ
レクタ２８に入力される。セレクタ２８の入力側には、
前記属性情報以外に注目画素値抽出部２５と平均値算出
部２６と０値レジスタ２７の出力部が入力される。この
注目画素値抽出部２５では、３×３の画素で構成される
窓の中心に位置する注目画素の値のみが抽出される。ま
た、平均値算出部２６では、窓内の注目画素以外の各画
素による重み付き平均値が算出される。また、０値レジ
スタ２７は、０データを保持している。
【００２３】前記パターンテーブル２４の属性の判定結
果が“血流情報”である時、注目画素値抽出部２５が選
択され、注目画素値抽出部２５の出力結果がセレクタ２
８から出力される。また、パターンテーブル２４の判定
結果が“信号の欠落”である時、平均値算出部２６が選
択され、平均値算出部２６の出力結果がセレクタ２８か
ら出力される。パターンテーブル２４の判定結果が“ノ
イズ”である時、０値レジスタ２７が選択され、０値レ
ジスタ２７の出力結果がセレクタ２８から出力される。
果が“血流情報”である時、注目画素値抽出部２５が選
択され、注目画素値抽出部２５の出力結果がセレクタ２
８から出力される。また、パターンテーブル２４の判定
結果が“信号の欠落”である時、平均値算出部２６が選
択され、平均値算出部２６の出力結果がセレクタ２８か
ら出力される。パターンテーブル２４の判定結果が“ノ
イズ”である時、０値レジスタ２７が選択され、０値レ
ジスタ２７の出力結果がセレクタ２８から出力される。
【００２４】以上より、原画像において、注目画素が、
“血流情報”と識別される場合は、注目画素の値をその
ままに画素の信号とする。また、“信号の欠落”と識別
される場合は、窓内の注目画素以外の８個の画素の重み
付き平均値を画素の信号とし、信号の欠落を補填する。
“ノイズ”と識別される場合は、０が出力され、ノイズ
が除去される。このようにして、血流イメージの画質の
改善を実現することが可能である。
“血流情報”と識別される場合は、注目画素の値をその
ままに画素の信号とする。また、“信号の欠落”と識別
される場合は、窓内の注目画素以外の８個の画素の重み
付き平均値を画素の信号とし、信号の欠落を補填する。
“ノイズ”と識別される場合は、０が出力され、ノイズ
が除去される。このようにして、血流イメージの画質の
改善を実現することが可能である。
【００２５】更に、このセレクタの出力を、メディアン
フィルタ２９に入力し、ノイズの除去と信号欠落の効果
を更に上げて、血流イメージの画質を改善することがで
きる。このため、メディアンフィルタ２９の出力結果を
画像表示回路４８の入力信号とすることにより、更に画
質の改善された血流情報の画像が得られる。
フィルタ２９に入力し、ノイズの除去と信号欠落の効果
を更に上げて、血流イメージの画質を改善することがで
きる。このため、メディアンフィルタ２９の出力結果を
画像表示回路４８の入力信号とすることにより、更に画
質の改善された血流情報の画像が得られる。
【００２６】以上、本発明の一実施例について説明した
が、原画像上におけるノイズと信号の欠落の多少に応じ
て、ｉ×ｊの窓の画素を、例えば、３×５と変えたり、
図３の（注）に示す二値パターン認識表の内容を変えた
り、画像の平滑化処理をするためのフィルタをメディア
ンフィルタ以外の一次元又は二次元フィルタに変えた
り、フィルタを使用しないなど、種々の仕様により実現
可能である。また、本実施例では、注目画素の属性の識
別について、窓内の“１”の値を持つ画素の数を用いた
が、“０”の値を持つ画素の数によっても同じ結果を得
ることが可能であることは言うまでもない。更にまた、
符号付きの血流平均速度の原画像について、窓内の画素
値を符号に応じて、二値ではなく、“−１”、“０”、
“１”の三値で三値化して、注目画素の値、窓内に
“１”又は“−１”の値を持つ画素の数、方向から注目
画素の属性を識別する変化も利用できる。
が、原画像上におけるノイズと信号の欠落の多少に応じ
て、ｉ×ｊの窓の画素を、例えば、３×５と変えたり、
図３の（注）に示す二値パターン認識表の内容を変えた
り、画像の平滑化処理をするためのフィルタをメディア
ンフィルタ以外の一次元又は二次元フィルタに変えた
り、フィルタを使用しないなど、種々の仕様により実現
可能である。また、本実施例では、注目画素の属性の識
別について、窓内の“１”の値を持つ画素の数を用いた
が、“０”の値を持つ画素の数によっても同じ結果を得
ることが可能であることは言うまでもない。更にまた、
符号付きの血流平均速度の原画像について、窓内の画素
値を符号に応じて、二値ではなく、“−１”、“０”、
“１”の三値で三値化して、注目画素の値、窓内に
“１”又は“−１”の値を持つ画素の数、方向から注目
画素の属性を識別する変化も利用できる。
【００２７】
【発明の効果】以上説明したように、本発明によれば、
各画素の値を忠実に伝達又はノイズの除去又は欠落情報
の補填を行うことにより、画質の良い血流情報を得る超
音波血流観測装置を提供できる。
各画素の値を忠実に伝達又はノイズの除去又は欠落情報
の補填を行うことにより、画質の良い血流情報を得る超
音波血流観測装置を提供できる。
【００２８】また、上記注目画素の最終結果を更にフィ
ルタリングすることにより各画素の値を平滑化し、更に
画質の良い血流情報を得る超音波血流観測装置を提供で
きる。
ルタリングすることにより各画素の値を平滑化し、更に
画質の良い血流情報を得る超音波血流観測装置を提供で
きる。
【図１】本発明の一実施例に係る超音波血流観測装置の
概略構成図である。
概略構成図である。
【図２】本発明の一実施例に係る画像処理回路の構成図
である。
である。
【図３】本発明の一実施例に係る窓内の“１”の数及び
注目画素の値と注目画素の属性との関係図である。
注目画素の値と注目画素の属性との関係図である。
【図４】従来の超音波血流情報装置のドプラ信号処理系
の概略構成図である。
の概略構成図である。
２０ 画像メモリ ２１ アドレス発生部 ２２ 絶対値化回路 ２３ 二値化回路 ２４ パターンテーブル ２５ 注目画素値抽出部 ２６ 平均値算出部 ２７ ０値レジスタ ２８ セレクタ ２９ メディアンフィルタ
フロントページの続き (56)参考文献 特開 平４−161146（ＪＰ，Ａ) 特開 昭56−91738（ＪＰ，Ａ) 特開 昭61−201372（ＪＰ，Ａ) 特開 平１−303132（ＪＰ，Ａ) 特開 昭62−270139（ＪＰ，Ａ)
Claims (4)
- 【請求項１】 超音波ビームをスキャンしながら、超音
波パルスを生体内に送信する送信装置と、前記送信された 超音波パルスによる生体からのエコーを
受信する受信装置と、前記受信装置からの受信信号により血流情報を得るドプ
ラ処理手段と、 前記血流情報を超音波画像として表示する手段と、 を含む 超音波血流観測装置において、前記超音波画像 を構成する各画素を注目画素とし、その
注目画素を中心とする複数の画素で構成される窓を設定
し、その窓内の各画素の画素値を抽出する窓抽出手段
と、前記 窓抽出手段により抽出された前記窓内の各画素の画
素値を多値化する多値化手段と、 前記多値化された注目画素の画素値に応じて及び前記多
値化された前記窓内における各画素値の個数に応じて、
前記注目画素が実信号であるか、実信号の欠落である
か、又はノイズであるかを判別する判別手段と、 前記実信号であると判別された場合には前記注目画素の
更新画素値としてそのままの画素値を出力し、前記実信
号の欠落であると判別された場合には前記窓内の他の画
素の画素値を利用して前記注目画素の更新画素値を決定
し、前記ノイズであると判別された場合には前記注目画
素の更新画素値として所定値を出力する画素値選択手段
と、 を含むことを特徴とする超音波血流観測装置。 - 【請求項２】 請求項１記載の装置において、 前記画素値選択手段から出力された前記注目画素の更新
画素値を線形又は非線形フィルタによりフィルタリング
する手段 を備えることを特徴とする超音波血流観測装
置。 - 【請求項３】 超音波ビームをスキャンしながら、超音
波パルスを生体内に送信する送信装置と、 前記送信された超音波パルスによる生体からのエコーを
受信する受信装置と、 前記受信装置からの受信信号により血流情報を得るドプ
ラ処理手段と、 前記血流情報を超音波画像として表示する手段と、 を含む超音波血流観測装置において、 前記超音波画像を構成する各画素を注目画素とし、その
注目画素を中心とする複数の画素で構成される窓を設定
し、その窓内の各画素の画素値を抽出する窓抽出手段
と、 前記窓抽出手段により抽出された前記窓内の各画素の画
素値を二値化する二値化手段と、 前記二値化された注目画素の画素値が１か又は０かに応
じて及び前記二値化された前記窓内の画素値１又は画素
値０の個数に応じて、前記注目画素が実信号であるか、
実信号の欠落であるか、又はノイズであるかを判別する
判別手段と、 前記実信号であると判別された場合には前記注目画素の
更新画素値としてそのままの画素値を出力し、前記実信
号の欠落であると判別された場合には前記窓内の他の画
素の画素値を利用して前記注目画素の更新画素値を決定
し、前記ノイズであると判別された場合には前記注目画
素の更新画素値として所定値を出力する画素値選択手段
と、 を含むことを特徴とする超音波血流観測装置。 - 【請求項４】 請求項３記載の装置において、 前記窓は３画素×３画素の大きさを有し、 前記判別手段は、 前記二値化された注目画素の画素値が１の場合であって
前記窓内の画素値１の個数が３以上の場合に、又は、前
記二値化された注目画素の画素値が１の場合であって前
記窓内の画素値１が２個ありかつ縦又は横に画素値１が
連続する場合に、前記注目画素の画素値が実信号である
と判別し、 前記二値化された注目画素の画素値が０の場合であって
前記窓内の画素値１の個数が５以上である場合に、前記
注目画素の画素値が実信号の欠落であると判別し、 上記条件以外の場合に前記注目画素の画素値がノイズで
あると判別することを特徴とする超音波血流観測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5255953A JP2723453B2 (ja) | 1993-10-13 | 1993-10-13 | 超音波血流観測装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5255953A JP2723453B2 (ja) | 1993-10-13 | 1993-10-13 | 超音波血流観測装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07108043A JPH07108043A (ja) | 1995-04-25 |
| JP2723453B2 true JP2723453B2 (ja) | 1998-03-09 |
Family
ID=17285877
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5255953A Expired - Fee Related JP2723453B2 (ja) | 1993-10-13 | 1993-10-13 | 超音波血流観測装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2723453B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4528529B2 (ja) * | 2004-01-20 | 2010-08-18 | 株式会社東芝 | 超音波診断装置及び超音波画像データ処理方法 |
| JP4616614B2 (ja) * | 2004-10-21 | 2011-01-19 | アロカ株式会社 | 超音波診断装置 |
| KR101055581B1 (ko) * | 2008-03-20 | 2011-08-23 | 삼성메디슨 주식회사 | 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법 |
| CN112263279B (zh) * | 2020-11-18 | 2023-08-25 | 飞依诺科技股份有限公司 | 高帧率平面波远场噪声抑制方法、设备及存储介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5691738A (en) * | 1979-12-27 | 1981-07-24 | Tokyo Shibaura Electric Co | Ultrasonic diagnosis apparatus |
| JPS61201372A (ja) * | 1985-03-02 | 1986-09-06 | Toshiba Corp | 画像処理装置 |
| JP2639391B2 (ja) * | 1986-03-06 | 1997-08-13 | ヒューレット・パッカード・カンパニー | 流体マップの画像形成装置 |
| JPH01303132A (ja) * | 1988-05-31 | 1989-12-07 | Yokogawa Medical Syst Ltd | 画像処理装置 |
| JP2950600B2 (ja) * | 1990-10-26 | 1999-09-20 | 株式会社東芝 | カラードプラ超音波診断装置 |
-
1993
- 1993-10-13 JP JP5255953A patent/JP2723453B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07108043A (ja) | 1995-04-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2016905B1 (en) | Ultrasound diagnostic apparatus | |
| US6277075B1 (en) | Method and apparatus for visualization of motion in ultrasound flow imaging using continuous data acquisition | |
| EP1358507B1 (en) | Adaptive image processing for spatial compounding | |
| US8684934B2 (en) | Adaptively performing clutter filtering in an ultrasound system | |
| JP3406106B2 (ja) | 超音波画像表示方法および超音波診断装置 | |
| KR100742467B1 (ko) | 촬상 시스템, 촬상 방법 및 혈액 움직임 촬상 시스템 | |
| US20010031921A1 (en) | Ultrasonic image processing method and examination system for displaying an ultrasonic color-coded image sequence of an object having moving parts | |
| US5879302A (en) | Imaging method and apparatus using MIP projection of doppler power values | |
| US9039620B2 (en) | Ultrasound diagnostic apparatus | |
| JP3397748B2 (ja) | カラードップラー映像システムにおけるカラー映像表示方法及び装置 | |
| EP1796037B1 (en) | Image processing system and method for controlling gains for color flow images | |
| EP0406909A2 (en) | Ultrasonic diagnosing apparatus | |
| JP2000149015A (ja) | 画像のエッジ強調方法及びイメ―ジング装置 | |
| JP4347860B2 (ja) | 超音波診断装置 | |
| JP2723453B2 (ja) | 超音波血流観測装置 | |
| JP2003529422A (ja) | Ｂカラー優先閾値計算を用いた超音波システム | |
| US6048313A (en) | Method and apparatus for fractal-based enhancement of ultrasound imaging | |
| US6045504A (en) | Method and apparatus for polynomial approximation of nonlinear operations in medical ultrasound imaging | |
| JPH0751134B2 (ja) | 超音波診断画像処理装置 | |
| JP3267739B2 (ja) | 超音波カラードプラ診断装置 | |
| JP3514544B2 (ja) | 超音波診断装置 | |
| JP2723458B2 (ja) | 超音波ドプラ診断装置 | |
| JP4616614B2 (ja) | 超音波診断装置 | |
| JP3043862B2 (ja) | 超音波診断装置 | |
| JP3569330B2 (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101128 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121128 Year of fee payment: 15 |
|
| LAPS | Cancellation because of no payment of annual fees |