(57)【特許請求の範囲】
【請求項１】
原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成し、前記原画像信号および前記複数の非鮮鋭マスク画像信号に基づいて、または前記複数の非鮮鋭マスク画像信号に基づいて前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、該各帯域制限画像信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、該各変換画像信号を積算して得た積算信号を前記原画像信号に加算することにより該原画像信号の所定の周波数成分が強調された処理済画像信号を得る画像処理方法において、
前記処理済画像信号の所望の周波数応答特性を指定し、
該指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定して該変換関数を定義することを特徴とする画像処理方法。
【請求項２】
原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数の非鮮鋭マスク画像信号を作成する非鮮鋭マスク画像信号作成手段と、前記原画像信号および前記複数の非鮮鋭マスク画像信号に基づいて、または前記複数の非鮮鋭マスク画像信号に基づいて前記原画像信号の複数の周波数帯域ごとの信号を表す複数の帯域制限画像信号を作成し、該各帯域制限原画像信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、該各変換画像信号を積算して得た積算信号を前記原画像信号に加算することにより該原画像信号の所定の周波数成分が強調された処理済画像信号を得る非線形処理手段とからなる画像処理装置において、
前記処理済画像信号の所望の周波数応答特性を指定する周波数応答特性入力手段と、
前記指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定して該変換関数を定義する変換関数定義手段とを有することを特徴とする画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像信号の所定の周波数成分を強調するための画像処理方法および装置に関し、特に詳しくは所望の画像処理結果を得るための画像処理パラメータの自動設定に関する。
【０００２】
【従来の技術】
従来、本出願人により、非鮮鋭マスク画像信号（以下、ボケ画像信号という）を用いて周波数強調処理を行って放射線画像の診断性能を向上させる数々の画像処理方法および装置が提案されている（特開昭55-163472 号、同55-87953号など）。ここで、ボケ画像信号とは、画素数は原画像信号と同じであるが原画像信号よりも鮮鋭度が低い画像を表す画像信号であり、原画像信号の所定の周波数以上の高周波成分が除去された周波数応答特性を有する信号である。
【０００３】
上記周波数強調処理は、原画像信号Ｓorg からボケ画像信号Ｓusを引いたものに強調係数βを乗じたものを、原画像信号Ｓorgに加算することにより、原画像信号の所定の空間周波数成分を強調するものである。これを式で表すと下記の式（１）のようになる。
【０００４】
Ｓproc＝Ｓorg＋β×（Ｓorg−Ｓus）  …（１）
（Ｓproc：周波数強調処理された信号、Ｓorg ：原画像信号、Ｓus：ボケ画像信号、β：強調係数）
この際、上記処理では信号を加算することによりアーチファクトが発生してしまうことがあるが、これは、原画像信号Ｓorgに加算する加算信号の周波数応答特性を調整することにより解決でき、そのための具体的な調整方法として、次のような方法が提案されている。（特願平8-182155号など）。
【０００５】
この方法とは、まず鮮鋭度の異なる、すなわち周波数応答特性の異なる複数のボケ画像信号を作成し、そのボケ画像信号および原画像信号の中の２つの信号の差分をとることにより、原画像信号の、ある限られた周波数帯域の周波数成分を表す複数の帯域制限画像信号（以下、バンドパス信号という）を作成し、さらにそのバンドパス信号をそれぞれ異なる変換関数によって所望の大きさとなるように抑制してから、その複数の抑制されたバンドパス信号を積算することにより上記加算信号を作成するものである。この処理を式として表すと、例えば、下記の式（２）
0003696339.tif 000002
（但し、Ｓproc：高周波成分が強調された画像信号
Ｓorg ：原画像信号
Ｓusk(k=1〜N）：ボケ画像信号
ｆk(k=1〜N）：各バンドパス信号を変換する変換関数
β（Ｓorg）：原画像信号に基づいて定められる強調係数）
のようになる。
【０００６】
この方法によれば、バンドパス信号を変換する変換関数ｆ1〜ｆNの定義を変えることによって原画像信号Ｓorgに加算する加算信号の周波数応答特性を調整することができる。したがって、アーチファクトを防止するためには、変換関数ｆ1〜ｆNを、アーチファクトが発生しないような周波数応答特性を有する画像信号が作成されるような関数として定義すればよい。さらには、アーチファクトの防止という目的に限らず、各変換関数の定義次第で、所望の周波数応答特性を有する処理済画像信号を得ることができる。
【０００７】
【発明が解決しようとする課題】
しかしながら、上記方法に基づいて、実際に変換関数を定義しようとした場合、各変換関数をどのように定義すれば所望の結果が得られるかは、上記（２）のような式から簡単にわかるというものではない。そのため、従来は、強調処理された画像を観察しながら、変換関数を定義するパラメータ（例えば傾きなど）を少しずつ調整して、それぞれの変換関数を定義していた。これは、経験の少ないユーザにとっては非常に困難な作業であり、また熟練したユーザにとっても、時間のかかる作業であった。
【０００８】
本発明は上記問題に鑑みて、経験の少ないユーザでも簡単に、上記変換関数の定義パラメータを調整して、所望の周波数応答特性を有する処理済画像信号を得ることができる画像処理方法および装置を提供することを目的とする。
【０００９】
【課題を解決するための手段】
本発明の画像処理方法は、原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成し、原画像信号およびその複数のボケ画像信号に基づいて、あるいは複数のボケ画像信号に基づいて原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号を作成し、その各バンドパス信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、その各変換画像信号を積算して得た積算信号を原画像信号に加算することにより原画像信号の所定の周波数成分が強調された処理済画像信号を得る画像処理方法であって、前記処理済画像信号の所望の周波数応答特性を指定し、指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義することを特徴とするものである。
【００１０】
また、本発明の画像処理装置は、上記画像処理方法に基づいて画像処理を行う装置であって、原画像を表す原画像信号に基づいて互いに周波数応答特性が異なる複数のボケ画像信号を作成するボケ画像信号作成手段と、原画像信号およびその複数のボケ画像信号に基づいて、あるいは複数のボケ画像信号に基づいて原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号を作成し、その各バンドパス信号を複数の変換関数に基づいて変換して複数の変換画像信号を作成し、その各変換画像信号を積算して得た積算信号を原画像信号に加算することにより原画像信号の所定の周波数成分が強調された処理済画像信号を得る処理手段とからなる画像処理装置であり、前記処理済画像信号の所望の周波数応答特性を指定する周波数応答特性入力手段と、指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義する変換関数定義手段とを有することを特徴とするものである。
【００１１】
ここで、「ボケ画像信号」とは、画素数は原画像信号と同じであるが原画像信号よりも鮮鋭度が低い画像を表す画像信号である。ボケ画像信号は、まず原画像信号の画素に対して所定間隔ごとに所定のフィルタリング処理を施すことによって画素を間引きし、そのようにして得た画像信号に対して同様のフィルタリング処理を繰り返してさらに画素数を少なくした画像信号を複数作成し、そのそれぞれに対して、所定の補間方法により原画像と画素数が同じになるように補間処理を施すことにより作成する。ここで、フィルタリングや補間については、一般に広く使用されている種々の方法を適用することができる。
【００１２】
また、「原画像信号の複数の周波数帯域ごとの信号を表す複数のバンドパス信号」は、例えば隣接する周波数帯域のボケ画像信号同士で差分をとって作成してもよいし、原画像信号と各ボケ画像信号の差分をとって作成してもよい。あるいは原画像信号とボケ画像信号の他の組み合わせで差分をとって作成することもできる。なお、変換画像信号の作成、積算信号の作成、積算信号の原画像信号への加算は、例えば上記（２）式で表すことができる。
【００１３】
また、「所定の周波数成分を強調する」とは、例えば画像のエッジ部を強調するために、高周波成分を強調するといったことを意味する。
【００１４】
また、「処理済画像信号の所望の周波数応答特性を指定する」方法としては、例えば、操作画面上でユーザに対していくつかの周波数を示し、そのそれぞれに対して所望の周波数応答特性を数値として入力させてもよいし、あるいは周波数強調された画像信号の周波数応答特性のグラフを画面上に表示し、ユーザがマウス等のポインティングデバイスを用いてそのグラフを所望の形状となるように変形するようにしてもよい。
【００１５】
また、「指定された周波数応答特性に基づいて前記変換関数の定義パラメータを決定してその変換関数を定義する」とは、例えば、前記バンドパス信号、入力された所望の周波数応答特性、および変換関数の定義パラメータの関係式を、バンドパス信号と周波数応答特性を既知の値とし、変換関数の定義パラメータを変数とする連立方程式として解くことによりパラメータを求めるなどの方法を意味する。
【００１６】
【発明の効果】
本発明の画像処理方法および装置は、指定された周波数応答特性から、画像処理で用いられる変換関数の定義パラメータを逆算して求め、そのパラメータにより定義される変換関数を用いて画像処理を行うため、所望の周波数応答特性を指定するだけで、指定したとおりの周波数応答特性の画像信号を得る画像処理を行うことができ、パラメータ決定のための試行錯誤が不要となるばかりでなく、経験の少ないユーザでも簡単に所望の画像信号が得られるような画像処理を行うことができる。
【００１７】
【発明の実施の形態】
以下、本発明の画像処理方法および装置の一実施の形態について、図面を参照して詳細に説明する。以下に示す画像処理装置は、蓄積性蛍光体シートに記録された人体の放射線画像を読み取って得た画像信号に対して、その画像が診断に適した画像となるように、ボケ画像信号を使用して周波数強調処理を施すものであり、処理された画像信号は主としてフィルムに記録され、診断に用いられる。
【００１８】
図１はこの画像処理装置の概略を表す図である。画像処理装置１は、ボケ画像信号を作成するボケ画像信号作成手段２と、特定の周波数を強調するための周波数強調処理を行う非線形処理手段３を有する。さらに、この画像処理装置１は、変換関数定義手段４と、周波数応答特性入力手段５とを備える。変換関数定義手段４は、非線形処理手段３が変換処理に使用する変換関数を定義する手段であり、例えば関数の傾きなどのパラメータを決定して変換関数を定義する。周波数応答特性入力手段５は、強調処理された画像信号の所望の周波数応答特性をユーザが指定できるようにするための手段であり、変換関数定義手段４は周波数応答特性入力手段５によりユーザが指定した周波数応答特性に基づいて、パラメータを決定する。
【００１９】
ここで、まずボケ画像信号の作成処理について詳細に説明する。図２はボケ画像信号作成処理の概要を示すブロック図である。図２に示されるように、図１のボケ画像信号作成手段２は、まずフィルタリング処理手段10により、原画像信号Ｓorgに対し、原画像の画素のｘ方向およびｙ方向に対してフィルタリング処理を施して原画像信号よりも解像度が低い画像信号Ｂ1（以下、低解像度画像信号という）を作成し、次にこの低解像度画像信号Ｂ1に対して同様のフィルタリング処理を施してこの低解像度画像信号Ｂ1よりもさらに解像度が低い低解像度画像信号Ｂ2を作成し、以降順次同様のフィルタリング処理を重ねていくものである。そして、補間処理手段11により、このフィルタリング処理の各段において得られる低解像度画像信号Ｂkに対して、それぞれ補間拡大処理を施して、鮮鋭度の異なる複数のボケ画像信号Ｓus1〜ＳusNを得るものである。
【００２０】
本実施の形態においては、上記フィルタリング処理のフィルタとして、一次元ガウス分布に略対応したフィルタを使用する。すなわちフィルタのフィルタ係数を、ガウス信号に関する下記の式（３）
【００２１】
【数１】
0003696339.tif 000003
【００２２】
にしたがって定める。これは、ガウス信号は周波数空間および実空間の双方において、局在性がよいためであり、例えば上記（３）式においてσ＝１とした場合の５×１の１次元フィルタは図３に示すようなものとなる。
【００２３】
フィルタリング処理は、図４に示すように、原画像信号Ｓorg に対して、あるいは低解像度画像信号に対して１画素おきに行う。このような１画素おきのフィルタリング処理をｘ方向、ｙ方向に行うことにより、低解像度画像信号Ｂ1の画素数は原画像の１／４となり、フィルタリング処理により得られる低解像度画像信号に対して繰り返しこのフィルタリング処理を施すことにより、得られるｎ個の低解像度画像信号Ｂk（ｋ＝１〜ｎ）は、それぞれ、画素数が原画像信号の１／２2kの画像信号となる。
【００２４】
次に、このようにして得られた低解像度画像信号Ｂkに対して施される補間拡大処理について説明する。補間演算の方法としては、Ｂスプラインによる方法など種々の方法が挙げられるが、本実施の形態においては、上記フィルタリング処理においてガウス信号に基づくローパスフィルタを用いているため、補間演算についてもガウス信号を用いるものとする。具体的には、下記の式（４）
【００２５】
【数２】
0003696339.tif 000004
【００２６】
において、σ＝２k-1 と近似したものを用いる。
【００２７】
画像信号Ｂ1を補間する際には、ｋ＝１であるためσ＝１となる。この場合、補間処理を行うためのフィルタは、図５に示すように５×１の一次元フィルタとなる。この補間処理は、まず低解像度画像信号Ｂ1に対して１画素おきに値が０の画素を１つずつ補間することにより低解像度画像信号Ｂ1を原画像と同一のサイズに拡大し、次に、この補間された低解像度画像信号Ｂ1に対して上述した図５に示す一次元フィルタによりフィルタリング処理を施すことにより行われる。
【００２８】
同様に、この補間拡大処理を全ての低解像度画像信号Ｂkに対して行う。低解像度画像信号Ｂkを補間する際には、上記式（４）に基づいて、３×２k−１の長さのフィルタを作成し、画像信号Ｂkの各画素の間に値が０の画素を２k−１個ずつ補間することにより、原画像と同一サイズに拡大し、この値が０の画素が補間された画像信号Ｂkに対して３×２k−１の長さのフィルタにより、フィルタリング処理を施すことにより補間拡大する。
【００２９】
次に、上記のようにして作成されたボケ画像信号を用いて行われる非線形処理について説明する。以下に示す処理は、特定の周波数成分を強調するための処理であり、原画像信号およびボケ画像信号は、上記（２）式にしたがって処理されるものとする。
【００３０】
図６はこのような処理を行う周波数強調処理装置の一例を示す図である。この周波数強調処理装置では、減算器21として示されるように、原画像信号Ｓorgとボケ画像信号Ｓuskのうち隣接する２つの信号に基づいてバンドパス信号が作成される。図７は、このようにして作成されたバンドパス信号の周波数応答特性の一例を示す図であり、６個のバンドパス信号の特性が示されている。各バンドパス信号は、変換器22においてそれぞれ異なる変換関数f1〜fNにより所望の大きさとなるように抑制され、さらに上記（２）式にしたがって、その複数の抑制されたバンドパス信号が演算器23において積算され、さらに原画像信号に加算されて、処理済画像信号Ｓprocが生成される。
【００３１】
ここで、２種類の変換関数と、それにより得られる処理済画像信号の例を示す。図８は変換関数の一例を示す図であり、No.1が図６のバンドパス信号（Ｓorg−Ｓus1）を変換する変換関数、No.2が図６のバンドパス信号（Ｓus1−Ｓus2）を変換する変換関数、以下No.が大きい変換関数ほど低い周波数帯域のバンドパス信号を変換する関数を表している。この６つの変換関数により、図７に示される６個のバンドパス信号をそれぞれ変換し、前記（２）式にしたがって処理した結果得られる画像信号の入力値が小さな時の周波数特性は、図９に示されるようなものとなる。これに対し、変換関数として、図10に示されるような関数を用いた場合（但し図８同様、No.が大きい変換関数ほど低い周波数帯域のバンドパス信号を変換する関数を表している）、図７のバンドパス信号はそれぞれ図11に示されるような周波数応答特性の信号に変換され、最終的な処理済画像信号としては、図12に示されるな周波数応答特性の画像信号が得られる。
【００３２】
この例から明らかであるように、変換関数を調整することにより所望の異なる周波数応答特性の画像信号を得ることができる。しかし、図８と図９、そして図10と図12を比べてわかるように、図８の変換関数を図10の変換関数に変えたら、得られる画像信号の周波数応答特性が図９のようなものから図12のようなものになるということは、必ずしも容易にわかることではない。すなわち、従来は、試行錯誤しながら変換関数の傾きなどのパラメータを調整することによって、所望の周波数応答特性が得られるような関数を見つけなければならなかった。
【００３３】
本発明の画像処理方法および装置は、この問題を解決するものであるが、以下そのための手段である、図１の周波数応答特性指定手段５と変換関数定義手段４について説明する。ここで、周波数応答特性指定手段５は、ディスプレイなどの表示装置を介してユーザに周波数応答特性の入力を促して、その入力からユーザが所望する周波数応答特性を認識するものであり、具体的には、そのような処理を行うソフトウェア、入力装置、表示装置などを意味する。
【００３４】
本実施の形態では、図13に示されるような操作画面が装置に付属するディスプレイに表示される。このとき画面上には、その時点で設定されている変換関数に基づいて入力信号を処理した場合に得られる信号の周波数応答特性曲線が表示されている（図13の（１））。その特性曲線上には、移動可能な６個の指示点が表示されている。すなわち、その画像処理がｎ個のバンドパス信号を作成して非線形処理を行うものである場合、ｎ個の指示点が表示されるものとする。
【００３５】
ユーザは、各指示点を、例えばマウスなどのポインティングデバイスを用いて所望の周波数応答に対応する点に移動する（図13の（２））。但し、表示されているままの周波数応答でよければ、必ずしも移動する必要はない。移動操作後の点（移動されなかった点を含む）は、その点が示す周波数およびその周波数における所望の周波数応答として認識され、これに基づいて特性曲線が計算しなおされ、移動操作後の点を通るような新しい特性曲線が画面上に表示される（図13の（３））。ユーザは、確認画面において、新たに表示された特性曲線が所望のものであれば、画像処理を行うことを指示し、所望のものでなければ、再度指示点の移動操作を行うことができる（図13の（４））。
【００３６】
但し、周波数応答特性指定手段５の形態は上記実施の形態に限られるものではなく、例えばいくつかの周波数を順次画面に示し、ユーザがその周波数に対する所望の周波数応答特性を数値として順次入力するような形態など、他の様々な形態が考えられる。
【００３７】
次に、変換関数定義手段４により行われる処理について説明する。本実施の形態では、上記のように６つのバンドパス信号を用いて処理を行っているため、上記のようなユーザインタフェースにより入力された所望の周波数応答特性は、以下のような連立方程式（５）
a1 = p1×S11＋p2×S21＋p3×S31＋p4×S41＋p5×S51＋p6×S61
a2 = p1×S12＋p2×S22＋p3×S32＋p4×S42＋p5×S52＋p6×S62
a3 = p1×S13＋p2×S23＋p3×S33＋p4×S43＋p5×S53＋p6×S63
a4 = p1×S14＋p2×S24＋p3×S34＋p4×S44＋p5×S54＋p6×S64
a5 = p1×S15＋p2×S25＋p3×S35＋p4×S45＋p5×S55＋p6×S65
a6 = p1×S16＋p2×S26＋p3×S36＋p4×S46＋p5×S56＋p6×S66
p1〜p6 ≧ 0                                        …（５）
（但し、a1〜a6 ： 所望の値として指定された周波数応答
S11〜S66 ： 原画像信号に基づいて作成されたバンドパス信号の周波数応答特性
p1〜p6 ： 変換関数の０付近の傾きを表すパラメータ）
に、a1〜a6として当てはめらる。但し、方程式の次元は、バンドパス信号の数、すなわち指定された値の数に基づいて決まるものであり、上記連立方程式は一例にすぎない。
【００３８】
変換関数定義手段４は、上記（５）式のような連立方程式を、指定された周波数応答特性a1〜a6およびバンドパス信号S11〜S66を既定値とし、パラメータp1〜p6を変数として解くことによりこのパラメータを求め、これにより変換関数を定義するものである。なお、上記（５）式において、バンドパス信号の周波数応答特性S11〜S66として表されている値は、S11〜S16（S21〜S26,…S61〜S66も同様）が１つのバンドパス信号の所定の６つの周波数の応答特性に対応するものであり、これらはボケ画像信号作成において用いられるフィルタのフィルタ係数と、補間係数から、フーリエ変換など、周波数解析において通常用いられている変換方法により求めることができる。このための計算は、上記連立方程式を解く度に行ってもよいし、あらかじめ計算しておいてその結果をメモリに記憶しておいてもよい。
【００３９】
また、本実施の形態では、変換関数のパラメータp1〜p6は、それぞれ各バンドパス信号を変換する関数の０付近の傾きを表すものとしている。これは、実際に変換処理に使用される範囲が０付近であるからである。なお、その関数が非線形関数であっても、０付近の所定の範囲を線形関数として近似し、その傾きをパラメータとして上記方程式に当てはめればよい。ただし、これは一例にすぎず、変換関数を定義するためのパラメータは必ずしも傾きに限られるものではない。
【００４０】
以上のようにして定義された変換関数に基づいて、非線形処理手段３により画像信号を処理すれば、ユーザが指定したとおりの周波数応答特性を有する処理済画像信号を、直ちに得ることができる。つまり、本発明により、変換関数を定義するための従来のような試行錯誤が必要なくなり、経験の少ないユーザでも、画像処理装置の機能を十分に生かして所望の画像処理を行い、所望の処理済画像信号を得ることができる。これは、経験の少ないユーザばかりでなく、経験豊富なユーザにとっても作業時間の短縮という利益をもたらすものであり、実用上の効果は極めて大きい。
【図面の簡単な説明】
【図１】本発明の画像処理装置の概略を示す図
【図２】ボケ画像信号作成処理の概要を示すブロック図
【図３】フィルタリング処理に使用されるフィルタの一例を示す図
【図４】低解像度画像信号作成処理の詳細を示す図
【図５】補間拡大処理に使用されるフィルタの一例を示す図
【図６】周波数強調処理装置の一例を示す図
【図７】バンドパス信号の周波数応答特性の一例を示す図
【図８】変換関数の一例を示す図
【図９】図８の変換関数を変換に使用した場合に得られる画像信号の周波数応答特性を示す図
【図１０】変換関数の他の例を示す図
【図１１】図１０の変換関数により変換された図６のバンドパス信号を示す図
【図１２】図１０の変換関数を変換に使用した場合に得られる画像信号の周波数応答特性を示す図
【図１３】周波数応答特性を指定するためのユーザインタフェースの一例を示す図
【符号の説明】
１  画像処理装置
２  ボケ画像信号（非鮮鋭マスク画像信号）作成手段
３  非線形処理手段
４  変換関数定義手段
５  周波数応答特性入力手段
10  フィルタリング処理手段
11  補間処理手段
21  減算器
22  変換器
23  演算器
【図１】
0003696339.tif 000005
【図２】
0003696339.tif 000006
【図３】
0003696339.tif 000007
【図４】
0003696339.tif 000008
【図５】
0003696339.tif 000009
【図６】
0003696339.tif 000010
【図７】
0003696339.tif 000011
【図８】
0003696339.tif 000012
【図９】
0003696339.tif 000013
【図１０】
0003696339.tif 000014
【図１１】
0003696339.tif 000015
【図１２】
0003696339.tif 000016
【図１３】
0003696339.tif 000017