JP3002937B2 - エネルギーサブトラクション画像処理方法 - Google Patents

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、放射線画像のエネルギ
ーサブトラクションで得られた画像のノイズを低減する
ようにしたエネルギーサブトラクション画像処理方法に
関するものである。

【０００２】

【従来の技術】従来より、放射線画像のサブトラクショ
ン処理が公知となっている。この放射線画像のサブトラ
クションとは、異なった条件で撮影した複数（基本的に
は２つ）の放射線画像を読み取ってデジタル画像信号を
得た後、これらのデジタル画像信号を両画像の各画素を
対応させて減算処理（サブトラクト）し、放射線画像中
の特定の構造物を抽出させる差信号を得る方法であり、
このようにして得た差信号を用いれば、特定構造物のみ
が抽出された放射線画像を再生することができる。

【０００３】このサブトラクション処理には、基本的に
次の２つの方法がある。すなわち、(1) 造影剤注入によ
り特定の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像信号を
引き算することによって特定の構造物を抽出するいわゆ
る時間サブトラクション処理と、(2) 同一の被写体に対
して相異なるエネルギー分布を有する放射線を照射し、
あるいは被写体透過後の放射線をエネルギー分布状態を
変えて複数の放射線検出手段に照射して、それにより特
定の構造物が異なる複数の放射線画像を記録し、その後
これらの放射線画像を担持する画像信号間で適当な重み
づけをした上で引き算を行なって、特定の構造物の画像
を抽出するいわゆるエネルギーサブトラクション処理で
ある。

【０００４】一方、ある種の蛍光体に放射線（Ｘ線、α
線、β線、γ線、紫外線、電子線等）を照射すると、こ
の放射線のエネルギーの一部がその蛍光体中に蓄積さ
れ、その後その蛍光体に可視光等の励起光を照射する
と、蓄積されたエネルギーに応じて蛍光体が輝尽発光を
示す。このような性質を示す蛍光体を蓄積性蛍光体（輝
尽性蛍光体）と言う。

【０００５】この蓄積性蛍光体を利用して、人体等の被
写体の放射線画像情報を一旦蓄積性蛍光体の層を有する
シート（以下、蓄積性蛍光体シートと称する）に記録
し、これを励起光で走査して輝尽発光させ、この輝尽発
光光を光電的に読み取って画像信号を得、この画像信号
を処理して診断適性の良い被写体の放射線画像を得る方
法が提案されている（例えば特開昭５５−１２４２９
号、同５５−１１６３４０号、同５５−１６３４７２
号、同５６−１１３９５号、同５６−１０４６４５号な
ど）。

【０００６】このような蓄積性蛍光体シートを利用する
放射線画像情報記録再生システムにおいては、該シート
に記録されている放射線画像情報が直接電気的画像信号
の形で読み取られるから、このシステムによれば、上述
のサブトラクション処理を容易に行なうことが可能とな
る。また、その他、Ｘ線写真フィルム等に放射線画像を
撮影した場合でも、そのフィルム等からマイクロフォト
メータ等によって放射線画像を読み取って画像信号を得
ることができるので、上記のサブトラクション処理を実
施可能である。

【０００７】

【発明が解決しようとする課題】従来、上記エネルギー
サブトラクションで得られた画像においては、撮影時の
到達放射線量が少なかった低濃度部のノイズが多いとい
う問題が認められている。そこで、従来からよくなされ
ているようにエネルギーサブトラクション画像に対して
全体的に粒状性改善の処理を施すと、低濃度部において
高周波のノイズは低減するものの、低周波のノイズが残
るために、いわゆる「もこもこ」とした感じのアーチフ
ァクト（偽画像）が発生してしまう。

【０００８】上に述べた低濃度部は、例えば人体の胸部
画像における椎体等であるので、このアーチファクトが
実際の診断に悪影響を及ぼすことは比較的少ないが、そ
れでもこのようなアーチファクトが発生しているエネル
ギーサブトラクション画像は非常に見苦しくて、悪印象
を与えるものとなってしまう。

【０００９】本発明は上記の事情に鑑みてなされたもの
であり、エネルギーサブトラクション画像の低濃度部の
ノイズを低減し、そしてアーチファクトが発生すること
も防止できるエネルギーサブトラクション画像処理方法
を提供することを目的とするものである。

【００１０】

【課題を解決するための手段】本発明によるエネルギー
サブトラクション画像処理方法は、前述したように同一
被写体の少なくとも一部の画像情報が互いに異なる放射
線画像が記録された複数の記録媒体から各々放射線画像
情報を読み取って複数組の原画像信号を得、これら複数
組の原画像信号について、相対応する画素についての信
号間で減算を行なって、上記被写体中の特定構造物の画
像を形成する差信号を得るエネルギーサブトラクション
において、上記原画像信号を平滑化処理にかけて平滑化
画像信号を得、上記原画像信号と差信号とを、相対応す
る各画素についての信号毎に、その画素に関する上記平
滑化画像信号が示す画像濃度が低くなるにつれて原画像
信号への重み付けを大きくして加重平均し、この加重平
均によって得られた合成画像信号に基づいて上記被写体
の放射線画像を再生することを特徴とするものである。

【００１１】なお上記の「画像濃度が低くなるにつれて
原画像信号への重み付けを大きくして」ということは、
画像濃度の全範囲において画像濃度低下に伴なって重み
付けを単調増大させることのみならず、画像濃度が低下
しても重み付けを変えない領域を一部設定する場合も含
むものとする。

【００１２】

【作用および発明の効果】上記原画像信号と差信号との
加重平均によって得られる合成画像信号が担持する画像
は、いわばエネルギーサブトラクション画像に部分的に
原画像が、あるいは原画像に近い画像が組み合わされた
ものとなる。そしてこの合成画像信号は、原画像信号へ
の重み付けが大きいほど（つまり、差信号への重み付け
が小さいほど）、より原画像信号に近いものとなる。先
に述べたエネルギーサブトラクション画像の低濃度部の
ノイズは、サブトラクション処理を行なうために発生す
るものであって本来原画像には無いものであるから、こ
のような低濃度部に関する合成画像信号を上記のように
して原画像信号に近いものとすれば、この合成画像信号
により再生される画像において上記ノイズが低減するこ
とになる。

【００１３】そして、画像濃度が低くなるにつれて原画
像信号への重み付けを大きくしているので、濃度勾配が
有る画像部分の相隣接する画素において、合成画像信号
が原画像信号に近いものから急に差信号に近いものとな
るように変化することがない。したがって、そのような
合成画像信号の急激な変化のために、再生画像において
濃度段差が生じることも防止できる。

【００１４】また原画像においては、被写体構造物のエ
ッジ部分等で濃度が急激に変化することもあるが、本発
明方法においてはこのような原画像を担う信号ではな
く、原画像信号を平滑化処理した信号が示す画像濃度に
基づいて加重平均の重み付けを変えているので、上記エ
ッジ部分等において合成画像信号が原画像信号に近いも
のから急に差信号に近いものとなるように変化すること
もない。そこで、そのような合成画像信号の急激な変化
のために、再生画像において濃度段差が生じることも防
止できる。

【００１５】以上説明した通り本発明方法は、エネルギ
ーサブトラクション画像のノイズが発生しやすい部分
を、本質的にそのようなノイズを含まない原画像と、あ
るいはそれに近い画像と置き換えるという処理を行なっ
ているので、エネルギーサブトラクションに起因する低
周波のノイズが残ってアーチファクトとなることも確実
に防止できる。

【００１６】

【実施例】以下、図面を参照して、本発明の実施例につ
いて説明する。なおここでは、前述した蓄積性蛍光体シ
ートを用いる例について説明する。

【００１７】まず本発明の第１実施例について説明す
る。図２は、Ｘ線撮影装置の概略図である。このＸ線撮
影装置１のＸ線管２から発せられたＸ線３が被写体（人
体の胸部）４に照射される。被写体４を透過したＸ線３
ａは第一の蓄積性蛍光体シート５に照射され、このＸ線
３ａのうち比較的低エネルギーのＸ線が該第一の蓄積性
蛍光体シート５に蓄積され、これにより該シート５に被
写体４のＸ線画像が蓄積記録される。シート５を透過し
たＸ線３ｂはさらに低エネルギーのＸ線をカットするフ
ィルタ６を透過し、該フィルタ６を透過した高エネルギ
ーＸ線３ｃが第二の蓄積性蛍光体シート７に照射され
る。これにより該シート７にも被写体４のＸ線画像が蓄
積記録される。被写体４には、サブトラクション処理を
行なうにあたって２つのＸ線画像の位置合わせを行なう
ための基準となる２つのマーク８が付されている。

【００１８】なお、上記Ｘ線撮影装置は一回の撮影で２
枚のシート５，７にＸ線画像を蓄積記録するものである
が、時間的に相前後した２つのタイミングでそれぞれ１
枚ずつ撮影を行なってもよい。

【００１９】図３は、Ｘ線画像読取装置と本発明のエネ
ルギーサブトラクション画像処理方法を実施するための
画像処理表示装置の斜視図である。図２に示したＸ線撮
影装置１で撮影が行なわれた後、第一および第二の蓄積
性蛍光体シート５，７が一枚ずつＸ線画像読取装置10の
所定位置にセットされる。ここでは、第一の蓄積性蛍光
体シート５に蓄積記録された第一のＸ線画像（低圧画
像）の読取りについて説明する。

【００２０】所定位置にセットされた蓄積性蛍光体シー
ト５は、図示しない駆動手段によって駆動されるエンド
レスベルト等のシート搬送手段15により、励起光副走査
のために矢印Ｙ方向に搬送される。一方、レーザ光源16
から発せられた励起光としての光ビーム17は、モータ18
により駆動され矢印Ｚ方向に高速回転する回転多面鏡19
によって反射偏向され、ｆθレンズ等の集束レンズ20を
通過した後、ミラー21により光路を変えて蓄積性蛍光体
シート５に入射し、副走査の方向（矢印Ｙ方向）と略直
角な矢印Ｘ方向に主走査する。

【００２１】蓄積性蛍光体シート５の、光ビーム17が照
射された箇所からは、蓄積記録されているＸ線画像情報
に応じた光量の輝尽発光光22が発せられ、この輝尽発光
光22は光ガイド23によって導かれ、フォトマルチプライ
ヤ（光電子増倍管）24によって光電的に検出される。光
ガイド23はアクリル板等の導光性材料を成形して作られ
たものであり、直線状をなす入射端面23ａが蓄積性蛍光
体シート５上の主走査線に沿って延びるように配され、
円環状に形成された射出端面23ｂにフォトマルチプライ
ヤ24の受光面が結合されている。入射端面23ａから光ガ
イド23内に入射した輝尽発光光22は、該光ガイド23の内
部を全反射を繰り返して進み、射出端面23ｂから射出し
てフォトマルチプライヤ24に受光され、放射線画像を表
わす輝尽発光光22がフォトマルチプライヤ24によって電
気信号に変換される。

【００２２】フォトマルチプライヤ24から出力されたア
ナログの光検出信号Ｓは、ログアンプ25で対数的に増幅
された後、Ａ／Ｄ変換器26によってＡ／Ｄ変換され、デ
ィジタルの画像信号Ｓ０が得られる。この画像信号Ｓ０
は第一の蓄積性蛍光体シート５に蓄積記録されていた第
一のＸ線画像を表わすものであり、ここでは第一の画像
信号Ｓ０₁ と称する。この第一の画像信号Ｓ０₁ は画像
処理表示装置30内の内部メモリに一旦記憶される。

【００２３】この画像処理表示装置30は、種々の指示を
入力するキーボード31、指示のための補助情報や画像信
号に基づく可視画像を表示するＣＲＴディスプレイ32、
補助記憶媒体としてのフロッピィディスクを駆動するフ
ロッピィディスク駆動装置33、およびＣＰＵや内部メモ
リが内蔵された本体部34を備えている。

【００２４】次に上記と同様にして、第二の蓄積性蛍光
体シート７に蓄積記録された第二のＸ線画像（高圧画
像）を表わす第二の画像信号Ｓ０₂ が得られ、この第二
の画像信号Ｓ０₂ も画像処理表示装置30内の内部メモリ
に一旦記憶される。

【００２５】図１は、画像処理表示装置30内の内部メモ
リに記憶された第一および第二のＸ線画像を表わす２つ
の画像信号Ｓ０₁ ，Ｓ０₂ に基づいて、該画像処理表示
装置30内で行なわれる処理の流れを表わした図である。
画像処理表示装置30の内部メモリに記憶された、第一お
よび第二のＸ線画像信号Ｓ０₁ ，Ｓ０₂ は、この図１に
示す第一のＸ線画像41，第二のＸ線画像42を担持する信
号である。第一のＸ線画像41は比較的低エネルギーのＸ
線による低圧画像であり、第二のＸ線画像42は比較的高
エネルギーのＸ線による高圧画像であるが、互いに軟部
および骨部の濃度は異なるものの両者ともこれら軟部お
よび骨部の双方を記録している。

【００２６】これら第一および第二のＸ線画像信号Ｓ０
₁ ，Ｓ０₂ が画像処理表示装置30の内部メモリから読み
出されると、先ずこれら２つの画像信号Ｓ０₁ ，Ｓ０₂
がそれぞれ担持する各Ｘ線画像41，42の相対的な位置合
わせが画像信号上で行なわれる（特開昭５８−１６３３
３８号公報参照）。この位置合わせは、図２に示した２
つのマーク８が重なるように、２つのＸ線画像を相対的
に直線移動および回転移動させて行なわれる。

【００２７】この後、サブトラクション処理が行なわれ
る。ここでＸ線の吸収係数μを、被写体の軟部と骨部、
および低エネルギーＸ線と高エネルギーＸ線とに分けて
次のように定める。

【００２８】 μ_L ^T ：低エネルギーＸ線による軟部の吸収係数 μ_H ^T ：高エネルギーＸ線による軟部の吸収係数 μ_L ^B ：低エネルギーＸ線による骨部の吸収係数 μ_H ^B ：高エネルギーＸ線による骨部の吸収係数 このとき、２つの画像信号Ｓ０₁ ，Ｓ０₂ の互いに対応
する画素毎に、式

【００２９】

【数１】

【００３０】に従って重み付け引き算を行なえば、骨部
の陰影が抽出された骨部画像を表わす骨部画像信号が求
められる。

【００３１】また、式

【００３２】

【数２】

【００３３】に従って重み付け引き算を行なうことによ
り、軟部画像43を表わす軟部画像信号Ｓ２を求めること
ができる。本実施例では上記２つの演算のうち、(2) 式
の演算のみを行なって、軟部画像信号Ｓ２を求める。

【００３４】さらに、式 Ｓ０＝（Ｓ０₁ ＋Ｓ０₂ ）／２ …(3) に従って互いに対応する各画素毎に加算処理を行なうこ
とにより、２つのＸ線画像41，42の重ね合わせ画像44が
生成される。この重ね合わせ画像44も軟部および骨部の
双方が記録された原画像である。

【００３５】本発明では、この重ね合わせ画像44に代え
てＸ線画像（低圧画像）41もしくはＸ線画像（高圧画
像）42を用いることも可能であるが、重ね合わせ画像44
は２つのＸ線画像41，42を重ね合わせたものであるため
これら各Ｘ線画像のいずれと比べてもノイズ成分が低減
されており、したがってその後の処理に有利となる。た
だし、重ね合わせ画像44を用いる場合は計算量が増大す
るので、それを避けたい場合は低圧画像41（これは粒状
性の点で高圧画像42よりも優れている）を用いるのが好
ましい。

【００３６】次に上記重ね合わせ画像信号Ｓ０を平滑化
処理にかけて、平滑化画像信号Ｓ０’を得る。この平滑
化処理方法としては、例えば各画素に対し該画素を中心
とした所定領域内の各画素に関する画像信号の平均を求
め、この平均値を中心の画素の画像信号とする単純な平
均化処理（ボケマスク処理）方法、上記所定領域内の画
像信号の中央値（メジアン）を中心の画素の画像信号と
するメジアンフィルタを用いる方法、上記所定領域内を
さらに複数の小領域に分け、各小領域毎に分散を求めて
分散の最も小さい小領域の平均値を中央の画素の画像信
号の値とするエッジ保存フィルタ（Ｖ−フィルタ）を用
いる方法、画像信号をフーリエ変換し、ノイズ成分に対
応する高空間周波数成分を取り除いた後に逆フーリエ変
換する方法等を用いることができる。

【００３７】それとともに、上記重ね合わせ画像信号Ｓ
０および軟部画像信号Ｓ２の、互いに対応する画素につ
いての信号間で、式 Ｎ＝Ｔ（Ｌｕ）・Ｓ０＋｛１−Ｔ（Ｌｕ）｝・Ｓ２ …(4) に従って加重平均処理がなされる。ここでＴ（Ｌｕ）
は、サブトラクション処理前の原画像信号の一種である
重ね合わせ画像信号Ｓ０に対する重み係数、１−Ｔ（Ｌ
ｕ）は軟部画像信号Ｓ２に対する重み係数である。

【００３８】上記重み係数Ｔ（Ｌｕ）は、前述した平滑
化画像信号Ｓ０’が示す画素値（画像濃度を示すディジ
タル値）Ｌｕに基づいて、例えば図４に示すような特性
に従って決定される。つまりこの重み係数Ｔ（Ｌｕ）
は、画素値Ｌｕが最大値Ｌ３から所定値Ｌ２までの間に
ある場合はゼロとされ、画素値Ｌｕが所定値Ｌ２を下回
る場合はその値の低下つまり画像濃度の低下につれてよ
り大きく設定され、画素値Ｌｕが最小値Ｌ１のとき最大
の１とされるものである。

【００３９】なお具体的には、以上のような画素値Ｌｕ
と重み係数Ｔ（Ｌｕ）との関係をテーブルの形で記憶手
段に記憶しておき、平滑化画像信号Ｓ０’が示す画素値
Ｌｕに従ってテーブル参照することにより重み係数Ｔ
（Ｌｕ）を決定するのが便利である。

【００４０】以上の加重平均処理で得られた合成画像信
号Ｎに基づいて、画像処理表示装置30のＣＲＴディスプ
レイ32において該信号Ｎが担持する放射線画像が再生表
示される。この合成画像45は、サブトラクション処理の
ためにノイズが発生しやすい軟部画像43の椎体部Ｑ等の
低濃度部が、濃度が低いほど重ね合わせ画像44への比率
を高くして、該重ね合わせ画像44と軟部画像43との中間
的な画像と（あるいは重ね合わせ画像44そのものと）置
き換えられた形のものとなる。このような合成画像45に
おいては、濃度段差を生じることなく上記ノイズの発生
が抑えられるようになる。その理由は、先に詳しく説明
した通りである。

【００４１】次に図５を参照して、本発明による第２実
施例について説明する。なお本実施例および後述の第３
実施例は、サブトラクション画像である軟部画像の求め
方が第１実施例と異なるものであり、その他の処理は本
質的に第１実施例と同様に行なえばよいので、以下は軟
部画像を求める点について説明する。またこの図５にお
いて、図１中と同等の要素には図１に付した番号，記号
を付し、図１を用いて説明した部分については重複説明
を避けるために、ここでの説明は省略する。

【００４２】図５に示される通りこの第２実施例におい
ては、前述の(1) 式に従ってサブトラクション処理がな
され、骨部の陰影が抽出された骨部画像46を表わす骨部
画像信号Ｓ１が求められる。またそれとともに、第１実
施例におけるのと同様にして、２つのＸ線画像41，42を
重ね合わせた重ね合わせ画像44を表わす重ね合わせ画像
信号Ｓ０が求められる。

【００４３】次に、骨部画像信号Ｓ１に平滑化処理が施
され、平滑化骨部画像47を担持する平滑化骨部画像信号
Ｓ１’が求められる。この平滑化処理方法としては、先
に説明した重ね合わせ画像信号Ｓ０に対する平滑化処理
と同様に、平均化処理（ボケマスク処理）方法、メジア
ンフィルタを用いる方法、エッジ保存フィルタ（Ｖ−フ
ィルタ）を用いる方法、フーリエ変換および逆フーリエ
変換する方法等を用いることができるが、例えば特開平
３−２８５４７５号公報に示されるヒストグラム適応フ
ィルタを用いた平滑化処理を用いると、画像情報として
必要なエッジ（互いに異なる２つの組織の陰影の境界を
定めるステップ状の濃度変化）を保存したままかつアー
チファクトなしにノイズを除去することができ、また簡
単な演算で短時間にノイズを除去することができる。

【００４４】次に、各画素毎に重ね合わせ画像信号Ｓ０
と平滑化画像信号Ｓ１’との間で重みづけ引き算、すな
わちサブトラクション処理がなされる。それにより、画
像情報としては前記(2) 式で表わされる軟部画像信号Ｓ
２と略同一の情報を担持する一方、この(2) 式で表わさ
れる軟部画像信号Ｓ２よりもノイズ成分が低減された処
理済軟部画像信号Ｓ２’が求められる。

【００４５】この処理済軟部画像43’を表わす画像信号
Ｓ２’を、図１に示した軟部画像信号Ｓ２に代えて、重
ね合わせ画像信号Ｓ０との加重平均処理にかければ、第
１実施例におけるのと同様に、低濃度部のノイズが低減
された軟部画像を担持する合成画像信号が得られる。

【００４６】なおこの実施例では、２回目のサブトラク
ション処理のために重ね合わせ画像信号Ｓ０が求められ
るので、この重ね合わせ画像信号Ｓ０を上記加重平均処
理および、重み係数Ｔ（Ｌｕ）を求める処理にそのまま
利用することができる。

【００４７】次に図６を参照して、本発明の第３実施例
について説明する。この実施例ではまず、Ｘ線画像41，
42をそれぞれ担持する画像信号Ｓ０₁ ，Ｓ０₂ から、前
述の(1) および(2) 式に基づいて、骨部画像46を示す骨
部画像信号Ｓ１と、軟部画像43を示す軟部画像信号Ｓ２
が求められる。

【００４８】次に、骨部画像信号Ｓ１を第２実施例と同
様にして平滑化処理することにより、骨部画像46に含ま
れるノイズ成分が低減された平滑化骨部画像信号Ｓ１’
が求められ、その後各画素毎に骨部画像信号Ｓ１から平
滑化骨部画像信号Ｓ１’を引き算することにより、ノイ
ズ成分のみが抽出されたノイズ画像48を示すノイズ信号
Ｓ_N が求められる。

【００４９】ここで平滑化骨部画像信号Ｓ１’において
は、骨部画像のエッジの情報はたとえノイズ成分と同程
度の高空間周波数であっても保存されているため、骨部
画像信号Ｓ１と平滑化骨部画像信号Ｓ１’との差を求め
ることによりエッジの情報はきれいにキャンセルされ、
したがってエッジの情報を失わせるような平滑化処理を
行なった場合と比べ、ノイズ信号Ｓ_N はより純粋に骨部
画像のノイズ成分のみを担持した信号となる。

【００５０】次にこのようにして求められたノイズ信号
Ｓ_N と、軟部画像43を表わす軟部画像信号Ｓ２とが各画
素毎に重み付け加算され、これにより、画像情報として
は上記軟部画像43と略同一の情報を担持する一方、該軟
部画像43よりもノイズ成分が低減された処理済軟部画像
43’を示す処理済軟部画像信号Ｓ２’が求められる。

【００５１】この処理済軟部画像43’を表わす画像信号
Ｓ２’を、図１に示した軟部画像信号Ｓ２に代えて、重
ね合わせ画像信号Ｓ０との加重平均処理にかければ、第
１実施例におけるのと同様に、低濃度部のノイズが低減
された軟部画像を担持する合成画像信号が得られる。

【００５２】以上、本発明における加重平均処理にかけ
る前のサブトラクション画像を、予め全体的なノイズ低
減処理にかけるようにした２つの実施例について説明し
たが、このようなノイズ低減処理は例えば前記特開平３
−２８５４７５号公報等によりそれ以外のものが種々公
知となっており、本発明においてはそれら公知のノイズ
低減処理も適宜併せて実行可能である。

【００５３】また、以上の各実施例はいずれも人体の胸
部のＸ線画像に基づいて軟部画像もしくは骨部画像を求
める例であるが、本発明は軟部画像もしくは骨部画像を
求める場合のみを対象とするものではなく、例えば、乳
腺が強調された画像もしくは悪性腫瘍が強調された画像
等にも適用可能であり、一般に被写体中の互いに異なる
２つの組織がそれぞれ強調もしくは抽出された２つの画
像のうちの一方、もしくは双方を求める際に広く適用す
ることができるものである。

【００５４】さらに、上記実施例は、蓄積性蛍光体シー
トを用いた例であるが、本発明は蓄積性蛍光体シートを
用いる場合に限られるものではなく、Ｘ線フイルム（撮
影に際して一般に増感スクリーンと組合わされる）等を
用いる場合にも適用することができる。

【図面の簡単な説明】

【図１】本発明によるサブトラクション画像処理の流れ
の一例を示す概略図

【図２】Ｘ線撮影装置の一例を示す概略図

【図３】本発明方法を実施する装置の一例を示す斜視図

【図４】平滑化画像信号が示す画素値Ｌｕと、重み係数
Ｔ（Ｌｕ）との関係の一例を示すグラフ

【図５】本発明方法におけるサブトラクション処理の流
れの一例を示す概略図

【図６】本発明方法におけるサブトラクション処理の流
れの別の例を示す概略図

【符号の説明】

１ Ｘ線撮影装置 ２ Ｘ線管 ３，３ａ，３ｂ，３ｃ Ｘ線 ４ 被写体 ５ 第一の蓄積性蛍光体シート ６ フィルタ ７ 第二の蓄積性蛍光体シート ８ マーク 10 Ｘ線画像読取装置 15 シート搬送手段 16 レーザ光源 17 光ビーム 19 回転多面鏡 22 輝尽発光光 23 光ガイド 24 フォトマルチプライヤ 25 ログアンプ 26 Ａ／Ｄ変換器 30 画像処理表示装置 41 Ｘ線画像（低圧画像） 42 Ｘ線画像（高圧画像） 43 軟部画像 43’ ノイズ成分低減処理済軟部画像 44 重ね合わせ画像 45 合成画像 46 骨部画像 47 平滑化骨部画像 48 ノイズ画像 Ｓ０₁ ，Ｓ０₂ Ｘ線画像信号 Ｓ０ 重ね合わせ画像信号 Ｓ０’ 平滑化重ね合わせ画像信号 Ｓ１ 骨部画像信号 Ｓ１’ 平滑化骨部画像信号 Ｓ２ 軟部画像信号 Ｓ２’ ノイズ成分低減処理済軟部画像信号 Ｎ 合成画像信号

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.⁷，ＤＢ名) G06T 1/00 G06T 5/00 H04N 1/409 A61B 6/00 - 6/14

Claims (1)

(57)【特許請求の範囲】
1. 【請求項１】 同一被写体の少なくとも一部の画像情報
   が互いに異なる放射線画像が記録された複数の記録媒体
   から各々放射線画像情報を読み取って複数組の原画像信
   号を得、 これら複数組の原画像信号について、相対応する画素に
   ついての信号間で減算を行なって、前記被写体中の特定
   構造物の画像を形成する差信号を得るエネルギーサブト
   ラクションにおいて、 前記原画像信号を平滑化処理にかけて平滑化画像信号を
   得、 前記原画像信号と差信号とを、相対応する各画素につい
   ての信号毎に、その画素に関する前記平滑化画像信号が
   示す画像濃度が低くなるにつれて原画像信号への重み付
   けを大きくして加重平均し、 この加重平均によって得られた合成画像信号に基づいて
   前記被写体の放射線画像を再生することを特徴とするエ
   ネルギーサブトラクション画像処理方法。