"放射線を放出する放射線源と、前記放射線源に対峙して配置され、かつ、前記放射線を入射したときに当該放射線に対応したデジタル電気量の2次元データをフレーム単位で出力する検出器と、前記放射線源と前記検出器との間に提供される撮像空間において当該放射線源と当該検出器とを結ぶラインが位置的に変化するように、当該放射線源及び当該検出器、当該放射線源、又は、前記対象物の何れかを、当該放射線源、当該検出器、及び当該対象物のうちの残りの要素に対して相対的に移動させる移動手段と、前記移動手段により前記放射線源及び前記検出器、当該検出器、又は、前記対象物を移動させている間に、前記検出器から出力される前記データをフレーム単位で収集するデータ収集手段と、を備え、前記データ収集手段により収集された前記データを用いて前記対象物の撮像部位の3次元画像を生成する放射線撮像装置において、前記撮像空間に配置されるとともに、この配置によって当該撮像空間内の予め定めた断層面に位置付けられ且つ既知の位置情報を前記放射線で画像化可能なマーカを有するファントムと、前記ファントム装置を前記撮像空間に配置した状態で、前記放射線源から放出された放射線に応じて前記データ収集手段により収集されたデータに基づいて画像を作成する画像作成手段と、前記マーカの既知の位置情報と前記画像から得られた前記マーカの位置の情報とに基づき、前記放射線源と前記検出器の間の距離情報及び前記検出器に対する前記放射線源の高さ情報を演算する第1の演算手段と、前記第1の演算手段の演算結果と前記データに基づいて、前記ラインの位置の変化量を加味した、前記撮像空間における前記放射線源、前記検出器、及び前記断層面の位置関係を規定するパラメータを演算する第2の演算手段と、を備えたことを特徴とする放射線撮像装置。前記第2の演算手段により演算されるパラメータをキャリブレーションデータとして記憶する記憶手段を備えたことを特徴とする放射線撮像装置。前記放射線により前記対象物の撮像部位を撮像したときに、前記キャリブレーションデータを前記記憶手段から読み出し、当該キャリブレーションデータを用いて当該撮像部位の実在位置を3次元的に演算する第3の演算手段を備えたことを特徴とする放射線撮像装置。前記第2の演算手段は、前記ラインの位置の変化量として、当該ラインが移動することに伴う当該ライン上の所定位置の移動量を加味して前記パラメータを演算する手段であることを特徴とする請求項1~3の何れか一項に記載の放射線撮像装置。前記移動手段は、前記放射線源と前記検出器の対を前記対象物の周りに互いに異なる、円軌道及び楕円軌道を含む、曲率を有する軌道並びに角速度で回転させる手段であり、前記第2の演算手段は、前記放射線源と前記検出器が回転するときに共有する回転中心の位置を前記ライン上の所定位置として用い、この回転中心の位置の変化量を加味して前記パラメータを演算する手段であることを特徴とする請求項4に記載の放射線撮像装置。前記放射線源は前記放射線としてのX線を曝射するX線管であり、前記検出器は前記X線を検出する検出器であり、前記移動手段は、前記X線管と前記検出器の対を前記対象物の周りに互いに異なる、円軌道及び楕円軌道を含む、曲率を有する軌道並びに角速度で回転させ、且つ、当該X線管と当該検出器とを常に互いに正対させるように回転させる手段であり、前記画像作成手段は、前記データ収集手段により収集されたデータをトモシンセシス法に基づいて前記画像としてのパノラマ画像を作成する手段である、ことを特徴とする請求項1~4の何れか一項に記載の放射線撮像装置。前記軌道は円軌道であり、前記所定位置の移動量は、前記放射線源と前記検出器が互いに異なる円軌道を回転中心を共有して回転するときの当該回転中心の位置の変化量である、ことを特徴とする請求項6に記載の放射線撮像装置。前記第1の演算手段は、前記距離情報として、前記X線管と前記回転中心との間の第1の距離及び当該回転中心と前記検出器との間の第2距離を演算するとともに、前記高さ情報として、前記高さ方向における前記検出器の基準位置からの前記X線管の高さを演算する手段であることを特徴とする請求項7に記載の放射線撮像装置。前記第2の演算手段は、前記パラメータとして、前記データのフレーム番号の入力とする関数であって、前記検出器の検出面の移動距離と当該フレーム番号の変化率、前記X線管から放射されるX線束の前記対象物への投影角度、当該投影角度と前記フレーム番号の変化率、前記回転中心から前記撮像部位までの第3の距離、前記回転中心の位置の変化量、及び当該回転中心の位置座標を演算する手段であることを特徴とする請求項8に記載の放射線撮像装置。前記ファントムは、ベースと、前記断層面としての基準断層面を前記ベースへ投影して生成される基準面軌道、及び、当該基準面軌道から離間し且つ当該基準面軌道に併走する他の軌道のそれぞれに沿って、軌道毎に複数、立設される支柱と、前記複数のそれぞれ支柱に配設され、X線透過率が少なくとも当該支柱のX線透過率とは異なるマーカと、を備えたことを特徴とする請求項1~9の何れか一項に記載の放射線撮像装置。前記ファントムが有する支柱は、前記撮像部位である歯列を辿った略馬蹄形の基準断層面を前記ベースへ投影して生成される基準面軌道、及び、当該基準面軌道から離間し且つ当該基準面軌道に併走する略馬蹄形の他の軌道のそれぞれに沿って、軌道毎に複数、立設されていることを特徴とする請求項10に記載の放射線撮像装置。前記マーカは、予め設定されている前記X線の照射角度が前記位置の変化を与えない角度である基準角度における前記第1の距離、前記第2の距離、前記高さ、及び前記第3の距離の情報を得るために前記基準面軌道及び前記他の軌道のそれぞれに位置させ、高さ方向の画像の拡大率を測定するための第1のマーカを含んでいることを特徴とする請求項10又は11に記載の放射線撮像装置。前記マーカは、前記基準角度以外の照射角度における前記第3の距離の情報を得るために前記基準面軌道及び前記他の軌道の何れか一方に位置させ、高さ方向の画像の拡大率を測定するための第2のマーカを更に含んでいることを特徴とする請求項12に記載の放射線撮像装置。前記マーカは、前記投影角度の実際値の情報を得るために、予め設定されている前記X線の照射角度毎に、前記基準面軌道及び前記他の軌道のそれぞれに配置された第3のマーカを更に含んでいることを特徴とする請求項13に記載の放射線撮像装置。前記第1及び第2のマーカは、高さ方向において既知の2点の位置情報を与えるマーカであり、前記第3のマーカは前記高さ方向に沿い且つ前記基準面軌道及び前記他の軌道において当該高さ方向の位置が互いに異なり且つライン状の又はラインを想定可能なマーカであることを特徴とする請求項14に記載の放射線撮像装置。前記第1の演算手段は、前記X線の照射の基準角度において前記画像から、前記データのフレーム番号の入力とする関数であって、前記検出器の検出面の移動距離と当該フレーム番号の変化率と、前記X線管から放射されるX線束の前記対象物への投影角度と前記フレーム番号の変化率とを算出する第1の変化率算出手段と、この第1の変化率算出手段の算出結果と前記第1のマーカの拡大率とを用いて前記第1の距離、前記第2の距離、及び前記高さを算出する第1の情報算出手段とを備えたことを特徴とする請求項12~15の何れか一項に記載の放射線撮像装置。前記第2の演算手段は、前記第3のマーカの前記画像への画像化情報を用いて前記X線の照射角度の実際値を前記予め定めた照射角度毎に算出する角度算出手段と、前記X線の照射角度の実際値毎に前記検出器の検出面の移動距離と当該フレーム番号の変化率と、前記X線管から放射されるX線束の前記対象物への投影角度と前記フレーム番号の変化率とを算出する第2の変化率算出手段と、この第2の変化率算出手段の算出結果と前記第2のマーカの拡大率とを用いて前記X線の照射角度の実際値毎に、前記回転中心から前記撮像部位までの第3の距離、前記回転中心の位置の変化量、及び当該回転中心の位置座標を算出する第2の情報算出手段とを備えたことを特徴とする請求項16に記載の放射線撮像装置。前記第2の演算手段は、前記予め設定したX線の照射角度毎に前記検出器の検出面の移動距離と当該フレーム番号の変化率と、前記X線管から放射されるX線束の前記対象物への投影角度と前記フレーム番号の変化率とを算出する第2の変化率算出手段と、この第2の変化率算出手段の算出結果と前記第2のマーカの拡大率とを用いて前記X線の照射角度の実際値毎に、前記回転中心から前記撮像部位までの第3の距離、前記回転中心の位置の変化量、及び当該回転中心の位置座標を算出する第2の情報算出手段とを備えたことを特徴とする請求項16に記載の放射線撮像装置。前記記憶手段は、前記X線の照射角度の実際値毎に又は予め設定された前記X線の照射角度毎に、前記検出器の検出面の移動距離と前記フレーム番号の変化率、前記投影角度と前記フレーム番号の変化率、前記回転中心から前記撮像部位までの第3の距離、前記回転中心の位置の変化量、及び当該回転中心の位置座標をルックアップテーブルとして有することを特徴とする請求項17又は18に記載の放射線撮像装置。前記第3の演算手段は、前記基準断層面とは異なる複数の断層面を前記第1の距離、前記第2の距離、前記投影角度のフレーム番号の変化率、前記第3の距離及び前記回転中心の位置の変化量を加味して設定する断層面設定手段と、前記画像のデータ用いて当該断層面設定手段が設定した断層面それぞれの画像を再構成する再構成手段と、この再構成した複数の断層像を用いて前記撮像部位の実在位置を、前記照射角度の実際値及び前記回転中心の位置座標を用いて、常に前記X線源を睨む方向に3次元的に同定する同定手段を備えたことを特徴とする請求項17~19の何れか一項に記載の放射線撮像装置。X線を放出するX線源と前記X線を電気信号として検出する検出器とを対象物を挟んで相互に対向させ、当該X線源と当該検出器を当該対象物の周りに回転させながら、前記X線を当該検出器によりデジタル電気量のフレームデータとして検出し、前記フレームデータに基づいて前記対象物の断層面のパノラマ画像を作成するパノラマ撮像装置における、前記X線源と前記検出器との間の空間に配置されるファントムであって、ベースと、前記断層面としての基準断層面を前記ベースへ投影して生成される基準面軌道、及び、当該基準面軌道から離間し且つ当該基準面軌道に併走する他の軌道のそれぞれに沿って、軌道毎に複数、立設される支柱と、前記複数のそれぞれ支柱に配設され、X線透過率が少なくとも当該支柱のX線透過率とは異なるマーカと、を備えたことを特徴とするファントム。上述した特許文献7に記載のパノラマ撮像装置の場合、再構成される複数の断層面は、各回転角度の位置に在る検出器とX線管とを結ぶライン上のある位置に存在すると仮定し、この仮定の元にトモシンセス法を用いて各断層面のパノラマ画像を作成している。このため、別の断層面のパノラマ画像を作成するなど、断層面を変更したときに、画像の拡大率が変わるため、この変化に因る歪が画像の縦方向歯列の上下方向に生じる。画像はデジタル化されているので、基準断層面に歯列が正確に位置決めされていれば、縦横の両方向に共に歪のない画像は作ることは可能ではある。しかし、この位置決めの条件から外れた場合、必ず画像に歪が発生する。また、歯列が基準断層面に沿って位置していない場合、再構成されたパノラマ画像には横方向のボケも生じる。その場合、画像処理によって焦点の合った画像が作成されれば、横方向の歪ボケは少なくなる。しかし、その場合でも、縦方向の歪はシフトアンドアッド量には関与しないため、パノラマ画像に縦方向の歪が残る。画像に歪があると、例えば画像上の2点の距離が実体に比して正しく描出されない。むろん、その2点間の距離も正確に計測できないなど、パノラマ画像を用いた計測能力は低い。上述した縦方向の歪の原因は以下のようである。データ収集時においてX線管と検出器とが互いに異なる半径で且つ互いに正対しながら歯列の周りを回転する。この回転の間、X線管及び検出器の対の歯列に対する回転角度、すなわちX線管から検出器に照射されるX線束の照射角度を順次変化させてX線スキャンが実行される。このとき、回転角度が変化するにつれて、X線管と検出器の回転の中心回転中心:X線管及び検出器は共に、この同一の回転中心を中心に回転するの位置が歯列に近づいたり、離れたりする。この回転中心の変化によって、歯列の縦方向の拡大率が歯列の歯列方向の各位置で異なり、これがパノラマ画像の縦方向の歪発生の原因となっていた。勿論、断層面を変えたときにもこの歪は生じる。このため、特許文献7に記載のパノラマ撮像装置によって生成されたパノラマ画像は定量的な計測には不向きである。また、サブトラクションなど時系列の変化を見ることも難しい。したがって、その臨床的な用途は限られている。このような理由が、従来のパノラマ撮像装置を真の口内撮影の代替手段としては使えず、歯科用CTには及ばない、一つの理由である。ところで、パノラマ撮像装置は、各メーカ間で違いがあるのは勿論のこと、同じメーカの製造であっても、その装置毎に機械的な動作にバラツキを有している。特に、X線管と検出器を回転させる機構では、そのようなバラツキの影響は誤差として、再構成されるパノラマ画像に影響する。このため、X線管と検出器との間の位置関係、X線管及び検出器の対の回転中心の移動状況、その移動の速度、更には、X線の照射方向などのファクタが設計通りかどうかなどについて、装置毎にチェックし、そのバラツキの情報を持っておくことが必要になる。この情報を得るということは、撮像空間の構造基準断層面に対する位置関係を3次元的に把握することを意味する。この情報は装置毎に把握して、得られた情報をパノラマ再構成に反映させるべきであるが、従来、そのような必要性がないと同時に、そのような仕組みも手法も無いというのが現状であった。本発明は、上述の従来の状況に鑑みてなされたもので、撮像空間内の予め定めた断層面に対する、放射線の照射角度毎の撮像系の構造、すなわち撮像空間の3次元構造をファントムを用いたスキャンで的確に把握し、その構造情報を用いて歪がより少なく、撮影対象の3次元的な実位置をより精度良く反映した画像を提供することができる放射線撮像装置及び同装置に用いるファントムを提供することを、その目的とする。本発明は、特に、上述した構造情報を用いて、トモシンセス法の基で再構成される画像上の縦方向の歪を排除または軽減し、実体物の3次元的な位置をより精度良く反映した断層像を提供することができる放射線撮像装置及び同装置に用いるファントムを提供することを、更なる目的とする。上述した目的を達成するために、本発明は、その1つの態様として、放射線を放出する放射線源と、前記放射線源に対峙して配置され、かつ、前記放射線を入射したときに当該放射線に対応したデジタル電気量の2次元データをフレーム単位で出力する検出器と、前記放射線源と前記検出器との間に提供される撮像空間において当該放射線源と当該検出器とを結ぶラインが位置的に変化するように、当該放射線源及び当該検出器、当該放射線源、又は、前記対象物の何れかを、当該放射線源、当該検出器、及び当該対象物のうちの残りの要素に対して相対的に移動させる移動手段と、前記移動手段により前記放射線源及び前記検出器、当該検出器、又は、前記対象物を移動させている間に、前記検出器から出力される前記データをフレーム単位で収集するデータ収集手段と、を備え、前記データ収集手段により収集された前記データを用いて前記対象物の撮像部位の3次元画像を生成する放射線撮像装置において、前記撮像空間に配置されるとともに、この配置によって当該撮像空間内の予め定めた断層面に位置付けられ且つ既知の位置情報を前記放射線で画像化可能なマーカを有するファントムと、前記ファントム装置を前記撮像空間に配置した状態で、前記放射線源から放出された放射線に応じて前記データ収集手段により収集されたデータに基づいて画像を作成する画像作成手段と、前記マーカの既知の位置情報と前記画像から得られた前記マーカの位置の情報とに基づき、前記放射線源と前記検出器の間の距離情報及び前記検出器に対する前記放射線源の高さ情報を演算する第1の演算手段と、前記第1の演算手段の演算結果と前記データに基づいて、前記ラインの位置の変化量を加味した、前記撮像空間における前記放射線源、前記検出器、及び前記断層面の位置関係を規定するパラメータを演算する第2の演算手段と、を備えたことを特徴とする放射線撮像装置が提供される。また、前述した目的を達成するため、本発明は別の態様として、X線を放出するX線源と前記X線を電気信号として検出する検出器とを対象物を挟んで相互に対向させ、当該X線源と当該検出器を当該対象物の周りに回転させながら、前記X線を当該検出器によりデジタル電気量のフレームデータとして検出し、前記フレームデータに基づいて前記対象物の断層面のパノラマ画像を作成するパノラマ撮像装置における、前記X線源と前記検出器との間の空間に配置されるファントムであって、ベースと、前記断層面としての基準断層面を前記ベースへ投影して生成される基準面軌道、及び、当該基準面軌道から離間し且つ当該基準面軌道に併走する他の軌道のそれぞれに沿って、軌道毎に複数、立設される支柱と、前記複数のそれぞれ支柱に配設され、X線透過率が少なくとも当該支柱のX線透過率とは異なるマーカと、を備えたことを特徴とするファントムが提供される。以上のように、本発明によれば、放射線の照射角度毎の撮像系の構造、すなわち撮像空間の3次元構造をファントムを用いたスキャンで的確に把握し、その構造情報を用いて歪がより少なく、撮影対象の3次元的な実位置をより精度良く反映した画像を提供することができる放射線撮像装置及び同装置に用いるファントムを提供することができる。以下、添付図面を参照して、本発明の実施形態を説明する。図1~43を参照して、本発明に係る放射線撮像装置及び同装置に用いるファントムの1つの実施形態を説明する。本実施形態では、この放射線撮像装置は、歯科用のX線パノラマ撮像装置として実施されている。以下、このパノラマ撮像装置を詳述する。図1に、かかるパノラマ撮像装置1の外観を示す。このパノラマ撮像装置1は、被検体の顎部をX線でスキャンし、そのデジタル量のX線透過データから顎部の歯列の実際の3次元的位置実在位置を同定し、かつ、その歯列の縦方向の拡大率の変化違いを補償した歯列のパノラマ画像をトモシンセシス法tomosynthesisの元で生成再構成することができる。このパノラマ撮像装置1の構成の概要を説明する。図1に示すように、このパノラマ撮像装置1は、被検体患者Pからデータを例えば被検体Pの立位や車椅子に座った姿勢で収集する筐体11と、この筐体11が行うデータの収集を制御し、その収集したデータを取り込んでパノラマ画像を作成し、かつ、操作者医師、技師などとの間でインターラクティブに又は自動的にパノラマ画像の後処理を行うための、コンピュータで構成される制御演算装置12とを備える。筐体11は、スタンド部13と、このスタンド部13に対して上下動可能な撮影部14とを備える。撮影部14は、スタンド部13の支柱に所定範囲で上下動可能に取り付けられている。ここで、説明の便宜のため、パノラマ撮像装置については、スタンド部13の長手方向、すなわち上下方向をZ軸とするXYZ直交座標系を設定する。なお、後述する2次元のパノラマ画像については、その横軸方向をj軸、縦軸方向をi軸=Z軸と表記する。撮影部14は、側面からみて、略コ字状を成す上下動ユニット23と、この上下動ユニット23に回転回動可能に支持された回転ユニット24とを備える。上下動ユニット23は、スタンド部13に設置された、図示しない上下駆動機構例えば、モータ及びラックピニオンを介して、高さ方向の所定範囲に渡ってZ軸方向縦軸方向に移動可能になっている。この移動のための指令が、制御演算装置12から上記上下動駆動機構に出される。上下動ユニット23は、前述したように、その一方の側面からみて略コ字状を成し、上下それぞれの側の上側アーム23A及び下側アーム23Bと、その上側、下側アーム23A,23Bを繋ぐ縦アーム23Cとが一体に形成されている。縦アーム23Cが、前述したスタンド部13に上下動可能に支持されている。上側アーム23Aの内部には、回転駆動用の回転駆動機構30A例えば、電動モータ及び減速ギヤなどが設置されている。この回転駆動機構30Aは、制御演算装置12から回転駆動用の指令を受ける。回転駆動機構30Aの出力軸、すなわち電動モータの回転軸は、上側アーム23Aから下側Z軸方向下側に突出するように配置されており、この回転軸に、回転ユニット24が回転可能に結合されている。つまり、回転ユニット24は、上下動ユニット23に垂下されており、回転駆動機構30Aの駆動に付勢されて回転する。また、回転駆動機構30Aは移動機構30Bに連結している。この移動機構30Bは図示しない電動モータ、ギヤなどから構成されている。この移動機構30Bも、制御演算装置12から回転駆動用の指令を受けて動作し、回転駆動機構30A、すなわち回転ユニット24をXY面に沿って移動可能に構成されている。これにより、後述するX線管及び検出器の対の回転中心の移動の軌跡を、XY面に沿った所定範囲において2次元的に一定軌道に沿って移動させることができる。一方、下側アーム23Bは、その先端部にチンレスト25が形成されている。このチンレスト25には、ヘッドレスト26が着脱自在に取り付けられる。このため、被検体Pは、バイトブロック単にバイトとも呼ばれるを咥え、顎をチンレスト25に載せ、且つ、額をヘッドレストに押し当てる。これにより、被検体Pの口腔部が後述する撮像空間内の所定位置に固定される。回転ユニット24は、その使用状態において、その一方の側面からみて略コ字状に形成され、その開放端側を下側に向けて回転自在に上側アーム23Aのモータ出力軸に取り付けられている。詳しくは、横方向、すなわちXY平面内で略平行に回転回動する横アーム24Aと、この横アーム24Aの両端部から下方Z軸方向に伸びた左右の縦アーム第1の縦アーム、第2の縦アーム24B,24Cとを一体に備える。回転ユニット24も制御演算装置12の制御下で駆動及び動作するようになっている。第1の縦アーム24Bの内部の下端部に放射線放出源としてのX線管31が装備されている。このX線管31は、例えば回転陽極X線管で構成されており、そのターゲット陽極からX線を第2の縦アーム24Cに向けて放射状に放射させる。このターゲットに衝突させる電子線の焦点は、例えば径0.5mm~1mm程度と小さく、したがって、このX線管31は点状のX線源として機能する。X線管31の前面の所定位置には、スリット状のコリメータ33が取り付けられている。このコリメータ33により、検出器32に入射させるX線を、その検出面すなわち実際の収集用の窓例えば5.0mm幅の窓に絞ることができる。一方、第2の縦アーム24Cの内部の下端部に放射線検出手段としての、X線検出素子を2次元状例えば、64×1500のマトリクス状に配置したデジタル形X線検出器32が装備されており、その検出面に入射するX線を検出する。この検出器32は、一例として、CdTeで作られた、縦長形の検出面例えば、横6.4mm×縦150mmを有している。なお、本実施形態はトモシンセシス法を採用しているため、検出器32はその横幅方向にも複数のX線検出素子を持つことが必須である。この検出器32は、その縦方向をZ軸方向に一致させて縦方向に配置される。この検出器32の検出時の横方向の有効幅は、前述したコリメータ33によって例えば約5.0mmに設定される。この検出器32は、例えば300fpsのフレームレート1フレームは、例えば、64×1500画素で入射X線を、当該X線の量に応じたデジタル電気量の画像データとして収集することができる。以下、この収集データをフレームデータと呼ぶ。撮影時には、X線管31及び検出器32は、被検体Pの口腔部を挟んで互いに斜めに対向して又は互いに正対して位置し、口腔部の周りを回転する。X線管31及び検出器32が互いに正対してとは、X線管31から照射されたX線ビームの中心軸XY面に投影したときの中心軸が検出器32の検出面に直交する状態を言う。また、X線管31及び検出器32が互いに斜めに対向してとは、上記X線ビームの中心軸が検出器32の検出面に90°以外の角度0°<角度<90°で入射する状態を言う。したがって、トモシンセス法に拠るスキャンを実行するときに必要なX線管31及び検出器32の回転移動のさせ方には、このように様々な態様を採ることができる。ただし、本実施形態においては、X線管31及び検出器32は、被検体Pの口腔部を挟んで常に互いに正対するように位置し、それらの対毎、一体に口腔部の周りを回転するように駆動される。ただし、この回転は単純な円を画く回転ではない。つまり、本実施形態の場合、X線管31及び検出器32の対は、その移動する回転中心RC図2,3参照を中心に回転駆動される。後述するが、回転中心RCからX線管31及び検出器32までの距離は縦方向の拡大率を考慮して決められている。回転中心RCは、図2に示す如く、歯列の奥側から前歯に向けてほぼ直線状に往復する軌道A、又は、同様に歯列の奥側から前歯に向けて直線的に進み、その後、鋭角に折り返して直線的に戻るという、少し角度を付けた3角形状の軌道Bを描くように回転駆動される。このため、X線管31及び検出器32は、角速度を適宜に変えながら回転している。なお、回転中心RCの軌道としては、一方の臼歯部の側から前歯に向けて円弧状に進み、その後、同様に円弧状に反対側の臼歯部の側戻るという山形状の軌道Cもある。この軌道Cは、特に、標準的な形状及びサイズの歯列に沿った断層面以下、3D基準断層面にX線焦点を合わせ且つその3D基準断層面を機械的に追従するように予め設計された軌道である。この3D基準断層面SSにX線焦点を追従させる場合も、X線管31及び検出器32は角速度を適宜に変えながら回転している。このように、X線管31及び検出器32が回転しながらX線によるスキャンが実行されるので、X線管31及び検出器32の回転軌跡によって囲まれる内部の空間は撮像空間を形成している。3D基準断層面をZ軸方向から見たときのXY面上の軌跡は、上述したように、略馬蹄形を成すもので、図2に一例を示す。この3D基準断層面の軌跡は、例えば文献RMolteniAuniversaltestphantomfordentalpanoramicradiographyMedicaMudivol36no31991によっても知られている。X線管31、3D基準断層面SS、検出器32、及び、回転中心RCの位置を貫く回転軸AXzの幾何学的な位置関係は図3に示すようになる。ただし、回転中心RCの位置は、図2に点線A、B、及びCで示したように、X線照射角度に応じて変わる。3D基準断層面SSは検出器32の入射口X線検出面Ldet:図6参照に平行であり、Z軸方向に沿った湾曲した断面であって2次元に展開したときには細長い矩形状の断面として設定されている。図4に、このパノラマ撮像装置の制御及び処理のための電気的なブロック図を示す。同図に示す如く、X線管31は高電圧発生器41及び通信ライン42を介して制御演算装置12に接続され、検出器32は通信ライン43を介して制御演算装置12に接続されている。高電圧発生器41は、スタンド部13、上下動ユニット23、又は回転ユニット24に備えられ、制御演算装置12からの制御信号により、X線管31に対する管電流及び管電圧などのX線曝射条件、並びに、曝射タイミングのシーケンスに応じて制御される。制御演算装置12は、例えば大量の画像データを扱うため、大容量の画像データを格納可能な、例えばパーソナルコンピュータで構成される。つまり、制御演算装置12は、その主要な構成要素して、内部バス50を介して相互に通信可能に接続されたインターフェース51,52,62、バッファメモリ53、画像メモリ54、フレームメモリ55、画像プロセッサ56、コントローラCPU57、及びD/A変換器59を備える。コントローラ57には操作器58が通信可能に接続され、また、D/A変換器59はモニタ60にも接続されている。このうち、インターフェース51,52はそれぞれ高電圧発生器41、検出器32に接続されており、コントローラ57と高電圧発生器41、検出器32との間で交わされる制御情報や収集データの通信を媒介する。また、別のインターフェース62は、内部バス50と通信ラインとを結ぶもので、コントローラ57が外部の装置と通信可能になっている。これにより、コントローラ57は、外部に在る口内X線撮影装置により撮影された口内画像をも取り込めるとともに、本撮影装置で撮影したパノラマ画像を例えばDICOMDigitalImagingandCommunicationsinMedicine規格により外部のサーバに送出できるようになっている。バッファメモリ53は、インターフェース52を介して受信した、検出器32からのデジタル量のフレームデータを一時的に記憶する。また、画像プロセッサ56は、コントローラ57の制御下に置かれ、撮像空間の構造解析、キャリブレーションのための処理、3D基準断層面のパノラマ画像の作成、及びそのパノラマ画像の後利用のための処理を操作者との間でインターラクティブに実行する機能を有する。この機能を実現するためのプログラムは、ROM61に予め格納されている。このため、このROM61は、本発明に係るプログラムを格納する記録媒体として機能する。なお、このプログラムは予めROM61に格納しておいてもよいが、場合によっては、外部システムから通信回線や持ち運び可能なメモリを介して、図示しないRAMなどの記録媒体にインストールするようにしてもよい。上述した3D基準断層面は、本実施形態では、装置側で予め用意されているものである。なお、3D基準断層面は、装置側で予め用意された複数の断層面から撮影前に選択するようにしてもよい。つまり、3D基準断層面としての固定した断面であることには変わりは無いが、かかる選択動作によって、3D基準断層面の位置を、歯列の奥行き前後方向の一定範囲で変更可能にしてもよい。画像プロセッサ56により処理される又は処理途中のフレームデータ、画像データ、及び、キャリブレーション用のルックアップテーブルLUTは、画像メモリ54に読出し書込み可能に格納される。画像メモリ54には、例えばハードディスクなどの大容量の記録媒体不揮発性且つ読出し書込み可能が使用される。また、フレームメモリ55は、再構成されたパノラマ画像データ、後処理されるパノラマ画像データなどを表示するために使用される。フレームメモリ55に記憶される画像データは、所定周期でD/A変換器59に呼び出されてアナログ信号に変換され、モニタ60の画面に表示される。コントローラ57は、ROM61に予め格納されている制御及び処理の全体を担うプログラムに沿って、装置の構成要素の全体の動作を制御する。かかるプログラムは、操作者からそれぞれに制御項目についてインターラクティブに操作情報を受け付けるように設定されている。このため、コントローラ57は、後述するように、フレームデータの収集スキャンなどを実行可能に構成されている。このため、患者は、図1に示すように、立位又は座位などの姿勢でチンレスト25の位置に顎を置いてバイトブロックを咥えるともに、ヘッドレスト26に額を押し当てる。これにより、患者の頭部顎部の位置が回転ユニット24の回転空間のほぼ中央部で固定される。この状態で、コントローラ57の制御の元、回転ユニット24が患者の顎部の周りをXY面に沿って回転する。なお、バイトブロックを噛むことによる位置決めの代わりに、患者は綿のようなものを噛み、イヤーロッドで耳を固定し、後は正面の鏡で回転が曲がっていないか確認することで、位置決めするようにしてもよい。この回転の間に、コントローラ57からの制御の元で、高電圧発生器41が所定周期のパルスモードで曝射用の高電圧指定された管電圧及び管電流をX線管31に供給させ、X線管31をパルスモードで駆動させる。これにより、X線管31から所定周期でパルス状のX線が曝射される。このパルス駆動には、半波整流した駆動信号を使う場合もあるし、インバータ回路を用いたDC駆動方式の駆動信号を使う場合もある。このX線は、撮影位置に位置する患者の顎部歯列部分を透過して検出器32に入射する。検出器32は、前述したように、非常に高速のフレームレート例えば300fpsで入射X線を検出し、対応する電気量の2次元のデジタルデータ例えば64×1500画素をフレーム単位で順次出力する。このフレームデータは、通信ライン43を介して、制御演算装置12のインターフェース52を介してバッファメモリ53に一時的に保管される。この一時保管されたフレームデータは、その後、画像メモリ53に転送されて保管される。このため、画像プロセッサ56は、画像メモリ53に保管されたフレームデータを用いて3D基準断層面SSに撮像焦点を当てた断層像を基準面パノラマ画像として再構成する。つまり、この基準面パノラマ画像は、3D基準断層面SSに沿って歯列が存在していると仮定したときのパノラマ画像であると定義される。また、この画像プロセッサ56は、この基準面パノラマ画像を用いて3次元3D画像及び3Dオートフォーカス画像を作成するなどの処理を行う。この処理の概要を図5に示す。なお、上記3D基準画像は、3D基準断層面SSに沿って歯列が存在していると仮定したとき3次元画像として定義される。また、3Dオートフォーカス画像は、3D基準画像からフレームデータ又は基準面パノラマ画像のデータを用いて歯列を自動的に最適焦点化した表面画像として定義される。つまり、この3Dオートフォーカス画像は、ボケが少なく、かつ、歯列の実在位置及びその実際のサイズを精度良く表現した最適焦点画像である。とくに、3Dオートフォーカス画像は、被検体個々によって歯列の実在位置が異なるという事実を考慮した画像である。つまり、被検体個々の歯列は3D基準断層面SS図6参照に沿っていることは無く、3D基準断層面SSから部分的に又は全体的にずれていたり、その面から傾いていたりする。このため、3Dオートフォーカス画像は、被検体個々の歯列の実際の3次元位置形状を自動的に且つ精度良く同定するとともに、その同定結果から実際の歯列形状を自動的に描出することで作成される。X線管31X線管焦点は点状と見做されるから照射されたX線は、コリメータ33を通して照射される。このX線X線束は、被検体Pの口腔部を透過して、Z軸方向に一定の長さを有し且つ横幅のある縦長の検出器32により検出される。このため、X線の照射方向は図3,6に示すようにオブリークになる。したがって、歯の実際の大きさとその歯の陰影が検出器32の検出面Ldetに作る投影像の大きさとの比本実施例では、この比を拡大率というは、回転中心RCの位置に応じて変化する。この拡大率は縦方向及び横方向つまり、通常、歯の縦方向及び横方向になるのそれぞれに存在する。図6の縦方向の例で言えば、歯の実際の高さP1realと検出面Ldet上の高さP1detとの比が回転中心RCの位置に応じて変わる。この回転中心RCの位置は、図2に例示した如く、1回のスキャンデータ収集の間に変化するように、X線管31及び検出器32の対の回転軌道が予め設定されている。この理由は以下のようである。図6に示すように、X線管31と検出器32との間の距離RsRdは一定に保持され、かつ、回転中心RCからX線管31及び検出器32に至る距離Rs,Rd本実施形態ではRs>Rdも一定に保持される。その一方で、3D基準断層面SSに焦点を合わせたスキャンを行うため、1回のスキャンデータ収集の間に、回転中心RCの位置の軌道は、略馬蹄形状に湾曲している歯列に対して、一例として、前述のように鋭角な3角形状の軌道B図2参照を採用している。具体的には、回転中心RCから3D基準断層面SSまでの距離Dと、検出器32から3D基準断層面SSまでの距離Rd-Dとがスキャンが進むにつれて変化する。これに応じて、回転中心RCは歯列に近づいたり遠ざかったりするので、X線管31も歯列に近づいたり遠ざかったりする。X線管31のX線焦点は点状と見做されるので、高さについて言えば、同一高さの歯であっても、X線管31が歯列に近いほど検出面Ldetへの投影像は大きくなる。すなわち、拡大率は大きい。図2の例で言えば前歯部をスキャンするときの方が臼歯部奥歯側をスキャンするときに比べて、回転中心RCが歯列に近くなり、その分、拡大率は大きくなる。例えば、前歯部をスキャンするX線照射方向θ=0°のときの距離d1は、臼歯部をスキャンするX線照射方向θ=60°、75°のときの距離d2、d3に対して、d1<d2、d1<d3、d2<d3の関係にある。図2に示す回転中心RCの軌跡はあくまで一例であるが、この回転中心RCが歯列に近づいて遠ざかることは、3D基準断層面SSに焦点を合わせるようにX線をスキャンさせるパノラマ撮像装置の場合、通常、当てはまる設計事項である。このように拡大率は歯列の位置によって変わるので、この拡大率の影響を補正したパノラマ画像をしない限り、口腔部の構造や時系列的な変化を精度良くかつ定量的に解析することはできない。さらに、被検体の実際の歯列は、その全体にせよ部分的にせよ、3D基準断層面SSの位置には無いことが殆どである。このため、拡大率の影響を回避する場合、歯列の3D基準断層面からのずれも合わせて考慮しなければならない。しかしながら、従来のパノラマ画像は、上述した拡大率に因る問題及び実際の歯列のずれを考慮しないで作成されている。このため、従来のパノラマ画像を用いた歯列の定量的な構造解析は非常に困難である。したがって、被検体毎の様々な形状の歯列や様々な位置にある歯列を対象とする場合であっても、また、同一歯列の中の歯の位置の如何に関わらず、3次元の実在位置を精度良く同定することができるパノラマ画像が望まれていた。この点、本実施例に係るパノラマ撮像装置は、同一の歯列であっても拡大率が歯の位置毎に異なることに因る画像の歪みを解消しつつ、歯列の実際の3次元位置形状を含むを自動的に且つ精度良く同定することを特徴の一つとしている。これにより、極めて位置形状の同定精度の高い3次元パノラマ画像を提供することができる。本実施例では、断層面の画像を得るためにトモシンセシス法tomosynthesisを用いている。つまり、スキャンによって一定レートで収集される多数のフレームデータは、2次元のパノラマ画像の写像位置に応じて互いにシフトさせつつ相互に加算処理シフトアッドされる。このため、本実施例で言う最適焦点とは、焦点が一番合っている、焦点ボケが少ないという意味であり、注目する部位がそれ以外の部位よりも解像度が良い、又は、画像の全体の解像度がより高いことを言う。基準面パノラマ画像が作成されると、そのデータは画像メモリ54に保管されるともに、モニタ60に適宜な態様で表示される。このうち、表示態様などについて、操作器58から与える操作者の意思が反映される。撮像空間を規定するパラメータのキャリブレーション撮像を説明する前に、図7~図20を用いて、ファントムを使った、撮像空間における基準断層面に対する撮像系の3次元構造を示す幾何学的なパラメータの値や変化量を推定する処理、すなわちキャリブレーションを説明する。このキャリブレーションの結果は、画像再構成に反映させられるとともに、必要に応じて、撮像空間の構造解析や設計に用いられる。このキャリブレーションに伴う処理は、コントローラ57及び画像プロセッサ56が協働して実行される。キャリブレーション専用のプロセッサを設けてもよい。このキャリブレーションには、本実施例では、被検体Pの歯列を模したファントムを用いることを特徴する。ファントム図7に、このファントム101の一部破断した外観を示す。このファントム101は、かかるキャリブレーションに必要なパラメータの測定を1個で間に合わせることができる万能型ファントムである。なお、本発明に係るファントムは、必ずしも、この万能型ファントムに限定されるものではなく、後述するように3D画像再構成に必要なパラメータをキャリブレーションを実施できるものであれば、その形態は様々に変形可能であることは勿論である。この変形例の幾つかも後に説明される。この万能型ファントム101は、透明な樹脂製の板状のベース111及び天板112と、このベース111及び天板112に挟持された複数の支柱113とを備える。これらの支柱113113 ́には、後述するようにX線透過率が樹脂材とは異なる金属製のマーカを備えている。なお、樹脂の種類は例えばアクリル樹脂であるが、X線透過率がマーカのそれと異なるものであればよい。また、樹脂が透明としたのは、マーカが見易いからである。支柱113113 ́のそれぞれは、その上下端がベース111及び天板112に差し込まれて固定されている。以下、これを詳述する。ベース111は、図7,8に示すように、四角い板状を成し、透明な樹脂部材で製造される。このベース111の上面には、3次元の基準断層面SSをXY面に投影したときの基準面軌道ORsと、この基準面軌道ORsよりも所定距離DS、例えば20mm外側に当該基準面軌道ORsに例えば平行に引かれた外側面軌道ORouterとが設定されている。これらの軌道ORsORouterは、オペレータに分かり易いように、ベース111の面上に実際に線として描出してもよいし、仮想的なものであってもよい。このベース111の上面には、それらの両方の軌道ORsORouterと基準断層面SSに焦点が合うようにX線管31及び検出器32を回転移動させるときのX線照射角度θそれぞれとの交点に四角い植設穴111Aが形成される。なお、上記両軌道ORs、ORouter間の距離DSは必ずしも20mmに設定する必要はなく、限られたサイズ関係の中で後述するパラメータをより精度良く演算できる値であればよい。複数の支柱113は、それぞれ、図9,10に示すように、アクリルなどの樹脂製の角柱として形成されている。各支柱113は、一定長さの角柱状の支柱本体113Aと、その上下端それぞれに一体に突設された四角柱状の突起113Bとを備える。支柱本体113Aは、その長手方向に直交する断面サイズが例えば5mm×5mmであり、その長さが92mmである。各突起113Bのサイズは、支柱本体113Aのそれよりも小さい断面であって、例えば高さ5mm程度の長さになっている。各支柱本体113Aの一面には、キャリブレーション用の第1、第2、及び第3のマーカ114,115、及び116が配設されている。これらのマーカ114,115、及び116は全て、アルミ製又は黄銅製の小径のロッドであり、その径は例えば0.6mmである。このうち、第1及び第2のマーカ114、115は、支柱本体113Aの上端、下端からそれぞれ所定距離、例えば10mm、15mmだけ離れた位置に横向きに配設されている。支柱本体11Aの表面に例えば径0.6mmの断面半円状の削りを入れて、その削り部分に小径ロッドとしての第1及び第2のマーカ114、115を固設する。さらに、図9に示すように、第3のマーカ116は、支柱本体113Aの上端から例えば30mm離れた位置を中心に縦方向に沿って固設されている。この第3のマーカ116は一定の長さを有し、その長さは例えば20mmである。この第3のマーカ116は、上述した第1及び第2のマーカ114,115と同様の方法で植設されている。なお、上述した支柱113及びマーカ位置の寸法はあくまで例示であって、その他の適宜な寸法で設計することができる。以上、基準面軌道ORsに沿って配設されるファントム113を図7について説明した。一方、外側面軌道ORouterに沿って配設されるファントム113 ́は、図10に示したように構成されている。ここで興味深い特徴は、図32に示したファントム113を上下逆さまにしたものが図10のファントム113 ́になっているということである。このため、各ファントム113 ́にも第2、第1のマーカ115,114が上下端寄りに横向きに位置し、且つ第1のマーカ114寄りの位置で縦向きに位置している。マーカの植設方法も全く同じであるので、ファントム101を組み立てるときに、基準面軌道ORsと外側面軌道ORouterとの間で逆さまに向きを変えればよく、パーツの共通化を図り、製造コストを下げることができる。勿論、互いの逆さまの向きを混同しないようにX線透過に影響を与えない上下端の目印を付けたり、植設用の突起113B及び植設穴111Aの形状を、上下のベース111及び天板112の間で違えるなどの変形を施してもよい。上述のように、第1及び第2のマーカ114,115と第3のマーカ116の植設方向及び長さは互いに異なる。その理由は、キャリブレーションにおいて異なるパラメータの測定が必要であり、そのパラメータの属性に合わせた形状の異なる種類のマーカが必要であることによる。このように、本実施例では、1つのファントム101に全ての必要な種類のマーカを効率的に且つ無駄無く配したことも特徴の1つである。これ故に、パラメータの種類に合わせた複数のファントムを使用しなくても良いという効果がある。第1及び第2のマーカ114,115は、後述するが、撮像空間に存在するX線管31、検出器32、回転中心RC、及び3D基準断層面SSの間の距離関係の情報及びX線管31の検出器32に対する高さ位置の情報を得るためのマーカである。これに対し、第3のマーカ116は、後述するゲインと呼ばれる量=ΔX/ΔFi、および、X線照射角度θそれぞれに対する実際の投影角度を測定するためのマーカである。すなわち、基準面軌道ORs及び外側面軌道ORouterに在る第1、第2、及び第3のマーカ114,115、及び116が基準面パノラマ画像及び外側面パノラマ画像に写り込む。これを例えばX線照射角度θ=75°についてみた場合、例えば基準面パノラマ画像には図11Aに示すように描出される。つまり、図11Bに示す幾何学的関係から、基準面パノラマ画像には、上から、基準面軌道ORsに在る第1のマーカ114ORs、外側面軌道ORouterに在る第2のマーカ115ORouter、基準面軌道ORsに在る第3のマーカ116ORs、外側面軌道ORouterに在る第3のマーカ116ORouter、基準面軌道ORsに在る第2のマーカ115ORs、及び、外側面軌道ORouterに在る第1のマーカ114ORouterの順で並んで黒く描出される。逆に言えば、このような並びで描出されるように、外側面軌道ORouterの基準面軌道ORsに対する離間距離及び各マーカの縦方向の位置を設定してある。ただし、外側面軌道ORouterに在るマーカ114ORouter,115ORouter,116ORouterの画像は、基準面軌道ORsに在るそれらの画像よりはボケている。なお、シフトアンドアッド量を変更して外側断層面に焦点を合わせてパノラマ画像を再構成すれば、そのボケの程度、つまり、最適焦点の画像か否かの関係は逆になる。パノラマ画像上では、4つのマーカ114ORs、115ORouter、115ORs、及び114ORouterの画像は横向の黒線として描出され、X線管31、検出器32、回転中心RC、及び基準断層面SSの間の距離関係のパラメータ、並びに、X線管31の検出器32に対する高さ位置のパラメータを測定するために使用される。また、2つのマーカ116ORs及び116ORouterの画像は、縦向きの黒線として描出され、後述するゲインと呼ばれる量=ΔX/ΔFi、および、X線照射角度θそれぞれに対する実際の投影角度を測定するために使用される。X線の照射角度が設計値又は想定値からずれている場合、実際の投影角度もそれらの値からずれるので、2つのマーカ116ORs及び116ORouterの縦向きの黒線の位置は一致せず、横向方向にずれて描出される。このずれを演算することで実際の投影角度のずれを測定できる。このように、ファントム101は、撮像空間における撮像系の距離及び高さに関して一度のスキャンで必要充分な位置情報を与えることができる。このため、このファントム101、異なる種類のパラメータを単独で測定可能な汎用性を発揮する。なお、天板112は必ずしも設けなくてもよい。しかし、ベース111に植設した複数の支柱113には、そのマーカ114,115,116が高精度な空間位置を保つことが要求される。このため、設置や保管のときに、支柱113が傾いたり、ずれたり、損傷したりすることを防止するためには、天板112は設ける方がよい。天板112とベース111との間に、両板を支えるためだけの樹脂性の支柱を設けてもよい。<再構成の原理>ここで、パノラマ撮像装置における再構成の基本原理を数式的に説明する。図12Aは、略馬蹄形の歯列の周りを、互いに正対し且つそれぞれ異なる曲線軌道TSTDに沿って回転移動するX線管31及び検出器32の様子を示す。X線管31が一方の軌道TSに沿って回転し、検出器32が他方の軌道TDに沿って回転する。つまり、X線管31及び検出器32は対となって回転するが、その対の回転の中心回転中心RCも移動する。同図Aの例では、回転中心RCは図2の山形状の軌道Cに沿って移動する場合で示しているが、図2の軌道Aや軌道Bであっても同様の原理が成り立つ。いま、図12Aに示すように、X線管31及び検出器32の対の回転中心RCが位置Oに在り、X線管31の焦点位置をS1、検出器32の幅方向の中心位置をC1、X線管31の回転半径以下、X線管回転中心距離をRS、検出器32の回転半径以下、検出器回転中心距離をRD、回転中心RCの位置Oから歯列のある点Qまでの距離以下、回転中心基準断層面距離をD、及び、回転中心RCが描く軌道をTOとする。この距離RsとRdは固定値である。さらに、図12Bは、X線管31の焦点位置SがS1からS2に回転移動するときの状態を示すもので、この移動により、回転中心RCの軌道TOは半径αの円上を角速度ωで動いて、回転中心RCがOS1からOS2に移動するものとする。この場合、その回転中心RCの2つの位置OS1、OS2が両方の焦点位置SがS1、S2に作る角度θ1はθ1=wtt:時間となる。一方、かかる回転移動によって、歯列の点Qの検出器32への投影点はそれまでのPAS1からPAS2に変化する。このとき、検出器32の幅方向の中心位置はそれぞれC1、C2である。この図12Bの幾何学的関係を、回転中心RCの位置がOS1からOS2への移動とその軌道TOとの関係に着目して抜き出すと、図13Aのように表される。この両位置OS1、OS2間の距離は微小であるから、角度θ1と半径αを使ってθ1αと表すことができる。この結果、回転中心位置OS1、OS2、X線焦点位置S2、及び再構成位置Qは全て距離の関係で表すことができ、図13Aに示すようになる。つまり、X線管位置S2と回転中心位置OS2との間の距離がX線管回転中心間距離Rs、及び、回転中心位置OS2と検出器中心位置C2との間の距離が検出器回転中心間距離Rdであるから、回転中心位置OS1とOS2の間の距離がαθ1、回転中心位置OS2と再構成位置Qとの間の距離がD-αθ1、再構成位置Qから線分OS2-C2に垂直に下ろした線分がD-αθ1sinθ1、さらに、その垂線の交点Bと回転中心位置OS2との間に距離がD-αθ1cosθ1となる。本実施態様では、撮像空間における撮像系の幾何学的な位置関係の解析構造解析や、撮像空間における歯列の実体位置の抽出を行う3D画像再構成オートフォーカスと呼ぶに必要なパラメータをキャリブレーションするための演算に回転中心位置OS1、OS2間の距離αを考慮する点がポイントである。ゲインの演算上述した図13Bに示す距離関係を用いてゲインと呼ばれる量=ΔX/ΔFiを求める。図13Bに示す幾何学的関係から、xRssinθ1...1の関係が成り立つ。αθ1を補正項M=αθ1と捉えると、θ1及びxは微小なので、近似式としてΔxΔθ...2が成り立つ。検出器32が出力するフレームデータをFiとおくと、ΔxΔθ...3なので、ΔxΔFi...4となる。この4式の左辺Δx/ΔFiはゲインシフトアンドアッド量の変化率と呼ばれる。このゲインΔX/ΔFiは、複数のフレームデータを互いにシフトして加算するトモシンセス法つまり、シフトアンドアッド演算におけるシフトアンドアッド量の変化率を示している。また、4式の右辺中のRSRDは、検出器とX線管との間の距離検出器X線管間距離を表し、RSD-Mは回転中心位置OS1、OS2間の移動距離αθ1の分だけ補正された、X線管と焦点との間の距離焦点位置X線管間距離を表している。また、D-Mは、かかる移動距離αθの分だけ差し引いた、新しい回転中心の位置と再構成点Qとの間の距離を表している。この結果、ゲインΔX/ΔFiの曲線以下、単にゲイン曲線というは、検出器X線管間距離RSRD、焦点位置X線管間距離RSD-M、回転中心再構成点間距離D-M、及び、フレームデータFiと回転角度θとの関係を表す角速度曲線Δθ/ΔFi図17参照に基づいて演算できる。このゲイン曲線を積分し、前歯の中心位置を画像の中心位置になるようにすれば、回転角度のそれぞれにて、回転中心RCから距離Dの位置で焦点の合ったパノラマ画像を再構成することができる。なお、特開2007-136163に示されているように、上述したゲインΔX/ΔFiの大小は通常の電気回路などのそれとは概念が異なり、ゲインΔX/ΔFiが大きいほど、フレームデータ同士を相互に加算するときのフレームデータの重ね合わせ量シフト量は小さくなる。反対に、ゲインΔX/ΔFiが小さくなるほど、その重ね合わせ量は大きくなる。本実施例では、撮像空間の構造解析やキャリブレーションに必要なパラメータを上述した4式のゲイン式に基づき、且つ、キャリブレーションファントムを使って求める。このため、先にキャリブレーションファントムの構成及び機能を説明する。パラメータの演算次に、図14に基づいて、撮像空間の構造解析やキャリブレーションに必要なパラメータを測定するための演算について説明する。つまり、ここで挙げるパラメータとしては、構造解析において、X線管回転中心距離RS、X線管回転中心距離RS、X線管31の検出器32に対するZ軸方向の高さB1、及び、キャリブレーションにおいて、ゲインΔx/ΔFi、X線照射角度θ、角速度曲線Δθ/ΔFi、回転中心基準断層面距離D、補正項M、移動する回転中心RCのXY面上の座標CX,CYである。この内、キャリブレーション用のパラメータΔx/ΔFi、θ、Δθ/ΔFi、D、M、CX,CYは、入力値FiのルックアップテーブルLUTとして記憶更新される。これらのパラメータを演算するための処理として、処理1ファントムの設置とキャリブレーション用のX線撮影スキャン、処理2ゲインΔx/ΔFiのプロファイルの演算、処理3X線照射角度θのずれ実際の投影角度θ ́の演算処理4角速度曲線θ=fFi:ΔθΔFiの演算、処理5パラメータRs,Rd,B1の演算、処理6パラメータΔx/ΔFi、θ、Δθ/ΔFi、D、M、CX,CYの演算、更新:つまりキャリブレーション及び、処理7歯列の実体位置を抽出した3D再構成が挙げられる。これを処理は、コントローラ57と画像プロセッサ56が協働して実行する、図14に示すフローチャートの中で実行される。このフローチャートに沿って説明する。<処理1>コントローラ57はオペレータにファントム101を、パノラマ撮像装置の撮像空間の所定位置に設置するように画面や音声で指示するステップS1。この所定位置とは、撮像時に患者Pが顎を載せるチンレスト25の位置である。次いで、コントローラ47はキャリブレーション用スキャンの実行を操作器58で指示するステップS2。この指示に応答して、コントローラ57はROM61に予め格納しているキャリブレーション用スキャンのプログラムをそのワークメモリに読み出す。コントローラ57は、このプログラムを実行することで、コリメータ33が付いたX線管31と検出器32とをファントムの周りを回転させる。この回転の間に、X線管31の点状のX線焦点から例えばパルスX線が曝射され。このパルスX線はコリメータ33によりコリメートされてファン状のX線ビームとなる。このX線ビームがファントムを透過し、検出器32の検出面に入射する。これにより、検出器32はファントムを透過してX線を検出して、それに対応するデジタル電気量のフレームデータを一定時間毎例えば300fpsに出力する。X線管31と検出器32は、単純にファントムの周りを回るのではなく、図12Aに示すように、X線管31と検出器32とが常に正対しつつも、その両者を結ぶ線分上の回転中心RCの位置がファントムの前側に接近した後、離れる軌道を追従するように回転する。つまり、実際のスキャン時には、歯列の前歯付近に進むほど、回転中心RCが歯列に接近し、回転中心RCの位置がずれていく。この移動を許容するように、X線管31及び検出器32の回転位置、角速度がそれぞれ個別に制御される。検出器32から出力されたフレームデータはバッファメモリ53に一時保管される。画像プロセッサ53は、そのフレームデータを用いてトモシンセス法の元に基準断層面SSの基準面パノラマ画像を再構成するステップS3。<処理2>次いで、画像プロセッサ56はゲインΔX/ΔFiを求めるステップS4。まず、再構成された基準面パノラマ画像上で、基準断層面SSに沿った軌道を辿る基準面位置にX線照射角度θ毎に配置された、ファントム101の支柱のマーカが中心に描出されているフレームデータの番号Fi0を決める。この決定は、オペレータが基準面パノラマ画像を目視しながら決めればよい。なお、この基準面パノラマ画像には、基準断層面SSよりも20mm外側の断層面の軌道を辿る外側面位置にX線照射角度θ毎に配置した支柱のファントムも当然に写り込んでいる。次いで、基準面位置の支柱それぞれのファントムに対して、一番焦点が合うフレームデータFiの重ね合わせ量シフトアンドアッド量Xを求める。これもオペレータが基準面パノラマ画像を観察しながら操作器58を操作し、中心フレームデータFi0の両サイドにあるフレームデータFiを重ねてその画像のボケを観察するという方法を試行錯誤で繰り返して決める。この結果、基準断層面SSに沿った各支柱のマーカに対応した中心フレームデータFi0とその最適な重ね合わせ量Xが決まったので、それらのデータを滑らかに繋ぎ合せて重ね合わせ量のプロファイルPxを求める。このプロファイルPxからX線照射角度θの設定値毎のゲインΔX/ΔFiを求める。なお、重ね合わせ量を横軸にとり、マーカ像のエッジ統計量例えば半値幅を縦軸にとったグラフを生成し、このグラフのエッジ統計量がピークとなる点を推定するようにしてもよい。この推定値から最適な重ね合わせ量を演算すればよい。これによれば、基準面パノラマ画像上で、写りこんでいるマーカ像の位置を点ROIなどで指定すれば、その指定位置における最適な重ね合わせ量をほぼ自動的に演算することができる。次いで、画像プロセッサ56は、コントローラ57からの指示に基づいてキャリブレーションの程度の指示を受ける。本実施形態では、X線照射角度θについてキャリブレーションを実施せずに、システムが予め有しているX線照射角度θの設計値をそのまま採用する簡易型のキャリブレーションと、ファントム101をスキャンして得たパノラマ画像からX線照射角度θもキャリブレーションする詳細型のキャリブレーションとが予め用意されている。このため、コントローラ57は、例えばモニタ60への画像表示を通じて、簡易型キャリブレーションを行うのか、詳細型キャリブレーションを行うのか、オペレータから事前に情報を得る。したがって、画像処理プロセッサ56はコントローラ57からの指示を受けて、キャリブレーションが簡易型であるのか詳細型であるのかを判断するステップS5。画像プロセッサ56が簡易型キャリブレーションを行うと判断した場合、例えば図2に示すように予め決まっている照射角度θ=0°、±15°、±30°、...の値をそのまま読み出して設定するステップS6。これに対し、詳細型キャリブレーションを行うと判断した場合、パノラマ画像からX線照射角度θのずれ、即ち実際の照射角度θ ́を演算する。<処理3>次に、X線照射角度θに対する実際の照射角度投影角度θ ́のずれ量θshiftを演算するステップS7。この演算には、再構成された基準面パノラマ画像上で、基準断層面SSよりも20mm外側の断層面に沿った、外側面位置に沿ってX線照射角度θ毎に配置された支柱のマーカについて、ステップS4と同様に、X線照射角度θ毎のゲインΔX/ΔFi、及び、そのゲインΔX/ΔFiのプロファイルを生成する。このプロファイルのデータを用いて、基準断層面SSよりも20mm外側の断面の外側面パノラマ画像を再構成する。この外側面パノラマ画像において、かかる外側面位置に在る支柱113 ́それぞれのファントムの横方向2次元基準面パノラマ画像上での横方向の物理的な中心位置を決める。この決定もオペレータがパノラマ画像を目視しながら行う。ステップS4において基準面位置に在る支柱それぞれのマーカが使用している中心フレームデータの番号Fi0は既に決定している。そこで、当該中心フレームデータに対する外側面パノラマ画像におけるマーカの横方向2次元基準面パノラマ画像上での横方向の位置図15A参照、及び、外側面位置にあるマーカ116の外側面パノラマ画像における横方向の位置図15B参照から、両者間の画像上のずれ量Pshiftを演算する。このずれ量Pshiftを実長のずれ量Lに変換する図15C参照。このずれ量Lと両軌道ORs、ORouter間の既知の距離DS実施形態では20mmとを用いて、実際の照射角度θ ́のずれ量θshift=arctanの演算を角度θ=0°、±15°、±30°、...毎に行う。これにより、所定値刻みのX線照射角度θそれぞれに対する実際の照射角度θ ́のずれ量Pshiftを求めることができる。このずれ量Pshiftの例を図16に示す。<処理4>次に、画像プロセッサ56は、投影角曲線θ=fFi、すなわち角速度曲線ΔθΔFiを演算するステップS8。詳細型キャリブレーションの場合、既に、各X線照射角度θからの実際の照射角度θ ́のずれ量θshiftが求められている。このため、このずれ量θshiftから基準断層面SSに在るマーカの照射角度θ ́をそれぞれ求める。簡易型キャリブレーションの場合は、簡便的に採用した設計値θがそのまま使用される。一方、前述したステップS4において、基準断層面SSに在る、X線照射角度θ毎のマーカの中心フレームデータFi0の番号は求められている。したがって、照射角度の実際値θ ́又は照射角度の設計値θのそれぞれのフレームデータFiを繋ぎ合わせてスムージングを掛けることで、投影角曲線θ=fFiを求める。この投影角曲線θ=fFiの一例を図17に示す。同図において、投影角曲線θ ́=fFiは、元の投影角曲線θ=fFiから、実際の照射角度θ ́の分だけ補正された曲線を示す。<処理5:X線照射角度θ=0度の位置における定数パラメータの演算>次に、画像プロセッサ56は、X線ビームの照射角度θ=0°のときのX線管回転中心間距離RS、検出器回転中心間距離Rd、及び、X線管の焦点位置の高さ情報B1を定数パラメータとしてパノラマ画像から演算するステップS9。図18に示すように、X線管31及び検出器32が対向して配置されており、その間に、回転中心RC及び基準断層面SSが位置しているとする。基準断層面SSの位置に上下に67mm、互いに離間した2つのマーカ114,115が存在しているとする。X線管31のX線焦点は点光源と見做せるほど小さい焦点例えば0.5mmの径。また、X線の照射角度θは0°であるとする。つまり、コリメータ33で絞られたX線ビームは基準断層面SSに存在すると仮定される歯列の前歯の中心部に照射される。このX線ビームは、2つのマーカ114,115をオブリークに透過して検出器32の検出面の高さB2,B3の位置にそれらの投影点を作る。つまり、マーカ114,115の位置は縦方向Z軸方向に拡大されて画像として投影点B2,B3を作る。なお、検出器32の検出面の最下限の位置を座標0の原点として設定しており、この点を通る水平面XY座標面から起算してX線焦点位置の高さをB1としている。したがって、検出器32の検出面には、座標0の原点、X線焦点高さB1、及び、マーカ114,115の投影高さB2,B3が下から順に並ぶ。この図18に模式的に示すX線照射角度θ=0°における幾何学関係に、前記ゲインの式4:ΔxΔFiを適用する。X線照射角度θ=0°のときは、図13A,Bから分かるように、補正項M=0と見做せる。このため、4式は、ΔxΔFiD...5と表すことができる。マーカ114,115に対する画像の拡大率の計算からB2/67=KD...6が成り立つ。ここで、拡大率KDは、検出器32により、基準断層面の位置であってX線照射角度θ=0°の位置にあるマーカ114,115の投影点B2,B3の位置を検出すれば判るので、既知の値である。同様に、B2/67=KD20...7が成り立ち、拡大率KD20は、検出器32により、基準断層面よりも20mm外側の位置であってX線照射角度θ=0°の位置にあるファントムの投影点B2,B3の位置を検出すれば判るので、既知の値である。このため、上記6、7式=KD...8=KD20...9において、X=RsRd,Y=RsD...10とおくと、X/Y=KD...11X/Y20=KD20...12となり、この2つの式からX,Yの値を求めておく。一方、前記5式は、8式を使えば、ΔxΔFiKDD...13と書き換えることができる。この13式において、回転中心基準断層面間距離D以外の項は既知であるので、13式からX線照射角度θ=0°における回転中心基準断層面間距離Dが分かる。距離Dが既知になると、X,Yの値が既に分っているので、10式を使って、X線照射角度θ=0°におけるX線管回転中心間距離RS、及び、検出器回転中心間距離Rdがそれぞれ求められる。この距離D,Rs,Rdが求められると、図18において幾何学的に成立する、B2D20-B1/H=KD20...14B2D-B1/H=KD...15の2式を解くことで、X線管31の上下方向Z軸方向の位置B1、及び、下側のファントムのX線管31からの高さHを求める。処理6:X線照射角度θ=0度以外の角度位置における、フレームデータFiを入力とする関数パラメータの演算このときには、照射角度θ毎のX線管31、検出器32、回転中心RC、及びファントムマーカの幾何学的な位置関係は、図19のように表される。X線照射角度θが0度以外の角度であっても、前述した式6及び8は成立している。このため、これらの式に基づいて、各角度θにおけるマーカ114,115が位置B2,B3に作る投影像B3、B2を求めることで、各照射角度θにおける回転中心基準断層面距離Dが演算されるステップS10。この距離Dが分かれば、既に既知となっているX線照射角度θ又はその実際値θ ́を用いて回転中心RCの位置座標CX,CYも演算されるステップS11。さらに、X線照射角度θが0°以外の角度のときには、図13A,Bから分かるように、補正項Msup10の考慮が必要である。このため、4式:ΔxΔFiを用いる必要がある。既に、各X線照射角度θにおける補正項M以外の項は演算されているので、それらを4式に当てはめて補正項Mが演算されるステップS12。このようにステップS10~S12を通して、関数パラメータΔxΔFi、θ、ΔθΔFi、D、M、CX,CYが演算される。次いで、画像プロセッサ56は、画像メモリ54に書き込まれている、それらの関数パラメータを今回求めた新しい値で更新するステップS13。これにより、3D画像再構成に必要なパラメータがキャリブレーションされたことになる。以上の構造解析及びキャリブレーションのための演算が終わると、画像プロセッサ56は、演算した定数パラメータRd,Rs,B1及び関数パラメータΔxΔFi、θ、ΔθΔFi、D、M、CX,CYが印刷したり表示したりして出力するか否かを、オペレータの操作情報から判断するステップS14。そのような出力が必要な場合、画像プロセッサ56がそれらの値を印刷又は表示するステップS15。さらに、かかるパラメータの出力が終わるか、又は、そのような出力が不要な場合、処理はコントローラ57に渡され、コンピュータ57がオペレータとの間で、インターラクティブに患者の撮像を行うか否かを判断するステップS16。撮像が不要な場合は、一連の処理を終わる。これにより、撮像空間の構造解析及び簡易型又は詳細型のキャリブレーションが完了する。一方、患者の顎部の撮像を行う場合は、後述するように、撮像空間における歯列の実体位置を正確に把握した3D再構成が実行される。これは、図20に概説するように、3D基準断層面SSから、X線管31を睨むX線の斜めの照射方向に沿って投影が行われ、歯列などの撮像対象実体物の3次元位置が高精度に同定される。以下、この位置同定の処理を含む撮像を説明する。<被検体の撮像>図5に戻って、コントローラ57及び画像プロセッサ56が協働して実行される撮像のための処理を説明する。この処理には、上述したように、スキャンによりデータ収集、プレ処理としての基準面パノラマ画像の再構成、並びに、メインの処理としての3次元オートフォーカス画像3次元表面画像の作成及びその3次元オートフォーカス画像を用いた各種態様に応じた表示や計測などが含まれる。データ収集及び基準面パノラマ画像の再構成まず、コントローラ57は、被検体Pの位置決めなど撮影の準備が済むと、操作器58を介して与えられる操作者の指示に応答し、データ収集のためのスキャンを指令する図5、ステップS1。これにより、回転駆動機構30A、移動機構30B、及び、高電圧発生器41が予め設定されている制御シーケンスに沿って駆動するように指令される。このため、X線管31及び検出器32の対を被検体Pの顎部の周囲に回転させながら、その回転動作の間に、X線管31にパルス状又は連続波のX線を所定周期で又は連続的に曝射させる。このとき、X線管31及び検出器32の対は、前述したようにキャリブレーションされた3D基準断層面SS図6参照を焦点化するように予め設定されている駆動条件に基づいて回転駆動される。この結果、X線管31から曝射されたX線は被検体Pを透過して検出器32により検出される。したがって、前述したように、検出器32から例えば300fpsのレートでX線透過量を反映したデジタル量のフレームデータ画素データが出力される。このフレームデータはバッファメモリ53に一時保管される。このスキャンの指令が済むと、処理の指示は画像プロセッサ56に渡される。画像プロセッサ56は、3D基準断層面SSをX線照射方向のフレーム番号Fi毎に、ルックアップテーブルLUTから照射角度、角速度、回転中心基準断層面間距離D、及び補正項Mの最新値を読み出して、当該3D基準断層面SSを補正する。これにより、この断層面SSが部分的にその前後方向に位置変更されてスムージングされるステップS2A。次いで、この補正された3D基準断層面SSの空間位置に対応したトモシンセシス法に基づくシフトアッドに拠り基準面パノラマ画像PIstを再構成するとともに、その再構成した画像の各画素値を記憶するステップS2B。なお、この再構成処理において、従来と同様に、前歯部の中心で縦横の拡大率の比が同じになるように係数を掛ける処理も実行される。この再構成の仕方は公知ではあるが、若干説明しておく。この再構成に使用するフレームデータのセットは、例えば図21に示すパノラマ画像の横方向の写像位置とその写像位置の画像を作成するために相互加算するフレームデータのセットとの関係を示す写像特性から求められる。この写像特性を示す曲線は、フレームデータ方向横軸において両サイドの臼歯部に応じて傾斜が急な両曲線部分と前歯部に応じて傾斜が臼歯部のそれよりも緩やかな曲線部分とから成っている。この投影特性上で、図示の如く、パノラマ画像の横方向における所望の写像位置を指定する。これに応じて、その写像位置の画像を作成するために使用するフレームデータのセットとそのシフト量重ね合わせの程度:つまり傾斜度が求められる。そこで、それらのフレームデータ画素値をその指定したシフト量を以ってシフトさせながら相互に加算して、指定した写像位置範囲の縦方向の画像データを求める。パノラマ画像の横方向の全範囲に亘って、上記写像位置の指定とシフトアッドを行うことにより、3D基準断層面SSに焦点を当てたときの基準面パノラマ画像PIstが再構成される。画像プロセッサ56は次いで、この基準面パノラマ画像PIstをモニタ60に表示させるステップS3。この基準面パノラマ画像PIstの例を図22に模式的に示す。この基準面パノラマ画像PIstは、フレームデータをシフトさせながら相互に加算した画像であるので、矩形状の2次元画像である。拡大率について言えば、前歯部の中心で縦横の拡大率の比が同じになるように係数を掛ける処理を行っているので、従来と同様に、拡大率に因る前歯部の縦横の画像歪はある程度改善されている。しかし、臼歯部に進むにつれて歯の縦横比は崩れてくる。つまり、臼歯部の歯は実寸より縮んで描出される。従来は、多くの場合、このような歪が在るパノラマ画像で我慢していた。基準面パノラマ画像上でのROI設定次いで、画像プロセッサ56は操作者から操作器58を使って基準面パノラマ画像PIstにROI関心領域が設定するか否かを判断するステップS4。ここで設定するROIは、読影者が特に関心を寄せる例えば矩形状の部分領域である。勿論、ROIは必ずしも矩形でなくてもよい。なお、このROIは、後述するオートフォーカスにより作成したパノラマ画像について設定してもよく、この処理も後述される。このステップS4の判断がYESとなると、画像プロセッサ56は操作者の操作情報に基づいて基準面パノラマ画像PIstにROIを設定するステップS5。次いで、ROIにより設定された部分領域の部分画像を切り出し、その部分画像を例えば拡大して表示するステップS6。この部分画像は、例えば図23に示すように、元の基準面パノラマ画像PIstに重畳して表示される。また、この1つ以上の部分画像を上歯、下歯の歯列の模式的に表すようにブロックを所定順に並べた、いわゆるテンプレートに収めるように表示してもよい。次いで、画像プロセッサ56は処理を終了させるか否かを判断する。この判断は操作者からの所定の操作情報が有るか否かによるステップS7。未だ処理を終了させないと判断した場合ステップS7、NO、ステップS4まで戻って上述した処理を繰り返す。一方、処理終了の判断ができた場合、図5に示す処理を終了させる。一方、画像プロセッサ56は、ステップS4の判断でNOとなる場合、すなわちROIを設定しないと判断した場合、次の判断に移行する。つまり、メインの処理としての3Dオートフォーカス画像を作成するか否かを、操作者の操作情報から判断するステップS8。この作成も行わないと判断した場合ステップS8、NO、ステップS7に戻って処理終了か否かを前述と同様に判断する。最適焦点の断面位置の特定これに対して、3Dオートフォーカス画像を作成すると判断した場合ステップS8、YES、ステップS9のサブルーチン処理に移行する。このステップS9で実行される処理は、本発明の特徴の一つを成すもので、回転中心RCの位置変化を考慮し、且つ、各画素から常にX線管31のX線焦点を睨んだオブリークな投影方向DRxに沿って、歯列の縦方向の歪みを補正しながら行なう、自動的な歯列の実存位置形状の同定処理である。この実在位置形状の同定のためのサブルーチン処理を図24に示す。まず、画像プロセッサ56は、基準面パノラマ画像PIst矩形を3D基準断層面SS湾曲面に平行な湾曲面に座標変換して3Dパノラマ画像を一度、作成する。そして、ルックアップテーブルLUTからフレーム番号Fi毎の照射角度θおよび回転中心の位置座標CX,CYの最新値を読み出し、この位置座標CX,CYからX線管回転中心間距離RSだけ方向を延長し、X線照射角度θ毎にX線管31の位置を演算する。そして、作成した3Dパノラマ画像の画素それぞれから常にX線管31のX線焦点を睨んだオブリークな投影方向DRxを決める。その上で、各投影方向DRxに沿って3D基準断層面SSに、断層面変更の演算によりフレームデータを求め、これを座標変換することで投影し、その湾曲した3D基準断層面SSの投影画像を作成するステップS51。この投影像の画素値は画像メモリ54に保管される。ここで行われる投影は、図25に説明するように、回転中心RCRC1、RC2の位置、すなわちX線管31の位置に向けたオブリークな投影方向DRxに沿って行われる。図25の例で言えば、3Dパノラマ画像上の高さ方向Z軸方向における同じ位置Pnの画素であっても、X線管31の位置の違いによって3D基準断層面SSの画像上の異なる位置SS1、SS2に投影される。この投影処理により作成される投影画像を3D基準画像PIrefと呼ぶことにする。この3D基準画像PIrefは、基準面パノラマ画像PIstの部位毎に、前述した縦方向の拡大率を考慮した斜め方向の投影によって作成されている。前歯部の歯の拡大率が大であったものが、その拡大は上述の投影により実サイズに是正され、一方、臼歯部の歯の拡大率が小であったものが、その拡大も上述の投影よりに実サイズに是正される。このため、3D基準画像PIrefは歯の実寸で表示された画像であり、スキャン中に回転中心RCが移動することによる拡大率の大小による歪が除去された画像である。ただし、この3D基準画像PIrefは歯列が3D基準断層面SSに沿って存在すると仮定したときの画像でもある。被検体Pの実際の歯は3D基準断層面SSに沿っていることは稀であるので、後述する更なる実在位置の同定処理が必要になる。画像プロセッサ56は、その3D基準画像PIrefをモニタ60に表示させ、操作者の参照に供するステップS52。この様子を図26に示す。この後、画像プロセッサ56は、3D基準断層面SSに、その面に平行な複数の湾曲した断層面を付加するステップS53。この様子を図27に示す。同図には、3D基準断層面SSの投影方向DRx歯列の奥行き方向の前後それぞれに複数の断層面が付加されている。この複数の断層面のそれぞれも、3D基準断層面SSが回転中心基準断層面間距離D及び補正項Mで補正された分だけ、その面の前後方向の位置が部分的に補正された断層面である。一例として、3D基準断層面SSの前側に複数の断層面SFm~SF1を間隔D1例えば0.5mmで設定し、その後側に複数の断層面SR1~SRnを間隔D2例えば0.5mmで設定している。間隔D1、D2は同じであっても、互いに相違していてもよい。また、付加する断層面は、3D基準断層面SSの前後に1枚ずつm、n=1であってもよいし、前後の何れかに1枚又は複数枚であってもよい。なお、この仮想的に付加する断層面SFm~SF1、SR1~SRnの位置データは、3D基準断層面SSの位置データと共に予めROM61に格納されているので、これを画像プロセッサ56のワークエリアに読み出すことで、かかる付加が実行される。断層面SFm~SF1、SS、SR1~SRnの高さは投影方向DRxの最大の傾きと歯列の高さとを考慮して適宜に設定されている。また、同定処理の都度、付加する断層面の位置間隔D1、D2及び枚数をインターラクティブに変更するようにしてもよい。次いで、画像プロセッサ56は、ステップS51で行ったと同様に、回転中心RCの位置座標CX,CYの変化に応じた投影方向DRxを求め、その投影方向DRXに沿って、基準面パノラマ画像PIstを、付加した断層面SFm~SF1、SR1~SRnそれぞれに、断層面変更の演算によりフレームデータを求めて、これを座標変換することで投影するステップS54。この結果、付加断層面SFm~SF1、SR1~SRnそれぞれの投影画像が作成される。これらの投影像の画素値は画像メモリ54に保管される。ここで作成される投影画像を3D付加画像PIsfm...PIsf1PIsr1...PIsrnと呼ぶ。これらの3D付加画像PIsfm...PIsf1PIsr1...PIsrnも、それぞれ、基準面パノラマ画像PIstの位置毎に、前述した縦方向の拡大率を考慮した斜め方向の投影によって作成されている。これを図28の例で言えば、3Dパノラマ画像上の高さ方向Z軸方向における同じ位置Pnの画素であっても、X線管31の位置の違いによって3D付加画像PIsfm...PIsf1PIsr1...PIsrnそれぞれの上で異なる位置に投影される。このため、これらの3D付加画像PIsfm...PIsf1PIsr1...PIsrnも歯の実寸で表示された画像であり、スキャン中に回転中心RCが移動することによる拡大率の大小による歪が除去された画像である。ただし、これらの3D付加画像PIsfm...PIsf1PIsr1...PIsrnは歯列がそれぞれの付加断層面SFm~SF1、SR1~SRnに沿って存在すると仮定したときの画像でもある。なお、この作成された複数枚の3D付加画像PIsfm...PIsf1PIsr1...PIsrnはそのまま3次元画像として、又は、座標変換した上で長方形状の2次元画像としてモニタ60に表示させるようにしてもよい。この後、画像プロセッサ56は3D基準画像PIref、すなわち3D基準断層面SSにおける初期位置Pを指定するステップS55:図29A参照。これが済むと、3D基準画像PIrefにおいて、指定した位置Px、y、zを中心とする一定長さの線分Lcを指定するステップS56:図29B参照。この線分Lcは2n個n=1,2,3、...;例えば128分の画素に相当する長さを有する。なお、線分Lcは湾曲する3D基準断層面SSの一部に沿って湾曲していてもよいし、直線と見做せる範囲で設定してもよい。次いで、画像プロセッサ56は、指定された線分Lcx、y、zの画像上の上下に複数本の同一長さの線分Laddを仮想的に付加するステップS57:図29C参照。さらに、上述した線分L及び複数の線分Laddのそれぞれを構成する2n個分の画素それぞれの画素値Pijを画像メモリ54から読み出し、これを各線分に割り当てるステップS58。この画素値Pijは、前述したステップS51,S54で既に取得して保管していた値である。次いで、複数の線分L及びLaddの対応する画素の画素値Pij同士を加算して、線分Lcx、y、zを構成する周波数解析用の2n個の画素値Pijを求めるステップS59:図29D参照。この加算より、線分Lx、y、zの元の画素値に統計的ノイズが混入している場合でも、その画素値の変化について後述する周波数解析を行なうときの統計的ノイズを低減させることができる。次いで、画像プロセッサ56は、付加した3D付加画像PIsf1...PIsf1PIsr1...PIsrnのそれぞれにおいて、上述の3D基準画像PIref上で現在指定されている線分Lcが、現在指定されている位置P及びその長さ、並びに、スキャン中のX線管31の回転位置が分っているので、線分Lcの両端とX線管31との結んでできる、Z軸方向から見たときに扇状となるX線照射範囲RAを演算できる。このため、位置Pが指定されれば、そのX線照射範囲RAに位置する線分Lfm~Lf1、Lr1~Lrnの位置を特定できる。なお、3D基準画像PIref上に位置Pを指定するステップS60の処理は全部の位置指定が終わるまで繰り返される。このため、実効的には、仮想した断層面SFm~SF1、SS、SR1~SRnを、位置が遠近するX線管31から照射されたX線は範囲H1~H2Z軸方向の範囲で扇形に透過していることになる図29F。このため、断層面SFm~SF1、SS、SR1~SRnそのものを、その高さがスキャン方向毎に変わり且つ互いに平行な略馬蹄形の断面にとして設定してもよい。上述のように線分Lfm~Lf1、Lr1~Lrnが決まると、画像プロセッサ56は、それらの線分の画素値Pijを画像メモリ54から読み出すステップS61。図29Eに示すように、X線管31は点源であるから、X線照射範囲RAは扇状Z軸方向から見たときにになっている。このため、線分Lfm~Lf1、Lr1~Lrnそれぞれの画素数は2n個からずれてしまっている。そこで、画像プロセッサ56は、付加した線分Lfm~Lf1、Lr1~Lrnの画素数が基準となる線分Lc=Pを投影方向DRxに貫く複数の断層面SFm~SF1、SS、SR1~SRn位置を横軸としたプロファイルとして演算するステップS65。このプロファイルの一例を図31に示す。同図において断面位置とは、複数の断層面SF1~SFm、SS、FR1~FRnの投影方向DRx歯列の奥行き方向の位置である。図32には、物質がエナメル質、海綿骨、空気、バイトブロックである場合の複数種のプロファイルPR1,PR2,PR3PR4の典型的なパターンが例示されている。仮に、現在指定している位置Pを通る投影方向DRxの何れかの位置にエナメル質の物質、すなわち歯が存在している場合、そのプロファイルPR1はシャープなピークを有する。また、かかる投影方向DRxに海綿骨が存在している場合、そのプロファイルPR2はなだらかな凸曲線となる。同様に、かかる投影方向DRxに空気しか存在している場合、そのプロファイルPR3は特定のピークを持たない傾向を示す曲線となる。さらに、かかる投影方向DRxにバイトブロックが存在している場合、そのプロファイルPR4は、2つのシャープなピークを有する。このうち、投影方向DRxの内側X線管の側に相当するピークがエナメル質の物質に対するピークを示し、外側検出器の側に相当するピークがバイトブロックに対するピークを示す。図32に示すプロファイルPR1~PR4のパターンを示すデータは、参照プロファイルとして、例えばROM61に参照テーブルとして予め記憶されている。そこで、画像プロセッサ56は、かかる参照テーブルを用いて、現在指定している位置Pを通る投影方向DRxにおける、歯に対する最適焦点の位置を特定するステップS66。つまり、前のステップS65で求めたプロファイルが参照プロファイルPR1~PR4の何れに該当するのか、パターン認識の手法で判断する。まず、求めたプロファイルが参照プロファイルPR2、PR4である場合には処理の対象から外す。一方、求めたプロファイルが参照プロファイルPR1エナメル質に該当する場合、そのピークを呈する断面位置、すなわち、複数の断層面SF1~SFm、SS、FR1~FRnのうちの何れかの位置が最適焦点であるとして特定する。さらに、求めたプロファイルが参照プロファイルPR4に該当する場合、その内側X線管の側にピークを呈する断面位置エナメル質の位置、すなわち、複数の断層面SFm~SF1、SS、FR1~FRnのうちの何れかの位置が最適焦点であるとして特定する。これらの位置の特定処理により、いま指定している位置Pに在る歯の部分が、実際は、奥行き方向のどの位置に在るかを決めたことになる。つまり、3D基準断層面SS上に沿った3D基準画像PIrefに描出された歯の部分は実際には、その断層面SSの前側に在るかもしれないし、後側に在るかもしれない。この実在位置が上述の特定処理により正確に決定される。別の言い方をすれば、3D基準断層面SS上に在ると仮定して描出された3D基準画像PIrefの歯の部分が、上述の特定処理により、実在する位置にシフトされると言える。この結果、図33~図36に示すように、位置Pの1回の指定毎に、3D基準断層面SS3D基準画像PIrefにおける位置P1がP1realまたはP2がP2realにシフトされる。とくに、複数の付加断層面SFm~SF1、FR1~FRnに設定する線分Lfm~Lf1、Lr1~Lrnの位置が投影方向DRxのオブリーク角度θを考慮して設定されている。このため、シフトされる位置P1realは、オブリーク角度θが小さい場合図34A、図35A参照よりも大きい場合図34B、図35B参照の方が低くなる。したがって、このシフト位置P1realは、オブリークなX線照射角度θ、すなわち拡大率の大小による歪みが補償されている。なお、図36に示すように、歯が3D基準断層面SSに沿って実在する場合、P1=P1realとなって、歯が位置するものと仮定していた3D基準断層面SSが実在位置として決まる。この場合はシフト量=0のシフトが実行されたことになる。画像プロセッサ56は、ステップS65において、これらの特定した、歯の実在位置を示すデータを位置P毎に、そのワークエリアに記憶する。このようにして、3D基準画像PIrefすなわち3D基準断層面SSで現在指定されている位置Pか否かで判定することで行う。この判断がNOとなって、全ての判断位置Pについて特定処理が完了していない場合、画像プロセッサ56は、その判断位置Pを1つ分シフトさせステップS68、その処理を前述したステップS55に戻し、上述した一連の特定処理を繰り返す。なお、図37に示すように、複数の判断位置Pは3D基準画像PIrefすなわち3D基準断層面SSに沿って所定間隔を以って2次元的に予め配置されている。同図の例では、3D基準画像PIrefの縦軸方向i及び横軸方向jに沿って縦横同一の所定間隔dを空けて配置されている。ただし、この所定間隔dは縦軸方向i及び横軸方向jそれぞれにて互いに相違させてもよい。ステップS68の処理におけるシフトの方向は、3D基準画像PIrefに沿った縦、横、及び斜めの何れの方向であってもよい。図37に示すように、3D基準画像PIrefの縦軸方向iに沿ってシフトさせた後、横軸方向jにシフトしてまた縦軸方向iに沿ってシフトさせることを規則正しく繰り返してもよい図の符号SCを参照。その逆に、横軸方向jにシフトさせて後、縦軸方向iにシフトさせることを繰り返してよい。さらに、斜め方向にシフトさせてもよい。その一方で、複数の判断位置Pの全てにおいて上述した一連の判断が終了すると、上述した繰り返し判断の中で前述したステップS67における判断がYESとなる。つまり、3D基準断層面SSの奥行き方向における判断位置P毎に最適焦点の断面位置の検出最適焦点位置の有無の判断を含むの処理が完了したことになる。この場合、最適焦点の断面位置の結合処理に移行する。最適焦点の断面位置を結合する処理上述したステップS67の判断がYESとなると、画像プロセッサ56はステップS65において特定し記憶していた最適焦点の断面位置を表すデータを読み出すステップS69。この断面位置のデータは、それぞれの判断位置Pを通る投影方向DRxの位置である。この様子を図38に模式的に示す。同図において、黒丸は3D基準画像PIref3D基準断層面SSの判断位置Pを示す。ここで、湾曲した3D基準画像PIrefの縦方向及び横方向をと表す。図38において、白丸で示す如く、例えば、i,j=0,0の判断位置Px00、y00、z00に対する最適焦点断面位置が内側X線管の側に1つ寄った断層面SR1の位置であり、その隣のi,j=0,1の判断位置Px01、y01、z01に対する最適焦点断面位置が内側さらに1つ寄った断層面SR2の位置であり、その隣のi,j=0,2の判断位置Px02、y02、z02に対する最適焦点断面位置が内側さらに1つ寄った断層面SR3の位置であり、といった具合になる。なお、図38は、図を見易くするため、Z軸方向縦方向の1つの位置におけるステップS68を示しているが、このZ軸方向の他の位置それぞれについてもステップS68の処理が実行される。次いで、画像プロセッサ56はノイズの除去を行うステップS70。図38の例で例えば、画像の縦横方向の位置i,j=0,3の判断位置Px03、y03、z03に対する最適焦点断面位置が外側検出器の側にm個も寄った断層面SFmの位置であるとする。このような場合、画像プロセッサ56は、断面位置同士の差分を例えば閾値判断に掛けてノイズであり異常であると見做す。この場合、隣接する断面同士の位置のデータを滑らかに繋がるように例えば平滑化し、その平滑化した新たな位置データに置換する、又は、選択的に検出器の外側に近いデータを優先させる、などの処理を行う。なお、このような置換による補償を行わずに、単に、異常データを処理対象から外すようにしてもよい。この異常データの排除にZ軸方向のデータの異常を加味することも当然可能である。この後、画像プロセッサ56は、ノイズ除去した位置すなわちエナメル質の位置を結合し、この結合した位置のデータを3次元的にスムージングして、エナメル質の部分の形状をトレースした表面画像を作成するステップS71。さらに、この画像プロセッサ56は、この表面画像を、その部位全てが自動的に最適焦点処理に付された3次元パノラマ画像、すなわち3Dオートフォーカス画像PIfocusとしてモニタ60に所定のビュー角度で表示させるステップS72。これにより、図39に示すように、所定のビュー角度で見た、被検体Pの口腔部の歯列の構造体が最も明瞭に見える輪郭に沿ってできる3Dオートフォーカス画像PIfocusを提供できる。同図において、湾曲している馬蹄形の範囲Sは、3Dオートフォーカス画像PIfocusを表示するための範囲であり、実線部分が歯列の実際の位置及び形状を表している。A-A ́線及びB-B ́線で示す如く、歯茎歯槽骨の部分や下顎洞、顎関節、頚動脈などは、歯主にエナメル質の端部から一定距離にした断層距離をキープし、断層面を作り3D断層面投影する方法も可能である。この場合は、これらの部位が最適焦点になっていることは保証できないが、3Dのパノラマ画像としては、違和感を覚えない画像として再構成可能である。勿論、これらの部位も最適焦点面の計算に工夫を加え、そのまま計算し用いる方法も、診断の目的によっては有り得るのは言うまでもない。このように、3Dオートフォーカス画像PIfocusは、歯列に沿って湾曲しながらも、その表面はでこぼこしており、このでこぼこにより個々の歯の実際の位置及びその形状輪郭を画素の濃淡で表している。その他の部分も違和感のない画像として表現できる。このように個々の被検体Pの歯列の実在位置形状を表す3Dオートフォーカス画像PIfocusが作成される。種々の表示処理この後、画像プロセッサ56は、その3Dオートフォーカス画像PIfocusを他の態様で観察する機会を操作者に与える。つまり、画像プロセッサ56は、操作者から操作情報に基づいて、その3Dオートフォーカス画像PIfocusを他の態様でインターラクティブに表示するか否かを判断する。その一例として、画像プロセッサ56は、3Dオートフォーカス画像3次元パノラマ画像PIfocusの部分領域を観察するか否かを判断する図5、ステップS10。このステップS10の判断がYESになると、さらに、その部分領域の観察を3D基準断層面SSで行うのか、又は、基準面パノラマ画像の矩形面2次元で行うのか、操作者からの情報を基づいて判断するステップS11。このステップS11において3D基準断層面SSを使用すると判断されると、画像プロセッサ56は、3Dオートフォーカス画像PIfocusを3D基準断層面SSに、個々の画素を通る投影方向DRxに沿って再投影するステップS12。この再投影の様子を図40に示す。この再投影は例えば3D基準断層面の一画素を、対応する3次元の画素をサブピクセルで区切り再投影するサブピクセル法により実行される。この3D基準断層面SSへの再投影像は、3D参照画像PIproj3Dとして、モニタ60に表示されるステップS13。この3D参照画像PIproj3Dの一例を図41に示す。一方、ステップS11において基準面パノラマ画像PIstの矩形面を使用すると判断されると、画像プロセッサ56は3Dオートフォーカス画像PIfocusをその矩形面、つまり基準面パノラマ画像の面に再投影するステップS14。この再投影も、標準パノラマ画像面の一画素を、対応する3次元の画素をサブピクセルで区切り再投影するいわゆるサブピクセル法により実行されるのは言うまでもない。この再投影の概念を図42に示す。この再投影像は、2D参照画像PIproj2Dとして、モニタ60に表示されるステップS15。この2D参照画像PIproj2Dの一例を図43に示す。そこで、操作者は、この3D参照画像PIproj3Dまたは2D参照画像PIproj2Dに所望の、例えば矩形のROI関心領域を設定するステップS16:図41及び図43を参照。このROIにより指定された部分領域の画像は例えば拡大され、例えば現在表示されている3D参照画像PIproj3Dまたは2D参照画像PIproj2Dに重畳表示されるステップS17。勿論、この表示は、パノラマ画像とは別個の単独画像であってもよいし、同パノラマ画像との分割表示であってもよいし、歯列を模した複数のブロックから成るテンプレートの1つに収めた表示であってもよい。この後、画像プロセッサ64はかかる一連の処理を終了するか否かを操作情報から判断しステップS18、この判断がYESの場合は処理を前述したステップS7に戻す。これに対し、NOの場合は処理をステップS10に戻して上述した処理を繰り返す。その一方で、前述したステップS10において部分画像の観察をしないと判断する場合、画像プロセッサ56は、現在表示されている3Dオートフォーカス画像PIfocusを回転、移動、及び/又は拡大縮小して表示するか否かをインターラクティブに判断するステップS19。この判断がYESのなる場合、指令情報に応じて3Dオートフォーカス画像PIfocusを回転、移動、及び/又は拡大縮小し、その画像を表示するステップS20,S21。この後、処理はステップS81に移され、前述と同様の処理を繰り返す。勿論、表示態様の種類は上述したものに限定されず、例えばカラー化など、その他の様々な態様を採り得る。操作者が処理の終了を指示している場合、画像プロセッサ64はステップS18、S7を経て、かかる処理を終了させる。なお、上述したステップS16の設定処理を行った後、ステップS17の表示処理を行わずに、ステップS19の処理に移行するようにしてもよい。その場合、設定したROIは、回転、移動、拡大縮小した画像と共にステップS21において表示される。以上のように、本実施形態によれば、パノラマ撮像空間の構造を3次元的に把握することで、投影方向が3次元的に表現できる。従って、パノラマ画像の焦点が合っている限りは、3次元表現された画像に歪が生じず、正確なパノラマ撮影画像を構築することができる。このことにより、パノラマ画像をより、位置決めの良否に関わらず安定に表示でき、かつパノラマ画像全体で鮮明な画像を作るようなこともできる。また、本実施形態に係るファントムを用いて撮像空間の構造解析及び3D画像再構成に必要なパラメータをキャリブレーションすることができる。キャリブレーションも、X線照射角度θを対象としない簡易型から、X線照射角度θを考慮した詳細型のキャリブレーションまで選択的に実行できる。このため、装置毎に、そして、医療施設に据え付けた後にも、必要性に応じて適宜なタイミングでキャリブレーションを簡単に行える。しかも、この構造解析やキャリブレーションはX線管31及び検出器32の対の回転中心RCの位置変化αθ1を考慮している。3D画像再構成は、その位置変化αθ1を補償したパラメータをも参酌して実行される。したがって、装置毎のパラメータの個体差や経時的なパラメータの変化が確実に補償されるので、歯列の実体位置が正確に描出された、つまり、3次元の距離の精度が高い、3次元表面画像が得られる。一方、基準断層面の軌道が不明のパノラマ撮像装置に対しても、ファントム101を用いた撮像を通して、その撮像空間を特徴付ける定数パラメータ及び関数パラメータを計測することもできる。ファントムの変形例なお、本発明に係る放射線撮像装置で採用可能なファントムは上述したものに限定されず、上述以外にも様々な構成で実施することができる。以下、この変形例を説明する。前述したファントム101には、基準面軌道ORsとは異なる他の軌道として、この基準面軌道ORsよりも外側の軌道である外側面軌道ORouterを設けるようにしたが、基準面軌道ORsとよりも内側に所定距離だけ離間させて設けた軌道であってもよい。また、ファントム101は、異なる2つの軌道ORs、ORouterに植設する支柱113、113 ́は、生産性を向上させる上で有利なように、共に同じ形状及び寸法のものにし、組み立て時に一方を逆さまに設置するだけの構造にした。しかしながら、X線照射角度θ=0以外の角度に設置する支柱113、113 ́のうち、一方、例えば外側面軌道ORouterに設置するファントム113 ́は、ゲインΔX/ΔFiおよびX線照射角度θを測定できれば足りるので、第1及び第2のマーカ114,115を外すようにしてもよい。この簡略化した一例を図44に示す。さらに、支柱113,113 ́については、前述した実施形態では細長い角柱状に形成したが、細長い丸棒201であってもよい図45A参照。また、支柱113、113 ́の代わりに、紐状の懸垂体202を用意し、その懸垂体の途中にマーカを設置する構造であってもよい図45B参照。第1、第2、及び第3のマーカ114,115,116の夫々についても、必ずしもライン状でなくてよく、複数の点がライン状に並んだ点列のマーカであってもよい図45C参照。特に、第1及び第2のマーカ114,115は上下方向の位置が分かればよいので、1つの点状のマーカにしてもよい図45D参照。さらには、第1、第2、及び第3のマーカ114,115,116を1つのマーカ体で構成することもできる。例えば図32において、第1のマーカ114と第2のマーカ115を上下に繋ぐ細長いマーカ板203を設け、そのマーカ板203の第3のマーカ116に相当する部分にスリットSを形成する図45E参照。このマーカ板の上下端位置が第1、第2のマーカの機能を果たし、且つ、スリットSが第3のマーカ116の機能を果たすので、このマーカ板を支柱113の所定位置に貼り付けるようにしてもよい。さらに、マーカ114~116はその他のファントムの部位とはX線透過率が異なればよいので、例えば支柱113,113 ́を金属性とし、マーカに相当する部分をX線透過率が異なる物質単なる凹みも含むで形成することもできる。このように、必要な位置や距離の情報を得るためにファントム101に与えられるべき構造の変形例は数多く存在することは言うまでもない。なお、すべての情報を一度に測定可能な構造のほか、必要な情報を個別に分けて収集できる複数のファントムを用意してもよく、そのようなファントムも本発明の主旨に含まれることは言うまでもない。例えば、前述した万能型ファントム101の基準面軌道ORsのみに前述した一組の支柱体113を植設した第1のファントムと、外側面軌道ORouter或いは、図示していなが内側面軌道のみに前述した一組の支柱体113 ́を植設した第2のファントムとを用意することができる。この場合、第1及び第2のファントムについて2回のスキャンを実行し、それぞれのスキャン結果から2つのパノラマ画像を再構成し、前述した測定を行なうことができる。さらに、これらの第1及び第2のファントムは、物理的に前述した1つの万能型ファントム101に結合可能な構成にすることもできる。この場合には、第1及び第2のファントム自体は個別に保管し、測定時には1つのファントムとして使用できる。なお、本発明は上述した実施形態のものに限定されず、本発明の要旨を逸脱しない範囲でさらに適宜に変形可能なものであり、それらも本発明に含まれることは言うまでもない。例えば、X線管31と検出器32を被検体Pの口腔部の周りに回転移動させるときに、X線管31と検出器32とが互いに斜めに対向しながらX線によるスキャンが行われるように、回転ユニット24の構造を変更してもよいし、かかる互いに斜めに対向したスキャンと実施形態で説明した互いに正対したスキャンとを適宜に組み合わせてもよい。この組み合わせは、口腔部の歯列のどの部位をスキャンするかに応じて設定することができる。これにより、常に歯列の各部位に極力、直交した状態でX線をスキャンさせ、歯の重なりや頸椎の写りこみなどを極力、回避した角度でX線照射を行い、アーチファクトなどのより少ないパノラマ画像を再構成させることができる。ところで、本発明に係る放射線撮像装置は、歯科用のパノラマ撮像装置に実施するものに限定されず、トモシンセシス法を用いて対象物の実体形状位置を3次元的に把握する装置にも実施することができる。そのような応用として、例えば医療用としては、トモシンセシス法を用いたマンモグラフィ、肺がん検査用スキャナへの用途がある。"