JP3425780B2 - 放射線画像のダイナミックレンジ圧縮方法 - Google Patents
放射線画像のダイナミックレンジ圧縮方法Info
- Publication number
- JP3425780B2 JP3425780B2 JP07649593A JP7649593A JP3425780B2 JP 3425780 B2 JP3425780 B2 JP 3425780B2 JP 07649593 A JP07649593 A JP 07649593A JP 7649593 A JP7649593 A JP 7649593A JP 3425780 B2 JP3425780 B2 JP 3425780B2
- Authority
- JP
- Japan
- Prior art keywords
- sus
- image signal
- dynamic range
- image
- signal
- 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.)
- Expired - Fee Related
Links
- 238000007906 compression Methods 0.000 title claims description 30
- 230000006835 compression Effects 0.000 title claims description 29
- 238000000034 method Methods 0.000 title claims description 20
- 238000012937 correction Methods 0.000 claims description 40
- 230000005855 radiation Effects 0.000 claims description 35
- 238000012545 processing Methods 0.000 claims description 18
- 230000007423 decrease Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 description 25
- 230000006870 function Effects 0.000 description 21
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 5
- 238000012886 linear function Methods 0.000 description 5
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 4
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000004020 luminiscence type Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000007740 vapor deposition Methods 0.000 description 1
- 238000012800 visualization Methods 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/73—Deblurring; Sharpening
- G06T5/75—Unsharp masking
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Compression Of Band Width Or Redundancy In Fax (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Picture Signal Circuits (AREA)
- Closed-Circuit Television Systems (AREA)
- Facsimile Image Signal Circuits (AREA)
Description
【０００１】
【産業上の利用分野】本発明は放射線画像のダイナミッ
クレンジ圧縮方法に関し、詳しくは、オリジナル画像信
号を処理して該オリジナル画像よりもダイナミックレン
ジの狭い画像を担持する画像信号を得る放射線画像のダ
イナミックレンジ圧縮方法の改善技術に関する。
クレンジ圧縮方法に関し、詳しくは、オリジナル画像信
号を処理して該オリジナル画像よりもダイナミックレン
ジの狭い画像を担持する画像信号を得る放射線画像のダ
イナミックレンジ圧縮方法の改善技術に関する。
【０００２】
【従来の技術】従来、放射線画像において、画像領域内
の微細構造の観察適正を確保しつつ、濃度域を圧縮する
方法として、例えば特開平３−２２２５７７号公報に開
示されるようなものがある。
の微細構造の観察適正を確保しつつ、濃度域を圧縮する
方法として、例えば特開平３−２２２５７７号公報に開
示されるようなものがある。
【０００３】前記特開平３−２２２５７７号公報に開示
される圧縮方法は、各画素点に対応して該各画素点を含
む所定マスク領域内のオリジナル画像信号Ｓorgを平均
化することで非鮮鋭マスク信号（ボケマスク信号）Ｓus
を求め、この非鮮鋭マスク信号Ｓusが増大するにつれて
単調減少する関数をｆ（Ｓus）としたときに、ダイナミ
ックレンジが圧縮された処理済み画像信号Ｓprocを、Ｓ
proc＝Ｓorg ＋ｆ（Ｓus）として得るものである。
される圧縮方法は、各画素点に対応して該各画素点を含
む所定マスク領域内のオリジナル画像信号Ｓorgを平均
化することで非鮮鋭マスク信号（ボケマスク信号）Ｓus
を求め、この非鮮鋭マスク信号Ｓusが増大するにつれて
単調減少する関数をｆ（Ｓus）としたときに、ダイナミ
ックレンジが圧縮された処理済み画像信号Ｓprocを、Ｓ
proc＝Ｓorg ＋ｆ（Ｓus）として得るものである。
【０００４】
【発明が解決しようとする課題】ところで、上記のよう
なダイナミックレンジ圧縮においては、最終的にハード
コピーを得たときなどの可視化状態で、読影に最適なコ
ントラスト特性を関心領域において得ることが目的であ
るが、従来のダイナミックレンジ圧縮方法によると、入
力されたディジタル画像信号（オリジナル画像信号）に
基づいて補正度合いを決定していたために以下のような
問題が発生していた。
なダイナミックレンジ圧縮においては、最終的にハード
コピーを得たときなどの可視化状態で、読影に最適なコ
ントラスト特性を関心領域において得ることが目的であ
るが、従来のダイナミックレンジ圧縮方法によると、入
力されたディジタル画像信号（オリジナル画像信号）に
基づいて補正度合いを決定していたために以下のような
問題が発生していた。
【０００５】即ち、入力されたディジタル信号値に対し
て一定の補正関数（例えば１次関数）を設定しても、ハ
ードコピーを得るときの記録材料等における階調特性に
よっては、最終的に可視化されたときにダイナミックレ
ンジ圧縮の補正度合いが変わってしまう（１次関数でな
くなってしまう）ことがあり、最終的な表示状態（ハー
ドコピー上）に対して最適なダイナミックレンジ圧縮処
理を、安定的に施すことができなかった。
て一定の補正関数（例えば１次関数）を設定しても、ハ
ードコピーを得るときの記録材料等における階調特性に
よっては、最終的に可視化されたときにダイナミックレ
ンジ圧縮の補正度合いが変わってしまう（１次関数でな
くなってしまう）ことがあり、最終的な表示状態（ハー
ドコピー上）に対して最適なダイナミックレンジ圧縮処
理を、安定的に施すことができなかった。
【０００６】本発明は上記問題点に鑑みなされたもので
あり、非鮮鋭マスク信号を用いたダイナミックレンジの
圧縮方法において、最終的な表示状態（ハードコピー
上）に対して最適なダイナミックレンジ圧縮処理を、安
定的に施すことができるようにすることを目的とする。
あり、非鮮鋭マスク信号を用いたダイナミックレンジの
圧縮方法において、最終的な表示状態（ハードコピー
上）に対して最適なダイナミックレンジ圧縮処理を、安
定的に施すことができるようにすることを目的とする。
【０００７】
【課題を解決するための手段】そのため本発明にかかる
放射線画像のダイナミックレンジ圧縮方法では、ダイナ
ミックレンジが圧縮された処理済み画像信号Ｓprocを得
る画像処理を行った後に、該処理済み画像信号Ｓprocを
表示装置に出力する構成において、各画素点に対応して
非鮮鋭マスク信号Ｓusを求める一方、前記オリジナル画
像信号Ｓorgに対応する前記表示装置における表示濃度
をＤorg、前記処理済み画像信号Ｓprocに対応する前記
表示装置における表示濃度をＤproc、前記非鮮鋭マスク
信号Ｓusに対応する前記表示装置における表示濃度をＤ
usとし、前記表示濃度Ｄusの関数である濃度補正値をｆ
２（Ｄus）としたときに、Ｄproc＝Ｄorg＋ｆ２（Ｄu
s）なる関係を満たすように前記処理済み画像信号Ｓpro
cを得るようにした。
放射線画像のダイナミックレンジ圧縮方法では、ダイナ
ミックレンジが圧縮された処理済み画像信号Ｓprocを得
る画像処理を行った後に、該処理済み画像信号Ｓprocを
表示装置に出力する構成において、各画素点に対応して
非鮮鋭マスク信号Ｓusを求める一方、前記オリジナル画
像信号Ｓorgに対応する前記表示装置における表示濃度
をＤorg、前記処理済み画像信号Ｓprocに対応する前記
表示装置における表示濃度をＤproc、前記非鮮鋭マスク
信号Ｓusに対応する前記表示装置における表示濃度をＤ
usとし、前記表示濃度Ｄusの関数である濃度補正値をｆ
２（Ｄus）としたときに、Ｄproc＝Ｄorg＋ｆ２（Ｄu
s）なる関係を満たすように前記処理済み画像信号Ｓpro
cを得るようにした。
【０００８】また、非鮮鋭マスク信号Ｓusに応じた補正
値ｆ１（Ｓus）によってオリジナル画像信号Ｓorgを補
正することで処理済み画像信号Ｓprocを得る放射線画像
のダイナミックレンジ圧縮方法において、前記非鮮鋭マ
スク信号Ｓusの変化に対する前記補正値ｆ１（Ｓus）の
変化割合の絶対値が、オリジナル画像信号Ｓorgを増大
補正する補正値ｆ１（Ｓus）については、非鮮鋭マスク
信号Ｓusの減少に応じて徐々に増大変化し、オリジナル
画像信号Ｓorgを減少補正する補正値ｆ１（Ｓus）につ
いては、非鮮鋭マスク信号Ｓusの増大に応じて徐々に増
大変化するようにした。
値ｆ１（Ｓus）によってオリジナル画像信号Ｓorgを補
正することで処理済み画像信号Ｓprocを得る放射線画像
のダイナミックレンジ圧縮方法において、前記非鮮鋭マ
スク信号Ｓusの変化に対する前記補正値ｆ１（Ｓus）の
変化割合の絶対値が、オリジナル画像信号Ｓorgを増大
補正する補正値ｆ１（Ｓus）については、非鮮鋭マスク
信号Ｓusの減少に応じて徐々に増大変化し、オリジナル
画像信号Ｓorgを減少補正する補正値ｆ１（Ｓus）につ
いては、非鮮鋭マスク信号Ｓusの増大に応じて徐々に増
大変化するようにした。
【０００９】
【作用】かかる構成によると、ダイナミックレンジ圧縮
処理を施した処理済みの画像信号Ｓprocを表示装置に出
力する構成において、かかる表示装置における表示濃度
上で所期のダイナミックレンジ圧縮が行われるように、
表示濃度ベースで要求されるダイナミックレンジ圧縮を
設定し、かかる要求に見合った処理済み画像信号Ｓproc
を得る構成とした。
処理を施した処理済みの画像信号Ｓprocを表示装置に出
力する構成において、かかる表示装置における表示濃度
上で所期のダイナミックレンジ圧縮が行われるように、
表示濃度ベースで要求されるダイナミックレンジ圧縮を
設定し、かかる要求に見合った処理済み画像信号Ｓproc
を得る構成とした。
【００１０】また、前記表示装置における階調特性は、
一般に低濃度及び高濃度側ほどコントラストが低下する
から、かかる階調特性に見合った補正特性とすることが
望まれ、そのために、非鮮鋭マスク信号Ｓusの変化に対
して一定割合で補正値ｆ１（Ｓus）が変化するのではな
く、低濃度又は高濃度側に行くに従ってより大きな変化
割合で補正値ｆ１（Ｓus）を変化させるようにした。
一般に低濃度及び高濃度側ほどコントラストが低下する
から、かかる階調特性に見合った補正特性とすることが
望まれ、そのために、非鮮鋭マスク信号Ｓusの変化に対
して一定割合で補正値ｆ１（Ｓus）が変化するのではな
く、低濃度又は高濃度側に行くに従ってより大きな変化
割合で補正値ｆ１（Ｓus）を変化させるようにした。
【００１１】
【実施例】以下に本発明の実施例を説明する。一実施例
を示す図１は、本発明にかかる放射線画像のダイナミッ
クレンジ圧縮方法が適用された画像処理装置を含む放射
線画像読取装置を示すものであり、医療用としての人体
の撮影を行う例を示す。
を示す図１は、本発明にかかる放射線画像のダイナミッ
クレンジ圧縮方法が適用された画像処理装置を含む放射
線画像読取装置を示すものであり、医療用としての人体
の撮影を行う例を示す。
【００１２】ここで、放射線発生源１は、放射線制御装
置２によって制御されて、被写体（人体胸部等）Ｍに向
けて放射線（一般的にはＸ線）を照射する。記録読取装
置３は、被写体Ｍを挟んで放射線源１と対向する面に変
換パネル４を備え、該変換パネル４は放射線源１からの
照射放射線量に対する被写体Ｍの放射線透過率分布に従
ったエネルギーを輝尽層に蓄積し、そこに被写体Ｍの潜
像を形成する。
置２によって制御されて、被写体（人体胸部等）Ｍに向
けて放射線（一般的にはＸ線）を照射する。記録読取装
置３は、被写体Ｍを挟んで放射線源１と対向する面に変
換パネル４を備え、該変換パネル４は放射線源１からの
照射放射線量に対する被写体Ｍの放射線透過率分布に従
ったエネルギーを輝尽層に蓄積し、そこに被写体Ｍの潜
像を形成する。
【００１３】前記変換パネル４は、支持体上に輝尽層
を、輝尽性蛍光体の気相堆積、或いは輝尽性蛍光体塗料
塗布によって設けてあり、該輝尽層は環境による悪影響
及び損傷を遮断するために保護部材によって遮蔽若しく
は被覆される。該輝尽性蛍光体材料としては、例えば、
特開昭６１−７２０９１号公報、或いは、特開昭５９−
７５２００号公報に開示されるような材料が使われる。
を、輝尽性蛍光体の気相堆積、或いは輝尽性蛍光体塗料
塗布によって設けてあり、該輝尽層は環境による悪影響
及び損傷を遮断するために保護部材によって遮蔽若しく
は被覆される。該輝尽性蛍光体材料としては、例えば、
特開昭６１−７２０９１号公報、或いは、特開昭５９−
７５２００号公報に開示されるような材料が使われる。
【００１４】光ビーム発生部（ガスレーザ，固体レー
ザ，半導体レーザ等）５は、出射強度が制御された光ビ
ームを発生し、その光ビームは種々の光学系を経由して
走査器６に到達し、そこで偏向を受け、更に、反射鏡７
で光路を偏向させて、変換パネル４に輝尽励起走査光と
して導かれる。
ザ，半導体レーザ等）５は、出射強度が制御された光ビ
ームを発生し、その光ビームは種々の光学系を経由して
走査器６に到達し、そこで偏向を受け、更に、反射鏡７
で光路を偏向させて、変換パネル４に輝尽励起走査光と
して導かれる。
【００１５】集光体８は、輝尽励起光が走査される変換
パネル４に近接して光ファイバである集光端が位置さ
れ、上記光ビームで走査された変換パネル４からの潜像
エネルギーに比例した発光強度の輝尽発光を受光する。
パネル４に近接して光ファイバである集光端が位置さ
れ、上記光ビームで走査された変換パネル４からの潜像
エネルギーに比例した発光強度の輝尽発光を受光する。
【００１６】９は、集光体８から導入された光から輝尽
発光波長領域の光のみを通過させるフィルタであり、該
フィルタ９を通過した光は、フォトマル10に入射して、
その入射光に対応した電流信号に光電変換される。
発光波長領域の光のみを通過させるフィルタであり、該
フィルタ９を通過した光は、フォトマル10に入射して、
その入射光に対応した電流信号に光電変換される。
【００１７】フォトマル10からの出力電流は、電流／電
圧変換器11で電圧信号に変換され、増幅器12で増幅され
た後、Ａ／Ｄ変換器13で画素毎のディジタルデータから
なる放射線画像信号に変換される。そして、このディジ
タル放射線画像信号（オリジナル画像信号Ｓorg）は、
マイクロコンピュータを内蔵した画像処理装置14に順次
出力される。
圧変換器11で電圧信号に変換され、増幅器12で増幅され
た後、Ａ／Ｄ変換器13で画素毎のディジタルデータから
なる放射線画像信号に変換される。そして、このディジ
タル放射線画像信号（オリジナル画像信号Ｓorg）は、
マイクロコンピュータを内蔵した画像処理装置14に順次
出力される。
【００１８】15は画像信号を記憶させておくための画像
メモリ（磁気ディスク装置）である。また、16は画像処
理装置14から直接又は前記画像メモリ15から読み出され
た放射線画像信号をプリンタ17に伝送するためのインタ
ーフェイスである。
メモリ（磁気ディスク装置）である。また、16は画像処
理装置14から直接又は前記画像メモリ15から読み出され
た放射線画像信号をプリンタ17に伝送するためのインタ
ーフェイスである。
【００１９】18は読取ゲイン調整回路であり、この読取
ゲイン調整回路18により光ビーム発生部５の光ビーム強
度調整、フォトマル用高圧電源19の電源電圧調整による
フォトマル10のゲイン調整、電流／電圧変換器11と増幅
器12のゲイン調整、及びＡ／Ｄ変換器13の入力ダイナミ
ックレンジの調整が行われ、放射線画像信号の読取ゲイ
ンが総合的に調整される。
ゲイン調整回路18により光ビーム発生部５の光ビーム強
度調整、フォトマル用高圧電源19の電源電圧調整による
フォトマル10のゲイン調整、電流／電圧変換器11と増幅
器12のゲイン調整、及びＡ／Ｄ変換器13の入力ダイナミ
ックレンジの調整が行われ、放射線画像信号の読取ゲイ
ンが総合的に調整される。
【００２０】尚、前記画像処理装置14に出力させるオリ
ジナル放射線画像信号Ｓorgの取得方法を、輝尽性蛍光
体を励起光で走査して発光せしめた輝尽発光を光電変換
して得る方法に限定するものではなく、例えば放射線フ
ィルムの画像を、光電変換によって読み取る方法や、被
写体を透過した放射線を蛍光体に照射されて蛍光に変換
し、該蛍光を光電変換して読み取る方法などであっても
良い。
ジナル放射線画像信号Ｓorgの取得方法を、輝尽性蛍光
体を励起光で走査して発光せしめた輝尽発光を光電変換
して得る方法に限定するものではなく、例えば放射線フ
ィルムの画像を、光電変換によって読み取る方法や、被
写体を透過した放射線を蛍光体に照射されて蛍光に変換
し、該蛍光を光電変換して読み取る方法などであっても
良い。
【００２１】オリジナル放射線画像信号Ｓorgは、検出
された放射線の強度に比例する形でも良いし、検出され
た放射線の強度の対数に比例する形でも良いが、後者の
方が好ましい。
された放射線の強度に比例する形でも良いし、検出され
た放射線の強度の対数に比例する形でも良いが、後者の
方が好ましい。
【００２２】ここで、前記画像処理装置14には、入力さ
れるオリジナル画像信号Ｓorgのダイナミックレンジを
圧縮して、オリジナル画像よりもダイナミックレンジの
狭い画像を担持する処理済みの画像信号Ｓprocを得る画
像処理機能が備えられており、かかるダイナミックレン
ジ圧縮のための画像処理は、以下の式に従って行われ
る。
れるオリジナル画像信号Ｓorgのダイナミックレンジを
圧縮して、オリジナル画像よりもダイナミックレンジの
狭い画像を担持する処理済みの画像信号Ｓprocを得る画
像処理機能が備えられており、かかるダイナミックレン
ジ圧縮のための画像処理は、以下の式に従って行われ
る。
【００２３】Ｓproc＝Ｓorg＋ｆ１（Ｓus）
【００２４】上式で、Ｓusは、各画素点に対応して該各
画素点を含む所定マスク領域内のオリジナル画像信号Ｓ
org を平均化することにより求められる非鮮鋭マスク信
号である。また、オリジナル画像信号Ｓorgに加算され
るｆ１（Ｓus）は、非鮮鋭マスク信号Ｓusの関数として
求められる補正値であり、非鮮鋭マスク信号Ｓusの増大
に応じて単調減少する。
画素点を含む所定マスク領域内のオリジナル画像信号Ｓ
org を平均化することにより求められる非鮮鋭マスク信
号である。また、オリジナル画像信号Ｓorgに加算され
るｆ１（Ｓus）は、非鮮鋭マスク信号Ｓusの関数として
求められる補正値であり、非鮮鋭マスク信号Ｓusの増大
に応じて単調減少する。
【００２５】尚、「単調減少」とは、非鮮鋭マスク信号
Ｓusの増大に伴ってｆ１（Ｓus）が必ず減少するもので
ある必要はなく、部分的にはＳusが変化してもｆ１（Ｓ
us）が変化しない領域が存在しても良いこととする。
Ｓusの増大に伴ってｆ１（Ｓus）が必ず減少するもので
ある必要はなく、部分的にはＳusが変化してもｆ１（Ｓ
us）が変化しない領域が存在しても良いこととする。
【００２６】ところで、ダイナミックレンジの圧縮にお
いては、読影に最適な濃度再現が行われることが望ま
れ、係る濃度再現特性とは無関係にダイナミックレンジ
圧縮を行った場合には、最終的に表示を行わせたときに
所期の圧縮効果が得られない場合がある。
いては、読影に最適な濃度再現が行われることが望ま
れ、係る濃度再現特性とは無関係にダイナミックレンジ
圧縮を行った場合には、最終的に表示を行わせたときに
所期の圧縮効果が得られない場合がある。
【００２７】そこで、以下に濃度の再現特性を考慮して
ダイナミックレンジ圧縮を行う実施例を、図２のフロー
チャートに従って説明する。
ダイナミックレンジ圧縮を行う実施例を、図２のフロー
チャートに従って説明する。
【００２８】まず、オリジナル画像信号Ｓorgに基づき
画像内の関心領域のヒストグラム解析を行い、予め設定
されている標準の階調変換テーブルＬＵＴを前記ヒスト
グラム解析の結果に応じて回転させたり、平行移動させ
ることによって、放射線画像毎に最適の階調変換テーブ
ルＬＵＴを決定する（Ｓ１）。
画像内の関心領域のヒストグラム解析を行い、予め設定
されている標準の階調変換テーブルＬＵＴを前記ヒスト
グラム解析の結果に応じて回転させたり、平行移動させ
ることによって、放射線画像毎に最適の階調変換テーブ
ルＬＵＴを決定する（Ｓ１）。
【００２９】次に、前記ヒストグラム解析結果（例えば
関心領域内における最大値，最小値，中央値など）に基
づいて、ダイナミックレンジ圧縮の領域を設定すること
になる非鮮鋭マスク信号Ｓusの基準値Ｓus１を決定する
（Ｓ２）。そして、前記階調変換テーブルＬＵＴを用い
て前記基準値Ｓus１に対応する表示濃度Ｄus１を求める
（Ｓ３）。
関心領域内における最大値，最小値，中央値など）に基
づいて、ダイナミックレンジ圧縮の領域を設定すること
になる非鮮鋭マスク信号Ｓusの基準値Ｓus１を決定する
（Ｓ２）。そして、前記階調変換テーブルＬＵＴを用い
て前記基準値Ｓus１に対応する表示濃度Ｄus１を求める
（Ｓ３）。
【００３０】次に、前記表示濃度Ｄus１を基準として、
以下の１次関数で表される表示濃度補正関数ｆ２（Ｄu
s）を決定する（Ｓ４及び図３（ｂ）参照）。
以下の１次関数で表される表示濃度補正関数ｆ２（Ｄu
s）を決定する（Ｓ４及び図３（ｂ）参照）。
【００３１】ｆ２（Ｄus）＝β（Ｄus１−Ｄus）…（Ｄ
us≦Ｄus１）、ｆ２（Ｄus）＝０…（Ｄus＞Ｄus１）
us≦Ｄus１）、ｆ２（Ｄus）＝０…（Ｄus＞Ｄus１）
【００３２】ここで、オリジナル画像信号Ｓorgに対応
する表示濃度をＤorg、前記処理済み画像信号Ｓprocに
対応する表示濃度をＤprocとしたときに、Ｄproc＝Ｄor
g＋ｆ２（Ｄus）を満足させる表示濃度Ｄorgと表示濃度
Ｄprocとの関係を求める（図３（ａ）における第３象
限）。前記表示濃度Ｄprocは、表示濃度ベースで、濃度
補正値ｆ２（Ｄus）によってオリジナル画像に対応する
表示濃度Ｄorgが低濃度側で圧縮されたデータである。
する表示濃度をＤorg、前記処理済み画像信号Ｓprocに
対応する表示濃度をＤprocとしたときに、Ｄproc＝Ｄor
g＋ｆ２（Ｄus）を満足させる表示濃度Ｄorgと表示濃度
Ｄprocとの関係を求める（図３（ａ）における第３象
限）。前記表示濃度Ｄprocは、表示濃度ベースで、濃度
補正値ｆ２（Ｄus）によってオリジナル画像に対応する
表示濃度Ｄorgが低濃度側で圧縮されたデータである。
【００３３】そして、最終的には、ダイナミックレンジ
圧縮された処理済み画像信号Ｓprocを、階調変換テーブ
ルＬＵＴで変換して表示させたときに、前記表示濃度Ｄ
procに合致する特性となれば良いことになる。
圧縮された処理済み画像信号Ｓprocを、階調変換テーブ
ルＬＵＴで変換して表示させたときに、前記表示濃度Ｄ
procに合致する特性となれば良いことになる。
【００３４】そこで、前記処理済み画像Ｓprocに対応す
る表示濃度Ｄproc及びオリジナル画像Ｓorgに対応する
表示濃度Ｄorgを、それぞれ前記階調変換テーブルＬＵ
Ｔによる逆変換によって信号値に変換し、所定の階調変
換テーブルＬＵＴにおいて、Ｄproc＝Ｄorg ＋ｆ２（Ｄ
us）なる関係を満足させるオリジナル画像信号Ｓorgと
処理済み画像信号Ｓprocとの関係を求める（図３（ａ）
における第１象限）。
る表示濃度Ｄproc及びオリジナル画像Ｓorgに対応する
表示濃度Ｄorgを、それぞれ前記階調変換テーブルＬＵ
Ｔによる逆変換によって信号値に変換し、所定の階調変
換テーブルＬＵＴにおいて、Ｄproc＝Ｄorg ＋ｆ２（Ｄ
us）なる関係を満足させるオリジナル画像信号Ｓorgと
処理済み画像信号Ｓprocとの関係を求める（図３（ａ）
における第１象限）。
【００３５】そして、前記表示濃度から逆変換して求め
た信号の関係（図３（ａ）における第１象限）におい
て、Ｓproc＝Ｓorg＋ｆ１（Ｓus）を満たす信号補正関
数ｆ１（Ｓus）を算出する（Ｓ５及び図３（ｃ）参
照）。
た信号の関係（図３（ａ）における第１象限）におい
て、Ｓproc＝Ｓorg＋ｆ１（Ｓus）を満たす信号補正関
数ｆ１（Ｓus）を算出する（Ｓ５及び図３（ｃ）参
照）。
【００３６】次いで、オリジナル画像信号Ｓorgから各
画素点毎の非鮮鋭マスク信号Ｓusを算出し（Ｓ６）、該
算出した非鮮鋭マスク信号Ｓusを上記のようにして求め
た補正関数ｆ１（Ｓus）によって補正値に変換し、該補
正値ｆ１（Ｓus）によってオリジナル画像信号Ｓorgを
補正することで、処理済み画像信号Ｓprocを得る（Ｓ
７）。
画素点毎の非鮮鋭マスク信号Ｓusを算出し（Ｓ６）、該
算出した非鮮鋭マスク信号Ｓusを上記のようにして求め
た補正関数ｆ１（Ｓus）によって補正値に変換し、該補
正値ｆ１（Ｓus）によってオリジナル画像信号Ｓorgを
補正することで、処理済み画像信号Ｓprocを得る（Ｓ
７）。
【００３７】前記補正関数ｆ１（Ｓus）によってダイナ
ミックレンジが圧縮された処理済み画像信号Ｓprocは、
前記階調変換テーブルＬＵＴによって変換された後（Ｓ
８）、プリンタ（表示装置）に出力されてハードコピー
が作成される（Ｓ９）。
ミックレンジが圧縮された処理済み画像信号Ｓprocは、
前記階調変換テーブルＬＵＴによって変換された後（Ｓ
８）、プリンタ（表示装置）に出力されてハードコピー
が作成される（Ｓ９）。
【００３８】尚、表示濃度上で、図３（ｂ）に示すよう
な１次関数で表される補正特性を設定し、図３（ａ）の
第２，第４象限に示すような階調変換テーブルＬＵＴを
用いる場合には、図３（ｃ）に示すように、オリジナル
画像信号Ｓorg の低濃度側を補正するための補正関数ｆ
１（Ｓus）としては、非鮮鋭マスク信号Ｓusが基準点Ｓ
us１から減少するに従ってその傾きを増すような特性と
なる。これは、前記階調変換テーブルＬＵＴが、一般的
に低濃度及び高濃度側で傾きが徐々に減少する特性によ
るものである。
な１次関数で表される補正特性を設定し、図３（ａ）の
第２，第４象限に示すような階調変換テーブルＬＵＴを
用いる場合には、図３（ｃ）に示すように、オリジナル
画像信号Ｓorg の低濃度側を補正するための補正関数ｆ
１（Ｓus）としては、非鮮鋭マスク信号Ｓusが基準点Ｓ
us１から減少するに従ってその傾きを増すような特性と
なる。これは、前記階調変換テーブルＬＵＴが、一般的
に低濃度及び高濃度側で傾きが徐々に減少する特性によ
るものである。
【００３９】上記実施例のように、処理済み画像信号Ｓ
procが出力されるプリンタやＣＲＴなどの表示装置の濃
度再現特性を考慮したダイナミックレンジ圧縮を行え
ば、画像信号上では所期の圧縮が行われているのに、表
示画像上では充分な圧縮効果を得ることができないとい
うことがなく、最終的な表示状態で最適となるダイナミ
ックレンジ圧縮を行わせることができる。
procが出力されるプリンタやＣＲＴなどの表示装置の濃
度再現特性を考慮したダイナミックレンジ圧縮を行え
ば、画像信号上では所期の圧縮が行われているのに、表
示画像上では充分な圧縮効果を得ることができないとい
うことがなく、最終的な表示状態で最適となるダイナミ
ックレンジ圧縮を行わせることができる。
【００４０】上記実施例では、表示濃度の関係を階調変
換テーブルＬＵＴを用いて画像信号の関係に逆変換し、
表示濃度上で所期のダイナミックレンジ圧縮がなされる
ように、信号補正関数ｆ１（Ｓus）を設定したが、階調
変換テーブルＬＵＴの特性を見込んで簡易的に信号補正
関数ｆ１（Ｓus）を設定させる構成としても良い。
換テーブルＬＵＴを用いて画像信号の関係に逆変換し、
表示濃度上で所期のダイナミックレンジ圧縮がなされる
ように、信号補正関数ｆ１（Ｓus）を設定したが、階調
変換テーブルＬＵＴの特性を見込んで簡易的に信号補正
関数ｆ１（Ｓus）を設定させる構成としても良い。
【００４１】即ち、一般的な階調変換テーブルＬＵＴの
特性では、表示濃度上で１次関数で表される特性の補正
を実現させるためには、図３（ｃ）に示すように、信号
補正関数f1(Sus）としては非鮮鋭マスク信号Ｓusが高濃
度側或いは低濃度側に変化するに従って傾きを増すよう
な特性が要求される。
特性では、表示濃度上で１次関数で表される特性の補正
を実現させるためには、図３（ｃ）に示すように、信号
補正関数f1(Sus）としては非鮮鋭マスク信号Ｓusが高濃
度側或いは低濃度側に変化するに従って傾きを増すよう
な特性が要求される。
【００４２】従って、Ｓproc＝Ｓorg＋ｆ１（Ｓus）な
る演算を行ってダイナミックレンジの圧縮を行う構成に
おいて、予め前述のような要求を見込んで、図４に示す
ように、例えば高濃度側の圧縮を行う場合には、非鮮鋭
マスク信号Ｓusが増大するに従ってその減少割合を増す
ような補正関数ｆ１（Ｓus）を設定させれば良い。前記
図４に示した補正関数ｆ１（Ｓus）を式で示すと以下の
ようになる。
る演算を行ってダイナミックレンジの圧縮を行う構成に
おいて、予め前述のような要求を見込んで、図４に示す
ように、例えば高濃度側の圧縮を行う場合には、非鮮鋭
マスク信号Ｓusが増大するに従ってその減少割合を増す
ような補正関数ｆ１（Ｓus）を設定させれば良い。前記
図４に示した補正関数ｆ１（Ｓus）を式で示すと以下の
ようになる。
【００４３】ｆ１（Ｓus）＝β（Ｓus１−Ｓus）2…
（Ｓus≧Ｓus１）、ｆ１（Ｓus）＝０…（Ｓus＜Ｓus
１）
（Ｓus≧Ｓus１）、ｆ１（Ｓus）＝０…（Ｓus＜Ｓus
１）
【００４４】ここで、前記基準値Ｓus１は、オリジナル
画像のヒストグラム解析に基づいて決定することが好ま
しく、例えば関心領域内の最大信号値と最小信号値との
平均値とする。
画像のヒストグラム解析に基づいて決定することが好ま
しく、例えば関心領域内の最大信号値と最小信号値との
平均値とする。
【００４５】また、係数βは、関心領域内の最大信号値
をＳmaxとしたときに、β＝ｋ／（Ｓus１−Ｓmax）とし
て求めた。前記ｋは定数であり、0.5≦ｋ≦1.5程度とす
ることが好ましい。但し、前記βを演算するときに用い
る基準値は、最大値Ｓmax に限定されるものではない。
をＳmaxとしたときに、β＝ｋ／（Ｓus１−Ｓmax）とし
て求めた。前記ｋは定数であり、0.5≦ｋ≦1.5程度とす
ることが好ましい。但し、前記βを演算するときに用い
る基準値は、最大値Ｓmax に限定されるものではない。
【００４６】上記のように、実際に階調変換テーブルＬ
ＵＴから表示濃度を逆変換して、信号補正関数を設定す
るのではなく、階調変換テーブルＬＵＴの特性を見込ん
だ関数形に設定する構成とすれば、精度は低下するもの
の信号補正関数ｆ１（Ｓus）を簡便に設定できる。
ＵＴから表示濃度を逆変換して、信号補正関数を設定す
るのではなく、階調変換テーブルＬＵＴの特性を見込ん
だ関数形に設定する構成とすれば、精度は低下するもの
の信号補正関数ｆ１（Ｓus）を簡便に設定できる。
【００４７】尚、上記では、高濃度側の圧縮を階調変換
テーベルＬＵＴの特性に略適合させるべく設定した例を
示したが、低濃度側でも同様に設定することができ、Ｓ
proc＝Ｓorg＋ｆ１（Ｓus）なる演算を行ってダイナミ
ックレンジの圧縮を行う構成では、非鮮鋭マスク信号Ｓ
usの減少の応じて補正値ｆ１（Ｓus）の増大割合が増大
するような関数を設定させれば良い。
テーベルＬＵＴの特性に略適合させるべく設定した例を
示したが、低濃度側でも同様に設定することができ、Ｓ
proc＝Ｓorg＋ｆ１（Ｓus）なる演算を行ってダイナミ
ックレンジの圧縮を行う構成では、非鮮鋭マスク信号Ｓ
usの減少の応じて補正値ｆ１（Ｓus）の増大割合が増大
するような関数を設定させれば良い。
【００４８】即ち、前記非鮮鋭マスク信号Ｓusの変化に
対する前記補正値ｆ１（Ｓus）の変化割合の絶対値が、
オリジナル画像信号Ｓorgを増大補正する補正値ｆ１
（Ｓus）については、非鮮鋭マスク信号Ｓusの減少に応
じて徐々に増大変化し、オリジナル画像信号Ｓorgを減
少補正する補正値ｆ１（Ｓus）については、非鮮鋭マス
ク信号Ｓusの増大に応じて徐々に増大変化する関数を設
定すれば良い。
対する前記補正値ｆ１（Ｓus）の変化割合の絶対値が、
オリジナル画像信号Ｓorgを増大補正する補正値ｆ１
（Ｓus）については、非鮮鋭マスク信号Ｓusの減少に応
じて徐々に増大変化し、オリジナル画像信号Ｓorgを減
少補正する補正値ｆ１（Ｓus）については、非鮮鋭マス
ク信号Ｓusの増大に応じて徐々に増大変化する関数を設
定すれば良い。
【００４９】
【発明の効果】以上説明したように本発明にかかる放射
線画像のダイナミックレンジ圧縮方法によると、表示濃
度上で所期のダイナミックレンジ圧縮特性を設定し、こ
れから圧縮処理された画像信号を得るようにしたので、
最終的にダイナミックレンジが圧縮された放射線画像を
表示させたときに、所期の圧縮効果を確実に得られるよ
うになる。
線画像のダイナミックレンジ圧縮方法によると、表示濃
度上で所期のダイナミックレンジ圧縮特性を設定し、こ
れから圧縮処理された画像信号を得るようにしたので、
最終的にダイナミックレンジが圧縮された放射線画像を
表示させたときに、所期の圧縮効果を確実に得られるよ
うになる。
【００５０】また、表示における階調変換特性の傾向に
合わせて、補正値の関数を設定させるようにすること
で、簡便的に略表示画像上での要求に見合ったダイナミ
ックレンジの圧縮が行えるようになる。
合わせて、補正値の関数を設定させるようにすること
で、簡便的に略表示画像上での要求に見合ったダイナミ
ックレンジの圧縮が行えるようになる。
【図１】本発明が適用される放射線画像読取処理装置を
示すブロック図。
示すブロック図。
【図２】表示濃度に応じた圧縮処理の流れを示すフロー
チャート。
チャート。
【図３】表示濃度に応じた圧縮処理の特性を示す線図。
【図４】表示濃度特性を見込んだ補正値の関数形を示す
線図。
線図。
１…放射線発生源
３…記録読取装置
14…画像処理装置
15…画像メモリ
16…インタフェイス
17…プリンタ
─────────────────────────────────────────────────────
フロントページの続き
(51)Int.Cl.7 識別記号 ＦＩ
Ｈ０４Ｎ 7/18 Ｈ０４Ｎ 1/40 １０１Ｅ
(58)調査した分野(Int.Cl.7，ＤＢ名)
A61B 6/00 - 6/14
H04N 1/407
H04N 1/41
Claims (2)
- 【請求項１】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Ｓorgを処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Ｓprocを得る画像処
理を行った後に、該処理済み画像信号Ｓprocを表示装置
に出力する放射線画像のダイナミックレンジ圧縮方法で
あって、 各画素点に対応して非鮮鋭マスク信号Ｓusを求める一
方、前記オリジナル画像信号Ｓorgに対応する前記表示
装置における表示濃度をＤorg、前記処理済み画像信号
Ｓprocに対応する前記表示装置における表示濃度をＤpr
oc、前記非鮮鋭マスク信号Ｓusに対応する前記表示装置
における表示濃度をＤusとし、前記表示濃度Ｄusの関数
である濃度補正値をｆ２（Ｄus）としたときに、 Ｄproc＝Ｄorg ＋ｆ２（Ｄus） なる関係を満たすように前記処理済み画像信号Ｓprocを
得ることを特徴とする放射線画像のダイナミックレンジ
圧縮方法。 - 【請求項２】被写体を透過した放射線画像情報に基づく
オリジナル画像を表すオリジナル画像信号Ｓorgを処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Ｓprocを得る放射線
画像のダイナミックレンジ圧縮方法であって、各画素点
に対応して非鮮鋭マスク信号Ｓusを求め、該非鮮鋭マス
ク信号Ｓusの関数である補正値ｆ１（Ｓus）によってオ
リジナル画像信号Ｓorgを補正することで前記処理済み
画像信号Ｓprocを得る放射線画像のダイナミックレンジ
圧縮方法において、 前記非鮮鋭マスク信号Ｓusの変化に対する前記補正値ｆ
１（Ｓus）の変化割合の絶対値が、オリジナル画像信号
Ｓorgを増大補正する補正値ｆ１（Ｓus）については、
非鮮鋭マスク信号Ｓusの減少に応じて徐々に増大変化
し、オリジナル画像信号Ｓorgを減少補正する補正値ｆ
１（Ｓus）については、非鮮鋭マスク信号Ｓusの増大に
応じて徐々に増大変化することを特徴とする放射線画像
のダイナミックレンジ圧縮方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07649593A JP3425780B2 (ja) | 1993-03-30 | 1993-04-02 | 放射線画像のダイナミックレンジ圧縮方法 |
| US08/215,261 US5471987A (en) | 1993-03-30 | 1994-03-21 | Method of compressing a dynamic range for a radiation image |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7203193 | 1993-03-30 | ||
| JP5-72031 | 1993-03-30 | ||
| JP07649593A JP3425780B2 (ja) | 1993-03-30 | 1993-04-02 | 放射線画像のダイナミックレンジ圧縮方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003074985A Division JP3637049B2 (ja) | 1993-03-30 | 2003-03-19 | 放射線画像のダイナミックレンジ圧縮方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06339025A JPH06339025A (ja) | 1994-12-06 |
| JP3425780B2 true JP3425780B2 (ja) | 2003-07-14 |
Family
ID=26413157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP07649593A Expired - Fee Related JP3425780B2 (ja) | 1993-03-30 | 1993-04-02 | 放射線画像のダイナミックレンジ圧縮方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3425780B2 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1156451B1 (en) * | 1995-09-29 | 2004-06-02 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| US6226618B1 (en) * | 1998-08-13 | 2001-05-01 | International Business Machines Corporation | Electronic content delivery system |
| JP3937660B2 (ja) * | 1998-09-18 | 2007-06-27 | キヤノン株式会社 | 画像処理方法、及び、コンピュータ読み取り可能な記憶媒体 |
| US6721441B1 (en) * | 1999-12-30 | 2004-04-13 | General Electric Company | Extended dynamic range system for digital X-ray imaging detectors |
| FR2856170B1 (fr) * | 2003-06-10 | 2005-08-26 | Biospace Instr | Procede d'imagerie radiographique pour la reconstruction tridimensionnelle, dispositif et programme d'ordinateur pour mettre en oeuvre ce procede |
| WO2009104459A1 (ja) * | 2008-02-19 | 2009-08-27 | コニカミノルタエムジー株式会社 | 小規模施設用診断支援装置及びプログラム |
| JP5401945B2 (ja) * | 2008-11-25 | 2014-01-29 | カシオ電子工業株式会社 | 印刷装置、印刷方法及びそのプログラム |
-
1993
- 1993-04-02 JP JP07649593A patent/JP3425780B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06339025A (ja) | 1994-12-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3467285B2 (ja) | 放射線画像処理方法 | |
| US5454044A (en) | Apparatus for enhancing image data using a monotonously decreasing function | |
| US5319719A (en) | Processing apparatus for radiographic image signals | |
| EP0313042B1 (en) | Method and apparatus for radiation image processing and x-ray image processing | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP3425780B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP2509503B2 (ja) | 画像処理方法及び装置 | |
| JPH06292008A (ja) | 放射線画像のダイナミックレンジ圧縮処理装置 | |
| JP3455566B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP3434532B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP3637049B2 (ja) | 放射線画像のダイナミックレンジ圧縮方法 | |
| JP3563509B2 (ja) | 放射線画像読み取り方法及び放射線画像読み取り装置 | |
| JP2852794B2 (ja) | デジタル放射線画像信号の処理装置 | |
| JP3586274B2 (ja) | 放射線画像のダイナミックレンジ圧縮処理装置及び方法 | |
| JPH07248557A (ja) | 放射線画像の処理方法 | |
| JP4258092B2 (ja) | 画像処理装置および画像処理方法 | |
| JP2530223B2 (ja) | 放射線撮影の補正方法 | |
| JP2852799B2 (ja) | デジタル放射線画像信号の処理装置 | |
| JPH05192318A (ja) | Ｘ線画像作成装置 | |
| JP3165461B2 (ja) | 画像処理方法、画像読取方法および画像再生方法 | |
| JP3241488B2 (ja) | 画像処理方法 | |
| JPS6175966A (ja) | 放射線画像情報処理方法および装置 | |
| JPH02166882A (ja) | 画像情報読取装置 | |
| JPH0786925B2 (ja) | 放射線画像処理方法および装置 | |
| JPH0281276A (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: 20090509 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100509 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |