JPH04271677A - X線記録のダイナミック圧縮方法及びその装置 - Google Patents
X線記録のダイナミック圧縮方法及びその装置Info
- Publication number
- JPH04271677A JPH04271677A JP3272850A JP27285091A JPH04271677A JP H04271677 A JPH04271677 A JP H04271677A JP 3272850 A JP3272850 A JP 3272850A JP 27285091 A JP27285091 A JP 27285091A JP H04271677 A JPH04271677 A JP H04271677A
- Authority
- JP
- Japan
- Prior art keywords
- value
- image
- low
- image value
- pass
- 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 abstract description 29
- 230000006835 compression Effects 0.000 title claims description 21
- 238000007906 compression Methods 0.000 title claims description 21
- 230000009466 transformation Effects 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 12
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 2
- 238000003745 diagnosis Methods 0.000 abstract description 21
- 238000006243 chemical reaction Methods 0.000 abstract description 8
- 238000001914 filtration Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 17
- 210000004072 lung Anatomy 0.000 description 12
- 238000012545 processing Methods 0.000 description 12
- 230000008859 change Effects 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 210000000056 organ Anatomy 0.000 description 5
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 229910052711 selenium Inorganic materials 0.000 description 4
- 239000011669 selenium Substances 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/409—Edge or detail enhancement; Noise or error suppression
- H04N1/4092—Edge or detail enhancement
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/407—Control or modification of tonal gradation or of extreme levels, e.g. background level
- H04N1/4072—Control or modification of tonal gradation or of extreme levels, e.g. background level dependent on the contents of the original
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
め要約のデータは記録されません。
Description
が画素に割当てられるX線記録のダイナミック圧縮方法
及びその方法を実施する装置に係る。
密度差を有する。例えば、肺の写真において、脊柱及び
心臓はX線に対し肺組織よりかなり大きい密度を有する
。即ち、それらはX線をかなり大きい程度に吸収する。 対象におけるこれらの密度差はX線画像検出器により1
0,000の係数までの強度差に変換さるべきX線三次
元画像になる。従来のX線フィルムは局部コントラスト
を失うことなくそのようなダイナミックレンジに対処で
きない。しかし、人は肺写真に対して所定の限られたダ
イナミックレンジでのみ関心を持つので、例えば比較的
弱く吸収する肺組織において、フィルムはこのダイナミ
ックサブレンジがフィルムダイナミックを充分に用いる
よう普通露出され、現像される。他のダイナミックレン
ジは失なわれた細部コントラストがあるよう露出角又は
露出不足でフィルムに示される。現在の画像検出器はよ
り大きいダイナミックレンジを有し、これによりそれら
は従来のフィルムで記録されえないダイナミック境界範
囲内に画像として用いられうるデータを供給することが
できる。これらの検出器は、例えば光導電層を含み、望
ましくはセレンを含み、X線三次元画像を電荷パターン
に変換し、画素を電荷センサーにより走査され、ディジ
タル入力画像値に順次に変換しうる。潜像としてX線三
次元画像を蓄積する記憶発光体で被われた板は、レーザ
ビームで照射されると、同等のダイナミックレンジを有
するX線放射の強度に比例した発光を生ずる。光センサ
を用いると、発光は電気信号に変換されえ、次にディジ
タル化されうる。
レンジを含む入力画像値がその画像を発生する表示ユニ
ットに更なる処理なしに印加される場合、「従来の」X
線記録の文字に対応する画像文字が用いられるべき時に
、入力画像値に含まれる大部分の情報は失なわれる。 全ての画像領域にけおるコントラストを減ずることによ
り診断に適切でない領域での情報のこの損失を避けるこ
とは事実可能であるが、診断に重要である密度範囲でア
クティブであるこのコントラスト減少は画像文字を従来
の記録からはずれさせ;多数の医者はかかる低コントラ
スト(「かすかな」)画像を拒絶にする。
放射が制御可能な減衰素子により検査の領域を貫通する
形に変化され、検査領域後で画像検出器前の強度差は非
常に大きく平滑化される方法を開示する。その場合、従
来のX線フィルムを用いて、一方で従来のX線画像の文
字(平均密度範囲での)を保ち、他方で暗い画像部分の
細部コントラストの表示を可能にするX線画像を得るこ
とも又可能である。
果が放射強度の位置依存変化なしに得られ、X線記録の
入力画像値が記録された対象の密度差により決定された
ダイナミックレンジを有するように前記のタイプの方法
を実施することを目的とする。
的は、低域通過画像値は入力画像値から形成され、各画
素に対して、出力画像値は関連入力画像値と関連低通過
画像値に依存する等化値の重畳により形成され、出力画
像値は小さい低通過画像値で入力画像値より大きく、大
きい低通過画像値で入力画像値より小さいことにより達
成される。
造の画像情報が除去される、即ち高い空間周波数を有す
る成分は抑圧される入力画像値とは異なる。現在の低域
通過画像値に依存する等化値の重畳により、低い空間周
波数を有する大きい領域(粗)構造のダイナミックが大
きく制限されることが達成される。出力画像値が表示ユ
ニット(例えばレーザ画像器)に印加される場合、入力
画像値を処理すると、従来の記録文字が得られ、等化値
の適切な変化で、従来の画像の文字が診断に関連する密
度範囲で得られ、一方他の密度範囲で、細部構造が減じ
られた大領域コントラストで得られるように、その転送
特性が調整される。視聴者のブラインドをする入力画像
の明るい領域は作られた出力画像中実質的な平均密度を
有し、細部コントラストがより容易に検出可能である領
域に転送される。本発明は、前記の画像検出器からの入
力画像の使用に制限されず、他の方法で作られたX線記
録、例えば画像増培管−テレビジョンチェーンで作られ
た画像で使用するのに適している。
通過画像値がX線記録のディジタル入力画像値から形成
される方法を既に示していることに注目さるべきである
。これらの低域通過画像値は、高域通過画像が差として
得られるよう入力画像値から減算される。この高域画像
は、任意に1より大きい係数で逓倍された後、元の画像
に加算される「アンシャープマスキング」方法といわれ
るこの効果は、縁部が減衰されることであり;ダイナミ
ック減少を目的とせずまたこれを達成しない。本発明の
更なる展開は、等化値の低域通過画像値への依存度を示
す等化関数はルックアップテーブルに蓄積され、各低域
通過画像値に対して、等化値がルックアップテーブルか
ら読出されることを提供とする。これは、所定の低域通
過画像値に対して、関連等化値をX線記録の個々の画素
に非常に素早く割り当てることを可能にする。
は望ましくは相互作用的に利用者により設定されえ、等
化関数はプリセットパラメータに応じて計算され、ルッ
クアップテーブルに蓄積されることを提供する。これは
モニターのパラメータの変化の効果をチェックすること
を可能とする。利用者は、最適画像印象を得るまでパラ
メータを相互作用的に変化させうる。従って、データを
集め或いはいわゆる投写放射線写真で現在まで慣例でな
いモニターの表示スクリーンに表示されるX線画像を基
に診断をすることが可能である。
るセレン又は記憶発生層は、それで作られた画像信号か
種々の画素でX線放射の強度に直線的に依存する性質を
有する。このタイプの検出器に適している本発明の更な
る展開は、画像信号の入力画像値が対数変換により得ら
れることを提供する。強度に直線的に依存する画像信号
から得られ、対数変換により得られた入力画像値は関連
画素の対象の密度に略直線的に依存する。本発明による
方法は、強度に直線的に依存する入力画像値に対する密
度に依存する入力画像値に対し、良い結果を与える。
X線記録の位置依存放射強度を電気信号に変換する画像
検出器と、b)信号から入力画像値を得るための手段と
、c)記憶容量がX線記録を蓄積するのに必要とされる
容量より大きい記憶装置と、d)入力画像値から低域通
過映像値を発生する手段と、e)低域通過画像値から等
化画像値を生じるための手段と、f)等化画像値を出力
画像値の再生用入力画像値と重畳する手段と、g)出力
画像値に応じてX線記録を表示する表示ユニットとから
なることを特徴とする。
。
X線ビーム2を出射する。患者の後ろに作られたX線三
次元画像は画像検出器4,5により時間による変化が対
象の後ろのX線放射の空間変化に対応する電気信号に変
換される。
在であり、光導電体、望ましくはセレンで被覆された円
筒形ドラム4を含む。X線ビーム2は放射扇形を有さな
ければならず、その面に回転軸40が位置する。X線記
録中、ドラムは軸40の回りに回転し、同時に扇形放射
ビーム2は、患者3に対して動かされ、これによりX線
記録中患者の2次元領域がドラム面に表示される。回転
軸40に並行に変位自在であるよう配置される適切なセ
ンサ5により電気信号に変換されるX線放射は、ドラム
4に充電パターンを発生する。この画像検出器の更なる
詳細は欧州特許PS第219897号を引用する。
うる。画像検出器は、時間による変化が対象の後のX線
放射の強度の空間的変化に依存する電気信号(又は複数
の電気信号)を発生することが必要なだけである。
れ、その後ディジタルアナログ変換器7で例えば16ビ
ットワイドデータワードのシーケンスに変換される。こ
れらのデータは、記憶容量がX線記録の記憶に必要とさ
れるより大きい半導体記憶器81と協働する画像処理コ
ンピュータ80を含む画像処理器8に印加される。画像
処理器8は一連の処理ステップ、なかんずくアナログデ
ィジタル変換器により供給される信号のダイナミック圧
縮を実行する。その結果、暗又は明画像部分での細部コ
ントラストを損ずることなく、表示ユニット9にX線記
録を表示することが可能である。表示ユニットはモニタ
ー又は画像処理器の出力信号からフィルム画像を発生す
るレーザ画像器でよい。
テップの時間経過を示す。ディジタルアナログ変換の後
、その目的の為、検出器特定誤差を除去しなければなら
ない画像補正方法(ブロック81)が先ず実行される。 かかる誤差は、センサ5によるセレンドラムの面の電荷
の走査中、自己放電が生じるという事実から起こる。セ
ンサがドラム面を同時に異なる場所で走査する複数の素
子からなる場合、更なる誤差は感度の差及び感度変動で
生じる。かかる画像誤差用補正方法はそれ自体公知であ
り、更なる説明は必要ない。他のタイプの画像検出器が
用いられる場合、他の補正方法が用いられなければなら
ない。この補正方法により理想的に検出誤差がないデー
タセットSが得られる。
に関してより詳細に説明されるダイナミック圧縮(ブロ
ック82)を受ける。ダイナミックが上記方法で圧縮さ
れた画像データAは次に更なる処理ステップ(ブロック
83)例えばエッジ強調方法か表示ユニット9への適合
を受けうる。
に示す如く、先ず対数変換(ブロック821)を含み、
それによりX線の強度に直線的に依存する補正画像デー
タSは異なる画素で対象の密度に少なくとも実質的に直
線的に依存するデータEに変換れる。画像検出器が密度
に直線的に関連するデータ自体を供給する場合、この変
換は省きうる。
E=Emax * log (Sma
x (S)/log (Smax )) (1)
ここでSmax は画像検出器で供給される最も高い可
能値、例えば30,000である。Emax は密度画
像での最も高い可能値に対応し、非常に高い値、コンピ
ュータで表わされうる整数、例えば32,000に設定
される。式1によって、値E=0はS=Smax に対
して得られ、値E=Emax が最も小さい値S=1に
対して得られる。変換パラメータのこの選択で、ダイナ
ミック圧縮は結果的に起こらない。効果的に、対数変換
はアドレスとして画像データSを解釈するルックアップ
テーブルの助けでなされ、その状態下で式1による値E
が記憶される。対数変換により値Sから発生した値は以
後入力画像値という。これらの映像値は矢印で示す如く
外部記憶例えばディスク記憶に印加されうる。
画像データが発生される前、即ち、アナログ領域でなさ
れうる。対数増幅器はこの目的の為用いられうる;ブロ
ック821は省かれうる。
ク822)である。ここで、入力画像値Eの和は、X線
記録の望ましくは方形サブ領域、零空間に位置する画素
にその和が割り当てられる公知の方法でなされる。作ら
れた和は零空間の画像値の数で分割され、零空間の中央
に低域通過画像値Lとして画素が割当てられる。この値
は結果的に零空間における画像値及び画素の算術平均を
示す。この方法は他の画素に対して繰り返され、線方向
に隣接する画素用低域通過画像値は、一方で1例の画素
が零空間に加算され、他方で一列の画素が除かれる公知
の方法で得られる。零空間の寸法は、医療診断に関する
細部構造よりかなり大きく、例えば3×3cmの領域に
対応しなければならない。
L+H
(2)が成り立つよう関連する画素及びその環境で
のX線記録の強度分布に依存する値Hで形成される低域
通過画像値Lと異なる。値Hは入力画像値Eのより高い
空間周波数成分を示し、又、全ての画素の値Hから、細
部構造だけを含む高域通過画像が夫々形成される。
のすぐ後に続き、それにより単に低域通過画像のダイナ
ミックは少なくともサブ領域に圧縮される。各画素に対
して、低域通過画像値Lに依存する等化値Cが決められ
、入力画像値Eに重畳され、等化値は、出力画像におい
て、診断に対して重要な領域で大きい領域構造のダイナ
ミックが最小可能範囲で影響され、一方診断に対して第
2の関係である別な密度領域で圧縮されるように選択さ
れる。等化機能、即ち、等化値Cの低域通過画像値Lへ
の依存関係、C値に対するスケールはL値に対するスケ
ールの2倍であることを示す詳細の例示を図4に示す。
度範囲は値L1及びL2の間に位置するとする。これら
の値は、肺の複数のX線写真のヒントグラム分析により
例えば、L1=5,00及びL2=10,000に決定
される。L1及びL2間の密度範囲において、圧縮は出
来る限り小さくなければならない。これは、等化値C=
0が診断に対して第1の関係の領域に位置することで達
成される。
することにより形成される場合)、より低い限界値L1
以下で正でなければならない。等化値は、限界値L1及
び現在の値Lの間の差がより高くなるに従いより高くな
ければならない。しかし、差L1−Lは、密度比が反転
されるのを防ぐようLと関連した値Cより小さい必要は
ない。実線により図1に示され、点L1で終る曲線は値
L上の値Cの適切な依存性を示す。それは下式で表わさ
れる: C=x1* L
1* (L1−L/L1)y1 (
3)ここで、x1は1より大きくない係数であり、一方
y1は1より小さくない指数である。曲線は値x1=0
.5及びy1=1.2に基づく。
Cはこの密度範囲でダイナミック圧縮を行うよう負でな
ければならない。この範囲でも、曲線は全く単調な(降
下)形を有さなければならず、LとL2の間の差はLの
現在の値に割当てられたCの値より小さい必要はない。 L=L2から始まる実線は下式に定義される:
C=x2* (L2 −Emax )* ((L−L2
)/(Emax −L2))y2
(4)ここで、x2は1以上でなく、y2は1
より小さくない。図4に示す曲線は値x2=0.5及び
y2=1.2に基づく。
ット値L1及びL2は診断に対して第1の関係の密度領
域を限定する。しかし、自動露出タイマーにより発生さ
れた同じ露出時間を用いる肺写真に対して、この仮定が
有効なだけである。
に適合される場合、異なる露出時間が用いられる肺写真
を用いることも可能である。単に肺X線写真用線量の変
化は、全体密度値が常に同じ量だけ変化するよう作用す
るので、この量のみが適合をするためにL1及びL2の
プリセット値に重畳されるべきで、これによりL1とL
2の間の差はこの適合の後でさえ変化しない;しかし、
L1及びL2を不変のままにし、その代わり低域通過画
像値について反転記号を有する値を重畳することも可能
である。患者の異なる露出時間又は異なる厚さに適合す
るために、夫々加算又は減算されなければならない量は
関連した記録の画素で異なる密度値の周波数を示すヒス
トグラムを分析することにより容易に得られうる。実際
、そのようなヒストグラム分析はそれらが他の理由、例
えば画像値を画像表示ユニットの伝達特性に適合させる
為に常に必要とするルールとして、追加コスト及び設計
努力を意味しない。
=1が選択される場合、変化Cが3つの直線圧分で形成
される破線により示される曲線が得られる。曲線のL1
及びL2での2つの曲りはやや欠点である。一点鎖線に
より示される如く、Cの適切な変化は診断に重要である
領域の(L0での)中央を通る負の傾斜を有する直線に
より達成されうる。この直線は、式(3)及び(4)に
おいて、x1=x2=0.5;y1=y2=1及びL1
=L2=L0が選択される場合に得られる。これは下式
になる: C=x* (L0−L)
(5)ここで例えばx=0.5。これはL1と
しての間の診断に関して密度範囲での大きい領域コント
ラストのかかる影響に帰するがこの影響は比較的小さい
。
A=E+C
(6)によりそれに重畳された後、出力画像値
から形成される出力画像になる、即ち肺写真に最適であ
る。診断に関連した肺組織は従来の写真で慣例の大きい
領域で、細部コントラストで表示される。他方で、大き
い領域で荒いコントラストはより明るく、より暗い画像
部分で減少され、一方出力画像の細部コントラスト(H
)はこれらの領域で保たれる。出力画像値Aが成分Hを
含む式(2)により細部構造を示す入力映像値Eからの
合計で形成される事実からこれは説明されうる。
依存性は、全てのX線記録に対して等しく適してはいな
い。骨の画像に対して、例えば高密度領域は特別の診断
有意性があり、一方より低密度を有する領域は診断に対
する重要性は低い。この場合、ダイナミック圧縮はより
低い密度範囲で作用されなければならず、一方より高い
密度範囲は変わらないままにあるべきである。このため
に、Cは(図4の左でプロットされた曲線部分に似てい
るが、より高い密度値で始まる)、より低い密度値に向
け増加しなければならず、一方より高い密度の範囲のC
は、基本的に変化せず、従来の映像文字が保たれること
を提供するため、ゼロに実質的に等しくなければならな
い。
る程度の依存関係を異なる器官の記録に対し与えること
が必要である。しかし、全ての場合において、CはLの
機能として、単調な変化を有し、特に式(6)によりE
及びCが共に加算される場合、変化を単調に減少させ、
その代わりCがEから減算される場合変化を単調に増加
する。
域通過画像値と関連した等化値C、即ち図3のブロック
823が決定される。そうするための最も単純な方法は
、低域通過画像値Lで形成される各入力アドレスに対し
て、それに割当てられる等化値−画像値Cが蓄積される
ルックアップテーブルによる。異なる器官のX線記録に
対して、多分異なる露出時間を有する写真に対して、又
は異なる厚さの同じ器官の患者の写真に対しても、値C
のLへの依存関係の異なる程度が利用されなければなら
ないので、最近撮影された器官で生じるCの変化は先ず
、、できれば先行する計算(ブロック824)の後記録
からロードされなければならない。
力画像値及び式(6)による関連した等化値Cは各画素
に対して重畳される。しかし、値Cが軸C=0に関して
鏡像反転される方法で変化する場合、値CをEから減算
することが基本的に可能となる。出力画像値Aで形成さ
れた出力画像では、大きな領域コントラストは、診断に
対して単に2次的関心である領域で少なくとも減少され
る。細部コントラストがこれらの領域で変わらないまま
であるので、それらは出力画像ではさらに識別できる。
置換えられる場合、F及びuを保持する式:F=(1−
u)* E+u* L
(7)はプリセット可能な値であり、uが1より大き
い場合細部コントラストは診断に対して基本的関係の密
度範囲の両側の領域で改善される。u=1に対して、F
はLと等しくなり、細部コントラストは変わらないまま
である。 Eは実質的にLに対応するので(高域通過用値H及び細
部構造は夫々普通低域通過画像値Lとわずかに比較され
る)、Fは実質的にLに対応し、u=1に対してそれと
同一である。
と、低域通過画像値Lと、出力画像値Aが利用出来なけ
ればならないが、3つの(半導体)画像記憶を要する全
てのこれらの値を同時に記憶することは必要ではない。 例えば、全ての値Lを有する低域通過画像が計算され、
蓄積された場合、入力画像値E又は低域画像値LはCが
決められ、Eに加算された後、出力画像値Aにより重ね
書きされえ;これは従って2つの画像記憶だけを必要と
する。低域通過画像に対し必要でない値E又は低域通過
画像の前に計算された低域通過画像値が完全に計算され
た場合、1つの画像を蓄積するのに必要とされるより多
くない記憶容量で十分である。
ック821)を相殺する指数変換(826)に追加的に
従ってよい;しかし指数変換は省かれてもよい。
はハードコピーユニット(レーザ画像器)と呼ばれるも
のによりX線記録の供給に基づく。しかし、その代わり
にモニターの表示スクリーンにX線記録の表示のダイナ
ミック圧縮をすることも可能である。モニターはX線記
録をフィルムにロードするユニットより更に小さいダイ
ナミックレンジを有する。圧縮方法の適切な実行におい
て、即ち、所望の計算方法が十分に高速でなされる場合
、検査員は相互作用にてダイナミック圧縮を実行し、そ
れによりX線記録は最適にモニターのダイナミックレン
ジに適合されうる。利用者は結果を得るがモニターの表
示スクリーンに生じるX線画像に基づいて診断を行う。 図5は適切な表示装置、即ちモニターに加えて所望の計
算及び役立つ記憶容量を有する装置により実現されうる
概略ブロック回路図を基にしたこの相互作用方法を示す
。
21で再び対数変換を受ける。画像データSは直接検出
器1…5(図1参照)から得られうるが、或いはそれら
は外部記憶から減じられる。対数変換の後に得られた入
力画像値Eは第1の画像記憶に入れられる。
れる相互作用されたダイナミック圧縮の後、これらの値
は画像ハーフトーン変換(ブロック830)を受け、デ
ィジタルアナログ変換器831を介してモニターに印加
される。画像ハーフトーン変換を得るための最も単純な
方法は各出力画像値Aにディジタルアナログ変換の後モ
ニター900の表示スクリーンに画像ハーフトーンを生
じるディジタル値を割当てるルックアップテーブルによ
る。
ダイナミック圧縮なしの出力画像値Aとの間の関係を示
す。第1の出力画像値以下で、画素は完全に黒であり、
第2の出力画像値以上で、画素はその最大の白を有する
。これらの値の間で図6の(A)に示す画像ハーフトー
ンへの値Aの変換が生じる。
ボードでよい入力ユニット851を用いて、プリセット
曲線のユニットコントラスト及び輝度は第1の動作モー
ドで変換されうる。輝度の変化は曲線の水平ずれを生じ
、コントラストの変化はこの曲線の傾斜に影響する。 利用者は、出力画像値Aの範囲が診断に関することが黒
と白の間の転移領域に位置するよう輝度及びコントラス
ト値を選択する。この場合、値A=L0は診断に対して
基本的関係の振幅範囲に位置する。
略的に示す画像値処理動作は表示装置に慣例の処理動作
に対応する。
理ステップ及び成分を夫々用いて実行される。図3を参
照して説明した如く、入力画像値Eは、一定の比較的に
大きい零空間を有する低域通過濾波を受ける。従って作
られた低域通過画像値Lは第2の画像記憶828に蓄積
される。第2の画像記憶から読取られた低域通過画像値
Lは、各画素に対して、出力画像値A(x,y)を得る
ために、下記動作 A(x,y)=E(x,y)
=C(L(x,y)) (8)がなされるよ
う、等化関数C(L)が蓄積されるルックアップテーブ
ル823により等化値に常に割当てられる。従って、こ
の動作は各画素に対して(記憶827,828の値のE
(x,y)及びL(x,y)への)2つの記憶アクセス
、ルックアップテーブル動作及び加算動作(ブロック8
25での)を必要とする。最近の慣用計算パワーの助け
により、これらの動作は非常に素早くなされうる。ダイ
ナミック圧縮は利用者の特別な要望に相互作用的に適合
されうる。その結果、入力ユニット851は(コントラ
スト及び輝度設定モード以外に)更に動作するモードで
点線ボックスと協働する。等化関数が基本的に3つの直
線圧分からなるとすると;L1以下の低域通過画像値用
第1の直線圧分は角度β1で定義され、その下に、この
直線圧分が水平に延びる。角度β1は45°より小さく
なければならず、即ち数値で示される傾斜m1=tan
(β1)は1より小さくなければならない。L1からL
2への第2の直線区分は診断に対して最も関係のある入
力値に割当てられる。この範囲の位置は、一方でX線記
録の露出時間に依存し、他方でX線記録で表示された器
官又は人体範囲に依存する。輝度が相互作用的に設定れ
るので、診断に対する第1義の振幅範囲内に位置する値
L0が既に存在する。値a及びbのみが依然プリセット
されなければならず、この範囲の最後の値はそれにより
下記の如く定義されうる:L1 =L0−a
(9)L2 =L
0+b (1
0)L2 以上の低域通過画像値用第3の直線圧分は夫
々角度β2又は傾斜m2=tamβ2で定義され、ここ
で数値で示されるm2は1より小さくなければならない
。
けをセットしなければならない。これらのパラメータは
、例えばスクリーンに図形表示された等化関数を基にし
て表示スクリーンにマウスを「クリック」することによ
り互いに独立に調整されうる。しかし、その代わりa=
b及びm1=m2を設定することにより4つのパラメー
タを2つに削減することが可能である。残る2つのパラ
メータはマウスによる双次元変位、例えばパラメータa
又はbを係るx−方向での変位によりプリセットされ、
一方y方向の変位はパラメータm1及びm2を変える。 これは(異なる動作モードで)x方向の変位が輝度を制
御し、y方向の変位がコントラストを制御することによ
る方法に全く似ている。
でなされる。第1の直線圧分に対して、下式、C=m1
* (L−L1)・L1 (1
1)は有効であり、L1の計算は式(9)によりL0及
びaからなされる。
して、C=0であり、第3の直線圧分に対して下記が成
り立つ: C=m2* (L−L2)ここでL2=L0+b。
わば実時間でなされるので、利用者が相互作用的にパラ
メータを変化させる時、等化関数は非常に素早く計算さ
れうる。
れてもよい更なる「角度」である。従って、等化関数C
(L)を示すデータワードシーケンスを(一次元)低域
通過濾波するようにすることが推奨される。そのため、
低域通過フィルタブロック852が設けられ、低域通過
零空間の寸法は診断に関連する振幅範囲の幅(a+b)
から適切に導かれ、例えばこの値の1/10に設定され
る。例えばL1=5000及びL2=10,000であ
る場合、等化関数の平滑化は、等化関数C(L)の50
0連続値の平均がとれらるようになされる。計算され平
滑された関数はルックアップテーブルにロードされ、そ
れから新たなプリセットパラメータにより圧縮がなされ
る。必要により、これらのパラメータは再び変換等され
うる。
ある。
。
す複数の曲線を示す。
図である。
縮で使用されるモニターの特性曲線と、相互作用的に調
整された等化関数を示す。
理器 9 表示ユニット 40 回転軸 80 画像処理コンピュータ 81 半導体記憶 83 処理ステップ 821 対数変換 822 低域通過濾波動作 823 ルックアップテーブル 825 加算動作 826 指数変換 827,828 記憶 830 画像ハーフトーン変換 851 入力ユニット 900 モニター
Claims (9)
- 【請求項1】 低域通過画像値(L)は入力画像値(
E)から形成され、各画素に対して、出力画像値(A)
は関連入力画像値(E)と関連低通過画像値(L)に依
存する等化値(C)の重畳により形成され、出力画像値
は小さい低通過画像値(L)で入力画像値より大きく、
大きい低通過画像値(L)で入力画像値(E)より小さ
いことを特徴とする、その画素にディジタル入力画像値
が割当てられるX線記録におけるダイナミック圧縮方法
。 - 【請求項2】 等化値(C)の低域通過画像値(L)
への依存度を示す等化関数(C,L)はルックアップテ
ーブル(823)に蓄積され、各低域通過画像値(L)
に対して、等化値(C)がルックアップテーブルから読
出されることを特徴とする請求項1記載の方法。 - 【請求項3】 等化関数(C(L))のパラメータ(
m1,m2,L1,L2)は望ましくは相互作用的に利
用者により設定されえ、等化関数はプリセットパラメー
タに応じて計算(852)され、ルックアップテーブル
(823)に蓄積されることを特徴とする請求項2記載
の方法。 - 【請求項4】 等化機能のパラメータのプリセットに
対して、低域通過画像値の範囲(L1−L2)は、等化
関数(C(1))が値0を有する範囲内で設定され、等
化関数の上向き傾斜(m1,m2)はこの範囲外に設定
されることを特徴とする請求項1記載の方法。 - 【請求項5】 等化関数は平滑化手順を受けることを
特徴とする請求項1乃至4のうちいずれか一項記載の方
法。 - 【請求項6】 入力画像値(E)は対数変換(821
)により画像信号(S)から得られることを特徴とする
、X線放射の強度に直線的に依存する画像信号を発生す
る画像検出器を有する請求項1記載の方法。 - 【請求項7】 出力画像値(A)は指数変換を受ける
ことを特徴とする請求項1乃至6のうちいずれか一項記
載の方法。 - 【請求項8】 a)X線記録の位置依存放射強度を電
気信号に変換する画像検出器(4,5)と、b)信号か
ら入力画像値(E)を得るための手段(7,8;821
)と、 c)記憶容量がX線記録を蓄積するのに必要とされる容
量より大きい記憶装置(81)と、 d)入力画像値から低域通過画像値を発生する手段(8
;822)と、 e)低域通過画像値から等化画像値を得るための手段(
8;823,824)と、 f)等化画像値を出力画像値(A)の再生用入力画像値
と重畳する手段(8,825)と、 g)出力画像値(A)に応じてX線記録を表示する表示
ユニット(9)とからなることを特徴とする請求項1記
載の方法を実施する装置。 - 【請求項9】 出力画像値を得る手段は、出力画像値
(L)に割当てられたアドレス下で、等化画像値(C)
が蓄積されるルックアップテーブル(823)を含むこ
とを特徴とする請求項8記載の装置。
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE19904033716 DE4033716A1 (de) | 1990-10-24 | 1990-10-24 | Verfahren zur dynamikkompression in roentgenaufnahmen und vorrichtung zur durchfuehrung des verfahrens |
| DE19914122264 DE4122264A1 (de) | 1991-07-05 | 1991-07-05 | Verfahren zur dynamikkompression in roentgenaufnahmen und vorrichtung zur durchfuehrung des verfahrens |
| DE41222644 | 1991-07-05 | ||
| DE40337162 | 1991-07-05 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04271677A true JPH04271677A (ja) | 1992-09-28 |
| JP3188491B2 JP3188491B2 (ja) | 2001-07-16 |
Family
ID=25897931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP27285091A Expired - Fee Related JP3188491B2 (ja) | 1990-10-24 | 1991-10-21 | X線記録のダイナミック圧縮方法及びその装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5357549A (ja) |
| EP (1) | EP0482712B1 (ja) |
| JP (1) | JP3188491B2 (ja) |
| DE (1) | DE59108446D1 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000209501A (ja) * | 1998-11-24 | 2000-07-28 | Ge Medical Syst Sa | 器官の厚さを補償する方法 |
| US6163029A (en) * | 1997-09-22 | 2000-12-19 | Kabushiki Kaisha Toshiba | Radiation detector, radiation detecting method and X-ray diagnosing apparatus with same radiation detector |
| JP2001111891A (ja) * | 1999-10-05 | 2001-04-20 | Toshiba Corp | 画像処理装置およびこれを用いたx線診断装置 |
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルx線イメージング検出器のダイナミックレンジ拡張システム |
| JP2002314881A (ja) * | 2001-04-13 | 2002-10-25 | Toshiba Corp | X線診断装置 |
| JP2010069319A (ja) * | 2009-12-17 | 2010-04-02 | Toshiba Corp | X線診断装置 |
Families Citing this family (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3467285B2 (ja) * | 1993-04-02 | 2003-11-17 | コニカミノルタホールディングス株式会社 | 放射線画像処理方法 |
| DE4415990A1 (de) | 1994-05-06 | 1995-11-23 | Philips Patentverwaltung | Verfahren zur Wiedergabe insbesondere einer digitalen Röntgenaufnahme als sichtbares Bild sowie Anordnung zur Durchführung des Verfahrens |
| JP3568279B2 (ja) * | 1995-06-30 | 2004-09-22 | 富士写真フイルム株式会社 | 画像再生方法および装置 |
| EP0813336B1 (en) * | 1996-06-12 | 2007-08-08 | FUJIFILM Corporation | Image processing method and apparatus |
| US5835618A (en) * | 1996-09-27 | 1998-11-10 | Siemens Corporate Research, Inc. | Uniform and non-uniform dynamic range remapping for optimum image display |
| US6782137B1 (en) * | 1999-11-24 | 2004-08-24 | General Electric Company | Digital image display improvement system and method |
| JP3893827B2 (ja) * | 1999-11-26 | 2007-03-14 | コニカミノルタホールディングス株式会社 | X線画像撮影システム |
| FR2803070B1 (fr) * | 1999-12-28 | 2002-06-07 | Ge Medical Syst Sa | Procede et systeme de gestion de la dynamique d'une image radiologique numerisee |
| US6633684B1 (en) | 2000-07-07 | 2003-10-14 | Athentech Technologies Corp. | Distortion-free image contrast enhancement |
| US6956977B2 (en) * | 2001-08-08 | 2005-10-18 | General Electric Company | Methods for improving contrast based dynamic range management |
| DE60137514D1 (de) * | 2001-10-31 | 2009-03-12 | Sony Deutschland Gmbh | Verfahren zum Verbessern der Qualität eines Bildes |
| FR2847698B1 (fr) * | 2002-11-27 | 2005-05-06 | Ge Med Sys Global Tech Co Llc | Procede de gestion de la dynamique d'une image radiologique numerique |
| US7218763B2 (en) * | 2003-02-27 | 2007-05-15 | Eastman Kodak Company | Method for automated window-level settings for magnetic resonance images |
| US20050111756A1 (en) * | 2003-11-25 | 2005-05-26 | Turner Robert W. | System and method for generating coherent data sets of images from various sources |
| DE102005003227A1 (de) * | 2005-01-24 | 2006-08-03 | Siemens Ag | Verfahren zur Reduzierung von Cupping-Artefakten in ConeBeam-CT-Bilddatensätzen |
| US8150110B2 (en) * | 2006-11-22 | 2012-04-03 | Carestream Health, Inc. | ROI-based rendering for diagnostic image consistency |
| US20080187194A1 (en) * | 2007-02-05 | 2008-08-07 | Zhang Daoxian H | Cad image normalization |
| US7844097B2 (en) * | 2007-12-03 | 2010-11-30 | Samplify Systems, Inc. | Compression and decompression of computed tomography data |
| US7916830B2 (en) * | 2008-09-11 | 2011-03-29 | Samplify Systems, Inc. | Edge detection for computed tomography projection data compression |
| US7852977B2 (en) * | 2008-09-11 | 2010-12-14 | Samplify Systems, Inc. | Adaptive compression of computed tomography projection data |
| US8588485B2 (en) | 2008-11-25 | 2013-11-19 | Carestream Health, Inc. | Rendering for improved diagnostic image consistency |
| US8045811B2 (en) * | 2008-11-26 | 2011-10-25 | Samplify Systems, Inc. | Compression and storage of projection data in a computed tomography system |
| US8290295B2 (en) * | 2009-03-03 | 2012-10-16 | Microsoft Corporation | Multi-modal tone-mapping of images |
| US8433120B2 (en) * | 2010-02-25 | 2013-04-30 | Carestream Health, Inc. | Method for image processing of mammographic images |
| US8457418B2 (en) * | 2011-08-02 | 2013-06-04 | Raytheon Company | Local area contrast enhancement |
| JP2013059511A (ja) * | 2011-09-14 | 2013-04-04 | Fujifilm Corp | 画像処理装置 |
| US9324137B2 (en) | 2012-10-24 | 2016-04-26 | Marvell World Trade Ltd. | Low-frequency compression of high dynamic range images |
| US11257261B2 (en) | 2015-07-23 | 2022-02-22 | Koninklijke Philips N.V. | Computed tomography visualization adjustment |
| WO2017148502A1 (en) * | 2016-02-29 | 2017-09-08 | Brainlab Ag | Automatic detection of an artifact in patient image data |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4317179A (en) * | 1978-12-26 | 1982-02-23 | Fuji Photo Film Co., Ltd. | Method and apparatus for processing a radiographic image |
| JPS5588740A (en) * | 1978-12-26 | 1980-07-04 | Fuji Photo Film Co Ltd | Method of treating gradation of radiation picture of breast and its device |
| NL8400845A (nl) * | 1984-03-16 | 1985-10-16 | Optische Ind De Oude Delft Nv | Inrichting voor spleetradiografie. |
| US4839805A (en) * | 1983-11-17 | 1989-06-13 | General Electric Company | Dual control of image level and window parameters of a display and the like |
| FR2580827B1 (fr) * | 1985-04-19 | 1987-05-22 | Thomson Cgr | Installation de radiologie |
| FR2581775B1 (fr) * | 1985-05-10 | 1990-03-23 | Thomson Cgr | Procede de reglage de l'echelle des gris dans la reproduction d'images numerisees et systeme mettant en oeuvre ce procede |
| DE3534768A1 (de) * | 1985-09-30 | 1987-04-02 | Philips Patentverwaltung | Anordnung zum erzeugen von roentgenaufnahmen mittels eines fotoleiters |
| EP0252327B1 (en) * | 1986-06-13 | 1991-02-20 | Fuji Photo Film Co., Ltd. | Radiation image read-out method and apparatus |
| DE3629409C2 (de) * | 1986-08-29 | 1994-02-17 | Agfa Gevaert Ag | Verfahren und Vorrichtung zur Kontrastbewertung bei der elektronischen Bildverarbeitung |
| US5042077A (en) * | 1987-10-02 | 1991-08-20 | General Electric Company | Method of highlighting subtle contrast in graphical images |
| US4918534A (en) * | 1988-04-22 | 1990-04-17 | The University Of Chicago | Optical image processing method and system to perform unsharp masking on images detected by an I.I./TV system |
| JP2663189B2 (ja) * | 1990-01-29 | 1997-10-15 | 富士写真フイルム株式会社 | 画像のダイナミックレンジ圧縮処理方法 |
-
1991
- 1991-10-21 JP JP27285091A patent/JP3188491B2/ja not_active Expired - Fee Related
- 1991-10-22 EP EP91202726A patent/EP0482712B1/de not_active Expired - Lifetime
- 1991-10-22 DE DE59108446T patent/DE59108446D1/de not_active Expired - Fee Related
- 1991-10-22 US US07/781,534 patent/US5357549A/en not_active Expired - Lifetime
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6163029A (en) * | 1997-09-22 | 2000-12-19 | Kabushiki Kaisha Toshiba | Radiation detector, radiation detecting method and X-ray diagnosing apparatus with same radiation detector |
| JP2000209501A (ja) * | 1998-11-24 | 2000-07-28 | Ge Medical Syst Sa | 器官の厚さを補償する方法 |
| JP2001111891A (ja) * | 1999-10-05 | 2001-04-20 | Toshiba Corp | 画像処理装置およびこれを用いたx線診断装置 |
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルx線イメージング検出器のダイナミックレンジ拡張システム |
| JP2002314881A (ja) * | 2001-04-13 | 2002-10-25 | Toshiba Corp | X線診断装置 |
| JP2010069319A (ja) * | 2009-12-17 | 2010-04-02 | Toshiba Corp | X線診断装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5357549A (en) | 1994-10-18 |
| EP0482712A2 (de) | 1992-04-29 |
| DE59108446D1 (de) | 1997-02-13 |
| EP0482712A3 (en) | 1992-09-16 |
| EP0482712B1 (de) | 1997-01-02 |
| JP3188491B2 (ja) | 2001-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3188491B2 (ja) | X線記録のダイナミック圧縮方法及びその装置 | |
| US5224177A (en) | High quality film image correction and duplication method and system | |
| JP3402630B2 (ja) | コントラストを強化する方法と装置 | |
| JP2551552B2 (ja) | 原画像の細部を強調した画像を表示もしくは記録する方法 | |
| Dobbins 3rd et al. | Threshold perception performance with computed and screen-film radiography: implications for chest radiography. | |
| US5550888A (en) | Method of displaying notably a digital X-ray image as a visible image, and device for carrying out the method | |
| CN116018611B (zh) | 使用深度卷积网络的噪声抑制 | |
| US6718056B1 (en) | Method of automatic determination of the contrast and brightness of a digital radiographic image | |
| US5714764A (en) | Method for detecting prospective abnormal patterns | |
| CA2349224A1 (en) | Method for automatic detection of region of interest for digital x-ray detectors using a filtered histogram | |
| JP2003283927A (ja) | 画像のコントラストを強調する方法 | |
| Kimme-Smith et al. | Effects of reduced exposure on computed radiography: comparison of nodule detection accuracy with conventional and asymmetric screen-film radiographs of a chest phantom. | |
| JP2952519B2 (ja) | 放射線画像の階調変換装置 | |
| JPS6253178B2 (ja) | ||
| Uffmann et al. | Computed radiography and direct radiography: influence of acquisition dose on the detection of simulated lung lesions | |
| Ishida | Image processing | |
| Cowen et al. | An image processing algorithm for PPCR imaging | |
| JPS5910839A (ja) | サブトラクシヨン画像の処理方法 | |
| Dobbins 3rd et al. | Variable compensation chest radiography performed with a computed radiography system: design considerations and initial clinical experience. | |
| JPS61169971A (ja) | X線画像処理装置 | |
| EP0527524B1 (en) | Image quality enhancing method | |
| Yarwood et al. | Image processing of digital chest ionograms | |
| Van Metter et al. | Exposure latitude for thoracic radiography | |
| Angelhed et al. | Optimization and evaluation of an image intensifier TV system for digital chest imaging | |
| JPS6262375B2 (ja) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080511 Year of fee payment: 7 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080511 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090511 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090511 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100511 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110511 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |