JP2002170114A - 画像処理方法および装置並びに記録媒体 - Google Patents

Abstract

(57)【要約】 【課題】 より細かな周波数帯域に帯域分割した帯域制
限画像信号を得て、より細かな周波数帯域での画像処理
を行う。 【解決手段】 ボケ画像信号作成手段１において、正方
状画素配列の原画像信号Ｓorgに正方市松縮小処理であ
る補間処理ａを施して縮小画像信号Ｓ１を得、これに市
松正方縮小処理である補間処理ｃを施して縮小画像信号
Ｓ２を得る。この処理を繰り返して複数の縮小画像信号
Ｓｋを得、縮小画像信号Ｓｋを原画像と同一サイズに補
間拡大してボケ画像信号Ｓusｋを得る。隣接するボケ画
像信号Ｓusｋ間および原画像信号Ｓorgとボケ画像信号
Ｓus１との差分値を帯域制限画像信号とする。帯域制限
画像信号のうち少なくとも１つの信号について変換処理
を施し、変換された帯域制限画像信号を積算して積算信
号を作成する。そして、積算信号に所定の強調係数を乗
じて原画像信号Ｓorgと加算することにより、原画像の
高周波成分が強調された処理済み画像信号Ｓprocを得
る。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は画像信号に対して、
この画像信号の高周波成分を強調する非鮮鋭マスク処理
（いわゆるボケマスク処理）を行う画像処理方法および
装置並びに画像処理方法をコンピュータに実行させるた
めのプログラムを記録したコンピュータ読取り可能な記
録媒体に関するものである。

【０００２】

【従来の技術】画像を表す画像信号を得、この画像信号
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば、放射線画像の診
断性能を向上させるために、画像信号に対して非鮮鋭マ
スク（ボケマスク）を用いて周波数強調処理を施す方法
が本出願人により提案されている（特開昭55-163472
号、特開昭55-87953号等）。この周波数強調処理は、読
み取った原画像信号Ｓorgに、この原画像信号Ｓorgから
非鮮鋭マスク画像信号Ｓusを減算することにより原画像
の高周波成分に関する信号を得、この高周波成分に関す
る信号に強調係数βを乗じたものを原画像信号Ｓorgに
加える処理を施すもので、これにより、画像において所
定の空間周波数成分を強調するようにしたものである。
これを式で表すと下記の式（１）のようになる。

【０００３】 Ｓproc＝Ｓorg＋β×（Ｓorg−Ｓus） （１） （Ｓproc：周波数強調処理された信号、Ｓorg：原画像
信号、Ｓus：非鮮鋭マスク画像信号、β：強調係数） ここで、非鮮鋭マスク画像信号Ｓusは、例えば画像を構
成する画素の１画素おきに、各画素の周囲のＭ×Ｎの範
囲内の原画像信号Ｓorgについて、 Ｓus＝ΣＳorg／（Ｍ×Ｎ） （２） なる演算を行うことにより求められる。

【０００４】また、特開平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-75395号に記載された処理に用いられ
る非鮮鋭マスク画像信号は、まず原画像信号の画素に対
して所定間隔毎に所定のフィルタリング処理を施すこと
によって画素数を縦横１／２となるように間引きし、そ
のようにして得た画像信号に対して同様のフィルタリン
グ処理を繰り返してさらに画素数を１／２^2kとなるよう
に段階的に少なくした縮小画像信号を複数作成し、縮小
画像信号のそれぞれに対して、所定の補間方法により原
画像と画素数が同一となるように補間処理を施すことに
より作成される。したがって、非鮮鋭マスク画像信号は
画素数は原画像信号と同一であるが、原画像信号よりも
鮮鋭度が低い画像を表す画像信号となる。

【０００６】また、帯域制限画像信号は、例えば隣接す
る周波数帯域の非鮮鋭マスク画像信号同士で差分を求め
る、あるいは原画像信号と各非鮮鋭マスク画像信号の差
分を求めることにより作成される。したがって、帯域制
限画像信号は、それにより表される画像の画素数が原画
像と同一となり、原画像信号の周波数帯域毎の周波数応
答特性を表す信号となる。

【０００７】

【発明が解決しようとする課題】上記帯域制限画像信号
は、原画像の１／２^2Nのサイズの画像を表すものである
ため、各帯域制限画像信号の周波数帯域は原画像信号の
ナイキスト周波数の１／２^N間隔の周波数帯域を表すも
のとなる。このため、ナイキスト周波数の１／２^Nの周
波数間隔で原画像の周波数処理を行うことが可能であ
る。その一方で、ナイキスト周波数の１／２^Nの周波数
帯域よりも細かい周波数帯域に帯域分割して処理を行う
ことが望まれている。

【０００８】本発明は上記事情に鑑みなされたものであ
り、より細かな周波数帯域に帯域分割した帯域制限画像
信号を得て、より細かな周波数帯域での画像処理を行う
ことができる画像処理方法および装置並びに画像処理方
法をコンピュータに実行させるためのプログラムを記録
したコンピュータ読取り可能な記録媒体を提供すること
を目的とするものである。

【０００９】

【課題を解決するための手段】本発明の画像処理方法
は、原画像を表す原画像信号に、該原画像の高周波成分
に関する信号を加算することにより、該原画像の高周波
成分を強調する画像処理方法において、前記原画像信号
に対して正方市松縮小処理および市松正方縮小処理に基
づく処理を施して、互いに周波数応答特性が異なる複数
の非鮮鋭マスク画像信号を作成し、前記原画像信号およ
び前記複数の非鮮鋭マスク画像信号、または前記複数の
非鮮鋭マスク画像信号に基づいて、前記原画像の複数の
周波数帯域毎の画像を表す複数の帯域制限画像信号を作
成し、該複数の帯域制限画像信号のうち少なくとも１つ
の帯域制限画像信号に対して、所定の変換処理を施して
少なくとも１つの変換画像信号を作成し、該変換画像信
号に基づいて前記原画像信号に加算する前記高周波成分
に関する信号を得ることを特徴とするものである。

【００１０】「正方市松縮小処理」とは、正方状画素配
列の画像信号に対して、隣接する画素列の各画素を交互
に間引く、あるいは上下左右に隣接する画素の画素値に
基づくフィルタリング処理を行う等により、画素数を減
少させつつ画素配列を正方状から市松状に変更する処理
のことをいう。この正方市松縮小処理により、処理後の
画像信号は市松状画素配列となり、画素数が処理前の正
方状画素配列の画像信号の略１／２となる。

【００１１】「市松正方縮小処理」とは、市松状画素配
列の画像信号に対して、隣接する画素列の各画素を交互
に間引く、あるいは対角方向に隣接する画素の画素値に
基づくフィルタリング処理を行う等により、画素値を減
少させつつ画素配列を市松状から正方状に変更する処理
のことをいう。この市松正方縮小処理により、処理後の
画像信号は正方状配列となり、画素数が処理前の市松状
配列の画像信号の略１／２となる。

【００１２】ここで、「正方状画素配列」とは、画像信
号により表される画像の各画素が図４７に示すように、
最も近接する画素が縦横等間隔に配列されてなるものを
いう。

【００１３】「市松状画素配列」とは、画像信号により
表される画像の各画素が図４８に示すように、最も近接
する画素が斜め方向に等間隔に配列され、ある画素列の
各画素間にその画素列に隣接する画素列の各画素が位置
せしめられるような配列のことをいう。

【００１４】「正方市松縮小処理および市松正方縮小処
理に基づく処理」とは、下記の処理のことをいう。ま
ず、原画像信号が正方状画素配列の場合、原画像信号に
対して正方市松縮小処理を施して、市松状画素配列の縮
小画像信号を作成し、さらにこの縮小画像信号に対して
市松正方縮小処理を施して正方状画素配列の縮小画像信
号を作成する。そして、得られた縮小画像信号に対して
正方市松縮小処理および市松正方縮小処理を交互に施す
ことにより、画素数を段階的に少なくした複数の縮小画
像信号を得る。一方、原画像信号が市松状画素配列の場
合、原画像信号に対して市松正方縮小処理を施して、市
松状画素配列の縮小画像信号を作成し、さらにこの縮小
画像信号に対して正方市松縮小処理を施して市松状画素
配列の縮小画像信号を作成する。そして、得られた縮小
画像信号に対して市松正方縮小処理および正方市松縮小
処理を交互に施すことにより、画素数を段階的に少なく
した複数の縮小画像信号を得る。そして、各縮小画像信
号に対して、スプライン補間演算などの補間方法によ
り、原画像と画素数が同一となるように補間演算を施す
処理である。

【００１５】なお、補間演算を行う際には、縮小画像信
号の画素配列に応じて補間係数が配置された補間フィル
タを使用する。例えば、縮小画像信号が正方状画素配列
である場合には、正方状に補間係数が配置された補間フ
ィルタを使用し、縮小画像信号が市松状画素配列である
場合には、市松状に補間係数が配置された補間フィルタ
を使用する。また、補間フィルタとしては例えばＢスプ
ライン補間演算を行う補間フィルタを使用することがで
きる。この場合、市松状に補間係数が配列された補間フ
ィルタは、正方状に補間係数が配列された補間フィルタ
の補間係数を４５度回転させた補間係数を有するものを
用いることができる。

【００１６】そして、「正方市松縮小処理および市松正
方縮小処理に基づく処理」により、互いに周波数応答特
性が異なる複数の非鮮鋭マスク画像信号が作成される。

【００１７】「変換画像信号に基づいて原画像信号に加
算する高周波成分に関する信号を得る」とは、１つの帯
域制限画像信号に対してのみ所定の変換処理を施した場
合は、１つの変換画像信号のみが得られるため、この１
つの変換画像信号を高周波成分に関する信号とすること
をいう。一方、２以上の帯域制限画像信号に対して所定
の変換処理を施した場合には、２以上の変換画像信号が
得られるため、これらの変換画像信号を加算することに
より、高周波成分に関する信号を得ることをいう。

【００１８】なお、本発明による画像処理方法において
は、前記所定の変換処理は、前記帯域制限画像信号の信
号値の少なくとも一部を小さくする処理であることが好
ましい。

【００１９】また、前記所定の変換処理は、前記帯域制
限画像信号の絶対値が所定の閾値より大きい場合に、該
信号値の絶対値を小さくする処理としてもよい。

【００２０】この際、前記帯域制限画像信号の信号値の
絶対値が前記所定の閾値よりも小さい他の閾値よりも小
さいほど、該信号値の絶対値を小さくする処理としても
よい。

【００２１】また、前記帯域制限画像信号の周波数帯域
に応じて該帯域制限画像信号の信号値の絶対値の大きさ
を変化させる処理としてもよい。

【００２２】また、前記所定の変換処理は、前記帯域制
限画像信号を、該帯域制限画像信号の周波数帯域に応じ
て異なる関数に基づいて、該帯域制限画像信号の信号値
の絶対値に基づいて定まる該絶対値以下の値となるよう
に変換する処理であってもよい。なお、複数の帯域制限
画像信号を変換する場合には、各帯域制限画像信号の周
波数帯域毎に異なる関数が用いられる。

【００２３】この際、前記関数は、前記帯域制限画像信
号の信号値の絶対値が所定の値より大きい場合は、前記
変換画像信号の信号値が略一定の値となるように前記帯
域制限画像信号を変換する関数であり、前記所定の値
は、前記帯域制限画像信号が高周波数帯域であるほど小
さい値であることが望ましい。

【００２４】さらに、前記関数は、前記帯域制限画像信
号が低周波帯域であるほど、該帯域制限画像信号におけ
る絶対値が０近傍の所定の範囲内の信号値を変換するこ
とにより得られる前記変換画像信号の絶対値がより小さ
い値となるように、前記帯域制限画像信号を変換する関
数であることが望ましい。

【００２５】なお、前記帯域制限画像信号の作成、前記
変換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算は、具体的には、下記の式にしたがって行って
もよい。

【００２６】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓo
rg,Ｓ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）：原画像信号に基づいて定められる強調係
数） あるいは、下記の式にしたがって行ってもよい。

【００２７】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓo
rg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ1
（Ｓorg−Ｓus1)＋ｆ2（Ｓorg−Ｓus2)＋…＋ｆk(Ｓorg
−Ｓusk)＋…＋ｆN（Ｓorg−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） なお、原画像信号が正方画素配列の場合、ｋが奇数の非
鮮鋭マスク画像信号は市松画素配列の縮小画像信号に対
して補間演算を施すことにより得られ、ｋが偶数の非鮮
鋭マスク画像信号は正方画素配列の縮小画像信号に対し
て補間演算を施すことにより得られる。一方、原画像信
号が市松画素配列の場合、ｋが奇数の非鮮鋭マスク画像
信号は正方画素配列の縮小画像信号に対して補間演算を
施すことにより得られ、ｋが偶数の非鮮鋭マスク画像信
号は市松画素配列の縮小画像信号に対して補間演算を施
すことにより得られる。

【００２８】また、前記帯域制限画像信号の作成、前記
変換画像信号の作成、前記高周波成分に関する信号の作
成、および該高周波成分に関する信号の前記原画像信号
への加算を、下記の式にしたがって行ってもよい。

【００２９】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓu
s1,Ｓus2,…ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝｛ｆ2（Ｓus1−Ｓus2)
＋ｆ3（Ｓus2−Ｓus3)＋…＋ｆk(Ｓusk-1−Ｓusk)＋…
＋ｆN（ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） あるいは、下記の式にしたがって行ってもよい。

【００３０】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓu
s1,Ｓus2,…ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ2（Ｓus1
−Ｓus2)＋ｆ3（Ｓus1−Ｓus3)＋…＋ｆk(Ｓus1−Ｓus
k)＋…＋ｆN（Ｓus1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらには、上記各式における所定の変換処理は、前記強
調係数に応じた変換処理としてもよい。この際、この強
調係数に応じた変換処理は、前記帯域制限画像信号の周
波数帯域に応じて異なる変換処理であって、該帯域制限
画像信号の周波数帯域が低周波帯域であるほど、該帯域
制限画像信号を大きく抑制する処理であり、前記強調係
数が大きいほど、高周波帯域の帯域制限画像信号に対す
る抑制の度合いと、低周波帯域の帯域制限画像信号に対
する抑制の度合いの差が大きいことが望ましい。

【００３１】さらに、前記所定の変換処理は、前記帯域
制限画像信号を、該帯域制限画像信号の信号値の絶対値
に基づいて定まる該絶対値以下の値となるように変換す
ることにより抑制画像信号を作成し、前記原画像信号お
よび前記複数の非鮮鋭マスク画像信号、または前記複数
の非鮮鋭マスク画像信号に基づいて、前記抑制画像信号
の作成に使用された帯域制限画像信号よりも低い周波数
帯域の信号を含む補助画像信号を作成し、該補助画像信
号を、該補助画像信号の信号値の絶対値が小さいほど１
に近く、大きいほど０に近い値となるように変換するこ
とにより前記抑制画像信号にそれぞれ対応する倍率信号
を作成し、前記抑制画像信号に、該抑制画像信号に対応
する前記倍率信号を乗ずることにより前記変換画像信号
を得る処理としてもよい。

【００３２】これは、具体的には、前記帯域制限画像信
号の作成、前記変換画像信号の作成、前記高周波成分に
関する信号の作成、および該高周波成分に関する信号の
前記原画像信号への加算を、下記の式にしたがって行っ
てもよい。

【００３３】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓo
rg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
Ｓus1)・ｇ（Ｓus1−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
（Ｓus2−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓu
sk−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN
−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） あるいは、下記の式にしたがって行ってもよい。

【００３４】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓo
rg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
Ｓus1)・ｇ（Ｓorg−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
（Ｓorg−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓo
rg−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg
−ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） にしたがって行ってもよい。

【００３５】さらに、前記所定の変換処理は、前記少な
くとも１つの帯域制限画像信号である被変換帯域制限画
像信号よりも低い周波数帯域の帯域制限画像信号である
低周波側帯域制限画像信号を、原点を通り該原点におけ
る傾きが略０で、処理される値が大きくなるにつれて該
傾きが漸増する非線形関数に基づいて変換することによ
り前記被変換帯域制限画像信号の補助画像信号を作成
し、該補助画像信号を前記被変換帯域制限画像信号に加
算することにより複合帯域制限画像信号を作成し、該複
合帯域制限画像信号を、該複合帯域制限画像信号の信号
値の絶対値に基づいて定まる該絶対値以下の値となるよ
うに変換することにより前記変換画像信号を得る処理と
してもよい。

【００３６】具体的には、前記帯域制限画像信号の作
成、前記変換画像信号の作成、前記高周波成分に関する
信号の作成、および該高周波成分に関する信号の前記原
画像信号への加算を、下記の式にしたがって行ってもよ
い。

【００３７】Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓo
rg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝［ｆ1｛（Ｓorg
−Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝＋ｆ2｛（Ｓus1−Ｓus
2)＋ｇ（Ｓus2−Ｓus3）｝＋…＋ｆk｛（Ｓusk-1−Ｓus
k)＋ｇ（Ｓusk−Ｓusk+1）｝＋…＋ｆN｛（ＳusN-1−Ｓ
usN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：複合帯域制限画像信号を変換するための
関数 ｇ：帯域制限画像信号を変換して補助画像信号を作成す
るための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） また、上記所定の変換処理は、全て、前記原画像を撮影
した際の撮影部位に応じて、前記帯域制限画像信号を変
換することにより行うことが望ましい。

【００３８】本発明による画像処理装置は、原画像を表
す原画像信号に、該原画像の高周波成分に関する信号を
加算することにより、該原画像の高周波成分を強調する
画像処理装置において、前記原画像信号に対して正方市
松縮小処理および市松正方縮小処理に基づく処理を施し
て、互いに周波数応答特性が異なる複数の非鮮鋭マスク
画像信号を作成する非鮮鋭マスク画像信号作成手段と、
前記原画像信号および前記複数の非鮮鋭マスク画像信
号、または前記複数の非鮮鋭マスク画像信号に基づい
て、前記原画像の複数の周波数帯域毎の画像を表す複数
の帯域制限画像信号を作成する帯域制限画像信号作成手
段と、該複数の帯域制限画像信号のうち少なくとも１つ
の帯域制限画像信号に対して、所定の変換処理を施して
少なくとも１つの変換画像信号を作成する変換手段と、
該変換画像信号に基づいて前記原画像信号に加算する前
記高周波成分に関する信号を得る周波数強調処理手段と
からなることを特徴とするものである。

【００３９】なお、本発明による画像処理装置において
は、変換手段において行われる上記所定の変換処理は、
本発明による画像処理方法と同様の処理とすればよい。

【００４０】また、本発明による画像処理装置において
は、前記変換手段が、前記帯域制限画像信号を、該帯域
制限画像信号の信号値の絶対値に基づいて定まる該絶対
値以下の値となるように変換することにより抑制画像信
号を作成する抑制画像信号作成手段と、前記原画像信号
および前記複数の非鮮鋭マスク画像信号、または前記複
数の非鮮鋭マスク画像信号に基づいて、前記抑制画像信
号の作成に使用された帯域制限画像信号よりも低い周波
数帯域の信号を含む補助画像信号を作成する補助画像信
号作成手段と、該補助画像信号を、該補助画像信号の信
号値の絶対値が小さいほど１に近く、大きいほど０に近
い値となるように変換することにより前記抑制画像信号
にそれぞれ対応する倍率信号を作成する倍率信号作成手
段と、前記抑制画像信号に、該抑制画像信号に対応する
前記倍率信号を乗ずることにより前記変換画像信号を作
成する乗算手段とを備えてなるものとしてもよい。

【００４１】さらに、本発明による画像処理装置におい
ては、前記変換手段が、前記少なくとも１つの帯域制限
画像信号である被変換帯域制限画像信号よりも低い周波
数帯域の帯域制限画像信号である低周波側帯域制限画像
信号を、原点を通り該原点における傾きが略０で、処理
される値が大きくなるにつれて該傾きが漸増する非線形
関数に基づいて変換することにより前記被変換帯域制限
画像信号の補助画像信号を作成する補助画像信号作成手
段と、該補助画像信号を前記被変換帯域制限画像信号に
加算することにより複合帯域制限画像信号を作成する複
合帯域制限画像信号作成手段と、該複合帯域制限画像信
号を、該複合帯域制限画像信号の信号値の絶対値に基づ
いて定まる該絶対値以下の値となるように変換すること
により前記変換画像信号を作成する変換画像信号作成手
段とを備えてなるものとしてもよい。

【００４２】なお、本発明による画像処理方法をコンピ
ュータに実行させるためのプログラムとしてコンピュー
タ読取り可能な記録媒体に記録して提供してもよい。

【００４３】

【発明の効果】本発明による画像処理方法および装置に
よれば、原画像信号に対して正方市松縮小処理および市
松正方縮小処理に基づく処理が施されて、互いに周波数
応答特性が異なる複数の非鮮鋭マスク画像信号が作成さ
れ、この非鮮鋭マスク画像信号から前記原画像の複数の
周波数帯域毎の画像を表す複数の帯域制限画像信号が作
成される。ここで、帯域制限画像信号は、原画像の濃度
変化が比較的小さいいわゆる平坦部においては、各周波
数帯域において値は小さくなるものである。これに対し
て、濃度が急激に変化するエッジ部近傍においては、帯
域制限画像信号が比較的低周波数帯域である場合、すな
わち非鮮鋭マスク画像信号を得る際の非鮮鋭マスクのサ
イズが比較的大きい場合は、そのエッジ部近傍の画素に
ついては、非鮮鋭マスクにエッジ部が含まれてしまうた
め、帯域制限画像信号がエッジ部の影響を受けて信号値
の絶対値が比較的大きなものとなる。このように、本来
エッジ部でない部分がエッジ部の濃度値の影響を受ける
ことにより、画像処理を行うことにより得られた画像の
エッジ部にオーバーシュート、アンダーシュートなどの
アーチファクトが発生してしまう。

【００４４】本発明はこの点に鑑みてなされたものであ
り、帯域制限画像信号のうち少なくとも１つの信号につ
いて、例えば帯域制限画像信号の少なくとも一部を小さ
くするような変換処理を施して変換画像信号を作成し、
この変換画像信号に基づいて原画像信号に加算するため
の高周波成分に関する信号を得るようにしたものであ
る。このため、比較的信号値の絶対値の大きい帯域制限
画像信号は、原画像信号に加算するための高周波成分に
関する信号に対する影響力が小さくなり実質的の非鮮鋭
マスクのサイズを小さくしたものと同様の信号となる。
そしてこれにより、濃度が急激に変化するエッジ部近傍
においても、アーチファクトの原因となる信号は影響力
が弱められるため、処理を施すことにより得られる画像
をアーチファクトのない良好なものとすることができ
る。

【００４５】一方、上述した特開平10-75395号に記載さ
れた処理においては、非鮮鋭マスク画像信号を作成する
際の画素の間引きは１／２とされるため、帯域制限画像
信号は、原画像信号のナイキスト周波数の１／２^N〜１
／２^N+1の周波数帯域の画像を表すものとなっている。
一方、本発明において得られる帯域制限画像信号は、原
画像信号のナイキスト周波数の１／２^N〜１／２^kおよび
１／２^k〜１／２^N+1（Ｎ＜ｋ＜Ｎ＋１）の周波数帯域の
画像を表すものとなっている。例えば、最高周波数帯域
および最高周波数帯域より１段階低周波数帯域の帯域制
限画像信号を第１および第２の帯域制限画像信号とする
と、第２の帯域制限画像信号は原画像の１／２²サイズ
の画像を表すものであるが、第１の帯域制限画像信号
は、原画像と原画像の１／２²サイズの画像との中間サ
イズの画像を表すものとなっている。よって、原画像信
号のナイキスト周波数をｆｓ／２とすると、第１の帯域
制限画像信号はｆｓ／２〜ｆｓ／ｋ、第２の帯域制限画
像信号はｆｓ／ｋ〜ｆｓ／４（２＜ｋ＜４）の周波数帯
域の画像を表すものとなる。このため、本発明によれ
ば、従来の方法と比較して、原画像信号をより細かい周
波数帯域に帯域分割した帯域制限画像信号を得ることが
できる。

【００４６】したがって、少なくとも１つの帯域制限画
像信号に対して、所定の変換処理を施すことにより、よ
り細かな周波数帯域の間隔で変換処理を施すことができ
るため、より詳細に画質を制御することができる。した
がって、得られた変換画像信号に基づいて高周波成分に
関する信号を得、これを原画像信号に加算することによ
り、より高画質の処理済み画像を再生可能な処理済み画
像信号を得ることができる。

【００４７】また、上記特開平10-75395号に記載された
処理に用いられる帯域制限画像信号は、正方状画素配列
の原画像信号を正方状に帯域分割して正方状配列の非鮮
鋭マスク画像信号を得、この非鮮鋭マスク画像信号から
帯域制限画像信号を得るものであるため、とくに所定の
変換処理において急峻なエッジを抑制したい場合、縦横
方向すなわち正方状に画素が並ぶ方向のエッジについて
は良好にエッジを抑制できる。しかしながら、帯域制限
画像信号は斜め方向すなわち市松状に画素が並ぶ方向に
帯域分割されていないため、斜め方向については良好に
エッジを抑制することができなかった。これに対し、本
発明により得られる帯域制限画像信号は、原画像が正方
状および市松状に帯域分割されているため、斜め方向に
帯域分割された市松状画素配列の帯域制限画像信号を用
いれば、斜め方向についても良好にエッジを抑制するこ
とができることとなる。したがって、エッジの方向によ
らずエッジが過度に強調されることを防止することがで
きる。

【００４８】さらに、原画像信号を取得する画像取得装
置によっては、原画像の鮮鋭度が縦横方向と斜め方向と
で鮮鋭度が異なる場合がある。例えば、デジタルカメラ
においては、撮影時のノイズを除去するために光学的な
ローパスフィルタを使用しているが、このローパスフィ
ルタの特性が得られた画像信号を再生することにより得
られる再生画像の縦横方向および斜め方向で異なる場合
がある。ここで、画像取得装置の特性が、原画像の縦横
方向の鮮鋭度が斜め方向よりも高いものとすると、上記
特開平10-75395号に記載された方法では、帯域制限画像
信号は縦横方向にのみ帯域分割されているため、斜め方
向についてのみ鮮鋭度を高くすることができなかった。
これに対して、本発明は帯域制限画像信号は斜め方向に
も帯域分割されているため、斜め方向についてのみ鮮鋭
度を高くすることができる。したがって、画像取得装置
の特性に拘わらず、高画質の画像を再現可能な処理済み
画像信号を得ることができる。

【００４９】なお、上記所定の変換処理として、少なく
とも１つの帯域制限画像信号の信号値の絶対値が所定の
閾値より大きい場合に、帯域制限画像信号の信号値の絶
対値を小さくする処理とすることにより、とくに大きな
信号値の影響力を弱めることができる。

【００５０】また、少なくとも１つの帯域制限画像信号
の信号値の絶対値が所定の閾値よりも小さい他の閾値よ
りも小さいほどこの帯域制限画像信号の信号値の絶対値
を小さくする処理とすれば、画像中においてノイズと見
なせるような信号値の絶対値の小さい成分のレスポンス
を小さくすることができ、これにより得られる画像のノ
イズを低減することができる。

【００５１】また、帯域制限画像信号の周波数帯域に応
じて帯域制限画像信号の信号値の絶対値の大きさを変化
させる処理とすることにより、周波数帯域に応じた処理
を行うことも可能となる。

【００５２】また、所定の閾値よりも大きいか否かによ
らず、周波数帯域に応じて異なる関数によって変換処理
を施せば、より適切な処理を施すことが可能となり、こ
れにより、最終的に得られる処理済み画像信号の全体の
周波数応答特性を自在にコントロールすることもでき
る。これは、上記アーチファクトのみならず、周波数帯
域の境界に発生する階段状のアーチファクトをも抑制で
きるという効果を有する。

【００５３】あるいは、所定の変換処理を施す帯域制限
画像信号よりも低い周波数帯域の画像信号を用いて、２
種類の信号を加算あるいは乗算するなどして変換画像信
号を作成することにより、上記階段状のアーチファクト
をさらに抑制し、より滑らかな画像信号を作成して良好
な処理済み画像を得ることができる。

【００５４】さらに、原画像を得た際の撮影部位に応じ
てこのような変換を行うことにより、各撮影部位に適し
た高周波成分の強調を行うことができる。

【００５５】

【発明の実施の形態】以下、図面を参照して本発明の実
施形態について説明する。図１は本発明の第１の実施形
態による画像処理装置の構成を示す概略ブロックであ
る。図１に示すように第１の実施形態による画像処理装
置は、入力された原画像信号Ｓorgに対して正方市松縮
小処理および市松正方縮小処理に基づく処理を施すこと
により、互いに周波数応答特性が異なる多重解像度のボ
ケ画像信号（非鮮鋭マスク画像信号）Ｓusｋ（ｋ＝１〜
Ｎ）を作成するボケ画像信号作成手段１と、ボケ画像信
号作成手段１において作成されたボケ画像信号Ｓusｋに
基づいて複数の帯域制限画像信号を作成する帯域制限画
像信号作成手段２と、帯域制限画像信号作成手段２にお
いて作成された複数の帯域制限画像信号のうち少なくと
も１つの帯域制限画像信号について、帯域制限画像信号
の少なくとも一部を小さくするような変換処理を施して
変換画像信号を得る変換手段３と、変換画像信号を積算
して積算信号（高周波成分に関する信号）を作成する積
算手段４と、積算信号に強調係数を乗じて原画像信号Ｓ
orgと加算することにより、原画像の高周波成分が強調
された処理済み画像信号Ｓprocを得る周波数強調処理手
段５とを備える。

【００５６】はじめに、ボケ画像信号作成手段１におい
て行われる処理について説明する。図２はボケ画像信号
作成手段１において行われる処理を示す概略ブロック図
である。ここで、原画像信号Ｓorgにより表される原画
像は、図３に示すように正方状に画素が配列されてなる
ものとする。なお、図３および以降の図において○の位
置が信号値を有する画素の位置を表すものとする。

【００５７】まず、原画像信号Ｓorgに対して正方市松
縮小処理である補間処理ａが施される。この補間処理ａ
は、図４に示すフィルタにより原画像信号Ｓorgの各画
素を１画素おきにフィルタリング処理することにより行
われる。なお、このフィルタリング処理は、隣接する画
素列においては画素を１画素ずつずらして行われる。こ
れにより、図５に示す○位置の画素にのみフィルタリン
グ処理が施され、○位置にのみ信号値を有する、画素が
市松状に配列された縮小画像を表す縮小画像信号Ｓ１が
得られる。なお、フィルタリング処理に代えて単純に１
画素おきに画素を間引く処理を行ってもよい。

【００５８】この縮小画像信号Ｓ１は原画像信号Ｓorg
と比較して画素配列が市松状となり、画素数が略１／２
となっている。なお、「略１／２」としたのは、補間処
理ａは画素数を１／２とする処理であるが、画素列にお
ける画素数が奇数である場合には、画素数を１／２とし
ても１画素余りが生じ、完全に１／２とはならないため
である。

【００５９】さらに、上記補間処理ａにより求められた
縮小画像信号Ｓ１に対して市松正方縮小処理である補間
処理ｃが施される。この補間処理ｃは、図６に示すフィ
ルタにより縮小画像信号Ｓ１の各画素列を１列おきにフ
ィルタリング処理することにより行われる。これによ
り、図７に示す○位置の画素にのみフィルタリング処理
が施され、○位置にのみ信号値を有する、画素が正方状
に配列された縮小画像を表す縮小画像信号Ｓ２が得られ
る。なお、フィルタリング処理に代えて単純に画素を間
引く処理を行ってもよい。

【００６０】この縮小画像信号Ｓ２は縮小画像信号Ｓ１
と比較して画素配列が正方状となり、画素数が略１／２
となっている。また、この縮小画像信号Ｓ２は原画像信
号Ｓorgと比較して縦横の画素数が１／２すなわち画素
数が１／４となっている。したがって、縮小画像信号Ｓ
２により表される画像のサイズは原画像の１／２²とな
っている。

【００６１】さらに、上記補間処理ｃにより求められた
縮小画像信号Ｓ２に対して正方市松縮小処理である補間
処理ａを施すことにより、図８に示す○位置の画素にの
みフィルタリング処理が施され、○位置にのみ信号値を
有する、画素が市松状に配列された縮小画像を表す縮小
画像信号Ｓ３が得られる。

【００６２】そして補間処理ａ，ｃを繰り返すことによ
り、Ｎ個の縮小画像信号Ｓｋ（ｋ＝１〜Ｎ）が得られ
る。ここで、ｋが奇数の場合、縮小画像信号Ｓｋにより
表される縮小画像は市松状の画素配列となる。一方、ｋ
が偶数の場合、縮小画像信号Ｓｋにより表される縮小画
像は正方状の画素配列となる。また、最低周波数帯域の
縮小画像信号ＢＮを補間処理ｃまで行うことにより得た
場合、その縮小画像信号ＳＮにより表される縮小画像は
市松状の画素配列となる。一方、最低周波数帯域の縮小
画像信号ＳＮを補間処理ａまで行うことにより得た場
合、その縮小画像信号ＳＮにより表される縮小画像は正
方状の画素配列となる。

【００６３】また、ｋが偶数の縮小画像信号Ｓｋは、原
画像に対して１／２^kサイズの画像を表すものとなって
いる。一方、ｋが奇数の縮小画像信号は、原画像に対す
る１／２^2kサイズと１／２^2(k+1)サイズとの中間サイズ
の画像を表すものとなっている。ここで、縮小画像信号
Ｓｋの周波数応答特性を図９に示す。図９に示すよう
に、縮小画像信号Ｓｋのレスポンスはｋが大きいほど高
周波成分が除去されたものとなっている（但し、図９に
おいてはｋ＝１〜３としている）。なお、図９において
ｆｓ／２は原画像信号Ｓorgのナイキスト周波数であ
る。

【００６４】次いで、このようにして得られた縮小画像
信号Ｓｋに対して、補間処理ｂ，ｄが施される。ここ
で、補間処理ｂは補間処理ａにより得られた市松状画素
配列の縮小画像信号Ｓｋ（ｋ：奇数）に対して、補間処
理ｄは補間処理ｃにより得られた正方状画素配列の縮小
画像信号Ｓｋ（ｋ：偶数）に対して施される。まず、簡
便のため補間処理ｄについて説明する。この補間処理ｄ
は、正方状画素配列の縮小画像信号Ｓｋに対してＢスプ
ライン補間演算を行う補間フィルタによるフィルタリン
グ処理を施して、原画像と同一サイズのボケ画像を表す
ボケ画像信号Ｓusｋを得る処理である。

【００６５】以下、Ｂスプライン補間演算処理について
説明する。Ｂスプライン補間演算は、比較的鮮鋭度の低
い滑らかな２次画像を再生可能な補間画像信号を得るた
めの補間演算方法である。なお、ここでは簡便のため１
次元のＢスプラインについて説明する。このＢスプライ
ン補間演算は、元のサンプル点（画素）を通ることは必
要とされない代わりに、第１階微分係数および第２階微
分係数（ｆ″（Ｘ）で表す）が各区間間で連続すること
が必要とされる。

【００６６】すなわち、 ｆ_k（ｘ）＝Ａ_kｘ³＋Ｂ_kｘ²＋Ｃ_kｘ＋Ｄ_k （４） において（式（４）においてＢ_kは便宜上用いられる係
数であり、フィルタリング処理により得られる画像とは
異なる。）、 ｆ_k′（Ｘ_k）＝ｆ_k-1′（Ｘ_k） （５） ｆ_k′（Ｘ_k+1）＝ｆ_k+1′（Ｘ_k+1） （６） ｆ_k″（Ｘ_k）＝ｆ_k-1″（Ｘ_k） （７） ｆ_k″（Ｘ_k+1）＝ｆ_k+1″（Ｘ_k+1） （８） が条件となる。但し、画素Ｘ_kにおける第１階微分係数
が、その画素Ｘ_kの前後の画素であるＸ_k-1とＸ_k+1とに
ついて、これらの画像信号Ｙ_k-1、Ｙ_k+1の勾配（Ｙ _k+1
−Ｙ_k-1）／（Ｘ_k+1−Ｘ_k-1）に一致することが条件で
あるから、下記式（９）を満たす必要がある。

【００６７】 ｆ_k′（Ｘ_k）＝（Ｙ_k+1−Ｙ_k-1）／（Ｘ_k+1−Ｘ_k-1） （９） 同様に、画素Ｘ_k+1における第１階微分係数が、その画
素Ｘ_k+1の前後の画素であるＸ_kとＸ_k+2とについて、こ
れらの画像信号Ｙ_k、Ｙ_k+2の勾配（Ｙ_k+2−Ｙ_k）／（Ｘ
_k+2−Ｘ_k）に一致することが条件であるから、下記式
（１０）を満たす必要がある。

【００６８】 ｆ_k′（Ｘ_k+1）＝（Ｙ_k+2−Ｙ_k）／（Ｘ_k+2−Ｘ_k） （１０） また関数ｆ（Ｘ）は一般に下記式（１１）に示すもので
近似される。

【００６９】 ｆ（Ｘ）＝ｆ(0)＋ｆ′(0)Ｘ＋｛ｆ″(0)/２｝Ｘ² （１１） ここで、各区間Ｘ_k-2〜Ｘ_k-1，Ｘ_k-1〜Ｘ_k，Ｘ_k〜
Ｘ_k+1，Ｘ_k+1〜Ｘ_k+2の間隔（格子間隔という）を１と
し、画素Ｘ_kからの画素Ｘ_k+1方向への補間点Ｘ_pの位置
をｔ（０≦ｔ≦１）とすれば、式（５）〜（８）および
（１１）より、 ｆ_k′（０）＝Ｃ_k＝（Ｙ_k+1−Ｙ_k-1）／２ ｆ_k′（１）＝３Ａ_k＋２Ｂ_k＋Ｃ_k＝（Ｙ_k+2−Ｙ_k）／２ ｆ_k″（０）＝Ｙ_k+1−２Ｙ_k＋Ｙ_k-1＝２Ｂ したがって、 Ａ_k＝（Ｙ_k+2−３Ｙ_k+1＋３Ｙ_k−Ｙ_k-1）／６ Ｂ_k＝（Ｙ_k+1−２Ｙ_k＋Ｙ_k-1）／２ Ｃ_k＝（Ｙ_k+1−Ｙ_k-1）／２ ここで、Ｄ_kは未知のため、 Ｄ_k＝（Ｄ₁Ｙ_k+2＋Ｄ₂Ｙ_k+1＋Ｄ₃Ｙ_k＋Ｄ₄Ｙ_k-1）／６ とおく。また、スプライン補間関数ｆ_k（ｘ）は上述の
通り、Ｘ＝ｔなる変数変換をしているため、 ｆ_k（ｘ）＝ｆ_k（ｔ） となる。よって、 ｆ_k（ｔ）＝｛（Ｙ_k+2−３Ｙ_k+1＋３Ｙ_k−Ｙ_k-1）／
６｝ｔ³＋｛（Ｙ_k+1−２Ｙ_k＋Ｙ_k-1）／２｝ｔ²＋
｛（Ｙ_k+1−Ｙ_k-1）／２｝ｔ＋（Ｄ₁Ｙ_k+2＋Ｄ₂Ｙ_k+1＋
Ｄ₃Ｙ_k＋Ｄ₄Ｙ_k-1）／６ となり、これを画像信号Ｙ_k-1，Ｙ_k，Ｙ_k+1，Ｙ_k+2につ
いて整理すると、下記式（１２）で表すことができる。

【００７０】 ｆ_k（ｔ）＝｛（−ｔ³＋３ｔ²−３ｔ＋Ｄ₄）／６｝Ｙ_k-1 ＋｛（３ｔ³−６ｔ²＋Ｄ₃）／６｝Ｙ_k ＋｛（−３ｔ³＋３ｔ²＋３ｔ＋Ｄ₂）／６｝Ｙ_k+1 ＋｛（ｔ³＋Ｄ₁）／６｝Ｙ_k+2 （１２） ここで、ｔ＝１とおけば、 ｆ_k（１）＝｛（Ｄ₄−１）／６｝Ｙ_k-1＋｛（Ｄ₃−３）
／６｝Ｙ_k＋｛（Ｄ₂＋３）／６｝Ｙ_k+1＋｛（Ｄ₁＋１）
／６｝Ｙ_k+2 次に区間Ｘ_k+1〜Ｘ_k+2についての式（９）は、 ｆ_k+1（ｔ）＝｛（−ｔ³＋３ｔ²−３ｔ＋Ｄ₄）／６｝Ｙ_k ＋｛（３ｔ³−６ｔ²＋Ｄ₃）／６｝Ｙ_k+1 ＋｛（−３ｔ³＋３ｔ²＋３ｔ＋Ｄ₂）／６｝Ｙ_k+2 ＋｛（ｔ³＋Ｄ₁）／６｝Ｙ_k+3 （１３） ここで、ｔ＝０とおけば、 ｆ_k+1（０）＝（Ｄ₄／６）Ｙ_k＋（Ｄ₃／６）Ｙ_k+1＋
（Ｄ₂／６）Ｙ_k+2＋（Ｄ₁／６）Ｙ_k+3 連続性の条件（ｆ_k（１）＝ｆ_k+1（０））、および各フ
ィルタリング処理画像信号に対応する係数同士が等しい
という条件により、Ｄ₄−１＝０，Ｄ₃−３＝Ｄ ₄，Ｄ₂＋
３＝Ｄ₃，Ｄ₁＋１＝Ｄ₂，Ｄ₁＝０、となり、したがっ
て、 Ｄ_k＝（Ｙ_k+1＋４Ｙ_k＋Ｙ_k-1）／６ となる。よって、 Ｙ_p＝ｆ_k（ｔ）＝｛（−ｔ³＋３ｔ²−３ｔ＋１）／６｝Ｙ_k-1 ＋｛（３ｔ³−６ｔ²＋４）／６｝Ｙ_k ＋｛（−３ｔ³＋３ｔ²＋３ｔ＋１）／６｝Ｙ_k+1 ＋｛ｔ³／６｝Ｙ_k+2 （ １４） したがって、フィルタリング処理画像信号Ｙ_k-1、Ｙ_k、
Ｙ_k+1、Ｙ_k+2にそれぞれ対応する補間係数ｂ_k-1、ｂ_k、
ｂ_k+1、ｂ_k+2は、 ｂ_k-1＝（−ｔ³＋３ｔ²−３ｔ＋１）／６ ｂ_k ＝（３ｔ³−６ｔ²＋４）／６ ｂ_k+1＝（−３ｔ³＋３ｔ²＋３ｔ＋１）／６ ｂ_k+2＝ｔ³／６ （１５） のように、ｔの関数として求められることとなる。

【００７１】そして、上記式（１５）により１次元の補
間係数を求め、１次元の補間係数を２次元に展開して補
間フィルタを作成し、この補間フィルタにより縮小画像
信号Ｓｋに対してフィルタリング処理を施すことによ
り、×位置における画素値を求め、原画像と同一サイズ
のボケ画像信号を得ることができる。

【００７２】ここで、信号値が既知の画素位置における
補間信号値を求めるための１次元のＢスプライン補間演
算の補間係数は、ｔ＝０であることから１／６，４／
６，１／６となるため、これを２次元に展開すると図１
０に示すように３×３の補間係数の行列からなる補間フ
ィルタとなる。一方、信号値が既知の画素位置の間の位
置における補間信号値を求めるための１次元のＢスプラ
イン補間演算の補間係数は、ｔ＝１／２であることか
ら、１／４８，４７／４８，４７／４８，１／４８とな
るため、これを２次元に展開すると図１１に示すように
３×３の補間係数の行列からなる補間フィルタとなる。

【００７３】したがって、図７に示す縮小画像信号Ｓ２
に対して、○位置においては図１０に示す補間フィルタ
により、×位置においては図１１に示す補間フィルタに
よりフィルタリング処理を施すことにより、○位置およ
び×位置の双方の信号値を求めることができ、これによ
り原画像と同一画素数を有するボケ画像を表すボケ画像
信号Ｓus２を得ることができる。なお、図１０および図
１１に示す補間フィルタによりフィルタリング処理を行
う際には、補間係数と対応する位置に信号値を有する画
素が位置するよう、補間係数の間に０の値を挿入すれば
よい。

【００７４】一方、補間処理ｂは市松状画素配列の縮小
画像信号Ｓｋ（ｋ：奇数）に対して補間処理を施すもの
である。ここで、正方状画素配列は画素が縦横方向に最
も小さい間隔で均等に並ぶが、市松状画素配列は画素が
斜め方向に最も小さい間隔で均等に並ぶものである。よ
って、補間処理ｂにおいて○位置における補間信号値を
求める補間フィルタは図１２に示すように、図１０に示
す補間フィルタの補間係数を４５度回転させたものを用
いる。また、×位置における補間信号値を求める補間フ
ィルタは図１３に示すように、図１１に示す補間フィル
タの補間係数を４５度回転させたものとすればよい。な
お、図１２および図１３に示す補間フィルタによりフィ
ルタリング処理を行う際には、補間係数と対応する位置
に信号値を有する画素が位置するよう、補間係数の間に
０の値を挿入すればよい。

【００７５】このように、Ｂスプライン補間演算を行う
補間フィルタによるフィルタリング処理を各縮小画像信
号Ｓｋに対して施すことにより、各縮小画像信号に対応
したボケ画像信号Ｓusｋを得ることができる。

【００７６】このようにして得られたボケ画像信号Ｓus
ｋの周波数応答特性を図１４に示す。図１４に示すよう
に、ボケ画像信号Ｓusｋのｋの値が大きくなるほど、原
画像信号Ｓorgの高周波成分が除去された信号となる。

【００７７】図１５は上記図２のボケ画像信号作成手段
１を含む本発明の第１の実施形態による画像処理装置の
構成を詳細に示すブロック図である。図１５に示すよう
に、ボケ画像信号作成手段１において作成された各ボケ
画像信号Ｓusｋは、次に、帯域制限画像信号作成手段２
および変換手段３において処理される。まず、原画像信
号Ｓorgおよびボケ画像信号作成手段１において作成さ
れた複数のボケ画像信号Ｓusｋに基づいて帯域制限画像
信号Ｂｋ（ｋ＝１〜Ｎ）が作成されるが、この帯域制限
画像信号Ｂｋは減算器２１により、互いに隣接する周波
数帯域同士のボケ画像信号Ｓusｋの減算を行うことによ
り得られる。すなわち、Ｓorg−Ｓus１、Ｓus１−Ｓus
２、…ＳusN-1−ＳusNを順次計算することにより、複数
の帯域制限画像信号Ｂｋが求められる。この帯域制限画
像信号Ｂｋの周波数応答特性を図１６に示す。図１６に
示すように、帯域制限画像信号Ｂｋ（ここではｋ＝１〜
３）はｋの値が大きくなるほど、原画像信号Ｓorgの低
周波数帯域を表す信号となる。

【００７８】次いで、変換手段３において、このように
して求められた帯域制限画像信号Ｂｋをこの帯域制限画
像信号Ｂｋの信号値に応じて変換する。この変換は変換
器２２において、例えば図１７に示すような関数ｆを用
いて行う。なお、この関数ｆは上記Ｂスプライン補間演
算の説明に用いた関数ｆ_k（Ｘ）とは異なる。この関数
ｆは、帯域制限画像信号Ｂｋの信号値の絶対値が閾値Ｔ
ｈ１よりも小さい場合は傾きが１であり、閾値Ｔｈ１よ
りも大きい場合は傾きが１よりも小さくなるような非線
形関数である。この関数ｆは、全ての帯域制限画像信号
Ｂｋにおいて同一のものであってもよいが、帯域制限画
像信号Ｂｋの周波数帯域毎に異なるものであってもよ
い。なお、図１５では帯域制限画像信号Ｂｋの周波数帯
域毎に異なる関数ｆｋ（ｋ＝１〜Ｎ）を使用するものと
している。

【００７９】このような関数ｆにより帯域制限画像信号
Ｂｋを変換することにより得られた変換信号Ｂｋ′は、
上述した積算手段４および周波数強調処理手段５を内包
する演算器２３に入力される。この演算器２３において
は下記のような処理が行われる。まず、変換画像信号Ｂ
ｋ′が積算されて積算信号が得られる。そしてこの積算
信号が求められると、周波数強調処理手段５において原
画像信号Ｓorgの値に応じた強調係数βが乗じられ、さ
らにこの強調係数βが乗じられた積算信号が原画像信号
Ｓorgと加算されて処理済み画像信号Ｓprocが得られ
る。

【００８０】以上の帯域制限画像信号作成手段２、変換
手段３、積算手段４および周波数強調処理手段５におい
て行われる処理を下記の式（１６）に示す。

【００８１】 Ｓ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）：原画像信号に基づいて定められる強調係
数） このようにして得られた処理済み画像信号Ｓprocは、例
えば図１８に示すような周波数応答特性を有するものと
なる。すなわち、上述した帯域制限画像信号Ｂｋは、原
画像の濃度変化が比較的小さいいわゆる平坦部において
は、各周波数帯域において信号値の絶対値は小さくなる
ものである。これに対して、濃度が急激に変化するエッ
ジ部近傍においては、帯域制限画像信号Ｂｋが比較的低
周波数帯域である場合、すなわちボケ画像信号Ｓusｋを
得る際のマスクのサイズが比較的大きい場合は、図４９
に示すように、そのエッジ部近傍の画素についておかれ
たマスクにエッジ部が含まれてしまうため、帯域制限画
像信号Ｂｋがエッジ部の影響を受けて信号値の絶対値が
比較的大きなものとなる。このように、本来エッジ部で
ない部分がエッジ部の濃度値の影響を受けることによ
り、画像処理を行うことにより得られた画像のエッジ部
にはオーバーシュート、アンダーシュートなどのアーチ
ファクトが発生してしまう。

【００８２】そこで、帯域制限画像信号Ｂｋの信号値の
絶対値が閾値Ｔｈ１よりも大きい場合に、上述した関数
ｆによりこの信号値の絶対値が小さくなるように帯域制
限画像信号Ｂｋして変換画像信号Ｂｋ′を得、この変換
画像信号Ｂｋ′を積算し、さらに所定の強調係数で強調
することにより、原画像信号Ｓorgに加算するための高
周波成分に関する信号を得るようにしたものである。

【００８３】これにより、図１８に示すように、エッジ
部が存在しない平坦部においては処理済み画像信号Ｓpr
ocの周波数応答特性は実線で示すようなものとなるが、
エッジ部近傍の領域においては、処理済み画像信号Ｓpr
ocは図１８の破線に示すように比較的低い周波数帯域の
レスポンスが低下されたような周波数応答特性を有する
ものとなる。これは、エッジ部近傍の領域においては、
ボケ画像信号（式（１）におけるＳus）を得る際のマス
クが実際のマスクよりも小さくされたことと同一の効果
を奏するものである。

【００８４】したがって、エッジ部近傍の領域に対応す
る比較的信号値の絶対値の大きい帯域制限画像信号Ｂｋ
は、原画像信号Ｓorgに加算するための高周波成分に関
する信号に対する影響力が小さくなる。このため、濃度
が急激に変化するエッジ部近傍においても、アーチファ
クトの原因となる信号は影響力が弱められることとな
り、処理を施すことにより得られる画像をアーチファク
トのない良好な画像とすることができる。

【００８５】図１９は上記特開平10-75395号に記載され
た処理に用いられる帯域制限画像信号と本実施形態の手
法により得られた帯域制限画像信号Ｂｋとの周波数応答
特性の相違を説明するための図であり、（ａ）は上記特
開平10-75395号に記載された処理に用いられる帯域制限
画像信号（以下従来の帯域制限画像信号とする）の周波
数応答特性を、（ｂ）は本実施形態の手法により得られ
た帯域制限画像信号Ｂｋの周波数応答特性を示す図であ
る。なお、図１９においてｆｓ／２は正方状画素配列の
原画像信号Ｓorgのナイキスト周波数であり、実線が各
帯域制限画像信号が再現可能な周波数帯域の区分すなわ
ち帯域制限画像信号の最大境界を表している。

【００８６】上記従来の帯域制限画像信号は、原画像信
号Ｓorgにより表される原画像の１／２^2Nのサイズの画
像を表すものとなることから、原画像信号Ｓ０は矩形状
に周波数帯域が帯域分割され、図１９（ａ）に示すよう
に、各帯域制限画像信号の周波数応答特性が矩形状とな
る。

【００８７】一方、本実施形態により得られる帯域制限
画像信号Ｂｋは、原画像の１／２^2Nのサイズの細部画像
とともに、１／２^2Nサイズと１／２^2(N+1)サイズとの中
間のサイズの細部画像を表すものとなる。このため、原
画像信号Ｓ０は矩形状および菱形状に交互に周波数帯域
が帯域分割され、図１９（ｂ）に示すように、各帯域制
限画像信号も矩形状および菱形状の周波数応答特性が交
互に現れるものとなる。

【００８８】ここで、本実施形態により得られる帯域制
限画像信号Ｂｋのうち、ｋが偶数となり矩形状の周波数
応答特性を有する帯域制限画像信号Ｂｋの周波数帯域の
最大境界は、従来の帯域制限画像信号の周波数帯域の最
大境界と同一となる。一方、ｋが奇数となり菱形状の周
波数応答特性となる帯域制限画像信号Ｂｋの周波数帯域
の最大境界は、従来の帯域制限画像信号の最大境界の間
に位置するものとなる。このため、本実施形態によれ
ば、従来の帯域制限画像信号と比較して、原画像信号Ｓ
０をより細かい周波数帯域に帯域分割した帯域制限画像
信号Ｂｋを得ることができる。具体的には、従来の帯域
制限画像信号は、原画像信号Ｓorgのナイキスト周波数
ｆｓ／２の１／２^N〜１／２^N+1の周波数帯域を表すもの
であるのに対し、本実施形態においては、従来の帯域制
限画像信号の１つの周波数帯域内に、原画像信号Ｓorg
のナイキスト周波数ｆｓ／２の１／２^N〜１／２^Mおよび
１／２^M〜１／２^N+1（Ｎ＜Ｍ＜Ｎ＋１）の２つの周波数
帯域の帯域制限画像信号が含まれるものとなる。

【００８９】したがって、本実施形態により得られた帯
域制限画像信号Ｂｋに対しては、より細かな周波数帯域
の間隔で上記変換処理を施すことができるため、より詳
細に画質を制御することができる。その結果、処理が施
された変換画像信号Ｂｋ′を積算して得られた積算画像
信号を原画像信号Ｓorgに加算することにより得られる
処理済み画像信号Ｓprocを、より高画質の処理済み画像
を再生可能なものとすることができる。

【００９０】また、上記従来の帯域制限画像信号は、原
画像信号Ｓ０を矩形状に帯域分割して得るものであり、
とくに上記図１７に示すような関数ｆを用いて急峻なエ
ッジを抑制したい場合、縦横方向すなわち正方状に画素
が並ぶ方向のエッジについては良好にエッジを抑制でき
る。しかしながら、従来の帯域制限画像信号は菱形状に
帯域分割されていないため、斜め方向については良好に
エッジを抑制することができず、その結果斜め方向につ
いてはエッジが過度に強調されてしまっていた。これに
対し、本実施形態により得られる帯域制限画像信号Ｂｋ
は、原画像が正方状および市松状に交互に帯域分割され
ているため、斜め方向に帯域分割された市松状配列の帯
域制限画像信号を用いれば、斜め方向についても良好に
エッジを抑制することができることとなる。したがっ
て、エッジの方向によらずエッジが過度に強調されるこ
とを防止することができる。

【００９１】なお、本実施形態においては、上記式（１
６）により処理済み画像信号Ｓprocを得るようにしてい
るが、下記式（１７）により処理済み画像信号Ｓprocを
得るようにしてもよい。式（１６）と式（１７）とで
は、帯域制限画像信号Ｂｋを得る際に、式（１６）にお
いては隣接する周波数帯域同士で減算を行っているが、
式（１７）においては、全ての周波数帯域のボケ画像信
号Ｓusｋと、原画像信号Ｓorgとで減算処理を行ってい
る点が異なるものである。式（１７）により得られる処
理済み画像信号Ｓprocの周波数応答特性を図２０に示
す。図２０に示すようにエッジ部が存在しない平坦部に
おいては、処理済み画像信号Ｓprocの周波数応答特性は
実線で示すようなものとなるが、エッジ部近傍の領域に
おいては、処理済み画像信号Ｓprocは図２０の破線に示
すように比較的低い周波数帯域のレスポンスが大きく低
下されたような周波数応答特性を有するものとなる。こ
れは、エッジ部近傍の領域においては、ボケ画像信号Ｓ
usｋを得る際のマスクが実際のマスクよりも小さくされ
たことと同一の効果を奏するものである。

【００９２】ここで、図１８と比較すると、図２０の場
合はエッジ部近傍の領域において、全周波数帯域に亘っ
てレスポンスが低下したものとなっている。このため、
式（１７）にしたがって処理済み画像信号Ｓprocを求め
る方が平坦部のレスポンスは低下されることなくエッジ
部近傍の領域のみのレスポンスが低下されるため、より
好ましいものである。

【００９３】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝（1/N)・｛ｆ1（Ｓorg−Ｓus1)＋ｆ2（Ｓorg−Ｓus2)＋… ＋ｆk(Ｓorg−Ｓusk)＋…＋ｆN（Ｓorg−ＳusN)｝ （１７） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：ボケ画像信号 ｆk(k=1〜N)：各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらに、上述した実施形態においては、図１７に示すよ
うに、閾値Ｔｈ１よりも信号値の絶対値が大きい場合に
この絶対値を小さくするような関数を用いて帯域制限画
像信号Ｂｋを変換しているが、例えば図２１に示すよう
に、帯域制限画像信号Ｂｋの信号値の絶対値が閾値Ｔｈ
１よりも大きい場合に、この絶対値を小さくするととも
に、閾値Ｔｈ２よりも小さい場合に、絶対値を小さくす
るような関数を用いるようにしてもよい。

【００９４】このように、帯域制限画像信号Ｂｋの絶対
値が閾値Ｔｈ１よりも小さい閾値Ｔｈ２よりも小さいほ
ど、この帯域制限画像信号Ｂｋの絶対値を小さくするよ
うに、帯域制限画像信号Ｂｋを変換することにより、画
像中においてノイズと見なせるような信号値の絶対値の
小さい成分のレスポンスを小さくすることができ、これ
によりノイズが低減された処理済み画像を表す処理済み
画像信号Ｓprocを得ることができる。

【００９５】また、下記式（１８）により処理済み画像
信号Ｓprocを得るようにしてもよい。上述した式（１
７）においては、帯域制限画像信号Ｂｋを得るために、
原画像信号Ｓorgからボケ画像信号Ｓusｋを減算するよ
うにしているが、式（１８）においては、原画像信号Ｓ
orgを使用することなく、ボケ画像信号Ｓus１からボケ
画像信号Ｓusｋ（ｋ＝２〜Ｎ）を減算して帯域制限画像
信号Ｂｋを作成しているものである。式（１８）におい
て得られる帯域制限画像信号Ｂｋの周波数応答特性を図
２２に示す。図２２に示すように式（１８）において得
られる帯域制限画像信号Ｂｋの周波数応答特性は、原画
像信号Ｓorgを使用していないことから、最高周波数帯
域に対応する周波数帯域が除去されたものとなってい
る。

【００９６】式（１７）により処理を行うものにおいて
は、原画像信号Ｓorgを用いていることから、画像中の
ノイズと見なせる高周波成分をも強調してしまう場合が
あるため、得られる処理済み画像にノイズが目立つもの
となってしまうことがある。これに対して式（１８）に
より処理を行うものにおいては、原画像信号Ｓorgを用
いて帯域制限画像信号Ｂｋを作成していないため、ノイ
ズと見なせる最高周波数帯域の帯域制限画像信号Ｂ１が
除去されたものとなる。したがって、上記変換処理を行
っても、ノイズと見なせる成分が強調されることがなく
なりより高画質の処理済み画像を再現可能な処理済み画
像信号Ｓprocを得ることができる。

【００９７】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN） ＝（1/N)・｛ｆ2（Ｓus1−Ｓus2)＋ｆ3（Ｓus1−Ｓus3)＋… ＋ｆk(Ｓus1−Ｓusk)＋…＋ｆN（Ｓus1−ＳusN)｝ （１８） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：ボケ画像信号 ｆk(k=2〜N)：各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらに、下記の式（１９）により処理を行うようにして
もよい。上述した式（１６）においては、原画像信号Ｓ
orgを用いて作成した帯域制限画像信号Ｂｋを変換処理
に使用しているのに対して、式（１９）においては、原
画像信号Ｓorgを用いた帯域制限画像信号Ｂ１（Ｓorg−
Ｓus1）を変換処理に使用しないものである。これによ
り図２２に示すように、最高周波数帯域の帯域制限画像
信号Ｂ１が除去されることとなるため、式（１８）によ
り処理を行うものと同様に、処理済み画像信号Ｓprocの
高周波成分が除去されてノイズが強調されることがなく
なり、より高画質の処理済み画像を再生可能な処理済み
画像信号Ｓprocを得ることができる。

【００９８】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ2（Ｓus1−Ｓus2)＋ｆ3（Ｓus2−Ｓus3)＋… ＋ｆk(Ｓusk-1−Ｓusk)＋…＋ｆN（ＳusN-1−ＳusN)｝ （１９） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：ボケ画像信号 ｆk(k=2〜N)：各帯域制限画像信号を変換する関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） なお、原画像が医療用放射線画像の場合、撮影部位に応
じて診断に必要な周波数帯域が異なるものである。例え
ば、肺と胃の二重造影の画像において、肺は比較的低周
波成分を強調するのが好ましく、胃は胃壁のヒダの観察
を行うため比較的高周波成分を強調することが好まし
い。一方、骨や人工骨等の金属を含む画像においては、
これらのエッジ部分が強調されすぎることによるアーチ
ファクトを防止する必要があるが、例えばマンモグラフ
ィのように骨などのエッジ部分を含まない画像に対して
は、アーチファクトは発生しにくく、骨を含む画像と同
様に周波数成分の強調を行うと、観察に必要な部分が逆
に見にくくなるという問題もある。したがって、原画像
を得る際の撮影部位、あるいは帯域制限画像信号Ｂｋの
周波数帯域に応じて変換処理を行う関数ｆを変化させる
ことが望ましい。

【００９９】例えば、図２３に示すように、骨を含む画
像においては関数ｆをＢのようにして帯域制限画像信号
Ｂｋの高周波成分を大きく抑制してエッジ部分にアーチ
ファクトを発生しにくくすることが好ましい。逆に、マ
ンモグラフィのように骨を含まない画像においては関数
ｆをＡのようにしてＢよりも高周波成分を強調し、略全
周波数帯域に亘って帯域制限画像信号Ｂｋの信号値の絶
対値が大きくなるように変換処理を行うことが好まし
い。

【０１００】また、胸部画像の場合、全周波数帯域に亘
って帯域制限画像信号Ｂｋが強調されることが好ましい
ため、帯域制限画像信号Ｂｋの周波数帯域に応じて図２
４に示す関数ｆｋのように変化させることが好ましい。
一方、胃の二重造影像の場合は帯域制限画像信号Ｂｋの
高周波成分を強調し、低周波成分を抑制した方が好まし
いため、帯域制限画像信号Ｂｋの周波数帯域に応じて図
２５に示す関数ｆｋのように変化させることが好まし
い。ここで、図中の「高」「低」はそれぞれ、高周波帯
域の帯域制限画像信号Ｂｋの変換に使用される関数、低
周波帯域の帯域制限画像信号Ｂｋの変換に使用される関
数を示している。

【０１０１】帯域制限画像信号Ｂｋを図２４に示す関数
ｆｋにしたがって変換することにより得られた処理済み
画像信号Ｓprocの周波数応答特性を図２６に、帯域制限
画像信号Ｂｋを図２５に示す関数ｆｋにしたがって変換
することにより得られた処理済み画像信号Ｓprocの周波
数応答特性を図２７にそれぞれ示す。図２４に示す関数
ｆｋを使用した場合、図２６に示すように略全周波数帯
域に亘ってレスポンスが強調されているのに対して、図
２５に示す関数ｆｋを使用した場合、高周波帯域のレス
ポンスが他の周波数帯域よりも大きく強調されているこ
とが分かる。

【０１０２】このように、胸部画像の場合に全周波数帯
域に亘ってレスポンスが強調されるように変換処理を施
すことにより、肺野がより観察し易いものとなる。一
方、胃の二重造影像において高周波帯域のレスポンスが
強調されるように変換処理を施すことにより、胃壁のヒ
ダがより観察し易いものとなる。そしてこのように、帯
域制限画像信号Ｂｋの周波数帯域に応じてあるいは撮影
部位に応じて帯域制限画像信号Ｂｋの信号値の絶対値を
変化させるように処理を行うことにより、撮影部位ある
いは周波数帯域に応じたより観察に適した画像を得るこ
とができることとなる。

【０１０３】このように、帯域制限画像信号Ｂｋを変換
する関数ｆｋを、帯域制限画像信号Ｂｋの周波数帯域毎
に異ならしめることにより、得られる処理済み画像信号
Ｓprocの周波数応答特性を任意の周波数応答特性となる
ように調整可能である。したがって、処理対象である画
像に求められる条件に応じて、適切な関数ｆｋの組み合
わせを設定することにより、処理済み画像信号Ｓprocの
周波数応答特性を調整することが可能となる。

【０１０４】以下、帯域制限画像信号Ｂｋを変換する関
数ｆｋを周波数帯域毎に異なるものとすることにより得
られる効果について詳細に説明する。図２８は、関数ｆ
ｋを全て同一の関数とし、全周波数帯域の帯域制限画像
信号Ｂｋを同一の関数を用いた変換処理により変換する
方法の問題点を示す図であり、画像の濃度が急激に変化
しているエッジ部の近傍における処理を段階的に示した
ものである。（ａ）はエッジ部近傍における階段状に信
号値が変化する原画像信号Ｓorgとその原画像信号Ｓorg
に基づいて作成されたボケ画像信号Ｓusｋを、（ｂ）は
（ａ）に対応する帯域制限画像信号Ｂｋを、（ｃ）は変
換画像信号Ｂｋ′を、（ｄ）は変換画像信号Ｂｋ′を積
算して得た積算信号をそれぞれ示す。

【０１０５】図２８（ｄ）に示すように、変換画像信号
Ｂｋ′を積算して得た積算信号は周波数帯域の境界部に
おいて不自然なつなぎ目ができてしまい、これがすじ状
のアーチファクトの原因となる。これを防止するために
は、この境界部ができるだけ自然につながるよう考慮し
て変換画像信号Ｂｋ′を作成しなければならない。しか
しながら、関数ｆｋが一意に定められている場合には各
帯域制限画像信号Ｂｋに対して境界部の影響を考慮した
任意の変換処理を行うことができず、結果としてすじ状
のアーチファクトの発生を防止することができなかっ
た。

【０１０６】したがって、変換処理に使用する関数ｆｋ
を帯域制限画像信号Ｂｋの周波数帯域毎に異なるものと
し、これらの関数ｆｋを周波数帯域の境界部を考慮して
設定することにより、すじ状のアーチファクトの発生を
防止するものである。

【０１０７】次に、このような帯域制限画像信号Ｂｋの
周波数帯域毎に異なる関数ｆｋについて例を示して説明
する。図２９は、変換画像信号Ｂｋ′の作成に使用され
る関数ｆｋの例を示す図であり、横軸に処理される帯域
制限画像信号Ｂｋの信号値を、縦軸に変換画像信号Ｂ
ｋ′の信号値をそれぞれ示している。これらの関数ｆｋ
は、帯域制限画像信号Ｂｋをその帯域制限画像信号Ｂｋ
の信号値の絶対値に基づいて定まる、その絶対値以下の
値となるように変換するもので、帯域制限画像信号Ｂｋ
の信号値の絶対値が所定の値より大きい帯域制限画像信
号Ｂｋについては、変換画像信号Ｂｋ′の値が略一定の
値となるように変換を行い、関数ｆｋが変換する帯域制
限画像信号Ｂｋの周波数帯域が高周波帯域であるほど、
上記所定の値が小さい値であることを特徴とするもので
ある。

【０１０８】言い換えれば、これらの関数ｆｋはそれぞ
れ、原点を通り、関数ｆｋの傾きがその関数ｆｋにより
変換される帯域制限画像信号Ｂｋの信号値にかかわらず
１以下であり、その関数ｆｋにより変換される帯域制限
画像信号Ｂｋの信号値の絶対値の増加にともない、その
関数ｆｋの傾きが０になるかまたは０に収束し、さら
に、その傾きが０または０に近い所定の値となるときの
帯域制限画像信号Ｂｋの信号値の絶対値が、高周波帯域
を変換する関数ｆｋであるほど小さい値であることを特
徴とする。

【０１０９】図２９に示す関数ｆｋは、振幅の大きな帯
域制限画像信号Ｂｋを抑制するような変換を行うもので
あり、周波数帯域の高い帯域制限画像信号Ｂｋの抑制の
度合いを、周波数帯域の低い帯域制限画像信号Ｂｋより
も強くするものであるが、これは実際の放射線画像のエ
ッジに含まれている高周波成分が、低周波成分に比べて
その振幅が小さいということを考慮したものである。す
なわち図３０に示すように、実際の放射線画像において
は、かなり急峻なエッジでさえも、図３０（ａ）のよう
に正確な階段状にはなっておらず、同図（ｂ）のように
高周波成分になるほどその振幅が小さくなっていること
が多い。このため、各周波数成分の振幅に合わせて、周
波数の高い帯域制限画像信号Ｂｋほど小さい振幅から抑
制を行うことが望ましく、図２９に示す関数ｆｋにより
それを実現することができるものである。

【０１１０】次に図３１に示す関数について説明する。
図３１に示す関数ｆｋは、帯域制限画像信号Ｂｋを、帯
域制限画像信号Ｂｋの信号値の絶対値に基づいて定ま
る、その絶対値以下の値となるように変換を行うもの
で、変換する帯域制限画像信号Ｂｋの周波数帯域が低周
波数帯域であるほど、信号値の絶対値が０近傍の所定の
範囲内の値である帯域制限画像信号Ｂｋを変換すること
により得られる変換画像信号Ｂｋ′の絶対値が小さい値
となることを特徴とするものである。

【０１１１】言い換えれば、これらの関数ｆｋはそれぞ
れ、原点を通り、関数ｆｋの傾きがその関数ｆｋにより
変換される帯域制限画像信号Ｂｋの信号値にかかわらず
１以下であり、その関数ｆｋの０近傍における傾きが、
低周波帯域の帯域制限画像信号Ｂｋを変換する関数ｆｋ
であるほど小さいことを特徴とするものである。

【０１１２】図３１に示す関数ｆｋは、変換画像信号Ｂ
ｋ′を積算して得られる積算信号（図２８（ｄ））を、
原画像信号Ｓorgに加えた場合、原画像信号Ｓorgと加算
された信号とのつなぎ目、すなわち信号の立ち上がりを
より自然なものとするという効果がある。

【０１１３】次に図３２に示す関数について説明する。
図３２に示す関数ｆｋは、図２９および図３１に示す関
数の双方の特徴を兼ね備えたものであり、双方の関数の
効果を得ることができるものである。

【０１１４】以上のように、帯域制限画像信号Ｂｋを変
換する関数ｆｋをその周波数帯域毎に目的に応じて異な
るものとすることにより、処理済み画像信号Ｓprocの周
波数応答特性を任意のものとすることが可能となる。図
３３はその効果の一例を示す図である。図３３におい
て、（ａ−１），（ｂ−１），（ｃ−１）は周波数帯域
に応じて異なる関数ｆｋにより帯域制限画像信号Ｂｋを
変換することにより得られた変換画像信号Ｂｋ′の周波
数帯域毎の周波数応答特性を示し、（ａ−２），（ｂ−
２），（ｃ−２）は（ａ−１），（ｂ−１），（ｃ−
１）にそれぞれ対応する処理済み画像信号Ｓprocの周波
数応答特性を示す。（ａ−１），（ａ−２）は全ての周
波数帯域において傾き１の関数を設定した場合、（ｂ−
１），（ｂ−２）は低周波帯域ほど傾きを小さく設定し
た場合、（ｃ−１），（ｃ−２）は特定の周波数帯域の
み関数の傾きを１とし、それ以外の周波数帯域の関数の
傾きを１より小さくした場合の周波数応答特性をそれぞ
れ示すものである。

【０１１５】図３３（ａ−１），（ａ−２）に示すよう
に、全ての周波数帯域において傾き１の関数ｆｋを使用
した場合、処理済み画像信号Ｓprocの周波数応答特性は
全ての周波数帯域に亘って略一定のレスポンスとなる。
これに対して、図３３（ｂ−１），（ｂ−２）に示すよ
うに、低周波数帯域ほど関数ｆｋの傾きを小さく設定し
た場合、処理済み画像信号Ｓprocの周波数応答特性は高
周波数帯域ほどレスポンスを高くすることができる。さ
らに、図３３（ｃ−１），（ｃ−２）に示すように、特
定の周波数帯域のみ関数ｆｋの傾きを１としそれ以外の
周波数帯域の関数ｆｋの傾きを１より小さくした場合、
処理済み画像信号Ｓprocにおける特定の周波数帯域のレ
スポンスを高くすることができる。

【０１１６】また、変換手段３により行う変換処理を強
調係数βに応じて異なる処理としてもよい。以下、その
効果を帯域制限画像信号Ｂｋの周波数帯域に応じて異な
る関数ｆｋを用いる場合を例に説明する。

【０１１７】このような関数としては、強調係数βが比
較的小さい場合には図３４に示すような関数ｆｋが、ま
た強調係数βが大きい場合には図３５に示すような関数
ｆｋが挙げられる。ここで、図中の「高」「低」はそれ
ぞれ、高周波帯域の帯域制限画像信号の変換に使用され
る関数、低周波帯域の帯域制限画像信号の変換に使用さ
れる関数を示している。

【０１１８】これらの関数ｆｋはいずれも帯域制限画像
信号Ｂｋの信号値をその信号値の絶対値よりも小さくな
るように抑制するものであるが、図３４および図３５を
比較するとわかるように、強調係数βが大きいほど、高
周波帯域の帯域制限画像信号Ｂｋを変換する関数ｆｋに
よる信号値の抑制の度合いと、低周波帯域の帯域制限画
像信号Ｂｋを変換する関数ｆｋによる信号値の抑制の度
合いとの差が大きくなるように関数ｆｋの組み合わせを
定義している。具体的には、高周波帯域の帯域制限画像
信号Ｂｋを変換する関数ｆｋは強調係数βに拘わらず一
定とし、低周波帯域の帯域制限画像信号Ｂｋを変換する
関数ｆｋの抑制の度合いを強調係数βが大きくなるほど
強くすることにより、抑制度合いの差を広げている。な
お、低周波帯域の帯域制限画像信号Ｂｋを変換する関数
ｆｋを一定として高周波帯域の帯域制限画像信号Ｂｋを
変換する関数ｆｋの抑制の度合いを緩めてもよく、また
全ての関数ｆｋを抑制の度合いの差が広がるように変化
させてもよい。

【０１１９】図３６は図３４に示す関数ｆｋを用いて処
理を行うことにより得られた処理済み画像信号Ｓprocの
周波数応答特性を示す図、図３７は図３５に示す関数ｆ
ｋを用いて処理を行うことにより得られた処理済み画像
信号Ｓprocの周波数応答特性を示す図である。図３６お
よび図３７に示すように、強調係数βを大きくした場合
には、高周波成分の強調の度合いが強くなり、低周波成
分の強調の度合いは強調係数βが小さい場合とそれほど
変わらないようになっていることがわかる。

【０１２０】強調係数βによらずに関数ｆｋの組み合わ
せを決定する場合には、強調係数βを大きくすると全て
の周波数帯域において強調の度合いが強くなる。一般に
強調係数βを大きくするのは、情報量が少ない高周波成
分を強調するためであるが、これでは、低周波成分まで
もが強調されてしまうことになり、その結果処理済み画
像においてコントラストがつきすぎ、アーチファクトが
発生する場合がある。したがって、図３５に示すよう
に、強調係数βを大きくしても低周波帯域の帯域制限画
像信号Ｂｋの強調の度合いをあまり変化させないように
することにより、アーチファクトの発生を防止すること
ができる。

【０１２１】なお、上記図３４および図３５に示す関数
ｆｋの説明においては、強調係数βが小さい場合および
大きい場合用に１組ずつ関数ｆｋの組み合わせを示した
が、これは関数ｆｋの組み合わせが２組であることを意
味するものではなく、強調係数βの変化に応じて関数ｆ
ｋの特徴をどのように変えればよいかを示すための例に
すぎない。したがって、強調係数βの変化の度合いを何
段階かのレベルに分類し、各レベルに応じて関数ｆｋの
組み合わせを決定するといったことも可能であり、その
ような段階をより細かく設定することにより、より高画
質な画像が得られることはいうまでもない。

【０１２２】また、上記図３４および図３５に示す関数
ｆｋの説明においては、強調係数βを大きくするのは高
周波成分のみを強調したいためであるとの解釈に基づい
て関数ｆｋを決定しているが、強調係数βを変化させる
目的はこれに限られるものではなく、関数ｆｋとしては
目的に応じて様々なものを適用可能である。

【０１２３】さらに、図３４および図３５に示す関数ｆ
ｋの説明においては、関数ｆｋは各周波数帯域毎に異な
る複数の関数ｆｋの組としているが、達成したい目的に
よっては、全ての周波数帯域に同一の関数ｆｋを用いる
ことも可能である。

【０１２４】次に、本発明の第２の実施形態について説
明する。なお、第２の実施形態におけるボケ画像信号作
成手段１、帯域制限画像信号作成手段２、積算手段４、
周波数強調処理手段５の処理は上記第１の実施形態と同
一であるため、ここでは説明を省略し、変換手段３の処
理についてのみ説明する。

【０１２５】図３８は第２の実施形態による画像処理装
置の構成を詳細に示すブロック図である。上記第１の実
施形態と同様、原画像信号Ｓorgおよびボケ画像信号作
成手段１において作成された複数のボケ画像信号Ｓusｋ
に基づいて帯域制限画像信号Ｂｋが作成される。この帯
域制限画像信号Ｂｋは減算器２１により、互いに隣接す
る周波数帯域の２つのボケ画像信号Ｓusｋ（但し原画像
信号ＳorgについてはＳorgとＳus１)の減算を行うこと
により得られる。すなわち、Ｓorg−Ｓus１、Ｓus１−
Ｓus２、…ＳusN-1−ＳusNを順次計算することにより、
複数の帯域制限画像信号Ｂｋ（ｋ＝１〜Ｎ）が求められ
る。第２の実施形態では、例えば帯域制限画像信号Ｓus
１−Ｓus２に対し、これに対応する補助画像信号をＳus
２−Ｓus３としている。したがって、帯域制限画像信号
Ｂｋを作成する手段と補助画像信号を作成する手段とは
実質的に同一の手段となる。すなわち、作成された帯域
制限画像信号Ｂｋは、帯域制限画像信号Ｂｋとして処理
されると同時に、隣接する帯域制限画像信号Ｂｋ−１に
対応する補助画像信号としても処理されるものとする。

【０１２６】上記のようにして求められた帯域制限画像
信号Ｂｋは、変換手段３において変換される。この変換
処理は、各帯域制限画像信号Ｂｋに対して、変換器２２
および変換器２４において行われる。すなわち、変換器
２２においては、各帯域制限画像信号Ｂｋに対して周波
数帯域毎に異なる関数ｆｋにより変換処理が施されて抑
制画像信号Ｒｋ（ｋ＝１〜Ｎ）が得られ、変換器２４に
おいては、関数ｇにより変換処理が施されて倍率信号Ｍ
ｋ（ｋ＝１〜Ｎ）が得られる。そして、乗算器２５にお
いて抑制画像信号Ｒｋと倍率信号Ｍｋとを乗算すること
により変換画像信号Ｂｋ′が得られる。この際、乗算
は、図３８に示すように、例えば帯域制限画像信号Ｓus
１−Ｓus２を変換して得た抑制画像信号Ｒ１と、補助画
像信号Ｓus２−Ｓus３を変換して得た倍率信号Ｍ１とい
うように、隣接した周波数帯域の信号同士で行われる。
ここで、関数ｆｋにより変換を行う変換器２２が抑制画
像信号作成手段に相当し、関数ｇにより変換を行う変換
器２４が倍率信号作成手段に相当し、それらの変換器２
２，２４の出力信号を乗ずる乗算器２５が乗算手段に相
当する。

【０１２７】関数ｆｋは、周波数帯域毎に異なる関数と
したが、全て同じ関数であってもよく、画像処理の目的
に応じて任意に設定可能なものとする。第２の実施形態
では、この関数ｆｋとして、例えば上記図２９に示すよ
うに、帯域制限画像信号Ｂｋをその帯域制限画像信号Ｂ
ｋの信号値の絶対値よりも小さくなるように抑制する、
帯域制限画像信号Ｂｋの周波数帯域毎に異なる関数ｆｋ
が使用される。

【０１２８】関数ｇとしては、例えば図３９に示すよう
なものが用いられる。図３９に示す関数ｇは、この関数
ｇによって補助画像信号を変換する場合、補助画像信号
の信号値の絶対値が小さいときには、変換後の値として
１に近い値が得られ、絶対値が大きいほど０に近い値が
得られることを示している。なお、Ｋは変換後の値が０
となる値のうち最小の値を示している。

【０１２９】上記変換手段３により得られた変換画像信
号Ｂｋ′は、演算器２３に入力される。演算器２３は、
積算手段４および周波数強調処理手段５を内包するもの
である。演算器２３においては、複数の変換画像信号Ｂ
ｋ′が積算され、この積算により得られた積算信号に原
画像信号Ｓorgに応じて決定される強調係数βが乗じら
れ、さらにこの強調係数βが乗じられた積算信号が原画
像信号Ｓorgと加算されて処理済み画像信号Ｓprocが得
られる。この処理を、上述した他の処理とともに数式と
して表すと、下記の式（２０）に示すものとなる。

【０１３０】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg−Ｓus1)・ｇ（Ｓus1−Ｓus2） ＋ｆ2(Ｓus1−Ｓus2)・ｇ（Ｓus2−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓusk−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN−ＳusN+1）｝ （２０） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：ボケ画像信号 ｆk(k=1〜N）：各帯域制限画像信号を変換して抑制信号
を作成するための関数 ｇ：各補助画像信号を変換して倍率信号を作成するため
の関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） 図４０は、第２の実施形態において、画像の濃度が急激
に変化しているエッジ部の近傍の信号値を処理した際の
効果を説明するための図である。（１）はエッジ部近傍
における階段状に信号値が変化する原画像信号Ｓorgと
その原画像信号Ｓorgに基づいて作成されたボケ画像信
号Ｓusｋを、（２）は（１）に対応する帯域制限画像信
号Ｂｋを、（３）は（２）の帯域制限画像信号Ｂｋより
１つ低い周波数帯域の信号を含む補助画像信号を、
（４）は（３）の補助画像信号を関数ｇにより処理する
ことにより得られる倍率信号Ｍｋを示し、（５）は
（２）の信号を関数ｆｋで処理した抑制画像信号Ｒｋに
（４）の倍率信号Ｍｋを乗算することにより得られた変
換画像信号Ｂｋ′を示している。なお（３）に示されて
いる値Ｋは、図３９に示した値Ｋであり、（３）の補助
画像信号がＫとなるときに（４）の倍率信号Ｍｋが０と
なることが図示されている。

【０１３１】図４０に示すように、（２）の信号を単に
絶対値が小さくなるように変換した場合には、変換画像
信号のピークの形状は滑らかになるものの、立ち上がり
部分は急峻なままとなる。これに対し、（５）に示す変
換画像信号Ｂｋ′は立ち上がり部分が滑らかになってい
る。このように変換画像信号Ｂｋ′の立ち上がり部分を
滑らかにすることにより、各変換画像信号Ｂｋ′を積算
して得られる積算信号の周波数帯域の境界に階段状に生
じるアーチファクトを防止することができる。

【０１３２】なお、関数ｇにより処理される補助画像信
号は、上述したように帯域制限画像信号Ｂｋと同一のも
のでなくてもよく、原画像信号Ｓorgとボケ画像信号Ｓu
sとの差分値（すなわちＳorg−Ｓusｋ）を補助画像信号
として用いてもよい。このような補助画像信号を用いた
処理を下記の式（２１）に示す。

【０１３３】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝｛ｆ1(Ｓorg−Ｓus1)・ｇ（Ｓorg−Ｓus2） ＋ｆ2(Ｓus1−Ｓus2)・ｇ（Ｓorg−Ｓus3）＋… ＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓorg−Ｓusk+1）＋… ＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg−ＳusN+1）｝ （２１） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：ボケ画像信号 ｆk(k=1〜N）：各帯域制限画像信号を変換して抑制信号
を作成するための関数 ｇ：各補助画像信号を変換して倍率信号を作成するため
の関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） 図４１は、式（２１）を適用してエッジ部近傍の画像信
号を処理した際の効果を説明するための図であり、図４
０と同様に、（１）はエッジ部近傍における階段状に信
号値が変化する原画像信号Ｓorgとその原画像信号Ｓorg
に基づいて作成されたボケ画像信号Ｓusｋを、（２）は
（１）に対応する帯域制限画像信号Ｂｋを、（３）は
（２）の帯域制限画像信号Ｂｋに対応する補助画像信号
を、（４）は（３）の補助画像信号を関数ｇにより処理
することにより得られる倍率信号Ｍｋを示し、（５）は
（２）の信号を関数ｆｋで処理した抑制画像信号Ｒｋに
（４）の倍率信号Ｍｋを乗算することにより得られた変
換画像信号Ｂｋ′を示している。

【０１３４】図４１（５）に示すように、関数ｇにより
処理する値としてＳorg−Ｓusｋを用いる場合、コント
ラストの大きなエッジについては、変換画像信号Ｂｋ′
の信号値は小さくなり、コントラストの小さなエッジに
ついては変換画像信号Ｂｋ′の信号値はもとの帯域制限
画像信号Ｂｋに近い大きさとなる。

【０１３５】この変換画像信号Ｂｋ′は積算された後
に、原画像信号Ｓorgに加算される。この場合、処理済
み画像信号Ｓprocにおいて、コントラストの大きなエッ
ジはほとんど強調されないのに対し、コントラストの小
さなエッジはコントラストの大きなエッジに比べると相
対的に強く強調されるという効果が得られる。

【０１３６】以上、２種類の式（２０）、（２１）につ
いて説明したが、関数ｆｋおよびｇについて、あるいは
関数ｇにより処理する帯域制限画像信号Ｂｋの作成方法
については、様々な変更が可能である。

【０１３７】次に、本発明の第３の実施形態について説
明する。但し、第３の実施形態についても、ボケ画像信
号作成手段１、帯域制限画像信号作成手段２、積算手段
４、周波数強調処理手段５の処理は上記第１の実施形態
と同一であるため説明を省略し、変換手段３の処理につ
いてのみ説明する。

【０１３８】図４２は本発明の第３の実施形態による画
像処理装置の構成を詳細に示すブロック図である。上記
第１および第２の実施形態と同様、原画像信号Ｓorgお
よびボケ画像信号作成手段１において作成された複数の
ボケ画像信号Ｓusｋに基づいて帯域制限画像信号Ｂｋが
作成される。この帯域制限画像信号Ｂｋは減算器２１に
より、互いに隣接する周波数帯域の２つのボケ画像信号
Ｓusｋ（但し原画像信号ＳorgについてはＳorgとＳus
１)の減算を行うことにより得られる。すなわち、Ｓorg
−Ｓus１、Ｓus１−Ｓus２、…ＳusN-1−ＳusNを順次計
算することにより、複数の帯域制限画像信号Ｂｋが求め
られる。

【０１３９】上記のようにして求められた帯域制限画像
信号Ｂｋは、変換手段３において変換される。この変換
処理では、図４２に示すように、処理を施すべき周波数
帯域の帯域制限画像信号を被変換帯域制限画像信号Ｂｋ
とし、被変換帯域制限画像信号Ｂｋの周波数帯域よりも
１つ低い周波数帯域の帯域制限画像信号Ｂｋ＋１を変換
器２４において関数ｇを用いて変換することにより補助
画像信号Ａｋを得、補助画像信号Ａｋと被変換帯域制限
画像信号Ｂｋとを加算して複合帯域制限画像信号Ｃｋを
得、さらに複合帯域制限画像信号Ｃｋを変換器２２にお
いて関数ｆｋを用いて変換して変換画像信号Ｂｋ′を得
るものである。ここで、関数ｇにより変換を行う変換器
２４が、補助画像信号作成手段に相当し、加算器２６が
複合帯域制限画像信号作成手段に相当する。

【０１４０】第３の実施形態において使用される関数ｇ
としては、例えば図４３に示すものが用いられる。図４
３に示す関数ｇは原点を通り原点における傾きが略０で
あり、処理される信号値が大きくなるにつれて傾きが漸
増し、最終的に傾きが略１となるような関数である。す
なわち、信号値の大きさが小さいときには強く抑制し、
信号値が大きくなるにつれて抑制の度合いを緩めるとい
うものである。上述の変換処理において、この関数ｇ
の、原点から徐々に傾きが増加する部分は、補助画像信
号Ａｋの立ち上がりの部分の波形に影響を与える。すな
わち、この関数ｇにより変換処理を行うことにより、急
峻な立ち上がり部分を滑らかにすることができる。ここ
で、実際の帯域制限画像信号Ｂｋの大きさには限度があ
るため、関数ｇは原点付近で傾きが０から徐々に増加す
るような関数であれば、どのようなものでもよく例え
ば、図４４に示すものであってもよい。また、第３の実
施形態では、信号の増幅は行わないため、傾きは最大１
としているが、第３の実施形態の効果は関数の傾きが原
点付近で漸増することにより得られるものであるため、
この条件が満たされていれば、必ずしも傾きを最大１に
限定する必要はない。

【０１４１】関数ｆｋは、全て同じ関数であってもよい
し、異なる関数であってもよく、画像処理の目的に応じ
て任意に設定可能なものとする。第３の実施形態では、
この関数ｆｋとして、例えば図２９のような関数ｆｋを
使用する。

【０１４２】上記変換手段３により得られた変換画像信
号Ｂｋ′は、演算器２３に入力される。演算器２３は、
積算手段４および周波数強調処理手段５を内包するもの
である。演算器２３においては、複数の変換画像信号Ｂ
ｋ′が積算され、この積算により得られた積算信号に原
画像信号Ｓorgに応じて決定される強調係数βが乗じら
れ、さらにこの強調係数βが乗じられた積算信号が原画
像信号Ｓorgと加算されて処理済み画像信号Ｓprocが得
られる。この処理を、上述した他の処理とともに数式と
して表すと、下記の式（２２）に示すものとなる。

【０１４３】 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN） ＝［ｆ1｛（Ｓorg−Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝ ＋ｆ2｛（Ｓus1−Ｓus2)＋ｇ（Ｓus2−Ｓus3）｝＋… ＋ｆk｛（Ｓusk-1−Ｓusk)＋ｇ（Ｓusk−Ｓusk+1）｝＋… ＋ｆN｛（ＳusN-1−ＳusN)＋ｇ（ＳusN−ＳusN+1）｝］ （２２） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：ボケ画像信号 ｆk(k=1〜N）：各複合帯域制限画像信号を変換するため
の関数 ｇ：各帯域制限画像信号を変換して補助画像信号を作成
するための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） 図４５は、第３の実施形態において画像の濃度が急激に
変化しているエッジ部の近傍の信号値を処理した際の効
果を説明するための図である。（１）はエッジ部近傍に
おける階段状に信号値が変化する原画像信号Ｓorgとそ
の原画像信号Ｓorgに基づいて作成されたボケ画像信号
Ｓusｋを、（２）および（３）は（１）の信号にかかる
帯域制限画像信号Ｂｋであって、（２）を被変換帯域制
限画像信号とした場合の低周波側帯域制限画像信号が
（３）である。（４）は（３）の低周波側帯域制限画像
信号を関数ｇにより変換することにより得られる補助画
像信号Ａｋを示し、（５）は（２）の被変換帯域制限画
像信号Ｂｋと（４）の補助画像信号Ａｋを加算すること
により得られる複合帯域制限画像信号Ｃｋを示し、
（６）は（５）の複合帯域制限画像信号Ｃｋを関数ｆｋ
により変換することにより得られる変換画像信号Ｂｋ′
を示している。図４５において、（２）の被変換帯域制
限画像信号Ｂｋを単に絶対値が小さくなるように変換し
た場合には、変換画像信号のピークの形状は滑らかにな
るものの、立ち上がり部分は急峻なままとなる。これに
対し、（５）に示される変換画像信号Ｂｋ′は立ち上が
り部分が滑らかになっている。このように変換画像信号
Ｂｋ′の立ち上がり部分を滑らかにすることにより、変
換画像信号Ｂｋ′を積算して得られる積算信号の周波数
帯域の境界に階段状に生じるアーチファクトを防止する
ことができる。

【０１４４】以上、式（２２）の処理について説明した
が、第３の実施形態においても、例えば、関数ｆｋおよ
びｇについて、様々な変更が可能である。

【０１４５】なお、上記第１の実施形態において、帯域
制限画像信号Ｂｋに対して撮影部位毎に異なる変換処理
を施すことにより得られる効果について説明したが、こ
の撮影部位毎に処理を変えるという思想は上記全ての実
施形態に適用できるものである。

【０１４６】また、上記第１から第３の実施形態におい
ては、正方状画素配列の原画像信号Ｓorgからボケ画像
信号Ｓusｋおよび帯域制限画像信号Ｂｋを作成して処理
済み画像信号Ｓprocを得ているが、市松状画素配列の原
画像信号Ｓorgに対して同様の処理を行うことも可能で
ある。この場合、ボケ画像信号作成手段１においては、
図４６に示すように、原画像信号Ｓorgに対して市松正
方縮小処理である補間処理ｃが施されて、正方状画素配
列の縮小画像を表す縮小画像信号Ｓ１が得られ、さらに
縮小画像信号Ｓ１に対して正方市松縮小処理である補間
処理ａが施されて、市松状画素配列の縮小画像を表す縮
小画像信号Ｓ２が得られる。

【０１４７】そして補間処理ｃ，ａを繰り返すことによ
り、Ｎ個の縮小画像信号Ｓｋ（ｋ＝１〜Ｎ）が得られ
る。ここで、ｋが奇数の場合、縮小画像信号Ｓｋにより
表される縮小画像は正方状の画素配列となる。一方、ｋ
が偶数の場合、縮小画像信号Ｓｋにより表される縮小画
像は市松状の画素配列となる。また、最低周波数帯域の
縮小画像信号ＢＮを補間処理ａまで行うことにより得た
場合、その縮小画像信号ＳＮにより表される縮小画像は
正方状の画素配列となる。一方、最低周波数帯域の縮小
画像信号ＳＮを補間処理ｃまで行うことにより得た場
合、その縮小画像信号ＳＮにより表される縮小画像は市
松状の画素配列となる。

【０１４８】そして、このようにして得られた縮小画像
信号Ｓｋに対して、補間処理ｄ（ｋ：奇数）、補間処理
ｂ（ｋ：偶数）を施して、原画像と同一サイズのボケ画
像を表すボケ画像信号Ｓusｋを得る。以下、上記各実施
形態と同一の処理を行うことにより、処理済み画像信号
Ｓprocを得ることができる。

【０１４９】なお、上記各実施形態においては、全ての
周波数帯域の帯域制限画像信号Ｂｋに対して関数ｆｋに
よる変換処理を施しているが、任意の１以上の周波数帯
域の帯域制限画像信号Ｂｋにのみ変換処理を施してもよ
い。この場合、変換処理により得られた変換画像信号Ｂ
ｋ′のみが積算され、さらに原画像信号Ｓorgに加算さ
れて、処理済み画像信号Ｓprocが得られることとなる。

【図面の簡単な説明】

【図１】本発明の実施形態による画像処理装置の構成を
示す概略ブロック図

【図２】ボケ画像信号作成手段において行われる処理を
示す概略ブロック図

【図３】原画像信号の画素配列を示す図

【図４】補間処理ａにおいて用いられるフィルタを示す
図

【図５】補間処理ａにより得られる縮小画像信号の画素
配列を示す図

【図６】補間処理ｃにおいて用いられるフィルタを示す
図

【図７】補間処理ｃにより得られる縮小画像信号の画素
配列を示す図

【図８】補間処理ｃにより得られる縮小画像信号にさら
に補間処理ａを施すことにより得られる縮小画像信号の
画素配列を示す図

【図９】縮小画像信号の周波数応答特性を表す図

【図１０】補間処理ｄにおいて用いられるフィルタを示
す図（その１）

【図１１】補間処理ｄにおいて用いられるフィルタを示
す図（その２）

【図１２】補間処理ｂにおいて用いられるフィルタを示
す図（その１）

【図１３】補間処理ｂにおいて用いられるフィルタを示
す図（その２）

【図１４】ボケ画像信号の周波数応答特性を表す図

【図１５】本発明の第１の実施形態による画像処理装置
の構成を詳細に示すブロック図

【図１６】帯域制限画像信号の周波数応答特性を表す図

【図１７】変換手段において帯域制限画像信号を変換す
る関数の一例を示す図

【図１８】処理済み画像信号の周波数応答特性の例を示
す図

【図１９】従来の方法により得られた帯域制限画像信号
と本実施形態により得られた帯域制限画像信号との周波
数応答特性の相違を説明するための図

【図２０】処理済み画像信号の周波数応答特性の他の例
を示す図

【図２１】変換手段において帯域制限画像信号を変換す
る関数の他の例を示す図

【図２２】帯域制限画像信号の周波数応答特性の他の例
を示す図

【図２３】撮影部位に応じて帯域制限画像信号を変換す
る関数の例を示す図

【図２４】周波数帯域に応じて帯域制限画像信号を変換
する関数の例を示す図（その１）

【図２５】周波数帯域に応じて帯域制限画像信号を変換
する関数の例を示す図（その２）

【図２６】図２４に示す関数により得られる処理済み画
像信号の周波数応答特性を示す図

【図２７】図２５に示す関数により得られる処理済み画
像信号の周波数応答特性を示す図

【図２８】全ての帯域制限画像信号を１種類の関数で変
換した際の問題点を示す図

【図２９】周波数帯域に応じて帯域制限画像信号を変換
する関数の例を示す図（その３）

【図３０】理想のエッジと実際のエッジの信号を示す図

【図３１】周波数帯域に応じて帯域制限画像信号を変換
する関数の例を示す図（その４）

【図３２】図２９および図３１の関数の特徴を兼ね備え
た関数の一例を示す図

【図３３】周波数帯域毎の周波数応答特性と全体の周波
数応答特性とを示す図

【図３４】強調係数が小さい場合に各帯域制限画像信号
を変換する関数の一例を示す図

【図３５】強調係数が大きい場合に各帯域制限画像信号
を変換する関数の一例を示す図

【図３６】図３４に示す関数を変換に用いることにより
得られる処理済み画像信号の周波数応答特性を示す図

【図３７】図３５に示す関数を変換に用いることにより
得られる処理済み画像信号の周波数応答特性を示す図

【図３８】本発明の第２の実施形態による画像処理装置
の構成を詳細に示すブロック図

【図３９】倍率信号作成手段において使用される関数の
一例を示す図

【図４０】第２の実施形態の効果を説明するための図
（その１）

【図４１】第２の実施形態の効果を説明するための図
（その２）

【図４２】本発明の第３の実施形態による画像処理装置
の構成を詳細に示すブロック図

【図４３】各帯域制限画像信号を変換して補助画像信号
を作成するための関数の例を示す図（その１）

【図４４】各帯域制限画像信号を変換して補助画像信号
を作成するための関数の例を示す図（その２）

【図４５】第３の実施形態の効果を説明するための図

【図４６】市松状画素配列の原画像信号からのボケ画像
信号の作成を説明するための図

【図４７】正方状の画素配列を示す図

【図４８】市松状の画素配列を示す図

【図４９】エッジ部近傍のボケ画像信号の作成処理を示
す図

【符号の説明】

１ ボケ画像信号作成手段 ２ 帯域制限画像信号作成手段 ３ 変換手段 ４ 積算手段 ５ 周波数強調処理手段 10 フィルタリング処理手段 11 補間演算処理手段 21 減算器 22 変換器 23 演算器 24 変換器 25 乗算手段 26 加算器

Claims (53)

【特許請求の範囲】
1. 【請求項１】 原画像を表す原画像信号に、該原画像
   の高周波成分に関する信号を加算することにより、該原
   画像の高周波成分を強調する画像処理方法において、 前記原画像信号に対して正方市松縮小処理および市松正
   方縮小処理に基づく処理を施して、互いに周波数応答特
   性が異なる複数の非鮮鋭マスク画像信号を作成し、 前記原画像信号および前記複数の非鮮鋭マスク画像信
   号、または前記複数の非鮮鋭マスク画像信号に基づい
   て、前記原画像の複数の周波数帯域毎の画像を表す複数
   の帯域制限画像信号を作成し、 該複数の帯域制限画像信号のうち少なくとも１つの帯域
   制限画像信号に対して、所定の変換処理を施して少なく
   とも１つの変換画像信号を作成し、 該変換画像信号に基づいて前記原画像信号に加算する前
   記高周波成分に関する信号を得ることを特徴とする画像
   処理方法。
2. 【請求項２】 前記所定の変換処理は、前記帯域制限
   画像信号の信号値の少なくとも一部を小さくする処理で
   あることを特徴とする請求項１記載の画像処理方法。
3. 【請求項３】 前記所定の変換処理は、前記帯域制限
   画像信号の信号値の絶対値が所定の閾値より大きい場合
   に、該信号値の絶対値を小さくする処理であることを特
   徴とする請求項１または２記載の画像処理方法。
4. 【請求項４】 前記所定の変換処理は、さらに前記帯
   域制限画像信号の信号値の絶対値が前記所定の閾値より
   も小さい他の閾値よりも小さいほど、該信号値の絶対値
   を小さくする処理であることを特徴とする請求項３記載
   の画像処理方法。
5. 【請求項５】 前記所定の変換処理は、前記帯域制限
   画像信号の周波数帯域に応じて、該帯域制限画像信号の
   信号値の絶対値の大きさを変化させる処理であることを
   特徴とする請求項３または４記載の画像処理方法。
6. 【請求項６】 前記所定の変換処理は、前記帯域制限
   画像信号を、該帯域制限画像信号の周波数帯域に応じて
   異なる関数に基づいて、該帯域制限画像信号の信号値の
   絶対値に基づいて定まる該絶対値以下の値となるように
   変換する処理であることを特徴とする請求項１または２
   記載の画像処理方法。
7. 【請求項７】 前記関数は、前記帯域制限画像信号の
   信号値の絶対値が所定の値より大きい場合は、前記変換
   画像信号の信号値が略一定の値となるように前記帯域制
   限画像信号を変換する関数であり、前記所定の値は、前
   記帯域制限画像信号が高周波数帯域であるほど小さい値
   であることを特徴とする請求項６記載の画像処理方法。
8. 【請求項８】 前記関数は、前記帯域制限画像信号が
   低周波帯域であるほど、該帯域制限画像信号における絶
   対値が０近傍の所定の範囲内の信号値を変換することに
   より得られる前記変換画像信号の絶対値がより小さい値
   となるように、前記帯域制限画像信号を変換する関数で
   あることを特徴とする請求項６または７記載の画像処理
   方法。
9. 【請求項９】 前記帯域制限画像信号の作成、前記変
   換画像信号の作成、前記高周波成分に関する信号の作
   成、および該高周波成分に関する信号の前記原画像信号
   への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
   s2,…Ｓ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. 【請求項１０】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ1
    （Ｓorg−Ｓus1)＋ｆ2（Ｓorg−Ｓus2)＋…＋ｆk(Ｓorg
    −Ｓusk)＋…＋ｆN（Ｓorg−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１から８のい
    ずれか１項記載の画像処理方法。
11. 【請求項１１】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝｛ｆ2（Ｓus1−Ｓus2)
    ＋ｆ3（Ｓus2−Ｓus3)＋…＋ｆk(Ｓusk-1−Ｓusk)＋…
    ＋ｆN（ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１から８のい
    ずれか１項記載の画像処理方法。
12. 【請求項１２】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ2（Ｓus1
    −Ｓus2)＋ｆ3（Ｓus1−Ｓus3)＋…＋ｆk(Ｓus1−Ｓus
    k)＋…＋ｆN（Ｓus1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１から８のい
    ずれか１項記載の画像処理方法。
13. 【請求項１３】 前記所定の変換処理は、前記強調係
    数に応じた変換処理であることを特徴とする請求項９か
    ら１２のいずれか１項記載の画像処理方法。
14. 【請求項１４】 前記強調係数に応じた変換処理が、
    前記帯域制限画像信号の周波数帯域に応じて異なる変換
    処理であって、該帯域制限画像信号の周波数帯域が低周
    波帯域であるほど、該帯域制限画像信号を大きく抑制す
    る処理であり、 前記強調係数が大きいほど、高周波帯域の帯域制限画像
    信号に対する抑制の度合いと、低周波帯域の帯域制限画
    像信号に対する抑制の度合いとの差が大きいことを特徴
    とする請求項１３記載の画像処理方法。
15. 【請求項１５】 前記所定の変換処理は、前記帯域制
    限画像信号を、該帯域制限画像信号の信号値の絶対値に
    基づいて定まる該絶対値以下の値となるように変換する
    ことにより抑制画像信号を作成し、 前記原画像信号および前記複数の非鮮鋭マスク画像信
    号、または前記複数の非鮮鋭マスク画像信号に基づい
    て、前記抑制画像信号の作成に使用された帯域制限画像
    信号よりも低い周波数帯域の信号を含む補助画像信号を
    作成し、 該補助画像信号を、該補助画像信号の信号値の絶対値が
    小さいほど１に近く、大きいほど０に近い値となるよう
    に変換することにより前記抑制画像信号にそれぞれ対応
    する倍率信号を作成し、 前記抑制画像信号に、該抑制画像信号に対応する前記倍
    率信号を乗ずることにより前記変換画像信号を作成する
    処理であることを特徴とする請求項１または２記載の画
    像処理方法。
16. 【請求項１６】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓus1−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓus2−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓu
    sk−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１５記載の画
    像処理方法。
17. 【請求項１７】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓorg−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓorg−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓo
    rg−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１５記載の画
    像処理方法。
18. 【請求項１８】 前記所定の変換処理は、前記少なく
    とも１つの帯域制限画像信号である被変換帯域制限画像
    信号よりも低い周波数帯域の帯域制限画像信号である低
    周波側帯域制限画像信号を、原点を通り該原点における
    傾きが略０で、処理される値が大きくなるにつれて該傾
    きが漸増する非線形関数に基づいて変換することにより
    前記被変換帯域制限画像信号の補助画像信号を作成し、 該補助画像信号を前記被変換帯域制限画像信号に加算す
    ることにより複合帯域制限画像信号を作成し、 該複合帯域制限画像信号を、該複合帯域制限画像信号の
    信号値の絶対値に基づいて定まる該絶対値以下の値とな
    るように変換することにより前記変換画像信号を作成す
    る処理であることを特徴とする請求項１または２記載の
    画像処理方法。
19. 【請求項１９】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝［ｆ1｛（Ｓorg
    −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝＋ｆ2｛（Ｓus1−Ｓus
    2)＋ｇ（Ｓus2−Ｓus3）｝＋…＋ｆk｛（Ｓusk-1−Ｓus
    k)＋ｇ（Ｓusk−Ｓusk+1）｝＋…＋ｆN｛（ＳusN-1−Ｓ
    usN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：複合帯域制限画像信号を変換するための
    関数 ｇ：帯域制限画像信号を変換して補助画像信号を作成す
    るための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項１８記載の画
    像処理方法。
20. 【請求項２０】 前記所定の変換処理は、前記原画像
    を撮影した際の撮影部位に応じて、前記帯域制限画像信
    号を変換する処理であることを特徴とする請求項１から
    １９のいずれか１項記載の画像処理方法。
21. 【請求項２１】 原画像を表す原画像信号に、該原画
    像の高周波成分に関する信号を加算することにより、該
    原画像の高周波成分を強調する画像処理装置において、 前記原画像信号に対して正方市松縮小処理および市松正
    方縮小処理に基づく処理を施して、互いに周波数応答特
    性が異なる複数の非鮮鋭マスク画像信号を作成する非鮮
    鋭マスク画像信号作成手段と、 前記原画像信号および前記複数の非鮮鋭マスク画像信
    号、または前記複数の非鮮鋭マスク画像信号に基づい
    て、前記原画像の複数の周波数帯域毎の画像を表す複数
    の帯域制限画像信号を作成する帯域制限画像信号作成手
    段と、 該複数の帯域制限画像信号のうち少なくとも１つの帯域
    制限画像信号に対して、所定の変換処理を施して少なく
    とも１つの変換画像信号を作成する変換手段と、 該変換画像信号に基づいて前記原画像信号に加算する前
    記高周波成分に関する信号を得る周波数強調処理手段と
    からなることを特徴とする画像処理装置。
22. 【請求項２２】 前記所定の変換処理は、前記帯域制
    限画像信号の信号値の少なくとも一部を小さくする処理
    であることを特徴とする請求項２１記載の画像処理装
    置。
23. 【請求項２３】 前記所定の変換処理は、前記帯域制
    限画像信号の信号値の絶対値が所定の閾値より大きい場
    合に、該信号値の絶対値を小さくする処理であることを
    特徴とする請求項２１または２２記載の画像処理装置。
24. 【請求項２４】 前記所定の変換処理は、さらに前記
    帯域制限画像信号の信号値の絶対値が前記所定の閾値よ
    りも小さい他の閾値よりも小さいほど、該信号値の絶対
    値を小さくする処理であることを特徴とする請求項２３
    記載の画像処理装置。
25. 【請求項２５】 前記所定の変換処理は、前記帯域制
    限画像信号の周波数帯域に応じて、該帯域制限画像信号
    の信号値の絶対値の大きさを変化させる処理であること
    を特徴とする請求項２３または２４記載の画像処理装
    置。
26. 【請求項２６】 前記所定の変換処理は、前記帯域制
    限画像信号を、該帯域制限画像信号の周波数帯域に応じ
    て異なる関数に基づいて、該帯域制限画像信号の信号値
    の絶対値に基づいて定まる該絶対値以下の値となるよう
    に変換する処理であることを特徴とする請求項２１また
    は２２記載の画像処理装置。
27. 【請求項２７】 前記関数は、前記帯域制限画像信号
    の信号値の絶対値が所定の値より大きい場合は、前記変
    換画像信号の信号値が略一定の値となるように前記帯域
    制限画像信号を変換する関数であり、前記所定の値は、
    前記帯域制限画像信号が高周波数帯域であるほど小さい
    値であることを特徴とする請求項２６記載の画像処理装
    置。
28. 【請求項２８】 前記関数は、前記帯域制限画像信号
    が低周波帯域であるほど、該帯域制限画像信号における
    絶対値が０近傍の所定の範囲内の信号値を変換すること
    により得られる前記変換画像信号の絶対値がより小さい
    値となるように、前記帯域制限画像信号を変換する関数
    であることを特徴とする請求項２６または２７記載の画
    像処理装置。
29. 【請求項２９】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…Ｓ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）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項２１から２８
    のいずれか１項記載の画像処理装置。
30. 【請求項３０】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ1
    （Ｓorg−Ｓus1)＋ｆ2（Ｓorg−Ｓus2)＋…＋ｆk(Ｓorg
    −Ｓusk)＋…＋ｆN（Ｓorg−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項２１から２８
    のいずれか１項記載の画像処理装置。
31. 【請求項３１】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝｛ｆ2（Ｓus1−Ｓus2)
    ＋ｆ3（Ｓus2−Ｓus3)＋…＋ｆk(Ｓusk-1−Ｓusk)＋…
    ＋ｆN（ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項２１から２８
    のいずれか１項記載の画像処理装置。
32. 【請求項３２】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ2（Ｓus1
    −Ｓus2)＋ｆ3（Ｓus1−Ｓus3)＋…＋ｆk(Ｓus1−Ｓus
    k)＋…＋ｆN（Ｓus1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項２１から２８
    のいずれか１項記載の画像処理装置。
33. 【請求項３３】 前記所定の変換処理は、前記強調係
    数に応じた変換処理であることを特徴とする請求項２９
    から３２のいずれか１項記載の画像処理装置。
34. 【請求項３４】 前記強調係数に応じた変換処理が、
    前記帯域制限画像信号の周波数帯域に応じて異なる変換
    処理であって、該帯域制限画像信号の周波数帯域が低周
    波帯域であるほど、該帯域制限画像信号を大きく抑制す
    る処理であり、 前記強調係数が大きいほど、高周波帯域の帯域制限画像
    信号に対する抑制の度合いと、低周波帯域の帯域制限画
    像信号に対する抑制の度合いとの差が大きいことを特徴
    とする請求項３３記載の画像処理装置。
35. 【請求項３５】 前記変換手段が、前記帯域制限画像
    信号を、該帯域制限画像信号の信号値の絶対値に基づい
    て定まる該絶対値以下の値となるように変換することに
    より抑制画像信号を作成する抑制画像信号作成手段と、 前記原画像信号および前記複数の非鮮鋭マスク画像信
    号、または前記複数の非鮮鋭マスク画像信号に基づい
    て、前記抑制画像信号の作成に使用された帯域制限画像
    信号よりも低い周波数帯域の信号を含む補助画像信号を
    作成する補助画像信号作成手段と、 該補助画像信号を、該補助画像信号の信号値の絶対値が
    小さいほど１に近く、大きいほど０に近い値となるよう
    に変換することにより前記抑制画像信号にそれぞれ対応
    する倍率信号を作成する倍率信号作成手段と、 前記抑制画像信号に、該抑制画像信号に対応する前記倍
    率信号を乗ずることにより前記変換画像信号を作成する
    乗算手段とを備えてなることを特徴とする請求項２１ま
    たは２２記載の画像処理装置。
36. 【請求項３６】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓus1−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓus2−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓu
    sk−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項３５記載の画
    像処理装置。
37. 【請求項３７】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓorg−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓorg−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓo
    rg−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項３５記載の画
    像処理装置。
38. 【請求項３８】 前記変換手段が、前記少なくとも１
    つの帯域制限画像信号である被変換帯域制限画像信号よ
    りも低い周波数帯域の帯域制限画像信号である低周波側
    帯域制限画像信号を、原点を通り該原点における傾きが
    略０で、処理される値が大きくなるにつれて該傾きが漸
    増する非線形関数に基づいて変換することにより前記被
    変換帯域制限画像信号の補助画像信号を作成する補助画
    像信号作成手段と、 該補助画像信号を前記被変換帯域制限画像信号に加算す
    ることにより複合帯域制限画像信号を作成する複合帯域
    制限画像信号作成手段と、 該複合帯域制限画像信号を、該複合帯域制限画像信号の
    信号値の絶対値に基づいて定まる該絶対値以下の値とな
    るように変換することにより前記変換画像信号を作成す
    る変換画像信号作成手段とを備えたことを特徴とする請
    求項２１または２２記載の画像処理装置。
39. 【請求項３９】 前記帯域制限画像信号の作成、前記
    変換画像信号の作成、前記高周波成分に関する信号の作
    成、および該高周波成分に関する信号の前記原画像信号
    への加算を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝［ｆ1｛（Ｓorg
    −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝＋ｆ2｛（Ｓus1−Ｓus
    2)＋ｇ（Ｓus2−Ｓus3）｝＋…＋ｆk｛（Ｓusk-1−Ｓus
    k)＋ｇ（Ｓusk−Ｓusk+1）｝＋…＋ｆN｛（ＳusN-1−Ｓ
    usN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：複合帯域制限画像信号を変換するための
    関数 ｇ：帯域制限画像信号を変換して補助画像信号を作成す
    るための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項３８記載の画
    像処理装置。
40. 【請求項４０】 前記所定の変換処理は、前記原画像
    を撮影した際の撮影部位に応じて、前記帯域制限画像信
    号を変換する処理であることを特徴とする請求項２１か
    ら３９のいずれか１項記載の画像処理装置。
41. 【請求項４１】 原画像を表す原画像信号に、該原画
    像の高周波成分に関する信号を加算することにより、該
    原画像の高周波成分を強調する画像処理方法をコンピュ
    ータに実行させるためのプログラムを記録したコンピュ
    ータ読取り可能な記録媒体において、 前記プログラムは、前記原画像信号に対して正方市松縮
    小処理および市松正方縮小処理に基づく処理を施して、
    互いに周波数応答特性が異なる複数の非鮮鋭マスク画像
    信号を作成する手順と、 前記原画像信号および前記複数の非鮮鋭マスク画像信
    号、または前記複数の非鮮鋭マスク画像信号に基づい
    て、前記原画像の複数の周波数帯域毎の画像を表す複数
    の帯域制限画像信号を作成する手順と、 該複数の帯域制限画像信号のうち少なくとも１つの帯域
    制限画像信号に対して、所定の変換処理を施して少なく
    とも１つの変換画像信号を作成する手順と、 該変換画像信号に基づいて前記原画像信号に加算する前
    記高周波成分に関する信号を得る手順とを有することを
    特徴とするコンピュータ読取り可能な記録媒体。
42. 【請求項４２】 前記所定の変換処理は、前記帯域制
    限画像信号の信号値の少なくとも一部を小さくする処理
    であることを特徴とする請求項４１記載のコンピュータ
    読取り可能な記録媒体。
43. 【請求項４３】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…Ｓ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）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４１または４
    ２記載のコンピュータ読取り可能な記録媒体。
44. 【請求項４４】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ1
    （Ｓorg−Ｓus1)＋ｆ2（Ｓorg−Ｓus2)＋…＋ｆk(Ｓorg
    −Ｓusk)＋…＋ｆN（Ｓorg−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４１または４
    ２記載のコンピュータ読取り可能な記録媒体。
45. 【請求項４５】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝｛ｆ2（Ｓus1−Ｓus2)
    ＋ｆ3（Ｓus2−Ｓus3)＋…＋ｆk(Ｓusk-1−Ｓusk)＋…
    ＋ｆN（ＳusN-1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４１または４
    ２記載のコンピュータ読取り可能な記録媒体。
46. 【請求項４６】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓus1,Ｓus2,…
    ＳusN） Ｆusm（Ｓus1,Ｓus2,…ＳusN）＝（1/N)・｛ｆ2（Ｓus1
    −Ｓus2)＋ｆ3（Ｓus1−Ｓus3)＋…＋ｆk(Ｓus1−Ｓus
    k)＋…＋ｆN（Ｓus1−ＳusN)｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：所定の変換処理を行う関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４１または４
    ２記載のコンピュータ読取り可能な記録媒体。
47. 【請求項４７】 前記所定の変換処理は、前記強調係
    数に応じた変換処理を行う手順を有することを特徴とす
    る請求項４３から４６のいずれか１項記載のコンピュー
    タ読取り可能な記録媒体。
48. 【請求項４８】 前記所定の変換処理は、前記帯域制
    限画像信号を、該帯域制限画像信号の信号値の絶対値に
    基づいて定まる該絶対値以下の値となるように変換する
    ことにより抑制画像信号を作成する手順と、 前記原画像信号および前記複数の非鮮鋭マスク画像信
    号、または前記複数の非鮮鋭マスク画像信号に基づい
    て、前記抑制画像信号の作成に使用された帯域制限画像
    信号よりも低い周波数帯域の信号を含む補助画像信号を
    作成する手順と、 該補助画像信号を、該補助画像信号の信号値の絶対値が
    小さいほど１に近く、大きいほど０に近い値となるよう
    に変換することにより前記抑制画像信号にそれぞれ対応
    する倍率信号を作成する手順と、 前記抑制画像信号に、該抑制画像信号に対応する前記倍
    率信号を乗ずることにより前記変換画像信号を作成する
    手順とを有することを特徴とする請求項４１または４２
    記載のコンピュータ読取り可能な記録媒体。
49. 【請求項４９】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓus1−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓus2−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓu
    sk−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（ＳusN
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４８記載のコ
    ンピュータ読取り可能な記録媒体。
50. 【請求項５０】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝｛ｆ1(Ｓorg−
    Ｓus1)・ｇ（Ｓorg−Ｓus2）＋ｆ2(Ｓus1−Ｓus2)・ｇ
    （Ｓorg−Ｓus3）＋…＋ｆk(Ｓusk-1−Ｓusk)・ｇ（Ｓo
    rg−Ｓusk+1）＋…＋ｆN(ＳusN-1−ＳusN)・ｇ（Ｓorg
    −ＳusN+1）｝ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：抑制画像信号を作成するための関数 ｇ：補助画像信号を変換して倍率信号を作成するための
    関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項４８記載のコ
    ンピュータ読取り可能な記録媒体。
51. 【請求項５１】 前記所定の変換処理は、前記少なく
    とも１つの帯域制限画像信号である被変換帯域制限画像
    信号よりも低い周波数帯域の帯域制限画像信号である低
    周波側帯域制限画像信号を、原点を通り該原点における
    傾きが略０で、処理される値が大きくなるにつれて該傾
    きが漸増する非線形関数に基づいて変換することにより
    前記被変換帯域制限画像信号の補助画像信号を作成する
    手順と、 該補助画像信号を前記被変換帯域制限画像信号に加算す
    ることにより複合帯域制限画像信号を作成する手順と、 該複合帯域制限画像信号を、該複合帯域制限画像信号の
    信号値の絶対値に基づいて定まる該絶対値以下の値とな
    るように変換することにより前記変換画像信号を作成す
    る手順とを有することを特徴とする請求項４１または４
    ２記載のコンピュータ読取り可能な記録媒体。
52. 【請求項５２】 前記帯域制限画像信号を作成する手
    順、前記変換画像信号を作成する手順、前記高周波成分
    に関する信号を作成する手順、および該高周波成分に関
    する信号を前記原画像信号へ加算する手順を、下記の式 Ｓproc＝Ｓorg＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓu
    s2,…ＳusN） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）＝［ｆ1｛（Ｓorg
    −Ｓus1)＋ｇ（Ｓus1−Ｓus2）｝＋ｆ2｛（Ｓus1−Ｓus
    2)＋ｇ（Ｓus2−Ｓus3）｝＋…＋ｆk｛（Ｓusk-1−Ｓus
    k)＋ｇ（Ｓusk−Ｓusk+1）｝＋…＋ｆN｛（ＳusN-1−Ｓ
    usN)＋ｇ（ＳusN−ＳusN+1）｝］ （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜N+1）：非鮮鋭マスク画像信号 ｆk(k=1〜N）：複合帯域制限画像信号を変換するための
    関数 ｇ：帯域制限画像信号を変換して補助画像信号を作成す
    るための関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
    数） にしたがって行うことを特徴とする請求項５１記載のコ
    ンピュータ読取り可能な記録媒体。
53. 【請求項５３】 前記所定の変換処理は、前記原画像
    を撮影した際の撮影部位に応じて、前記帯域制限画像信
    号を変換する手順を有することを特徴とする請求項４１
    から５２のいずれか１項記載のコンピュータ読取り可能
    な記録媒体。