JPH02156939A - エコーグラフ信号処理装置 - Google Patents

エコーグラフ信号処理装置

Info

Publication number
JPH02156939A
JPH02156939A JP1253803A JP25380389A JPH02156939A JP H02156939 A JPH02156939 A JP H02156939A JP 1253803 A JP1253803 A JP 1253803A JP 25380389 A JP25380389 A JP 25380389A JP H02156939 A JPH02156939 A JP H02156939A
Authority
JP
Japan
Prior art keywords
filter
signal
high frequency
circuit
sampling
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
Application number
JP1253803A
Other languages
English (en)
Inventor
Odile Bonnefous
オディル ボンヌフ
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.)
Koninklijke Philips NV
Original Assignee
Philips Gloeilampenfabrieken NV
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 Philips Gloeilampenfabrieken NV filed Critical Philips Gloeilampenfabrieken NV
Publication of JPH02156939A publication Critical patent/JPH02156939A/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/5206Two-dimensional coordinated display of distance and direction; B-scan display

Abstract

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

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明はエコーグラフ信号S。(1)を処理する装置に
関するものである。
(従来の技術) 本発明は特に器官の映像を形成する医療用エコーグラフ
の分野に有利に用いることができる。
通常、エコーグラフ信号処理装置により解決すべき技術
的な問題は検査される媒体のできるだけ正確な像を得る
必要があること、即ち、輪郭およびこれを含む反射壁を
正確に得る必要があることである。
かかる−船釣な技術問題に対する解決策は既に種々提案
されている。特にエコーグラフ像の形成に関する種々の
解析は、ジャーナル ドフィジック アプリケ18(1
983)の第527−556頁にエム フィフクが発表
した論文“超音波映像”に記載されている。これらの技
術により満足な結果がしばしば得られるも、これら既知
の解決策にはある制限が課せられている。実際上、器官
の輪郭の測定および低輝度物体の検出は、特に、圧電ト
ランスデユーサにより放出された超音波のコヒーレント
特性に関し分析された媒体に含まれるマルチスキャッタ
点によって生じるエコーの建設的および破懐的干渉によ
り妨害されるようになる。この寄生現象は、光学分野に
おいて“スペックル”として既知であり、かつ、レーザ
放出でしばしば遭遇する“先位”が発生する現象に類似
している。超音波分野におけるこの干渉雑音の主欠点は
エコーグラフ像の品質の劣化、特に低輝度区域における
可視性の損失に起因するとともに慣例の輪郭検出処理を
無効にする広周波数帯域および強力な可変性に起、因す
る。
(発明が解決しようとする課題) これがため、“スペックル”の低減はエコーグラフ像の
品質を改善する決定的なファクタであることは明ら・か
である。この目的のため、次に示す種々の方法が提案さ
れている。
“スペックル”非相関を達成するために、アクイジショ
ン前に信号処理を含む方法。これは、例えば、空間構成
である(IEEE )ランザクチオンソニック ウルト
ラソニック、SO251−6(1978)に記載された
C、 B、バークハートの論文参照)。一般に、この方
法によれば品質の良好なエコーグラフ像を得ることがで
きる。しかし、この場合には一層複雑な電子処理装置を
用いる必要のある欠点がある。
更に、アクイジション時間が比較的長く、像を実時間で
形成するのが困難となる。
また、例えば、スペックルによる変動を有効に減少する
が輪郭の拡大により像にブルーアを発生する低域通過フ
ィルタを用いて像を形成した後に適用される平滑方法。
従って、本発明は上述した技術的な問題を解決し、信号
を後処理するも個別の像の輪郭を保持しながら、“スペ
ックル”を充分に減少し得るようにした上述した種類の
エコーグラフ信号処理装置を提供することを目的とする
(課面を解決するための手段) 本発明エコーグラフ信号処理装置は、エコーグラフ信号
S0(t)のN個の個別の値C−ok(ここにkは1゜
2、・・・、 N)と、p個の縦続接続の分解装置とを
具え、各分解装置は、 スケーリングフィルタと称されるフィルタh(t)によ
り入力信号Sj−1(t)の相関を行う低域通過フィル
タと、前記相関の結果のkに対し2つの値のうちの1つ
のみを保持して全体でN/2j値自ばここにに=t、 
1+2j、 1+2x2j、 1+3x2j、 ・)が
順序jの値S0(t)低周波成分と称される値Sj(t
)を表わすデシメータを具える低周波伝送チャネルと、
ウェーブレットフィルタと称されるフィルタg(1)に
より入力信号Sj−1(t)の相関を行う高周波フィル
タと、この相関の結果のkに対し2つの値のうちの1つ
のみを保持して全体でN/2 ’値Cj、(ここにに=
1.1+2j、 l+2x 2j、 l+3X 2j、
・・・)がメモリに記憶され順序jの値S。(1)の高
周波成分と称されるり。
(1)を表わし、これら低周波フィルタおよび高周波フ
ィルタを直交ミラーフィルタとする高周波減少チャネル
とを具え、 他に、p個の縦続接続の再構成装置を具え、各再構成装
置は C′j−/L−1と順序jの再構成された信号であるs
′j(t)の値を表わす所定のスレシホルドAj、C’
ikとを比較し、全部のkに対しIC’ +b/11t
hl<A+または、逆の場合にN/2jとなる際、Wj
kのN/2 ’の値を発生する比較回路と、前記比較回
路により供給される信号の2つの順次の値間に零を挿入
するサンプリングを2倍とする回路と、この2倍回路に
より供給される信号のコンボリューションを行って、前
記ウェーブレットフィルタg(t)でサンプリングを2
倍にする高周波コンボリューション回路とを順次に具え
る高周波エンハンスメントチャネルと、C′jkの2つ
の連続値間に零を挿入するサンプリングを2倍にする回
路と、この2倍回路により供給される信号のコンボリュ
ーションを行い前記スケーリングフィルタh(t)でサ
ンプリングを2倍にする低周波コンボリューション回路
と、前記高周波コンボリューション回路及び低周波コン
ボリューション回路の出力信号の和を形成することによ
り順序j−1の再構成信号S′j−t(oを供給する加
算器とを有する実再構成チャネルとを具えることを特徴
とする。
直交ミラーフィルタは、例えば、S、G、マレ−の論文
“多重精細度信号分解の理論二ウェーブレット表示”か
ら既知である。周波数fの領域におけるかかるフィルタ
の2つをH(f )およびG(f)とすると、これらフ
ィルタはじじきの関係を満足する。
G(f)=exp(−2yr f)H(f+Fe/2)
       (1)H2(f)+G2(f)=1  
           (2)ここにFeはサンプリン
グ周波数である。
式(2)から明らかなように、直交ミラーフィルタはこ
れらが完全にエネルギーに関しては相補的であり、フィ
ルタにより処理された後にも全信号が保持されていると
云う興味深い特性を有する。これがため、本発明信号処
理装置に関しては概算なくしかも情報損失な(極めて厳
密に作動することは明らかである。
時間領域・では、フィルタg(t)およびh(t)が上
記フィルタG(f)およびH(f)に相当する。フィル
タH(f)を低域通過フィルタとする場合にはこれが測
定された物体の寸法を表わし、従って、I(f)および
h(t)はスケーリングフィルタと称される。逆に、フ
ィルタG(f)を高域通過フィルタとする場合には、こ
れが測定された像の高周波成分を表わし、従ってG(f
)およびg(t)はウェーブレットフィルタと称される
ようになる。
スケーリングフィルタ及びウェーブレフトフィルタの作
動を要約すると次の通りである。解析すべきエコーグラ
フ信号5o(i)はN値C6,lで所定のサンプリング
ステップによりサンプリングする。また、フィルタg(
t)およびh(t)も所定数のn値(n<N)で同一の
ステップを用いてサンプリングする。
次いで上式(1)を次・に示すように書換える。
g+(−1)’−’h+−+ 低周波フィルタおよびデシメータはスケ−りングフィル
タと相俟って相関操作を行い、次いで得られた2つの値
のうちの一方を抑圧する。この作動を次式で示す。
CIK=ΣlcO+ 1+2j[hl この場合の加算は1のn値として行う。
一般に、順序jに対し次式が成立する。
C+g=Σ1ci−1n 1+2+1hここにCIkは
分解順序jのスケーリングパラメータである。同様に、
高周波フィルタおよびデシメータの全効果を順序jのウ
ェーブレットW1.のパラメータを計算し、これを次式
に示す。
IFIK=Σ+Cl−1・1+2にglこれがため、各
々が分解装置を通過する際CIKで表わされる信号Sj
(t)の最高周波数成分が減少し、抽出された周波数は
毎回低くなる。
スケーリングファクタC11lおよびウェーブレットパ
ラメータfHはフィルタh(t)およびg(t)との相
関による量の場合とは相違するものとして解釈すること
ができる。実際上、上述した直交ミラーフィルタ理論お
よびアイ ドーベチイによる論文“オルソノーマルベー
シスオブコンパクトリ サポーテッドウエーブレッド”
から明らかなように、スケーリングフィルタh(t)及
びウェーブレットフィルタg(t)に対し、これらフィ
ルタg(t)及びh(t)に関連するスケーリング関数
およびウェーブレット関数と称される次式で示す関連す
る関数φ(1)およびψ(1)が存在する。
φ(t)=  xh(2jt) ψ(t)= g (2t) Xφ(2t)ここに×はコ
ンボリューション積を示す。
スケーリング関数φ(1)およびウェーブレッド関数ψ
(1)に基づき、オルソノーマライズドベースは次式で
示す関数φ1に(1)およびψ+K(t)に対し再構成
することができる。
φIK(t)・1/□Iφ((t−2jk)/2j) 
 およびψ++1(tM/v”l”ψ((t−2jk)
/2j)これがため、次式で示す関数C11lがベース
関数φIK(1)の信号S0(t)の投影であることは
明らかである。
c+x=<S0(t)、  φ+x(t)〉=fSo(
t)φ+x(t)dt上式は次式で示すように書換える
ことができる。
C+に=<Sj−1(t)、φ五に(t)>=f Sj
−1(t)φ、に(t)dt同様に、次式で示す関数W
IKが関数ψ+g(t)の信号S0(t)の投影である
ことは明らかである。
W+に=<S0(t)、  ψ +K(t)>=J S
0(t)ψ +a(t)dt及び W、g(Sl−1(t)、   ψ+x(t)>=fs
+−1(t)ψ+x(t)ctt信号の再構成状態は分
解状態に対し逆となる。
各再構成装置では、再構成される信号は・分解中メモリ
に記憶される高周波成分W、Kにより増強されるように
なる。しかし、このエンハンスメントは次に示すような
条件となる。即ち、各kに対し比Ic’ jK/WjK
 lは順序j−1の再構成信号を表わすスレシホルドA
1.C’lxよりも小さくする必要がある。
所定再構成状態において高周波の増大が行われない場合
にはC/、はC4とは相違する。その理由は上記条件が
満足されなくなるからである。換言すれば、各高周波増
大チャネルによって、全部のkに対し逆の場合にIcI
K/WIKI<A赤、LK”0であるとjIK’”Wi
llで与えられるパラメータl、に/を供給するように
なる。周波数のこの条件再構成は本質的なものとなる。
実・除土この場合、エコーグラフ信号で行われる処理作
動の関数として所望のように個別に調整され得るp個の
スレシホルドA、が存在する。
例えば、高周波デイテール(順序j=1)の除去および
輪郭(例えば順序4)の増強に対しては低いスレシホル
ドA1および高いスレシホルドA4を固定する。これが
ため、その構造が一般に順序j=1に相当するスペック
ルは充分選択的に減少し得るようになる。“スペックル
“に対してはスレシホルドA1を推定することができる
。実際上、比I C’ IK/WIKは信号対雑音比に
よって解析することができ、スペックルの信号対雑音比
が1.91となることは理論上から明らかである。従っ
て、スペックルを除去するためには1.91または2の
スレシホルドA1を適用し得るようにする。一般に、エ
コーグラフ像の物体は次に示すよ、うにその順序jの関
数として分類することができる。
スペックル   順序   j=1 壁部            1 組織            2 嚢胞            2.3 解剖物体          3.4 嚢胞のような大きな物体を解析する必要がある場合には
、例えば、A+”2. A2”3. A3”4のように
、低いスレシホルドを低い順序で置き、高いスレシホル
ドを高い順序で置(ようにする。しかし、壁部を解析す
る必要がある場合には順序1の高いスレシホルドを適用
し得るようにする。
2倍サンプリングおよびスケーリングフィルタh(t)
による低周波コンボリューションの接合効果は係数C’
l+l+IFを次式で示すkの群に従って規定された係
数C/、に変換することによって表わすことができる。
C/、□、=ΣIC’ i+l+ K+1h2j   
k偶数CZ2に+−ΣIc’ ++l+ K+I+1h
2j+l   k奇数同様に、2倍サンプリングおよび
ウェーブレットフィルタg(t)による高周波コンボリ
ューションの効果は次式で表わすことができる。
W′j2K”Σ+W1やl+に+Igg+   k偶数
W′j2に+1=ΣtW1やI+に+lやtgz++t
   k奇数これがた・め、順序jの再構成信号s′j
(t)はベース関数φ+x(t)およびψ+x(t)で
次式に示すように展開させることができる。
S’1(t)=Σxc’I+1.gφ++t、g(t)
+ΣKW’lやINKψIやt、x(t)後者の式によ
ってスケーリングベース関数およびウェーブレッドベー
ス関数に物理的な量を与えることができ、従って、再構
成信号を分解する機能を得ることができる。この重要な
特徴によって、構成を正確に記載することなく、直交ミ
ラーフィルタのかかる特性以外に今まで課せられていな
かったフィルタg(t)およびh(t)の選択基準を確
立することができる。初期エコーグラフ信号S0(t)
の出現が既知である場合には、この信号S0(t)の包
絡線にできるだけ近似したスケーリング関数φ(1)を
選択する。ウェーブレット関数ψ(1)は次式で示す直
交関数前を適用することによってこれから導出すること
ができる。
くφ(t)、ψ(t)>=J’ φ(1)ψ(t)dt
=0ここにφ(1)及びψ(1)は関数φ+x(t)及
びψ+x(t)の構成に対するベースとする。フィルタ
g(t)およびh(t)の形状はφ(1)及びψ(1)
と、g(t)及びh(t)とを組合わせることによって
上述した式から抽出することができる。
(実施例) 図面につき本発明の詳細な説明する。
第1図のブロック図の第1部分に示す本発明エコーグラ
フ信号S0(t)の処理装置は入力信号S0(t)のN
個の個別の値C0K(kは1,2.・・・、N)を供給
するサンプリング装置100を具える。次いで、サンプ
リングされた信号50(t)はp個の縦続接続された分
解装置2001によって処理する。各分解装置は2つの
チャネル、即ち、低周波伝送チャネルおよび高周波縮減
チャネルを具え、p個の順次の低周波伝送チャネルを互
いに直列に接続する。各低周波伝送チャネルは低周波フ
ィルタ2011を具え、このフィルタはスケーリングフ
ィルタと称されるフィルタh(t)と相俟って、順序j
−1の信号S0(t)の低周波成分である入力信号5l
−1(t)の相関を行う。上述したように、フィルタh
(t)は好適に選択されたスケーリング関数・φ(1)
により構成するとともに多数のn点(nは4〜20とす
る)で同一のサンプリングステップにより信号S0(t
)としてサンプリングし、従ってフィルタh(t)を容
易に実現することができる。
従って、フィルタ201.によって係数CIKを出力し
、そのデシメータ202.によって2つの値のうちの1
つのみを保持する。これがため、分解ユニット200+
の低周波伝送チャネルによって順序jの信号S。
(1)の低周波成分である5i(t)を表わすN/2 
’係数CjKを供給し、ここにkは1.1+2j、 1
+2x2j、 1+3x2j、・・・とする。これがた
め、p個の分解ユニットの数は有限となり、多くとも1
ogN/log2となる。これから明らかなようにpは
信号の後の再構成に悪影響を与えることな(、任意の数
に限定することができる。
同様に、高周波減衰チャネルはウェーブレッドフィルタ
g(t)による入力信号5t−1(t)の相関を行う高
周波フィルタを具え、フィルタh(t)およびg(t)
を直交ミラーフィルタとする。デシメータ204Iによ
り相関結果の2つの値のうちの1つのみを保持し、従っ
て全体でn/2 ’の値WjKによって順序jの信号S
。(1)の高周波成分であるDi(t)を表わす。一方
の分解ユニットから他方の分解ユニットを用いる係数C
l11に対し、後者の使用に対しメモリ300に係数f
lKを記憶する。この分解を順序pまで継続する。
第2図は本発明処理装置の第2部分を示す。この第2部
分は信号の再構成に関連する。分解装置に対し対称とな
るように、この処理装置にはp個の縦続接続の再構成ユ
ニット400+を設ける。各再i成エニット400Iは
2つのチャネル、即ち、高周波エンハンスメントチャネ
ルおよび実再構成チャネルを具え、p個の連続再構成チ
ャネルは互いに直列に接続する。各高周波エンハンスメ
ントチャネルは:!t I C’ IX/Wlに1を所
定のスレシホルドA、と比較する。ここにC/)には順
序jの再構成された入力信号S’j(t)の値を示す。
全部のkに対しlc’lK/WIKIがA、よりも小さ
い場合には、比較回路4011によってメモリ300に
記憶されたN/2 ’個の係数Wいを出力する。逆比較
の結果を得る場合には比較回路4011によって零値の
みを発生する。比較回路4011により供給される信号
の2つの連続値間に零を挿入することによりサンプリン
グを2倍にする回路402、によってサンプルの数をN
/2jからN/2 ’ −’に変化させるようにする。
最後に高周波コンボリューション回路403Iによって
ウェーブレットフィルタg(t)によりサンプリングを
2倍にする回路402jで供給される信号のコンボリュ
ーションを達成する。
又、実再構成チャネルはCIKの2つの連続値間に零を
挿入するサンプリングを2倍にする回路404Iを具え
る。このサンプリングを2倍にする回路4041からの
信号を低周波コンボリューション回路405、によって
スケーリングフィルタh(t)でコンボリューションす
る。最後に、加算器406.によって、高周波および低
周波コンボリューション回路403i、405、の出力
信号の和を形成することにより順序j−1の再構成され
た信号s′j−1(t)を供給する。各再構成回路40
0.により再注入された周波数は毎回高いか又は低くな
る。従って、信号S′。(1)が得られるまで再構成処
理を継続する。全部の係数WIKが係数C′、により再
構成されると C’ i K ” Ci K * およびS′。(t)
 = S0(t)が得られるようになる。
第3aおよび3b図はエコーグラフィック信号の分解お
よび再構成に用い得るスケーリングおよびウェーブレッ
ト関数φ(1)およびψ(1)の例を示す。これら関数
は1.ドウベッチイスによる論文(池の多くの例をも含
む)から導出することができる。
第4図は第3図の関数φ(1)およびψ(1)によるエ
コーグラフィック信号S0(t)に対し順序j=3まで
の連続再構成・ステップを示す。信号Sj(t)の高周
波数は信号D′j(t)に保持されるようになる。
第5図はスレシホルドAI”2. A2”3. A3”
4が供給される場合に第4図の信号Sj(t)およびD
j(t)を基にした再構成の結果S′。(1)を太い実
線で示す。初期信号30(t)は細い実線で示す。図面
から明らかなように、信号S′。(1)はスペックルを
含まないが、例えば壁部を通過する高周波数は含む。
【図面の簡単な説明】
第1図は本発明信号処理装置の第1部分の構成を示すブ
ロック図、 第2図は同じく信号処理装置の第2部分の構成を示すブ
ロック図、 第3aおよび3b図はスケーリング関数φ(1)および
ウェーブレット関数ψ(1)の例を示す説明図、第4図
は第3図の関数により信号の順序3までの連続分解状態
を示す説明図、 第5図は第4図の信号S0(t)の再構成の例を示す説
明図である。 100 ・・・ 200 ・・・ 201 ・・・ 202.204 300  ・・・ 400 ・・・ 401 ・・・ 402 ・・・ サンプリング装置 分解ユニット 低周波フィルタ ・・・ デシメータ メモリ 再構成ユニット 比較回路 サンプリング2倍回路 高周波コンボリューション回路 404.405 低周波コンボリューション回路 加算器

Claims (1)

  1. 【特許請求の範囲】 1、エコーグラフ信号S_0(t)のN個の個別の値C
    _0_k(ここにには1、2、・・・、N)と、p個の
    縦続接続の分解装置(200_j、j=1、2、・・・
    、p)とを具え、各分解装置(200j)は、 スケーリングフィルタと称されるフィルタ h(t)により入力信号S_j_−_1(t)の相関を
    行う低域通過フィルタ(201_j)と、前記相関の結
    果のkに対し2つの値のうちの1つのみを保持して全体
    でN/2^j値C_j_k(ここにk=1、1+2^j
    、1+2×2^j、1+3×2^j、・・・)が順序j
    の値S_0(t)低周波成分と称される値S_j(t)
    を表わすデシメータ(202_j)を具える低周波伝送
    チャネルと、 ウェーブレットフィルタと称されるフィル タg(t)により入力信号S_j_−_1(t)の相関
    を行う高周波フィルタ(203_j)と、この相関の結
    果のkに対し2つの値のうちの1つのみを保持して全体
    でN/2^j値C_j_k(ここにk=1、1+2^j
    、1+2×2^j、1+3×2^j、・・・)がメモリ
    (300)に記憶され順序jの値S_0(t)の高周波
    成分と称されるD_j(t)を表わし、これら低周波フ
    ィルタ(201_j)および高周波フィルタ(203_
    j)を直交ミラーフィルタとする高周波減少チャネルと
    を具え、 他に、p個の縦続接続の再構成装置(400_j)を具
    え、各再構成装置(400_j)は |C′_j_k/W_j_k|と順序jの再構成された
    信号であるS′_j(t)の値を表わす所定のスレシホ
    ルドA_j、C′_j_kとを比較し、全部のkに対し
    |C′_j_k/W_j_k|<A_jまたは、逆の場
    合にN/2^jとなる際、W_j_kのN/2^jの値
    を発生する比較回路(401_j)と、前記比較回路(
    401_j)により供給される信号の2つの順次の値間
    に零を挿入するサンプリングを2倍とする回路(402
    _j)と、この2倍回路(402_j)により供給され
    る信号のコンボリューションを行って、前記ウェーブレ
    ットフィルタg(t)でサンプリングを2倍にする高周
    波コンボリューション回路(403_j)とを順次に具
    える高周波エンハンスメントチャネルと、C′_j_k
    の2つの連続値間に零を挿入するサンプリングを2倍に
    する回路(404_j)と、この2倍回路(404_j
    )により供給される信号のコンボリューションを行い前
    記スケーリングフィルタh(t)でサンプリングを2倍
    にする低周波コンボリューション回路(405_j)と
    、前記高周波コンボリューション回路および低周波コン
    ボリューション回路(403_j、405_j)の出力
    信号の和を形成することにより順序j−1の再構成信号
    S′_j_−_1(t)を供給する加算器(406_j
    )とを有する実再構成チャネルとを具えることを特徴と
    するエコーグラフ信号処理装置。
JP1253803A 1988-09-30 1989-09-30 エコーグラフ信号処理装置 Pending JPH02156939A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR8812804A FR2637400B1 (fr) 1988-09-30 1988-09-30 Dispositif de traitement ameliore d'un signal echographique
FR8812804 1988-09-30

Publications (1)

Publication Number Publication Date
JPH02156939A true JPH02156939A (ja) 1990-06-15

Family

ID=9370567

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1253803A Pending JPH02156939A (ja) 1988-09-30 1989-09-30 エコーグラフ信号処理装置

Country Status (5)

Country Link
US (1) US5000183A (ja)
EP (1) EP0361606B1 (ja)
JP (1) JPH02156939A (ja)
DE (1) DE68920102T2 (ja)
FR (1) FR2637400B1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05288727A (ja) * 1992-04-10 1993-11-02 Nippon Butsuri Tanko Kk ウェーブレット変換を用いた打音調査方法
JP2006116307A (ja) * 2004-10-18 2006-05-11 Medison Co Ltd 後処理を通じた2次元超音波映像の画質改善方法及び装置

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5090412A (en) * 1987-08-31 1992-02-25 Yokogawa Medical Systems, Limited Ultrasonic diagnosis apparatus
FR2650071B1 (fr) * 1989-07-20 1991-09-27 Asulab Sa Procede de traitement d'un signal electrique
EP0456028B1 (en) * 1990-04-27 1994-12-28 Mitsubishi Denki Kabushiki Kaisha Inspection apparatus
US5384725A (en) * 1990-05-18 1995-01-24 Yale University Method and apparatus for encoding and decoding using wavelet-packets
WO1991018362A1 (en) * 1990-05-18 1991-11-28 Yale University Method and apparatus for machine processing using wavelet/wavelet-packet bases
US5105814A (en) * 1990-08-15 1992-04-21 Hewlett-Packard Company Method of transforming a multi-beam ultrasonic image
JPH07168809A (ja) * 1993-03-30 1995-07-04 Klics Ltd ウェーブレット変換方法及びウェーブレット変換回路
US5546477A (en) * 1993-03-30 1996-08-13 Klics, Inc. Data compression and decompression
US5481269A (en) * 1994-05-27 1996-01-02 Westinghouse Electric Corp. General frame wavelet classifier
US5748786A (en) 1994-09-21 1998-05-05 Ricoh Company, Ltd. Apparatus for compression using reversible embedded wavelets
JP3302229B2 (ja) 1994-09-20 2002-07-15 株式会社リコー 符号化方法、符号化/復号方法及び復号方法
US6549666B1 (en) * 1994-09-21 2003-04-15 Ricoh Company, Ltd Reversible embedded wavelet system implementation
US6195465B1 (en) 1994-09-21 2001-02-27 Ricoh Company, Ltd. Method and apparatus for compression using reversible wavelet transforms and an embedded codestream
US6229927B1 (en) 1994-09-21 2001-05-08 Ricoh Company, Ltd. Reversible embedded wavelet system implementation
US5881176A (en) 1994-09-21 1999-03-09 Ricoh Corporation Compression and decompression with wavelet style and binary style including quantization by device-dependent parser
US6873734B1 (en) 1994-09-21 2005-03-29 Ricoh Company Ltd Method and apparatus for compression using reversible wavelet transforms and an embedded codestream
US5966465A (en) * 1994-09-21 1999-10-12 Ricoh Corporation Compression/decompression using reversible embedded wavelets
US5576548A (en) * 1995-06-05 1996-11-19 University Of South Florida Nuclear imaging enhancer
US6059727A (en) * 1995-06-15 2000-05-09 The Regents Of The University Of Michigan Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound scan data
US5592171A (en) * 1995-08-17 1997-01-07 The United States Of America As Represented By The Secretary Of Commerce Wind profiling radar
US5999656A (en) * 1997-01-17 1999-12-07 Ricoh Co., Ltd. Overlapped reversible transforms for unified lossless/lossy compression
US7108663B2 (en) * 1997-02-06 2006-09-19 Exogen, Inc. Method and apparatus for cartilage growth stimulation
US7789841B2 (en) * 1997-02-06 2010-09-07 Exogen, Inc. Method and apparatus for connective tissue treatment
IT1293120B1 (it) * 1997-06-13 1999-02-11 Esaote Spa Metodo per aumentare le potenzialita' diagnostiche dei sistemi ecografici mediante l'uso di mappe spettrali in tempo reale e
US6044172A (en) * 1997-12-22 2000-03-28 Ricoh Company Ltd. Method and apparatus for reversible color conversion
AU763551B2 (en) 1998-05-06 2003-07-24 Exogen, Inc. Ultrasound bandages
US6585647B1 (en) 1998-07-21 2003-07-01 Alan A. Winder Method and means for synthetic structural imaging and volume estimation of biological tissue organs
US6213946B1 (en) * 1998-12-24 2001-04-10 Agilent Technologies, Inc. Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging
US6799141B1 (en) 1999-06-09 2004-09-28 Beamcontrol Aps Method for determining the channel gain between emitters and receivers
US6314452B1 (en) 1999-08-31 2001-11-06 Rtimage, Ltd. System and method for transmitting a digital image over a communication network
DE19954088C2 (de) * 1999-11-10 2002-01-17 Siemens Ag Verfahren zur Filterung einer digitalen Signalfolge
US20010047516A1 (en) * 2000-02-01 2001-11-29 Compaq Computer Corporation System for time shifting live streamed video-audio distributed via the internet
AU2001253563A1 (en) * 2000-04-18 2001-10-30 Rtimage Inc. System and method for the lossless progressive streaming of images over a communication network
WO2002040099A2 (en) * 2000-10-25 2002-05-23 Exogen, Inc. Transducer mounting assembly
WO2002048966A1 (en) * 2000-12-14 2002-06-20 Rtimage Inc. Three-dimensional image streaming system and method for medical images
US6898323B2 (en) 2001-02-15 2005-05-24 Ricoh Company, Ltd. Memory usage scheme for performing wavelet processing
US6950558B2 (en) * 2001-03-30 2005-09-27 Ricoh Co., Ltd. Method and apparatus for block sequential processing
US7006697B1 (en) 2001-03-30 2006-02-28 Ricoh Co., Ltd. Parallel block MQ arithmetic image compression of wavelet transform coefficients
US6859563B2 (en) 2001-03-30 2005-02-22 Ricoh Co., Ltd. Method and apparatus for decoding information using late contexts
US6895120B2 (en) 2001-03-30 2005-05-17 Ricoh Co., Ltd. 5,3 wavelet filter having three high pair and low pair filter elements with two pairs of cascaded delays
US7062101B2 (en) 2001-03-30 2006-06-13 Ricoh Co., Ltd. Method and apparatus for storing bitplanes of coefficients in a reduced size memory
JP2004528921A (ja) * 2001-05-22 2004-09-24 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ ノイズの多い超音波データにおけるコヒーレント反射体の検出
US7581027B2 (en) * 2001-06-27 2009-08-25 Ricoh Co., Ltd. JPEG 2000 for efficent imaging in a client/server environment
US7280252B1 (en) 2001-12-19 2007-10-09 Ricoh Co., Ltd. Error diffusion of multiresolutional representations
US6763322B2 (en) 2002-01-09 2004-07-13 General Electric Company Method for enhancement in screening throughput
US7095907B1 (en) 2002-01-10 2006-08-22 Ricoh Co., Ltd. Content and display device dependent creation of smaller representation of images
US7120305B2 (en) * 2002-04-16 2006-10-10 Ricoh, Co., Ltd. Adaptive nonlinear image enlargement using wavelet transform coefficients
US7653255B2 (en) 2004-06-02 2010-01-26 Adobe Systems Incorporated Image region of interest encoding
AU2005205820B2 (en) * 2004-09-04 2011-04-14 Smith & Nephew Plc Ultrasound device and method of use
US7639886B1 (en) 2004-10-04 2009-12-29 Adobe Systems Incorporated Determining scalar quantizers for a signal based on a target distortion
US9307958B2 (en) 2010-08-05 2016-04-12 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and an ultrasonic image processing apparatus
US8839672B2 (en) * 2010-10-19 2014-09-23 Board Of Regents, The University Of Texas System Combined ultrasound and photoacoustic imaging of metal objects
CN112148059B (zh) * 2020-10-12 2022-07-05 四川科陆新能电气有限公司 一种用于光伏发电站的mppt最大功率跟踪方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2556844B1 (fr) * 1983-12-14 1987-11-13 Labo Electronique Physique Appareil d'examen de milieux par echographie ultrasonore
GB8431374D0 (en) * 1984-12-12 1985-01-23 Bamber J C Adaptive filtering
NO862342L (no) * 1985-06-20 1986-12-22 Kontron Holding Ag Fremgangsmaate og apparat for fremskaffelse av et bilde.
US4817015A (en) * 1985-11-18 1989-03-28 The United States Government As Represented By The Secretary Of The Health And Human Services High speed texture discriminator for ultrasonic imaging
US4886069A (en) * 1987-12-21 1989-12-12 General Electric Company Method of, and apparatus for, obtaining a plurality of different return energy imaging beams responsive to a single excitation event
FR2632071B1 (fr) * 1988-05-31 1991-03-15 Labo Electronique Physique Echographe ultrasonore a bruit d'interference reduit

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05288727A (ja) * 1992-04-10 1993-11-02 Nippon Butsuri Tanko Kk ウェーブレット変換を用いた打音調査方法
JP2006116307A (ja) * 2004-10-18 2006-05-11 Medison Co Ltd 後処理を通じた2次元超音波映像の画質改善方法及び装置

Also Published As

Publication number Publication date
DE68920102T2 (de) 1995-06-29
EP0361606A1 (fr) 1990-04-04
DE68920102D1 (de) 1995-02-02
US5000183A (en) 1991-03-19
FR2637400A1 (fr) 1990-04-06
EP0361606B1 (fr) 1994-12-21
FR2637400B1 (fr) 1990-11-09

Similar Documents

Publication Publication Date Title
JPH02156939A (ja) エコーグラフ信号処理装置
Zhang et al. Doppler ultrasound signal denoising based on wavelet frames
Zong et al. Speckle reduction and contrast enhancement of echocardiograms via multiscale nonlinear processing
Michailovich et al. Despeckling of medical ultrasound images
Mateo et al. Finding out general tendencies in speckle noise reduction in ultrasound images
Choi et al. Analysis of wavelet-domain Wiener filters
JPH02280739A (ja) 超音波エコーグラフ結像装置
Bahoura et al. FPGA-implementation of wavelet-based denoising technique to remove power-line interference from ECG signal
Portilla et al. Image denoising via adjustment of wavelet coefficient magnitude correlation
Khan et al. Speckle suppression in medical ultrasound images through Schur decomposition
JPH0225785A (ja) 超音波エコーグラフ装置
Bhuiyan et al. New spatially adaptive wavelet-based method for the despeckling of medical ultrasound images
CN109724693A (zh) 一种基于平稳小波的融合光谱去噪方法
Elsharif et al. Wavelet decomposition–based speckle reduction method for ultrasound images by using speckle-reducing anisotropic diffusion and hybrid median
Izzetoglu et al. Wavelet domain least squares deconvolution for ultrasonic backscattered signals
Weiss et al. Wavelet-based denoising of underwater acoustic signals
Cardoso et al. Optimal wavelet estimation for data compression and noise suppression of ultrasonic NDE signals
Michailovich et al. Shift-invariant, DWT-based" projection" method for estimation of ultrasound pulse power spectrum
Bylund et al. Interactive 3D filter design for ultrasound artifact reduction
Leftwich et al. Fourier and wavelet analysis of sperm whale acoustic data in the northern Gulf of Mexico
Mohammed et al. Shift-invariant wavelet packet for signal de-noising in ultrasonic testing
Cena et al. Wavelet shrinkage of ultrasound data
Cena et al. De-noising ultrasound images with wavelets
Leftwich Denoising and Deconvolving Sperm Whale Data in the Northern Gulf of Mexico using Fourier and Wavelet Techniques
Feng et al. Implementation of de-noise DWT chip based on adaptive soft-threshold