放射線源と二次元検出器を用いたトモシンセシス撮影で得られた画像を処理する画像処理装置であって、トモシンセシス撮影で前記二次元検出器が出力する複数の投影データを取得する取得手段と、トモシンセシス撮影で得られた複数の投影データを、前記放射線源の照射中心方向と垂直となるように仮想的に設定された仮想CT検出面での仮想投影データに変換することなく、解析的再構成処理を実行し被検体の断層画像を得る再構成手段と、を有することを特徴とする画像処理装置。前記再構成手段は、前記投影データに再構成フィルターを合成し、逆投影処理を行うことを特徴とする請求項1に記載の画像処理装置。前記再構成手段は、前記複数の投影データの各画素値に対応する二次元検出器上の位置と、前記放射線源との配置関係に基づいて再構成処理を行うことを特徴とする請求項1または2に記載の画像処理装置。前記再構成手段は、前記複数の投影データと再構成フィルターと前記配置関係とに基づいて逆投影処理を行うことを特徴とする請求項3に記載の画像処理装置。前記再構成手段は、前記配置関係から定まる第1の係数を乗じつつ前記投影データと前記再構成フィルターとを合成し、さらに該合成されたデータに前記配置関係から定まる第2の係数を乗じて逆投影処理を行うことを特徴とする請求項4に記載の画像処理装置。前記再構成手段は、仮想検出器の座標系への投影によって、投影データにおいて近接する位置の画素値同士を加算し、前記仮想検出器における仮想画素を補間する処理を行わないことを特徴とする請求項1乃至5のいずれか1項に記載の画像処理装置。前記再構成手段は、前記二次元検出器で得られた投影データに再構成アルゴリズムを直接適用することを特徴とする請求項1乃至6のいずれか1項に記載の画像処理装置。前記再構成手段は、前記二次元検出器と平行な平面上から各再構成点の画素値を直接再構成することを特徴とする請求項1乃至7のいずれか1項に記載の画像処理装置。前記再構成手段は、FeldkampのコーンビームCT再構成アルゴリズムにおけるフィルター畳み込みの軸を前記二次元検出器と平行な平面の軸に変形した再構成アルゴリズムに基づき前記投影データを再構成することを特徴とする請求項1乃至8のいずれか1項に記載の画像処理装置。トモシンセシス撮影とは、放射線源と二次元の撮像面を有する二次元検出器とによる撮影であって、前記放射線源の焦点位置と前記撮像面の中心位置との距離を変えながら前記放射線源および前記二次元検出器の少なくともいずれか一方を移動させつつ前記放射線源が放射線を複数回照射し、各照射に応じて二次元検出器で得られる複数の投影データを得る撮影方法であることを特徴とする請求項1乃至9のいずれか1項に記載の画像処理装置。前記再構成された前記被検体の断層画像を表示部に表示させる表示制御手段を更に有することを特徴とする請求項1乃至10のいずれか1項に記載の画像処理装置。以下の式に基づき演算する演算手段を更に有することを特徴とする請求項1乃至11のいずれか1項に記載の画像処理装置。前記演算手段は、前記式を離散化して演算することを特徴とする請求項12に記載の画像処理装置。放射線源と二次元検出器を用いたトモシンセシス撮影で得られた画像を処理する画像処理方法であって、トモシンセシス撮影で前記二次元検出器が出力する複数の投影データを取得する取得ステップと、トモシンセシス撮影で得られた複数の投影データを、前記放射線源の照射中心方向と垂直となるように仮想的に設定されたCT検出面での仮想投影データに変換することなく解析的再構成処理により再構成処理し被検体の断層画像を再構成する再構成ステップと、を有することを特徴とする画像処理方法。前記再構成ステップでは、前記投影データに再構成フィルターを合成し逆投影処理を行うことを特徴とする請求項14に記載の画像処理方法。前記再構成ステップでは、FeldkampのコーンビームCT再構成アルゴリズムにおける前記投影データと再構成フィルターとの合成軸を前記二次元検出器に平行な軸に変形した再構成アルゴリズムにより再構成処理を行うことを特徴とする請求項14に記載の画像処理方法。放射線源と二次元検出器を用いたトモシンセシス撮影で得られた画像を処理する画像処理をコンピュータに実行させるためのプログラムを記憶したコンピュータによって読み取り可能な記憶媒体であって、トモシンセシス撮影で前記二次元検出器が出力する複数の投影データを取得する処理と、トモシンセシス撮影で得られた複数の投影データを、前記放射線源の照射中心方向と垂直となるように仮想的に設定されたCT検出面での仮想投影データに変換することなく解析的再構成処理を実行し被検体の断層画像を再構成する処理と、をコンピュータに実行させるためのプログラムを記憶したコンピュータによって読み取り可能な記憶媒体。しかしながら、図7に示すようにトモシンセシスで均等に配置されている検出器7001で得られた画素を仮想CT検出器7002の配置に幾何変換しようとすると不均等の配置になる。このため各点の画素値を周辺の画素値で補間して生成する。この補間操作は空間的なローパスフィルタに相当するので、この幾何変換時に高周波情報が落ちてしまい、結果的にコーンビームCTアルゴリズムで再構成を行って得られる断層画像の空間分解能が低下してしまう。さらに、特許文献1の手法では、コーンビームCTの画像を幾何変換して、保持するメモリ空間が必要になる。また、幾何変換及び補間という余分な処理を行うため、処理時間が短い点にメリットがあるトモシンセシスにあって、再構成にかかる処理時間が長くなってしまう。そこで発明の一形態に係る画像処理装置は、放射線源と二次元検出器を用いたトモシンセシス撮影で得られた画像を処理する画像処理装置であって、トモシンセシス撮影で前記二次元検出器が出力する複数の投影データを取得する取得手段と、トモシンセシス撮影で得られた複数の投影データを、前記放射線源の照射中心方向と垂直となるように仮想的に設定された仮想CT検出面での仮想投影データに変換することなく、逆投影処理し被検体の断層画像の解析的再構成処理を実行する再構成手段と、を有することを特徴とする。本発明によれば、トモシンセシスで得られた投影画像をコーンビームCTの配置に幾何変換を行わずに、直接フィルタードバックプロジェクションを行い、断層像を求めることが可能になる。これにより従来の技術より、空間分解能の高く、処理時間の短い、トモシンセシスによる断層像を提供することができる。本発明の実施形態に係る画像処理装置は、放射線源と二次元検出器を用いたトモシンセシス撮影で得られた画像を処理する。以下、本発明の実施形態に係る画像処理装置以下、断層画像診断装置、画像処理方法以下、断層画像生成方法について添付図面を参照して説明する。図1は本発明の実施形態に係る断層画像診断装置の機能構成を示す図である。断層画像診断装置100は複数の照射角度からコーンビーム状にX線を照射できるX線管101、被検体102を寝かせる寝台103、X線を受像してX線画像を取得するX線検出器106を備えている。X線検出器106はここでは二次元の撮像面を有する二次元検出器である。X線管101とX線管から照射されるX線を検出するX線検出器106とは、被検体を挟むように対向する位置に配置されている。機構制御部105はX線管101とX線検出器106の位置の制御を行う。この断層画像診断装置100は単純撮影、長尺撮影に加えてトモシンセシス撮影が可能である。ここで単純撮影とは被検体102にX線を照射し1枚のX線画像を得る撮影方法である。長尺撮影とは、全身や全脊椎、下肢全長などの大きな被写体を複数の撮影に分けて部分ごとに撮影する撮影方法である。機構制御部105によりX線管101およびX線検出器を撮影部位に沿って移動させつつ複数の撮影を行う。長尺撮影で得られた画像を繋ぎ合わせて1枚の被写体画像を得る。トモシンセシス撮影では、X線管101の焦点位置とX線検出器106の撮像面の中心位置との距離を変えながらX線管101およびX線検出器106の少なくともいずれか一方を平行移動させる。X線管101がX線を複数回照射し、各照射に応じてX線検出器106で得られる複数の投影データを得る撮影方法である。X線管またはX線検出器106の移動は投影データから被検体102の撮影部位の断層画像を再構成する。撮影制御部104はX線検出器106を電気的に制御して、X線画像を取得する。X線発生装置107はX線管101を電気的に制御して、所定の条件でX線を発生させる。X線撮影システム制御部108は、機構制御部105と撮影制御部104とを制御して、複数のX線照射角度からのX線画像を取得する。このX線画像は、被検体102の撮影部位の投影データである。さらにX線撮影システム制御部108には、画像処理部109、画像保存部112が備えられ、1又は複数の情報処理装置コンピュータが内蔵される。コンピュータには、例えば、CPU等の主制御部、ROMReadOnlyMemory、RAMRandomAccessMemory等の記憶部が具備される。また、コンピュータには、GPUGraphicsProcessingUnit等のグラフィック制御部、ネットワークカード等の通信部、キーボード、ディスプレイ又はタッチパネル等の入出力部等が具備されていてもよい。なお、これらの各構成部は、バス等により接続され、主制御部が記憶部に記憶されたプログラムを実行することで制御される。画像処理部109はX線管101とX線検出器106を用いたトモシンセシス撮影で得られた画像投影データを処理する。X線撮影システム制御部108の指示に従い、取得したX線画像を再構成し、断層画像を生成する。このために、画像処理部109は前処理部113、係数計算部114、重畳積分部115、重み計算部116、バックプロジェクション部117を備える。前処理部113にはX線撮影システム制御部108が撮影制御部104を介して、X線検出器106から取得した様々なX線照射角度からの複数のX線画像以下、投影画像または投影データと呼ぶが入力される。ここで、投影画像には欠陥補正やゲイン補正などが行われ、対数変換log変換が行われる。これによりX線検出器106やX線管101に起因するX線照射ムラや画素欠陥などが補正される。係数計算部114は、X線検出器106上の検出点とX線管101との幾何配置から決まる係数の計算を行う。ここで幾何配置とは、X線管101とX線検出器106の相対的かつ物理的な位置関係をいい、X線管101の焦点と、X線検出器106の撮像面と各画素の位置との関係である。撮影のたびX線検出器106とX線管101との幾何配置は異なるので、1回のX線照射ごと、画素の位置ごとに係数が異なる。重畳積分部115は係数計算部114で計算された係数と、X線検出器106の検出点の画素値との積に、再構成用のフィルター関数を重畳積分する。この重畳積分はX線検出器106の検出面に平行な座標軸上で行われる。再構成用のフィルター関数は、一般的な再構成に用いられるrampフィルターやsheppLoganフィルターなどを用いればよい。この結果、投影画像に再構成フィルター処理が行われたフィルター処理された画像フィルタード画像が生成される。重み計算部116は再構成点とX線管101との相対的な位置関係を示す幾何配置から決まる重み係数の計算を行う。再構成点とは、フィルター処理された投影画像を再構成して生成する際に、投影画像の画素の位置を示すアイソセンタを原点とした3次元座標上の点である。アイソセンタとはX線管の照射方向を変えたときに基準軸ビーム中心、照射中心が交わる回転中心のことである。重み計算部116は、X線管の照射方向を変えたときにビーム中心が交わる回転中心を原点として、フィルター処理された投影画像の画素の位置を示す3次元座標上の点と、X線管との相対的な位置関係を示す幾何配置により求められる重み係数を計算する。バックプロジェクション部117は重畳積分部115で計算されたフィルタード画像に、重み計算部116で計算された再構成点とX線管との幾何配置から決まる重みを乗算しながら、逆投影処理バックプロジェクションを行う。これにより、フィルタードバックプロジェクションによる再構成処理が行われ、所望の被検体の断層画像を再構成することができる。尚、重みの乗算は、逆投影処理バックプロジェクションにおいて必須のものではないが、被検体の断層画像投影画像の再構成を、より精度よく行うためには、重みの乗算を行うことが好ましい。断層画像の再構成に用いる具体的な係数、計算理論式については後述する。通常のCTなどの断層画像診断装置と断層画像診断装置100の異なる点は、一般撮影装置や透視撮影装置などを応用して、断層画像が撮影される点である。このために、X線管101は被検体102の周りを180度より小さい角度、例えば±40度程度の角度で撮影を行い、X線検出器106は水平方向にスライドするか、もしくは固定状態で撮影を行う。これにより、大掛かりなCT装置を用いなくても、所定の範囲で照射角が変更可能な一般撮影装置で断層像を撮影することができ、断層画像診断装置の大幅なコストダウンができる。また、短い撮影時間、開放的な空間で撮影ができるため、被検体の負担も小さくできる。次に、図2を用いて、図1に示す断層画像診断装置100における断層画像生成処理の流れの一例について説明する。まず、S201では投影画像の取得が行われる。これは、X線管101のX線照射角度を-40度~40度まで変えながら、被検体102をX線で撮影することにより行われる。撮影枚数は任意の数が可能であるが、例えば角度を1度ずつ変えながら、80枚の投影画像を15FPSで撮影すると6秒程度で画像の収集ができる。X線の撮影条件も任意の条件が設定可能であるが、胸などの撮影では100kV、1mAs程度で行えばよい。また、X線検出器106とX線管101との間の距離は透視撮影装置や一般撮影装置の設定範囲にように100cm~150cm程度に設定される。X線管101は円弧状の軌道を描くことが好ましいが、透視撮影装置や一般撮影装置では機構的に円弧状の軌道を描くことが難しい場合が多い。この場合、寝台103の長軸方向にX線管101を移動しながら、X線照射角βを変えて撮影すれば良い。このときのX線照射角と移動したX線管101の位置関係は、Dtanβで与えられる。ここで、Dはβ0のときのX線管101の焦点とアイソセンタとの距離である。また、X線検出器106をX線管101に対して平行移動させる。このときの平行移動量はPtanβで与えられ、PはアイソセンタとX線検出器106の中心までの距離である。このようにX線検出器106を平行移動させれば、X線管101のX線照射方向が変わっても、基準軸はX線検出器106の中心を常に通るようになる。透視撮影装置や一般撮影装置によっては、X線検出器106を平行移動させる機構が無い場合もある。この場合はアイソセンタをX線検出器106の特定の位置、たとえば中心位置と一致させれば、X線検出器106を平行移動させなくてもX線検出器106を固定した状態でトモシンセシス撮影が可能である。ただし、X線検出器106を移動する機構が無いと、X線照射角度βが大きくなるにつれ、X線照射範囲がX線検出器106からはみ出るため、有効視野FOVが欠け、再構成できる断層画像の範囲が小さくなる。S201で得られた一連の投影画像は画像処理部109に入力され、まず前処理S202が行われる。ここでは、X線検出器106の製造過程で生じた、欠陥画素の補正や、X線管101に起因する照射ムラを補正する処理が行われる。これらの処理はX線検出器で一般的に行われる処理を用いればよい。また、前処理部では式1で示される対数変換log変換が行われる。ここで、Iは投影画像の画素値であり、logは自然対数である。この処理により、投影画像の画素値はX線減弱係数の足し合わせたものになる。このX線減弱係数の加法性に基づいてX線画像の再構成が行われる。次に係数計算部114では、X線検出器106上の検出点とX線管101の相対的な位置関係を示す幾何配置から決まる係数計算S203を行う。具体的にはこの係数は式2であらわされる。各変数の関係は図3の再構成座標系により表わされる。3次元座標軸x、y、zは再構成座標空間を示し、アイソセンタを原点としている。xz平面はX線検出器106の検出面に平行な面で、アイソセンタ301を通る面である。また、y軸はX線検出器106の検出面に垂直な法線となっている。xtとztはX線検出器106上の点とX線管101の焦点302を結ぶ直線303が、xz平面と交わる点のx座標、z座標である。y軸とX線管101の基準軸の成す角βが、X線照射角度投影角度である。式2は直線303と直線304の成す角の余弦値を表している。直線304は直線303がxz平面と交わる点からz軸に降ろした垂線が交わる点と焦点302を結ぶ線である。なお、この式2は発明を実施するために具体的に係数を数式で表したものである。したがって、式2に相当する係数の計算は他の数学的方法で計算しても良く、式2による計算に本発明を限定するものではない。重畳積分部115は式2の係数とX線検出器106上の対応する点との積を、再構成用のフィルター関数と重畳積分して、フィルター処理を行う。具体的には式3で表わされる。hxt’xtはrampフィルターやsheppLoganなどの再構成用フィルター関数である。qが得られる。ここで、式2と式3の示す重畳積分処理を直感的に説明するために、図4と図5を用いて2次元の再構成の説明を行う。本来、フィルタードバックプロジェクションで数学的に厳密な再構成を行うためには、図4に示すようにX線を線状にコリメートして、回転をしながら、平行スキャンを繰り返して、t軸上投影データpを得る必要がある。このようにして得られた投影データは式4を用いることで再構成が可能である。式4はCT再構成の原理式の一形態であるRadonTransformを同値変形して得られる式である。ここで、htt’は再構成用のフィルター関数である。しかしながら、上述の方法では平行スキャンと微小回転を繰り返すため、データの取得に時間がかかり、被検体の被曝量も極めて大きくなってしまう。したがって、今日では医療用のX線CTに用いられることはほとんど無い。そこで、ファンビームやコーンビームを用いたCTにより、平行スキャンを行うことなく、2次元や3次元の断面像を生成することが今日では盛んに行われている。ファンビームやコーンビームのフィルタードバックプロジェクションによる再構成は式4をファンビームやコーンビームに適した形に変換し、直接再構成ができるようにして行われる。そこで本発明ではファンビームやコーンビームと同様にトモシンセシスに合わせた形に式4を変換することで、特許文献1のようにコーンビームCTへの幾何変換及びそれに伴う補間処理をすることなく、直接に断層像を得る技術を提供する。本実施形態で提案する式3は式4の投影データに式2を乗じて、重畳積分軸をt’からxtに置き換えた形になっている。前述のように、式2は直線303と直線304の成す角の余弦値を表している。これは直感的に図5を用いて理解できる。図5aの被写体501は厚さ1、X線減弱係数αの被検体である。これに、図4のように平行なビームでスキャンをした場合、投影画像は502のように均一な分布になる。X線の物質内の減弱は式5で表わされるため、式1の対数変換log変換を行えば減弱係数の分布が得られる。一方で、図5bのように平行でないファンビームを用いて得られる投影画像は503のように円弧状の分布になる。これは、中心のビーム505に比較して、周辺のビーム504は被写体501の内部を1cosφ長い経路通るためである。したがって、X線の物質内の減弱は式6のようになる。これから正しい減弱係数の分布を得るためには、式1の対数変換log変換を行ったあとに、cosφを乗じてやれば良い。図5bは直感的な描像を示すために、ファンビームを用いた簡単な例を示したが、コーンビームやトモシンセシスの場合も同様な考えが成り立つ。実際には式4をトモシンセシスの直接再構成に適した形にする過程で数学的に式3ただしzt0の場合が導出されることになる。以上2次元のRadonTransformから得られる2次元のCT再構成理論およびそのトモシンセシス再構成アルゴリズムへの変換を説明したが、これを3次元に拡張することでトモシンセシスの再構成式が得られる。三次元への拡張は、理論に従い3次元のRadonTransformに従い行ってもよいし、FeldkampのコーンビームCT再構成アルゴリズム同様にコーンビームを複数のファンビームの集まりと考えることで2次元式から3次元式を導出しても良い。S203とS205が全ラインに対して処理されていない場合S206-No、次のラインに処理が進められS204、再度、係数計算が得られる。S203とS205を全ラインに処理することで得られた2次元フィルタード画像Gを後述するS208でバックプロジェクションすることで、断層像を得ることができる。S208でバックプロジェクションする際にS207で計算される重みを乗じながらバックプロジェクションを行う。S207において、重み係数の計算を行なう。この重みは再構成点とX線管のトモシンセシス幾何配置により決まり、具体的には式7で表わされる。図6にバックプロジェクション逆投影処理の概略図を示す。602はトモシンセシスの再構成で生成される断層面を示し、r→はアイソセンタを原点とした602上の再構成点601を示す3次元ベクトル、603は平行移動するX線管の焦点、604は603と共に移動するX線検出器、y^はX線管の基準軸ビーム中心に沿った単位ベクトルである。式7で計算される重み係数を逆投影処理バックプロジェクション時に乗ずることで、X線管101のX線照射角度βと平行移動あるいは固定するX線検出器106の検出面に平行な座標軸上xtに関する再構成式が得られる。これにより、トモシンセシス画像データからの直接再構成が可能になる。最終的にトモシンセシスの直接再構成はS208のバックプロジェクションを行うことで達成され、断層像を得ることができる。S208のバックプロジェクションは式7で計算される重みをフィルタード画像の対応する画素値に乗じながら、X線管の照射角度範囲で積分を行う。具体的には式8で表される。ここで、βmは最大X線照射角、-βmは最小のX線照射角度であり、fr→は断層像の画素値である。すなわち、トモシンセシスのバックプロジェクションは断層位置r→とX線管の焦点603を結ぶ直線がX線検出器604と交差する点の画素値を全投影角度βに対して総和することで行われる。ただし、X線検出器604の画素値はS205、式3でフィルタリングされたものである。また、総和時にはS207、式7で計算されるトモシンセシス幾何配置で決定される重みを乗じながら行う。この式は、FeldkampのコーンビームCT再構成アルゴリズムにおけるフィルター畳み込みの軸を二次元検出器と平行な平面の軸に変形した再構成アルゴリズムに基づき投影データを再構成するための式である。この再構成アルゴリズムは、FeldkampのコーンビームCT再構成アルゴリズムにおける投影データと再構成フィルターとの合成軸を二次元検出器に平行な軸に変形した再構成アルゴリズムに基づくものである。また、8式は、X線検出器106と平行な平面上から各再構成点の画素値を直接再構成する式である。X線検出器106で得られた投影データに補間処理をせずに再構成アルゴリズムを直接適用する式である。これにより、特許文献1に記載されるようなコーンビームCT幾何配置への幾何変換、及びそれに起因する補間を行わずに直接トモシンセシス再構成を行うことが可能になる。上述の式8では、式3で得られるフィルタード画像Gを、トモシンセシスの幾何配置で決定される係数を乗じながら逆投影処理する。式8を用いることで、特許文献1のように実際の検出器で得られた投影データを変換し仮想CT検出器における仮想投影データを得ることなく、直接に実投影データを再構成し断層画像を得ることができる。たとえば胸部を撮影部位としてトモシンセシス撮影した際に、一般撮影では確認しがたい胸骨の奥にある構造を断層画像として画像化することができる。これにより、変換に伴う隣接ピクセルの加算による補間処理を行わず再構成が可能となるため、処理負荷を抑えつつ画質を向上させた断層画像を得ることができる。以上が本発明の代表的な実施形態であるが、本発明は、上記及び図面に示す実施形態に限定することなく、その要旨を変更しない範囲内で適宜変形して実施できるものである。例えば、本発明は、例えば、システム、装置、方法、プログラム若しくは記憶媒体等としての実施態様を採ることもできる。具体的には、複数の機器から構成されるシステムに適用してもよいし、また、一つの機器からなる装置に適用してもよい。その他の実施例また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。なお、上述の例に加えて、得られた断層画像を表示する表示部、及び表示制御部を備えていてもよい。上述の例では、理論式8に基づく計算を行うこととしたが、デジタル計算機での処理に伴う理論式8の離散化を行っても良いことはもちろんである。この場合、式3はフィルター合成の処理は畳み込み演算となる。なお、式8他の式をデジタル計算機で処理する場合において、デジタル値で計算処理を行うことに起因する誤差や、計算量その他の理由により生じる近似誤差は許容される。上述の例では画像処理部109が再構成処理を行う例を示したが、これを単体の画像処理装置が実行することとしても良いし、複数の装置からなるシステムで機能を分散させて処理することとしてもよい。上述の例では、コーンビーム状のX線を発生させるX線源を用いる例を示したが、これに限らず、絞りの形状等によって四角錘その他の形状であってもよい。上述の例では、X線撮影を例に示したが、その他の放射線撮影にも用いることができる。上述の実施例の再構成処理は、数学的な解析処理を伴ういわゆる解析的再構成の処理であり、単純加算法やシフト加算法とは異なる処理である。以上のとおり、上述の実施例では、画像処理部109の前処理部113または不図示の取得部がトモシンセシス撮影で二次元検出器が出力する複数の投影データを取得する。係数計算部114、重畳積分部115、重み計算部116及びバックプロジェクション部117逆投影処理部から再構成処理部は構成される。再構成処理部はトモシンセシス撮影で得られた複数の投影データをX線管の照射中心方向と垂直となるように仮想的に設定された仮想CT検出面での仮想投影データに変換することなく逆投影処理し被検体の断層画像を再構成する。すなわち、再構成処理部はトモシンセシス撮影で得られた複数の投影データを、仮想的に設定された仮想CT検出面での仮想投影データに変換することなく、解析的再構成処理を実行し被検体の断層画像を得ることができる。再構成処理部は、複数の投影データの各画素値に対応する二次元検出器上の位置と、放射線源との配置関係に基づいて再構成処理を行う。再構成処理部は、は、複数の投影データと再構成フィルターと配置関係とに基づいて逆投影処理を行う。再構成処理部は、配置関係から定まる第1の係数を乗じつつ投影データと再構成フィルターとを合成し、さらに合成されたデータに配置関係から定まる第2の係数を乗じて逆投影処理を行う。再構成処理部は、仮想検出器の座標系への投影によって、投影データにおいて近接する位置の画素値同士を加算し、仮想検出器における仮想画素を補間する処理を行わない。再構成処理部は、二次元検出器放射線検出器で得られた投影データに再構成アルゴリズムを直接適用する。再構成処理部は、二次元検出器放射線検出器と平行な平面上から各再構成点の画素値を直接再構成する。再構成処理部は、FeldkampのコーンビームCT再構成アルゴリズムにおけるフィルター畳み込みの軸を二次元検出器と平行な平面の軸に変形した再構成アルゴリズムに基づき投影データを再構成する。これにより、変換に伴う隣接ピクセルの加算による補間処理を行わず直接再構成が可能となるため、処理負荷を抑えつつ画質を向上させた断層画像を得ることができる。また、重畳積分部115は、投影データに再構成フィルターを合成し、バックプロジェクション部117が逆投影処理する。これによりトモシンセシスで得られた投影画像をコーンビームCTの配置に幾何変換を行わずに、直接フィルタードバックプロジェクションを行い、断層像を求めることが可能になる。これにより従来の技術より、空間分解能が高く、処理時間の短い、トモシンセシスによる断層像を提供することができる。Feldkamp法など理論式を変換して得られる再構成式を用いるため、CT撮影による再構成に用いるrampフィルター、shoppLoganフィルターその他の再構成フィルターをそのまま利用することができる。もちろん、トモシンセシス撮影に伴う照射角不足を鑑みたトモシンセシス用のフィルターの適用も妨げない。