JPS59214433A - 人工信号エツジ識別・選択信号処理操作を行なう多重測定雑音低減装置 - Google Patents
人工信号エツジ識別・選択信号処理操作を行なう多重測定雑音低減装置Info
- Publication number
- JPS59214433A JPS59214433A JP59032254A JP3225484A JPS59214433A JP S59214433 A JPS59214433 A JP S59214433A JP 59032254 A JP59032254 A JP 59032254A JP 3225484 A JP3225484 A JP 3225484A JP S59214433 A JPS59214433 A JP S59214433A
- Authority
- JP
- Japan
- Prior art keywords
- video signal
- processed video
- signal
- processing
- processed
- 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
- 238000005259 measurement Methods 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 9
- 229910052740 iodine Inorganic materials 0.000 description 9
- 239000011630 iodine Substances 0.000 description 9
- 230000009467 reduction Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 210000004872 soft tissue Anatomy 0.000 description 6
- 210000004204 blood vessel Anatomy 0.000 description 5
- 238000003708 edge detection Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 5
- 239000002872 contrast media Substances 0.000 description 4
- 239000003795 chemical substances by application Substances 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000011410 subtraction method Methods 0.000 description 3
- 238000001467 acupuncture Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000002594 fluoroscopy Methods 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 241000025470 Clusia rosea Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 238000002798 spectrophotometry method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/30—Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from X-rays
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/30—Transforming light or analogous information into electric information
- H04N5/32—Transforming X-rays
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
本発明は、大雑把に言って、X線システムのような結像
システムに関するものであり、一層詳しくは、複数の測
定値から導いた映像の処理に関するものである。
システムに関するものであり、一層詳しくは、複数の測
定値から導いた映像の処理に関するものである。
多くの結像応用分野では、映像は複数の測定値の重みつ
き合計として構成される。最も基本的な例としては、X
線結像分野における最新の開発例がある。ここで、測定
値というのは、種々のエネルギあるいは時間またはそれ
ら両方についてのものである。これらの測定値は組合わ
されて特定の具体物を表わす選択像を与える。
き合計として構成される。最も基本的な例としては、X
線結像分野における最新の開発例がある。ここで、測定
値というのは、種々のエネルギあるいは時間またはそれ
ら両方についてのものである。これらの測定値は組合わ
されて特定の具体物を表わす選択像を与える。
一例として、rscienceJ 1981年の第21
4巻761〜765頁に掲載されているC、A、Mis
tretta、’A、、B、Crummy共’4rDi
agnosis of Car’diovascu
la、r Disease byDigita
l 5ubtraction Angiogr
’aphyJに記載されているようなディジタル式X線
撮影法を用いる一時減法技術がある。ここでは、血管に
ヨウ素を注入する前後に11だ4111定植を減算して
血管だけの映像を得ている。別の例としては、rMed
ical PhysicsJ 1981年第8巻第6
59〜667頁に掲載されているり、A、Lehman
nその他の共著rGeneralized Imag
eCombination in Dual K
VP Digital R,adiographV
Jに記載れているようなエネルギ選択式X線撮影法があ
る。ここでは、種々のエネルギについて得た測定値を組
合わせて特定の物体を強調あるいは除去している。第3
の例としては、rPhysics in Medi
cine & Biology]1’976年第2
1巻第733〜744頁に掲載されているR、E、Al
varez、A、Macovski共著rEner、g
y−3ele、ctive Reconstru
ctions 1nX−ray Compu
terized T。
4巻761〜765頁に掲載されているC、A、Mis
tretta、’A、、B、Crummy共’4rDi
agnosis of Car’diovascu
la、r Disease byDigita
l 5ubtraction Angiogr
’aphyJに記載されているようなディジタル式X線
撮影法を用いる一時減法技術がある。ここでは、血管に
ヨウ素を注入する前後に11だ4111定植を減算して
血管だけの映像を得ている。別の例としては、rMed
ical PhysicsJ 1981年第8巻第6
59〜667頁に掲載されているり、A、Lehman
nその他の共著rGeneralized Imag
eCombination in Dual K
VP Digital R,adiographV
Jに記載れているようなエネルギ選択式X線撮影法があ
る。ここでは、種々のエネルギについて得た測定値を組
合わせて特定の物体を強調あるいは除去している。第3
の例としては、rPhysics in Medi
cine & Biology]1’976年第2
1巻第733〜744頁に掲載されているR、E、Al
varez、A、Macovski共著rEner、g
y−3ele、ctive Reconstru
ctions 1nX−ray Compu
terized T。
m’ography」に記載されているようなエネルギ
選択式コンピュータ断層撮影法かある。ここでは、細組
かの測定を2種のエネルギ・スペクトルについて行ない
、無ひずみ′復元を杓なって各区域の原子番号および密
度を描く。第4の例として、1981年5月5日出罪の
米国出j卯6−260694やrRadi o log
yJ 1981年第141巻第828〜831頁に掲載
されているW 、 R、B r o d y著rHyb
rid 5ubtraction for Im
proved Intravenous Arte
riography」に記載されているハイブリット減
法がある。ここでは、ヨウ素の注入前後に2重のエネル
ギ測定を行なう。各2重エネルギ対を用いて軟質組織を
除去する。こうして得た減算血管映像は軟質組織の動き
がない。
選択式コンピュータ断層撮影法かある。ここでは、細組
かの測定を2種のエネルギ・スペクトルについて行ない
、無ひずみ′復元を杓なって各区域の原子番号および密
度を描く。第4の例として、1981年5月5日出罪の
米国出j卯6−260694やrRadi o log
yJ 1981年第141巻第828〜831頁に掲載
されているW 、 R、B r o d y著rHyb
rid 5ubtraction for Im
proved Intravenous Arte
riography」に記載されているハイブリット減
法がある。ここでは、ヨウ素の注入前後に2重のエネル
ギ測定を行なう。各2重エネルギ対を用いて軟質組織を
除去する。こうして得た減算血管映像は軟質組織の動き
がない。
いずれにしても、多数の測定値を組合わせて特定の物体
の性り′↓を選択する場合、結果的に得た5NR(信”
i M ’A音比)は非選択映像に比べて少ない。した
かって、特定の物体を選択することによって疾病過程の
可視化を改善するこれら種々の処理技術はSNHの低減
についての欠点を持つ。
の性り′↓を選択する場合、結果的に得た5NR(信”
i M ’A音比)は非選択映像に比べて少ない。した
かって、特定の物体を選択することによって疾病過程の
可視化を改善するこれら種々の処理技術はSNHの低減
についての欠点を持つ。
このSNHの低減は疾病部位を可視化する能力と対立す
るものである。
るものである。
ハイブリッド減法その他の多重測定システムのSNRを
改善するシステムがA 、 M a c o v s
klによって開発されたが、これは1982年9月7
Bに出願された米国#!f許出願415,334rMu
ltjple MeasurementNoise
Reducing SystemJに開ノJ\され
ている。種々の条件下で成る対象物について一連の測定
が行なわれる。測定値は組合わされてこの対象物の成る
特定の局面を選択する。
改善するシステムがA 、 M a c o v s
klによって開発されたが、これは1982年9月7
Bに出願された米国#!f許出願415,334rMu
ltjple MeasurementNoise
Reducing SystemJに開ノJ\され
ている。種々の条件下で成る対象物について一連の測定
が行なわれる。測定値は組合わされてこの対象物の成る
特定の局面を選択する。
次に、同じ測定値を種々の重みにおいて組合わせ、所望
の選択を行なうことなく低雑音映像を得る。この改良選
択映像は、それの低周波数成分と低雑音映像の高周波数
成分とを用いて選択映像を低H音映像と組合わせること
によって形成される。高周波数成分の適切な振幅を確保
するため、これらの成分に、たとえば1選択映像の導関
数対低雑音映像の導関数の比で重みをつける。この場合
、非線形限界値を用いる。この比は高周波数成分の正し
い振幅の指標となるが、当てにはならない。M a c
o v s k iシステムの限界の1つは、低雑音
信号による異質信号(ここでは「人工に号」と呼ぶ)の
導入である。
の選択を行なうことなく低雑音映像を得る。この改良選
択映像は、それの低周波数成分と低雑音映像の高周波数
成分とを用いて選択映像を低H音映像と組合わせること
によって形成される。高周波数成分の適切な振幅を確保
するため、これらの成分に、たとえば1選択映像の導関
数対低雑音映像の導関数の比で重みをつける。この場合
、非線形限界値を用いる。この比は高周波数成分の正し
い振幅の指標となるが、当てにはならない。M a c
o v s k iシステムの限界の1つは、低雑音
信号による異質信号(ここでは「人工に号」と呼ぶ)の
導入である。
本発明はM a c o v s k iシステムの改
良をI]的とする。この改良により、M a c o
v s k iシステムでは抑Rnlできない1−1立
った人下信じエツジを除去できる。簡単に言えば、映像
信号内にある人工信号エツジの位置が識別され、ディス
プレイ信号がM a c o v s k i雑音低減
システム(NR3)信号から初期選択信号Sあるいはそ
れから導いた信号に切換えられる。好ましい実施例Cは
、人工信号エツジは、水平、垂直方向において複数の微
分演算器に人工映像信号を通して空間的な勾配ヘクトル
の2つの成分を得、次に各ポイントで勾配の大きさを測
定することによって識別される。この大きさを次に限界
値検出器に与える。いくつかのエツジか関係するので、
比較的大きな限界値を用いる。この限界値を越えたとき
、したかって、NR3信吟信号立った人工信号エツジが
あることが示されるとき、選択信号Sあるいはそれから
4かれて信号をティスプレィに与える。
良をI]的とする。この改良により、M a c o
v s k iシステムでは抑Rnlできない1−1立
った人下信じエツジを除去できる。簡単に言えば、映像
信号内にある人工信号エツジの位置が識別され、ディス
プレイ信号がM a c o v s k i雑音低減
システム(NR3)信号から初期選択信号Sあるいはそ
れから導いた信号に切換えられる。好ましい実施例Cは
、人工信号エツジは、水平、垂直方向において複数の微
分演算器に人工映像信号を通して空間的な勾配ヘクトル
の2つの成分を得、次に各ポイントで勾配の大きさを測
定することによって識別される。この大きさを次に限界
値検出器に与える。いくつかのエツジか関係するので、
比較的大きな限界値を用いる。この限界値を越えたとき
、したかって、NR3信吟信号立った人工信号エツジが
あることが示されるとき、選択信号Sあるいはそれから
4かれて信号をティスプレィに与える。
したがって、本発明の目的は目立った人工信号を除去す
る多重測定雑音低減システムである。
る多重測定雑音低減システムである。
本発明の別の目的は多重測定雑音低減システムにおける
映像を改善する方法である。
映像を改善する方法である。
本発明のまた別の目的は映像信号内の人工信号のエツジ
位置を検出し、検出位価で人工信号を持たない信号に切
換えることである。
位置を検出し、検出位価で人工信号を持たない信号に切
換えることである。
本発明およびその目的、特徴は、図面に関連した以下の
説明から最も容易に明らかとなろう。
説明から最も容易に明らかとなろう。
図面を参照して、第1図は可変電圧電源3によって電力
の供給を受けるX線v8を包含するX線システムのブロ
ックダイアグラムである。X線は対象物、通常は人体を
通過して伝えられる。伝達されたX線はX線検出器4、
たとえば、蛍光倍増管、TV左カメラ−次元あるいは二
次元列のシンチレータなどに受は取られる。映像信号5
は記憶システム6に与えられ、ここには複数の像9.1
0.11.12などを格納し得る。
の供給を受けるX線v8を包含するX線システムのブロ
ックダイアグラムである。X線は対象物、通常は人体を
通過して伝えられる。伝達されたX線はX線検出器4、
たとえば、蛍光倍増管、TV左カメラ−次元あるいは二
次元列のシンチレータなどに受は取られる。映像信号5
は記憶システム6に与えられ、ここには複数の像9.1
0.11.12などを格納し得る。
一時減法では、第1の像9が格納される。注射器2を用
いてコントラスト剤を注入後、第2の像10が格納され
る。これらの像を続く処理システムで使用する。エネル
ギ選択システムでは、像9.10.11.12などが種
々のアノード電圧3に対応する種々の値のビーム・工2
、ルキで格納される。2重エネルキ・システムては、2
種類の電圧を用いる。さらに、種々のX線ビーム◆フィ
ルタ(図示せず)を加えることもある。
いてコントラスト剤を注入後、第2の像10が格納され
る。これらの像を続く処理システムで使用する。エネル
ギ選択システムでは、像9.10.11.12などが種
々のアノード電圧3に対応する種々の値のビーム・工2
、ルキで格納される。2重エネルキ・システムては、2
種類の電圧を用いる。さらに、種々のX線ビーム◆フィ
ルタ(図示せず)を加えることもある。
先に述へたB r o 、d yのハイプリント鍼法シ
ステムでは、2種類の電圧でのalll yかコンi・
ラスト剤の注入前に行なわれる。これらの測定値は9、
工0として格納される。コントラスト 後、同じ2種類の電圧で2回の付加的な測定が行なわれ
、その測定値が11.12として格納される。これら4
つの測定値を次に処理して注射器2からのコントラス]
・剤によって不透明化された血清の低竹音映像を得る。
ステムでは、2種類の電圧でのalll yかコンi・
ラスト剤の注入前に行なわれる。これらの測定値は9、
工0として格納される。コントラスト 後、同じ2種類の電圧で2回の付加的な測定が行なわれ
、その測定値が11.12として格納される。これら4
つの測定値を次に処理して注射器2からのコントラス]
・剤によって不透明化された血清の低竹音映像を得る。
第21/1を参照して、線9.10.11.12は、成
る対象物の一連のX線Will定価を表わしており、ト
述のMacovskiシステムに従って処理される。こ
れらの信号は、たとえば、透視システムにおけるTVカ
メラの出力からイ14ることかできる。この透視システ
ムでは、一連の映像か、たとえは、磁気ディスクあるい
はンリンドステート・フレーム・メモリに形成、格納さ
れる。これらの(A ’4は、異なった光子エネルギ中
スペクI・ルによってmるか、あるいはコントラスト剤
の注入の前後に111だX線映像を表わすことができる
。
る対象物の一連のX線Will定価を表わしており、ト
述のMacovskiシステムに従って処理される。こ
れらの信号は、たとえば、透視システムにおけるTVカ
メラの出力からイ14ることかできる。この透視システ
ムでは、一連の映像か、たとえは、磁気ディスクあるい
はンリンドステート・フレーム・メモリに形成、格納さ
れる。これらの(A ’4は、異なった光子エネルギ中
スペクI・ルによってmるか、あるいはコントラスト剤
の注入の前後に111だX線映像を表わすことができる
。
般に、多重測定を行なって特定の物体の隔N1あるいは
除−1(を行なう。成る簡単な例では、ヨウ素コントラ
スト剤の注入の前後に得た測定値を鍼qすることによっ
て血管映像の隔離を行なう。他の例では1種々のX線光
子エネルギ・スペクトルて測定を行ない、これらの測定
値を組合わせて特定の物体を除去あるいは強調する。
除−1(を行なう。成る簡単な例では、ヨウ素コントラ
スト剤の注入の前後に得た測定値を鍼qすることによっ
て血管映像の隔離を行なう。他の例では1種々のX線光
子エネルギ・スペクトルて測定を行ない、これらの測定
値を組合わせて特定の物体を除去あるいは強調する。
重みつき加(至)器13が′A11l定イ1ηを選ひ、
所望物体の性質を含む選択信号17を発生する。しかし
ながら、この信号17は比較的貧弱なS N Rを持っ
ていることが多い。重みつき合ら1は、通常、信号を縮
小する負の重みを伴なう。しかしなから、各独立しγこ
測定値の雑音可変条件数が加わって大きな可変条件数、
したかつて、低いSNRを与える。このSNRは、帯域
幅を減する低域フィルタ15を使用してかなり改善でき
る。しかしながら、この雑音低減のために解像度が低減
するという犠牲を払わなければならない。本発明の主目
的は、雑音を低減しながら、高周波情tV源を用いて解
像度を回復させることにある。
所望物体の性質を含む選択信号17を発生する。しかし
ながら、この信号17は比較的貧弱なS N Rを持っ
ていることが多い。重みつき合ら1は、通常、信号を縮
小する負の重みを伴なう。しかしなから、各独立しγこ
測定値の雑音可変条件数が加わって大きな可変条件数、
したかつて、低いSNRを与える。このSNRは、帯域
幅を減する低域フィルタ15を使用してかなり改善でき
る。しかしながら、この雑音低減のために解像度が低減
するという犠牲を払わなければならない。本発明の主目
的は、雑音を低減しながら、高周波情tV源を用いて解
像度を回復させることにある。
信号9.10.11.12は重みつき加3つ構造14に
も与えられる。これらの重みは映像の成る特定の成分に
ついての信号対雑音比を最大にするように選択し得る。
も与えられる。これらの重みは映像の成る特定の成分に
ついての信号対雑音比を最大にするように選択し得る。
たとえば、血管結像の場合、これらの重みは、ヨウ素コ
ントラスト剤についての信は対剥1音比を最大にするよ
うに選定される。
ントラスト剤についての信は対剥1音比を最大にするよ
うに選定される。
これにより、重みつき加算構造13についての差冗かは
っきりする。重みはヨウ素を隔離し、種々の組織成分を
除去するように選定される。理想的には、SNRを最大
にする重みは、間に入った組織によって種々の映像部分
で異なる。したがって、これらの千みは信号を操作する
につれて変化することになるニジかしながら、簡略にす
るために1屯みは一″定とし、平均SNRを最大にして
性能面での悪化を比較的少なくしてもよい。
っきりする。重みはヨウ素を隔離し、種々の組織成分を
除去するように選定される。理想的には、SNRを最大
にする重みは、間に入った組織によって種々の映像部分
で異なる。したがって、これらの千みは信号を操作する
につれて変化することになるニジかしながら、簡略にす
るために1屯みは一″定とし、平均SNRを最大にして
性能面での悪化を比較的少なくしてもよい。
高いSNR重みつき合81.信号18が高周波数成分を
抽出する高域フィルタ16にケーえられる。この高域フ
ィルタの好ましい実施例は低域フィルタ15の相補要素
である。こうして、正規化伝送(iTt号の合旧は当該
周波数の合計に等しくなる。これらのフィルタは、いず
れも、周波数領域において実施されるか、あるいはたた
み込み関数として実施され得る。
抽出する高域フィルタ16にケーえられる。この高域フ
ィルタの好ましい実施例は低域フィルタ15の相補要素
である。こうして、正規化伝送(iTt号の合旧は当該
周波数の合計に等しくなる。これらのフィルタは、いず
れも、周波数領域において実施されるか、あるいはたた
み込み関数として実施され得る。
合成器21が低周波数映像すなわち信号19を高周波数
映像すなわち信号2oと組合わせて処理映像22を得、
これをディスプレイ23に表示する。信号19は濾波作
用のために低雑音なので、処理信号22は所望信号のう
ちの低雑音のものとなり、完全な帯域幅すなわち解像度
を有する。
映像すなわち信号2oと組合わせて処理映像22を得、
これをディスプレイ23に表示する。信号19は濾波作
用のために低雑音なので、処理信号22は所望信号のう
ちの低雑音のものとなり、完全な帯域幅すなわち解像度
を有する。
映像内の人工信号の程度は合成器21の性質に依存する
。簡略化のために、合成器21をり1なる加算器として
もよい。ここで、所望選択映像は必要な低雑音および高
解像度を持つことになる。しかしなから、信号20は望
ましくない構造の高周波数成分を含んでいる。たとえば
、血信映像において、信号20は、骨や種々の軟質組織
臓器、運動人工信号のような望ましくない構造のエツジ
を含んでいる。しかしながら、これらは工・ンジだけな
ので、熟練した放射線技師であれは、それを無視し、血
管に注意を集中することができる。
。簡略化のために、合成器21をり1なる加算器として
もよい。ここで、所望選択映像は必要な低雑音および高
解像度を持つことになる。しかしなから、信号20は望
ましくない構造の高周波数成分を含んでいる。たとえば
、血信映像において、信号20は、骨や種々の軟質組織
臓器、運動人工信号のような望ましくない構造のエツジ
を含んでいる。しかしながら、これらは工・ンジだけな
ので、熟練した放射線技師であれは、それを無視し、血
管に注意を集中することができる。
このシステムは、特に、先に述へたハイブリッド減法シ
ステムに応用できる。加算器13が4つの信号に重みを
つけてヨウ素以外の何かを除去する、と共に、軟質組織
運動天下信号を減算すると仮定する。たとえは、加算器
14は、軟質組織運動人1゛信ぢについてに関係なく低
雑音一時城法を行なうようにイ菖けに重みをつけること
ができる。これは、ヨウ素の注入前には低エネルギ信号
に対して大きな正の重みを使用し、ヨウ素成分の少ない
高エネルギ侶弓−に対してはもっと小さい正の重みを使
用することによって達成できる。これらの重みはヨウ素
社入後に得た測定値については逆になり、ヨウ素以外の
何かについて一時減算を行なう。
ステムに応用できる。加算器13が4つの信号に重みを
つけてヨウ素以外の何かを除去する、と共に、軟質組織
運動天下信号を減算すると仮定する。たとえは、加算器
14は、軟質組織運動人1゛信ぢについてに関係なく低
雑音一時城法を行なうようにイ菖けに重みをつけること
ができる。これは、ヨウ素の注入前には低エネルギ信号
に対して大きな正の重みを使用し、ヨウ素成分の少ない
高エネルギ侶弓−に対してはもっと小さい正の重みを使
用することによって達成できる。これらの重みはヨウ素
社入後に得た測定値については逆になり、ヨウ素以外の
何かについて一時減算を行なう。
この場合、運動がないときには、19.20を加算する
と、完全な解像度を持つ所望の隔離血管信号が形成され
ることになる。軟質組織の存在下では、信号19はハイ
ブリッド減法について先に述へた用途において述へたよ
うに影響を受けず、連h″じし、血管のみとなる。しか
しなから、低雑音は連動人工信号を含むことになる。し
たがって、それらか加算器21において加算されたとき
、これら人に信けのエツジは存在することになる。
と、完全な解像度を持つ所望の隔離血管信号が形成され
ることになる。軟質組織の存在下では、信号19はハイ
ブリッド減法について先に述へた用途において述へたよ
うに影響を受けず、連h″じし、血管のみとなる。しか
しなから、低雑音は連動人工信号を含むことになる。し
たがって、それらか加算器21において加算されたとき
、これら人に信けのエツジは存在することになる。
本発明によれば、M a c o v s k iの多
重測定雑音低減システムを改造して、M a c o
v s k iた人工信号を、それが存在する映像区域
でパラメータ選択信号あるいはそれから心いた信号に切
換えることによって除去できる。
重測定雑音低減システムを改造して、M a c o
v s k iた人工信号を、それが存在する映像区域
でパラメータ選択信号あるいはそれから心いた信号に切
換えることによって除去できる。
本発明による装器が第3図のブロックダイアグラムに示
しである。選択信号(S)17および低雑音(M)18
は全体的に30で示すNRSフィルタ兼合成器に加えら
れる。これは第2図の機能要素15.16.21に和名
する。低雑音4を号Mまたは差(M−3)である信号D
、あるいは人工信号エツジを含む任意他の信号がエツジ
検出演算器32に与えられる。この演q器は大下佑号工
ンジを識別12、それに応答して信号Eを発生する。
しである。選択信号(S)17および低雑音(M)18
は全体的に30で示すNRSフィルタ兼合成器に加えら
れる。これは第2図の機能要素15.16.21に和名
する。低雑音4を号Mまたは差(M−3)である信号D
、あるいは人工信号エツジを含む任意他の信号がエツジ
検出演算器32に与えられる。この演q器は大下佑号工
ンジを識別12、それに応答して信号Eを発生する。
この信号Eはスイッチ36に与えられ、それtこよって
、人工信号エツジの検出時にNRSフィルタ兼合成器3
0からの信号か無視される。選択信号Sあるいは空間フ
ィルタ34から導かれた空間波波された信号が出力映像
信号Sとして′j−えられる。したがって、エツジ検出
演算器によって識別された人工信号エツジの区域におい
ては、NRSフィルタ兼合成器30で発生した人工信号
がこのNRSフィルタ〕1艷合成2:)を分離すること
によって除去される。
、人工信号エツジの検出時にNRSフィルタ兼合成器3
0からの信号か無視される。選択信号Sあるいは空間フ
ィルタ34から導かれた空間波波された信号が出力映像
信号Sとして′j−えられる。したがって、エツジ検出
演算器によって識別された人工信号エツジの区域におい
ては、NRSフィルタ兼合成器30で発生した人工信号
がこのNRSフィルタ〕1艷合成2:)を分離すること
によって除去される。
第4図は第3図のエツジ検出演算器32の成る実施例の
機能ブロックダイアグラムであり、複数の微分前算器(
あるいは、他のエツジ・フィルタ)40.41.42を
包含する。各微分演算器は水・l・、垂直方向において
人力映像信号りまたはMを受けて空間勾配ペター・ルの
2つの成分を与える。たとえは、ポイント演算器44に
よって導かれた映像信号の各ポイントでの勾配ヘクトル
の犬きざを次に限界値と比較する。たとえば、この分野
で公知のテンブレー1・・フィルタのような他の工、ジ
・マツプ技術を用いてもエツジの向き、強1隻を検出す
ることができる。
機能ブロックダイアグラムであり、複数の微分前算器(
あるいは、他のエツジ・フィルタ)40.41.42を
包含する。各微分演算器は水・l・、垂直方向において
人力映像信号りまたはMを受けて空間勾配ペター・ルの
2つの成分を与える。たとえは、ポイント演算器44に
よって導かれた映像信号の各ポイントでの勾配ヘクトル
の犬きざを次に限界値と比較する。たとえば、この分野
で公知のテンブレー1・・フィルタのような他の工、ジ
・マツプ技術を用いてもエツジの向き、強1隻を検出す
ることができる。
人丁信θの存在を検出し、選択的に雑音低減システム(
NR3)を甥視することによって目立つた人丁信シ)を
除去する改良多重測定雑音低減システムについて説明し
てきたが、空間演算器を使用して人工信号のエツジを検
出し、NR3信号とより雑音の高い、すなわちややひす
んだ、人工信号なしの信号との間で選択を行なうことに
よって、映像を改善できる。本発明を成る特定の実施例
について説−明してきたが、この説明はほんの例示であ
り、発明を限定する意図はない。特許請求の範囲に定義
した発明の真の精神、範囲から逸脱することなく種ノ2
の変更、応用が行なえることは当業者にとって明らかで
あろう。
NR3)を甥視することによって目立つた人丁信シ)を
除去する改良多重測定雑音低減システムについて説明し
てきたが、空間演算器を使用して人工信号のエツジを検
出し、NR3信号とより雑音の高い、すなわちややひす
んだ、人工信号なしの信号との間で選択を行なうことに
よって、映像を改善できる。本発明を成る特定の実施例
について説−明してきたが、この説明はほんの例示であ
り、発明を限定する意図はない。特許請求の範囲に定義
した発明の真の精神、範囲から逸脱することなく種ノ2
の変更、応用が行なえることは当業者にとって明らかで
あろう。
第1図は本発明においてイ1効な信号をイ1するための
X線システムのブロックダイアクラムである。 第2図は従来技術によるX線システムのブロンフタイア
グラムである。 第3図は本発明によるX IQ信号処理システムの機能
ブロンフタイアグラムである。 第4図は第3図のシ7テムにイーJ効なエツジ検出回路
の機能ブロックダイアグラムである。 図面において、1・・・対象物、2・・・注射器、3・
・・′屯源、4・・・X線検出器、5・・・映像信す、
6争φ・記憶システム、8・Φ・X線′ご、9.10.
11.12会中拳像、13・・・屯みつき加q器、14
・・・重みつき加算構造、15・・・低域フィルタ、1
7・・・彦択信け、18・・・高SNR重みつき合計信
号、20・・・高周波数映像、21・・・合成器、22
・拳・処理映像、30・・・NRSフィルタ兼合成器、
32・・・エツジ検出演算器、36・・・スインチ、4
0.41.42・・・微分演算器、44・・・ポイント
演算器 手続補正書(方式) 特許庁長官 志 賀 学 殿 3、補正をする者 事件との関係 出願人 4、代理人
X線システムのブロックダイアクラムである。 第2図は従来技術によるX線システムのブロンフタイア
グラムである。 第3図は本発明によるX IQ信号処理システムの機能
ブロンフタイアグラムである。 第4図は第3図のシ7テムにイーJ効なエツジ検出回路
の機能ブロックダイアグラムである。 図面において、1・・・対象物、2・・・注射器、3・
・・′屯源、4・・・X線検出器、5・・・映像信す、
6争φ・記憶システム、8・Φ・X線′ご、9.10.
11.12会中拳像、13・・・屯みつき加q器、14
・・・重みつき加算構造、15・・・低域フィルタ、1
7・・・彦択信け、18・・・高SNR重みつき合計信
号、20・・・高周波数映像、21・・・合成器、22
・拳・処理映像、30・・・NRSフィルタ兼合成器、
32・・・エツジ検出演算器、36・・・スインチ、4
0.41.42・・・微分演算器、44・・・ポイント
演算器 手続補正書(方式) 特許庁長官 志 賀 学 殿 3、補正をする者 事件との関係 出願人 4、代理人
Claims (7)
- (1)、 複数の測定値から導いた所望パラメータを表
わす勾象物の処理映像信号内の異質映像信号を低減する
ための、結像システムで用いる装置であって、 11J記複数の測定値を処理し、前記所望パラメータを
表わす第1の処理映像信号を発生する第1の処理手段と
、 +iii記+h数の測定値を処理し、前記第1処理映像
信号より大きな信−り対雑音比を有する第2の処理映像
信号を発生する第2の処理手段と、前記第1処理映像信
号を受は取り、それを低域濾波する第1のフィルタ手段
と、 前記第2処理映像信号を受は取り、それを高域濾波する
第2のフィルタ手段と、 前記第1、@2の処理映像信号を濾波した後それらを組
合わせて組合わせ映像信号を発生する手段と、 前記第2処理映像信号内の異質映像のエツジの空間的位
置を決定する手段と、 異質映像のエツジの検出に応答して前記組合わせ映像信
号から前記第1処理映像信号に基つく信号に切換える手
段と を包含することを特徴とする袋打。 - (2)、 特許請求の範囲第1項記載の装置において、
エツジの空間的位置を決定する前記手段が、前記第1、
第2の処理映像信号間の差の空間的勾配信号ベクトルの
成分を決定する微分手段を包含することを特徴とする装
置市−0 - (3)、 特許請求の範囲第1項記載の装置において、
前記第1、第2のフィルタ手段が周波数相補型であるこ
とを特徴とする装置。 - (4)、 特許請求の範囲第1項記戦の装置において、
前記81処理映像信号か前記測定値の第1の重みつき合
計から導かれ、前記第2処理映像信号が前記測定値の第
2の重みつき合計から導かれたものであることを特徴と
する装置。 - (5)、 特許請求の範囲第4項記載の装置において、
「11j記結像システムがX線システムを包含し、tj
9記測定値が種、々のX線エネルギについてのものであ
ることを特徴とする装置。 - (6)、 特許請求の範囲第4項記載の装置において、
前記結像システムがX線システムを包含し、前記測定値
が種々の時間についてのものであることを特徴とする装
置。 - (7)、 複数の測定を行なうX線結像システムにおい
て処理映像信号内の異質信号を低減する方法であって、 複数の4IJ定f1ηを処理して対象物の所望パラメー
タを表わす第1の処理映像信号を発生させる段階と、 前記複数の測定値を処理して前記第1処理映像信号より
も大きい信号対雑音比を有する第2の処理映像信号を発
生し、人工映像を導き出す段階と。 これら人工映像のエツジの空間的位置を決定する段階と
、 前記人工映像の空間的位置で前記第2処理映像信号から
前記第1処理映像信号に切換える段階と からなることを特徴とする方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US468656 | 1983-02-22 | ||
US06/468,656 US4499493A (en) | 1983-02-22 | 1983-02-22 | Multiple measurement noise reducing system using artifact edge identification and selective signal processing |
Publications (1)
Publication Number | Publication Date |
---|---|
JPS59214433A true JPS59214433A (ja) | 1984-12-04 |
Family
ID=23860696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP59032254A Pending JPS59214433A (ja) | 1983-02-22 | 1984-02-22 | 人工信号エツジ識別・選択信号処理操作を行なう多重測定雑音低減装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US4499493A (ja) |
EP (1) | EP0117154A3 (ja) |
JP (1) | JPS59214433A (ja) |
AU (1) | AU558490B2 (ja) |
CA (1) | CA1214286A (ja) |
IL (1) | IL71012A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63147437A (ja) * | 1986-11-26 | 1988-06-20 | ピカー インターナシヨナル インコーポレイテツド | 画像化装置 |
JPH05174143A (ja) * | 1991-06-10 | 1993-07-13 | General Electric Co <Ge> | 空間フィルタ |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4541106A (en) * | 1984-02-22 | 1985-09-10 | General Electric Company | Dual energy rapid switching imaging system |
US4729100A (en) * | 1984-08-28 | 1988-03-01 | Kabushiki Kaisha Toshiba | CT System which convolutes projection data with a frequency varying filter function |
GB2170373B (en) * | 1984-12-28 | 1989-03-15 | Canon Kk | Image processing apparatus |
US4678996A (en) * | 1985-05-07 | 1987-07-07 | Picker International, Inc. | Magnetic resonance imaging method |
GB2181330B (en) * | 1985-09-26 | 1990-05-09 | Toshiba Kk | X-ray inspection apparatus |
US4827528A (en) * | 1985-11-15 | 1989-05-02 | Stanford University | Error-minimizing noise-reduction system |
US4761819A (en) * | 1987-02-27 | 1988-08-02 | Picker International, Inc. | Adaptive noise reduction filter for reconstructed images |
JP2536094B2 (ja) * | 1988-10-20 | 1996-09-18 | 富士ゼロックス株式会社 | カラ―画像処理装置 |
US4918543A (en) * | 1989-01-24 | 1990-04-17 | Eastman Kodak Company | Apparatus for thresholding an image signal |
DE69026159T2 (de) * | 1989-08-02 | 1996-10-17 | Canon Kk | Bildverarbeitungsgerät |
JP2993014B2 (ja) * | 1989-08-30 | 1999-12-20 | 富士ゼロックス株式会社 | 画像処理装置の画質制御方式 |
JPH03185574A (ja) * | 1989-12-15 | 1991-08-13 | Fuji Xerox Co Ltd | 多値画像フィルタ処理装置 |
US5273040A (en) * | 1991-11-14 | 1993-12-28 | Picker International, Inc. | Measurement of vetricle volumes with cardiac MRI |
KR950003035B1 (ko) * | 1992-04-17 | 1995-03-29 | 주식회사 금성사 | 영상신호의 천이영역 강조장치 및 방법 |
US5391199A (en) * | 1993-07-20 | 1995-02-21 | Biosense, Inc. | Apparatus and method for treating cardiac arrhythmias |
US5738096A (en) * | 1993-07-20 | 1998-04-14 | Biosense, Inc. | Cardiac electromechanics |
JP3372611B2 (ja) * | 1993-10-18 | 2003-02-04 | キヤノン株式会社 | 映像伝送システム、映像処理装置及び映像処理方法 |
US5868673A (en) * | 1995-03-28 | 1999-02-09 | Sonometrics Corporation | System for carrying out surgery, biopsy and ablation of a tumor or other physical anomaly |
US5797849A (en) * | 1995-03-28 | 1998-08-25 | Sonometrics Corporation | Method for carrying out a medical procedure using a three-dimensional tracking and imaging system |
US6246898B1 (en) | 1995-03-28 | 2001-06-12 | Sonometrics Corporation | Method for carrying out a medical procedure using a three-dimensional tracking and imaging system |
US5817022A (en) * | 1995-03-28 | 1998-10-06 | Sonometrics Corporation | System for displaying a 2-D ultrasound image within a 3-D viewing environment |
US5795298A (en) * | 1995-03-28 | 1998-08-18 | Sonometrics Corporation | System for sharing electrocardiogram electrodes and transducers |
US5830144A (en) * | 1995-03-28 | 1998-11-03 | Vesely; Ivan | Tracking data sheath |
US5718241A (en) * | 1995-06-07 | 1998-02-17 | Biosense, Inc. | Apparatus and method for treating cardiac arrhythmias with no discrete target |
US5954665A (en) * | 1995-06-07 | 1999-09-21 | Biosense, Inc. | Cardiac ablation catheter using correlation measure |
JPH09261481A (ja) * | 1996-03-22 | 1997-10-03 | Toshiba Corp | 画像形成装置 |
US6125213A (en) * | 1997-02-17 | 2000-09-26 | Canon Kabushiki Kaisha | Image processing method, an image processing apparatus, and a storage medium readable by a computer |
US6019725A (en) * | 1997-03-07 | 2000-02-01 | Sonometrics Corporation | Three-dimensional tracking and imaging system |
US6490474B1 (en) | 1997-08-01 | 2002-12-03 | Cardiac Pathways Corporation | System and method for electrode localization using ultrasound |
US6950689B1 (en) * | 1998-08-03 | 2005-09-27 | Boston Scientific Scimed, Inc. | Dynamically alterable three-dimensional graphical model of a body region |
US6917697B2 (en) * | 2001-05-08 | 2005-07-12 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus to automatically determine tissue cancellation parameters in X-ray dual energy imaging |
US20050261571A1 (en) * | 2004-05-21 | 2005-11-24 | Willis Nathaniel P | 3-D ultrasound navigation during radio-frequency ablation |
US7590303B2 (en) * | 2005-09-29 | 2009-09-15 | Samsung Electronics Co., Ltd. | Image enhancement method using local illumination correction |
JP5106978B2 (ja) * | 2007-10-15 | 2012-12-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4445226A (en) * | 1981-05-05 | 1984-04-24 | The Board Of Trustees Of The Leland Stanford Junior University | Multiple-energy X-ray subtraction imaging system |
NL8105223A (nl) * | 1981-11-18 | 1983-06-16 | Philips Nv | Inrichting voor verschilbeeldbepaling. |
US4463375A (en) * | 1982-09-07 | 1984-07-31 | The Board Of Trustees Of The Leland Standford Junior University | Multiple-measurement noise-reducing system |
-
1983
- 1983-02-22 US US06/468,656 patent/US4499493A/en not_active Expired - Lifetime
-
1984
- 1984-02-17 CA CA000447734A patent/CA1214286A/en not_active Expired
- 1984-02-20 AU AU24730/84A patent/AU558490B2/en not_active Ceased
- 1984-02-21 EP EP84301118A patent/EP0117154A3/en not_active Withdrawn
- 1984-02-21 IL IL71012A patent/IL71012A/xx unknown
- 1984-02-22 JP JP59032254A patent/JPS59214433A/ja active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63147437A (ja) * | 1986-11-26 | 1988-06-20 | ピカー インターナシヨナル インコーポレイテツド | 画像化装置 |
JPH05174143A (ja) * | 1991-06-10 | 1993-07-13 | General Electric Co <Ge> | 空間フィルタ |
Also Published As
Publication number | Publication date |
---|---|
IL71012A (en) | 1988-11-15 |
AU2473084A (en) | 1984-08-30 |
IL71012A0 (en) | 1984-05-31 |
EP0117154A2 (en) | 1984-08-29 |
AU558490B2 (en) | 1987-01-29 |
EP0117154A3 (en) | 1987-02-04 |
CA1214286A (en) | 1986-11-18 |
US4499493A (en) | 1985-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPS59214433A (ja) | 人工信号エツジ識別・選択信号処理操作を行なう多重測定雑音低減装置 | |
US4503461A (en) | Multiple measurement noise reducing system using space-variant filters | |
Dong et al. | Combined iterative reconstruction and image‐domain decomposition for dual energy CT using total‐variation regularization | |
US20200118306A1 (en) | Method for processing unmatched low-dose x-ray computed tomography image using neural network and apparatus therefor | |
US4463375A (en) | Multiple-measurement noise-reducing system | |
US7599575B2 (en) | Method and apparatus for reducing noise in a sequence of fluoroscopic images | |
JPH11155834A (ja) | 医用イメージング・システム及び医用画像を形成する方法 | |
US6760401B2 (en) | Apparatus and method for processing of digital images | |
JP2001238870A (ja) | 器官の厚みを補償する方法および装置 | |
JPH02278974A (ja) | X線画像処理装置 | |
CN111598989A (zh) | 一种图像渲染参数设置方法、装置、电子设备及存储介质 | |
JPH11313802A (ja) | 構造の二次元画像化用の医療診断方法 | |
CN111915538A (zh) | 一种用于数字血管减影的图像增强方法及系统 | |
Stahl et al. | Noise-resistant weak-structure enhancement for digital radiography | |
Rehm et al. | Design and testing of artifact-suppressed adaptive histogram equalization: a contrast-enhancement technique for display of digital chest radiographs | |
Macovski et al. | Measurement-dependent filtering: a novel approach to improved SNR | |
Zain et al. | Enhancement of bone fracture image using filtering techniques | |
Sprawls | Optimizing medical image contrast, detail and noise in the digital era | |
Gang et al. | Three-material decomposition using a dual-layer flat panel detector in the presence of soft tissue motion | |
Yu et al. | Large-factor Micro-CT super-resolution of bone microstructure | |
Sawant et al. | Hybrid filtering techniques for medical image denoising | |
JPH03263982A (ja) | エネルギーサブトラクション画像の表示方法および装置 | |
KR20100121273A (ko) | 엑스선 영상 처리 장치 및 방법 | |
Park et al. | Hybrid filter based on neural networks for removing quantum noise in low-dose medical X-ray CT images | |
Qian et al. | Sparse-view CT reconstruction based on improved re-sidual network |