JP2002074356A - 画像処理方法および装置並びに記録媒体 - Google Patents
画像処理方法および装置並びに記録媒体Info
- Publication number
- JP2002074356A JP2002074356A JP2000257395A JP2000257395A JP2002074356A JP 2002074356 A JP2002074356 A JP 2002074356A JP 2000257395 A JP2000257395 A JP 2000257395A JP 2000257395 A JP2000257395 A JP 2000257395A JP 2002074356 A JP2002074356 A JP 2002074356A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- image signal
- image
- original image
- frequency
- 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.)
- Withdrawn
Links
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
(57)【要約】
【課題】 周波数強調処理を行う際に、高周波のノイズ
を除去できるとともに、エッジ等のコントラストの大き
な情報の劣化を防止する。 【解決手段】 原画像信号Ｓorgからボケ画像信号Ｓus
を求め、原画像信号Ｓorgからボケ画像信号Ｓusを減算
して、信号Ｓorg−Ｓusを得る。信号Ｓorg−Ｓusを、値
が大きいほど変換後の値が小さくなる非線形の変換関数
ｆにより変換して変換信号ｆ（Ｓorg−Ｓus）を得、こ
れに強調係数β（Ｓorg）を乗算して、原画像信号Ｓorg
から減算する。すなわち、 Ｓproc＝Ｓorg−β（Ｓorg）×ｆ（Ｓorg−Ｓus） の演算を行う。
を除去できるとともに、エッジ等のコントラストの大き
な情報の劣化を防止する。 【解決手段】 原画像信号Ｓorgからボケ画像信号Ｓus
を求め、原画像信号Ｓorgからボケ画像信号Ｓusを減算
して、信号Ｓorg−Ｓusを得る。信号Ｓorg−Ｓusを、値
が大きいほど変換後の値が小さくなる非線形の変換関数
ｆにより変換して変換信号ｆ（Ｓorg−Ｓus）を得、こ
れに強調係数β（Ｓorg）を乗算して、原画像信号Ｓorg
から減算する。すなわち、 Ｓproc＝Ｓorg−β（Ｓorg）×ｆ（Ｓorg−Ｓus） の演算を行う。
Description
【０００１】
【発明の属する技術分野】本発明は、画像信号に対して
所定の周波数成分を強調する周波数強調処理を施す画像
処理方法および装置並びに画像処理方法をコンピュータ
に実行させるためのプログラムを記録したコンピュータ
読取り可能な記録媒体に関するものである。
所定の周波数成分を強調する周波数強調処理を施す画像
処理方法および装置並びに画像処理方法をコンピュータ
に実行させるためのプログラムを記録したコンピュータ
読取り可能な記録媒体に関するものである。
【０００２】
【従来の技術】本出願人により、非鮮鋭マスク画像信号
（以下、ボケ画像信号という）を用いて周波数強調処理
あるいはダイナミックレンジ圧縮処理等を行って放射線
画像の診断性能を向上させる数々の画像処理方法および
装置が提案されている（特開昭55-163472号、同55-8795
3号、特開平3-222577号、同10-75395号、同10-75364
号、同10-171983号等）。例えば周波数強調処理は、原
画像信号Ｓorgからボケ画像信号Ｓusを引いたものに強
調係数βを乗じたものを、原画像信号Ｓorgに加算する
ことにより、原画像信号の所定の空間周波数成分を強調
するものである。これを式で表すと下記の式（１）のよ
うになる。
（以下、ボケ画像信号という）を用いて周波数強調処理
あるいはダイナミックレンジ圧縮処理等を行って放射線
画像の診断性能を向上させる数々の画像処理方法および
装置が提案されている（特開昭55-163472号、同55-8795
3号、特開平3-222577号、同10-75395号、同10-75364
号、同10-171983号等）。例えば周波数強調処理は、原
画像信号Ｓorgからボケ画像信号Ｓusを引いたものに強
調係数βを乗じたものを、原画像信号Ｓorgに加算する
ことにより、原画像信号の所定の空間周波数成分を強調
するものである。これを式で表すと下記の式（１）のよ
うになる。
【０００３】 Ｓproc＝Ｓorg＋β×（Ｓorg−Ｓus） （１） （Ｓproc：周波数強調処理された信号、Ｓorg：原画像
信号、Ｓus：ボケ画像信号、β：強調係数） また、式（１）に示す強調係数βの値を負とする周波数
強調処理方法も本出願人により提案されている（特開平
1-106277号、同1-239681号）。このように、強調係数β
の値を負とすることにより、原画像から高周波のノイズ
を除去することができる。
信号、Ｓus：ボケ画像信号、β：強調係数） また、式（１）に示す強調係数βの値を負とする周波数
強調処理方法も本出願人により提案されている（特開平
1-106277号、同1-239681号）。このように、強調係数β
の値を負とすることにより、原画像から高周波のノイズ
を除去することができる。
【０００４】また、特開平10-75395号には、原画像信号
に加算する加算信号の周波数応答特性を調整することに
より、周波数強調処理された信号に対してアーチファク
トが発生することを防止する方法が提案されている。こ
の方法とは、まず鮮鋭度の異なる、すなわち周波数応答
特性の異なる複数のボケ画像信号を作成し、そのボケ画
像信号および原画像信号の中の２つの信号の差分をとる
ことにより、原画像信号の、ある限られた周波数帯域の
周波数成分を表す複数の帯域制限画像信号を作成し、さ
らにその帯域制限画像信号をそれぞれ異なる変換関数に
よって所望の大きさとなるように変換してから、その複
数の抑制された帯域制限画像信号を積算することにより
上記加算信号を作成するものである。この処理は例えば
下記の式（２）により表すことができる。
に加算する加算信号の周波数応答特性を調整することに
より、周波数強調処理された信号に対してアーチファク
トが発生することを防止する方法が提案されている。こ
の方法とは、まず鮮鋭度の異なる、すなわち周波数応答
特性の異なる複数のボケ画像信号を作成し、そのボケ画
像信号および原画像信号の中の２つの信号の差分をとる
ことにより、原画像信号の、ある限られた周波数帯域の
周波数成分を表す複数の帯域制限画像信号を作成し、さ
らにその帯域制限画像信号をそれぞれ異なる変換関数に
よって所望の大きさとなるように変換してから、その複
数の抑制された帯域制限画像信号を積算することにより
上記加算信号を作成するものである。この処理は例えば
下記の式（２）により表すことができる。
【０００５】 Ｓ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-75364号には、ダイナミックレンジ圧
縮処理を施す場合において、処理が施された信号に対し
てアーチファクトが発生することを防止する方法が提案
されている。この方法とは、上記特開平10-75395号に記
載されたように、複数の帯域制限画像信号を作成し、こ
の帯域制限画像信号に基づいて原画像信号の低周波成分
に関する信号（低周波成分信号）を得、原画像信号にこ
の低周波成分に関する信号を加算することによりダイナ
ミックレンジ圧縮処理を施すようにしたものである。こ
の処理は例えば下記の式（３）により表すことができ
る。
数） さらに、特開平10-75364号には、ダイナミックレンジ圧
縮処理を施す場合において、処理が施された信号に対し
てアーチファクトが発生することを防止する方法が提案
されている。この方法とは、上記特開平10-75395号に記
載されたように、複数の帯域制限画像信号を作成し、こ
の帯域制限画像信号に基づいて原画像信号の低周波成分
に関する信号（低周波成分信号）を得、原画像信号にこ
の低周波成分に関する信号を加算することによりダイナ
ミックレンジ圧縮処理を施すようにしたものである。こ
の処理は例えば下記の式（３）により表すことができ
る。
【０００６】 Ｓproc＝Ｓorg ＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）） （３ ） Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）＝｛ｆd1(Ｓorg−
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） さらにまた、特開平10-171983号には、周波数強調処理
とダイナミックレンジ圧縮処理とを同時に施す場合にお
いて、処理が施された信号に対してアーチファクトが発
生することを防止する方法が提案されている。この方法
とは、上記特開平10-75395号に記載されたように、複数
の帯域制限画像信号を作成し、この帯域制限画像信号に
基づいて原画像信号の高周波成分に関する信号（高周波
成分信号）および低周波成分信号を得、原画像信号にこ
れらの高周波成分に関する信号および低周波成分に関す
る信号を加算することにより周波数強調処理を施すとと
もにダイナミックレンジ圧縮処理を施すようにしたもの
である。この処理は例えば下記の式（４）により表すこ
とができる。
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） さらにまた、特開平10-171983号には、周波数強調処理
とダイナミックレンジ圧縮処理とを同時に施す場合にお
いて、処理が施された信号に対してアーチファクトが発
生することを防止する方法が提案されている。この方法
とは、上記特開平10-75395号に記載されたように、複数
の帯域制限画像信号を作成し、この帯域制限画像信号に
基づいて原画像信号の高周波成分に関する信号（高周波
成分信号）および低周波成分信号を得、原画像信号にこ
れらの高周波成分に関する信号および低周波成分に関す
る信号を加算することにより周波数強調処理を施すとと
もにダイナミックレンジ圧縮処理を施すようにしたもの
である。この処理は例えば下記の式（４）により表すこ
とができる。
【０００７】 Ｓproc＝Ｓorg ＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）） （４ ） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）＝｛ｆu1(Ｓorg−
Ｓus1)＋ｆu2(Ｓus1 −Ｓus2)＋…＋ｆuk(Ｓusk-1−Ｓu
sk )＋…＋ｆun(Ｓusn-1−Ｓusn)｝ Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）＝｛ｆd1(Ｓorg−
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆuk(k=1〜n）：高周波成分信号を得るために使用する
変換関数 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） これらの周波数強調処理やダイナミックレンジ圧縮処理
（以下変換処理とする）においては、帯域制限画像信号
を変換する変換関数等の定義を変更することによって原
画像信号に加算する加算信号の周波数応答特性を調整す
ることができる。このため、各変換関数の定義次第で、
アーチファクトの発生防止等所望の周波数応答特性を有
する処理済み画像信号を得ることができる。
Ｓus1)＋ｆu2(Ｓus1 −Ｓus2)＋…＋ｆuk(Ｓusk-1−Ｓu
sk )＋…＋ｆun(Ｓusn-1−Ｓusn)｝ Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）＝｛ｆd1(Ｓorg−
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆuk(k=1〜n）：高周波成分信号を得るために使用する
変換関数 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） これらの周波数強調処理やダイナミックレンジ圧縮処理
（以下変換処理とする）においては、帯域制限画像信号
を変換する変換関数等の定義を変更することによって原
画像信号に加算する加算信号の周波数応答特性を調整す
ることができる。このため、各変換関数の定義次第で、
アーチファクトの発生防止等所望の周波数応答特性を有
する処理済み画像信号を得ることができる。
【０００８】
【発明が解決しようとする課題】上述した特開平1-1062
77号等に記載された、周波数強調処理の際に強調係数を
負とする方法は、原画像から高周波のノイズを除去する
ことができるものの、画像中に含まれるコントラストが
大きなエッジなどの高周波情報も低減してしまうため、
得られる画像がぼけたものとなってしまう。
77号等に記載された、周波数強調処理の際に強調係数を
負とする方法は、原画像から高周波のノイズを除去する
ことができるものの、画像中に含まれるコントラストが
大きなエッジなどの高周波情報も低減してしまうため、
得られる画像がぼけたものとなってしまう。
【０００９】本発明は上記事情に鑑みなされたものであ
り、周波数強調処理を行う際に、高周波のノイズを除去
できるとともに、エッジ等のコントラストの大きな情報
の劣化を防止できる画像処理方法および装置並びに画像
処理方法をコンピュータに実行させるためのプログラム
を記録したコンピュータ読取り可能な記録媒体を提供す
ることを目的とするものである。
り、周波数強調処理を行う際に、高周波のノイズを除去
できるとともに、エッジ等のコントラストの大きな情報
の劣化を防止できる画像処理方法および装置並びに画像
処理方法をコンピュータに実行させるためのプログラム
を記録したコンピュータ読取り可能な記録媒体を提供す
ることを目的とするものである。
【００１０】
【課題を解決するための手段】本発明による第１の画像
処理方法は、原画像を表す原画像信号に対して、周波数
強調処理を施す画像処理方法において、前記原画像信号
から前記原画像の高周波成分を表す高周波信号を作成
し、該高周波信号を、信号値が大きいほどその値を抑制
する変換関数により変換して変換信号を得、該変換信号
を前記原画像信号から減算して処理済み画像信号を得る
ことを特徴とするものである。
処理方法は、原画像を表す原画像信号に対して、周波数
強調処理を施す画像処理方法において、前記原画像信号
から前記原画像の高周波成分を表す高周波信号を作成
し、該高周波信号を、信号値が大きいほどその値を抑制
する変換関数により変換して変換信号を得、該変換信号
を前記原画像信号から減算して処理済み画像信号を得る
ことを特徴とするものである。
【００１１】「変換関数」は、信号値が大きいほどその
値を抑制するような関数であり、具体的には特開平10-7
8395号に記載されたような非線形関数を用いることが好
ましい。
値を抑制するような関数であり、具体的には特開平10-7
8395号に記載されたような非線形関数を用いることが好
ましい。
【００１２】「変換信号を原画像信号から減算する」と
は、原画像信号から変換信号をそのまま減算するもので
あってもよく、変換信号に強調係数を乗算し、強調係数
が乗算された変換信号を原画像信号から減算するもので
あってもよい。
は、原画像信号から変換信号をそのまま減算するもので
あってもよく、変換信号に強調係数を乗算し、強調係数
が乗算された変換信号を原画像信号から減算するもので
あってもよい。
【００１３】本発明による第２の画像処理方法は、原画
像を表す原画像信号に対して、周波数強調処理を施す画
像処理方法において、前記原画像信号から複数の帯域制
限画像信号を作成し、該各帯域制限画像信号を、信号値
が大きいほどその値を抑制する変換関数により変換して
複数の変換帯域制限画像信号を得、該複数の変換帯域制
限画像信号を加算して加算変換信号を得、該加算変換信
号を前記原画像信号から減算して処理済み画像信号を得
ることを特徴とするものである。
像を表す原画像信号に対して、周波数強調処理を施す画
像処理方法において、前記原画像信号から複数の帯域制
限画像信号を作成し、該各帯域制限画像信号を、信号値
が大きいほどその値を抑制する変換関数により変換して
複数の変換帯域制限画像信号を得、該複数の変換帯域制
限画像信号を加算して加算変換信号を得、該加算変換信
号を前記原画像信号から減算して処理済み画像信号を得
ることを特徴とするものである。
【００１４】「原画像信号を多重解像度変換することに
より帯域制限画像信号を作成」するに際しては、ラプラ
シアンピラミッドの手法によるラプラシアンピラミッド
分解により、あるいはウェーブレット変換により原画像
信号を複数の周波数帯域毎の周波数応答特性を表す信号
に変換する方法などを用いることができる。
より帯域制限画像信号を作成」するに際しては、ラプラ
シアンピラミッドの手法によるラプラシアンピラミッド
分解により、あるいはウェーブレット変換により原画像
信号を複数の周波数帯域毎の周波数応答特性を表す信号
に変換する方法などを用いることができる。
【００１５】ラプラシアンピラミッドの手法は例えば特
開平5-244508号、特開平6-96200 、特開平6-301766号等
に記載されており、このラプラシアンピラミッドは、原
画像に対してガウス関数で近似されたようなマスクによ
りマスク処理を施した後、画像をサブサンプリングして
画素数を間引いて半分にすることにより、原画像の１／
４のサイズのボケ画像を得、このボケ画像のサンプリン
グされた画素に値が０の画素を補間して元の大きさの画
像に戻し、この画像に対してさらに上述したマスクによ
りマスク処理を施し、このマスク処理が施されたボケ画
像を原画像から減算して、原画像信号のある限られた周
波数帯域の周波数成分を表す、すなわち原画像の複数の
周波数帯域毎の周波数応答特性を表す帯域制限画像信号
を得るものである。この処理を得られたボケ画像に対し
て繰り返すことにより原画像の１／２2Nのサイズを有す
る帯域制限画像信号をＮ個作成するものである。なお、
最低周波数帯域のボケ画像は原画像の低周波成分を表す
ものとなる。
開平5-244508号、特開平6-96200 、特開平6-301766号等
に記載されており、このラプラシアンピラミッドは、原
画像に対してガウス関数で近似されたようなマスクによ
りマスク処理を施した後、画像をサブサンプリングして
画素数を間引いて半分にすることにより、原画像の１／
４のサイズのボケ画像を得、このボケ画像のサンプリン
グされた画素に値が０の画素を補間して元の大きさの画
像に戻し、この画像に対してさらに上述したマスクによ
りマスク処理を施し、このマスク処理が施されたボケ画
像を原画像から減算して、原画像信号のある限られた周
波数帯域の周波数成分を表す、すなわち原画像の複数の
周波数帯域毎の周波数応答特性を表す帯域制限画像信号
を得るものである。この処理を得られたボケ画像に対し
て繰り返すことにより原画像の１／２2Nのサイズを有す
る帯域制限画像信号をＮ個作成するものである。なお、
最低周波数帯域のボケ画像は原画像の低周波成分を表す
ものとなる。
【００１６】なお、ラプラシアンピラミッド分解やウェ
ーブレット変換などにより原画像信号を多重解像度変換
した場合、最低周波数帯域の信号は原画像を縮小した低
周波の情報を表すものであり、原画像信号の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号ではないた
め、本発明においては処理には使用しないか、または値
を０として使用する。
ーブレット変換などにより原画像信号を多重解像度変換
した場合、最低周波数帯域の信号は原画像を縮小した低
周波の情報を表すものであり、原画像信号の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号ではないた
め、本発明においては処理には使用しないか、または値
を０として使用する。
【００１７】「加算変換信号を原画像信号から減算す
る」とは、原画像信号から加算変換信号をそのまま減算
するものであってもよく、加算変換信号に強調係数を乗
算し、強調係数が乗算された加算変換信号を原画像信号
から減算するものであってもよい。
る」とは、原画像信号から加算変換信号をそのまま減算
するものであってもよく、加算変換信号に強調係数を乗
算し、強調係数が乗算された加算変換信号を原画像信号
から減算するものであってもよい。
【００１８】なお、本発明による第１および第２の画像
処理方法においては、前記処理済み画像信号に対して、
さらにダイナミックレンジ圧縮処理を施してダイナミッ
クレンジ圧縮処理済み画像信号を得るようにしてもよ
く、さらに、前記ダイナミックレンジ圧縮処理済み画像
信号に対して、階調変換処理を施して階調変換処理済み
画像信号を得るようにしてもよい。
処理方法においては、前記処理済み画像信号に対して、
さらにダイナミックレンジ圧縮処理を施してダイナミッ
クレンジ圧縮処理済み画像信号を得るようにしてもよ
く、さらに、前記ダイナミックレンジ圧縮処理済み画像
信号に対して、階調変換処理を施して階調変換処理済み
画像信号を得るようにしてもよい。
【００１９】本発明による第１の画像処理装置は、原画
像を表す原画像信号に対して、周波数強調処理を施す画
像処理装置において、前記原画像信号から前記原画像の
高周波成分を表す高周波信号を作成する高周波信号作成
手段と、該高周波信号を、信号値が大きいほどその値を
抑制する変換関数により変換して変換信号を得る変換手
段と、該変換信号を前記原画像信号から減算して処理済
み画像信号を得る減算手段とを備えたことを特徴とする
ものである。
像を表す原画像信号に対して、周波数強調処理を施す画
像処理装置において、前記原画像信号から前記原画像の
高周波成分を表す高周波信号を作成する高周波信号作成
手段と、該高周波信号を、信号値が大きいほどその値を
抑制する変換関数により変換して変換信号を得る変換手
段と、該変換信号を前記原画像信号から減算して処理済
み画像信号を得る減算手段とを備えたことを特徴とする
ものである。
【００２０】本発明による第２の画像処理装置は、原画
像を表す原画像信号に対して、周波数強調処理を施す画
像処理装置において、前記原画像信号から複数の帯域制
限画像信号を作成する帯域制限画像信号作成手段と、該
各帯域制限画像信号を、信号値が大きいほどその値を抑
制する変換関数により変換して複数の変換帯域制限画像
信号を得る変換手段と、該複数の変換帯域制限画像信号
を加算して加算変換信号を得る加算手段と、該加算変換
信号を前記原画像信号から減算して処理済み画像信号を
得る減算手段とを備えたことを特徴とするものである。
像を表す原画像信号に対して、周波数強調処理を施す画
像処理装置において、前記原画像信号から複数の帯域制
限画像信号を作成する帯域制限画像信号作成手段と、該
各帯域制限画像信号を、信号値が大きいほどその値を抑
制する変換関数により変換して複数の変換帯域制限画像
信号を得る変換手段と、該複数の変換帯域制限画像信号
を加算して加算変換信号を得る加算手段と、該加算変換
信号を前記原画像信号から減算して処理済み画像信号を
得る減算手段とを備えたことを特徴とするものである。
【００２１】なお、本発明による第１および第２の画像
処理装置においては、前記処理済み画像信号に対して、
ダイナミックレンジ圧縮処理を施してダイナミックレン
ジ圧縮処理済み画像信号を得るダイナミックレンジ圧縮
処理手段をさらに備えるようにしてもよく、前記ダイナ
ミックレンジ圧縮処理済み画像信号に対して、階調変換
処理を施して階調変換処理済み画像信号を得る階調処理
手段をさらに備えるようにしてもよい。
処理装置においては、前記処理済み画像信号に対して、
ダイナミックレンジ圧縮処理を施してダイナミックレン
ジ圧縮処理済み画像信号を得るダイナミックレンジ圧縮
処理手段をさらに備えるようにしてもよく、前記ダイナ
ミックレンジ圧縮処理済み画像信号に対して、階調変換
処理を施して階調変換処理済み画像信号を得る階調処理
手段をさらに備えるようにしてもよい。
【００２２】なお、本発明による第１および第２の画像
処理方法をコンピュータに実行させるためのプログラム
として、コンピュータ読取り可能な記録媒体に記録して
提供してもよい。
処理方法をコンピュータに実行させるためのプログラム
として、コンピュータ読取り可能な記録媒体に記録して
提供してもよい。
【００２３】
【発明の効果】本発明による第１の画像処理方法および
装置によれば、原画像の高周波成分を表す高周波信号を
作成し、この高周波信号を信号値が大きいほどその値を
抑制する変換関数により変換して変換信号を得、この変
換信号を原画像信号から減算して処理済み画像信号が得
られる。
装置によれば、原画像の高周波成分を表す高周波信号を
作成し、この高周波信号を信号値が大きいほどその値を
抑制する変換関数により変換して変換信号を得、この変
換信号を原画像信号から減算して処理済み画像信号が得
られる。
【００２４】また、本発明による第２の画像処理方法お
よび装置によれば、原画像から帯域制限画像信号を作成
し、この帯域制限画像信号を信号値が大きいほどその値
を抑制する変換関数により変換して変換帯域制限画像信
号を得る。そして、変換帯域制限画像信号を加算して加
算変換信号を得、この加算変換信号を原画像信号から減
算して処理済み画像信号が得られる。
よび装置によれば、原画像から帯域制限画像信号を作成
し、この帯域制限画像信号を信号値が大きいほどその値
を抑制する変換関数により変換して変換帯域制限画像信
号を得る。そして、変換帯域制限画像信号を加算して加
算変換信号を得、この加算変換信号を原画像信号から減
算して処理済み画像信号が得られる。
【００２５】ここで、原画像中に含まれるコントラスト
が大きな部分は、高周波信号または帯域制限画像信号に
おいて比較的大きな信号値を有するものとなる。本発明
により、高周波信号または帯域制限画像信号を変換関数
により変換すると、信号値が大きいほどその値が抑制さ
れるため、原画像中に含まれるコントラストが大きい部
分についてもその信号値は抑制されることとなる。した
がって、変換信号または加算変換信号を原画像信号から
減算することにより、原画像から高周波のノイズを除去
でき、かつコントラストが大きい部分についても情報が
低減されるのを防止することができる。
が大きな部分は、高周波信号または帯域制限画像信号に
おいて比較的大きな信号値を有するものとなる。本発明
により、高周波信号または帯域制限画像信号を変換関数
により変換すると、信号値が大きいほどその値が抑制さ
れるため、原画像中に含まれるコントラストが大きい部
分についてもその信号値は抑制されることとなる。した
がって、変換信号または加算変換信号を原画像信号から
減算することにより、原画像から高周波のノイズを除去
でき、かつコントラストが大きい部分についても情報が
低減されるのを防止することができる。
【００２６】また、処理済み画像信号に対してさらにダ
イナミックレンジ圧縮処理を施すことにより、さらに低
濃度域および高濃度域においてつぶれのない画像を得る
ことができる。
イナミックレンジ圧縮処理を施すことにより、さらに低
濃度域および高濃度域においてつぶれのない画像を得る
ことができる。
【００２７】さらに、ダイナミックレンジ圧縮処理済み
画像信号に対して階調変換処理を施すことにより、コン
トラストも良好な画像を得ることができる。
画像信号に対して階調変換処理を施すことにより、コン
トラストも良好な画像を得ることができる。
【００２８】
【発明の実施の形態】以下、本発明の実施形態につい
て、図面を参照して詳細に説明する。以下に示す画像処
理装置は、蓄積性蛍光体シートに記録された人体の放射
線画像を読み取って得た画像信号に対して、その画像が
診断に適した画像となるように、ボケ画像信号を使用し
て周波数強調処理を施すものであり、処理された画像信
号は主としてフィルムに記録され、診断に用いられる。
て、図面を参照して詳細に説明する。以下に示す画像処
理装置は、蓄積性蛍光体シートに記録された人体の放射
線画像を読み取って得た画像信号に対して、その画像が
診断に適した画像となるように、ボケ画像信号を使用し
て周波数強調処理を施すものであり、処理された画像信
号は主としてフィルムに記録され、診断に用いられる。
【００２９】図１は本発明の第１の実施形態による画像
処理装置の構成を示す概略ブロック図である。画像処理
装置１は、読取装置等において得られた原画像信号Ｓor
gから原画像のボケ画像を表すボケ画像信号Ｓusを作成
するボケ画像信号作成手段２と、原画像信号Ｓorgから
ボケ画像信号Ｓusを減算して原画像の高周波成分を表す
高周波信号Ｓorg−Ｓusを得る減算器３と、高周波信号
Ｓorg−Ｓusを非線形の変換関数ｆにより変換して変換
信号ｆ（Ｓorg−Ｓus）を得る変換器４と、変換信号ｆ
（Ｓorg−Ｓus）に正の強調係数β（Ｓorg）を乗じて強
調変換信号β（Ｓorg）・ｆ（Ｓorg−Ｓus）を得、さら
に原画像信号Ｓorgから強調変換信号β（Ｓorg）・ｆ
（Ｓorg−Ｓus）を減算して処理済み画像信号Ｓprocを
得る演算器５とを備える。なお、ボケ画像信号作成手段
２および減算器３が本発明による高周波信号作成手段、
変換器４が本発明による変換手段、演算器５が本発明に
よる減算手段と観念される。
処理装置の構成を示す概略ブロック図である。画像処理
装置１は、読取装置等において得られた原画像信号Ｓor
gから原画像のボケ画像を表すボケ画像信号Ｓusを作成
するボケ画像信号作成手段２と、原画像信号Ｓorgから
ボケ画像信号Ｓusを減算して原画像の高周波成分を表す
高周波信号Ｓorg−Ｓusを得る減算器３と、高周波信号
Ｓorg−Ｓusを非線形の変換関数ｆにより変換して変換
信号ｆ（Ｓorg−Ｓus）を得る変換器４と、変換信号ｆ
（Ｓorg−Ｓus）に正の強調係数β（Ｓorg）を乗じて強
調変換信号β（Ｓorg）・ｆ（Ｓorg−Ｓus）を得、さら
に原画像信号Ｓorgから強調変換信号β（Ｓorg）・ｆ
（Ｓorg−Ｓus）を減算して処理済み画像信号Ｓprocを
得る演算器５とを備える。なお、ボケ画像信号作成手段
２および減算器３が本発明による高周波信号作成手段、
変換器４が本発明による変換手段、演算器５が本発明に
よる減算手段と観念される。
【００３０】ボケ画像信号作成手段２では、原画像信号
Ｓorgに基づいて、各画素の周囲のＮ×Ｎ個の走査点の
信号値を平均化してボケ画像信号Ｓusが求められる。
Ｓorgに基づいて、各画素の周囲のＮ×Ｎ個の走査点の
信号値を平均化してボケ画像信号Ｓusが求められる。
【００３１】変換器４において用いられる変換関数ｆの
例を図２に示す。この変換関数ｆは、信号Ｓorg−Ｓus
の絶対値が閾値Ｔｈ１よりも小さい場合は傾きが１であ
り、閾値Ｔｈ１よりも大きい場合は傾きが１よりも小さ
くなるような非線形関数である。
例を図２に示す。この変換関数ｆは、信号Ｓorg−Ｓus
の絶対値が閾値Ｔｈ１よりも小さい場合は傾きが１であ
り、閾値Ｔｈ１よりも大きい場合は傾きが１よりも小さ
くなるような非線形関数である。
【００３２】演算器５には、原画像信号Ｓorgと強調係
数β（Ｓorg）との対応テーブルが記憶されており、原
画像信号Ｓorgの各画素値に対応する強調係数β（Ｓor
g）が変換信号ｆ（Ｓorg−Ｓus）に乗算される。
数β（Ｓorg）との対応テーブルが記憶されており、原
画像信号Ｓorgの各画素値に対応する強調係数β（Ｓor
g）が変換信号ｆ（Ｓorg−Ｓus）に乗算される。
【００３３】なお、第１の実施形態による画像処理装置
１において行われる処理は下記の式（５）に示すものと
なる。
１において行われる処理は下記の式（５）に示すものと
なる。
【００３４】 Ｓproc＝Ｓorg−β（Ｓorg）×ｆ（Ｓorg−Ｓus） （５） （Ｓproc：処理済み画像信号、Ｓorg：原画像信号、Ｓu
s：ボケ画像信号、β（Ｓorg）：強調係数、ｆ：変換関
数） 次いで、第１の実施形態の動作について説明する。図３
は第１の実施形態の動作を示すフローチャート、図４は
第１の実施形態の動作を説明するための図である。ま
ず、原画像信号Ｓorgからボケ画像信号Ｓusが作成され
る（ステップＳ１）。ここで、原画像信号Ｓorgのエッ
ジ部分におけるプロファイルを図４（ａ）に示す。この
ような原画像信号Ｓorgから得られるボケ画像信号Ｓus
のプロファイルは図４（ｂ）に示すものとなる。図４
（ｂ）に示すようにボケ画像信号Ｓusにおいては、原画
像信号Ｓorgにおけるエッジ部分がぼけたものとなる。
s：ボケ画像信号、β（Ｓorg）：強調係数、ｆ：変換関
数） 次いで、第１の実施形態の動作について説明する。図３
は第１の実施形態の動作を示すフローチャート、図４は
第１の実施形態の動作を説明するための図である。ま
ず、原画像信号Ｓorgからボケ画像信号Ｓusが作成され
る（ステップＳ１）。ここで、原画像信号Ｓorgのエッ
ジ部分におけるプロファイルを図４（ａ）に示す。この
ような原画像信号Ｓorgから得られるボケ画像信号Ｓus
のプロファイルは図４（ｂ）に示すものとなる。図４
（ｂ）に示すようにボケ画像信号Ｓusにおいては、原画
像信号Ｓorgにおけるエッジ部分がぼけたものとなる。
【００３５】次いで、減算器３において原画像信号Ｓor
gからボケ画像信号Ｓusが減算されて高周波信号Ｓorg−
Ｓusが得られる（ステップＳ２）。高周波信号Ｓorg−
Ｓusのプロファイルを図４（ｃ）に示す。図４（ｃ）に
示すように、高周波信号Ｓorg−Ｓusは、エッジ部分に
比較的大きな値を有するプロファイルとなる。そして、
変換器４において、高周波信号Ｓorg−Ｓusが図３に示
す変換関数ｆにより変換されて変換信号ｆ（Ｓorg−Ｓu
s）が得られる（ステップＳ３）。ここで、変換関数ｆ
は、高周波信号Ｓorg−Ｓusの絶対値が閾値Ｔｈ１より
も小さい場合は傾きが１であり、閾値Ｔｈ１よりも大き
い場合は傾きが１よりも小さくなるような非線形関数で
ある。したがって、変換信号ｆ（Ｓorg−Ｓus）は図４
（ｄ）に示すように、高周波信号Ｓorg−Ｓusにおいて
閾値Ｔｈ１を越える部分（原画像のエッジに対応）の信
号値が抑制されたプロファイルを有するものとなる。
gからボケ画像信号Ｓusが減算されて高周波信号Ｓorg−
Ｓusが得られる（ステップＳ２）。高周波信号Ｓorg−
Ｓusのプロファイルを図４（ｃ）に示す。図４（ｃ）に
示すように、高周波信号Ｓorg−Ｓusは、エッジ部分に
比較的大きな値を有するプロファイルとなる。そして、
変換器４において、高周波信号Ｓorg−Ｓusが図３に示
す変換関数ｆにより変換されて変換信号ｆ（Ｓorg−Ｓu
s）が得られる（ステップＳ３）。ここで、変換関数ｆ
は、高周波信号Ｓorg−Ｓusの絶対値が閾値Ｔｈ１より
も小さい場合は傾きが１であり、閾値Ｔｈ１よりも大き
い場合は傾きが１よりも小さくなるような非線形関数で
ある。したがって、変換信号ｆ（Ｓorg−Ｓus）は図４
（ｄ）に示すように、高周波信号Ｓorg−Ｓusにおいて
閾値Ｔｈ１を越える部分（原画像のエッジに対応）の信
号値が抑制されたプロファイルを有するものとなる。
【００３６】そして、演算器５において、変換信号ｆ
（Ｓorg−Ｓus）に強調係数β（Ｓorg）が乗算されると
ともに、原画像信号Ｓorgから減算されて処理済み画像
信号Ｓprocが得られる。処理済み画像信号Ｓprocのプロ
ファイルを図４（ｅ）に示す。図４（ｅ）に示すよう
に、処理済み画像信号Ｓprocにおいては、信号値が一定
の平坦部における高周波ノイズが除去されているととも
に、エッジ部分のコントラストも保持されたものとな
る。
（Ｓorg−Ｓus）に強調係数β（Ｓorg）が乗算されると
ともに、原画像信号Ｓorgから減算されて処理済み画像
信号Ｓprocが得られる。処理済み画像信号Ｓprocのプロ
ファイルを図４（ｅ）に示す。図４（ｅ）に示すよう
に、処理済み画像信号Ｓprocにおいては、信号値が一定
の平坦部における高周波ノイズが除去されているととも
に、エッジ部分のコントラストも保持されたものとな
る。
【００３７】このように、第１の実施形態においては、
高周波信号Ｓorg−Ｓusを図３に示す非線形関数ｆによ
り変換した後、強調係数β（Ｓorg）を乗算して原画像
信号Ｓorgから減算することにより、処理済み画像信号
Ｓprocを得るようにしたため、原画像から高周波のノイ
ズを除去でき、かつコントラストが大きい部分について
も情報の低減のない処理済み画像信号Ｓprocを得ること
ができる。
高周波信号Ｓorg−Ｓusを図３に示す非線形関数ｆによ
り変換した後、強調係数β（Ｓorg）を乗算して原画像
信号Ｓorgから減算することにより、処理済み画像信号
Ｓprocを得るようにしたため、原画像から高周波のノイ
ズを除去でき、かつコントラストが大きい部分について
も情報の低減のない処理済み画像信号Ｓprocを得ること
ができる。
【００３８】次いで、本発明の第２の実施形態について
説明する。図５は本発明の第２の実施形態による画像処
理装置の構成を示す概略ブロック図である。図５に示す
ように、本発明の第２の実施形態による画像処理装置１
０は、原画像信号Ｓorgの複数の周波数帯域毎の周波数
応答特性を表す帯域制限画像信号を作成する帯域制限画
像信号作成手段１２と、帯域制限画像信号に基づいて原
画像信号Ｓorgに対して特定の周波数を強調するための
周波数強調処理を行って処理済み画像信号Ｓprocを得る
周波数強調処理手段１３とを有する。
説明する。図５は本発明の第２の実施形態による画像処
理装置の構成を示す概略ブロック図である。図５に示す
ように、本発明の第２の実施形態による画像処理装置１
０は、原画像信号Ｓorgの複数の周波数帯域毎の周波数
応答特性を表す帯域制限画像信号を作成する帯域制限画
像信号作成手段１２と、帯域制限画像信号に基づいて原
画像信号Ｓorgに対して特定の周波数を強調するための
周波数強調処理を行って処理済み画像信号Ｓprocを得る
周波数強調処理手段１３とを有する。
【００３９】まず帯域制限画像信号の作成処理について
詳細に説明する。図６は帯域制限画像信号作成処理の概
要を示すブロック図、図７は帯域制限画像信号作成処理
を模式的に示す図である。なお、本実施形態において
は、例えば特開平5-244508号に記載されたラプラシアン
ピラミッドの手法により帯域制限画像信号を作成するも
のとする。まず図６に示すように、図５の帯域制限画像
信号作成手段１２は、フィルタリング処理手段２０にお
いて原画像信号Ｓorgに対し、原画像信号Ｓorgを得た際
の主走査方向（Ｘ方向）および副走査方向（Ｙ方向）
（図８参照）に対してフィルタリング処理を施して原画
像信号Ｓorgよりも解像度が低い画像信号Ｌ1（以下、低
解像度画像信号という）を作成し、次にこの低解像度画
像信号Ｌ1に対して同様のフィルタリング処理を施して
この低解像度画像信号Ｌ1よりもさらに解像度が低い低
解像度画像信号Ｌ2を作成し、以降順次同様のフィルタ
リング処理を繰り返して各解像度毎の低解像度画像信号
Ｌk（ｋ＝１〜ｎ）を得るものである。そして、補間処
理手段２１において、このフィルタリング処理の各段に
おいて得られる低解像度画像信号Ｌkに対して、それぞ
れ２倍の画素数となるように補間処理を施して、鮮鋭度
が異なる複数のボケ画像信号Ｓus1〜Ｓusｎ（以下Ｓusk
（ｋ＝１〜ｎ）で代表させる）を得る。この後、減算器
２２により互いに対応する画素数を有する低解像度画像
信号Ｌk-1とボケ画像信号Ｓuskおよび原画像信号Ｓorg
とボケ画像信号Ｓus１との差分を求め、これを帯域制限
画像信号Ｂkとする。
詳細に説明する。図６は帯域制限画像信号作成処理の概
要を示すブロック図、図７は帯域制限画像信号作成処理
を模式的に示す図である。なお、本実施形態において
は、例えば特開平5-244508号に記載されたラプラシアン
ピラミッドの手法により帯域制限画像信号を作成するも
のとする。まず図６に示すように、図５の帯域制限画像
信号作成手段１２は、フィルタリング処理手段２０にお
いて原画像信号Ｓorgに対し、原画像信号Ｓorgを得た際
の主走査方向（Ｘ方向）および副走査方向（Ｙ方向）
（図８参照）に対してフィルタリング処理を施して原画
像信号Ｓorgよりも解像度が低い画像信号Ｌ1（以下、低
解像度画像信号という）を作成し、次にこの低解像度画
像信号Ｌ1に対して同様のフィルタリング処理を施して
この低解像度画像信号Ｌ1よりもさらに解像度が低い低
解像度画像信号Ｌ2を作成し、以降順次同様のフィルタ
リング処理を繰り返して各解像度毎の低解像度画像信号
Ｌk（ｋ＝１〜ｎ）を得るものである。そして、補間処
理手段２１において、このフィルタリング処理の各段に
おいて得られる低解像度画像信号Ｌkに対して、それぞ
れ２倍の画素数となるように補間処理を施して、鮮鋭度
が異なる複数のボケ画像信号Ｓus1〜Ｓusｎ（以下Ｓusk
（ｋ＝１〜ｎ）で代表させる）を得る。この後、減算器
２２により互いに対応する画素数を有する低解像度画像
信号Ｌk-1とボケ画像信号Ｓuskおよび原画像信号Ｓorg
とボケ画像信号Ｓus１との差分を求め、これを帯域制限
画像信号Ｂkとする。
【００４０】第２の実施形態においては、上記フィルタ
リング処理のフィルタとして、１次元ガウス分布に略対
応したフィルタを使用する。すなわちフィルタのフィル
タ係数を、ガウス信号に関する下記の式（６）にしたが
って定める。
リング処理のフィルタとして、１次元ガウス分布に略対
応したフィルタを使用する。すなわちフィルタのフィル
タ係数を、ガウス信号に関する下記の式（６）にしたが
って定める。
【００４１】
【数１】 このように、１次元ガウス分布に略対応したフィルタを
使用するのは、ガウス信号は周波数空間および実空間の
双方において局在性がよいためである。なお、上記
（６）式においてσ＝１とした場合の５×１の１次元フ
ィルタは図９に示すようなものとなる。
使用するのは、ガウス信号は周波数空間および実空間の
双方において局在性がよいためである。なお、上記
（６）式においてσ＝１とした場合の５×１の１次元フ
ィルタは図９に示すようなものとなる。
【００４２】フィルタリング処理は、図１０に示すよう
に、原画像信号Ｓorgに対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をＸ方向、Ｙ方向に行うことによ
り、低解像度画像信号Ｌ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｌk（ｋ＝
１〜ｎ）は、それぞれ画素数が原画像信号Ｓorgの１／
２2kの画像信号となる。
に、原画像信号Ｓorgに対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をＸ方向、Ｙ方向に行うことによ
り、低解像度画像信号Ｌ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｌk（ｋ＝
１〜ｎ）は、それぞれ画素数が原画像信号Ｓorgの１／
２2kの画像信号となる。
【００４３】次に、このようにして得られた低解像度画
像信号Ｌkに対して施される補間処理について説明す
る。補間処理を行うための補間演算の方法としては、Ｂ
スプラインによる方法等種々の方法が挙げられるが、本
実施形態においては、上記フィルタリング処理において
ガウス信号に基づくローパスフィルタを用いているた
め、補間演算についてもガウス信号を用いるものとす
る。具体的には、下記の式（７）において、σ＝２k-1
と近似したものを用いる。
像信号Ｌkに対して施される補間処理について説明す
る。補間処理を行うための補間演算の方法としては、Ｂ
スプラインによる方法等種々の方法が挙げられるが、本
実施形態においては、上記フィルタリング処理において
ガウス信号に基づくローパスフィルタを用いているた
め、補間演算についてもガウス信号を用いるものとす
る。具体的には、下記の式（７）において、σ＝２k-1
と近似したものを用いる。
【００４４】
【数２】 例えば低解像度画像信号Ｌ1を補間する際には、ｋ＝１
であるためσ＝１となる。この場合、補間処理を行うた
めのフィルタは、図１１に示すように５×１の１次元フ
ィルタとなる。この補間処理は、まず低解像度画像信号
Ｌ1に対して１画素おきに値が０の画素を１つずつ補間
することにより低解像度画像信号Ｌ1を原画像と同一の
画素数となるように拡大し、次に、この補間された低解
像度画像信号Ｌ1に対して上述した図１１に示す１次元
フィルタによりフィルタリング処理を施すことにより行
われる。
であるためσ＝１となる。この場合、補間処理を行うた
めのフィルタは、図１１に示すように５×１の１次元フ
ィルタとなる。この補間処理は、まず低解像度画像信号
Ｌ1に対して１画素おきに値が０の画素を１つずつ補間
することにより低解像度画像信号Ｌ1を原画像と同一の
画素数となるように拡大し、次に、この補間された低解
像度画像信号Ｌ1に対して上述した図１１に示す１次元
フィルタによりフィルタリング処理を施すことにより行
われる。
【００４５】同様に、この補間拡大処理を全ての低解像
度画像信号Ｌkに対して行う。低解像度画像信号Ｌkを補
間する際には、上記式（７）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｌkの各
画素の間に値が０の画素を１つずつ補間することによ
り、１段階高解像度の低解像度画像信号Ｌk-1と同一画
素数となるように拡大し、この値が０の画素が補間され
た低解像度画像信号Ｌkに対して３×２k−１の長さのフ
ィルタにより、フィルタリング処理を施すことにより補
間拡大してボケ画像信号Ｓuskを得る。
度画像信号Ｌkに対して行う。低解像度画像信号Ｌkを補
間する際には、上記式（７）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｌkの各
画素の間に値が０の画素を１つずつ補間することによ
り、１段階高解像度の低解像度画像信号Ｌk-1と同一画
素数となるように拡大し、この値が０の画素が補間され
た低解像度画像信号Ｌkに対して３×２k−１の長さのフ
ィルタにより、フィルタリング処理を施すことにより補
間拡大してボケ画像信号Ｓuskを得る。
【００４６】次に、上記のようにして作成されたボケ画
像信号Ｓuskが、対応する画素数を有する低解像度画像
信号Ｌk-1から減算されて、帯域制限画像信号Ｂk（ｋ＝
１〜ｎ）が得られる。なお、帯域制限画像信号Ｂkは下
記の式（８）に示すものとなる。
像信号Ｓuskが、対応する画素数を有する低解像度画像
信号Ｌk-1から減算されて、帯域制限画像信号Ｂk（ｋ＝
１〜ｎ）が得られる。なお、帯域制限画像信号Ｂkは下
記の式（８）に示すものとなる。
【００４７】Ｂ1＝Ｓorg−Ｓus１ Ｂ2＝Ｌ1−Ｓus２ Ｂ3＝Ｌ2−Ｓus３ ・ （８） ・ Ｂk＝Ｌk−Ｓusｋ 具体的には、図７に示すように５段階の低解像度画像信
号Ｌ1〜Ｌ5が得られた場合、まず最低解像度の低解像度
画像信号Ｌ5に対して補間処理を施して、低解像度画像
信号Ｌ4と同一画素数を有するボケ画像信号Ｓus５を作
成する。そして、低解像度画像信号Ｌ4からボケ画像信
号Ｓus５を減算して帯域制限画像信号Ｂ5を得る。以下
順次Ｌ3−Ｓus４、Ｌ2−Ｓus３、Ｌ1−Ｓus２、Ｓorg−
Ｓus１の演算を行って、帯域制限画像信号Ｂ1〜Ｂ5を得
る。ここで、最低解像度の低解像度画像信号Ｌk（Ｌ5）
は、原画像を縮小した低周波の情報を表すものである
が、これ以降の演算において使用することはない。
号Ｌ1〜Ｌ5が得られた場合、まず最低解像度の低解像度
画像信号Ｌ5に対して補間処理を施して、低解像度画像
信号Ｌ4と同一画素数を有するボケ画像信号Ｓus５を作
成する。そして、低解像度画像信号Ｌ4からボケ画像信
号Ｓus５を減算して帯域制限画像信号Ｂ5を得る。以下
順次Ｌ3−Ｓus４、Ｌ2−Ｓus３、Ｌ1−Ｓus２、Ｓorg−
Ｓus１の演算を行って、帯域制限画像信号Ｂ1〜Ｂ5を得
る。ここで、最低解像度の低解像度画像信号Ｌk（Ｌ5）
は、原画像を縮小した低周波の情報を表すものである
が、これ以降の演算において使用することはない。
【００４８】次に、上述したように算出された帯域制限
画像信号Ｂkを用いて行われる変換処理について説明す
る。図１２は周波数強調処理手段１３の構成を帯域制限
画像信号作成手段１２とともに示す概略ブロック図、図
１３は周波数強調処理を模式的に示す図である。図１２
に示すように、帯域制限画像信号作成手段１２において
作成された帯域制限画像信号Ｂkが、変換器３２におい
て変換関数ｆ1〜ｆnにより所望の大きさとなるように抑
制されて、変換済み信号（変換帯域制限画像信号）ｆk
Ｂk（ｋ＝１〜ｎ）が得られる。この変換関数ｆ1として
は、第１の実施形態と同様に図２に示す特性を有するも
のとを使用すればよい。なお、変換関数は、各帯域制限
画像信号において同一のものであってもよいが、各信号
毎に異なるものであってもよい。また、変換器３２が本
発明による変換手段と観念される。
画像信号Ｂkを用いて行われる変換処理について説明す
る。図１２は周波数強調処理手段１３の構成を帯域制限
画像信号作成手段１２とともに示す概略ブロック図、図
１３は周波数強調処理を模式的に示す図である。図１２
に示すように、帯域制限画像信号作成手段１２において
作成された帯域制限画像信号Ｂkが、変換器３２におい
て変換関数ｆ1〜ｆnにより所望の大きさとなるように抑
制されて、変換済み信号（変換帯域制限画像信号）ｆk
Ｂk（ｋ＝１〜ｎ）が得られる。この変換関数ｆ1として
は、第１の実施形態と同様に図２に示す特性を有するも
のとを使用すればよい。なお、変換関数は、各帯域制限
画像信号において同一のものであってもよいが、各信号
毎に異なるものであってもよい。また、変換器３２が本
発明による変換手段と観念される。
【００４９】そして、この変換済み信号ｆkＢkのうち、
最低解像度の変換済み信号ｆnＢnは高周波信号Ｓｎとさ
れるとともに、高周波信号Ｓｎが１段階高解像度の変換
済み信号ｆn-1Ｂn-1と同一画素数となるように、補間処
理手段３３において上記補間処理手段２１と同様に補間
処理がなされて、拡大高周波信号Ｓｎ′が得られる。こ
の後、変換済み信号ｆn-1Ｂn-1と拡大高周波信号Ｓｎ′
とが加算器３４において加算されて、高周波信号Ｓｎ−
１が得られる。そして、高周波信号Ｓｋ−１の補間拡大
による拡大高周波信号Ｓｋ−１′の取得、および拡大高
周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加算
による高周波信号Ｓｋ−２の取得を繰り返し行って、最
高解像度の高周波信号（加算変換信号）Ｓ１を得る。な
お、補間処理手段３３および加算機３４が、本発明によ
る加算手段と観念される。
最低解像度の変換済み信号ｆnＢnは高周波信号Ｓｎとさ
れるとともに、高周波信号Ｓｎが１段階高解像度の変換
済み信号ｆn-1Ｂn-1と同一画素数となるように、補間処
理手段３３において上記補間処理手段２１と同様に補間
処理がなされて、拡大高周波信号Ｓｎ′が得られる。こ
の後、変換済み信号ｆn-1Ｂn-1と拡大高周波信号Ｓｎ′
とが加算器３４において加算されて、高周波信号Ｓｎ−
１が得られる。そして、高周波信号Ｓｋ−１の補間拡大
による拡大高周波信号Ｓｋ−１′の取得、および拡大高
周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加算
による高周波信号Ｓｋ−２の取得を繰り返し行って、最
高解像度の高周波信号（加算変換信号）Ｓ１を得る。な
お、補間処理手段３３および加算機３４が、本発明によ
る加算手段と観念される。
【００５０】具体的には図１３に示すように、５段階の
帯域制限画像信号Ｂ1〜Ｂ5が得られた場合には、まず変
換済み信号ｆ1Ｂ1〜ｆ5Ｂ5が得られ、最低解像度の変換
済み信号ｆ5Ｂ5が高周波信号Ｓ５とされる。そして、高
周波信号Ｓ５に対して１段階高解像度の変換済み信号ｆ
4Ｂ4と同一画素数となるように補間処理が施されて拡大
高周波信号Ｓ５′が得られる。そして、変換済み信号ｆ
4Ｂ4と拡大高周波信号Ｓ５′とが加算されて、高周波信
号Ｓ４が得られる。以下同様にして高周波信号Ｓ３，Ｓ
２が得られ、最終的に最高解像度の高周波信号Ｓ１が得
られる。
帯域制限画像信号Ｂ1〜Ｂ5が得られた場合には、まず変
換済み信号ｆ1Ｂ1〜ｆ5Ｂ5が得られ、最低解像度の変換
済み信号ｆ5Ｂ5が高周波信号Ｓ５とされる。そして、高
周波信号Ｓ５に対して１段階高解像度の変換済み信号ｆ
4Ｂ4と同一画素数となるように補間処理が施されて拡大
高周波信号Ｓ５′が得られる。そして、変換済み信号ｆ
4Ｂ4と拡大高周波信号Ｓ５′とが加算されて、高周波信
号Ｓ４が得られる。以下同様にして高周波信号Ｓ３，Ｓ
２が得られ、最終的に最高解像度の高周波信号Ｓ１が得
られる。
【００５１】このようにして最高解像度の高周波信号Ｓ
１が得られると、演算器３５において、下記の式（９）
に示すように、強調係数β（Ｓorg）が乗じられて強調
変換信号β（Ｓorg）・Ｓ１が得られ、この強調変換信
号β（Ｓorg）・Ｓ１が原画像信号Ｓorgから減算されて
処理済み画像信号Ｓprocが得られる。なお、演算器３５
が、本発明による減算手段と観念される。
１が得られると、演算器３５において、下記の式（９）
に示すように、強調係数β（Ｓorg）が乗じられて強調
変換信号β（Ｓorg）・Ｓ１が得られ、この強調変換信
号β（Ｓorg）・Ｓ１が原画像信号Ｓorgから減算されて
処理済み画像信号Ｓprocが得られる。なお、演算器３５
が、本発明による減算手段と観念される。
【００５２】 Ｓproc＝Ｓorg −β（Ｓorg ）・Ｓ１ （９） （但し、Ｓproc：処理済み画像信号 Ｓorg ：原画像信号 β（Ｓorg）：強調係数） 次いで、第２の実施形態の動作について説明する。図１
４は第２の実施形態の動作を示すフローチャート、図１
５および図１６は第２の実施形態の動作を説明するため
の図である。なお、動作説明図においては、３段階のボ
ケ画像信号が作成されるものとして説明する。まず、読
取装置等から原画像信号Ｓorgが画像処理装置１０に入
力される（ステップＳ１１）。原画像信号Ｓorgは帯域
制限画像信号作成手段１２に入力されてここで原画像信
号Ｓorgの周波数帯域毎の周波数応答特性を表す帯域制
限画像信号Ｂkが作成される（ステップＳ１２）。ここ
で、原画像信号Ｓorgのエッジ部分におけるプロファイ
ルを図１５（ａ）に示す。このような原画像信号Ｓorg
から得られるボケ画像信号Ｓus１，Ｓus２，Ｓus３のプ
ロファイルはそれぞれ図１５（ｂ）、（ｃ），（ｄ）に
示すものとなる。１５（ａ）〜（ｃ）に示すようにボケ
画像信号Ｓus１〜Ｓus３においては、原画像信号Ｓorg
におけるエッジ部分がぼけたものとなる。原画像信号Ｓ
orgおよびボケ画像信号Ｓus１〜Ｓus３から得られる帯
域制限画像信号Ｂ1〜Ｂ3のプロファイルを図１５
（ｅ），（ｆ），（ｇ）に示す。ここで、Ｂ1＝Ｓorg−
Ｓus１、Ｂ2＝Ｓus１−Ｓus２、Ｂ3＝Ｓus２−Ｓus３で
ある。図１５（ｅ）〜（ｇ）に示すように、帯域制限画
像信号Ｂ1〜Ｂ3は、原画像のエッジに対応する部分に比
較的大きな値を有するプロファイルとなる。
４は第２の実施形態の動作を示すフローチャート、図１
５および図１６は第２の実施形態の動作を説明するため
の図である。なお、動作説明図においては、３段階のボ
ケ画像信号が作成されるものとして説明する。まず、読
取装置等から原画像信号Ｓorgが画像処理装置１０に入
力される（ステップＳ１１）。原画像信号Ｓorgは帯域
制限画像信号作成手段１２に入力されてここで原画像信
号Ｓorgの周波数帯域毎の周波数応答特性を表す帯域制
限画像信号Ｂkが作成される（ステップＳ１２）。ここ
で、原画像信号Ｓorgのエッジ部分におけるプロファイ
ルを図１５（ａ）に示す。このような原画像信号Ｓorg
から得られるボケ画像信号Ｓus１，Ｓus２，Ｓus３のプ
ロファイルはそれぞれ図１５（ｂ）、（ｃ），（ｄ）に
示すものとなる。１５（ａ）〜（ｃ）に示すようにボケ
画像信号Ｓus１〜Ｓus３においては、原画像信号Ｓorg
におけるエッジ部分がぼけたものとなる。原画像信号Ｓ
orgおよびボケ画像信号Ｓus１〜Ｓus３から得られる帯
域制限画像信号Ｂ1〜Ｂ3のプロファイルを図１５
（ｅ），（ｆ），（ｇ）に示す。ここで、Ｂ1＝Ｓorg−
Ｓus１、Ｂ2＝Ｓus１−Ｓus２、Ｂ3＝Ｓus２−Ｓus３で
ある。図１５（ｅ）〜（ｇ）に示すように、帯域制限画
像信号Ｂ1〜Ｂ3は、原画像のエッジに対応する部分に比
較的大きな値を有するプロファイルとなる。
【００５３】帯域制限画像信号Ｂkは、上記図２に示す
ような変換関数により変換されて変換済み信号ｆkＢkが
得られ（ステップＳ１３）、さらに変換済み信号ｆkＢk
の１段階高周波数帯域への補間処理による高周波信号Ｓ
ｋの取得、および高周波信号Ｓｋと対応する周波数帯域
の変換済み信号ｆkＢkとの加算による高周波信号Ｓｋ−
１の取得を最高周波数帯域の変換済み信号ｆ1Ｂ1まで繰
り返し行って、高周波信号Ｓ１を得る（ステップＳ１
４）。変換済み信号ｆＢ1〜ｆＢ3のプロファイルを図１
６（ａ），（ｂ），（ｃ）に示す。また、高周波信号Ｓ
１のプロファイルを図１６（ｄ）に示す。図１６（ａ）
〜（ｃ）に示すように、変換済み信号ｆＢ 1〜ｆＢ3は、
帯域制限画像信号Ｂ1〜Ｂ3において原画像のエッジに対
応する部分の信号が抑制されたプロファイルを有するも
のとなる。
ような変換関数により変換されて変換済み信号ｆkＢkが
得られ（ステップＳ１３）、さらに変換済み信号ｆkＢk
の１段階高周波数帯域への補間処理による高周波信号Ｓ
ｋの取得、および高周波信号Ｓｋと対応する周波数帯域
の変換済み信号ｆkＢkとの加算による高周波信号Ｓｋ−
１の取得を最高周波数帯域の変換済み信号ｆ1Ｂ1まで繰
り返し行って、高周波信号Ｓ１を得る（ステップＳ１
４）。変換済み信号ｆＢ1〜ｆＢ3のプロファイルを図１
６（ａ），（ｂ），（ｃ）に示す。また、高周波信号Ｓ
１のプロファイルを図１６（ｄ）に示す。図１６（ａ）
〜（ｃ）に示すように、変換済み信号ｆＢ 1〜ｆＢ3は、
帯域制限画像信号Ｂ1〜Ｂ3において原画像のエッジに対
応する部分の信号が抑制されたプロファイルを有するも
のとなる。
【００５４】そして高周波信号Ｓ１を用いて上記式
（９）に示す演算を行って処理済み画像信号Ｓprocを得
て（ステップＳ１５）、処理を終了する。処理済み画像
信号Ｓprocのプロファイルを図１６（ｅ）に示す。図１
６（ｅ）に示すように、処理済み画像信号Ｓprocにおい
ては、信号値が一定の平坦部における高周波ノイズが除
去されているとともに、エッジ部分のコントラストも保
持されたものとなる。
（９）に示す演算を行って処理済み画像信号Ｓprocを得
て（ステップＳ１５）、処理を終了する。処理済み画像
信号Ｓprocのプロファイルを図１６（ｅ）に示す。図１
６（ｅ）に示すように、処理済み画像信号Ｓprocにおい
ては、信号値が一定の平坦部における高周波ノイズが除
去されているとともに、エッジ部分のコントラストも保
持されたものとなる。
【００５５】なお、上記第１および第２の実施形態にお
いては、原画像信号Ｓorgに対して周波数強調処理を施
しているが、周波数強調処理が施された処理済み画像信
号Ｓorgに対してさらにダイナミックレンジ圧縮処理を
施してもよい。以下、これを第３の実施形態として説明
する。図１７は第３の実施形態の構成を示す概略ブロッ
ク図である。図１７に示すように、本発明の第３の実施
形態による画像処理装置４０は、上記第１または第２の
実施形態と同様に原画像信号Ｓorgに対して周波数強調
処理を施して周波数強調処理が施された処理済み画像信
号Ｓprocを得る周波数強調処理手段４１と、処理済み画
像信号Ｓprocに対してダイナミックレンジ圧縮処理を施
してダイナミックレンジ圧縮処理済み画像信号Ｓproc１
を得るダイナミックレンジ圧縮処理手段４２とを備え
る。
いては、原画像信号Ｓorgに対して周波数強調処理を施
しているが、周波数強調処理が施された処理済み画像信
号Ｓorgに対してさらにダイナミックレンジ圧縮処理を
施してもよい。以下、これを第３の実施形態として説明
する。図１７は第３の実施形態の構成を示す概略ブロッ
ク図である。図１７に示すように、本発明の第３の実施
形態による画像処理装置４０は、上記第１または第２の
実施形態と同様に原画像信号Ｓorgに対して周波数強調
処理を施して周波数強調処理が施された処理済み画像信
号Ｓprocを得る周波数強調処理手段４１と、処理済み画
像信号Ｓprocに対してダイナミックレンジ圧縮処理を施
してダイナミックレンジ圧縮処理済み画像信号Ｓproc１
を得るダイナミックレンジ圧縮処理手段４２とを備え
る。
【００５６】ダイナミックレンジ圧縮処理手段４２にお
いては、下記の式（１０）に示すように処理済み画像信
号Ｓprocに対してダイナミックレンジ圧縮処理が施され
る。
いては、下記の式（１０）に示すように処理済み画像信
号Ｓprocに対してダイナミックレンジ圧縮処理が施され
る。
【００５７】 Ｓproc１＝Ｓproc＋Ｄ（Ｓusa） （１０） 但し、Ｓusa：ボケ画像信号 Ｄ：ダイナミックレンジ圧縮係数 なお、ボケ画像信号Ｓusaは第１の実施形態において求
めたボケ画像信号Ｓusと同一のものであってもよい。
めたボケ画像信号Ｓusと同一のものであってもよい。
【００５８】なお、下記の式（１１）に示すように、処
理済み画像信号Ｓprocに対してダイナミックレンジ圧縮
処理を施すようにしてもよい。
理済み画像信号Ｓprocに対してダイナミックレンジ圧縮
処理を施すようにしてもよい。
【００５９】 Ｓproc１＝Ｓproc＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）） （１１） Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）＝｛ｆd1(Ｓorg−
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓusk(k=1〜n）：第２の実施形態と同様にして
求めたボケ画像信号 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） さらに、ダイナミックレンジ圧縮処理済み画像信号Ｓpr
oc１に対して階調変換処理を施すようにしてもよい。以
下、これを第４の実施形態として説明する。図１８は第
４の実施形態の構成を示すブロック図である。図１８に
示すように、本発明の第４の実施形態による画像処理装
置４０′は、第３の実施形態と同様に周波数強調処理手
段４１およびダイナミックレンジ圧縮処理手段４２を備
え、さらにダイナミックレンジ圧縮処理手段４２におい
て得られたダイナミックレンジ圧縮処理済み画像信号Ｓ
proc１に対して階調変換処理を施して階調変換処理済み
画像信号Ｓproc２を得る階調処理手段４３を備えたもの
である。階調処理手段４３においては、下記の式（１
２）に示すように、ダイナミックレンジ圧縮処理済み画
像信号Ｓproc１に対して階調変換処理が施される。
Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋…＋ｆdk(Ｓusk-1−Ｓu
sk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓusk(k=1〜n）：第２の実施形態と同様にして
求めたボケ画像信号 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 （ＤはＳorg−Ｆdrcを変換する関数）） さらに、ダイナミックレンジ圧縮処理済み画像信号Ｓpr
oc１に対して階調変換処理を施すようにしてもよい。以
下、これを第４の実施形態として説明する。図１８は第
４の実施形態の構成を示すブロック図である。図１８に
示すように、本発明の第４の実施形態による画像処理装
置４０′は、第３の実施形態と同様に周波数強調処理手
段４１およびダイナミックレンジ圧縮処理手段４２を備
え、さらにダイナミックレンジ圧縮処理手段４２におい
て得られたダイナミックレンジ圧縮処理済み画像信号Ｓ
proc１に対して階調変換処理を施して階調変換処理済み
画像信号Ｓproc２を得る階調処理手段４３を備えたもの
である。階調処理手段４３においては、下記の式（１
２）に示すように、ダイナミックレンジ圧縮処理済み画
像信号Ｓproc１に対して階調変換処理が施される。
【００６０】Ｓproc２＝γ（Ｓproc１） （１２） 但し、γ：階調変換関数 なお、第４の実施形態においては、階調処理手段４３に
おいて行われる階調変換処理の程度に応じて、ダイナミ
ックレンジ圧縮処理手段４２において行われるダイナミ
ックレンジ圧縮処理の程度を変更することが好ましい。
すなわち、階調変換処理により得られる処理済み画像信
号Ｓproc２のコントラストが大きくなるほど、高濃度域
および低濃度域の画像が潰れてしまうため、潰れが生じ
ないように、式（１０）あるいは式（１１）におけるダ
イナミックレンジ圧縮係数の程度を変更することが好ま
しい。
おいて行われる階調変換処理の程度に応じて、ダイナミ
ックレンジ圧縮処理手段４２において行われるダイナミ
ックレンジ圧縮処理の程度を変更することが好ましい。
すなわち、階調変換処理により得られる処理済み画像信
号Ｓproc２のコントラストが大きくなるほど、高濃度域
および低濃度域の画像が潰れてしまうため、潰れが生じ
ないように、式（１０）あるいは式（１１）におけるダ
イナミックレンジ圧縮係数の程度を変更することが好ま
しい。
【図１】本発明の第１の実施形態による画像処理装置の
構成を示す概略ブロック図
構成を示す概略ブロック図
【図２】変換関数の一例を示す図
【図３】第１の実施形態の動作を示すフローチャート
【図４】第１の実施形態の動作を説明するための図
【図５】本発明の第２の実施形態による画像処理装置の
構成を示す概略ブロック図
構成を示す概略ブロック図
【図６】帯域制限画像信号作成処理の概要を示すブロッ
ク図
ク図
【図７】帯域制限画像信号作成処理を模式的に示す図
【図８】原画像の主走査方向および副走査方向を示す図
【図９】フィルタリング処理に使用されるフィルタの一
例を示す図
例を示す図
【図１０】低解像度画像信号作成処理の詳細を示す図
【図１１】補間処理に使用されるフィルタの一例を示す
図
図
【図１２】周波数強調処理手段の構成を帯域制限画像信
号作成手段とともに示す概略ブロック図
号作成手段とともに示す概略ブロック図
【図１３】変換処理を模式的に示す図
【図１４】第２の実施形態の動作を示すフローチャート
【図１５】第２の実施形態の動作を説明するための図
（その１）
（その１）
【図１６】第２の実施形態の動作を説明するための図
（その２）
（その２）
【図１７】本発明の第３の実施形態による画像処理装置
の構成を示す概略ブロック図
の構成を示す概略ブロック図
【図１８】本発明の第４の実施形態による画像処理装置
の構成を示す概略ブロック図
の構成を示す概略ブロック図
１，１０，４０ 画像処理装置 ２ ボケ画像信号作成手段 ３，２２ 減算器 ４，３２ 変換器 ５，３５ 演算器 １２ 帯域制限画像信号作成手段 １３ 周波数強調処理手段 ２０ フィルタリング処理手段 ２１ 補間処理手段 ３４ 加算器
Claims (15)
- 【請求項１】 原画像を表す原画像信号に対して、周
波数強調処理を施す画像処理方法において、 前記原画像信号から前記原画像の高周波成分を表す高周
波信号を作成し、 該高周波信号を、信号値が大きいほどその値を抑制する
変換関数により変換して変換信号を得、 該変換信号を前記原画像信号から減算して処理済み画像
信号を得ることを特徴とする画像処理方法。 - 【請求項２】 原画像を表す原画像信号に対して、周
波数強調処理を施す画像処理方法において、 前記原画像信号から複数の帯域制限画像信号を作成し、 該各帯域制限画像信号を、信号値が大きいほどその値を
抑制する変換関数により変換して複数の変換帯域制限画
像信号を得、 該複数の変換帯域制限画像信号を加算して加算変換信号
を得、 該加算変換信号を前記原画像信号から減算して処理済み
画像信号を得ることを特徴とする画像処理方法。 - 【請求項３】 前記処理済み画像信号に対して、さら
にダイナミックレンジ圧縮処理を施してダイナミックレ
ンジ圧縮処理済み画像信号を得ることを特徴とする請求
項１または２記載の画像処理方法。 - 【請求項４】 前記ダイナミックレンジ圧縮処理済み
画像信号に対して、さらに階調変換処理を施して階調変
換処理済み画像信号を得ることを特徴とする請求項３記
載の画像処理方法。 - 【請求項５】 前記所定の変換関数は、非線形関数で
あることを特徴とする請求項１から４のいずれか１項記
載の画像処理方法。 - 【請求項６】 原画像を表す原画像信号に対して、周
波数強調処理を施す画像処理装置において、 前記原画像信号から前記原画像の高周波成分を表す高周
波信号を作成する高周波信号作成手段と、 該高周波信号を、信号値が大きいほどその値を抑制する
変換関数により変換して変換信号を得る変換手段と、 該変換信号を前記原画像信号から減算して処理済み画像
信号を得る減算手段とを備えたことを特徴とする画像処
理装置。 - 【請求項７】 原画像を表す原画像信号に対して、周
波数強調処理を施す画像処理装置において、 前記原画像信号から複数の帯域制限画像信号を作成する
帯域制限画像信号作成手段と、 該各帯域制限画像信号を、信号値が大きいほどその値を
抑制する変換関数により変換して複数の変換帯域制限画
像信号を得る変換手段と、 該複数の変換帯域制限画像信号を加算して加算変換信号
を得る加算手段と、 該加算変換信号を前記原画像信号から減算して処理済み
画像信号を得る減算手段とを備えたことを特徴とする画
像処理装置。 - 【請求項８】 前記処理済み画像信号に対して、ダイ
ナミックレンジ圧縮処理を施してダイナミックレンジ圧
縮処理済み画像信号を得るダイナミックレンジ圧縮処理
手段をさらに備えたことを特徴とする請求項６または７
記載の画像処理装置。 - 【請求項９】 前記ダイナミックレンジ圧縮処理済み
画像信号に対して、階調変換処理を施して階調変換処理
済み画像信号を得る階調処理手段をさらに備えたことを
特徴とする請求項８記載の画像処理装置。 - 【請求項１０】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項６から９のいずれか１項
記載の画像処理装置。 - 【請求項１１】 原画像を表す原画像信号に対して、
周波数強調処理を施す画像処理方法をコンピュータに実
行させるためのプログラムを記録したコンピュータ読取
り可能な記録媒体において、 前記プログラムは、前記原画像信号から前記原画像の高
周波成分を表す高周波信号を作成する手順と、 該高周波信号を、信号値が大きいほどその値を抑制する
変換関数により変換して変換信号を得る手順と、 該変換信号を前記原画像信号から減算して処理済み画像
信号を得る手順とを有することを特徴とするコンピュー
タ読取り可能な記録媒体。 - 【請求項１２】 原画像を表す原画像信号に対して、
周波数強調処理を施す画像処理方法をコンピュータに実
行させるためのプログラムを記録したコンピュータ読取
り可能な記録媒体において、 前記プログラムは、前記原画像信号から複数の帯域制限
画像信号を作成する手順と、 該各帯域制限画像信号を、信号値が大きいほどその値を
抑制する変換関数により変換して複数の変換帯域制限画
像信号を得る手順と、 該複数の変換帯域制限画像信号を加算して加算変換信号
を得る手順と、 該加算変換信号を前記原画像信号から減算して処理済み
画像信号を得る手順とを有することを特徴とするコンピ
ュータ読取り可能な記録媒体。 - 【請求項１３】 前記処理済み画像信号に対して、さ
らにダイナミックレンジ圧縮処理を施してダイナミック
レンジ圧縮処理済み画像信号を得る手順を有することを
特徴とする請求項１１または１２記載のコンピュータ読
取り可能な記録媒体。 - 【請求項１４】 前記ダイナミックレンジ圧縮処理済
み画像信号に対して、さらに階調変換処理を施して階調
変換処理済み画像信号を得る手順を有することを特徴と
する請求項１３記載のコンピュータ読取り可能な記録媒
体。 - 【請求項１５】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項１１から１４のいずれか
１項記載のコンピュータ読取り可能な記録媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000257395A JP2002074356A (ja) | 2000-08-28 | 2000-08-28 | 画像処理方法および装置並びに記録媒体 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000257395A JP2002074356A (ja) | 2000-08-28 | 2000-08-28 | 画像処理方法および装置並びに記録媒体 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2002074356A true JP2002074356A (ja) | 2002-03-15 |
Family
ID=18745858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000257395A Withdrawn JP2002074356A (ja) | 2000-08-28 | 2000-08-28 | 画像処理方法および装置並びに記録媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2002074356A (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006106919A1 (ja) * | 2005-03-31 | 2006-10-12 | Nikon Corporation | 画像処理方法 |
| JP2007234034A (ja) * | 2007-03-30 | 2007-09-13 | Matsushita Electric Ind Co Ltd | 画像処理装置及び画像処理方法 |
| US7298917B2 (en) | 2002-11-11 | 2007-11-20 | Minolta Co., Ltd. | Image processing program product and device for executing Retinex processing |
| JP2009284550A (ja) * | 2003-07-01 | 2009-12-03 | Xerox Corp | 走査されたドキュメントのためのデジタル・スクリーニング解除技術 |
| US7630573B2 (en) | 2004-08-03 | 2009-12-08 | Fujifilm Corporation | Noise reduction apparatus and method |
| US7791656B2 (en) | 2005-08-16 | 2010-09-07 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
| US8244034B2 (en) | 2006-03-31 | 2012-08-14 | Nikon Corporation | Image processing method |
| US8369640B2 (en) | 2007-02-02 | 2013-02-05 | Nikon Corporation | Image processing for removing a noise component contained in an image |
| JP2020198540A (ja) * | 2019-06-03 | 2020-12-10 | 株式会社シグマ | 画像信号処理方法 |
| JP7550110B2 (ja) | 2021-06-18 | 2024-09-12 | 株式会社日立ハイテク | 放射線画像処理装置、画像処理方法、及びシステム |
-
2000
- 2000-08-28 JP JP2000257395A patent/JP2002074356A/ja not_active Withdrawn
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7298917B2 (en) | 2002-11-11 | 2007-11-20 | Minolta Co., Ltd. | Image processing program product and device for executing Retinex processing |
| JP2009284550A (ja) * | 2003-07-01 | 2009-12-03 | Xerox Corp | 走査されたドキュメントのためのデジタル・スクリーニング解除技術 |
| US7630573B2 (en) | 2004-08-03 | 2009-12-08 | Fujifilm Corporation | Noise reduction apparatus and method |
| WO2006106919A1 (ja) * | 2005-03-31 | 2006-10-12 | Nikon Corporation | 画像処理方法 |
| US8059910B2 (en) | 2005-03-31 | 2011-11-15 | Nikon Corporation | Image processing method for removing noise from an image |
| JPWO2006106919A1 (ja) * | 2005-03-31 | 2008-09-11 | 株式会社ニコン | 画像処理方法 |
| JP4535125B2 (ja) * | 2005-03-31 | 2010-09-01 | 株式会社ニコン | 画像処理方法 |
| US7990443B2 (en) | 2005-08-16 | 2011-08-02 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
| US7791656B2 (en) | 2005-08-16 | 2010-09-07 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
| US8244034B2 (en) | 2006-03-31 | 2012-08-14 | Nikon Corporation | Image processing method |
| US8369640B2 (en) | 2007-02-02 | 2013-02-05 | Nikon Corporation | Image processing for removing a noise component contained in an image |
| US8675982B2 (en) | 2007-02-02 | 2014-03-18 | Nikon Corporation | Image processing for removing a noise component contained in an image |
| JP2007234034A (ja) * | 2007-03-30 | 2007-09-13 | Matsushita Electric Ind Co Ltd | 画像処理装置及び画像処理方法 |
| JP2020198540A (ja) * | 2019-06-03 | 2020-12-10 | 株式会社シグマ | 画像信号処理方法 |
| JP7423033B2 (ja) | 2019-06-03 | 2024-01-29 | 株式会社シグマ | 画像信号処理方法 |
| JP7550110B2 (ja) | 2021-06-18 | 2024-09-12 | 株式会社日立ハイテク | 放射線画像処理装置、画像処理方法、及びシステム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3700804B2 (ja) | 画像処理方法および装置 | |
| JP3816151B2 (ja) | 画像処理方法および装置 | |
| US7280703B2 (en) | Method of spatially filtering a digital image using chrominance information | |
| JP4004562B2 (ja) | 画像処理方法および装置 | |
| JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP2001167264A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP2008511048A (ja) | 画像処理方法、及び画像処理用コンピュータソフトウェア | |
| JPH0944656A (ja) | 画像処理方法および装置 | |
| JP3783116B2 (ja) | 放射線画像強調処理方法および装置 | |
| JP2007018379A (ja) | 画像処理方法及び画像処理装置 | |
| JPH0944655A (ja) | 画像処理方法および装置 | |
| JPH0998293A (ja) | 画像処理方法および装置 | |
| JP2002183727A (ja) | 画像処理装置 | |
| JP4202395B2 (ja) | 画像変換方法、変換画像生成方法、および画像補正装置 | |
| JP2002319020A (ja) | 画像処理装置および画像処理方法 | |
| JP3669452B2 (ja) | 画像処理方法および装置 | |
| JPH1063838A (ja) | 画像処理方法および装置 | |
| JP4014671B2 (ja) | 多重解像度変換方法および装置 | |
| JP2000307868A (ja) | 画像処理方法および装置 | |
| JP2001056856A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2007042124A (ja) | ノイズ除去方法および装置並びにプログラム | |
| Cherifi et al. | Color contrast enhancement method using steerable pyramid transform |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20071106 |