JPH06292013A - 放射線画像のダイナミックレンジ圧縮方法 - Google Patents
放射線画像のダイナミックレンジ圧縮方法Info
- Publication number
- JPH06292013A JPH06292013A JP5075900A JP7590093A JPH06292013A JP H06292013 A JPH06292013 A JP H06292013A JP 5075900 A JP5075900 A JP 5075900A JP 7590093 A JP7590093 A JP 7590093A JP H06292013 A JPH06292013 A JP H06292013A
- Authority
- JP
- Japan
- Prior art keywords
- sus
- image
- signal
- value
- dynamic range
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 40
- 238000007906 compression Methods 0.000 title claims abstract description 39
- 230000006835 compression Effects 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims description 34
- 238000012937 correction Methods 0.000 claims abstract description 99
- 230000007423 decrease Effects 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000003247 decreasing effect Effects 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 49
- 238000004364 calculation method Methods 0.000 claims description 3
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 abstract description 18
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 abstract description 18
- 210000004072 lung Anatomy 0.000 description 18
- 238000010586 diagram Methods 0.000 description 16
- 230000008859 change Effects 0.000 description 13
- 238000006243 chemical reaction Methods 0.000 description 7
- 210000000038 chest Anatomy 0.000 description 6
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000001370 mediastinum Anatomy 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000007740 vapor deposition Methods 0.000 description 1
Landscapes
- Facsimile Image Signal Circuits (AREA)
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
いて、素抜け部や照射野絞りの部分が必要以上に補正さ
れることを回避する。 【構成】非鮮鋭マスク信号Susの関数f1(Sus)を補
正値とし、該補正値f1(Sus)をオリジナル画像信号
をSorg に加算することにより、ダイナミックレンジの
圧縮を行う。ここで、前記関数f1(Sus)の関数形
が、2つの基準値Sus1,Sus2で挟まれる領域ではマ
スク信号Susの増加に伴って単調減少し、基準値Sus1
を下回る領域ではマスク信号Susの減少に伴って単調減
少し、基準値Sus2を上回る領域ではマスク信号Susの
増加に伴って単調増加する。これにより、高濃度及び低
濃度側での補正度合いが抑制され、素抜け部や照射野絞
りと関心領域との信号差が保持される。
Description
クレンジ圧縮方法に関し、詳しくは、オリジナル画像信
号を処理して該オリジナル画像よりもダイナミックレン
ジの狭い画像を担持する画像信号を得る放射線画像のダ
イナミックレンジ圧縮方法の改善技術に関する。
細構造の観察適正を確保しつつ、濃度域を圧縮する方法
として、特開平3−222577号公報に開示されるよ
うなものがある。前記特開平3−222577号公報に
開示される圧縮方法は、各画素点に対応して該各画素点
を含む所定マスク領域内のオリジナル画像信号Sorg を
平均化することで非鮮鋭マスク信号(ボケマスク信号)
Susを求め、この非鮮鋭マスク信号が増大するにつれて
単調減少する関数をf(Sus)としたときに、処理済み
画像信号Sprocを、 Sproc=Sorg +f(Sus) として得るものである。
端に信号値が異なる場合、ダイナミックレンジの圧縮を
行わずに、前記信号値の違いをそのまま保持させたい場
合がある。しかしながら、従来のダイナミックレンジ圧
縮方法によると、非鮮鋭マスク信号が小さいときほど大
きな増大補正を施し、また、非鮮鋭マスク信号が大きい
ときほど大きな減少補正を施す構成であったから、例え
ば高濃度領域の圧縮を行う場合、信号値を補正したくな
い素抜け部分の信号値を必要以上に下げてしまったり、
低濃度領域の圧縮においては、照射野絞りを行って極端
に信号値が低い部分を、無用に上げてしまうという問題
があった。
極端に信号値が異なる部分があると、その部分に圧縮処
理を行ってしまうという問題がある。例えば、胸部正面
画像で低濃度部を圧縮する場合に、肺野部(高濃度部)
にペースメーカなどの金属があると、このペースメーカ
が存在する部分は信号値が小さくなるため、補正値を信
号値の関数で設定すると望まない圧縮が行われ、ペース
メーカの信号レベルを肺野の信号レベルに近づけてしま
い、信号値の違いを保持できないという問題が生じる。
あり、素抜け部や照射野絞りの部分に対して必要以上の
圧縮補正がなされることを回避でき、また、圧縮を行い
たくない領域に例えば信号値の極端に大きい又は小さい
部分があっても、かかる極端な信号差を有する部分に影
響されて圧縮処理がなされてしまうことを回避できるダ
イナミックレンジの圧縮方法を提供することを目的とす
る。
放射線画像のダイナミックレンジ圧縮方法は、被写体を
透過した放射線画像情報に基づくオリジナル画像を表す
オリジナル画像信号Sorg を処理して前記オリジナル画
像よりもダイナミックレンジの狭い画像を担持する処理
済み画像信号Sprocを得るもので、各画素点に対応して
非鮮鋭マスク信号Susを求め、該非鮮鋭マスク信号Sus
の関数f1(Sus)を補正値としたときに、 Sproc=Sorg +f1(Sus) なる演算を行って前記処理済み画像信号Sprocを得る放
射線画像のダイナミックレンジ圧縮方法において、前記
補正値f1(Sus)が、前記非鮮鋭マスク信号Susの所
定中間値SusAで正の最大値をとり、前記所定中間値S
usAを上回る領域では非鮮鋭マスク信号Susの増大に応
じて単調減少し、前記所定中間値SusAを下回る領域で
は前記正の最大値に固定されるか又は非鮮鋭マスク信号
Susの減少に応じて単調減少するようにした。
鮮鋭マスク信号Susの所定中間値SusBで負の最小値を
とり、前記所定中間値SusBを下回る領域では非鮮鋭マ
スク信号Susの増大に応じて単調減少し、前記所定中間
値SusBを上回る領域では前記負の最小値に固定される
か又は非鮮鋭マスク信号Susの増大に応じて単調増加す
るようにした。
せ、前記非鮮鋭マスク信号Susの2つの異なる中間値を
SusA,SusB(SusA<SusB)としたときに、前記
補正値f1(Sus)が、前記中間値SusAで正の最大値
をとり、かつ、前記中間値SusBで負の最小値をとる一
方、前記2つの中間値SusA,SusBで挟まれる信号領
域では非鮮鋭マスク信号Susの増大に応じて単調減少
し、中間値SusAを下回る領域では前記正の最大値に固
定されるか又は非鮮鋭マスク信号Susの減少に応じて単
調減少し、前記中間値SusBを上回る領域では前記負の
最小値に固定されるか又は非鮮鋭マスク信号Susの増大
に応じて単調増加するようにしても良い。
基づくオリジナル画像を表すオリジナル画像信号Sorg
を処理して前記オリジナル画像よりもダイナミックレン
ジの狭い画像を担持する処理済み画像信号Sprocを得る
放射線画像のダイナミックレンジ圧縮方法であって、補
正関数をFとしたときに、 Sproc=Sorg +F なる演算を行ってダイナミックレンジが圧縮された画像
信号Sprocを得る放射線画像のダイナミックレンジ圧縮
方法において、前記補正関数Fとして予め異なる複数種
の関数を備え、画像内の複数領域に対して異なる補正関
数Fを定めるようにした。
写体の輪郭と被写体の内部構造の輪郭と照射野絞りの境
界とのうちの少なくとも1つによって区分することが好
ましい。また、前記複数の補正関数Fが、各画素点に対
応する非鮮鋭マスク信号Susの増大に応じて単調減少す
る関数f1(Sus)と、前記非鮮鋭マスク信号Susに因
らない関数f2とを少なくとも含むようにすると良い。
域において、補正値f1(Sus)による増大補正量が増
加することがないから、信号値の低い部分での必要以上
の圧縮を抑止できる。同様に、中間値SusBを上回る領
域においては、補正値f1(Sus)による減少補正量が
増加することがないから、信号値の高い部分での必要以
上の圧縮を抑止できる。
を境とする補正値の変化特性を同時に設定することで、
信号値の高い部分と低い部分との両方で過圧縮を回避し
つつ、圧縮を行わせることが可能となる。また、同一画
像内で1つの補正関数のみによってダイナミックレンジ
を圧縮させるのではなく、領域を被写体の輪郭,被写体
の内部構造の輪郭,照射野絞りの境界によって区分し、
各領域毎に異なる関数を定めて圧縮処理を行わせること
によって、不本意な補正が行われることを回避する。
て、各画素点毎の非鮮鋭マスク信号Susを用いて補正値
を設定する構成の場合には、信号値の極端に異なる異物
が存在すると、信号値に基づく補正値の設定では、前記
異物の影響を受けることになってしまうので、前記非鮮
鋭マスク信号Susを用いない補正値を画像領域毎に使い
分けられるようにした。
を示す図1は、本発明にかかる放射線画像のダイナミッ
クレンジ圧縮方法が適用された画像処理装置を含む放射
線画像情報読取装置を示すものであり、医療用としての
人体の胸部放射線撮影を行う例を示す。
置2によって制御されて、被写体(人体胸部等)Mに向
けて放射線(一般的にはX線)を照射する。記録読取装
置3は、被写体Mを挟んで放射線源1と対向する面に変
換パネル4を備え、該変換パネル4は放射線源1からの
照射放射線量に対する被写体Mの放射線透過率分布に従
ったエネルギーを輝尽層に蓄積し、そこに被写体Mの潜
像を形成する。
を、輝尽性蛍光体の気相堆積、或いは輝尽性蛍光体塗料
塗布によって設けてあり、該輝尽層は環境による悪影響
及び損傷を遮断するために保護部材によって遮蔽若しく
は被覆される。該輝尽性蛍光体材料としては、例えば、
特開昭61−72091号公報、或いは、特開昭59−
75200号公報に開示されるような材料が使われる。
ザ,半導体レーザ等)5は、出射強度が制御された光ビ
ームを発生し、その光ビームは種々の光学系を経由して
走査器6に到達し、そこで偏向を受け、更に、反射鏡7
で光路を偏向させて、変換パネル4に輝尽励起走査光と
して導かれる。集光体8は、輝尽励起光が走査される変
換パネル4に近接して光ファイバである集光端が位置さ
れ、上記光ビームで走査された変換パネル4からの潜像
エネルギーに比例した発光強度の輝尽発光を受光する。
9は、集光体8から導入された光から輝尽発光波長領域
の光のみを通過させるフィルタであり、該フィルタ9を
通過した光は、フォトマル10に入射して、その入射光に
対応した電流信号に光電変換される。
圧変換器11で電圧信号に変換され、増幅器12で増幅され
た後、A/D変換器13で画素毎のディジタルデータから
なる放射線画像信号に変換される。そして、このディジ
タル放射線画像信号(オリジナル画像信号Sorg )は、
マイクロコンピュータを内蔵した画像処理装置14に順次
出力される。
メモリ(磁気ディスク装置)である。また、16は画像処
理装置14から直接又は前記画像メモリ15から読み出され
た放射線画像信号をプリンタ17に伝送するためのインタ
ーフェイスである。18は読取ゲイン調整回路であり、こ
の読取ゲイン調整回路18により光ビーム発生部5の光ビ
ーム強度調整、フォトマル用高圧電源19の電源電圧調整
によるフォトマル10のゲイン調整、電流/電圧変換器11
と増幅器12のゲイン調整、及びA/D変換器13の入力ダ
イナミックレンジの調整が行われ、放射線画像信号の読
取ゲインが総合的に調整される。
ジナル放射線画像信号Sorg の取得方法を、輝尽性蛍光
体を励起光で走査して発光せしめた輝尽発光を光電変化
して得る方法に限定するものではなく、例えば放射線フ
ィルムの画像を、光電変換によって読み取る方法や、被
写体を透過した放射線を蛍光体に照射されて蛍光に変換
し、該蛍光を光電変換して読み取る方法などであっても
良い。
された放射線の強度に比例する形でも良いし、検出され
た放射線の強度の対数に比例する形でも良いが、後者の
方が好ましい。ここで、前記画像処理装置14には、入力
されるオリジナル画像信号Sorg のダイナミックレンジ
を圧縮して、オリジナル画像よりもダイナミックレンジ
の狭い画像を担持する処理済みの画像信号Sprocを得る
画像処理機能が備えられており、かかるダイナミックレ
ンジ圧縮のための画像処理は、以下の式に従って行われ
る。
所定マスク領域内のオリジナル画像信号Sorg を平均化
することにより求められる非鮮鋭マスク信号である。但
し、前記非鮮鋭マスク信号Susは、マスク領域内の平均
化によって求める方法に限定されるものではなく、例え
ばメジアン値を用いる構成であっても良い。
れるf1(Sus)は、非鮮鋭マスク信号Susの関数とし
て求められる補正値であり、例えば図2に示すような関
数形に設定される。図2に示す関数形を式で示すと、以
下のようになる。
補正値変化の基準となる基準非鮮鋭マスク信号Susi
(i=0〜5)は、オリジナル画像の関心領域内の最大
信号値,最小信号値をSmax,Smin としたときに、 Susi=ki・Smax +(1−ki)・Smin なる式で求めるようにした。そして、前記係数ki(i
=0〜5)は、例えばk0=−0.8 、k1=−0.2 、k
2=0.5 、k3=1.0 、k4=1.5 、k5=1.8とし
た。ここで、前記係数k0,k1,k4,k5について
は、これらによって決定されるSus0,Sus1,Sus
4,Sus5が、関心領域の信号範囲の外側に設定される
ようにしてある。
信号範囲Sus1〜Sus4については、非鮮鋭マスク信号
Susの増加に伴って単調減少し、基準値Sus1(中間値
SusA)で正の最大値をとった補正値f1(Sus)は、
基準値Sus1を下回る信号領域では、逆に非鮮鋭マスク
信号Susの減少に伴って単調減少する。更に、基準値S
us4(中間値SusB)で負の最小値をとった補正値f1
(Sus)は、基準値Sus4を上回る信号領域では、非鮮
鋭マスク信号Susの増加に伴って単調増加する。
非鮮鋭マスク信号Susの増大に伴ってf1(Sus)が必
ず減少・増加するものである必要はなく、部分的にはS
usが変化してもf1(Sus)が変化しない領域が存在し
ても良いこととする。ここで、前記オリジナル画像信号
Sorg のヒストグラム解析結果が、図3(a)に示すよ
うな特性を有し、関心領域である人体の部分に対して、
信号値が極端に大きな素抜けの部分と、信号値が極端に
小さな照射野絞りの部分とが存在するものとする。
Sorg を補正するための補正値f1(Sus)が、図4に
示すように、基準値Sus2を下回る信号領域で、非鮮鋭
マスク信号Susの減少に応じて必ず増大し、基準値Sus
3を上回る信号領域で、非鮮鋭マスク信号Susの増加に
応じて必ず減少する特性であったとする。この場合、信
号値が大きくなるほど信号Sorg を大きく減少補正し、
また、信号値が小さくなるほど信号Sorg を大きく増大
補正することになるから、素抜けの部分や照射野絞りの
部分の信号値が、関心領域の信号値に大きく近づくこと
になって(図3(b)参照)、関心領域と素抜け或いは
照射野絞りとの信号値の違いを明確に保持することがで
きなくなって好ましくない。
における補正特性では、関心領域における高濃度部及び
低濃度部について圧縮を行わせるが、関心領域を外れた
信号領域では、少なくとも関心領域内における最大補正
度合いを大きく上回る補正が施されることがないから、
素抜けや照射野絞りの部分に相当する信号が必要以上に
補正されてしまうことを回避でき、以て、関心領域と素
抜け或いは照射野絞りとの信号差の違いを保持しつつ関
心領域に対する圧縮を行える(図3(c)参照)。
1(中間値SusA)を下回る領域ではマスク信号Susの
減少に応じて単調減少し、基準値Sus4(中間値Sus
B)を上回る領域ではマスク信号Susの増加に応じて単
調増加する特性としてあるが、基準値Sus1,Sus4か
ら外れた低濃度,高濃度領域で少なくとも補正度合いが
増大変化を示さなければ効果を発揮することになるか
ら、関数形は図2に示す特性に限定されるものではな
く、以下に示すような種々の関数形を用いることができ
る。
us1(中間値SusA)で正の最大値をとった補正値f1
(Sus)は、基準値Sus1を下回る信号領域では、前記
正の最大値に固定され、また、基準値Sus4(中間値S
usB)で負の最小値をとった補正値f1(Sus)は、基
準値Sus4を上回る信号領域では、前記負の最小値に固
定される。
せて、低濃度領域と高濃度領域とを同時に圧縮する構成
としても良い。また、図7及び図8に示す関数形では、
基準値Sus1を下回る信号領域では、補正値f1(Su
s)がマスク信号Susの減少に応じて必ず減少変化し、
また、基準値Sus4を上回る信号領域では、補正値f1
(Sus)がマスク信号の増加に応じて必ず増大変化する
特性としてある。即ち、図2に示す関数形では、基準値
Sus1を下回る信号領域及び基準値Sus4を上回る信号
領域に、マスク信号Susが変化しても補正値f1(Su
s)が変化しない領域が存在したが、図7及び図8に示
す関数形では、変化を示さない領域が存在しないパター
ンを示す。
せて、低濃度領域と高濃度領域とを同時に圧縮する構成
としても良い。また、前記図2に示す関数形では、低濃
度領域と高濃度領域とを同時に圧縮する構成としたが、
図9及び図10に示すように、低濃度領域と高濃度領域と
に関数形を分けて、いずれか一方の濃度領域のみを圧縮
する構成としても良い。
用い、高濃度領域では図8の関数形を用いる構成とする
など、低濃度側と高濃度側とで異なる特性の関数形を組
み合わせて構成しても良い。ところで、上記実施例で
は、基準値Sus1(中間値SusA)及び基準値Sus4
(中間値SusB)を関心領域の信号範囲の外側に設定さ
せるようにしたが、前記基準値Sus1,Sus4(>Sus
1)を関心領域の信号範囲内に設定する構成としても良
く、かかる設定とした実施例を以下に示す。
る式に基づくオリジナル画像信号Sorg の補正におい
て、前記補正値f1(Sus)を図11に示すような関数形
としてある。前記図11に示した関数形を式に示すと、以
下のようになる。
補正値変化の基準となる基準非鮮鋭マスク信号Susi
(i=1〜4)(Sus2=Sus3)は、オリジナル画像
の関心領域内の最大信号値,最小信号値をSmax,Smin
としたときに、 Susi=ki・Smax +(1−ki)・Smin なる式で求めるようにした。そして、前記係数ki(i
=1〜4)は、例えばk1=0.35、k2=k3=0.5 、
k4=0.65とした。また、補正値f1(Sus)の係数β
1 ,β2 =0.8 とした。
4については、これらによって決定されるSus1,Sus
2(Sus3),Sus4が、関心領域の信号範囲の内側に
設定されるようにしてある。上記特性の補正値f1(S
us)を用いてオリジナル画像信号Sorg を補正すれば、
例えば図12に示すような特性のオリジナル画像信号Sor
g を、その低濃度部(信号値の低い側)及び高濃度部の
コントラストを維持しながら、中濃度領域の圧縮を行い
(図13参照)、全体としてダイナミックレンジを圧縮す
ることができる。
値Sus1(中間値SusA)及び負の最小値をとる基準値
Sus4(中間値SusB)を関心領域の信号範囲の内側に
設定する構成においても、前記図7〜図10に示した種々
の関数形を用いることができる。ところで、図14に示す
ような胸部正面画像で、例えば肺の部分に金属製のペー
スメーカがあると、この部分は極端に信号値が低下する
ことになり、信号値が低い側(低濃度側)を、非鮮鋭マ
スク信号Susに基づいて補正する場合には、本来肺野の
部分は補正を行わせたくないのに、前記ペースメーカの
部分については信号値を上げる補正を行うことになって
しまう。
を行わせたくない領域に、信号値が極端に異なる異物が
存在しても、かかる異物に影響された非鮮鋭マスク信号
Susによって補正が加えられてしまうことが防止できる
構成としてある。まず、図14に示すようなオリジナル画
像(胸部正面画像)において、補正特性を使い分ける領
域区分を行う。図14に示す例では、周縁部の照射野絞り
領域(領域1)、素抜け領域(領域2)、肺野領域(領
域3)、その他の被写体領域(領域4)の4領域に分け
ている。
CRT等に表示させ、オペレータがかかる表示を見なが
ら領域を指定する情報を入力する構成であっても良い
し、また、オリジナル画像信号Sorg から画像信号解析
により被写体の輪郭,肺野(内部構造)の輪郭,照射野
絞りの境界を抽出して自動設定させる構成であっても良
い。
信号Sprocを得るものとした場合、領域1,2,3で
は、非鮮鋭マスク信号Susに因らずに補正関数Fを0
(一定値関数f2)とし、肺野を除く関心領域である領
域4では、非鮮鋭マスク信号Susの関数である補正値f
1(Sus)を用いて圧縮を行わせる。
うな関数形に設定してあり、かかる関数形は下式に示さ
れる。
us1,Sus2は、図16(a)に示すように、非鮮鋭マス
ク信号Susの肺野(及び肺領域に設けられたペースメー
カ)を横断する部分のプロファイルをとったときに、基
準値Sus2については、素抜け部よりも低く、かつ、肺
野部の信号レベルよりも高く設定してあり、また、基準
値Sus1については、肺野領域とその他の被写体領域と
を分ける境界レベル付近に設定してある。
スク信号Susの肺野を横断する部分のプロファイルに基
づいて、図15に示すような特性の補正値f1(Sus)を
画像の全領域に適用したとすると、前記図16(a)のプ
ロファイルに対応する補正値f1(Sus)のプロファイ
ルが図16(b)に示すように設定されることになる。即
ち、ペースメーカがなければ、肺野の信号は基準値Sus
1,Sus2で挟まれる信号範囲内に収まるので、補正値
が設定されることがない(補正値=0となる)が、ペー
スメーカが存在するとかかるペースメーカの部分で非鮮
鋭マスク信号Susが基準値Sus1を下回ることによっ
て、ペースメーカの部分の信号を上げる補正値が設定さ
れる。そして、ペースメーカの部分の信号が持ち上げら
れると、肺野部分との信号差が縮まり、読影に不適当と
なってしまう。
初から補正対象領域外としてあるから、たとえペースメ
ーカの部分で高い非鮮鋭マスク信号Susが設定されて
も、これに応じてペースメーカの部分の信号に補正が加
えられることがなく、ペースメーカと肺野領域との信号
差を保持できる。即ち、図16(b)に示す補正値f1
(Sus)のうち、領域1,2,3に該当する部分を0に
した図16(c)に示すような補正値f1(Sus)が設定
されることになるから、ペースメーカに影響されて大き
な非鮮鋭マスク信号Susが演算されても、これに基づい
て補正が行われることはない。
部分も補正対象外としてあるから、素抜けや照射野絞り
の部分に補正が加えられて、人体の部分の信号値に近づ
くことが回避でき、素抜けや照射野絞りの部分と人体の
部分との信号差を保持して、診断性能に優れた画像を提
供できる。尚、前記照射野絞り領域である領域1におい
て、前記補正関数Fを非常に大きな値をとる一定関数と
し、照射野絞りの部分を素抜け部と同じ位の信号値にし
てしまい、照射野絞りの部分も含めて表示させるとき
に、照射野絞りと素抜け部とが略一様に高濃度に表示さ
れ、人体部分を際立たせるようにしても良い。
せながら、補正特性を使い分ける領域を指定するとき
に、補正値f1(Sus)による補正特性を決定する基準
値Sus1,Sus2や係数β1 ,β2 のデータを入力させ
たり、また、補正値f1(Sus)を用いない領域に対し
て与える一定補正値(一定値関数f2)を可変入力させ
る構成としても良く、更には、前記補正データの入力結
果に基づいて得られた処理済み画像信号Sprocを表示さ
せるようにしても良い。
クサイズ或いは周波数特性は画像の診断性を左右する重
要なパラメータである。ダイナミックレンジ圧縮処理で
は、被写体の大まかな構造の変化(肺野部,縦隔部など
の滑らかな信号差)に対応する超低周波数成分のみを非
鮮鋭マスク信号Susとして抽出し、Susに基づいて補正
値f1(Sus)を設定することによって、細かい構造物
の変化(骨,血管など)を維持しつつ、全体の濃度範囲
を圧縮することができるのである。
号Susが被写体の大まかな変化に対応する超低周波数成
分だけでなく、細かい構造物の変化に相当する周波数成
分をも含んでしまい、非鮮鋭マスク信号Susに基づく補
正値を加算することによって細かい構造物の変化が打ち
消されて骨や血管などのコントラストが低下してしま
う。
化が急激な部分での非鮮鋭画像のエッジ切れが悪くな
り、圧縮を行いたい領域と行いたくない領域との境界付
近で望まない圧縮が施されてしまう。また、マスクサイ
ズを更に大きくし過ぎると、被写体の大まかな変化に相
当する周波数成分までも失ってしまう(極端な場合全く
平坦な画像になってしまう)ので、非鮮鋭マスク信号S
usに基づく補正値を加算してもダイナミックレンジ圧縮
効果が得られなくなる。
果、マスクサイズの大きさは等身大画像上の長さで10mm
から60mmが好ましく、より好ましくは15mmから30mmであ
り、最も好ましくは20mmから30mmであることを見出し
た。マスクサイズが10mmよりも小さいと、細かい構造物
の変化に対応する周波数成分が急激に増加するので、こ
のようなマスクサイズで求めた非鮮鋭マスク信号Susに
基づいて補正値を設定すると、著しく診断性能が低下し
てしまう。また、特に、胸部画像や腹部画像において
は、マスクサイズを15mm以上にすれば、Susは大動脈な
どの太い血管に対応する周波数成分を持たなくなり、マ
スクサイズを20mm以上にすればSusは肋骨などに対応す
る比較的低周波数であるが、コントラストを下げたくな
い周波数成分を含まなくなるので、診断性能の高い画像
が得られる。
長さと長辺の長さの平均値、正方形なら一辺の長さ、円
ならば直径、楕円ならば長径と短径の平均値を指す。ま
た、マスクサイズのかわりに非鮮鋭マスクのもつ周波数
特性で記述すると、非鮮鋭マスクの変調伝達関数が0.01
サイクル/mmのとき0.5 以上かつ0.06サイクル/mmのと
き0.5 以下であることが好ましく、より好ましくは0.02
サイクル/mmのとき0.5 以上かつ0.04サイクル/mmのと
き0.5 以下、更に好ましくは0.02サイクル/mmのとき0.
5 以上かつ0.03サイクル/mmのとき0.5 以下である。
usの関数である補正値f1(Sus)の絶対値の最大値は
被写体の関心領域のダイナミックレンジの1/8から1
/2であることが好ましい。例えば、被写体の関心領域
のダイナミックレンジが2桁であるとき、圧縮補正量の
絶対値の最大値は1/4桁から1桁であることが好まし
い。
−Sus)の様に非鮮鋭マスク信号Susの一次関数で表さ
れるとき、その傾きであって圧縮度合いを決めるβの好
ましい範囲は0.2 〜1.0 であり、より好ましくは0.4 〜
0.8 である。補正量が小さすぎるとダイナミックレンジ
圧縮効果が現れず、一方、補正量が大きすぎるとオリジ
ナル画像における領域毎の濃度の大小関係が逆転して
(例えば肺野の平均濃度よりも縦隔部の平均濃度の方が
高くなって)しまい、診断に耐えない画像になってしま
う。例えば、上記一次関数の傾きβを1より大きくした
場合このような問題が起こる。
リジナル画像信号Sorg を処理して前記オリジナル画像
よりもダイナミックレンジの狭い画像を担持する処理済
み画像信号Sprocを得る放射線画像のダイナミックレン
ジ圧縮方法において、ダイナミックレンジの圧縮処理を
施したくない照射野絞り,素抜け,異物の領域が補正さ
れて、関心領域に対する信号差を保持できなくなること
を防止でき、以て、読影性に優れた放射線画像を提供で
きるという効果がある。
示すブロック図。
トグラム線図。
ンを示す線図。
ンを示す線図。
ンを示す線図。
ンを示す線図。
ンを示す線図。
ンを示す線図。
す線図。
示す線図。
図。
プロファイル線図。
Claims (6)
- 【請求項1】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Sorg を処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Sprocを得る放射線
画像のダイナミックレンジ圧縮方法であって、 各画素点に対応して非鮮鋭マスク信号Susを求め、該非
鮮鋭マスク信号Susの関数f1(Sus)を補正値とした
ときに、 Sproc=Sorg +f1(Sus) なる演算を行って前記処理済み画像信号Sprocを得る放
射線画像のダイナミックレンジ圧縮方法において、 前記補正値f1(Sus)が、前記非鮮鋭マスク信号Sus
の所定中間値SusAで正の最大値をとり、前記所定中間
値SusAを上回る領域では非鮮鋭マスク信号Susの増大
に応じて単調減少し、前記所定中間値SusAを下回る領
域では前記正の最大値に固定されるか又は非鮮鋭マスク
信号Susの減少に応じて単調減少することを特徴とする
放射線画像のダイナミックレンジ圧縮方法。 - 【請求項2】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Sorg を処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Sprocを得る放射線
画像のダイナミックレンジ圧縮方法であって、 各画素点に対応して非鮮鋭マスク信号Susを求め、該非
鮮鋭マスク信号Susの関数f1(Sus)を補正値とした
ときに、 Sproc=Sorg +f1(Sus) なる演算を行って前記処理済み画像信号Sprocを得る放
射線画像のダイナミックレンジ圧縮方法において、 前記補正値f1(Sus)が、前記非鮮鋭マスク信号Sus
の所定中間値SusBで負の最小値をとり、前記所定中間
値SusBを下回る領域では非鮮鋭マスク信号Susの増大
に応じて単調減少し、前記所定中間値SusBを上回る領
域では前記負の最小値に固定されるか又は非鮮鋭マスク
信号Susの増大に応じて単調増加することを特徴とする
放射線画像のダイナミックレンジ圧縮方法。 - 【請求項3】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Sorg を処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Sprocを得る放射線
画像のダイナミックレンジ圧縮方法であって、 各画素点に対応して非鮮鋭マスク信号Susを求め、該非
鮮鋭マスク信号Susの関数f1(Sus)を補正値とした
ときに、 Sproc=Sorg +f1(Sus) なる演算を行って前記処理済み画像信号Sprocを得る放
射線画像のダイナミックレンジ圧縮方法において、 前記非鮮鋭マスク信号Susの2つの異なる中間値をSus
A,SusB(SusA<SusB)としたときに、前記補正
値f1(Sus)が、前記中間値SusAで正の最大値をと
り、かつ、前記中間値SusBで負の最小値をとる一方、
前記2つの中間値SusA,SusBで挟まれる信号領域で
は非鮮鋭マスク信号Susの増大に応じて単調減少し、中
間値SusAを下回る領域では前記正の最大値に固定され
るか又は非鮮鋭マスク信号Susの減少に応じて単調減少
し、前記中間値SusBを上回る領域では前記負の最小値
に固定されるか又は非鮮鋭マスク信号Susの増大に応じ
て単調増加することを特徴とする放射線画像のダイナミ
ックレンジ圧縮方法。 - 【請求項4】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Sorg を処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Sprocを得る放射線
画像のダイナミックレンジ圧縮方法であって、 補正関数をFとしたときに、 Sproc=Sorg +F なる演算を行ってダイナミックレンジが圧縮された画像
信号Sprocを得る放射線画像のダイナミックレンジ圧縮
方法において、 前記補正関数Fとして予め異なる複数種の関数を備え、
画像内の複数領域に対して異なる補正関数Fを定めるこ
とを特徴とする放射線画像のダイナミックレンジ圧縮方
法。 - 【請求項5】前記画像内の複数領域が前記被写体の輪郭
と被写体の内部構造の輪郭と照射野絞りの境界とのうち
の少なくとも1つによって区分されることを特徴とする
請求項4記載の放射線画像のダイナミックレンジ圧縮方
法。 - 【請求項6】前記複数の補正関数Fが、各画素点に対応
する非鮮鋭マスク信号Susの増大に応じて単調減少する
関数f1(Sus)と、前記非鮮鋭マスク信号Susに因ら
ない関数f2とを少なくとも含むことを特徴とする請求
項4又は5のいずれかに記載の放射線画像のダイナミッ
クレンジ圧縮方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07590093A JP3434532B2 (ja) | 1993-04-01 | 1993-04-01 | 放射線画像のダイナミックレンジ圧縮方法 |
| US08/215,261 US5471987A (en) | 1993-03-30 | 1994-03-21 | Method of compressing a dynamic range for a radiation image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07590093A JP3434532B2 (ja) | 1993-04-01 | 1993-04-01 | 放射線画像のダイナミックレンジ圧縮方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06292013A true JPH06292013A (ja) | 1994-10-18 |
| JP3434532B2 JP3434532B2 (ja) | 2003-08-11 |
Family
ID=13589676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP07590093A Expired - Lifetime JP3434532B2 (ja) | 1993-03-30 | 1993-04-01 | 放射線画像のダイナミックレンジ圧縮方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3434532B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0766202A3 (en) * | 1995-09-29 | 1998-01-14 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| US5796870A (en) * | 1995-04-20 | 1998-08-18 | Fuji Photo Film Co., Ltd. | Method and apparatus for compressing dynamic ranges of images |
| JP2000155838A (ja) * | 1998-09-18 | 2000-06-06 | Canon Inc | 画像処理装置、方法及びコンピュ―タ読み取り可能な記憶媒体 |
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルx線イメージング検出器のダイナミックレンジ拡張システム |
| US6993183B2 (en) | 2001-03-13 | 2006-01-31 | Canon Kabushiki Kaisha | Image processing apparatus, image processing system, image processing method, and computer-readable storage medium |
-
1993
- 1993-04-01 JP JP07590093A patent/JP3434532B2/ja not_active Expired - Lifetime
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5796870A (en) * | 1995-04-20 | 1998-08-18 | Fuji Photo Film Co., Ltd. | Method and apparatus for compressing dynamic ranges of images |
| EP0766202A3 (en) * | 1995-09-29 | 1998-01-14 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| EP1156451A3 (en) * | 1995-09-29 | 2002-01-23 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| JP2000155838A (ja) * | 1998-09-18 | 2000-06-06 | Canon Inc | 画像処理装置、方法及びコンピュ―タ読み取り可能な記憶媒体 |
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルx線イメージング検出器のダイナミックレンジ拡張システム |
| US6993183B2 (en) | 2001-03-13 | 2006-01-31 | Canon Kabushiki Kaisha | Image processing apparatus, image processing system, image processing method, and computer-readable storage medium |
| US7310444B2 (en) | 2001-03-13 | 2007-12-18 | Canon Kabushiki Kaisha | Image processing apparatus, image processing system, image processing method, and computer-readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3434532B2 (ja) | 2003-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5471987A (en) | Method of compressing a dynamic range for a radiation image | |
| JP3467285B2 (ja) | 放射線画像処理方法 | |
| EP0440166B1 (en) | Method for compressing dynamic ranges of images | |
| EP0032237B1 (en) | Method of and apparatus for processing a radiographic image | |
| US5319719A (en) | Processing apparatus for radiographic image signals | |
| US4804842A (en) | Radiation image read-out method and apparatus | |
| US5067163A (en) | Method for determining a desired image signal range from an image having a single background | |
| EP0442468A1 (en) | Method for forming energy substraction images | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| US6125166A (en) | Method of forming energy subtraction images | |
| JP2000033082A (ja) | 放射線画像の画像処理装置 | |
| US5291403A (en) | Method of and apparatus for processing radiation image | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| US5301107A (en) | Method for forming energy subtraction images | |
| JPS63183435A (ja) | 画像処理条件決定方法 | |
| JPH06292008A (ja) | 放射線画像のダイナミックレンジ圧縮処理装置 | |
| JP3434532B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP3814864B2 (ja) | 放射線画像の画像処理条件決定装置及び画像処理装置並びに放射線画像の画像処理条件決定方法及び画像処理方法 | |
| JP3455566B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JPH06339025A (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| US7483556B2 (en) | Energy subtraction processing method and apparatus | |
| JP3637049B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP2852794B2 (ja) | デジタル放射線画像信号の処理装置 | |
| JP4258092B2 (ja) | 画像処理装置および画像処理方法 | |
| JP3185105B2 (ja) | 放射線画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| 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: 20080530 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090530 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090530 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100530 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110530 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110530 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120530 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120530 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130530 Year of fee payment: 10 |