JPH1063837A - 画像処理方法および装置 - Google Patents
画像処理方法および装置Info
- Publication number
- JPH1063837A JPH1063837A JP8217333A JP21733396A JPH1063837A JP H1063837 A JPH1063837 A JP H1063837A JP 8217333 A JP8217333 A JP 8217333A JP 21733396 A JP21733396 A JP 21733396A JP H1063837 A JPH1063837 A JP H1063837A
- Authority
- JP
- Japan
- Prior art keywords
- image signal
- low
- pixels
- original image
- filter
- 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
- 238000000034 method Methods 0.000 title claims description 71
- 238000001914 filtration Methods 0.000 claims abstract description 69
- 238000003672 processing method Methods 0.000 claims abstract description 16
- 238000007906 compression Methods 0.000 claims description 8
- 230000006835 compression Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 19
- 230000000694 effects Effects 0.000 description 6
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 5
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 230000004807 localization Effects 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Picture Signal Circuits (AREA)
Abstract
(57)【要約】
【課題】 画像信号に対して所定のフィルタリング処理
を施すことにより非鮮鋭マスク画像信号を作成し、その
非鮮鋭マスク画像信号を用いて周波数強調などを目的と
する非線形処理を行う画像処理方法および装置におい
て、画像信号の読取密度に拘わらず同じ周波数特性の非
鮮鋭マスク画像信号を作成する。 【解決手段】 読取密度情報入力手段６により画像信号
の読取密度を得、その情報に基づいてフィルタ係数選択
手段５がフィルタ係数リスト７からフィルタ係数を選択
し、選択されたフィルタ係数のフィルタによってフィル
タリング手段２がフィルタリング処理を行う。
を施すことにより非鮮鋭マスク画像信号を作成し、その
非鮮鋭マスク画像信号を用いて周波数強調などを目的と
する非線形処理を行う画像処理方法および装置におい
て、画像信号の読取密度に拘わらず同じ周波数特性の非
鮮鋭マスク画像信号を作成する。 【解決手段】 読取密度情報入力手段６により画像信号
の読取密度を得、その情報に基づいてフィルタ係数選択
手段５がフィルタ係数リスト７からフィルタ係数を選択
し、選択されたフィルタ係数のフィルタによってフィル
タリング手段２がフィルタリング処理を行う。
Description
【０００１】
【発明の属する技術分野】本発明は、画像信号に対して
所定のフィルタリング処理を施すことにより非鮮鋭マス
ク画像信号を作成し、その非鮮鋭マスク画像信号を用い
て周波数強調などを目的とする非線形処理を行う画像処
理方法および装置に関し、特に詳しくは、非鮮鋭マスク
画像信号の作成に使用されるフィルタの特性に関するも
のである。
所定のフィルタリング処理を施すことにより非鮮鋭マス
ク画像信号を作成し、その非鮮鋭マスク画像信号を用い
て周波数強調などを目的とする非線形処理を行う画像処
理方法および装置に関し、特に詳しくは、非鮮鋭マスク
画像信号の作成に使用されるフィルタの特性に関するも
のである。
【０００２】
【従来の技術】従来、本出願人により、非鮮鋭マスク画
像信号（以下、ボケ画像信号という）を用いて周波数強
調処理あるいはダイナミックレンジ圧縮処理などを行っ
て放射線画像の診断性能を向上させる数々の画像処理方
法および装置が提案されている（特開昭55-163472 号、
同55-87953号、特開平3-222577号、特願平8-172498号、
同8-182155号、同8-182156号など）。ここで、ボケ画像
信号とは、画素数は原画像信号と同じであるが原画像信
号よりも鮮鋭度が低い画像を表す画像信号であり、ボケ
画像信号の作成は、原画像信号の画素に対して所定間隔
ごとに所定のフィルタリング処理を施すことによって画
素を間引きし、その後間引きした数分の画素を所定の補
間方法によって補間することにより行われる。
像信号（以下、ボケ画像信号という）を用いて周波数強
調処理あるいはダイナミックレンジ圧縮処理などを行っ
て放射線画像の診断性能を向上させる数々の画像処理方
法および装置が提案されている（特開昭55-163472 号、
同55-87953号、特開平3-222577号、特願平8-172498号、
同8-182155号、同8-182156号など）。ここで、ボケ画像
信号とは、画素数は原画像信号と同じであるが原画像信
号よりも鮮鋭度が低い画像を表す画像信号であり、ボケ
画像信号の作成は、原画像信号の画素に対して所定間隔
ごとに所定のフィルタリング処理を施すことによって画
素を間引きし、その後間引きした数分の画素を所定の補
間方法によって補間することにより行われる。
【０００３】上記所定のフィルタリング処理としては、
ローパスフィルタにより原画像信号の高周波成分を取り
除くような処理、具体的にはフィルタ内の画素値の平均
値や、荷重平均値を求めるといった処理が行われてい
る。このフィルタリング処理では、原画像信号に対して
フィルタリング処理を施して得られた画素数の少ない信
号に対してさらにこのフィルタリング処理を施して、フ
ィルタリングの各段階で得られた画素数の少ない画像信
号をそれぞれ原画像信号と同じ画素数となるように補間
し、複数の異なるボケ画像信号を作成することが行われ
ている。
ローパスフィルタにより原画像信号の高周波成分を取り
除くような処理、具体的にはフィルタ内の画素値の平均
値や、荷重平均値を求めるといった処理が行われてい
る。このフィルタリング処理では、原画像信号に対して
フィルタリング処理を施して得られた画素数の少ない信
号に対してさらにこのフィルタリング処理を施して、フ
ィルタリングの各段階で得られた画素数の少ない画像信
号をそれぞれ原画像信号と同じ画素数となるように補間
し、複数の異なるボケ画像信号を作成することが行われ
ている。
【０００４】このようにして作成されたボケ画像信号
は、原画像信号の所定の周波数帯域の成分を表すものと
なる。また、異なるボケ画像信号同士の差分をとること
により、さらに限られた周波数帯域の周波数成分を表す
信号を得ることもできる。これは、例えば周波数強調処
理やダイナミックレンジ圧縮処理などにおいて、原画像
信号の特定の周波数帯域の成分を操作する際に用いられ
ており、例えば、特定の周波数帯域の信号を必要に応じ
て抑制して原画像信号に加算することによって周波数強
調を行うなどの非線形処理が行われている。
は、原画像信号の所定の周波数帯域の成分を表すものと
なる。また、異なるボケ画像信号同士の差分をとること
により、さらに限られた周波数帯域の周波数成分を表す
信号を得ることもできる。これは、例えば周波数強調処
理やダイナミックレンジ圧縮処理などにおいて、原画像
信号の特定の周波数帯域の成分を操作する際に用いられ
ており、例えば、特定の周波数帯域の信号を必要に応じ
て抑制して原画像信号に加算することによって周波数強
調を行うなどの非線形処理が行われている。
【０００５】なお、ボケ画像信号は、上記のように原画
像信号に基づいて作成されるものであるが、この原画像
信号は、読取装置により原画像信号を所定の読取密度で
読み取ってデジタル信号としたものである。一般に、デ
ジタル化された画像信号を例えばプリント出力として再
生する場合、読取密度によって決定されるある周波数
（ナイキスト周波数）以下の周波数成分は正しく再生さ
れることが知られている。つまり、読取密度は再生時に
必要となる画質レベルを考慮して決定されるものである
ため、必ずしも一定であるとは限らないものである。
像信号に基づいて作成されるものであるが、この原画像
信号は、読取装置により原画像信号を所定の読取密度で
読み取ってデジタル信号としたものである。一般に、デ
ジタル化された画像信号を例えばプリント出力として再
生する場合、読取密度によって決定されるある周波数
（ナイキスト周波数）以下の周波数成分は正しく再生さ
れることが知られている。つまり、読取密度は再生時に
必要となる画質レベルを考慮して決定されるものである
ため、必ずしも一定であるとは限らないものである。
【０００６】例えば、上記引用例に示されている放射線
画像読取再生システムでは、蓄積性蛍光体シートに記録
された人体の放射線画像をレーザビーム走査によりデジ
タル画像信号として読み取っているが、読取密度は蓄積
性蛍光体シートの大きさによって異なるものとしてい
る。
画像読取再生システムでは、蓄積性蛍光体シートに記録
された人体の放射線画像をレーザビーム走査によりデジ
タル画像信号として読み取っているが、読取密度は蓄積
性蛍光体シートの大きさによって異なるものとしてい
る。
【０００７】
【発明が解決しようとする課題】ここで、読取密度が異
なる画像信号、すなわちナイキスト周波数が異なる画像
信号に対し、同一のローパスフィルタによるフィルタリ
ング処理、および同一の補間方法による補間処理を施し
た場合、得られるボケ画像信号の周波数特性は読取密度
によって異なることになる。これにより、例えば１つの
原画像を２種類の読取密度で読み取って２種類の原画像
信号を得た場合、同じボケ画像信号を用いて周波数強調
処理あるいはダイナミックレンジ圧縮処理を行っても、
強調される周波数帯域あるいは圧縮される周波数帯域は
２種類の原画像信号でそれぞれ異なってしまうという問
題が生ずる。
なる画像信号、すなわちナイキスト周波数が異なる画像
信号に対し、同一のローパスフィルタによるフィルタリ
ング処理、および同一の補間方法による補間処理を施し
た場合、得られるボケ画像信号の周波数特性は読取密度
によって異なることになる。これにより、例えば１つの
原画像を２種類の読取密度で読み取って２種類の原画像
信号を得た場合、同じボケ画像信号を用いて周波数強調
処理あるいはダイナミックレンジ圧縮処理を行っても、
強調される周波数帯域あるいは圧縮される周波数帯域は
２種類の原画像信号でそれぞれ異なってしまうという問
題が生ずる。
【０００８】本発明は、上記問題に鑑みて、読取密度に
拘わらず同じ周波数特性のボケ画像信号を作成すること
ができ、これにより所望の非線形処理、例えば特定の周
波数成分を強調するような処理を、常に同じように行う
ことができる画像処理方法および装置を提供することを
目的とするものである。
拘わらず同じ周波数特性のボケ画像信号を作成すること
ができ、これにより所望の非線形処理、例えば特定の周
波数成分を強調するような処理を、常に同じように行う
ことができる画像処理方法および装置を提供することを
目的とするものである。
【０００９】
【課題を解決するための手段】本発明の画像処理方法
は、原画像を所定の読取密度で読み取って複数の画素に
より構成される原画像信号を得、該原画像信号に対して
所定のフィルタにより所定のフィルタリング処理を施し
て前記原画像信号よりも画素数が少ない低解像度画像信
号を作成し、該低解像度画像信号に対して前記フィルタ
リング処理をさらに施して、前記低解像度画像信号より
もさらに画素数が少ない低解像度画像信号を作成し、前
記低解像度画像信号に対するフィルタリング処理を繰り
返し行うことによって、互いに解像度が異なる複数の低
解像度画像信号を作成し、該各低解像度画像信号を、該
低解像度画像信号の画素数が前記原画像信号の画素数と
同じになるようにそれぞれ補間拡大することにより、前
記原画像信号の複数の異なるボケ画像信号を作成し、そ
のボケ画像信号および前記原画像信号に基づいて所定の
非線形処理を行う画像処理方法において、前記所定のフ
ィルタとして、前記読取密度に基づいて定められたフィ
ルタ係数を有するフィルタを使用することを特徴とする
ものである。
は、原画像を所定の読取密度で読み取って複数の画素に
より構成される原画像信号を得、該原画像信号に対して
所定のフィルタにより所定のフィルタリング処理を施し
て前記原画像信号よりも画素数が少ない低解像度画像信
号を作成し、該低解像度画像信号に対して前記フィルタ
リング処理をさらに施して、前記低解像度画像信号より
もさらに画素数が少ない低解像度画像信号を作成し、前
記低解像度画像信号に対するフィルタリング処理を繰り
返し行うことによって、互いに解像度が異なる複数の低
解像度画像信号を作成し、該各低解像度画像信号を、該
低解像度画像信号の画素数が前記原画像信号の画素数と
同じになるようにそれぞれ補間拡大することにより、前
記原画像信号の複数の異なるボケ画像信号を作成し、そ
のボケ画像信号および前記原画像信号に基づいて所定の
非線形処理を行う画像処理方法において、前記所定のフ
ィルタとして、前記読取密度に基づいて定められたフィ
ルタ係数を有するフィルタを使用することを特徴とする
ものである。
【００１０】ここで、前記非線形処理とは、具体的に
は、原画像信号に含まれる特定の周波数成分を強調する
周波数強調処理や、原画像の最高濃度と最低濃度との差
すなわちダイナミックレンジを狭めるように高濃度域も
しくは低濃度域あるいは画像全体のコントラストを下げ
るダイナミックレンジ圧縮処理などが挙げられる。ま
た、上記方法により作成されるボケ画像信号は、必ずし
もそれを用いて非線形処理を行うためだけに作成される
ものではなく、原画像よりも鮮鋭度の低い画像としてそ
れ自体を使用することも可能である。
は、原画像信号に含まれる特定の周波数成分を強調する
周波数強調処理や、原画像の最高濃度と最低濃度との差
すなわちダイナミックレンジを狭めるように高濃度域も
しくは低濃度域あるいは画像全体のコントラストを下げ
るダイナミックレンジ圧縮処理などが挙げられる。ま
た、上記方法により作成されるボケ画像信号は、必ずし
もそれを用いて非線形処理を行うためだけに作成される
ものではなく、原画像よりも鮮鋭度の低い画像としてそ
れ自体を使用することも可能である。
【００１１】なお、「所定のフィルタ」とは、例えば上
述のようにフィルタ内の画素値の平均値を求めるような
ローパスフィルタであり、所定の周波数以上の高周波成
分を除去するようなフィルタである。また「所定のフィ
ルタリング処理」とは、原画像信号を構成する画素の、
例えば１画素おきに、フィルタによる演算を行うことで
あり、フィルタリングの結果得られる画像信号の画素数
が、原画像信号の画素数よりも少なくなるように、すな
わち画素を間引くような処理である。また、「補間拡
大」の処理の演算方法としては、一般に知られている、
例えばＢスプラインによる方法など、種々の方法を適用
することができる。
述のようにフィルタ内の画素値の平均値を求めるような
ローパスフィルタであり、所定の周波数以上の高周波成
分を除去するようなフィルタである。また「所定のフィ
ルタリング処理」とは、原画像信号を構成する画素の、
例えば１画素おきに、フィルタによる演算を行うことで
あり、フィルタリングの結果得られる画像信号の画素数
が、原画像信号の画素数よりも少なくなるように、すな
わち画素を間引くような処理である。また、「補間拡
大」の処理の演算方法としては、一般に知られている、
例えばＢスプラインによる方法など、種々の方法を適用
することができる。
【００１２】さらに、本発明の画像処理装置は、上記画
像処理方法に基づいて画像処理を実施する装置であり、
原画像を所定の読取密度で読み取って得られる複数の画
素により構成される原画像信号に対して所定のフィルタ
により所定のフィルタリング処理を施して前記原画像信
号よりも画素数が少ない低解像度画像信号を作成し、該
低解像度画像信号に対して前記フィルタリング処理をさ
らに施して、前記低解像度画像信号よりもさらに画素数
が少ない低解像度画像信号を作成し、前記低解像度画像
信号に対するフィルタリング処理を繰り返し行うことに
よって、互いに解像度が異なる複数の低解像度画像信号
を作成するフィルタリング手段と、前記各低解像度画像
信号を、該低解像度画像信号の画素数が前記原画像信号
の画素数と同じになるようにそれぞれ補間拡大すること
により、前記原画像信号の複数の異なる非鮮鋭マスク画
像信号を作成する補間拡大手段と、前記非鮮鋭マスク画
像信号および前記原画像信号に基づいて所定の非線形処
理を行う非線形処理手段と、前記読取密度に関する情報
を入力するための読取密度情報入力手段と、前記情報と
して入力された読取密度に基づいて、前記所定のフィル
タのフィルタ係数を予め保持された複数のフィルタ係数
の中から選択するフィルタ係数選択手段を有することを
特徴とするものである。
像処理方法に基づいて画像処理を実施する装置であり、
原画像を所定の読取密度で読み取って得られる複数の画
素により構成される原画像信号に対して所定のフィルタ
により所定のフィルタリング処理を施して前記原画像信
号よりも画素数が少ない低解像度画像信号を作成し、該
低解像度画像信号に対して前記フィルタリング処理をさ
らに施して、前記低解像度画像信号よりもさらに画素数
が少ない低解像度画像信号を作成し、前記低解像度画像
信号に対するフィルタリング処理を繰り返し行うことに
よって、互いに解像度が異なる複数の低解像度画像信号
を作成するフィルタリング手段と、前記各低解像度画像
信号を、該低解像度画像信号の画素数が前記原画像信号
の画素数と同じになるようにそれぞれ補間拡大すること
により、前記原画像信号の複数の異なる非鮮鋭マスク画
像信号を作成する補間拡大手段と、前記非鮮鋭マスク画
像信号および前記原画像信号に基づいて所定の非線形処
理を行う非線形処理手段と、前記読取密度に関する情報
を入力するための読取密度情報入力手段と、前記情報と
して入力された読取密度に基づいて、前記所定のフィル
タのフィルタ係数を予め保持された複数のフィルタ係数
の中から選択するフィルタ係数選択手段を有することを
特徴とするものである。
【００１３】
【発明の効果】本発明の画像処理方法および装置によれ
ば、画像信号に対して所定のフィルタリング処理を施す
ことによりボケ画像信号を作成し、そのボケ画像信号を
用いて周波数強調などを目的とする非線形処理を行う場
合に、ボケ画像信号の作成に使用するフィルタの特性
を、原画像信号を読み取った際の読取密度に応じて設定
することにより、読取密度が異なっても、同じ周波数特
性のボケ画像信号を作成することができ、これにより読
取密度の異なる原画像信号に対して同様の非線形処理を
施すことができ、原画像信号の読取密度に影響されない
一定の画像処理結果が得られるようになる。
ば、画像信号に対して所定のフィルタリング処理を施す
ことによりボケ画像信号を作成し、そのボケ画像信号を
用いて周波数強調などを目的とする非線形処理を行う場
合に、ボケ画像信号の作成に使用するフィルタの特性
を、原画像信号を読み取った際の読取密度に応じて設定
することにより、読取密度が異なっても、同じ周波数特
性のボケ画像信号を作成することができ、これにより読
取密度の異なる原画像信号に対して同様の非線形処理を
施すことができ、原画像信号の読取密度に影響されない
一定の画像処理結果が得られるようになる。
【００１４】
【発明の実施の形態】以下、本発明の画像処理方法およ
び装置の一実施の形態について、図面を参照して詳細に
説明する。以下に示す画像処理装置は、蓄積性蛍光体シ
ートに記録された人体の放射線画像を読み取って得た画
像信号に対して、その画像が診断に適した画像となるよ
うに、ボケ画像信号を使用して周波数強調処理あるいは
ダイナミックレンジ圧縮処理などを施すものであり、処
理された画像信号は主としてフィルムに記録され、診断
に用いられる。
び装置の一実施の形態について、図面を参照して詳細に
説明する。以下に示す画像処理装置は、蓄積性蛍光体シ
ートに記録された人体の放射線画像を読み取って得た画
像信号に対して、その画像が診断に適した画像となるよ
うに、ボケ画像信号を使用して周波数強調処理あるいは
ダイナミックレンジ圧縮処理などを施すものであり、処
理された画像信号は主としてフィルムに記録され、診断
に用いられる。
【００１５】図１はこの画像処理装置の概略を表す図で
ある。画像処理装置１は、ボケ画像信号を作成するため
の手段として、フィルタリング処理を行うフィルタリン
グ手段２と、補間係数リスト８として記憶された補間係
数を使用して補間拡大処理を行う補間拡大手段３とを有
する。さらに、画像処理の目的に応じて、例えば周波数
強調処理などを行う手段として、非線形処理手段４を有
する。この際、画像処理装置１はフィルタ係数リスト７
として数種類の係数を記憶しており、フィルタ係数選択
手段５がその中の１種類を選択し、フィルタリング手段
２はその選択されたフィルタ係数を使用する。この選択
は、処理する画像信号の読取密度に基づいて行われるも
のであり、この読取密度情報を得るための手段として読
取密度情報入力手段６が備えられている。
ある。画像処理装置１は、ボケ画像信号を作成するため
の手段として、フィルタリング処理を行うフィルタリン
グ手段２と、補間係数リスト８として記憶された補間係
数を使用して補間拡大処理を行う補間拡大手段３とを有
する。さらに、画像処理の目的に応じて、例えば周波数
強調処理などを行う手段として、非線形処理手段４を有
する。この際、画像処理装置１はフィルタ係数リスト７
として数種類の係数を記憶しており、フィルタ係数選択
手段５がその中の１種類を選択し、フィルタリング手段
２はその選択されたフィルタ係数を使用する。この選択
は、処理する画像信号の読取密度に基づいて行われるも
のであり、この読取密度情報を得るための手段として読
取密度情報入力手段６が備えられている。
【００１６】読取密度情報入力手段６による入力は、ユ
ーザがキーボードから数値として入力してもよいし、操
作画面に数種類の密度を表示してユーザに選択させると
いった形態でもよい。あるいは、読取装置側で、読取密
度情報を画像信号に添付しておき、画像処理装置１が入
力された画像信号ごとに添付された情報を認識するよう
にしてもよく、画像処理装置１、詳しくはフィルタ係数
選択手段５が読取密度を認識することができれば、どの
ような形態であってもよい。
ーザがキーボードから数値として入力してもよいし、操
作画面に数種類の密度を表示してユーザに選択させると
いった形態でもよい。あるいは、読取装置側で、読取密
度情報を画像信号に添付しておき、画像処理装置１が入
力された画像信号ごとに添付された情報を認識するよう
にしてもよく、画像処理装置１、詳しくはフィルタ係数
選択手段５が読取密度を認識することができれば、どの
ような形態であってもよい。
【００１７】ここで、まずボケ画像信号の作成、すなわ
ちフィルタリング手段２および補間拡大手段３の処理に
ついて詳細に説明する。図２はボケ画像信号作成処理の
概要を示すブロック図である。図２に示されるように、
フィルタリング手段２は、まず原画像信号Ｓorgに対
し、原画像の画素のｘ方向およびｙ方向に対してフィル
タリング処理を施して原画像信号よりも解像度が低い第
１の低解像度画像信号Ｂ1を作成し、次に第１の低解像
度画像信号Ｂ1に対して同様のフィルタリング処理を施
して第１の低解像度画像信号よりもさらに解像度が低い
第２の低解像度画像信号Ｂ2を作成し、以降順次同様の
フィルタリング処理を重ねていくものである。そして、
補間拡大手段３は、前記フィルタリング処理の各段にお
いて得られる低解像度画像信号に対して、それぞれ補間
拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号
Ｓus1〜ＳusNを得るものである。
ちフィルタリング手段２および補間拡大手段３の処理に
ついて詳細に説明する。図２はボケ画像信号作成処理の
概要を示すブロック図である。図２に示されるように、
フィルタリング手段２は、まず原画像信号Ｓorgに対
し、原画像の画素のｘ方向およびｙ方向に対してフィル
タリング処理を施して原画像信号よりも解像度が低い第
１の低解像度画像信号Ｂ1を作成し、次に第１の低解像
度画像信号Ｂ1に対して同様のフィルタリング処理を施
して第１の低解像度画像信号よりもさらに解像度が低い
第２の低解像度画像信号Ｂ2を作成し、以降順次同様の
フィルタリング処理を重ねていくものである。そして、
補間拡大手段３は、前記フィルタリング処理の各段にお
いて得られる低解像度画像信号に対して、それぞれ補間
拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号
Ｓus1〜ＳusNを得るものである。
【００１８】本実施の形態においては、上記フィルタリ
ング処理のフィルタとして、一次元ガウス分布に略対応
したフィルタを使用する。すなわちフィルタのフィルタ
係数を、ガウス信号に関する下記の式（１）
ング処理のフィルタとして、一次元ガウス分布に略対応
したフィルタを使用する。すなわちフィルタのフィルタ
係数を、ガウス信号に関する下記の式（１）
【００１９】
【数１】
【００２０】にしたがって定める。これは、ガウス信号
は周波数空間および実空間の双方において、局在性がよ
いためである。
は周波数空間および実空間の双方において、局在性がよ
いためである。
【００２１】フィルタリング処理は、図３に示すよう
に、原画像信号Ｓorg に対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をｘ方向、ｙ方向に行うことによ
り、低解像度画像信号Ｂ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｂk （ｋ＝
１〜ｎ）は、それぞれ、画素数が原画像信号の１／２2k
の画像信号となる。
に、原画像信号Ｓorg に対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をｘ方向、ｙ方向に行うことによ
り、低解像度画像信号Ｂ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｂk （ｋ＝
１〜ｎ）は、それぞれ、画素数が原画像信号の１／２2k
の画像信号となる。
【００２２】図４は、このようにして生成された低解像
度画像信号Ｂk の周波数レスポンスを示すものであり、
横軸中に示されるＮ／２はナイキスト周波数を表してい
る。この図のように、低解像度画像信号Ｂkのレスポン
スはｋが大きいほど高周波成分が除去されたものとなっ
ている（なお、図４においてはｋ＝１〜３としてい
る）。
度画像信号Ｂk の周波数レスポンスを示すものであり、
横軸中に示されるＮ／２はナイキスト周波数を表してい
る。この図のように、低解像度画像信号Ｂkのレスポン
スはｋが大きいほど高周波成分が除去されたものとなっ
ている（なお、図４においてはｋ＝１〜３としてい
る）。
【００２３】次に、このようにして得られた低解像度画
像信号Ｂkに対して施される補間拡大処理について説明
する。補間演算の方法としては、Ｂスプラインによる方
法など種々の方法が挙げられるが、本実施の形態におい
ては、上記フィルタリング処理においてガウス信号に基
づくローパスフィルタを用いているため、補間演算につ
いてもガウス信号を用いるものとする。具体的には、下
記の式（２）
像信号Ｂkに対して施される補間拡大処理について説明
する。補間演算の方法としては、Ｂスプラインによる方
法など種々の方法が挙げられるが、本実施の形態におい
ては、上記フィルタリング処理においてガウス信号に基
づくローパスフィルタを用いているため、補間演算につ
いてもガウス信号を用いるものとする。具体的には、下
記の式（２）
【００２４】
【数２】
【００２５】において、σ＝２k-1 と近似したものを用
いる。
いる。
【００２６】低解像度画像信号Ｂ1を補間する際には、
ｋ＝１であるためσ＝１となる。この場合、補間処理を
行うためのフィルタは、図５に示すように５×１の一次
元フィルタとなる。この補間処理は、まず低解像度画像
信号Ｂ1に対して１画素おきに値が０の画素を１つずつ
補間することにより低解像度画像信号Ｂ1を原画像と同
一のサイズに拡大し、次に、この補間された低解像度画
像信号Ｂ1に対して上述した図５に示す一次元フィルタ
によりフィルタリング処理を施すことにより行われる。
ｋ＝１であるためσ＝１となる。この場合、補間処理を
行うためのフィルタは、図５に示すように５×１の一次
元フィルタとなる。この補間処理は、まず低解像度画像
信号Ｂ1に対して１画素おきに値が０の画素を１つずつ
補間することにより低解像度画像信号Ｂ1を原画像と同
一のサイズに拡大し、次に、この補間された低解像度画
像信号Ｂ1に対して上述した図５に示す一次元フィルタ
によりフィルタリング処理を施すことにより行われる。
【００２７】図６は、値が０の画素が補間された低解像
度画像信号Ｂ1を一次元的に示すものである。この図か
らわかるように、図５のフィルタは５×１のフィルタで
あるが、このフィルタによるフィルタリング処理は、実
質的には２×１のフィルタ（0.5，0.5）および３×１の
フィルタ（0.1，0.8，0.1）の２種類のフィルタによる
フィルタリング処理と等価なものである。以上のような
補間拡大処理により、原画像信号Ｓorgと同一画素数、
すなわち原画像と同一サイズのボケ画像信号Ｓus1が得
られる。
度画像信号Ｂ1を一次元的に示すものである。この図か
らわかるように、図５のフィルタは５×１のフィルタで
あるが、このフィルタによるフィルタリング処理は、実
質的には２×１のフィルタ（0.5，0.5）および３×１の
フィルタ（0.1，0.8，0.1）の２種類のフィルタによる
フィルタリング処理と等価なものである。以上のような
補間拡大処理により、原画像信号Ｓorgと同一画素数、
すなわち原画像と同一サイズのボケ画像信号Ｓus1が得
られる。
【００２８】次に、低解像度画像信号Ｂ2に対して補間
拡大処理が施される。低解像度画像信号Ｂ2の補間にお
いては、ｋ＝２であるため、σ＝２となる。この場合、
補間処理を行うためのフィルタは、図７に示すように１
１×１の一次元フィルタとなる。この補間処理は、まず
低解像度画像信号Ｂ2に対して１画素おきに値が０の画
素を３つずつ補間することにより低解像度画像信号Ｂ2
を原画像と同一のサイズに拡大し、次に、この補間され
た低解像度画像信号Ｂ2に対して上述した図７に示す一
次元フィルタによりフィルタリング処理を施すことによ
り行われる。
拡大処理が施される。低解像度画像信号Ｂ2の補間にお
いては、ｋ＝２であるため、σ＝２となる。この場合、
補間処理を行うためのフィルタは、図７に示すように１
１×１の一次元フィルタとなる。この補間処理は、まず
低解像度画像信号Ｂ2に対して１画素おきに値が０の画
素を３つずつ補間することにより低解像度画像信号Ｂ2
を原画像と同一のサイズに拡大し、次に、この補間され
た低解像度画像信号Ｂ2に対して上述した図７に示す一
次元フィルタによりフィルタリング処理を施すことによ
り行われる。
【００２９】図８は、値が０の画素が補間された低解像
度画像信号Ｂ2を一次元的に示すものである。この図か
らわかるように、図７のフィルタは１１×１のフィルタ
であるが、このフィルタによるフィルタリング処理は、
実質的には２×１のフィルタ（0.5，0.5）および３×１
のフィルタ（0.3，0.65，0.05）、（0.13，0.74，0.1
3）および（0.05，0.65，0.3）の４種類のフィルタによ
るフィルタリング処理と等価なものである。以上のよう
な補間拡大処理により、原画像信号Ｓorgと同一画素
数、すなわち原画像と同一サイズのボケ画像信号Ｓus2
が得られる。
度画像信号Ｂ2を一次元的に示すものである。この図か
らわかるように、図７のフィルタは１１×１のフィルタ
であるが、このフィルタによるフィルタリング処理は、
実質的には２×１のフィルタ（0.5，0.5）および３×１
のフィルタ（0.3，0.65，0.05）、（0.13，0.74，0.1
3）および（0.05，0.65，0.3）の４種類のフィルタによ
るフィルタリング処理と等価なものである。以上のよう
な補間拡大処理により、原画像信号Ｓorgと同一画素
数、すなわち原画像と同一サイズのボケ画像信号Ｓus2
が得られる。
【００３０】同様に、この補間拡大処理を全ての低解像
度画像信号Ｂkに対して行う。低解像度画像信号Ｂkを補
間する際には、上記式（２）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｂkの各
画素の間に値が０の画素を２k−１個ずつ補間すること
により、原画像と同一サイズに拡大し、この値が０の画
素が補間された低解像度画像信号Ｂkに対して３×２k−
１の長さのフィルタにより、フィルタリング処理を施す
ことにより補間拡大する。
度画像信号Ｂkに対して行う。低解像度画像信号Ｂkを補
間する際には、上記式（２）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｂkの各
画素の間に値が０の画素を２k−１個ずつ補間すること
により、原画像と同一サイズに拡大し、この値が０の画
素が補間された低解像度画像信号Ｂkに対して３×２k−
１の長さのフィルタにより、フィルタリング処理を施す
ことにより補間拡大する。
【００３１】ここで、この３×２k−１の長さのフィル
タによるフィルタリング処理を施すことは、２k周期で
長さが２または３のフィルタによりフィルタリング処理
を施すことと等価なものとなる。つまり、フィルタ長が
長くなっても、実質的には長さが２または３のフィルタ
によりフィルタリング処理を施していることと同様であ
るため、演算量はそれほど多くはならず、したがって、
高速なボケ画像信号作成処理を行うことが可能である。
タによるフィルタリング処理を施すことは、２k周期で
長さが２または３のフィルタによりフィルタリング処理
を施すことと等価なものとなる。つまり、フィルタ長が
長くなっても、実質的には長さが２または３のフィルタ
によりフィルタリング処理を施していることと同様であ
るため、演算量はそれほど多くはならず、したがって、
高速なボケ画像信号作成処理を行うことが可能である。
【００３２】なお、この補間拡大処理におけるフィルタ
リング処理のフィルタは、必ずしも一次元フィルタに限
られるものではなく、２次元フィルタにより行ってもよ
い。この場合も同様に、フィルタサイズが大きくなって
も、実質的には２×２または３×３のフィルタによるフ
ィルタリング処理と等価になるため、高速な演算処理を
行うことができる。また、ガウス信号のフィルタの代わ
りに、Ｂスプライン補間演算により上記補間拡大処理を
行っても良く、その他種々の補間演算方法を適用するこ
とができる。
リング処理のフィルタは、必ずしも一次元フィルタに限
られるものではなく、２次元フィルタにより行ってもよ
い。この場合も同様に、フィルタサイズが大きくなって
も、実質的には２×２または３×３のフィルタによるフ
ィルタリング処理と等価になるため、高速な演算処理を
行うことができる。また、ガウス信号のフィルタの代わ
りに、Ｂスプライン補間演算により上記補間拡大処理を
行っても良く、その他種々の補間演算方法を適用するこ
とができる。
【００３３】以上、ボケ画像信号の作成について説明し
たが、次に、本発明が解決しようとする問題点およびそ
の解決手段について、例を示して説明する。従来、上記
のようなボケ画像信号作成処理では、原画像信号が取得
される際の読取密度は常に一定であることが前提となっ
ており、ボケ画像信号を作成するにあたり、この読取密
度については何ら考慮がなされていなかった。しかしな
がら、実際の画像処理システムにおいては、種々の画像
信号が入力されるため、全ての画像信号の読取密度が同
じであるとは限らない。例えば、本実施の形態の放射線
画像処理システムでは、この読取密度は蓄積性蛍光体シ
ートの大きさによって異なっており、半切／大角サイズ
で５本／mm、４つ切サイズで6.7本／mm、６つ切サイズ
で10本／mmとなっている。
たが、次に、本発明が解決しようとする問題点およびそ
の解決手段について、例を示して説明する。従来、上記
のようなボケ画像信号作成処理では、原画像信号が取得
される際の読取密度は常に一定であることが前提となっ
ており、ボケ画像信号を作成するにあたり、この読取密
度については何ら考慮がなされていなかった。しかしな
がら、実際の画像処理システムにおいては、種々の画像
信号が入力されるため、全ての画像信号の読取密度が同
じであるとは限らない。例えば、本実施の形態の放射線
画像処理システムでは、この読取密度は蓄積性蛍光体シ
ートの大きさによって異なっており、半切／大角サイズ
で５本／mm、４つ切サイズで6.7本／mm、６つ切サイズ
で10本／mmとなっている。
【００３４】しかしながら、読取密度が異なる画像信
号、すなわちナイキスト周波数が異なる画像信号に対
し、同一のローパスフィルタによるフィルタリング処
理、および同一の補間方法による補間処理を施した場
合、得られるボケ画像信号の周波数特性は読取密度によ
って異なるものとなる。例えば、読取密度が10本／mm、
５本／mm、6.7本／mmのそれぞれの原画像信号に対し
て、図９に示す一次元フィルタを使用したフィルタリン
グ処理を施して作成したボケ画像信号の周波数特性を、
読取密度が10本／mmの場合を図１０に、５本／mmの場合
を図１１に、6.7本／mmの場合を図１２にそれぞれ示
す。
号、すなわちナイキスト周波数が異なる画像信号に対
し、同一のローパスフィルタによるフィルタリング処
理、および同一の補間方法による補間処理を施した場
合、得られるボケ画像信号の周波数特性は読取密度によ
って異なるものとなる。例えば、読取密度が10本／mm、
５本／mm、6.7本／mmのそれぞれの原画像信号に対し
て、図９に示す一次元フィルタを使用したフィルタリン
グ処理を施して作成したボケ画像信号の周波数特性を、
読取密度が10本／mmの場合を図１０に、５本／mmの場合
を図１１に、6.7本／mmの場合を図１２にそれぞれ示
す。
【００３５】ここで、上記のようなボケ画像信号作成処
理では、ｋ段目の低解像度画像信号Ｂkのナイキスト周
波数はｋ-１段目の低解像度画像信号Ｂk-1のナイキスト
周波数の１／２になる。そのため、読取密度が10本／mm
と５本／mmのようにナイキスト周波数が１／２の関係に
なっている場合、読取密度10本／mmのときのボケ画像信
号Ｓuskと読取密度５本／mmのときのボケ画像信号Ｓusk
-1の周波数特性は等しくなる。これは、図１０と図１１
の比較からも明らかである。これは、読取密度が10本／
mmの画像信号と５本／mmの画像信号からは、同一の周波
数特性を持つボケ画像信号を作成することができるとい
うことであり、したがって、同一のボケ画像信号を使用
して同一の非線形処理を行うことができるということを
表している。
理では、ｋ段目の低解像度画像信号Ｂkのナイキスト周
波数はｋ-１段目の低解像度画像信号Ｂk-1のナイキスト
周波数の１／２になる。そのため、読取密度が10本／mm
と５本／mmのようにナイキスト周波数が１／２の関係に
なっている場合、読取密度10本／mmのときのボケ画像信
号Ｓuskと読取密度５本／mmのときのボケ画像信号Ｓusk
-1の周波数特性は等しくなる。これは、図１０と図１１
の比較からも明らかである。これは、読取密度が10本／
mmの画像信号と５本／mmの画像信号からは、同一の周波
数特性を持つボケ画像信号を作成することができるとい
うことであり、したがって、同一のボケ画像信号を使用
して同一の非線形処理を行うことができるということを
表している。
【００３６】これに対し、図１２に示される読取密度6.
7本／mmのときのボケ画像信号の周波数特性は、図１０
および図１１とは全く異なるものであり、したがってこ
のようなボケ画像信号を用いて上記読取密度が10本／mm
と５本／mmの場合と同様な非線形処理を施しても、10本
／mmと５本／mmの場合と同様の効果、例えば特定の周波
数成分の強調などの効果は得られない。
7本／mmのときのボケ画像信号の周波数特性は、図１０
および図１１とは全く異なるものであり、したがってこ
のようなボケ画像信号を用いて上記読取密度が10本／mm
と５本／mmの場合と同様な非線形処理を施しても、10本
／mmと５本／mmの場合と同様の効果、例えば特定の周波
数成分の強調などの効果は得られない。
【００３７】そこで、本発明の画像処理方法および装置
では、図９のフィルタの他に、図１３に示すような別の
フィルタをも用意しておき、読取密度情報に基づいて、
読取密度が10本／mmと５本／mmの場合には図９のフィル
タを、読取密度が6.7本／mmの場合には図１３のフィル
タを選択的に使用する。図１４は、読取密度が6.7本／m
mの場合に図１３のフィルタを使用して作成したボケ画
像信号の周波数特性を示す図であるが、図11の読取密度
が５本／mmの場合に図９のフィルタを使用して作成した
ボケ画像信号の周波数特性とほぼ等しくなっていること
が、図１１と図１４の比較によりわかる。すなわち、読
取密度が6.7本／mmの場合でも、図１３のフィルタを使
用してボケ画像信号を作成すれば、読取密度が10本／mm
と５本／mmの場合と同様な非線形処理により、10本／mm
と５本／mmの場合と同様の効果を得ることが可能とな
る。
では、図９のフィルタの他に、図１３に示すような別の
フィルタをも用意しておき、読取密度情報に基づいて、
読取密度が10本／mmと５本／mmの場合には図９のフィル
タを、読取密度が6.7本／mmの場合には図１３のフィル
タを選択的に使用する。図１４は、読取密度が6.7本／m
mの場合に図１３のフィルタを使用して作成したボケ画
像信号の周波数特性を示す図であるが、図11の読取密度
が５本／mmの場合に図９のフィルタを使用して作成した
ボケ画像信号の周波数特性とほぼ等しくなっていること
が、図１１と図１４の比較によりわかる。すなわち、読
取密度が6.7本／mmの場合でも、図１３のフィルタを使
用してボケ画像信号を作成すれば、読取密度が10本／mm
と５本／mmの場合と同様な非線形処理により、10本／mm
と５本／mmの場合と同様の効果を得ることが可能とな
る。
【００３８】以上、２種類のフィルタを原画像信号の読
取密度に応じて切り替える方法および装置について説明
したが、画像処理装置に入力される画像信号の読取密度
の種類の数に応じて何種類ものフィルタを用意してお
き、選択的に使用することができることはいうまでもな
い。また、本実施の形態では１次元フィルタを使用して
いるが、２次元フィルタでも同様の効果が得られること
は明らかである。
取密度に応じて切り替える方法および装置について説明
したが、画像処理装置に入力される画像信号の読取密度
の種類の数に応じて何種類ものフィルタを用意してお
き、選択的に使用することができることはいうまでもな
い。また、本実施の形態では１次元フィルタを使用して
いるが、２次元フィルタでも同様の効果が得られること
は明らかである。
【００３９】最後に、上記のようにして作成されたボケ
画像信号を用いて行う非線形処理の一例を以下に示す。
以下に示す処理は、特願平8-182155号などに示される周
波数強調処理であり、図15はこの周波数強調処理装置の
一例を示す図である。図15に示すように、原画像信号Ｓ
orgと、フィルタリング手段２および補間拡大手段３に
より作成されたボケ画像信号Ｓuskについて、減算器21
によりそれらの信号の差分がとられ、原画像信号の、限
られた周波数帯域の成分である帯域制限画像信号（Ｓor
g-Ｓus1、Ｓus1-Ｓus2など）が作成される。
画像信号を用いて行う非線形処理の一例を以下に示す。
以下に示す処理は、特願平8-182155号などに示される周
波数強調処理であり、図15はこの周波数強調処理装置の
一例を示す図である。図15に示すように、原画像信号Ｓ
orgと、フィルタリング手段２および補間拡大手段３に
より作成されたボケ画像信号Ｓuskについて、減算器21
によりそれらの信号の差分がとられ、原画像信号の、限
られた周波数帯域の成分である帯域制限画像信号（Ｓor
g-Ｓus1、Ｓus1-Ｓus2など）が作成される。
【００４０】そして、この各帯域制限画像信号が変換器
22によりそれぞれ異なる関数f1〜fNにより変換される。
この関数としては例えば図16に示されるような関数が使
用される。図中の「高」「低」はそれぞれ、高周波帯域
の帯域制限画像信号の変換に使用される関数、低周波帯
域の帯域制限画像信号の変換に使用される関数を示して
おり、高周波帯域を処理する関数ほど小さい値で傾き０
に収束するような関数であることを示している。
22によりそれぞれ異なる関数f1〜fNにより変換される。
この関数としては例えば図16に示されるような関数が使
用される。図中の「高」「低」はそれぞれ、高周波帯域
の帯域制限画像信号の変換に使用される関数、低周波帯
域の帯域制限画像信号の変換に使用される関数を示して
おり、高周波帯域を処理する関数ほど小さい値で傾き０
に収束するような関数であることを示している。
【００４１】次に、演算器23において所定の演算式に基
づいて原画像信号Ｓorgに変換された信号を加算するこ
とによって、それぞれの周波数成分を目的に応じた度合
いで強調する。これらの処理は、下記の式（３） Ｓproc＝Ｓorg ＋β（Ｓorg ）×Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)…（３） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） として表すことができる。
づいて原画像信号Ｓorgに変換された信号を加算するこ
とによって、それぞれの周波数成分を目的に応じた度合
いで強調する。これらの処理は、下記の式（３） Ｓproc＝Ｓorg ＋β（Ｓorg ）×Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆN(ＳusN-1−ＳusN)…（３） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：前記各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） として表すことができる。
【００４２】なお、この周波数強調処理は、ボケ画像信
号を用いて、強調のために原画像信号に加算する加算信
号を作成するものであり、周波数強調処理を行うことに
よってアーチファクトが発生しないようにするために、
この加算信号を構成する各周波数帯域の信号がそれぞれ
所望の信号となるよう、周波数帯域ごとに異なる変換処
理を施すことを特徴としている。所望の信号を作成する
ためには、その信号の基礎となるボケ画像信号が所望の
特性を有する信号であることが好ましい。つまり、従来
の、読取密度によってボケ画像信号の周波数特性が異な
ってしまう従来の画像処理方法では、このような周波数
強調処理の効果を十分に得ることができなかった。本発
明の画像処理方法および装置により、原画像信号の読取
密度に拘わらず、上記のような周波数強調処理、あるい
はその他の例えばダイナミックレンジ圧縮処理におい
て、その画像処理の効果を十分に得ることが可能とな
る。
号を用いて、強調のために原画像信号に加算する加算信
号を作成するものであり、周波数強調処理を行うことに
よってアーチファクトが発生しないようにするために、
この加算信号を構成する各周波数帯域の信号がそれぞれ
所望の信号となるよう、周波数帯域ごとに異なる変換処
理を施すことを特徴としている。所望の信号を作成する
ためには、その信号の基礎となるボケ画像信号が所望の
特性を有する信号であることが好ましい。つまり、従来
の、読取密度によってボケ画像信号の周波数特性が異な
ってしまう従来の画像処理方法では、このような周波数
強調処理の効果を十分に得ることができなかった。本発
明の画像処理方法および装置により、原画像信号の読取
密度に拘わらず、上記のような周波数強調処理、あるい
はその他の例えばダイナミックレンジ圧縮処理におい
て、その画像処理の効果を十分に得ることが可能とな
る。
【図１】本発明の画像処理装置の概略を示す図
【図２】ボケ画像信号作成処理の概要を示すブロック図
【図３】低解像度画像信号作成処理の詳細を示す図
【図４】低解像度画像信号の周波数応答特性を示す図
【図５】補間拡大処理に使用されるフィルタの一例を示
す図
す図
【図６】補間拡大処理におけるフィルタリング処理を説
明するための図
明するための図
【図７】補間拡大処理に使用されるフィルタの他の例を
示す図
示す図
【図８】補間拡大処理におけるフィルタリング処理を説
明するための図
明するための図
【図９】低解像度画像信号作成のためのフィルタリング
処理に使用されるフィルタを示す図
処理に使用されるフィルタを示す図
【図１０】読取密度が10本／mmの原画像信号に基づい
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
【図１１】読取密度が５本／mmの原画像信号に基づい
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
【図１２】読取密度が6.7本／mmの原画像信号に基づい
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
て、図９のフィルタを使用して作成されたボケ画像信号
の周波数応答特性を示す図
【図１３】低解像度画像信号作成のためのフィルタリン
グ処理に使用される他のフィルタを示す図
グ処理に使用される他のフィルタを示す図
【図１４】読取密度が6.7本／mmの原画像信号に基づい
て、図１３のフィルタを使用して作成されたボケ画像信
号の周波数応答特性を示す図
て、図１３のフィルタを使用して作成されたボケ画像信
号の周波数応答特性を示す図
【図１５】周波数強調処理装置の概要を示す図
【図１６】変換関数の一例を示す図
１ 画像処理装置 ２ フィルタリング手段 ３ 補間拡大手段 ４ 非線形処理手段 ５ フィルタ係数選択手段 ６ 読取密度情報入力手段 ７ フィルタ係数リスト ８ 補間係数リスト 21 減算器 22 変換器 23 演算器
Claims (8)
- 【請求項１】 原画像を所定の読取密度で読み取って複
数の画素により構成される原画像信号を得、該原画像信
号に対して所定のフィルタにより所定のフィルタリング
処理を施して前記原画像信号よりも画素数が少ない低解
像度画像信号を作成し、該低解像度画像信号に対して前
記フィルタリング処理をさらに施して、前記低解像度画
像信号よりもさらに画素数が少ない低解像度画像信号を
作成し、前記低解像度画像信号に対するフィルタリング
処理を繰り返し行うことによって、互いに解像度が異な
る複数の低解像度画像信号を作成し、該各低解像度画像
信号を、該低解像度画像信号の画素数が前記原画像信号
の画素数と同じになるようにそれぞれ補間拡大すること
により、前記原画像信号の複数の異なる非鮮鋭マスク画
像信号を作成し、該非鮮鋭マスク画像信号および前記原
画像信号に基づいて所定の非線形処理を行う画像処理方
法において、 前記所定のフィルタとして、前記読取密度に基づいて定
められたフィルタ係数を有するフィルタを使用すること
を特徴とする画像処理方法。 - 【請求項２】 前記非線形処理が、周波数強調処理であ
ることを特徴とする請求項１記載の画像処理方法。 - 【請求項３】 前記非線形処理が、ダイナミックレンジ
圧縮処理であることを特徴とする請求項１記載の画像処
理方法。 - 【請求項４】 原画像を所定の読取密度で読み取って複
数の画素により構成される原画像信号を得、該原画像信
号に対して所定のフィルタにより所定のフィルタリング
処理を施して前記原画像信号よりも画素数が少ない低解
像度画像信号を作成し、該低解像度画像信号に対して前
記フィルタリング処理をさらに施して、前記低解像度画
像信号よりもさらに画素数が少ない低解像度画像信号を
作成し、前記低解像度画像信号に対するフィルタリング
処理を繰り返し行うことによって、互いに解像度が異な
る複数の低解像度画像信号を作成し、該各低解像度画像
信号を、該低解像度画像信号の画素数が前記原画像信号
の画素数と同じになるようにそれぞれ補間拡大すること
により、前記原画像信号の複数の異なる非鮮鋭マスク画
像信号を作成する非鮮鋭マスク画像信号の作成方法にお
いて、 前記所定のフィルタとして、前記読取密度に基づいて定
められたフィルタ係数を有するフィルタを使用すること
を特徴とする非鮮鋭マスク画像信号の作成方法。 - 【請求項５】 原画像を所定の読取密度で読み取って得
られる複数の画素により構成される原画像信号に対して
所定のフィルタにより所定のフィルタリング処理を施し
て前記原画像信号よりも画素数が少ない低解像度画像信
号を作成し、該低解像度画像信号に対して前記フィルタ
リング処理をさらに施して、前記低解像度画像信号より
もさらに画素数が少ない低解像度画像信号を作成し、前
記低解像度画像信号に対するフィルタリング処理を繰り
返し行うことによって、互いに解像度が異なる複数の低
解像度画像信号を作成するフィルタリング手段と、前記
各低解像度画像信号を、該低解像度画像信号の画素数が
前記原画像信号の画素数と同じになるようにそれぞれ補
間拡大することにより、前記原画像信号の複数の異なる
非鮮鋭マスク画像信号を作成する補間拡大手段と、前記
非鮮鋭マスク画像信号および前記原画像信号に基づいて
所定の非線形処理を行う非線形処理手段とからなる画像
処理装置において、 前記読取密度に関する情報を入力するための読取密度情
報入力手段と、 前記情報として入力された読取密度に基づいて、前記所
定のフィルタのフィルタ係数を予め保持された複数のフ
ィルタ係数の中から選択するフィルタ係数選択手段を有
することを特徴とする画像処理装置。 - 【請求項６】 前記非線形処理手段が周波数強調処理を
行う手段であることを特徴とする請求項５記載の画像処
理装置。 - 【請求項７】 前記非線形処理手段がダイナミックレン
ジ圧縮処理を行う手段であることを特徴とする請求項５
記載の画像処理装置。 - 【請求項８】 原画像を所定の読取密度で読み取って得
られる複数の画素により構成される原画像信号に対して
所定のフィルタにより所定のフィルタリング処理を施し
て前記原画像信号よりも画素数が少ない低解像度画像信
号を作成し、該低解像度画像信号に対して前記フィルタ
リング処理をさらに施して、前記低解像度画像信号より
もさらに画素数が少ない低解像度画像信号を作成し、前
記低解像度画像信号に対するフィルタリング処理を繰り
返し行うことによって、互いに解像度が異なる複数の低
解像度画像信号を作成するフィルタリング手段と、前記
各低解像度画像信号を、該低解像度画像信号の画素数が
前記原画像信号の画素数と同じになるようにそれぞれ補
間拡大することにより、前記原画像信号の複数の異なる
非鮮鋭マスク画像信号を作成する補間拡大手段とからな
る非鮮鋭マスク画像信号作成装置において、 前記読取密度に関する情報を入力するための読取密度情
報入力手段と、 前記情報として入力された読取密度に基づいて、前記所
定のフィルタのフィルタ係数を予め保持された複数のフ
ィルタ係数の中から選択するフィルタ係数選択手段を有
することを特徴とする非鮮鋭マスク画像信号作成装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21733396A JP3738791B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
| US08/914,144 US5949915A (en) | 1996-08-19 | 1997-08-19 | Image processing method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21733396A JP3738791B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1063837A true JPH1063837A (ja) | 1998-03-06 |
| JP3738791B2 JP3738791B2 (ja) | 2006-01-25 |
Family
ID=16702542
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP21733396A Expired - Fee Related JP3738791B2 (ja) | 1996-08-19 | 1996-08-19 | 画像処理方法および装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5949915A (ja) |
| JP (1) | JP3738791B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000268171A (ja) * | 1999-01-14 | 2000-09-29 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| JP2007110450A (ja) * | 2005-10-13 | 2007-04-26 | Canon Inc | 画像処理装置 |
| JP2009070407A (ja) * | 1999-01-26 | 2009-04-02 | Fujifilm Corp | 画像処理方法および装置並びに記録媒体 |
| JP2011243233A (ja) * | 2011-09-09 | 2011-12-01 | Olympus Corp | 画像処理装置 |
| US9311695B2 (en) | 2010-03-26 | 2016-04-12 | Shimadzu Corporation | Image processing method and radiographic apparatus using the same |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AUPP779998A0 (en) * | 1998-12-18 | 1999-01-21 | Canon Kabushiki Kaisha | Continuous kernel image interpolation |
| DE60022237T2 (de) | 1999-01-14 | 2006-07-06 | Fuji Photo Film Co., Ltd., Minami-Ashigara | Bidverarbeitungsverfahren und -System und Aufzeichnungsmedium zur Durchfürung des Verfahrens |
| US7177481B2 (en) * | 2000-12-19 | 2007-02-13 | Konica Corporation | Multiresolution unsharp image processing apparatus |
| US7668390B2 (en) * | 2002-11-06 | 2010-02-23 | Digivision, Inc. | Systems and methods for image enhancement in multiple dimensions |
| JP2005040506A (ja) * | 2003-07-25 | 2005-02-17 | Morita Mfg Co Ltd | Ｘ線画像処理方法及び装置 |
| EP1583030A1 (en) * | 2004-03-31 | 2005-10-05 | Fujitsu Limited | Image magnification device and image magnification method |
| JP4999392B2 (ja) * | 2006-07-28 | 2012-08-15 | キヤノン株式会社 | 画像処理装置及びその制御方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
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 |
| FR2445536B1 (fr) * | 1978-12-26 | 1989-12-22 | Fuji Photo Film Co Ltd | Procede et appareil pour traiter une image radiographique |
| US5051902A (en) * | 1987-10-20 | 1991-09-24 | Fuji Photo Film Co., Ltd. | Method and apparatus for radiation image processing and x-ray image processing, including spatial frequency filtering to improve graininess |
| JP2663189B2 (ja) * | 1990-01-29 | 1997-10-15 | 富士写真フイルム株式会社 | 画像のダイナミックレンジ圧縮処理方法 |
-
1996
- 1996-08-19 JP JP21733396A patent/JP3738791B2/ja not_active Expired - Fee Related
-
1997
- 1997-08-19 US US08/914,144 patent/US5949915A/en not_active Expired - Lifetime
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000268171A (ja) * | 1999-01-14 | 2000-09-29 | Fuji Photo Film Co Ltd | 画像処理方法および装置並びに記録媒体 |
| JP2009070407A (ja) * | 1999-01-26 | 2009-04-02 | Fujifilm Corp | 画像処理方法および装置並びに記録媒体 |
| JP4598115B2 (ja) * | 1999-01-26 | 2010-12-15 | 富士フイルム株式会社 | 画像処理方法および装置並びに記録媒体 |
| JP2007110450A (ja) * | 2005-10-13 | 2007-04-26 | Canon Inc | 画像処理装置 |
| JP4649309B2 (ja) * | 2005-10-13 | 2011-03-09 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
| US9311695B2 (en) | 2010-03-26 | 2016-04-12 | Shimadzu Corporation | Image processing method and radiographic apparatus using the same |
| JP2011243233A (ja) * | 2011-09-09 | 2011-12-01 | Olympus Corp | 画像処理装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3738791B2 (ja) | 2006-01-25 |
| US5949915A (en) | 1999-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3882651B2 (ja) | 画像処理装置およびプログラム | |
| JP3816151B2 (ja) | 画像処理方法および装置 | |
| JP3700804B2 (ja) | 画像処理方法および装置 | |
| US20070086671A1 (en) | Image processing apparatus | |
| JPH10509824A (ja) | 画像データの適応的補間のためのシステム及び方法 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| JP3999432B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP6174903B2 (ja) | 空間・階調超解像装置及びプログラム | |
| JP3696339B2 (ja) | 画像処理方法および装置 | |
| JP3783116B2 (ja) | 放射線画像強調処理方法および装置 | |
| JPH0998293A (ja) | 画像処理方法および装置 | |
| JP2000209432A (ja) | 画像処理方法 | |
| JP4598115B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2002183727A (ja) | 画像処理装置 | |
| JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2007215188A (ja) | 画像変換方法、変換画像生成方法、および画像補正装置 | |
| JP4124380B2 (ja) | 画像処理再生装置 | |
| JP3669452B2 (ja) | 画像処理方法および装置 | |
| JP4267159B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPH1185971A (ja) | 階調変換方法および装置 | |
| JP2001056856A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP4213337B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2000268171A (ja) | 画像処理方法および装置並びに記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050701 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050705 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050905 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20051025 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20051025 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081111 Year of fee payment: 3 |
|
| 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: 20081111 Year of fee payment: 3 |
|
| 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: 20091111 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091111 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101111 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111111 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121111 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121111 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131111 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |