JPS60169980A - 画像のコントラスト劣化を修正する方法及びその装置 - Google Patents
画像のコントラスト劣化を修正する方法及びその装置Info
- Publication number
- JPS60169980A JPS60169980A JP59239290A JP23929084A JPS60169980A JP S60169980 A JPS60169980 A JP S60169980A JP 59239290 A JP59239290 A JP 59239290A JP 23929084 A JP23929084 A JP 23929084A JP S60169980 A JPS60169980 A JP S60169980A
- Authority
- JP
- Japan
- Prior art keywords
- contrast
- image
- scale
- manipulating
- modulation transfer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 32
- 230000006866 deterioration Effects 0.000 title description 2
- 210000004204 blood vessel Anatomy 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 13
- 230000015556 catabolic process Effects 0.000 claims description 10
- 238000006731 degradation reaction Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 6
- 230000002411 adverse Effects 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 238000002583 angiography Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 238000012986 modification Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 description 9
- 239000002872 contrast media Substances 0.000 description 6
- 238000012937 correction Methods 0.000 description 4
- 238000005316 response function Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- 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
- H04N5/321—Transforming X-rays with video transmission of fluoroscopic images
- H04N5/325—Image enhancement, e.g. by subtraction techniques using polyenergetic X-rays
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/60—Circuit arrangements for obtaining a series of X-ray photographs or for X-ray cinematography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10121—Fluoroscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Closed-Circuit Television Systems (AREA)
- Controls And Circuits For Display Device (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明はコンピュータ促進画像化に関し、更に詳しくは
コントラスト(対比)劣化を修正するため減算画像化に
おけるスケール操作を使用1゛る方法及び装置に関する
。
コントラスト(対比)劣化を修正するため減算画像化に
おけるスケール操作を使用1゛る方法及び装置に関する
。
(従来の技術)
放射線が小孔或いは狭いスリットを通過することにより
或いは小さなピン或いは狭いバーによって減衰すると、
その結果としての画像は゛拡散し″もどの対象よりも大
きく或いは幅が広くなり、画像ど背景との間のコントラ
ストが減少り”る。これ番よエネルギー回折及び散乱し
、検知手段が精密でないこと、電子ノイズ及びディジタ
ル誤差等による複合的な結果である。
或いは小さなピン或いは狭いバーによって減衰すると、
その結果としての画像は゛拡散し″もどの対象よりも大
きく或いは幅が広くなり、画像ど背景との間のコントラ
ストが減少り”る。これ番よエネルギー回折及び散乱し
、検知手段が精密でないこと、電子ノイズ及びディジタ
ル誤差等による複合的な結果である。
゛拡散″の正確な作用形fぶは゛変調伝達作用(mod
ulation trausfer function
)”或いはMTFとして知られる。画像化装置は、画像
中の特定の拡散槻栴の夫々のMTFの組合せの結果であ
る全体的なMTFである。要するに小孔。
ulation trausfer function
)”或いはMTFとして知られる。画像化装置は、画像
中の特定の拡散槻栴の夫々のMTFの組合せの結果であ
る全体的なMTFである。要するに小孔。
スリット、ピン及びバーを通る画像は拡散し実際の対象
と比較しそのコントラストが減衰することになる。MT
Fは通常直線体に対する固定した追加分、及び対象のコ
ントラストからの付随的な演紳を限定する。従ってスリ
ット、バー。
と比較しそのコントラストが減衰することになる。MT
Fは通常直線体に対する固定した追加分、及び対象のコ
ントラストからの付随的な演紳を限定する。従ってスリ
ット、バー。
小孔或いはピンが平均的な寸法よりも十分に大きい時に
は、コントラストは実際の目的のためには現実の対象の
それと変るものではない。対象が小さい場合には画像の
平均的なコントラストは対象のコントラストと幅の両方
に比例する。
は、コントラストは実際の目的のためには現実の対象の
それと変るものではない。対象が小さい場合には画像の
平均的なコントラストは対象のコントラストと幅の両方
に比例する。
MIF拡散を修正するための多くの試みがこれまでなさ
れてきた。例えば複雑な“復元°′フィルタ法であって
、これは使用された拡散画像を再度鮮明にするものであ
る。今日使用されている復元フィルタは時間がかかり、
゛ノイズ″が入る。従ってこの方法は既に画像中に入っ
ているノイズを更に強めるものであり、従って装置の効
率と画像の品質を低下しがちである。
れてきた。例えば複雑な“復元°′フィルタ法であって
、これは使用された拡散画像を再度鮮明にするものであ
る。今日使用されている復元フィルタは時間がかかり、
゛ノイズ″が入る。従ってこの方法は既に画像中に入っ
ているノイズを更に強めるものであり、従って装置の効
率と画像の品質を低下しがちである。
画像品質向上技術についての多くの研究のうち、対象を
僅かに拡大することはMTFによる拡散により生ずるコ
ントラスト劣化の問題はどには問題ではない。
僅かに拡大することはMTFによる拡散により生ずるコ
ントラスト劣化の問題はどには問題ではない。
例えば血管系が研究されるディジタル螢光透視診断技術
及びこれに関連する研究において、血管画像が僅かに拡
大する場合における血管の検出性については殆ど問題と
ならない。これに対し、コントラストが減衰する場合に
は小さな血管を検出すること或いは識別することが殆ど
不可能になる。
及びこれに関連する研究において、血管画像が僅かに拡
大する場合における血管の検出性については殆ど問題と
ならない。これに対し、コントラストが減衰する場合に
は小さな血管を検出すること或いは識別することが殆ど
不可能になる。
MTFによる拡散問題についての一つの解決策は小さな
血管を見るための異なる“ウィンドウ(window)
’″(グレイスケール延び)を使用することである。
血管を見るための異なる“ウィンドウ(window)
’″(グレイスケール延び)を使用することである。
この゛解決策″は結果的により大きなコントラストスケ
ールを付与するものであり、これは大きな血管における
活性的な細部の損失を伴わざるをえない。またこの゛解
決策を採用する場合“正しい″ウィンドウを手動により
発見しなければならず、これは時間を消費する作業であ
る。MTF拡散問題に対する別の解決策は非直線スケー
ルを使用するものである。
ールを付与するものであり、これは大きな血管における
活性的な細部の損失を伴わざるをえない。またこの゛解
決策を採用する場合“正しい″ウィンドウを手動により
発見しなければならず、これは時間を消費する作業であ
る。MTF拡散問題に対する別の解決策は非直線スケー
ルを使用するものである。
例えば二つの別々の直線部分を含む二重ウィンドウが大
ぎい血管と小さい血管のそれぞれをカバーするようにす
る。他の非直線スケール(例えば指数スケール及び対数
スケール)も採用されたがこれらは11論的背景もなく
、関連するコントラストを歪め、細部の損失を伴いI1
1度計算を不可能にするものである。
ぎい血管と小さい血管のそれぞれをカバーするようにす
る。他の非直線スケール(例えば指数スケール及び対数
スケール)も採用されたがこれらは11論的背景もなく
、関連するコントラストを歪め、細部の損失を伴いI1
1度計算を不可能にするものである。
〈発明が解決しようとする問題点)
従って画像中の対象のMTFによる拡散に起因する、画
像のコントラスト劣化を効果的に修正する方法は長く要
望されていたものである。
像のコントラスト劣化を効果的に修正する方法は長く要
望されていたものである。
(問題点を解決するための手段)
従って本発明の目的は対象のMTFにより生ずるコント
ラスト劣化を修正するための減算画像化に使用される方
法をif供することである。
ラスト劣化を修正するための減算画像化に使用される方
法をif供することである。
このような方法は、大きな対象に悪い影響を及ぼすこと
なく小さな対象の画像品質を向上させるしのであり、修
正は大きい対象及び小さい対象の双方を見る場合にも一
つのウィンドウを使用して効果的なものでなければなら
ない。
なく小さな対象の画像品質を向上させるしのであり、修
正は大きい対象及び小さい対象の双方を見る場合にも一
つのウィンドウを使用して効果的なものでなければなら
ない。
本発明の一つの局面においてはこの方法は下記段階:
画像化される対象に均質な背景を付与し;変調伝達作用
効果で得られるような実際の画像のコン1−ラストを決
定し;そして 決定した実際の画像コントラストのスケールを操作して
、変調伝達作用劣化により影響されない理想]ントラス
トに近い決定した実際の画像コン1〜ラストを生じせし
める段1t!!i;とを含む。
効果で得られるような実際の画像のコン1−ラストを決
定し;そして 決定した実際の画像コントラストのスケールを操作して
、変調伝達作用劣化により影響されない理想]ントラス
トに近い決定した実際の画像コン1〜ラストを生じせし
める段1t!!i;とを含む。
本発明の一つの特徴は、理想的なコントラスト函数の操
作からあらかじめ作成した表を使用してスケールを操作
することである。
作からあらかじめ作成した表を使用してスケールを操作
することである。
スケールを操作する段階は、本発明の範囲を逸脱するこ
となく種々の方法により達成しうる。
となく種々の方法により達成しうる。
例えばその一つの方法は、測定コントラス]〜函数と理
想的なコントラスト函数との比較によりあらかじめ作成
した表を使用して、スケールを操作するものである。
想的なコントラスト函数との比較によりあらかじめ作成
した表を使用して、スケールを操作するものである。
表を得る1=めに必要な測定は通常厄介で時間のかかる
ものである。そのため最少用の測定のみが必要な方法を
使用することが望ましい。測定する吊を減少させるため
、例えば研究中の面上の血管のみを画像が含むようにし
コントラスト材料の密度が一定であるなどのある仮定が
成りたつなら、理論的に函数の傾斜を決定することがで
きる。適当な仮定が成り立つなら、応動函数の形を得る
ために必要な測定は、正常化及びスケーリングにも使用
されるMTF測定のみである。真の値に対する測定値の
傾斜はMTF測定から容易に得られる。
ものである。そのため最少用の測定のみが必要な方法を
使用することが望ましい。測定する吊を減少させるため
、例えば研究中の面上の血管のみを画像が含むようにし
コントラスト材料の密度が一定であるなどのある仮定が
成りたつなら、理論的に函数の傾斜を決定することがで
きる。適当な仮定が成り立つなら、応動函数の形を得る
ために必要な測定は、正常化及びスケーリングにも使用
されるMTF測定のみである。真の値に対する測定値の
傾斜はMTF測定から容易に得られる。
別の操作方法は下記段階:
既知の正確なIVjTFと同じFWHMを有する三角形
函数である標本の形状を有する大体のM1’ Fを使用
し、そして、 所定の対象寸法に対してM T F拡散によるコントラ
スト損失を大体集める段階;とを含んでなる。この三角
形函数により以下のコン1−ラストと寸法との函数が得
られる: C=A (d ’ −d ’ /3) 、ここでd丘1
゜そして C=Δ(d−d 2/3) 、ここでd7.1一方、理
想的なコントラスト函数は: C= A’波d、ここでd≧1.そして。
函数である標本の形状を有する大体のM1’ Fを使用
し、そして、 所定の対象寸法に対してM T F拡散によるコントラ
スト損失を大体集める段階;とを含んでなる。この三角
形函数により以下のコン1−ラストと寸法との函数が得
られる: C=A (d ’ −d ’ /3) 、ここでd丘1
゜そして C=Δ(d−d 2/3) 、ここでd7.1一方、理
想的なコントラスト函数は: C= A’波d、ここでd≧1.そして。
2C=A”d ” 、 ここTd ’−1式中へ=β/
r。
r。
d =F?F、/r 。
そして、βは比例定数;
夕はコントラスト材料の密度;
rはMTFの最大(FWHM=)高さの半分である全幅
;そして、 Rは血管側部の寸法(方形円筒であるとする)。
;そして、 Rは血管側部の寸法(方形円筒であるとする)。
本発明の上述特徴及びその他の特徴ならびに目的は、添
附図面と関連する以下の本発明の説明からより明らかと
なろう。
附図面と関連する以下の本発明の説明からより明らかと
なろう。
(作用及び実施例)
第1図に示したように、この方法はディジタル螢光透視
診断装置(DF)に使用されるような減算画像化装置に
特に適している。第1図の装@11は典型的なりF型の
装置である。画像取得ユニット12は、コントラスト材
料が加えられる前或いは後に投影される画像を取得する
。
診断装置(DF)に使用されるような減算画像化装置に
特に適している。第1図の装@11は典型的なりF型の
装置である。画像取得ユニット12は、コントラスト材
料が加えられる前或いは後に投影される画像を取得する
。
減算ユニット13はコントラスト材料を加える前に得ら
れる画像をコントラスト材料を加えた後に得られる画像
から減算し、均質な背景を有する対象の画像を供するよ
うにしたものである。
れる画像をコントラスト材料を加えた後に得られる画像
から減算し、均質な背景を有する対象の画像を供するよ
うにしたものである。
減算ユニット13から出力される画像データはディジタ
ル形態の場合好ましいものではあるが不可欠ではない。
ル形態の場合好ましいものではあるが不可欠ではない。
なお本明細書で言及する画像は基本的には小さな対象物
、即ち最終画像におりるコントラストを減するMIF効
果を有する例えば血管であることに注意されるべきであ
る。
、即ち最終画像におりるコントラストを減するMIF効
果を有する例えば血管であることに注意されるべきであ
る。
血管を画像化する際通常生ずるコントラスト現象を修正
りる手段が設けられている。より詳しくいうと、変調伝
達作用特性に由来するコントラスト損失を復元すること
によりコントラストを向上させるため、モニタ14上の
画像のスケールを操作する手段が設けられている。減算
ユニット13からデータをディジタル−アナログコンバ
ータ16によりアナログビデオ信号に変換する時には、
モニタ上に画像が表われる。この装置はプロセッサ17
の制御下にあり、プロセッサはなかんずくコントラスト
制御のためのスケールを操作する。第1図に示すように
演算ユニツl−18及びテーブル19は、この目的のた
めプロセッサの制御下にこの装置に使用ξれることが好
ましい。
りる手段が設けられている。より詳しくいうと、変調伝
達作用特性に由来するコントラスト損失を復元すること
によりコントラストを向上させるため、モニタ14上の
画像のスケールを操作する手段が設けられている。減算
ユニット13からデータをディジタル−アナログコンバ
ータ16によりアナログビデオ信号に変換する時には、
モニタ上に画像が表われる。この装置はプロセッサ17
の制御下にあり、プロセッサはなかんずくコントラスト
制御のためのスケールを操作する。第1図に示すように
演算ユニツl−18及びテーブル19は、この目的のた
めプロセッサの制御下にこの装置に使用ξれることが好
ましい。
第3図及び4図は考えられるスケール操作の相異なる形
式を示す。実際に使用する形式は使用する装置及び得ら
れる結果に依存する。例えば目的が小さな血管を示す画
像を得たい場合異なるスケールの操作が使用され、次に
大きな血管を画像化するのに使用される。スケール操作
のすべてにおいて、対象が比較的均質な*景の視野に入
ることがm要である。
式を示す。実際に使用する形式は使用する装置及び得ら
れる結果に依存する。例えば目的が小さな血管を示す画
像を得たい場合異なるスケールの操作が使用され、次に
大きな血管を画像化するのに使用される。スケール操作
のすべてにおいて、対象が比較的均質な*景の視野に入
ることがm要である。
以下に記載する操作のすべては所望用法の血管画像を向
上させる一方コントラスト、従って異なる寸法の血管の
可視麿に悪影響を与えるものではない。
上させる一方コントラスト、従って異なる寸法の血管の
可視麿に悪影響を与えるものではない。
背景説明どして第2(a)図は点源からはなれた拡散M
IFを示す図である。軸は点源の線画像の位置からの距
離に対する強度を表わす。
IFを示す図である。軸は点源の線画像の位置からの距
離に対する強度を表わす。
理想的な画像は点(p)或いは寸法を有さない線Jに集
中される。実際の画像はMTFにより拡散される。図示
の曲線は単一のピークを有する対称的な函数、即ちゴー
シアン(Q aussian )函数である。本発明の
範囲を逸脱することなくゴーシアン函数に加え多くの形
式の函数、例えば三角、矩形ステップ函数及びメサ函数
を使用Jることができる。
中される。実際の画像はMTFにより拡散される。図示
の曲線は単一のピークを有する対称的な函数、即ちゴー
シアン(Q aussian )函数である。本発明の
範囲を逸脱することなくゴーシアン函数に加え多くの形
式の函数、例えば三角、矩形ステップ函数及びメサ函数
を使用Jることができる。
第2(b)図は血管の側部の直線寸法Rについて測定さ
れたコントラストを真のコントラストで除痺したコント
ラストの応動函数を示し、ここで血管は矩形の横断面を
有する中空管であって、−側部が視線と平行であるもの
とする。
れたコントラストを真のコントラストで除痺したコント
ラストの応動函数を示し、ここで血管は矩形の横断面を
有する中空管であって、−側部が視線と平行であるもの
とする。
血管の側部の寸法が増大すると測定値は真の値に近づく
。小さな血管については応動函数は漸近的に側部に依存
して直線に近づく。
。小さな血管については応動函数は漸近的に側部に依存
して直線に近づく。
第2(C)図は真のコントラスト値の函数としての同じ
測定〈実際)コントラストを示す。
測定〈実際)コントラストを示す。
この図からどの所望点においてもコントラストに対する
必要な修正を推論覆ることは容易である。
必要な修正を推論覆ることは容易である。
本発明の好ましい処理においてはMTF函数は三角であ
ると考える。ゴーシアン或いはその他の形であると考え
ることはもとより本発明の範囲を逸脱するものではない
。血管が円筒形というよりも方形の中空管であり側部が
視線に対して平行であるとJoると、ディジタル減界血
管造影法においては真のコントラストは下記式=C″β
I)Rで表わされ、 式中pはコントラスト材料の密度 Rは管側部の寸法、そして βは比例定数である。
ると考える。ゴーシアン或いはその他の形であると考え
ることはもとより本発明の範囲を逸脱するものではない
。血管が円筒形というよりも方形の中空管であり側部が
視線に対して平行であるとJoると、ディジタル減界血
管造影法においては真のコントラストは下記式=C″β
I)Rで表わされ、 式中pはコントラスト材料の密度 Rは管側部の寸法、そして βは比例定数である。
MTFは応動函数をほぼ下記式:
測定されたコントラスト/実際のコントラス1式中r=
MTFのF W HM 。
MTFのF W HM 。
はMTFの形状におよぼすある範囲に依存する比例定数
である。
である。
小さな血管においては測定コントラスト/実際づ″通り
)であり、測定されたコントラストは聞のエネルギーを
示すことをあられしている。
)であり、測定されたコントラストは聞のエネルギーを
示すことをあられしている。
A−βpr、 d−R7rrという表記をすると、MT
F効果のない理想的な状態では、 C−Δ−86(すべてのRに対して)。
F効果のない理想的な状態では、 C−Δ−86(すべてのRに対して)。
しかしながら実際にはMTF効果そして三角MTFによ
り、下記の式があてはまる:C=A (d −d 2/
3)、d″51の場合C=A (’d 2/2−d ’
/3)、d、=1の場合。
り、下記の式があてはまる:C=A (d −d 2/
3)、d″51の場合C=A (’d 2/2−d ’
/3)、d、=1の場合。
理想的な函数及び実際的な函数は第2(c)図に示され
ている。MTFのrが知られていると演免ユニットはあ
らゆる測定値に対してコントラストの理論的な理想値を
計nすることができる。
ている。MTFのrが知られていると演免ユニットはあ
らゆる測定値に対してコントラストの理論的な理想値を
計nすることができる。
しかしながら、大ぎい画像の夫々のビクセル値を修正す
ることはあまりにも時間がかかる。
ることはあまりにも時間がかかる。
従って好ましい方法は対照表を計算することであり、こ
こでは入力はいずれかの測定値であり、」一連出力はお
おまかな理想値である。このような対照表は画像処理の
分野ではよく知られている。
こでは入力はいずれかの測定値であり、」一連出力はお
おまかな理想値である。このような対照表は画像処理の
分野ではよく知られている。
rが知られていないとすると、その値はraであると考
えられ、修正方法及び手段に使用される。
えられ、修正方法及び手段に使用される。
r変化の代りに「aを使用すると得られるコントラスト
値は: d71の場合 従ってr a d rである時は得られる結果は理想値
に近い。
値は: d71の場合 従ってr a d rである時は得られる結果は理想値
に近い。
小さな血管画像のコントラストを改善づるのに有用な近
似値は下記関係式を利用するものである。
似値は下記関係式を利用するものである。
この近似式から得られる結果は変数A及びdを用いると
次のように表わされる: C” =、p (1−d /6)、d71(7)jJJ
合、及び、 C”=A、d(1−d/6) d丘1の場合。
次のように表わされる: C” =、p (1−d /6)、d71(7)jJJ
合、及び、 C”=A、d(1−d/6) d丘1の場合。
これは小さな血管(d<1)の場合良好な近似値である
。この関係は第3図に示されている。
。この関係は第3図に示されている。
タベての値の範囲からのある゛″ウインドウ″かかわる
小さな血管について有用な別の近似式%式%) 1゛興味のウィンドウ(window of 1nte
rest)”を選びウィンドウ中の曲線を理想的な直線
に近づけるよう変換する。
小さな血管について有用な別の近似式%式%) 1゛興味のウィンドウ(window of 1nte
rest)”を選びウィンドウ中の曲線を理想的な直線
に近づけるよう変換する。
例えばウィンドウが閾値Tと上側値Uの間にあるとづる
“延びた1′スケールは王から0及びUから上側数値W
の変換であるから: C’ =W (C−T)/ (U−T)逆のスケールは C” =W−C’ −W (U−C)/ (U−T)で
あり、式中TはWにまたUはOに変えられる。
“延びた1′スケールは王から0及びUから上側数値W
の変換であるから: C’ =W (C−T)/ (U−T)逆のスケールは C” =W−C’ −W (U−C)/ (U−T)で
あり、式中TはWにまたUはOに変えられる。
゛逆′”の操作はこの補完的な操作である二〇”=VV
−C”=C’ 。
−C”=C’ 。
特別なウィンドウは]=Oでスタートする。すると:
C”−W/u−C(Cは範囲o−u、cLは0−W)。
1“逆″のスケールは
C” =W (U−C) /U (この場合も同じ範囲
)。
)。
この点にJ3いて、二乗法スケールが使用されている。
C” =1./WC2=W(U−C)2 /U2(同じ
範囲) 再度逆にすると: C” −W−W (U−C)/U (同じ帽1゜延びの
ない特別の場合:W−C C” −U−(U−C)2/U−20−C” /U(C
とc′f は共にO−Uの範囲)。
範囲) 再度逆にすると: C” −W−W (U−C)/U (同じ帽1゜延びの
ない特別の場合:W−C C” −U−(U−C)2/U−20−C” /U(C
とc′f は共にO−Uの範囲)。
Cが測定されたコントラストであるとすると、A及びd
を使用してコントラストは次のように書くことができよ
う: C”w2Ad ’ −2Ad ” /3.d <1.そ
して、 07郊2Ad −(2A/3−A’ /U)d ” 。
を使用してコントラストは次のように書くことができよ
う: C”w2Ad ’ −2Ad ” /3.d <1.そ
して、 07郊2Ad −(2A/3−A’ /U)d ” 。
d>1゜
これはd−0について二乗法(放物線)として出発し、
次にほぼd−64及びd−1,5(第4図参照)で理想
値と交叉する。
次にほぼd−64及びd−1,5(第4図参照)で理想
値と交叉する。
この近似(+tjは従って従前のものより広い゛ウィン
ドウ″に適している。しかしながらこのスクールは二つ
の値の回答を生じさゼるからC〉Uの場合にはならない
。注意が払われる時にはこれは、画像において値Cより
も高いUを選ぶために極めて有用である。コントラス1
へ材料が黒にみえるよう実際にはDSA画像は多くの場
合逆スケールで得られるのであるのでこれは特に有用で
ある。
ドウ″に適している。しかしながらこのスクールは二つ
の値の回答を生じさゼるからC〉Uの場合にはならない
。注意が払われる時にはこれは、画像において値Cより
も高いUを選ぶために極めて有用である。コントラス1
へ材料が黒にみえるよう実際にはDSA画像は多くの場
合逆スケールで得られるのであるのでこれは特に有用で
ある。
本発明は特定の装置及び方法について記載し!こけれど
も、記載は例を説明するものであり、本発明の範囲を限
定するものでないことが理解されるべきである。特に他
の近似法も本発明の範囲に含まれることが指摘されるべ
きである。
も、記載は例を説明するものであり、本発明の範囲を限
定するものでないことが理解されるべきである。特に他
の近似法も本発明の範囲に含まれることが指摘されるべ
きである。
第1図は、本発明によるコントラスト劣化修正方法を採
用する画像化装置のブロックダイアグラム: 第2(a)図は典型的なMTFのグラフ、第2(a)及
び2(b)図はそれぞれMTFに基づく応動曲線の相異
なる態様を示すグラフ;第3図はあるコントラスト劣化
修正方法を用いた結果を示すグラフ;そして、 第、4図は別のコントラスト劣化修正方法を用いた結果
を示すグラフである。 図中、番号11はディジタル螢光透視診断装置(DF)
、12は画像取得ユニツ1〜.13は減算ユニット、1
6はディジタル−アナログコンバータ、17はプロセッ
サ、18は演算ユニット、19はテーブルである。 特許出願人 工ルシント リミテッド 手続補正書(方式) 昭和 6CM#3 月 ノδ 日 特許庁長官 志 賀 学 1、へ (特許庁審査官 殿) 1、 事件の表示 昭和59年 特許願第239290号 2、発明の名称 画像のコントラスト劣化を修正する方法及び装置3、補
正をする者 事件との関係 特許 出願人 氏名(名称) エルシント リミテッド4、代理人 住所 東京都文京区白山5丁目14番7号早川ビル 電
話東京+14610531番(代表)氏名 (Sの)弁
理士 早 川 政 名6、 補正の対象 −
用する画像化装置のブロックダイアグラム: 第2(a)図は典型的なMTFのグラフ、第2(a)及
び2(b)図はそれぞれMTFに基づく応動曲線の相異
なる態様を示すグラフ;第3図はあるコントラスト劣化
修正方法を用いた結果を示すグラフ;そして、 第、4図は別のコントラスト劣化修正方法を用いた結果
を示すグラフである。 図中、番号11はディジタル螢光透視診断装置(DF)
、12は画像取得ユニツ1〜.13は減算ユニット、1
6はディジタル−アナログコンバータ、17はプロセッ
サ、18は演算ユニット、19はテーブルである。 特許出願人 工ルシント リミテッド 手続補正書(方式) 昭和 6CM#3 月 ノδ 日 特許庁長官 志 賀 学 1、へ (特許庁審査官 殿) 1、 事件の表示 昭和59年 特許願第239290号 2、発明の名称 画像のコントラスト劣化を修正する方法及び装置3、補
正をする者 事件との関係 特許 出願人 氏名(名称) エルシント リミテッド4、代理人 住所 東京都文京区白山5丁目14番7号早川ビル 電
話東京+14610531番(代表)氏名 (Sの)弁
理士 早 川 政 名6、 補正の対象 −
Claims (1)
- 【特許請求の範囲】 (1)画像化における変調伝達作用効果によって生ずる
コントラスト劣化を修正する方法であって、大きい対象
の画像に悪い影響を与えることなく小さな対象の画像品
質を向上させ、かつ大きい対象及び小さい対象の双方を
見るための単一のウィンドウのみを使用する時にも修正
が効果的である方法において、下記段階: 画像化される対象に均質な背景を付与し;全体のシステ
ム変調伝達作用効果で得られる画像コントラストを決定
し;そして、 決定された画像コン1〜ラストのスケールを操作して、
変調伝達作用コントラスト減衰により影響されることな
い理想的なコントラストに、決定されたコントラストを
近づけるようにする段階; を含む方法。 (2)スケールを操作する段階が、あらかじめ用意され
た対照表を使用する段階を含む特許請求の範囲第1項記
載の方法。 (3)スケールを操作すφ段階が、測定コン1−ラスト
函数に対して理想コントラスト函数を比較することによ
り前記対照表を準備することを含む特許請求の範囲第2
項記載の方法。 く4)スケールを操作する段階が、既知FWHMを使用
することを含む特許請求の範囲第1項記載の方法。 (5)スケールを操作する段階が、近似的なFWHMを
使用することを含む特許請求の範囲第1項記載の方法。 (6)三角MTFを使用づ°ることを含む特許請求の範
囲第1〜6項のいずれか1項記載の方法。 (7)方形断面を有する血管のような長い細い対象を画
像中に含む特許請求の範囲第1〜6項のいずれか1項記
載の方法。 (8)均質な背景を、付与する段階が、ディジタル減算
血管造影法を採用することを含む特許請求の範囲第1項
記載の方法。 (9)スケールを操作する段階が、逆方向に二乗するこ
とを含む特許請求の範囲第1項記載の方法。 (10)画像化における変調伝達作用効果によって生ず
るコントラスト劣化を修正する装置であって、大きい対
象の画像に悪い影響を与えることなく小さな対象の画像
品質を向上させ、かつ大きい対象及び小さい対象の双方
を見るための単一のウィンドウのみを使用する時にも修
正が効果的である装置において、下記手段: 画像化される対象に均質な背景を付与する手段: 全体のシステム変調伝達作用効果で得られる画像コント
ラストを決定する手段;及び決定された画像コントラス
トのスケールを操作して、変調伝達作用コントラスト減
衰により影響されることのない理想的なコントラストに
、決定されたコントラストを近づけるようにする手段: を含む装置。 (11)スケールを操作する手段が、あらかじめ用意さ
れた対照表を使用する手段を含む特許請求の範囲第10
項記載の装置。 (12)スケールを操作する手段が前記あらかじめ□用
意された対照表を準備する手段を含み、前記対照表を準
備する手段が測定コン(〜ラス]−函数に対して理想コ
ントラスト函数を比較する手段を含む特許請求の範囲第
11項記載の装置。 (13)スケールを操作する手段が、既知FWH’Mを
決定し使用する手段を含む特許請求の範囲第10項記載
の装置。 (14)スケールを操作する手段が、F W I−I
Mを使用し近似させる手段を含む特許請求の範囲第(1
5)三角形MTFを使用する手段を含む特許請求の範囲
第10−14項のいずれか1項記載の装置。 ゛ (16)Zj形断面を有する例えば血管のような長く細
い対象を画像中に含む特許請求の範囲第10−15項の
いずれか1項記載の装置。 〈11)均質な背景を付与する手段が、ディジタル減詐
血管造影手段を含む特許請求の範囲第10項記載の装置 (18)スケールを操作する手段が、逆方向に二乗づる
手段を含む特許請求の範囲第10項記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IL70214 | 1983-11-13 | ||
IL70214A IL70214A (en) | 1983-11-13 | 1983-11-13 | Image contrast enhancement arrangement |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS60169980A true JPS60169980A (ja) | 1985-09-03 |
JP2575612B2 JP2575612B2 (ja) | 1997-01-29 |
Family
ID=11054673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP59239290A Expired - Fee Related JP2575612B2 (ja) | 1983-11-13 | 1984-11-13 | 画像のコントラスト劣化を修正する方法及びその装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US4907288A (ja) |
JP (1) | JP2575612B2 (ja) |
DE (1) | DE3441162C2 (ja) |
FR (1) | FR2555003B1 (ja) |
IL (1) | IL70214A (ja) |
NL (1) | NL193182C (ja) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2595153B1 (fr) * | 1986-02-28 | 1990-12-07 | Thomson Cgr | Systeme d'imagerie numerique a reglage de l'echelle des gris, notamment pour la visualisation des vaisseaux sanguins |
JPH01248187A (ja) * | 1988-03-30 | 1989-10-03 | Toshiba Corp | ディスプレイシステム |
GB8813514D0 (en) * | 1988-06-08 | 1988-07-13 | Quantel Ltd | Image enhancement |
FR2633746A1 (fr) * | 1988-07-01 | 1990-01-05 | Labo Electronique Physique | Methode numerique de modification d'image et systeme de mise en oeuvre |
NL8900217A (nl) * | 1989-01-30 | 1990-08-16 | Philips Nv | Werkwijze voor het bepalen van een modulatie overdrachtfunctie van een digitaal beeldvormend systeem. |
FR2648302A1 (fr) * | 1989-06-09 | 1990-12-14 | Peugeot | Procede de traitement d'un signal composite video d'une scene a examiner |
DE69032837T2 (de) * | 1989-09-22 | 1999-05-12 | Fuji Photo Film Co., Ltd., Minami-Ashigara, Kanagawa | Verfahren und Gerät zur Energiesubtraktionsverarbeitung, Verfahren und Gerät zur Superpositionsverarbeitung, und Gerät zur Strahlungsbildauslesung |
CA2028092C (en) * | 1989-10-23 | 1996-08-06 | Yuji Tanaka | Binary image processing apparatus |
NL9001265A (nl) * | 1990-06-05 | 1992-01-02 | Philips Nv | Werkwijze voor bepaling van een modulatie-overdrachtfunctie van een digitaal beeldvormend systeem. |
US5359675A (en) * | 1990-07-26 | 1994-10-25 | Ronald Siwoff | Video spectacles |
US5125046A (en) * | 1990-07-26 | 1992-06-23 | Ronald Siwoff | Digitally enhanced imager for the visually impaired |
FR2666426B1 (fr) * | 1990-08-31 | 1994-08-19 | Gen Electric Cgr | Procede de correction des mesures de densite optique effectuees sur un film radiographique. |
FR2670979A1 (fr) * | 1990-12-21 | 1992-06-26 | Philips Electronique Lab | Procede de segmentation binaire locale d'images numerisees, par seuillage d'histogramme. |
US5271064A (en) * | 1991-06-14 | 1993-12-14 | University Of Cincinnati | Apparatus and method for smoothing regions and enhancing edges in gray scale images |
US6466689B1 (en) * | 1991-11-22 | 2002-10-15 | Arch Development Corp. | Method and system for digital radiography |
US5297550A (en) * | 1992-08-06 | 1994-03-29 | Picker International, Inc. | Background darkening of magnetic resonance angiographic images |
US5297551A (en) * | 1992-08-06 | 1994-03-29 | Picker International, Inc. | Weighted ray projection imaging for MR angiography |
GB9420654D0 (en) * | 1994-10-13 | 1994-11-30 | Kodak Ltd | A method for improving the contrast in photographic film materials |
DE69623884T2 (de) * | 1995-07-25 | 2003-05-28 | Kabushiki Kaisha Toyota Chuo Kenkyusho, Aichi | Bildverarbeitungsvorrichtung |
JPH09261490A (ja) * | 1996-03-22 | 1997-10-03 | Minolta Co Ltd | 画像形成装置 |
US6864888B1 (en) * | 1999-02-25 | 2005-03-08 | Lockheed Martin Corporation | Variable acuity rendering for a graphic image processing system |
GB9904692D0 (en) * | 1999-03-01 | 1999-04-21 | Isis Innovation | X-ray image processing |
GB0028491D0 (en) * | 2000-11-22 | 2001-01-10 | Isis Innovation | Detection of features in images |
JP2002271611A (ja) * | 2001-03-14 | 2002-09-20 | Fujitsu Ltd | 画像処理装置 |
US20030160800A1 (en) | 2002-02-22 | 2003-08-28 | Agfa-Gevaert | Multiscale gradation processing method |
TWI237497B (en) * | 2004-05-07 | 2005-08-01 | Quanta Comp Inc | High speed and high accuracy image adjust apparatus and method |
US20060056722A1 (en) * | 2004-09-14 | 2006-03-16 | Nathan Moroney | Edge preserving method and apparatus for image processing |
DE102004049227B4 (de) * | 2004-10-08 | 2007-03-01 | Yxlon International Security Gmbh | Verfahren zur Bestimmung der Lageänderung eines Objekts in einem Gepäckstück |
JP2009211313A (ja) * | 2008-03-03 | 2009-09-17 | Fujitsu Ltd | 画像処理装置、画像処理方法および画像処理プログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5323236A (en) * | 1976-08-16 | 1978-03-03 | Fuji Electric Co Ltd | Optical discrimination equipment |
JPS5834677A (ja) * | 1981-08-25 | 1983-03-01 | Matsushita Graphic Commun Syst Inc | 画情報読取装置のmtf補正方式 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4193089A (en) * | 1976-03-05 | 1980-03-11 | Picker Corporation | Television radiation imaging system and method |
US4156914A (en) * | 1977-08-18 | 1979-05-29 | Baird Corporation | Computer image display system and processor therefor |
US4315318A (en) * | 1978-12-26 | 1982-02-09 | Fuji Photo Film Co., Ltd. | Method and apparatus for processing a radiation image |
US4317179A (en) * | 1978-12-26 | 1982-02-23 | Fuji Photo Film Co., Ltd. | Method and apparatus for processing a radiographic image |
JPS5691735A (en) * | 1979-12-25 | 1981-07-24 | Fuji Photo Film Co Ltd | Method and apparatus for treating xxray image |
SE426281B (sv) * | 1981-05-05 | 1982-12-20 | Ericsson Telefon Ab L M | Forfarande for visning pa en bildskerm av utsignaler fran en stralningskenslig detektor |
US4394688A (en) * | 1981-08-25 | 1983-07-19 | Hamamatsu Systems, Inc. | Video system having an adjustable digital gamma correction for contrast enhancement |
US4517607A (en) * | 1981-11-09 | 1985-05-14 | Ricoh Company, Ltd. | Method of and apparatus for compensating image in image reproduction system |
US4438495A (en) * | 1981-11-13 | 1984-03-20 | General Electric Company | Tomography window-level gamma functions |
US4456926A (en) * | 1982-06-21 | 1984-06-26 | Thomson-Csf Broadcast, Inc. | Enhancement of fluroscopically generated images |
JPS5928144A (ja) * | 1982-08-09 | 1984-02-14 | Fuji Photo Film Co Ltd | 放射線画像再生装置 |
US4559557A (en) * | 1984-06-01 | 1985-12-17 | General Electric Company | Region-of-interest digital subtraction angiography |
-
1983
- 1983-11-13 IL IL70214A patent/IL70214A/xx unknown
-
1984
- 1984-11-10 DE DE3441162A patent/DE3441162C2/de not_active Expired - Lifetime
- 1984-11-12 FR FR848417206A patent/FR2555003B1/fr not_active Expired - Lifetime
- 1984-11-13 JP JP59239290A patent/JP2575612B2/ja not_active Expired - Fee Related
- 1984-11-13 NL NL8403468A patent/NL193182C/nl not_active IP Right Cessation
-
1987
- 1987-09-16 US US07/098,742 patent/US4907288A/en not_active Expired - Lifetime
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5323236A (en) * | 1976-08-16 | 1978-03-03 | Fuji Electric Co Ltd | Optical discrimination equipment |
JPS5834677A (ja) * | 1981-08-25 | 1983-03-01 | Matsushita Graphic Commun Syst Inc | 画情報読取装置のmtf補正方式 |
Also Published As
Publication number | Publication date |
---|---|
US4907288A (en) | 1990-03-06 |
FR2555003B1 (fr) | 1991-11-29 |
JP2575612B2 (ja) | 1997-01-29 |
FR2555003A1 (fr) | 1985-05-17 |
DE3441162C2 (de) | 1994-01-13 |
IL70214A (en) | 1987-10-20 |
NL8403468A (nl) | 1985-06-03 |
DE3441162A1 (de) | 1985-07-25 |
IL70214A0 (en) | 1984-02-29 |
NL193182B (nl) | 1998-10-01 |
NL193182C (nl) | 1999-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPS60169980A (ja) | 画像のコントラスト劣化を修正する方法及びその装置 | |
Seibert et al. | Removal of image intensifier veiling glare by mathematical deconvolution techniques | |
JPH0233190B2 (ja) | ||
US4903205A (en) | Method and apparatus for displaying radiation image, and method and apparatus for calculating unsharp mask signal used for the same | |
EP0905650B1 (en) | Image enhancement processing apparatus and method | |
JP2849964B2 (ja) | 画像処理方法および装置 | |
JP3540914B2 (ja) | X線撮影装置 | |
Rehm et al. | Design and testing of artifact-suppressed adaptive histogram equalization: a contrast-enhancement technique for display of digital chest radiographs | |
JPH01259842A (ja) | 放射線診断装置 | |
JPS58201169A (ja) | 画像デ−タ処理装置 | |
JP4584550B2 (ja) | X線計測装置 | |
JPH0326597B2 (ja) | ||
JP3945976B2 (ja) | 画像表示方法及び装置 | |
US5604780A (en) | Method for processing a radiation image | |
Landini et al. | Notes on the implementation of the mass–radius method of fractal dimension estimation | |
US6781603B2 (en) | Image display method and apparatus | |
JPH0654984B2 (ja) | X線診断装置 | |
JPH1084504A (ja) | X線画像表示装置 | |
Angelhed et al. | Optimization and evaluation of an image intensifier TV system for digital chest imaging | |
Rajapakshe et al. | Quality control for on-line portal imaging systems | |
JP2696997B2 (ja) | デジタル画像データ変換装置 | |
JP2791558B2 (ja) | データ圧縮方法 | |
Katsulai et al. | A qualitative evaluation of image fidelity in image processing | |
Webber et al. | Digital subtraction of radiographs as an aid for detecting pathologic changes in periodontal bone | |
JP2000298076A (ja) | 鮮鋭度評価方法およびこれを用いる鮮鋭度評価装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |