JP4250372B2

JP4250372B2 - 放射線撮影装置、放射線撮影方法、放射線画像用システム、プログラム、及びコンピュータ可読記憶媒体

Info

Publication number: JP4250372B2
Application number: JP2002101208A
Authority: JP
Inventors: 仁司井上
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2002-04-03
Filing date: 2002-04-03
Publication date: 2009-04-08
Anticipated expiration: 2022-04-03
Also published as: JP2003298939A

Description

【０００１】
【発明の属する技術分野】
本発明は医療分野における放射線撮影装置、放射線撮影方法、放射線画像用システム、プログラム、及びコンピュータ可読記憶媒体に関する。
【０００２】
【従来の技術】
Ｘ線に代表される放射線の物質透過能力を用いて、その透過強度分布を画像化する技術は、近代医療技術発展の基本となるものである。Ｘ線発見以来、その強度分布の画像化は、Ｘ線強度分布を蛍光体により可視光に変換した後、銀塩フィルムで潜像を作り現像するという方法が取られてきた。近年、Ｘ線画像をデジタル化する際は輝尽性蛍光体を用い、Ｘ線照射による輝尽性蛍光体上の蓄積エネルギ分布としての潜像をレーザ光で励起して読み出し、デジタル画像化する、いわゆるイメージングプレートを用いる方法が一般化してきた。さらに、半導体技術の進歩により人体の大きさをカバーできる大判の固体撮像素子、いわゆるフラットパネルディテクタも開発され、潜像をつくることなく直接Ｘ線画像をデジタル化し、効率のよい診断が行えるようになって来た。
【０００３】
また一方、光電子増倍管（イメージインテンシファイア）に代表される高感度の撮像素子により微弱なＸ線による蛍光を画像化し、人体内部の動態を観察することも可能であり、一般に用いられてきている。そして、最新のフラットパネルディテクタはそのイメージインテンシファイアにも匹敵する感度を持ち、人体の広範囲における動態を撮影することが可能になってきている。
【０００４】
医療用のＸ線撮影で最も有効であるのは人体の胸部撮影である。腹部を含む胸部の広範囲を撮影すれば、肺疾患を含む多くの疾病の発見に役立つため、通常の健康診断では胸部Ｘ線撮影は不可欠なものになっている。また、近年健康診断のために撮影された膨大な量の胸部Ｘ線画像を効率よく診断するため、胸部デジタルＸ線画像に対し、計算機を用いて画像解析を行い、医師の初期診断を補助するいわゆる計算機支援診断（Ｃｏｍｐｕｔｅｒ−ＡｉｄｅｄＤｉａｇｎｏｓｉｓ，ＣＡＤ）も実用化しつつある。
【０００５】
そして、近年、人体胸部全体を撮影できるような大きさのフラットパネルディテクタを用いて、人体の呼吸動態を観察することにより、さらに診断精度の向上を図ることができるようになりつつあり、また、その膨大な画像データを効率よく処理する動態ＣＡＤ技術も開発されつつある。
【０００６】
【発明が解決しようとする課題】
呼吸動態のＣＡＤは、主に肺野部分の呼吸による動態を観察するものであるが、一般撮影で撮影されるのは肺野部分だけでなく、それを覆う肋骨などが同時に写り込む。肋骨などの呼吸に伴う動きは肺野部分に比べれば少ないものであるが、当然呼吸とともに移動する。この肋骨の移動は目的とする肺野の動きとは必ずしも一致しないため、計算機により肺野（肺門理）の移動過程を追跡する場合に、肋骨の動きはいわゆるノイズとなる。
【０００７】
そこで、本発明は、エネルギサブトラクション処理の適用により、安定した肺野部動態の把握が可能な放射線撮影装置、放射線撮影方法、放射線画像用システム、プログラム、及びコンピュータ可読記憶媒体を提供することを目的とする。
【０００８】
【課題を解決するための手段】
上記目的を達成するための放射線撮影装置は、第１のエネルギのＸ線で撮像された複数の画像データと、第１のエネルギと異なる第２のエネルギで撮像された複数の画像データとを、対応付ける第１の手段と、
前記第１の手段で対応付けられた画像データ間でエネルギサブトラクション演算を行う第２の手段と、
を備えることを特徴とする。
【０００９】
上記目的を達成するためのプログラムは、コンピュータを、
第１のエネルギのＸ線で撮像された複数の画像データと、第１のエネルギと異なる第２のエネルギで撮像された複数の画像データとを、対応付ける第１の手段と、
前記第１の手段で対応付けられた画像データ間でエネルギサブトラクション演算を行う第２の手段として機能させることを特徴とする。
【００１０】
上記目的を達成するための放射線撮影方法は、第１の手段が、第１のエネルギのＸ線で撮像された複数の画像データと、第１のエネルギと異なる第２のエネルギで撮像された複数の画像データとを、対応付ける工程と、
第２の手段が、前記第１の手段で対応付けられた画像データ間でエネルギサブトラクション演算を行う工程と、
を備えることを特徴とする。
【００２０】
【発明の実施の形態】
本実施形態では、２つのエネルギのＸ線を用いて撮影することにより、呼吸動態を表す複数フレームで構成される２組の動画像を得、エネルギサブトラクション技術を用いて、当該２組の動画像から、肋骨などの骨部組織が除去された、肺野のみの動きを表す１組の動画像を作成することにより、上述の課題を解決した。
【００２１】
エネルギサブトラクションは、人体中のカルシューム成分の多い骨部と、それ以外のいわゆる軟部組織ではＸ線に対する透過率が異なり、またＸ線エネルギ変化に対する透過率の変化量も異なる性質を利用して、２つ以上のエネルギのＸ線で撮影された複数の画像間の演算により骨部もしくは軟部を分離する技術である。
【００２２】
たとえば、あるＸ線エネルギＥ１の時の、骨部のＸ線透過率をμ_１厚みをＴ_１、軟部のＸ線透過率をμ_２厚みをＴ_２とし、入射Ｘ線量をＨ１とすると、透過するＸ線量Ｈ２は
Ｈ２＝Ｈ１ｅｘｐ（−Ｔ_１μ_１−Ｔ_２μ_２）・・・・（式１）
また、あるＸ線エネルギＥ２の時の、骨部のＸ線透過率をμ_３厚みをＴ_１、軟部のＸ線透過率をμ_４厚みをＴ_２とし、入射Ｘ線量をＨ３とすると、透過するＸ線量Ｈ４は
Ｈ４＝Ｈ３ｅｘｐ（−Ｔ_１μ_３−Ｔ_２μ_４）・・・・（式２）
これらＨ２、Ｈ４は画素値として現れるが、この（式１）と（式２）それぞれの対数を計算すると以下になる。
【００２３】
ｌｎ（Ｈ２）＝ｌｎ（Ｈ１）−Ｔ_１μ_１−Ｔ_２μ_２・・・・（式３）
ｌｎ（Ｈ４）＝ｌｎ（Ｈ３）−Ｔ_１μ_３−Ｔ_２μ_４・・・・（式４）
ここで、（式３）と（式４）から骨部の情報であるＴ_１と軟部の情報であるＴ_２をそれぞれ独立して求めることができる。
【００２４】
たとえば、（式３）をμ_４倍、（式４）をμ_２倍して、差分をとると軟部組織の情報Ｔ_２は消去され、出力される画像情報は骨部のみになる。一般にこの両者にかける係数は物理的に求まりにくく、最適な値を試行錯誤して求める場合が多い。すなわち、まとめると、対数変換された２つのエネルギで取得された画像Ｉ１とＩ２があった場合、あらたな画像Ｉ３は
Ｉ３＝Ｋ１・Ｉ１−Ｋ２・Ｉ２・・・・（式５）
で計算でき、Ｋ１、Ｋ２を適当に選択すると、骨部もしくは軟部組織のみの画像を得ることができる。このＫ１、Ｋ２は撮影条件、被写体の厚みなどに左右されるため、画像内でも一定でない場合もある。
【００２５】
また、エネルギサブトラクションの特徴の一つとして、良性腫瘍と悪性腫瘍との判別が行える可能性があることを挙げることができる。良性腫瘍と悪性腫瘍とではＸ線エネルギの変化に対する放射線透過率の変化率が異なる場合があるからである。
【００２６】
以下、本発明の好ましい実施形態を、添付図面を参照しながら詳細に説明する。
【００２７】
（第１の実施形態）
図１は本発明を実施した第１の形態を模式的に示すブロック図であり、１はＸ線発生装置（Ｘ線源）を表し、図で示す破線の矢印の方向へＸ線を放射する。このＸ線発生装置は外部からのコントロール入力に従って少なくとも２種類のエネルギを持つＸ線を発することが可能であり、瞬時にそのＸ線管球電圧を切り替えることが可能に構成されている。２は被写体である人体（患者）を示し、この場合胸部の撮影をするため背面よりＸ線を入射し胸部の画像を撮影する。３はＸ線強度分布を画像化するフラットパネルディテクタであり、受像面に画像を構成する複数の画素に対応する複数の受像素子（単に画素ともいう）がマトリックス状に配置されている。通常このマトリックス状の画素は１００μｍ〜２００μｍピッチの等間隔に配置さている。４は人体へ放射されるＸ線量が最適なＸ線量になるようにコントロールするため、人体を透過したＸ線量をモニタするフォトタイマと呼ばれる装置である。フラットパネルディテクタから出力される各画素値は初期段階ではアナログ電圧信号であるため、５で示すＡ／Ｄ変換器で数値であるデジタル値に変換される。通常、少なくともこのＡ／Ｄ変換器はフラットパネルディテクタと同筐体に内蔵され、外見的にはフラットパネルディテクタから直接デジタル値（画像データ）が出力されるように認識される。６および７は画像データを一旦格納するバッファメモリであり、それぞれいわゆるダブルバッファとして働き、片側が読み込み中はもう一方は読み出しを行う構成で、読み出し処理の連続性を保つ働きをスイッチ１６を用いて行う。９は差分器であるが、あらかじめメモリ８に記憶されたＸ線を放射せずにセンサから取得された画像データ（オフセット画像）を実際の被写体の画像データから差し引く機能を有する。具体的にはスイッチ１７をＢ側に倒し、Ｘ線を放射せずに画像データをメモリ８に格納し、実際の使用時にはスイッチ１７をＡ側に倒して用いる。１０は参照テーブル（ＬｏｏｋＵｐＴａｂｌｅ）であり、画像データ値を値変換する機能を有する。具体的にはこの参照テーブルは入力値をその対数値に比例した値に変換するように設定される。１１のブロックは差分器であり、あらかじめ被写体を置かずに、Ｘ線のみを放射して取得され、対数値に変換された画像データを１２−１もしくは１２−２のメモリに格納し、当該画像データを実際の被写体の画像から差し引き、フラットパネルディテクタの画素ごとのゲインのばらつきを補正するためのものである。具体的には被写体が無い状態で放射線を放射し、スイッチ１７をＡ側、スイッチ１８をＢ側に倒して、メモリ１２−１もしくは１２−２へゲインばらつきを表す画像データを格納し、実際の被写体画像の場合にはスイッチ１７、スイッチ１８それぞれＡ側に倒して用いる。１２−１、１２−２の２つのメモリは、エネルギの異なるＸ線によるゲインばらつきのそれぞれに対応している。すなわち、エネルギＥ１のＸ線によるゲインばらつきをメモリ１２−１に格納し、エネルギＥ２のＸ線によるゲインばらつきをメモリ１２−２に格納する。一般にフラットパネルディテクタのＸ線に対する感度はＸ線のエネルギに依存するものであるため、このように２つのゲイン補正用メモリを設けている。そして、使用するＸ線のエネルギに応じて、スイッチ２１をＡ側もしくはＢ側に設定して両者のメモリの出力を使い分ける。尚、この両者がほぼ同じである場合、メモリ１２−１及び１２−２を一つのメモリでまかなうことも可能である。１９は欠陥画素を補正する機能を有するブロックであり、２０で示すメモリにあらかじめ記憶された、使用されるフラットパネルセンサの欠陥画素について、そのデータを回りの正常な画素データから推測し、欠陥画素補正を行うためのものである。この補正には一般に周りの正常な画素データ値の平均値が用いられる。このようにして各種の補正がなされた画像は一旦画像メモリ１３に格納され、ついで１４のファイリング装置に記録される。このようにして得られた複数の画像データは２２で示されるプログラムおよび計算機システムからなる画像解析・表示装置で読み出され、後述される画像解析がなされる。画像解析・表示装置２２は、放射線撮影装置内に設けられることも、まったく別の外部の計算機システムとして構成されることも可能である。画像解析・表示装置２２は、一般の計算機システムを用いて構成することができる。
【００２８】
１５で示されるブロックは撮影の制御を行うコントローラ（制御機構）を表し、所定のタイミングでフラットパネルディテクタ３を駆動すると共に、Ｘ線発生装置１へＸ線パルスの放射タイミングのトリガおよびＸ線のエネルギ制御信号を出力する。
【００２９】
図２は撮影順序を図示するものであり、左のＡの列は被写体である患者の動作、中央のＢの列は操作者である放射線技師の動作、Ｃの列はＸ線撮影装置のモードを表す。最初の（Ａ１）の時点で患者は操作者の指示に従って、撮影台の前（図１の２の位置）に立つ。次の（Ｂ１）の時点で操作者は息を吸い込むように指示し、ついで（Ｂ２）の時点から息をゆっくり吐き出すように指示を出す。患者はその指示に従って息を吸い込んだ（Ａ２）後、ゆっくりと息を吐き出す（Ａ３）が、操作者は図１の撮影装置を操作し、患者の呼吸動態を連続的に撮影する（Ｃ１）。この撮影間隔は秒あたり３〜１０画像程度になる。操作者は患者の様子を見ながら適当な時間（数秒）経過したあとに、今度はゆっくりと息を吸い込むように指示を出す（Ｂ３）。この時点でも、連続的なＸ線撮影は継続している。操作者は患者の様子をみながら、患者が息を吸いきった時点で、撮影終了を告げる（Ｂ４）。そして、呼吸動態を表す連続的な画像データの収集が終了する（Ｃ２）。
【００３０】
図３はこのときの状態を模式的に表すタイミングチャートであり、上段がＸ線パルス、下段がフラットパネルディテクタ３を含むセンサシステムの動きを表す。患者が操作者の指示に従って呼気もしくは吸気を行う間、Ａ１及びＡ２で示すＸ線パルスが発せられる。このＡ１とＡ２に対応するＸ線パルスは異なるエネルギ（たとえば管球電圧１２０ＫＶｐと１００ＫＶｐ）を持つものである。このＡ１のパルスとＡ２のパルスの間隔Ｔ１は約５０ｍＳｅｃと比較的短く設定する。このＴ１の期間にＸ線の強度分布をフラットパネルディテクタへ蓄積し、読み出す動作を行う。この蓄積に要する時間は３〜４ｍＳｅｃ程度で、残りの時間でフラットパネルディテクタからデータを読み取る。すなわちフラットパネルディテクタは、Ｘ線パルスが発せられている間は画像情報を蓄積し、次にそれを読み取る。このＡ１もしくはＡ２のパルス幅は、規定の一定値を用いてもよく、また、図１にあるフォトタイマ４を用いて制御することも可能である。この場合、フォトタイマ４で計測される総Ｘ線量（積分値）が所定の設定値になった時点で、コントローラ１５はＸ線放射停止信号をＸ線発生装置１に送出し、Ｘ線放射を止める（パルスを停止する）。また、図３に示されるように、次のサイクルの撮影までの時間間隔Ｔ２を約３００ｍＳｅｃ程度の比較的長い時間に設定する。通常、ゆっくり呼吸する場合に必要な時間サンプリング間隔は秒あたり３〜１０フレーム程度であり、ここでは患者に比較的ゆっくり呼吸してもらい、秒あたり３フレーム程度の撮影を行うものとしている。
【００３１】
上述の操作により、深呼吸をゆっくりするだけの従来の健康診断とあまり変わらない動作を患者に行ってもらうことで、呼吸動態撮影が行える。
【００３２】
本実施形態の呼気・吸気の順、回数は上述のものに限られるものではない。
【００３３】
このような操作により、患者の呼吸の各位相に対応した複数のフレームからなる動画像が得られる。
【００３４】
以下、図１の画像解析・表示装置２２で行われる動作を説明する。
【００３５】
一般に呼吸による肺野部の動きは肺の動きがほとんどであり、肋骨はほとんど移動しない。したがって、図３のＴ１（５０ｍＳｅｃ）間隔で撮影された２枚の画像においては肺の動きは多少あるものの、肋骨はほとんど移動していないと考えられる。ここでは、前述の（式５）に従って、最適な係数Ｋ１、Ｋ２を用いて、それぞれの呼吸位相における２組の異なるエネルギで撮影された画像から軟部組織である肺のみを抽出する。この場合、図１の参照テーブル１０によりすでに対数値の画像になっているため、そのままで係数をかけて差分を取る演算を行うことが可能である。このとき肺自体には多少の動きがあるが、それほど問題にはならない。これは、３００ｍＳｅｃの時間間隔でサンプリングされる中５０ｍＳｅｃ隔たった２つの画像を平均化したような結果になるためである。これは時間アパーチャすなわちＸ線の蓄積時間が５０ｍＳｅｃになったようなものである。もし問題がある場合には、画像間で近似するブロックを見出すことにより、画像間の移動ベクトルを見つけ、その後ワーピング処理を行って画像間の動きを補償してから上述のエネルギサブトラクションを行えばよい。
【００３６】
このようにして求められた、軟部組織のみの動き画像を観察するか、もしくは、当該動画像フレーム間の動きベクトルを演算することにより、肋骨等に邪魔されない呼吸動態の診断が容易に行える。
【００３７】
（第２の実施形態）
第２の実施形態でも第１の実施形態と同様の図１の構成を用いるが、患者の呼吸動態画像中から肺が同一の呼吸位相状態の画像を見出し、エネルギサブトラクション演算を行う点が第１の実施形態とは異なる。
【００３８】
図４は、一例として、呼吸動態を観察するためのＦ１〜Ｆ１６の１６枚の画像を示し、Ｆ１が最大吸気、Ｆ８およびＦ９が最大呼気の状態である。この息を吐きながらの状態と息を吸い込みながらの２つの状態から肺が同一状態となる２つの画像を見出せれば、その２画像間でエネルギサブトラクションが可能になる。同図でＦ１〜Ｆ８の８枚の区間はエネルギＥ１で撮影され、Ｆ９〜Ｆ１６の８枚の区間はエネルギＥ２で撮影される。Ｆ１及びＦ１６からＦ１’で示されるエネルギサブトラクション画像としての軟部組織画像を得る。同様に、Ｆ２とＦ１５からＦ２’、Ｆ３とＦ１４からＦ３’、Ｆ４とＦ１３からＦ４’、Ｆ５とＦ１２からＦ５’、Ｆ６とＦ１１からＦ６’、Ｆ７とＦ１０からＦ７’、Ｆ８とＦ９からＦ８’の、肺の各状態に対応した８枚のエネルギサブトラクション画像を得る。最大呼気であるＦ８とＦ９の期間は、患者にその状態を少し長めに保持するように指示することが望ましい。
【００３９】
肺が同一の状態（単に呼吸の位相ともいう）を的確に見出すことは一般に困難であるが、次のような各種手法を提案することができる。患者に伸縮可能なベローズを通して呼吸してもらい、このベローズの長さを基準として、同一呼吸位相を見出す方法、患者の体軸方向の人体インピーダンスを連続的にモニタし、インピーダンス値に基づいて、同一呼吸位相を見出す方法等が考えられる。
【００４０】
図５は図１の構成に、２３で示されるベローズ等を追加したもので、当該ベローズを介して患者に呼吸してもらい、その長さ情報を不図示の装置で読み取り、同一長さ時点を同一呼吸位相として、２つの時点でＸ線撮影が行われるように構成されるシステムを模式的に示すものである。
【００４１】
図６は２４で示される人体のインピーダンス測定装置により、人体のインピーダンスを連続的にモニタしながら、同一インピーダンス値であるときを同一呼吸位相として、２つの時点でＸ線撮影が行われるように構成されるシステムを模式的に示している。
【００４２】
本実施形態では、呼気状態と吸気状態をそれぞれ異なるエネルギで撮影したが、患者の呼気状態及び吸気状態を２回ずつ撮影する場合には、呼気状態及び吸気状態をそれぞれ２つのエネルギで撮影可能なため、呼気状態もしくは吸気状態の画像どうしでエネルギサブトラクション演算を行うことが可能となる。
【００４３】
また画像間に多少のずれがある場合には、前述のワーピング処理により画像間のずれを補償することも可能である。
【００４４】
（他の実施形態）
本発明を、複数の装置（例えば、１以上の画像処理装置、１以上のインターフェース、１以上の放射線撮影装置、及び１以上のＸ線発生装置、等）から構成されるシステム、又は画像処理装置及び放射線撮影装置を統合した構成等のいずれにも適用可能であることは言うまでもない。
【００４５】
また、本発明の目的は、上述の実施の形態の機能を実現するためのソフトウェアのプログラムコードを記憶した記憶媒体を、システム或いは装置に供給し、そのシステム或いは装置のコンピュータ（ＣＰＵ又はＭＰＵ等）が記憶媒体に格納されたプログラムコードを読み出して実行することによっても、達成されることは言うまでもない。
【００４６】
この場合、記憶媒体から読み出されたプログラムコード自体が上述の実施の形態の機能を実現することとなり、当該プログラムコード、及び当該プログラムコードを記憶した記憶媒体は本発明を構成することとなる。
【００４７】
プログラムコードを供給するための記憶媒体としては、ＲＯＭ、フロッピー（Ｒ）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード等を用いることができる。
【００４８】
また、コンピュータが読み出したプログラムコードを実行することにより、上述の実施の形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ等が実際の処理の一部又は全部を行い、その処理によって上述の実施の形態の機能が実現される場合も本発明を構成することは言うまでもない。
【００４９】
さらに、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された拡張機能ボードやコンピュータに接続された機能拡張ユニット等に備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニット等に備わるＣＰＵなどが実際の処理の一部又は全部を行い、その処理によって上述の実施の形態の機能が実現される場合も本発明を構成することは言うまでもない。
【００５０】
以上説明したように、本実施形態によれば、エネルギサブトラクション処理を適用して骨部情報を除去することにより、より安定した肺野部の動態を観察・解析できるようになった。
【００５１】
【発明の効果】
本発明によれば、エネルギサブトラクション処理の適用により、安定した肺野部動態の把握が可能な放射線撮影装置、放射線撮影方法、放射線画像用システム、プログラム、及びコンピュータ可読記憶媒体を提供することができる。
【図面の簡単な説明】
【図１】ブロック図
【図２】撮影シーケンスの一例を表す図
【図３】撮影に係るタイミングチャート
【図４】同一呼吸位相の２画像間でエネルギサブトラクションを行う模式図
【図５】呼吸検知をベローズを用いて行う場合のブロック図
【図６】呼吸検知を人体インピーダンス測定により行う場合の模式図
【符号の説明】
３フラットパネルディテクタ
１５コントローラ

Claims

第１のエネルギのＸ線で撮像された複数の画像データと、第１のエネルギと異なる第２のエネルギで撮像された複数の画像データとを、呼吸位相に応じて対応付ける第１の手段と、
前記第１の手段で対応付けられた画像データ間でエネルギサブトラクション演算を行う第２の手段とを備え、
呼気状態と吸気状態とで異なるエネルギを用いて撮影を行うように構成されていることを特徴とする放射線撮影装置。
前記第１の手段は、撮影時間に応じて対応付けることを特徴とする請求項１に記載の放射線撮影装置。
前記複数の画像は、呼吸状態にある被写体を時系列に撮像してえられた画像データであることを特徴とする請求項１又は２に記載の放射線撮影装置。
前記エネルギサブトラクション演算の対象となる画像間の動きを補償する動き補償手段とを備えることを特徴とする請求項１に記載の放射線撮影装置。
前記呼吸状態を検出する呼吸状態検出手段を更に備え前記呼吸状態検出手段の出力に基づいて前記呼吸位相が求められることを特徴とする請求項１に記載の放射線撮影装置。
コンピュータを、
第１のエネルギのＸ線で撮像された複数の画像データと、第１のエネルギと異なる第２のエネルギで撮像された複数の画像データとを、呼吸位相に応じて対応付ける第１の手段と、
前記第１の手段で対応付けられた画像データ間でエネルギサブトラクション演算を行う第２の手段として機能させ、
呼気状態と吸気状態とで異なるエネルギを用いて撮影を行うように機能させることを特徴とするプログラム。