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
Application number
JP3272850A
Other languages
English (en)
Other versions
JP3188491B2 (ja
Inventor
Hanns-Ingo Maack
ハンス−インゴ マック
Ulrich Neitzel
ウルリッヒ ニーゼル
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
Priority claimed from DE19904033716 external-priority patent/DE4033716A1/de
Priority claimed from DE19914122264 external-priority patent/DE4122264A1/de
Application filed by Philips Gloeilampenfabrieken NV filed Critical Philips Gloeilampenfabrieken NV
Publication of JPH04271677A publication Critical patent/JPH04271677A/ja
Application granted granted Critical
Publication of JP3188491B2 publication Critical patent/JP3188491B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/409Edge or detail enhancement; Noise or error suppression
    • H04N1/4092Edge or detail enhancement
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/407Control or modification of tonal gradation or of extreme levels, e.g. background level
    • H04N1/4072Control 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

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

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、ディジタル入力画像値
が画素に割当てられるX線記録のダイナミック圧縮方法
及びその方法を実施する装置に係る。
【0002】X線記録に示される対象領域は屡々大きい
密度差を有する。例えば、肺の写真において、脊柱及び
心臓はX線に対し肺組織よりかなり大きい密度を有する
。即ち、それらはX線をかなり大きい程度に吸収する。 対象におけるこれらの密度差はX線画像検出器により1
0,000の係数までの強度差に変換さるべきX線三次
元画像になる。従来のX線フィルムは局部コントラスト
を失うことなくそのようなダイナミックレンジに対処で
きない。しかし、人は肺写真に対して所定の限られたダ
イナミックレンジでのみ関心を持つので、例えば比較的
弱く吸収する肺組織において、フィルムはこのダイナミ
ックサブレンジがフィルムダイナミックを充分に用いる
よう普通露出され、現像される。他のダイナミックレン
ジは失なわれた細部コントラストがあるよう露出角又は
露出不足でフィルムに示される。現在の画像検出器はよ
り大きいダイナミックレンジを有し、これによりそれら
は従来のフィルムで記録されえないダイナミック境界範
囲内に画像として用いられうるデータを供給することが
できる。これらの検出器は、例えば光導電層を含み、望
ましくはセレンを含み、X線三次元画像を電荷パターン
に変換し、画素を電荷センサーにより走査され、ディジ
タル入力画像値に順次に変換しうる。潜像としてX線三
次元画像を蓄積する記憶発光体で被われた板は、レーザ
ビームで照射されると、同等のダイナミックレンジを有
するX線放射の強度に比例した発光を生ずる。光センサ
を用いると、発光は電気信号に変換されえ、次にディジ
タル化されうる。
【0003】一方、かく得られ、対象の全ダイナミック
レンジを含む入力画像値がその画像を発生する表示ユニ
ットに更なる処理なしに印加される場合、「従来の」X
線記録の文字に対応する画像文字が用いられるべき時に
、入力画像値に含まれる大部分の情報は失なわれる。 全ての画像領域にけおるコントラストを減ずることによ
り診断に適切でない領域での情報のこの損失を避けるこ
とは事実可能であるが、診断に重要である密度範囲でア
クティブであるこのコントラスト減少は画像文字を従来
の記録からはずれさせ;多数の医者はかかる低コントラ
スト(「かすかな」)画像を拒絶にする。
【0004】
【従来の技術】欧州特許出願第158382号は、X線
放射が制御可能な減衰素子により検査の領域を貫通する
形に変化され、検査領域後で画像検出器前の強度差は非
常に大きく平滑化される方法を開示する。その場合、従
来のX線フィルムを用いて、一方で従来のX線画像の文
字(平均密度範囲での)を保ち、他方で暗い画像部分の
細部コントラストの表示を可能にするX線画像を得るこ
とも又可能である。
【0005】
【発明が解決しようとする課題】本発明は比較可能な結
果が放射強度の位置依存変化なしに得られ、X線記録の
入力画像値が記録された対象の密度差により決定された
ダイナミックレンジを有するように前記のタイプの方法
を実施することを目的とする。
【0006】
【課題を解決するための手段】本発明によれば、この目
的は、低域通過画像値は入力画像値から形成され、各画
素に対して、出力画像値は関連入力画像値と関連低通過
画像値に依存する等化値の重畳により形成され、出力画
像値は小さい低通過画像値で入力画像値より大きく、大
きい低通過画像値で入力画像値より小さいことにより達
成される。
【0007】低域通過映像値は、X線記録を含む細部構
造の画像情報が除去される、即ち高い空間周波数を有す
る成分は抑圧される入力画像値とは異なる。現在の低域
通過画像値に依存する等化値の重畳により、低い空間周
波数を有する大きい領域(粗)構造のダイナミックが大
きく制限されることが達成される。出力画像値が表示ユ
ニット(例えばレーザ画像器)に印加される場合、入力
画像値を処理すると、従来の記録文字が得られ、等化値
の適切な変化で、従来の画像の文字が診断に関連する密
度範囲で得られ、一方他の密度範囲で、細部構造が減じ
られた大領域コントラストで得られるように、その転送
特性が調整される。視聴者のブラインドをする入力画像
の明るい領域は作られた出力画像中実質的な平均密度を
有し、細部コントラストがより容易に検出可能である領
域に転送される。本発明は、前記の画像検出器からの入
力画像の使用に制限されず、他の方法で作られたX線記
録、例えば画像増培管−テレビジョンチェーンで作られ
た画像で使用するのに適している。
【0008】西独特許明細書第2952422号は低域
通過画像値がX線記録のディジタル入力画像値から形成
される方法を既に示していることに注目さるべきである
。これらの低域通過画像値は、高域通過画像が差として
得られるよう入力画像値から減算される。この高域画像
は、任意に1より大きい係数で逓倍された後、元の画像
に加算される「アンシャープマスキング」方法といわれ
るこの効果は、縁部が減衰されることであり;ダイナミ
ック減少を目的とせずまたこれを達成しない。本発明の
更なる展開は、等化値の低域通過画像値への依存度を示
す等化関数はルックアップテーブルに蓄積され、各低域
通過画像値に対して、等化値がルックアップテーブルか
ら読出されることを提供とする。これは、所定の低域通
過画像値に対して、関連等化値をX線記録の個々の画素
に非常に素早く割り当てることを可能にする。
【0009】更なる実施例では、等化関数のパラメータ
は望ましくは相互作用的に利用者により設定されえ、等
化関数はプリセットパラメータに応じて計算され、ルッ
クアップテーブルに蓄積されることを提供する。これは
モニターのパラメータの変化の効果をチェックすること
を可能とする。利用者は、最適画像印象を得るまでパラ
メータを相互作用的に変化させうる。従って、データを
集め或いはいわゆる投写放射線写真で現在まで慣例でな
いモニターの表示スクリーンに表示されるX線画像を基
に診断をすることが可能である。
【0010】X線記録をディジタル電気信号に変換でき
るセレン又は記憶発生層は、それで作られた画像信号か
種々の画素でX線放射の強度に直線的に依存する性質を
有する。このタイプの検出器に適している本発明の更な
る展開は、画像信号の入力画像値が対数変換により得ら
れることを提供する。強度に直線的に依存する画像信号
から得られ、対数変換により得られた入力画像値は関連
画素の対象の密度に略直線的に依存する。本発明による
方法は、強度に直線的に依存する入力画像値に対する密
度に依存する入力画像値に対し、良い結果を与える。
【0011】上記の方法を実施するための装置は、a)
X線記録の位置依存放射強度を電気信号に変換する画像
検出器と、b)信号から入力画像値を得るための手段と
、c)記憶容量がX線記録を蓄積するのに必要とされる
容量より大きい記憶装置と、d)入力画像値から低域通
過映像値を発生する手段と、e)低域通過画像値から等
化画像値を生じるための手段と、f)等化画像値を出力
画像値の再生用入力画像値と重畳する手段と、g)出力
画像値に応じてX線記録を表示する表示ユニットとから
なることを特徴とする。
【0012】
【実施例】本発明を以下図面を参照して詳細に説明する
【0013】図1において、X線放射器は患者3を通る
X線ビーム2を出射する。患者の後ろに作られたX線三
次元画像は画像検出器4,5により時間による変化が対
象の後ろのX線放射の空間変化に対応する電気信号に変
換される。
【0014】画像検出器は例えば軸40の回りを回転自
在であり、光導電体、望ましくはセレンで被覆された円
筒形ドラム4を含む。X線ビーム2は放射扇形を有さな
ければならず、その面に回転軸40が位置する。X線記
録中、ドラムは軸40の回りに回転し、同時に扇形放射
ビーム2は、患者3に対して動かされ、これによりX線
記録中患者の2次元領域がドラム面に表示される。回転
軸40に並行に変位自在であるよう配置される適切なセ
ンサ5により電気信号に変換されるX線放射は、ドラム
4に充電パターンを発生する。この画像検出器の更なる
詳細は欧州特許PS第219897号を引用する。
【0015】他のタイプの画像検出器が代りに用いられ
うる。画像検出器は、時間による変化が対象の後のX線
放射の強度の空間的変化に依存する電気信号(又は複数
の電気信号)を発生することが必要なだけである。
【0016】センサ5の出力信号は、増幅器6に印加さ
れ、その後ディジタルアナログ変換器7で例えば16ビ
ットワイドデータワードのシーケンスに変換される。こ
れらのデータは、記憶容量がX線記録の記憶に必要とさ
れるより大きい半導体記憶器81と協働する画像処理コ
ンピュータ80を含む画像処理器8に印加される。画像
処理器8は一連の処理ステップ、なかんずくアナログデ
ィジタル変換器により供給される信号のダイナミック圧
縮を実行する。その結果、暗又は明画像部分での細部コ
ントラストを損ずることなく、表示ユニット9にX線記
録を表示することが可能である。表示ユニットはモニタ
ー又は画像処理器の出力信号からフィルム画像を発生す
るレーザ画像器でよい。
【0017】図2は画像処理器8での異なる画像処理ス
テップの時間経過を示す。ディジタルアナログ変換の後
、その目的の為、検出器特定誤差を除去しなければなら
ない画像補正方法(ブロック81)が先ず実行される。 かかる誤差は、センサ5によるセレンドラムの面の電荷
の走査中、自己放電が生じるという事実から起こる。セ
ンサがドラム面を同時に異なる場所で走査する複数の素
子からなる場合、更なる誤差は感度の差及び感度変動で
生じる。かかる画像誤差用補正方法はそれ自体公知であ
り、更なる説明は必要ない。他のタイプの画像検出器が
用いられる場合、他の補正方法が用いられなければなら
ない。この補正方法により理想的に検出誤差がないデー
タセットSが得られる。
【0018】すぐその後で、得られる画像データは図3
に関してより詳細に説明されるダイナミック圧縮(ブロ
ック82)を受ける。ダイナミックが上記方法で圧縮さ
れた画像データAは次に更なる処理ステップ(ブロック
83)例えばエッジ強調方法か表示ユニット9への適合
を受けうる。
【0019】ブロック82で示される処理動作は、図3
に示す如く、先ず対数変換(ブロック821)を含み、
それによりX線の強度に直線的に依存する補正画像デー
タSは異なる画素で対象の密度に少なくとも実質的に直
線的に依存するデータEに変換れる。画像検出器が密度
に直線的に関連するデータ自体を供給する場合、この変
換は省きうる。
【0020】対数変換は下式によりなされる:    
        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から発生した値は以
後入力画像値という。これらの映像値は矢印で示す如く
外部記憶例えばディスク記憶に印加されうる。
【0021】根本的に、対数変換は、或いはディジタル
画像データが発生される前、即ち、アナログ領域でなさ
れうる。対数増幅器はこの目的の為用いられうる;ブロ
ック821は省かれうる。
【0022】次のステップは低域通過濾波動作(ブロッ
ク822)である。ここで、入力画像値Eの和は、X線
記録の望ましくは方形サブ領域、零空間に位置する画素
にその和が割り当てられる公知の方法でなされる。作ら
れた和は零空間の画像値の数で分割され、零空間の中央
に低域通過画像値Lとして画素が割当てられる。この値
は結果的に零空間における画像値及び画素の算術平均を
示す。この方法は他の画素に対して繰り返され、線方向
に隣接する画素用低域通過画像値は、一方で1例の画素
が零空間に加算され、他方で一列の画素が除かれる公知
の方法で得られる。零空間の寸法は、医療診断に関する
細部構造よりかなり大きく、例えば3×3cmの領域に
対応しなければならない。
【0023】各画素に対して、入力画像値Eは式:E=
L+H                      
  (2)が成り立つよう関連する画素及びその環境で
のX線記録の強度分布に依存する値Hで形成される低域
通過画像値Lと異なる。値Hは入力画像値Eのより高い
空間周波数成分を示し、又、全ての画素の値Hから、細
部構造だけを含む高域通過画像が夫々形成される。
【0024】これは動作(ブロック823,…825)
のすぐ後に続き、それにより単に低域通過画像のダイナ
ミックは少なくともサブ領域に圧縮される。各画素に対
して、低域通過画像値Lに依存する等化値Cが決められ
、入力画像値Eに重畳され、等化値は、出力画像におい
て、診断に対して重要な領域で大きい領域構造のダイナ
ミックが最小可能範囲で影響され、一方診断に対して第
2の関係である別な密度領域で圧縮されるように選択さ
れる。等化機能、即ち、等化値Cの低域通過画像値Lへ
の依存関係、C値に対するスケールはL値に対するスケ
ールの2倍であることを示す詳細の例示を図4に示す。
【0025】例えば、肺写真に対して、診断に関する密
度範囲は値L1及びL2の間に位置するとする。これら
の値は、肺の複数のX線写真のヒントグラム分析により
例えば、L1=5,00及びL2=10,000に決定
される。L1及びL2間の密度範囲において、圧縮は出
来る限り小さくなければならない。これは、等化値C=
0が診断に対して第1の関係の領域に位置することで達
成される。
【0026】値Cは(出力画像値がE及びCを共に加算
することにより形成される場合)、より低い限界値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に基づく。
【0027】L2を越える低域通過画像値Lに対して、
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に基づく。
【0028】前記において、各肺写真に対して、プリセ
ット値L1及びL2は診断に対して第1の関係の密度領
域を限定する。しかし、自動露出タイマーにより発生さ
れた同じ露出時間を用いる肺写真に対して、この仮定が
有効なだけである。
【0029】しかし、値L1及びL2が現在の露出時間
に適合される場合、異なる露出時間が用いられる肺写真
を用いることも可能である。単に肺X線写真用線量の変
化は、全体密度値が常に同じ量だけ変化するよう作用す
るので、この量のみが適合をするためにL1及びL2の
プリセット値に重畳されるべきで、これによりL1とL
2の間の差はこの適合の後でさえ変化しない;しかし、
L1及びL2を不変のままにし、その代わり低域通過画
像値について反転記号を有する値を重畳することも可能
である。患者の異なる露出時間又は異なる厚さに適合す
るために、夫々加算又は減算されなければならない量は
関連した記録の画素で異なる密度値の周波数を示すヒス
トグラムを分析することにより容易に得られうる。実際
、そのようなヒストグラム分析はそれらが他の理由、例
えば画像値を画像表示ユニットの伝達特性に適合させる
為に常に必要とするルールとして、追加コスト及び設計
努力を意味しない。
【0030】式(3)及び(4)において、y1=y2
=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と
しての間の診断に関して密度範囲での大きい領域コント
ラストのかかる影響に帰するがこの影響は比較的小さい
【0031】等化値の上記依存性は、入力画像値が式:
A=E+C                    
    (6)によりそれに重畳された後、出力画像値
から形成される出力画像になる、即ち肺写真に最適であ
る。診断に関連した肺組織は従来の写真で慣例の大きい
領域で、細部コントラストで表示される。他方で、大き
い領域で荒いコントラストはより明るく、より暗い画像
部分で減少され、一方出力画像の細部コントラスト(H
)はこれらの領域で保たれる。出力画像値Aが成分Hを
含む式(2)により細部構造を示す入力映像値Eからの
合計で形成される事実からこれは説明されうる。
【0032】図4に示す等化値Cの低域通過画像値Lの
依存性は、全てのX線記録に対して等しく適してはいな
い。骨の画像に対して、例えば高密度領域は特別の診断
有意性があり、一方より低密度を有する領域は診断に対
する重要性は低い。この場合、ダイナミック圧縮はより
低い密度範囲で作用されなければならず、一方より高い
密度範囲は変わらないままにあるべきである。このため
に、Cは(図4の左でプロットされた曲線部分に似てい
るが、より高い密度値で始まる)、より低い密度値に向
け増加しなければならず、一方より高い密度の範囲のC
は、基本的に変化せず、従来の映像文字が保たれること
を提供するため、ゼロに実質的に等しくなければならな
い。
【0033】従って、等化値の低域通過画像値への異な
る程度の依存関係を異なる器官の記録に対し与えること
が必要である。しかし、全ての場合において、CはLの
機能として、単調な変化を有し、特に式(6)によりE
及びCが共に加算される場合、変化を単調に減少させ、
その代わりCがEから減算される場合変化を単調に増加
する。
【0034】ダイナミック圧縮を行なうよう、現在の低
域通過画像値と関連した等化値C、即ち図3のブロック
823が決定される。そうするための最も単純な方法は
、低域通過画像値Lで形成される各入力アドレスに対し
て、それに割当てられる等化値−画像値Cが蓄積される
ルックアップテーブルによる。異なる器官のX線記録に
対して、多分異なる露出時間を有する写真に対して、又
は異なる厚さの同じ器官の患者の写真に対しても、値C
のLへの依存関係の異なる程度が利用されなければなら
ないので、最近撮影された器官で生じるCの変化は先ず
、、できれば先行する計算(ブロック824)の後記録
からロードされなければならない。
【0035】続いて、出力画像値を形成するために、入
力画像値及び式(6)による関連した等化値Cは各画素
に対して重畳される。しかし、値Cが軸C=0に関して
鏡像反転される方法で変化する場合、値CをEから減算
することが基本的に可能となる。出力画像値Aで形成さ
れた出力画像では、大きな領域コントラストは、診断に
対して単に2次的関心である領域で少なくとも減少され
る。細部コントラストがこれらの領域で変わらないまま
であるので、それらは出力画像ではさらに識別できる。
【0036】式(3)から(5)において、値LがFで
置換えられる場合、F及びuを保持する式:F=(1−
u)* E+u* L               
 (7)はプリセット可能な値であり、uが1より大き
い場合細部コントラストは診断に対して基本的関係の密
度範囲の両側の領域で改善される。u=1に対して、F
はLと等しくなり、細部コントラストは変わらないまま
である。 Eは実質的にLに対応するので(高域通過用値H及び細
部構造は夫々普通低域通過画像値Lとわずかに比較され
る)、Fは実質的にLに対応し、u=1に対してそれと
同一である。
【0037】X線記録の各画素に対して、入力画像値E
と、低域通過画像値Lと、出力画像値Aが利用出来なけ
ればならないが、3つの(半導体)画像記憶を要する全
てのこれらの値を同時に記憶することは必要ではない。 例えば、全ての値Lを有する低域通過画像が計算され、
蓄積された場合、入力画像値E又は低域画像値LはCが
決められ、Eに加算された後、出力画像値Aにより重ね
書きされえ;これは従って2つの画像記憶だけを必要と
する。低域通過画像に対し必要でない値E又は低域通過
画像の前に計算された低域通過画像値が完全に計算され
た場合、1つの画像を蓄積するのに必要とされるより多
くない記憶容量で十分である。
【0038】必要により、出力画像値は対数変換(ブロ
ック821)を相殺する指数変換(826)に追加的に
従ってよい;しかし指数変換は省かれてもよい。
【0039】本発明によるダイナミック圧縮の前記説明
はハードコピーユニット(レーザ画像器)と呼ばれるも
のによりX線記録の供給に基づく。しかし、その代わり
にモニターの表示スクリーンにX線記録の表示のダイナ
ミック圧縮をすることも可能である。モニターはX線記
録をフィルムにロードするユニットより更に小さいダイ
ナミックレンジを有する。圧縮方法の適切な実行におい
て、即ち、所望の計算方法が十分に高速でなされる場合
、検査員は相互作用にてダイナミック圧縮を実行し、そ
れによりX線記録は最適にモニターのダイナミックレン
ジに適合されうる。利用者は結果を得るがモニターの表
示スクリーンに生じるX線画像に基づいて診断を行う。 図5は適切な表示装置、即ちモニターに加えて所望の計
算及び役立つ記憶容量を有する装置により実現されうる
概略ブロック回路図を基にしたこの相互作用方法を示す
【0040】図5により補正画像データSはブロック8
21で再び対数変換を受ける。画像データSは直接検出
器1…5(図1参照)から得られうるが、或いはそれら
は外部記憶から減じられる。対数変換の後に得られた入
力画像値Eは第1の画像記憶に入れられる。
【0041】更に示されるべき、出力画像値Aが発生さ
れる相互作用されたダイナミック圧縮の後、これらの値
は画像ハーフトーン変換(ブロック830)を受け、デ
ィジタルアナログ変換器831を介してモニターに印加
される。画像ハーフトーン変換を得るための最も単純な
方法は各出力画像値Aにディジタルアナログ変換の後モ
ニター900の表示スクリーンに画像ハーフトーンを生
じるディジタル値を割当てるルックアップテーブルによ
る。
【0042】図6の(A)はモニターの画素の輝度Hと
ダイナミック圧縮なしの出力画像値Aとの間の関係を示
す。第1の出力画像値以下で、画素は完全に黒であり、
第2の出力画像値以上で、画素はその最大の白を有する
。これらの値の間で図6の(A)に示す画像ハーフトー
ンへの値Aの変換が生じる。
【0043】トラックボール、いわゆるマウス又はキー
ボードでよい入力ユニット851を用いて、プリセット
曲線のユニットコントラスト及び輝度は第1の動作モー
ドで変換されうる。輝度の変化は曲線の水平ずれを生じ
、コントラストの変化はこの曲線の傾斜に影響する。 利用者は、出力画像値Aの範囲が診断に関することが黒
と白の間の転移領域に位置するよう輝度及びコントラス
ト値を選択する。この場合、値A=L0は診断に対して
基本的関係の振幅範囲に位置する。
【0044】ここに記載される限りにおいて、図5に概
略的に示す画像値処理動作は表示装置に慣例の処理動作
に対応する。
【0045】ダイナミック圧縮は破線ボックスで示す処
理ステップ及び成分を夫々用いて実行される。図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より小さくなければならない
【0046】利用者は従って値a,b,m1及びm2だ
けをセットしなければならない。これらのパラメータは
、例えばスクリーンに図形表示された等化関数を基にし
て表示スクリーンにマウスを「クリック」することによ
り互いに独立に調整されうる。しかし、その代わりa=
b及びm1=m2を設定することにより4つのパラメー
タを2つに削減することが可能である。残る2つのパラ
メータはマウスによる双次元変位、例えばパラメータa
又はbを係るx−方向での変位によりプリセットされ、
一方y方向の変位はパラメータm1及びm2を変える。 これは(異なる動作モードで)x方向の変位が輝度を制
御し、y方向の変位がコントラストを制御することによ
る方法に全く似ている。
【0047】等化関数を罫線することはブロック852
でなされる。第1の直線圧分に対して、下式、C=m1
* (L−L1)・L1            (1
1)は有効であり、L1の計算は式(9)によりL0及
びaからなされる。
【0048】L1からL2の範囲の第2の直線圧分に対
して、C=0であり、第3の直線圧分に対して下記が成
り立つ: C=m2* (L−L2)ここでL2=L0+b。
【0049】直線区分だけが含まれ、即ち等化関数がい
わば実時間でなされるので、利用者が相互作用的にパラ
メータを変化させる時、等化関数は非常に素早く計算さ
れうる。
【0050】計算の後、等化関数は多くの場合に妨害さ
れてもよい更なる「角度」である。従って、等化関数C
(L)を示すデータワードシーケンスを(一次元)低域
通過濾波するようにすることが推奨される。そのため、
低域通過フィルタブロック852が設けられ、低域通過
零空間の寸法は診断に関連する振幅範囲の幅(a+b)
から適切に導かれ、例えばこの値の1/10に設定され
る。例えばL1=5000及びL2=10,000であ
る場合、等化関数の平滑化は、等化関数C(L)の50
0連続値の平均がとれらるようになされる。計算され平
滑された関数はルックアップテーブルにロードされ、そ
れから新たなプリセットパラメータにより圧縮がなされ
る。必要により、これらのパラメータは再び変換等され
うる。
【図面の簡単な説明】
【図1】本発明を実施する装置のブロック回路系統図で
ある。
【図2】画像処理中になされる動作の時間的順序を示す
【図3】ダイナミック圧縮のフローチャートを示す。
【図4】出力画像値対低域通過画像値の適切な変化を示
す複数の曲線を示す。
【図5】相互作用ダイナミック圧縮のブロック回路系統
図である。
【図6】(A),(B)は夫々相互作用ダイナミック圧
縮で使用されるモニターの特性曲線と、相互作用的に調
整された等化関数を示す。
【符号の説明】
2  X線ビーム 3  患者 4,5  画像検出器 6  増幅器 7,831  ディジタルアナログ変換器8  画像処
理器 9  表示ユニット 40  回転軸 80  画像処理コンピュータ 81  半導体記憶 83  処理ステップ 821  対数変換 822  低域通過濾波動作 823  ルックアップテーブル 825  加算動作 826  指数変換 827,828  記憶 830  画像ハーフトーン変換 851  入力ユニット 900  モニター

Claims (9)

    【特許請求の範囲】
  1. 【請求項1】  低域通過画像値(L)は入力画像値(
    E)から形成され、各画素に対して、出力画像値(A)
    は関連入力画像値(E)と関連低通過画像値(L)に依
    存する等化値(C)の重畳により形成され、出力画像値
    は小さい低通過画像値(L)で入力画像値より大きく、
    大きい低通過画像値(L)で入力画像値(E)より小さ
    いことを特徴とする、その画素にディジタル入力画像値
    が割当てられるX線記録におけるダイナミック圧縮方法
  2. 【請求項2】  等化値(C)の低域通過画像値(L)
    への依存度を示す等化関数(C,L)はルックアップテ
    ーブル(823)に蓄積され、各低域通過画像値(L)
    に対して、等化値(C)がルックアップテーブルから読
    出されることを特徴とする請求項1記載の方法。
  3. 【請求項3】  等化関数(C(L))のパラメータ(
    m1,m2,L1,L2)は望ましくは相互作用的に利
    用者により設定されえ、等化関数はプリセットパラメー
    タに応じて計算(852)され、ルックアップテーブル
    (823)に蓄積されることを特徴とする請求項2記載
    の方法。
  4. 【請求項4】  等化機能のパラメータのプリセットに
    対して、低域通過画像値の範囲(L1−L2)は、等化
    関数(C(1))が値0を有する範囲内で設定され、等
    化関数の上向き傾斜(m1,m2)はこの範囲外に設定
    されることを特徴とする請求項1記載の方法。
  5. 【請求項5】  等化関数は平滑化手順を受けることを
    特徴とする請求項1乃至4のうちいずれか一項記載の方
    法。
  6. 【請求項6】  入力画像値(E)は対数変換(821
    )により画像信号(S)から得られることを特徴とする
    、X線放射の強度に直線的に依存する画像信号を発生す
    る画像検出器を有する請求項1記載の方法。
  7. 【請求項7】  出力画像値(A)は指数変換を受ける
    ことを特徴とする請求項1乃至6のうちいずれか一項記
    載の方法。
  8. 【請求項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. 【請求項9】  出力画像値を得る手段は、出力画像値
    (L)に割当てられたアドレス下で、等化画像値(C)
    が蓄積されるルックアップテーブル(823)を含むこ
    とを特徴とする請求項8記載の装置。
JP27285091A 1990-10-24 1991-10-21 X線記録のダイナミック圧縮方法及びその装置 Expired - Fee Related JP3188491B2 (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 富士写真フイルム株式会社 画像のダイナミックレンジ圧縮処理方法

Cited By (6)

* Cited by examiner, † Cited by third party
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