X線を照射するX線源と、前記X線を検出するX線検出器と、被検体を戴置する寝台と、前記X線源と前記X線検出器を搭載して前記被検体の周囲を回転する回転円盤と、前記X線検出器によって検出された前記X線が変換された実投影データを記憶する記憶装置と、を備えるX線CT装置であって、前記実投影データに対して第1ノイズ低減処理を施したデータを用いて第1画像を再構成し、更に、前記実投影データに対して前記第1ノイズ低減処理よりもノイズ低減度合が低い第2ノイズ低減処理を施したデータ、又は前記実投影データを用いて第2画像を再構成する画像再構成部と、前記第1画像及び前記第2画像の両方に含まれる所定の構造物について、前記第1画像から前記構造物の境界である第1境界を抽出し、更に、前記第2画像から前記構造物の境界である第2境界を抽出する境界抽出部と、前記第1境界及び前記第2境界の間の領域について、少なくとも前記第1画像及び前記第2画像を用いて合成処理を行い、合成画像を生成する画像合成部と、を備えるX線CT装置。前記第1画像に基づいて方向依存フィルタの方向を決定し、前記方向依存フィルタを前記第2画像に適用して平滑化処理を行い、第3画像を生成する画像平滑部、を更に備え、前記画像合成部は、更に、前記第3画像を用いて合成処理を行う請求項1に記載のX線CT装置。前記画像平滑部は、前記第1画像の注目画素に隣接する複数の方向の隣接画素に対して、それぞれ前記注目画素と前記隣接画素との差の絶対値を算出し、算出された値が最も小さい前記隣接画素の方向を前記方向依存フィルタの方向として決定する請求項2に記載のX線CT装置。前記画像再構成部は、前記第1ノイズ低減処理及び前記第2ノイズ低減処理として、前記X線検出器の検出器出力に応じたフォトン重み、及び罰則項を評価関数に含む逐次近似法によって前記実投影データを修正する請求項1に記載のX線CT装置。前記第1画像に基づいて、前記第1画像及び前記第2画像をどのように合成するかを示す第1合成マップを生成する第1合成マップ生成部と、少なくとも前記第2画像に基づいて前記第1合成マップを修正し、前記第1画像、前記第2画像、及び前記第3画像をどのように合成するかを示す第2合成マップを生成する第2合成マップ生成部と、を更に備え、前記画像合成部は、前記第2合成マップに従って合成処理を行う請求項2に記載のX線CT装置。前記第2合成マップによって示される前記第1画像を用いる領域と前記第3画像を用いる領域との境界周辺の画素について、前記合成画像の値を、前記第1画像の値と前記第3画像の値を重み付加算した値に変換する境界領域重み付加算部、を更に備える請求項5に記載のX線CT装置。前記第1画像と前記合成画像の全面において、操作者によって指定された重みの値によって重み付加算した画像を生成する画像重み付加算部、を更に備える請求項1に記載のX線CT装置。前記第1画像におけるエア領域及び前記エア領域に接する領域を示すエアマップを生成するエアマップ生成部、を更に備え、前記第1合成マップ生成部は、更に、前記エアマップに基づいて前記第1合成マップを生成し、前記第2合成マップ生成部は、更に、前記エアマップに基づいて前記第1合成マップを修正し、前記第2合成マップを生成する請求項5に記載のX線CT装置。X線を照射するX線源と、前記X線を検出するX線検出器と、被検体を戴置する寝台と、前記X線源と前記X線検出器を搭載して前記被検体の周囲を回転する回転円盤と、前記X線検出器によって検出された前記X線が変換された実投影データを記憶する記憶装置と、を備えるX線CT装置が実行する画像処理方法であって、前記実投影データに対して第1ノイズ低減処理を施したデータを用いて第1画像を再構成し、更に、前記実投影データに対して前記第1ノイズ低減処理よりもノイズ低減度合が低い第2ノイズ低減処理を施したデータ、又は前記実投影データを用いて第2画像を再構成するステップと、前記第1画像及び前記第2画像の両方に含まれる所定の構造物について、前記第1画像から前記構造物の境界である第1境界を抽出し、更に、前記第2画像から前記構造物の境界である第2境界を抽出するステップと、前記第1境界及び前記第2境界の間の領域について、少なくとも前記第1画像及び前記第2画像を用いて合成処理を行い、合成画像を生成するステップと、を含む画像処理方法。しかしながら、特許文献1に記載の方法では、強いストリークアーチファクトや強いノイズに対しては、構造物とノイズが同等に扱われ、低減しない場合があるという問題がある。本発明は、前述した問題点に鑑みてなされたものであり、その目的とすることは、修正された投影データを用いて再構成された画像SS再構成画像のノイズ低減効果を保持しつつ、構造物の輪郭に画質の隔たりや凹凸が発生させずに、ぼけてしまった構造物のコントラストを回復させることが可能なX線CT装置等を提供することである。前述した目的を達成するために第1の発明は、X線を照射するX線源と、前記X線を検出するX線検出器と、被検体を戴置する寝台と、前記X線源と前記X線検出器を搭載して前記被検体の周囲を回転する回転円盤と、前記X線検出器によって検出された前記X線が変換された実投影データを記憶する記憶装置と、を備えるX線CT装置であって、前記実投影データに対して第1ノイズ低減処理を施したデータを用いて第1画像を再構成し、更に、前記実投影データに対して前記第1ノイズ低減処理よりもノイズ低減度合が低い第2ノイズ低減処理を施したデータ、又は前記実投影データを用いて第2画像を再構成する画像再構成部と、前記第1画像及び前記第2画像の両方に含まれる所定の構造物について、前記第1画像から前記構造物の境界である第1境界を抽出し、更に、前記第2画像から前記構造物の境界である第2境界を抽出する境界抽出部と、前記第1境界及び前記第2境界の間の領域について、少なくとも前記第1画像及び前記第2画像を用いて合成処理を行い、合成画像を生成する画像合成部と、を備えるX線CT装置である。第2の発明は、X線を照射するX線源と、前記X線を検出するX線検出器と、被検体を戴置する寝台と、前記X線源と前記X線検出器を搭載して前記被検体の周囲を回転する回転円盤と、前記X線検出器によって検出された前記X線が変換された実投影データを記憶する記憶装置と、を備えるX線CT装置が実行する画像処理方法であって、前記実投影データに対して第1ノイズ低減処理を施したデータを用いて第1画像を再構成し、更に、前記実投影データに対して前記第1ノイズ低減処理よりもノイズ低減度合が低い第2ノイズ低減処理を施したデータ、又は前記実投影データを用いて第2画像を再構成するステップと、前記第1画像及び前記第2画像の両方に含まれる所定の構造物について、前記第1画像から前記構造物の境界である第1境界を抽出し、更に、前記第2画像から前記構造物の境界である第2境界を抽出するステップと、前記第1境界及び前記第2境界の間の領域について、少なくとも前記第1画像及び前記第2画像を用いて合成処理を行い、合成画像を生成するステップと、を含む画像処理方法である。本発明により、修正された投影データを用いて再構成された画像SS再構成画像のノイズ低減効果を保持しつつ、構造物の輪郭に画質の隔たりや凹凸が発生させずに、ぼけてしまった構造物のコントラストを回復させることができる。また、コントラストの回復度合を調整することができる。以下図面に基づいて、本発明の実施形態を詳細に説明する。最初に、図1、図2を参照しながら、X線CT装置1の構成を説明する。図1に示すように、X線CT装置1は、X線管11や検出器12が搭載されるスキャナ2、被検体10を載置する寝台4、検出器12から得られるデータの処理を行う演算装置5、マウス、トラックボール、キーボード、タッチパネルなどの入力装置6、及び再構成画像などを表示する表示装置7などを含む。操作者は、入力装置6を介して、撮影条件や再構成パラメータなどを入力する。撮影条件は、例えば、寝台送り速度、管電流、管電圧、スライス位置などである。また、再構成パラメータは、例えば、関心領域、再構成画像サイズ、逆投影位相幅、再構成フィルタ関数などである。図2に示すように、X線CT装置1は、大きく分けて、スキャナ2、操作ユニット3、寝台4から構成される。スキャナ2は、X線管11X線源、検出器12X線検出器、コリメータ13、駆動装置14、中央制御装置15、X線制御装置16、高電圧発生装置17、スキャナ制御装置18、寝台制御装置19、寝台移動計測装置20、コリメータ制御装置21、プリアンプ22、A/Dコンバータ23などから構成されている。中央制御装置15は、操作ユニット3における入力装置6から撮影条件や再構成パラメータを入力し、撮影に必要な制御信号を、コリメータ制御装置21、X線制御装置16、スキャナ制御装置18、寝台制御装置19に送信する。コリメータ制御装置21は、制御信号に基づいてコリメータ13の位置を制御する。撮影スタート信号を受けて撮影が開始されると、X線制御装置16は、制御信号に基づいて高電圧発生装置17を制御する。高電圧発生装置17は、X線管11に管電圧、管電流を印加する。X線管11では、印加された管電圧に応じたエネルギーの電子が陰極から放出され、放出された電子がターゲットに衝突することによって電子エネルギーに応じたエネルギーのX線が被検体10に照射される。また、スキャナ制御装置18は、制御信号に基づいて駆動装置14を制御する。駆動装置14は、X線管11、検出器12、プリアンプ22等が搭載されているガントリ部を被検体10の周りに周回させる。寝台制御装置19は、制御信号に基づいて寝台4を制御する。X線管11から照射されるX線は、コリメータ13によって照射領域が制限され、被検体10内の各組織においてX線減弱係数に応じて吸収減衰され、被検体10を通過し、X線管11に対向する位置に配置された検出器12によって検出される。検出器12は、2次元方向チャネル方向およびこれに直交する列方向に配置された複数の検出素子によって構成される。各検出素子によって受光したX線は、実投影データに変換される。すなわち、検出器12によって検出されるX線は、電流に変換され、プリアンプ22によって増幅され、A/Dコンバータ23によってデジタルデータに変換され、LOG変換され、キャリブレーションが行われて実投影データとして演算装置5に入力される。このとき、互いに対向するX線管11と検出器12が、被検体10の周囲を回転するので、実投影データは、回転方向の離散的なX線管位置および対応する検出器位置において収集される。各々のX線管位置における実投影データの取得単位が、ビューと呼ばれている。演算装置5は、再構成演算装置31、画像処理装置32等から構成される。また、入出力装置9は、入力装置6、表示装置7、記憶装置8記憶部等から構成される。再構成演算装置31は、実投影データを用いて画像再構成処理を行い、再構成画像を生成する。再構成演算装置31は、各ビューの実投影データに再構成フィルタを重畳してフィルタ補正投影データを生成し、フィルタ補正投影データに対して、ビュー方向に重み以降、ビュー方向重みという。を加重して逆投影処理を行うことによって、被検体10内部のX線減弱係数の分布図として非破壊的に断層像を画像化する。再構成演算装置31は、生成される再構成画像を記憶装置8に保存する。また、再構成演算装置31は、表示装置7にCT画像として表示する。あるいは、画像処理装置32が、記憶装置8に保存される再構成画像に対して画像処理を行い、表示装置7にCT画像として表示する。X線CT装置1は、2次元方向に検出素子が配列された検出器12を用いるマルチスライスCT、検出素子が1列すなわち1次元方向チャネル方向のみに配列された検出器12を用いるシングルスライスCTに大別される。マルチスライスCTでは、検出器12に合わせてX線源であるX線管11から円錐状、もしくは角錐状に広がるX線ビームが照射される。シングルスライスCTでは、X線管11から扇状に広がるX線ビームが照射される。通常、X線CT装置1による撮影では、ガントリ部が、寝台4に載置された被検体10の周りを周回しながら、X線の照射が行われる。はじめに、図3を参照しながら、本実施形態の概要を述べる。本実施形態では、演算装置5が、実投影データに対してノイズ低減処理を施したデータを用いてSS再構成画像50を再構成し、更に、実投影データに対して、SS再構成画像50の再構成時よりもノイズ低減度合が低いノイズ低減処理を施したデータ、又は実投影データを用いてNSS再構成画像60を再構成する。例えば、演算装置5は、逐次近似投影データフィルタを用いて実投影データを修正し、修正されたデータを用いて再構成された画像をSS再構成画像50とし、逐次近似投影データフィルタを用いた修正を行っていないデータ=実投影データを用いて再構成された画像をNSS再構成画像60とする。逐次近似投影データフィルタを用いた実投影データの修正処理としては、例えば、X線検出器の検出器出力に応じたフォトン重み、及び罰則項を評価関数に含む逐次近似法によって実投影データを修正する処理が挙げられる。例えば、非特許文献1に記載の次式を、以下の説明のように用いる。1、、i、、Iは、チャネルおよび列、ビューの通し番号である。piおよびyiは、i番目の検出素子における暫定投影データおよび実投影データである。p=p1、、pi、、pIは、暫定投影データのベクトルである。κiは、i番目の検出素子の近接素子の集合である。νikは、i番目の検出素子とk番目の検出素子の相関性を決定する定数である。非特許文献1では、チャネルおよび列、ビュー方向でνikを経験的に決定している。Ψpi-pkは、i番目の暫定投影データとk番目の暫定投影データ間のコントラストを変数とするポテンシャル関数である。非特許文献1では、ポテンシャル関数として2次関数を使用している。diは、i番目の実投影データと順投影データの差分値に加重される重み係数以降、フォトン重みよいう。である。X線CT装置の画像再構成では、i番目の検出素子における検出器出力を反映した値である。1式から種々の数値解析法により導出される更新式において、更新の過程における暫定投影データは、第一項により実投影データおよびフォトン重みによる制約を受けながら、第二項により近接の暫定投影データ間のコントラストに応じて平滑化される。よって、第一項と第二項のバランスを調整する役目を果たす1式の任意の定数βは、暫定投影データの平滑化の度合いを調整するパラメータである。逐次近似再構成および逐次近似投影データフィルタにおいて、フォトン重みは重要な役割を持つ。説明を分かり易くするため、ここではビュー方向およびチャネル方向の2次元として考える。水平方向のビューX線管位置が被写体の左または右になるビューでは、人体の体幅方向にX線が透過する為、鉛直方向のビューX線管位置が被写体の上又は下になるビューと比較して人体における透過長が長く、検出器出力が小さなチャネルが多い。これにより、断層像において主に水平方向にストリーク状のアーチファクトが生ずる。1式のフォトン重みは検出器出力に応じた値となるため、鉛直方向のビューのチャネルより水平方向のビューのチャネルで小さな値となる。これに応じて、鉛直方向のビューより水平方向のビューにおいて、1式の第一項による制約は相対的に弱くなり、第二項による平滑化効果を強く受ける。図3に示すSS再構成画像50及びNSS再構成画像60は、同一の構造物40例えば、骨などの部分を拡大した模式図である。斜線が付された領域が構造物40の本来の位置を示している。SS再構成画像50は、構造物40のエッジ輪郭がぼけてしまうので、エッジ抽出処理境界抽出処理を行うと、本来の構造物40よりも広がったエッジが抽出される。図3に示すSS再構成画像50では、実線にて図示されているSSエッジ51が、エッジ抽出処理によって抽出されたエッジを示している。一方、NSS再構成画像60は、構造物40のエッジ輪郭が鮮明ではあるが、ストリークアーチファクト62に代表されるノイズを含んでいる。NSS再構成画像60に対してエッジ抽出処理境界抽出処理を行うと、本来の構造物40に近いエッジが抽出される。図3に示すNSS再構成画像60では、実線にて図示されているNSSエッジ61が、エッジ抽出処理によって抽出されたエッジを示している。しかし、NSS再構成画像60ではストリークアーチファクト62が含まれており、構造物40のエッジ周辺の画質に問題がある。そこで、本実施形態では、演算装置5が、SS再構成画像50及びNSS再構成画像60の両方に含まれる構造物40について、SS再構成画像50から構造物40のエッジ境界を抽出し、更に、NSS再構成画像60から構造物40のエッジ境界を抽出する。そして、演算装置5が、SS再構成画像50から抽出されたSSエッジ51及びNSS再構成画像60から抽出されたNSSエッジ61の間の領域について、少なくともSS再構成画像50及びNSS再構成画像60の合成処理を行い、合成画像70を生成する。このように生成された合成画像70は、SS再構成画像50のノイズ低減効果を保持しつつ、構造物40の輪郭に画質の隔たりや凹凸が発生させずに、ぼけてしまった構造物40のコントラストが回復されている。尚、後述するように、SSエッジ51の外側、及びNSSエッジ61の内側については、演算装置5は、SS再構成画像50、NSS再構成画像60、或いはNSS再構成画像60を平滑化した画像、のいずれかのうち、適切な画像のCT値を用いて合成画像70を生成する。また、演算装置5は、SSエッジ51及びNSSエッジ61の間の領域の一部についても、NSS再構成画像60を平滑化した画像のCT値を用いて合成画像70を生成する。次に、図4を参照しながら、本実施形態における画像処理の特徴とその効果を説明する。本実施形態における画像処理は、図4に示すように、10のステップに分かれている。S401では、演算装置5は、円形などの再構成FOV範囲を示す再構成領域マップを生成する。S402では、演算装置5は、エアの領域を示すマップ=エアマップを生成する。S403では、演算装置5は、SS再構成画像のエッジ強度を示すエッジデータとそれらを2値化して得られた領域を表すエッジマップ=SSエッジマップを生成する。S404では、演算装置5は、SSエッジマップから、SS再構成画像にNSS再構成画像を合成する境界を示す合成マップ=SS合成マップを生成する。S405では、演算装置5は、NSS再構成画像のエッジ領域を表すエッジマップ=NSSエッジマップを生成する。S406では、演算装置5は、NSS再構成画像に対してエッジ保存型の平滑化処理が施された画像=SMSMoothed-NSS再構成画像を生成する。S407では、演算装置5は、NSSエッジマップを用いてSS合成マップを修正し、NSS合成マップを生成する。S408では、演算装置5は、NSS合成マップをもとに、SS再構成画像とNSS再構成画像とSM-NSS再構成画像が合成された画像=合成画像を生成する。S409では、演算装置5は、合成画像中にあるSS再構成画像を用いた領域とSM-NSS再構成画像を用いた領域の境界周辺の画素について、合成画像の値を、SS再構成画像の値とSM-NSS再構成画像の値を重み付加算した値に変換する。S410では、演算装置5は、SS再構成画像と合成画像の全面において、操作者によって指定された重みの値によって重み付加算した画像を生成する。本実施形態における画像処理の特徴とその効果は、少なくとも以下の5つが挙げられる。1つ目は、S404によって、SS再構成画像のエッジマップから、SS再構成画像にNSS再構成画像を合成する境界を示すSS合成マップを生成し、S407によって、NSSエッジマップを用いてSS合成マップの修正を行い、NSS合成マップを生成することである。これは、はじめにSS再構成画像を用いることで、ストリークアーチファクトの影響を少なくした境界抽出を行い、その境界の内側でNSS再構成画像を用いたぼけの影響の少ない境界を抽出することに相当する。そのため、ストリーク状のアーチファクトの影響が少なく、ぼけの影響の少ない境界を捕らえられるようになる。2つ目は、S404及びS407によって、エッジ強度だけでなく、CT値の参照とエア領域に隣接するかどうかの判定をすることである。これにより、ボケの多い骨などのハイコントラストな構造物のみにコントラスト回復を行い、逐次近似投影データフィルタによるノイズ低減効果を保持できる。3つ目は、S408によって、NSS再構成画像中の構造物の輪郭に対し、SS再構成画像の方向を用いて方向依存フィルタを施すことである。これにより、NSS再構成画像の構造物の輪郭のボケが少なく、その輪郭にある凹凸が軽減される。4つ目は、S409によって、合成画像中にあるSS再構成画像を用いた領域とSM-NSS再構成画像を用いた領域の境界周辺の値を、SS再構成画像の値とSM-NSS再構成画像の値を重み付加算した値に変換することである。これは、構造物の境界位置周辺にSS再構成画像とSM-NSS再構成画像の重み付き加算を行うことに相当する。そのため、合成画像の構造物の輪郭の隔たりや凹凸が軽減される。5つ目は、S410によって、SS再構成画像と合成画像の全面において、操作者によって指定された重みの値によって重み付加算した画像を生成することである。これによって、操作者が、合成画像に近い画像を生成するのか、或いはSS再構成画像に近い画像を生成するのかについて、調整できるようになる。すなわち、コントラスト回復の度合を調整できるようになる。以下では、図5~図19を参照しながら、S401~S410の詳細を説明する。<S401>S401の再構成範囲抽出処理では、演算装置5は、2次元の画像平面の中で、再構成画像が出力されているFOV範囲を示す再構成領域マップを生成する。通常は円形の範囲となる。<S402>S402のエア領域判定処理では、演算装置5は、SS再構成画像を入力し、エア領域となる画素、及びエア領域と構造物の境界となる中間的なCT値をもつ画素を含むエアマップを生成する。図5に示すように、エア領域判定処理のフローは、4つのステップに分けられる。S501では、演算装置5は、エアマップの全ての画素を0に初期化する。S502では、演算装置5は、エアマップの全ての画素を1つずつ走査していき、SS再構成画像におけるCT値が所定のパラメータDAIRCTVAL1以下の画素と同じ位置のエアマップの画素に1を設定する。S503では、演算装置5は、S502の結果、依然としてエアマップの値が0の画素に対して、エアマップの値が1の画素と隣接する位置、かつSS再構成画像における同じ位置の画素のCT値が所定のパラメータDAIRCTVAL2以下、の条件を満たすか否かを判定し、条件を満たすエアマップの画素に2を設定する。S504では、演算装置5は、S503の結果、依然としてエアマップの値が0の画素に対して、エアマップの値が2の画素と隣接する位置、かつSS再構成画像における同じ位置の画素のCT値が所定のパラメータDAIRCTVAL2以下、の条件を満たすか否かを判定し、条件を満たすエアマップの画素に3を設定する。図5に示す処理によって、確実に、エア領域となる画素、及びエア領域に隣接する画素を抽出し、エアマップを生成することができる。そして、このように生成されたエアマップによって、エア領域と接する体表や肺野内血管のエッジと、所望の構造物のエッジとを区別することができる。<S403>S403のSSエッジマップ生成処理では、演算装置5は、SS再構成画像と再構成領域マップとエアマップを入力し、エッジが高い領域を示すSSエッジマップを生成する。図6に示すように、SSエッジマップ生成処理のフローは、2つのステップに分けられる。S601のエッジ抽出処理では、演算装置5は、Sobelフィルタの値を算出し、これらをエッジ強度としたエッジデータを生成する。Sobelフィルタの値Gは、各画素に対して、図7に示す水平方向の係数行列を乗算した値Gh、及び垂直方向の係数行列を乗算した値Gvを用いて、以下の式となる。尚、S601の処理は、エッジの強度が算出される他の方法を用いてもよい。S602のエッジ領域2値化処理では、演算装置5は、S601によって生成されたエッジデータ中にあるエッジ強度が高い画素が2値化されたSSエッジマップを生成する。このエッジ領域2値化処理は、図8に示すように、3つのステップに分けられる。S801では、演算装置5は、再構成領域マップの値が1、かつエアマップの値が0の画素=エアではない画素の位置におけるエッジデータの値のヒストグラムを作成する。S802では、演算装置5は、S801によって作成されたヒストグラムにおいて、エッジ強度が高い方から頻度を積分し、全頻度×DEDGERATEとなる閾値thHistを算出する。DEDGERATEは所定のパラメータである。S803では、演算装置5は、S601によって生成されたエッジデータにおいて、S802によって算出された閾値thHist以上となるSSエッジマップの画素に1を設定する。例えば、DEDGERATEが0.4であれば、SSエッジマップにおいて、エッジ強度の値が高い方から全体の40に入る画素にエッジ領域であることを示す1が設定される。<S404>S404のSS再構成画像によるSS合成マップ生成処理では、演算装置5は、SS再構成画像とSSエッジデータとSSエッジマップとエアマップを入力し、SS再構成画像とNSS再構成画像をどの位置に合成するかを表すSS合成マップを生成する。S404の処理では、演算装置5は、1画素毎に、図9に示す条件判定を行う。はじめに、演算装置5は、SS再構成画像のCT値が所定のパラメータCTVTH2以上であるかを判定しS901、CTVTH2以上であればS901のYes、SS合成マップに1を設定するS902。このように判定された画素は、CT値が高い画像である。そして、このような画素については、基本的にNSS再構成画像のCT値を用いて合成画像が生成されることになる。例えば、図3の例における構造物40が、CT値が高い骨のような場合、演算装置5は、図3に示す合成画像70の斜線領域については、基本的には、NSS再構成画像のCT値を用いて合成画像を生成する。但し、後述するS407の処理によって再判定が行われ、NSS再構成画像のCT値を用いない場合もある。SS再構成画像のCT値がCTVTH2よりも小さい場合S901のNo、演算装置5は、SSエッジマップの値が1であるかを判定しS903、SSエッジマップの値が1でない場合S903のNo、SS合成マップに5を設定するS904。このように判定された画素は、エッジではない画素である。SSエッジマップの値が1である場合S903のYes、演算装置5は、SS再構成画像のCT値が所定のパラメータCTVTH以上であるかを判定しS905、CTVTH以上であればS905のYes、SS合成マップに1を設定するS906。このように判定された画素は、所望の構造物、或いは所望の構造物のエッジに相当する画素である。SS再構成画像のCT値がCTVTHより小さい場合S905のNo、演算装置5は、エアマップの値が0より大きいかを判定するS907。エアマップの値が0の場合、すなわちエア領域ではない場合S907のNo、演算装置5は、SSエッジデータの値が所定のパラメータEDVTH以上であるかを判定しS908、EDVTH以上であればS907のYes、SS合成マップに2を設定するS909。このように判定された画素は、CT値が高くなく、かつエッジの強度が高い画素である。そして、このような画素については、基本的にNSS再構成画像のCT値を用いて合成画像が生成されることになる。但し、後述するS407の処理によって再判定が行われ、NSS再構成画像のCT値を用いない場合もある。SSエッジデータの値がEDVTHより小さい場合S908のNo、演算装置5は、SS合成マップに3を設定するS910。一方、エアマップの値が1以上の場合、すなわちエア領域と判定された画素の場合S907のYes、演算装置5は、SSエッジデータの値がEDVTH以上であるかを判定しS911、EDVTH以上である場合S911のYes、SS合成マップに4を設定するS912。このように判定された画素は、基本的には、体表や肺野内血管のエッジに相当する画素である。そして、このような画素については、エッジを保存するために、NSS再構成画像のCT値を用いて合成画像が生成されることになる。SSエッジデータの値がEDVTHより小さい場合S911のNo、演算装置5は、SS合成マップに2を設定するS913。このように判定された画素は、エア領域の境界においてエッジ強度が低い画素である。そして、このような画素については、後述するS407において再判定が行われる。<S405>S405のNSSエッジマップ生成処理では、演算装置5は、S403のSSエッジマップ生成処理と同じ処理を行い、NSSエッジデータとNSSエッジマップを生成する。<S406>S406のNSS再構成画像のエッジ保存型平滑化処理では、演算装置5は、図10に示す条件判定処理により、特定の領域においてSS再構成画像に基づいて方向が決定された方向依存フィルタが施された画像=SMNSS再構成画像を生成する。はじめに、図10に示す条件判定処理について説明する。次に、図11に示す方向依存フィルタ処理について説明する。図10に示す条件判定処理では、演算装置5は、エア領域と接する位置にあるか否かを判定するS1001。エア領域に接する位置にある場合S1001のYes、演算装置5は、何も処理を施さないS1002。これは、エア領域に接する画素については、体表や肺野内血管に接する画素であるため、平滑化を行わない。エア領域と接する位置にない場合S1001のNo、演算装置5は、NSS再構成画像のCT値が所定のパラメータDSMCTVALより大きい、或いは、NSS再構成画像のエッジ強度が所定のパラメータDSMEVALより大きいかを判定するS1003。NSS再構成画像のCT値がDSMCTVALよりも大きい、或いはNSS再構成画像のエッジ強度が所定のパラメータDSMEVALよりも大きい場合S1003のYes、演算装置5は、SS再構成画像のCT値がDSMCTVALよりも大きい、かつ、NSS再構成画像のエッジ強度がDSMEVALよりも小さいかを判定するS1004。SS再構成画像のCT値がDSMCTVALよりも大きく、かつ、NSS再構成画像のエッジ強度がDSMEVALよりも小さい場合S1004のYes、演算装置5は、何も処理を施さないS1005。このように判定された画素は、骨の内部や臓器の輪郭の画素であるため、平滑化を行わない。SS再構成画像の値がDSMCTVAL以下、或いはNSS再構成画像のエッジ強度がDSMEVAL以上である場合S1004のNo、演算装置5は、SS再構成画像に基づいて方向が決定された方向依存フィルタ処理を施すS1006。このように判定された画素は、骨の輪郭の画素などであるため、そのままNSS再構成画像のCT値を用いると、ストリークアーチファクトの影響によって、骨などの輪郭に凹凸が発生してしまう可能性がある。従って、方向依存フィルタ処理を施す。一方、NSS再構成画像のCT値がDSMCTVAL以下、かつNSS再構成画像のエッジ強度がDSMEVAL以下の場合S1003のNo、演算装置5は、NSS再構成画像のCT値がDSMCTVALMTHより小さいかを判定するS1007。NSS再構成画像のCT値がDSMCTVALMTHよりも小さい場合S1007のYes、演算装置5は、何も処理を施さないS1008。このように判定された画素は、主にエア領域の画素であるため、平滑化を行わない。NSS再構成画像のCT値がDSMCTVALMTH以上である場合S1007のNo、演算装置5は、SS再構成画像に基づいて方向が決定された方向依存フィルタ処理を施すS1009。このように判定された画素は後述する重み付加算領域となるため、ノイズによる凹凸を軽減しておく必要がある。従って、方向依存フィルタ処理を施す。図11に示す方向依存フィルタ処理では、演算装置5は、ストリークアーチファクトが軽減されたSS再構成画像のエッジの方向を利用して、NSS再構成画像の平滑化を行う。方向依存フィルタ処理のフローは、図11に示すように、2つのステップからなる。S1101のSS再構成画像を用いた平滑化方向の決定処理では、演算装置5は、例えば、左右方向、上下方向、右上がり方向、及び左上がり方向の4方向の隣接画素それぞれ、2つずつの画素との差の絶対値を算出する。そして、S1102において、演算装置5は、4つの値の中で最も小さい値となる方向の画素を用いて平滑化を行う。具体的な算出式は、以下の通りである。まず、演算装置5は、式3~式6によって、注目画素と各方向の隣接画素との差の絶対値を算出する。但し、X軸の負方向が左方向、X軸の正方向が右方向、Y軸の負方向が上方向、Y軸の正方向が下方向とし、Pif:SS画像のCT値、F1if:左右方向の隣接画素との差の絶対値、F2if:上下方向の隣接画素との差の絶対値、F3if:右上がり方向の隣接画素との差の絶対値、F4if:左上がり方向の隣接画素との差の絶対値、である。次に、MinF1ifであれば、演算装置5は、式7に示す左右方向の隣接画素を用いた方向依存フィルタによって、NSS画像の平滑化処理を行う。但し、Pf:NSS画像のCT値、Psm:平滑化後のCT値=SM-NSS再構成画像のCT値、である。また、MinF2ifであれば、演算装置5は、式8に示す上下方向の隣接画素を用いた方向依存フィルタによって、NSS画像の平滑化処理を行う。また、MinF3ifであれば、演算装置5は、式9に示す右上がり方向の隣接画素を用いた方向依存フィルタによって、NSS画像の平滑化処理を行う。また、MinF4ifであれば、演算装置5は、式10に示す左上がり方向の隣接画素を用いた方向依存フィルタによって、NSS画像の平滑化処理を行う。前述の式7~式10では、単純平均によるフィルタとしたが、加重平均によるフィルタとしても良い。このように、SS再構成画像の情報に基づいて決定した方向を利用した方向依存フィルタを用いて、NSS再構成画像を平滑化することによって、NSS再構成画像の構造物の輪郭のボケが少なく、その輪郭にある凹凸例えば、ストリークアーチファクトなどが軽減された、SM-NSS再構成画像を生成することができる。但し、本実施の形態における平滑化処理はこれに限られるものではなく、既存の平滑化フィルタを用いてNSS再構成画像を平滑化しても良い。<S407>S407のNSSエッジマップによるSS合成マップ修正処理では、演算装置5は、図12の条件判定に基づいて、S404によって生成されたSS合成マップを修正し、NSS合成マップを生成する。S407の処理では、演算装置5は、1画素毎に、図12に示す条件判定を行う。はじめに、演算装置5は、SS合成マップが4であるかを判定するS1201。SS合成マップが4である場合S1201のYes、演算装置5は、NSS合成マップに4をそのまま設定する。SS合成マップが4でない場合S1201のNo、演算装置5は、SS合成マップが2以下であるかを判定するS1203。SS合成マップが2より大きい場合S1203のNo、演算装置5は、NSS合成マップに5を設定する。SS合成マップが2以下の場合S1203のNo、演算装置5は、SS再構成画像のCT値がCTVTH2NSS以上であるかを判定しS1205、CTVTH2NSS以上であればS1205のYes、NSS合成マップに1を設定するS1206。このように判定された画素は、CT値が高い画素である。SS再構成画像のCT値がCTVTH2NSSよりも小さい場合S1205のNo、演算装置5は、NSSエッジマップの値が1であるかを判定しS1207、1でない場合S1207のNo、NSS合成マップに5を設定するS1208。このように判定された画素は、エッジではない画素である。NSSエッジマップの値が1である場合S1207のYes、演算装置5は、SMNSS再構成画像のCT値がCTVTHNSS以上であるかを判定しS1209、CTVTHNSS以上であればS1209のYes、NSS合成マップに1を設定するS1210。SMNSS再構成画像のCT値がCTVTHNSSより小さい場合S1209のNo、演算装置5は、NSSエッジデータの値がEDVTHNSS以上であるかを判定しS1211、EDVTHNSSより小さい場合S1211のNo、NSS合成マップに3を設定するS1212。NSSエッジデータの値がEDVTHNSS以上の場合S1211のYes、演算装置5は、エアマップの値が0より大きいかを判定するS1213。エアマップの値が0以下の場合、すなわちエア領域に接しない場合S1213のNo、演算装置5は、NSS合成マップに2を設定するS1214。エアマップの値が1以上の場合S1213のYes、演算装置5は、NSS合成マップに4を設定するS1215。このように判定された画素は、体表や肺野内血管のエッジに相当する画素である。<S408>S408の合成処理では、演算装置5は、NSS合成マップに従って、SS再構成画像、NSS再構成画像、及びSM-NSS再構成画像の中から適切な画像を用いて合成画像を生成する。図12の判定結果と合成処理に用いられる画像の組合せを纏めると、以下の通りである。1NSS合成マップの値が1CT値が高い画素=所望の構造物に相当する画素SM-NSS再構成画像2NSS合成マップの値が2所望の構造物のエッジに相当する画素SM-NSS再構成画像3NSS合成マップの値が3所望の構造物でもなく、エッジでもない画素SS再構成画像4NSS合成マップの値が4体表や肺野内血管のエッジに相当する画素NSS再構成画像5NSS合成マップの値が5所望の構造物でもなく、エッジでもない画素SS再構成画像<S409>S409の境界領域重み付加算処理では、演算装置5は、構造物の外側にSS再構成画像とNSS再構成画像の値の重み付加算によって算出する領域重み付加算領域を設けることによって、境界の隔たりと構造物外側のぼけを軽減する。演算装置5は、構造物の境界からスタートし、NSS再構成画像とSS再構成画像のCT値差が閾値よりも小さくなる境界=CT値差境界を抽出する。そして、演算装置5は、合成画像における構造物の境界からCT値差境界までのCT値を、構造物の境界からの距離とCT値差境界からの距離を用いて、NSS再構成画像のCT値とSS再構成画像のCT値の重み付け加算値とする。S409の処理は、図13に示すように、4つのステップに分けられる。S1301では、演算装置5は、SS再構成画像とNSS再構成画像の差分の絶対値を算出し、その値がDIFMAPTH以上の画素に1を設定し、それ以外の画素に0を設定することによって、差分大マップを生成する。S1302では、演算装置5は、SMNSS再構成画像領域との境界からの距離を表す距離画像を生成する。合成画像の例を図14、差分大マップの例を図15に示し、これらを入力とした場合のS1302におけるSMNSS再構成画像領域との境界からの距離画像の生成処理過程を図16に示す。プロセス1では、演算装置5は、SMNSS再構成画像利用領域の輪郭となる画素に-1、SMNSS再構成画像利用領域の輪郭よりも内側の画素に-2、NSS再構成画像利用領域の画素に-10を設定する。プロセス2では、演算装置5は、SS再構成画像領域においてSMNSS再構成画像領域と隣接する画素に1を設定する。これが重み付加算領域とSMNSS再構成画像領域との境界であり、1は境界からの距離となる。尚、図16の例では、隣接する画素は、注目画素の上、下、右、左、右上、右下、左上、左下の8つの画素である。プロセス3では、演算装置5は、1が設定された画素から差分大マップ上でフラグが立っている領域の範囲内、かつ所定の距離DMAXDIST以内を表す距離画像を生成する。図16の例では、DMAXDISTを6としている。距離画像の生成の方法は様々あるが、例えば次のように行う。はじめに、演算装置5は、現時点で0が設定されている画素のうち、1が設定されている画素と隣接する画素に2を設定する。次に、演算装置5は、現時点で0が設定されている画素のうち、2が設定された画素に隣接する画素に3を設定する。演算装置5は、このような処理を、設定される値がDMAXDISTになるまで繰り返す。このようにして得られたプロセス3の画像が、SMNSS再構成画像領域との境界からの距離画像となる。S1303では、演算装置5は、重み付加算領域におけるSS再構成画像領域との境界からの距離を表す距離画像を生成する。この処理過程を図17及び図19に示す。プロセス1では、演算装置5は、SMNSS再構成画像領域との境界からの距離画像において、距離値が割り当てられた画素と割り当てられていない画素の境界に1を設定する。つまり、1が設定された画素群が、重み付加算領域とSS再構成画像領域の境界となる。プロセス2では、演算装置5は、SS再構成画像領域との境界からSMNSS再構成画像領域の境界までの距離を表す距離画像を生成する。演算装置5は、SMNSS再構成画像領域との境界からの距離画像の1以上の画素に対して、基本的にはSMNSS再構成画像領域との境界からの距離画像生成処理と同じ処理を行う。これによって、SS再構成画像領域との境界からの距離画像が生成される。但し、プロセス2におけるSS再構成画像領域との境界からの距離画像生成処理において、前述の算出処理では距離が算出できない場合がある。例えば、合成画像、差分大マップ、及びSMNSS再構成画像領域との境界からの距離画像が図18に示す例の場合、SMNSS再構成画像領域との境界からの距離画像において、距離値が割り当てられた画素と割り当てられていない画素の境界、すなわち重み付加算領域とSS再構成画像領域の境界がない。このような場合、前述の算出処理では、SS再構成画像領域との境界からの距離を算出できない。これに対して、図19に示すプロセス3では、演算装置5は、SMNSS再構成画像領域との境界からの距離画像において、隣接する8近傍のすべての画素の値以上の画素値をもつ画素図19における斜線部の画素を抽出し、SS再構成画像との境界からの距離画像中の該当画素に1を設定する。プロセス4では、演算装置5は、SMNSS再構成画像境界からの距離画像において1以上の画素に対し、SMNSS再構成画像領域との境界からの距離画像生成と同じ処理を行うことによって、SS再構成画像領域からの距離画像を生成する。このようにして得られたプロセス4の画像が、SS再構成画像領域との境界からの距離画像となる。S1304では、演算装置5は、SMNSS再構成画像領域との境界からの距離画像の値とSS再構成画像領域との境界からの距離画像の値を重み係数に用いて、SMNSS再構成画像の値とSS再構成画像の値との重み付け和を算出する。算出式の例を以下に示す。但し、ds:SM-NSS画像領域との境界からの距離、di:SS画像領域との境界からの距離、Ps:SM-NSS画像のCT値、Pi:SS画像のCT値、Pm:重み付加算した画素値、である。<S410>S410の全面重み付加算処理では、演算装置5は、SS再構成画像と合成画像の全面において、操作者によって指定された重みの値によって重み付加算した画像を生成する。算出式の例を以下に示す。但し、w:操作者が指定した重み、Pa:合成画像の画素値、Pb:SS画像の画素値、Pc:重み付加算した画素値、である。重みwは、00~10までの値である。wを大きい値とするほど合成画像に近い画像が生成され、小さい値とするほどSS再構成画像に近い画像が生成される。これによって、コントラスト回復の度合を調整できるようになる。なお、本実施形態では逐次近似投影データフィルタを適用した画像と適用していない画像を用いたが、非特許文献2に提案されているような逐次近似画像再構成によって生成された画像をSS再構成画像とし、逐次近似再構成を行う前の画像や従来の再構成法にて生成された画像をNSS再構成画像として、同じように適用することが可能である。以上、添付図面を参照しながら、本発明に係るX線CT装置等の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。