JPH04318775A - エネルギーサブトラクション画像生成方法 - Google Patents

エネルギーサブトラクション画像生成方法

Info

Publication number
JPH04318775A
JPH04318775A JP3086881A JP8688191A JPH04318775A JP H04318775 A JPH04318775 A JP H04318775A JP 3086881 A JP3086881 A JP 3086881A JP 8688191 A JP8688191 A JP 8688191A JP H04318775 A JPH04318775 A JP H04318775A
Authority
JP
Japan
Prior art keywords
image
energy
subtraction
radiation
image data
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.)
Granted
Application number
JP3086881A
Other languages
English (en)
Other versions
JP2640582B2 (ja
Inventor
志村一男
Kazuo Shimura
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP3086881A priority Critical patent/JP2640582B2/ja
Priority to US07/866,469 priority patent/US5301107A/en
Publication of JPH04318775A publication Critical patent/JPH04318775A/ja
Application granted granted Critical
Publication of JP2640582B2 publication Critical patent/JP2640582B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Image Analysis (AREA)
  • Conversion Of X-Rays Into Visible Images (AREA)
  • Radiography Using Non-Light Waves (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は放射線画像のエネルギー
サブトラクション画像のノイズを低減させ、観察性能の
優れた画像を得るエネルギーサブトラクション画像生成
方法に関するものである。
【0002】
【従来の技術】記録された放射線画像を読み取って画像
データを得、この画像データに適切な画像処理を施した
後、画像を再生記録することが種々の分野で行われてい
る。たとえば、後の画像処理に適合するように設計され
たガンマ値の低いX線フイルムを用いてX線画像を記録
し、このX線画像が記録されたフイルムからX線画像を
読み取って電気信号に変換し、この電気信号(画像デー
タ)に画像処理を施した後コピー写真等に可視像として
再生することにより、コントラスト,シャープネス,粒
状性等の画質性能の良好な再生画像を得ることの出来る
システムが開発されている(特公昭61−5193号公
報参照)。
【0003】また本出願人により、放射線(X線,α線
,β線,γ線,電子線,紫外線等)を照射するとこの放
射線エネルギーの一部が蓄積され、その後可視光等の励
起光を照射すると蓄積されたエネルギーに応じた光量の
輝尽発光光を放射する蓄積性蛍光体(輝尽性蛍光体)を
利用して、人体等の被写体の放射線画像を一旦シート状
の蓄積性蛍光体に撮影記録し、蓄積性蛍光体シートをレ
ーザ光等の励起光で走査して輝尽発光光を生ぜしめ、得
られた輝尽発光光を光電的に読み取って画像信号を得、
この画像信号に基づいて被写体の放射線画像を写真感光
材料等の記録材料、CRT等に可視像として出力させる
放射線記録再生システムがすでに提案されている(特開
昭55−12429号,同56−11395号,同55
−0163472号,同56−164645 号,同5
5−116340 号等)。
【0004】このシステムは、従来の銀塩写真を用いる
放射線写真システムと比較して極めて広い放射線露光域
にわたって画像を記録し得るという実用的な利点を有し
ている。すなわち、放射線露光量に対する、蓄積後に励
起によって発光する輝尽発光光の光量が極めて広い範囲
に渡って比例することが認められており、従って種々の
撮影条件により放射線露光量がかなり大幅に変動しても
、蓄積性蛍光体シートより放射される輝尽発光光を読取
りゲインを適当な値に設定して光電変換手段により読み
取って電気信号(画像データ)に変換し、この画像デー
タを用いて写真感光材料、CRT等の表示装置に放射線
画像を可視像として出力することによって、放射線露光
量の変動に影響されない放射線画像を得ることができる
。
【0005】上記のようにX線フイルムや蓄積性蛍光体
シートを用いるシステムにおいて、記録された複数の放
射線画像を読み取って複数の画像データを得た後、これ
らの画像データに基づいて放射線画像のサブトラクショ
ン処理を施すことがある。
【0006】ここで、放射線画像のサブトラクション処
理とは、互いに異なった条件で撮影された複数の放射線
画像の差に対応する画像を得る処理をいい、具体的には
これら複数の放射線画像を所定のサンプリング間隔で読
み取って各放射線画像に対応する複数のディジタルの画
像信号を得、これら複数のディジタルの画像信号の各対
応するサンプリング点毎に減算処理を施すことにより、
放射線画像中の特定の被写体部分のみを強調または抽出
した放射線画像を得る処理をいう。
【0007】このサブトラクション処理には基本的には
次の二つがある。すなわち、造影剤の注入により被写体
の特定の部分(たとえば人体を被写体としたときの血管
等)が強調された放射線画像から造影剤が注入されてい
ない放射線画像を引き算(サブトラクト)することによ
って被写体の特定の部分(たとえば血管等)を抽出する
いわゆる時間サブトラクションと、被写体の特定の部分
が互いに異なるエネルギーを有する放射線に対して異な
る放射線吸収率を有することを利用して、同一の被写体
に対して互いに異なるエネルギーを有する放射線を照射
してこれら互いに異なるエネルギーを有する各放射線に
よる複数の放射線画像を得、これら複数の放射線画像を
適当に重み付けしてその差を演算することによって被写
体の特定部分を抽出するいわゆるエネルギーサブトラク
ションとがある。本出願人も蓄積性蛍光体シートを用い
たエネルギーサブトラクションについて種々提案してい
る(特開昭59−83486号公報,特開昭60−22
5541号公報参照)。
【0008】しかし上記エネルギーサブトラクション処
理後の画像は、処理前の複数の放射線画像(以下エネル
ギーサブトラクション処理前の放射線画像を「原画像」
と称する。)を減算処理することにより得られた画像で
あるため、原画像と比べるとS/N比が低下して見にく
い画像となってしまうという問題点がある。
【0009】例えば人体の胸部のように軟部および骨部
から構成された被写体に互いにエネルギーの異なる放射
線を照射して複数の放射線画像を得、これら複数の放射
線画像を読み取ってこれら複数の放射線画像のそれぞれ
を表わす複数の画像データを得、これら複数の画像デー
タに基づいてエネルギーサブトラクション処理を行なっ
て被写体の主として軟部が記録された軟部画像を表わす
軟部画像データもしくは被写体の主として骨部が記録さ
れた骨部画像を表わす骨部画像データを求め、求められ
た軟部画像もしくは骨部画像を観察の対象とする場合が
ある。この軟部画像,骨部画像はそれぞれ骨部,軟部の
陰影が消去された画像であるため、骨部もしくは軟部に
隠れてしまっていた陰影や骨部もしくは軟部の影響で見
にくくなってしまっていた陰影を浮かび上がらせること
ができ、所定の観察目的に良くマッチングする場合があ
る。しかし前述したように、これら軟部画像,骨部画像
はサブトラクション処理により得られた画像であるため
、原画像と比べノイズ成分が強調され、この点からはむ
しろ観察適性が悪化してしまっていた。
【0010】そこで本出願人は、サブトラクション処理
前の原画像とほぼ同程度にまでノイズを低減した観察適
性の優れたサブトラクション画像を生成する方法を提供
する方法として、◆エネルギーの異なる(高圧と低圧)
放射線を、該エネルギーの異なる放射線に対して放射線
吸収率の異なる複数の組織(例えば骨部と軟部組織)を
含む被写体に照射して得られた複数の放射線画像(高圧
画像と低圧画像)のそれぞれを表わす複数の原画像デー
タ(高圧画像データと低圧画像データ)に基づいて、前
記被写体中の主として第一の組織(例えば骨部)が記録
された第一の画像を表わす第一の画像データ(骨部画像
データ)を求め、◆前記第一の画像データ(骨部画像デ
ータ)を処理することにより前記第一の画像のノイズ成
分が低減又は除去された第一の平滑化画像を表わす第一
の平滑化画像データを求め、◆前記原画像データから前
記平滑化画像データを減算処理することにより、前記被
写体の主として第二の組織(例えば軟部組織)が記録さ
れた第二の画像を表わす第二の画像データ(軟部組織画
像データ)を求める方法を提案した。(特願平2−94
653 号)この方法は、その実施に際しさらに細かな
ステップに分解したり順序を変更したりして、表面的に
は各種変形した方法を採用してこの方法と実質的同一の
方法を実現することができる。
【0011】例えば、表現上は異なるものの、この方法
と実質同一の方法の例として、エネルギーの異なる(高
圧と低圧)放射線を該エネルギーの異なる放射線に対し
て放射線吸収率の異なる複数の組織(例えば骨部と軟部
組織)を含む被写体に照射して得られた複数の放射線画
像(高圧画像と低圧画像)のそれぞれを表わす複数の原
画像データ(高圧画像データと低圧画像データ)に基づ
いて、被写体中の主として第一の組織(例えば骨部)が
記録された第一の画像を表わす第一の画像データ(骨部
画像データ)と、前記被写体中の主として第二の組織(
例えば軟部組織)が記録された第二の画像を表わす第二
の画像データ(軟部組織画像データ)とを求め、前記第
一の画像データ(骨部画像データ)を処理して前記第一
の画像の主に低空間周波数成分となる前記第一の組織成
分が低減又は除去されたノイズ画像を表わすノイズ画像
データを求め、前記第二の画像データ(軟部組織画像デ
ータ)と前記ノイズ画像データとを加算処理することに
より前記被写体中の主として第二の組織(軟部組織)が
記録された新たな第二の画像を表わす新たな第二の画像
データを求めるという態様が挙げられる。この態様は、
上記第一のエネルギーサブトラクション画像生成方法と
実質的に同一の演算を行なうものである。
【0012】また、上記方法と同類の第二のエネルギー
サブトラクション画像生成方法として、◆互いにエネル
ギーの異なる複数の放射線を該複数の放射線に対して互
いに放射線吸収率の異なる複数の組織から構成される被
写体に照射して得られた複数の放射線画像のそれぞれを
表わす複数の原画像データに基づいて、前記被写体中の
主として第一の組織が記録された第一の画像を表わす第
一の画像データを求める第一の処理を行なった後、◆前
記第一の画像データを処理することにより前記第一の画
像のノイズ成分が低減された第一の平滑化画像を表わす
第一の平滑化画像データを求めて、前記原画像データか
ら該第一の平滑化画像データを減算処理することにより
、前記被写体の主として第二の組織が記録された第二の
画像を表わす第二の画像データを求める第二の処理を行
ない、該第二の処理の後、◆前記第二の画像データを処
理することにより前記第二の画像のノイズ成分が低減さ
れた第二の平滑化画像を表わす第二の平滑化画像データ
を求めて、前記原画像データから該第二の平滑化画像デ
ータを減算処理することにより、前記被写体の主として
前記第一の組織が記録された新たな第一の画像を表わす
新たな第一の画像データを求める第三の処理を行なう方
法が挙げられる。
【0013】ここで前記第二のエネルギーサブトラクシ
ョン画像生成方法における前記第二の処理、第三の処理
を繰り返し行ない、これによりさらに画質性能の良好な
画像を得ることも可能である。即ち、上記方法の第三の
エネルギーサブトラクション画像生成方法は、前記第二
のエネルギーサブトラクション画像生成方法における各
処理を行なった後、前記第三の処理により求められた前
記新たな第一の画像データを前記第二の処理における前
記第一の画像データとして再度前記第二の処理を行なう
ことにより、前記被写体の主として前記第二の組織が記
録された新たな第二の画像を表わす新たな第二の画像デ
ータを求める新たな第二の処理と、該新たな第二の画像
データを前記第三の処理における前記第二の画像データ
として再度前記第三の処理を行なうことにより、前記被
写体の主として前記第一の組織が記録された新たな第一
の画像を表わす新たな第一の画像データを求める新たな
第三の処理とを一回または複数回繰り返すことを特徴と
するものである。
【0014】また、前記第二もしくは第三のエネルギー
サブトラクション画像生成方法を適用して、最終的に被
写体の第二の組織が記録された第二の画像を表わす第二
の画像データを求めることも可能である。即ち、この方
法の第四のエネルギーサブトラクション画像生成方法は
、前記第二もしくは第三のエネルギーサブトラクション
画像生成方法における処理を行なった後、前記第三の処
理もしくは前記新たな第三の処理により求められた前記
新たな第一の画像データを前記第二の処理もしくは前記
新たな第二の処理における前記第一の画像データとして
再度前記第二の処理もしくは前記新たな第二の処理を行
なうことにより、前記被写体の主として前記第二の組織
が記録された新たな第二の画像を表わす新たな第二の画
像データを求めることを特徴とするものである。
【0015】ここで上記第二〜第四のエネルギーサブト
ラクション画像生成方法には、上記第一のエネルギーサ
ブトラクション画像生成方法と同様なステップが含まれ
ており、したがって上記第一のエネルギーサブトラクシ
ョン画像生成方法について説明したと同様に、上記第二
〜第四のエネルギーサブトラクション画像生成方法も実
質同一の各種の態様を包含する概念として把握されるも
のである。また、実質同一を含めた上記各方法を含みさ
えすれば上記方法の概念に含まれることは当然であり、
例えばこの方法を実施する前に他の方法によるノイズ低
減化処理等のステップが含まれていてもよく、この方法
を実施した後にさらにノイズを低減化させるための他の
ステップが含まれていてもよい。
【0016】なお、上記各エネルギーサブトラクション
画像生成方法における「第一の画像」(上記「新たな第
一の画像」を含む)と上記「第二の画像」(上記「新た
な第二の画像」を含む)はエネルギーサブトラクション
処理により得られた、同一被写体の互いに異なる組織の
陰影が強調もしくは抽出された2つの画像をいい、特定
のものに限定されるものではないが、例えば前述した軟
部画像と骨部画像や、人体の乳房を被写体とした場合の
、乳腺が強調された画像と悪性腫瘍が強調された画像等
をいう。
【0017】上記各方法は、エネルギーサブトラクショ
ン処理により得られた画像は、サブトラクション処理(
減算処理)を行なったためにS/N比が低下したもので
あることに着目して提案されたものである。
【0018】即ち、上記方法の第一のエネルギーサブト
ラクション画像生成方法は、サブトラクション処理によ
り得られる2つの画像(第一の画像および第二の画像)
のうちの一方の第一の画像を求め、この画像にノイズ部
分が低減または除去された第一の平滑化画像を求め、原
画像からこの第一の平滑化画像を減算処理するようにし
たため、もとの原画像と同程度までノイズが低減され観
察適性の優れた第二の画像が生成される。
【0019】ここで、高画質の第二の画像を得るために
は、上記第一の平滑化画像を求めるに際して被写体の第
一の組織の陰影は保存したままノイズ成分のみを削除す
る必要がある。ところが第一の組織の陰影とノイズ成分
とはその空間周波数成分の一部が互いに重なっており、
このため極力ノイズ成分のみを取り除くような非線形フ
ィルタを用いたとしても第一の組織の陰影とノイズ成分
との完全な分離にはおのずから限界がある。そこで上記
第二〜第四のエネルギーサブトラクション画像生成方法
は一回のノイズ低減化処理でノイズを完全に分離するこ
とを放棄し、ノイズ低減化処理を繰り返し行なうことに
よりノイズが低減化された観察適性の優れた画像を生成
するものである。
【0020】即ち、上記の第二のエネルギーサブトラク
ション画像生成方法は、第一の画像データを処理するこ
とによりノイズ成分を低減し、その後第二の画像データ
を求めて該第二の画像データを処理することによりさら
にノイズ成分を低減することにより新たな第一の画像デ
ータを求めるようにしたものであり、上記2回のノイズ
低減化処理においてそれぞれ得意とするノイズ成分の低
減化を行なうことができ、上記第一のエネルギーサブト
ラクション画像生成方法よりもさらにノイズが低減され
た、さらに観察適性の優れた画像が生成される。
【0021】また、上記方法の第三のエネルギーサブト
ラクション画像生成方法は、上記第二のエネルギーサブ
トラクション画像生成方法を繰り返すことによりさらに
ノイズ低減化を図るものであり、各ノイズ低減化処理に
おいて互いに得意とするノイズ低減化処理を分担させる
ことができ、さらにノイズが低減された画像が生成され
る。
【0022】また、上記方法の第四のエネルギーサブト
ラクション画像生成方法は、上記第二もしくは第三のエ
ネルギーサブトラクション画像生成方法を行なった後、
この第二もしくは第三のエネルギーサブトラクション画
像生成方法により得られた新たな第一の画像データにつ
いてノイズ低減化処理を行なって原画像との減算処理を
行なうようにしたため、ノイズ成分の低減化された新た
な第二の画像が生成される。
【0023】これらの方法によるエネルギーサブトラク
ションを総称して、以下、ここでは粒状改善エネルギー
サブトラクション画像生成方法あるいは粒状改善エネル
ギーサブトラクション(または粒状改善エネサブ)とい
うことにする。
【0024】
【発明が解決しようとする課題】ところで、上記の各方
法により粒状改善エネルギーサブトラクションを行なう
ときは、放射線エネルギーを分離させるために例えばX
線フイルムや蓄積性蛍光体シート等の記録シートを2枚
使用する場合、2枚の記録シートの間に低エネルギーの
放射線をカットするフィルタを介在させて放射線を照射
し、高エネルギーの放射線のみが下の記録シートに達す
るようにしている。この放射線をカットするフィルタと
しては、銅板等の金属のフィルタを使用することが多い
が、その代りに蓄積性蛍光体シートをフィルタとして使
用することもできる。すなわち、蓄積性蛍光体シートも
多少の低エネルギー放射線をカットする作用を有してい
るので、蓄積性蛍光体シートをフィルタの代用として用
いることができる。
【0025】このとき、2枚の記録シートの間に介在さ
れる蓄積性蛍光体シートは、単なるフィルタでなく、そ
れ自身も画像データを蓄積することができるので、その
蓄積した画像データを利用すればさらにエネルギーサブ
トラクション画像の粒状を改善することができることが
想像される。
【0026】本発明は、蓄積性蛍光体シートを低エネル
ギー放射線をカットするフィルタとして使用するエネル
ギーサブトラクション方法において、該フィルタとして
使用される蓄積性蛍光体シートに蓄積した画像データを
利用して粒状をさらに改善したエネルギーサブトラクシ
ョン画像生成方法を提供することを目的とするものであ
る。
【0027】
【課題を解決するための手段】本発明によるエネルギー
サブトラクション画像生成方法は、蓄積性蛍光体シート
を低エネルギー放射線をカットするフィルタとして使用
し、このフィルタを挟む2つの記録シートに記録された
放射線画像から得られた画像データの間で減算処理をし
てサブトラクション画像を得る粒状改善エネサブにおい
て、この蓄積性蛍光体シートに蓄積された画像データの
高周波成分を前記エネルギーサブトラクション画像に加
算して粒状をさらに改善したエネルギーサブトラクショ
ン画像を得ることを特徴とするものである。
【0028】すなわち、本発明によるエネルギーサブト
ラクション画像生成方法は、放射線エネルギーを分離さ
せるために上下2枚の記録シートの間に低エネルギーの
放射線をカットするフィルタを介在させて放射線を照射
し、高エネルギーの放射線のみが下の記録シートに達す
るようにし、該両記録シートに記録された放射線画像を
読み取って得られた画像データの間で減算処理をしてサ
ブトラクション画像を得るエネルギーサブトラクション
方法であって、高周波成分については減算処理をせず、
低周波成分についてのみ減算処理をして粒状を改善した
エネルギーサブトラクション画像を得るエネルギーサブ
トラクション方法において、◆蓄積性蛍光体シートを低
エネルギー放射線をカットするフィルタとして使用し、
この蓄積性蛍光体シートに蓄積された画像データの高周
波成分を前記エネルギーサブトラクション画像に加算し
て粒状を改善したエネルギーサブトラクション画像を得
ることを特徴とするものである。
【0029】前記蓄積性蛍光体シートに蓄積された画像
データの高周波成分をエネルギーサブトラクション画像
に加算するには、エネルギーサブトラクションにより得
られた画像には高周波を低減する周波数フィルタをかけ
、蓄積性蛍光体シートに蓄積された画像データには高周
波透過フィルタをかけて、両フィルタを通した画像デー
タを加算すればよい。
【0030】このとき、両フィルタのレスポンスの和が
、全周波数領域に亘って1となるようにすれば、エネル
ギーサブトラクション画像の周波数特性はそのまま維持
される。
【0031】
【作用および効果】本発明によるエネルギーサブトラク
ション画像生成方法は、高周波成分についてはサブトラ
クションの減算処理をせず低周波成分についてのみ減算
処理をして粒状を改善したエネルギーサブトラクション
画像を得る粒状改善エネルギーサブトラクション方法に
おいて、蓄積性蛍光体シートを低エネルギー放射線をカ
ットするフィルタとして使用し、この蓄積性蛍光体シー
トに蓄積された画像データの高周波成分を前記エネルギ
ーサブトラクション画像に加算するものであるから、エ
ネルギーサブトラクション画像の粒状をさらに改善する
ことができる。
【0032】
【実施例】以下、図面を参照して、本発明の基本的プロ
セスについて説明する。
【0033】図1は本発明の基本的プロセスを示すブロ
ック図である。図1に示すように、本発明によるエネル
ギーサブトラクション画像生成方法は、低エネルギー放
射線をカットするフィルタとして使用された蓄積性蛍光
体シート6(図10参照)を挟む2つの記録シート5、
7に記録された放射線画像、すなわち高圧画像5Aと低
圧画像5Bから得られた画像データ(高圧画像データと
低圧画像データ)の間で減算処理をしてサブトラクショ
ン画像11を得る粒状改善エネルギーサブトラクション
において、この蓄積性蛍光体シート6に蓄積された画像
データ(中圧画像データ)6Aの高周波成分を前記粒状
改善エネサブ画像に加算して粒状をさらに改善したエネ
ルギーサブトラクション画像を得るものである。
【0034】すなわち、本発明によるエネルギーサブト
ラクション画像生成方法は、図10に示すように、放射
線エネルギーを分離させるために上下2枚の記録シート
5、7の間に低エネルギーの放射線をカットするフィル
タ6を介在させて放射線3を照射し、高エネルギーの放
射線のみが下の記録シート7に達するようにし、両記録
シート5、7に記録された放射線画像を読み取って得ら
れた画像データの間で減算処理をしてサブトラクション
画像を得るエネルギーサブトラクション方法であって、
高周波成分については減算処理をせず、低周波成分につ
いてのみ減算処理をして粒状を改善したエネルギーサブ
トラクション画像を得る粒状改善エネルギーサブトラク
ション方法において、蓄積性蛍光体シートを低エネルギ
ー放射線をカットするフィルタ6として使用し、この蓄
積性蛍光体シート6に蓄積された画像データ(中圧画像
データ)6Aの高周波成分を前記エネルギーサブトラク
ション画像に加算して粒状を改善したエネルギーサブト
ラクション画像を得るものである。
【0035】前記蓄積性蛍光体シート6に蓄積された画
像データ(中圧画像データ)の高周波成分をエネルギー
サブトラクション画像に加算するには、図1に示すよう
に、エネルギーサブトラクションにより得られた画像1
1には高周波を低減する高域低減フィルタ12をかけ、
蓄積性蛍光体シート6に蓄積された画像データ6Aには
高域透過フィルタ13をかけて、両フィルタ12、13
を通した画像データを加算すればよい。
【0036】このとき、両フィルタ12、13のレスポ
ンスの和が、図1の下に示すように、全周波数領域に亘
って1となるようにすれば、エネルギーサブトラクショ
ン画像の周波数特性はそのまま維持される。
【0037】以下、本発明の前提となる粒状改善エネル
ギーサブトラクションの例について詳細に説明する。な
お、ここでは前述した蓄積性蛍光体シートを用いる例に
ついて説明する。
【0038】図10は、X線撮影装置の概略図である。 このX線撮影装置1のX線管2から発せられたX線3に
より被写体(人体の胸部)4が照射される。被写体4を
透過したX線3aは第一の蓄積性蛍光体シート5に照射
され、X線3aのエネルギーのうち比較的低エネルギー
のX線が該第一の蓄積性蛍光体シート5に蓄積され、こ
れにより該シート5に被写体4のX線画像が蓄積記録さ
れる。シート5を透過したX線3bはさらに低エネルギ
ーのX線をカットするフィルタ6を透過し、該フィルタ
6を透過した高エネルギーX線3cが第二の蓄積性蛍光
体シート7に照射される。これにより該シート7にも被
写体4のX線画像が蓄積記録される。被写体4には、サ
ブトラクション処理を行なうにあたって2つのX線画像
の位置合わせを行なうための基準となる2つのマーク8
が付されている。
【0039】なお、上記X線撮影装置は一回の撮影で2
枚のシート5,7にX線画像を蓄積記録するものである
が、時間的に相前後した2つのタイミングでそれぞれ1
枚ずつ撮影を行なってもよい。
【0040】図11は、X線画像読取装置と粒状改善エ
ネルギーサブトラクションを実施するための画像処理表
示装置の斜視図である。図10に示すX線撮影装置1で
撮影が行なわれた後、第一および第二の蓄積性蛍光体シ
ート5,7が一枚ずつX線画像読取装置10の所定位置
にセットされる。ここでは、第一の蓄積性蛍光体シート
5に蓄積記録された第一のX線画像の読取りの場合につ
いて説明する。
【0041】所定位置にセットされた、第一のX線画像
が蓄積記録された蓄積性蛍光体シート5は、図示しない
駆動手段により駆動されるエンドレスベルト等のシート
搬送手段15により、矢印Y方向に搬送(副走査)され
る。一方、レーザ光源16から発せられた光ビーム17
はモータ18により駆動され矢印Z方向に高速回転する
回転多面鏡19によって反射偏向され、fθレンズ等の
集束レンズ20を通過した後、ミラー21により光路を
かえて蓄積性蛍光体シート14に入射し、副走査の方向
(矢印Y方向)と略直角な矢印X方向に主走査する。蓄
積性蛍光体シート14の、光ビーム17が照射された箇
所からは、蓄積記録されているX線画像情報に応じた光
量の輝尽発光光22が発せられ、この輝尽発光光22は
光ガイド23によって導かれ、フォトマルチプライヤ(
光電子増倍管)24によって光電的に検出される。光ガ
イド23はアクリル板等の導光性材料を成形して作られ
たものであり、直線状をなす入射端面23a が蓄積性
蛍光体シート14上の主走査線にそって延びるように配
され、円環状に形成された射出端面23bにフォトマル
チプライヤ24の受光面が結合されている。入射端面2
3a から光ガイド23内に入射した輝尽発光光22は
、該光ガイド23の内部を全反射を繰り返して進み、射
出端面23b から射出してフォトマルチプライヤ24
に受光され、放射線画像を表わす輝尽発光光22がフォ
トマルチプライヤ24によって電気信号に変換される。
【0042】フォトマルチプライヤ24から出力された
アナログ信号Sは、ログアンプ25で対数的に増幅され
た後、A/D変換器26に入力され、サンプリングされ
て、ディジタルの画像信号S0が得られる。この画像信
号S0は第一の蓄積性蛍光体シート5に蓄積記録された
第一のX線画像を表わすものであり、ここでは第一の画
像信号S01 と呼ぶ。この第一の画像信号S01 は
画像処理表示装置30内の内部メモリに一旦記憶される
。
【0043】この画像処理表示装置30には、種々の指
示を入力するキーボード31、指示のための補助情報や
画像信号に基づく可視画像を表示するCRTディスプレ
イ32、補助記憶媒体としてのフロッピィディスクが装
填され駆動されるフロッピィディスク駆動装置33、お
よびCPUや内部メモリが内蔵された本体部34が備え
られている。
【0044】次に上記と同様にして、第二の蓄積性蛍光
体シート7に蓄積記録された第二のX線画像を表わす第
二の画像信号S02 が得られ、この第二の画像信号S
02 も画像処理表示装置30内の内部メモリに一旦記
憶される。
【0045】図2は、画像処理表示装置内の内部メモリ
に記憶された第一および第二のX線画像を表わす2つの
画像信号S01 ,S02 に基づいて、該画像処理表
示装置内で行なわれる処理の流れの一例を表わした図で
ある。
【0046】画像処理表示装置内の内部メモリに記憶さ
れた、第一および第二のX線画像信号S01 ,S02
 は、図2に示すそれぞれ第一のX線画像41,第二の
X線画像42を担持する信号である。第一のX線画像4
1は比較的低エネルギーX線による画像であり、第二の
X線画像42は比較的高エネルギーX線による画像であ
るが、互いに軟部および骨部の濃度は異なるものの両者
ともこれら軟部および骨部の双方が記録された原画像で
ある。
【0047】これら第一および第二のX線画像信号S0
1 ,S02 は図11に示す画像処理表示装置30内
の内部メモリから読み出され、先ずこれら2つの画像信
号S01 ,S02 がそれぞれ担持する各X線画像4
1,42の相対的な位置合わせが画像信号上で行なわれ
る(特開昭58−163338 号公報参照)。この位
置合わせは、図2に示す2つのマーク8が重なるように
2つのX線画像を相対的に直線的な移動および回転移動
を行なうことにより行なわれる。
【0048】この後、サブトラクション処理が行なわれ
る。
【0049】ここでX線の吸収係数μを、被写体の軟部
と骨部、および低エネルギーX線と高エネルギーX線と
に分けて次のように定める。
【0050】μL T :低エネルギーX線による軟部
の吸収係数 μH T :高エネルギーX線による軟部の吸収係数μ
L B :低エネルギーX線による骨部の吸収係数μH
 B :高エネルギーX線による骨部の吸収係数このと
き、2つの画像信号S01 ,S02 の互いに対応す
る画素毎に、式
【0051】
【数1】
【0052】(ただしCはバイアス成分を表わす)◆に
従って重み付け引き算を行なうことによって、骨部の陰
影が抽出された骨部画像43(図2参照)を表わす骨部
画像信号S1が求められる。
【0053】また、式
【0054】
【数2】
【0055】(ただしC′はバイアス成分を表わす)◆
に従って重み付け引き算を行なうことにより軟部画像を
表わす軟部画像信号S2を求めることができるが、本例
ではこの演算は不必要である。
【0056】さらに、式
【0057】
【数3】
【0058】に従って互いに対応する各画素毎に加算処
理を行なうことにより2つのX線画像41,42の重ね
合わせ画像44が生成される。この重ね合わせ画像44
も軟部および骨部の双方が記録された原画像である。こ
の重ね合わせ画像44に代えてX線画像41もしくはX
線画像42を用いることも可能であるが、重ね合わせ画
像44は2つのX線画像41,42を重ね合わたもので
あるためこれら各X線画像のいずれと比べてもノイズ成
分が低減されており、したがってその後の処理に有利と
なる。
【0059】次に骨部画像信号S1を処理することによ
り、骨部画像43に含まれるノイズ成分の抽出が行なわ
れる。
【0060】図3は、骨部画像および骨部画像信号を処
理して求めた画像の、空間周波数fに対するスペクトル
を表わした図である。図3に示すグラフ51が骨部画像
43のスペクトルを表わしており、ノイズ成分53が含
まれている。
【0061】ここで、先ず、骨部画像信号S1に平滑化
処理が行なわれる。この平滑化処理方法としては、例え
ば各画素に対し該画素を中心とした所定領域内の各画素
に対応する画像信号の平均を求め、この平均値を中心の
画素の画像信号とする単純な平均化処理方法、上記所定
領域内の画像信号の中央値(メジアン)を中心の画素の
画像信号とするというメジアンフィルタを用いる方法、
上記所定領域内をさらに複数の小領域に分け、各小領域
毎に分散を求めて分散の最も小さい小領域の平均値を中
央の画素の画像信号の値とするエッジ保存フィルタ(V
−フィルタ)を用いる方法、画像信号をフーリエ変換し
、ノイズ成分に対応する高空間周波数成分を取り除いた
後逆フーリエ変換する方法等を用いることができるが、
上記ボケマスク処理方法はエッジがぼけるという欠点を
有し、また上記メジアンフィルタを用いる方法は画素を
入れ替えることになるので等高線状のアーチファクトが
発生する場合があり、さらに上記エッジ保存フィルタを
用いるとハニカム状のアーチファクトが発する場合があ
り、フーリエ変換する方法は演算に時間がかかるという
問題がある。そこで本例では、上記いずれの方法とも異
なる以下に示すヒストグラム適応フィルタを用いた平滑
化が行なわれる。この方法を用いると画像情報として必
要なエッジ(互いに異なる2つの組織の陰影の境界を定
めるステップ状の濃度変化)を保存したままかつ上記ア
ーチファクトなしにノイズを除去することができ、また
簡単な演算で短時間にノイズを除去することができると
いう長所を有する。
【0062】まず骨部画像の各画素毎に該画素を中心と
した所定領域内の多数の画素の画像信号S1のヒストグ
ラムを作成する。
【0063】図4(A),図4(B)は、上記のように
して求めた、ある画素(画像信号S1′)を中心とした
所定領域内の多数の画素に対応する画像信号S1の出現
頻度をプロットした、互いに異なる2つのヒストグラム
を表わした図、図5は、画像信号S1と中央の画素の画
像信号S1′との差を変数とした関数の一例を表わした
図である。
【0064】図4(A),図4(B)に示すようなヒス
トグラムを表わす関数を一般にh(S1)で表わし、絶
対値|S1−S1′|が増加するに従って単調減少する
、例えば図5に示すような関数をf(S1−S1′)と
する。このとき、式
【0065】
【数4】
【0066】に従って処理後の頻度を表わす関数g(S
1)を求める。この関数g(S1)は、関数h(S1)
が図4(A)に示すように複数の山を有する場合は中央
の画素の画像信号S1′が属する山のみを抽出する作用
を有する。
【0067】上記(4) 式に従って関数g(S1)を
求めた後、該関数g(S1)で重みづけをした画像信号
S1の平均的な値S1*を求める。即ち、具体的には例
えば関数g(S1)の一次モーメントが次式に従って求
められる。
【0068】
【数5】
【0069】骨部画像の各画素をそれぞれ中心の画素と
して上記(4),(5) 式に従う処理が行なわれ、こ
れにより平滑化画像信号S1*(簡単のため、各画素に
対応する画像信号と画像全体を表わす画像信号とで同一
の記号を用いている。)が求められる。この平滑化画像
信号S1*は図2のグラフ52に示すように、主として
もとの骨部画像信号S1の高空間周波数成分を取り除い
た信号であるが、エッジ近傍の画素については図4(A
)に示すようにその画素の属する山のみを抽出した後の
平均的な値を求めた信号であるため、もとの骨部画像中
のエッジはぼけることなく保存されている。
【0070】次に各画素毎に重ね合わせ画像44を表わ
す重ね合わせ画像信号S0(上記(3)式参照)から平
滑化画像信号S1*を重みづけ引き算、即ち
【0071
】
【数6】
【0072】(ただしC″はバイアス分を表わす。)◆
を行なうことにより、画像情報としては上記(2) 式
で表わされる軟部画像と略同一の情報を担持するととも
に上記(2) 式で表わされる軟部画像よりもノイズ成
分が低減された処理済軟部画像46(図2参照)が求め
られる。
【0073】(6) 式に従って求められた画像信号S
2′は画像処理表示装置30のCRTディスプレイ32
に送られ、この画像信号S2′に基づく可視画像がCR
Tディスプレイ32に再生表示される。
【0074】なお、上記例は骨部画像信号S1を平滑化
して原画像から引くことにより軟部画像信号S2′を求
める例であるが、骨部画像を観察対象とする場合は、上
記(2) 式に基づいて軟部画像信号S2を求め、この
軟部画像信号S2を平滑化して原画像から引くことによ
りノイズ成分が低減された骨部画像を求めればよい。
【0075】次に図2を参照して説明した上記例と実質
同一の処理である、他の例について説明する。図6はこ
の実質同一の例を説明するために画像処理表示装置内の
内部メモリに記憶された第一および第二のX線画像を表
わす2つの画像信号S01 ,S02 に基づいて、画
像処理表示装置内で行なわれる処理の流れの他の例を示
した図である。図2と同一の要素には図2に付した番号
,記号を付し、図2を用いて説明した箇所については重
複説明を避けるために、説明は省略する。
【0076】2つのX線画像41,42から上記(1)
 式,(2)式に基づいて骨部画像43(骨部画像信号
S1)と軟部画像47(軟部画像信号S2)が求められ
る。
【0077】次に前述した例と同様にして骨部画像信号
S1を上記(4) 式,(5)式に基づいて処理するこ
とにより、骨部画像43に含まれるノイズ成分が低減化
された平滑化画像信号S1*が求められ、その後各画素
毎に骨部画像信号S1から平滑化画像信号S1*を引き
算することにより、ノイズ成分のみが抽出されたノイズ
画像48(ノイズ信号SN )が求められる。
【0078】
【数7】
【0079】このノイズ信号SN は図3のグラフ53
に示すように骨部画像のノイズ成分を抽出した信号であ
る。 ここで平滑化画像信号S1*は骨部画像のエッジの情報
はたとえノイズ成分と同程度の高空間周波数であっても
保存されているため、上記(7)式に従って骨部画像信
号S1と平滑化画像信号S1*との差を求めることによ
りエッジの情報はきれいにキャンセルされ、したがって
エッジの情報を失わせるような平滑化処理を行なった場
合と比べ、ノイズ信号SN はより純粋に骨部画像のノ
イズ成分のみを担持した信号となる。
【0080】次にこのようにして求められたノイズ信号
SN と軟部画像47(図6参照)を表わす軟部画像信
号S2とが各画素毎に重み付け足し算され、これにより
画像情報とては上記軟部画像47と略同一の情報を担持
するとともに該軟部画像47よりもノイズ成分が低減さ
れた処理済軟部画像46が求められる。本例ではこの重
み付け足し算は、式
【0081】
【数8】
【0082】に従って行なわれ、これによりノイズ成分
の一層の低減が図られる。
【0083】なお、詳細な説明は省略するが、図2を参
照して説明した例と図6を参照して説明した例とは実質
的に同一のものである。(詳細については特願平2−9
4653 号参照)図7は、粒状改善エネサブの他の例
の処理の流れを表わした図、図8は、図7に示す各画像
の所定の一方向についてのプロファイルを模式的に示し
た図である。図7において、図2もしくは図6と対応す
る要素にはこれら図2,図6と同一の番号,記号を付し
重複説明は省略する。
【0084】図8(a),(b) は、それぞれ、X線
画像(原画像)41,42を模式的に表わした図であり
、X線画像41,42上の所定の一方向(x方向)に沿
った画像信号S01 ,S02 の値をプロットしたも
のであり、これらの画像信号S01 ,S02 には互
いにその値は異なるものの一様な軟部(図に斜線を施し
た部分)を表わす信号成分とステップ状に変化した骨部
を表わす信号成分とが重畳され、かつランダムなノイズ
成分が重畳されている。
【0085】2つのX線画像(原画像)41,42を表
わすこの2つの画像信号S01 ,S02 に基づいて
上記(2) 式に基づいて重みづけ減算処理(記号−で
表わす)を行なうことにより軟部画像47を表わす軟部
画像信号S2が求められ、また2つの画像信号S01 
,S02 に基づいて上記(3) 式に基づいて加算処
理(記号+で表わす)を行なうことにより、重ね合わせ
画像44を表わす重ね合わせ画像信号S0が求められる
。
【0086】図8(c) は、重ね合わせ画像信号S0
を模式的に表わした図であり、図8(a),(b) と
同様に軟部を表わす一様な信号成分(図に斜線を施した
部分)とステップ的に変化した骨部を表わす信号成分と
、さらにランダムなノイズ成分とが重畳されているが、
このノイズ成分は図8(a),(b) に示す2つのX
線画像41,42と比べ低減化されている。また図8(
d) は、上記(2) 式に基づいて求められた軟部画
像信号S2を表わした図である。一様な軟部を表わす信
号成分のみが抽出されているが、ランダムなノイズ成分
は上記2つのX線画像41,42(図8(a),(b)
 )のいずれよりも増加している。
【0087】また本例では求める必要はないが、仮に上
記(1) 式に基づいて骨部画像信号S1を求めたとし
た場合の該骨部画像信号S1を表わした図である。ステ
ップ状に変化した骨部を表わす信号成分が抽出されてい
るが、軟部画像信号S2(図8(d) )と同様に、ラ
ンダムなノイズ成分は上記2つのX線画像41,42(
図8(a),(b) )のいずれよりも増加している。
【0088】ここで軟部画像47(軟部画像信号S2,
図8(d) )に平滑化処理51(図7参照)が施され
、平滑化軟部画像61を表わす平滑化軟部画像信号S2
*(図8(f) )が求められる。この平滑化処理51
では、軟部画像47の例えば1.0 サイクル/mm以
上の高空間周波数成分がカットされる。
【0089】次に重ね合わせ画像信号S0から平滑化軟
部画像信号S2*が重みづけ減算され、これにより骨部
画像62を表わす骨部画像信号S1′が求められる。こ
の骨部画像信号S1′は図8(g) に示されるように
、骨部画像信号S1(図8(e) )と比べランダムな
ノイズ成分が低減化されているが、軟部画像47を平滑
化処理した影響が表われ、軟部画像の高空間周波数成分
が若干混入している。
【0090】次に上記のようにして求められた骨部画像
信号S1′に平滑化処理52が施される。ここで施され
る平滑化処理52では、骨部画像62の例えば0.5 
サイクル/mm以上の空間周波数帯にある低コントラス
トの陰影(骨部画像信号S1′の変化の小さいもの)の
みがカットされる。この処理方法としては、例えば所定
の画素P0 に対して0.5 サイクル/mmに対応す
る面積のウィンドウを考え、このウィンドウ内の各画素
にそれぞれ対応する各信号S1′のうち、所定の画素P
0 に対応する信号S10 ′の値±所定値内にある信
号S1′の平均値を求めて該平均値を所定の画素P0 
の新たな信号S10 ′とするフィルタを用いて骨部画
像62上を走査する方法等が採用される。この平滑化処
理52により、平滑化骨部画像63を表わす平滑化骨部
画像信号S1′*が求められる。この平滑化骨部画像信
号S1′*は、図8(i) に示すようにノイズ成分お
よび混入した軟部画像の高周波成分は低減されているも
のの立ち上がり部分も鈍ってしまっている。
【0091】次に、重ね合わせ画像信号S0から平滑化
骨部画像信号S1′*が重みづけ引き算され、軟部画像
64を表わす軟部画像信号S2′が求められる。この軟
部画像64は図8(h) に示すように、軟部画像47
(図8(d) )よりもノイズ成分は低減されているが
、平滑化骨部画像信号S1′*(図8(i) )の立上
り部分が鈍っている分だけ、その部分の骨部画像の情報
がノイズとして重畳されている。ただし、ランダムなノ
イズ部分、およびノイズとしての骨部画像の情報はかな
り小さく、したがってこの段階で一連の処理を停止し、
軟部画像信号S2′を画像処理表示装置30のCRTデ
ィスプレイ32(図11参照)に送ってこの軟部画像信
号S2′に基づく可視画像をCRTディスプレイに再生
表示し、観察するようにしてもよい。ただし本例では、
上記と同様な処理をさらに繰り返し、さらに画質の改善
が図られている。
【0092】軟部画像64を表わす軟部画像信号S2′
を求めた後、該軟部画像信号S2′に平滑化処理53が
施され、平滑化軟部画像65を表わす平滑化軟部画像信
号S2′*(図8(j)  )が求められる。この平滑
化処理53としては例えば1.5 サイクル/mm以上
の空間周波数成分をカットする処理が施される。
【0093】この平滑化軟部画像信号S2′*は重ね合
わせ画像信号S0から重みづけ減算処理され、骨部画像
66を表わす骨部画像信号S1″が求められる。この骨
部画像66は、図8(k) に示すように、骨部画像6
2(図8(g) )と比べランダムノイズおよびノイズ
として混入する軟部画像の情報も減少している。骨部画
像を観察対象とする場合はこの骨部画像信号S1″に基
づく可視画像をCRTディスプレイ32上に再生表示し
てもよい。
【0094】本例では、上記のように求められた骨部画
像信号S1″にさらに平滑化処理54が施されて平滑化
骨部画像67を表わす平滑化骨部画像信号S1″*(図
8(m) )が求められる。この平滑化処理54として
は例えば1.0 サイクル/mm以上の低コントラスト
成分のカットが行なわれる。
【0095】次に、重ね合わせ画像信号S0からこの平
滑化骨部画像信号S1″*が重みづけ引き算され、軟部
画像信号S2″が求められる。この軟部画像信号S2″
は図8(l) に示すように、前回求めた軟部画像信号
S2′(図8(h) )と比べランダムノイズおよびノ
イズとしての骨部画像の情報の双方ともさらに低減化さ
れた信号となる。
【0096】このようにして平滑化処理と重ね合わせ画
像(原画像)との重みづけ引き算を繰り返すことにより
、ノイズが順次低減化された骨部画像と軟部画像とを交
互に得ることができる。
【0097】図9は、図7を参照して説明した例と実質
同一の他の処理の流れを表わした図である。図7等と同
一の要素には該図7等と同一の番号,記号を付し説明は
省略する。図9に示す処理は、図7に示す骨部画像62
を求めるまでの処理(図2を参照して説明した処理(但
し図2とは骨部画像と軟部画像が入れ替っている))を
図2を参照して説明した処理(但し図2とは骨部画像と
軟部画像とが入れ替っている)に置き換えたものであっ
て、前述したようにこれらは互いに実質同一の処理であ
る。
【0098】なお、図9に示した処理では最初の段階の
み、図2を参照して説明した処理方法に置き換えたが、
この置き換えは繰り返し行なわれる処理の任意の段階に
ついて行なうことができいずれも実質同一の処理であり
、本発明にはこれら任意の1つもしくは複数の段階につ
いて変更された全ての実質同一の処理態様が包含される
ものである。
【0099】以上の各例はいずれも人体の胸部のX線画
像に基づいて軟部画像もしくは骨部画像を求める例であ
るが、また、本発明は軟部画像もしくは骨部画像を求め
るものに限られるものではなく、例えば、乳腺が強調さ
れた画像もしくは悪性腫瘍が強調された画像でもよく、
一般に被写体中の互いに異なる2つの組織がそれぞれ強
調もしくは抽出された2つの画像のうちの一方もしくは
双方を求める際に広く適用することができるものである
。
【0100】さらに、上記各例は、蓄積性蛍光体シート
を用いた例であるが、本発明は蓄積性蛍光体シートを用
いたものに限られるものではなくX線フイルム(撮影に
際して一般に増感スクリーンと組合わされる)等を用い
たものにも適用することができる。
【図面の簡単な説明】
【図1】本発明の基本的プロセスを示すブロック図
【図
2】画像処理表示装置内で行なわれる処理の流れの一例
を表わした図、
【図3】骨部画像および骨部画像信号を処理して求めた
画像の、空間周波数スペクトルを表わした図
【図4】あ
る画素を中心とした所定領域内の多数の画素に対応する
画像信号の出現頻度をプロットした、異なる2つのヒス
トグラムを表わした図
【図5】画像信号S1と所定領域の中心の画素の画像信
号S1′との差を変数とした関数の一例を表わした図
【
図6】画像処理表示装置内で行なわれる、図2に示した
処理と実質同一の他の処理の流れを表わした図
【図7】
本発明に適用される粒状改善エネサブの他の例の処理の
流れを表わした図
【図8】図7に示す各画像の所定方向についてのプロフ
ァイルを模式的に表わした図
【図9】図7に示した処理と実質同一の他の処理の流れ
を表わした図
【図10】X線撮影装置の概略図
【図11】X線画像読取装置と、本発明の前提となる粒
状改善エネルギーサブトラクション画像生成方法を実施
する画像処理表示装置の斜視図
【符号の説明】
1    X線撮影装置 2    X線管 3,3a,3b,3c    X線 4    被写体 5    第一の蓄積性蛍光体シート 5A    高圧画像 6    フィルタ(蓄積性蛍光体シートシート)6A
    中圧画像 7    第二の蓄積性蛍光体シート 7A    低圧画像 8    マーク 11    粒状改善エネルギーサブトラクション画像
(データ) 12    周波数フィルタ(高域低減フィルタ)13
    周波数フィルタ(高域透過フィルタ)16  
  レーザ光源 19    回転多面鏡 22    輝尽発光光 23    光ガイド 24    フォトマルチプライヤ 25    ログアンプ 26    A/D変換器 30    画像処理表示装置 41,42    X線画像(原画像)43,62,6
6    骨部画像 44    重ね合わせ画像(原画像)45,63,6
7    平滑化骨部画像46,47,64    軟
部画像 48    ノイズ画像 61,65    平滑化軟部画像

Claims (3)

    【特許請求の範囲】
  1. 【請求項1】  放射線エネルギーを分離させるために
    上下2枚の記録シートの間に低エネルギーの放射線をカ
    ットするフィルタを介在させて放射線を照射し、高エネ
    ルギーの放射線のみが下の記録シートに達するようにし
    、該両記録シートに記録された放射線画像を読み取って
    得られた画像データの間で減算処理をしてサブトラクシ
    ョン画像を得るエネルギーサブトラクション方法であっ
    て、高周波成分については減算処理をせず、低周波成分
    についてのみ減算処理をして粒状を改善するエネルギー
    サブトラクション方法において、蓄積性蛍光体シートを
    前記低エネルギー放射線をカットするフィルタとして使
    用し、この蓄積性蛍光体シートに蓄積された画像データ
    の高周波成分を前記エネルギーサブトラクション画像に
    加算して粒状を改善したことを特徴とするエネルギーサ
    ブトラクション画像生成方法。
  2. 【請求項2】  前記蓄積性蛍光体シートに蓄積された
    画像データの高周波成分をエネルギーサブトラクション
    画像に加算する際、前記エネルギーサブトラクション画
    像のデータには高周波を低減する周波数フィルタをかけ
    、前記蓄積性蛍光体シートに蓄積された画像データには
    高周波透過フィルタをかけて、両フィルタを通した画像
    データを加算することを特徴とする請求項1記載のエネ
    ルギーサブトラクション画像生成方法。
  3. 【請求項3】  前記両フィルタのレスポンスの和が全
    周波数領域に亘って1であることを特徴とするエネルギ
    ーサブトラクション画像生成方法。
JP3086881A 1991-04-18 1991-04-18 エネルギーサブトラクション画像生成方法 Expired - Lifetime JP2640582B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP3086881A JP2640582B2 (ja) 1991-04-18 1991-04-18 エネルギーサブトラクション画像生成方法
US07/866,469 US5301107A (en) 1991-04-18 1992-04-10 Method for forming energy subtraction images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3086881A JP2640582B2 (ja) 1991-04-18 1991-04-18 エネルギーサブトラクション画像生成方法

Publications (2)

Publication Number Publication Date
JPH04318775A true JPH04318775A (ja) 1992-11-10
JP2640582B2 JP2640582B2 (ja) 1997-08-13

Family

ID=13899180

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3086881A Expired - Lifetime JP2640582B2 (ja) 1991-04-18 1991-04-18 エネルギーサブトラクション画像生成方法

Country Status (2)

Country Link
US (1) US5301107A (ja)
JP (1) JP2640582B2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002197440A (ja) * 2000-12-27 2002-07-12 Fuji Photo Film Co Ltd 放射線画像データ処理方法および装置
JP2009061253A (ja) * 2007-08-15 2009-03-26 Fujifilm Corp 画像成分分離装置、方法、およびプログラム
JP2018038646A (ja) * 2016-09-08 2018-03-15 富士フイルム株式会社 画像処理装置、方法およびプログラム
JP2018134205A (ja) * 2017-02-21 2018-08-30 富士フイルム株式会社 エネルギーサブトラクション処理装置、方法およびプログラム
JP2019520885A (ja) * 2016-06-07 2019-07-25 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. スペクトル及び非スペクトル物質分解を組み合わせることによる定量的撮像の精度及び分解能の向上
WO2020012766A1 (ja) * 2018-07-09 2020-01-16 キヤノン株式会社 画像処理装置およびその制御方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3467285B2 (ja) * 1993-04-02 2003-11-17 コニカミノルタホールディングス株式会社 放射線画像処理方法
US6125166A (en) * 1998-01-13 2000-09-26 Fuji Photo Film Co., Ltd. Method of forming energy subtraction images
DE10109586A1 (de) * 2001-02-28 2002-09-05 Philips Corp Intellectual Pty Verfahren und Vorrichtung zur Bildverarbeitung von Röntgenaufnahmen
US6917697B2 (en) * 2001-05-08 2005-07-12 Ge Medical Systems Global Technology Company, Llc Method and apparatus to automatically determine tissue cancellation parameters in X-ray dual energy imaging
JP4008328B2 (ja) * 2002-09-30 2007-11-14 富士フイルム株式会社 放射線画像情報読取方法および装置
JP5106978B2 (ja) * 2007-10-15 2012-12-26 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー X線ct装置
GB201813372D0 (en) * 2018-08-16 2018-10-03 Rolls Royce Plc Machining error compensation

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5512429A (en) * 1978-07-12 1980-01-29 Fuji Photo Film Co Ltd Radioactive image reader
US4315318A (en) * 1978-12-26 1982-02-09 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiation image
JPS55116340A (en) * 1979-02-28 1980-09-06 Fuji Photo Film Co Ltd Method and device for processing gradation of radiation picture
JPS5611395A (en) * 1979-07-11 1981-02-04 Fuji Photo Film Co Ltd Radiation image writeereading device
JPS56104645A (en) * 1979-12-25 1981-08-20 Fuji Photo Film Co Ltd Radiation picture treating method and its device
JPS58163338A (ja) * 1982-03-20 1983-09-28 富士写真フイルム株式会社 放射線画像のサプトラクシヨン処理方法
JPS5983486A (ja) * 1982-11-04 1984-05-14 Fuji Photo Film Co Ltd 放射線画像のエネルギ−・サブトラクシヨン方法およびその方法に用いられる蓄積性螢光体シ−ト、蓄積性螢光体シ−ト積層体並びに蓄積性螢光体シ−トフイルタ積層体
JPS60225541A (ja) * 1984-04-24 1985-11-09 富士写真フイルム株式会社 エネルギ−サブトラクシヨン用高速撮影装置
JPS615193A (ja) * 1984-06-16 1986-01-10 星野 謙三 破砕法
US5049746A (en) * 1989-10-19 1991-09-17 Fuji Photo Film Co., Ltd. Method and apparatus for displaying energy subtraction images

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002197440A (ja) * 2000-12-27 2002-07-12 Fuji Photo Film Co Ltd 放射線画像データ処理方法および装置
JP2009061253A (ja) * 2007-08-15 2009-03-26 Fujifilm Corp 画像成分分離装置、方法、およびプログラム
JP2019520885A (ja) * 2016-06-07 2019-07-25 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. スペクトル及び非スペクトル物質分解を組み合わせることによる定量的撮像の精度及び分解能の向上
JP2018038646A (ja) * 2016-09-08 2018-03-15 富士フイルム株式会社 画像処理装置、方法およびプログラム
JP2018134205A (ja) * 2017-02-21 2018-08-30 富士フイルム株式会社 エネルギーサブトラクション処理装置、方法およびプログラム
WO2020012766A1 (ja) * 2018-07-09 2020-01-16 キヤノン株式会社 画像処理装置およびその制御方法
JP2020005918A (ja) * 2018-07-09 2020-01-16 キヤノン株式会社 画像処理装置およびその制御方法
US11756176B2 (en) 2018-07-09 2023-09-12 Canon Kabushiki Kaisha Image processing apparatus and control method for same

Also Published As

Publication number Publication date
US5301107A (en) 1994-04-05
JP2640582B2 (ja) 1997-08-13

Similar Documents

Publication Publication Date Title
EP0031952B1 (en) Method of and apparatus for processing a radiation image
JP3002937B2 (ja) エネルギーサブトラクション画像処理方法
US5402338A (en) Method for forming energy subtraction images
EP0689794B1 (en) Method for forming energy subtraction images
JP2640582B2 (ja) エネルギーサブトラクション画像生成方法
JP3901770B2 (ja) エネルギーサブトラクション処理方法および装置
JP2627097B2 (ja) エネルギーサブトラクション画像生成方法および装置
US6125166A (en) Method of forming energy subtraction images
US5291403A (en) Method of and apparatus for processing radiation image
JP2729872B2 (ja) エネルギーサブトラクション画像生成方法
JP2981693B2 (ja) エネルギーサブトラクション画像生成方法
JPH0312785A (ja) 画像処理方法
JP2651625B2 (ja) エネルギーサブトラクション画像生成方法
JP2000060834A (ja) エネルギ―サブトラクション画像生成方法
JP2981694B2 (ja) エネルギーサブトラクション画像生成方法
JPH03222076A (ja) エッジ判定方法および装置
JPH11205682A (ja) エネルギーサブトラクション画像生成方法
JPH02193270A (ja) 放射線画像信号生成方法
JP2739372B2 (ja) 画像平滑化処理方法および装置
JPH08294054A (ja) エネルギーサブトラクション処理方法および装置
JP2867096B2 (ja) 画像平滑化処理方法
JPH0738758A (ja) 放射線画像のダイナミックレンジ圧縮方法
JP2002359781A (ja) エネルギーサブトラクション処理方法および装置
JP2631032B2 (ja) 放射線画像のエネルギーサブトラクション方法および装置
JPH02272681A (ja) 放射線画像処理装置

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 19970311

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080502

Year of fee payment: 11

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090502

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090502

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100502

Year of fee payment: 13

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110502

Year of fee payment: 14

EXPY Cancellation because of completion of term