(57)【特許請求の範囲】
【請求項１】
同一の被写体を撮影したエネルギー状態の互いに異なる複数の放射線画像を表す複数の放射線画像データについて、対応する画素ごとに前記放射線画像データの、所定のパラメータに基づいたサブトラクション処理を行なうことにより、前記放射線画像中の所望の組織の陰影を抽出または強調するエネルギーサブトラクション処理方法において、
前記サブトラクション処理が、下記演算式（１′）により規定されるエネルギーサブトラクション処理であり、
前記所定のパラメータのうち、前記エネルギーサブトラクション処理によって得られるエネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを除いた他のパラメータを、前記エネルギーサブトラクション処理の対象となる２つの画像データのうち一方または両者の加算である加算画像データを平滑化して得られた平滑化画像データに基づいて算出し、
前記濃度シフト量に関するパラメータを前記他のパラメータに基づいて算出することを特徴とするエネルギーサブトラクション処理方法。
Ｓproc＝Ｋａ・Ｈ−Ｋｂ・Ｌ＋Ｋｃ  （１′）
ただし、Ｓprocはサブトラクション画像を表すサブトラクション画像データ、Ｈはサブトラクション処理の対象となる２つの画像データのうち高圧側の画像データ、Ｌは低圧側の画像データ、Ｋａ，Ｋｂ，Ｋｃは前記所定のパラメータをそれぞれ意味する。
【請求項２】
前記エネルギーサブトラクション処理が下記演算式（１）により規定されるものであり、前記濃度シフト量に関するパラメータが下記演算式（２）により算出されるものであることを特徴とする請求項１記載のエネルギーサブトラクション処理方法。
Ｓproc＝Ｋａ（Ｌｕ）・Ｈ−Ｋｂ（Ｌｕ）・Ｌ＋Ｋｃ  （１）
Ｋｃ＝（１−Ｋａ（Ｌｕ）＋Ｋｂ（Ｌｕ））×Ｌｕ    （２）
ただし、Ｌｕは低圧側の画像データを平滑化して得られた平滑化低圧画像データ、Ｋａ（Ｌｕ），Ｋｂ（Ｌｕ）は前記他のパラメータ、Ｋｃは前記濃度シフト量に関するパラメータをそれぞれ意味する。
【請求項３】
前記他のパラメータを、前記放射線画像中の所望とする組織の陰影を最適な階調で再生し得る、前記各放射線画像に対する画像処理の処理条件の比に応じて変更することを特徴とする請求項１または２記載のエネルギーサブトラクション処理方法。
【請求項４】
前記処理条件が、前記放射線画像データの入出力比を表すスケールファクター値またはラチチュード値であることを特徴とする請求項３記載のエネルギーサブトラクション処理方法。
【請求項５】
同一の被写体を撮影したエネルギー状態の互いに異なる複数の放射線画像を表す複数の放射線画像データについて、対応する画素ごとに前記放射線画像データの、所定のパラメータに基づいたサブトラクション処理を行なうことにより、前記放射線画像中の所望の組織の陰影を抽出または強調するエネルギーサブトラクション処理装置において、
前記サブトラクション処理が、下記演算式（１′）により規定されるエネルギーサブトラクション処理であり、
前記所定のパラメータのうち、前記エネルギーサブトラクション処理によって得られるエネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを除いた他のパラメータを、前記エネルギーサブトラクション処理の対象となる２つの画像データのうち一方または両者の加算である加算画像データを平滑化して得られた平滑化画像データに基づいて算出する重み付けパラメータ算出手段と、
前記濃度シフト量に関するパラメータを前記他のパラメータに基づいて算出する濃度パラメータ算出手段とを備えたことを特徴とするエネルギーサブトラクション処理装置。
Ｓproc＝Ｋａ・Ｈ−Ｋｂ・Ｌ＋Ｋｃ  （１′）
ただし、Ｓprocはサブトラクション画像を表すサブトラクション画像データ、Ｈはサブトラクション処理の対象となる２つの画像データのうち高圧側の画像データ、Ｌは低圧側の画像データ、Ｋａ，Ｋｂ，Ｋｃは前記所定のパラメータをそれぞれ意味する。
【請求項６】
前記エネルギーサブトラクション処理が下記演算式（１）により規定されるものであり、前記濃度パラメータ算出手段が下記演算式（２）にしたがって前記濃度シフト量に関するパラメータを算出するものであることを特徴とする請求項５記載のエネルギーサブトラクション処理装置。
Ｓproc＝Ｋａ（Ｌｕ）・Ｈ−Ｋｂ（Ｌｕ）・Ｌ＋Ｋｃ  （１）
Ｋｃ＝（１−Ｋａ（Ｌｕ）＋Ｋｂ（Ｌｕ））×Ｌｕ    （２）
ただし、Ｌｕは低圧側の画像データを平滑化して得られた平滑化低圧画像データ、Ｋａ（Ｌｕ），Ｋｂ（Ｌｕ）は前記他のパラメータ、Ｋｃは前記濃度シフト量に関するパラメータをそれぞれ意味する。
【請求項７】
前記他のパラメータを、前記放射線画像中の所望とする組織の陰影を最適な階調で再生し得る、前記各放射線画像に対する画像処理の処理条件の比に応じて変更する重み付けパラメータ変更手段をさらに備えたことを特徴とする請求項５または６記載のエネルギーサブトラクション処理装置。
【請求項８】
前記パラメータ変更手段が、前記放射線画像データの入出力比を表すスケールファクター値またはラチチュード値を前記処理条件として、前記他のパラメータを変更するものであることを特徴とする請求項７記載のエネルギーサブトラクション処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明はエネルギーサブトラクション処理方法および装置に関し、詳細にはエネルギーサブトラクション処理のためのパラメータ変更の改良に関するものである。
【０００２】
【従来の技術】
従来より、放射線写真フイルムに記録された放射線画像を光電的に読み取って画像データを得、この画像データに適切な画像処理を施した後、画像を再生記録することが種々の分野で行われている。たとえば、後の画像処理に適合するように設計されたガンマ値の低いＸ線フィルムを用いてＸ線画像を記録し、このＸ線画像が記録されたフィルムからＸ線画像を光電的に読み取って電気信号を得、この画像信号をデジタル画像データに変換し、この画像データに画像処理を施した後、コピー写真等に可視像として再生することにより、コントラスト、シャープネス、粒状性等の画質性能の良好な再生画像を得ることが行われている（特公昭61-5193 号参照）。
【０００３】
また本願出願人により、放射線（Ｘ線，α線，β線，γ線，電子線，紫外線等）を照射すると、この放射線エネルギーの一部が蓄積され、その後可視光等の励起光を照射すると蓄積されたエネルギーに応じて輝尽発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、人体等の被写体の放射線画像情報を一旦シート状の蓄積性蛍光体に記録し、この蓄積性蛍光体シートをレーザー光等の励起光で走査して輝尽発光光を生ぜしめ、得られた輝尽発光光を光電的に読み取って画像データを得、この画像データに基づき被写体の放射線画像を写真感光材料等の記録材料、ＣＲＴ等に可視像として出力させる放射線画像記録再生システムがすでに提案されている（特開昭55-12429号，同56-11395号，同55-163472 号，同56-104645 号，同55- 116340号等）。このシステムは、従来の銀塩写真を用いる放射線写真システムと比較して極めて広い放射線露出域にわたって画像を記録しうるという実用的な利点を有している。
【０００４】
上記システムにおいて、蓄積性蛍光体シートに照射された放射線の線量等に応じて最適な読取条件で読み取って画像信号を得る前に、予め低レベルの光ビームにより蓄積性蛍光体シートを走査してこのシートに記録された放射線画像の概略を読み取る先読みを行い、この先読みにより得られた先読画像信号を分析し、その後上記シートに高レベルの光ビームを照射して走査し、この放射線画像に最適な読取条件で読み取って画像信号を得る本読みを行うように構成されたシステムが知られている。
【０００５】
ここで読取条件とは、読取りにおける輝尽発光光の光量と読取装置の出力との関係に影響を与える各種の条件を総称するものであり、例えば入出力の関係を定める読取ゲイン、スケールファクタ、あるいは読取りにおける励起光のパワー等を意味するものである。
【０００６】
また、この先読みを行わずに読取条件を設定するための方法が本願出願人により提案されている（特開昭55-48672号，同55-50180号，同56-11348号等）。この方法は、蓄積性蛍光体シートに放射線を照射したときにシートから発せられる瞬時発光光を検光検出器により検出し、この瞬時発光光から放射線画像の特性あるいは蓄積性蛍光体シートに蓄積記録された放射線の蓄積量等の情報を得て、読取条件を決定する方法である。
【０００７】
さらに、この先読みを行うシステムか先読みを行わないシステムかによらず、得られた画像信号（先読画像信号を含む）を分析し、画像信号に画像処理を施す際の最適な画像処理条件を決定するようにしたシステムもある。ここで画像処理条件とは、画像信号に基づく再生画像の階調や感度等に影響を及ぼす処理を該画像信号に施す際の各種の条件を総称するものであり、前述した先読みを行なわないシステムにおいては前述した読取条件である読取ゲイン、スケールファクタをも含むものである。この画像信号に基づいて最適な画像処理条件を決定する方法は、蓄積性蛍光体シートを用いるシステムに限られず、たとえば従来のＸ線フィルム等の記録シートに記録された放射線画像から画像信号を得るシステムにも適用されている（以下、本明細書中、このような読取条件および／または画像処理条件を決定するシステムをＥＤＲ（Exposure Data Recognizer）処理システムまたはＥＤＲ処理手段ということがある）。
【０００８】
また上記のようにＸ線フイルムや蓄積性蛍光体シート等を用いるシステムにおいて、記録された複数の放射線画像を読み取って複数の画像データを得た後、これらの画像データに基づいて上記放射線画像のサブトラクション処理を施すことがある。
【０００９】
ここで、放射線画像のサブトラクション処理とは、互いに異なった条件で撮影された複数の放射線画像の差に対応する画像を得る処理をいい、具体的にはこれら複数の放射線画像を所定のサンプリング間隔で読み取って各放射線画像に対応する複数のデジタルの画像データを得、これら複数の画像データの各対応するサンプリング点（画素）毎に減算処理を施すことにより、放射線画像中の特定の被写体部分（以下、組織または構造物等の陰影とも称する）のみを強調または抽出した放射線画像を得る処理をいう。
【００１０】
このサブトラクション処理には基本的には次の二つがある。すなわち、造影剤の注入により被写体の特定の部分（たとえば人体を被写体としたときの血管等の陰影）が強調された放射線画像から造影剤が注入されていない放射線画像を引き算（サブトラクト）することによって被写体の特定の部分（血管等の陰影）を抽出するいわゆる時間サブトラクションと、被写体の特定の部分が互いに異なるエネルギーを有する放射線に対して異なる放射線吸収率を有することを利用して、同一の被写体に対して互いに異なるエネルギーを有する放射線を照射してこれら互いに異なるエネルギーを有する各放射線による複数の放射線画像を得、これら複数の放射線画像を適当に重み付けしてその差を演算すること（下記式（３）参照）によって被写体の特定部分を抽出するいわゆるエネルギーサブトラクションとがある。本出願人も蓄積性蛍光体シートを用いたエネルギーサブトラクションについて提案している（特開昭59-83486号，特開昭60-225541号）。
【００１１】
Ｓproc＝Ｋａ・Ｈ−Ｋｂ・Ｌ＋Ｋｃ              （３）
ただし、Ｓprocはサブトラクション処理により得られるサブトラクション画像データ、Ｋａ，Ｋｂは重み付け係数、Ｋｃはバイアス成分（以下、Ｋａ，Ｋｂ，Ｋｃをまとめてサブトラクション処理のパラメータという）、Ｈはいわゆる高圧側の画像データ、Ｌはいわゆる低圧側の画像データをそれぞれ意味する。
【００１２】
【発明が解決しようとする課題】
ところで、上記エネルギーサブトラクション処理に関し、被写体の放射線撮影を行なう際、所定のエネルギー分布を有する放射線を被写体に照射すると、例えば被写体の体厚に応じて放射線透過率が異なり、しかも低エネルギーの放射線ほど透過率が低いため、放射線が被写体を透過するにしたがって放射線のエネルギー分布が全体として高エネルギー側に偏るいわゆるビームハードニングの現象が生じることが知られている。
【００１３】
例えば骨部の骨塩定量をエネルギーサブトラクション処理によって行なう場合、骨部の厚さが同一であってもこの骨部を取り囲む軟部組織の厚さが厚い場合はビームハードニングによる影響が大きいため骨部濃度が薄くなり、一方、軟部組織の厚さが薄い場合はビームハードニングによる影響が小さいため骨部濃度が濃くなる。
【００１４】
このようにビームハードニングの影響の大小によって、抽出した組織の濃度が変動したのでは的確な診断に悪影響を与えることになる。特に過去に得られたエネルギーサブトラクション画像と現在得られたエネルギーサブトラクション画像とを比較する場合に、両画像間でビームハードニングの影響度合に差があれば、現在も病巣部が存在するにも拘らず、治癒したとの誤った判断がなされる虞れもある。
【００１５】
またエネルギーサブトラクション画像は、放射線の被写体内部における散乱の影響を受けやすく、散乱の度合に応じてエネルギーサブトラクション画像の画質が変化する。したがってエネルギーサブトラクション画像の画質を向上させるためには放射線の被写体内部における散乱の影響も考慮する必要がある。
【００１６】
そこで本出願人は、ビームハードニングを低減すべく、被写体の体厚やこの体厚と相関を有する所定の画像処理条件に応じて、上記サブトラクション処理の各パラメータＫａ，Ｋｂ，Ｋｃをそれぞれ独立して変更するエネルギーサブトラクション処理方法を提案している（特願平8-278236号（特開平10-118056号））。
【００１７】
しかし、上述したパラメータのうちバイアス成分であるパラメータＫｃは、エネルギーサブトラクション画像全体の濃度シフト量を規定するパラメータであるため、他のパラメータから独立して変更すると、エネルギーサブトラクション画像全体の濃度がビームハードニングの大小によって変化してしまうおそれがある。
【００１８】
本発明は上記事情に鑑みなされたものであって、放射線のビームハードニング、散乱の影響を低減しつつ、得られるエネルギーサブトラクション画像の濃度を安定させることができるエネルギーサブトラクション処理方法および装置を提供することを目的とするものである。
【００１９】
【課題を解決するための手段】
本発明のエネルギーサブトラクション処理方法および装置は、所定のパラメータに基づいたサブトラクション処理を行なうときの、エネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを、他のパラメータに基づいて算出することにより、放射線のビームハードニングや散乱の影響を低減しつつ、エネルギーサブトラクション画像の濃度を安定させるものである。
【００２０】
すなわち本発明のエネルギーサブトラクション処理方法は、同一の被写体を撮影したエネルギー状態の互いに異なる複数の放射線画像を表す複数の放射線画像データについて、対応する画素ごとに前記放射線画像データの、所定のパラメータに基づいたサブトラクション処理を行なうことにより、前記放射線画像中の所望の組織の陰影を抽出または強調するエネルギーサブトラクション処理方法において、
前記所定のパラメータのうち、エネルギーサブトラクション処理によって得られるエネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを、前記所定のパラメータのうち前記濃度シフト量に関するパラメータを除いた他のパラメータに基づいて算出することを特徴とするものである。
【００２１】
ここで、所定のパラメータに基づいたエネルギーサブトラクション処理とは、例えば下記式に示すものを意味する。
【００２２】
Ｓproc＝Ｋａ・Ｈ−Ｋｂ・Ｌ＋Ｋｃ                        （１′）
ただし、Ｓprocはサブトラクション画像を表すサブトラクション画像データ、Ｈはサブトラクション処理の対象となる２つの画像データのうち高圧側の画像データ、Ｌは低圧側の画像データ、Ｋａ，Ｋｂは上記他のパラメータ（重み付けのパラメータ）、Ｋｃは濃度シフト量に関するパラメータをそれぞれ意味する。
【００２３】
高圧側の画像データ、低圧側の画像データとは、それぞれ、他方に対して相対的に高エネルギー成分が強調された（または、低エネルギー成分が抑制された）画像データ、低エネルギー成分が強調された（または高エネルギー成分が抑制された）画像データを意味するものであり、例えば、エネルギー分離フィルターを介して２枚の蓄積性蛍光体シートを重ねるなどして、１回の放射線照射でエネルギー分布が互いに異なるこれら２つの画像データをそれぞれの蓄積性蛍光体シートに略同時に記録する、いわゆる１ショットエネルギーサブトラクション撮影においては、被写体から遠い側に配された蓄積性蛍光体シートからは高圧側の画像データが得られ、被写体に近い側に配された蓄積性蛍光体シートからは低圧側の画像データが得られる。
【００２４】
ここで上記式（１′）に示したエネルギーサブトラクション処理のパラメータのうち上記他のパラメータＫａ，Ｋｂは、エネルギーサブトラクション処理の対象となる２つの画像データＨ，Ｌのうち一方（ＨまたはＬ）または両者を加算して得られた加算画像データ（αＨ＋βＬ）を、平滑化して得られた平滑化低圧画像データＸｕに基づいて算出されたもの（Ｋａ（Ｘｕ），Ｋｂ（Ｘｕ））とするのが好ましい。特にパラメータＫａ，Ｋｂを、低圧側の画像データＬを平滑化して得られた平滑化低圧画像データＬｕに依存するように算出されたもの（Ｋａ（Ｌｕ），Ｋｂ（Ｌｕ））とするのが好ましく、式（１′）は下記式（１）のように示される。
【００２５】
Ｓproc＝Ｋａ（Ｌｕ）・Ｈ−Ｋｂ（Ｌｕ）・Ｌ＋Ｋｃ  （１）
ただし、Ｓprocはサブトラクション画像を表すサブトラクション画像データ、Ｈはサブトラクション処理の対象となる２つの画像データのうち高圧側の画像データ、Ｌは低圧側の画像データ、Ｋａ（Ｌｕ），Ｋｂ（Ｌｕ）は平滑化低圧画像データについての重み付けのパラメータ、Ｋｃは濃度シフト量に関するパラメータをそれぞれ意味する。
【００２６】
濃度シフト量に関するパラメータを除いた他のパラメータに基づいて算出する場合の演算としては、例えば下記式（２）を適用するのが好ましい。
【００２７】
Ｋｃ＝（１−Ｋａ（Ｌｕ）＋Ｋｂ（Ｌｕ））×Ｌｕ    （２）
ただし、Ｌｕは平滑化低圧画像データを意味する。
【００２８】
上記式（２）は、平滑化低圧画像データＬｕに基づいて式（１′）に示したパラメータＫａ，Ｋｂを求め（すなわち、Ｋａ（Ｌｕ），Ｋｂ（Ｌｕ）を求める）、サブトラクション画像Ｓprocとして軟部組織を適用した場合（軟部組織を強調・抽出することを意味する）に、式（１）において、Ｓproc＝Ｌｕ＝Ｈ＝Ｌなる近似を適用することにより導出される。なお、この近似を適用する場合は、平滑化する以前の低圧側画像データＬを、所望とする組織の陰影を最適な階調で再生し得る画像処理条件にしたがった規格化処理を予め施しておくことが好ましく、高圧側画像データについても同様に規格化処理を施しておくことが好ましい。また平滑化して得られた平滑化低圧画像データＬｕについてさらに規格化処理を施してもよい。高圧側画像データＨ、低圧側画像データＬ、および平滑化低圧画像データＬｕの濃度およびコントラストが安定し、近似するのに適切なものとなるからである。
【００２９】
パラメータＫａ（Ｌｕ），Ｋｂ（Ｌｕ）は、そのような低圧画像についての画像処理条件と高圧画像についての画像処理条件との比に応じて変更するのも好適であり、エネルギーサブトラクション処理の対象となる２つの画像の例えばコントラストに関する条件の比、具体的には上述したＥＤＲ処理等の規格化処理において決定される読取条件および／または画像処理条件の比などである。このような読取条件または画像処理条件の比としては、その処理による放射線画像データの入出力比を表すスケールファクター値（Ｇｐ値）の比またはラチチュード値（Ｌ値＝４／Ｇｐ）の比を適用するのが好ましい。
【００３０】
具体的には、スケールファクタ（Ｇｐ１；高圧側，Ｇｐ２；低圧側）の比（＝Ｇｐ２／Ｇｐ１）が大きくなるにしたがって、パラメータＫａ（Ｌｕ），Ｋｂ（Ｌｕ）を大きくすればよく、またはラチチュード（Ｌ１；高圧側，Ｌ２；低圧側）の比（＝Ｌ２／Ｌ１）が小さくなるにしたがって、パラメータＫａ（Ｌｕ），Ｋｂ（Ｌｕ）を大きくすればよい。
【００３１】
なお上述した本発明のエネルギーサブトラクション処理方法においては、エネルギーサブトラクション処理の各パラメータを、サブトラクション処理の対象画素の全てにおいて同一の値を採る必要はなく、対象画素ごとに変化するものとするのが好ましい。
【００３２】
本発明のエネルギーサブトラクション処理装置は、上記本発明のエネルギーサブトラクション処理方法を実施するための装置であって、同一の被写体を撮影したエネルギー状態の互いに異なる複数の放射線画像を表す複数の放射線画像データについて、対応する画素ごとに前記放射線画像データの、所定のパラメータに基づいたサブトラクション処理を行なうことにより、前記放射線画像中の所望の組織の陰影を抽出または強調するエネルギーサブトラクション処理装置において、
前記所定のパラメータのうち、エネルギーサブトラクション処理によって得られるエネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを、前記所定のパラメータのうち前記濃度シフト量に関するパラメータを除いた他のパラメータに基づいて算出する濃度パラメータ算出手段を備えたことを特徴とするものである。
【００３３】
ここで、エネルギーサブトラクション処理の対象となる２つの画像データのうち一方を平滑化して得られた平滑化画像データに基づいて、または２つの画像データの加算である加算画像データを平滑化して得られた平滑化画像データに基づいて、上記他のパラメータを算出する重み付けパラメータ算出手段を備えたものとするのが好ましい。特に、エネルギーサブトラクション処理を上記演算式（１）により規定されるものとし、濃度パラメータ算出手段を上記演算式（２）にしたがって濃度シフト量に関するパラメータを算出するものとするのが好ましい。
【００３４】
また、上記他のパラメータを、放射線画像中の所望とする組織の陰影を最適な階調で再生し得る、各放射線画像に対する画像処理の処理条件の比に応じて変更する重み付けパラメータ変更手段をさらに備えたものとするのが好ましく、このパラメータ変更手段を、放射線画像データの入出力比を表すスケールファクター値またはラチチュード値を上記画像処理条件として適用して、他のパラメータを変更するものとするのがより好ましい。
【００３５】
【発明の効果】
本発明のエネルギーサブトラクション処理方法、装置によれば、所定のパラメータに基づいたサブトラクション処理を行なうときの、エネルギーサブトラクション画像全体の濃度シフト量に関するパラメータを、他のパラメータに基づいて算出することにより、放射線のビームハードニングや散乱の影響を低減しつつ、エネルギーサブトラクション画像の濃度を安定させることができる。
【００３６】
【発明の実施の形態】
以下、本発明のエネルギーサブトラクション処理方法、装置の具体的な実施の形態について、図面を用いて説明する。
【００３７】
図１は本発明のエネルギーサブトラクション処理装置の一実施形態の構成を示すブロック図、図２は放射線撮影装置の概略図、図３は放射線画像読取装置と、画像処理装置および図１に示したエネルギーサブトラクション処理装置を内包する画像処理表示装置を示す斜視図である。
【００３８】
図示のエネルギーサブトラクション処理装置１０は、エネルギー状態が相異なる２つの放射線画像Ｐ１，Ｐ２（Ｐ１；高エネルギー成分が相対的に強調された画像（以下、高圧画像という），Ｐ２；低エネルギー成分が相対的に強調された画像（以下、低圧画像という））をそれぞれ読み取って得られた２つの放射線画像データＳ１，Ｓ２（Ｓ１；高圧画像を読み取って得られた高圧画像データ，Ｓ２；低圧画像を読み取って得られた低圧画像データ）の規格化処理された画像データＳ１′，Ｓ２′（Ｓ１′；高圧画像データを規格化処理して得られた規格化高圧画像データ，Ｓ２′；低圧画像データを規格化処理して得られた規格化低圧画像データ）を記憶する内部メモリ４と、外部の画像処理装置２０から入力された、低圧画像データＳ２を規格化した後（規格化低圧画像データＳ２′）に平滑化処理して得られた平滑化低圧画像データＳLu′に基づいて、エネルギーサブトラクションのパラメータＫａ(Lu′)，Ｋｂ(Lu′)を算出する重み付けパラメータ算出手段１と、同じく画像処理装置２０から入力された高圧画像データＳ１についてのスケールファクターＧｐ１と低圧側画像データＳ２についてのスケールファクターＧｐ２との比（＝Ｇｐ２／Ｇｐ１）が大きくなるにしたがって、重み付けパラメータＫａ(Lu′)，Ｋｂ(Lu′)をそれぞれＫａ(Lu′)′，Ｋｂ(Lu′)′に変更する重み付けパラメータ変更手段２と、変更された重み付けパラメータＫａ(Lu′)′，Ｋｂ(Lu′)′および平滑化低圧画像データＳLu′に基づいて濃度シフトパラメータＫｃを算出する濃度パラメータ算出手段３と、パラメータＫａ(Lu′)′，Ｋｂ(Lu′)′およびＫｃに基づいて、下記式（１″）にしたがって内部メモリ４に記憶された２つの規格化画像データＳ１′，Ｓ２′に対してサブトラクション処理を行なうエネルギーサブトラクション処理手段５とを備えた構成である。
【００３９】
Ｓproc＝Ｋａ(Lu′)′・Ｓ１′−Ｋｂ(Lu′)′・Ｓ２′＋Ｋｃ  （１″）
ただし、Ｓprocはサブトラクション処理により得られるサブトラクション画像データを意味する。
【００４０】
ここで外部の画像処理装置２０は、入力された各画像データＳ１，Ｓ２についてそれぞれ、画像データのレベルごとの頻度の分布を示すヒストグラムを作成し、このヒストグラムの形状に基づいて所望とする画像部分に対応する画像データを認識し、その画像部分を最適な階調の可視像として再生し得るように画像データＳ１，Ｓ２に対して規格化処理（ＥＤＲ処理）を行うとともに、規格化して得られた規格化低圧画像データＳ２′に対してさらに平滑化処理を施して平滑化低圧画像データＳLu′を算出するものであり、具体的には、高圧画像データＳ１についてその頻度分布を表すヒストグラムを作成し、このヒストグラムに基づいてスケールファクターＧｐ１，ゲインＳｋ１を求め、これらスケールファクターＧｐ１，ゲインＳｋ１に基づいて高圧画像データＳ１に対して規格化処理を施して規格化高圧画像データＳ１′を得、低圧画像データＳ２についてヒストグラムを作成し、このヒストグラムに基づいてスケールファクターＧｐ２，ゲインＳｋ２を求め、これらスケールファクターＧｐ２，ゲインＳｋ２に基づいて低圧画像データＳ２に対して規格化処理を施して規格化低圧画像データＳ２′を得、規格化低圧画像データＳ２′を平滑化して平滑化低圧画像データＳLu′を得る処理をなす。
【００４１】
濃度パラメータ算出手段３は詳しくは、重み付けパラメータＫａ(Lu′)′，Ｋｂ(Lu′)′および平滑化低圧画像データＳLu′に基づき、下記式（２′）にしたがって式（１″）における濃度シフトパラメータＫｃを算出する作用をなす。
【００４２】
Ｋｃ＝（１−Ｋａ(Lu′)′＋Ｋｂ(Lu′)′）×ＳLu′    （２′）
放射線撮影装置３０の放射線管３１から発せられた放射線Ｘが被写体（人体の胸部）４０に照射される。被写体４０を透過した放射線Ｘａは、被写体４０に近い側に配された第２の蓄積性蛍光体シート５２に照射され、第２のシート５２に被写体４０の放射線画像Ｐ２が蓄積記録される。なおこの放射線画像Ｐ２は、後述する高圧画像Ｐ１に対して相対的に低エネルギー成分を多く含むものであるため、低圧画像と称する。
【００４３】
第２のシート５２を透過した放射線Ｘｂはさらに放射線の低エネルギー成分を抑制するフィルタ５５に到達し、放射線Ｘｂのエネルギーのうち低エネルギー成分が抑制された成分（相対的に高エネルギー成分が強調された成分）Ｘｃがフィルタ５５を透過し、この相対的に高エネルギー成分が強調された放射線Ｘｃが第１の蓄積性蛍光体シート５１に到達し、これにより第１のシート５１に被写体４０の高圧画像Ｐ１が蓄積記録される。
【００４４】
被写体４０には、サブトラクション処理を行なうにあたって２つの放射線画像Ｐ１，Ｐ２の位置合わせを行なうための基準となる２つのマーク４１が付されている。なお上記放射線撮影装置３０は１回の撮影（１ショット）で２つのシート５１，５２にそれぞれ放射線画像Ｐ１，Ｐ２を略同時に蓄積記録するものであるが、中間層にエネルギー分離フィルター５５を有し両面に蓄積性蛍光体層を備えた単一の蓄積性蛍光体シートを用いて撮影を行なうものであってもよいし、時間的に相前後した２つのタイミングでそれぞれ１枚ずつのシートに撮影を行なうもの（２ショット）であってもよい。
【００４５】
図２に示した放射線撮影装置３０で撮影が行なわれた後、第１および第２の蓄積性蛍光体シート５１，５２が一枚ずつ放射線画像読取装置６０の所定位置にセットされる。ここでは、第１の蓄積性蛍光体シート５１に蓄積記録された高圧画像Ｐ１の読取りの場合について説明する。
【００４６】
所定位置にセットされた、高圧画像Ｐ１が蓄積記録された蓄積性蛍光体シート５１は、図示しない駆動手段により駆動されるエンドレスベルト等のシート搬送手段６１により、矢印Ｙ方向に搬送（副走査）される。一方、レーザ光源６２から発せられた光ビームＬはモータ６３により駆動され矢印Ｒ方向に高速回転する回転多面鏡６４によって反射偏向され、ｆθレンズ等の集束レンズ６５を通過した後、ミラー６６により光路を変えて蓄積性蛍光体シート５１に入射し、副走査の方向（矢印Ｙ方向）に略直交する矢印Ｘ方向に主走査する。蓄積性蛍光体シート５１の、光ビームＬが照射された箇所からは、蓄積記録されている高圧画像Ｐ１に応じた光量の輝尽発光光Ｍ１が発せられ、この輝尽発光光Ｍ１は光ガイド６７によって導かれ、フォトマルチプライヤ（光電子増倍管）６８によって光電的に検出される。
【００４７】
フォトマルチプライヤ６８から出力されたアナログ信号Ａ１は、ログアンプ６９で対数的に増幅された後、Ａ／Ｄ変換器７０に入力され、サンプリングされてデジタルの高圧画像データＳ１が得られる。この高圧画像データＳ１は画像処理装置２０内の内部メモリ（図示せず）に一旦記憶される。
【００４８】
画像処理表示装置８０は、画像処理装置２０およびサブトラクション処理装置１０を内包するものであり、画像信号に対して各種の画像処理を施すものであり、種々の指示等を入力するキーボード８１、指示のための補助情報や画像データに基づく可視画像を表示するＣＲＴディスプレイ８２、並びに画像処理装置２０、サブトラクション処理装置１０およびその他の演算機能を備えた本体部８３が備えられている。
【００４９】
次に上記と同様にして、第２の蓄積性蛍光体シート５２に蓄積記録された低圧画像Ｐ２を表わす低圧画像データＳ２が得られ、この低圧画像データＳ２も画像処理装置２０内の内部メモリに一旦記憶される。
【００５０】
画像処理装置２０は内部メモリに記憶された各画像データＳ１，Ｓ２についてそれぞれ規格化処理を施して規格化画像データＳ１′，Ｓ２′を得、その過程でスケールファクターＧｐ１，Ｇｐ２を得る。さらに画像処理装置２０は、規格化低圧画像データＳ２′に対して平滑化処理を施して平滑化画像データＳLu′を算出し、この平滑化低圧画像データＳLu′を重み付けパラメータ算出手段１に入力する。また、規格化画像データＳ１′，Ｓ２′は、サブトラクション処理装置１０内の内部メモリ３に一旦記憶され、各スケールファクターＧｐ１，Ｇｐ２は重み付けパラメータ変更手段２に入力される。
【００５１】
重み付けパラメータ算出手段１は、入力された平滑化低圧画像データＳLu′に基づいてエネルギーサブトラクションのパラメータＫａ(Lu′)およびＫｂ(Lu′)を算出する。パラメータ変更手段２は入力されたスケールファクターＧｐ１，Ｇｐ２の比（＝Ｇｐ２／Ｇｐ１）を算出し、重み付けパラメータ算出手段１から入力されたＫａ(Lu′)，Ｋｂ(Lu′)を、それぞれＫａ(Lu′)′，Ｋｂ(Lu′)′に変更する。そして変更された重み付けパラメータＫａ(Lu′)′，Ｋｂ(Lu′)′は、濃度パラメータ算出手段３およびサブトラクション処理手段５に入力される。
【００５２】
濃度パラメータ算出手段３は、入力された重み付けパラメータＫａ(Lu′)′，Ｋｂ(Lu′)′および規格化された平滑化低圧画像データＳLu′に基づいて、式（２′）にしたがって濃度シフトパラメータＫｃを算出し、得られた濃度シフトパラメータＫｃをサブトラクション処理手段５に入力する。
【００５３】
サブトラクション処理手段５には、内部メモリ３から規格化高圧画像データＳ１′および規格化低圧画像データＳ２′が入力され、重み付けパラメータ変更手段２から入力された新たなパラメータＫａ(Lu′)′およびＫｂ(Lu′)′、並びに濃度パラメータ算出手段３から入力された濃度シフトパラメータＫｃを用いた上記式（１″）にしたがって、サブトラクション画像データＳprocを算出する。
【００５４】
このように算出されたサブトラクション画像データＳprocは、濃度シフト量に関するパラメータＫｃを、他のパラメータＫａ(Lu′)′，Ｋｂ(Lu′)′に基づいて算出することにより、放射線のビームハードニングや散乱の影響を低減しつつ、エネルギーサブトラクション画像Ｐproc（エネルギーサブトラクション画像データＳprocが表すエネルギーサブトラクション画像）の濃度を安定させることができる。
【００５５】
そして得られたサブトラクション画像データＳprocは、ＣＲＴ８２に入力され、ビームハードニングの影響が十分に低減されるとともに安定的に所定の濃度範囲になる、診断性能の高い可視像であるサブトラクション画像Ｐprocとして再生され、画像読影および診断に供される。
【００５６】
なお、本実施形態のエネルギーサブトラクション処理装置は、パラメータ変更手段２が、スケールファクターの比（＝Ｇｐ２／Ｇｐ１）に応じて、パラメータＫａ(Lu′)，Ｋｂ(Lu′)を変更するものとしたが、本発明のエネルギーサブトラクション処理方法、装置は、この態様に限定されるものではなく、高圧画像データＳ１の画像処理条件（特にコントラストに係るもの）と低圧画像データＳ２の画像処理条件（特にコントラストに係るもの）との比に応じてパラメータを変更するものであれば、如何なるものであってもよく、例えば、ラチチュード値（Ｌ値＝４／Ｇｐ）の比に応じてパラメータを変更するものなどを適用することもできる。ラチチュード値の比に応じてパラメータを変更する場合は、このラチチュード値の比（＝Ｌ２／Ｌ１：Ｌ１；高圧側画像データＳ１についてのラチチュード，Ｌ２；低圧側画像データＳ２についてのラチチュード）が小さくなるにしたがって、パラメータＫａ(Lu′)，Ｋｂ(Lu′)が大きくなるように変更すればよい。
【００５７】
なお、処理対象の複数の放射線画像データを得る方式は、上記のような蓄積性蛍光体シートを走査する方式に限られるものではなく、固体検出器を用いる方式などでもよい。
【図面の簡単な説明】
【図１】本発明のエネルギーサブトラクション処理装置の一実施形態の構成を示すブロック図
【図２】放射線撮影装置の概略図
【図３】放射線画像読取装置と画像処理表示装置を示す図
【符号の説明】
１    重み付けパラメータ算出手段
２    重み付けパラメータ変更手段
３    濃度パラメータ算出手段
４    内部メモリ
５    サブトラクション処理手段
10    エネルギーサブトラクション処理装置
20    画像処理装置
30    放射線撮影装置
40    被写体
51,52 蓄積性蛍光体シート
82    ＣＲＴ
【図１】
0003775588.tif 000002
【図２】
0003775588.tif 000003
【図３】
0003775588.tif 000004