(57)【特許請求の範囲】
【請求項１】
原画像を表す原画像信号に、該原画像の高周波成分に関する信号を加算して該原画像の高周波成分を強調するとともに該原画像の低周波成分に関する信号を加算して該原画像のダイナミックレンジを圧縮する画像処理方法において、
前記原画像信号に基づいて、互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成し、
前記原画像信号および前記複数の非鮮鋭マスク画像信号、または前記複数の非鮮鋭マスク画像信号に基づいて、前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、
該各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第１の変換処理を施して複数の第１の変換画像信号を作成し、該各変換画像信号を積算することにより前記原画像信号の前記高周波成分に関する信号を得るとともに、前記各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第２の変換処理を施して複数の第２の変換画像信号を作成し、該各変換画像信号を積算することにより積算信号を作成し、該積算信号を前記原画像信号から減ずることにより差分信号を作成し、該差分信号を変換することにより前記原画像信号の前記低周波成分に関する信号を得、
前記原画像信号に前記高周波成分に関する信号および前記低周波成分に関する信号を加算することにより、該原画像信号の高周波成分を強調するとともにダイナミックレンジを圧縮することを特徴とする画像処理方法。
【請求項２】
前記帯域制限画像信号の作成、前記第１および第２の変換画像信号の作成、前記高周波成分に関する信号の作成、前記低周波成分に関する信号の作成、および前記原画像信号への前記高周波成分に関する信号と低周波成分に関する信号の加算を、下記の式
0003700804.tif 000006
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：非鮮鋭マスク画像信号
ｆuk(k=1〜N）：前記第１の変換処理に使用する関数
ｆdk(k=1〜N）：前記第２の変換処理に使用する関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを変換する関数））
にしたがって行うことを特徴とする請求項１記載の画像処理方法。
【請求項３】
原画像を表す原画像信号に、該原画像の高周波成分に関する信号を加算して該原画像の高周波成分を強調するとともに該原画像の低周波成分に関する信号を加算して該原画像のダイナミックレンジを圧縮する画像処理装置において、
前記原画像信号に基づいて、互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成する非鮮鋭マスク画像信号作成手段と、
前記原画像信号および前記複数の非鮮鋭マスク画像信号、または前記複数の非鮮鋭マスク画像信号に基づいて、前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成する帯域制限画像信号作成手段と、
前記各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第１の変換処理を施して複数の第１の変換画像信号を作成する第１の変換手段と、
前記第１の変換手段により作成された各変換画像信号を積算することにより第１の積算信号を作成する第１の積算手段と、
前記第１の積算信号に基づいて前記原画像信号の前記高周波成分に関する信号を得る周波数強調処理手段と、
前記各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第２の変換処理を施して複数の第２の変換画像信号を作成する第２の変換手段と、
前記第２の変換手段により作成された各変換画像信号を積算することにより第２の積算信号を作成する第２の積算手段と、
前記第２の積算信号を前記原画像信号から減ずることにより差分信号を作成する減算手段と、
前記差分信号を変換することにより前記原画像信号の前記低周波成分に関する信号を得るダイナミックレンジ圧縮処理手段と、
前記原画像信号に前記高周波成分に関する信号および前記低周波成分に関する信号を加算して、該原画像信号の高周波成分を強調するとともにダイナミックレンジを圧縮する加算手段とを備えてなることを特徴とする画像処理装置。
【請求項４】
前記帯域制限画像信号の作成、前記第１および第２の変換画像信号の作成、前記高周波成分に関する信号の作成、前記低周波成分に関する信号の作成、および前記原画像信号への前記高周波成分に関する信号と低周波成分に関する信号の加算を、下記の式
0003700804.tif 000007
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：非鮮鋭マスク画像信号
ｆuk(k=1〜N）：前記第１の変換処理に使用する関数
ｆdk(k=1〜N）：前記第２の変換処理に使用する関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを変換する関数））
にしたがって行うことを特徴とする請求項３記載の画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像信号に対し、その画像信号の非鮮鋭マスクを利用した周波数強調処理と、同じく非鮮鋭マスクを利用したダイナミックレンジ圧縮処理の両方を同時に施す画像処理方法および装置に関するものである。
【０００２】
【従来の技術】
従来、本出願人により、放射線画像を表す画像信号に対し非鮮鋭マスク画像信号（以下、ボケ画像信号という）を用いた周波数強調処理あるいはダイナミックレンジ圧縮処理などを行って放射線画像の診断性能を向上させる数々の画像処理方法および装置が提案されている（特開昭55-163472号、同55-87953号、特開平3-222577号など）。
【０００３】
また、特開平6-292009号には、１つの画像信号に対して周波数強調処理とダイナミックレンジ圧縮処理の両方を施す画像処理方法および装置が開示されている。この方法および装置は、ダイナミックレンジ圧縮によって劣化したコントラスト分解能を補正すること、およびダイナミックレンジ圧縮によって「コントラストが低下したように見える」という感覚的なものを補正することを目的として周波数強調処理を行うものであり、より詳細には下記の式（１）
Ｓproc＝Ｓorg＋ｆ1（Ｓus1）＋ｆ2（Ｓorg−Ｓus2）  …  （１）
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓus1，Ｓus2：ボケ画像信号
ｆ1：単調減少関数
ｆ2：単調増加関数）
に従って処理を行っている。
【０００４】
【発明が解決しようとする課題】
しかしながら、上記特開平6-292009号の方法では、信号値の変化が急激なエッジなどの部分において、作成されるボケ画像信号がエッジの影響を受けてしまい、このボケ画像信号により処理が行われた画像に、オーバーシュート、アンダーシュートあるいは偽輪郭などのアーチファクトが生じてしまうという問題がある。
【０００５】
本発明はこの問題に鑑みて、アーチファクトによる画質劣化を生じさせずに周波数強調処理とダイナミックレンジ圧縮処理の両方の処理を画像信号に対して施す画像処理方法および装置を提供することを目的とするものである。
【０００６】
【課題を解決するための手段】
本発明の画像処理方法および装置は、原画像を表す原画像信号に、該原画像の高周波成分に関する信号を加算して該原画像の高周波成分を強調するとともに該原画像の低周波成分に関する信号を加算して該原画像のダイナミックレンジを圧縮する画像処理方法において、
前記原画像信号に基づいて、互いに周波数応答特性が異なる複数のボケ画像信号を作成し、
前記原画像信号および前記複数のボケ画像信号、または前記複数のボケ画像信号に基づいて、前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、
該各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第１の変換処理を施して複数の第１の変換画像信号を作成し、該各変換画像信号を積算することにより前記原画像信号の前記高周波成分に関する信号を得るとともに、前記各帯域制限画像信号のうち少なくとも１つの信号について、該帯域制限画像信号の少なくとも一部を小さくするような第２の変換処理を施して複数の第２の変換画像信号を作成し、該各変換画像信号を積算することにより積算信号を作成し、該積算信号を前記原画像信号から減ずることにより差分信号を作成し、該差分信号を変換することにより前記原画像信号の前記低周波成分に関する信号を得、
前記原画像信号に前記高周波成分に関する信号および前記低周波成分に関する信号を加算することにより、該原画像信号の高周波成分を強調するとともにダイナミックレンジを圧縮することを特徴とするものである。
【０００７】
ここで、ボケ画像信号とは、画素数は原画像信号と同じであるが原画像信号よりも鮮鋭度が低い画像を表す画像信号であり、ボケ画像信号の作成は、原画像信号の画素に対して所定間隔ごとに所定のフィルタリング処理を施すことによって画素を間引きし、その後間引きした数分の画素を所定の補間方法によって補間することにより行われる。
【０００８】
ここで、前記帯域制限画像信号の作成、前記第１および第２の変換画像信号の作成、前記高周波成分に関する信号の作成、前記低周波成分に関する信号の作成、および前記原画像信号への前記高周波成分に関する信号と低周波成分に関する信号の加算は、下記の式
0003700804.tif 000002
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆuk(k=1〜N）：前記第１の変換処理に使用する関数
ｆdk(k=1〜N）：前記第２の変換処理に使用する関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを変換する関数））
にしたがって行うことが望ましい。
【０００９】
【発明の効果】
本発明の画像処理方法および装置によれば、原画像信号のボケ画像信号を作成し、そのボケ画像信号を使用して原画像信号の周波数強調処理とダイナミックレンジ圧縮処理の両方の処理を行う際に、互いに周波数応答特性が異なる複数のボケ画像信号を作成してボケ画像信号同士の差をとることにより周波数帯域ごとの信号を作成し、周波数強調処理やダイナミックレンジ圧縮処理に使用する所定の周波数成分の信号を、その周波数帯域ごとの信号をそれぞれ変換してから積算することによって得るので、アーチファクトの原因となる周波数成分の信号を効率よく抑制することができ、高画質の処理済画像信号を得ることができる。
【００１０】
【発明の実施の形態】
以下、本発明の画像処理方法および装置の一実施の形態について、図面を参照して詳細に説明する。以下に示す画像処理装置は、蓄積性蛍光体シートに記録された人体の放射線画像を読み取って得た画像信号に対して、その画像が診断に適した画像となるように、ボケ画像信号を使用して画像処理を施すものであり、処理された画像信号は主としてフィルムに記録され、診断に用いられる。
【００１１】
図１に示すように本発明による画像処理装置は、入力された原画像信号Ｓorgに基づいて互いに周波数応答特性が異なる多重解像度のボケ画像信号Ｓusk （ｋ＝１〜ｎ）を作成するボケ画像信号作成手段１と、ボケ画像信号作成手段１において作成されたボケ画像信号Ｓuskに基づいて複数の帯域制限画像信号を作成する帯域制限画像信号作成手段２と、帯域制限画像信号作成手段２において作成された複数の帯域制限画像信号のうち少なくとも１つの信号について、その帯域制限画像信号の少なくとも一部を小さくするような変換処理を施す第１の変換手段3aと、第１の変換手段3aにより変換された帯域制限画像信号を積算して積算信号を作成する第１の積算手段4aと、第１の積算手段4aにより作成された積算信号に所定の強調係数を乗ずる周波数強調処理手段５と、帯域制限画像信号作成手段２において作成された複数の帯域制限画像信号のうち少なくとも１つの信号について、その帯域制限画像信号の少なくとも一部を小さくするような変換処理を施す第２の変換手段3bと、第２の変換手段3bにより変換された帯域制限画像信号を積算して積算信号を作成する第２の積算手段4bと、第２の積算手段4bにより作成された積算信号を原画像信号Ｓorgから減算することにより差分信号を得る減算手段６と、減算手段６により得られた差分信号を変換するダイナミックレンジ圧縮処理手段７と、周波数強調処理手段５およびダイナミックレンジ圧縮処理手段７により作成された信号を原画像信号Ｓorgに加算する手段８とからなる。
【００１２】
はじめに、ボケ画像信号作成手段１において行われる処理について説明する。図２はボケ画像信号作成手段１の処理を示すブロック図である。図２に示すように原画像を表すデジタルの原画像信号Ｓorgがフィルタリング処理手段10においてローパスフィルタによりフィルタリングされる。このローパスフィルタとしては、例えば図３に示すような５×１のグリッド状の１次元ガウス分布に略対応したフィルタＦを用いることができる。このフィルタＦは下記の式（２）
【００１３】
【数１】
0003700804.tif 000003
【００１４】
において、σ＝１としたものである。ここで、フィルタＦとしてガウス信号を用いるのは、ガウス信号は周波数空間および実空間の双方において、局在性がよいためである。
【００１５】
そして、このフィルタＦにより原画像の画素のｘ方向およびｙ方向に対してフィルタリング処理を施すことにより、原画像信号Ｓorg 全体に対してフィルタリング処理が施される。
【００１６】
フィルタリング処理手段10においては、このようなフィルタＦにより以下のようにしてフィルタリング処理が施される。図４はフィルタリング処理の詳細を表す図である。図４に示すように、原画像信号Ｓorgに対し、図３に示すフィルタＦにより１画素おきにフィルタリング処理が施される。そしてこのフィルタリング処理により、フィルタリング処理画像信号Ｂ1が得られる。このフィルタリング処理画像信号Ｂ1は、原画像に対するサイズが１／４（ｘ方向、ｙ方向にそれぞれ１／２）となっている。次いで、再度フィルタＦによりフィルタリング処理画像信号Ｂ1に対して１画素おきにフィルタリング処理が施される。そしてこのようなフィルタＦによるフィルタリング処理が繰り返し施されることにより、ｎ個のフィルタリング処理画像信号Ｂk （k=1〜n）が得られる。このフィルタリング処理画像信号Ｂkは、原画像に対してサイズが１／２2kとなっている。この際、フィルタリング処理画像信号Ｂk の周波数特性は図５に示されるようなものとなる。図５に示されるようにフィルタリング処理画像信号Ｂk のレスポンスはｋが大きいほど高周波成分が除去されたものとなっている（但し、図５においてはk=1〜3としている）。
【００１７】
なお、上記実施の形態においては、図３に示す１次元フィルタＦにより、原画像のｘ方向およびｙ方向に対してフィルタリング処理が施されるようになっているが、図６に示すような５×５の２次元フィルタにより原画像信号Ｓorg およびフィルタリング処理画像信号に対して一度にフィルタリング処理が施されるようにしてもよい。
【００１８】
次いで、このようにして得られたフィルタリング処理画像信号Ｂkに対して図２に示される補間演算処理手段11において補間演算処理が施され、これにより原画像と同一の大きさの多重解像度のボケ画像が得られる。以下この補間演算処理について説明する。
【００１９】
補間演算の方法としては、Ｂスプラインによる方法など種々の方法が挙げられるが、本発明による実施の形態においては、ローパスフィルタとしてガウス信号に基づくフィルタＦを用いているため、補間演算を行うための補間係数としてもガウス信号を用いるものとする。ここで、ガウス信号を用いた補間係数とは、下記の式（３）
【００２０】
【数２】
0003700804.tif 000004
【００２１】
において、σ＝２k-1と近似したものを用いる。
【００２２】
フィルタリング処理画像信号Ｂ1を補間する際には、ｋ＝１であるためσ＝１となる。上記式（３）において、σ＝１としたときの補間を行うためのフィルタは、図７に示すように５×１の１次元フィルタとなる。まずフィルタリング処理画像信号Ｂ1に対して１画素おきに値が０の画素を１つずつ補間することによりフィルタリング処理画像信号Ｂ1 が原画像と同一のサイズに拡大される。このように値が０の画素が補間されたフィルタリング処理画像信号Ｂ1を１次元的に図８に示す。そして、この補間されたフィルタリング処理画像信号Ｂ1に対して上述した図７に示すフィルタＦ1によりフィルタリング処理が施される。
【００２３】
ここで、図７に示すフィルタＦ1は５×１のフィルタであるが、図８に示すようにフィルタリング処理画像信号Ｂ1は１画素おきに値が０の画素が補間されている。したがって、フィルタリング処理画像Ｂ1に対してフィルタＦ1により施されるフィルタリング処理は、実質的には２×１のフィルタ（0.5, 0.5）および３×１のフィルタ（0.1, 0.8, 0.1）の２種類のフィルタにより、施されるフィルタリング処理と等価なものとなる。そしてこのフィルタリング処理により、原画像信号Ｓorg と同一データ数、すなわち原画像と同一サイズのボケ画像の信号Ｓus1が得られる。
【００２４】
次いで、フィルタリング処理画像信号Ｂ2に対してフィルタリング処理が施される。フィルタリング処理画像信号Ｂ2を補間する際には、ｋ＝２であるため、σ＝２となる。上記式（３）において、σ＝２としたときの補間を行うためのフィルタは、図９に示すように１１×１の１次元フィルタとなる。そして、まずフィルタリング処理画像信号Ｂ2に対して図11のように１画素おきに値が０の画素を３つずつ補間することによりフィルタリング処理画像信号Ｂ2 が原画像と同一のサイズに拡大される。値が０の画素が補間されたこのフィルタリング処理画像信号Ｂ2に対して上述した図９に示すフィルタＦ2によるフィルタリング処理が施される。
【００２５】
ここで、図９に示すフィルタＦ2は１１×１のフィルタであるが、図10に示すようにフィルタリング処理画像信号Ｂ2は１画素おきに値が０の画素が３つずつ補間されている。したがって、フィルタリング処理画像信号Ｂ2に対してフィルタＦ2により施されるフィルタリング処理は、実質的には２×１のフィルタ（0.5, 0.5）および３×１のフィルタ（0.3, 0.65, 0.05）、（0.13, 0.74, 0.13）および（0.05, 0.65, 0.3）の４種類のフィルタにより施されるフィルタリング処理と等価なものとなる。そしてこのフィルタリング処理により、原画像信号Ｓorg と同一データ数のボケ画像信号Ｓus2が得られる。
【００２６】
そしてこのようなフィルタリング処理が全てのフィルタリング処理画像信号Ｂkに対して行われる。フィルタリング処理画像信号Ｂkを補間する際には、上記式（３）に基づいて、３×２k −１の長さのフィルタを作成し、フィルタリング処理画像信号Ｂk の各画素の間に値が０の画素を２k −１個ずつ補間することにより、原画像と同一サイズに拡大する。この値が０の画素が補間されたフィルタリング処理画像信号Ｂkに対して３×２k −１の長さのフィルタにより、フィルタリング処理が行われる。
【００２７】
ここで、この３×２k −１の長さのフィルタにより施されるフィルタリング処理は、２k 周期で長さが２または３のフィルタにより施されるフィルタリング処理と等価なものとなる。そしてこのフィルタリング処理により、ｎ個のボケ画像信号Ｓusk が得られる。このボケ画像信号Ｓuskを可視像として表すと、結果としてそれぞれ解像度が異なる、すなわち周波数応答特性が異なる多重解像度のボケ画像となる。このように、フィルタは長くなるものの、実質的には長さが２または３のフィルタによりフィルタリング処理を施していることと同様であるため、演算量は、フィルタが長くなろうともそれほど多くはならないものである。このため、演算量を減らして、多重解像度のボケ画像信号Ｓuskを高速に作成することが可能となる。
【００２８】
なお、本実施の形態においては、長さが３×２k −１の１次元のフィルタにより画像のｘ方向およびｙ方向にフィルタリング処理が施されるようになっているが、予め２次元状のフィルタを作成しておき、このフィルタによりフィルタリング処理画像に対するフィルタリング処理を施すことにより、ボケ画像信号Ｓusk が得られるようにしてもよい。この場合、フィルタリング処理画像に対して、補間演算を行うためのフィルタリング処理に使用するフィルタは、（３×２k −１）×（３×２k −１）のフィルタとなるが、上述した１次元のフィルタを用いる場合と同様に、これは２k 周期で２×２または３×３のフィルタにより施されるフィルタリング処理と等価なものとなり、上述した１次元フィルタを用いる場合と同様に、フィルタのサイズが大きくなっても、フィルタリング処理を行うための演算量は実質的にはそれほど大きくならないものである。
【００２９】
このようにして得られたボケ画像信号Ｓusk の周波数特性を図11に示す。図11に示すように、ボケ画像信号Ｓusk のｋの値が大きくなるほど、原画像信号Ｓorg の高周波成分が除去された信号となる。
【００３０】
図12は、上記図２のボケ画像信号作成手段を含む本発明の画像処理装置の一実施の形態について、その全体の構成を示す図であるが、この図に示されるように、ボケ画像信号作成手段１において生成された各ボケ画像信号は、次に、帯域制限画像信号作成手段２および変換手段3a、3bにおいて処理される。図12に示すように、まず、原画像信号Ｓorgおよびボケ画像信号作成手段１において作成された複数のボケ画像信号Ｓuskに基づいて帯域制限画像信号が作成されるが、この帯域制限画像信号は減算器21により、原画像信号Ｓorgとボケ画像信号Ｓus1、あるいは互いに隣接する周波数帯域同士のボケ画像信号Ｓuskの減算を行うことにより得られる。すなわち、Ｓorg−Ｓus1、Ｓus1−Ｓus2、…ＳusN-1−ＳusNを順次計算することにより、複数の帯域制限画像信号が求められる。この帯域制限画像信号の周波数特性を図13に示す。図13に示すように、帯域制限画像信号はボケ画像信号Ｓusk のｋの値が大きくなるほど、原画像信号Ｓorg の低周波数成分の帯域を表す信号となる。
【００３１】
このようにして求められた帯域制限画像信号は、図12に示されるように第１の変換手段3aおよび第２の変換手段3bにそれぞれ入力され、各変換手段の変換器22により処理される。
【００３２】
変換手段3aの変換器22aによる変換は、例えば図14、図15に示される関数、あるいはこれらを組み合わせた関数を使用して行う。ここで各帯域制限画像信号を変換する関数fu1〜fuNは全て同じ関数であってもよいが、できれば処理済画像信号Ｓproc の周波数特性が所望の周波数特性となるように各関数を個別に定義することが望ましい。
【００３３】
例えば図14の関数は、振幅の大きな帯域制限画像信号を抑制するような変換を行うものであり、周波数帯域の高い帯域制限画像信号の抑制の度合いを、周波数帯域の低い帯域制限画像信号よりも強くするものであるが、これは実際の放射線画像のエッジに含まれている高周波成分が、低周波成分に比べてその振幅が小さいということを考慮したものである。実際の放射線画像においては、かなり急峻なエッジでさえも正確な階段状にはなっておらず、高周波成分になるほどその振幅が小さくなっていることが多い。このため、各周波数成分の振幅に合わせて、周波数の高い帯域制限画像信号ほど小さい振幅から抑制を行うことが望ましく、本関数によりそれを実現することができる。
【００３４】
また図15の関数は、帯域制限画像信号を、帯域制限画像信号の絶対値に基づいて決まる、その絶対値以下の値となるように変換を行うもので、この関数が低周波帯域を処理する関数であるほど、帯域制限画像信号の絶対値が０近傍の所定の範囲内の値である帯域制限画像信号を変換した際に得られる変換画像信号の絶対値が小さい値であることを特徴とするものである。言い換えれば、これらの関数はそれぞれ、原点を通り、関数の傾きがその関数により処理される値にかかわらず１以下であり、その関数の０近傍における傾きが、低周波帯域を処理する関数であるほど小さいことを特徴とするものである。これらの関数は、変換画像信号を積算して得られる信号を、原画像信号Ｓorgに加えた場合、原画像信号Ｓorg と加算された信号とのつなぎ目、すなわち信号の立ち上がりをより自然なものとするという効果がある。
【００３５】
同様に、変換手段3bの変換器22bによる変換は、例えば図16に示される関数、前述の図15の関数、あるいは図16と図15の関数を組み合わせたものを使用して行う。変換手段3aの場合と同様、各帯域制限画像信号を変換する関数fd1〜fdNは全て同じ関数であってもよいが、できれば処理済画像信号Ｓproc の周波数特性が所望の周波数特性となるように各関数を個別に定義することが望ましい。
【００３６】
上記各変換手段3a、3bにより出力された変換画像信号は、それぞれ演算器23a、23bに入力される。ここで演算器23aは周波数強調処理に必要な信号を作成するための演算を行い、演算器23bはダイナミックレンジ圧縮処理に必要な信号を作成するための演算を行うものである。
【００３７】
演算器23aは、上述した第１の積算手段4aおよび周波数強調処理手段５を内包するものである。この演算器23aにおいては次のような処理が行われる。まず、上述したように変換手段3aにより変換された帯域制限画像信号は積算される。そしてこの積算信号が求められると、周波数強調処理手段５において原画像信号Ｓorg の値に応じた強調度βが乗じられる。
【００３８】
一方、演算器23bは、上述した第２の積算手段4b、減算手段６およびダイナミックレンジ圧縮処理手段７を内包するものである。この演算器23bにおいては次のような処理が行われる。まず、上述したように変換手段3aにより変換された帯域制限画像信号は積算される。そして求められた積算信号は原画像信号Ｓorgから減じられる。さらにその減算により作成された差分信号が所定の関数により変換され、ダイナミックレンジ圧縮係数が得られる。
【００３９】
演算器23aおよび23bにより求められた信号は、それぞれ加算手段８により原画像信号Ｓorgに加算され処理済画像信号Ｓprocが得られる。
【００４０】
以上の処理は、下記の式（４）
0003700804.tif 000005
（但し、Ｓproc：処理済画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆuk(k=1〜N）：前記第１の変換処理に使用する関数
ｆdk(k=1〜N）：前記第２の変換処理に使用する関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数
Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定められるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを変換する関数））
により表すことができる。
【００４１】
すなわち、本発明によれば、上述のように周波数帯域ごとの信号を作成して、変換処理を周波数帯域ごとに行うため、従来方法に比べて画像のエッジ部におけるオーバーシュート、アンダーシュートあるいは偽輪郭などのアーチファクトが発生しにくく、特に、各関数を調整を行いながら個別に設定するようにすれば、より高い画質の処理済画像を得ることが可能となる。
【００４２】
なお、上述の実施の形態は周波数強調処理に使用する信号とダイナミックレンジ圧縮処理に使用する信号をそれぞれ原画像信号Ｓorgに基づいて作成して最後に加算を行うものであるが、本発明は、原画像信号Ｓorgに対して周波数強調処理あるいはダイナミックレンジ圧縮処理のいずれか一方の処理を先に行い、その結果得られた信号に対して他の一方の処理を行うようにする場合も含むものとする。
【００４３】
但し、一般に放射線画像の低濃度部は撮影時の撮影線量が少ない部分であるためノイズが比較的多く含まれている。これに対し、周波数強調処理では高濃度部ほど強調の度合いが高くなる濃度依存の強調処理を行っているため、原画像信号に対して直接周波数強調処理を行えばこの低濃度部は強調されず、ノイズが強調されることがない。
【００４４】
一方、例えば原画像信号に対しダイナミックレンジ圧縮処理を施してから周波数強調処理を行うと、低濃度部はダイナミックレンジ圧縮処理により濃度が高められることになる。つまり、強調処理においてその部分が強調されることとなり、同時にその部分に含まれるノイズが強調されてしまうことになる。
【００４５】
したがって、望ましくは上記実施の形態のように原画像信号Ｓorgに基づいて各信号を作成するのがよい。また処理時間の短縮という意味でも、２つの処理は並列に行うことが望ましい。
【図面の簡単な説明】
【図１】本発明による画像処理装置の概念を表す図
【図２】ボケ画像信号作成手段の詳細を表す図
【図３】ボケ画像信号作成手段において用いられるフィルタ（１次元）を表す図
【図４】フィルタリング処理の詳細を表す図
【図５】フィルタリング処理画像信号の周波数特性を表す図
【図６】フィルタリング処理手段において用いられるフィルタ（２次元）を表す図
【図７】フィルタリング処理画像信号Ｂ1 の補間演算に用いられるフィルタを表す図
【図８】補間演算の詳細を表す図
【図９】フィルタリング処理画像信号Ｂ2 の補間演算に用いられるフィルタを表す図
【図１０】補間演算の詳細を表す図
【図１１】ボケ画像信号の周波数特性を表す図
【図１２】本発明の画像処理装置の一実施の形態における全体構成を示す図
【図１３】帯域制限画像信号の周波数特性を表す図
【図１４】第１の変換手段における帯域制限画像信号の変換処理の一例を表す図
【図１５】第１の変換手段における帯域制限画像信号の変換処理の他の例を表す図
【図１６】第２の変換手段における帯域制限画像信号の変換処理の一例を表す図
【符号の説明】
１  ボケ画像信号作成手段
２  帯域制限画像信号作成手段
３  変換手段
４  積算手段
５  周波数強調処理手段
６  減算手段
７  ダイナミックレンジ圧縮処理手段
８  加算手段
10  フィルタリング処理手段
11  補間処理手段
21  減算器
22  変換器
23  演算器
【図１】
0003700804.tif 000008
【図２】
0003700804.tif 000009
【図３】
0003700804.tif 000010
【図４】
0003700804.tif 000011
【図５】
0003700804.tif 000012
【図６】
0003700804.tif 000013
【図７】
0003700804.tif 000014
【図８】
0003700804.tif 000015
【図９】
0003700804.tif 000016
【図１０】
0003700804.tif 000017
【図１１】
0003700804.tif 000018
【図１２】
0003700804.tif 000019
【図１３】
0003700804.tif 000020
【図１４】
0003700804.tif 000021
【図１５】
0003700804.tif 000022
【図１６】
0003700804.tif 000023