JPH0944645A - 画像処理方法および装置 - Google Patents

Abstract

(57)【要約】 【課題】 グリッド像を含む画像に対して多重解像度変
換を施して画像処理を行う画像処理方法において、干渉
縞に基づくモアレが発生することのない画像を得る。 【解決手段】 画像入力手段31から入力された、グリッ
ド像を含む画像を表す画像信号Ｓを多重解像度分解処理
手段32においてラプラシアンピラミッド等の方法により
多重解像度の画像に分解する。この際、原画像に対して
フィルタリング処理を施すためのフィルタを、グリッド
像を除去する特性を有するフィルタとする。分解された
画像の所望とする周波数帯域の強調度をより低周波数帯
域の画像の信号値に基づいて決定し、強調処理手段33に
おいて画像に強調処理を施す。強調処理が施された画像
とそれ以外の画像とを復元処理手段34において復元し、
処理済画像信号Ｓ′を得る。処理済画像信号Ｓ′は画像
出力手段35において可視像として再生される。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、原画像の所定の周
波数帯域の画像に画像処理を施す画像処理方法および装
置に関し、特にグリッドを使用して撮影を行うことによ
り得られるグリッド像が記録された画像に対して処理を
行う画像処理方法および装置に関するものである。

【０００２】

【従来の技術】画像を表す画像信号を得、この画像信号
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば放射線画像の診断
性能を向上させるために、画像信号に対してボケマスク
処理等の周波数強調処理を施す方法が本出願人により提
案されている（特開昭55-163772 等）。この周波数処理
は、原画像を表す画像信号からボケマスク信号を減算し
たものに強調度を掛けたものを加える処理を施すもの
で、これにより画像において所定の空間周波数成分を強
調するようにしたものである。

【０００３】また、画像信号に対して周波数処理を施す
別の方法として、フーリエ変換、ウェーブレット変換、
サブバンド変換等により画像を多重解像度画像に変換す
ることにより画像を表す画像信号を複数の周波数帯域の
信号に分解し、この分解された信号のうち、所望とする
周波数帯域の信号に対して強調等の所定の画像処理を施
す方法が提案されている。

【０００４】また、近年画像処理の分野において、画像
を多重解像度に変換する新規な方法としてラプラシアン
ピラミッドなる方法が提案されている（例えば特開平6-
301766号）。このラプラシアンピラミッドは、原画像に
対してガウス関数で近似されたようなマスクによりマス
ク処理を施した後、画像をサブサンプリングして画素数
を間引いて半分にすることにより、原画像の１／４のサ
イズのボケ画像を得、このボケ画像のサンプリングされ
た画素に値が０の画素を補間して元の大きさの画像に戻
し、この画像に対してさらに上述したマスクによりマス
ク処理を施してボケ画像を得、このボケ画像を原画像か
ら減算して原画像の所定の周波数帯域の細部画像を得る
ものである。この処理を得られたボケ画像に対して繰り
返すことにより原画像の１／２^2Nの大きさのボケ画像を
Ｎ個作成するものである。ここで、ガウス関数で近似さ
れたようなマスクによりマスク処理を施した画像に対し
てサンプリングを行っているため、実際にはガウシアン
フィルタを用いているが、ラプラシアンフィルタをかけ
た場合と同様の処理画像が得られる。そしてこのように
原画像サイズの画像から順に１／２^2Nの大きさの低周波
数帯域の画像が得られるため、この処理の結果得られた
画像はラプラシアンピラミッドと呼ばれる。

【０００５】なお、このラプラシアンピラミッドについ
ては、Burt P.J.,“Fast Filter Transforms for Image
Processing ”，Computer Graphics and Image Proces
sing16 巻、20〜51頁、1981年；Crowley J.L.,Stern R.
M.,“Fast Computation ofthe Difference of Low・Pass
Transform”IEEETrans.on Pattern Analysis and Mach
ine Intelligence、６巻、２号、1984年３月、Mallat
S.G.,“A Theory for Multiresolution Signal Decompo
sition ；The Wavelet Representation”IEEE Trans.on
Pattern Analysis and Machine Intelligence 、11
巻、７号、1989年７月；Ebrahimi T.,Kunt M.,“Image
compression by Gabor Expansion”，Optical Engineer
ing,30巻、７号、873 〜880 頁、1991年７月、およびPi
eter Vuylsteke,Emile Schoeters，“Multiscale Image
Contrast Amplification ”SPIEVol.2167 Image Proce
ssing(1994),pp551 〜560 に詳細が記載されている。

【０００６】そしてこのようにして得られたラプラシア
ンピラミッドの全ての周波数帯域の画像に対して、画像
の値を強調するような処理を施し、この強調処理が施さ
れた各周波数帯域の画像を逆変換して処理が施された画
像を得る方法が上記特開平6-301766号に記載されてい
る。このように処理が施された画像は、各周波数帯域に
おいて画像が強調されているため、実質的に上述したボ
ケマスク処理において複数のサイズのマスクによりボケ
マスク処理を施したような画像となっている。

【０００７】一方、放射線画像をＸ線フイルムや蓄積性
蛍光体シート等の記録シートに撮影記録することが行わ
れているが、この撮影記録の際に被写体で散乱された放
射線が記録シートに照射されないように4.0 本／mm程度
の細かなピッチで放射線の透過しないたとえば鉛と透過
しやすいたとえばアルミニウムや木材等が交互に配置さ
れたグリッドを被写体と記録シートとの間に配置して撮
影を行うことがある。グリッドを用いて撮影を行うと被
写体により散乱された放射線が記録シートに照射されに
くいため、被写体の放射線画像のコントラストを向上さ
せることができるが、一方記録シートには被写体像とと
もに細かな縞模様状のグリッド像が記録される。

【０００８】

【発明が解決しようとする課題】ところで、放射線画像
が記録された記録シートから画像信号を得る放射線画像
読取装置は、通常放射線画像情報を担持した光を光電的
に読み取り、画像情報として必要な最高の空間周波数
（以下、この必要な最高の空間周波数をｆssで表す。）
に対応したサンプリング間隔Δｘ＝１／（２・ｆss）で
サンプリングしてデジタルの画像信号を得るように構成
される。このようにして得られた画像信号には、被写体
の放射線画像を表す有用な情報だけでなく、たとえグリ
ッド像の空間周波数が必要な最高の空間周波数ｆssより
高周波であってもグリッド像に起因するノイズも含まれ
ることがある。

【０００９】図14は記録シートに記録された放射線画像
（グリッド像を含む）の、該グリッド像の縞模様と直交
する方向の空間周波数特性を表すグラフを表す図であ
る。

【００１０】ここでは4.0 本／mmのグリッドを用いて撮
影を行い、グリッド像の空間周波数が４cycle ／mmにあ
るものとし、被写体の放射線画像の再現に必要な最高の
空間周波数ｆssは2.5cycle／mmであるとして説明する。

【００１１】図15は被写体の放射線画像の再現に必要な
ｆss＝2.5(cycle ／mm）以下の空間周波数帯の情報を得
るために、ｆss＝2.5(cycle ／mm）に対応したサンプリ
ング間隔Δｘ＝１／（２・ｆss）＝0.2 （mm）、すなわ
ち１mmにつき５回サンプリングを行った場合のノイズの
混入の状態を説明するための図である。

【００１２】この場合図に実線で示すグラフ（図14のグ
ラフと同一）がｆss＝2.5(cycle ／mm) で折り返した位
置にいわゆるエリアジングと呼ばれるノイズが混入す
る。したがってこの場合グリッド像の空間周波数４cycl
e ／mmのエリアジングは１cycle ／mmに生じる。

【００１３】図16は、上記サンプリング間隔Δｘ＝１／
（２・ｆss）＝0.2 （mm）でサンプリングして得られた
画像信号が担持する放射線画像の空間周波数特性を示し
ている。

【００１４】１cycle ／mmの位置に上記グリッド像に対
応したノイズが混入し、可視画像を再生記録した際、該
可視画像に１cycle ／mmの縞模様として現われる結果と
なる。放射線画像を観察するに際し、グリッド像がそれ
ほど目立たない空間周波数帯にある場合であっても、サ
ンプリングして画像信号を得ると目立ちやすい空間周波
数帯に縞模様が現われる。そしてこの縞模様により干渉
縞が生じ、モアレが生じてしまう。そしてこのモアレに
よりこの画像データに基づいて再生した可視画像が非常
に見にくいものとなってしまうことがあった。

【００１５】また、上述した特開平6-301766号に記載さ
れた方法において、グリッド像を含む画像を多重解像度
に変換すると、画像を多重解像度に変換するためのフィ
ルタとしてガウス関数を用いているため、画像が低周波
数帯域のものとなるにしたがって上述したエリアジング
が繰り返し起こり、モアレが非常に目立つものとなって
しまう。そしてこのようなモアレが含まれる画像に対し
て画像処理を施し、画像処理が施された画像を逆多重解
像度変換することにより復元しても、復元された処理済
画像はモアレによるアーチファクトが多く含まれる見に
くい画像となってしまう。

【００１６】本発明は上記事情に鑑み、グリッド像を含
む画像に対して多重解像度変換を施して画像処理を行う
画像処理方法において、干渉縞に基づくモアレが発生す
ることのない画像を得ることができる画像処理方法およ
び装置を提供することを目的とするものである。

【００１７】

【課題を解決するための手段】本発明による画像処理方
法および装置は、グリッド像が記録された原画像に対し
て、所定のフィルタによるフィルタリング処理を繰り返
し施すことにより該原画像を多重解像度空間に変換して
得られた複数の周波数帯域ごとの画像のうち所望の周波
数帯域の画像に対して画像処理を施し、該画像処理が施
された周波数帯域の画像および他の周波数帯域の画像を
逆多重解像度空間に変換することにより処理済画像を得
る画像処理方法および装置において、少なくとも最初に
前記フィルタリング処理を施す際の前記所定のフィルタ
を、前記グリッド像を除去するグリッド像除去フィルタ
としたことを特徴とするものである。

【００１８】

【発明の効果】本発明による画像処理方法および装置
は、原画像をフィルタリング処理して多重解像度に変換
する際に、少なくとも最初にフィルタリング処理を施す
フィルタをグリッドを除去する特性を有するフィルタと
したため、画像を多重解像度複数に変換しても、変換さ
れた画像が全ての周波数帯域においてグリッド像が除去
されたものとなり、したがって、所望とする周波数帯域
の画像に対して画像処理を施して逆多重解像度変換を施
すことにより得られる処理済画像信号をグリッド像が消
去された、干渉縞によるモアレのない見易い画像とする
ことができる。

【００１９】

【発明の実施の形態】以下図面を参照して本発明の実施
の形態について説明する。

【００２０】図１は、放射線画像撮影装置の一例の概略
を示した図である。ここでは記録シートとして前述した
蓄積性蛍光体シートを用いる例について説明する。

【００２１】放射線源１から放射された放射線２は、被
写体３を経由しさらにグリッド４を経由して、蓄積性蛍
光体シート11を照射する。グリッド４は４本／mmのピッ
チで鉛4aとアルミニウム4bとが交互に配置されている。
放射線２は、鉛4aには遮ぎられ、アルミニウム4bは透過
してシート11を照射する。このためシート11には被写体
像とともに４本／mmの縞模様状のグリッド像が蓄積記録
される。被写体３内で散乱された放射線2aはグリッド４
に斜めに入射するため該グリッドに遮られ、または該グ
リッド４により反射され、シート11には照射されず、し
たがってシート11には散乱放射線の照射の少ない鮮明な
放射線画像が蓄積記録される。

【００２２】図２は、グリッドを使用して撮影を行うこ
とにより蓄積性蛍光体シート11に蓄積記録された、被写
体像（図の斜線部）と縞模様状のグリッド像（図の縦
縞）を示した図である。このようにシート11には被写体
像５と縞模様６とが重畳された放射線画像が記録され
る。

【００２３】図３は、本発明の画像処理方法の一実施の
形態を使用した、放射線画像読取装置の一実施例の斜視
図である。

【００２４】読取手段の一例である読取部10の所定位置
にセットされた放射線画像が記録された蓄積性蛍光体シ
ート11は、図示しない駆動手段により駆動されるエンド
レスベルト等のシート搬送手段15により、矢印Ｙ方向に
搬送（副走査）される。一方、レーザー光源16から発せ
られた光ビーム17はモータ24により駆動され矢印方向に
高速回転する回転多面鏡18によって反射偏向され、ｆθ
レンズ等の集束レンズ19を通過した後、ミラー20により
光路を変えて前記シート11に入射し副走査の方向（矢印
Ｙ方向）と略垂直な矢印Ｘ方向に主走査する。光ビーム
17が照射されたシート11の箇所からは、蓄積記録されて
いる放射線画像情報に応じた光量の輝尽発光光21が発散
され、この輝尽発光光21は光ガイド22によって導かれ、
フォトマルチプライヤ（光電子増倍管）23によって光電
的に検出される。上記光ガイド22はアクリル板等の導光
性材料を成形して作られたものであり、直線状をなす入
射端面22a が蓄積性蛍光体シート11上の主走査線に沿っ
て延びるように配され、円環状に形成された射出端面22
b にフォトマルチプライヤ23の受光面が結合されてい
る。入射端面22a から光ガイド22内に入射した輝尽発光
光21は、該光ガイド22の内部を全反射を繰り返して進
み、射出端面22b から射出してフォトマルチプライヤ23
に受光され、放射線画像を表す輝尽発光光21がフォトマ
ルチプライヤ23によって電気信号に変換される。

【００２５】フォトマルチプライヤ23から出力されたア
ナログ出力信号Ｓ_A には、良好な放射線可視画像を再生
出力するために必要な、所望とする範囲の空間周波数帯
のうちの最高の第１の空間周波数ｆss＝2.5cycle／mmよ
り高い空間周波数帯域の情報、特に図２に示すグリッド
像６に関する情報も含まれている。このグリッド像６に
関する情報は可視画像を観察する際にその可視画像を見
にくくする原因のひとつとなるものであり、取り除く必
要がある情報である。

【００２６】アナログ出力信号Ｓ_A はログアンプ26で対
数的に増幅された後、Ａ／Ｄ変換器28において所定のサ
ンプリング間隔でサンプリングされてディジタル化さ
れ、ディジタルの画像信号Ｓが得られる。この画像信号
Ｓは一旦記憶部29に記憶される。

【００２７】図４は、蓄積性蛍光体シート11上の各サン
プリング点を示した図である。図の横方向（ｘ方向）は
図３のＸ方向（主走査方向）に対応し、縦方向（ｙ方
向）は図３のＹ方向（副走査方向）に対応している。図
の・印が初期画像信号における各サンプリング点、○印
が再サンプリング後の各サンプリング点を表している。
本実施の形態においては、この各サンプリング点はｘ方
向（主走査方向）およびｙ方向（副走査方向）双方と
も、蓄積性蛍光体シート11上に記録された放射線画像を
読み取るために必要な所望とする空間周波数帯の最高の
空間周波数ｆss＝2.5cycle／mmに対応するサンプリング
間隔２・Δ＝１／（２・ｆss）＝0.2 （mm）の半分のサ
ンプリング間隔Δ＝0.1 （mm）、すなわち上記最高の空
間周波数ｆss＝2.5(cycle ／mm）の２倍の、空間周波数
ｆsw＝5.0(cycle ／mm）に対応するサンプリング間隔Δ
＝１／（２・ｆsw）＝0.1 （mm）でサンプリングされ
る。このサンプリング間隔の調整は、ｘ方向（主走査方
向）についてはＡ／Ｄ変換器28（図３参照）のサンプリ
ング時間間隔を調整すること等により行われ、ｙ方向
（副走査方向）については、シート搬送手段15（図３参
照）によるシート11の搬送速度を調整すること等により
行われる。

【００２８】このようにして得られた図４の・点で示す
サンプリング間隔（0.1 mm間隔）の画像信号Ｓは、ｆsw
＝5.0(cycle ／mm）以下の情報を担持しており、したが
って図２に示す縞模様６の情報(4.0cycle ／mm）も担持
している。またこの実施の形態においては縞模様６のエ
リアジングは発生しない。

【００２９】この画像信号Ｓは、図３に示す記憶部29に
一旦記憶された後、画像信号処理手段30に入力され、以
下のようにして画像処理が施される。

【００３０】図５は本発明による画像処理方法を実施す
るための装置の概略を表すブロック図である。図５に示
すように本発明による画像処理方法を実施するための装
置は、装置に画像を入力するための上述した図３に示す
画像読取手段のような画像入力手段31と、入力された画
像に対して多重解像度分解処理を施す多重解像度分解処
理手段32と、多重解像度分解処理手段32において複数の
周波数帯域に分解された画像のうち、所定の周波数帯域
の画像に対して後述するような強調処理を施す強調処理
手段33と、強調処理手段33により強調処理が施された周
波数帯域の画像および他の周波数帯域の画像を復元して
処理済画像を得るための復元処理手段34と、復元処理手
段34により復元された処理済画像を可視像として再生す
るための画像出力手段35とからなるものである。

【００３１】次いで本発明による画像処理方法の作用に
ついて説明する。図６は図５における多重解像度画像分
解処理手段32において行われる処理を説明するためのブ
ロック図である。なお、本実施の形態においてはラプラ
シアンピラミッドの手法により画像信号Ｓを多重解像度
画像に分解するものとする。図３に示すように原画像を
表すデジタルの画像信号Ｓが多重解像度分解処理手段32
に入力されると、まず、以下のようなフィルタリング処
理が施される。

【００３２】すなわち、図６に示すフィルタリング手段
40において、図７に示すようなフィルタＬ_G により、画
像信号Ｓに対してフィルタリング処理を施す。すなわ
ち、図４に示す画素値ａij（i,j ＝1,2 …）においてフ
ィルタリング処理後の画素値ａ′ijが

【００３３】

【数１】

【００３４】となるようにフィルタリング処理が施され
る。ここで、図７に示すフィルタＬ_Gは、図８に示すよ
うに４cycle ／mmの伝達特性が０となる特性を有するも
のであり、したがって図７に示すフィルタＬ_G によりフ
ィルタリング処理を施すことにより、画像信号Ｓは縞模
様６の情報を担持しないものとなる。なお、このフィル
タＬ_G は、特開平3-114039号に開示されているものであ
る。

【００３５】このようなフィルタＬ_G によりフィルタリ
ングされた画像信号Ｓはフィルタリング手段40において
１画素おきにサンプリングされ、縞模様除去画像ｇ₁ が
得られる。この縞模様除去画像ｇ₁ は、原画像の１／４
の大きさとなっている。次いで補間手段41において、こ
の縞模様除去画像ｇ₁ のサンプリングされた間隔に値が
０の画素が補間される。この補間は縞模様除去画像ｇ₁
の一列毎および一行毎に値が０の行および列を挿入する
ことにより行う。このように値が０の画素が補間された
縞模様除去画像ｇ₁ は一画素おきに値が０の画素が挿入
されているため、信号値の変化が滑らかではないものと
なっている。

【００３６】そしてこのようにして補間が行われた後、
さらにこの補間がなされた縞模様除去画像ｇ₁ に対して
フィルタＬにより再度フィルタリング処理を施し、低解
像度近似画像ｇ₁ ′を得る。この低解像度近似画像
ｇ₁ ′は上述した補間がなされた縞模様除去画像ｇ₁ と
比較して信号値の変化が滑らかなものとなっている。ま
た原画像と比較して周波数帯域的には半分より高い高周
波数が消えたような画像となっている。これは画像の大
きさを１／４にして一画素おきに値が０の画素を補間
し、さらにフィルタＬによりフィルタリング処理を施し
ているため、空間周波数が半分よりも高い周波数帯域の
画像がぼかされたようになっているからである。

【００３７】次いで減算器42において、原画像から低解
像度近似画像ｇ₁ ′の減算が行われ、細部画像ｂ₀ が得
られる。この減算は原画像と低解像度近似画像ｇ ₁′と
の相対応する画素についての信号間で行われる。ここ
で、低解像度近似画像ｇ₁ ′は上述したように原画像の
空間周波数のうち半分より高い周波数帯域の画像がぼけ
たようになっているため、細部画像ｂ₀ は原画像のうち
半分より上の周波数帯域のみを表す画像となっている。
すなわち、図９に示すように細部画像ｂ₀ は原画像のナ
イキスト周波数ＮのうちＮ／２〜Ｎの周波数帯域の画像
を表すものとなっている。

【００３８】次いで、縞模様除去画像ｇ₁ はフィルタリ
ング手段40に入力され、フィルタリング手段40において
ローパスフィルタＬによりフィルタリングされる。この
ローパスフィルタは例えば図10に示すように５×５のグ
リッド上の二次元ガウス分布に略対応している。このロ
ーパスフィルタＬは以降の全ての解像度の画像に対して
適用される。

【００３９】このようなローパスフィルタによりフィル
タリングされた縞模様除去画像ｇ₁は、フィルタリング
手段40において１画素おきにサンプリングされ、低解像
度近似画像ｇ₂ が得られる。この低解像度近似画像ｇ₂
は、縞模様除去画像ｇ₁ の１／４すなわち原画像の１／
16の大きさとなっている。次いで補間手段41において、
この低解像度近似画像ｇ₂ のサンプリングされた間隔に
値が０の画素が補間される。この補間は低解像度近似画
像ｇ₂ の一列毎および一行毎に値が０の行および列を挿
入することにより行う。このように値が０の画素が補間
された低解像度近似画像ｇ₂ はぼけてはいるものの一画
素おきに値が０の画素が挿入されているため、信号値の
変化が滑らかではないものとなっている。

【００４０】そしてこのようにして補間が行われた後、
さらにこの補間がなされた低解像度近似画像ｇ₂ に対し
て図10に示すローパスフィルタＬにより再度フィルタリ
ング処理を施し、低解像度近似画像ｇ₂ ′を得る。この
低解像度近似画像ｇ₂ ′は上述した補間がなされた低解
像度近似画像ｇ₂ と比較して信号値の変化が滑らかなも
のとなっている。また縞模様除去画像ｇ₁ と比較して周
波数帯域的には半分より高い周波数帯域の画像が消えた
ようになっている。

【００４１】次いで減算器42において、縞模様除去画像
ｇ₁ から低解像度近似画像ｇ₂ ′の減算が行われ、細部
画像ｂ₁ が得られる。この減算は縞模様除去画像ｇ₁ と
低解像度近似画像ｇ₂ ′との相対応する画素についての
信号間で行われる。ここで、低解像度近似画像ｇ₂ ′は
上述したように縞模様除去画像ｇ₁ の空間周波数のうち
半分より高い周波数帯域の画像がぼけたようになってい
るため、細部画像ｂ₁は縞模様除去画像ｇ₁ のうち半分
より上の周波数帯域のみを表す画像となっている。すな
わち、図９に示すように細部画像ｂ₁ は縞模様除去画像
ｇ₁ のうちの半分より上の周波数帯域のみ、すなわち原
画像のナイキスト周波数ＮのうちＮ／４〜Ｎ／２の周波
数帯域の画像を表すものとなっている。このようにガウ
ス分布のローパスフィルタＬによりフィルタリング処理
を施して細部画像を得るようにしているが、フィルタリ
ング処理が施された画像を低解像度近似画像から減算し
ていることから、実質的にはラプラシアンフィルタによ
りフィルタリング処理を施した場合と同様の結果とな
る。

【００４２】そして上述した処理をフィルタリング手段
40によりフィルタリングされかつサンプリングされた低
解像度近似画像ｇ_k （ｋ＝１〜Ｎ）に対して順次繰り返
し行い、図９に示すようにｎ個の細部画像ｂ_k （ｋ＝１
〜ｎ）および低解像度近似画像の残留画像ｇ_L を得る。
ここで、細部画像ｂ_k は、ｂ₀ から順に解像度が低くな
る、すなわち画像の空間周波数が低くなるものであり、
原画像のナイキスト周波数Ｎに対して、細部画像ｂ_k は
Ｎ／２^k+1 〜Ｎ／２^k の周波数帯域を表し、画像の大き
さが原画像の１／２^2k倍となっている。すなわち、最も
解像度が高い細部画像ｂ₀ は原画像と同じ大きさである
が、細部画像ｂ₀ の次に高解像度の細部画像ｂ₁ 原画像
の大きさの１／４となっている。このように、細部画像
が原画像と同一の大きさのものから順次小さくなり、ま
た細部画像はラプラシアンフィルタを施したものと実質
的に同一の画像であることから、本実施例による多重解
像度変換はラプラシアンピラミッドと呼ばれるものであ
る。また、残留画像ｇ_L は原画像の非常に解像度が低い
近似画像であると見なすことができ、極端な場合は、残
留画像ｇ_L は原画像の平均値を表す１つだけの画像から
なるものとなる。そしてこのようにして得られた細部画
像ｂ_k および残留画像ｇ_L は図示しないメモリに記憶さ
れる。

【００４３】次いでこのようにして得られた細部画像ｂ
_k のうち、所望とする周波数帯域の細部画像ｂ_k に対し
て強調処理手段33において強調処理が施される。この強
調処理は、例えば、ある周波数帯域の細部画像ｂ_k に対
してある強調係数を乗じることにより行われる。

【００４４】次いで、強調処理が施された周波数帯域の
細部画像ｂ_k および他の周波数帯域の細部画像を逆変換
する。この逆変換の処理は復元処理手段34において以下
のようにして行われる。

【００４５】図11は細部画像の逆変換の詳細を表す図で
ある。まず、残留画像ｇ_L が補間手段44により各画素の
間が補間されて元の大きさの４倍の大きさの画像ｇ_L ′
とされる。次に加算器45においてその補間された画像ｇ
_L ′と最も低解像度の細部画像ｂ_n-1 の相対応する画素
同志で加算を行い、加算画像（ｇ_L ′＋ｂ_n-1 ）を得
る。次いでこの加算画像（ｇ_L ′＋ｂ_n-1 ）は補間手段
44に入力され、この補間手段44において各画素の間が補
間されて元の大きさの４倍の大きさの画像ｂ_n-1´とさ
れる。

【００４６】次いでこの画像ｂ_n-1 ′は、加算器45にお
いて細部画像ｂ_n-1の一段階高解像度の画像ｂ_n-2 と相
対応する画素同志の加算が行われ、加算された加算信号
（ｂ_n-1 ′＋ｂ_n-2 ）は補間手段44において各画素の間
隔が補間され、細部画像ｂ_n-2 の４倍の大きさの画像ｂ
_n-2 とされる。

【００４７】以上の処理を繰り返し、強調画像ｂ_kpにつ
いても同様の処理を施す。すなわち、強調画像ｂ_kpと上
述した処理が施された一段階低解像度の画像ｂ_k-1 ′と
の加算が加算器45において行われ、さらに加算信号（ｂ
_kp＋ｂ_k-1 ′）に対して補間手段44において各画素の間
が補間され、補間信号ｂ_kp′を得る。そしてこの処理を
より高周波の細部画像に対して順次行い、最終的に加算
器45において補間画像ｂ₁ ′と最高解像度の細部画像ｂ
₀ との加算が行われ、処理済画像信号Ｓ′を得る。

【００４８】このようにして得られた処理済画像信号
Ｓ′は画像出力手段35に入力され、可視像として表示さ
れる。この画像出力手段５はＣＲＴ等のディスプレイ手
段でもよいし、感光フィルムに光走査記録を行う記録装
置であってもよいし、あるいはそのために画像信号を一
旦光ディスク、磁気ディスク等の画像ファイルに記憶さ
せる装置であってもよい。

【００４９】なお、逆変換フィルタには画像ｇ_n ′を求
めるために用いたローパスフィルタ（図６の補間手段41
で用いられたフィルタ）と同じものを用いることにより
原画像を完全に復元することができる。また上記実施の
形態以外にも、変換時に画像ｇ₁ ′を求めるフィルタと
して、フィルタＬ_G を用いた場合、補間画像ｂ₁ ′を求
めるフィルタとしてはＬ_G を用いればよい。

【００５０】また、上述した実施の形態においては、処
理後の信号がグリッドを除去するためのフィルタとして
図８に示すような周波数特性を有するフィルタＬ_G （図
７参照）を用いるようにしているが、これに限定される
ものではなく、例えば図13に示すような周波数特性を有
するフィルタＬ_G ′（図12参照）を用いるようにしても
よい。

【図面の簡単な説明】

【図１】グリッドを用いて放射線画像の撮影を行う撮影
装置を表す図

【図２】グリッドを使用することにより得られる放射線
画像を表す図

【図３】放射線画像読取装置を表す図

【図４】蓄積性蛍光体シート上のサンプリング点を表す
図

【図５】本発明による画像処理方法を説明するための図

【図６】多重解像度分解処理手段の詳細を表す図

【図７】グリッド像除去フィルタを表す図

【図８】グリッド像除去フィルタの周波数特性を表す図

【図９】ラプラシアンピラミッド処理が施された複数の
周波数帯域ごとの細部画像を表す図

【図１０】ローパスフィルタを表す図

【図１１】復元処理手段の詳細を表す図

【図１２】他のグリッド像除去フィルタを表す図

【図１３】他のグリッド像除去フィルタの周波数特性を
表す図

【図１４】グリッド像の縞模様と直交する方向の空間周
波数特性を表す図

【図１５】図１４のグラフと2.5cycle／mmに対応するサ
ンプリング間隔でサンプリングされた画像信号のエリア
ジングを表す図

【図１６】図１５と同一のサンプリング間隔でサンプリ
ングされた画像信号が担持する放射線画像の空間周波数
特性を表す図

【符号の説明】

11 蓄積性蛍光体シート 17 光ビーム 23 フォトマルチプライヤ 26 ログアンプ 27 アナログフィルタ 28 Ａ／Ｄ変換器 29 記憶部 30 画像処理手段 31 画像入力手段 32 多重解像度分解処理手段 33 強調処理手段 34 復元処理手段 35 画像出力手段 40 フィルタリング手段 41 補間手段 42 減算器 44 補間手段 45 加算器

Claims (2)

【特許請求の範囲】
1. 【請求項１】 グリッドを使用して撮影を行うことによ
   り得られるグリッド像が記録された原画像に対して、所
   定のフィルタによるフィルタリング処理を繰り返し施す
   ことにより該原画像を多重解像度空間に変換して複数の
   周波数帯域ごとの画像を得、 該複数の周波数帯域ごとの画像のうち所望の周波数帯域
   の画像に対して画像処理を施し、 該画像処理が施された周波数帯域の画像および他の周波
   数帯域の画像を逆多重解像度空間に変換することにより
   処理済画像を得る画像処理方法であって、 前記原画像に対して少なくとも最初に前記フィルタリン
   グ処理を施す際の前記所定のフィルタが、前記グリッド
   像を除去するグリッド像除去フィルタであることを特徴
   とする画像処理方法。
2. 【請求項２】 グリッドを使用して撮影を行うことによ
   り得られるグリッド像が記録された原画像に対して、所
   定のフィルタによるフィルタリング処理を繰り返し施す
   ことにより該原画像を多重解像度空間に変換して複数の
   周波数帯域ごとの画像を得る多重解像度変換手段と、 該複数の周波数帯域ごとの画像のうち所望の周波数帯域
   の画像に対して画像処理を施す画像処理手段と、 該画像処理が施された周波数帯域の画像および他の周波
   数帯域の画像を逆多重解像度空間に変換することにより
   処理済画像を得る逆変換手段とを備えた画像処理装置で
   あって、 前記原画像に対して少なくとも最初に前記フィルタリン
   グ処理を施す際の前記所定のフィルタが、前記グリッド
   像を除去するグリッド像除去フィルタであることを特徴
   とする画像処理装置。