JPH0225785A - 超音波エコーグラフ装置 - Google Patents
超音波エコーグラフ装置Info
- Publication number
- JPH0225785A JPH0225785A JP1132933A JP13293389A JPH0225785A JP H0225785 A JPH0225785 A JP H0225785A JP 1132933 A JP1132933 A JP 1132933A JP 13293389 A JP13293389 A JP 13293389A JP H0225785 A JPH0225785 A JP H0225785A
- Authority
- JP
- Japan
- Prior art keywords
- echographic
- ultrasonic
- component
- receiving
- region
- 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.)
- Pending
Links
- 230000000877 morphologic effect Effects 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000003628 erosive effect Effects 0.000 claims description 23
- 238000002604 ultrasonography Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 4
- 230000000737 periodic effect Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 13
- 230000033001 locomotion Effects 0.000 description 5
- 206010011732 Cyst Diseases 0.000 description 3
- 208000031513 cyst Diseases 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000007373 indentation Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 206010039729 Scotoma Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000000007 visual effect 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/5206—Two-dimensional coordinated display of distance and direction; B-scan display
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明は、少なくとも1個の圧電変換器と、超音波ビー
ム伝達ステージと、及び前記変換器へ帰還したエコーグ
ラフ信号eを受信し且つ処理するステージとを具える超
音波エコーグラフ装置に関するものである。
ム伝達ステージと、及び前記変換器へ帰還したエコーグ
ラフ信号eを受信し且つ処理するステージとを具える超
音波エコーグラフ装置に関するものである。
(従来の技術とその問題点)
本発明は、特に器官の映像の形成用の医療エコーグラフ
の分野で実際に魅力的に用いれる。
の分野で実際に魅力的に用いれる。
前述の種類の何らかの超音波エコーグラフ装置によって
解決されるべき一般技術問題は、できる限り正確な映像
が、その輪郭に関してと同時にそこに含まれる鏡状壁に
関しても、試験される対象のものとして得られねばなら
ぬのが一般的である。
解決されるべき一般技術問題は、できる限り正確な映像
が、その輪郭に関してと同時にそこに含まれる鏡状壁に
関しても、試験される対象のものとして得られねばなら
ぬのが一般的である。
種々の解答がこの一般技術問題について既に提案されて
きた。特に雑誌フィジーク・アップリケ(Physiq
ue App1iqu6e)18(1983)の527
〜556頁にに発表されたエム・フィフク(M、Fin
k)により論文“イマジュリー・ウルトラソノール”
(Imagerieultrasonore)が参照
され、その論文はエコーグラフ映像の形成について種々
の局面を論じ尽すほどに分析している。それらがしばし
ば満足な結果を提供するとしても、すべての既知の解答
はそれにも拘らず所定の限界を有しない。器官の輪郭の
決定と低コントラストの対象の探知とは、研究されてい
る対象中の複数の散乱点により生じる建設的及び破壊的
な反響妨害、及び圧電変換器により伝達される超音波の
コヒーレント性による反響妨害によって妨害される。音
響分野におけるこの寄生現象は、光学分野において“ス
ペックル”として知られ且つレーザー放射でしばしば遭
遇する“ライト・グレインズを引き起こす現象と類似す
る。
きた。特に雑誌フィジーク・アップリケ(Physiq
ue App1iqu6e)18(1983)の527
〜556頁にに発表されたエム・フィフク(M、Fin
k)により論文“イマジュリー・ウルトラソノール”
(Imagerieultrasonore)が参照
され、その論文はエコーグラフ映像の形成について種々
の局面を論じ尽すほどに分析している。それらがしばし
ば満足な結果を提供するとしても、すべての既知の解答
はそれにも拘らず所定の限界を有しない。器官の輪郭の
決定と低コントラストの対象の探知とは、研究されてい
る対象中の複数の散乱点により生じる建設的及び破壊的
な反響妨害、及び圧電変換器により伝達される超音波の
コヒーレント性による反響妨害によって妨害される。音
響分野におけるこの寄生現象は、光学分野において“ス
ペックル”として知られ且つレーザー放射でしばしば遭
遇する“ライト・グレインズを引き起こす現象と類似す
る。
超音波分野におけるこの妨害雑音の第一の欠点は、一方
ではエコーグラフ映像の質の劣化、特に低コントラスト
の場合における詳細部分の消失であり、他方では広い周
波数帯域と普通の輪郭探知法を非能率にする強い分散と
である。
ではエコーグラフ映像の質の劣化、特に低コントラスト
の場合における詳細部分の消失であり、他方では広い周
波数帯域と普通の輪郭探知法を非能率にする強い分散と
である。
従って“スペックル″の低減がエコーグラフ映像の改善
には決定的である。これを達成するために種々の方法が
提案されてきた。
には決定的である。これを達成するために種々の方法が
提案されてきた。
一方では、スペックルの関係除去を得るために捕捉に先
だった信号の処理から成る方法がある。
だった信号の処理から成る方法がある。
これは、例えば空間構成である〔米国電気・電子通信学
会(IBBB)雑誌、音響光学、超音波学(Sonic
s Ultrasonics) 、5025巻1−6号
(1978)中のシー・ビー・パークハード(C,B、
Burkhardt)による論文参照〕。この方法は
一般に高品質のエコーグラフ映像を提供する。それは、
然し乍ら、複雑な電子的処理装置の使用を必要とすると
いう欠点を有する。更に、捕捉時間が比較的長く、だか
ら映像のリアルタイム形式が妨げられる。
会(IBBB)雑誌、音響光学、超音波学(Sonic
s Ultrasonics) 、5025巻1−6号
(1978)中のシー・ビー・パークハード(C,B、
Burkhardt)による論文参照〕。この方法は
一般に高品質のエコーグラフ映像を提供する。それは、
然し乍ら、複雑な電子的処理装置の使用を必要とすると
いう欠点を有する。更に、捕捉時間が比較的長く、だか
ら映像のリアルタイム形式が妨げられる。
他方では、例えばスペックルによる変動を有効に低減し
、然し輪郭を拡張することにより映像内へ汚れを導入す
る低域濾波器によって、映像の形成の後に実行される平
滑化方法がある。
、然し輪郭を拡張することにより映像内へ汚れを導入す
る低域濾波器によって、映像の形成の後に実行される平
滑化方法がある。
従って、本発明により解決されるべき技術問題は、少な
くとも1個の圧電変換器と、超音波ビーム伝達ステージ
と、及び前記変換器へ帰還したエコーグラフ信号eを受
信し且つ処理するステージとを具える超音波エコーグラ
フ装置を提案することであり、その超音波エコーグラフ
装置は“信号後処理”形の装置であるが、明確な映像輪
郭を維持し乍らスペックル効果を実質的に低減しなけれ
ばならない。
くとも1個の圧電変換器と、超音波ビーム伝達ステージ
と、及び前記変換器へ帰還したエコーグラフ信号eを受
信し且つ処理するステージとを具える超音波エコーグラ
フ装置を提案することであり、その超音波エコーグラフ
装置は“信号後処理”形の装置であるが、明確な映像輪
郭を維持し乍らスペックル効果を実質的に低減しなけれ
ばならない。
(課題を解決するための手段)
本発明によって、この問題は、前記受信し且つ処理する
ステージが、領域り内に構成要素gを有し且つ前記エコ
ーグラフ信号へ適用される基本的な侵食及び拡張手段に
よって実現される、少なくとも1個のグレイスケール形
態学的フィルターを具えることにより解決される。
ステージが、領域り内に構成要素gを有し且つ前記エコ
ーグラフ信号へ適用される基本的な侵食及び拡張手段に
よって実現される、少なくとも1個のグレイスケール形
態学的フィルターを具えることにより解決される。
数学的形態学とその映像処理への応用は、ジエイ・セラ
(J、 5erra)による1982年の学術図書“映
像分析と数学的形態学”から、及びパターン分析及び機
械知能についての米国電気・電子学会(IEBB)雑誌
PAMI−9巻第4号(1987年7月)の532〜5
50頁から知られる。これらの文書は基本的に2進映像
の形態学的分析に関し、且つ侵食及び拡張動作を介して
、不規則なものが除去され且つそれで正確な測定が実行
され得る映像がいかにして構成され得るかを開示してい
る。数学的形態学は構成の量を計ることを基本的に目指
し、且つ映像はその映像とその映像を横切って移動され
る構成要素として参照“される1個の小さい限られた要
素との間の関係によって具体化する。これらの関係は基
本的には、欠陥(def 1ciency)により更に
構成された映像を図式的に算出する(正確な性質は以下
に詳細に説明する)侵食(erosion)動作(e
egとして表示される)と、逆に過剰(excess)
により更に構成された映像を算出する拡張(dilat
ion)動作(e■gとして表示される)とから生じる
。厳密に言えば、侵食と拡張とは等次元でないから形態
学的フィルターではない。従って、形態学的動作の侵食
及び拡張に基づくものは閉鎖(closing)及び開
放(opening)動作として参照し編集される。
(J、 5erra)による1982年の学術図書“映
像分析と数学的形態学”から、及びパターン分析及び機
械知能についての米国電気・電子学会(IEBB)雑誌
PAMI−9巻第4号(1987年7月)の532〜5
50頁から知られる。これらの文書は基本的に2進映像
の形態学的分析に関し、且つ侵食及び拡張動作を介して
、不規則なものが除去され且つそれで正確な測定が実行
され得る映像がいかにして構成され得るかを開示してい
る。数学的形態学は構成の量を計ることを基本的に目指
し、且つ映像はその映像とその映像を横切って移動され
る構成要素として参照“される1個の小さい限られた要
素との間の関係によって具体化する。これらの関係は基
本的には、欠陥(def 1ciency)により更に
構成された映像を図式的に算出する(正確な性質は以下
に詳細に説明する)侵食(erosion)動作(e
egとして表示される)と、逆に過剰(excess)
により更に構成された映像を算出する拡張(dilat
ion)動作(e■gとして表示される)とから生じる
。厳密に言えば、侵食と拡張とは等次元でないから形態
学的フィルターではない。従って、形態学的動作の侵食
及び拡張に基づくものは閉鎖(closing)及び開
放(opening)動作として参照し編集される。
閉鎖は侵食に追従される拡張であり、一方開放は拡張に
追従される侵食である。閉鎖と開放とは映像を再構成し
、その映像は更に構成されており且つその寸法は最初の
映像の寸法とほぼ等しいが、然し乍ら、閉鎖はやや大き
くなり、開放はやや小さくなる。映像が2進ではなくて
アナログである普通のエコーグラフでは、この説明で後
に定義するグレースケール形態学的フィルターを使用す
る。
追従される侵食である。閉鎖と開放とは映像を再構成し
、その映像は更に構成されており且つその寸法は最初の
映像の寸法とほぼ等しいが、然し乍ら、閉鎖はやや大き
くなり、開放はやや小さくなる。映像が2進ではなくて
アナログである普通のエコーグラフでは、この説明で後
に定義するグレースケール形態学的フィルターを使用す
る。
すべてのこれらの動作において、構成する要素の領域り
の大きさが非常に重要であることは明らかになるだろう
。実際に、大きい構成要素は映像の不規則なものが一層
有効に吸収でき、従って実質的に“スペックル”効果を
低減する。然し乍ら、それは同時に考慮された装置の分
解能の劣化をきたす。従って、信号変動の平滑化と小さ
い細部の保全との間の適度な妥協を見出すことが魅力的
である。本発明による特別の変形では、構成要素gの領
域りの寸法がエコーグラフ装置の分解能をなるべく超過
しない方がよい。従って、医療分野では決定的である装
置の分解能が維持される。−層正確に言えば、本出願人
は魅力的な形態は構成要素の寸法がエコーグラフ装置の
分解能の半分にほぼ等しい形態であることを実証した。
の大きさが非常に重要であることは明らかになるだろう
。実際に、大きい構成要素は映像の不規則なものが一層
有効に吸収でき、従って実質的に“スペックル”効果を
低減する。然し乍ら、それは同時に考慮された装置の分
解能の劣化をきたす。従って、信号変動の平滑化と小さ
い細部の保全との間の適度な妥協を見出すことが魅力的
である。本発明による特別の変形では、構成要素gの領
域りの寸法がエコーグラフ装置の分解能をなるべく超過
しない方がよい。従って、医療分野では決定的である装
置の分解能が維持される。−層正確に言えば、本出願人
は魅力的な形態は構成要素の寸法がエコーグラフ装置の
分解能の半分にほぼ等しい形態であることを実証した。
(実施例)
以下、例をあげて、添付の図面を参照しながら本発明の
詳細な説明する。
詳細な説明する。
第1図は超音波エコーグラフ装置のブロック線図を示し
、圧電変換器100と、超音波ビーム伝達ステージ20
0と、及び圧電変換器100へ帰還したエコーグラフ信
号eを受信する受信及び処理ステージ300と、同時に
機械的な変換器走査制御装置101 とを具える。この
変換器の代わりに、電子的走査制御装置を伴う変換器の
列を使用することもできる。この伝達ステージ200は
、例えば32MHzの周波数を有する発振器から成る普
通のシーケンサ−と、及び5 MHzの繰り返えし周波
数で発生器を制御する分周器とを具え、その発生器の電
気的励起信号がこれらの信号を超音波パルス信号の周期
的な列に変換する圧電変換器100へ印加される。
、圧電変換器100と、超音波ビーム伝達ステージ20
0と、及び圧電変換器100へ帰還したエコーグラフ信
号eを受信する受信及び処理ステージ300と、同時に
機械的な変換器走査制御装置101 とを具える。この
変換器の代わりに、電子的走査制御装置を伴う変換器の
列を使用することもできる。この伝達ステージ200は
、例えば32MHzの周波数を有する発振器から成る普
通のシーケンサ−と、及び5 MHzの繰り返えし周波
数で発生器を制御する分周器とを具え、その発生器の電
気的励起信号がこれらの信号を超音波パルス信号の周期
的な列に変換する圧電変換器100へ印加される。
この伝達ステージ200と受信及び処理ステージ300
とを分離するためのT/Rスイッチ102が圧電変換器
100と前記ステージ200.300との間に挿入され
、受信中の回路が伝達信号により妨害されることを防止
する。
とを分離するためのT/Rスイッチ102が圧電変換器
100と前記ステージ200.300との間に挿入され
、受信中の回路が伝達信号により妨害されることを防止
する。
受信及び処理ステージ300 は、T/Rスイッチ10
2の出力端子へ接続された、深さの関数としての利得補
償を有する高周波増幅器301と、対数圧縮増幅器30
4と、記憶及び走査変換装置305と、及び表示装置3
06とを具える。第1図かられかるように、受信及び処
理ステージは領域り内に構成要素gを有するグレイスケ
ール形態学的フィルター302゜303をも具え、エコ
ーグラフ信号eへ印加される侵食手段302e、 30
3e及び拡張手段302d、 303dによって実現さ
れる。
2の出力端子へ接続された、深さの関数としての利得補
償を有する高周波増幅器301と、対数圧縮増幅器30
4と、記憶及び走査変換装置305と、及び表示装置3
06とを具える。第1図かられかるように、受信及び処
理ステージは領域り内に構成要素gを有するグレイスケ
ール形態学的フィルター302゜303をも具え、エコ
ーグラフ信号eへ印加される侵食手段302e、 30
3e及び拡張手段302d、 303dによって実現さ
れる。
拡張手段302d、 303dは2元信号e (x、
y)を下記関係式による付加的定義で与えられる拡張さ
れた信号D(x、y)に変形する。
y)を下記関係式による付加的定義で与えられる拡張さ
れた信号D(x、y)に変形する。
x’ 、 y’がDに属するとき
D(xI y) = (e+g) (x、 y)=Ma
x [:e(x+x’ 、 y+y’ )+g(x’
、 y’ ) )逆に、信号e(x、y)に対して侵食
手段302e、 303eは下式で定義される侵食され
た信号ε(x、 y)を生産する。
x [:e(x+x’ 、 y+y’ )+g(x’
、 y’ ) )逆に、信号e(x、y)に対して侵食
手段302e、 303eは下式で定義される侵食され
た信号ε(x、 y)を生産する。
x/ 、 y/がDに属するとき
E(x、 y) = (eeg) (x、 y)=Mi
n (:e(x+x’ 、 y+y’ )−g(x’
、 y’ ) )第2図の最初のエコーグラフ信号e(
x)に対するこれら2つの動作を図解し、本例では、こ
の最初のaに示したエコーグラフ信号e(x)は散乱環
境内の嚢胞を表現するエコーの包路線である。スペック
ルによる信号の変動は散乱背景中と同様に嚢胞の領域中
にも明らかに見える。第2図b−fは“平らな”要素を
参照され且つ次式によって定義される構成要素g(x)
によりこの信号に実行された種々の形態学的作用を図解
する。
n (:e(x+x’ 、 y+y’ )−g(x’
、 y’ ) )第2図の最初のエコーグラフ信号e(
x)に対するこれら2つの動作を図解し、本例では、こ
の最初のaに示したエコーグラフ信号e(x)は散乱環
境内の嚢胞を表現するエコーの包路線である。スペック
ルによる信号の変動は散乱背景中と同様に嚢胞の領域中
にも明らかに見える。第2図b−fは“平らな”要素を
参照され且つ次式によって定義される構成要素g(x)
によりこの信号に実行された種々の形態学的作用を図解
する。
領域りを通じてg (x) =0
この構成要素と領域りとは第2図aに示される。
与えられたものから不正な描出へ導き得る信号への大き
い歪を防止するために、構成要素の領域りはなるべく適
当な寸法、即ち多くてもそのエコーグラフ装置の分解能
と等しく寸法状めされる。第2図に示した例では、その
寸法は軸方向分解能の約半分に等しい。
い歪を防止するために、構成要素の領域りはなるべく適
当な寸法、即ち多くてもそのエコーグラフ装置の分解能
と等しく寸法状めされる。第2図に示した例では、その
寸法は軸方向分解能の約半分に等しい。
信号e(×)の拡張D (x)と侵食B(x)とが第2
図のbとCとにそれぞれ示される。これらの図から、侵
食(拡張)は次式と一致する構成要素の領域りを適して
局所極小(極大)の近くで各試料を単純に置き換えるこ
とがわかる。
図のbとCとにそれぞれ示される。これらの図から、侵
食(拡張)は次式と一致する構成要素の領域りを適して
局所極小(極大)の近くで各試料を単純に置き換えるこ
とがわかる。
X′がDに属するとき
D(x)=Max e(x+x’ )
X′がDに属するとき
E!(x)=Min e(x+x’ )これらの侵食及
び拡張動作が形態学的閉鎖及び開放動作の定義を可能に
する。この信号e (x)の閉鎖F (x)は次式によ
って与えられる。
び拡張動作が形態学的閉鎖及び開放動作の定義を可能に
する。この信号e (x)の閉鎖F (x)は次式によ
って与えられる。
F(x)= [:(60g)eg〕(x)従って、それ
が拡張へ適用された侵食に関係する。
が拡張へ適用された侵食に関係する。
これが第1図に示した形態学的フィルター302により
実行される動作であり、これは侵食手段302eによっ
て追従される拡張要素302dを経由して実現される。
実行される動作であり、これは侵食手段302eによっ
て追従される拡張要素302dを経由して実現される。
逆に開放0(x)は次式によって与えられる。
0(x)= ((eeg)○g〕(x)この開放動作は
従って拡張によって追従される侵食から成る。拡張手段
303dによって追従される侵素手段303eによって
実現される第1図の形態学的フィルター303が開放す
る形態学的フィルターである。
従って拡張によって追従される侵食から成る。拡張手段
303dによって追従される侵素手段303eによって
実現される第1図の形態学的フィルター303が開放す
る形態学的フィルターである。
第2図d及びeが第2図aに示したエコーグラフ信号の
それぞれ閉鎖及び開放を示す。図解的に言えば、開放(
又は閉鎖)が、それらの絶対的高さ(深さ)と無関係に
、最も鋭い先端(窪み)を切り捨てるその能力によって
特徴付けられる。最後に、第2図fは、開放する形態学
的フィルター303により追従される閉鎖する形態学的
フィルター302によって、第2図aのエコーグラフ信
号e (c)に実行された閉鎮−開放動作F [:0(
x)]の結果を示す。
それぞれ閉鎖及び開放を示す。図解的に言えば、開放(
又は閉鎖)が、それらの絶対的高さ(深さ)と無関係に
、最も鋭い先端(窪み)を切り捨てるその能力によって
特徴付けられる。最後に、第2図fは、開放する形態学
的フィルター303により追従される閉鎖する形態学的
フィルター302によって、第2図aのエコーグラフ信
号e (c)に実行された閉鎮−開放動作F [:0(
x)]の結果を示す。
第2図かられかるように、嚢胞と散乱背景との間の境界
はぼけさせられず、且つ映像の認識はその輪郭の質に強
く依存することが知られていることが強調されるべきで
ある。これが、例えば低域濾波器を超える、形態学的フ
ィルターの実質的な利点を表現する。
はぼけさせられず、且つ映像の認識はその輪郭の質に強
く依存することが知られていることが強調されるべきで
ある。これが、例えば低域濾波器を超える、形態学的フ
ィルターの実質的な利点を表現する。
異なる映像に実行された試験は、閉鎖する形態学的フィ
ルター302と開放する形態学的フィルター303とが
、映像品質改善の観点からスペックル低減に最も適して
いることを、出願人に立証した。
ルター302と開放する形態学的フィルター303とが
、映像品質改善の観点からスペックル低減に最も適して
いることを、出願人に立証した。
第2図dとeとの比較がこれか2つのフィルターの形式
間の相違を示す。
間の相違を示す。
閉鎖する形態学的フィルターは映像中に暗点を発生する
破壊的な妨害による変動を低減する。然し乍ら、連続的
な妨害に関連する先端は残存する。
破壊的な妨害による変動を低減する。然し乍ら、連続的
な妨害に関連する先端は残存する。
閉鎖−開放フィルターが後者の先端を切り揄てるが、構
成要素の領域の寸法よりも小さい、血管の壁のような細
かい詳細は見難い。従って、細かい構成が散乱環境内に
存在する場合には、閉鎖フィルターが一層有効であり、
一方、閉鎖−開放フィルターは、例えば輪郭検出工程の
最初の段階用に一層有益であることがわかる。
成要素の領域の寸法よりも小さい、血管の壁のような細
かい詳細は見難い。従って、細かい構成が散乱環境内に
存在する場合には、閉鎖フィルターが一層有効であり、
一方、閉鎖−開放フィルターは、例えば輪郭検出工程の
最初の段階用に一層有益であることがわかる。
第2図かられかるように、“平らな”構成要素の使用が
端縁と厳しく切り捨てられた窪み又は先端とを浸み込ま
せることによる不快な映像表現へ導く。2つの方法がこ
の望ましくない効果を除去するために提案され得る。第
1の方法は領域りにおいてg≠0 によって規定される立体的な要素の形をした構成要素の
使用である。
端縁と厳しく切り捨てられた窪み又は先端とを浸み込ま
せることによる不快な映像表現へ導く。2つの方法がこ
の望ましくない効果を除去するために提案され得る。第
1の方法は領域りにおいてg≠0 によって規定される立体的な要素の形をした構成要素の
使用である。
第3図は凸面の構成要素g (x)によって得られた、
第2図aのエコーグラフ信号の閉鎖されたF’ (x)
を示す。
第2図aのエコーグラフ信号の閉鎖されたF’ (x)
を示す。
第2の方法は、立体的構成要素の場合に、次式のごとき
侵食及び拡大動作の掛算的定義の使用である。
侵食及び拡大動作の掛算的定義の使用である。
X ′、 y/がDに属するとき
(e()g) (x、 y)=Min e(x+x’
、 y+y’ ) / g(x’ 、 y’ )x′、
y′がDに属するとき (e■g)(x、y)=Max e(x+x’ 、y+
y’ ) X g(x’ 、y’ )第4図は前述の掛
算的動作によって得られた、第2図aに示されたエコー
グラフ信号の閉鎖されたF’ (x)を示す。
、 y+y’ ) / g(x’ 、 y’ )x′、
y′がDに属するとき (e■g)(x、y)=Max e(x+x’ 、y+
y’ ) X g(x’ 、y’ )第4図は前述の掛
算的動作によって得られた、第2図aに示されたエコー
グラフ信号の閉鎖されたF’ (x)を示す。
いずれの場合にも浸み込み移行が一層進行的になること
がわかる。理論上は、第2の方法(掛算的方法として参
照された)は頂上と谷間とが平均して鋭くない低コント
ラスト領域でより良い平滑化を提供する。然し乍ら、構
成要素が適当に選ばれた場合には、実際には2つの方法
はほぼ同じ結果へ導くことを出願人は見出した。他方、
付加的な方法がもっと容易に実行され得ることは注意さ
れるべきである。
がわかる。理論上は、第2の方法(掛算的方法として参
照された)は頂上と谷間とが平均して鋭くない低コント
ラスト領域でより良い平滑化を提供する。然し乍ら、構
成要素が適当に選ばれた場合には、実際には2つの方法
はほぼ同じ結果へ導くことを出願人は見出した。他方、
付加的な方法がもっと容易に実行され得ることは注意さ
れるべきである。
1次元映像e(x)に関して発展されたすべての思考が
、2次元エコーグラフ映像e(x、y)にも適用できる
。然し乍ら、後者の場合には、幾つかの要求が構成要素
の形状に関して課せられる。
、2次元エコーグラフ映像e(x、y)にも適用できる
。然し乍ら、後者の場合には、幾つかの要求が構成要素
の形状に関して課せられる。
最も単純な構成要素はその領域りが矩形であり(正方形
状も含んで)、その矩形寸法はエコーグラフ装置の軸方
向及び横方向の分解能に適応される。軸方向分解能は横
方向分解能よりも一般に高い。
状も含んで)、その矩形寸法はエコーグラフ装置の軸方
向及び横方向の分解能に適応される。軸方向分解能は横
方向分解能よりも一般に高い。
矩形構成要素による映像の不快な表現を避けるために、
楕円形であり且つその長軸(短軸)が横方向分解能(軸
方向分解能)の半分に事実上等しい領域りを有する平ら
な構成要素の使用が有利である。
楕円形であり且つその長軸(短軸)が横方向分解能(軸
方向分解能)の半分に事実上等しい領域りを有する平ら
な構成要素の使用が有利である。
“平らな”構成要素を使用する代わりに、例えば2つの
別のガウス関数 g(x、 y) =kgax(x)gLat(y)ここ
で、gaM(x) =exp(−x”/ σ、2)gt
at (y) =exp(−y”/σ【2)の形で書か
れ得る立体的構成要素を使用することも可能である。
別のガウス関数 g(x、 y) =kgax(x)gLat(y)ここ
で、gaM(x) =exp(−x”/ σ、2)gt
at (y) =exp(−y”/σ【2)の形で書か
れ得る立体的構成要素を使用することも可能である。
この場合には領域りはガウス関数の凸部分に制限される
。
。
部分的走査を伴うエコーグラフ装置の場合には、構成要
素の方向が全体映像を横切るビームの方向を追従しなけ
ればならぬことも注意されるべきである。
素の方向が全体映像を横切るビームの方向を追従しなけ
ればならぬことも注意されるべきである。
第1図は本発明による超音波エコーグラフのブロック線
図を示し、 第2図は平らな構成要素の場合におけるエコーグラフ信
号aの、拡張b1侵食C1閉鎖d1開放e1及び閉鎮−
開放fのそれぞれグラフを示し、第3図は立体的構成要
素の場合における第2図aのエコーグラフ信号の閉鎖を
示し、 第4図は掛算的形態学の場合における第2図aのエコー
グラフ信号の閉鎖を示す。 100・・・圧電変換器 101・・・走査制御装置 102 ・・・T/Rスイッチ 200・・・伝達ステージ 300・・・受信及び処理ステージ 301・・・開放する形態学的フィルター302d・・
・拡張手段 302e・・・侵食手段 303・・・閉鎖する形態学的フィルター303d・・
・拡張手段 303e・・・侵食手段 304・・・対数圧縮増幅器 305・・・記憶及び走査変換器 306・・・表示装置 e・・・エコーグラフ信号 g・・・構成要素
図を示し、 第2図は平らな構成要素の場合におけるエコーグラフ信
号aの、拡張b1侵食C1閉鎖d1開放e1及び閉鎮−
開放fのそれぞれグラフを示し、第3図は立体的構成要
素の場合における第2図aのエコーグラフ信号の閉鎖を
示し、 第4図は掛算的形態学の場合における第2図aのエコー
グラフ信号の閉鎖を示す。 100・・・圧電変換器 101・・・走査制御装置 102 ・・・T/Rスイッチ 200・・・伝達ステージ 300・・・受信及び処理ステージ 301・・・開放する形態学的フィルター302d・・
・拡張手段 302e・・・侵食手段 303・・・閉鎖する形態学的フィルター303d・・
・拡張手段 303e・・・侵食手段 304・・・対数圧縮増幅器 305・・・記憶及び走査変換器 306・・・表示装置 e・・・エコーグラフ信号 g・・・構成要素
Claims (1)
- 【特許請求の範囲】 1、少なくとも1個の圧電変換器(100)と、超音波
ビーム伝達ステージ(200)と、及び前記圧電変換器
(100)へ帰還したエコーグラフ信号eを受信し且つ
処理する受信及び処理ステージ(300)とを具える超
音波エコーグラフ装置において、 前記受信及び処理ステージ(300)が、領域D内に構
成要素gを有し且つ前記エコーグラフ信号へ適用される
基本的な侵食手段(302e、303e)及び拡張手段
(302d、303d)によって実現される、少なくと
も1個のグレイスケール形態学的フィルター(302、
303)を具えることを特徴とする超音波エコーグラフ
装置。 2、構成要素gの領域Dの寸法が最大でも前記エコーグ
ラフ装置の分解能の半分に等しいことを特徴とする請求
項1記載の超音波エコーグラフ装置。 3、構成要素gの領域Dの寸法が前記エコーグラフ装置
の分解能の半分に実質的に等しいことを特徴とする請求
項2記載の超音波エコーグラフ装置。 4、前記形態的フィルター(302)が侵食手段(30
2e)により追従される拡張手段(302d)によって
実現される閉鎖するフィルターであることを特徴とする
請求項1〜3のいずれか1項記載の超音波エコーグラフ
装置。 5、前記閉鎖する形態的フィルター(302)が、拡大
手段(303d)により追従される侵食手段(303e
)によって実現される開放する形態学的フィルター(3
03)により追従されることを特徴とする請求項4記載
の超音波エコーグラフ装置。 6、構成要素gが領域Dでg=0で定義される平らな要
素であることを特徴とする請求項1〜5のいずれか1項
記載の超音波エコーグラフ装置。 7、構成要素gが領域Dでg≠0で定義される立体的要
素であることを特徴とする請求項1〜5のいずれか1項
記載の超音波エコーグラフ装置。 8、拡張手段(302d、303d)及び侵食手段(3
02e、303e)がそれぞれ次の動作により定義され
ることを特徴とする請求項6又は7記載の超音波エコー
グラフ装置。 x′、y′がDに属するとき (e+g)(x、y)=Max〔e(x+x′、y+y
′)+g(x′、y′)〕x′、y′がDに属するとき (e−g)(x、y)=Min〔e(x+x′、y+y
′)−g(x′、y′)〕9、拡張手段(302d、3
03d)及び侵食手段(302e、303e)がそれぞ
れ次の動作により定義されることを特徴とする請求項7
記載の超音波エコーグラフ装置。 x′、y′がDに属するとき (e+g)(x、y)=Max e(x+x′、y+y
′)×g(x′、y′)x′、y′がDに属するとき (e−g)(x、y)=Min e(x+x′、y+y
′)/g(x′、y′)10、領域Dが矩形であること
を特徴とする請求項1〜9のいずれか1項記載のエコー
グラフ装置。 11、領域Dが楕円形であることを特徴とする請求項1
〜9のいずれか1項記載のエコーグラフ装置。 12、構成要素が g(x、y)=kg_a_x(x)g_t_a_t(y
)ここでg_a_x(x)=exp(−x^2/σ_a
^2)及びg_t_a_t(x)=exp(−y^2/
σ_t^2)で与えられ、領域Dは前記ガウス関数g_
a_x(x)及びg_t_a_t(y)の凸部に制限さ
れていることを特徴とする請求項9記載のエコーグラフ
装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR8807210A FR2632071B1 (fr) | 1988-05-31 | 1988-05-31 | Echographe ultrasonore a bruit d'interference reduit |
FR8807210 | 1988-05-31 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH0225785A true JPH0225785A (ja) | 1990-01-29 |
Family
ID=9366763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP1132933A Pending JPH0225785A (ja) | 1988-05-31 | 1989-05-29 | 超音波エコーグラフ装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US4951676A (ja) |
EP (1) | EP0345852A1 (ja) |
JP (1) | JPH0225785A (ja) |
FR (1) | FR2632071B1 (ja) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3854303T2 (de) * | 1987-08-31 | 1995-12-21 | Yokogawa Medical Syst | Ultraschall-diagnosegerät. |
FR2637400B1 (fr) * | 1988-09-30 | 1990-11-09 | Labo Electronique Physique | Dispositif de traitement ameliore d'un signal echographique |
EP0596560A1 (fr) * | 1992-11-04 | 1994-05-11 | Laboratoires D'electronique Philips | Appareil d'examen de milieux par echographie ultrasonore |
US6155978A (en) * | 1998-12-09 | 2000-12-05 | General Electric Company | Three-dimensional imaging by projecting morphologically filtered pixel data |
US8021301B2 (en) * | 2003-12-26 | 2011-09-20 | Fujifilm Corporation | Ultrasonic image processing apparatus, ultrasonic image processing method and ultrasonic image processing program |
US7601119B2 (en) * | 2006-04-25 | 2009-10-13 | Hrayr Kamig Shahinian | Remote manipulator with eyeballs |
US8323182B2 (en) | 2007-12-18 | 2012-12-04 | Manohara Harish M | Endoscope and system and method of operation thereof |
EP2498667A4 (en) | 2009-11-13 | 2017-12-27 | California Institute of Technology | Stereo imaging miniature endoscope with single imaging chip and conjugated multi-bandpass filters |
FR2974220B1 (fr) | 2011-04-18 | 2013-04-12 | Michelin Soc Tech | Analyse de l'image numerique de la surface interne d'un pneumatique - traitement des points de fausse mesure |
FR2974218A1 (fr) | 2011-04-18 | 2012-10-19 | Michelin Soc Tech | Analyse de l'image numerique de la surface d'un pneumatique - traitement des points de non mesure |
FR2974219A1 (fr) * | 2011-04-18 | 2012-10-19 | Michelin Soc Tech | Analyse de l'image numerique de la surface externe d'un pneumatique - traitement des points de fausse mesure |
US9456735B2 (en) | 2012-09-27 | 2016-10-04 | Shahinian Karnig Hrayr | Multi-angle rear-viewing endoscope and method of operation thereof |
US9295375B2 (en) | 2012-09-27 | 2016-03-29 | Hrayr Karnig Shahinian | Programmable spectral source and design tool for 3D imaging using complementary bandpass filters |
JP5777650B2 (ja) * | 2013-01-29 | 2015-09-09 | 富士フイルム株式会社 | 超音波診断装置および超音波画像生成方法 |
US9861261B2 (en) | 2014-03-14 | 2018-01-09 | Hrayr Karnig Shahinian | Endoscope system and method of operation thereof |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8431374D0 (en) * | 1984-12-12 | 1985-01-23 | Bamber J C | Adaptive filtering |
US4697594A (en) * | 1985-08-21 | 1987-10-06 | North American Philips Corporation | Displaying a single parameter image |
FR2589022B1 (fr) * | 1985-10-18 | 1988-05-27 | Thomson Csf | Procede et dispositif de generation d'images a partir de signaux ultra-sonores obtenus par echographie |
-
1988
- 1988-05-31 FR FR8807210A patent/FR2632071B1/fr not_active Expired - Lifetime
-
1989
- 1989-05-04 US US07/347,553 patent/US4951676A/en not_active Expired - Fee Related
- 1989-05-25 EP EP89201334A patent/EP0345852A1/fr not_active Ceased
- 1989-05-29 JP JP1132933A patent/JPH0225785A/ja active Pending
Also Published As
Publication number | Publication date |
---|---|
US4951676A (en) | 1990-08-28 |
FR2632071B1 (fr) | 1991-03-15 |
FR2632071A1 (fr) | 1989-12-01 |
EP0345852A1 (fr) | 1989-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0225785A (ja) | 超音波エコーグラフ装置 | |
Bamber et al. | Adaptive filtering for reduction of speckle in ultrasonic pulse-echo images | |
Michailovich et al. | Despeckling of medical ultrasound images | |
Karaman et al. | An adaptive speckle suppression filter for medical ultrasonic imaging | |
DE60028952T2 (de) | Verfahren und vorrichtung zur bilderzeugung durch verwendung von scherwellen | |
TW475982B (en) | Acoustic microscope | |
US5619998A (en) | Enhanced method for reducing ultrasound speckle noise using wavelet transform | |
US7811233B2 (en) | Finite amplitude distortion-based inhomogeneous pulse echo ultrasonic imaging | |
JP4022393B2 (ja) | 超音波診断装置 | |
JPH02280739A (ja) | 超音波エコーグラフ結像装置 | |
Koo et al. | Speckle reduction with edge preservation in medical ultrasonic images using a homogeneous region growing mean filter (HRGMF) | |
KR20140040679A (ko) | 향상된 초음파 이미징 시스템의 스펙클 저감/억제를 위한 향상된 초음파 이미징 방법/기술 | |
JPH11313822A (ja) | 3次元イメ―ジング・システム並びに走査平面の動きを追跡する方法及び装置 | |
JPH07506915A (ja) | 画像処理方法及び装置 | |
US20090209858A1 (en) | System and method for ultrasonic image processing | |
Tasnim et al. | Study of speckle noise reduction from ultrasound B-mode images using different filtering techniques | |
Bamber et al. | Texture analysis and speckle reduction in medical echography | |
Argenti et al. | Speckle suppression in ultrasonic images based on undecimated wavelets | |
JPH08511184A (ja) | 超音波走査のための方法及び装置 | |
Loupas et al. | Noise reduction in ultrasonic images by digital filtering | |
US5471878A (en) | Method of suppressing grain noise in ultrasonic waveform data by 3D filtering | |
US5524623A (en) | Adaptive artifact suppression for ultrasound imaging | |
Parker et al. | Analysis of B-scan speckle reduction by resolution limited filtering | |
EP0493909A2 (en) | Image processing system for producing an attenuation map of a scanned image | |
JP7052588B2 (ja) | 超音波診断装置、超音波画像表示方法及びプログラム |