JPH03276264A - エネルギーサブトラクション画像表示装置 - Google Patents

Abstract

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

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明は、同一被写体のエネルギーサブトラクション画
像のそれぞれを表わす二つの画像データに基づいてエネ
ルギーサブトラクション画像を求め、該エネルギーサブ
トラクション画像を可視画像として表示するエネルギー
サブトラクション画像表示装置に関するものである。

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

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

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

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

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

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

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

このエネルギーサブトラクションの例としては、例えば
人体の胸部のように軟部および骨部から構成された被写
体に互いにエネルギーの異なる放射線を照射して複数の
放射線画像を得、これら複数の放射線画像を読み取って
これら複数の放射線画像のそれぞれを表わす複数の画像
データを得、これら複数の画像データに基づいてエネル
ギーサブトラクション処理を行なって被写体の主として
軟部が記録された軟部画像を表わす軟部画像データもし
くは被写体の主として骨部が記録された骨部画像を表わ
す骨部画像データを求め、求められた軟部画像もしくは
骨部画像を観察の対象とする場合がある。この軟部画像
、骨部画像はそれぞれ骨部、軟部の陰影が消去された画
像であるため、骨部もしくは軟部に隠れてしまっていた
陰影や骨部もしくは軟部の影響で見にくくなってしまっ
ていた陰影を浮かび上がらせることができる。

またエネルギーサブトラクションの他の例としては、例
えば人体の乳房を被写体とした場合の、乳腺の陰影が強
調された画像、悪性腫瘍が強調された画像等も観察対象
とされている。

（発明が解決しようとする課題） 上記のようなエネルギーサブトラクション画像を得る際
に、該エネルギーサブトラクション画像のコントラスト
や画像全体の濃度（ＣＲＴデイスプレィに画像を表示す
る場合の、表示された画像の輝度を含む）が適切なもの
となるようにサブトラクション演算を行なう必要があり
、これまでは適切なコントラスト、濃度レベルを試行錯
誤的に調節していた。

ところがエネルギーサブトラクション画像を観察するに
際し、サブトラクション処理前の複数の放射線画像の重
みづけを順次変更してサブトラクション処理を行ない、
この結果得られた順次変更されたエネルギーサブトラク
ション画像をＣＲＴデイスプレィ等に順次表示して観察
を行なう場合がある。

この場合、各重みづけに対応してその都度コントラスト
や濃度レベルを適切に調節しないと観察しやすい画像と
ならず、重みづけ変更する都度試行錯誤的に調整するこ
ととするときわめて頻雑であり使い勝手の悪い画像表示
装置となってしまうという問題がある。

本発明は、上記事情に鑑み、重みづけを順次変更したエ
ネルギーサブトラクション画像を求める際に適切なコン
トラスト、濃度が自動的に設定され、したがって観察し
やすい使い勝手のよいエネルギーサブトラクション画像
表示装置を提供することを目的とするものである。

（課題を解決するための手段） 本発明の第一のエネルギーサブトラクション画像表示装
置は、 同一被写体のエネルギーサブトラクション画像のそれぞ
れを表わす二つの画像データをＳ１，Ｓ２、第一および
第二の重み係数をそれぞれＷ１，Ｗ２、バイアス分を０
５サブトラクシジン画像データを８３としたとき、次の
（１）式 ％式％（１） に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
て、順次変更されたサブトラクション画像データＳ３を
求める画像演算部と、該順次変更されたサブトラクショ
ン画像データＳ３に基づく順次変更されたサブトラクシ
ョン画像データＳ３に基づく順次変更されたサブトラク
ション画像表示装置において、 前記画像演算部が、 次の（２１式 ％式％（２１ に基づいて差分画像データＳＯを求める差分画像演算手
段と、 前記差分画像データＳＯのばらつきの程度を指標する第
一の特徴量と前記差分画像データＳＯの平均的な値を指
標する第二の特徴量とを求める特徴量演算手段と、 前記第一の特徴量に基づいて前記第二の重み係数Ｗ２を
求める重み係数演算手段と、 前記第二の特徴量に基づいて前記バイアス分Ｃを求める
バイアス演算手段と、 前記重み係数演算手段で求められた第二の重み係数Ｗ２
と前記バイアス演算手段で求められたバイアス分Ｃとを
用いて、上記（１）式に従ってサブトラクション画像デ
ータＳ３を求めるサブトラクション画像演算手段とを備
えたことを特徴とするものである。

また本発明の第二のエネルギーサブトラクション画像表
示装置は、上記第一のエネルギーサブトラクション画像
表示装置において、 各画素毎の前記差分画像データＳＯの最大値および最小
値をそれぞれＳ　０ｍａｘ　、　　Ｓ　０ｍ１ｎとし、
二つの定数をＣ１，Ｃ２としたとき、 前記第一の特徴量および前記第二の特徴量が、それぞれ
Ｓ　０ＩＩａｘ　−Ｓ　Ｏａ＋ｉｎ　、　　（Ｓ　０ｍ
ａｘ　＋　Ｓ　０ｉｉｎ）／２であり、 前記重み係数演算手段が、前記第二の重み係数Ｗ２を次
の（３）式 ％式％）（３） に従って求めるものであり、かつ 前記バイアス演算手段が前記バイアス分Ｃを次の（４）
式 ％式％） に従って求めるものであることを特徴とするものである
。

また本発明の第三のエネルギーサブトラクション画像表
示装置は、上記第一のエネルギーサブトラクション画像
表示装置において、 前記重み係数演算手段が、前記第一の特徴量とともに前
記重み係数演算手段が、前記第一の重み係数Ｗ２を求め
るものであり、 前記バイアス演算手段が、前記第二の特徴量とともに前
記重み係数比αに基づいて前記バイアス分Ｃを求めるも
のであることを特徴とするものである。

また、本発明の第四のエネルギーサブトラクション画像
表示装置は、上記第三のエネルギーサブトラクション画
像表示装置において、 各画素毎の前記差分画像データＳＯの最大値および最小
値をそれぞれＳ　０ｍａｘ　、　　Ｓ　Ｏｍｊｎ　％前
記重み係数比αを変数とした二つの関数をＦ（α）。

Ｇ（α）としたとき、 前記第一の特徴量および前記第二の特徴量がそれぞれＳ
　０ｍａｘ　−Ｓ　０ｍ１ｎ　、　　（Ｓ　０ｍａｘ　
＋　Ｓ　Ｏａ＋ｉｎ）／２であり、 前記重み係数比αを変数とした二つの関数を次の（５）
式 ％式％） に従って求めるものであり、かつ 前記バイアス演算手段が、前記バイアス分Ｃを次の（６
）式 ％式％） （６） に従って求めるものであることを特徴とするものである
。

また、本発明の第五のエネルギーサブトラクション画像
表示装置は、 同一被写体のエネルギーサブトラクション画像のそれぞ
れを表わす二つの画像データをＳ１，Ｓ２、第一および
第二の重み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ
、サブトラクション画像データを８３としたとき、次の
（７）式 ％式％（７） に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
て、順次変更されたサブトラクション画像データＳ３を
求める画像演算部と、該順次変更されたサブトラクショ
ン画像データＳ３に基づく順次変更されたサブトラクシ
ョン画像データＳ３に基づく順次変更されたサブトラク
ション画像表示装置において、 前記画像演算部が、 前記重み係数比αを変数とした前記第二の重み係数Ｗ２
を表わす関数Ｗ２　（αを変数とした前記重み係数比α
を変数とした前記バイアス分Ｃを表わす関数Ｃ（α）と
を記憶しておく記憶手段と、前記記憶手段から読み出さ
れた第二の重み係数Ｗ２とバイアス分Ｃとを用いて、上
記（７）式に従ってサブトラクション画像データＳ３を
求めるサブトラクション画像演算手段とを備えたことを
特徴どするものである。

また、本発明の第六のエネルギーサブトラクション画像
表示装置は、 同一被写体のエネルギーサブトラクション画像のそれぞ
れを表わす二つの画像データをＳ１，Ｓ２、第一および
第二の重み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ
、サブトラクション画像データを８３としたとき、次の
（８）式 ％式％（８１ に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
て、順次変更されたサブトラクション画像データＳ３を
求める画像演算部と、該順次変更されたサブトラクショ
ン画像データｓ３に基づく順次変更されたサブトラクシ
ョン画像データＳ３に基づく順次変更されたサブトラク
ション画像表示装置において、 前記画像演算部が、 次の（９）式 ％式％（９） に基づいて差分画像データＳＯを求める差分画像演算手
段と、 前記差分画像データＳＯのばらつきの程度を指標する第
一の特徴量を求める特徴量演算手段と、前記第一の特徴
量に基づいて前記第二の重み係数Ｗ２を求める重み係数
演算手段と、 前記重み係数比αを変数とした前記バイアス分Ｃを表わ
す関数Ｃ（α）を記憶しておく記憶手段と、 前記重み係数演算手段で求められた第二の重み係数Ｗ２
と前記記憶手段から読み出されたバイアス分Ｃとを用い
て、上記（８）式に従ってサブトラクション画像データ
Ｓ３を求めるサブトラクション画像演算手段とを備えた
ことを特徴とするものである。

さらに、本発明の第七のエネルギーサブトラクション画
像表示装置は、 同一被写体のエネルギーサブトラクション画像のそれぞ
れを表わす二つの画像データをＳ１，Ｓ２、第一および
第二の重み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ
，サブトラクション画像データを８３としたとき、次の
（ｌＯ）式％式％（） に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
て、順次変更されたサブトラクション画像データＳ３を
求める画像演算部と、該順次変更されたサブトラクショ
ン画像データＳ３に基づく順次変更されたサブトラクシ
ョン画像データＳ３に基づく順次変更されたサブトラク
ション画像表示装置において、 前記画像演算部が、 次の第（１１）式 ％式％（１１） に基づいて差分画像データＳＯを求める差分画像演算手
段と、 前記差分画像データＳＯの平均的な値を指標する第二の
特徴量を求める特徴量演算手段と、前記第二の特徴量に
基づいて前記バイアス分Ｃを求めるバイアス演算手段と
、 前記重み係数比αを変数とした前記第二の重み係数Ｗ２
を表わす関数Ｗ２　（α）を記憶しておく記憶手段と、 前記記憶手段から読み出された第二の重み係数Ｗ２と前
記バイアス演算手段で求められたバイアス分Ｃとを用い
て、上記（ｌＯ）式に従ってサブトラクション画像デー
タＳ３を求めるサブトラクション画像演算手段とを備え
たことを特徴とするものである。

（作　　用） 本発明の第一のエネルギーサブトラクション画像表示装
置は、上記（２）式に基づいて差分画像データＳＯを求
めた際に該差分画像データｓｏのばらつきの程度と平均
的な値とをそれぞれ指標する第一の特徴量と第二の特徴
量とを求め、上記第一の特徴量に基づいて第二の重み係
数Ｗ２を求めるとともに上記第二の特徴量に基づいてバ
イアス分Ｃを求めるようにしたものである。

ここで、上記第一の特徴量は差分画像データＳＯのばら
つきの程度を指標するものであるため、これに基づいて
適切なコントラストを求めることが可能となる。ここで
は上記第一の特徴量に基づいて第二の重み係数Ｗ２を求
めているがこのことがコントラストを求めることに対応
している。即ち、重み係数比α＝Ｗ１／Ｗ２であること
から第一の重み係数Ｗ１はＷｌ−α・Ｗ２となる。ここ
でαの値はこれを先ず定めた後に演算を行なうものであ
るから上記（１）式に基づく各サブトラクシジン演算処
理を行なう際はαは既に設定されており、したがって上
記第二の重み係数Ｗ２を定めると一義的に上記第一の重
み係数Ｗ１も定まることを意味し、上記（１）式におい
てバイアス分Ｃを除くＳ３’　−Ｗｌ−Ｓｌ−Ｗ２・８
２　　　・・・（１２）の演算が可能となり、したがっ
てサブトラクション画像のコントラストが定まることと
なる。この第一の特徴量としては、本発明の第二のエネ
ルギーサブトラクション画像表示装置における。第一の
特徴量Ｓ　０ｓａｘ　−Ｓ　０Ｉｌｉｎであってもよい
が、これに限られるものではなく、例えば差分画像デー
タＳＯの分散値等該差分画像データＳＯのばらつきの程
度を指標するものであればよい。

また上記第二の特徴量は差分画像データＳＯの平均的な
値を指標するものであるため、これに基づいて適切な平
均濃度を求めることができる。ここではバイアス分Ｃを
定めているため、これによりサブトラクション画像全体
の濃度レベルが定まることになる。

この第二の特徴量としては本発明の第二のエネルギーサ
ブトラクション画像表示装置における第二の特徴量（Ｓ
　０ｔｉａｘ　＋　Ｓ　Ｃ１＋ｉｎ　）　／　２を採用
してもよいが、これに限られるものではなく差分画像デ
ータＳＯの平均値、中央値（メジアン値）等差分画像デ
ータの平均的な値を指標するものであればよい。

本発明の第一のエネルギーサブトラクション画像表示装
置は、重み係数比α＝Ｗ１／Ｗ２を順次変換して上記の
ようにして第一および第二の重み係数Ｗ１，Ｗ２．およ
びバイアス分Ｃを求めて上記（１）式に従ってサブトラ
クション画像データＳ３を求めるようにしたものである
ため、順次求められたサブトラクション画像データＳ３
は、常に適切なコントラストと適切な濃度レベルを有す
るサブトラクション画像を得ることかできる。

また本発明の第二のエネルギーサブトラクション画像表
示装置は、上記第一のエネルギーサブトラクション画像
表示装置の実施態様のひとつであり、上記第一の特徴量
、第二の特徴量としてそれぞれ５Ｏ１１ａｘ　−８０ｍ
１ｎ　、　　（Ｓ０ｍａｘ　＋ＳＯａ＋ｉｎ　）／２を
選択し、上記第二の重み係数Ｗ２．バイアス分Ｃを求め
るについてそれぞれ上記（３）式、（４）式を採用した
ものであり、これにより適切なコントラストと適切な濃
度レベルを有するサブトラクション画像を表示すること
ができる。

次に本発明の第三のエネルギーサブトラクション画像表
示装置の意義について、前述した胸部Ｘ線画像のエネル
ギーサブトラクションを例に説明する。

複数の胸部Ｘ線画像をサブトラクション処理することに
より得られた、被写体の軟部のみが抽出された軟部画像
は、元来ダイナミックレンジの広い画像であるため、こ
の軟部画像を例えばＣＲＴデイスプレィに可視画像とし
て表示する場合、その軟部画像のダイナミックレンジを
ＣＲＴデイスプレィに表示可能な最大のダイナミックレ
ンジ程度にまで広げることが好ましく、一方被写体の骨
部のみが抽出された骨部画像は元来ダイナミックレンジ
の狭い画像であるため、この骨部画像をＣＲＴデイスプ
レィに可視画像として表示する場合、その骨部画像のダ
イナミックレンジをＣＲＴデイスプレィに表示可能な最
大のダイナミックレンジ程度まで広げてしまうとコント
ラストが強すぎる画像となってしまい、かえって観察適
性の低下を招くことになる場合がある。また画像を表示
する場合に、例えば観察者の習慣や好み等により、画像
全体をより白っぽく（濃度を低く）表示したりより黒っ
ぽく　（濃度を高く）表示したりすることを望む場合が
ある。

本発明の第三のエネルギーサブトラクション画像表示装
置は、上記第一のエネルギーサブトラクション画像表示
装置において、第一の特徴量とともに重み係数比αに基
づいて第二の重み係数Ｗ２を求め、かつ第二の特徴量と
ともに重み係数比αに基づいてバイアス分Ｃを求めて上
記（１）式に従ってサブトラクション画像データＳ３を
求めるようにしたため、上記第一のエネルギーサブトラ
クション画像表示装置の作用効果はそのまま保存される
とともに、重み係数比αに応じてダイナミックレンジお
よび濃度レベルが適切に変化した画像を表示することが
できる。

また本発明の第四のエネルギーサブトラクション画像表
示装置は、上記第三のエネルギーサブトラクション画像
表示装置の実施態様のひとってあり、上記第一の特徴量
、第二の特徴量としてそれぞれＳ　０ａａｘ　−５０ｍ
１ｎ　、　　（Ｓ　Ｏａ＋ａｘ　＋　Ｓ　０ｍ１ｎ　）
／２を選択し、上記第二の重み係数Ｗ２．バイアス分Ｃ
を求めるについてそれぞれ第一の特徴量と重み係数比α
との関数である（５）式、第二の特微量と重み係数比α
との関数である（６）式を採用したものであり、これに
より画像が順次表示されるにつれてダイナミックレンジ
、濃度レベルを適切に順次変化させることができる。尚
、上記（５）式、（６）式に挿入されている関数Ｆ（α
）、Ｇ（α）は、特に特定の形状の関数に限定されるも
のではなく、順次変化させて表示する際の両極に位置す
る２つの画像としてどのような画像を選択するかという
ことや観察者の好み等により適宜定められるものである
。

以上の第一〜第四のエネルギーサブトラクション画像表
示装置は、差分画像データＳＯを求め、この差分画像デ
ータＳＯに基づいてその後の処理を行なうように構成さ
れた装置であるが、差分画像データＳＯを求めることは
本発明の目的を達成するために必ずしも必要なことでは
ない。

すなわち、例えば人体の胸部Ｘ線画像のみを対象とする
装置等限定された用途にのみ用いるエネルギーサブトラ
クション画像表示装置においては、重み係数比αに対す
る差分画像データＳＯのばらつきの概略値と概略平均値
の値、したがって第一の特徴量、第二の特徴量の概略値
を経験的、実験的に求めておくことかできる。したがっ
て各サブトラクション画像データＳ３を求める過程で差
分画像データＳＯを求め、この求めた差分画像データＳ
Ｏに基づいて第一および第二の特徴量を求めることに代
え、経験的、実験的に求められた第一および第二の特徴
量を採用することが可能である。

本発明の第五のエネルギーサブトラクション画像表示装
置は上記観点に基づくものであり、上記第一〜第四のエ
ネルギーサブトラクション画像表示装置が差分画像デー
タＳＯを求めるものであることに代え、重み係数比αを
変数とした第二の重み係数Ｗ２を表わす関数Ｗ２　（α
）と重み係数比αを変数としたバイアス分Ｃを表わす関
数Ｃ（α）とをあらかじめ求めて記憶しておいて、この
記憶された第二の重み係数Ｗ２とバイアス分Ｃとを用い
て上記（７）式（（１）式と同一）に従ってサブトラク
ション画像データＳ３を求めるようにしたため、上記第
一のエネルギーサブトラクション画像表示装置はど厳密
ではないか実用上十分な精度で適切なコントラスト、濃
度レベルを備えたサブトラクション画像を表示すること
ができる。またこれらの関数Ｗ２（α）、Ｃ（α）の関
数形により上記第三のエネルギーサブ］・ラクション画
像表示装置のように表示される画像が順次変更されるに
従ってコントラスト、濃度レベルを順次適切に変化させ
ることも可能となる。

上記第一〜第四のエネルギーサブトラクション画像表示
装置は、差分画像データＳＯを求めることにより第二の
重み係数Ｗ２とバイアス分Ｃとを求めるものであり、上
記第五のエネルギーサブトラクション画像表示装置は第
二の重み係数Ｗ２とバイアス分Ｃとを重み係数比αの関
数Ｗ２　（α）。

Ｃ（α）としてＳ己憶しておくものであり、いずれも本
発明の目的を達成し得るものであることから、これらの
組み合わせも可能となる。即ち、本発明の第六のエネル
ギーサブトラクション画像表示装置は、差分画像データ
ＳＯを求めることにより第二の重み係数Ｗ２を求めると
ともにバイアス分Ｃについては重み係数比αの関数Ｃ（
α）として記憶しておくものであり、また本発明の第七
のエネルギーサブトラクション画像表示装置は、第二の
重み係数Ｗ２は重み係数比αの関数Ｃ（α）として記憶
しておきバイアス分Ｃは差分画像データＳＯを求めるこ
とにより求めるようにしたものであり、いずれの場合も
上記第一〜第五のエネルギーサブトラクション画像表示
装置と同様に適切なコントラスト、濃度レベルを有する
サブトラクション画像を表示することができる。

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

第１図は、Ｘ線撮影装置の概略図である。

このＸ線撮影装置１のＸ線管２から発せられたＸ線３に
より被写体（人体の胸部）４が照射される。被写体４を
透過したＸ線３ａは第一の蓄積性蛍光体シート５に照射
され、Ｘ線３ａのエネルギーのうち比較的低エネルギー
のＸ線が該第一の、蓄積性蛍光体シート５に蓄積され、
これにより該シート５に被写体４のＸ線画像が蓄積記録
される。シート５を透過したＸ線３ｂはさらに低エネル
ギーのＸ線をカットするフィルタ６を透過し、該フィル
タ６を透過した高エネルギーＸ線３Ｃが第二の蓄積性蛍
光体シート７に照射される。これにより該シート７にも
被写体４のＸ線画像が蓄積記録される。

被写体４には、サブトラクション処理を行なうにあたっ
て２つのＸ線画像の位置合わせを行なうための基準とな
る２つのマーク８が付されている。

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

第２図は、Ｘ線画像読取装置と本発明のエネルギーサブ
トラクション画像表示装置の一例を内包した画像処理表
示装置の斜視図である。

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

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

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

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

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

第３図は、画像処理表示装置内の内部メモリに記憶され
た第一および第二のＸ線画像を表わす２つの画像信号Ｓ
１，Ｓ２に基づいて、該画像処理表示装置内で行なわれ
る処理の流れを表わした図である。

画像処理表示装置内の内部メモリに記憶された、第一お
よび第二のＸ線画像信号Ｓ１，Ｓ２は、第３図に示すそ
れぞれ第一のＸ線画像４１．第二のＸ線画像４２を担持
する信号である。第一のＸ線画像４１は比較的低エネル
ギーＸ線による画像であり、第二のＸ線画像４２は比較
的高エネルギーＸ線による画像であるが、互いに軟部お
よび骨部の濃度は異なるものの両者ともこれら軟部およ
び骨部の双方が記録された画像である。

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

次に位置合わせ後の２つのＸ線画像の各画素毎に式 ％式％（１３） ただし、Ｗ１，Ｗ２は重み係数、 Ｃはバイアス分を表わす に従って演算を行なってサブトラクション画像信号Ｓ３
が求められる。

ここで上記（１３）式を変形すると、 Ｓ３−Ｗ２・（Ｗ１／Ｗ２・５ｌ−３２）＋Ｃ・・・（
１４） となりＷｌ／Ｗ２をαと置くと、 Ｓ３−Ｗ２・　（α・５ｌ−５２）　十Ｃ・・・（１５
）と書き表わすことができる。ここでαは２つの重み係
数Ｗ１，Ｗ２の比Ｗｌ／Ｗ２であり、ここではこのαを
重み係数比と呼ぶこととする。

ここで上記（１５）式（上記（１３）式と同等）に従っ
て演算を行なってサブトラクション画像信号Ｓ３を求め
るにあたり、上記（１５）式のかっこ内、即ちＳＯ−α
・５ｌ−Ｓ２　　　　　　　　　・・・（１６）の演算
を各画素毎に行なって差分画像信号ＳＯか求められる。

尚本実施例では、画像処理装置３０における、上記（Ｉ
Ｂ）式に従って演算を行なう機能が本発明にいう差分画
像演算手段の一例と観念される。この差分画像信号ＳＯ
は各画素毎に求められたものであり、各画素毎に一つの
信号が対応している。

次に、これら各画素に対応する多数の差分画像信号ＳＯ
のうちの最大値５ＯＩＩａｘ、最小値５０ｍ１ｎを求め
る。ここでは求められた最大値ＳＯｏ＋ａｘ。

最小値５０ｍ１ｎ、、さらに重み係数比αの関数ｆ（α
）１ｇ（α）を用いて、（１５）式の重み係数Ｗ２と、
バイアス分Ｃが求められる。　第４Ａ図、第４Ｂ図は、
重み係数比αの関数ｆ（α）９ｇ（α）の−例を表わし
たグラフである。

ここではこれら第４Ａ図、第４Ｂ図に示した関数ｆ（α
）１ｇ（α）を用いて、式 ％式％） （１７） ） （１８） を用いて重み係数Ｗ２．バイアス分Ｃが求められる。こ
こで、（１７）式中の１０２３はＣＲＴデイスプレィ３
２（第２図参照）に可視画像を表示する際の濃度（輝度
）分解能を表わしており、ｆ（α）−１の場合、画像上
の最大値５ＯＩＩａｘを有する部分が表示可能な最大濃
度（最小輝度）で表示され、画像上の最小値ＳＯｓ＋ｉ
ｎを有する部分が表示可能な最小濃度（最大輝度）で表
示されることを意味している。

このように上記（１７）、　（ｉｌｌ）式で求められた
重み係数Ｗ２とバイアス分Ｃを用い上記（１５）式に基
づいてサブトラクション画像信号Ｓ３が求められる。

求めたサブトラクション画像信号Ｓ３はＣＲＴデイスプ
レィ３２に送られ、ＣＲＴデイスプレィ３２にはこのサ
ブトラクション画像信号Ｓ３に基づく可視画像が表示さ
れる。

ここで、α＋１は被写体中の主として骨部が記録された
骨部画像４３（第３図参照）を表わしαを１から２に近
づけると骨部と軟部との両者が記録された画像４４とな
り、α−２では被写体の軟部のみが記録された軟部画像
４５となる。

骨部画像４３、即ちα崎１の場合、第４Ａ図に示すよう
にｆ　（α）　＋０．３である。前述したように骨部画
像は元来ダイナミックレンジの狭い画像であ、るため、
ｆ（α）−１、即ち骨部画像のダイナミックレンジをＣ
ＲＴデイスプレィ３２に表示可能な最大のダイナミック
レンジに調整するとコントラストが強すぎ、かえって観
察適性の低下を招くため、ｆ　（α）−０，３と比較的
小さな値としたものである。また、骨部画像４３の場合
ｇ（α）′−。

４３０と比較的小さい値を有している。この場合（Ｉ８
）式よりバイアス分Ｃか小さい値となり、全体として白
っぽい画像としてＣＲＴデイスプレィ３２に表示される
ことになる。

一方、軟部画像４５、即ちα＋２の場合、ｆ（α）−０
，８６，ｇ　（α）−ｅｏｏとなる。したがって軟部画
像４５の場合は比較的高コントラストかつ全体の濃度（
輝度）レベルが中程度の可視画像としてＣＲＴデイスプ
レィ３２に表示される。

尚、上記実施例では重み係数比α＝Ｗ１／Ｗ２の値に応
じて値の変化する関数ｆ（α）１ｇ（α）を用いて重み
係数Ｗ２．バイアス分Ｃを求めたが、重み係数Ｗ２．バ
イアス分Ｃは必ずしも重み係数比αの関数である必要は
なく、骨部画像から軟部画像に至る一連の画像のコント
ラストや平均的な濃度（輝度）を固定しておいてもよい
。

また上記実施例では、重み係数Ｗ２を求めるに際し、（
１Ｂ）式に従って差分画像信号ＳＯを求め、その最大値
ＳＯａ＋ａｘと最小値５０ｍ１ｒ＋の差Ｓ０ｍａｘ−Ｓ
Ｏｗｉｎを用い（Ｉ７）式に従って重み係数Ｗ２を求め
たが、この差Ｓ　Ｏｉａ、ｘ　−Ｓ　Ｏ■ｉｎに代えて
差分画像信号ＳＯのばらつきの程度を指標する他の特微
量、例えば差分画像信号ＳＯの分散値等を用いてもよい
ものである。またこれと同様に、バイアス分Ｃを求める
に際し、差分画像信号ＳＯの最大値５ＯＩＩｌａｘと最
小値ＳＯａ＋ｉｏとの中央値（Ｓ０ｍａｘ　＋　Ｓ　０
ｍ１ｎ　）　／　２を用い（１８）式に従ってバイアス
分Ｃを求めたが、この中央値（Ｓ０ｍａｘ＋Ｓ　Ｏ層ｉ
ｎ　）　／　２に代えて差分画像信号ＳＯの平均的な値
を指標する他の特微量、例えば差分画像信号ＳＯの平均
値、メジアン値、出現頻度が最大となる差分画像信号Ｓ
Ｏの値等を用いてもよいものである。

さらに、上記実施例では重み係数Ｗ２．バイアス分Ｃを
求めるに際し、（１６）式に従って差分画像信号ＳＯを
求めたが、例えば胸部Ｘ線画像に限定した場合、重み係
数比αの各値に対する差分画像信号ＳＯのばらつきの程
度、平均的な値は統計的に求めることができるため、差
分画像信号ＳＯを求めることに代え、あらかじめ重み係
数比αの関数しての重み係数Ｗ２　（α）、バイアス分
Ｃ（α）を求めて記憶しておくようにしてもよい。

第５Ａ図、第５Ｂ図は、重み係数比αの関数としての重
み係数Ｗ２　（α）、バイアス分Ｃ（α）の−例を表わ
したグラフである。

これらのグラフをあらかじめ記憶しておき、重み係数比
αの各値に従って重み係数Ｗ２　（α）。

バイアス分Ｃ（α）が求められ、求められた重み係数Ｗ
２（α）、バイアス分Ｃ（α）を用いて、上記（１５）
式に従ってサブトラクション画像信号Ｓ３を求め、ＣＲ
Ｔデイスプレィ３２にこのサブトラクション画像信号Ｓ
３に基づく可視画像を再生表示するようにしてもよい。

尚、重み係数Ｗ２　（α）、バイアス分Ｃ（α）をあら
かじめ記憶しておくように構成する場合において、胸部
ＸＩＪｔ画像のみでなく、腹部画像等種々の被写体部分
を取扱う場合には、これら各被写体部分毎にそれぞれ適
切な重み係数Ｗ２　（α）。

バイアス分Ｃ（α）をあらかじめ記憶しておくようにし
てもよい。

また上記各実施例では、重み係数Ｗ２とバイアス分Ｃの
両者とも差分画像信号ＳＯを求めてこれに基づいて求め
るが、もしくは重み係数Ｗ２とバイアス分Ｃの両者とも
重み係数比αの関数としてあらかじめ求めて記憶してお
くかのいずれであるが、重み係数Ｗ２とバイアス分Ｃの
うちの一方は差分画像信号ＳＯに基づいて求め、他方は
重み係数比αの関数として求めておくようにしてもよい
。

さらに、上記各実施例は人体の胸部を被写体とした二つ
のＸ線画像に基づいて骨部画像から軟部画像に至る一連
の画像を表示する例であるが、本発明は胸部画像に限ら
れるものではなく、同一被写体のエネルギーサブトラク
ション画像のサブトラクション処理を行なって該処理後
の画像を表示する場合に広く適用できるものである。

（発明の効果） 以上詳細に説明したように、本発明のエネルギーサブト
ラクション画像表示装置は、差分画像データＳＯの、ば
らつきの程度および平均的な値を指標するそれぞれ第一
および第二の特徴量に基づいて重み係数およびバイアス
分を求めるが、もしくはこれらの重み係数およびバイア
ス分を重み係数比αの関数として記憶しておくようにし
たため、重みづけを順次変更したサブトラクション画像
を求める際に適切なコントラスト、濃度（輝度）が自動
的に設定され、観察適性の優れた可視画像が表示され、
しかも使い勝手のよい装置となる。

【図面の簡単な説明】

第１図は、Ｘ線撮影装置の概略図、 第２図は、Ｘ線画像読取装置と本発明のエネルギーサブ
トラクション画像表示装置の一例を内包した画像処理表
示装置の斜視図、 第３図は、画像処理表示装置内の内部メモリに記憶され
た第一および第二のＸ線画像を表わす２つの画像信号Ｓ
１，Ｓ２に基づいて、該画像処理表示装置内で行なわれ
る処理の流れを表わした図、第４Ａ図、第４Ｂ図は、重
み係数比αの関数ｆ（α）１ｇ（α）の−例を表わした
グラフ、第５Ａ図、第５Ｂ図は、重み係数比αの関数と
しての重み係数Ｗ２　（α）、バイアス分Ｃ（α）の−
例を表わしたグラフである。 １・・・Ｘ線撮影装置　　３・・・Ｘ線４・・・被写体 ５．７・・・蓄積性蛍光体シート ６・・・フィルタ　　　　１０・・・Ｘ線画像読取装置
１７・・・光ビーム　　　　１９・・・回転多面鏡２２
・・・輝尽発光光　　　２３・・・光ガイド２４・・・
フォトマルチプライヤ ３０・・・画像処理表示装置 ４１・・・第一のＸ線画像　４２・・・第二のＸ線画像
４３・・・骨部画像 ４４・・・骨部と軟部の双方か記録された画像４５・・
・軟部画像 ０ 第２図 ’３ｃ 第４Ａ図 第４Ｂ図 第５Ａ図 第５Ｂ図

Claims (7)

【特許請求の範囲】
1. （１）同一被写体の二つの放射線画像のそれぞれを表わ
   す二つの画像データをＳ１，Ｓ２、第一および第二の重
   み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ、サブト
   ラクション画像データをＳ３としたとき、次の第一式 Ｓ３＝Ｗ１・Ｓ１−Ｗ２・Ｓ２＋Ｃ に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
   て、順次変更されたサブトラクション画像データＳ３を
   求める画像演算部と、該順次変更されたサブトラクショ
   ン画像データＳ３に基づく順次変更されたサブトラクシ
   ョン画像を表示する表示部とを備えたエネルギーサブト
   ラクション画像表示装置において、 前記画像演算部が、 次の第二式 Ｓ０＝α・Ｓ１−Ｓ２ に基づいて差分画像データＳ０を求める差分画像演算手
   段と、 前記差分画像データＳ０のばらつきの程度を指標する第
   一の特徴量と前記差分画像データＳ０の平均的な値を指
   標する第二の特徴量とを求める特徴量演算手段と、 前記第一の特徴量に基づいて前記第二の重み係数Ｗ２を
   求める重み係数演算手段と、 前記第二の特徴量に基づいて前記バイアス分Ｃを求める
   バイアス演算手段と、 前記重み係数演算手段で求められた第二の重み係数Ｗ２
   と前記バイアス演算手段で求められたバイアス分Ｃとを
   用いて、前記第一式に従ってサブトラクション画像デー
   タＳ３を求めるサブトラクション画像演算手段とを備え
   たことを特徴とするエネルギーサブトラクション画像表
   示装置。
2. （２）各画素毎の前記差分画像データＳ０の最大値およ
   び最小値をそれぞれＳ０ｍａｘ，Ｓ０ｍｉｎとし、二つ
   の定数をＣ１，Ｃ２としたとき、 前記第一の特徴量および前記第二の特徴量が、それぞれ
   Ｓ０ｍａｘ−Ｓ０ｍｉｎ，（Ｓ０ｍａｘ＋Ｓ０ｍｉｎ）
   ／２であり、 前記重み係数演算手段が、前記第二の重み係数Ｗ２を次
   の第三式 Ｗ２＝Ｃ１／（Ｓ０ｍａｘ−Ｓ０ｍｉｎ） に従って求めるものであり、かつ 前記バイアス演算手段が前記バイアス分Ｃを次の第四式 Ｃ＝Ｃ２−（Ｓ０ｍａｘ＋Ｓ０ｍｉｎ）／２に従って求
   めるものであることを特徴とする請求項１記載のエネル
   ギーサブトラクション画像表示装置。
3. （３）前記重み係数演算手段が、前記第一の特徴量とと
   もに前記重み係数比αに基づいて前記第二の重み係数Ｗ
   ２を求めるものであり、 前記バイアス演算手段が前記第二の特徴量とともに前記
   重み係数比αに基づいて前記バイアス分Ｃを求めるもの
   であることを特徴とする請求項１記載のエネルギーサブ
   トラクション画像表示装置。
4. （４）各画素毎の前記差分画像データＳ０の最大値およ
   び最小値をそれぞれＳ０ｍａｘ，Ｓ０ｍｉｎ、前記重み
   係数比αを変数とした二つの関数をＦ（α），Ｇ（α）
   としたとき、 前記第一の特徴量および前記第二の特徴量がそれぞれＳ
   ０ｍａｘ−Ｓ０ｍｉｎ，（Ｓ０ｍａｘ＋Ｓ０ｍｉｎ）／
   ２であり、 前記重み係数演算手段が、前記第二の重み係数Ｗ２を次
   の第五式 Ｗ２＝Ｆ（α）／（Ｓ０ｍａｘ−Ｓ０ｍｉｎ）に従って
   求めるものであり、かつ 前記バイアス演算手段が、前記バイアス分Ｃを次の第六
   式 Ｃ＝Ｇ（α）−（Ｓ０ｍａｘ＋Ｓ０ｍｉｎ）／２に従っ
   て求めるものであることを特徴とする請求項３記載のエ
   ネルギーサブトラクション画像表示装置。
5. （５）同一被写体の二つの放射線画像のそれぞれを表わ
   す二つの画像データをＳ１，Ｓ２、第一および第二の重
   み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ、サブト
   ラクション画像データをＳ３としたとき、次の第七式 Ｓ３＝Ｗ１・Ｓ１−Ｗ２・Ｓ２＋Ｃ に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
   て、順次変更されたサブトラクション画像データＳ３を
   求める画像演算部と、該順次変更されたサブトラクショ
   ン画像データＳ３に基づく順次変更されたサブトラクシ
   ョン画像を表示する表示部とを備えたエネルギーサブト
   ラクション画像表示装置において、 前記画像演算部が、 前記重み係数比αを変数とした前記第二の重み係数Ｗ２
   を表わす関数Ｗ２（α）と、前記重み係数比αを変数と
   した前記バイアス分Ｃを表わす関数Ｃ（α）とを記憶し
   ておく記憶手段と、 前記記憶手段から読み出された第二の重み係数Ｗ２とバ
   イアス分Ｃとを用いて、前記第七式に従ってサブトラク
   ション画像データＳ３を求めるサブトラクション画像演
   算手段とを備えたことを特徴とするエネルギーサブトラ
   クション画像表示装置。
6. （６）同一被写体の二つの放射線画像のそれぞれを表わ
   す二つの画像データをＳ１，Ｓ２、第一および第二の重
   み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ、サブト
   ラクション画像データをＳ３としたとき、次の第八式 Ｓ３＝Ｗ１・Ｓ１−Ｗ２・Ｓ２＋Ｃ に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
   て、順次変更されたサブトラクション画像データＳ３を
   求める画像演算部と、該順次変更されたサブトラクショ
   ン画像データＳ３に基づく順次変更されたサブトラクシ
   ョン画像を表示する表示部とを備えたエネルギーサブト
   ラクション画像表示装置において、 前記画像演算部が、 次の第九式 Ｓ０＝α・Ｓ１−Ｓ２ に基づいて差分画像データＳ０を求める差分画像演算手
   段と、 前記差分画像データＳ０のばらつきの程度を指標する第
   一の特徴量を求める特徴量演算手段と、前記第一の特徴
   量に基づいて前記第二の重み係数Ｗ２を求める重み係数
   演算手段と、 前記重み係数比αを変数とした前記バイアス分Ｃを表わ
   す関数Ｃ（α）を記憶しておく記憶手段と、 前記重み係数演算手段で求められた第二の重み係数Ｗ２
   と前記記憶手段から読み出されたバイアス分Ｃとを用い
   て、前記第八式に従ってサブトラクション画像データＳ
   ３を求めるサブトラクション画像演算手段とを備えたこ
   とを特徴とするエネルギーサブトラクション画像表示装
   置。
7. （７）同一被写体の二つの放射線画像のそれぞれを表わ
   す二つの画像データをＳ１，Ｓ２、第一および第二の重
   み係数をそれぞれＷ１，Ｗ２、バイアス分をＣ、サブト
   ラクション画像データをＳ３としたとき、次の第十式 Ｓ３＝Ｗ１・Ｓ１−Ｗ２・Ｓ２＋Ｃ に基づき、かつ重み係数比α＝Ｗ１／Ｗ２を順次変更し
   て、順次変更されたサブトラクション画像データＳ３を
   求める画像演算部と、該順次変更されたサブトラクショ
   ン画像データＳ３に基づく順次変更されたサブトラクシ
   ョン画像を表示する表示部とを備えたエネルギーサブト
   ラクション画像表示装置において、 前記画像演算部が、 次の第十一式 Ｓ０＝α・Ｓ１−Ｓ２ に基づいて差分画像データＳ０を求める差分画像演算手
   段と、 前記差分画像データＳ０の平均的な値を指標する第二の
   特徴量を求める特徴量演算手段と、前記第二の特徴量に
   基づいて前記バイアス分Ｃを求めるバイアス演算手段と
   、 前記重み係数比αを変数とした前記第二の重み係数Ｗ２
   を表わす関数Ｗ２（α）を記憶しておく記憶手段と、 前記記憶手段から読み出された第二の重み係数Ｗ２と前
   記バイアス演算手段で求められたバイアス分Ｃとを用い
   て、前記第十式に従ってサブトラクション画像データＳ
   ３を求めるサブトラクション画像演算手段とを備えたこ
   とを特徴とするエネルギーサブトラクション画像表示装
   置。