JPH06303430A - 画像のダイナミックレンジ圧縮処理方法 - Google Patents
画像のダイナミックレンジ圧縮処理方法Info
- Publication number
- JPH06303430A JPH06303430A JP5091276A JP9127693A JPH06303430A JP H06303430 A JPH06303430 A JP H06303430A JP 5091276 A JP5091276 A JP 5091276A JP 9127693 A JP9127693 A JP 9127693A JP H06303430 A JPH06303430 A JP H06303430A
- Authority
- JP
- Japan
- Prior art keywords
- sus
- image
- signal
- contrast
- function
- 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
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/407—Control or modification of tonal gradation or of extreme levels, e.g. background level
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
- G06T9/005—Statistical coding, e.g. Huffman, run length coding
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
(57)【要約】
【目的】 オリジナル画像の画像信号のダイナミックレ
ンジを圧縮して処理済画像信号を得るダイナミックレン
ジ圧縮処理方法において、オリジナル画像のコントラス
トが通常の画像のコントラストと異なる場合であって
も、ダイナミックレンジを適切に圧縮できるようにす
る。 【構成】 オリジナル画像信号Ｓorg のボケマスク信号
Ｓusを求め、このボケマスク信号Ｓusの値を変数とする
単調減少関数ｆ1 （Ｓus），ｆ2 （Ｓus）を用いて、 Ｓproc＝Ｓorg ＋ｆ1 ，2 （Ｓus） に従って、オリジナル画像信号Ｓorg よりダイナミック
レンジの狭い画像を担持する処理済画像信号Ｓprocを得
る。この際オリジナル画像信号Ｓorg のコントラストの
大小に応じて、関数ｆ1 （Ｓus），ｆ2 （Ｓus）の適用
信号範囲および／または値の大きさを関数ｆ1 ′ｆ2 ′
あるいは関数ｆ1 ″ｆ2 ″のように変化させる。
ンジを圧縮して処理済画像信号を得るダイナミックレン
ジ圧縮処理方法において、オリジナル画像のコントラス
トが通常の画像のコントラストと異なる場合であって
も、ダイナミックレンジを適切に圧縮できるようにす
る。 【構成】 オリジナル画像信号Ｓorg のボケマスク信号
Ｓusを求め、このボケマスク信号Ｓusの値を変数とする
単調減少関数ｆ1 （Ｓus），ｆ2 （Ｓus）を用いて、 Ｓproc＝Ｓorg ＋ｆ1 ，2 （Ｓus） に従って、オリジナル画像信号Ｓorg よりダイナミック
レンジの狭い画像を担持する処理済画像信号Ｓprocを得
る。この際オリジナル画像信号Ｓorg のコントラストの
大小に応じて、関数ｆ1 （Ｓus），ｆ2 （Ｓus）の適用
信号範囲および／または値の大きさを関数ｆ1 ′ｆ2 ′
あるいは関数ｆ1 ″ｆ2 ″のように変化させる。
Description
【０００１】
【産業上の利用分野】本発明は、オリジナル画像を表す
オリジナル画像信号を処理してオリジナル画像よりもダ
イナミックレンジの狭い画像を担持する処理済画像信号
を求める画像のダイナミックレンジ圧縮処理方法に関す
るものである。
オリジナル画像信号を処理してオリジナル画像よりもダ
イナミックレンジの狭い画像を担持する処理済画像信号
を求める画像のダイナミックレンジ圧縮処理方法に関す
るものである。
【０００２】
【従来の技術】記録された画像を読み取って画像信号を
得、この画像信号に適切な画像処理を施した後、画像を
再生表示することは種々の分野で行われている。たとえ
ば、後の画像処理に適合するように設計されたガンマ値
の低いＸ線フイルムを用いてＸ線画像を記録し、このＸ
線画像が記録されたフイルムからＸ線画像を読み取って
電気信号（画像信号）に変換し、この画像信号に画像処
理を施した後コピー写真等に可視像として再生すること
により、コントラスト，シャープネス，粒状性等の画質
性能の良好な再生画像を得ることが行われている（特公
昭61-5193 号公報参照）。
得、この画像信号に適切な画像処理を施した後、画像を
再生表示することは種々の分野で行われている。たとえ
ば、後の画像処理に適合するように設計されたガンマ値
の低いＸ線フイルムを用いてＸ線画像を記録し、このＸ
線画像が記録されたフイルムからＸ線画像を読み取って
電気信号（画像信号）に変換し、この画像信号に画像処
理を施した後コピー写真等に可視像として再生すること
により、コントラスト，シャープネス，粒状性等の画質
性能の良好な再生画像を得ることが行われている（特公
昭61-5193 号公報参照）。
【０００３】また、本願出願人により、放射線（Ｘ線，
α線，β線，γ線，電子線，紫外線等）を照射するとこ
の放射線エネルギーの一部が蓄積され、その後可視光等
の励起光を照射すると蓄積されたエネルギーに応じて輝
尽発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用し
て、人体等の被写体の放射線画像情報をシート状の蓄積
性蛍光体に一旦記録し、この蓄積性蛍光体シートをレー
ザー光等の励起光で走査して輝尽発光光を生ぜしめ、得
られた輝尽発光光を光電的に読み取って画像信号を得、
この画像信号に基づき被写体の放射線画像を写真感光材
料等の記録材料、ＣＲＴデイスプレイ装置等に可視画像
として出力させる放射線画像記録再生システムがすでに
提案されている（特開昭55-12429号，同56-11395号，同
55-163472 号，同56-104645 号，同55- 116340号等）。
α線，β線，γ線，電子線，紫外線等）を照射するとこ
の放射線エネルギーの一部が蓄積され、その後可視光等
の励起光を照射すると蓄積されたエネルギーに応じて輝
尽発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用し
て、人体等の被写体の放射線画像情報をシート状の蓄積
性蛍光体に一旦記録し、この蓄積性蛍光体シートをレー
ザー光等の励起光で走査して輝尽発光光を生ぜしめ、得
られた輝尽発光光を光電的に読み取って画像信号を得、
この画像信号に基づき被写体の放射線画像を写真感光材
料等の記録材料、ＣＲＴデイスプレイ装置等に可視画像
として出力させる放射線画像記録再生システムがすでに
提案されている（特開昭55-12429号，同56-11395号，同
55-163472 号，同56-104645 号，同55- 116340号等）。
【０００４】このシステムは、従来の銀塩写真を用いる
放射線写真システムと比較して極めて広い放射線露出域
にわたって画像を記録しうるという実用的な利点を有し
ている。すなわち、蓄積性蛍光体においては、放射線露
光量に対して蓄積後に励起によって輝尽発光する発光光
の光量が極めて広い範囲にわたって比例することが認め
られており、従って種々の撮影条件により放射線露光量
がかなり大幅に変動しても、蓄積性蛍光体シートより放
射される輝尽発光光の光量を読取ゲインを適当な値に設
定して光電変換手段により読み取って電気信号（画像信
号）に変換し、この画像信号を用いて写真感光材料等の
記録材料、ＣＲＴデイスプレイ装置等の表示装置に放射
線画像を可視画像として出力させることによって、放射
線露光量の変動に影響されない放射線画像を得ることが
できる。
放射線写真システムと比較して極めて広い放射線露出域
にわたって画像を記録しうるという実用的な利点を有し
ている。すなわち、蓄積性蛍光体においては、放射線露
光量に対して蓄積後に励起によって輝尽発光する発光光
の光量が極めて広い範囲にわたって比例することが認め
られており、従って種々の撮影条件により放射線露光量
がかなり大幅に変動しても、蓄積性蛍光体シートより放
射される輝尽発光光の光量を読取ゲインを適当な値に設
定して光電変換手段により読み取って電気信号（画像信
号）に変換し、この画像信号を用いて写真感光材料等の
記録材料、ＣＲＴデイスプレイ装置等の表示装置に放射
線画像を可視画像として出力させることによって、放射
線露光量の変動に影響されない放射線画像を得ることが
できる。
【０００５】上記のように画像信号を得てこの画像信号
に基づいて可視画像を得るシステムにおいて、観察対象
領域が低濃度域から高濃度域の広い濃度域に渡っている
場合や、どの濃度域を観察対象とするかが不明な場合等
に、オリジナル画像の高濃度域が観察に適した濃度とな
るように画像処理を施した画像とオリジナル画像の低濃
度域が観察に適した濃度となるように画像処理を施した
画像との双方を並べて再生することが行われている。
に基づいて可視画像を得るシステムにおいて、観察対象
領域が低濃度域から高濃度域の広い濃度域に渡っている
場合や、どの濃度域を観察対象とするかが不明な場合等
に、オリジナル画像の高濃度域が観察に適した濃度とな
るように画像処理を施した画像とオリジナル画像の低濃
度域が観察に適した濃度となるように画像処理を施した
画像との双方を並べて再生することが行われている。
【０００６】しかしながら、複数の画像を並べて再生す
ると各画像の寸法が小さくなり見ずらいという問題があ
る。
ると各画像の寸法が小さくなり見ずらいという問題があ
る。
【０００７】また、一つの画像で観察対象濃度域を広げ
るために、最高濃度と最低濃度との差すなわちダイナミ
ックレンジを狭めるように高濃度域もしくは低濃度域も
しくは画像全体のコントラストを下げることも行われて
いる。
るために、最高濃度と最低濃度との差すなわちダイナミ
ックレンジを狭めるように高濃度域もしくは低濃度域も
しくは画像全体のコントラストを下げることも行われて
いる。
【０００８】しかし、コントラストを下げると、そのコ
ントラストを下げた領域内の微細構造も見にくくなると
いう問題が生じる。
ントラストを下げた領域内の微細構造も見にくくなると
いう問題が生じる。
【０００９】以下この問題点について詳述する。
【００１０】図15は、オリジナル画像上のある方向（ｘ
方向）に沿うオリジナル画像信号Ｓorg の値の一例をプ
ロットした図である。オリジナル画像はｘ方向に沿って
全体としては階段状の濃度（画像信号Ｓorg ）を有し、
かつ各階段内に比較的高空間周波数の微細構造を有して
いる。
方向）に沿うオリジナル画像信号Ｓorg の値の一例をプ
ロットした図である。オリジナル画像はｘ方向に沿って
全体としては階段状の濃度（画像信号Ｓorg ）を有し、
かつ各階段内に比較的高空間周波数の微細構造を有して
いる。
【００１１】図16(a) ，(b) は高濃度域のコントラスト
を下げることによりダイナミックレンジを圧縮した場合
の、それぞれ画像信号の変換グラフの一例を示した図、
および画像上のｘ方向に沿う変換後の画像信号Ｓorg ′
の値を示した図である。
を下げることによりダイナミックレンジを圧縮した場合
の、それぞれ画像信号の変換グラフの一例を示した図、
および画像上のｘ方向に沿う変換後の画像信号Ｓorg ′
の値を示した図である。
【００１２】図15に示した値を有する画像信号Ｓorg
が、図16(a) に示すグラフＡに沿って画像信号Ｓorg ′
に変換されると、図16(b) に示すように高濃度域の濃度
が低下するが、これとともに高濃度域の階段内の微細構
造もコントラストを失う結果となり、このことは観察対
象としている、全体として高濃度の領域内の微細構造が
非常に観察しにくいものとなってしまうことを意味して
いる。
が、図16(a) に示すグラフＡに沿って画像信号Ｓorg ′
に変換されると、図16(b) に示すように高濃度域の濃度
が低下するが、これとともに高濃度域の階段内の微細構
造もコントラストを失う結果となり、このことは観察対
象としている、全体として高濃度の領域内の微細構造が
非常に観察しにくいものとなってしまうことを意味して
いる。
【００１３】図17(a) ，(b) は低濃度域のコントラスト
を下げることによりダイナミックレンジを圧縮した場合
の、それぞれ画像信号の変換グラフの一例を示した図、
および画像上のｘ方向に沿う変換後の画像信号Ｓorg ′
の値を示した図である。
を下げることによりダイナミックレンジを圧縮した場合
の、それぞれ画像信号の変換グラフの一例を示した図、
および画像上のｘ方向に沿う変換後の画像信号Ｓorg ′
の値を示した図である。
【００１４】図15に示した値を有する画像信号Ｓorg
が、図17(a) に示すグラフＢに沿って画像信号Ｓorg ′
に変換されると、図17(b) に示すように今度は低濃度域
内の微細構造が非常に観察しにくいものとなってしまう
結果となる。
が、図17(a) に示すグラフＢに沿って画像信号Ｓorg ′
に変換されると、図17(b) に示すように今度は低濃度域
内の微細構造が非常に観察しにくいものとなってしまう
結果となる。
【００１５】そこで、画像のダイナミックレンジを圧縮
するための方法が提案されている（日本放射線技術学会
雑誌 第45巻第８号 1989年８月 1030頁 阿南充洋ほ
か）。この方法は、人体の胸部Ｘ線画像について、 (1) ボケマスク信号Ｓusを求める (2) ボケマスク信号Ｓusの値を２倍し、最大値（1023）
以上となった場合は最大値（1023）でクリップする。
するための方法が提案されている（日本放射線技術学会
雑誌 第45巻第８号 1989年８月 1030頁 阿南充洋ほ
か）。この方法は、人体の胸部Ｘ線画像について、 (1) ボケマスク信号Ｓusを求める (2) ボケマスク信号Ｓusの値を２倍し、最大値（1023）
以上となった場合は最大値（1023）でクリップする。
【００１６】 ａ＝ ２×Ｓus （２×Ｓus＜1023） 1023 （２×Ｓus≧1023） (3) 反転画像を作成する。
【００１７】ｂ＝1023−ａ (4) オリジナル画像に反転画像×係数αを加える。
【００１８】ｃ＝Ｓorg ＋α・ｂ （α＝0.3 ） のステップにより、画像のダイナミックレンジを圧縮す
る手法が記載されている。この手法は本発明と同様に画
像のダイナミックレンジを圧縮するとともに各領域内の
細部構造のコントラストは維持されるという効果を有し
ている。
る手法が記載されている。この手法は本発明と同様に画
像のダイナミックレンジを圧縮するとともに各領域内の
細部構造のコントラストは維持されるという効果を有し
ている。
【００１９】しかしこの手法は、低濃度域のダイナミッ
クレンジのみを圧縮するものであって、高濃度域のダイ
ナミックレンジを圧縮することはできず、例えば前述し
た四肢骨画像等の観察には適さないという欠点がある。
また、この手法は処理された画像に偽輪郭を生じる可能
性があるため、ダイナミックレンジは圧縮されたものの
今度はこの偽輪郭のために画像が見にくくなるおそれが
あるという欠点もある。 そこで本願出願人により、オ
リジナル画像を表すオリジナル画像信号Ｓorgを平均化
することによりボケマスク信号Ｓusを求め、このボケマ
スク信号Ｓusの値が増大するにつれて単調減少する関数
をｆ1 （Ｓus）としたとき、式 Ｓproc＝Ｓorg ＋ｆ1 （Ｓus） …(1) に従ってオリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済画像信号Ｓprocを求める画像の
ダイナミックレンジ圧縮処理方法が提案されている（特
開平3-222577号公報）。この方法によれば、低濃度領域
のみならず高濃度領域のダイナミックレンジをも圧縮で
き、また、関数ｆ1 （Ｓus）の微係数を連続させれば、
偽輪郭が生じる問題もなくなり、濃度域を圧縮すること
により画像内の適性観察領域を広げるとともに、各領域
内の微細構造の観察適性をも確保することができる。
クレンジのみを圧縮するものであって、高濃度域のダイ
ナミックレンジを圧縮することはできず、例えば前述し
た四肢骨画像等の観察には適さないという欠点がある。
また、この手法は処理された画像に偽輪郭を生じる可能
性があるため、ダイナミックレンジは圧縮されたものの
今度はこの偽輪郭のために画像が見にくくなるおそれが
あるという欠点もある。 そこで本願出願人により、オ
リジナル画像を表すオリジナル画像信号Ｓorgを平均化
することによりボケマスク信号Ｓusを求め、このボケマ
スク信号Ｓusの値が増大するにつれて単調減少する関数
をｆ1 （Ｓus）としたとき、式 Ｓproc＝Ｓorg ＋ｆ1 （Ｓus） …(1) に従ってオリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済画像信号Ｓprocを求める画像の
ダイナミックレンジ圧縮処理方法が提案されている（特
開平3-222577号公報）。この方法によれば、低濃度領域
のみならず高濃度領域のダイナミックレンジをも圧縮で
き、また、関数ｆ1 （Ｓus）の微係数を連続させれば、
偽輪郭が生じる問題もなくなり、濃度域を圧縮すること
により画像内の適性観察領域を広げるとともに、各領域
内の微細構造の観察適性をも確保することができる。
【００２０】
【発明が解決しようとする課題】しかしながら、上述し
た特開平3-222577号公報に開示された方法においては、
処理済画像信号をＳorg を求めるための関数ｆ1 （Ｓu
s）はいかなる画像信号であっても固定的であり、画像
のコントラストを考慮しない方法であるため、通常処理
されるの画像よりもダイナミックレンジが広く、コント
ラストが大きい画像の画像信号については、適切なダイ
ナミックレンジの圧縮を行うことができず、この画像信
号の高濃度領域あるいは低濃度領域もしくは画像全体の
コントラストを適切に下げることができなかった。
た特開平3-222577号公報に開示された方法においては、
処理済画像信号をＳorg を求めるための関数ｆ1 （Ｓu
s）はいかなる画像信号であっても固定的であり、画像
のコントラストを考慮しない方法であるため、通常処理
されるの画像よりもダイナミックレンジが広く、コント
ラストが大きい画像の画像信号については、適切なダイ
ナミックレンジの圧縮を行うことができず、この画像信
号の高濃度領域あるいは低濃度領域もしくは画像全体の
コントラストを適切に下げることができなかった。
【００２１】本発明は上記事情に鑑み、コントラストが
通常と異なる画像であっても、適切にダイナミックレン
ジの圧縮を行うことができる画像のダイナミックレンジ
圧縮処理方法を提供することを目的とするものである。
通常と異なる画像であっても、適切にダイナミックレン
ジの圧縮を行うことができる画像のダイナミックレンジ
圧縮処理方法を提供することを目的とするものである。
【００２２】
【課題を解決するための手段】本発明よる第１の画像の
ダイナミックレンジ圧縮処理方法は、オリジナル画像を
表すオリジナル画像信号Ｓorg の各画素点に対応して該
各画素点の周囲の所定範囲内の前記オリジナル画像信号
Ｓorg を平均化することによりボケマスク信号Ｓusを求
め、該ボケマスク信号Ｓusの値が増大するにつれて単調
減少する関数をｆ1 （Ｓus）としたとき、式 Ｓproc＝Ｓorg ＋ｆ1 （Ｓus） に従って前記オリジナル画像信号を処理して前記オリジ
ナル画像よりもダイナミックレンジの狭い画像を担持す
る処理済画像信号Ｓprocを求める画像のダイナミックレ
ンジ圧縮処理方法において、前記オリジナル画像信号に
基づいて該オリジナル画像信号のコントラストを求め、
該コントラストに応じて前記関数ｆ1 （Ｓus）の適用信
号範囲および／または値の大きさを決定することを特徴
とするものである。
ダイナミックレンジ圧縮処理方法は、オリジナル画像を
表すオリジナル画像信号Ｓorg の各画素点に対応して該
各画素点の周囲の所定範囲内の前記オリジナル画像信号
Ｓorg を平均化することによりボケマスク信号Ｓusを求
め、該ボケマスク信号Ｓusの値が増大するにつれて単調
減少する関数をｆ1 （Ｓus）としたとき、式 Ｓproc＝Ｓorg ＋ｆ1 （Ｓus） に従って前記オリジナル画像信号を処理して前記オリジ
ナル画像よりもダイナミックレンジの狭い画像を担持す
る処理済画像信号Ｓprocを求める画像のダイナミックレ
ンジ圧縮処理方法において、前記オリジナル画像信号に
基づいて該オリジナル画像信号のコントラストを求め、
該コントラストに応じて前記関数ｆ1 （Ｓus）の適用信
号範囲および／または値の大きさを決定することを特徴
とするものである。
【００２３】ここで、関数ｆ1 （Ｓus）の適用信号範囲
とは、画像信号のうち関数ｆ1 （Ｓus）により処理をう
ける濃度領域の範囲をいい、値の大きさとは、信号を強
調する強調度の大きさをいう。
とは、画像信号のうち関数ｆ1 （Ｓus）により処理をう
ける濃度領域の範囲をいい、値の大きさとは、信号を強
調する強調度の大きさをいう。
【００２４】また、本発明による第２の画像のダイナミ
ックレンジ圧縮処理方法は、前記単調減少する関数ｆ1
（Ｓus）が、少なくとも一部に曲線を含むことを特徴と
するものである。
ックレンジ圧縮処理方法は、前記単調減少する関数ｆ1
（Ｓus）が、少なくとも一部に曲線を含むことを特徴と
するものである。
【００２５】ここで、単調減少するとは、ボケマスク信
号Ｓusの値の増大にしたがって、関数ｆ1 （Ｓus）が必
ず減少するものである必要はなく、部分的にはボケマス
ク信号Ｓusが変化しても関数ｆ1 （Ｓus）が変化しない
領域が存在してもよい。
号Ｓusの値の増大にしたがって、関数ｆ1 （Ｓus）が必
ず減少するものである必要はなく、部分的にはボケマス
ク信号Ｓusが変化しても関数ｆ1 （Ｓus）が変化しない
領域が存在してもよい。
【００２６】
【作用】本発明による画像のダイナミックレンジ圧縮処
理方法は、オリジナル画像信号のコントラストの大きさ
に応じてダイナミックレンジの圧縮を行うための関数ｆ
1 （Ｓus）の適用信号範囲および／または値の大きさを
決定するようにしたため、通常処理される画像よりもコ
ントラストが広い画像の場合は関数ｆ1 （Ｓus）の適用
信号範囲を広くするおよび／または値の大きさを大きく
することにより、またコントラストが狭い画像の場合は
関数ｆ1 （Ｓus）の適用信号範囲を狭くするおよび／ま
たは値の大きさを小さくすることにより、画像の濃度域
は適性観察領域に圧縮されることとなる。したがって、
コントラストの大小に応じて、観察に的した処理済画像
を得ることができる。
理方法は、オリジナル画像信号のコントラストの大きさ
に応じてダイナミックレンジの圧縮を行うための関数ｆ
1 （Ｓus）の適用信号範囲および／または値の大きさを
決定するようにしたため、通常処理される画像よりもコ
ントラストが広い画像の場合は関数ｆ1 （Ｓus）の適用
信号範囲を広くするおよび／または値の大きさを大きく
することにより、またコントラストが狭い画像の場合は
関数ｆ1 （Ｓus）の適用信号範囲を狭くするおよび／ま
たは値の大きさを小さくすることにより、画像の濃度域
は適性観察領域に圧縮されることとなる。したがって、
コントラストの大小に応じて、観察に的した処理済画像
を得ることができる。
【００２７】また、前述した関数ｆ1 （Ｓus）の少なく
とも一部に曲線を含ませるようにすれば、処理済画像信
号Ｓprocが担持する画像に儀輪郭が生じることが防止さ
れ、より観察しやすい画像となる。
とも一部に曲線を含ませるようにすれば、処理済画像信
号Ｓprocが担持する画像に儀輪郭が生じることが防止さ
れ、より観察しやすい画像となる。
【００２８】
【実施例】以下、図面を参照して、本発明の実施例につ
いて説明する。なお、ここでは前述した蓄積性蛍光体シ
ートを用いた例について説明する。
いて説明する。なお、ここでは前述した蓄積性蛍光体シ
ートを用いた例について説明する。
【００２９】図１は、Ｘ線撮影装置の一例の概略図であ
る。
る。
【００３０】Ｘ線撮影装置10のＸ線源11からＸ線12が被
写体13に向けて照射され、被写体13を透過したＸ線が蓄
積性蛍光体シート14に照射されることにより、被写体13
の透過Ｘ線画像が蓄積性蛍光体シート14に蓄積記録され
る。
写体13に向けて照射され、被写体13を透過したＸ線が蓄
積性蛍光体シート14に照射されることにより、被写体13
の透過Ｘ線画像が蓄積性蛍光体シート14に蓄積記録され
る。
【００３１】図２は、Ｘ線画像読取装置の一例を表した
斜視図である。
斜視図である。
【００３２】図２に示したＸ線撮影装置で撮影が行わ
れ、Ｘ線画像が記録された蓄積性蛍光体シート14が読取
部20の所定位置にセットされる。
れ、Ｘ線画像が記録された蓄積性蛍光体シート14が読取
部20の所定位置にセットされる。
【００３３】蓄積性蛍光体シート14が読取部20の所定位
置にセットされると、このシート14はモータ21により駆
動されるエンドレスベルト22により、矢印Ｙ方向に搬送
（副走査）される。一方、レーザー光源23から発せられ
た光ビーム24はモータ25により駆動され矢印方向に高速
回転する回転多面鏡26によって反射偏向され、ｆθレン
ズ等の集束レンズ27を通過した後、ミラー28により光路
を変えて前記シート14に入射し副走査の方向（矢印Ｙ方
向）と略垂直な矢印Ｘ方向に主走査する。シート14の励
起光24が照射された箇所からは、蓄積記録されているＸ
線画像情報に応じた光量の輝尽発光光29が発散され、こ
の輝尽発光光29は光ガイド30によって導かれ、フォトマ
ルチプライヤ（光電子増倍管）31によって光電的に検出
される。上記光ガイド30はアクリル板等の導光性材料を
成形して作られたものであり、直線状をなす入射端面30
a が蓄積性蛍光体シート14上の主走査線に沿って延びる
ように配され、円環状に形成された射出端面30b にフォ
トマルチプライヤ31の受光面が結合されている。入射端
面30a から光ガイド30内に入射した輝尽発光光29は、該
光ガイド30の内部を全反射を繰り返して進み、射出端面
30b から射出してフォトマルチプライヤ31に受光され、
Ｘ線画像を表す輝尽発光光29がフォトマルチプライヤ31
によって電気信号に変換される。
置にセットされると、このシート14はモータ21により駆
動されるエンドレスベルト22により、矢印Ｙ方向に搬送
（副走査）される。一方、レーザー光源23から発せられ
た光ビーム24はモータ25により駆動され矢印方向に高速
回転する回転多面鏡26によって反射偏向され、ｆθレン
ズ等の集束レンズ27を通過した後、ミラー28により光路
を変えて前記シート14に入射し副走査の方向（矢印Ｙ方
向）と略垂直な矢印Ｘ方向に主走査する。シート14の励
起光24が照射された箇所からは、蓄積記録されているＸ
線画像情報に応じた光量の輝尽発光光29が発散され、こ
の輝尽発光光29は光ガイド30によって導かれ、フォトマ
ルチプライヤ（光電子増倍管）31によって光電的に検出
される。上記光ガイド30はアクリル板等の導光性材料を
成形して作られたものであり、直線状をなす入射端面30
a が蓄積性蛍光体シート14上の主走査線に沿って延びる
ように配され、円環状に形成された射出端面30b にフォ
トマルチプライヤ31の受光面が結合されている。入射端
面30a から光ガイド30内に入射した輝尽発光光29は、該
光ガイド30の内部を全反射を繰り返して進み、射出端面
30b から射出してフォトマルチプライヤ31に受光され、
Ｘ線画像を表す輝尽発光光29がフォトマルチプライヤ31
によって電気信号に変換される。
【００３４】フォトマルチプライヤ31から出力されたア
ナログ出力信号Ｓ０は対数増幅器32で対数的に増幅さ
れ、Ａ／Ｄ変換器33でディジタル化され、これによりオ
リジナル画像信号Ｓorg が得られ、画像処理表示部40に
入力される。この画像処理表示部40は、可視画像を再生
表示するためのＣＲＴディスプレイ41、ＣＰＵ，内部メ
モリ，インターフェイス等が内蔵された本体部42、フロ
ッピィディスクが装填され駆動されるフロッピィディス
クドライブ部43、およびこのＸ線画像読取装置に必要な
情報を入力するためのキーボード44から構成されてい
る。
ナログ出力信号Ｓ０は対数増幅器32で対数的に増幅さ
れ、Ａ／Ｄ変換器33でディジタル化され、これによりオ
リジナル画像信号Ｓorg が得られ、画像処理表示部40に
入力される。この画像処理表示部40は、可視画像を再生
表示するためのＣＲＴディスプレイ41、ＣＰＵ，内部メ
モリ，インターフェイス等が内蔵された本体部42、フロ
ッピィディスクが装填され駆動されるフロッピィディス
クドライブ部43、およびこのＸ線画像読取装置に必要な
情報を入力するためのキーボード44から構成されてい
る。
【００３５】この画像処理表示部40に被写体13のオリジ
ナルＸ線画像を担持するオリジナル画像信号Ｓorg が入
力されると、このオリジナル画像信号Ｓorg を以下のよ
うに変換することにより画像のダイナミックレンジが圧
縮される。
ナルＸ線画像を担持するオリジナル画像信号Ｓorg が入
力されると、このオリジナル画像信号Ｓorg を以下のよ
うに変換することにより画像のダイナミックレンジが圧
縮される。
【００３６】図３は、オリジナル画像上の各画素点と該
各画素点に対応するオリジナル画像信号Ｓorg を表した
図である。図にＳij等で示した記号が対応する各画素点
のオリジナル画像信号Ｓorg を表している。ここで中央
の白丸で示した画素のボケマスク信号Ｓusijが
各画素点に対応するオリジナル画像信号Ｓorg を表した
図である。図にＳij等で示した記号が対応する各画素点
のオリジナル画像信号Ｓorg を表している。ここで中央
の白丸で示した画素のボケマスク信号Ｓusijが
【００３７】
【数１】
【００３８】の演算により求められ、この演算を各画素
について行うことにより画像全体のボケマスク信号Ｓus
が求められる。なお、ｍ，ｎは、オリジナル画像信号Ｓ
org を得る際のサンプリング間隔やオリジナルＸ線画像
の性質等により適宜選択される値である。
について行うことにより画像全体のボケマスク信号Ｓus
が求められる。なお、ｍ，ｎは、オリジナル画像信号Ｓ
org を得る際のサンプリング間隔やオリジナルＸ線画像
の性質等により適宜選択される値である。
【００３９】これと同時に、オリジナル画像信号Ｓorg
のコントラストが求められる。図４は、オリジナル画像
信号Ｓorg のコントラストＧp をコントラスト変換処理
前の濃度とコントラスト変換処理後の濃度との関係にお
いて表す図である。図４に示すコントラストＧp によ
り、コントラスト変換処理前の画像は所定の濃度領域を
有するコントラスト変換処理後画像として再生される。
なお、コントラストＧpは高濃度側の点Ｏを固定し所定
範囲内に（Ｇpmin≦Ｇp ≦Ｇpmax）にあるように制限さ
れる。
のコントラストが求められる。図４は、オリジナル画像
信号Ｓorg のコントラストＧp をコントラスト変換処理
前の濃度とコントラスト変換処理後の濃度との関係にお
いて表す図である。図４に示すコントラストＧp によ
り、コントラスト変換処理前の画像は所定の濃度領域を
有するコントラスト変換処理後画像として再生される。
なお、コントラストＧpは高濃度側の点Ｏを固定し所定
範囲内に（Ｇpmin≦Ｇp ≦Ｇpmax）にあるように制限さ
れる。
【００４０】図５(a) ，(b) は、ボケマスク信号Ｓusを
変数とした単調減少関数の一例を表した図である。
変数とした単調減少関数の一例を表した図である。
【００４１】ボケマスク信号Ｓusは最高値1023を有し、
図５(a) の関数ｆ1 （Ｓus）は、ボケマスク信号Ｓusの
値が小さい領域で変化し、ボケマスク信号Ｓusの値がｄ
点より大きい領域ではｆ（Ｓus）＝０となる関数であ
る。図５(b) の関数ｆ2 （Ｓus）は途中のｅ点までは
零、それ以上は直線的に傾斜した関数形を有している。
この関数ｆ1 （Ｓus）および／またはｆ2 （Ｓus）を用
いて各画素点（ｉ，ｊ）について、 Ｓprocij＝Ｓij＋ｆk （Ｓusij）（ｋ＝１，２） …(3) の演算を行い、画像全体について処理済画像信号Ｓproc
が求められる。
図５(a) の関数ｆ1 （Ｓus）は、ボケマスク信号Ｓusの
値が小さい領域で変化し、ボケマスク信号Ｓusの値がｄ
点より大きい領域ではｆ（Ｓus）＝０となる関数であ
る。図５(b) の関数ｆ2 （Ｓus）は途中のｅ点までは
零、それ以上は直線的に傾斜した関数形を有している。
この関数ｆ1 （Ｓus）および／またはｆ2 （Ｓus）を用
いて各画素点（ｉ，ｊ）について、 Ｓprocij＝Ｓij＋ｆk （Ｓusij）（ｋ＝１，２） …(3) の演算を行い、画像全体について処理済画像信号Ｓproc
が求められる。
【００４２】ここで、図４に示すコントラストＧp は、
所定の範囲内（Ｇpmin≦Ｇp ≦Ｇpmax）にあるように制
限されるものであるが、画像によってはコントラストＧ
p がこの範囲外となるような場合がある。この場合は、
算出されたコントラストＧpに応じて図５に示す関数ｆ
1 （Ｓus）および／またはｆ2 （Ｓus）の適用信号範囲
および／または値の大きさを変化させて、ダイナミック
レンジ圧縮処理を行う。
所定の範囲内（Ｇpmin≦Ｇp ≦Ｇpmax）にあるように制
限されるものであるが、画像によってはコントラストＧ
p がこの範囲外となるような場合がある。この場合は、
算出されたコントラストＧpに応じて図５に示す関数ｆ
1 （Ｓus）および／またはｆ2 （Ｓus）の適用信号範囲
および／または値の大きさを変化させて、ダイナミック
レンジ圧縮処理を行う。
【００４３】例えば、図６に示すように求められたコン
トラストがコントラストＧp1に示すようにＧp1＞Ｇpmax
となった場合は、全体のコントラストは圧縮気味とな
り、またコントラスト変換処理後の濃度は高くなるた
め、関数ｆ1 （Ｓus）および／または関数ｆ2 （Ｓus）
を図７の関数ｆ1 ′，ｆ2 ′に示すように関数ｆ1 （Ｓ
us）の適用信号範囲を広くかつ値を小さくし、ｆ2 （Ｓ
us）の適用信号範囲を狭くすることにより、圧縮の程度
を小さくしてダイナミックレンジ圧縮処理を行う。
トラストがコントラストＧp1に示すようにＧp1＞Ｇpmax
となった場合は、全体のコントラストは圧縮気味とな
り、またコントラスト変換処理後の濃度は高くなるた
め、関数ｆ1 （Ｓus）および／または関数ｆ2 （Ｓus）
を図７の関数ｆ1 ′，ｆ2 ′に示すように関数ｆ1 （Ｓ
us）の適用信号範囲を広くかつ値を小さくし、ｆ2 （Ｓ
us）の適用信号範囲を狭くすることにより、圧縮の程度
を小さくしてダイナミックレンジ圧縮処理を行う。
【００４４】また、求められたコントラストがコントラ
ストＧp2に示すようにＧp2＜Ｇpminとなった場合は、全
体のコントラストは伸長気味となり、またコントラスト
変換処理後の濃度は低くなるため、関数ｆ1 （Ｓus）お
よび／または関数ｆ2 （Ｓus）を図７の関数ｆ1 ″
ｆ2 ″に示すように関数ｆ1 （Ｓus）の適用信号範囲を
狭くかつ値を大きくし、関数ｆ2 （Ｓus）の適用信号範
囲を広くすることにより、圧縮の程度を大きくしてダイ
ナミックレンジ圧縮処理を行う。
ストＧp2に示すようにＧp2＜Ｇpminとなった場合は、全
体のコントラストは伸長気味となり、またコントラスト
変換処理後の濃度は低くなるため、関数ｆ1 （Ｓus）お
よび／または関数ｆ2 （Ｓus）を図７の関数ｆ1 ″
ｆ2 ″に示すように関数ｆ1 （Ｓus）の適用信号範囲を
狭くかつ値を大きくし、関数ｆ2 （Ｓus）の適用信号範
囲を広くすることにより、圧縮の程度を大きくしてダイ
ナミックレンジ圧縮処理を行う。
【００４５】このようにオリジナル画像信号Ｓorg から
得られたコントラストＧp の大きさに応じて、関数ｆ1
（Ｓus），ｆ2 （Ｓus）の適用信号範囲および／または
値の大きさを変化させることにより、オリジナル画像の
最高濃度と最低濃度との差に応じたダイナミックレンジ
の圧縮処理を行うことができる。
得られたコントラストＧp の大きさに応じて、関数ｆ1
（Ｓus），ｆ2 （Ｓus）の適用信号範囲および／または
値の大きさを変化させることにより、オリジナル画像の
最高濃度と最低濃度との差に応じたダイナミックレンジ
の圧縮処理を行うことができる。
【００４６】一方、コントラストＧp が図８に示すよう
に低濃度側の点Ｏ2 を固定し、所定範囲内（Ｇpmin≦Ｇ
p ≦Ｇpmax）にあるように制限される場合において、コ
ントラストがコントラストＧp3に示すようにＧp3＞Ｇpm
axとなった場合は、全体のコントラストは圧縮気味とな
り、またコントラスト変換処理後の濃度は低くなるた
め、関数ｆ1 （Ｓus）および／またはｆ2 （Ｓus）を図
９の関数ｆ1 ′，ｆ2 ′に示すように、関数ｆ1 （Ｓu
s）の適用信号範囲を狭くし、関数ｆ2 （Ｓus）の適用
信号範囲を広くかつ値を小さくして、圧縮の程度を低く
してダイナミックレンジ圧縮処理を行う。
に低濃度側の点Ｏ2 を固定し、所定範囲内（Ｇpmin≦Ｇ
p ≦Ｇpmax）にあるように制限される場合において、コ
ントラストがコントラストＧp3に示すようにＧp3＞Ｇpm
axとなった場合は、全体のコントラストは圧縮気味とな
り、またコントラスト変換処理後の濃度は低くなるた
め、関数ｆ1 （Ｓus）および／またはｆ2 （Ｓus）を図
９の関数ｆ1 ′，ｆ2 ′に示すように、関数ｆ1 （Ｓu
s）の適用信号範囲を狭くし、関数ｆ2 （Ｓus）の適用
信号範囲を広くかつ値を小さくして、圧縮の程度を低く
してダイナミックレンジ圧縮処理を行う。
【００４７】また、求められたコントラストがコントラ
ストＧp4に示すように、Ｇp4＜Ｇpminとなった場合は、
全体のコントラストは伸長気味となり、またコントラス
ト変換処理後の濃度は高くなるため、関数ｆ1 （Ｓus）
および／または関数ｆ2 （Ｓus）を図９の関数ｆ1 ″，
ｆ2 ″に示すように、関数ｆ1 （Ｓus）の適用信号範囲
を広くし、関数ｆ2 （Ｓus）の適用信号範囲を狭くかつ
値を大きくして、圧縮の程度を大きくしてダイナミック
レンジ圧縮処理を行えばよい。
ストＧp4に示すように、Ｇp4＜Ｇpminとなった場合は、
全体のコントラストは伸長気味となり、またコントラス
ト変換処理後の濃度は高くなるため、関数ｆ1 （Ｓus）
および／または関数ｆ2 （Ｓus）を図９の関数ｆ1 ″，
ｆ2 ″に示すように、関数ｆ1 （Ｓus）の適用信号範囲
を広くし、関数ｆ2 （Ｓus）の適用信号範囲を狭くかつ
値を大きくして、圧縮の程度を大きくしてダイナミック
レンジ圧縮処理を行えばよい。
【００４８】さらに、コントラストＧp が図10に示すよ
うに濃度の略中間地点Ｏ3 を固定し、所定範囲内（Ｇpm
in≦Ｇp ≦Ｇpmax）にあるように制限される場合は、低
濃度側の濃度を上げるための関数ｆ1 （Ｓus）は図７
(a) に示すものを用い、高濃度側の濃度を下げるための
関数ｆ2 （Ｓus）は図９(b) に示すものを用い、コント
ラストＧp の大小に応じて各関数ｆ1 （Ｓus），ｆ
2 （Ｓus）を関数ｆ1 ′，ｆ1 ″，ｆ2 ′，ｆ2 ″に示
すように変化させればよい。
うに濃度の略中間地点Ｏ3 を固定し、所定範囲内（Ｇpm
in≦Ｇp ≦Ｇpmax）にあるように制限される場合は、低
濃度側の濃度を上げるための関数ｆ1 （Ｓus）は図７
(a) に示すものを用い、高濃度側の濃度を下げるための
関数ｆ2 （Ｓus）は図９(b) に示すものを用い、コント
ラストＧp の大小に応じて各関数ｆ1 （Ｓus），ｆ
2 （Ｓus）を関数ｆ1 ′，ｆ1 ″，ｆ2 ′，ｆ2 ″に示
すように変化させればよい。
【００４９】なお、上述した実施例においては、ダイナ
ミックレンジ圧縮処理を行うための関数として、関数ｆ
1 （Ｓus），ｆ2 （Ｓus）の２つの関数を用いている
が、いずれか一方の関数のみを用いるようにしてもよい
ものである。
ミックレンジ圧縮処理を行うための関数として、関数ｆ
1 （Ｓus），ｆ2 （Ｓus）の２つの関数を用いている
が、いずれか一方の関数のみを用いるようにしてもよい
ものである。
【００５０】また、上述した実施例においてダイナミッ
クレンジ圧縮処理に用いる関数は、上述した関数ｆ
1 （Ｓus），ｆ2 （Ｓus）に限定されるものではなく、
例えば図11に示す関数ｆ3 （Ｓus）のように一部に曲線
を含むものであってもよく、図12に示す関数ｆ4 （Ｓu
s）のようにｆ4 （Ｓus）＝０とならない部分を含むも
のであってもよい。さらに図13に示す関数ｆ4 （Ｓus）
のように低濃度側と高濃度側とで圧縮の程度が高なる関
数を用いるようにしてもよい。さらに、図14に示す関数
ｆ6 （Ｓus）のように、ｆ6 （Ｓus）＝０となる点ｄに
おいて折れ線とならずに微係数が連続する関数としても
よい。点ｄが折れ線となる場合、この点ｄに相当する濃
度領域に、オリジナル画像には何ら特別の輪郭は存在し
ないにも拘らず、処理済画像信号Ｓprocに基づいて再生
表示された可視画像に偽輪郭が生じる場合がある。そこ
で、点ｄの部分を微係数が連続する関数とすることによ
り偽輪郭が生じることも防止される。
クレンジ圧縮処理に用いる関数は、上述した関数ｆ
1 （Ｓus），ｆ2 （Ｓus）に限定されるものではなく、
例えば図11に示す関数ｆ3 （Ｓus）のように一部に曲線
を含むものであってもよく、図12に示す関数ｆ4 （Ｓu
s）のようにｆ4 （Ｓus）＝０とならない部分を含むも
のであってもよい。さらに図13に示す関数ｆ4 （Ｓus）
のように低濃度側と高濃度側とで圧縮の程度が高なる関
数を用いるようにしてもよい。さらに、図14に示す関数
ｆ6 （Ｓus）のように、ｆ6 （Ｓus）＝０となる点ｄに
おいて折れ線とならずに微係数が連続する関数としても
よい。点ｄが折れ線となる場合、この点ｄに相当する濃
度領域に、オリジナル画像には何ら特別の輪郭は存在し
ないにも拘らず、処理済画像信号Ｓprocに基づいて再生
表示された可視画像に偽輪郭が生じる場合がある。そこ
で、点ｄの部分を微係数が連続する関数とすることによ
り偽輪郭が生じることも防止される。
【００５１】上記各実施例は蓄積性蛍光体シートに記録
されたＸ線画像を読み取って画像信号を得るシステムで
あるが、本発明は蓄積性蛍光体シートを用いるシステム
に限られるものではなく、Ｘ線フイルム等に記録された
Ｘ線画像、その他記録シートに記録された一般の画像等
を読み取って画像信号を得るシステムに広く適用できる
ものである。
されたＸ線画像を読み取って画像信号を得るシステムで
あるが、本発明は蓄積性蛍光体シートを用いるシステム
に限られるものではなく、Ｘ線フイルム等に記録された
Ｘ線画像、その他記録シートに記録された一般の画像等
を読み取って画像信号を得るシステムに広く適用できる
ものである。
【００５２】
【発明の効果】以上詳細に説明したように、本発明によ
るダイナミックレンジ圧縮処理方法は、オリジナル画像
のコントラストに応じて圧縮処理を行うための関数の適
用信号範囲および／または値の大きさを決定するように
したため、オリジナル画像の大小に拘らず、高濃度領
域、低濃度領域あるいは画像全体のダイナミックレンジ
を圧縮でき、画像の適性観察領域を広げることができ
る。
るダイナミックレンジ圧縮処理方法は、オリジナル画像
のコントラストに応じて圧縮処理を行うための関数の適
用信号範囲および／または値の大きさを決定するように
したため、オリジナル画像の大小に拘らず、高濃度領
域、低濃度領域あるいは画像全体のダイナミックレンジ
を圧縮でき、画像の適性観察領域を広げることができ
る。
【図１】Ｘ線撮影の装置の一例の概略図
【図２】Ｘ線画像読取装置の一例を表した斜視図
【図３】オリジナル画像上の各画素点とこの各画素点に
対応するオリジナル画像信号Ｓorg を表す図
対応するオリジナル画像信号Ｓorg を表す図
【図４】オリジナル画像信号のコントラストの例を表す
図
図
【図５】ボケマスク信号を変数とした単調減少関数の例
を表す図
を表す図
【図６】オリジナル画像信号のコントラストの他の例を
表す図
表す図
【図７】ボケマスク信号を変数とした単調減少関数の他
の例を表す図
の例を表す図
【図８】オリジナル画像信号のコントラストのさらに他
の例を表す図
の例を表す図
【図９】ボケマスク信号を変数とした単調減少関数のさ
らに他の例を表す図
らに他の例を表す図
【図１０】オリジナル画像信号のコントラストのさらに
他の例を表す図
他の例を表す図
【図１１】ボケマスク信号を変数とした単調減少関数の
さらに他の例を表す図
さらに他の例を表す図
【図１２】ボケマスク信号を変数とした単調減少関数の
さらに他の例を表す図
さらに他の例を表す図
【図１３】ボケマスク信号を変数とした単調減少関数の
さらに他の例を表す図
さらに他の例を表す図
【図１４】ボケマスク信号を変数とした単調減少関数の
さらに他の例を表す図
さらに他の例を表す図
【図１５】オリジナル画像上のある方向（ｘ方向）に沿
うオリジナル画像信号の値の一例をプロットした図
うオリジナル画像信号の値の一例をプロットした図
【図１６】高濃度域のコントラストを下げることにより
ダイナミックレンジを圧縮した場合の、それぞれ画像信
号の変換グラフの一例を示した図、および画像上のｘ方
向に沿う変換後の画像信号の値を表す図
ダイナミックレンジを圧縮した場合の、それぞれ画像信
号の変換グラフの一例を示した図、および画像上のｘ方
向に沿う変換後の画像信号の値を表す図
【図１７】低濃度域のコントラストを下げることにより
ダイナミックレンジを圧縮した場合の、それぞれ画像信
号の変換グラフの一例を示した図、および画像上のｘ方
向に沿う変換後の画像信号の値を表す図
ダイナミックレンジを圧縮した場合の、それぞれ画像信
号の変換グラフの一例を示した図、および画像上のｘ方
向に沿う変換後の画像信号の値を表す図
10 Ｘ線撮影装置 14 蓄積性蛍光体シート 20 読取部 40 画像処理表示部
Claims (2)
- 【請求項１】 オリジナル画像を表すオリジナル画像信
号Ｓorg の各画素点に対応して該各画素点の周囲の所定
範囲内の前記オリジナル画像信号Ｓorg を平均化するこ
とによりボケマスク信号Ｓusを求め、該ボケマスク信号
Ｓusの値が増大するにつれて単調減少する関数をｆ
1 （Ｓus）としたとき、式 Ｓproc＝Ｓorg ＋ｆ1 （Ｓus） に従って前記オリジナル画像信号を処理して前記オリジ
ナル画像よりもダイナミックレンジの狭い画像を担持す
る処理済画像信号Ｓprocを求める画像のダイナミックレ
ンジ圧縮処理方法において、 前記オリジナル画像信号に基づいて該オリジナル画像信
号のコントラストを求め、 該コントラストに応じて前記関数ｆ1 （Ｓus）の適用信
号範囲および／または値の大きさを決定することを特徴
とする画像のダイナミックレンジ圧縮処理方法。 - 【請求項２】 前記単調減少する関数ｆ1 （Ｓus）が、
少なくとも一部に曲線を含むことを特徴とする請求項１
記載の画像のダイナミックレンジ圧縮処理方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP09127693A JP3196864B2 (ja) | 1993-04-19 | 1993-04-19 | 画像のダイナミックレンジ圧縮処理方法 |
| US08/229,217 US5454053A (en) | 1993-04-19 | 1994-04-18 | Method for compressing dynamic ranges of images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP09127693A JP3196864B2 (ja) | 1993-04-19 | 1993-04-19 | 画像のダイナミックレンジ圧縮処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06303430A true JPH06303430A (ja) | 1994-10-28 |
| JP3196864B2 JP3196864B2 (ja) | 2001-08-06 |
Family
ID=14021940
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP09127693A Expired - Lifetime JP3196864B2 (ja) | 1993-04-19 | 1993-04-19 | 画像のダイナミックレンジ圧縮処理方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5454053A (ja) |
| JP (1) | JP3196864B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルｘ線イメージング検出器のダイナミックレンジ拡張システム |
| WO2006129529A1 (ja) | 2005-06-02 | 2006-12-07 | Konica Minolta Holdings, Inc. | 画像処理方法及び画像処理装置 |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4404047C1 (de) * | 1994-02-09 | 1995-08-31 | Daimler Benz Aerospace Ag | Verfahren zur Beeinflussung der Bildqualität von durch eine elektronische Bildaufnahmevorrichtung gelieferten Bildern |
| US5761333A (en) * | 1995-01-31 | 1998-06-02 | General Electric Company | Contrast enhancement for CT systems |
| JPH08307683A (ja) * | 1995-04-28 | 1996-11-22 | Minolta Co Ltd | 画像処理装置 |
| EP1168243B1 (en) * | 1995-09-29 | 2004-06-09 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| JP3380660B2 (ja) * | 1995-10-25 | 2003-02-24 | 大日本スクリーン製造株式会社 | 画像処理装置 |
| EP0813336B1 (en) * | 1996-06-12 | 2007-08-08 | FUJIFILM Corporation | Image processing method and apparatus |
| JP3363735B2 (ja) * | 1996-06-26 | 2003-01-08 | 松下電器産業株式会社 | Ｘ線画像装置 |
| JP4124380B2 (ja) * | 1996-08-26 | 2008-07-23 | 富士フイルム株式会社 | 画像処理再生装置 |
| US5835618A (en) * | 1996-09-27 | 1998-11-10 | Siemens Corporate Research, Inc. | Uniform and non-uniform dynamic range remapping for optimum image display |
| US6330362B1 (en) | 1996-11-12 | 2001-12-11 | Texas Instruments Incorporated | Compression for multi-level screened images |
| DE19722761B4 (de) * | 1997-06-02 | 2006-06-01 | Hell Gravure Systems Gmbh | Verfahren zur Aufbereitung von Bilddaten |
| EP1667066B1 (en) * | 2003-09-11 | 2020-01-08 | Panasonic Intellectual Property Corporation of America | Visual processing device, visual processing method, visual processing program, integrated circuit, display device, imaging device, and mobile information terminal |
| EP1667064A4 (en) * | 2003-09-11 | 2009-06-10 | Panasonic Corp | IMAGE PROCESSING DEVICE, METHOD AND PROGRAM AND SEMICONDUCTOR ELEMENT |
| EP1667065B1 (en) * | 2003-09-11 | 2018-06-06 | Panasonic Intellectual Property Corporation of America | Visual processing apparatus, visual processing method, visual processing program, and semiconductor device |
| WO2005027040A1 (ja) * | 2003-09-11 | 2005-03-24 | Matsushita Electric Industrial Co., Ltd. | 画像処理装置、画像処理方法および画像処理プログラム |
| US9829984B2 (en) * | 2013-05-23 | 2017-11-28 | Fastvdo Llc | Motion-assisted visual language for human computer interfaces |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4315318A (en) * | 1978-12-26 | 1982-02-09 | Fuji Photo Film Co., Ltd. | Method and apparatus for processing a radiation image |
| JPS55116339A (en) * | 1979-02-28 | 1980-09-06 | Fuji Photo Film Co Ltd | Method and device for processing gradation of radiation picture of chest |
| US4571635A (en) * | 1984-02-17 | 1986-02-18 | Minnesota Mining And Manufacturing Company | Method of image enhancement by raster scanning |
| JP2663189B2 (ja) | 1990-01-29 | 1997-10-15 | 富士写真フイルム株式会社 | 画像のダイナミックレンジ圧縮処理方法 |
-
1993
- 1993-04-19 JP JP09127693A patent/JP3196864B2/ja not_active Expired - Lifetime
-
1994
- 1994-04-18 US US08/229,217 patent/US5454053A/en not_active Expired - Lifetime
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001314396A (ja) * | 1999-12-30 | 2001-11-13 | General Electric Co <Ge> | デジタルｘ線イメージング検出器のダイナミックレンジ拡張システム |
| WO2006129529A1 (ja) | 2005-06-02 | 2006-12-07 | Konica Minolta Holdings, Inc. | 画像処理方法及び画像処理装置 |
| US8131101B2 (en) | 2005-06-02 | 2012-03-06 | Konica Minolta Holdings, Inc. | Image processing method and image processing apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| US5454053A (en) | 1995-09-26 |
| JP3196864B2 (ja) | 2001-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2663189B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| US4747052A (en) | Radiation image processing | |
| US5067163A (en) | Method for determining a desired image signal range from an image having a single background | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| JP3599435B2 (ja) | 画像データ補間演算方法および装置 | |
| JPH06304159A (ja) | 画像の撮影体位認識方法 | |
| JPH0792830B2 (ja) | 放射線画像情報読取方法 | |
| JPH06274614A (ja) | 画像処理方法 | |
| JPH06274615A (ja) | 画像処理方法 | |
| US5832142A (en) | Interpolating operation method and apparatus for image signals | |
| JP2952418B2 (ja) | 放射線画像信号生成方法 | |
| JP3720081B2 (ja) | 階調補正方法および装置 | |
| JPH01237887A (ja) | 放射線画像の粒状雑音低減方法 | |
| JPH11205682A (ja) | エネルギーサブトラクション画像生成方法 | |
| JP3399684B2 (ja) | 画像処理方法および装置 | |
| JP2739385B2 (ja) | 放射線画像読取条件及び／又は画像処理条件決定方法および装置 | |
| JP2663189C (ja) | ||
| JP3241488B2 (ja) | 画像処理方法 | |
| JP3165530B2 (ja) | 放射線画像読取条件決定方法および装置 | |
| JPH0786925B2 (ja) | 放射線画像処理方法および装置 | |
| JP3165461B2 (ja) | 画像処理方法、画像読取方法および画像再生方法 | |
| JPH02272681A (ja) | 放射線画像処理装置 | |
| JP2631032B2 (ja) | 放射線画像のエネルギーサブトラクション方法および装置 | |
| JPH0281276A (ja) | 画像処理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20010522 |
|
| 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: 20080608 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: 20080608 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: 20080608 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090608 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090608 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100608 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100608 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110608 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110608 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120608 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120608 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130608 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term |