JP3399684B2 - 画像処理方法および装置 - Google Patents
画像処理方法および装置Info
- Publication number
- JP3399684B2 JP3399684B2 JP04518895A JP4518895A JP3399684B2 JP 3399684 B2 JP3399684 B2 JP 3399684B2 JP 04518895 A JP04518895 A JP 04518895A JP 4518895 A JP4518895 A JP 4518895A JP 3399684 B2 JP3399684 B2 JP 3399684B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- image data
- directions
- qorg
- pixel
- 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
- 238000003672 processing method Methods 0.000 title claims description 16
- 238000000034 method Methods 0.000 claims description 30
- 238000001914 filtration Methods 0.000 claims description 23
- 230000004044 response Effects 0.000 claims description 16
- 230000005855 radiation Effects 0.000 description 21
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 14
- 238000003745 diagnosis Methods 0.000 description 6
- 230000005284 excitation Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/409—Edge or detail enhancement; Noise or error suppression
- H04N1/4092—Edge or detail enhancement
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Radiography Using Non-Light Waves (AREA)
- Facsimile Image Signal Circuits (AREA)
Description
【０００１】
【産業上の利用分野】本発明は画像処理方法、詳しくは
周波数強調処理に関し、さらに詳細には、再生画像の
Ｘ、Ｙ方向の鮮鋭度を互いに略等しく揃えられるように
した周波数強調処理に関するものである。
周波数強調処理に関し、さらに詳細には、再生画像の
Ｘ、Ｙ方向の鮮鋭度を互いに略等しく揃えられるように
した周波数強調処理に関するものである。
【０００２】また本発明は、このような周波数強調処理
を実施するための画像処理装置に関するものである。
を実施するための画像処理装置に関するものである。
【０００３】
【従来の技術】ある種の蛍光体に放射線（Ｘ線、α線、
β線、γ線、電子線、紫外線等）を照射すると、この放
射線のエネルギーの一部がその蛍光体中に蓄積され、そ
の後該蛍光体に可視光等の励起光を照射すると、蓄積さ
れたエネルギーに応じて蛍光体が輝尽発光を示す。この
ような性質を示す蛍光体を蓄積性蛍光体と言う。
β線、γ線、電子線、紫外線等）を照射すると、この放
射線のエネルギーの一部がその蛍光体中に蓄積され、そ
の後該蛍光体に可視光等の励起光を照射すると、蓄積さ
れたエネルギーに応じて蛍光体が輝尽発光を示す。この
ような性質を示す蛍光体を蓄積性蛍光体と言う。
【０００４】この蓄積性蛍光体を利用して、人体等の被
写体の放射線画像情報を一旦蓄積性蛍光体のシート（以
下、蓄積性蛍光体シートと称する）に記録し、該シート
を励起光で走査して輝尽発光させ、この輝尽発光光を光
電的に読み取って画像信号を得、この画像信号を処理し
て診断適性の良い被写体の放射線画像を得る放射線画像
情報記録再生システムが提案されている（例えば特開昭
５５―１２４２９号、同５５―１１６３４０号、同５５
―１６３４７２号、同５６―１１３９５号、同５６―１
０４６４５号など）。
写体の放射線画像情報を一旦蓄積性蛍光体のシート（以
下、蓄積性蛍光体シートと称する）に記録し、該シート
を励起光で走査して輝尽発光させ、この輝尽発光光を光
電的に読み取って画像信号を得、この画像信号を処理し
て診断適性の良い被写体の放射線画像を得る放射線画像
情報記録再生システムが提案されている（例えば特開昭
５５―１２４２９号、同５５―１１６３４０号、同５５
―１６３４７２号、同５６―１１３９５号、同５６―１
０４６４５号など）。
【０００５】ところで、このような放射線画像情報記録
再生システムにあっては、再生された放射線画像におい
て診断上重要な空間周波数領域が明確に示されず、コン
トラスト、鮮鋭度が低くて診断適性が劣ることがあっ
た。このような事情に鑑みて本出願人は既に、診断上有
効な超低空間周波数以上の周波数成分を強調することに
より、診断適性を向上させうる画像処理方法を提案した
（特開昭５５―１６３４７２号公報参照）。
再生システムにあっては、再生された放射線画像におい
て診断上重要な空間周波数領域が明確に示されず、コン
トラスト、鮮鋭度が低くて診断適性が劣ることがあっ
た。このような事情に鑑みて本出願人は既に、診断上有
効な超低空間周波数以上の周波数成分を強調することに
より、診断適性を向上させうる画像処理方法を提案した
（特開昭５５―１６３４７２号公報参照）。
【０００６】この画像処理方法は、Ｘ、Ｙ方向に並ぶ画
素のそれぞれについての画素値を示しているデジタル画
像データを、各画素についての原画像データをＱorg 、
この画素に関する超低空間周波数に対応するボケマスク
（非鮮鋭マスク）データをＱus、強調係数をβとしたと
き、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけるようにしたものである。
素のそれぞれについての画素値を示しているデジタル画
像データを、各画素についての原画像データをＱorg 、
この画素に関する超低空間周波数に対応するボケマスク
（非鮮鋭マスク）データをＱus、強調係数をβとしたと
き、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけるようにしたものである。
【０００７】ここで、超低空間周波数に対応するボケマ
スクデータＱusとは、オリジナル画像を超低空間周波数
成分より低い周波数成分しか含まないようにぼかした非
鮮鋭画像の各画素の濃度に対応するデータを指す。
スクデータＱusとは、オリジナル画像を超低空間周波数
成分より低い周波数成分しか含まないようにぼかした非
鮮鋭画像の各画素の濃度に対応するデータを指す。
【０００８】
【発明が解決しようとする課題】上記の演算処理で得ら
れた新たな画像データＱ’に基づいて放射線画像を再生
すれば、診断上有効な超低空間周波数以上の周波数成分
が強調され、診断適性に優れた再生放射線画像が得られ
るようになる。
れた新たな画像データＱ’に基づいて放射線画像を再生
すれば、診断上有効な超低空間周波数以上の周波数成分
が強調され、診断適性に優れた再生放射線画像が得られ
るようになる。
【０００９】しかしその場合、画像再生装置がＸ、Ｙ走
査により画像を再生するもの、例えばＣＲＴ表示装置、
サーマルヘッドを用いる感熱記録装置、光走査記録装置
等であると、Ｘ、Ｙ方向の鮮鋭度が互いに異なってしま
うという不具合が従来より認められていた。この不具合
は、画像再生装置のＸ、Ｙ方向のレスポンスが互いに異
なっていることに起因していると考えられる。
査により画像を再生するもの、例えばＣＲＴ表示装置、
サーマルヘッドを用いる感熱記録装置、光走査記録装置
等であると、Ｘ、Ｙ方向の鮮鋭度が互いに異なってしま
うという不具合が従来より認められていた。この不具合
は、画像再生装置のＸ、Ｙ方向のレスポンスが互いに異
なっていることに起因していると考えられる。
【００１０】再生画像の鮮鋭度は、上述した演算処理に
おける強調係数βを変えることによって調節可能である
ので、強調係数βを比較的大きく設定して、レスポンス
が低い方向の鮮鋭度も十分となるようにすることも可能
であるが、そのようにすると、レスポンスが高い方向に
ついて過度の周波数強調がなされてアーチファクト（偽
画像）が発生することがある。
おける強調係数βを変えることによって調節可能である
ので、強調係数βを比較的大きく設定して、レスポンス
が低い方向の鮮鋭度も十分となるようにすることも可能
であるが、そのようにすると、レスポンスが高い方向に
ついて過度の周波数強調がなされてアーチファクト（偽
画像）が発生することがある。
【００１１】本発明は上記の事情に鑑みてなされたもの
であり、先に説明した周波数強調処理を行なった上で、
画像再生装置において再生される画像のＸ、Ｙ方向の鮮
鋭度を略一定に揃えることができる画像処理方法を提供
することを目的とするものである。
であり、先に説明した周波数強調処理を行なった上で、
画像再生装置において再生される画像のＸ、Ｙ方向の鮮
鋭度を略一定に揃えることができる画像処理方法を提供
することを目的とするものである。
【００１２】また本発明は、そのような画像処理方法を
実施することができる画像処理装置を提供することを目
的とするものである。
実施することができる画像処理装置を提供することを目
的とするものである。
【００１３】
【課題を解決するための手段】本発明による一つの画像
処理方法は、前述したようにＸ、Ｙ方向に並ぶ画素のそ
れぞれについての画素値を示しているデジタル画像デー
タを、各画素についての原画像データをＱorg 、この画
素に関する超低空間周波数に対応するボケマスクデータ
をＱus、強調係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） ……（１） なる演算処理にかけて、この演算処理後の画像データ
を、Ｘ、Ｙ走査により画像を再生する画像再生装置に送
る画像処理方法において、上記演算処理を、画像再生装
置のＸ、Ｙ方向のレスポンスの差に応じて、該Ｘ、Ｙ方
向について互いに異なる特性の下に行なうことにより、
この画像再生装置で再生される画像のＸ、Ｙ方向の最終
的な鮮鋭度を互いに略等しく揃えることを特徴とするも
のである。
処理方法は、前述したようにＸ、Ｙ方向に並ぶ画素のそ
れぞれについての画素値を示しているデジタル画像デー
タを、各画素についての原画像データをＱorg 、この画
素に関する超低空間周波数に対応するボケマスクデータ
をＱus、強調係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） ……（１） なる演算処理にかけて、この演算処理後の画像データ
を、Ｘ、Ｙ走査により画像を再生する画像再生装置に送
る画像処理方法において、上記演算処理を、画像再生装
置のＸ、Ｙ方向のレスポンスの差に応じて、該Ｘ、Ｙ方
向について互いに異なる特性の下に行なうことにより、
この画像再生装置で再生される画像のＸ、Ｙ方向の最終
的な鮮鋭度を互いに略等しく揃えることを特徴とするも
のである。
【００１４】なお、上記の演算処理をＸ、Ｙ方向につい
て互いに異なる特性の下に行なうためには、より具体的
に、上記演算処理をまずＸ、Ｙの一方向について強調係
数をβ1 に設定して行ない、次いでこの演算処理後の画
像を原画像データとして、Ｘ、Ｙの他方向について強調
係数をβ2 （β1 ≠β2 ）に設定して上記演算処理を行
なえばよい。さらには、このように１次元の演算処理を
２回行なう代わりに、それと等価の２次元フィルタリン
グ処理を行なうようにしてもよい。
て互いに異なる特性の下に行なうためには、より具体的
に、上記演算処理をまずＸ、Ｙの一方向について強調係
数をβ1 に設定して行ない、次いでこの演算処理後の画
像を原画像データとして、Ｘ、Ｙの他方向について強調
係数をβ2 （β1 ≠β2 ）に設定して上記演算処理を行
なえばよい。さらには、このように１次元の演算処理を
２回行なう代わりに、それと等価の２次元フィルタリン
グ処理を行なうようにしてもよい。
【００１５】また、本発明による別の画像処理方法は、
前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、前記
（１）式で示される演算処理にかけて、この演算処理後
の画像データを、Ｘ、Ｙ走査により画像を再生する画像
再生装置に送る画像処理方法において、上記演算処理
を、画像再生装置のＸ、Ｙ方向のレスポンスの差に応じ
て、該Ｘ、Ｙ方向のうちの一方向のみについて行なうこ
とにより、この画像再生装置で再生される画像のＸ、Ｙ
方向の最終的な鮮鋭度を互いに略等しく揃えることを特
徴とするものである。
前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、前記
（１）式で示される演算処理にかけて、この演算処理後
の画像データを、Ｘ、Ｙ走査により画像を再生する画像
再生装置に送る画像処理方法において、上記演算処理
を、画像再生装置のＸ、Ｙ方向のレスポンスの差に応じ
て、該Ｘ、Ｙ方向のうちの一方向のみについて行なうこ
とにより、この画像再生装置で再生される画像のＸ、Ｙ
方向の最終的な鮮鋭度を互いに略等しく揃えることを特
徴とするものである。
【００１６】一方、本発明による一つの画像処理装置
は、前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれに
ついての画素値を示しているデジタル画像データを、前
記（１）式で示される演算処理にかけて、この演算処理
後の画像データを、Ｘ、Ｙ走査により画像を再生する画
像再生装置に送る画像処理装置において、上記演算処理
を、Ｘ、Ｙ方向について互いに異なる特性の下に行ない
得るように構成されたことを特徴とするものである。
は、前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれに
ついての画素値を示しているデジタル画像データを、前
記（１）式で示される演算処理にかけて、この演算処理
後の画像データを、Ｘ、Ｙ走査により画像を再生する画
像再生装置に送る画像処理装置において、上記演算処理
を、Ｘ、Ｙ方向について互いに異なる特性の下に行ない
得るように構成されたことを特徴とするものである。
【００１７】なお、上記の演算処理をＸ、Ｙ方向につい
て互いに異なる特性の下に行なうためには、より具体的
に、上記演算処理をまずＸ、Ｙの一方向について強調係
数をβ1 に設定して行ない、次いでこの演算処理後の画
像を原画像データとして、Ｘ、Ｙの他方向について強調
係数をβ2 （β1 ≠β2 ）に設定して上記演算処理を行
うように構成すればよい。さらには、このように１次元
の演算処理を２回行なう代わりに、それと等価の２次元
フィルタリング処理を行なうように構成してもよい。
て互いに異なる特性の下に行なうためには、より具体的
に、上記演算処理をまずＸ、Ｙの一方向について強調係
数をβ1 に設定して行ない、次いでこの演算処理後の画
像を原画像データとして、Ｘ、Ｙの他方向について強調
係数をβ2 （β1 ≠β2 ）に設定して上記演算処理を行
うように構成すればよい。さらには、このように１次元
の演算処理を２回行なう代わりに、それと等価の２次元
フィルタリング処理を行なうように構成してもよい。
【００１８】また、本発明による別の画像処理装置は、
前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、前記
（１）式で示される演算処理にかけて、この演算処理後
の画像データを、Ｘ、Ｙ走査により画像を再生する画像
再生装置に送る画像処理装置において、上記演算処理
を、Ｘ、Ｙ方向のうちの一方向のみについて行ない得る
ように構成されたことを特徴とするものである。
前述したようにＸ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、前記
（１）式で示される演算処理にかけて、この演算処理後
の画像データを、Ｘ、Ｙ走査により画像を再生する画像
再生装置に送る画像処理装置において、上記演算処理
を、Ｘ、Ｙ方向のうちの一方向のみについて行ない得る
ように構成されたことを特徴とするものである。
【００１９】
【作用および発明の効果】以上の通り本発明によれば、
演算処理をＸ、Ｙ方向について互いに異なる特性の下に
行なうことにより、あるいはＸ、Ｙ方向のうちの一方向
のみについて行なうことにより、画像再生装置のＸ、Ｙ
方向のレスポンスの差を補償して、再生画像のＸ、Ｙ方
向の最終的な鮮鋭度を互いに略等しく揃えることができ
る。そこで本発明によれば、Ｘ、Ｙ方向のうちの一方向
の鮮鋭度が不足することを防止でき、また、Ｘ、Ｙ方向
のうちの一方向について過度の周波数強調処理をかけて
アーチファクトを発生させてしまうことも防止できる。
演算処理をＸ、Ｙ方向について互いに異なる特性の下に
行なうことにより、あるいはＸ、Ｙ方向のうちの一方向
のみについて行なうことにより、画像再生装置のＸ、Ｙ
方向のレスポンスの差を補償して、再生画像のＸ、Ｙ方
向の最終的な鮮鋭度を互いに略等しく揃えることができ
る。そこで本発明によれば、Ｘ、Ｙ方向のうちの一方向
の鮮鋭度が不足することを防止でき、また、Ｘ、Ｙ方向
のうちの一方向について過度の周波数強調処理をかけて
アーチファクトを発生させてしまうことも防止できる。
【００２０】なお本発明は、前述した蓄積性蛍光体シー
トを用いる放射線画像情報記録再生システムにおいて適
用されれば、診断適性が極めて高い放射線画像を再生で
きるので特に好ましいが、このような放射線画像情報記
録再生システム以外に適用された場合でも、上記の通り
の作用、効果を奏することは勿論である。
トを用いる放射線画像情報記録再生システムにおいて適
用されれば、診断適性が極めて高い放射線画像を再生で
きるので特に好ましいが、このような放射線画像情報記
録再生システム以外に適用された場合でも、上記の通り
の作用、効果を奏することは勿論である。
【００２１】
【実施例】以下、図面に示す実施例に基づいて本発明を
詳細に説明する。図１は本発明の方法により周波数強調
処理を行なう放射線画像情報記録再生システムの一例を
示すものである。このシステムは、前述の蓄積性蛍光体
シートに記録された放射線画像情報を読み取り、その読
み取った放射線画像を一例として感熱記録装置により再
生記録するものである。
詳細に説明する。図１は本発明の方法により周波数強調
処理を行なう放射線画像情報記録再生システムの一例を
示すものである。このシステムは、前述の蓄積性蛍光体
シートに記録された放射線画像情報を読み取り、その読
み取った放射線画像を一例として感熱記録装置により再
生記録するものである。
【００２２】例えばＸ線等の放射線が人体等の被写体を
介して照射されることによりこの被写体の透過放射線画
像情報を蓄積記録した蓄積性蛍光体シート10は、エンド
レスベルト等のシート搬送手段11により、副走査のため
に矢印Ｙ方向に搬送される。半導体レーザ等の励起光源
12から射出された励起光（読取光）としてのレーザビー
ム13は、高速回転する回転多面鏡14によって反射偏向さ
れ、通常ｆ・θレンズからなる走査レンズ18によって集
束され、ミラー19で反射して蓄積性蛍光体シート10上を
上記副走査方向Ｙと略直角な矢印Ｘ方向に主走査する。
介して照射されることによりこの被写体の透過放射線画
像情報を蓄積記録した蓄積性蛍光体シート10は、エンド
レスベルト等のシート搬送手段11により、副走査のため
に矢印Ｙ方向に搬送される。半導体レーザ等の励起光源
12から射出された励起光（読取光）としてのレーザビー
ム13は、高速回転する回転多面鏡14によって反射偏向さ
れ、通常ｆ・θレンズからなる走査レンズ18によって集
束され、ミラー19で反射して蓄積性蛍光体シート10上を
上記副走査方向Ｙと略直角な矢印Ｘ方向に主走査する。
【００２３】こうしてレーザビーム13が照射されたシー
ト10の箇所からは、蓄積記録されている放射線画像情報
に応じた光量の輝尽発光光15が発散され、この輝尽発光
光15は集光体16によって集光され、光検出器としてのフ
ォトマルチプライヤー（光電子増倍管）17によって光電
的に検出される。
ト10の箇所からは、蓄積記録されている放射線画像情報
に応じた光量の輝尽発光光15が発散され、この輝尽発光
光15は集光体16によって集光され、光検出器としてのフ
ォトマルチプライヤー（光電子増倍管）17によって光電
的に検出される。
【００２４】上記集光体16はアクリル板等の導光性材料
を成形して作られたものであり、直線状をなす入射端面
16ａが蓄積性蛍光体シート10上のビーム走査線に沿って
延びるように配され、円環状に形成された出射端面16ｂ
に上記フォトマルチプライヤー17の受光面が結合されて
いる。上記入射端面16ａから集光体16内に入射した輝尽
発光光15は、該集光体16の内部を全反射を繰り返して進
み、出射端面16ｂから出射してフォトマルチプライヤー
17に受光され、前記放射線画像情報を担持する輝尽発光
光15の光量がこのフォトマルチプライヤー17によって検
出される。
を成形して作られたものであり、直線状をなす入射端面
16ａが蓄積性蛍光体シート10上のビーム走査線に沿って
延びるように配され、円環状に形成された出射端面16ｂ
に上記フォトマルチプライヤー17の受光面が結合されて
いる。上記入射端面16ａから集光体16内に入射した輝尽
発光光15は、該集光体16の内部を全反射を繰り返して進
み、出射端面16ｂから出射してフォトマルチプライヤー
17に受光され、前記放射線画像情報を担持する輝尽発光
光15の光量がこのフォトマルチプライヤー17によって検
出される。
【００２５】フォトマルチプライヤー17のアナログ出力
信号（画像信号）Ｓは対数増幅器20によって増幅され、
Ａ／Ｄ変換器21において所定の収録スケールファクター
でデジタル化される。こうして得られた、２次元画像を
担持するデジタルの画像データＤは、画像処理装置27に
おいて階調処理、周波数強調処理等の画像処理を受け
る。この画像処理を受けた画像データＤ’は、サーマル
ヘッド28ａを用いる感熱記録装置28に入力され、該感熱
記録装置28において画像データＤ’が示す２次元画像が
ハードコピー29に再生記録される。
信号（画像信号）Ｓは対数増幅器20によって増幅され、
Ａ／Ｄ変換器21において所定の収録スケールファクター
でデジタル化される。こうして得られた、２次元画像を
担持するデジタルの画像データＤは、画像処理装置27に
おいて階調処理、周波数強調処理等の画像処理を受け
る。この画像処理を受けた画像データＤ’は、サーマル
ヘッド28ａを用いる感熱記録装置28に入力され、該感熱
記録装置28において画像データＤ’が示す２次元画像が
ハードコピー29に再生記録される。
【００２６】次に、上記画像処理装置27における周波数
強調処理について詳しく説明する。この画像処理装置27
は、１画像分の画素値を示しているデジタル画像データ
を一時的に格納しておくバッファメモリ23、後述する強
調係数βを記憶しておくメモリ24、およびフィルタリン
グ処理部25等からなり、まずＸ方向（主走査方向）に延
びる１走査ライン毎に、そのラインを構成している各画
素に関する原画像データＱorg に対して、下記の演算処
理 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） に相当するフィルタリング処理を行なう。なおβは強調
係数、Ｑusは各画素に関する超低空間周波数に対応する
ボケマスクデータで、マスクサイズをＮ、そのマスクサ
イズ内の画素の画素値の合計値をΣとすると、Ｑus＝Σ
／Ｎである。
強調処理について詳しく説明する。この画像処理装置27
は、１画像分の画素値を示しているデジタル画像データ
を一時的に格納しておくバッファメモリ23、後述する強
調係数βを記憶しておくメモリ24、およびフィルタリン
グ処理部25等からなり、まずＸ方向（主走査方向）に延
びる１走査ライン毎に、そのラインを構成している各画
素に関する原画像データＱorg に対して、下記の演算処
理 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） に相当するフィルタリング処理を行なう。なおβは強調
係数、Ｑusは各画素に関する超低空間周波数に対応する
ボケマスクデータで、マスクサイズをＮ、そのマスクサ
イズ内の画素の画素値の合計値をΣとすると、Ｑus＝Σ
／Ｎである。
【００２７】例えば、Ｎ＝３とし、Ｘ方向に並ぶ３つの
画素の原画像データをｘ1 、ｘ2 、ｘ3 とすると、中央
の画素の原画像データｘ2 については、 Ｑ’＝ｘ2 ＋β｛ｘ2 −（ｘ1 ＋ｘ2 ＋ｘ3 ）／３｝ ＝（−β／３）ｘ1 ＋（１＋２β／３）ｘ2 ＋（−β／３）ｘ3 であるから、 ｘ2 ’＝ａｘ1 ＋ｂｘ2 ＋ｃｘ3 （ｂ＝１＋２β／３、ａ＝ｃ＝−β／３） ……（２） と置き代えるフィルタリング処理がなされる。
画素の原画像データをｘ1 、ｘ2 、ｘ3 とすると、中央
の画素の原画像データｘ2 については、 Ｑ’＝ｘ2 ＋β｛ｘ2 −（ｘ1 ＋ｘ2 ＋ｘ3 ）／３｝ ＝（−β／３）ｘ1 ＋（１＋２β／３）ｘ2 ＋（−β／３）ｘ3 であるから、 ｘ2 ’＝ａｘ1 ＋ｂｘ2 ＋ｃｘ3 （ｂ＝１＋２β／３、ａ＝ｃ＝−β／３） ……（２） と置き代えるフィルタリング処理がなされる。
【００２８】そして画像処理装置27は、以上の通りのフ
ィルタリング処理を受けた後の画像データを原画像デー
タとして、今度はＹ方向（副走査方向）に延びる１本の
画素列毎に、上記（２）式で示されるフィルタリング処
理を行なう。
ィルタリング処理を受けた後の画像データを原画像デー
タとして、今度はＹ方向（副走査方向）に延びる１本の
画素列毎に、上記（２）式で示されるフィルタリング処
理を行なう。
【００２９】ここで、１回目のフィルタリング処理と２
回目のフィルタリング処理とでは、強調係数βがそれぞ
れβ1 、β2 と、互いに異なる値に設定される。これら
の強調係数β1 、β2 は、感熱記録装置28のＸ、Ｙ方向
のレスポンスの差に応じて定められるものである。すな
わち、感熱記録装置28のＸ、Ｙ方向のレスポンス特性が
各々図２の曲線ａ、ｂであるとし、Ｘ、Ｙ方向で主に鮮
鋭度を揃えたい空間周波数をαとし、この空間周波数α
に対応するＸ、Ｙ方向のＭＴＦをそれぞれｍ1、ｍ2 と
すると、強調係数β1 、β2 はそれぞれ図３に示すよう
に互いに相違して、ｍ1 ／ｍ2 ＝（１＋β2 ）／（１＋
β1 ）の関係を満たす値に設定される。
回目のフィルタリング処理とでは、強調係数βがそれぞ
れβ1 、β2 と、互いに異なる値に設定される。これら
の強調係数β1 、β2 は、感熱記録装置28のＸ、Ｙ方向
のレスポンスの差に応じて定められるものである。すな
わち、感熱記録装置28のＸ、Ｙ方向のレスポンス特性が
各々図２の曲線ａ、ｂであるとし、Ｘ、Ｙ方向で主に鮮
鋭度を揃えたい空間周波数をαとし、この空間周波数α
に対応するＸ、Ｙ方向のＭＴＦをそれぞれｍ1、ｍ2 と
すると、強調係数β1 、β2 はそれぞれ図３に示すよう
に互いに相違して、ｍ1 ／ｍ2 ＝（１＋β2 ）／（１＋
β1 ）の関係を満たす値に設定される。
【００３０】上記のようなフィルタリング処理を行なっ
て、超低空間周波数以上の周波数成分を強調することに
より、感熱記録装置28によりハードコピー29として再生
された放射線画像においては、診断に重要な周波数領域
が大幅に強調され、鮮鋭度が向上して該放射線画像の診
断適性が向上する。以上の処理によるメリットは、前記
特開昭５５―１６３４７２号公報にさらに詳しく記載さ
れている。
て、超低空間周波数以上の周波数成分を強調することに
より、感熱記録装置28によりハードコピー29として再生
された放射線画像においては、診断に重要な周波数領域
が大幅に強調され、鮮鋭度が向上して該放射線画像の診
断適性が向上する。以上の処理によるメリットは、前記
特開昭５５―１６３４７２号公報にさらに詳しく記載さ
れている。
【００３１】また、１回目、２回目のフィルタリング処
理における強調係数β1 、β2 を上述のように設定する
と、感熱記録装置28のＸ、Ｙ方向のレスポンスの差が補
償されて、この感熱記録装置28で記録される画像は、主
に空間周波数αにおいてＸ、Ｙ方向の鮮鋭度が互いに略
等しいものとなる。
理における強調係数β1 、β2 を上述のように設定する
と、感熱記録装置28のＸ、Ｙ方向のレスポンスの差が補
償されて、この感熱記録装置28で記録される画像は、主
に空間周波数αにおいてＸ、Ｙ方向の鮮鋭度が互いに略
等しいものとなる。
【００３２】なお画像処理装置27のメモリ24には、主に
鮮鋭度を揃えたい空間周波数α毎の強調係数β1 、β2
の組合わせを記憶させておき、この空間周波数αに応じ
てそこから強調係数β1 、β2 を読み出して、フィルタ
リング処理部25において設定させるようにすると便利で
ある。
鮮鋭度を揃えたい空間周波数α毎の強調係数β1 、β2
の組合わせを記憶させておき、この空間周波数αに応じ
てそこから強調係数β1 、β2 を読み出して、フィルタ
リング処理部25において設定させるようにすると便利で
ある。
【００３３】以上の実施例では、１次元方向についての
フィルタリング処理を２回行なっているが、それと等価
の２次元フィルタリング処理を１回行なって、同様の効
果を得ることも可能である。以下、そのような２次元フ
ィルタリング処理の例を説明する。
フィルタリング処理を２回行なっているが、それと等価
の２次元フィルタリング処理を１回行なって、同様の効
果を得ることも可能である。以下、そのような２次元フ
ィルタリング処理の例を説明する。
【００３４】前述した（２）式のフィルタリング処理の
係数が、Ｘ方向については強調係数をβ1 とすることに
よりａ、ｂ、ｃとなり、Ｙ方向については強調係数をβ
2 とすることによりａ’、ｂ’、ｃ’となるものとす
る。
係数が、Ｘ方向については強調係数をβ1 とすることに
よりａ、ｂ、ｃとなり、Ｙ方向については強調係数をβ
2 とすることによりａ’、ｂ’、ｃ’となるものとす
る。
【００３５】この場合、まず原画像データに対してＸ方
向に強調係数β1 で前記のフィルタリング処理を行な
い、その処理後の画像データに対してＹ方向に強調係数
β2 で同様のフィルタリング処理を行なうのと等価な作
用を果たす２次元フィルタは、次の行列式で表される。
向に強調係数β1 で前記のフィルタリング処理を行な
い、その処理後の画像データに対してＹ方向に強調係数
β2 で同様のフィルタリング処理を行なうのと等価な作
用を果たす２次元フィルタは、次の行列式で表される。
【００３６】
【数１】
【００３７】上記の２次元フィルタの各要素を３×３画
素の原画像データに掛けて足すことにより、中央位置の
画素に関する変換画像データが求められる。すなわち、
３×３画素の原画像データを例えば ｘ1 ｘ2 ｘ3 ｘ1 ’ ｘ2 ’ ｘ3 ’ ｘ1 ” ｘ2 ” ｘ3 ” としたとき、中央位置の画素の原画像データｘ2 ’は、 Ｑ’＝ａ’ａｘ1 ＋ａ’ｂｘ2 ＋ａ’ｃｘ3 +b'ax1 ’＋ｂ’ｂｘ2 ’＋ｂ’ｃｘ3 ’ ＋ｃ’ａｘ1 ”＋ｃ’ｂｘ2 ”＋ｃ’ｃｘ3 ” なる画像データＱ’に変換される。
素の原画像データに掛けて足すことにより、中央位置の
画素に関する変換画像データが求められる。すなわち、
３×３画素の原画像データを例えば ｘ1 ｘ2 ｘ3 ｘ1 ’ ｘ2 ’ ｘ3 ’ ｘ1 ” ｘ2 ” ｘ3 ” としたとき、中央位置の画素の原画像データｘ2 ’は、 Ｑ’＝ａ’ａｘ1 ＋ａ’ｂｘ2 ＋ａ’ｃｘ3 +b'ax1 ’＋ｂ’ｂｘ2 ’＋ｂ’ｃｘ3 ’ ＋ｃ’ａｘ1 ”＋ｃ’ｂｘ2 ”＋ｃ’ｃｘ3 ” なる画像データＱ’に変換される。
【００３８】また以上の各実施例では、２次元フィルタ
リング処理をＸ方向、Ｙ方向で互いに特性を変えて行な
うことにより、感熱記録装置28のＸ、Ｙ方向のレスポン
スの差を補償するようにしているが、このフィルタリン
グ処理を上記レスポンスの差に応じてＸ、Ｙ方向のうち
の一方向のみについて行なうことにより、感熱記録装置
28のＸ、Ｙ方向のレスポンスの差を補償することも可能
である。
リング処理をＸ方向、Ｙ方向で互いに特性を変えて行な
うことにより、感熱記録装置28のＸ、Ｙ方向のレスポン
スの差を補償するようにしているが、このフィルタリン
グ処理を上記レスポンスの差に応じてＸ、Ｙ方向のうち
の一方向のみについて行なうことにより、感熱記録装置
28のＸ、Ｙ方向のレスポンスの差を補償することも可能
である。
【図面の簡単な説明】
【図１】本発明の方法を実施する装置の一例を示す概略
構成図
構成図
【図２】本発明に関わる、画像再生装置のＸ、Ｙ各方向
のレスポンス特性を示すグラフ
のレスポンス特性を示すグラフ
【図３】画像再生装置のＸ、Ｙ各方向のレスポンス特性
の差を補償する周波数強調係数βの値を説明するグラフ
の差を補償する周波数強調係数βの値を説明するグラフ
10 蓄積性蛍光体シート
11 シート搬送手段
12 励起光源
13 レーザビーム
14 回転多面鏡
15 輝尽発光光
16 集光体
17 フォトマルチプライヤー
23 バッファメモリ
24 メモリ
25 フィルタリング処理部
27 画像処理装置
28 感熱記録装置
29 ハードコピー
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 昭56−75140（ＪＰ，Ａ)
特開 平４−90679（ＪＰ，Ａ)
特開 昭62−252260（ＪＰ，Ａ)
特開 昭55−50180（ＪＰ，Ａ)
特開 昭63−103227（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G03B 42/02 - 42/04
H04N 1/40
Claims (8)
- 【請求項１】 Ｘ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、各画素
についての原画像データをＱorg 、この画素に関する超
低空間周波数に対応するボケマスクデータをＱus、強調
係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけて、 この演算処理後の画像データを、Ｘ、Ｙ走査により画像
を再生する画像再生装置に送る画像処理方法において、 前記演算処理を、前記画像再生装置のＸ、Ｙ方向のレス
ポンスの差に応じて、該Ｘ、Ｙ方向について互いに異な
る特性の下に行なうことにより、この画像再生装置で再
生される画像のＸ、Ｙ方向の最終的な鮮鋭度を互いに略
等しく揃えることを特徴とする画像処理方法。 - 【請求項２】 前記演算処理をまずＸ、Ｙの一方向につ
いて強調係数をβ1に設定して行ない、 次いでこの演算処理後の画像を原画像データとして、
Ｘ、Ｙの他方向について強調係数をβ2 （β1 ≠β2 ）
に設定して前記演算処理を行なうことを特徴とする請求
項１記載の画像処理方法。 - 【請求項３】 前記演算処理をまずＸ、Ｙの一方向につ
いて強調係数をβ1に設定して行ない、次いでこの演算
処理後の画像を原画像データとして、Ｘ、Ｙの他方向に
ついて強調係数をβ2 （β1 ≠β2 ）に設定して前記演
算処理を行なうことと等価の２次元フィルタリング処理
を、前記原画像データに対してかけることを特徴とする
請求項１記載の画像処理方法。 - 【請求項４】 Ｘ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、各画素
についての原画像データをＱorg 、この画素に関する超
低空間周波数に対応するボケマスクデータをＱus、強調
係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけて、 この演算処理後の画像データを、Ｘ、Ｙ走査により画像
を再生する画像再生装置に送る画像処理方法において、 前記演算処理を、前記画像再生装置のＸ、Ｙ方向のレス
ポンスの差に応じて、該Ｘ、Ｙ方向のうちの一方向のみ
について行なうことにより、この画像再生装置で再生さ
れる画像のＸ、Ｙ方向の最終的な鮮鋭度を互いに略等し
く揃えることを特徴とする画像処理方法。 - 【請求項５】 Ｘ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、各画素
についての原画像データをＱorg 、この画素に関する超
低空間周波数に対応するボケマスクデータをＱus、強調
係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけて、 この演算処理後の画像データを、Ｘ、Ｙ走査により画像
を再生する画像再生装置に送る画像処理装置において、 前記演算処理を、前記Ｘ、Ｙ方向について互いに異なる
特性の下に行ない得ることを特徴とする画像処理装置。 - 【請求項６】 前記演算処理をまずＸ、Ｙの一方向につ
いて強調係数をβ1に設定して行ない、次いでこの演算
処理後の画像を原画像データとして、Ｘ、Ｙの他方向に
ついて強調係数をβ2 （β1 ≠β2 ）に設定して前記演
算処理を行なうことを特徴とする請求項５記載の画像処
理装置。 - 【請求項７】 前記演算処理をまずＸ、Ｙの一方向につ
いて強調係数をβ1に設定して行ない、次いでこの演算
処理後の画像を原画像データとして、Ｘ、Ｙの他方向に
ついて強調係数をβ2 （β1 ≠β2 ）に設定して前記演
算処理を行なうことと等価の２次元フィルタリング処理
を、前記原画像データに対してかけることを特徴とする
請求項５記載の画像処理装置。 - 【請求項８】 Ｘ、Ｙ方向に並ぶ画素のそれぞれについ
ての画素値を示しているデジタル画像データを、各画素
についての原画像データをＱorg 、この画素に関する超
低空間周波数に対応するボケマスクデータをＱus、強調
係数をβとしたとき、 Ｑ’＝Ｑorg ＋β（Ｑorg −Ｑus） なる演算処理にかけて、 この演算処理後の画像データを、Ｘ、Ｙ走査により画像
を再生する画像再生装置に送る画像処理装置において、 前記演算処理を、前記Ｘ、Ｙ方向のうちの一方向のみに
ついて行ない得ることを特徴とする画像処理装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP04518895A JP3399684B2 (ja) | 1995-03-06 | 1995-03-06 | 画像処理方法および装置 |
| US08/611,488 US5857037A (en) | 1995-03-06 | 1996-03-05 | Image processing method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP04518895A JP3399684B2 (ja) | 1995-03-06 | 1995-03-06 | 画像処理方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08240870A JPH08240870A (ja) | 1996-09-17 |
| JP3399684B2 true JP3399684B2 (ja) | 2003-04-21 |
Family
ID=12712301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP04518895A Expired - Fee Related JP3399684B2 (ja) | 1995-03-06 | 1995-03-06 | 画像処理方法および装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5857037A (ja) |
| JP (1) | JP3399684B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10257328A (ja) * | 1997-03-10 | 1998-09-25 | Ricoh Co Ltd | 画像処理装置 |
| EP1217571B2 (en) * | 2000-12-21 | 2014-03-12 | Datalogic IP TECH S.r.l. | Method and device for focusing an electrical signal representative of an optical code |
| IE87469B1 (en) | 2016-10-06 | 2024-01-03 | Google Llc | Image processing neural networks with separable convolutional layers |
Family Cites Families (2)
| 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 |
| JPS61230571A (ja) * | 1985-04-05 | 1986-10-14 | Dainippon Screen Mfg Co Ltd | 画像走査装置 |
-
1995
- 1995-03-06 JP JP04518895A patent/JP3399684B2/ja not_active Expired - Fee Related
-
1996
- 1996-03-05 US US08/611,488 patent/US5857037A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US5857037A (en) | 1999-01-05 |
| JPH08240870A (ja) | 1996-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH03222577A (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| US6614044B2 (en) | Image signal generating method, apparatus and program | |
| JP3399684B2 (ja) | 画像処理方法および装置 | |
| JP2532940B2 (ja) | 偽画像信号の検出方法 | |
| JPH0617984B2 (ja) | 放射線画像情報読取方法 | |
| US7366358B2 (en) | Method of and system for generating image signal | |
| JP2952428B2 (ja) | 放射線画像のエネルギーサブトラクション方法および装置 | |
| JPS61170178A (ja) | 放射線画像情報の読取条件決定方法 | |
| JPH01238650A (ja) | 放射線画像の被写体認識方法 | |
| JPH0743765B2 (ja) | 放射線画像処理方法および装置 | |
| JP3165530B2 (ja) | 放射線画像読取条件決定方法および装置 | |
| US20020186873A1 (en) | Image reproduction system | |
| JP2867096B2 (ja) | 画像平滑化処理方法 | |
| JPH0314372A (ja) | 放射線画像の時間差サブトラクション方法および放射線画像処理装置 | |
| JP3241488B2 (ja) | 画像処理方法 | |
| JP3165461B2 (ja) | 画像処理方法、画像読取方法および画像再生方法 | |
| JPH0786925B2 (ja) | 放射線画像処理方法および装置 | |
| JPH06235982A (ja) | 放射線画像処理方法 | |
| JPH01194683A (ja) | 放射線画像処理方法 | |
| JPH02272681A (ja) | 放射線画像処理装置 | |
| JPH0743767B2 (ja) | 放射線画像処理方法および装置 | |
| JPH0743768B2 (ja) | Ｘ線画像処理方法および装置 | |
| JPH06304160A (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: 20030204 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| 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: 20080221 Year of fee payment: 5 |
|
| 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: 20080221 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090221 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100221 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100221 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110221 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |