JPH0944659A - エネルギーサブトラクション処理方法および装置 - Google Patents
エネルギーサブトラクション処理方法および装置Info
- Publication number
- JPH0944659A JPH0944659A JP7191524A JP19152495A JPH0944659A JP H0944659 A JPH0944659 A JP H0944659A JP 7191524 A JP7191524 A JP 7191524A JP 19152495 A JP19152495 A JP 19152495A JP H0944659 A JPH0944659 A JP H0944659A
- Authority
- JP
- Japan
- Prior art keywords
- image
- signal
- image signal
- frequency
- processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】
【課題】 エネルギーサブトラクション処理方法および
装置において、得られるサブトラクション画像の画質を
向上する。 【解決手段】 画像入力手段１から入力された２つの画
像信号ＳＯ1 ，ＳＯ2を多重解像度分解処理手段２にお
いてラプラシアンピラミッド等の方法により複数の周波
数帯域ごとの細部画像信号に分解し、記憶手段３に記憶
する。エネルギーサブトラクション処理手段４が記憶手
段３から各周波数帯域ごとの細部画像信号を読み出し、
各周波数帯域に応じて重み付け係数を変更してサブトラ
クション処理を行い差信号を得る。差信号は復元処理手
段５において復元され、エネルギーサブトラクション信
号Ｓ′を得る。エネルギーサブトラクション信号Ｓ′は
画像出力手段６において可視像と して再生される。
装置において、得られるサブトラクション画像の画質を
向上する。 【解決手段】 画像入力手段１から入力された２つの画
像信号ＳＯ1 ，ＳＯ2を多重解像度分解処理手段２にお
いてラプラシアンピラミッド等の方法により複数の周波
数帯域ごとの細部画像信号に分解し、記憶手段３に記憶
する。エネルギーサブトラクション処理手段４が記憶手
段３から各周波数帯域ごとの細部画像信号を読み出し、
各周波数帯域に応じて重み付け係数を変更してサブトラ
クション処理を行い差信号を得る。差信号は復元処理手
段５において復元され、エネルギーサブトラクション信
号Ｓ′を得る。エネルギーサブトラクション信号Ｓ′は
画像出力手段６において可視像と して再生される。
Description
【０００１】
【発明の属する技術分野】本発明は、同一被写体の放射
線画像情報を担持する複数の画像信号のサブトラクショ
ン処理を行うエネルギーサブトラクション処理方法およ
び装置に関するものである。
線画像情報を担持する複数の画像信号のサブトラクショ
ン処理を行うエネルギーサブトラクション処理方法およ
び装置に関するものである。
【０００２】
【従来の技術】画像を表す画像信号を得、この画像信号
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば放射線画像の診断
性能を向上させるために、画像信号に対してボケマスク
処理等の周波数強調処理を施す方法が本出願人により提
案されている（特開昭55-163772 等）。この周波数処理
は、原画像を表す画像信号からボケマスク信号を減算し
たものに強調度を乗じたものを加える処理を施すもの
で、これにより画像において所定の空間周波数成分を強
調するようにしたものである。
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば放射線画像の診断
性能を向上させるために、画像信号に対してボケマスク
処理等の周波数強調処理を施す方法が本出願人により提
案されている（特開昭55-163772 等）。この周波数処理
は、原画像を表す画像信号からボケマスク信号を減算し
たものに強調度を乗じたものを加える処理を施すもの
で、これにより画像において所定の空間周波数成分を強
調するようにしたものである。
【０００３】また、画像信号に対して周波数処理を施す
別の方法として、フーリエ変換、ウェーブレット変換、
サブバンド変換等により画像を多重解像度画像に変換す
ることにより画像を表す画像信号を複数の周波数帯域の
信号に分解し、この分解された信号のうち、所望とする
周波数帯域の信号に対して強調等の所定の画像処理を施
す方法が提案されている。
別の方法として、フーリエ変換、ウェーブレット変換、
サブバンド変換等により画像を多重解像度画像に変換す
ることにより画像を表す画像信号を複数の周波数帯域の
信号に分解し、この分解された信号のうち、所望とする
周波数帯域の信号に対して強調等の所定の画像処理を施
す方法が提案されている。
【０００４】また、近年画像処理の分野において、画像
を多重解像度空間に変換する新規な方法としてラプラシ
アンピラミッドなる方法が提案されている（例えば特開
平6-301766号）。このラプラシアンピラミッドは、原画
像に対してガウス関数で近似されたようなマスクにより
マスク処理を施した後、画像をサブサンプリングして画
素数を間引いて半分にすることにより、原画像の１／４
のサイズのボケ画像を得、このボケ画像のサンプリング
された画素に値が０の画素を補間して元の大きさの画像
に戻し、この画像に対してさらに上述したマスクにより
マスク処理を施してボケ画像を得、このボケ画像を原画
像から減算して原画像の所定の周波数帯域を表す細部画
像を得るものである。この処理を得られたボケ画像に対
して繰り返すことにより原画像の１／２2Nの大きさのボ
ケ画像をＮ個作成するものである。ここで、ガウス関数
で近似されたようなマスクによりマスク処理を施した画
像に対してサンプリングを行っているため、実際にはガ
ウシアンフィルタを用いているが、ラプラシアンフィル
タをかけた場合と同様の処理画像が得られる。そしてこ
のように原画像サイズの画像から順に１／２2Nの大きさ
の低周波数帯域の画像が得られるため、この処理の結果
得られた画像はラプラシアンピラミッドと呼ばれる。
を多重解像度空間に変換する新規な方法としてラプラシ
アンピラミッドなる方法が提案されている（例えば特開
平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 andMachi
ne Intelligence、６巻、２号、1984年３月、Mallat S.
G.,“A Theory forMultiresolution Signal Decomposit
ion ；The Wavelet Representation”IEEETrans.on Pat
tern Analysis and Machine Intelligence 、11巻、７
号、1989年７月；Ebrahimi T.,Kunt M.,“Image compre
ssion by Gabor Expansion”，Optical Engineering,30
巻、７号、873 〜880 頁、1991年７月、およびPieter V
uylsteke,Emile Schoeters，“Multiscale Image Contr
ast Amplification ”SPIEVol.2167 Image Processing
(1994),pp551 〜560 に詳細が記載されている。
ては、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 andMachi
ne Intelligence、６巻、２号、1984年３月、Mallat S.
G.,“A Theory forMultiresolution Signal Decomposit
ion ；The Wavelet Representation”IEEETrans.on Pat
tern Analysis and Machine Intelligence 、11巻、７
号、1989年７月；Ebrahimi T.,Kunt M.,“Image compre
ssion by Gabor Expansion”，Optical Engineering,30
巻、７号、873 〜880 頁、1991年７月、およびPieter V
uylsteke,Emile Schoeters，“Multiscale Image Contr
ast Amplification ”SPIEVol.2167 Image Processing
(1994),pp551 〜560 に詳細が記載されている。
【０００６】そしてこのようにして得られたラプラシア
ンピラミッドの全ての周波数帯域の画像に対して強調処
理を施し、この強調処理が施された各周波数帯域の画像
を逆変換して処理済画像を得る方法が上記特開平6-3017
66号に記載されている。このように処理が施された画像
は、各周波数帯域において画像が強調されているため、
実質的に上述したボケマスク処理において複数のサイズ
のマスクによりボケマスク処理を施したような画像とな
っている。
ンピラミッドの全ての周波数帯域の画像に対して強調処
理を施し、この強調処理が施された各周波数帯域の画像
を逆変換して処理済画像を得る方法が上記特開平6-3017
66号に記載されている。このように処理が施された画像
は、各周波数帯域において画像が強調されているため、
実質的に上述したボケマスク処理において複数のサイズ
のマスクによりボケマスク処理を施したような画像とな
っている。
【０００７】一方、従来より放射線画像のサブトラクシ
ョン処理が公知となっている。この放射線画像のサブト
ラクションとは、異なった条件で撮影した２つの放射線
画像を光電的に読み出してデジタル画像信号を得た後、
これらのデジタル画像信号を両画像の各画素を対応させ
て減算処理し、放射線画像中の特定の構造物を抽出させ
る差信号を得る方法であり、このようにして得た差信号
を用いれば、特定構造物のみが抽出された放射線画像を
再生することができる。
ョン処理が公知となっている。この放射線画像のサブト
ラクションとは、異なった条件で撮影した２つの放射線
画像を光電的に読み出してデジタル画像信号を得た後、
これらのデジタル画像信号を両画像の各画素を対応させ
て減算処理し、放射線画像中の特定の構造物を抽出させ
る差信号を得る方法であり、このようにして得た差信号
を用いれば、特定構造物のみが抽出された放射線画像を
再生することができる。
【０００８】このサブトラクション処理には、基本的に
次の２つの方法がある。すなわち、(1) 造影剤注入によ
り特定の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像信号を
引き算（サブトラクト）することによって特定の構造物
を抽出するいわゆる時間サブトラクション処理と、(2)
同一の被写体に対して相異なるエネルギー分布を有する
放射線を照射し、あるいは被写体透過後の放射線をエネ
ルギー分布状態を変えて２つの放射線検出手段に照射し
て、それにより特定の構造物が異なる画像を２つの放射
線画像間に存在せしめ、その後この２つの放射線画像の
画像信号間で適当な重み付けをした上で引き算（サブト
ラクト）を行って、特定の構造物の画像を抽出するいわ
ゆるエネルギーサブトラクション処理である。
次の２つの方法がある。すなわち、(1) 造影剤注入によ
り特定の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像信号を
引き算（サブトラクト）することによって特定の構造物
を抽出するいわゆる時間サブトラクション処理と、(2)
同一の被写体に対して相異なるエネルギー分布を有する
放射線を照射し、あるいは被写体透過後の放射線をエネ
ルギー分布状態を変えて２つの放射線検出手段に照射し
て、それにより特定の構造物が異なる画像を２つの放射
線画像間に存在せしめ、その後この２つの放射線画像の
画像信号間で適当な重み付けをした上で引き算（サブト
ラクト）を行って、特定の構造物の画像を抽出するいわ
ゆるエネルギーサブトラクション処理である。
【０００９】先に述べた蓄積性蛍光体シートを利用する
放射線画像情報記録再生システムにおいては、蓄積性蛍
光体シートに記録されている放射線画像情報が直接電気
的画像信号の形で読み取られるから、このシステムによ
れば、上述のようなサブトラクション処理を容易に行う
ことが可能となる。この蓄積性蛍光体シートを用いてエ
ネルギーサブトラクション処理を行うためには、例えば
２枚の蓄積性蛍光体シートに特定の構造物に対応する部
分の画像情報が異なるように画像記録（撮影）を行えば
よく、具体的には、エネルギー分布の異なる２種類の放
射線を用いて撮影を２回行う２ショット法と、例えば被
写体を透過した放射線を、重ねられた２枚の蓄積性蛍光
体シート（それらは互いに接していても、離れていても
よい）に同時に曝射することによって、両シートに互い
にエネルギー分布が異なる放射線を照射するようにした
１ショット法が知られている。
放射線画像情報記録再生システムにおいては、蓄積性蛍
光体シートに記録されている放射線画像情報が直接電気
的画像信号の形で読み取られるから、このシステムによ
れば、上述のようなサブトラクション処理を容易に行う
ことが可能となる。この蓄積性蛍光体シートを用いてエ
ネルギーサブトラクション処理を行うためには、例えば
２枚の蓄積性蛍光体シートに特定の構造物に対応する部
分の画像情報が異なるように画像記録（撮影）を行えば
よく、具体的には、エネルギー分布の異なる２種類の放
射線を用いて撮影を２回行う２ショット法と、例えば被
写体を透過した放射線を、重ねられた２枚の蓄積性蛍光
体シート（それらは互いに接していても、離れていても
よい）に同時に曝射することによって、両シートに互い
にエネルギー分布が異なる放射線を照射するようにした
１ショット法が知られている。
【００１０】
【発明が解決しようとする課題】このようなエネルギー
サブトラクション処理を施すことにより、画像中の特定
の構造物に関する情報のみを抽出することができるが、
差信号を得る際に画像の空間周波数に関わらず重み付け
がなされるため、サブトラクション処理により画像中の
ノイズなどの高周波数成分が強調されてしまい、得られ
るエネルギーサブトラクション画像の画質が低下し、観
察がしにくいものとなっていた。
サブトラクション処理を施すことにより、画像中の特定
の構造物に関する情報のみを抽出することができるが、
差信号を得る際に画像の空間周波数に関わらず重み付け
がなされるため、サブトラクション処理により画像中の
ノイズなどの高周波数成分が強調されてしまい、得られ
るエネルギーサブトラクション画像の画質が低下し、観
察がしにくいものとなっていた。
【００１１】本発明による発明は上記事情に鑑み、得ら
れるサブトラクション画像の画質を向上させることがで
きるエネルギーサブトラクション処理方法および装置を
提供することを目的とするものである。
れるサブトラクション画像の画質を向上させることがで
きるエネルギーサブトラクション処理方法および装置を
提供することを目的とするものである。
【００１２】
【課題を解決するための手段】本発明によるエネルギー
サブトラクション処理方法および装置は、同一被写体を
透過したそれぞれエネルギー分布が異なる放射線により
得られた少なくとも一部の画像情報が互いに異なる複数
の放射線画像を表す複数の画像信号のそれぞれを、多重
解像度空間に変換することにより得られる複数の周波数
帯域ごとの周波数信号に対して、相対応する画素につい
ての信号間で前記各周波数信号に周波数帯域に応じた重
み付け係数を乗じて減算を行って前記複数の周波数帯域
ごとの差信号を得、この差信号を逆変換することにより
エネルギーサブトラクション信号を得ることを特徴とす
るものである。
サブトラクション処理方法および装置は、同一被写体を
透過したそれぞれエネルギー分布が異なる放射線により
得られた少なくとも一部の画像情報が互いに異なる複数
の放射線画像を表す複数の画像信号のそれぞれを、多重
解像度空間に変換することにより得られる複数の周波数
帯域ごとの周波数信号に対して、相対応する画素につい
ての信号間で前記各周波数信号に周波数帯域に応じた重
み付け係数を乗じて減算を行って前記複数の周波数帯域
ごとの差信号を得、この差信号を逆変換することにより
エネルギーサブトラクション信号を得ることを特徴とす
るものである。
【００１３】
【発明の効果】本発明によるエネルギーサブトラクショ
ン処理方法および装置は、エネルギーサブトラクション
処理を施す複数の画像信号のそれぞれを多重解像度空間
に変換して複数の周波数帯域ごとの周波数信号を得、こ
の複数の周波数信号ごとの差信号を得る際の重み付け係
数を各周波数信号の周波数帯域に応じて、変更するよう
にしたものである。例えば、各周波数帯域のうち比較的
高周波数帯域の信号についてはノイズ成分を多く含んで
いるため、ノイズを低減させるような重み付け係数によ
り差信号を得るようにしたものである。このようにノイ
ズなどの不要な周波数成分を含む周波数帯域については
この不要周波数成分が目立たないような重み付け係数を
設定することができるため、各周波数帯域の差信号を逆
変換することにより得られるエネルギーサブトラクショ
ン信号を可視像として再生することにより、ノイズなど
の不要な周波数成分が目立たない、観察に適した良好な
エネルギーサブトラクション画像を得ることができる。
ン処理方法および装置は、エネルギーサブトラクション
処理を施す複数の画像信号のそれぞれを多重解像度空間
に変換して複数の周波数帯域ごとの周波数信号を得、こ
の複数の周波数信号ごとの差信号を得る際の重み付け係
数を各周波数信号の周波数帯域に応じて、変更するよう
にしたものである。例えば、各周波数帯域のうち比較的
高周波数帯域の信号についてはノイズ成分を多く含んで
いるため、ノイズを低減させるような重み付け係数によ
り差信号を得るようにしたものである。このようにノイ
ズなどの不要な周波数成分を含む周波数帯域については
この不要周波数成分が目立たないような重み付け係数を
設定することができるため、各周波数帯域の差信号を逆
変換することにより得られるエネルギーサブトラクショ
ン信号を可視像として再生することにより、ノイズなど
の不要な周波数成分が目立たない、観察に適した良好な
エネルギーサブトラクション画像を得ることができる。
【００１４】また、複数の周波数帯域ごとの周波数信号
は、強調処理などの種々の画像処理を施すために得ら
れ、記憶手段に記憶されて画像処理の際にこの記憶手段
から読み出されて画像処理に供されるものである。本発
明によるエネルギーサブトラクション処理方法および装
置はこの記憶手段に予め記憶された周波数信号を用いる
ことができるため、新たな画像を得る必要が無くなり、
エネルギーサブトラクション処理を行うための時間を短
縮し、高速に処理を行うことができる。
は、強調処理などの種々の画像処理を施すために得ら
れ、記憶手段に記憶されて画像処理の際にこの記憶手段
から読み出されて画像処理に供されるものである。本発
明によるエネルギーサブトラクション処理方法および装
置はこの記憶手段に予め記憶された周波数信号を用いる
ことができるため、新たな画像を得る必要が無くなり、
エネルギーサブトラクション処理を行うための時間を短
縮し、高速に処理を行うことができる。
【００１５】
【発明の実施の形態】以下図面を参照して本発明の実施
の形態について説明する。
の形態について説明する。
【００１６】図１は本発明によるエネルギーサブトラク
ション処理方法を実施するための装置の概略を表すブロ
ック図である。図１に示すように本発明によるエネルギ
ーサブトラクション処理方法を実施するための装置は、
装置にエネルギーサブトラクション処理を施す画像信号
を入力するための画像入力手段１と、入力された画像信
号に対して多重解像度分解処理を施して複数の周波数帯
域ごとの周波数信号を得る多重解像度分解処理手段２
と、多重解像度分解処理手段２において複数の周波数帯
域に分解された周波数信号を記憶する記憶手段３と、記
憶手段３に記憶された周波数信号を読み出してエネルギ
ーサブトラクション処理を施すエネルギーサブトラクシ
ョン処理手段４と、エネルギーサブトラクション処理手
段４によりエネルギーサブトラクション処理が施された
差信号Ｓsub を復元してエネルギーサブトラクション信
号を得るための復元処理手段５と、復元処理手段５によ
り復元されたエネルギーサブトラクション信号を可視像
として再生するための画像出力手段６と、エネルギーサ
ブトラクション処理手段４の他に周波数信号に対して強
調処理などの画像処理を施すための画像処理手段７とか
らなるものである。
ション処理方法を実施するための装置の概略を表すブロ
ック図である。図１に示すように本発明によるエネルギ
ーサブトラクション処理方法を実施するための装置は、
装置にエネルギーサブトラクション処理を施す画像信号
を入力するための画像入力手段１と、入力された画像信
号に対して多重解像度分解処理を施して複数の周波数帯
域ごとの周波数信号を得る多重解像度分解処理手段２
と、多重解像度分解処理手段２において複数の周波数帯
域に分解された周波数信号を記憶する記憶手段３と、記
憶手段３に記憶された周波数信号を読み出してエネルギ
ーサブトラクション処理を施すエネルギーサブトラクシ
ョン処理手段４と、エネルギーサブトラクション処理手
段４によりエネルギーサブトラクション処理が施された
差信号Ｓsub を復元してエネルギーサブトラクション信
号を得るための復元処理手段５と、復元処理手段５によ
り復元されたエネルギーサブトラクション信号を可視像
として再生するための画像出力手段６と、エネルギーサ
ブトラクション処理手段４の他に周波数信号に対して強
調処理などの画像処理を施すための画像処理手段７とか
らなるものである。
【００１７】図２は、本発明に用いられるエネルギーサ
ブトラクション処理を施すべき画像信号を得るための放
射線画像撮影装置を表す図である。
ブトラクション処理を施すべき画像信号を得るための放
射線画像撮影装置を表す図である。
【００１８】フィルタ46を間に挟み蓄積性蛍光体シート
45および47がシート47を下にして重ねられている。この
上には、被写体44を介して放射線43を発する放射線源42
が配置されている。以上のように、放射線撮影装置41が
構成されている。
45および47がシート47を下にして重ねられている。この
上には、被写体44を介して放射線43を発する放射線源42
が配置されている。以上のように、放射線撮影装置41が
構成されている。
【００１９】この放射線源42から発せられた放射線43が
被写体44に照射される。被写体44を透過した放射線43a
は第１の蓄積性蛍光体シート45に照射され、放射線43a
のエネルギーの一部が第１の蓄積性蛍光体シート45に記
録され、これによりシート45に被写体44の放射線画像
（高圧画像）が蓄積記録される。シート45を透過した放
射線43b はさらにフィルタ46を透過し、フィルタ46を透
過した放射線43c が第２の蓄積性蛍光体シート47に照射
される。これによりシート47にも被写体44の放射線画像
（低圧画像）が蓄積記録される。
被写体44に照射される。被写体44を透過した放射線43a
は第１の蓄積性蛍光体シート45に照射され、放射線43a
のエネルギーの一部が第１の蓄積性蛍光体シート45に記
録され、これによりシート45に被写体44の放射線画像
（高圧画像）が蓄積記録される。シート45を透過した放
射線43b はさらにフィルタ46を透過し、フィルタ46を透
過した放射線43c が第２の蓄積性蛍光体シート47に照射
される。これによりシート47にも被写体44の放射線画像
（低圧画像）が蓄積記録される。
【００２０】図３は、本発明に用いる放射線画像を読み
取る放射線画像読取装置の斜視図である。
取る放射線画像読取装置の斜視図である。
【００２１】図２に示す放射線撮影装置41で撮影が行わ
れた後、第１および第２の蓄積性蛍光体シート45，47が
１枚ずつ放射線画像読取装置の所定位置にセットされ
る。ここでは、第１の蓄積性蛍光体シート45に蓄積記録
された第１の放射線画像の読取りの場合について説明す
る。
れた後、第１および第２の蓄積性蛍光体シート45，47が
１枚ずつ放射線画像読取装置の所定位置にセットされ
る。ここでは、第１の蓄積性蛍光体シート45に蓄積記録
された第１の放射線画像の読取りの場合について説明す
る。
【００２２】所定位置にセットされた、第１の放射線画
像が蓄積記録された蓄積性蛍光体シート45は、図示しな
い駆動手段により駆動されるエンドレスベルト等のシー
ト搬送手段27により、矢印Ｙ方向に搬送（副走査）され
る。一方、レーザ光源16から発せられた光ビーム17はモ
ータ18により駆動され矢印Ｚ方向に高速回転する回転多
面鏡19によって反射偏向され、ｆθレンズ等の集束レン
ズ20を透過した後、ミラー21により光路を変えて蓄積性
蛍光体シート45に入射し、副走査の方向（矢印Ｙ方向）
と略垂直な矢印Ｘ方向に主走査する。蓄積性蛍光体シー
ト45の光ビーム17が照射されたか所からは、蓄積記録さ
れている放射線画像情報に応じた光量の輝尽発光光22が
発せられ、この輝尽発光光22は光ガイド23によって導か
れ、フォトマルチプライヤ（光電子増倍管）24によって
光電的に検出される。光ガイド23はアクリル板等の導光
性材料を成形して作られたものであり、直線状をなす入
射端面23a が蓄積性蛍光体シート45上の主走査線にそっ
て延びるように配され、円環状に形成された射出端面23
b にフォトマルチプライヤ24の受光面が結合されてい
る。入射端面23a から光ガイド23内に入射した輝尽発光
光22は、該光ガイド23の内部を全反射を繰り返して進
み、射出端面23b から射出してフォトマルチプライヤ24
に受光され、放射線画像を表す輝尽発光光22がフォトマ
ルチプライヤ24によって電気信号に変換される。フォト
マルチプライヤ24から出力されたアナログ信号Ｓは、ロ
グアンプ25で対数的に増幅された後、Ａ／Ｄ変換器26に
入力され、サンプリングされてデジタルの画像信号ＳＯ
が得られる。この画像信号ＳＯは第１の蓄積性蛍光体シ
ート45に蓄積記録された第１の放射線画像を表すもので
あり、第１の画像信号ＳＯ1 と呼ぶ。この第１の画像信
号ＳＯ1 は本発明によるエネルギーサブトラクション処
理方法を実施するための装置を内包する画像処理表示装
置30内の記憶手段に一旦記録される。
像が蓄積記録された蓄積性蛍光体シート45は、図示しな
い駆動手段により駆動されるエンドレスベルト等のシー
ト搬送手段27により、矢印Ｙ方向に搬送（副走査）され
る。一方、レーザ光源16から発せられた光ビーム17はモ
ータ18により駆動され矢印Ｚ方向に高速回転する回転多
面鏡19によって反射偏向され、ｆθレンズ等の集束レン
ズ20を透過した後、ミラー21により光路を変えて蓄積性
蛍光体シート45に入射し、副走査の方向（矢印Ｙ方向）
と略垂直な矢印Ｘ方向に主走査する。蓄積性蛍光体シー
ト45の光ビーム17が照射されたか所からは、蓄積記録さ
れている放射線画像情報に応じた光量の輝尽発光光22が
発せられ、この輝尽発光光22は光ガイド23によって導か
れ、フォトマルチプライヤ（光電子増倍管）24によって
光電的に検出される。光ガイド23はアクリル板等の導光
性材料を成形して作られたものであり、直線状をなす入
射端面23a が蓄積性蛍光体シート45上の主走査線にそっ
て延びるように配され、円環状に形成された射出端面23
b にフォトマルチプライヤ24の受光面が結合されてい
る。入射端面23a から光ガイド23内に入射した輝尽発光
光22は、該光ガイド23の内部を全反射を繰り返して進
み、射出端面23b から射出してフォトマルチプライヤ24
に受光され、放射線画像を表す輝尽発光光22がフォトマ
ルチプライヤ24によって電気信号に変換される。フォト
マルチプライヤ24から出力されたアナログ信号Ｓは、ロ
グアンプ25で対数的に増幅された後、Ａ／Ｄ変換器26に
入力され、サンプリングされてデジタルの画像信号ＳＯ
が得られる。この画像信号ＳＯは第１の蓄積性蛍光体シ
ート45に蓄積記録された第１の放射線画像を表すもので
あり、第１の画像信号ＳＯ1 と呼ぶ。この第１の画像信
号ＳＯ1 は本発明によるエネルギーサブトラクション処
理方法を実施するための装置を内包する画像処理表示装
置30内の記憶手段に一旦記録される。
【００２３】この画像処理表示装置30は、本発明による
エネルギーサブトラクション処理装置を内包するもので
あり、種々の指示を入力するキーボード31、指示のため
の補助情報や画像信号に基づく可視画像を表示するＣＲ
Ｔ３、補助記憶媒体としてのフロッピーディスクが装填
され駆動されるフロッピィディスク駆動装置33、および
ＣＰＵや内部メモリが内蔵された本体部34が備えられて
いる。
エネルギーサブトラクション処理装置を内包するもので
あり、種々の指示を入力するキーボード31、指示のため
の補助情報や画像信号に基づく可視画像を表示するＣＲ
Ｔ３、補助記憶媒体としてのフロッピーディスクが装填
され駆動されるフロッピィディスク駆動装置33、および
ＣＰＵや内部メモリが内蔵された本体部34が備えられて
いる。
【００２４】次に上記と同様にして、第２の蓄積性蛍光
体シート47に蓄積記録された第２の放射線画像を表す第
２の画像信号ＳＯ2 が得られ、この第２の画像信号ＳＯ
2 も画像処理表示装置30内の記憶手段に一旦記憶され
る。
体シート47に蓄積記録された第２の放射線画像を表す第
２の画像信号ＳＯ2 が得られ、この第２の画像信号ＳＯ
2 も画像処理表示装置30内の記憶手段に一旦記憶され
る。
【００２５】このようにして得られた第１および第２の
画像信号ＳＯ1 ，ＳＯ2 は図１に示す画像入力手段１に
より多重解像度分解処理手段２に入力される。多重解像
度分解処理手段２においては、第１および第２の画像信
号ＳＯ1 ，ＳＯ2 に対して以下のような処理が施され
る。図４は図１における多重解像度分解処理手段２にお
いて行われる処理を説明するためのブロック図である。
なお、本実施の形態においてはラプラシアンピラミッド
の手法により第１の画像信号ＳＯ1 を多重解像度画像に
分解するものとする。図４に示すように原画像を表すデ
ジタルの画像信号Ｓが多重解像度分解処理手段２に入力
されると、フィルタリング手段10においてローパスフィ
ルタによりフィルタリングされる。このローパスフィル
タは例えば図５に示すように５×５のグリッド上の二次
元ガウス分布に略対応している。このローパスフィルタ
は後述するように全ての解像度の画像に対して適用され
る。
画像信号ＳＯ1 ，ＳＯ2 は図１に示す画像入力手段１に
より多重解像度分解処理手段２に入力される。多重解像
度分解処理手段２においては、第１および第２の画像信
号ＳＯ1 ，ＳＯ2 に対して以下のような処理が施され
る。図４は図１における多重解像度分解処理手段２にお
いて行われる処理を説明するためのブロック図である。
なお、本実施の形態においてはラプラシアンピラミッド
の手法により第１の画像信号ＳＯ1 を多重解像度画像に
分解するものとする。図４に示すように原画像を表すデ
ジタルの画像信号Ｓが多重解像度分解処理手段２に入力
されると、フィルタリング手段10においてローパスフィ
ルタによりフィルタリングされる。このローパスフィル
タは例えば図５に示すように５×５のグリッド上の二次
元ガウス分布に略対応している。このローパスフィルタ
は後述するように全ての解像度の画像に対して適用され
る。
【００２６】このようなローパスフィルタによりフィル
タリングされた画像信号Ｓはフィルタリング手段10にお
いて１画素おきにサンプリングされ、低解像度近似画像
信号ｇｈ1 が得られる。この低解像度近似画像信号ｇｈ
1 により表される画像は、原画像の１／４の大きさとな
っている。次いで補間手段11において、この低解像度近
似画像信号ｇｈ1 のサンプリングされた間隔に値が０の
画素が補間される。この補間は低解像度近似画像信号ｇ
ｈ1 の一列ごとおよび一行ごとに値が０の行および列を
挿入することにより行う。このように値が０の画素が補
間された低解像度近似画像信号ｇｈ1 はぼけてはいるも
のの一画素おきに値が０の画素が挿入されているため、
信号値の変化が滑らかではないものとなっている。
タリングされた画像信号Ｓはフィルタリング手段10にお
いて１画素おきにサンプリングされ、低解像度近似画像
信号ｇｈ1 が得られる。この低解像度近似画像信号ｇｈ
1 により表される画像は、原画像の１／４の大きさとな
っている。次いで補間手段11において、この低解像度近
似画像信号ｇｈ1 のサンプリングされた間隔に値が０の
画素が補間される。この補間は低解像度近似画像信号ｇ
ｈ1 の一列ごとおよび一行ごとに値が０の行および列を
挿入することにより行う。このように値が０の画素が補
間された低解像度近似画像信号ｇｈ1 はぼけてはいるも
のの一画素おきに値が０の画素が挿入されているため、
信号値の変化が滑らかではないものとなっている。
【００２７】そしてこのようにして補間が行われた後、
さらにこの補間がなされた低解像度近似画像信号ｇｈ1
に対して図５に示すローパスフィルタにより再度フィル
タリング処理を施し、低解像度近似画像信号ｇｈ1 ′を
得る。この低解像度近似画像信号ｇｈ1 ′により表され
る画像は上述した補間がなされた低解像度近似画像信号
ｇｈ1 により表される画像と比較して信号値の変化が滑
らかなものとなっている。また原画像と比較して周波数
帯域的には半分より高い高周波数が消えたような画像と
なっている。これは画像の大きさを１／４にして一画素
おきに値が０の画素を補間し、さらに図５に示すローパ
スフィルタによりフィルタリング処理を施しているた
め、ガウス関数により空間周波数が半分よりも高い周波
数帯域の画像がぼかされたようになっているからであ
る。
さらにこの補間がなされた低解像度近似画像信号ｇｈ1
に対して図５に示すローパスフィルタにより再度フィル
タリング処理を施し、低解像度近似画像信号ｇｈ1 ′を
得る。この低解像度近似画像信号ｇｈ1 ′により表され
る画像は上述した補間がなされた低解像度近似画像信号
ｇｈ1 により表される画像と比較して信号値の変化が滑
らかなものとなっている。また原画像と比較して周波数
帯域的には半分より高い高周波数が消えたような画像と
なっている。これは画像の大きさを１／４にして一画素
おきに値が０の画素を補間し、さらに図５に示すローパ
スフィルタによりフィルタリング処理を施しているた
め、ガウス関数により空間周波数が半分よりも高い周波
数帯域の画像がぼかされたようになっているからであ
る。
【００２８】次いで減算器12において、原画像から低解
像度近似画像信号ｇｈ1 ′の減算が行われ、細部画像信
号ｂｈ0 が得られる。この減算は原画像と低解像度近似
画像信号ｇｈ 1′との相対応する画素についての信号間
で行われる。ここで、低解像度近似画像信号ｇｈ1 ′は
上述したように原画像の空間周波数のうち半分より高い
周波数帯域の画像がぼけたようになっているため、細部
画像信号ｂｈ0 は原画像のうち半分より上の周波数帯域
のみを表す画像となっている。すなわち、図６に示すよ
うに細部画像信号ｂｈ0 は原画像のナイキスト周波数Ｎ
のうちＮ／２〜Ｎの周波数帯域の画像を表すものとなっ
ている。
像度近似画像信号ｇｈ1 ′の減算が行われ、細部画像信
号ｂｈ0 が得られる。この減算は原画像と低解像度近似
画像信号ｇｈ 1′との相対応する画素についての信号間
で行われる。ここで、低解像度近似画像信号ｇｈ1 ′は
上述したように原画像の空間周波数のうち半分より高い
周波数帯域の画像がぼけたようになっているため、細部
画像信号ｂｈ0 は原画像のうち半分より上の周波数帯域
のみを表す画像となっている。すなわち、図６に示すよ
うに細部画像信号ｂｈ0 は原画像のナイキスト周波数Ｎ
のうちＮ／２〜Ｎの周波数帯域の画像を表すものとなっ
ている。
【００２９】次いで、低解像度近似画像信号ｇｈ1 はフ
ィルタリング手段10に入力され、図５に示すローパスフ
ィルタによりフィルタリング処理が施される。そしてフ
ィルタリング処理が施された低解像度近似画像信号ｇｈ
1 は、フィルタリング手段10において１画素おきにサン
プリングされ、低解像度近似画像信号ｇｈ2 が得られ
る。この低解像度近似画像信号ｇｈ2 により表される画
像は、低解像度近似画像信号ｇｈ1 により表される画像
の１／４すなわち原画像の１／16の大きさとなってい
る。次いで補間手段11において、この低解像度近似画像
信号ｇｈ2 のサンプリングされた間隔に値が０の画素が
補間される。この補間は低解像度近似画像信号ｇｈ2 の
一列ごとおよび一行ごとに値が０の行および列を挿入す
ることにより行う。このように値が０の画素が補間され
た低解像度近似画像信号ｇｈ2 はぼけてはいるものの一
画素おきに値が０の画素が挿入されているため、信号値
の変化が滑らかではないものとなっている。
ィルタリング手段10に入力され、図５に示すローパスフ
ィルタによりフィルタリング処理が施される。そしてフ
ィルタリング処理が施された低解像度近似画像信号ｇｈ
1 は、フィルタリング手段10において１画素おきにサン
プリングされ、低解像度近似画像信号ｇｈ2 が得られ
る。この低解像度近似画像信号ｇｈ2 により表される画
像は、低解像度近似画像信号ｇｈ1 により表される画像
の１／４すなわち原画像の１／16の大きさとなってい
る。次いで補間手段11において、この低解像度近似画像
信号ｇｈ2 のサンプリングされた間隔に値が０の画素が
補間される。この補間は低解像度近似画像信号ｇｈ2 の
一列ごとおよび一行ごとに値が０の行および列を挿入す
ることにより行う。このように値が０の画素が補間され
た低解像度近似画像信号ｇｈ2 はぼけてはいるものの一
画素おきに値が０の画素が挿入されているため、信号値
の変化が滑らかではないものとなっている。
【００３０】そしてこのようにして補間が行われた後、
さらにこの補間がなされた低解像度近似画像信号ｇｈ2
に対して図５に示すローパスフィルタにより再度フィル
タリング処理を施し、低解像度近似画像信号ｇｈ2 ′を
得る。この低解像度近似画像信号ｇｈ2 ′は上述した補
間がなされた低解像度近似画像信号ｇｈ2 と比較して信
号値の変化が滑らかなものとなっている。また低解像度
近似画像信号ｇｈ1 と比較して周波数帯域的には半分よ
り高い周波数帯域の画像が消えたようになっている。
さらにこの補間がなされた低解像度近似画像信号ｇｈ2
に対して図５に示すローパスフィルタにより再度フィル
タリング処理を施し、低解像度近似画像信号ｇｈ2 ′を
得る。この低解像度近似画像信号ｇｈ2 ′は上述した補
間がなされた低解像度近似画像信号ｇｈ2 と比較して信
号値の変化が滑らかなものとなっている。また低解像度
近似画像信号ｇｈ1 と比較して周波数帯域的には半分よ
り高い周波数帯域の画像が消えたようになっている。
【００３１】次いで減算器12において、低解像度近似画
像信号ｇｈ1 から低解像度近似画像信号ｇｈ2 ′の減算
が行われ、細部画像信号ｂｈ1 が得られる。この減算は
低解像度近似画像信号ｇｈ1 と低解像度近似画像信号ｇ
ｈ2 ′との相対応する画素についての信号間で行われ
る。ここで、低解像度近似画像信号ｇｈ2 ′は上述した
ように低解像度近似画像信号ｇｈ1 の空間周波数のうち
半分より高い周波数帯域の画像がぼけたようになってい
るため、細部画像信号ｂｈ1 は低解像度近似画像信号ｇ
ｈ1 のうち半分より上の周波数帯域のみを表す画像とな
っている。すなわち、図４に示すように細部画像信号ｂ
ｈ1 は低解像度近似画像信号ｇｈ1 のうちの半分より上
の周波数帯域のみ、すなわち原画像のナイキスト周波数
ＮのうちＮ／４〜Ｎ／２の周波数帯域の画像を表すもの
となっている。このようにガウス分布のローパスフィル
タによりフィルタリング処理を施して細部画像信号を得
るようにしているが、フィルタリング処理が施された画
像を低解像度近似画像信号から減算していることから、
実質的にはラプラシアンフィルタによりフィルタリング
処理を施した場合と同様の結果となる。
像信号ｇｈ1 から低解像度近似画像信号ｇｈ2 ′の減算
が行われ、細部画像信号ｂｈ1 が得られる。この減算は
低解像度近似画像信号ｇｈ1 と低解像度近似画像信号ｇ
ｈ2 ′との相対応する画素についての信号間で行われ
る。ここで、低解像度近似画像信号ｇｈ2 ′は上述した
ように低解像度近似画像信号ｇｈ1 の空間周波数のうち
半分より高い周波数帯域の画像がぼけたようになってい
るため、細部画像信号ｂｈ1 は低解像度近似画像信号ｇ
ｈ1 のうち半分より上の周波数帯域のみを表す画像とな
っている。すなわち、図４に示すように細部画像信号ｂ
ｈ1 は低解像度近似画像信号ｇｈ1 のうちの半分より上
の周波数帯域のみ、すなわち原画像のナイキスト周波数
ＮのうちＮ／４〜Ｎ／２の周波数帯域の画像を表すもの
となっている。このようにガウス分布のローパスフィル
タによりフィルタリング処理を施して細部画像信号を得
るようにしているが、フィルタリング処理が施された画
像を低解像度近似画像信号から減算していることから、
実質的にはラプラシアンフィルタによりフィルタリング
処理を施した場合と同様の結果となる。
【００３２】そして上述した処理をフィルタリング手段
10によりフィルタリングされかつサンプリングされた低
解像度近似画像信号ｇｈk （ｋ＝１〜Ｎ）に対して順次
繰り返し行い、図４に示すようにｎ個の細部画像信号ｂ
ｈk （ｋ＝１〜ｎ）および低解像度近似画像信号の残留
画像信号ｇｈL を得る。ここで、細部画像信号ｂｈ
kは、ｂｈ0 から順に解像度が低くなる、すなわち画像
の周波数帯域が低くなるものであり、原画像のナイキス
ト周波数Ｎに対して、細部画像信号ｂｈk はＮ／２k+1
〜Ｎ／２k の周波数帯域を表し、画像の大きさが原画像
の１／２2k倍となっている。すなわち、最も解像度が高
い細部画像信号ｂｈ0 により表される画像は原画像と同
じ大きさであるが、細部画像信号ｂｈ0 の次に高解像度
の細部画像信号ｂｈ1 により表される画像は原画像の大
きさの１／４となっている。このように、細部画像信号
により表される画像が原画像と同一の大きさのものから
順次小さくなり、また細部画像信号はラプラシアンフィ
ルタを施したものと実質的に同一の画像であることか
ら、本実施の形態による多重解像度変換はラプラシアン
ピラミッドと呼ばれるものである。また、残留画像信号
ｇｈL は原画像の非常に解像度が低い近似画像信号であ
ると見なすことができ、極端な場合は、残留画像信号ｇ
ｈL は原画像の平均値を表す１つだけの画像信号からな
るものとなる。そしてこのようにして得られた細部画像
信号ｂｈk および残留画像信号ｇｈL は図１に示す記憶
手段３に記憶される。
10によりフィルタリングされかつサンプリングされた低
解像度近似画像信号ｇｈk （ｋ＝１〜Ｎ）に対して順次
繰り返し行い、図４に示すようにｎ個の細部画像信号ｂ
ｈk （ｋ＝１〜ｎ）および低解像度近似画像信号の残留
画像信号ｇｈL を得る。ここで、細部画像信号ｂｈ
kは、ｂｈ0 から順に解像度が低くなる、すなわち画像
の周波数帯域が低くなるものであり、原画像のナイキス
ト周波数Ｎに対して、細部画像信号ｂｈk はＮ／２k+1
〜Ｎ／２k の周波数帯域を表し、画像の大きさが原画像
の１／２2k倍となっている。すなわち、最も解像度が高
い細部画像信号ｂｈ0 により表される画像は原画像と同
じ大きさであるが、細部画像信号ｂｈ0 の次に高解像度
の細部画像信号ｂｈ1 により表される画像は原画像の大
きさの１／４となっている。このように、細部画像信号
により表される画像が原画像と同一の大きさのものから
順次小さくなり、また細部画像信号はラプラシアンフィ
ルタを施したものと実質的に同一の画像であることか
ら、本実施の形態による多重解像度変換はラプラシアン
ピラミッドと呼ばれるものである。また、残留画像信号
ｇｈL は原画像の非常に解像度が低い近似画像信号であ
ると見なすことができ、極端な場合は、残留画像信号ｇ
ｈL は原画像の平均値を表す１つだけの画像信号からな
るものとなる。そしてこのようにして得られた細部画像
信号ｂｈk および残留画像信号ｇｈL は図１に示す記憶
手段３に記憶される。
【００３３】一方、第２の画像信号ＳＯ2 に対してもラ
プラシアンピラミッドの手法により多重解像度変換処理
が施され、図６に示すように第１の画像信号ＳＯ1 と同
様の周波数帯域の細部画像信号ｂｌk および残留画像信
号ｇｌL が得られ図１に示す記憶手段３に記憶される。
プラシアンピラミッドの手法により多重解像度変換処理
が施され、図６に示すように第１の画像信号ＳＯ1 と同
様の周波数帯域の細部画像信号ｂｌk および残留画像信
号ｇｌL が得られ図１に示す記憶手段３に記憶される。
【００３４】ついで、記憶手段３に記憶された細部画像
信号ｂｈk ，ｂｌk および残留画像信号ｇｈL ，ｇｌL
が読み出されてエネルギーサブトラクション処理手段４
に入力される。エネルギーサブトラクション処理手段４
においては細部画像信号ｂｈk ，ｂｌk および残留画像
信号ｇｈL ，ｇｌL に対して以下に示すようなサブトラ
クション処理が施される。すなわち、各周波数帯域の細
部画像信号ｂｈk ，ｂｌk および残留画像信号ｇｈL ，
ｇｌL に対して下記の式 Ｓsub k ＝Ｋｈk ×ｂｈk −Ｋｌk ×ｂｌk Ｓsubg＝ＫｈL ×ｇｈL −ＫｌL ×ｇｌL …(1) 但し、Ｋｈk ，Ｋｌk ：細部画像信号ｂｈk ，ｂｌk
の重み付け係数 ＫｈL ，ＫｌL ：残留画像信号ｇｈL ，ｇｌk の重み付
け係数 により重み付け減算が行われ、これにより各周波数帯域
ごとの２つの放射線画像の差の画像に対応する差信号Ｓ
sub ｂk ，Ｓsubgが得られる。この際、重み付け係数と
しては、以下の表１に示すように、高周波数帯域におい
ては、第１の細部画像信号ｂｈk の重み付け係数Ｋｈk
を０、第２の細部画像信号ｂｌk の重み付け係数を−
１として第１の細部画像信号ｂｈk については後述する
逆変換に用いないようにするものである。これは、エネ
ルギーサブトラクション処理は２つの画像信号の差分値
を求めるものであるため、画像中のノイズのようにラン
ダムに現れる成分が強調されてしまうものである。そし
てこのノイズ成分が強調されることを防止するために、
第１の細部画像信号ｂｈk の重み付け係数を０としてい
るものである。なお、他の周波数帯域の細部画像信号ｂ
ｈk ，ｂｌk および残留画像信号ｇｈL ，ｇｌk につい
ては、３：２の重み付け係数によりサブトラクション処
理が行われる。
信号ｂｈk ，ｂｌk および残留画像信号ｇｈL ，ｇｌL
が読み出されてエネルギーサブトラクション処理手段４
に入力される。エネルギーサブトラクション処理手段４
においては細部画像信号ｂｈk ，ｂｌk および残留画像
信号ｇｈL ，ｇｌL に対して以下に示すようなサブトラ
クション処理が施される。すなわち、各周波数帯域の細
部画像信号ｂｈk ，ｂｌk および残留画像信号ｇｈL ，
ｇｌL に対して下記の式 Ｓsub k ＝Ｋｈk ×ｂｈk −Ｋｌk ×ｂｌk Ｓsubg＝ＫｈL ×ｇｈL −ＫｌL ×ｇｌL …(1) 但し、Ｋｈk ，Ｋｌk ：細部画像信号ｂｈk ，ｂｌk
の重み付け係数 ＫｈL ，ＫｌL ：残留画像信号ｇｈL ，ｇｌk の重み付
け係数 により重み付け減算が行われ、これにより各周波数帯域
ごとの２つの放射線画像の差の画像に対応する差信号Ｓ
sub ｂk ，Ｓsubgが得られる。この際、重み付け係数と
しては、以下の表１に示すように、高周波数帯域におい
ては、第１の細部画像信号ｂｈk の重み付け係数Ｋｈk
を０、第２の細部画像信号ｂｌk の重み付け係数を−
１として第１の細部画像信号ｂｈk については後述する
逆変換に用いないようにするものである。これは、エネ
ルギーサブトラクション処理は２つの画像信号の差分値
を求めるものであるため、画像中のノイズのようにラン
ダムに現れる成分が強調されてしまうものである。そし
てこのノイズ成分が強調されることを防止するために、
第１の細部画像信号ｂｈk の重み付け係数を０としてい
るものである。なお、他の周波数帯域の細部画像信号ｂ
ｈk ，ｂｌk および残留画像信号ｇｈL ，ｇｌk につい
ては、３：２の重み付け係数によりサブトラクション処
理が行われる。
【００３５】 表 １ Ｋｈ0 ：Ｋｌ0 ＝０：−１ Ｋｈ1 ：Ｋｌ1 ＝０：−１ Ｋｈ2 ：Ｋｌ2 ＝３：２ ： ： ： ： Ｋｈn ：Ｋｌn ＝３：２ ＫｈL ：ＫｌL ＝３：２ 次いでこのようにして得られた差信号Ｓsub ｂk ，Ｓsu
bgを逆変換する。この逆変換の処理は復元処理手段５に
おいて以下のようにして行われる。
bgを逆変換する。この逆変換の処理は復元処理手段５に
おいて以下のようにして行われる。
【００３６】図７は差信号の逆変換の詳細を表す図であ
る。まず、残留画像信号の差信号Ｓsubgが補間手段14に
より各画素の間が補間されて元の大きさの４倍の大きさ
の画像を表す画像信号ｇL ′とされる。次に加算器15に
おいてその補間された画像信号ｇL ′と最も低解像度の
差信号Ｓsub ｂn-1 の相対応する画素同志で加算を行
い、加算画像信号（ｇL ′＋Ｓsub ｂn-1 ）を得る。次
いでこの加算画像信号（ｇL ′＋Ｓsub ｂn-1 ）は補間
手段14に入力され、この補間手段14において各画素の間
が補間されて元の大きさの４倍の大きさの画像を表す画
像信号ｂn-1 ′とされる。
る。まず、残留画像信号の差信号Ｓsubgが補間手段14に
より各画素の間が補間されて元の大きさの４倍の大きさ
の画像を表す画像信号ｇL ′とされる。次に加算器15に
おいてその補間された画像信号ｇL ′と最も低解像度の
差信号Ｓsub ｂn-1 の相対応する画素同志で加算を行
い、加算画像信号（ｇL ′＋Ｓsub ｂn-1 ）を得る。次
いでこの加算画像信号（ｇL ′＋Ｓsub ｂn-1 ）は補間
手段14に入力され、この補間手段14において各画素の間
が補間されて元の大きさの４倍の大きさの画像を表す画
像信号ｂn-1 ′とされる。
【００３７】次いでこの画像ｂn-1 ′は、加算器15にお
いて差信号Ｓsub ｂn-1 の一段階高解像度の差信号Ｓsu
b ｂn-2 と相対応する画素同志の加算が行われ、加算さ
れた加算画像信号（ｂn-1 ′＋Ｓsub ｂn-2 ）は補間手
段14において各画素の間隔が補間され、差信号Ｓsub ｂ
n-2 により表される画像の４倍の大きさの画像を表す画
像信号ｂn-2 とされる。
いて差信号Ｓsub ｂn-1 の一段階高解像度の差信号Ｓsu
b ｂn-2 と相対応する画素同志の加算が行われ、加算さ
れた加算画像信号（ｂn-1 ′＋Ｓsub ｂn-2 ）は補間手
段14において各画素の間隔が補間され、差信号Ｓsub ｂ
n-2 により表される画像の４倍の大きさの画像を表す画
像信号ｂn-2 とされる。
【００３８】以上の処理を高周波の細部画像信号に対し
て順次行い、最終的に加算器15において補間画像信号ｂ
1 ′と最高解像度の差信号Ｓsub ｂ0 との加算が行わ
れ、エネルギーサブトラクション信号Ｓ′を得る。
て順次行い、最終的に加算器15において補間画像信号ｂ
1 ′と最高解像度の差信号Ｓsub ｂ0 との加算が行わ
れ、エネルギーサブトラクション信号Ｓ′を得る。
【００３９】このようにして得られたエネルギーサブト
ラクション信号Ｓ′は画像出力手段６に入力され、可視
像として表示される。この画像出力手段６はＣＲＴ等の
ディスプレイ手段でもよいし、感光フィルムに光走査記
録を行う記録装置であってもよいし、あるいはそのため
に画像信号を一旦光ディスク、磁気ディスク等の画像フ
ァイルに記憶させる装置であってもよい。
ラクション信号Ｓ′は画像出力手段６に入力され、可視
像として表示される。この画像出力手段６はＣＲＴ等の
ディスプレイ手段でもよいし、感光フィルムに光走査記
録を行う記録装置であってもよいし、あるいはそのため
に画像信号を一旦光ディスク、磁気ディスク等の画像フ
ァイルに記憶させる装置であってもよい。
【００４０】このように、本発明によるエネルギーサブ
トラクション処理方法および装置は、多重解像度空間に
変換された複数の周波数帯域の細部画像信号に対して周
波数帯域に応じて重み付け係数を変更して差信号を得る
ようにしたため、上述した実施の形態のようにノイズを
多く含む高周波数帯域についてはノイズを強調しないよ
うな重み付け係数とすることができるため、エネルギー
サブトラクション信号を可視像として再生することによ
り得られるエネルギーサブトラクション画像は、ノイズ
が低減された粒状の良い画像となる。
トラクション処理方法および装置は、多重解像度空間に
変換された複数の周波数帯域の細部画像信号に対して周
波数帯域に応じて重み付け係数を変更して差信号を得る
ようにしたため、上述した実施の形態のようにノイズを
多く含む高周波数帯域についてはノイズを強調しないよ
うな重み付け係数とすることができるため、エネルギー
サブトラクション信号を可視像として再生することによ
り得られるエネルギーサブトラクション画像は、ノイズ
が低減された粒状の良い画像となる。
【００４１】また、複数の周波数帯域ごとの周波数信号
は、画像処理手段７において行われる強調処理などの種
々の画像処理を施すために得られ、記憶手段３に記憶さ
れてエネルギーサブトラクション処理や画像処理の際に
この記憶手段３から読み出されてエネルギーサブトラク
ション処理、画像処理に供されるものである。このよう
に記憶手段３に予め記憶された細部画像信号を用いるこ
とができるため、新たな画像を獲る必要が無くなり、エ
ネルギーサブトラクション処理を行うための時間を短縮
し、高速に処理を行うことができる。
は、画像処理手段７において行われる強調処理などの種
々の画像処理を施すために得られ、記憶手段３に記憶さ
れてエネルギーサブトラクション処理や画像処理の際に
この記憶手段３から読み出されてエネルギーサブトラク
ション処理、画像処理に供されるものである。このよう
に記憶手段３に予め記憶された細部画像信号を用いるこ
とができるため、新たな画像を獲る必要が無くなり、エ
ネルギーサブトラクション処理を行うための時間を短縮
し、高速に処理を行うことができる。
【００４２】また、上述した実施の形態においては、画
像を多重解像度画像に変換するためにラプラシアンピラ
ミッドの手法を用いているが、これに限定されるもので
はなく、例えばウェーブレット変換、あるいはサブバン
ド変換等他の方法により多重解像度画像に変換するよう
にしてもよいものである。
像を多重解像度画像に変換するためにラプラシアンピラ
ミッドの手法を用いているが、これに限定されるもので
はなく、例えばウェーブレット変換、あるいはサブバン
ド変換等他の方法により多重解像度画像に変換するよう
にしてもよいものである。
【００４３】ここで、ウェーブレット変換は、周波数解
析の方法として近年開発されたものであり、ステレオの
パターンマッチング、データ圧縮等に応用がなされてい
るものである（OLIVIER RIOUL and MARTIN VETTERLI;Wa
velets and Signal Processing,IEEE SP MAGAZINE,P.14
-38,OCTOBER 1991、Stephane Mallat;Zero-Crossingsof
a Wavelet Transform,IEEE TRANSACTIONS ON INFORM
ATION THEORY,VOL.37,NO.4,P.1019-1033,JULY 1991
）。
析の方法として近年開発されたものであり、ステレオの
パターンマッチング、データ圧縮等に応用がなされてい
るものである（OLIVIER RIOUL and MARTIN VETTERLI;Wa
velets and Signal Processing,IEEE SP MAGAZINE,P.14
-38,OCTOBER 1991、Stephane Mallat;Zero-Crossingsof
a Wavelet Transform,IEEE TRANSACTIONS ON INFORM
ATION THEORY,VOL.37,NO.4,P.1019-1033,JULY 1991
）。
【００４４】このウェーブレット変換は、
【００４５】
【数１】
【００４６】なる式において信号を複数の周波数帯域ご
との周波数信号に変換するものである。すなわち、関数
ｈの周期および縮率を変化させ、原信号を移動させるこ
とによりフィルタリング処理を行えば、細かな周波数か
ら粗い周波数までの所望とする周波数に適合した周波数
信号を作成することができる。
との周波数信号に変換するものである。すなわち、関数
ｈの周期および縮率を変化させ、原信号を移動させるこ
とによりフィルタリング処理を行えば、細かな周波数か
ら粗い周波数までの所望とする周波数に適合した周波数
信号を作成することができる。
【００４７】一方、サブバンド変換は、ウェーブレット
変換のように１種類のフィルタにより２つの周波数帯域
の画像を得るのみではなく、複数種類のフィルタを用い
て複数の周波数帯域の画像を一度に得ることをも含む変
換方法である。
変換のように１種類のフィルタにより２つの周波数帯域
の画像を得るのみではなく、複数種類のフィルタを用い
て複数の周波数帯域の画像を一度に得ることをも含む変
換方法である。
【００４８】そして、このようにウェーブレット変換あ
るいはサブバンド変換により得られた複数の周波数帯域
ごとの画像に対して上述したラプラシアンピラミッドの
場合と同様に、各周波数帯域の細部画像信号における重
み付け係数を変更してエネルギーサブトラクション処理
を行うことにより、ノイズなどの不要な周波数成分を強
調することを防止し、ノイズが目立たない観察読影に適
した良好なエネルギーサブトラクション画像を得ること
ができる。
るいはサブバンド変換により得られた複数の周波数帯域
ごとの画像に対して上述したラプラシアンピラミッドの
場合と同様に、各周波数帯域の細部画像信号における重
み付け係数を変更してエネルギーサブトラクション処理
を行うことにより、ノイズなどの不要な周波数成分を強
調することを防止し、ノイズが目立たない観察読影に適
した良好なエネルギーサブトラクション画像を得ること
ができる。
【図１】本発明によるエネルギーサブトラクション処理
方法を実施するための装置のブロック図
方法を実施するための装置のブロック図
【図２】放射線画像撮影装置を表す図
【図３】放射線画像読取装置の詳細を表す図
【図４】多重解像度分解処理手段の詳細を表す図
【図５】ローパスフィルタを表す図
【図６】ラプラシアンピラミッドによる処理が施された
複数の周波数帯域ごとの細部画像信号を表す図
複数の周波数帯域ごとの細部画像信号を表す図
【図７】復元処理手段の詳細を表す図
１ 画像入力手段 ２ 多重解像度分解処理手段 ３ 記憶手段 ４ エネルギーサブトラクション処理手段 ５ 復元処理手段 ６ 画像出力手段 ７ 画像処理手段 10 フィルタリング手段 11 補間手段 12 減算器 14 補間手段 15 加算器
Claims (2)
- 【請求項１】 同一被写体を透過したそれぞれエネルギ
ー分布が異なる放射線により得られた少なくとも一部の
画像情報が互いに異なる複数の放射線画像を表す複数の
画像信号のそれぞれを、多重解像度空間に変換すること
により得られる複数の周波数帯域ごとの周波数信号に対
して、相対応する画素についての信号間で前記各周波数
信号に周波数帯域に応じた重み付け係数を乗じて減算を
行って前記複数の周波数帯域ごとの差信号を得、 該複数の周波数帯域の差信号を逆変換することによりエ
ネルギーサブトラクション信号を得ることを特徴とする
エネルギーサブトラクション処理方法。 - 【請求項２】 同一被写体を透過したそれぞれエネルギ
ー分布が異なる放射線により得られた少なくとも一部の
画像情報が互いに異なる複数の放射線画像を表す複数の
画像信号のそれぞれを、多重解像度空間に変換すること
により得られる複数の周波数帯域ごとの周波数信号を記
憶する記憶手段と、 該複数の周波数帯域ごとの周波数信号に対して相対応す
る画素についての信号間で前記各周波数信号に周波数帯
域に応じた重み付け係数を乗じて減算を行って前記複数
の周波数帯域ごとの差信号を得るサブトラクション手段
と、 該複数の周波数帯域の差信号を逆変換することによりエ
ネルギーサブトラクション信号を得る逆変換手段とから
なることを特徴とするエネルギーサブトラクション処理
装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7191524A JPH0944659A (ja) | 1995-07-27 | 1995-07-27 | エネルギーサブトラクション処理方法および装置 |
| US08/686,516 US5694447A (en) | 1995-07-27 | 1996-07-26 | Energy subtraction processing method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7191524A JPH0944659A (ja) | 1995-07-27 | 1995-07-27 | エネルギーサブトラクション処理方法および装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0944659A true JPH0944659A (ja) | 1997-02-14 |
Family
ID=16276102
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7191524A Withdrawn JPH0944659A (ja) | 1995-07-27 | 1995-07-27 | エネルギーサブトラクション処理方法および装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5694447A (ja) |
| JP (1) | JPH0944659A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004166801A (ja) * | 2002-11-18 | 2004-06-17 | Kose Corp | 肌のつやの評価方法 |
| JP2025500595A (ja) * | 2022-04-11 | 2025-01-09 | ディアールテック コーポレーション | Ｘ線撮影装置およびこれを利用するｘ線撮影方法 |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6965697B1 (en) * | 1998-07-15 | 2005-11-15 | Sony Corporation | Coding apparatus and method, decoding apparatus and method, data processing system, storage medium, and signal |
| DE60022237T2 (de) * | 1999-01-14 | 2006-07-06 | Fuji Photo Film Co., Ltd., Minami-Ashigara | Bidverarbeitungsverfahren und -System und Aufzeichnungsmedium zur Durchfürung des Verfahrens |
| US7177481B2 (en) * | 2000-12-19 | 2007-02-13 | Konica Corporation | Multiresolution unsharp image processing apparatus |
| US7181086B2 (en) * | 2002-06-06 | 2007-02-20 | Eastman Kodak Company | Multiresolution method of spatially filtering a digital image |
| US7397964B2 (en) * | 2004-06-24 | 2008-07-08 | Apple Inc. | Gaussian blur approximation suitable for GPU |
| GB0716045D0 (en) * | 2007-08-17 | 2007-09-26 | Durham Scient Crystals Ltd | Method and apparatus for inspection of materials |
| DE102007040472A1 (de) * | 2007-08-28 | 2009-03-05 | Bernd Hansen | Trennvorrichtung |
| KR101486776B1 (ko) * | 2010-07-29 | 2015-01-29 | 삼성전자주식회사 | 영상 처리 방법 및 장치와 이를 채용한 의료영상시스템 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5049748A (en) * | 1989-10-19 | 1991-09-17 | Fuji Photo Film Co. Ltd. | Method and apparatus for forming energy subtraction images |
| JP2627097B2 (ja) * | 1990-04-04 | 1997-07-02 | 富士写真フイルム株式会社 | エネルギーサブトラクション画像生成方法および装置 |
| DE69331719T2 (de) * | 1992-06-19 | 2002-10-24 | Agfa-Gevaert, Mortsel | Verfahren und Vorrichtung zur Geräuschunterdrückung |
-
1995
- 1995-07-27 JP JP7191524A patent/JPH0944659A/ja not_active Withdrawn
-
1996
- 1996-07-26 US US08/686,516 patent/US5694447A/en not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004166801A (ja) * | 2002-11-18 | 2004-06-17 | Kose Corp | 肌のつやの評価方法 |
| JP2025500595A (ja) * | 2022-04-11 | 2025-01-09 | ディアールテック コーポレーション | Ｘ線撮影装置およびこれを利用するｘ線撮影方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5694447A (en) | 1997-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3683914B2 (ja) | ピラミッド的画像分解に基づいた放射線画像の多重処理法 | |
| EP1120743B1 (en) | Method and unit for suppressing a periodic pattern | |
| EP0677780B1 (en) | Image superposition processing method | |
| EP0756247B1 (en) | Image processing method and apparatus | |
| US5907642A (en) | Method and apparatus for enhancing images by emphasis processing of a multiresolution frequency band | |
| JPH0944659A (ja) | エネルギーサブトラクション処理方法および装置 | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| JP3675896B2 (ja) | 画像処理方法および装置 | |
| JPH0312785A (ja) | 画像処理方法 | |
| JP2640582B2 (ja) | エネルギーサブトラクション画像生成方法 | |
| JP2754068B2 (ja) | 放射線画像信号生成方法および放射線画像読取装置 | |
| JP3370797B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JP3205417B2 (ja) | 画像処理方法 | |
| US6614044B2 (en) | Image signal generating method, apparatus and program | |
| JPH06274614A (ja) | 画像処理方法 | |
| JPH0944651A (ja) | 画像処理方法および装置 | |
| JPH03285475A (ja) | エネルギーサブトラクション画像生成方法 | |
| US6587594B1 (en) | Moire eliminating filter and image processing method and apparatus using the filter | |
| JP4265869B2 (ja) | 周期的パターン抑制処理方法および装置 | |
| JPH09160149A (ja) | 画像処理方法および装置 | |
| JPH0944645A (ja) | 画像処理方法および装置 | |
| JP3432917B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JP3454318B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JPH0622218A (ja) | エネルギーサブトラクション画像生成方法 | |
| JPH09161061A (ja) | 画像処理方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20021001 |