JPH04154385A - 放射線画像のエネルギーサブトラクション方法および装置 - Google Patents

Abstract

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

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明は、放射線画像のエネルギーサブトラクション方
法および装置に関するものである。

（従来の技術） 記録された放射線画像を読み取って画像信号を得、この
画像信号に適切な画像処理を施した後、画像を再生記録
することが種々の分野で行われている。たとえば、後の
画像処理に適合するように設計されたガンマ値の低いＸ
線フィルムを用いてＸ線画像を記録し、このＸ線画像が
記録されたフィルムからＸ線画像を読み取って電気信号
に変換し、この電気信号（画像信号）に画像処理を施し
た後コピー写真等に可視像として再生することにより、
コントラスト　シャープネス、粒状性等の画質性能の良
好な再生画像を得ることの出来るシステムが開発されて
いる（特公昭８１−５１９３号公報参照）。

また本出願人により、放射線（Ｘ線、α線、β線、γ線
、電子線、紫外線等）を照射するとこの放射線エネルギ
ーの一部が蓄積され、その後可視光等の励起光を照射す
ると蓄積されたエネルギーに応じた光量の輝尽発光光を
放射する蓄積性蛍光体（輝尽性蛍光体）を利用して、人
体等の被写体の放射線画像を一部シート状の蓄積性蛍光
体に撮影記録し、蓄積性蛍光体シートをレーザ光等の励
起光で走査して輝尽発光光を生ぜしめ、得られた輝尽発
光光を光電的に読み取って画像信号を得、この画像信号
に基づいて被写体の放射線画像を写真感光材料等の記録
材料、ＣＲＴ等に可視像として出力させる放射線記録再
生システムがすでに提案されている（特開昭５５−１２
４２９号、同５Ｂ−１１３９５号。

同５５−１１３４０号、同５Ｂ−１６４６４５号、同５
５−１１６３４０号等）。

このシステムは、従来の銀塩写真を用いる放射線写真シ
ステムと比較して極めて広い放射線露出域にわたって画
像を記録し得るという実用的な利点を有している。すな
わち、放射線露光量に対する、蓄積後に励起によって発
光する輝尽発光光の光量が極めて広い範囲に渡って比例
することが認められており、従って種々の撮影条件によ
り放射線露光量がかなり大幅に変動しても、蓄積性蛍光
体シートより放射される輝尽発光光を読取りゲインを適
当な値に設定して光電変換手段により読み取って電気信
号（画像信号）に変換し、この画像信号を用いて写真感
光材料、ＣＲＴ等の表示装置に放射線画像を可視像とし
て出力することによって、放射線露光量の変動に影響さ
れない放射線画像を得ることができる。

上記のようにＸ線フィルムや蓄積性蛍光体シート等の記
録シートを用いるシステムにおいて、記録シートに記録
された複数の放射線画像を読み取って複数の画像信号を
得た後、これらの画像信号に基づいて上記放射線画像の
サブトラクション処理を施すことがある。

ここで、放射線画像のサブトラクション処理とは、互い
に異なった条件で撮影された複数の放射線画像の差に対
応する画像を得る処理をいい、具体的にはこれら複数の
放射線画像を所定のサンプリング間隔で読み取って各放
射線画像に対応する複数のディジタルの画像信号を得、
これら複数のディジタルの画像信号の各対応するサンプ
リング点毎に減算処理を施すことにより、放射線画像中
の特定の被写体部分のみを強調または抽出した放射線画
像を得る処理をいう。

このサブトラクション処理には基本的には次の二つがあ
る。すなわち、造影剤の注入により被写体の特定の部分
（たとえば人体を被写体としたときの血管等）が強調さ
れた放射線画像から造影剤が注入されていない放射線画
像を引き算（サブトラクト）することによって被写体の
特定の部分（たとえば血管等）を抽出するいわゆる時間
差サブトラクションと、被写体の特定の部分が互いに異
なるエネルギーを有する放射線に対して異なる放射線吸
収率を有することを利用して、同一の被写体に対して互
いに異なるエネルギーを有する放射線を照射してこれら
互いに異なるエネルギーを有する各放射線による複数の
放射線画像を得、これら複数の放射線画像を適当に重み
付けしてその差を演算することによって被写体の特定部
分を抽出するいわゆるエネルギーサブトラクションとが
ある。本出願人も蓄積性蛍光体シートを用いたエネルギ
ーサブトラクションについて提案している（特開昭５９
−８３４８８号公報、特開昭６０−２２５５４１号公報
参照）。

（発明が解決しようとする課題） 上記特開昭６０−２２５５４１号公報に記載された方法
は、互いにエネルギーの異なる放射線を用いた２回の放
射線撮影を行ない、その結果得られた２枚の放射線画像
を読み取って２つのディジタル画像信号を得、これらの
画像信号に基づいてサブトラクションを行なうものであ
る。しかしこの方法では２回の撮影の間に時間差がある
ため、その間に被写体が動いてしまい、サブトラクショ
ン処理後の画像信号に基づいて再生した可視画像に、こ
の動きに起因する、複数の画像の不一致による偽画像（
モーションアーティファクト）が発生し、この可視画像
の画質が著しく低下してしまうという問題がある。

また、上記特開昭５９−８３４８８号公報には、放射線
エネルギーにより吸収率の異なるフィルタを挾んだ２枚
の記録シートに被写体を透過した放射線を照射すること
等により１回の撮影でエネルギーサブトラクションを行
なうことのできる方法が提案されている。この方法を用
いると被写体の動きによる偽画像は生じない。

しかしこの方法を用いると、放射線画像が記録される２
枚の記録シートはフィルタを挾んで空間的に互いに異な
る位置に配置されることとなり、したがってこれら２枚
の記録シートの放射線源および被写体からの距離が異な
り、このためこれら２枚の記録シートに記録される放射
線画像の寸法が異なり、これが偽画像の原因となる。

この偽画像を防止するため、被写体とともに２つのマー
クを各記録シートに写し込み、これら各記録シートに記
録された放射線画像を読取って画像信号を得た後、この
画像信号上で各放射線画像に記録された２つのマークの
間の距離を互いに同一とするように放射線画像の一方も
しくは双方の寸法を補正し、これにより互いに同一寸法
の放射線画像を生成した後サブトラクション処理を行な
うことが提案されている（特開平２−１０４１７４号公
報。

特願平１−２４８９４５号参照）。

しかし記録シートに記録された放射線画像を読み取って
画像信号を得る際は、該放射線画像を例えば０．１ｍｍ
間隔等、所定の間隔毎に離散的にサンプリングしながら
読み取るのが通常であり、このためこのような離散的な
読み取りにより得られた画像信号に基づいて各放射線画
像に記録された２つのマークの間の距離を求めた場合、
その距離は上記サンプリングに伴う誤差を含んだものと
なる。

例えば、２枚の放射線画像の内の第一の放射線画像にお
けるマーク間の実際の距離ｄがｄ　−２００，０１１で
あり、上記２枚の放射線画像のうちの第二の放射線画像
のうちの第二の放射線画像におけるマーク間の距離ｄ′
が実際の寸法ではｄ’　−２００，２−一であり、これ
ら２つの放射線画像は互いに寸法が異なるにも拘らず、
これを０．１ｍｓ間隔でサンプリングした画像信号に基
づいて各マークの位置を求めた場合、各マーク毎に最大
−つのサンプリング間隔の半分、即ち０．０５ｉ■の誤
差をもって求められ、したがって上記第一の放射線画像
についてはｄ　−２００，１諺ｌと求められる場合があ
り、一方第二の放射線画像についてもｄ’　＝２００．
１ｍｍと求められる場合があり、この場合寸法の差異は
ないと判定されてしまうこととなり、寸法が互いに異な
ることにより生じる偽画像が改善されないまま残ってし
まうこととなる。

これを解決する方法の一つとしてサンプリング間隔を小
さくすることも考えられるが、この解決方法を採用する
と放射線画像−枚あたりの読み取り速度が遅くなり、も
しくは放射線画像１枚あたり同じ速度で読み取るために
はより高速の読み取り回路を必要とし、また読み取って
得た画像信号を蓄えておくためにより大容量のメモリを
必要とすることとなる。

本発明は、上記事情に鑑み、サンプリング間隔を小さく
することなく、複数の放射線画像の寸法の差をより高精
度に求めることのできる放射線画像のエネルギーサブト
ラクシジン方法、および該方法を採用した装置を提供す
ることを目的とするものである。

（課題を解決するための手段） 上記目的を達成するための本発明の放射線画像のエネル
ギーサブトラクション方法は、互いに異なるエネルギー
分布を有する放射線による互いに同一の被写体の複数の
放射線画像を記録シートに記録し、該記録シートから前
記複数の放射線画像を読み取って該複数の放射線画像の
それぞれを表わす複数の画像信号を得、該複数の画像信
号に基づいて該複数の画像信号がそれぞれ担持する前記
放射線画像の寸法が互いに同一となるように該寸法を補
正し、該補正後の画像信号に基づいてサブトラクション
処理を行なう放射線画像のエネルギーサブトラクション
方法において、前記放射線画像の記録に当たり、前記複
数の放射線画像上の互いに対応する位置に前記寸法の補
正を行なうためのマークが３つ以上ずつ記録された前記
複数の放射線画像を記録し、 前記寸法の補正に当たり、前記複数の放射線画像のそれ
ぞれに３つ以上ずつ記録された前記マークの前記放射線
画像上の位置関係に基づいて、前記寸法の補正を行なう
ことを特徴とするものである。

また本発明の放射線画像のエネルギーサブトラクション
装置は、 放射線を発する放射線源と、被写体が配置される被写体
配置部と、互いに異なるエネルギー分布を有する放射線
による複数の放射線画像を記録するための記録シートが
配置される、前記被写体配置部を挟んで前記放射線源と
対向する位置に設けられたシート保持部とからなる記録
ユニット、前記記録ユニットにおいて前記被写体の放射
線画像の記録の行なわれた前記記録シートから、前記被
写体の複数の放射線画像のそれぞれを表わす複数の画像
信号を得る読取りユニット、前記読取りユニットで得ら
れた前記複数の画像信号に基づいて該複数の画像信号が
それぞれ担持する前記放射線画像の寸法が互いに同一と
なるように該寸法を補正する寸法補正ユニット、および
前記寸法補正ユニットで寸法の補正の行なわれた後の放
射線画像を担持する画像信号に基づいてサブトラクショ
ン処理を行なう演算ユニットを備えた放射線画像のエネ
ルギーサブトラクション装置において、前記記録ユニッ
トが、前記複数の放射線画像上の互いに対応する位置に
前記寸法の補正を行なうためのマークが３つ以上ずつ記
録された前記複数の放射線画像を記録するものであり、 前記寸法補正ユニットが、前記複数の放射線画像のそれ
ぞれに３つ以上ずつ記録された前記マークの前記放射線
画像上の位置関係に基づいて、前記寸法の補正を行なう
ものであることを特徴とするものである。

（作　　用） 本発明の放射線画像のエネルギーサブトラクション方法
および装置では、各放射線画像に３つ以上のマークを記
録し、これら３つ以上のマークに基づいて寸法を補正す
るようにしたため、前述した特開平２−１０４１７４号
公報、特願平１−２４６９４５号に記載された、各放射
線画像につき２つのマークを記録してこれらのマークに
基づいて寸法を補正する方法と比べ、サンプリング間隔
を小さくすることなく、より正確に寸法の補正をするこ
とができ、したがって従来よりも確実に偽画像の発生を
防止することができる。

（実　施　例） 以下、図面を参照して、本発明の実施例について説明す
る。

第１図は、本発明のエネルギーサブトラクション装置の
記録ユニットの実施例であるＸ線撮影装置の概略図であ
る。

このＸ線撮影装置１のＸ線管２から発せられたＸ線３に
より被写体４が照射される。被写体４を透過したＸ線３
ａは第一の蓄積性蛍光体シート５に照射され、Ｘ線３ａ
のエネルギーの一部が該第−の蓄積性蛍光体シート５に
蓄積され、これにより該シート５に被写体４のＸ線画像
が蓄積記録される。

シート５を透過したＸ線３ｂはさらにフィルタ６を透過
し、該フィルタ６を透過したＸ線３Ｃが第二の蓄積性蛍
光体シート７に照射される。これにより該シート７にも
被写体４のＸ線画像が蓄積記録される。尚、被写体４に
は、各蓄積性蛍光体シート５．７に蓄積記録されたＸ線
画像の寸法の補正を行なうための４つのマーカ１１．１
２．１３．１４が付されている。

第２図は各蓄積性蛍光体シート５．７に蓄積記録された
Ｘ線画像を模式的に表わした図である。

各蓄積性蛍光体シート５，７のほぼ全面に被写体４の各
Ｘ線画像４ａ、　４ａ’が蓄積記録され、またこれらの
Ｘ線画像４ａ、　４ａ’　とともにマーカ１１．１２゜
１３．１４の像１１ａ、１２ａ、１３ａ、１４ａ　；　
ｌｌａ　’　、１２ａ’１３ａ　’　、１４ａ’　　（
以後この像をマークと称する。）が形成されている。こ
のマーク１１ａ、１２ａ、１３ａ、１４ａ　；Ｌｌａ　
’　、１２ａ’　、１３ａ’　、１４ａ’　は蓄積性蛍
光体シート５．７の互いに対応する位置に記録されてお
り、画像信号を得た後このマークｌｌａ、１２ａ、１３
ａ、１４ａ　；１１ａ　’　、１２ａ’　、１３ａ’　
、１４ａ’に基づいて２つのＸ線画像の相対的な位置合
わせが行なわれる。

ここで、第１図に示すように各蓄積性蛍光体シート５，
７はＸ線管２および被写体４からの距離が互いに異なっ
ているため、各蓄積性蛍光体シート５．７に蓄積記録さ
れたＸ線画像４ａ、　４ａ’の倍率およびマーク１１ａ
、１２ａ、１３ａ、１４ａ　；　ｌｌａ　’　、１２ａ
’、　１３ａ’　、　１４ａ’　の位置が互いに異なっ
ている。この倍率の相違が偽画像の原因の一つとなる。

第３図は、本発明のエネルギーサブトラクション装置の
読取ユニットの一実施例であるＸ線画像読取装置、およ
び寸法補正ユニットと演算ユニットの一実施例である画
像処理表示装置の斜視図である。

第１図に示すＸ線撮影装置１で撮影が行なわれた後、第
一および第二の蓄積性蛍光体シート５゜７が一枚ずつＸ
線画像読取装置１０の所定位置にセットされる。ここで
は、第一の蓄積性蛍光体シート５に蓄積記録された第一
のＸ線画像の読取りの場合について説明する。

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

フォトマルチプライヤ２４から出力されたアナログ信号
Ｓは、ログアンプ２５で対数的に増幅された後、Ａ／Ｄ
変換器２６に入力され、サンプリングされて、ディジタ
ルの画像信号ＳＯが得られる。この画像信号ＳＯは第一
の蓄積性蛍光体シートらに蓄積記録された第一のＸ線画
像を表わすものであり、第一の画像信号ＳＯ１と呼ぶ。

この第一の画像信号Ｓ０１は画像処理表示装置３０内の
内部メモリに一旦記憶される。

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

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

このようにしてサブトラクション演算を行なうべき２つ
の画像信号ＳＯ１，ＳＯｚが内部メモリに記憶されると
、これら２つの画像信号Ｓ０１゜Ｓ０２が読み出されて
、これら２つの画像信号ＳＯ□、ＳＯ７が担持する各Ｘ
線画像の寸法が同一ととなるように寸法の補正演算が行
なわれる。尚、この寸法の補正演算はＣＰＵ等の内蔵さ
れた本体部３４で実行されるが、この寸法の補正演算を
行なう機能を担うハードウェアとソフトウェアとの組合
せが、本発明の寸法補正ユニットの一例と観念される。

ここで、本実施例における、画像信号Ｓ０１゜ＳＯ□、
が表わす２つのＸ線画像の寸法補正演算、回転演算およ
び位置合わせ演算について説明する。

第一の画像信号Ｓ０１が担持する第一のＸ線画像上の各
サンプリング点の座標を（Ｘｌ、Ｙｓ　）、第二の画像
信号Ｓ０２が担持する第二のＸ線画像上の各サンプリン
グ点の座標を（ＸＺ、Ｙ２）とし、ａ、ｂ、ｃ、ｄを係
数としたとき、ヘルマート変換 ・・・・・・（１） に従って第一のＸ線画像の座標を変換することにより、
第一のＸ線画像と第二のＸ線画像とを重ね合わせること
がきる。ここで（１）式に基づく座標変換では、第一の
Ｘ線画像全体を均等に拡大もしくは縮小すること、該第
−のＸ線画像全体を回転すること、および該第−のＸ線
画像をＸ方向、Ｙ方向に平行移動すること、の全てが同
時に行なわれる。

次に（１）式に含まれる係数ａ、ｂ、ｃ、ｄの求め方に
ついて説明する。

（１）式は、 ＸＺ　−ａ　ＸＩ　＋　ｂ　Ｙｌ　＋　ｃ　　　　−・
”　（２）Ｙ２−−ｂＸ１＋ａＹ１＋ｄ　　　−−−−
−−（３）に分けられる。ここで第一のＸｌｉＩ画像上
のマーク１１ａ、　１２ａ、　１３ａ、　１４ａの各座
標をそれぞれ（Ｘ＋＋、Ｙ＋＋）、　　（Ｘ１２・　Ｙ
ｌ２）・　（Ｘｌ３・　Ｙｌ３）・　（Ｘ１４・Ｙ１４
）とし、第二のＸ線画像上のマーク１１ａ′、１２ａ’
　、１３ａ’　、１４ａ’の各座標をそれぞれ（ＸＺ。

Ｙ２１）、　　（Ｘ２２・Ｙｌ２）・　（Ｘ２３・Ｙｌ
３）・　（Ｘ２４１　　Ｙ２４）とする。このとき、（
２）式、（３）式より、 Ｘ、−ａＸ、十ｂＹ１１＋ｃ　　　”’−（２ａ）Ｘ２
□−ａＸ＋□＋ｂＹ、□＋ｃ　　　−＝−（２ｂ）ＸＺ
、　−ａ　Ｘ、３＋ｂ　Ｙ、十ｃ　　　−−−−（２ｃ
）Ｘ２４−ａＸ、４＋ｂＹ、４＋ｃ　　　・＝−＝（２
ｄ）Ｙ２．−−ｂＸ、、十ａＹ、＋ｄ　　・・・・・１
３ａ）Ｙ２２ｍ　　ｂｘ１２＋ａｙ１２＋ｄ　・＝−（
ａｂ）Ｙ、−−ｂＸ、、＋ａＹ、＋ｄ　　−＝−（３ｃ
）Ｙ２４−−ｂＸ、４＋ａＹ、４＋ｄ　　・−・・（３
ｄ）となる。ここで求めるべき係数はａ、ｂ、ｃ、ｄの
４つであるため、例えば上記（２ａ）、　　（２ｂ）　
；　（３ａ）。

（３ｂ）の４つの式、即ち各Ｘ線画像に２つずつマーり
を記録した場合でも係数ａ、ｂ、ｃ、ｄを求めることが
でき、これが従来の方法に相当するものであるが、これ
では前述したように、Ｘ線画像を読み取った際のサンプ
リング誤差により、偽画像を十分に防止するには不十分
な程度に大きな誤差をもった係数ａ、ｂ、ｃ、ｄが求め
られるおそれがある。そこでここでは最小二乗近似法に
より係数ａ、ｂ、ｃ、ｄが求められる。ここで（２ａ）
。

（２ｂ）、　　（２ｃ）、　　（２ｄ）　；　（３ａ）
、　　（３ｂ）、　　（３ｃ）、　　（３ｄ）式により
、 Ｓｌ　−Σ　（Ｘ２１−　（ａＸＨ＋ｂＹ＋＋＋ｃ）ｌ
　２・・・・・・（４） Ｓ２−ｍＩΣ　ｆＹ２＋−（ｂ　Ｘ１＋＋　ａ　ｙ、、
＋　ｄ）　）　　２・・・・・・（５） を求め、これら（４）、（５）式を最小にするように、
ａ　Ｓｌ　／　ａ　ａ　−０・・・・・（６）ａｓ、／
ａｂ−０−・・−・・（７） ａ　Ｓｌ　／　ａ　ｃ　−０・・・・・４ｇ）により、
ａ、ｂ、ｃを求め（この求められたａ。

ｂ、ｃをａｌ、ｂｌ、ｃｌとする）、またａ　Ｓ２　／
　ａ　ａ　−０−（９） ａｓ２／ａｂ−０−−−−・−（１０）ａ　Ｓ２　／　
ａ　ｄ　−０−・・・−・（１１）によりａ、ｂ、ｄを
求める（この求められたａ。

ｂ、ｄをａｇ、ｂｌ、　　ｄ２とすル）。コノヨウニし
て求められたａｌ　＊　　ｂｌ　、’１　　；　”Ｚ　
、ｂＺ　＋ｄ２を用い、第一のＸ線画像の座標をヘルマ
ート変換式 ％式％（１２） に従って変換することにより、第一のＸ線画像を第二の
Ｘ線画像に重ね合わせることができる。上記例では、各
Ｘ線画像にマークを４点記録し、これらのマークすべて
に基づいて最小二乗近似法によりヘルマート変換式の係
数ａ、　　ｂ、　　ｃ、　　ｄを求め、このようにして
求められた係数ａ１．Ｉｇ。

Ｃ１；ａｚ＋　　ｂｌ、ｄ２を用いて（１２）式にした
がってヘルマート変換を行なうようにしたため、従来の
方法と比べ、複数のＸ線画像をより高精度に重ね併せる
ことができ、従って偽画像の発生が一層確実に防止され
る。尚、上記実施例では、第一のＸ線画像の座標を変換
することにより該第−のＸ線画像を第二のＸ線画像に重
ね合わせるようにしたが、第二のＸ線画像の座標を変換
して該第二のＸ線画像を第一のＸ線画像に重ね合わせる
ようにしてもよく、第一のＸ線画像と第二のＸ線画像の
双方の座標を変換してこれら２つのＸ線画像を重ね合わ
せるようにしてもよいことはもちろんである。

以上のようにして２つのＸ線画像の寸法補正演算が行な
われた後、サブトラクション処理、即ち、２つの画像信
号ＳＯ！、ＳＯｚの上記補正演算後の画像信号をそれぞ
れＳＯｒ　　、ＳＯｚ’　としたとき、 ５ｌ−Ｗａ　−５Ｏ１−Ｗｂ−５Ｏ２’＋Ｃ・・・・・
・（１３） ただし、Ｗａ　、Ｗｂは重み付は係 数、Ｃはバイアス分を表わす。

に従って重み付は引き算が行なわれ、これにより、２つ
のＸ線画像の差の画像に対応する画像信号Ｓ１が生成さ
れる。この画像信号Ｓ１は画像処理表示装置３０のＣＲ
Ｔデイスプレィ３２に送られ、この画像信号Ｓ１に基づ
く可視画像（エネルギーサブトラクション画像）がＣＲ
Ｔデイスプレィ３２に再生表示される。尚、本体部３４
で実行される上記サブトラクション処理を行なう機能（
ハードウェアとソフトウェアとの組合せ）が本発明の演
算ユニットの一例と観念される。

尚、上記実施例ではヘルマート変換により、２つのＸ線
画像の倍率補正、回転補正、Ｘ方向、Ｙ方向の位置補正
（平行移動）の全てを同時に行なったが、これらを同時
に行なうことに代え、一つずつ順次に行なってもよいこ
とはいうまでもない。

以下に、これらを順次に行なう例について説明する。

先ず、第一のＸ線画像に対する第二のＸ線画像の倍率α
は、各Ｘ線画像にマークが４つ記録されているため、複
数衣められる。ここで第一のＸ線画像に記録されたマー
クｌｌａ、１２ａ、１３ａ、１４ａの座標を（Ｘ＋＋＋
　Ｙ＋＋）　　（ｉ　ｓ”ｌ、２４．４　）で表わし、
第二のＸ線画像に記録されたマークｌｌａ　　、１２ａ
’、１８ａ’　、１４ａ’の座標を（Ｘ　２１　、　Ｙ
　２υ　（ｉ　−１，２，３゜４）で表わすと、第一の
Ｘ線画像に記録された４つのマークｌｌａ、　１２ａ、
　１３ａ、　１４ａのうちの任意の２点をとり、この２
点と対応する第二のＸ線画像のマークｌｌａ　’　、１
２ａ’　、１３ａ’　、１４ａ’　ののうちの２点を用
いて、倍率αは、 α−ｉ（Ｘｌ−Ｘ２１）　２”（Ｙ２１−Ｙ２υ２）／
（（Ｘｚ−Ｘｕ）　２”（Ｙｕ−Ｙ＋＋）　２１但し、
ｉ、ｊ　−１，２，３，４；　ｉ　≠ｊ・・・・・・（
１４） として求められる。

この倍率αを用いて第一のＸ線画像をα倍に拡大するこ
とにより、第一のＸ線画像の寸法を第二のＸ線画像の寸
法に合わせることができる。

ここで各Ｘ線画像にマークが２つずつ記録されている場
合、倍率αは一義的に定まるが、これは前述した従来例
に相当し、偽画像を防止するには不十分な精度をもった
倍率αが求められるおそれがある。これに対し本実施例
では、（１３）式のｉ。

ｊに１〜４のうちの任意の数（ただしｉ≠ｊ）を当ては
めることにより１．６つの倍率α（α□。

α２．α３．α４．α５．α６とする）が求められるた
め、これら６つの倍率α１．α２．α３゜ａ４＋　　α
３．α、の平均値１を ａ　−（（２１＋（１２＋（２３＋ａ４　＋ａ５　＋ａ
６　）　／６・・・・・・（１５） として求め、第一のＸ線画像が１倍される。

また上記２つのＸ線画像の相対的な回転角（互いに対応
する線分、即ち（Ｘ＋＋、Ｙｒ＋）と（Ｘｚ。

Ｙ＋＋）とを結ぶ線分と（Ｘ２１．　Ｙ２１）と（Ｘｚ
ｌ。

Ｙ２．）とを結ぶ線分とのなす角をθとすると、このθ
は θ−ｔａｎ　’　ｌ（Ｙ　２．−Ｙ２υバＸ２ｌ−Ｘ２
．）１−ｊａｎ−１［（Ｙ　＋＋−Ｙ＋＋）／（Ｘ＋＋
−Ｘ＋＋）１・・・・・・（１６） と表わされ、第一のＸ線画像を角度θだけ回転させるこ
とにより、２つのＸ線画像の回転位置が合わせられる。

尚ここでも６つのθが求められ、回転角度を正確に求め
るため、それら６つのθの平均値θが求められ、この平
均角度θだけ回転される。

さらに２つのＸ線画像のＸ方向、Ｙ方向の位置ずれΔＸ
、ΔＹは、 ΔＸ−Ｘ２＋　　ＩＱ　’　（Ｘｚ　　Ｃｘ　）　’　
ＣＯ８θａ　・（Ｙ＋＋　　Ｃｙ　）　争ｓ１ｎ　Ｎ＋
Ｃｘ　ｌ・・・・・・（１７） ΔＹ−Ｙ２＋　　Ｉｎ　’　（Ｘｚ−Ｃｘ　）　’Ｓｉ
ｎθＵ　’　（Ｙ＋＋−Ｃｙ　）　・ｃｏｓθ＋ＣＹＩ
・・・・・・（ｌＢ） と表わされる。尚、（１７）、（１８）式において、Ｃ
ｘ。

ＣＹは回転中心のそれぞれＸ座標、Ｙ座標である。

ここでは（１７）式、（１８）式に従って４つのΔＸ　
ΔＹが求められるため、これらの平均値ΔＸ、ΔＹが求
められ、これらの平均値ΔＸ、ΔＹだけ第一のＸ線画像
のＸ方向、Ｙ方向の位置調整が行なわれる。

このようにして、第一のＸ線画像の寸法の調整、回転角
度調整、平行移動が行なわれ、これにより２つのＸ線画
像が高精度に重ね合わされ、偽画像の防止が図られる。

尚、上記実施例では、各Ｘ線画像にはマークが４つずつ
記録されたが、これらは４つに限られたものではなく、
３つ以上記録されていれば複数のＸ線画像を従来の方法
よりも高精度に重ね合わせることができ、従って従来の
方法よりも偽画像の防止が図られることとなる。

また上記各実施例は、２枚の蓄積性蛍光体シートを用い
た例であるが、蓄積性蛍光体層が複数層に形成された一
枚のシートを用いてもよいことはもちろんである。また
前述したように蓄積性蛍光体層は一層のみであってその
一層の表面付近、裏面付近に記録されたＸ線画像を読み
取るようにしてもよい。

また上記実施例は、蓄積性蛍光体シートを用いた例であ
るが、本発明の記録シートは蓄積性蛍光体シートに限ら
れるものではなくＸａフィルム（撮影に際して一般に増
感スクリーンと組合わされる）等であってもよい。

（発明の効果） 以上詳細に説明したように、本発明の放射線画像のエネ
ルギーサブトラクション方法および装置は、各放射線画
像の互いに対応する位置に３つ以上のマークを記録し、
これら３つ以上のマークの位置関係に基づいて寸法の補
正を行なうようにしたため、サブトラクション処理の行
なわれる複数の放射線画像を高精度に重ね合わせること
ができ、偽画像の発生が防止される。

【図面の簡単な説明】

第１図は、本発明のエネルギーサブトラクション装置の
記録ユニットの一実施例であるＸＩＩ撮影装置の概略図
、 第２図は、各蓄積性蛍光体シートに蓄積記録されたＸ線
画像を模式的に表わした図、 第３図は、本発明のエネルギーサブトラクション装置の
読取ユニットおよび演算ユニットの一実施例であるＸ線
画像読取装置および画像処理表示装置の斜視図である。 １・・・ＸＩＩ撮影装置　　　　２・・・Ｘ線管３、３
ａ、　３ｂ、　３ｃｍＸ線　　４・・・被写体５・・・
第一の蓄積性蛍光体シート ６・・・フィルタ ７・・・第二の蓄積性蛍光体シート １１．１２．１３，１４・・・マーカ １１ａ、１２ａ、１３ａ、１４ａ　；　ｌｌａ　’、１
２ａ’　　、１３ａ’　　、１４ａ’・・・マーク １９・・・回転多面鏡 ２３・・・光ガイド １６・・・レーザ光源 ２２・・・輝尽発光光 ２４・・・フォトマルチプライヤ ２５・・・ログアンプ ３０・・・画像処理表示装置 ２６・・・Ａ／Ｄ変換器 仔 ＞噸トー一

Claims (2)

【特許請求の範囲】
1. （１）互いに異なるエネルギー分布を有する放射線によ
   る互いに同一の被写体の複数の放射線画像を記録シート
   に記録し、該記録シートから前記複数の放射線画像を読
   み取って該複数の放射線画像のそれぞれを表わす複数の
   画像信号を得、該複数の画像信号に基づいて該複数の画
   像信号がそれぞれ担持する前記放射線画像の寸法が互い
   に同一となるように該寸法を補正し、該補正後の画像信
   号に基づいてサブトラクション処理を行なう放射線画像
   のエネルギーサブトラクション方法において、前記放射
   線画像の記録に当たり、前記複数の放射線画像上の互い
   に対応する位置に前記寸法の補正を行なうためのマーク
   が３つ以上ずつ記録された前記複数の放射線画像を記録
   し、 前記寸法の補正に当たり、前記複数の放射線画像のそれ
   ぞれに３つ以上ずつ記録された前記マークの前記放射線
   画像上の位置関係に基づいて、前記寸法の補正を行なう
   ことを特徴とする放射線画像のエネルギーサブトラクシ
   ョン方法。
2. （２）放射線を発する放射線源と、被写体が配置される
   被写体配置部と、互いに異なるエネルギー分布を有する
   放射線による複数の放射線画像を記録するための記録シ
   ートが配置される、前記被写体配置部を挟んで前記放射
   線源と対向する位置に設けられたシート保持部とからな
   る記録ユニット、前記記録ユニットにおいて前記被写体
   の放射線画像の記録の行なわれた前記記録シートから、
   前記被写体の複数の放射線画像のそれぞれを表わす複数
   の画像信号を得る読取りユニット、前記読取りユニット
   で得られた前記複数の画像信号に基づいて該複数の画像
   信号がそれぞれ担持する前記放射線画像の寸法が互いに
   同一となるように該寸法を補正する寸法補正ユニット、
   および前記寸法補正ユニットで寸法の補正の行なわれた
   後の放射線画像を担持する画像信号に基づいてサブトラ
   クション処理を行なう演算ユニットを備えた放射線画像
   のエネルギーサブトラクション装置において、前記記録
   ユニットが、前記複数の放射線画像上の互いに対応する
   位置に前記寸法の補正を行なうためのマークが３つ以上
   ずつ記録された前記複数の放射線画像を記録するもので
   あり、 前記寸法補正ユニットが、前記複数の放射線画像のそれ
   ぞれに３つ以上ずつ記録された前記マークの前記放射線
   画像上の位置関係に基づいて、前記寸法の補正を行なう
   ものであることを特徴とする放射線画像のエネルギーサ
   ブトラクション装置。