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

Abstract

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

Description

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

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

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

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

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

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

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

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

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

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

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

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

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

しかし撮影に際し、複数の記録シートを用いる場合にお
いて、これら複数の記録シートが互いに平行とならない
場合が生じることがあり、この場合上記特開平２−１０
４１７４号公報、特願平１−２４６９４５号に記載され
た方法では複数の放射線画像を正確に重畳させることが
できず、偽画像が発生してしまうおそれがある。

本発明は、上記事情に鑑み、撮影に際し複数の記録シー
トが互いに平行でなかった場合であっても複数の放射線
画像を正確に重ね合わせることができる放射線画像のエ
ネルギーサブトラクション方法、および該方法を採用し
た装置を提供することを目的とするものである。

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

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

（作　　用） 本発明の放射線画像のエネルギーサブトラクション方法
および装置では、各放射線画像に３つ以上のマークを記
録し、これら３つ以上のマークに基づいて放射線画像上
の互いに異なる２方向について互いに独立した補正量と
なるように寸法補正を行なうようにしたため、撮影に際
し複数枚の記録シートが互いに斜めに配置されていた場
合であってもこれら複数枚の記録シートに記録されてい
た放射線画像を正確に重畳させることができ、従って偽
画像の発生を防止することができる。

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

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

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

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

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

各蓄積性蛍光体シート５．７のほぼ全面に被写体４の各
Ｘ線画像４ａ、　４ａ’が蓄積記録され、またこれらの
Ｘ線画像４ａ、　４ａ’　とともにマーカ１１．１２゜
１３の像１１ａ、１２ａ、１３ａ　；　ｌｌａ　’　、
１２ａ’　、１３ａ’　　（以後この像をマークと称す
る。）が形成されている。

このマーク１１ａ、１２ａ、１３ａ　；　ｌｌａ　’　
、１２ａ’　、１３ａ’は蓄積性蛍光体シート５，７の
互いに対応する位置に記録されており画像信号を得た後
このマーク１１ａ、１２ａ、１３ａ　；　ｌｌａ　’　
、１２ａ’　、１３ａ’に基づいて２つのＸ線画像の相
対的な位置合わせが行なわれる。

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

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

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

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

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

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

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

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

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

次に（１）式に含まれる係数ａ、　　ｂ、　　ｃ、　　
ｄ、　　ｅ。

ｆの求め方について説明する。

（１）式は、 Ｘ　　−ａＸ　　＋ｂＹ　　＋ｅ　　　　−−−−・−
（２）Ｙ　　−ｃＸ　　＋ｄＹ　　＋　ｆ　　　　−−
−−−−（Ｓ）に分けられる。ここで第一のＸ線画像上
のマーク１１ａ、１２ａ、１３ａの各座標をそれぞれ（
Ｘｌ、、　Ｙｌｌ）　。

（Ｘ１２．　Ｙｌ２）　、　　（Ｘ１３．　Ｘ１３）と
し、第二のＸ線画像上のマークｌｌａ　’　、１２ａ’
　、１３ａ’の各座標をそれぞれ（Ｘ　２□、　Ｙ２、
）　、　　（Ｘ２２．　Ｙ２２）　、　　（Ｘ２３”　
２３）とする。このとき、（２）式、（３）式より、 Ｘ２ｌ−ａＸ１１＋ｂＹ１１＋ｅ　　・旧・・（２ａ）
Ｘ２２−ａＸ１２＋ｂＹ１２＋ｅ・・・・・・（２ｂ）
Ｘ２３−ａＸ１３＋ｂＹ１３＋ｅ・・・・・・（２ｃ）
Ｙｌ１−ＣＸ１１＋ｄＹｌ１＋ｆ・回・（３ａ）Ｙ２２
−ＣＸ１２＋ｄＹ１２＋ｆ　　　−旧−（３ｂ）Ｙ２３
−ＣＸ１３＋ｄＹｌ３＋ｆ　　　−旧−（３Ｃ＞となる
。ここで求めるべき係数はａ、ｂ、ｃ、ｄ。

ｅ、ｆの６つであるため、上記（２ａ）、　（２ｂ）、
　　（２ｃ）　；（３ａ）、　（３ｂ）、　（３ｃ）の
６本の式に基づいて求めることができ、 ””（Ｘ２１−Ｘ２２　）（Ｙ１１’１３）−（Ｘ２１
−Ｘ２３）（Ｙｌ、−Ｙ１□）ｌ　／ｌ（Ｘ１□−Ｘ１
□）（Ｙ１□−Ｙｌ３）−（Ｘ１□−Ｘ１３）（Ｙｌ、
−Ｙ１□）ｌ　　　　　−−−−−−（４）ｂ−（Ｘ２
１−Ｘ２２）　（Ｘ１１−Ｘ１３）−（Ｘ２１−Ｘ２Ｂ
）　（Ｘｌｌ−Ｘ１２）１／　（（Ｙｌｌ−Ｙｌ２）　
（Ｘ１１−Ｘ１３）−（Ｙｌｌ−Ｙｌ３）　（Ｘｌｌ−
Ｘ１２）１・・・・・・（５） ｃ−（Ｙｌ１−Ｙ２２）　（Ｙｌｌ−Ｙｌ３）−（Ｙｌ
１−Ｙ２３）　（Ｙｌｌ−Ｙｌ２）１／（（Ｘｌｌ−１
１２）　（Ｙｌｌ−Ｙｌ３）−（Ｘｌｌ−Ｘ１３）　（
Ｙｌｌ−Ｙｌ２）１・・・・・・（８） ｄ−（Ｙｌ１−Ｙ２２）　（ＸＩＬ−Ｘｌ３＞−（Ｙｌ
１−Ｙ２３）　（Ｘｌｌ−Ｘ１２））／（（Ｙｌｌ−Ｙ
ｌ２）　（Ｘ１１−Ｘ１３）−（Ｙｌｌ−Ｙｌ３）　（
Ｘｌｌ−Ｘ１２）１・・・・・・（７） ｅ−（（Ｘ２１ＸＩ２−ＸＩＩＸ２２）　（Ｙ１２Ｘ１
３−Ｙ１３Ｘ１２）−（Ｘ２２ＸＩ３−Ｘ２３Ｘ１２）
　（Ｙ１１Ｘ１２−Ｙ１２Ｘ１１）１／　（（Ｘ１２−
Ｘｌｌ）　（Ｙ１２Ｘ１３−Ｙ１３Ｘ１２）−（Ｘ１３
−Ｘ１２）　（Ｙ１１Ｘ１２−Ｙ１２Ｘ１１）１・・・
・・・（８） ｒ−（（Ｙ２１Ｘ１２−ＸＩ　ＬＹ２２）　（ＹＬ２Ｘ
１３−Ｙ１３Ｘ１２）−（Ｙ２２ＸＬ３−Ｘ２３Ｘ１２
）　（Ｙ１１Ｘ１２−Ｙ１２Ｘ１１）１／　（（Ｘ１２
−Ｘｌｌ）　（Ｙ１２Ｘ１３−Ｙ１３Ｘ１２＞−（Ｘ１
３−Ｘ１２）　ｍ１Ｘ１２−Ｙ１２Ｘｌ　１）　１・・
・・・・（９） となる。

このようにして（４） 〜（９） 式にしたがって求め られた係数ａ、ｂ、ｃ、ｄ、ｅ、ｆを用い、（１）式に
従って座標変換を行なうことにより、第一のＸ線画像を
第二のＸ線画像に重ね合わせることができる。

上記例では、各Ｘ線画像にマークを３点記録し、これら
のマークに基づいてアフィン変換を行なうようにしたた
め、Ｘ方向とＹ方向とで互いに独立に寸法補正が行なわ
れ、第１図に示すように第一の蓄積性蛍光体シート５と
第二の蓄積性蛍光体シート７とが互いに傾いた状態で撮
影が行なわれた場合であっても、複数の放射線画像を高
精度に重ね合わせることができ、偽画像の発生が防止さ
れる。

ここで上記実施例では各Ｘ線画像に３つずつマークが付
されていたが、このマークの数は３つに限られるもので
はなく、４つ以上ずつ付されていてもよい。この場合、
４つ以上ずつのマークのうちの任意の３つずつを用いて
（４）〜（９）式にしたがって演算を行なうことにより
、（１）式のアフィン変換の各係数ａ、　　ｂ、　　ｃ
、　　ｄ、　　ｅ、　　ｆのそれぞれとして複数の値が
求められるが、これら求められた複数の値を用いて平均
値ａ、ｂ、ｃ、ｄ、ｅ。

Ｔ求め、この平均値ａ、ｂ、τ、ｄ、ｅ、了を用いて（
１）式に従ってアフィン変換を行なうことにより、より
高精度の重ね合わせを行なうことができる。尚、上記平
均値に代えて中央値、　（最大値＋最小値）／２等を求
め、この求めた値を用いてアフィン変換を行なってもよ
いことはもちろんであり、また最小二乗近似法等を用い
て各係数値。

ａ、　ｂ、　　ｃ、　　ｄ、　　ｅ、　　ｆを求めても
よいことももちろんである。

また、上記実施例では、第一のＸ線画像の座標を変換す
ることにより該第−のＸ線画像を第二のＸ線画像に重ね
合わせるようにしたが、第二のＸ線画像の座標を変換し
て該第二のＸ線画像を第一のＸ線画像に重ね合わせるよ
うにしてもよく、第一のＸ線画像と第二のＸ線画像の双
方の座標を変換してこれら２つのＸ線画像を重ね合わせ
るようにしてもよいことはもちろんである。

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

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

尚、上記実施例ではアフィン変換により、２つのＸ線画
像のＸ方向とＹ方向とで互いに独立した寸法補正、回転
補正、Ｘ方向、Ｙ方向の位置補正（平行移動）の全てを
同時に行なったが、これらを同時に行なうことに代え、
一つずつ順次に行なってもよいことはいうまでもない。

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

先ず、第一のＸ線画像に対する第二のＸ線画像のＸ方向
の倍率α　は、マーク１２ａとマーク１３ａとを結ぶ線
分、およびマーク１２ａ　’　とマーク１３ａ　’　と
を結ぶ線分がＸ方向に延びる線分であることから、これ
らのマーク１２ａ　、　１３ａ　；　１２ａ　’１３ａ
′　を用い、これらのマーク１２ａ　、　１３ａ　；１
２ａ　’　１３ａ　’　の座標を（Ｘ１２．　Ｙ１２）
　、　　（Ｘ１３゜Ｙ　　）　　；　　（Ｘ２２．　Ｙ
２２）　、　　（Ｘ２３．　Ｙ２３）としたとき、 ・・・・・・（１１） として求められる。

また第一のＸ線画像に対する第二のＸ線画像のＹ方向の
倍率α　は、マークｌｌａとマーク１２ａとを結ぶ線分
、およびマークｌｌａ　’　とマーク１２ａとを結ぶ線
分がＹ方向に延びる線分であることから、これらのマー
クｌｌａ　、　１２ａ　　；　Ｌｌａ　’　１２ａ　’
　を用い、これらのマークｌｌａ　、　１２ａ　；　Ｌ
ｌａ　’　１２ａ　’の座標を（Ｘ　　　Ｙ　　）、　
　（Ｘ１２．　Ｙｌ。）　　；　　（Ｘ１１’　　　１
１ ・・・・・・（１２） として求められる。

これらの倍率α　、α　を用いて第一のＸ線画ｙ 像をＸ方向にα　倍、Ｙ方向にα　倍することにｙ より、第一のＸ線画像の縦横の寸法を第二のＸ線画像に
合わせることができる。尚、各Ｘ線画像に４つずつ以上
のマークを付し、上記倍率α　、α、の６値として複数
の値を求め、これら複数の値の平均値ａ、ａ　　等を求
め、この平均値等を用ｙ いて寸法補正してもよいことは前述した実施例と同様で
ある。

また上記２つのＸ線画像の相対的な回転角（互いに対応
する線分、即ち（Ｘ　　、　Ｙ　、）と（Ｘ１１Ｌｊ　
　　　ＩＪ Ｙ１ρとを結ぶ線分と（Ｘ２ｊ、Ｙ２ｊ）と（Ｘ２ｉ’
Ｙ２１）とを結ぶ線分とのなす角をθとすると、このθ
は θ−ｔａｎ　’　Ｉ（Ｙ２ｊ−Ｙ２ｉ）／（Ｘ２ｊ−Ｘ
２１）１−ｊａｎ　’　１（Ｙｌｊ−Ｙｌｉ）／（Ｘｌ
ｊ−Ｘ１１）１（ｔ、、ｒｌ、２．３；　ｉ≠ｊ） ・・・・・・（１３） と表わされ、第一のＸ線画像を角度θだけ回転させるこ
とにより、２つのＸ線画像の回転位置が合わせられる。

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

さらに２つのＸ線画像のＸ方向、Ｙ方向の位置ずれΔＸ
、ΔＹは、 Δｘ−ｘ−＋α・（Ｘ１１−Ｃｘ）・ＣＯ５θ−ａ　　
・（Ｙ　　−Ｃ）　・ｓｉｎθ＋Ｃｘ）ｘ　　　　ｌｉ
　　　Ｙ ・・・・・・（１４） ΔＹ−Ｙ−ｆα−（Ｘ１１−Ｃｘ）・ｓｉｎθｉ −ａ　　　（Ｙ　　−Ｃ）　・ｃｏｓθ＋ＣＹ）ｙ　　
　　　　ｌｉ　　　　Ｙ ・・・・・・（１５） と表わされる。尚、（１４）、（１５）式において、Ｃ
Ｘ。

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

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

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

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

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

【図面の簡単な説明】

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

Claims (2)

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