(57)【特許請求の範囲】
【請求項１】
入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成し、該各帯域制限画像信号に対応して設定された複数の変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、前記入力された画像信号に対して所定の処理を施して処理済画像信号を得る画像処理方法において、
基準解像度を有する基準解像度画像について前記変換処理を行なう際に用いられる、該基準解像度画像についての前記各帯域制限画像信号に対応する基準変換関数を用意し、
前記基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、前記処理用画像の解像度と前記基準解像度の違いに基づいて、前記処理用画像についての処理済画像信号の周波数応答特性が、前記基準解像度画像についての処理済画像信号の周波数応答特性と略同じになるように前記基準変換関数を補正し、この補正により得た関数を前記処理用画像についての変換関数として設定するものであり、
前記解像度が、画素密度を表すものであることを特徴とする画像処理方法。
【請求項２】
前記基準変換関数の補正を、前記処理用画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性が、該処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、前記基準解像度画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性とほぼ一致するように行なうことを特徴とする請求項１項記載の画像処理方法。
【請求項３】
前記処理用画像の解像度に関する情報を取得し、該情報に基づいて、前記処理用画像についての変換関数を設定することを特徴とする請求項１または２記載の画像処理方法。
【請求項４】
前記入力された画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、該複数のボケ画像信号および前記入力された画像信号に基づいて前記帯域制限画像信号を作成することを特徴とする請求項１から３いずれか１項記載の画像処理方法。
【請求項５】
前記基準変換関数は、非線形関数であることを特徴とする請求項１から４いずれか１項記載の画像処理方法。
【請求項６】
前記所定の処理は、周波数強調処理であることを特徴とする請求項１から５いずれか１項記載の画像処理方法。
【請求項７】
前記所定の処理は、ダイナミックレンジ圧縮処理であることを特徴とする請求項１から５いずれか１項記載の画像処理方法。
【請求項８】
入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成する画像信号作成手段と、該各帯域制限画像信号に対応する複数の変換関数を設定し、設定した変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、前記入力された画像信号に対して所定の処理を施して処理済画像信号を得る変換処理手段とを備えた画像処理装置において、
前記変換処理手段が、基準解像度を有する基準解像度画像について前記変換処理を行なう際に用いられる、該基準解像度画像についての前記各帯域制限画像信号に対応する基準変換関数を保持すると共に、前記基準解像度とは異なる解像度を有する処理用画像を表す画像信号が前記画像信号作成手段に入力されたとき、前記処理用画像の解像度と前記基準解像度の違いに基づいて、前記処理用画像についての処理済画像信号の周波数応答特性が、前記基準解像度画像についての処理済画像信号の周波数応答特性と略同じになるように前記基準変換関数を補正し、この補正により得た関数を前記処理用画像についての変換関数として設定するものであり、
前記解像度が、画素密度を表すものであることを特徴とする画像処理装置。
【請求項９】
前記変換処理手段が、前記基準変換関数の補正を、前記処理用画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性が、該処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、前記基準解像度画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性とほぼ一致するように行なうものであることを特徴とする請求項８記載の画像処理装置。
【請求項１０】
前記処理用画像の解像度に関する情報を取得する解像度情報取得手段をさらに備え、前記変換処理手段は、該解像度情報取得手段により取得された情報に基づいて、前記処理用画像についての変換関数を設定するものであることを特徴とする請求項８または９記載の画像処理装置。
【請求項１１】
前記画像信号作成手段は、前記入力された画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、該複数のボケ画像信号および前記入力された画像信号に基づいて前記帯域制限画像信号を作成するものであることを特徴とする請求項８から１０いずれか１項記載の画像処理装置。
【請求項１２】
前記基準変換関数は、非線形関数であることを特徴とする請求項８から１１いずれか１項記載の画像処理装置。
【請求項１３】
前記所定の処理は、周波数強調処理であることを特徴とする請求項８から１２いずれか１項記載の画像処理装置。
【請求項１４】
前記所定の処理は、ダイナミックレンジ圧縮処理であることを特徴とする請求項８から１２いずれか１項記載の画像処理装置。
【請求項１５】
入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成する手順と、該各帯域制限画像信号に対応して設定された複数の変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、前記入力された画像信号に対して所定の処理を施して処理済画像信号を得る手順とを有する画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読取り可能な記録媒体において、
基準解像度を有する基準解像度画像について前記変換処理を行なう際に用いられる、該基準解像度画像についての前記各帯域制限画像信号に対応する基準変換関数を用意する手順と、
前記基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、前記処理用画像の解像度と前記基準解像度の違いに基づいて、前記処理用画像についての処理済画像信号の周波数応答特性が、前記基準解像度画像についての処理済画像信号の周波数応答特性と略同じになるように前記基準変換関数を補正し、この補正により得た関数を前記処理用画像についての変換関数として設定する手順を有し、
前記解像度が、画素密度を表すものであることを特徴とするコンピュータ読取り可能な記録媒体。
【請求項１６】
前記変換関数を設定する手順が、前記基準変換関数の補正を、前記処理用画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性が、該処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、前記基準解像度画像を表す画像信号を処理して得た前記処理済画像信号により表される画像の周波数応答特性とほぼ一致するように行なう手順であることを特徴とする請求項１５項記載のコンピュータ読取り可能な記録媒体。
【請求項１７】
前記処理用画像の解像度に関する情報を取得する手順をさらに有し、前記変換関数を設定する手順が、該情報に基づいて、前記処理用画像についての変換関数を設定する手順であることを特徴とする請求項１５または１６記載のコンピュータ読取り可能な記録媒体。
【請求項１８】
前記複数の画像信号を作成する手順は、前記入力された画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、該複数のボケ画像信号および前記入力された画像信号に基づいて前記帯域制限画像信号を作成する手順であることを特徴とする請求項１５から１７いずれか１項記載のコンピュータ読取り可能な記録媒体。
【請求項１９】
前記基準変換関数は、非線形関数であることを特徴とする請求項１５から１８いずれか１項記載のコンピュータ読取り可能な記録媒体。
【請求項２０】
前記所定の処理は、周波数強調処理であることを特徴とする請求項１５から１９いずれか１項記載のコンピュータ読取り可能な記録媒体。
【請求項２１】
前記所定の処理は、ダイナミックレンジ圧縮処理であることを特徴とする請求項１５から１９いずれか１項記載のコンピュータ読取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像信号に対して所定の周波数成分を強調するなどの処理を行なうための画像処理方法および装置並びに画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読取り可能な記録媒体に関するものである。
【背景技術】
【０００２】
従来、本出願人により、非鮮鋭マスク画像信号（以下、ボケ画像信号という）を用いて周波数強調処理あるいはダイナミックレンジ圧縮処理などを行って放射線画像の診断性能を向上させる数々の画像処理方法および装置が提案されている（特開昭55-163472号、同55-87953号、特開平3-222577号、特開平10-75395号、同10-171983号など）。例えば周波数強調処理は、原画像信号Ｓorgからボケ画像信号Ｓusを引いたものに強調係数βを乗じたものを、原画像信号Ｓorgに加算することにより、原画像信号の所定の空間周波数成分を強調するものである。これを式で表すと下記の式（１）のようになる。
Ｓproc＝Ｓorg＋β×（Ｓorg−Ｓus） …（１）
（Ｓproc：周波数強調処理された信号、Ｓorg：原画像信号、Ｓus：ボケ画像信号、β：強調係数）
【０００３】
また、特開平10-75395号には、原画像信号に加算する加算信号の周波数応答特性を調整することにより、周波数強調処理された信号に対してアーチファクトが発生することを防止する方法が提案されている。この方法とは、まず鮮鋭度の異なる、すなわち周波数応答特性の異なる複数のボケ画像信号を作成し、そのボケ画像信号および原画像信号の中の２つの信号の差分をとることにより、原画像信号の、ある限られた周波数帯域の周波数成分を表す複数の帯域制限画像信号（バンドパス信号）を作成し、さらにその帯域制限画像信号をそれぞれ異なる変換関数によって所望の大きさとなるように変換してから、その複数の抑制された帯域制限画像信号を積算することにより上記加算信号を作成するものである。この処理は例えば下記の式（２）により表すことができる。
Ｓproc＝Ｓorg＋β（Ｓorg）×Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝ｆ１(Ｓorg−Ｓus1)＋ｆ２(Ｓus1−Ｓus2)＋…
＋ｆk(Ｓusk-1−Ｓusk)＋…＋ｆＮ(ＳusN-1−ＳusN)  …（２）
（但し、Ｓproc：高周波成分が強調された画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆｋ(k=1〜N）：各帯域制限画像信号を変換する変換関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数）
【０００４】
これらの周波数強調処理においては、帯域制限画像信号を変換する変換関数のパラメータを変更することによって原画像信号に加算する加算信号の周波数応答特性を調整することができる。このため、各変換関数の定義次第で、アーチファクトの発生防止など所望の周波数応答特性を有する処理済画像信号を得ることができる。
【０００５】
さらに、特開平 10-171983号には、周波数強調処理とダイナミックレンジ圧縮処理とを同時に施す場合において、処理が施された信号に対してアーチファクトが発生することを防止する方法が提案されている。この方法とは、上記特開平10-75395号に記載されたように、複数の帯域制限画像信号を作成し、この帯域制限画像信号に基づいて原画像信号の高周波成分に関する信号（高周波成分信号）および低周波成分に関する信号（低周波成分信号）を得、原画像信号にこれらの高周波成分に関する信号および低周波成分に関する信号を加算することにより周波数強調処理を施すとともにダイナミックレンジ圧縮処理を施すようにしたものである。この処理は例えば下記の式（３）により表すことができる。
Ｓproc＝Ｓorg ＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…ＳusN））…（３）
Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝｛ｆｕ１(Ｓorg−Ｓus1)＋ｆｕ２(Ｓus1 −Ｓus2)＋…
＋ｆｕｋ(Ｓusk-1−Ｓusk )＋…＋ｆｕＮ(ＳusN-1−ＳusN)｝
Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝｛ｆｄ１(Ｓorg−Ｓus1)＋ｆｄ２(Ｓus1 −Ｓus2)＋…
＋ｆｄｋ(Ｓusk-1−Ｓusk )＋…＋ｆｄＮ(ＳusN-1−ＳusN)｝
（但し、Ｓproc：処理済画像信号
Ｓorg：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆｕｋ(k=1〜N）：高周波成分信号を得るために使用する変換関数
ｆｄｋ(k=1〜N）：低周波成分信号を得るために使用する変換関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められる                                ダイナミックレンジ圧縮用の信号（ＤはＳorg−Ｆdrcを変換する関数））
【０００６】
一方、上述した周波数強調処理などに用いられるボケ画像信号は、原画像信号の画素に対して所定間隔毎に所定のフィルタリング処理を施すことによって画素を間引きし、その後間引きした数分の画素を所定の補間方法によって補間することにより作成される。このフィルタリング処理としては、ローパスフィルタにより原画像信号の高周波成分を取り除くような処理、具体的にはフィルタ内の画素値の平均値や、荷重平均値を求めるといった処理が行われている。上記特開平10-75395号などにおいてボケ画像信号を得るために行われるフィルタリング処理では、原画像信号に対してフィルタリング処理を施して得られた画素数の少ない信号に対してさらにこのフィルタリング処理を施して、フィルタリングの各段階で得られた画素数の少ない画像信号をそれぞれ原画像信号と同じ画素数となるように補間することにより、複数の異なるボケ画像信号を作成している。
【０００７】
ボケ画像信号は、上記のように原画像信号に基づいて作成されるものであるが、この原画像信号は、読取装置により原画像信号を所定の読取密度で読取って所定の画素密度を有する画像を再現可能なデジタル信号としたものである。一般に、デジタル化された画像信号を例えばプリント出力として再生する場合、画素密度によって決定されるある周波数（ナイキスト周波数）以下の周波数成分は正しく再生されることが知られている。すなわち、読取密度すなわち画素密度は再生時に必要となる画質レベルを考慮して決定されるため、必ずしも一定であるとは限らない。
【０００８】
例えば、放射線画像読取再生システムでは、蓄積性蛍光体シートに記録された人体の放射線画像をレーザビーム走査によりデジタル画像信号として読取っているが、読取密度すなわち画素密度は蓄積性蛍光体シートの大きさによって異なったり、ユーザの設定により任意の値に変更することができるものである。
【０００９】
ここで、読み取りの際の画素密度（読取密度）すなわちナイキスト周波数が異なり、該信号に基づいて得られる画像の解像度がそれぞれ異なる画像信号に対し、同一のローパスフィルタによるフィルタリング処理、および同一の補間方法による補間処理を施した場合、得られる帯域制限画像信号の周波数特性、具体的には帯域制限画像信号の周波数帯域は画素密度によって異なることとなる。これにより、例えば１つの原画像を２種類の読取密度で読み取って２種類の画素密度を有する原画像信号を得た場合、同じボケ画像信号を用いて帯域制限画像信号を得て周波数強調処理あるいはダイナミックレンジ圧縮処理を行っても、強調される周波数帯域あるいは圧縮される周波数帯域は２種類の原画像信号でそれぞれ異なってしまうという問題が生じる。
【００１０】
このため、原画像信号の画素密度に関する情報を得、その情報に基づいて複数のフィルタ係数リストからフィルタ係数を選択し、選択されたフィルタ係数のフィルタによって原画像信号に対してフィルタリング処理を行ってボケ画像信号を得るようにした方法が提案されている（特開平10-63837号）。ここで例えば読取密度が５本／ｍｍと6.7本／ｍｍの原画像信号に対して同一のローパスフィルタによりフィルタリング処理を施すことにより得られる帯域制限画像信号はその周波数帯域が異なるものとなるが、この方法によれば各原画像信号に対してそれぞれ異なるローパスフィルタによりフィルタリング処理を施すことにより、各原画像信号から得られる帯域制限画像信号の周波数帯域を略一致させることができる。したがって、画素密度に拘わらず同じ周波数特性のボケ画像信号を作成することができ、これにより同じ周波数特性の帯域制限画像信号を作成し、所望の処理、例えば上述したような周波数強調処理やダイナミックレンジ圧縮処理などを、常に同じように行なうことができることとなる。
【００１１】
一方、原画像信号を圧縮する（データ数を少なくする）形式としては、ＪＰＥＧ、ＧＩＦ、ＴＩＦＦなど種々の形式が存在するが、近年画像信号を解像度毎に階層的に分解し、各階層毎のデータ（階層データ）を符号化して圧縮する形式が提案されている。この圧縮形式は、具体的には原画像信号をウェーブレット変換などにより原画像の１／２ｎ倍の解像度を有する複数の解像度毎の階層画像信号に分解し、この分解された各解像度毎の階層画像信号を階層順に符号化して１つのファイルとして圧縮するものである。
【００１２】
この圧縮方式は以下のような特徴を有する。
（１）従来のＪＰＥＧで用いられているＤＣＴ方式のように、画像信号をブロック毎に処理していないため、ブロック歪みのようなアーチファクトが生じない。
（２）画像信号が階層的に符号化されているため、画像信号の転送の際に必要な解像度の情報のみを転送すればよく、効率的な画像転送が可能となる。
（３）画像信号が多重解像度あるいは多重分解能に分解されているため、周波数強調処理など種々の画像処理を比較的簡単に行なうことができる。
（４）多重解像度解析による空間と周波数との同時分解が可能であり、符号化効率に大きく影響を与える低周波数領域に対しては広い範囲で直交変換を行ない、高周波領域に対しては狭い範囲で直交変換が可能となるため、画像中のエッジ周辺部に量子化ノイズが発生しても、その空間的広がりを抑えることができる。このため、ノイズが知覚されにくい。
【００１３】
また、イーストマンコダック社が提案するＦｌａｓｈＰｉｘファイルのように、１つのファイル内に複数の性質の異なるデータを記憶することができるファイル形式が提案されているが、このようなＦｌａｓｈＰｉｘ規格のファイルにも、多重解像度に分解された階層画像信号を保管することも可能である。
【００１４】
このように、原画像信号を多重解像度に分解することにより、原画像の１／２ｎ 倍の解像度を有する複数の画像を表す階層画像信号から原画像信号を構成することができる。そして、このような多重解像度に分解された信号を用いて圧縮された画像は、その用途によって原画像と同じ解像度レベルで復元して使用されたり、あるいは途中の解像度レベルで復元して使用されたりする。例えば、プリンタのように高画質の画像を再生する必要がある場合には、最高解像度までの画像を表す階層画像信号に基づいて画像を復元することにより原画像と同一解像度を有する高画質の画像を再現可能な画像信号を得ることができる。また、ＣＲＴのようにプリンタほど高解像度の画像を再現する必要がない場合には、原画像よりも低解像度の画像を表す階層画像信号に基づいて画像信号を復元し、さらに必要であれば拡大、縮小することにより、原画像よりは解像度が低いもののＣＲＴの解像度に適した画像を再生することができる。
【発明の開示】
【発明が解決しようとする課題】
【００１５】
しかしながら、多重解像度に分解された信号を用いて圧縮された画像を、原画像と同じ解像度レベルで復元しないで途中の解像度レベルで復元して使用する場合、復元された画像の解像度は、原画像の解像度に対して１／２ｎ 倍の関係を有するものの、上記特開平10-63837号に記載のような、予め想定されている読取密度に応じた解像度とは異なるため、上記特開平10-63837号に記載の方法を用いることは困難であり、原画像信号について設定されるべき変換関数を用いて所望の変換処理を行なうこととなる。この場合、原画像信号と原画像信号よりも低い解像度の画像を表す階層画像信号とでは周波数応答特性が異なるため、上述した周波数強調処理を多重解像度に分解された画像信号における原画像よりも低解像度の画像を表す画像信号に施す場合に、原画像の解像度に対応する変換関数を低解像度の画像に適用したのでは、原画像信号に対して周波数強調処理を施した場合と異なる周波数応答特性を有する画像となるおそれがある。一方この場合においても、上記特開平10-63837号に記載の方法を適用するべく、各レベルの解像度に応じた変換関数を用意し、解像度に応じて変換関数を選択して周波数強調処理を施すことが考えられる。しかしながら、各解像度に応じて変換関数を用意したのでは、変換関数の数が膨大なものとなり、変換関数の管理が非常に煩わしいものとなる。これは原画像信号を多重解像度に分解した画像信号に対して周波数強調処理を施す場合のみならず、原画像よりも低い解像度を有する画像を再生するに際して周波数強調処理を施す場合にも生じる問題である。
【００１６】
また、原画像と同じ解像度レベルまで至らない途中の解像度レベルで復元することに限らず、原画像と同じ解像度レベルで復元した原画像と同サイズの画像を拡大し、この拡大により得た画像に対して原画像について設定されるべき変換関数を用いて強調処理などを施す場合においても同様の問題を生じる。例えばＣＲＴ上で画像を見て画像の一部分を詳細に読影する場合には、画像の該当する部分を拡大して表示したりあるいはフィルム出力することもある。この場合、同じ変換関数を用いて処理しているにも拘わらず、拡大画像の空間解像度がアップしているために、拡大後に強調処理などを施した画像の画質と、原画像サイズで通常通り強調処理などを施した画像の画質とが異なり、読影に影響をおよぼすことがある。
【００１７】
さらに、上述したような周波数強調処理などを行なう画像処理装置においては、放射線画像読取再生システムにおいて得られた、予め定められた読取密度（解像度の一指標）で読み取られた原画像信号のみならず、入力モダリティの種類などに応じて読取密度が異なる画像を表す原画像信号が入力される、つまり入力画像の解像度が多様化することもある。このような場合に予め定められた解像度（以下基準解像度とする）に対応する変換関数により、基準解像度とは異なる解像度の画像を表す画像信号に対して上述した周波数強調処理などを施すと、得られる画像が、基準解像度の画像を表す画像信号に対して処理を施した場合と異なる周波数応答特性を有するものとなるおそれがある。この場合、解像度に応じて変換関数を用意し、解像度に応じた変換関数を選択して周波数強調処理などを施すことにより、処理を施す画像の解像度に拘わらず同一の周波数応答特性を有する画像を再現可能な処理済画像信号を得ることができる。しかしながら、解像度に応じて変換関数を用意したのでは、変換関数の数が膨大なものとなり、変換関数の管理が非常に煩わしいものとなる。
【００１８】
本発明は上記事情に鑑みなされたものであり、処理する画像の解像度に拘わらず、周波数強調処理などを常に同じように行なうことができる画像処理方法および装置並びに画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読取り可能な記録媒体を提供することを目的とするものである。
【課題を解決するための手段】
【００１９】
本発明による第１の画像処理方法は、入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成し、該各帯域制限画像信号に対応して設定された複数の変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、入力された画像信号に対して所定の処理を施して処理済画像信号を得る画像処理方法において、
基準解像度を有する基準解像度画像について変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を用意し、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、基準変換関数に対応する各帯域制限画像信号の周波数帯域が処理用画像についての各帯域制限画像信号の周波数帯域に最も近くなるように基準変換関数を相対移動させ、この移動させた基準変換関数を処理用画像についての変換関数として設定することを特徴とするものである。
【００２０】
なお「解像度」という言葉は、一般的には非常に広い意味を有し、例えば鮮鋭度を意味するものとして使用されたり、画像（あるいは画像信号）の画素密度を意味するものとして使用されたりする場合もあるが、本願発明において用いられる解像度は鮮鋭度に関するものを意味するものとする。また、画像（あるいは画像信号）の総画素数を意味するものとして画像サイズという言葉を用いる。ここで「鮮鋭度に関するものを意味する」とは、例えば画素密度を増減（この場合画像サイズも増減する）して鮮鋭度を増減したり、あるいは画素密度を同じに維持しつつ再生する（できる）周波数範囲を拡大・縮小したり信号のレスポンスを変更して鮮鋭度を増減することのいずれをも含む意味である。したがって、例えば、基準解像度画像を画素間引きを用いたフィルタリング処理などによって低解像度画像にしたときには、画素密度が低下されており、低解像度画像は基準解像度画像よりも画像サイズが小さくなる。一方、前記画素間引きを用いたフィルタリング処理によって得た低解像度画像を補間拡大して得た鮮鋭度の低いボケ画像は、画素密度が基準解像度画像と同じで、画像サイズも同じものである。換言すれば、前記画素間引きを用いたフィルタリング処理によって得た低解像度画像は基準解像度画像よりも解像度の低い（小さい）画像であると共に画像サイズが小さい画像である一方、ボケ画像は基準解像度画像よりも解像度の低い（小さい）画像であると共に画像サイズが同じ画像である。
【００２１】
上記において「各帯域制限画像信号に対応して設定された複数の変換関数」は、各帯域制限画像信号が６つの周波数帯域からなるものである場合は、各周波数帯域に対応して６つの変換関数となる。
【００２２】
なお、本発明においては、基準変換関数としては、線形関数、非線型関数、定数など種々の態様のものを用いることができる（後述する第２の方法においても同様）。
【００２３】
また、「所定の処理」としては、例えば、上記式（２）に示すような原画像信号に含まれる特定の周波数成分を強調する周波数強調処理や、上記式（３）に示すような原画像の最高濃度と最低濃度との差すなわちダイナミックレンジを狭めるように高濃度域もしくは低濃度域あるいは高濃度域および低濃度域の双方のコントラストを下げるダイナミックレンジ圧縮処理などが挙げられる（後述する第２の方法においても同様）。
【００２４】
本発明による第１の画像処理方法において、「基準変換関数を相対移動」させて処理用画像についての変換関数を設定するに際しては、基準変換関数と帯域制限画像信号の組数、並びに基準解像度画像についての帯域制限画像信号と処理用画像についての帯域制限画像信号の共通する周波数範囲の数に基づいて、前記相対移動を行なうのが望ましく、具体的には、以下のようにするとよい。すなわち、処理用画像についての変換関数と該変換関数に対応する各帯域制限画像信号の組数を、処理用画像の解像度が基準解像度の２−（ｐ＋１）倍から２−ｐ倍（ｐは正の整数）の範囲内にあるときにはｐ個だけ減少させ、処理用画像の解像度が基準解像度の２ｑ−１ 倍から２ｑ 倍（ｑは正の整数）の範囲内にあるときにはｑ個だけ増加させるとよい。なお、第２の解像度が第１の解像度の２−１倍から１倍（上記においてｑ，ｐが０の場合に相当）の範囲内にあるときには前記組数を増減させない。
【００２５】
ここで処理用画像の解像度が基準解像度の２−（ｐ＋１）倍から２−ｐ倍の範囲内にあるときにｐ個だけ減少させるのは、解像度の関係が上述のような状態にあるときには、各帯域制限画像信号のうちの共通する周波数範囲の数がｐ個分だけ減少するため、処理用画像について強調処理など行なう際には、各帯域制限画像信号のうちの、基準解像度画像についての強調処理には使用されない低周波側の周波数帯域に対応するｐ個分の基準変換関数を用いる必要がなくなる（むしろ、用いると不具合を生じる虞れがある）からである。
【００２６】
一方、処理用画像の解像度が基準解像度の２ｑ−１ 倍から２ｑ 倍の範囲内にあるときには、帯域制限画像信号としては、基準解像度画像についての強調処理には使用される低周波側の周波数帯域成分が不足するため、単に基準変換関数を相対移動させるだけでは低周波側の変換関数に用いるべき帯域制限画像信号が存在せず、画質に影響を与える低周波側の特性を基準解像度画像についての強調処理などによる特性と同じにできないため、処理用画像について強調処理など行なう際には帯域制限画像信号の数をｑ個だけ増やす必要があるからである。そして、この場合、基準変換関数を相対移動させただけでは、処理用画像についての変換関数の数が生成される帯域制限画像信号の数よりも少なくなるので、生成される帯域制限画像信号のうち、基準解像度画像についての強調処理には使用されない高周波側の周波数帯域に対応する変換関数（不足する分）を新たに定義するのが望ましく、結果的に、処理用画像についての変換関数と該変換関数に対応する各帯域制限画像信号の組数を、ｑ個だけ増加させることになる。
【００２７】
また、上記第１の画像処理方法においては、処理用画像の解像度が基準解像度よりも低い場合は、複数の変換関数のうち処理用画像の解像度以下の周波数帯域に対応する変換関数を、前記変換処理を行なう変換関数として設定することもできる。
【００２８】
ここで、「処理用画像の解像度以下の周波数帯域に対応する変換関数」とは、例えば処理用画像の解像度が基準解像度の１／２ｋ 倍である場合には、各周波数帯域毎の変換関数のうちｋ個分低周波数帯域側の変換関数のことをいう。例えば基準解像度画像を表す画像信号に対して６つの変換関数ｆ１ 〜ｆ６ （ｆ１ が最高周波数帯域に対応）があり、処理用画像の解像度が基準解像度の１／２である場合には、最高周波数帯域に対応する変換関数以外の５つの変換関数ｆ２ 〜ｆ６ のことをいう。そして処理用画像を表す画像信号から得られる６つの帯域制限画像信号のうちの、最も低周波側のものを除く５つの帯域制限画像信号を５つの変換関数により変換して変換画像信号が得られることとなる。
【００２９】
また、上記第１の画像処理方法において、基準解像度画像の画像サイズとは異なる画像サイズの出力画像を得る際には、出力画像の画像サイズに近く（好ましくは最も近く）なるように基準解像度画像の画像サイズを２２ｋ（ｋは整数）倍した変倍画像を処理用画像として変換処理および所定の処理を施し、その後、処理済画像信号が表す画像サイズが２２ｋ倍の画像に対して、出力画像の画像サイズと等しくなるように変倍処理を施す、すなわち前記所定の処理が施された処理済画像信号に対して変倍処理を施して変倍処理済画像信号を得ることが望ましい。
【００３０】
「変倍画像を処理用画像として変換処理および所定の処理を施す」とは、変倍画像を表す変倍画像信号を前記処理用画像を表す画像信号として入力し、この入力された画像信号（変倍画像信号）から各帯域制限画像信号を作成し、設定された変換関数に基づいて変換処理を施し、さらに所定の処理を施すことを意味する。
【００３１】
「画像サイズが２２ｋ倍」とは、主走査方向および副走査方向に画素密度をそれぞれ２ｋ 倍することによって画像サイズを２２ｋ倍とすることを意味する。なお、ｋは整数、すなわちゼロ、あるいは正負いずれかの整数であり、基準解像度画像の画像サイズを２２ｋ倍するとは、変倍画像を、拡大（ｋ＞０）あるいは縮小（ｋ＜０）することだけでなく、基準解像度画像の画像サイズそのまま（ｋ＝０）とすることも含む意味である。
【００３２】
また、「出力画像の画像サイズと等しくなるように変倍処理を施す」場合においても、「変倍」は、画像サイズを、拡大（変倍率＞１）あるいは縮小（変倍率＜１）することだけでなく、画像サイズをそのまま（変倍率＝１）とすることを含む。例えば、出力画像の画像サイズが基準解像度画像の画像サイズの２２ｋであるときには、前記所定の処理を施して得られた処理済画像信号が表す画像の画像サイズそのものが出力画像の画像サイズと同じになるので、「出力画像の画像サイズと等しくなるように変倍処理を施す」に際しては、変倍率を１倍にすればよく、実際には処理済画像信号をそのまま変倍処理済画像信号とすればよい（つまり変倍処理を行なわなくてもよい）。
【００３３】
本発明による第２の画像処理方法は、入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成し、該各帯域制限画像信号に対応して設定された複数の変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、入力された画像信号に対して所定の処理を施して処理済画像信号を得る画像処理方法において、
基準解像度を有する基準解像度画像について前記変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を用意し、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、処理用画像の解像度と基準解像度の違いに基づいて、処理用画像についての処理済画像信号のレスポンスが、基準解像度画像についての処理済画像信号のレスポンスと略同じになるように基準変換関数を補正（新たな関数を算出）し、この補正により得た関数を処理用画像についての変換関数として設定することを特徴とするものである。
【００３４】
「処理用画像の解像度と基準解像度の違いに基づいて、処理用画像についての処理済画像信号のレスポンスが、基準解像度画像についての処理済画像信号のレスポンスと略同じになるように基準変換関数を補正」するに際しては、例えば、基準解像度画像を表す画像信号および処理用画像を表す画像信号から得られる、それぞれの各帯域制限画像信号、所望の周波数応答特性、および各画像それぞれに設定されるべき変換関数のパラメータの関係式を、各帯域制限画像信号と周波数応答特性を既知の値とし、変換関数のパラメータを変数とする連立方程式として解くことによりパラメータを求めるなどの方法を用いることができる。あるいは、処理用画像についての変換関数および処理用画像を表す画像信号から得られる複数の画像信号の周波数帯域から、内挿、外挿などの演算を行なうことにより算出する方法であってもよい。
【００３５】
上記第２の画像処理方法においては、基準変換関数の補正（新たな関数の算出）を、入力された処理用画像を表す画像信号を処理して得た処理済画像信号により表される画像の特性が、処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、基準解像度画像を表す画像信号を処理して得た処理済画像信号により表される画像の特性とほぼ一致するように行なうことが望ましい。
【００３６】
さらに、本発明による第１および第２の画像処理方法においては、処理用画像の解像度に関する解像度情報を取得し、該解像度情報に基づいて、処理用画像についての変換関数を設定することが好ましい。
【００３７】
なお、解像度情報を取得するに際しては、ユーザ（オペレータ）がキーボードから数値として入力してもよいし、操作画面に数種類の解像度情報を表示してユーザに選択させるといった形態でもよいし、あるいは処理用画像を表す画像信号に解像度情報を添付しておき、処理用画像を表す画像信号を処理する際に解像度情報を自動的に取得するようにしてもよく、解像度情報を取得することができれば、どのような形態であってもよい。
【００３８】
なお、「解像度情報」は、画像（あるいは画像信号）の鮮鋭度に関する情報である限りどのようなものを用いてもよく、例えば蓄積性蛍光体シートに記録された放射線画像を読み取る際の読取密度（例えばｄｐｉや本／mm）のみならず、例えば画像信号を得る際のサンプリング周波数であってもよい。
【００３９】
なお、本発明における各帯域制限画像信号は、それぞれが、入力された処理用画像を表す画像信号（原画像信号）の所定の周波数帯域成分を表すものである限りどのような方法を用いて生成してもよく、例えば、原画像信号をウェーブレット変換やラプラシアンピラミッドなどの手法により多重解像度展開して作成することができる。
【００４０】
また、この帯域制限画像信号は、入力された画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、該複数のボケ画像信号および入力された画像信号に基づいて作成することもできる。なお、上述したように原画像信号を多重解像度展開した場合には、各解像度毎の画像を表す階層画像信号を低解像度画像信号とし、低解像度画像信号を補間拡大することによりボケ画像信号を作成することもできる。
【００４１】
ここで「ボケ画像信号」とは、入力された処理用画像を表す画像信号（原画像信号）と画素数は同じであるが原画像信号よりも鮮鋭度（解像度）が低い画像を表す画像信号である。ボケ画像信号は、まず原画像信号の画素に対して所定間隔毎に所定のフィルタリング処理、線形補間処理、あるいは単純に画素数を減少させる処理などを施すことによって画素を間引きし、そのようにして得た画像信号に対して同様のフィルタリング処理を繰り返してさらに画素数を少なくして解像度が原画像（処理用画像）の１／２ｎ 倍となる低解像度画像を表す低解像度画像信号を複数作成し、そのそれぞれに対して所定の補間方法により原画像と画素数が同じになるように補間処理を施すことにより作成する。ここで、フィルタリング処理や補間は入力された処理用画像の解像度（例えば画素密度）に基づいて定められたフィルタ係数を有するフィルタにより行われるものであり、その具体的手法については、例えば上述した特開平10-63837号に記載されている方法など、一般に広く使用されている種々の方法を用いることができる。
【００４２】
この場合、「帯域制限画像信号」は、例えば隣接する周波数帯域のボケ画像信号同士で差分を求めて作成してもよいし、原画像信号と各ボケ画像信号の差分を求めて作成してもよい。あるいは原画像信号とボケ画像信号の他の組み合わせで差分を求めて作成することもできる。なお、変換画像信号の作成、および周波数強調処理やダイナミックレンジ圧縮処理は例えば上記（２）式や（３）式で表すことができる。
【００４３】
本発明による第１の画像処理装置は、上記第１の画像処理方法を実施する装置、すなわち、入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成する画像信号作成手段と、該各帯域制限画像信号に対応する複数の変換関数を設定し、設定した変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、前記入力された画像信号に対して所定の処理を施して処理済画像信号を得る変換処理手段とを備えた画像処理装置であって、
変換処理手段を、基準解像度を有する基準解像度画像について変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を保持すると共に、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が画像信号作成手段に入力されたとき、基準変換関数に対応する各帯域制限画像信号の周波数帯域が処理用画像についての各帯域制限画像信号の周波数帯域に最も近くなるように基準変換関数を相対移動させ、この相対移動させた基準変換関数を処理用画像についての変換関数として設定するものとしたことを特徴とするものである。
【００４４】
上記第１の画像処理装置においては、変換処理手段を、処理用画像についての変換関数と該変換関数に対応する各帯域制限画像信号の組数を、処理用画像の解像度が基準解像度の２−（ｐ＋１）倍から２−ｐ倍（ｐは正の整数）の範囲内にあるときにはｐ個だけ減少させ、処理用画像の解像度が基準解像度の２ｑ−１ 倍から２ｑ 倍（ｑは正の整数）の範囲内にあるときにはｑ個だけ増加させるものとすることが望ましい。
【００４５】
また上記第１の画像処理装置においては、変換処理手段を、処理用画像の解像度が基準解像度よりも低い場合は、複数の変換関数のうち処理用画像の解像度以下の周波数帯域に対応する変換関数を、変換処理を行なう変換関数として設定するものとすることもできる。
【００４６】
さらに、本発明による第１の画像処理装置においては、基準解像度画像の画像サイズとは異なる画像サイズの出力画像を得ることができるように、さらに、該出力画像の画像サイズに近い基準解像度画像の２２ｋ（ｋは整数）倍となる第１の倍率と、該第１の倍率で変倍された画像を該画像の画像サイズが前記出力画像の画像サイズと等しくなるように変倍するための第２の倍率とを算出する倍率算出手段と、第１の倍率により基準解像度画像を変倍処理して処理用画像を得る第１の変倍処理手段と、変倍処理を行なって得た処理用画像を表す画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成し、複数の変換関数のうち処理用画像の解像度以下の周波数帯域に対応する変換関数を、変換処理を行なう変換関数として設定して該変換処理を行なって変換画像信号を得、該変換画像信号に基づいて前記処理済画像信号を得るよう画像信号作成手段および変換処理手段を制御する制御手段と、該処理済画像信号を第２の倍率によりさらに変倍処理して出力画像を表す変倍処理済画像信号を得る第２の変倍処理手段とを備えたものとするとよい。
【００４７】
なお、倍率算出手段と第１の変倍処理手段で、出力画像の画像サイズに最も近い、基準解像度画像の画像サイズを２２ｋ倍した変倍画像を得、この変倍画像を処理用画像とする手段が構成される。また第２の変倍処理手段は、処理済画像信号が表す２２ｋ倍の画像に対して出力画像の画像サイズと等しくなるように変倍処理を施す手段として機能するものである。
【００４８】
本発明による第２の画像処理装置は、上記第２の画像処理方法を実施する装置、すなわち、入力された画像信号から互いに異なる帯域制限特性を有する複数の帯域制限画像信号を作成する画像信号作成手段と、該各帯域制限画像信号に対応する複数の変換関数を設定し、設定した変換関数に基づいて変換処理を施して変換画像信号を得、該変換画像信号に基づいて、前記入力された画像信号に対して所定の処理を施して処理済画像信号を得る変換処理手段とを備えた画像処理装置であって、
変換処理手段を、基準解像度を有する基準解像度画像について変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を保持すると共に、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が画像信号作成手段に入力されたとき、処理用画像の解像度と基準解像度の違いに基づいて、処理用画像についての処理済画像信号のレスポンスが、基準解像度画像についての処理済画像信号のレスポンスと略同じになるように基準変換関数を補正し、この補正により得た関数を処理用画像についての変換関数として設定するものとしたことを特徴とするものである。
【００４９】
本発明による第２の画像処理装置においては、変換処理手段を、処理用画像を表す画像信号を処理して得た処理済画像信号により表される画像の特性が、該処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、基準解像度画像を表す画像信号を処理して得た処理済画像信号により表される画像の特性とほぼ一致するように、基準変換関数の補正を行なうものとすることが望ましい。
【００５０】
本発明による上記第１および第２の画像処理装置においては、処理用画像の解像度に関する情報を取得する解像度情報取得手段をさらに備えたものとすると共に、変換処理手段を、該解像度情報取得手段により取得された情報に基づいて、処理用画像についての変換関数を設定するものとするのが望ましい。
【００５１】
また、本発明による上記第１および第２の画像処理装置における基準変換関数は非線形関数であってもよい。
【００５２】
さらに、本発明による上記第１および第２の画像処理装置における所定の処理は、周波数強調処理あるいはダイナミックレンジ圧縮処理であってもよい。
【００５３】
なお、本発明による画像処理方法をコンピュータに実行させるためのプログラムとして、コンピュータ読取り可能な記録媒体に記録して提供してもよい。
【発明の効果】
【００５４】
本発明による第１の画像処理方法および装置によれば、基準解像度を有する基準解像度画像について変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を用意し、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、基準変換関数に対応する各帯域制限画像信号の周波数帯域が処理用画像についての各帯域制限画像信号の周波数帯域に最も近くなるように基準変換関数を相対移動させ、この移動させた基準変換関数を処理用画像についての変換関数として設定しているので、基準解像度画像を表す入力画像信号（基準解像度画像信号）および処理用画像を表す入力画像信号（処理用画像信号）のそれぞれから得られる各帯域制限画像信号のうちの、同じか若しくは近傍の周波数範囲を示す帯域制限画像信号に対して、ほぼ同様の度合いで強調処理などが施されるように、処理用画像についての変換関数を設定することができる。したがって、処理用画像について所定の処理を施して得た処理済画像信号を再生することにより得られる出力画像においては、基準解像度画像（を表す画像信号）と略同じ周波数範囲である限り、その周波数応答特性を基準解像度画像信号から得られた処理済画像信号を再生することにより得られる出力画像の周波数応答特性と略一致させる、つまり、画像の解像度に影響されることなく、ほぼ同様の周波数応答特性を有する画像を再現可能な処理済画像信号を得ることができる。
【００５５】
また、処理用画像についての変換関数と該変換関数に対応する各帯域制限画像信号の組数を、処理用画像の解像度が基準解像度の２−（ｐ＋１）倍から２−ｐ倍（ｐは正の整数）の範囲内にあるときにはｐ個だけ減少させ、処理用画像の解像度が基準解像度の２ｑ−１ 倍から２ｑ 倍（ｑは正の整数）の範囲内にあるときにはｑ個だけ増加させたり、さらには処理用画像の解像度が基準解像度よりも低い場合は、複数の変換関数のうち処理用画像の解像度以下の周波数帯域に対応する変換関数を、変換処理を行なう変換関数として設定するようにすれば、処理用画像についての変換関数として、基準解像度画像について用いる基準変換関数と同一の変換関数を用いることができ、解像度毎に変換関数を用意する必要が無くなり、これにより本発明を適用したシステムの構成を簡易なものとすることができ、また変換関数の管理の煩わしさをなくすことができる。
【００５６】
ここで、基準解像度よりも低い解像度を有する処理用画像に対して周波数強調処理を施す場合に、基準解像度画像についての基準変換関数のうち処理用画像の解像度以下の周波数帯域に対応する変換関数を、変換処理を行なう変換関数として設定するということは、基準解像度画像信号から処理済画像信号を得る際に用いられる各周波数帯域毎の所定変換関数を処理用画像信号の周波数帯域に対応させ、対応させた変換関数により処理用画像信号から得られる互いに帯域制限特性が異なる複数の画像信号の変換を行なうようにすることである。処理用画像の解像度が基準解像度の１／２ｎ 倍である場合には、処理用画像信号から得られる各帯域制限画像信号のピーク周波数は基準解像度画像信号から得られる各帯域制限画像信号のピーク周波数と一致する部分を有するため、基準解像度画像に適用する基準変換関数を処理用画像信号から得られる複数の画像信号の周波数帯域に対応させて変換処理を行なうことにより、処理用画像信号から得られる変換画像信号の周波数応答特性を基準解像度画像信号から得られる変換画像信号の周波数応答特性と同一のものとすることができる。このため、処理用画像信号から得られる処理済画像信号を再生することにより得られる画像においては、その周波数応答特性を基準解像度画像信号から得られた処理済画像信号を再生することにより得られる画像の周波数応答特性と略一致させることができ、これにより解像度に影響されない一定の周波数応答特性を有する画像を再現可能な処理済画像信号を得ることができる。
【００５７】
さらに、本発明による第１の画像処理方法および装置において、基準解像度画像の画像サイズとは異なる画像サイズの出力画像を得る際に、出力画像の画像サイズに近くなるように基準解像度画像の画像サイズを２２ｋ倍した変倍画像を処理用画像として変換処理および所定の処理を施すと、第１の倍率は処理用画像の解像度を基準解像度の２ｋ 倍とするものともなり、上述のように、処理用画像信号から得られる各帯域制限画像信号のピーク周波数は基準解像度画像信号から得られる各帯域制限画像信号のピーク周波数と一致する部分を有することとなる。したがって、基準解像度画像に適用する基準変換関数を処理用画像信号から得られる各帯域制限画像信号の周波数帯域に対応させて変換処理を行なうことにより、処理用画像信号から得られる変換画像信号の周波数応答特性を基準解像度画像信号から得られる変換画像信号の周波数応答特性と略同一のものとすることができる。このため、処理用画像信号から得られる処理済画像信号が表す画像サイズが２２ｋ倍の画像に対して、所望とする出力画像の画像サイズと等しくなるように変倍処理を施すことにより得られる変倍処理済画像信号を再生して得た画像においては、その周波数応答特性を基準解像度画像信号から得られた処理済画像信号を再生することにより得られる画像の周波数応答特性と略一致させることができ、これにより解像度に影響されない一定の周波数応答特性を有する画像を再現可能な変倍処理済画像信号を得ることができる。
【００５８】
また、出力画像の画像サイズを基準解像度画像の画像サイズよりも小さくする場合において、前記第１の倍率が２２ｋ（ｋは負の整数）倍となるときには、前記処理用画像が基準解像度画像の画素数よりも少ない縮小画像となり、この縮小画像について強調処理などの所定の処理を施すことができるので、該所定の処理のための演算時間を短縮することができる。
【００５９】
また、本発明の第２の画像処理方法および装置によれば、基準解像度を有する基準解像度画像について前記変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数を用意し、基準解像度とは異なる解像度を有する処理用画像を表す画像信号が入力されたとき、処理用画像の解像度と基準解像度の違いに基づいて、処理用画像についての処理済画像信号のレスポンスが、基準解像度画像についての処理済画像信号のレスポンスと略同じになるように基準変換関数を補正し、この補正により得た関数を処理用画像についての変換関数として設定しているので、処理用画像信号から得られる処理済画像信号においてはその周波数応答特性を基準変換関数を用いて基準解像度画像信号を処理した場合に得られる処理済画像信号の周波数応答特性と略同一のものとすることができ、これにより得られる処理済画像信号を再生することにより得られる画像においては、その周波数応答特性を解像度に拘わらず略一定のものとすることができる。したがって、解像度に影響されない一定の周波数応答特性を有する画像を再現可能な処理済画像信号を得ることができる。
【００６０】
また、基準変換関数を補正しているので、基準変換関数に基づいて処理用画像信号に対して変換処理を行なう場合の変換関数を算出することができ、解像度毎に複数の変換関数を用意する必要が無くなり、これにより変換関数の管理の煩わしさをなくすことができる。
【００６１】
また、処理用画像信号を処理して得た処理済画像信号により表される画像の特性が、処理用画像の解像度に応じたナイキスト周波数の少なくとも１／５以下の周波数範囲において、基準解像度画像信号を処理して得た処理済画像信号により表される画像の特性とほぼ一致するように前記基準変換関数の補正を行なうと、画質に影響を与える度合いの大きい比較的低周波領域における周波数応答特性を、基準解像度画像信号を処理した場合に得られる周波数応答特性と略同一のものとすることができる。
【発明を実施するための最良の形態】
【００６２】
以下、図面を参照して本発明の実施形態について説明する。なお、本実施形態においては、蓄積性蛍光体シートに記録された人体の放射線画像を読取って得た画像信号に対して、その画像が診断に適した画像となるように、ボケ画像信号を使用して周波数強調処理を施すものであり、処理された画像信号は主としてフィルムに記録され、診断に用いられる。
【００６３】
図１は本発明の第１の実施形態による画像処理装置の構成を示す概略ブロック図である。画像処理装置１は、読取装置などにおいて得られ、後述するように原画像信号Ｓorg を多重解像度空間に変換して符号化することにより得られた原画像信号Ｓorg′ からボケ画像信号を作成するボケ画像信号作成手段２と、特定の周波数を強調するための周波数強調処理を行って中間処理済画像信号Ｓproc′を得る周波数強調処理手段３とを有する。
【００６４】
さらに、この画像処理装置１は、パラメータ設定手段４と拡大率入力手段５と拡大縮小手段６とを備える。パラメータ設定手段４は、周波数強調処理手段３が変換処理に使用する変換関数を設定する手段であり、拡大率入力手段５から入力された拡大率に関する情報Ｍに基づいて変換関数を設定する。拡大率入力手段５は原画像信号Ｓorg の拡大率情報Ｍを得るための手段であり、ユーザが拡大率をキーボードから数値として入力してもよく、操作画面に数種類の拡大率を表示してユーザに選択させるといった形態でもよい。拡大縮小手段６は周波数強調処理手段３において得られた中間処理済画像信号Ｓproc′を拡大率情報Ｍに基づいて拡大縮小して最終的な処理済画像信号Ｓprocを得るものである。
【００６５】
原画像信号Ｓorg は以下のようにして多重解像度空間に変換されて符号化される。まず図２（ａ）に示すように、原画像信号Ｓorg がウェーブレット変換されて複数の解像度毎の４つのデータＬＬ１、ＨＬ０、ＬＨ０およびＨＨ０に分解される。ここで、データＬＬ１は原画像の縦横を１／２に縮小した画像を表し、データＨＬ０、ＬＨ０およびＨＨ０はそれぞれ縦エッジ、横エッジおよび斜めエッジ成分の画像を表すものとなる。そして、４つのデータＬＬ１、ＨＬ０、ＬＨ０およびＨＨ０を逆ウェーブレット変換することにより原画像信号Ｓorgが得られる。次に、図２（ｂ）に示すようにデータＬＬ１をさらにウェーブレット変換して４つのデータＬＬ２、ＨＬ１、ＬＨ１およびＨＨ１を得る。ここで、データＬＬ２はデータＬＬ１の縦横をさらに１／２に縮小した画像を表すものとなり、データＨＬ１、ＬＨ１およびＨＨ１はそれぞれデータＬＬ１の縦エッジ、横エッジおよび斜めエッジ成分の画像を表すものとなる。そしてこれら４つのデータＬＬ２、ＨＬ１、ＬＨ１およびＨＨ１を逆ウェーブレット変換することにより原画像の１／２の解像度の画像を表す低解像度画像信号が得られる。さらに、ウェーブレット変換を行なう毎に得られるデータＬＬに対してウェーブレット変換を所望とする回数繰り返して、複数の解像度毎のデータを得る。その後、図２（ｃ）に示すように、各解像度毎のデータを符号化して原画像信号Ｓorg′ を得る。
【００６６】
なお、原画像信号Ｓorg′ において所望とする解像度までのデータのみを復号化し、かつ逆ウェーブレット変換を施すことにより、原画像の１／２ｋ （ｋ：所望とする解像度）倍の解像度且つ画像サイズが１／２２ｋ倍の画像を表す低解像度画像信号ＢＢｋ を復元することができる。この機能をなす手段として、ボケ画像信号作成手段２には、後述する第２の実施形態において説明する復号化手段３０（図１７参照）を複数段設けた構成と同じもの、および復元した所望とする解像度の低解像度画像信号ＢＢｋ をフィルタリング処理手段１０に入力する切替手段が設けられている（いずれも図１においては不図示）。なお、図１７に示すような複数の復号化手段３０を設けた構成とは異なり、データＬＬＫ＋１，ＨＬＫ，ＬＨＫおよびＨＨｋを逆ウェーブレット変換して直ちに原画像の１／２ｋ 倍の解像度の画像を表す低解像度画像信号ＢＢｋ を復元する構成とすることもできる。
【００６７】
次に、ボケ画像信号の作成処理について詳細に説明する。図３は第１の実施形態におけるボケ画像信号作成手段の構成を示す概略ブロック図である。
【００６８】
図３に示すように、第１の実施形態におけるボケ画像信号作成手段２は、まずフィルタリング処理手段１０に入力される原画像信号Ｓorg として原画像信号Ｓorg′を復元することにより得られた任意の解像度レベルの画像信号を使用し、フィルタリング処理手段１０により、この原画像信号Ｓorg に対し、原画像信号Ｓorg が表す原画像の画素のｘ方向およびｙ方向に対してフィルタリング処理を施して低解像度画像信号Ｂ１ を作成し、次にこの低解像度画像信号Ｂ１ に対して同様のフィルタリング処理を施してこの低解像度画像信号Ｂ１ よりもさらに解像度が低い低解像度画像信号Ｂ２ を作成し、以降順次同様のフィルタリング処理を重ねていくものである。そして、補間処理手段１１により、このフィルタリング処理の各段において得られる低解像度画像信号Ｂｋ に対して、それぞれ補間拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号Ｓus1 〜ＳusN （以下Ｓusk （ k＝１〜Ｎ）で代表させる）を得るものである。
【００６９】
本実施形態においては、上記フィルタリング処理のフィルタとして、１次元ガウス分布に略対応したフィルタを使用する。すなわちフィルタのフィルタ係数を、ガウス信号に関する下記の式（４）にしたがって定める。
【数１】
0004598115.tif 000002
【００７０】
これは、ガウス信号は周波数空間および実空間の双方において局在性がよいためであり、例えば上記（４）式においてσ＝１とした場合の５×１の１次元フィルタは図４に示すようなものとなる。
【００７１】
フィルタリング処理は、図５に示すように、原画像信号Ｓorg に対して、あるいは低解像度画像信号Ｂｋ に対して１画素おきに行なう。このような１画素おきのフィルタリング処理をｘ方向、ｙ方向に行なうことにより、低解像度画像信号Ｂ１ の画素数は原画像の１／４となり、フィルタリング処理により得られる低解像度画像信号に対して繰り返しこのフィルタリング処理を施すことにより、得られるｎ個の低解像度画像信号Ｂｋ （ｋ＝１〜ｎ）は、それぞれ、画素数が原画像信号Ｓorg の１／２２Ｋの画像信号となる。
【００７２】
次に、このようにして得られた低解像度画像信号Ｂｋ に対して施される補間拡大処理について説明する。補間演算の方法としては、Ｂスプラインによる方法など種々の方法が挙げられるが、本実施形態においては、補間演算についてガウス信号を用いるものとする。具体的には、下記の式（５）において、σ＝２ｋ−１ と近似したものを用いる。
【数２】
0004598115.tif 000003
【００７３】
低解像度画像信号Ｂ１ を補間する際には、ｋ＝１であるためσ＝１となる。この場合、補間処理を行なうためのフィルタは、図６に示すように５×１の１次元フィルタとなる。この補間処理は、まず低解像度画像信号Ｂ１ に対して１画素おきに値が０の画素を１つずつ補間することにより低解像度画像信号Ｂ１ を原画像と同一のサイズに拡大し、次に、この補間された低解像度画像信号Ｂ１ に対して上述した図６に示す１次元フィルタによりフィルタリング処理を施すことにより行われる。
【００７４】
同様に、この補間拡大処理を全ての低解像度画像信号Ｂｋ に対して行なう。低解像度画像信号Ｂｋ を補間する際には、上記式（５）に基づいて、３×２ｋ −１の長さのフィルタを作成し、画像信号Ｂｋ の各画素の間に値が０の画素を２ｋ −１個ずつ補間することにより、原画像と同一サイズに拡大し、この値が０の画素が補間された画像信号Ｂｋ に対して３×２ｋ −１の長さのフィルタにより、フィルタリング処理を施すことにより補間拡大する。
【００７５】
次に、上記のようにして作成されたボケ画像信号Ｓusk を用いて行われる周波数強調処理について説明する。図７は周波数強調処理を行なう装置の構成をボケ画像信号作成手段２とともに示す概略ブロック図である。図７に示すように、原画像信号Ｓorg と、フィルタリング処理手段１０および補間処理手段１１により作成されたボケ画像信号Ｓusk について、減算器２１によりそれらの信号の差分が求められ、原画像信号Ｓorg の限られた周波数帯域の成分である帯域制限画像信号（Ｓorg-Ｓus1、Ｓus1-Ｓus2など）が作成される。
【００７６】
一方、画像処理装置１に設けられているパラメータ設定手段４は、基準解像度を有する基準解像度画像について変換処理を行なう際に用いられる、該基準解像度画像についての各帯域制限画像信号に対応する基準変換関数f１〜fＮを不図示のメモリに保持している。そして、フィルタリング処理手段１０に入力される原画像信号Ｓorg として、基準解像度画像信号が入力されたときには、保持している基準変換関数f１〜fＮを各帯域制限画像信号に応じて設定する。
【００７７】
そして、フィルタリング処理手段１０に基準解像度画像信号が入力されたときには、各帯域制限画像信号が変換器２２においてそれぞれ異なる変換関数f１〜fＮにより所望の大きさとなるように変換され、さらに下記の式（２′）にしたがって、その複数の変換された帯域制限画像信号（本発明の変換画像信号）が演算器２３において積算されて周波数強調処理に必要な信号（積算信号）Ｆusm が得られ、さらに信号Ｆusm に所定の係数が乗じられた後原画像信号Ｓorg と加算されて、処理済画像信号Ｓproc′が生成される。これにより所望とする周波数成分を目的に応じた度合いで強調した処理済画像信号Ｓproc′を得ることができる。
Ｓproc′＝Ｓorg ＋β（Ｓorg ）×Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN）
Ｆusm （Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝ｆ１(Ｓorg −Ｓus1)＋ｆ２(Ｓus1 −Ｓus2)＋…
＋ｆｋ(Ｓusk-1−Ｓusk )＋…＋ｆＮ(ＳusN-1−ＳusN)…（２′）
（但し、Ｓproc′：高周波成分が強調された画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆｋ(k=1〜N）：各帯域制限画像信号を変換する変換関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
【００７８】
以上、処理済画像信号Ｓproc′の作成について説明したが、次に、第１の実施形態が解決しようとする問題点およびその解決手段について例を示して説明する。上述したように原画像と同一解像度の画像（原画像と同一の画像であっても異なる画像であってもよい）を表す処理済画像信号Ｓproc′を得る場合には、帯域制限画像信号の周波数帯域に応じた変換関数ｆｋ を用いて周波数強調処理を行えばよいが、原画像よりも低解像度の画像（原画像と同一の画像であっても異なる画像であってもよい）を表す処理済画像信号を得る場合に、その解像度に応じた変換関数を用意する必要があるため、変換関数の管理が非常に煩わしいものとなる。ここで、原画像信号Ｓorg を６つの周波数帯域に分離して周波数強調処理を行なう場合は、帯域制限画像信号の周波数特性は６つのピークを有するものとなる。
【００７９】
そして、例えば、基準解像度としての10本／ｍｍの読取密度で読み取られた基準解像度画像を表す原画像信号Ｓorg （以下10本／ｍｍの解像度を有する原画像信号ともいう；他の読取密度においても同様）から帯域制限画像信号を得る場合、原画像信号Ｓorg のナイキスト周波数は５cycle/mmであり、図８に示すように、最高周波数帯域の帯域制限画像信号のピーク周波数がこのナイキスト周波数５cycle/mmとなり、最高周波数帯域の次の周波数帯域の帯域制限画像信号のピーク周波数はナイキスト周波数の１／５の値を有する１．０cycle/mmとなり、以下低周波数帯域となるにつれて、０．５cycle/mm、０．２５cycle/mm、０．１２cycle/mm、０．０６cycle/mmのようにピーク周波数は前段のピーク周波数の１／２となる。また、基準解像度画像の１／２の解像度の画像（読取密度を５本／ｍｍとしてもよいし、基準解像度画像を間引きなどにより低解像度化してもよい）を表す画像信号を原画像信号Ｓorg としたときにおける、帯域制限画像信号のピーク周波数は、図９に示すように２．５cycle/mm、０．５cycle/mm、０．２５cycle/mm、０．１２５cycle/mmおよび０．０６cycle/mmとなる。
【００８０】
以下基準解像度画像の１／４、１／８、１／１６の解像度の画像を表す画像信号を原画像信号Ｓorg としたときにおける、帯域制限画像信号のピーク周波数は、図１０から図１２に示すようになる。つまり、帯域制限画像信号のうち最高周波数帯域の帯域制限画像信号のピーク周波数は解像度レベルに応じたナイキスト周波数となり、最高周波数帯域の次の周波数帯域の帯域制限画像信号のピーク周波数はナイキスト周波数の１／５となり、以下低周波数帯域となるにつれて、前段のピーク周波数の１／２となる。
【００８１】
次に、基準解像度画像の画像サイズとは異なる画像サイズ（以下所望画像サイズともいう）の出力画像を得る場合について説明する。
【００８２】
例えば、基準解像度画像よりも小さい所望画像サイズの画像を表す処理済画像信号を得る場合には、本発明に係る処理用画像を表す画像信号、すなわちフィルタリング処理手段１０に入力する原画像信号Ｓorg として所望画像サイズに最も近い、基準解像度画像の画像サイズを２２ｋ（ｋは負の整数）倍した変倍画像である低解像度画像を表す画像信号を使用する。以下、この画像サイズを２２ｋ倍した低解像度画像の解像度を所定解像度といい、この低解像度画像を表す画像信号を基準低解像度画像信号ともいう。
【００８３】
そして、この基準低解像度画像信号から低解像度ボケ画像信号を得、この低解像度ボケ画像信号に基づいて低解像度帯域制限画像信号を作成し、さらに低解像度帯域制限画像信号に対して帯域制限画像信号の周波数帯域に対応した変換関数により変換画像信号を得て、中間処理済画像信号Ｓproc′を得、さらに中間処理済画像信号Ｓproc′を拡大率情報Ｍに基づいて拡大または縮小して所望画像サイズの画像を表す処理済画像信号Ｓprocを得る。
【００８４】
より具体的な態様として、基準解像度画像の１／４の画像サイズで且つ１／２解像度の画像（１／２２ 縮小画像）を表す処理済画像信号Ｓprocを得る場合を図１３のフローチャートを参照して説明する。なお、ここでは原画像信号Ｓorg においては図８に示すように帯域制限画像信号は６つの周波数帯域からなるものとする。
【００８５】
まず、拡大率入力手段５からユーザが所望とする前記所望画像サイズに対応する拡大率（所望拡大率とする）を入力する（ステップＳ１）。これにより拡大率入力手段５から拡大率情報Ｍがボケ画像信号作成手段２、パラメータ設定手段４および拡大縮小手段６に入力される。ボケ画像信号作成手段２においては、拡大率情報Ｍに基づいて、所望拡大率に最も近い画像サイズが２２ｋ（ｋは負の整数）倍の画像を表す低解像度画像信号ＢＢｋ を復号化手段３０（図１７参照）により復元してフィルタリング処理手段１０に入力するとともに（ステップＳ２）、低解像度ボケ画像信号Ｓusk を作成する（ステップＳ３）。なお、低解像度画像信号ＢＢｋ が原画像信号Ｓorg に対応する基準低解像度画像信号となる。また、本実施形態においては、基準解像度画像の１／４の画像サイズに対応する拡大率情報Ｍが入力されることになるので、画像サイズが１／４倍の画像を表す低解像度画像信号ＢＢ１ を復元する。
【００８６】
ここで、フィルタリング処理手段１０に原画像信号Ｓorg として基準解像度画像信号が入力された場合に得られる６つのボケ画像信号（原画像信号含む）Ｓorg 〜Ｓus6 と、所定解像度（所望画像サイズ）に近い、１／２解像度画像（１／２２ 縮小画像）を表す低解像度画像信号ＢＢ１ から得られる低解像度ボケ画像信号Ｓorg 〜Ｓus5 との対応関係を、各ボケ画像信号の周波数範囲が同じものが同列となるように図１４に示す。なお、図１４には、１／４解像度画像（１／４２ 縮小画像）〜１／１６解像度画像（１／１６２ 縮小画像）までの原画像信号Ｓorg およびボケ画像信号Ｓusk の対応関係も示す。図１４に示すように、基準解像度画像におけるボケ画像信号Ｓus1 が、例えば１／２解像度画像における低解像度ボケ画像信号Ｓorg に対応し、以下順次図中右側に相対移動している。
【００８７】
次に、低解像度ボケ画像信号Ｓusk に基づいて低解像度帯域制限画像信号が作成される（ステップＳ４）。基準解像度画像信号から得られる６つの帯域制限画像信号と、各解像度画像を表す低解像度画像信号ＢＢｋ から得られる低解像度帯域制限画像信号との対応関係を、各帯域制限画像信号の周波数範囲が同じものが同列となるように図１５に示す。なお、フィルタリング処理手段１０に低解像度画像信号ＢＢｋ が入力された場合において生成される低解像度帯域制限画像信号のうちの０．０３cycle/mm以下の低周波数帯域のものは、基準解像度画像についての強調処理には使用されない周波数範囲の信号であるため、強調処理には使用しないようにする。図１４〜図１６における１／２〜１／１６解像度画像について、右側を途中までしか示していないのはこのためである。
【００８８】
一方、パラメータ設定手段４においては各低解像度帯域制限画像信号に対して施す変換処理の変換関数ｆｋ の設定が行われる（ステップＳ５）。ここで、基準解像度画像信号から得られる６つの帯域制限画像信号に対しては上述したように変換関数ｆ１〜ｆ６が用いられるが、１／２解像度画像を表す低解像度画像信号ＢＢ１ から得られる低解像度帯域制限画像信号に対しては、設定する変換関数に対応する低解像度帯域制限画像信号のピーク波長が、それぞれ基準解像度画像信号から得られる各帯域制限画像信号の波長と同じになるように基準変換関数ｆ１〜ｆ６をシフト（相対移動）させる。
【００８９】
各解像度画像における変換関数の対応関係を図１６に示す。例えば、図８、図９、および図１４から判るように、例えば基準解像度画像における帯域制限画像信号Ｓus1−Ｓus2のピーク波長が、１／２解像度画像における低解像度帯域制限画像信号Ｓorg−Ｓus1のピーク波長と同じになり、以下順次図１５中右側に相対移動している。したがって、図１６に示すように、処理用画像としての低解像度画像の解像度が基準解像度の１／２ｋ 倍となると、低解像度画像についての変換関数ｆＮ は基準変換関数ｆ１〜ｆ６のうちのｋ個分低周波数帯域側のものが用いられることとなる。したがって、基準解像度画像に対して周波数強調処理を施す場合の上記式（２′）における積算信号Ｆusm を式（６）に示すと、１／２解像度画像に対して周波数強調処理を施す場合の積算信号Ｆusm は式（７）に示すものとなる。
Ｆusm （Ｓorg,Ｓus1,Ｓus2,…Ｓus6）
＝ｆ１(Ｓorg −Ｓus1)＋ｆ２(Ｓus1 −Ｓus2)＋…
＋ｆ５(Ｓus4−Ｓus5 )＋ｆ６(Ｓus5−Ｓus6)…（６）
Ｆusm （Ｓorg,Ｓus1,Ｓus2,…Ｓus5）
＝ｆ２(Ｓorg −Ｓus1)＋ｆ３(Ｓus1 −Ｓus2)＋…
＋ｆ４(Ｓus3−Ｓus4 )＋ｆ５(Ｓus4−Ｓus5)…（７）
【００９０】
そして、このようにして低解像度帯域制限画像信号から低解像度変換画像信号が作成されて（ステップＳ６）低解像度積算信号Ｆusm が得られると、上記式（２′）に示すように、強調係数β（Ｓorg） が低解像度積算信号Ｆusm に乗算され、さらに原画像信号Ｓorg （ここでは基準低解像度画像信号）に加算されて中間処理済画像信号Ｓproc′が作成される（ステップＳ７）。この中間処理済画像信号Ｓproc′は、基準解像度画像の１／４の画像サイズで且つ１／２解像度の画像（１／２２ 縮小画像）を表すことになる。なお、式（２′）におけるβ（Ｓorg） は基準低解像度画像信号に基づいて設定されることになる。
【００９１】
そしてこのようにして得られた中間処理済画像信号Ｓproc′は拡大縮小手段６に入力される。拡大縮小手段６においては、拡大率入力手段５から入力された拡大率情報Ｍに基づいて、所望拡大率の画像を再現可能なように中間処理済画像信号Ｓproc′を拡大または縮小して最終的な処理済画像信号Ｓprocを得る（ステップＳ８）。ここで、本実施形態においては、拡大率入力手段５から入力された拡大率情報Ｍは１／４の画像サイズに対応するため、拡大縮小手段６においては１／２２ 縮小画像（１／２解像度画像）を表す中間処理済画像信号Ｓproc′を何ら拡大、縮小することなく最終的な処理済画像信号Ｓprocとする。
【００９２】
このように、本実施形態においては、基準解像度画像よりも画像サイズが小さい画像を表す処理済画像信号を得る、あるいは基準解像度画像よりも低解像度の画像を表す処理済画像信号を得る際に、基準解像度画像に対して周波数強調処理を行なうときに使用する基準変換関数を、基準解像度画像についての各帯域制限画像信号の周波数帯域が処理用画像についての各帯域制限画像信号の周波数帯域に最も近くなるように基準変換関数を相対移動させ、この相対移動させた基準変換関数を処理用画像についての変換関数として設定して変換画像信号を求めた後に周波数強調処理を施すようにしたので、得られる処理済画像信号を再生することにより得られる画像においては、その周波数応答特性を基準解像度画像に基づく処理済画像信号Ｓprocを再生することにより得られる画像の周波数応答特性と略一致させることができ、これにより解像度に影響されない一定の周波数応答特性を有する画像を再現可能な処理済画像信号Ｓprocを得ることができる。
【００９３】
また、基準解像度画像に対して周波数強調処理を行なう際に使用する変換関数と同一の変換関数を用いているため、解像度毎に変換関数を用意する必要が無くなり、これにより本発明を適用したシステムの構成を簡易なものとすることができ、また変換関数の管理の煩わしさをなくすことができる。
【００９４】
なお、画素数（画像サイズ）に拘わらず、基準解像度画像よりも低解像度の画像を表す処理済画像信号を得る場合においても、基準変換関数ｆ１〜ｆ６を同様にシフトさせるとよく、得られる効果も略同じである。ただし、拡大縮小手段６における処理では、鮮鋭度を向上させることができないので、最終的に得られる画像は、画像サイズが所望画像サイズになるものの、完全に所望解像度に復元できないこともある。
【００９５】
さらに、出力画像の画像サイズを基準解像度画像の画像サイズよりも小さくする場合において、フィルタリング処理手段１０に入力される画像の大きさが基準解像度画像の２２ｋ（ｋは負の整数）倍となるときには、処理用画像が基準解像度画像の画素数よりも少ない縮小画像となり、この縮小画像について強調処理などの所定の処理を施すことができるので、該所定の処理のための演算時間を短縮することができる。
【００９６】
なお、上記第１の実施形態においては、基準解像度画像信号あるいは基準低解像度画像信号からフィルタリングおよび補間拡大によりボケ画像信号を求め、基準解像度画像信号およびボケ画像信号から帯域制限画像信号を作成しているが、これに限定されるものではなく、例えば基準解像度画像信号あるいは基準低解像度画像信号をウェーブレット変換やラプラシアンピラミッドなどの手法により多重解像度に変換し、変換された各解像度毎の画像信号を帯域制限画像信号としてもよい。
【００９７】
また、上記第１の実施形態においては、多重解像度に変換されて符号化された原画像信号Ｓorg′ に基づいて所定解像度の画像を復元した後に周波数強調処理を施しているが、これに限定されるものではなく、何ら符号化されない画像信号、あるいはこの符号化されない画像信号とは異なる解像度の画像信号に対して周波数強調処理を施してもよい。
【００９８】
次いで、本発明の第２の実施形態について説明する。上記第１の実施形態においては、原画像信号Ｓorg を多重解像度空間に変換して符号化した原画像信号Ｓorg′を所望画像サイズに最も近い、１／２ｋ 解像度の画像（１／２２ｋ縮小画像）を表す低解像度画像信号ＢＢｋ まで復元し、この復元された低解像度画像信号ＢＢｋ をフィルタリングおよび補間拡大を順次繰り返して複数のボケ画像信号を作成していたが、第２の実施形態においては、最低解像度（最小画像サイズ）から所望画像サイズに最も近い、１／２ｋ 解像度の画像（１／２２ｋ縮小画像）までの各解像度レベル（各画像サイズ）の低解像度画像信号ＢＢｋ をそれぞれ補間拡大して複数のボケ画像信号を作成するものである。なお、第２の実施形態においては、ボケ画像信号作成手段２における処理が異なるのみであって、他の構成については第１の実施形態と同様であるため、詳細な説明は省略する。
【００９９】
ここで、第２の実施形態におけるボケ画像信号の作成処理について説明する。図１７は第２の実施形態におけるボケ画像信号作成手段２の構成を示す概略ブロック図である。図１７に示すように、まず復号化手段３０により原画像信号Ｓorg′ における最低解像度の画像を表すデータを復号化および逆ウェーブレット変換して原画像の１／２ｎ 倍の解像度の画像を表す低解像度画像信号ＢＢｎ を作成し、次にこの低解像度画像信号ＢＢｎ を用いて原画像の１／２ｎ−１ 倍の解像度の画像を表す低解像度画像信号ＢＢｎ−１ を作成し、以降順次復号化および逆ウェーブレット変換を行なって原画像の１／２ｋ （ｋ＝１〜ｎ）倍の解像度を表す低解像度画像信号ＢＢｋ および最高解像度の画像すなわち原画像（基準解像度画像）を表す原画像信号Ｓorg を作成する。なお、低解像度画像信号ＢＢｋ は原画像信号Ｓorg をウェーブレット変換することにより得られるデータＬＬｋに対応するものとなる。そして、補間処理手段３１により、この復号化処理の各段において得られる低解像度画像信号ＢＢｋ に対して、それぞれ上記第１の実施形態と同様の補間拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号Ｓus1 〜ＳusN （以下Ｓusk（k＝１〜Ｎ）で代表させる）を得る。
【０１００】
次に、このボケ画像信号Ｓusk から帯域制限画像信号を得る。そして、上記式（２′）に示すように周波数強調処理を施して処理済画像信号Ｓproc′を得る。
【０１０１】
ここで、第２の実施形態においても、第１の実施形態と同様にして基準解像度画像よりも小さい所望画像サイズの画像（基準解像度よりも低い所定解像度の画像と同じ意味）を表す処理済画像信号Ｓprocを得る場合について説明する。なお、第２の実施形態においても第１の実施形態と同様に生成される帯域制限画像信号の周波数帯域は６つであり、所望画像サイズが基準解像度画像の１／４サイズであるとする。
【０１０２】
まず、拡大率入力手段５からユーザが所望拡大率を入力する。これにより拡大率入力手段５から拡大率情報Ｍがボケ画像信号作成手段２、パラメータ設定手段４および拡大縮小手段６に入力される。ボケ画像信号作成手段２においては、拡大率情報Ｍに基づいて、最低解像度から所望拡大率に最も近い画像サイズが２２ｋ（ｋは負の整数）倍までの各画像サイズ（各解像度）の画像を表す低解像度画像信号ＢＢｋ を復元するとともに、低解像度ボケ画像信号Ｓusk を作成する。なお、低解像度画像信号ＢＢｋ が原画像信号Ｓorg に対応する基準低解像度画像信号となる。また、本実施形態においては、基準解像度画像の１／４の画像サイズに対応する拡大率情報Ｍが入力されることになるので、画像サイズが１／４倍の画像を表す低解像度画像信号ＢＢｎ 〜ＢＢ１ を復元する。
【０１０３】
ここで、生成される６つのボケ画像信号（原画像信号含む）Ｓorg 〜Ｓus6 と、所定解像度（所望画像サイズ）に最も近い、１／２解像度画像（１／２２ 縮小画像）を表す低解像度画像信号ＢＢ１ から得られる低解像度ボケ画像信号Ｓorg 〜Ｓus5 との対応関係、および基準解像度画像信号から得られる帯域制限画像信号と各解像度画像を表す低解像度画像信号ＢＢｋ から得られる帯域制限画像信号の関係は、上記図１４および図１５に示すものとなる。また、パラメータ設定手段４においては各低解像度帯域制限画像信号に対して施す変換処理の変換関数の設定が第１の実施形態と同様に行われ、解像度が１／２ｋ となると変換関数ｆＮ は、上記式（６）および式（７）に示すように、基準変換関数ｆ１ 〜ｆ６ のうちのｋ個分低周波数帯域側のものが用いられることとなる。
【０１０４】
そして、このようにして低解像度帯域制限画像信号から低解像度変換画像信号が作成されて低解像度積算信号が得られると、強調係数β（Ｓorg） が積算信号に乗算され、さらに所定解像度の画像に対応する原画像信号Ｓorg （基準低解像度画像信号）に加算されて中間処理済画像信号Ｓproc′が得られる。
【０１０５】
そしてこのようにして得られた中間処理済画像信号Ｓproc′は拡大縮小手段６に入力される。拡大縮小手段６においては、拡大率入力手段５から入力された拡大率情報Ｍに基づいて、所望拡大率の画像を再現可能なように中間処理済画像信号Ｓproc′を拡大または縮小して最終的な処理済画像信号Ｓprocを得る。ここで、第２の実施形態においても、拡大率入力手段５から入力された拡大率情報Ｍは１／４の画像サイズに対応するため、拡大縮小手段６においては１／２２ 縮小画像（１／２解像度画像）を表す中間処理済画像信号Ｓproc′を何ら拡大、縮小することなく最終的な処理済画像信号Ｓprocとする。
【０１０６】
このように、第２の実施形態においても、基準解像度画像よりも画像サイズが小さい画像を表す処理済画像信号を得る際に、基準解像度画像に対して周波数強調処理を行なう際に使用する基準変換関数を所望画像サイズ（所定解像度）に応じてシフト（相対移動）して帯域制限画像信号の変換を行なうようにしたため、得られる処理済画像信号を再生することにより得られる画像においては、その周波数応答特性を基準解像度画像に基づく処理済画像信号Ｓprocを再生することにより得られる画像の周波数応答特性と略一致させることができ、これにより解像度に影響されない一定の周波数応答特性を有する画像を再現可能な処理済画像信号Ｓprocを得ることができる。
【０１０７】
なお、上記第１および第２の実施形態においては、１／２縮小画像を得る場合について説明しているが、例えば１／３２ 縮小画像（１／３解像度画像）を得る場合には、ボケ画像信号作成手段２および周波数強調処理手段３において、１／３２ 縮小画像に最も近い１／４２ 縮小画像（１／４解像度画像）を表す低解像度画像信号に基づいて周波数強調処理が行われて中間処理済画像信号Ｓproc′が得られる。そして、この中間処理済画像信号Ｓproc′を拡大縮小手段６において、（４／３）２ 倍に拡大して最終的な処理済画像信号Ｓprocが得られることとなる。なお、所望画像サイズが基準解像度画像の１／２２ｋ倍の場合には、拡大縮小手段６を省略することができるため、装置の構成を簡易なものとすることができる。
【０１０８】
また、上記第１および第２の実施形態においては、式（２′）に基づいて帯域制限画像信号を変換関数ｆｋ により変換して変換画像信号を作成しているが、下記の式（８）に示すように、変換関数ｆｋ に代えて強調係数αｋ を用いて帯域制限画像信号を強調して変換画像信号を得てもよい。この場合、基準解像度画像用に設定される強調係数α１〜α６を、上記変換関数ｆｋ をシフトさせたと同じように所望画像サイズ（所定解像度）に応じてシフトして強調処理が行われることとなる。
Ｓproc＝Ｓorg＋β（Ｓorg）×Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝α１(Ｓorg−Ｓus1)＋α２(Ｓus1−Ｓus2)＋…
＋αｋ(Ｓusk-1−Ｓusk)＋…＋αＮ(ＳusN-1−ＳusN)  …（８）
【０１０９】
上記各実施形態は所望画像サイズ（所定解像度）が基準解像度画像よりも小さい場合についてのものであるが、逆に、所望画像サイズ（所定解像度）が基準解像度画像よりも大きい場合についても、上記各実施形態と略同様の方法を用いることができる。以下、この場合について説明する。
【０１１０】
医療診断などにおいて、ＣＲＴ上で画像を見て画像の一部分を詳細に読影する場合には、画像の該当する部分を拡大して表示したりあるいはフィルム出力することがある。この場合、従来の態様では、拡大処理した画像信号に対して基準解像度画像に対する変換関数を用いて処理して強調処理を施すことになる。一方、拡大処理に際しては、補間処理を施して２２ｎ倍拡大画像（ｎは正の整数）を得る態様と、多重解像度信号として、予め拡大画像用に用意された２ｎ 倍解像度レベルの信号を用いる態様とが考えられる。ここで、拡大画像を表す画像信号の周波数帯域は、前者の場合には画素数は２２ｎ倍になるものの最大でも基準解像度画像の周波数帯域と同じであり、後者の場合には２ｎ 倍解像度レベルの信号を用いている分だけ拡がる。
【０１１１】
これに対して、帯域制限画像信号の生成過程、変換処理、および強調処理としては、入力される画像信号の態様に拘わらず同じ処理がなされ、拡大処理に際して補間処理を施すか２ｎ 倍解像度レベルの信号を用いるかに拘わらず、画像サイズと帯域制限画像信号との間には、図１５から推測されるように、処理対象となる帯域制限画像信号と変換関数の組数が増えるため基準解像度に存在しない変換関数を定義する必要が生じる。
【０１１２】
一方、拡大画像における各帯域制限画像信号が担持する周波数成分は、拡大処理に際して２ｎ 倍解像度レベルの信号を用いた場合には、図１８に示す同列部分が同じ周波数範囲となり、各帯域制限画像信号のピーク波長は、例えば、基準解像度画像におけるＳorg−Ｓus1が、２倍解像度画像（２２ 倍拡大画像）におけるＳus1 −Ｓus2 、４倍解像度画像（４２ 倍拡大画像）におけるＳus2 −Ｓus3 と同じになる。そして、２倍解像度画像におけるＳorg−Ｓus1、４倍解像度画像におけるＳorg −Ｓus1 およびＳus1 −Ｓus2 が、基準解像度画像に存在しない高周波成分（高周波帯域制限画像信号）となる。そこで、拡大画像用に用意された２ｎ 倍解像度レベルの信号を用いて強調処理を行なうに際しては、図１９に示すように、基準解像度画像と同じ周波数成分に対しては基準解像度画像についてのものと同じ変換関数が用いられるようにすると共に、基準解像度画像に存在しない高周波成分に対応する変換関数として、基準解像度画像の最高周波数帯域で用いた変換関数ｆ１ をそのまま利用することとする。これにより、拡大画像を担持する処理済画像信号Ｓprocの応答を、視覚的に影響が大きな低周波側から基準解像度画像の最高周波数（ナイキスト周波数）までに亘って、基準解像度画像と略同じにすることができると共に、基準解像度画像に存在しない高周波成分についても、変換関数ｆ１ を用いることで、基準解像度画像についてのものとある程度の対応関係をもつ強調処理を施すことができる。
【０１１３】
他方、補間処理を施して２２ｎ倍拡大画像を得、この拡大画像を表す画像信号を用いて強調処理を行なう場合には、拡大画像における各帯域制限画像信号のピーク波長は、例えば基準解像度画像におけるＳorg−Ｓus1が２２ 倍拡大画像におけるＳus1 −Ｓus2 、４２ 倍拡大画像におけるＳus2 −Ｓus3 と同じになる。そして、２２ 倍拡大画像におけるＳorg−Ｓus1、４２ 倍拡大画像におけるＳorg −Ｓus1 およびＳus1 −Ｓus2 が基準解像度画像の最高周波数（ナイキスト周波数）と同じになる。そこで、補間処理を施して得た２２ｎ倍拡大画像を用いて強調処理などを行う場合には、基準解像度画像と同じ周波数成分に対しては基準解像度画像と同じ変換関数が用いられるように、図１９と同様の変換関数を用いる。つまり、２２ 倍拡大画像におけるＳus1 −Ｓus2 、４２ 倍拡大画像におけるＳus2 −Ｓus3 、および２２ 倍拡大画像におけるＳorg−Ｓus1、４２ 倍拡大画像におけるＳorg −Ｓus1 、Ｓus1 −Ｓus2 に対して変換関数ｆ１ を用いる。これにより、拡大画像を担持する処理済画像信号Ｓprocの応答を、視覚的に影響が大きな低周波側から基準解像度画像の最高周波数まで（本例では略全帯域となる）に亘って、基準解像度画像と略同じにすることができる。
【０１１４】
次に、本発明の第３の実施形態について説明する。図２０は本発明の第３の実施形態による画像処理装置の構成を示す概略ブロック図である。この図２０において、第１の実施形態の要素と同等の要素には同番号を付し、それらについての説明は特に必要のない限り省略する。
【０１１５】
第３の実施形態に係る画像処理装置１は、読取装置などにおいて得られた任意の解像度を有する原画像信号Ｓorg からボケ画像信号を作成するボケ画像信号作成手段２と、上述した式（２）に基づいて特定の周波数を強調するための周波数強調処理を行って処理済画像信号Ｓprocを得る、第１の実施形態にかかる周波数強調処理手段３と同じ機能をなす変換処理手段３とを有する。さらに、この第３の実施形態にかかる画像処理装置１は、変換関数算出手段４０と解像度情報入力手段５０とを備える。変換関数算出手段４０は、変換処理手段３が変換処理に使用する変換関数を算出する手段であり、後述するように解像度情報入力手段５０から入力される解像度情報に基づいて、基準解像度を有する基準解像度画像信号（以下基準原画像信号とする）Ｓorg０ に対して周波数強調処理を施す際に使用する変換関数を用いて、入力された画像信号が表す画像の解像度に応じた変換関数を算出するものである。解像度情報入力手段５０は原画像信号Ｓorg の解像度情報Ｍを得るための手段である。解像度情報入力手段５０による入力は、ユーザがキーボードから数値として入力してもよいし、操作画面に数種類の密度を表示してユーザに選択させるといった形態でもよい。あるいは読取装置側で、解像度情報Ｍを原画像信号Ｓorgに添付しておき、画像処理装置１が入力された画像信号毎に添付された解像度情報Ｍを認識するようにしてもよく、変換関数算出手段４０が解像度を認識することができれば、どのような形態であってもよい。
【０１１６】
なお、第３の実施形態の変換器２２において用いられる変換関数f１〜fＮは定数とする。ここで、第１の実施形態の図８に示すように、10本／ｍｍの読取密度で読み取られた原画像信号Ｓorg から得られる各帯域制限画像信号に対して設定される変換関数は下記の表１に示すような値を有するものとする。
【表１】
0004598115.tif 000004
【０１１７】
次に、第３の実施形態が解決しようとする問題点およびその解決手段について例を示して説明する。上述したように原画像信号Ｓorgの所望とする周波数成分を強調するためには、帯域制限画像信号の周波数帯域に応じた変換関数ｆｋを用いて周波数強調処理を行えばよい。しかしながら、10本／ｍｍの解像度を有する原画像信号Ｓorgに対して周波数強調処理を行なう際に使用する変換関数を、３００ｄｐｉの解像度を有する原画像信号Ｓorgに対して適用した場合、10本／ｍｍの解像度を有する原画像信号Ｓorgから得られる処理済画像信号Ｓprocと３００ｄｐｉの解像度を有する原画像信号Ｓorgから得られる処理済画像信号Ｓprocとにおいて、得られる画像の周波数応答特性が図２１に示すように異なるものとなってしまう。
【０１１８】
したがって、第３の実施形態においては、例えば10本／ｍｍの解像度を有する原画像信号Ｓorgを基準原画像信号Ｓorg０とし、基準原画像信号Ｓorg０に対して周波数強調処理を施す際に使用する基準変換関数を用いて、３００ｄｐｉの原画像信号Ｓorgに対して周波数強調処理を施す際に使用する変換関数を算出するようにしたものである。以下、変換関数の算出について説明する。
【０１１９】
図２２は３００ｄｐｉ解像度の原画像信号Ｓorgから得られる帯域制限画像信号の周波数応答特性を示す図である。図２２に示すように、３００ｄｐｉ解像度の原画像信号Ｓorgのナイキスト周波数は５．９cycle/mmであり、最高周波数帯域の帯域制限画像信号のピーク周波数がこのナイキスト周波数となる。そして、最高周波数帯域の次の周波数帯域の帯域制限画像信号のピーク周波数はナイキスト周波数の１／５の値である１．１８cycle/mmとなり、以下低周波数帯域となるにつれて、０．５９cycle/mm、０．３０cycle/mm、０．１５cycle/mm、０．０７cycle/mmのようにピーク周波数は前段の周波数帯域の１／２の関係となっている。この関係は図７に示した基準原画像信号Ｓorg０の帯域制限画像信号の周波数応答特性の関係と一致するものであり、この関係は、低解像度画像信号を作成する際のフィルタを変更しない限り原画像信号Ｓorgの解像度に拘わらず成立するものである。
【０１２０】
したがって、基準原画像信号Ｓorg０から得られる帯域制限画像信号のピーク周波数、３００ｄｐｉ解像度の原画像信号Ｓorgから得られる帯域制限画像信号のピーク周波数および基準変換関数に基づいて、内挿あるいは外挿の線形補間を行なうことにより、下記の表２に示すように３００ｄｐｉ画像に対する変換関数を求めることができる。具体的には、ピーク周波数が１．１８cycle/mmとなる帯域制限画像信号に対応する変換関数は、表１における５cycle/mmおよび１cycle/mmの変換関数の値から、ピーク周波数と変換関数の値との関係を下記の式（９）により求め、この式（９）に３００ｄｐｉ解像度のピーク周波数の値を代入してそのピーク周波数における変換関数の値を求める。そして、この線形補間を全ピーク周波数に対して行ない、各帯域制限画像信号に対応する変換関数を求める。なお、ここでは線形補間により求めているが、これに限定されるものではない。
【表２】
0004598115.tif 000005
【０１２１】
ｆｘ＝１／４０×Ｃ＋０．８７５ （９）
但し、ｆｘ：求めるべき変換関数の値
Ｃ：ピーク周波数
【０１２２】
そしてこのようにして求められた変換関数により変換処理手段３において、変換処理（周波数強調処理）が行われて処理済画像信号Ｓprocが得られる。ここで、上述したようにして求められた変換関数により処理が施された処理済画像信号Ｓprocと基準原画像信号Ｓorg０に対して基準変換関数により処理を施すことにより得られた処理済画像信号Ｓproc０とにおける周波数応答特性の関係を図２３に示す。図２３に示すように、第１の実施形態によれば処理を施す原画像信号の解像度に拘わらず略同一の周波数応答特性を有する処理済画像信号Ｓprocが得られることが分かる。
【０１２３】
次いで、第３の実施形態の動作について説明する。図２４は本実施形態の動作を示すフローチャートである。まず、読取装置などから原画像信号Ｓorgが画像処理装置１に入力される（ステップＳ１）。原画像信号Ｓorgはボケ画像信号作成手段２に入力されてここでボケ画像信号Ｓuskが作成される（ステップＳ２）。一方、解像度情報入力手段５０には、原画像信号Ｓorgの解像度情報Ｍが入力され（ステップＳ３）、入力された解像度情報Ｍに基づいて変換関数算出手段４０において上述したように解像度に応じた変換関数が算出される（ステップＳ４）。なお、ステップＳ１，２の処理よりもステップＳ３，４の処理を先に行ってもよいが、これらの処理を並列して行なうことにより演算時間を短縮することができる。変換処理手段３においては、ボケ画像信号Ｓuskに基づいて帯域制限画像信号が作成され、さらに変換関数算出手段４０において算出された変換関数に基づいて上記式（２）に示す周波数強調処理が行われ処理済画像信号Ｓprocが得られる（ステップＳ５）。
【０１２４】
この周波数強調処理は、ボケ画像信号を用いて強調のために原画像信号に加算する加算信号を作成するものであり、周波数強調処理を行なうことによってアーチファクトが発生しないようにするために、この加算信号を構成する各周波数帯域の信号がそれぞれ所望の信号となるよう、周波数帯域毎に異なる変換処理を施すものである。所望の信号を作成するためには、原画像信号Ｓorgの解像度に拘わらず所望の周波数応答特性を有する信号となることが好ましい。第３の実施形態によれば、解像度に拘わらず処理済画像信号Ｓprocの周波数応答特性が略同一のものとなるように変換関数を算出するようにしたため、原画像の解像度に拘わらず略同一の周波数応答特性を有する処理済画像信号Ｓprocを得ることができる。
【０１２５】
また、基準変換関数に基づいて解像度毎に変換関数を算出しているため、解像度毎に複数の変換関数を用意する必要が無くなり、これにより変換関数の管理の煩わしさをなくすことができる。
【０１２６】
なお、上記第３の実施形態においては、基準変換関数を10本／ｍｍの解像度に対応する１組のみ用いて任意の解像度に対応する変換関数を算出しているが、例えば利用頻度が高い数種類の基準解像度に対応する変換関数を基準変換関数として用意し、処理を施す画像の解像度が上記数種類の基準解像度のいずれに近い値を有するか否かを判断し、解像度が近い方の基準変換関数を用いて処理を施す画像の解像度に対応した変換関数を求めてもよい。例えば、２種類の基準解像度およびこれらに対応する変換関数を第１および第２の基準解像度、第１および第２の基準変換関数とし、さらに第１および第２の基準変換関数の使用範囲を図２５（ａ）に示すように設定した場合において、処理を施す原画像信号により表される画像の解像度が図中矢印Ａの部分である場合には第１の基準変換関数を使用し、矢印Ｂの部分である場合には第２の基準変換関数を使用する。また、第１および第２の基準変換関数の使用範囲を図２５（ｂ）に示すように設定した場合において、処理を施す原画像信号により表される画像の解像度が図中矢印Ｃの部分である場合には第１の基準変換関数を使用し、矢印Ｄの部分である場合には第２の基準変換関数を使用してもよい。
【０１２７】
また、上記第３の実施形態においては変換関数を定数としているが、所定の傾きを有する非線形関数を用いてもよい。以下、変換関数を非線形関数とした場合における解像度に応じた変換関数の算出について第４の実施形態として説明する。第４の実施形態においては、変換関数ｆｋとして下記の式（１０）に示す非線形関数が用いられる。
f(Sin)=Ｓout＝Sin×Ｙ×（exp(Ｘ/Sin)-1)/(exp(Ｘ/Sin)+1) …（１０）
【０１２８】
式（１０）において、Ｓinは入力信号、Ｓoutは出力信号、Ｘは非線形の程度すなわち抑制条件を決定するパラメータ、Ｙは関数全体の傾きすなわち周波数応答特性をコントロールするパラメータである。そして、パラメータＸ，Ｙを調整することにより帯域制限画像信号の周波数応答特性を変更することができる。第４の実施形態においては、パラメータＸ，Ｙとして、下記の表３に示すものを用いる。なお、表３においては、10本／ｍｍ、6.7本／ｍｍおよび５本／ｍｍの画素密度を有する原画像信号に対するパラメータを示すものであり、上の行ほど高周波帯域の帯域制限画像信号に対するパラメータとなっている。ここで、ナイキスト周波数が低くなるほど非線形処理に使用する帯域制限画像信号の数は図２６に示すように少なくなる。例えば10本／ｍｍの解像度の場合には６つの帯域制限画像信号を用いるが6.7本／ｍｍおよび５本／ｍｍの解像度の場合には５つの帯域制限画像信号のみが用いられる。これは、6.7本／ｍｍおよび５本／ｍｍの画素密度において０．０３cycle/ｍｍにピークを有する低周波数帯域に対しては、10本／ｍｍに相当する低周波数帯域が存在しないためである。また、ここでも基準原画像信号の解像度を10本／ｍｍとし、表３に示すパラメータＸ，Ｙの組からなる変換関数を基準変換関数とする。また、処理が施される原画像信号により表される画像の解像度を２００ｄｐｉとする。
【表３】
0004598115.tif 000006
【０１２９】
図２７は２００ｄｐｉ解像度の原画像信号Ｓorgから得られる帯域制限画像信号の周波数応答特性を概略的に示す図である。図２７に示すように、２００ｄｐｉ解像度の原画像信号Ｓorgのナイキスト周波数は約４．０（厳密には３．９３７）cycle/mmであり、最高周波数帯域の帯域制限画像信号のピーク周波数がこのナイキスト周波数となる。そして、最高周波数帯域の次の周波数帯域の帯域制限画像信号のピーク周波数はナイキスト周波数の１／５の値である０．８cycle/mmとなり、以下低周波数帯域となるにつれて、０．４cycle/mm、０．２cycle/mm、０．１cycle/mm、０．０５cycle/mmのようにピーク周波数は前段の周波数帯域の１／２の関係となっている。この関係は図８に示した10本／ｍｍの解像度を有する基準原画像信号Ｓorg０の帯域制限画像信号の周波数応答特性の関係と一致するものであり、この関係は、低解像度画像信号を作成する際のフィルタを変更しない限り原画像信号Ｓorgの解像度に拘わらず成立するものである。このような周波数応答特性を有する２００ｄｐｉ画像を表す原画像信号Ｓorgに対して基準変換関数により処理を施すことにより得られる処理済画像信号Ｓprocと基準原画像信号Ｓorg０に対して基準変換関数により処理を施すことにより得られる処理済画像信号Ｓproc０との周波数応答特性を図２８に示す。図２８に示すように処理済画像信号Ｓprocは処理済画像信号Ｓproc０と比較して低周波数帯域が強調された画像となってしまう。
【０１３０】
このため、上記第３の実施形態と同様にして、処理を施す原画像信号の解像度に応じて基準変換関数から処理に用いる変換関数を算出する。以下変換関数の算出について説明する。第４の実施形態においては、基準原画像信号から得られる帯域制限画像信号のピーク周波数、変換関数のパラメータ（ここでは表３におけるＹの値）および処理を施す２００ｄｐｉ画像を表す原画像信号Ｓorgから得られる帯域制限画像信号のピーク周波数が予め分かっており、さらには求める変換関数の数が６個であることも分かっている。そこで、まず下記の式（１２）により基準原画像信号の帯域制限画像信号の周波数応答特性、２００ｄｐｉ解像度の原画像信号に基づく帯域制限画像信号のピーク周波数、および基準変換関数のパラメータＹから２００ｄｐｉ解像度の原画像信号Ｓorgに基づく帯域制限画像信号のピーク周波数に対応する周波数でのレスポンスを求める。
R1=X1[F1]*Y1＋X2[F1]*Y2＋X3[F1]*Y3＋X4[F1]*Y4＋X5[F1]*Y5＋X6[F1]*Y6
R2=X1[F2]*Y1＋X2[F2]*Y2＋X3[F2]*Y3＋X4[F2]*Y4＋X5[F2]*Y5＋X6[F2]*Y6
R3=X1[F3]*Y1＋X2[F3]*Y2＋X3[F3]*Y3＋X4[F3]*Y4＋X5[F3]*Y5＋X6[F3]*Y6
R4=X1[F4]*Y1＋X2[F4]*Y2＋X3[F4]*Y3＋X4[F4]*Y4＋X5[F4]*Y5＋X6[F4]*Y6
R5=X1[F5]*Y1＋X2[F5]*Y2＋X3[F5]*Y3＋X4[F5]*Y4＋X5[F5]*Y5＋X6[F5]*Y6
R6=X1[F6]*Y1＋X2[F6]*Y2＋X3[F6]*Y3＋X4[F6]*Y4＋X5[F6]*Y5＋X6[F6]*Y6
…（１２）
但し、R1〜R6 ： ピーク周波数におけるレスポンス（２００ｄｐｉ画像）
X1[F1]〜X6[F6]：基準原画像信号に基づいて作成された帯域制限画像信号のF1〜F6の周波数におけるレスポンス
Ｙ1〜Ｙ6 ： 基準変換関数のパラメータ
【０１３１】
そしてこのようにしてレスポンスＲ１〜Ｒ２が求められると下記の連立方程式（１３）を解くことにより、２００ｄｐｉ解像度画像に対応する変換関数のパラメータＡ１〜Ａ６を求める。
R1=Z1[F1]*A1＋Z2[F1]*A2＋Z3[F1]*A3＋Z4[F1]*A4＋Z5[F1]*A5＋Z6[F1]*A6
R2=Z1[F2]*A1＋Z2[F2]*A2＋Z3[F2]*A3＋Z4[F2]*A4＋Z5[F2]*A5＋Z6[F2]*A6
R3=Z1[F3]*A1＋Z2[F3]*A2＋Z3[F3]*A3＋Z4[F3]*A4＋Z5[F3]*A5＋Z6[F3]*A6
R4=Z1[F4]*A1＋Z2[F4]*A2＋Z3[F4]*A3＋Z4[F4]*A4＋Z5[F4]*A5＋Z6[F4]*A6
R5=Z1[F5]*A1＋Z2[F5]*A2＋Z3[F5]*A3＋Z4[F5]*A4＋Z5[F5]*A5＋Z6[F5]*A6
R6=Z1[F6]*A1＋Z2[F6]*A2＋Z3[F6]*A3＋Z4[F6]*A4＋Z5[F6]*A5＋Z6[F6]*A6
…（１３）
但し、R1〜R6 ： 上記式（１２）により求めたレスポンス
Z1[F1]〜Z6[F6]：２００ｄｐｉ解像度の原画像信号に基づいて作成された帯域制限画像信号のF1〜F6の周波数におけるレスポンス
Ａ1〜Ａ6 ： 変換関数のパラメータ
【０１３２】
このようにして求められた変換関数を用いて上記第３の実施形態と同様にして２００ｄｐｉ解像度の画像信号に対して周波数強調処理を施すことにより、解像度に拘わらず処理済画像信号Ｓprocの周波数応答特性が略同一のものとすることができ、これにより原画像の解像度に拘わらず略同一の周波数応答特性を有する処理済画像信号Ｓprocを得ることができる。
【０１３３】
ここで、上記式（１３）におけるZ1[F1]〜Z6[F6]の算出方法について説明する。10本／ｍｍの解像度を有する基準原画像信号Ｓorg０の帯域制限画像信号は上記図８に示す特性を有するものとなる。図８に示す特性において実線で示す最高周波数帯域の特性は、図４に示す１次元フィルタのフィルタ係数をフーリエ変換した結果と、図６に示す１次元フィルタのフィルタ係数をフーリエ変換した結果とを乗算することによりボケ画像信号の周波数応答特性を得、さらにこの周波数応答特性を基準原画像信号Ｓorg０の周波数応答特性（全周波数帯域に亘って１の値を有する）から減算することにより求められる。さらに低周波数帯域の特性については、その周波数帯域に対応する低解像度画像信号を求めるフィルタのフィルタ係数およびボケ画像信号を求めるフィルタ係数に基づいて上記と同様にしてボケ画像信号の周波数応答特性を得、これを１つ高周波数帯域のボケ画像信号の周波数応答特性から減算することにより求められる。
【０１３４】
ここで、基準原画像信号Ｓorg０の最高周波数帯域の帯域制限画像信号のレスポンスをX1[Fi]（Fi：周波数）なる関数で表し、周波数を2048によりサンプリングした場合の帯域制限画像信号のレスポンスをX1[i]なる関数とすると、Fiとiとの関係は、
ｉ＝2047×Fi／fnq （１４）
但し、fnq：ナイキスト周波数
０≦ｉ≦2047
となる。また、この対応関係を図２９に示す。なお、ｉは関数X1[i]をテーブルとして扱った場合におけるアドレスとなる。このようにして帯域制限画像信号の全周波数帯域におけるレスポンスを関数X1[i]〜X6[i]として求め、これらの関数X1[i]〜X6[i]をテーブルとして扱うことにより、上記式（１２）および式（１３）におけるX1[F1]〜X6[F6]の値を式（１４）を用いて容易に求めることができる。
【０１３５】
このようにして求められた全周波数帯域のレスポンスの関数X1[i]〜X6[i]に基づいてレスポンスZ1[F1]〜Z6[F6]を求める。なお、ここでは関数X1[i]から任意の周波数FiにおけるレスポンスZ1[Fi]を求める場合について説明する。レスポンスZ1[Fi]は２００ｄｐｉ解像度の原画像信号に基づいて作成された帯域制限画像信号のうち最高周波数帯域の帯域制限画像信号の周波数Fiにおけるレスポンスである。２００ｄｐｉ解像度の原画像信号のナイキスト周波数は４cycle/mmであり、レスポンスZ1[Fi]の周波数を2048によりサンプリングした場合のレスポンスをZ1[i]なる関数で表すものとすると、周波数Fiとiとの関係はX1[i]と同様に式（１４）に示すものとなる。関数X1[i]と関数Z1[i]の関係を図３０に示す。概念的には関数X1[i]と関数Z1[i]とはナイキスト周波数が異なる別の関数であるが、関数X1[i]を用いてレスポンスZ1[Fi]を近似的に求めることができる。例えば、２cycle/mmのレスポンスZ1[Fi]を求める場合は、上記式（１４）により、
ｉ＝2047×２／４＝1024
となるため、関数X1[i]におけるアドレスが1024となる所の値を求め、それを２cycle/mmにおけるレスポンスZ1[Fi]として求める。
【０１３６】
このように、関数X1[i]〜X6[i]をテーブルとして持ち、周波数強調処理を行なう原画像信号のナイキスト周波数および上記式（１４）に基づいてアドレスｉを求めることにより、関数X1[i]〜X6[i]のテーブルを参照して近似的にレスポンスZ1[F1]〜Z6[F6]を求めることができる。
【０１３７】
なお、３００ｄｐｉ解像度の画像を表す原画像信号Ｓorgに対して周波数強調処理を施す場合には、上記式（１２）により基準原画像信号の帯域制限画像信号の周波数応答特性、３００ｄｐｉ解像度画像信号に基づく帯域制限画像信号のピーク周波数、および基準変換関数のパラメータＹから３００ｄｐｉ解像度の原画像信号に基づく帯域制限画像信号のピーク周波数に対応する周波数でのレスポンスを求める。そして、上記式（１３）により変換関数のパラメータを求めることができる。なお、３００ｄｐｉ解像度においてはナイキスト周波数は５．９cycle/mmであって基準原画像信号Ｓorg０のナイキスト周波数５cycle/mmを越えるものとなる。ここで、上記式（１２）におけるX1[F1]〜X6[F6]は、０〜５cycle/mmの範囲に対する関数値と考えているため、５．９cycle/mmに対応する値がなく、式（１２）におけるR1を求めることができない。したがって、３００ｄｐｉ解像度の画像においては、ピーク周波数として基準原画像信号から得られる帯域制限画像信号のピーク周波数を用いる。
【０１３８】
また、１００ｄｐｉ解像度の画像を表す原画像信号Ｓorgに対して周波数強調処理を施す場合について下記に説明する。１００ｄｐｉ解像度のナイキスト周波数は約２．０cycle/mm（厳密には１．９７cycle/mm）であるため、帯域制限画像信号の数は５つとなる。したがって、上記式（１２）、（１３）は下記の式（１５）、（１６）に示すように変更され、これにより変換関数の５つのパラメータが求められることとなる。
R1=X1[F1]*Y1＋X2[F1]*Y2＋X3[F1]*Y3＋X4[F1]*Y4＋X5[F1]*Y5＋X6[F1]*Y6
R2=X1[F2]*Y1＋X2[F2]*Y2＋X3[F2]*Y3＋X4[F2]*Y4＋X5[F2]*Y5＋X6[F2]*Y6
R3=X1[F3]*Y1＋X2[F3]*Y2＋X3[F3]*Y3＋X4[F3]*Y4＋X5[F3]*Y5＋X6[F3]*Y6
R4=X1[F4]*Y1＋X2[F4]*Y2＋X3[F4]*Y3＋X4[F4]*Y4＋X5[F4]*Y5＋X6[F4]*Y6
R5=X1[F5]*Y1＋X2[F5]*Y2＋X3[F5]*Y3＋X4[F5]*Y4＋X5[F5]*Y5＋X6[F5]*Y6
…（１５）
但し、R1〜R5 ： ピーク周波数におけるレスポンス（１００ｄｐｉ画像）
X1[F1]〜X6[F5]：原画像信号に基づいて作成された帯域制限画像信号のF1〜F5の周波数におけるレスポンス
Ｙ1〜Ｙ6 ： 基準変換関数のパラメータ
R1=Z1[F1]*A1＋Z2[F1]*A2＋Z3[F1]*A3＋Z4[F1]*A4＋Z5[F1]*A5
R2=Z1[F2]*A1＋Z2[F2]*A2＋Z3[F2]*A3＋Z4[F2]*A4＋Z5[F2]*A5
R3=Z1[F3]*A1＋Z2[F3]*A2＋Z3[F3]*A3＋Z4[F3]*A4＋Z5[F3]*A5
R4=Z1[F4]*A1＋Z2[F4]*A2＋Z3[F4]*A3＋Z4[F4]*A4＋Z5[F4]*A5
R5=Z1[F5]*A1＋Z2[F5]*A2＋Z3[F5]*A3＋Z4[F5]*A4＋Z5[F5]*A5
…（１６）
但し、R1〜R5 ：上記式（１５）により求めたレスポンス
Z1[F1]〜Z5[F5]：１００ｄｐｉ解像度の原画像信号に基づいて作成された帯域制限画像信号のF1〜F5の周波数におけるレスポンス
Ａ1〜Ａ5 ： 変換関数のパラメータ
【０１３９】
なお、上記第４および第４の実施形態においては、変換処理として上記式（２）に示す周波数強調処理を行っているが、周波数強調処理とダイナミックレンジ圧縮処理とを同時に行なう処理であってもよい。以下この処理を行なう装置について第５の実施形態として説明する。
【０１４０】
図３１は、本発明の第５の実施形態による変換処理の一例としての周波数強調処理およびダイナミックレンジ圧縮処理を行なう装置の構成をボケ画像信号作成手段２とともに示す概略ブロック図である。図３１に示すように、原画像信号Ｓorgと、フィルタリング処理手段１０および補間処理手段１１により作成されたボケ画像信号Ｓuskについて、減算器２１によりそれらの信号の差分が求められ、原画像信号の限られた周波数帯域の成分である帯域制限画像信号（Ｓorg-Ｓus1、Ｓus1-Ｓus2など）が作成される。このようにして求められた帯域制限画像信号は、図３１に示すように第１の変換手段３ａおよび第２の変換手段３ｂにそれぞれ入力され、各変換手段の変換器２２ａ，２２ｂにより処理される。
【０１４１】
第１の変換手段３ａの変換器２２ａによる変換は、上述したように原画像信号Ｓorgの解像度に応じて算出された変換関数を使用して行なう。具体的には、例えば図３２、図３３に示される変換関数、あるいはこれらを組み合わせた関数を基準変換関数とし、原画像信号Ｓorgの解像度に拘わらず処理済画像信号Ｓprocの周波数応答特性が略同一のものとなるように、基準となる変換関数に基づいて解像度に応じた変換関数を算出し、算出された変換関数を使用して行なう。
【０１４２】
ここで、図３２に示す変換関数は、振幅の大きな帯域制限画像信号を抑制するような変換を行なうものであり、周波数帯域の高い帯域制限画像信号の抑制の度合いを、周波数帯域の低い帯域制限画像信号よりも強くするものであるが、これは実際の放射線画像のエッジに含まれている高周波成分が、低周波成分に比べてその振幅が小さいということを考慮したものである。実際の放射線画像においては、かなり急峻なエッジでさえも正確な階段状にはなっておらず、高周波成分になるほどその振幅が小さくなっていることが多い。このため、各周波数成分の振幅に合わせて、周波数の高い帯域制限画像信号ほど小さい振幅から抑制を行なうことが望ましく、本関数によりそれを実現することができる。
【０１４３】
また図３３の関数は、帯域制限画像信号を、帯域制限画像信号の絶対値に基づいて決まる、その絶対値以下の値となるように変換を行なうもので、この関数が低周波帯域を処理する関数であるほど、帯域制限画像信号の絶対値が０近傍の所定の範囲内の値である帯域制限画像信号を変換した際に得られる変換画像信号の絶対値が小さい値であることを特徴とするものである。言い換えれば、これらの関数はそれぞれ、原点を通り、関数の傾きがその関数により処理される値に拘わらず１以下であり、その関数の０近傍における傾きが、低周波帯域を処理する関数であるほど小さいことを特徴とするものである。これらの関数は、変換画像信号を積算して得られる信号を、原画像信号Ｓorgに加えた場合、原画像信号Ｓorgと加算された信号とのつなぎ目、すなわち信号の立ち上がりをより自然なものとするという効果がある。
【０１４４】
同様に、第２の変換手段３ｂの変換器２２ｂによる変換は、例えば図３４に示される関数、前述の図３３の関数、あるいは図３４と図３３の関数を組み合わせたものを使用して行なう。上記各変換手段３ａ，３ｂにより出力された変換画像信号は、それぞれ演算器２３ａ，２３ｂに入力される。ここで演算器２３ａは周波数強調処理に必要な信号を作成するための演算を行ない、演算器２３ｂはダイナミックレンジ圧縮処理に必要な信号を作成するための演算を行なうものである。
【０１４５】
演算器２３ａは、上述した第３の実施形態と同様の周波数強調処理を行なうものである。この演算器２３ａにおいては次のような処理が行われる。まず、上述したように第１の変換手段３ａにより変換された帯域制限画像信号は積算される。そしてこの積算信号が求められると、さらに原画像信号Ｓorg の値に応じた強調度βが乗じられる。
【０１４６】
一方、演算器２３ｂは、ダイナミックレンジ圧縮処理を行なうものである。この演算器２３ｂにおいては次のような処理が行われる。まず、上述したように第２の変換手段３ｂにより変換された帯域制限画像信号は積算される。そして求められた積算信号は原画像信号Ｓorgから減じられる。さらにその減算により作成された差分信号が変換関数により変換され、ダイナミックレンジ圧縮係数が得られる。
【０１４７】
演算器２３ａ，２３ｂにより求められた信号は、それぞれ加算手段２８により原画像信号Ｓorgに加算され処理済画像信号Ｓprocが得られる。
【０１４８】
以上の処理は、従来技術で述べたと同様の下記の式（３）により表される。
Ｓproc＝Ｓorg ＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…ＳusN））…（３）
Ｆusm（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝｛ｆｕ１(Ｓorg−Ｓus1)＋ｆｕ２(Ｓus1 −Ｓus2)＋…
＋ｆｕｋ(Ｓusk-1−Ｓusk )＋…＋ｆｕＮ(ＳusN-1−ＳusN)｝
Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…ＳusN）
＝｛ｆｄ１(Ｓorg−Ｓus1)＋ｆｄ２(Ｓus1 −Ｓus2)＋…
＋ｆｄｋ(Ｓusk-1−Ｓusk )＋…＋ｆｄＮ(ＳusN-1−ＳusN)｝
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆｕｋ(k=1〜N）：第１の変換処理に使用する変換関数
ｆｄｋ(k=1〜N）：第２の変換処理に使用する変換関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを変換する関数））
【０１４９】
ここで、図３１に示す第５の実施形態においては、周波数強調処理に使用する信号とダイナミックレンジ圧縮処理に使用する信号をそれぞれ原画像信号Ｓorgに基づいて作成して最後に加算を行なうものであるが、原画像信号Ｓorgに対して周波数強調処理あるいはダイナミックレンジ圧縮処理のいずれか一方の処理を先に行ない、その結果得られた信号に対して他の一方の処理を行なうようにする場合も含むものとする。但し、一般に放射線画像の低濃度部は撮影時の撮影線量が少ない部分であるためノイズが比較的多く含まれている。これに対し、周波数強調処理では高濃度部ほど強調の度合いが高くなる濃度依存の強調処理を行っているため、原画像信号に対して直接周波数強調処理を行えばこの低濃度部は強調されず、ノイズが強調されることがない。一方、例えば原画像信号に対しダイナミックレンジ圧縮処理を施してから周波数強調処理を行なうと、低濃度部はダイナミックレンジ圧縮処理により濃度が高められることになる。つまり、強調処理においてその部分が強調されることとなり、同時にその部分に含まれるノイズが強調されてしまうことになる。したがって、望ましくは図３１に示す第５の実施形態のように原画像信号Ｓorgに基づいて各信号を作成するのがよい。また処理時間の短縮という意味でも、２つの処理は並列に行なうことが望ましい。
【０１５０】
このように、ダイナミックレンジ圧縮処理を行なう場合においても、原画像信号Ｓorgの解像度に応じて処理済画像信号Ｓprocの周波数応答特性が略同一のものとなるように変換関数を算出することにより、原画像信号Ｓorgの解像度に拘わらず略同一の周波数応答特性を有する処理済画像信号Ｓprocを得ることができる。
【０１５１】
なお、上記第３の実施形態において述べたと同様に、第４および第５の実施形態においても、10本／ｍｍの解像度に対応する１組の基準変換関数を用いるのみならず、表３右側に示すように、利用頻度が高い数種類の基準解像度に対応する変換関数を基準変換関数の組として用意し、基本解像度以外の所定解像度の画像を原画像とする場合には、この処理を施す原画像の所定解像度が上記数種類の基準解像度のいずれに近い値を有するか否かを判断して解像度が近い方を選択するなど、前記組から最適な組を選択し、選択した基準変換関数の組と基本解像度および所定解像度情報を用いて処理を施す画像の解像度に対応した変換関数を算出し、算出した変換関数を用いて変換処理を行なうようにしてもよい。このとき、複数の基本解像度から処理対象に適切な基本解像度を選択する方法としては、特開平10-63838号に記載の方法を用いたり、試行錯誤により設定してもよいし、あるいは、以下のような方法を用いることもできる。
【０１５２】
まず、上記表３に示すような基本解像度（ａ，ｂ，ｃ；ａ＜ｂ＜ｃ：ｃが最高解像度、上記例では、ａ＝５，ｂ＝6.7 ，ｃ＝10）の定義パラメータのテーブルを用意しておき、処理対象となる原画像解像度Ｒｓが基本解像度（ａ，ｂ，ｃ；ａ＜ｂ＜ｃ：ｃが最高解像度）と異なる場合には、基本解像度ａ，ｂ，ｃの内、最も近い基本解像度の定義パラメータを以下のルールにしたがって選択する。
Ｒｓ＞ｂ＋（ｃ−ｂ）／２のとき、ｃ用のパラメータ
ｂ＋（ｃ−ｂ）／２≧Ｒｓ＞ａ＋（ｂ−ａ）／２のとき、ｂ用のパラメータ
ｂ＋（ｃ−ｂ）／２≧Ｒｓのときａ用のパラメータ
この方法を用いると、変換関数の定義パラメータの決定を簡単に行なうことができる。
【０１５３】
また、このように、基本解像度が複数存在し且つ各基本解像度用パラメータがそれぞれ存在するシステムとする場合には、それらの中から１つの基本解像度と基本解像度パラメータを選択し、選択した基本解像度とパラメータおよび処理対象となる画像の解像度情報を処理対象となる画像に関連付けて保存するとよい。これにより、画像をフィルムやモニターに再出力するとき、前回出力した画像に対応する定義パラメータなどを記憶装置から読み出し、この読み出した定義パラメータを使用して処理済画像信号procを生成して、前回出力した画像と同等画質の画像を得ることができる。また、定義パラメータを決定する処理を省略できるので、短時間で画像が出力される。
【０１５４】
なお、上記第３および第４の実施形態において示した周波数強調処理、あるいは第５の実施形態において示したダイナミックレンジ圧縮処理を行なうに際しては、所定解像度画像信号を処理して得た処理済画像信号により表される画像の特性が、所定解像度を規定するナイキスト周波数の１／５（より好ましくは１／２）以下の周波数において、原画像信号（上記例では基準原画像信号Ｓorg０ など）を処理して得た処理済画像信号Ｓprocにより表される画像の特性と略一致するように変換関数を算出するのが好ましい。以下、この点について、強調特性と処理済画像信号が表す画像のレスポンスの関係を示した図３５を参照して説明する。
【０１５５】
例えば、基本解像度として、５本／ｍｍ，6.7 本／ｍｍ，10本／ｍｍが存在し、各々に対してそれぞれ変換関数群が定義されているシステムの場合において、例えば処理対象となる画像の解像度が４本／ｍｍである場合には、最も近い基本解像度は５本／ｍｍとなるので、５本／ｍｍ用の変換関数群を用いて、上述したようにして４本／ｍｍ用の変換関数群を計算することになる。
【０１５６】
ここで、各帯域制限画像信号を得る際には画像信号を１／２単位で縮小しているため、４本／ｍｍ用の画像信号における、各帯域制限画像信号のうちの最高周波数帯域信号の特性幅（例えばレスポンス０．８以上の略平坦部分）は、少なくとも1.0cycle/mm〜2.0cycle/mm位となり、各帯域制限画像信号のゲインを変換関数で制御するような変換処理を行なう場合において、例えば1.0 cycle/mm近傍の特性を基本解像度特性に一致させようとすれば2.0 cycle/mm近傍のレスポンスを持ち上げることができず2.0 cycle/mm近傍では基本解像度特性と一致しなくなり、逆に2.0 cycle/mm近傍の特性を基本解像度特性に一致させようとすれば、1.0 cycle/mm近傍のレスポンスも持ち上がるため1.0 cycle/mm近傍では基本解像度特性と一致しなくなるというように、1.0cycle/mm〜2.0cycle/mmの間のレスポンスを微調整することが難しい（できない）。
【０１５７】
したがって、基本解像度５本／ｍｍ用の変換関数群として、図３５（Ａ）に示すように、比較的低周波側から比較的高周波側に亘って強調するような特性となる変換関数群（パラメータ）の設定となっている場合において、４本／ｍｍ用の変換関数群を計算する際には、最高周波数帯域信号を除く低周波側の帯域信号に対する変換関数として、最高周波数帯域信号に対する変換関数によって1.0 cycle/mm近傍が持ち上がる分を補正するような値にすることができ、計算された変換関数群に基づく処理済画像信号Ｓprocのレスポンスを、図３５（Ｂ）に示すように、1.0cycle/mm〜2.0cycle/mm近傍においても、基本解像度特性と略一致させることができる。
【０１５８】
しかしながら、基本解像度５本／ｍｍ用の変換関数群として、図３５（Ｃ）に示すように、比較的高周波成分のみを強調するなどナイキスト周波数に近い周波数帯域で急激な特性変化を持つような特性となる変換関数群の設定となっている場合には、上述のように1.0cycle/mm〜2.0cycle/mmの間のレスポンスを微調整することが難しいという問題があるため、処理済画像信号Ｓprocの全周波数帯域に亘ってレスポンスを基本解像度特性と一致させるということはできない。この場合、図３５（Ｄ）に示すように、視覚的に影響が大きな（画質に影響を与える度合いの大きい）低周波側のレスポンスを重視し、例えばナイキスト周波数（本例では2.0 cycle/mm）の１／５以下の周波数において、処理済画像信号Ｓprocのレスポンスが基本解像度特性と一致するようにするのが好ましい。また、ナイキスト周波数の１／２以下の低周波帯域成分を重視して基本解像度特性に合わせるようにすれば、よりよい結果となる。
【０１５９】
また、上記第３から第５の実施形態においては、入力される原画像信号Ｓorgの解像度に応じて変換関数を算出しているが、入力する原画像信号Ｓorgのレスポンスに関する特性情報に基づいて変換関数を補正してもよい。すなわち、原画像信号Ｓorgを取得するための装置としては、蓄積性蛍光体シートから放射線画像を読み取る装置の他、半導体センサなどもあり、得られる原画像信号Ｓorgのレスポンスが装置に応じて異なるものである。このため、装置に応じた原画像信号Ｓorgのレスポンスに関する特性情報に基づいて変換関数を補正することにより、原画像信号Ｓorgを取得する装置に拘わらず、一定の周波数応答特性を有する処理済画像信号Ｓprocを得ることができる。例えば、ある装置において得られる原画像信号Ｓorgの２cycle/mmにおけるレスポンスが基準となるレスポンスから４０％劣っているような場合には、その旨を特性情報として本発明による画像処理装置に入力し、画像処理装置においては２cycle/mmのレスポンスが基準レスポンスとなるように変換関数を補正し、補正された変換関数により処理済画像信号Ｓprocを得る。
【０１６０】
さらに、上記第４および第５の実施形態においては、変換関数を非線形関数として帯域制限画像信号に対して変換処理を施しているが、非線形関数に限定されるものではなく、変換関数を線形関数や定数としてもよい。
【０１６１】
また、上記第３から第５の実施形態においては、原画像信号からフィルタリングおよび補間拡大によりボケ画像信号を求め、原画像信号およびボケ画像信号から帯域制限画像信号を作成しているが、これに限定されるものではなく、例えば原画像信号をウェーブレット変換やラプラシアンピラミッドなどの手法により多重解像度に変換し、変換された各解像度毎の画像信号からボケ画像信号を作成して帯域制限画像信号を作成してもよい。
【図面の簡単な説明】
【０１６２】
【図１】本発明の実施形態による画像処理装置の構成を示す概略ブロック図
【図２】原画像信号をウェーブレット変換して階層毎に符号化する状態を示す図
【図３】第１の実施形態におけるボケ画像信号作成手段の構成を示す概略ブロック図
【図４】フィルタリング処理に使用されるフィルタの一例を示す図
【図５】低解像度画像信号作成処理の詳細を示す図
【図６】補間拡大処理に使用されるフィルタの一例を示す図
【図７】周波数強調処理装置の一例を示す図
【図８】帯域制限画像信号の周波数応答特性を示す図（基準解像度）
【図９】帯域制限画像信号の周波数応答特性を示す図（１／２解像度）
【図１０】帯域制限画像信号の周波数応答特性を示す図（１／４解像度）
【図１１】帯域制限画像信号の周波数応答特性を示す図（１／８解像度）
【図１２】帯域制限画像信号の周波数応答特性を示す図（１／１６解像度）
【図１３】第１の実施形態の動作を示すフローチャート
【図１４】基準解像度画像から得られる６つのボケ画像信号と、最低解像度から所定解像度に最も近い周波数帯域までの低解像度画像から得られる低解像度ボケ画像信号との対応関係を示す図
【図１５】基準解像度画像から得られる６つの帯域制限画像信号と、各解像度画像に基づく原画像信号および低解像度ボケ画像信号から得られる低解像度帯域制限画像信号との対応関係を示す図
【図１６】基準解像度画像と各解像度画像とにおける変換関数の対応関係を示す図
【図１７】第２の実施形態におけるボケ画像信号作成手段の構成を示す概略ブロック図
【図１８】基準解像度の画像と拡大画像における、各帯域制限画像信号の対応関係を示す図
【図１９】基準解像度の画像と拡大画像における、各変換関数の対応関係を示す図
【図２０】本発明の第３の実施形態による画像処理装置の構成を示す概略ブロック図
【図２１】処理済画像信号の周波数応答特性を示す図
【図２２】帯域制限画像信号の周波数応答特性を示す図（３００ｄｐｉ）
【図２３】第３の実施形態により得られた処理済画像信号の周波数応答特性を示す図
【図２４】第３の実施形態の処理を示すフローチャート
【図２５】変換関数の選択状態を示す図
【図２６】解像度に応じた帯域制限画像信号の数を示す図
【図２７】帯域制限画像信号の周波数応答特性を示す概略図（２００ｄｐｉ）
【図２８】処理済画像信号の周波数応答特性を示す図
【図２９】レスポンスX1[F1]と関数X1[i]との対応関係を示す図
【図３０】関数X1[i]と関数Z1[i]との関係を示す図
【図３１】ダイナミックレンジ圧縮処理装置の一例を示す図
【図３２】第１の変換手段における変換関数の一例を表す図
【図３３】第１の変換手段における変換関数の他の例を表す図
【図３４】第２の変換手段における変換関数の一例を表す図
【図３５】強調特性と処理済画像信号が表す画像のレスポンスの関係を示す図
【符号の説明】
【０１６３】
１ 画像処理装置
２ ボケ画像信号作成手段
３ 周波数強調処理手段、変換処理手段
４ パラメータ設定手段
５ 拡大率入力手段
６ 拡大縮小手段
10 フィルタリング処理手段
11，31 補間処理手段
21 減算器
22 変換器
23 演算器
30 復号化手段
40 変換関数算出手段
50 解像度情報入力手段
【図１】
0004598115.tif 000007
【図２】
0004598115.tif 000008
【図３】
0004598115.tif 000009
【図４】
0004598115.tif 000010
【図５】
0004598115.tif 000011
【図６】
0004598115.tif 000012
【図７】
0004598115.tif 000013
【図８】
0004598115.tif 000014
【図９】
0004598115.tif 000015
【図１０】
0004598115.tif 000016
【図１１】
0004598115.tif 000017
【図１２】
0004598115.tif 000018
【図１３】
0004598115.tif 000019
【図１４】
0004598115.tif 000020
【図１５】
0004598115.tif 000021
【図１６】
0004598115.tif 000022
【図１７】
0004598115.tif 000023
【図１８】
0004598115.tif 000024
【図１９】
0004598115.tif 000025
【図２０】
0004598115.tif 000026
【図２１】
0004598115.tif 000027
【図２２】
0004598115.tif 000028
【図２３】
0004598115.tif 000029
【図２４】
0004598115.tif 000030
【図２５】
0004598115.tif 000031
【図２６】
0004598115.tif 000032
【図２７】
0004598115.tif 000033
【図２８】
0004598115.tif 000034
【図２９】
0004598115.tif 000035
【図３０】
0004598115.tif 000036
【図３１】
0004598115.tif 000037
【図３２】
0004598115.tif 000038
【図３３】
0004598115.tif 000039
【図３４】
0004598115.tif 000040
【図３５】
0004598115.tif 000041