(57)【特許請求の範囲】
【請求項１】
ポジトロン放射断層撮影（ＰＥＴ）において被検体を載せる天板である被検体天板が連続して移動する際の前記被検体天板のたわみを検出し、当該検出されたたわみに基づくたわみ情報を提供する測定サブシステムと、
前記測定サブシステムから前記たわみ情報を受信し、ポジトロン放射断層撮影（ＰＥＴ）スキャナから複数の同時発生事象に対応するＰＥＴ測定データを受信し、当該受信したたわみ情報及びＰＥＴ測定データを再構成サブシステムに伝達する収集サブシステムと、
を備え、
前記再構成サブシステムは、前記受信したたわみ情報を用いて、前記検出されたたわみに応じた変換行列を生成し、
ＰＥＴスキャン画像を再構成するのに使用される再構成空間に前記変換行列の転置を適用して補正済み再構成空間を生成し、
前記補正済み再構成空間に基づいて、前記ＰＥＴ測定データからＰＥＴスキャン画像を再構成するプロセッサを有する、ポジトロン放射断層撮影（ＰＥＴ）イメージングシステム。
【請求項２】
前記測定サブシステムは、前記被検体天板の所定の位置において検出した少なくとも１つのたわみ量に基づいて、前記被検体天板の前記たわみを推定する、請求項１に記載のＰＥＴイメージングシステム。
【請求項３】
前記測定サブシステムは、複数のブラッグ反射回折格子を含む１本の光ファイバから反射光を受信し、前記光ファイバから受信した前記反射光の波長に基づいて、前記被検体天板の縦方向の複数位置におけるたわみ量を決定する、請求項２に記載のＰＥＴイメージングシステム。
【請求項４】
前記測定サブシステムは、１つのブラッグ反射回折格子をそれぞれが含む複数本の光ファイバから反射光を受信し、前記複数本の光ファイバから受信した前記反射光の波長に基づいて、前記被検体天板の縦方向の複数位置におけるたわみ量を決定する、請求項２に記載のＰＥＴイメージングシステム。
【請求項５】
前記測定サブシステムは、１本の光ファイバ又は複数本の光ファイバから光を受信し、前記１本の光ファイバ又は複数本の光ファイバそれぞれから受信した光の、前記被検体天板における異なるたわみ量に関する少なくとも１つの特性を測定し、前記１本の光ファイバ又は複数本の光ファイバそれぞれから受信した前記光の少なくとも１つの特性と前記被検体天板の前記異なるたわみ量との間の対応関係を規定するテーブルを生成する、請求項２に記載のＰＥＴイメージングシステム。
【請求項６】
ポジトロン放射断層撮影（ＰＥＴ）イメージングシステムが行なうＰＥＴスキャン画像を再構成するためのポジトロン放射断層撮影（ＰＥＴ）イメージング方法であって、
測定サブシステムが、ポジトロン放射断層撮影（ＰＥＴ）において被検体を載せる天板である被検体天板が連続して移動する際の前記被検体天板のたわみを検出する検出ステップと、
前記測定サブシステムが、前記検出ステップにより検出されたたわみに基づいてたわみ情報を提供する提供ステップと、
収集サブシステムが、前記測定サブシステムから前記たわみ情報を受信し、ポジトロン放射断層撮影（ＰＥＴ）スキャナから複数の同時発生事象に対応するＰＥＴ測定データを受信する受信ステップと、
再構成サブシステムが有するプロセッサが、前記受信ステップにより受信されたたわみ情報及びＰＥＴ測定データを用いて、ＰＥＴスキャン画像を再構成する再構成ステップと、
を含み、
前記再構成ステップは、前記検出されたたわみに応じて変換行列を生成し、ＰＥＴスキャン画像を再構成するのに使用される再構成空間に前記変換行列の転置を適用して補正済み再構成空間を生成し、前記補正済み再構成空間に基づいて、前記ＰＥＴ測定データから前記ＰＥＴスキャン画像を再構成する、ポジトロン放射断層撮影（ＰＥＴ）イメージング方法。
【請求項７】
前記検出ステップは、前記被検体天板の所定の位置において検出した少なくとも１つのたわみ量に基づいて、前記被検体天板の前記たわみを推定する、請求項６に記載のＰＥＴイメージング方法。
【請求項８】
前記検出ステップは、
１つ又は複数のブラッグ反射回折格子を含む光ファイバから反射光を受信し、
前記光ファイバから受信した前記反射光の波長に基づいて、前記被検体天板の縦方向の複数位置におけるたわみ量を決定する、請求項７に記載のＰＥＴイメージング方法。
【請求項９】
前記検出ステップは、
１つのブラッグ反射回折格子をそれぞれが含む複数本の光ファイバから反射光を受信し、
前記複数本の光ファイバから受信した前記反射光の波長に基づいて、前記被検体天板の縦方向の複数位置におけるたわみ量を決定する、請求項７に記載のＰＥＴイメージング方法。
【請求項１０】
前記測定サブシステムが、
１本の光ファイバ又は複数本の光ファイバから光を受信し、
前記１本の光ファイバ又は複数本の光ファイバそれぞれから受信した光の、前記被検体天板における異なるたわみ量に関する少なくとも１つの特性を測定し、
前記１本の光ファイバ又は複数本の光ファイバそれぞれから受信した前記光の少なくとも１つの特性と前記被検体天板の前記異なるたわみ量との間の対応関係を規定するテーブルを生成する、ことを更に含む、請求項７に記載のＰＥＴイメージング方法。
【請求項１１】
核医学イメージングにおいて被検体を載せる天板である被検体天板が連続して移動する際の前記被検体天板のたわみを検出し、当該検出されたたわみに基づくたわみ情報を提供する測定サブシステムと、
前記測定サブシステムから前記たわみ情報を受信し、スキャナから複数の放射線測定データを受信し、当該受信したたわみ情報及び放射線測定データを再構成サブシステムに伝達する収集サブシステムと、
を備え、
前記再構成サブシステムは、前記受信したたわみ情報を用いて、前記検出されたたわみに応じた変換行列を生成し、
核医学画像を再構成するのに使用される再構成空間に前記変換行列の転置を適用して補正済み再構成空間を生成し、
前記補正済み再構成空間に基づいて、前記放射線測定データから核医学画像を再構成するプロセッサを有する、核医学イメージングシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、ポジトロン放射断層撮影（ＰＥＴ）イメージングシステム、方法及び核医学イメージングシステムに関する。
【背景技術】
【０００２】
ポジトロン放射断層撮影（ＰＥＴ：Positron  Emission  Tomography)イメージングは、医用イメージングの分野で利用が広がっている。ＰＥＴイメージングにおいては、最初に患者に放射性医薬品が投与される。かかる放射性医薬品は、多くの場合、注射で患者に注入されるが、吸入又は摂取させることも可能である。放射性医薬品を投与してしばらくすると、かかる薬品は、当該薬品の物理的性質及び生体分子的性質により、人体内の特定の部位に集積する。薬品の実際の空間分布、蓄積点若しくは蓄積領域の強度、及び投与から捕獲、そして最終的な排出に至るプロセスの動態は、全て、臨床的に重要な意味を持つ可能性のある要素である。このプロセスを通して、放射性医薬品に付着させたポジトロン放出体は、半減期、分岐比といった同位体の物理的性質に従ってポジトロン（陽電子）を放出する。
【０００３】
放出されたポジトロンが電子と衝突する際に消滅事象（annihilation  event）が起こり、その結果、ポジトロン及び電子は崩壊する。ほとんどの場合、消滅事象によって、実質的に１８０°離れて反対方向に移動する２つのガンマ線（５１１ｋｅＶ）が発生する。
【０００４】
これら２つのガンマ線を検出し、それらの部位を結ぶ線、即ち同時計数線（ＬＯＲ：Line  Of  Response）を引くことによって、高い確率で元の消滅位置を探し出すことができる。このプロセスは、相互作用が起こっている可能性のある線を識別するに過ぎないが、これらの線を数多く蓄積すれば、断層を再構成するプロセスを用いることによって、元の分布を推定することができる。シンチレータ結晶では２つのガンマ線の相互作用から２つのシンチレーション事象が生じており、シンチレーション事象が生じた部位の他に、正確なタイミング（数百ピコ秒以内）も利用可能であれば、飛行時間（ＴＯＦ：Time  Of  Flight）の計算により、上記の同時計数線に沿った、消滅事象が生じた確率の高い位置に関する更なる情報を加えることができる。スキャナが有するタイミング分解能の限界によって、この線に沿った位置判定の精度が決まることになる。更に、元のシンチレーション事象の部位を決定する際の限界によって、スキャナの最終的空間分解能が決まることになる。一方、同位体の特定の特性（例えばポジトロンのエネルギー）も、（ポジトロンの範囲及び２つのガンマ線の共直線性を介して）この特定の薬品の空間分解能を決める一因となる。
【０００５】
多数の事象を集めれば、断層再構成から、推定すべき被検体の画像に関する必要な情報を生成できる。該当の検出エレメントで検出された実質的に同時に起こっている２つの事象によって、同時計数線が形成される。それらの事象の幾何学的属性に従ってこの線をヒストグラム化して、再構成すべき投影図、即ちサイノグラムを規定することができる。事象は、画像に個別に加えることも可能である。
【０００６】
したがって、データを収集し画像を再構成する基本的要素は、システム−被検体開口（system-patient  aperture）を横切る線であるＬＯＲである。他にも、事象の部位に関する情報を得ることができる。第１に、点の再構成又は配置を行うシステムの能力は、有効視野（ＦＯＶ：Field  Of  View）全域で空間的に不変なのではなく、中央の方が良好なのであって、周辺に向かうに従って緩やかに劣化することがサンプリングや再構成で明らかになっている。この挙動の特徴を表すのに、典型的には、点拡がり関数（ＰＳＦ：point-spread-function）を使用する。そのＰＳＦを再構成プロセスに組み込むために、いくつかのツールが開発されている。第２に、ＬＯＲに沿ったどの場所の方が、事象が起こった確率がより高いかを決定するのに、飛行時間、又は、対の検出に伴って生じる各検出器へのガンマ線の到達時間差を使用することができる。
【０００７】
上記の検出プロセスは、多数の事象について反復する必要がある。イメージング作業を支持するのに計数（即ち事象の対）がどれだけ必要になるのかを決定するには、各イメージング事例を解析しなければならないものの、数億の計数を蓄積する必要があるというのが、全身検査としての「典型的な長さ１００ｃｍのＦＤＧ（フルオロデオキシグルコース：fluoro-deoxyglucose）の研究」での現状である。これだけの計数を蓄積するのに要する時間は、注入量及びスキャナの感度とスキャナの計数能力によって決まる。
【０００８】
ＰＥＴイメージングシステムは、被検体から飛散するガンマ線を検出するために、相互に対向して配設された検出器を使用する。典型的には、各角度から飛来するガンマ線を検出するために、環状に配置された検出器を使用する。したがって、ＰＥＴスキャナは、本質的に等方性になっているはずの照射をできるだけ多く捕獲できるように、典型的には実質的に円筒形になっている。検出器として部分円環を使用し、抜けている角度のものを捕獲するのに検出器を回転することも考えられないことではない。しかし、これらの手法を実施しても、結果は、スキャナ全体の感度に対して厳しいものとなる。１つの面に含まれる全てのガンマ線が検出器と相互作用を起こす機会を有する円筒形状では、軸方向の寸法が大きくなると、感度又は照射を捕獲する能力に非常に有利な効果が現れる。したがって、最良の構造は、全てのガンマ線が検出される可能性がある球体構造である。当然ながら、人体への適用では、球状構造は、極めて大きくなり、極めて高価にならざるをえない。したがって、現実的には、検出器の軸方向の長さが可変である円筒形状が、最新のＰＥＴスキャナの構造の基本である。
【０００９】
ひとたびＰＥＴスキャナ全体の形状が分かれば、もう一つの問題は、できるだけ多くのシンチレーション料をガンマ線経路に配置して、できるだけ多くのガンマ線を停止させて光に変換できるようにすることである。断層再構成の原理によって放射性同位体の時空間分布を再構成できるようにするために、検出された各事象のエネルギー（即ち発生した光の量）、部位、及びタイミングの特徴を決定する必要があるであろう。ほとんどの最新ＰＥＴスキャナは、数千の個別結晶で構成されている。それらの結晶はモジュールの形で配置されており、シンチレーション事象の位置を識別するのに用いられている。典型的には、結晶エレメントの断面は、概ね４ｍｍ×４ｍｍである。それよりも小さい又は大きい寸法で、正方形以外の断面も可能である。結晶の長さ又は深さがガンマ線を捕獲する確率を決定することになるのであるが、典型的には、１０〜３０ｍｍの範囲である。検出器モジュールが、スキャナの主たる構成部品である。
【００１０】
ＰＥＴイメージングは、高速且つ高輝度なシンチレーション結晶によるガンマ線から光への変換に依存している。シンチレータの中の相互作用位置を決定し、個別事象の時間ペアリングを行った後で、消滅プロセスの部位を再生することができる。これらの行為を行うには非常に高速のコンポーネント（検出器及び電子機器）が必要であり、卓越した信号対ノイズ比も必要である。高品質の電子機器を用いれば、信号対ノイズ比は主に、検出プロセスに関与する固有のポアソン統計によって決まる。より多くのフォトンを検出すれば信号対ノイズ比が向上し、したがって空間分解能及びタイミング分解能が一層高まることになる。検出プロセスでの重大な光損失は、検出器の設計及び電子機器の改善をもっても補償することはできない。捕集した光全体の割合（シンチレータで発生した量に対して）は、設計の効率を表すのにふさわしい指標である。光の捕集量を最大化しようとする人は誰でも、光センサをできるだけシンチレーション結晶に接近させるとともに、反射及び他のエッジ効果を回避しようとするであろう。このようなことを行えば、出来上がるものは、否応なく、結晶とセンサとの間の距離が短い検出器の大きなアレイになってしまうであろう。
【００１１】
上記のように、ＰＥＴイメージングシステムは単なる計数器ではなく、シンチレーション事象の存在を検出することに加え、その検出部位の識別も行う必要がある。各相互作用の部位を識別できるようにするための最も直接的な設計は、概念的には、恐らくシンチレータ結晶毎に独立したフォトセンサ及びデータ収集チャネルを有するようにすることであろう。共通のフォトセンサの物理的な大きさ、各データ収集チャネルに必要な電力、及びこれらの品目の関連コスト、といった制約のために、フォトセンサ数及び電子機器のチャネル数を削減する目的で、通常は何らかの多重化が採用される。２つの最も一般的な多重化の形態は、光多重（光分配）、又は、アナログ式の電子的多重（抵抗性電荷分配網）である。
【００１２】
ありとあらゆる事象を空間的及び時間的に適切に把握するために、かなりの取り組みがなされてきた。ある一連の追加的補正は、イメージングシステムが有する若干の非理想的状態を補償する。例えば、感度補正は個々の結晶の微小な違いに対処し、利得補正は光電子増倍管（ＰＭＴ：Photomultiplier  Tube）の若干の固有利得差を補償し、複合システム行列（complex  system  matrix）は環状検出器の結晶配置に存在する小さな間隙を補うことができる、といった具合である。しかし、これらの非理想的状態がイメージングシステムに及ぼす影響はどれも、起こり得る患者の動きの作用に比べれば、あまり重要ではない。
【００１３】
患者の動きは、指示及び固定ベルトによって、ある程度は統御することができる。しかし、被検体を載せる天板である被検体天板（patient  pallet）のたわみは、スキャナの設計に内在するものである。被検体天板の機械的特性（例えば、たるみ、たわみ、曲がり等）は、ひどすぎる場合は画像品質に影響しかねない。例えば、これらの機械的特性は、連続的で累積的な各ＰＥＴスキャンの間、又はＰＥＴ／ＣＴイメージングシステムのＣＴスキャナによるスキャンとＰＥＴスキャナによるスキャンとの間の垂直方向の位置ずれを招くおそれがある。被検体天板は、典型的にはたわみを最小限に抑えるように図られているものの、ほとんどのイメージングシステムで数ｍｍ（ミリメートル）の変位が起こり得る。この変位によって、スキャナの画像品質が最適なものから外れてしまうことは、考えられないことではない。
【００１４】
したがって、たわみの影響に対して採り得る１つの対処手法は、被検体天板の設計段階であらゆるたわみを補償することである。しかし、被検体天板のたわみを最小限に抑えるのに十分なだけの剛性を有するシステムを構築すれば、被検体天板のコストが上昇するとともに複雑性が増すことになるであろう。加えて、必要となる新たな機械部品を組み込むことができないおそれがあるため、既存の機械は、再設計された被検体天板を再装備するのに必ずしも適していないであろう。
【００１５】
画像を利用した測定及び補償は、たわみの影響に対するもう一つの対処手法である。例えば、ＣＴイメージングでは、被検体天板のたわみ量を推定し、したがって、必要な補正情報を全てイメージングシステムに提供するのに、画像そのものを使用することができる。但し、ＰＥＴイメージングシステムでは空間分解能が遥かに低く、せいぜい４〜５ｍｍである。更に、空間分解能は、ＰＥＴスキャナの中心から離れると低下する。被検体天板の構造を認識できる所では、空間分解能は１０ｍｍもの大きさになるであろう。軸から外れて放出された放射線から生じるこのような漸進的劣化は、例えば相互作用の視差及び深さに起因する。例えば、ＰＥＴスキャナの中心から放出された放射線は、結晶の４×４ｍｍの表面と相互に作用し合うことになるであろう。しかし、中心から外れている場合に、同じ放射線はより大きな結晶（例えば長さ１２ｍｍの結晶の斜め部分）を「確認」することになるであろう。かかる空間分解能は、イメージングがなされているＰＥＴ領域のたわみをＰＥＴスキャン画像誘導式に測定することに使用することを極めて困難なものにするであろう。加えて、被検体天板の表面に放射線源を付着させない限り、ＰＥＴイメージングシステムは被検体天板を認識できないであろう。
【００１６】
ＰＥＴ／ＣＴイメージングシステムにおいて、ＣＴイメージング空間でたわみを測定し、そのたわみをＰＥＴイメージング空間で流用することは、理論的に可能であろう。このＰＥＴイメージング空間は、典型的には２０〜１００ｃｍである。上記のように流用するには、諸々の変数が必要になるであろう。例えば患者の体重及び台上の患者の体重分布であり、この体重分布は、大きく変動するものである。しかしながら、たわみ値をＣＴフィールドからＰＥＴに適切に流用するのに必要な測定結果を得ることは、ＰＥＴの有効視野（ＦＯＶ）で、たわみそのものを測定するのと同程度に困難であると思われる。
【００１７】
図１は、ＰＥＴイメージングシステムの被検体天板の曲がりが及ぼす影響例を説明するための図である。図１は、被検体天板のたわみがＰＥＴのサジタル画像に及ぼす影響例を示している。同図の曲がり（即ち被検体天板のたわみ）は、たわみの影響がより明白になるように誇張されている。図１に例示しているように、被検体天板のたわみの度合いは、イメージングすべき領域によってまちまちである。すなわち、ＰＥＴイメージングシステムは、天板（テーブル）を寝台から移動させて、スキャン領域である「ベッド１」、「ベッド２」及び「ベッド３」において、例えばＰＥＴのステップアンドシュートスキャン（step  and  shoot  scan）を行なう。かかる場合、「ベッド１」、「ベッド２」及び「ベッド３」のＰＥＴサジタル（sagittal）画像を組み合わせた組み合わせ画像において、例えば、「ベッド１」では、天板の移動距離が長いために、天板は、深い曲がり角で沈下する。また、例えば、組み合わせ画像において、「ベッド３」では、天板の移動距離が短いために、天板は、浅い曲がり角でわずかに沈下する。すなわち、「ベッド１」、「ベッド２」及び「ベッド３」のＰＥＴ画像それぞれで対応する画素は、実際の被検体の同一部位を示す画素とは異なる。
【先行技術文献】
【特許文献】
【００１８】
【特許文献１】米国特許第７６９７７３８号明細書
【発明の概要】
【発明が解決しようとする課題】
【００１９】
本発明が解決しようとする課題は、天板のたわみを補償することができるポジトロン放射断層撮影（ＰＥＴ）イメージングシステム、方法及び核医学イメージングシステムを提供することである。
【課題を解決するための手段】
【００２０】
実施形態のポジトロン放射断層撮影（ＰＥＴ）イメージングシステムは、測定サブシステムと、収集サブシステムと、再構成サブシステムとを備える。測定サブシステムは、ポジトロン放射断層撮影（ＰＥＴ）において被検体を載せる天板である被検体天板が連続して移動する際の前記被検体天板のたわみを検出し、当該検出されたたわみに基づくたわみ情報を提供する。収集サブシステムは、前記測定サブシステムから前記たわみ情報を受信し、ポジトロン放射断層撮影（ＰＥＴ）スキャナから複数の同時発生事象に対応するＰＥＴ測定データを受信し、当該受信したたわみ情報及びＰＥＴ測定データを再構成サブシステムに伝達する。再構成サブシステムは、前記受信したたわみ情報を用いて、前記検出されたたわみに応じた変換行列を生成し、ＰＥＴスキャン画像を再構成するのに使用される再構成空間に前記変換行列の転置を適用して補正済み再構成空間を生成し、前記補正済み再構成空間に基づいて、前記ＰＥＴ測定データからＰＥＴスキャン画像を再構成するプロセッサを有する。
【図面の簡単な説明】
【００２１】
【図１】図１は、ＰＥＴイメージングシステムの被検体天板の曲がりが及ぼす影響例を説明するための図である。
【図２Ａ】図２Ａは、本実施形態のＰＥＴイメージングシステムの一例を示す図である。
【図２Ｂ】図２Ｂは、本実施形態のＰＥＴイメージングシステムの別の一例を示す図である。
【図３Ａ】図３Ａは、本実施形態のＰＥＴイメージングシステムの主要部を成すコンポーネント群のブロック図である。
【図３Ｂ】図３Ｂは、本実施形態のＰＥＴイメージングシステムに含まれるサブシステムのより詳細な模式図である。
【図４】図４は、本実施形態に係るＰＥＴイメージングシステムの計算ユニットの一例におけるブロック図である。
【図５】図５は、本実施形態で実行されるＰＥＴイメージング方法の各ステップを示すフローチャートである。
【図６Ａ】図６Ａは、被検体天板のたわみを補償する方法の一例を示すフローチャートである。
【図６Ｂ】図６Ｂは、被検体天板のたわみを補償する図６Ａの一例とは異なる方法の一例を示すフローチャートである。
【図７Ａ】図７Ａは、ブラッグ干渉回折格子を有する光ファイバ及び２芯前方時分割多重（ＦＴＤＭ）光ファイバをそれぞれ湾曲させたときに出力光信号に及ぶ影響例を示す図（１）である。
【図７Ｂ】図７Ｂは、ブラッグ干渉回折格子を有する光ファイバ及び２芯前方時分割多重（ＦＴＤＭ）光ファイバをそれぞれ湾曲させたときに出力光信号に及ぶ影響例を示す図（２）である。
【発明を実施するための形態】
【００２２】
以下、添付図面を参照して、ポジトロン放射断層撮影（ＰＥＴ）イメージングシステムの実施形態を詳細に説明する。
【００２３】
（実施形態）
本実施形態は、まるで完全に剛直なシステムで収集されたかのようなＰＥＴ画像を形成する目的で、ＰＥＴイメージングシステムの被検体天板のたわみを補償するためのシステム及び方法に関する。被検体天板のたわみの影響を補償することによって、患者の体重及び被検体天板上の体重分布に依存しない、高画質のＰＥＴスキャン画像を提供することができる。更に、たわみを補償することによって、ＰＥＴスキャン画像を、対応するＣＴ画像と位置合わせすること、並びに絶対基準点となって参照されるＰＥＴスキャン画像を形成することが、容易に行えるようになる。
【００２４】
したがって、以下で説明する実施形態には、被検体天板のたわみを患者毎に決定できるシステム及びそのたわみ情報を使用してＰＥＴ画像の一連の再構成処理におけるたわみを補償する方法の各説明が含まれている。
【００２５】
図２Ａは、本実施形態のＰＥＴイメージングシステムの一例を示す図である。図２Ａは、ＰＥＴイメージングシステム１０の一実施形態を示している。ＰＥＴイメージングシステム１０は、ＰＥＴスキャナ１４及びＣＴスキャナ１６（任意選択）が搭載された構造体１２を備える。ＣＴスキャナ１６は任意選択であり、ＰＥＴイメージングシステム１０の他の実施形態では除外している。ＰＥＴイメージングシステム１０は、更に、被検体天板１８と、被検体天板１８のたわみを検出するのに用いる光ファイバ２０と、を備える。光ファイバ２０を使用することによって、長手方向の測定システムを設けて天板のたわみを測定している。一実施形態では、光ファイバ２０を、被検体天板１８の側面に固定している。他の実施形態では、光ファイバ２０を、被検体天板１８の他の表面に取り付けたり、又は、被検体天板１８の中に組み込んだりする。光ファイバ２０は、被検体表面及び被検体天板機構（例えばローラー、支持体等）との干渉が回避されている限りは、被検体天板１８の上、又は、被検体天板１８の内部のどこに配置しても構わない。
【００２６】
更に、被検体天板１８のたわみを検出するのに、１本又は複数本の光ファイバ２０を使用することができる。或いは、光ファイバ２０に代えて、又は、光ファイバ２０に追加して、レーザ、ストリング若しくはバーコード付きエンコーダ、といった他のセンサのうちの１つ又はそれらの組み合わせを使用することができる。但し、他のセンサを使用するには、典型的には被検体天板の動きに垂直な平面で測定を行うことが必要になる。図２Ｂは、１つの例を示している。同図において、被検体天板１８の下に配設されたレーザ２２が、光ファイバ２０に追加されている。この例では、ＰＥＴスキャナ１４とＣＴスキャナ１６との間にレーザ２２が設置されている。但し、レーザの場所はこの場所に限定されず、構造体１２の外部に配置することも可能である。更に、他の実施形態では、レーザ２２は、被検体天板１８の側部（例えば図２Ｂの紙面に垂直な方向に）又はその上といった他の場所に配置されている。なお、図２Ａ及び図２Ｂでは、被検体天板１８が寝台装置から送り出されることで、下方向（図中の方向４２５を参照）に沈み込んでいる様相を示している。また、図２Ａ及び図２Ｂでは、被検体天板１８がたわみ無しで水平に移動した場合の位置を破線で示している。なお、図２Ａ及び図２Ｂは、被検体天板１８の側面に固定された光ファイバ２０を水平方向にて図示しているが、実際には、被検体天板１８に固定された光ファイバ２０は、被検体天板１８がたわむことにより、被検体天板１８と同様に曲がる。
【００２７】
図３Ａは、本実施形態のＰＥＴイメージングシステム１０の主要部を成すコンポーネント群のブロック図である。ＰＥＴイメージングシステム１０は、図３Ａに示すように、測定サブシステム５０と、収集サブシステム５２と、再構成サブシステム５４と、被検体天板１８と、ＰＥＴスキャナ１４と、ＣＴスキャナ１６（任意選択）と、を備える。測定サブシステム５０は、被検体を載せる天板である被検体天板１８のたわみを検出し、当該検出されたたわみに基づくたわみ情報を提供する。収集サブシステム５２は、測定サブシステム５０からたわみ情報を受信し、ＰＥＴスキャナ１４から複数の同時発生事象に対応するＰＥＴ測定データを受信し、当該受信したたわみ情報及びＰＥＴ測定データを再構成サブシステム５４に伝達する。再構成サブシステム５４は、受信したたわみ情報及びＰＥＴ測定データを用いて、ＰＥＴスキャン画像を再構成するプロセッサを有する。
【００２８】
図３Ｂは、本実施形態のＰＥＴイメージングシステム１０に含まれるサブシステムのより詳細な模式図である。図３Ｂに例示しているように、収集サブシステム５２は、ＰＥＴスキャナ１４、ＣＴスキャナ１６（任意選択）、及び測定サブシステム５０からデータを収集し、そのデータを再構成サブシステム５４に提供する。収集サブシステム５２は、１つ又は複数の有線式及び／若しくは無線式の通信方法によって、データの収集及び提供を行う。
【００２９】
測定サブシステム５０は、被検体天板１８のたわみを正確にリアルタイムに検出する。本実施形態では、１本又は複数本の光ファイバ２０と、上記１本又は複数本の光ファイバ２０が出力する光照射の少なくとも１つの特性を決定する少なくとも１つの測定ユニット３３２と、が測定サブシステム５０を形成している。上記少なくとも１つの測定ユニット３３２は、被検体天板１８のたわみに比例する上記少なくとも１つの特性（例えば光照射強度）を決定し、上記少なくとも１つの特性に関するたわみ情報を、収集サブシステム５２に含まれるデータ収集ユニット３２４に提供する。比例の仕方は、例えば単純なもの、線形のもの、又はモデルによるもっと複雑なものである。
【００３０】
たわみ情報は、収集サブシステム５２によって、（ステップアンドシュート型寝台動作に対しては）ばらばらの時間に、又は、（連続型寝台動作に対しては）一定時間間隔で取り込まれる。本実施形態では、収集サブシステム５２は、測定サブシステム５０からたわみ情報を受信すると、受信したたわみ情報を、例えば天板の公称高さ、天板のイメージング開口の中への伸び、他の生理学的信号（例えば心電図）等に関する諸環境情報のうちの１つ又はそれらの組み合わせと共に、再構成サブシステム５４へと提供する。或いは、上記１つ又は組み合わせの情報は、独立して再構成サブシステム５４に提供してもよい。再構成サブシステム５４は、取り込んだたわみ情報を用いて、被検体天板１８のたわみ量の全体又はその一部分を推定する。
【００３１】
図３Ｂに示す収集サブシステム５２は、更に、ＰＥＴスキャナ１４及びＣＴスキャナ１６からそれぞれデータを収集するデータ収集ユニット３２０及び３２２も備えているが、それらのデータ収集ユニットは、同一の収集サブシステム内に含まれている必要はない。例えば、他の実施形態において、データ収集ユニット３２４は別のサブシステムの一部分であってもよく、測定サブシステム５０若しくは再構成サブシステム５４の中に組み込まれていてもよい。図３Ｂに示す一例では、データ収集ユニット３２４は、測定サブシステム５０からたわみ情報を受信する。また、図３Ｂに示す一例では、データ収集ユニット３２０は、ＰＥＴスキャナ１４から複数の同時発生事象（すなわち、消滅事象により、略同時に生じるシンチレーション事象）に対応するＰＥＴ測定データを受信する。そして、収集サブシステム５２は、受信したたわみ情報及びＰＥＴ測定データを再構成サブシステム５４に伝達する。
【００３２】
一実施形態において、既存スキャナの改造に適応するために、天板のたわみを測定する測定サブシステム５０は、被検体天板１８だけに限定されている。実際の撮像開口内の付加的な機構は、実装するのが更に複雑であり、被検体天板１８の動きそのもの又はイメージング（もしＦＯＶ内に配置される場合に）と干渉する可能性もある。
【００３３】
更に、上記では、測定サブシステム５０が提供するたわみ情報を、少なくとも１つの決定された光照射特性に関連するものとして論じているが、別の実施形態では、１本又は複数本の光ファイバ２０からの光照射そのものを、たわみ情報として収集サブシステム５２に提供する。この場合、上記少なくとも１つの測定ユニット３３２は、測定サブシステム５０から除外される。更に他の実施形態では、図３Ｂに例示しているように、測定サブシステム５０は、上記１本又は複数本の光ファイバ２０と、上記少なくとも１つの測定ユニット３３２と、被検体天板１８の実際のたわみの全体又はその一部分を、上記少なくとも１つの決定された光照射特性に基づいて推定する計算ユニット３３０と、を備える。一実施形態では、計算ユニット３３０は、たわみテーブル３３４を用いて実際のたわみを推定する。このテーブルは、一実施形態では、校正プロセスの間に発生される。これについては後で論じる。上記更に他の実施形態において、測定サブシステム５０は、その推定したたわみを、たわみ情報として収集サブシステム５２に提供する。
【００３４】
開示の実施形態の測定サブシステム５０は更に、前に論じた他のセンサのうちの１つ又はそれらの組み合わせも具備することができる。或いは、上記１本又は複数本の光ファイバ２０を、上記他のセンサのうちの１つ又はそれらの組み合わせで置き換えてしまってもよい。加えて、前に言及したように、測定サブシステム５０は、有線式及び／又は無線式の通信方法によってたわみ情報（例えば、上記少なくとも１つの決定された光放射の特性）を提供することができる。更に、他の実施形態では、たわみ情報を再構成サブシステム５４に、収集サブシステム５２を介さず直接提供する。
【００３５】
測定サブシステム５０を実現するのに使用可能な１つのシステムは、ブラッグ干渉の原理を用いて被検体天板１８のたわみを検出する。一実施形態において、測定サブシステム５０は、被検体天板１８沿いに固定された１本又は複数本の光ファイバ２０を具備している。上記１本又は複数本の光ファイバ２０は、被検体天板１８そのものの上に、ＣＴイメージング又はＰＥＴイメージングに無視できる程度の効果しか生じない状態で配置することができる。実施形態の一例において、被検体天板１８にたわみが生じると、上記１本又は複数本の光ファイバ２０は、たわみ量に比例する少なくとも１つの特性を有する光照射を提供する。上記光照射の上記少なくとも１つの特性は、被検体天板１８の所定の場所におけるたわみ量を決定するのに用いられる。例えば、被検体天板１８の１箇所又は複数箇所の所定の場所におけるたわみ量は、少なくとも１つのブラッグ干渉回折格子（Bragg  interference  grating）を有する光ファイバに送入された光照射と、その光照射が反射されて同一光ファイバを通って戻って来た光照射との間の差を監視することによって決定される。別の実施形態では、被検体天板１８の所定の場所におけるたわみ量は、２芯（dual  core）のＦＴＤＭ（前方時分割多重：Forward  time  division  multiplexing）光ファイバの内側及び／又は外側のコアから出力される光照射の強度の変化と、内側のコア及び外側のコアから出力される光放射の到達時間差と、に基づいて決定される。
【００３６】
例えばブラッグ干渉を用いた測定サブシステム５０は、たわみに比例する信号を１本又は複数本の光ファイバ２０から生成することになるが、既知のたわみ（又は、例えばレーザ２２を使って、別のシステムで測定したたわみ）がブラッグシステムとたわみとの間の対応テーブル（たわみテーブル３３４）を作製するような校正段階が必要になるものと考えられる。次いで、計算ユニット３３０は、たわみテーブル３３４を測定サブシステム５０に保存し、収集サブシステム５２に報告するときに、たわみテーブル３３４を用いて上記の信号を実際のたわみ（たわみ情報）に変換する。たわみテーブルは、自動的に発生させても良いし、手動で発生させても良く、例えば製造業者が予め設定しておいてもよい。更に、収集サブシステム５２が被検体天板の高さ又は長手方向の位置の変化を報告するのと同時に、測定サブシステム５０からのたわみ情報が、収集されたＰＥＴ測定データ及び／又はＣＴ測定データに加えられるような形に、収集サブシステム５２を設計しておく必要がある。
【００３７】
このように、測定サブシステム５０は、被検体天板１８の所定の位置において検出した少なくとも１つのたわみ量に基づいて、被検体天板１８のたわみを推定する。一例を挙げると、測定サブシステム５０は、複数のブラッグ反射回折格子を含む１本の光ファイバから反射光を受信し、光ファイバから受信した反射波の波長に基づいて、被検体天板１８の縦方向（長手方向）の複数位置におけるたわみ量を決定する。或いは、他の一例では、測定サブシステム５０は、１つのブラッグ反射回折格子をそれぞれが含む複数本の光ファイバから反射光を受信し、複数本の光ファイバから受信した反射波の波長に基づいて、被検体天板１８の縦方向（長手方向）の複数位置におけるたわみ量を決定する。ここで、波長等の反射波の情報に基づいて実際のたわみ（たわみ情報）を推定するために、測定サブシステム５０は、１本の光ファイバ又は複数本の光ファイバから光を受信し、１本の光ファイバ又は複数本の光ファイバそれぞれから受信した光の、被検体天板１８における異なるたわみ量に関する少なくとも１つの特性を測定する。そして、測定サブシステム５０は、１本の光ファイバ又は複数本の光ファイバそれぞれから受信した光の少なくとも１つの特性と被検体天板１８の異なるたわみ量との間の対応関係を規定するテーブル（たわみテーブル３３４）を生成する。
【００３８】
図２Ａは、被検体天板１８に取り付けられた光ファイバ２０をどのように使用すれば、ブラッグ干渉の原理によって、さえぎらないように被検体開口全体を残した状態で、被検体天板１８のたわみを評価できるのかを図解している。図７Ａ及び図７Ｂは、ブラッグ干渉回折格子を有する光ファイバ及び２芯前方時分割多重（ＦＴＤＭ）光ファイバをそれぞれ湾曲させたときに出力光信号に及ぶ影響例を示す図である。図７Ａ及び７Ｂでは、光ファイバが曲がったときに、３つのブラッグ干渉回折格子を有する光ファイバと２芯ＦＴＤＭ光ファイバのそれぞれの出力光に及ぶ影響の実例を示している。
【００３９】
例えば、図７Ａでは、光ファイバが曲がったときに、入力光の波長の反射強度（実線）と出力光の波長の反射強度（破線）との関係を示している。図７Ａに示す一例では、ある曲がりでは、波長「λＢ」の入力光に対する出力光の反射強度分布が、波長「λＢ」の入力光の反射強度分布から変化しないことを示している。また、図７Ａに示す一例では、ある曲がりでは、波長「λＡ」の入力光に対する出力光の反射強度分布のピークが、波長「λＡ」の入力光の反射強度分布のピークより「ΔλＡ」増加していることを示している。また、図７Ａに示す一例では、ある曲がりでは、波長「λＣ」の入力光に対する出力光の反射強度分布のピークが、波長「λＣ」の入力光の反射強度分布のピークより「ΔλＣ」減少していることを示している。また、例えば、図７Ｂでは、光ファイバが曲がったときに生じる出力光の波長から変換された電圧値の時間変化曲線の一例を示している。図７Ｂの一例に示すように、曲がりが有る場合と、曲がりが無い場合とで、電圧値の時間変化曲線は、変化する。
【００４０】
図３Ｂに例示しているように、再構成サブシステム５４は、電子記憶装置３０４と、インタフェース３０６と、ディスプレイ３０８と、計算ユニット３０２とを備える。電子記憶装置３０４は、ＰＥＴスキャナ１４、ＣＴスキャナ１６、測定サブシステム５０から取り込んだデータ又は計算ユニット３０２が再構成したＰＥＴスキャン画像のうちの１つ又はそれらの組み合わせを保存する。これについては、後で更に考察する。更に、再構成サブシステム５４が被検体天板１８の実際のたわみを推定するときは、たわみテーブル３３４が電子記憶装置３０４に保存されている。インタフェース３０６は、計算ユニット３０２の設定及び／若しくは制御を行うために、且つ／又は、計算ユニット３０２に更に別の命令を与えるために使用される。更に、ディスプレイ３０８は、ＰＥＴイメージングシステム１０を操作するために、ユーザによって使用される。図４は、本実施形態に係るＰＥＴイメージングシステムの計算ユニット３０２の一例におけるブロック図である。計算ユニット３０２は、主記憶装置４４０、及び／又は、ＲＯＭ４５０に保存されているデータ及び命令を処理するプロセッサ４８０を具備している。プロセッサ４８０は、更に、ディスク４１０又はＣＤ−ＲＯＭ４２０に保存されている情報も処理することができる。本例のプロセッサ４８０は、米国インテル社によるＸｅｏｎプロセッサ（登録商標）又は米国ＡＭＤ社によるＯｐｔｅｒｏｎプロセッサ（登録商標）とすることができる。当業者であれば気付くように、プロセッサ４８０は、Ｐｅｎｔｉｕｍ（登録商標）プロセッサ、Ｃｏｒｅ２Ｄｕｏプロセッサ（登録商標）、及びこれらに類するものとすることも可能である。したがって、ガンマ線検出方法に対応する命令は、ディスク４１０、ＣＤ−ＲＯＭ４２０、主記憶装置４４０、又はＲＯＭ４５０のうちの任意のものに保存することができる。
【００４１】
計算ユニット３０２は更に、インターネット又は私的ネットワークといったネットワークとの接続を行うための、米国インテル社によるＩｎｔｅｌＥｔｈｅｒｎｅｔ（登録商標）Ｐｒｏネットワークインタフェースカード等のネットワークインタフェース４７５も具備することができる。ディスプレイ制御部４３０は、ディスプレイ３８５と接続するための、米国エヌビディア社のＮＶＩＤＩＡ  Ｇ−ＦｏｒｃｅＧＴＸグラフィックスアダプタ（登録商標）とすることができる。計算ユニット３０２は更に、キーボード２９５、ポインティングデバイス２８５、又はマイク、トラックボール、ジョイスティック、タッチスクリーン、及びこれらに類するインタフェース等の他の汎用インタフェース３０６と接続するためのＩ／Ｏインタフェース４９０も具備することができる。
【００４２】
ディスク制御部４５０は、ディスク４１０をバス４７０と相互接続する。このディスクは、ハードディスクドライブ若しくはフラッシュメモリドライブ、及びＣＤ−ＲＯＭ４２０若しくはＤＶＤドライブとすることができる。このバスは、ＩＳＡ、ＥＳＩＡ、ＶＥＳＡ、ＰＣＩ、又は同様のものとすることができ、計算ユニット３０２の全コンポーネントを相互接続する。計算ユニット３０２のコンポーネント群はその特徴が既によく知られているため、それらの全般的特徴及び機能については、簡潔にする目的で説明を省略する。もちろん、米国フリースケール社のＦｒｅｅｓｃａｌｅＣｏｌｄＦｉｒｅ、Ｉ．ＭＸ、ＡＲＭプロセッサ（登録商標）といった、当該技術分野で知られている他のプロセッサ及びハードウェアの販売業者及びタイプも、本発明に使用することができる。
【００４３】
本例の計算ユニット３０２は更に、ＦＰＧＡ、ＡＳＩＣ、マイクロコントローラ、ＰＬＤ、又は光ディスク等の他のコンピュータ可読媒体上に独立して実装することも可能である。加えて、本例の計算ユニット３０２は、ＰＣ等の計算装置のハードウェアプラットフォームであり、プロセッサ４８０は、例えばＩｎｔｅｌＰｅｎｔｉｕｍ（登録商標）Ｐｒｏｃｅｓｓｏｒ又は当該技術分野で知られている他の任意のプロセッサとすることができる。主記憶装置４４０、ＲＯＭ４５０、ディスク４１０、又はＣＤ−ＲＯＭ４２０のいずれかに保存されるコンピュータ可読命令は、プロセッサ４８０及びＭｉｃｒｏｓｏｆｔＶＩＳＴＡ（登録商標）、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ（登録商標）、ＬＩＮＵＸ（登録商標）、ＡｐｐｌｅＭＡＣ−ＯＳ（登録商標）、及び当業者に公知の他のシステム、といったオペレーティングシステムと一緒に実行する、ユーティリティプログラム、バックグラウンドデーモン、又はオペレーティングシステムの構成要素若しくはそれらの組み合わせとして提供することができる。
【００４４】
主記憶装置４４０及び／又はＲＯＭ４５０は、レジストリ機能及びこれらに類する計算ユニット３０２の機能をサポートしている。したがって、ＲＯＭ４５０はＰＲＯＭ等の読出専用メモリであるが、主記憶装置４４０は、ランダムアクセスメモリ（ＲＡＭ）、ＦＬＡＳＨメモリ、ＥＥＰＲＯＭメモリ、又はこれらに類する記憶装置とすることができる。主記憶装置４４０及びＲＯＭ４５０は既によく知られているため、簡潔にする目的で、それらの説明はこれ以上行わない。
【００４５】
計算ユニット３０２の例示的実施形態として図４の説明を行ってきたが、図４のコンポーネント群は、被検体天板１８のたわみを推定する測定サブシステム５０の計算ユニット３３０を実現するのに使用することも可能である。
【００４６】
図５は、本実施形態で実行されるＰＥＴイメージング方法の各ステップを示すフローチャートである。具体的には、図５は、被検体天板１８のたわみを補償しながらＰＥＴスキャン画像を再構成するための各ステップを示すフローチャートである。ステップ５１０において、被検体天板１８のたわみが検出される。被検体天板１８のたわみは、前に論じたように、測定サブシステム５０を用いて検出することができる。ステップ５２０において、測定サブシステム５０によって、検出されたたわみに基づいて、たわみ情報（例えば上記光照射、上記少なくとも１つの決定された照射特性、又は上記推定たわみ）が生成され、収集サブシステム５２に、たわみ情報が提供される。一実施形態では、そのたわみ情報は、被検体天板１８の所定の位置における少なくとも１つのたわみ量に対応している。ステップ５３０及びステップ５４０において、収集サブシステム５２が、たわみ情報及びＰＥＴスキャナ１４からのＰＥＴ測定データが受信する。一実施形態では、収集サブシステム５２は、たわみ情報を、再構成サブシステム５４へのデータの流れの中で、例えば寝台（天板）の公称高さ、天板のイメージング開口の中への伸び、他の生理学的信号（例えば心電図）等に関する環境情報と共に提供する。そのデータの流れは更に、ＰＥＴ測定データも含んでいる。環境情報を含んでいるため、再構成サブシステム５４は、ＰＥＴ測定データが収集される環境を表現する正確且つ「そのとき」のパラメータセットを用いてＰＥＴスキャン画像を再構成することができる。或いは、測定サブシステム５０によって、たわみ情報を再構成サブシステム５４に、直接伝達してもよい。更に、ステップ５５０において、再構成サブシステム５４が、受信したたわみ情報及びＰＥＴ測定データに基づくＰＥＴ走査画像を再構成する。
【００４７】
以下、被検体天板１８のたわみを補償するための方法の実施形態について説明する。被検体天板制御部によって監視され通知される被検体天板１８の水平方向の動き、及び被検体天板１８のうちのイメージングすべき被検体の領域に対応する部分の実際のたわみ量を提供する測定サブシステム５０によって、図２に示すＰＥＴイメージングシステム１０の完全な幾何学的形状が決まる。この関係が決まった後、再構成サブシステム５４の計算ユニット３０２によって、次の２つの手法のうちの一方を用いて、被検体天板１８のあらゆるたわみが補正される。すなわち、第１の手法では、図６Ａに例示しているように、（１）再構成画像に直線移動及び回転移動の補正要素を適用する。又は、第２の手法では、図６Ｂに例示しているように、（２）画像空間の変換を行う。これら２つの手法は、技術的には等価であるが、統計上の小さな不一致及び／又はアルゴリズム上の利点が、それらを別個のものにしている。
【００４８】
具体的には、第１の手法では、たわみが有る状態でＰＥＴスキャン画像を再構成し、その再構成画像に一定の変換を適用する。図６Ａは、被検体天板１８のたわみを補償する方法の一例を示すフローチャートである。図６Ａのステップ６１０に例示しているように、計算ユニット３０２のプロセッサ４８０は、検出した被検体天板１８のたわみに従って変換行列「Ｔ」を生成する。次いで、ステップ６２０において、プロセッサ４８０は、その変換行列を、被検体天板１８のたわみを補正しないで再構成された、すなわち、たわみを伴って再構成された仮の再構成ＰＥＴスキャン画像に適用して、補正済みの再構成ＰＥＴスキャン画像を生成する。被検体天板１８のたわみを補償する図６Ａの一例とは異なる方法の一例を示すフローチャートである。一方、第２の手法を用いると、図６Ｂのステップ６５０において、プロセッサ４８０は、検出した被検体天板１８のたわみに従って変換行列「Ｔ」を生成する。次いで、ステップ６６０において、プロセッサ４８０は、第１の手法で用いた変換行列「Ｔ」の転置を、ＰＥＴスキャン画像を再構成するのに使用される再構成空間（又はフィールド）に適用して補正済み再構成空間を生成する。したがって、再構成データは既に「正しい」空間に入っている。ステップ６７０において、プロセッサ４８０は、補正済みの再構成空間に基づいて、補正済みの再構成されたＰＥＴスキャン画像を生成する。上記の変換で使用した変換行列は、典型的には、３軸全てにおいて直線移動及び回転移動を含んでいる。或いは、無視できると分かっている変位量又は回転量は、ゼロに設定されている。
【００４９】
以下、数式を用いて、上記の手法について詳細に説明する。ＰＥＴ測定データの形成は、通常は、以下の式（１）のようにモデル化される。
【００５０】
ｇ＝Ｈｆ＋ｂ  ・・・（１）
【００５１】
式（１）において、「ｇ」は、Ｍの検出器ビン（bin）を有する列ベクトルである測定データであり、「ｆ」は、Ｎのボクセルを有する同じく列ベクトルである放射画像（emission  image）である。「Ｈ」は、「Ｍ×Ｎ」のエレメントのシステム行列であり、一定の放射ボクセルから放出されたポジトロンが特定の検出器ビンによって検出される確率である。システム行列「Ｈ」は、「あるＬＯＲ」上の「ある位置」に集積した薬品から放出されるガンマ線が、「ある投影方向」において検出される確率を示す検出確率を示す行列であり、ＰＥＴイメージングシステム１０において物理的に一意に特定される要素である。「ｂ」は、背景のゆらぎであって、Ｍのビンを有する列ベクトルであり、ランダムな散乱フォトンに起因しており、ランダム散布推定法から分かると考えられている。
【００５２】
ＯＳＥＭ（Ordered  Subset  MLEM）法等、ＰＥＴ画像の再構成に用いられる逐次近似法では、目標は、前述のモデルを用いて測定データ「ｇ」に最善に整合する推定値「ｆ」を見つけ出すことである。なお、上記の再構成画像「ｆ」は、「未補正」の被検体天板たわみ空間に対応している。ＰＥＴ画像再構成は、通常は１回に１つの場所で静止した寝台位置に対して行われ、したがって、「ｆ」は、ステップアンドシュート型寝台動作における１つの寝台位置又は連続型寝台運動に属する特定のインスタンスのいずれかに対応する。
【００５３】
「補正済み」被検体天板たわみ空間で、放射画像「ｆｃｏｒｒ」を得るには、以下の式（２）により未補正再構成画像を変換するのが１つのやり方である（第１の手法）。
【００５４】
ｆｃｏｒｒ＝Ｔｆ  ・・・（２）
【００５５】
式中、「Ｔ」はＮ×Ｎの変換行列であり、測定サブシステム５０（例えば光ファイバ測定サブシステム）から得られたアフィン変換と、ボクセルグリッド間に位置する変換座標の再サンプリングの双方を備えている。この手法では、変換は、単純で剛性な直線移動及び回転移動の行列に対応している。
【００５６】
被検体天板のたわみを補償する別のやり方（第２の手法）は、以下の式（３）に示すように、たわみをＰＥＴ測定データ形成プロセスの中でモデル化することである。
【００５７】
ｇ＝ＨＴ−1ｆｃｏｒｒ＋ｂ  ・・・（３）
【００５８】
式中、「Ｔ−1」はＮ×Ｎの変換行列であり、たわみの無い画像空間をたわみの有る画像空間（補正済み再構成空間）に変換する。ここで、以下の式（４）及び式（５）のように、「Ｔ−1」をシステム行列「Ｈ」と組み合わせて、補正済みシステム行列「Ｈｃｏｒｒ」を形成することができる。
【００５９】
Ｈｃｏｒｒ＝ＨＴ−１  ・・・（４）
【００６０】
ｇ＝Ｈｃｏｒｒｆｃｏｒｒ＋ｂ  ・・・（５）
【００６１】
すなわち、第２の手法では、式（４）で求まる補正済みシステム行列「Ｈｃｏｒｒ」を用いて、式（５）に基づく逐次近似法によりＰＥＴ画像を再構成する。したがって、再構成画像は、たわみの無い画像空間に存在することになる。この手法は、連続型寝台動作に対して特に有用である。
【００６２】
中に画像フレームが存在しているステップアンドシュート型収集法の場合にはどちらの方法も使用されており、実質的に等価である。しかし、連続型寝台運動では、入来するＰＥＴ測定データを、画像再構成の前に補正しなければならない。というのは、実際のたわみは絶えず変化しており、したがって、再構成空間で継続的に調整する必要があるからである。画像再構成後に補正することは、連続型寝台運動には適用できない。というのは、複数の時間点で、再構成画像が、異なる量のたわみを含むことになるはずだからである。
【００６３】
以上まとめると、計算ユニット３０２のプロセッサ４８０は、検出されたたわみに応じて変換行列を生成する。そして、第１の手法を行なう場合、プロセッサ４８０は、たわみを有して再構成された仮のＰＥＴスキャン画像に対して変換行列を適用して再構成されるＰＥＴスキャン画像を生成する。或いは、第２の手法を行なう場合、プロセッサ４８０は、ＰＥＴスキャン画像を再構成するのに使用される再構成空間に変換行列の転置を適用して補正済み再構成空間を生成し、補正済み再構成空間に基づいてＰＥＴ走査画像を再構成する。
【００６４】
なお、上記の実施形態で説明した方法は、ＰＥＴイメージングシステムだけでなく、天板のたわみが発生し、ＰＥＴ画像のように、再構成画像に天板が写りこまない他のイメージングシステムにおいても適用することができる。具体的には、上記の実施形態で説明した方法は、シングルフォトンエミッションＣＴ（ＳＰＥＣＴ:  Single  Photon  Emission  computed  Tomography）イメージングシステムにおいても適用することができる。すなわち、上記の実施形態で説明した方法は、核医学イメージングシステムに適用することができる。かかる場合、測定サブシステムは、被検体天板のたわみを検出し、当該検出されたたわみに基づくたわみ情報を提供する。そして、収集サブシステムは、測定サブシステムからたわみ情報を受信し、スキャナから複数の放射線測定データを受信し、当該受信したたわみ情報及び放射線測定データを再構成サブシステムに伝達する。再構成サブシステムは、受信したたわみ情報及び放射線測定データを用いて、核医学画像を再構成するプロセッサを有する。
【００６５】
以上、説明したとおり、本実施形態によれば、天板のたわみを補償することができる。
【００６６】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００６７】
１０  ＰＥＴイメージングシステム
１２  構造体
１４  ＰＥＴスキャナ
１６  ＣＴスキャナ
１８  被検体天板
２０  光ファイバ
５０  測定サブシステム
５２  収集サブシステム
５４  再構成サブシステム
【図１】
0005931380.tif 000002
【図２Ａ】
0005931380.tif 000003
【図２Ｂ】
0005931380.tif 000004
【図３Ａ】
0005931380.tif 000005
【図３Ｂ】
0005931380.tif 000006
【図４】
0005931380.tif 000007
【図５】
0005931380.tif 000008
【図６Ａ】
0005931380.tif 000009
【図６Ｂ】
0005931380.tif 000010
【図７Ａ】
0005931380.tif 000011
【図７Ｂ】
0005931380.tif 000012