"ヨウ素を含む造影剤が投与された被検体のX線CTデュアルエネルギー撮影により得られたデータに基づいて、前記被検体を水およびヨウ素の合成で表すモデルに基づく、前記被検体における仮想的な水の密度分布を表す水密度画像と、前記被検体における仮想的なヨウ素の密度分布を表すヨウ素密度画像とを生成する第1の生成手段と、前記水密度画像およびヨウ素密度画像に基づいて、前記被検体の非造影撮影により得られる画像に相当する仮想非造影画像を生成する第2の生成手段とを備えた画像生成装置。前記第2の生成手段は、前記仮想非造影画像の少なくとも一部を、仮想非造影画像=β×水密度画像α×ヨウ素密度画像-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式1なる数式1にしたがって生成する請求項1に記載の画像生成装置。前記第2の生成手段は、前記仮想非造影画像の少なくとも一部を、仮想非造影画像=β×{水密度画像α×ヨウ素密度画像α2×ヨウ素密度画像}-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式2なる数式2にしたがって生成する請求項1に記載の画像生成装置。前記仮想非造影画像の少なくとも一部は、前記被検体の軟部に相当する部分である請求項2または請求項3に記載の画像生成装置。前記第2の生成手段は、前記軟部に相当する部分を、前記データを用いて生成される画像の画素値に基づいて特定する請求項4に記載の画像生成装置。前記第2の生成手段は、前記被検体の軟部とは異なる物質に相当する部分を、前記データを基に、所定の実効X線エネルギーのX線を用いたX線CT撮影により得られる画像に相当する画像として生成する請求項4または請求項5に記載の画像生成装置。前記第2の生成手段は、前記被検体の軟部とは異なる物質に相当する部分を、仮想非造影画像=b×水密度画像f被検体中のヨウ素濃度に係る指標但し、1.0≧b≧0.9、f:前記指標をパラメータとする関数...数式3なる数式3にしたがって生成する請求項4または請求項5に記載の画像生成装置。前記第2の生成手段は、前記被検体の軟部とは異なる物質に相当する部分を、仮想非造影画像=b×水密度画像c但し、1.1≧b≧0.9、c:定数...数式4なる数式4にしたがって生成する請求項4または請求項5に記載の画像生成装置。前記数式1において、1.05≧β≧0.95、0.29≧α≧0.25、1050≧γ≧950なる条件を満たす請求項2から請求項8のいずれか一項に記載の画像生成装置。前記数式1におけるαは、操作者により指定される請求項2から請求項9のいずれか一項に記載の画像生成装置。前記第1の生成手段は、前記デュアルエネルギー撮影により得られた、第1のエネルギーのX線による第1の投影データと、前記第1のエネルギーとは異なる第2のエネルギーのX線による第2の投影データとを、ヨウ素の成分が消去されるように加重減算処理して得られたデータを画像再構成処理することにより、前記水密度画像を生成し、前記第1の投影データと前記第2の投影データとを、水の成分が消去されるように加重減算処理して得られたデータを画像再構成処理することにより、前記ヨウ素密度画像を生成する請求項1から請求項10のいずれか一項に記載の画像生成装置。前記データに基づいて、前記造影剤が投与された被検体を所望の実効X線エネルギーのX線を用いてX線CT撮影したときに得られる画像に相当する画像を生成する第3の生成手段をさらに備えた請求項1から請求項11のいずれか一項に記載の画像生成装置。請求項1から請求項12のいずれか一項に記載の画像生成装置を備えたX線CT装置。コンピュータを、請求項1から請求項12のいずれか一項に記載の画像生成装置として機能させるためのプログラム。単純画像と造影画像を取得するには、通常、被検体の単純撮影と造影撮影とを別々に行う必要があるため、被検体の被曝が一般の検査と比較して大きくなる。このような事情により、造影検査における被検体の低被曝化が望まれている。第1の観点の発明は、ヨウ素を含む造影剤が投与された被検体のX線CTデュアルエネルギー撮影により得られたデータdataに基づいて、前記被検体を水およびヨウ素の合成で表すモデルに基づく、前記被検体における仮想的な水の密度分布を表す水密度画像と、前記被検体における仮想的なヨウ素の密度分布を表すヨウ素密度画像とを生成する第1の生成手段と、前記水密度画像およびヨウ素密度画像に基づいて、前記被検体の非造影撮影により得られる画像に相当する仮想非造影画像を生成する第2の生成手段とを備えた画像生成装置を提供する。第2の観点の発明は、前記第2の生成手段が、前記仮想非造影画像の少なくとも一部を、仮想非造影画像=β×水密度画像α×ヨウ素密度画像-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式1なる数式1にしたがって生成する上記第1の観点の画像生成装置を提供する。第3の観点の発明は、前記第2の生成手段が、前記仮想非造影画像の少なくとも一部を、仮想非造影画像=β×{水密度画像α×ヨウ素密度画像α2×ヨウ素密度画像}-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式2なる数式2にしたがって生成する上記第1の観点の画像生成装置を提供する。第4の観点の発明は、前記仮想非造影画像の少なくとも一部が、前記被検体の軟部に相当する部分である上記第2の観点または第3の観点の画像生成装置を提供する。第5の観点の発明は、前記第2の生成手段が、前記軟部に相当する部分を、前記データを用いて生成される画像の画素値に基づいて特定する上記第4の観点の画像生成装置を提供する。第6の観点の発明は、前記第2の生成手段が、前記被検体の軟部とは異なる物質に相当する部分を、前記データを基に、所定の実効X線エネルギーのX線を用いたX線CT撮影により得られる画像に相当する画像として生成する上記第4の観点または第5の観点の画像生成装置を提供する。第7の観点の発明は、前記第2の生成手段が、前記被検体の軟部とは異なる物質に相当する部分を、仮想非造影画像=b×水密度画像f被検体中のヨウ素濃度に係る指標但し、1.0≧b≧0.9、f:前記指標をパラメータとする関数...数式3なる数式3にしたがって生成する上記第4の観点または第5の観点の画像生成装置を提供する。第8の観点の発明は、前記第2の生成手段が、前記被検体の軟部とは異なる物質に相当する部分を、仮想非造影画像=b×水密度画像c但し、1.1≧b≧0.9、c:定数...数式4なる数式4にしたがって生成する上記第4の観点または第5の観点の画像生成装置を提供する。第9の観点の発明は、前記数式1において、1.05≧β≧0.95、0.29≧α≧0.25、1050≧γ≧950なる条件を満たす上記第2の観点から第8の観点のいずれか一つの観点の画像生成装置を提供する。第10の観点の発明は、前記数式1におけるαが、操作者により指定される上記第2の観点から第9の観点のいずれか一つの観点の画像生成装置を提供する。第11の観点の発明は、前記第1の生成手段が、前記デュアルエネルギー撮影により得られた、第1のエネルギーのX線による第1の投影データと、前記第1のエネルギーとは異なる第2のエネルギーのX線による第2の投影データとを、ヨウ素の成分が消去されるように加重減算処理して得られたデータを画像再構成処理することにより、前記水密度画像を生成し、前記第1の投影データと前記第2の投影データとを、水の成分が消去されるように加重減算処理して得られたデータを画像再構成処理することにより、前記ヨウ素密度画像を生成する上記第1の観点から第10の観点のいずれか一つの観点の画像生成装置を提供する。第12の観点の発明は、前記データに基づいて、前記造影剤が投与された被検体を所望の実効X線エネルギーのX線を用いてX線CT撮影したときに得られる画像に相当する画像を生成する第3の生成手段をさらに備えた上記第1の観点から第11の観点のいずれか一つの観点の画像生成装置を提供する。第13の観点の発明は、上記第1の観点から第12の観点のいずれか一つの観点の画像生成装置を備えたX線CT装置を提供する。第14の観点の発明は、コンピュータcomputerを、上記第1の観点から第12の観点のいずれか一つの観点の画像生成装置として機能させるためのプログラムを提供する。なお、ここで言う密度とは、実際の密度ではなく仮想的な密度であり、所定の2種類の物質を基底物質とし、ある物質のX線吸収係数をその2種類の基底物質のX線吸収係数の加重加算値で表す場合に、その2種類の基底物質のX線吸収係数に乗ずる加重係数に相当するものである。したがって、例えば、水およびヨウ素を1組の基底物質としたときの水密度画像は、撮影対象の画像における各画素が表す物質のX線吸収係数を、水のX線吸収係数とヨウ素のX線吸収係数との加重加算値で表す場合における、水のX線吸収係数に乗ずる加重係数に相当する値の分布図となる。上記観点の発明によれば、ヨウ素を含む造影剤が投与された被検体のX線CTデュアルエネルギー撮影により、仮想非造影画像を高精度に得ることができる。その結果、造影画像に相当する画像と、非造影画像に相当する画像の両方を、1回のデュアルエネルギー撮影によって得ることができ、造影撮影および非造影撮影をそれぞれ行う必要がなく、造影検査における低被曝化を実現できる。以下、発明の実施形態について説明する。第一実施形態図1は、第一実施形態によるX線CT装置の構成を示す図である。X線CT装置100は、操作コンソールconsole1と、撮影テーブルtable10と、走査ガントリgantry20とを具備している。操作コンソール1は、操作者からの入力を受け付ける入力装置2と、被検体40をスキャンscanするための各部の制御や各種演算などを行う中央処理装置3と、走査ガントリ20で取得したデータを収集するデータ収集バッファbuffer5と、画像を表示するモニタmonitor6と、プログラムやデータなどを記憶する記憶装置7とを具備している。撮影テーブル10は、被検体40を載せて走査ガントリ20の開口部に搬入搬出するクレードルcradle12を具備している。クレードル12は、撮影テーブル10に内蔵するモータmotorで昇降および水平直線移動される。走査ガントリ20は、回転部15と、回転部15を回転可能に支持する支持部16とを有する。回転部15には、X線管21と、X線管21を制御するX線コントローラcontroller22と、X線管21から発生したX線81をコリメートcollimateして整形するコリメータcollimator23と、X線管21から照射され、被検体40を透過したX線81を検出するX線検出器24と、X線検出器24の出力を投影データに変換して収集するデータ収集装置DAS;DataAcquisitionSystem25と、X線コントローラ22,コリメータ23,DAS25の制御を行う回転部コントローラ26とが搭載される。支持部16は、制御信号などを操作コンソール1や撮影テーブル10と通信する制御コントローラ29を具備する。回転部15と支持部16とは、スリップリングslipring30を介して電気的に接続されている。これより、第一実施形態に係るX線CT装置における処理の流れについて説明する。図2は、第一実施形態に係るX線CT装置における処理の流れを示す図である。ステップstepS1では、撮影対象である被検体40をデュアルエネルギー撮影して、第1X線管電圧による複数ビューviewの投影データである第1管電圧投影データPv1と、第2X線管電圧による複数ビューの投影データである第2管電圧投影データPv2とを収集する。例えば、X線管電圧を1または数ビュー単位で第1X線管電圧と第2X線管電圧とに交互に切り換えながら被検体40のスキャンscanを行い、180°X線ビームbeamのファンfan角α°分以上の投影データを収集する。この場合、欠落したビューの投影データは、その欠落したビューに近接するビューの投影データを加重加算処理して求める。また例えば、まずX線管電圧を第1X線管電圧にしてスキャンを行って、180°ファン角α°分以上の投影データを収集し、次にX線管電圧を第2X線管電圧にしてスキャンを行って、180°ファン角α°分以上の投影データを収集する。なお、第1および第2X線管電圧は、例えば、80kVpと140kVpである。ステップS2では、被検体40を2種類の基底物質の合成によって表すという所定のモデルに基づき、水およびヨウ素を基底物質として、第1および第2管電圧投影データPv1,Pv2を基に、被検体40における仮想的な水の密度分布を表す水密度画像Gwと、被検体40における仮想的なヨウ素の密度分布を表すヨウ素密度画像Gioとを再構成する。図3および図4は、水およびヨウ素を基底物質として、第1および第2管電圧投影データを基に、ヨウ素密度画像および水密度画像を再構成する方法の一例を概念的に示す図である。例えば、図3に示すように、第1管電圧投影データPv1と第2管電圧投影データPv2とをヨウ素の成分が消去されるように加重減算処理して、水の成分を表す第1中間投影データPwaを生成する。また、第1管電圧投影データPv1と第2管電圧投影データPv2とを水の成分が消去されるように加重減算処理して、ヨウ素の成分を表す第2中間投影データPioを生成する。次に、第1中間投影データPwaを画像再構成処理して水密度画像Gwを再構成する。また、第2中間投影データPioを逆投影処理等により画像再構成処理してヨウ素密度画像Gioを再構成する。また例えば、図4に示すように、第1管電圧投影データPv1を画像再構成処理して第1管電圧画像Gv1を再構成する。また、第2管電圧投影データPv2を画像再構成処理して第2管電圧画像Gv2を再構成する。次に、第1管電圧画像Gv1と第2管電圧画像Gv2とをヨウ素の成分が消去されるように加重減算処理して、水密度画像Gwを再構成する。また、第1管電圧画像Gv1と第2管電圧画像Gv2とを水の成分が消去されるように加重減算処理して、ヨウ素密度画像Gioを再構成する。ステップS3では、水密度画像Gwおよびヨウ素密度画像Gioを基に、造影剤が投与された被検体40を所望の実効X線エネルギーE1のX線を用いてX線CT撮影したときに得られる画像に相当するモノクロマチック画像Gmを再構成する。このモノクロマチック画像Gmは、被検体40の造影画像に相当する画像である。実効X線エネルギーE1は、例えば、120keVである。図5は、水密度画像およびヨウ素密度画像を基にモノクロマチック画像を再構成する方法の一例を概念的に示す図である。例えば、図5に示すように、水密度画像Gwおよびヨウ素密度画像Gioを、実効X線エネルギーE1に対応した所定の加重係数により加重加算処理して、モノクロマチック画像Gmを再構成する。なお、各物質密度画像およびモノクロマチック画像の生成方法については、特許文献:特開2010-172590号公報や、米国特許文献:米国特許7724865号公報US2009/0052612等を参照されたい。ステップS4では、水密度画像Gwおよびヨウ素密度画像Gioを基に、造影剤が投与される前の被検体40をX線CT撮影したときに得られる画像に相当する仮想非造影画像Gncを、次の数式1にしたがって生成する。Gnc=βGwαGio-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式1すなわち、仮想非造影画像Gncの各画素値を、水密度画像Gwにおける対応画素の画素値に係数βを乗じた値と、ヨウ素密度画像Gioにおける対応画素の画素値に係数αを乗じた値と、オフセットoffset成分である-γとを加算して算出する。ここで、数式1の意味について説明する。造影剤が投与された被検体40のデュアルエネルギー撮影により得られたデータを基に水密度画像Gwおよびヨウ素密度画像Gioを生成すると、被検体40中のヨウ素の密度は、ヨウ素密度画像Gioに反映され、水密度画像Gioにはほとんど反映されないはずである。本実施形態で求めようとしている仮想非造影画像Gncは、被検体40を造影剤なしでX線CT撮影したときに得られる、いわゆる単純画像に相当する画像である。したがって、仮想非造影画像Gncは、基本的には、被検体40中のヨウ素の密度がほとんど反映されないはずの水密度画像Gwに近い画像であると考えることができる。ここで、濃度が異なる複数種類の造影剤ファントムphantomをデュアルエネルギー撮影して得られた、ヨウ素密度画像における画素値と、水密度画像における画素値との関係を図6に示す。図6のグラフから理解されるように、造影剤の濃度が変化すると、ヨウ素密度画像の画素値は当然変化するが、造影剤の濃度によって水密度画像の画素値も変化していることが分かる。これは、被検体40に投与する造影剤の濃度や量、造影剤の染まり具合被検体の個体差、造影剤投与からの撮影タイミングtimingなどに起因する、被検体40中のヨウ素の濃度によっては、ヨウ素密度画像Gioの画素値と共に、水密度画像Gwの画素値も変化することを意味している。つまり、水密度画像Gwそのままでは、単純画像に相当する画像として不安定であることを意味している。そこで、被検体40中のヨウ素の濃度変化に伴う水密度画像Gwの画素値の変化分を補うように、仮想非造影画像Gncを、水密度画像Gwにヨウ素密度画像Gioの成分を少しだけ加えるようにする。後は、仮想非造影画像Gncにおける画素値がCT値として適正な値を取るように、水密度画像Gwに乗ずる係数βと、ヨウ素密度画像Gioに乗ずる係数αと、オフセット分の-γとを調整する。数式1は、このような考え方によって導き出されたものである。ところで、オフセット分の-γは、仮想非造影画像Gncの画素値がCT値となり、撮影対象が空気であるときにそのCT値が-1000HUとなるようにするため、-1000が適当である。しかし、実際には、水は温度によって比重が変わり、そのX線吸収係数も変化する。また、X線CT装置のデータ収集系により得られるデータ値も、個体差による誤差が生じる。そのため、これらを考慮すると、経験的には、1100≧γ≧900の範囲が適当であり、実装レベルでは、1050≧γ≧950の範囲で調整される可能性が高い。また、係数β,αは、例えば、次のようにして求めることができる。まず、仮想非造影画像のCT値Gncを、下記のように定義する。Gnc=βα-γ...数式A1但し、OFFSETw:水密度画像Gwにおけるオフセット成分、OFFSETio:ヨウ素密度画像Gioにおけるオフセット成分次に、空気を撮影対象としてデュアルエネルギー撮影を行う。空気を撮影したときの理論上の水密度画像のCT値は0になるはずなので、このときの水密度画像の実測値をGwとしたとき、次式が成り立つ。-OFFSETw=0...数式A2∴OFFSETw=Gwこのように、空気を撮影対象としたデュアルエネルギー撮影から、オフセットOFFSETwを求めることができる。次に、水を撮影対象としてデュアルエネルギー撮影を行う。水を撮影したときの理論上のヨウ素密度画像のCT値は0になるはずなので、このときのヨウ素密度画像の実測値をGioとしたとき、次式が成り立つ。-OFFSETio=0...数式A3∴OFFSETio=Gioまた、このときの仮想非造影画像のCT値Gncは0になるはずであるから、さらに次式が成り立つ。Gnc=β-OFFSETwα-OFFSETio-γ=β-OFFSETw-γ=0...数式A4∴β=γ/-OFFSETwこのように、水を撮影対象としたデュアルエネルギー撮影から、オフセットOFFSETioと、係数βを求めることができる。なお、実際の実験結果によれば、係数βは、1.1≧β≧0.9の範囲で求められており、実装レベルでは、1.05≧β≧0.95の範囲で調整されるのが好ましい。次に、造影剤ファントムを撮影対象としてデュアルエネルギー撮影を行う。造影剤を撮影したときの仮想非造影画像GncのCT値は0になるはずなので、このときの水密度画像Gwおよびヨウ素密度画像の実測値をそれぞれGw,Gioとしたとき、次式が成り立つ。Gnc=β-OFFSETwα-OFFSETio-γ0=β-OFFSETwα-OFFSETio-γ-OFFSETwα-OFFSETio=γβα-OFFSETio=γβ--OFFSETw...数式A5∴α=γβ--OFFSETw-OFFSETioこのように、造影剤を撮影対象としたデュアルエネルギー撮影から、係数αを求めることができる。なお、造影剤の濃度を変えた複数パターンの実際の実験結果によれば、係数αは、0.30≧α≧0.24の範囲で求められており、実装レベルでは、0.28≧α≧0.26の範囲で調整されるのが好ましい。実際の係数αの決定については、造影剤の濃度別に得られた係数αの値を平均化して求めてもよい。あるいは、操作者に選択肢を与えて選択させてもよいし、操作者に推奨範囲を示して値を直接入力させてもよい。ちなみに、β,α,γを固定値とするなら、例えば、β=1,α=0.27,γ=1000が好ましい。なお、β,α,γを上記の値にそれぞれ決定し、本実施形態によって人工的に生成された仮想非造影画像と、実際に撮影された非造影画像単純画像とを比較してみたところ、被検体の肝臓やすい臓などにおいては、CT値が誤差-2~1以内という高い精度でほぼ一致することが確認されている。ステップS5では、生成したモノクロマチック画像Gmおよび仮想非造影画像Gncを表示する。これにより、造影剤ありとなしとでそれぞれ撮影することなく、1回のデュアルエネルギー撮影により、被検体40の造影画像と非造影画像とを生成することができ、被検体40の低被曝化を図ることができる。なお、上記の数式1の代わりに、例えば、次式のような多次項を含む数式を用いることもできる。Gnc=βGwαGioα2Gio-γ但し、1.1≧β≧0.9、0.30≧α≧0.24、1100≧γ≧900...数式2第二実施形態図7は、第二実施形態に係るX線CT装置における処理の流れを示す図である。ステップT1~T3は、ステップS1~S3と同様であるため、説明を省略する。ステップT4では、生成する仮想非造影画像Gncの画像空間において、注目画素を選択する。ステップT5では、モノクロマチック画像における上記注目画素の対応画素の画素値CT値Gmが、被検体40の軟部に相当する画素値範囲、例えば0~300HUに入っているか否かを判定する。ここで、対応画素の画素値Gmが、軟部に相当する画素値範囲に入っていると判定されたときには、注目画素は、造影剤が染まり易く造影剤の有無による変化が大きい軟部組織を表すものと判断し、ステップT6に進む。一方、対応画素の画素値Gmがこの画素値範囲に入っていないと判定されたときには、注目画素は、造影剤が染まり難く造影剤の有無による変化が少ない、軟部組織以外の物質を表すものと判断し、ステップT7に進む。ステップT6では、仮想非造影画像における注目画素の画素値Gncを、次の数式1にしたがって算出する。Gnc=βGwαGio-γ但し、1.05≧β≧0.95、0.29≧α≧0.25、1050≧γ≧950...数式1すなわち、注目画素の画素値Gncを、水密度画像Gwにおける対応画素の画素値Gwに係数βを乗じた値と、ヨウ素密度画像Gioにおける対応画素の画素値Gioに係数αを乗じた値と、-γとを加算して算出する。ステップT6における注目画素の画素値Gncの算出後は、ステップT8に進む。ステップT7では、注目画素の画素値Gncとして、モノクロマチック画像における対応画素の画素値Gmを採用する。軟部組織以外の物質では、造影剤が染まりにくく、造影剤の有無による変化が元々少ないので、このような物質を表す画像部分については、数式1や数式2のように、造影剤ありの軟部組織に焦点をあてて導かれた数式により生成される画像よりも、単純画像に相当する画像として一般的で実績のあるモノクロマチック画像の方がより単純画像に近くなることが期待される。そこで、軟部組織以外の物質を表す画素については、モノクロマチック画像における対応画素の画素値を採用することにする。ステップT8では、生成する仮想非造影画像の画像空間において、注目画素として次に選択すべき画素がまだあるか否かを判定する。注目画素として選択すべき画素がまだあると判定されたときには、ステップT4に戻って、新たな画素を注目画素に選択して処理を続ける。注目画素として選択すべき画素がもうないと判定されたときには、ステップT9に進む。ステップT9では、生成したモノクロマチック画像および仮想非造影画像を表示する。このような第二実施形態によれば、仮想非造影画像のうち、軟部組織以外の物質を表す部分をモノクロマチック画像とするので、造影剤の有無による変化が元々少ない部分については、単純画像に相当する画像としてより一般的で実績のある画像を採用することができ、仮想非造影画像を、実際の単純画像により近い画像とすることが期待できる。第三実施形態図8は、第三実施形態に係るX線CT装置における処理の流れを示す図である。ステップU1~U6,U8は、ステップT1~T6,T8と同様であるため、説明を省略する。ステップU7では、仮想非造影画像Gncにおける注目画素の画素値Gncを、次の数式3にしたがって算出する。Gnc=bGwfρ但し、1.0≧b≧0.9、ρ:被検体中のヨウ素濃度に係る指標、f:ρの関数...数式3すなわち、仮想非造影画像Gncにおける注目画素の画素値Gncを、水密度画像Gwにおける対応画素の画素値Gwに係数bを乗じた値に、被検体中のヨウ素濃度に係る指標ρをパラメータとする関数fの値を加算して算出する。軟部組織以外の物質では、造影剤が染まりにくく、造影剤の有無による変化が元々少ないので、このような物質を表す画像部分については、数式1や数式2のように、造影剤ありの軟部組織に焦点をあてて導かれた数式により生成される画像よりも、基本的に水密度画像Gwをベースにして、被検体中のヨウ素濃度によって定まるオフセット成分を加える方が、より単純画像に近くなることが期待される。そこで、軟部組織以外の物質を表す画素については、数式3に示すような、被検体中のヨウ素濃度に係る指標をパラメータとする関数値を含む数式によって生成される画像における対応画素の画素値を採用することにする。被検体中のヨウ素濃度に係る指標としては、例えば、被検体に投与する造影剤の濃度や、造影剤の投与開始からの撮影タイミングなどが考えられる。また、関数fとしては、例えば、f=nρn=定数を考えることができる。なお、数式3の代わりに、例えば、次式のような数式を用いることもできる。Gnc=bGwc但し、1.1≧b≧0.9、C=定数...数式4軟部組織以外の物質では、造影剤の有無による変化が元々少ないので、水密度画像Gwをベースbaseにしてオフセット成分を加えるだけでも、より単純画像に近くなる可能性がある。なお、上記の数式3,4では、係数bを、1.1≧b≧0.9としているが、これは、水密度画像をベースにするという考え方によれば、係数bは、1.1≧b≧0.9の範囲で調整されるのが好ましいという理由による。このように上記の各実施形態によれば、ヨウ素を含む造影剤が投与された被検体のX線CTデュアルエネルギー撮影により、仮想非造影画像を高精度に得ることができる。その結果、造影画像に相当するモノクロマチック画像と、非造影画像に相当する仮想非造影画像の両方を、1回のデュアルエネルギー撮影によって得ることができ、造影撮影および非造影撮影をそれぞれ行う必要がなく、造影検査における低被曝化を実現できる。以上、発明の実施形態について説明したが、発明の実施形態は、上記のものに限定されず、発明の趣旨を逸脱しない範囲内において種々の形態を取り得る。例えば、上記の各実施形態において、モノクロマチック画像は、水密度画像およびヨウ素密度画像を基に生成しているが、他の物質を基底物質とする物質密度画像を基に生成してもよい。また例えば、上記の各実施形態において、造影画像に相当する画像として、モノクロマチック画像を採用しているが、ヨウ素密度画像を採用してもよい。また例えば、上記の各実施形態において、仮想非造影画像における注目画素が被検体の軟部に相当するか否かの画素値判定に、モノクロマチック画像以外の画像、例えば第1X線管電圧画像や第2X線管電圧画像、水密度画像、ヨウ素密度画像などを用いてもよい。また例えば、上記の各実施形態はX線CT装置であるが、上記のような仮想非造影画像の生成機能を有する画像生成装置も発明の一実施形態である。また例えば、コンピュータを、上記の画像生成装置として機能させるためのプログラムもまた、発明の一実施形態である。"