JP3604470B2

JP3604470B2 - Positron CT apparatus and image reconstruction method thereof

Info

Publication number: JP3604470B2
Application number: JP25104195A
Authority: JP
Inventors: 真介森
Original assignee: Hamamatsu Photonics KK
Current assignee: Hamamatsu Photonics KK
Priority date: 1995-09-28
Filing date: 1995-09-28
Publication date: 2004-12-22
Anticipated expiration: 2015-09-28
Also published as: JPH0990043A

Description

【０００１】
【発明の属する技術分野】
本発明は、測定対象物に投入されたＲＩ線源により発生する電子・陽電子対消滅に伴って放出される光子対を検出することにより、その測定対象物内の物質分布を測定するポジトロンＣＴ装置およびその画像再構成方法に関するものである。
【０００２】
【従来の技術】
ポジトロンＣＴ装置（ＰｏｓｉｔｒｏｎＥｍｉｓｓｉｏｎＣｏｍｐｕｔｅｄ−Ｔｏｍｏｇｒａｐｈｙ；以下、ＰＥＴと呼ぶ）は、生体や疾病患の研究あるいは臨床検査等に応用され、体内に投入された陽電子放出核種（以下、ＲＩ線源と呼ぶ）の分布を画像化し、生体機能を見るための装置である。
【０００３】
ＲＩ線源は、神経伝達に関与するドーパミンや体内でのグルコース代謝に関係するＦＤＧ（^１８Ｆ−フルオロデオキシグルコース）等の生体内物質、或いは、例えば新規開発中の薬剤に、部分的に付加されて用いられる。ＰＥＴは、このような物質の生体内での分布、消費量あるいは時間的変化の様子を見ることができる。又、ＰＥＴは、脳血流量や酸素消費量などの生体の基礎代謝を測定することもできる。
【０００４】
このようなＰＥＴの検出部は、リング状に配置された多数の光子検出器（これをリングと呼ぶ）からなり、そのリング内に、ＲＩ線源を注入あるいは吸入された人体などの測定対象物が置かれる。測定対象物内のＲＩ線源から放出された陽電子は、直ちに近くの電子と結合して、それぞれ５１１ｋｅＶのエネルギを持つ１対の光子（ガンマ線）が互いに反対方向に放出される。リングによってこの１対の光子を同時計数することにより、電子・陽電子対消滅がどの直線上で発生したかを特定することができる。ＰＥＴは、このような同時計数情報を蓄積して画像再構成処理を行って、ＲＩ線源の分布画像を作成する。
【０００５】
このようなＰＥＴによる計測に要する時間は、ＲＩ線源の半減期に依存し、数十分ないし数時間に及ぶ場合があり、従来、測定対象物は、この計測時間内では位置および方位ともに固定されている必要があった。このように計測時間が長時間に及ぶ場合には、測定対象物である人体やその他の動物にとっては非常な苦痛であり、計測時間内に測定対象物が動いて位置や方位が変化することがある。計測中に測定対象物が少しでも動くと正確が測定ができず、画像再構成して得られた画像には、測定対象物の動きに起因するアーティファクトが発生し、正確な再構成画像を得ることができなかった。
【０００６】
また、仮りに測定対象物を固定することができたとしても、測定対象物は麻酔状態ではなく覚醒状態でその生理状態を測定する場合が多いので、長時間の固定に起因するストレスの為に生理状態が変化する。この変化は測定に悪影響を与え、やはり正確な測定はできない。
【０００７】
そこで、測定対象物を固定することなく、すなわち、測定対象物の位置および方位に関する動き（以下、体動と言う）を容認することによって生理状態を維持したまま測定し、そして、その測定中に体動を測定し、この体動情報を用いて同時計数情報の補正を行う技術が知られている（例えば、特開平２−２０９１３３号公報、特開平４−１２８６７９号公報）。
【０００８】
なお、リングの光子検出器の空間的な測定分解能を向上させるために、リングを回転運動やウォブリング運動させることが行われているが、逆に測定対象物が動くことによっても測定分解能が向上する。体動補正技術は、この点に関しても有効である。
【０００９】
【発明が解決しようとする課題】
上記の体動補正技術を採用することによって、測定対象物が動く場合であっても、以前に比べれば正確な再構成画像を得ることが可能となった。
【００１０】
一方、一般にリングを構成する多数の光子検出器間に感度ムラがあるので、以上のように測定対象物から発生する光子対を検出して蓄積された同時計数情報（エミッションデータ）そのままに基づいて画像再構成処理すると正確なＲＩ線源分布画像を得ることができない。そこで、較正用線源を用いて感度補正データ（ブランクデータ）を取得し、エミッションデータをブランクデータで除算することによって感度補正が行われる。
【００１１】
しかし、測定対象物が動く場合には、このように単純に感度補正を行うことは問題がある。すなわち、測定対象物にとっては同一位置で発生し同一方向に飛行する光子対であっても、測定対象物の位置と方位が異なれば、その光子対は、異なる検出感度を有する異なる光子検出器対で受光される。それにも拘らず、従来の感度補正では、エミッションデータを同一のブランクデータで感度補正するので、正確な感度補正はできない。そこで、動きのある測定対象物を正確に測定するために、体動補正技術とともに、この場合に適した感度補正技術が求められていた。
【００１２】
本発明は、上記問題点を解消する為になされたものであり、測定対象物が測定中に動く場合であっても、その体動補正を行うとともに、精度のよい感度補正を行って正確な再構成画像を得ることができるポジトロンＣＴ装置およびその画像再構成方法を提供することを目的とする。
【００１３】
【課題を解決するための手段】
請求項１に係るポジトロンＣＴ装置は、測定空間を囲んで所定軸の周囲に配列された複数の光子検出器からなるリングによって、測定空間における電子・陽電子対消滅に伴って発生する光子対を検出し、光子対のそれぞれの光子を検出した２つの光子検出器を結ぶ直線について測定空間に設定した極座標で表現した座標値に対応した番地に所定数を累積加算して蓄積された投影データに基づいて、電子・陽電子対消滅の発生頻度の空間分布を測定するポジトロンＣＴ装置であって、（１）電子・陽電子対消滅に伴う光子対を発生させる較正用線源を測定空間で所定軸を中心にしてリングに対して相対的に回転運動させる較正用線源回転手段と、（２）較正用線源回転手段に対して較正用線源の回転運動を指示し、較正用線源が回転運動している状態で、較正用線源から発生する光子対を検出して蓄積された投影データを基準感度補正データとして獲得する基準感度補正データ生成手段と、（３）測定空間に置かれた測定対象物のリングに対する相対的な位置および方位に応じた位置方位データを出力する位置方位測定手段と、（４）測定空間に較正用線源がない状態で、測定対象物から発生する光子対を検出し、位置方位データに基づいて相対的な位置および方位の変化を補償して測定時間中に蓄積された投影データを測定データとして獲得する測定手段と、（５）測定時間中に位置方位データを入力してその位置方位データが表す相対的な位置および方位を極座標上の座標値に写像し、測定時間中における測定対象物の極座標上での存在密度分布を生成する存在密度分布生成手段と、（６）基準感度補正データを存在密度分布で加重平均計算した結果に基づいて、感度補正データを生成する感度補正データ生成手段と、（７）測定データを感度補正データに基づいて感度補正する感度補正手段と、（８）感度補正手段によって感度補正された測定データに基づいて、測定対象物における電子・陽電子対消滅の発生頻度の空間分布を算出し画像再構成を行なう画像再構成手段と、を備えることを特徴とする。
【００１４】
本装置の作用は以下のとおりである。基準感度補正データ生成手段によって、較正用線源回転手段により回転させられている較正用線源から発生した光子対がリングによって検出され蓄積された投影データが基準感度補正データとして獲得される。位置方位測定手段によって、測定空間に置かれた測定対象物のリングに対する相対的な位置および方位が測定されて位置方位データが出力され、測定手段によって、測定対象物についての投影データが測定対象物の動きを補償して測定データとして獲得されるとともに、存在密度分布生成部によって、測定対象物の動き分布を表す存在密度分布が生成される。そして、感度補正データ生成手段によって測定対象物の動きを補償した感度補正データが生成され、感度補正手段によって測定データが正確に感度補正され、画像再構成手段によって正確な再構成画像を得る。
【００１５】
請求項２に係るポジトロンＣＴ装置では、存在密度分布生成手段は、測定時間中の各時刻における座標値を蓄積し、蓄積された座標値に基づいて測定時間経過後に存在密度分布を生成する。請求項３に係るポジトロンＣＴ装置では、存在密度分布生成手段は、測定時間中の各時刻において座標値に基づいて存在密度分布を逐次生成する。
【００１６】
請求項４に係るポジトロンＣＴ装置の画像再構成方法は、測定空間を囲んで所定軸の周囲に配列された複数の光子検出器からなるリングによって、測定空間における電子・陽電子対消滅に伴って発生する光子対を検出し、光子対のそれぞれの光子を検出した２つの光子検出器を結ぶ直線について測定空間に設定した極座標で表現した座標値に対応した番地に所定数を累積加算して蓄積された投影データに基づいて、電子・陽電子対消滅の発生頻度の空間分布を測定するポジトロンＣＴ装置の画像再構成方法であって、（１）電子・陽電子対消滅に伴う光子対を発生させる較正用線源を測定空間で所定軸を中心にしてリングに対して相対的に回転運動させた状態で、較正用線源から発生する光子対を検出して蓄積された投影データを基準感度補正データとして獲得する第１のステップと、（２）測定空間に較正用線源がない状態で、測定空間に置かれた測定対象物のリングに対する相対的な位置および方位を測定し、測定対象物から発生する光子対を検出して相対的な位置および方位の変化を補償して測定時間中に蓄積された投影データを測定データとして獲得するとともに、測定時間中に相対的な位置および方位を極座標上の座標値に写像し、測定時間中における測定対象物の極座標上での存在密度分布を生成する第２のステップと、（３）基準感度補正データを存在密度分布で加重平均計算した結果に基づいて、感度補正データを生成する第３のステップと、（４）測定データを感度補正データに基づいて感度補正する第４のステップと、（５）第４のステップによって感度補正された測定データに基づいて、測定対象物における電子・陽電子対消滅の発生頻度の空間分布を算出し画像再構成を行なう第５のステップと、を備えることを特徴とする。
【００１７】
本方法の作用は以下のとおりである。第１のステップにおいて、較正用線源回転手段により回転させられている較正用線源から発生した光子対がリングによって検出され蓄積された投影データが基準感度補正データとして獲得される。第２のステップにおいて、測定空間に置かれた測定対象物のリングに対する相対的な位置および方位に応じた位置方位データが獲得され、測定対象物の動き分布を表す存在密度分布が生成されるとともに、測定対象物についての投影データが測定対象物の動きを補償して測定データとして獲得される。そして、第３のステップにおいて測定対象物の動きを補償した感度補正データが生成され、第４のステップにおいて測定データが正確に感度補正され、第５のステップにおいて正確な再構成画像を得る。
【００１８】
請求項５に係るポジトロンＣＴ装置の画像再構成方法では、第２のステップは、測定時間中の各時刻における座標値を蓄積し、蓄積された座標値に基づいて測定時間経過後に存在密度分布を生成する。請求項６に係るポジトロンＣＴ装置の画像再構成方法では、第２のステップは、測定時間中の各時刻において座標値に基づいて存在密度分布を逐次生成する。
【００１９】
【発明の実施の形態】
以下、添付図面を参照して本発明の実施の形態を詳細に説明する。尚、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
【００２０】
実施形態の要部の説明に先立って、実施形態が適用されるＰＥＴの構成、動作、画像再構成方法および感度補正方法等について説明する。
【００２１】
先ず、２Ｄ−ＰＥＴの構成と動作について説明する。図８は、２Ｄ−ＰＥＴのシステム構成図である。
【００２２】
測定空間内にある測定対象物を囲んで中心軸の周囲にｎ個の光子検出器Ｄ_ｋ（ｋ＝１，２，３，．．．，ｎ）がリング状に配置されたリングは、光子検出器Ｄ_ｋ（ｋ＝１，２，３，．．．，ｎ）それぞれと同時計数回路３０とが信号線で接続されている。リング２０の何れかの光子検出器が光子を検出すると、その光子エネルギに応じた信号が信号線を通じて同時計数回路３０に送られ、同時計数回路３０は、リング２０のそれぞれの光子検出器から到達した信号に基づいて、リング２０の内の２つの光子検出器がＲＩ線源１０から放出された所定の光子エネルギ（５１１ｋｅＶ）を有する光子対を同時検出したことを認識して、その時のこれら２つの光子検出器それぞれを示す検出器識別信号（Ｉ，Ｊ）を出力する。
【００２３】
これら検出器識別信号（Ｉ，Ｊ）は、ｔ−θ変換部４０に入力され、光子対を検出した２つの光子検出器を結ぶ直線Ｌ１をリング２０内の測定空間における極座標で表現する場合に用いられる２つの変数（ｔ，θ）それぞれを座標軸とするｔ−θ平面上の写像位置（Ｔ，θ’）に変換される。ここで、Ｔは、この直線と座標原点との間の距離であり、θ’は、この直線に座標原点から下ろした垂線の座標系主軸に対する角度である。ｔ−θメモリ５０は、この（Ｔ，θ’）に対応する番地に記憶されている投影データに１を累積加算する。このようにして、ＲＩ線源１０で発生した１対の光子対についての同時計数情報が、ｔ−θメモリ５０に蓄積される。
【００２４】
このようにしてＲＩ線源１０で発生しリング２０で検出された多数の光子対の同時計数情報は、ｔ−θメモリ５０に投影データとして蓄積される。図９は、リングにおける光子対検出と画像再構成とを説明する図である。リング２０内の測定空間に設定した極座標系のθ’方向（図９の直線Ｌ１の方向）から見たＲＩ線源１０の分布を表す画像は、この蓄積された投影データの内の、極座標系でθ＝θ’上の各Ｔ点に対応する投影データから画像再構成部６０によって再構成され、この再構成画像は画像表示部７０で表示される。
【００２５】
しかし、測定空間で発生した全ての光子対が検出されるのではない。図１０は、リングにおける光子対検出の説明図である。ＲＩ線源１０からの光子対は、ＲＩ線源１０が存在するあらゆる位置から放出され、あらゆる方向に向かって飛行するが、その内で、リング２０面上で発生しリング２０面に沿った方向に飛行した光子対のみが、リング２０で検出され得る。例えば、図９で示すように、リング２０内に存在するＲＩ線源１０から発生する光子対であっても、リング２０面に沿った方向以外の方向Ｌ２ないしＬ４に飛行した光子対は検出されない。リング２０面に沿った方向Ｌ１に飛行した光子対は、何れかの光子検出器Ｉ，Ｊで検出され、ｔ−θ変換部４０で（Ｔ，θ’）に変換され、ｔ−θメモリ５０の対応する番地の投影データに１が加算される。
【００２６】
又、製作が容易であり安価であることから、光子検出器がブロック型構成とされる。このタイプのリングでは、図１１に示すように、複数のブロックＢ_ｐそれぞれは、多数の光子検出器Ｄ_ｐｑが並列配置されていて、これらのブロックＢ_ｐがリング状に配置される（ｐ＝１，２，３，．．．、ｑ＝１，２，３，．．．）。このブロック型ＰＥＴは、単なる光子検出器の配置が異なるだけであって、動作原理、投影データの蓄積および画像再構成等は、上述の説明と差異がない。
【００２７】
上述のＰＥＴでは、ＲＩ線源から放出されてあらゆる方向に飛行する光子対の内、リング面に沿った方向に飛行する光子対のみを検出するため、ＲＩ線源から放出される光子対を捕捉する確率が小さく、検出感度が低くて統計ノイズが大きいという問題がある。検出感度を向上させるために測定対象物に多量のＲＩ線源を注入することも考えられるが、測定対象物が生体である場合には限界がある。
【００２８】
そこで、このような場合には、３次元タイプのＰＥＴ（３Ｄ−ＰＥＴ）が用いられる。図１２は、３Ｄ−ＰＥＴのシステム構成図であり、図１３は、３Ｄ−ＰＥＴのリングの構成図である。３Ｄ−ＰＥＴの多層リング２１は、前述の図９あるいは図１０に示したのと同様の光子検出器の一層リングＲ_１，Ｒ_２，．．．．，Ｒ_ｍを多層配置したものであり、ＲＩ線源１０で放出され直線Ｌ１の方向に飛行した光子対は、互いに異なる単層リングＲ_ｐ、Ｒ_ｑ（ｐ≠ｑ）それぞれに属する２つの光子検出器によっても同時計数され得る。多層リング２１の内の或２つの光子検出器が、ＲＩ線源１０から放出されたエネルギ５１１ｋｅＶの光子対を同時検出すると、同時計数回路３１は、その２つの光子検出器それぞれを示す検出器識別信号（Ｉ，Ｊ）、および、その２つの光子検出器それぞれが属する２つの単層リング間の差信号ＲＤ（ＲｉｎｇＤｉｆｆｅｒｅｎｃｅ）を出力する。これら検出器識別信号（Ｉ，Ｊ）およびリング間差信号ＲＤは、ｘ−ｙ−θ−ψ変換部４１に入力され、光子対を検出した２つの光子検出器を結ぶ直線Ｌ１を多層リング２１内の測定空間における極座標で表現する場合に用いられる４つの変数（ｘ，ｙ，θ，ψ）それぞれを座標軸とするｘ−ｙ−θ−ψ空間上の写像位置（ｘ，ｙ，θ，ψ）に変換される。ここで、θとψは、直線Ｌ１の方向を表し、ｘとｙは、直線Ｌ１に垂直な投射平面（ＰｒｏｊｅｃｔｉｏｎＰｌａｎｅ）上の直交座標系による位置を表す。ｘ−ｙ−θ−ψメモリ５１は、この（ｘ，ｙ，θ，ψ）に対応する番地に記憶されている投影データに１を累積加算する。このようにして、ＲＩ線源１０で発生した１対の光子対についての同時計数情報が、ｘ−ｙ−θ−ψメモリ５１に投影データとして蓄積され、この投影データから画像再構成部６１によって再構成され、この再構成画像は画像表示部７１で表示される。
【００２９】
以上の何れのタイプのＰＥＴにおいても、光子検出器それぞれの検出感度は一定ではなくムラが多いので、光子検出器の感度補正を以下のようにして行う。すなわち、図１４に示すように、リング内にＲＩ線源を投与された測定対象物１０ａを置いて計測し、同時計数回路３２とｔ−θ変換部４２とを介して同時計数情報を蓄積することをエミッション計測と言い、これによってｔ−θメモリ５２に蓄積されたデータをエミッションデータと言い、Ｅ（ｔ，θ）で表す。また、これとは別に、図１５に示すように、リング内に測定対象物のない状態で、リング内でリング中心軸の周りに較正用のＲＩ線源１０ｂを回転させ、模擬的な平行光を光子検出器それぞれに入射させて計測（ブランク計測）を行う。このようにしてｔ−θメモリ５２に蓄積されたデータ（ブランクデータＢ（ｔ，θ））は、光子検出器対の検出感度ばらつきを表す。そして、エミッションデータＥ（ｔ，θ）をブランクデータＢ（ｔ，θ）で割ることにより、感度補正を行う。
【００３０】
さらに、測定対象物の光子吸収を以下のようにして補正する。図１６に示すように、リング内にエミッション計測時と同じ位置にＲＩ線源が投与されていない測定対象物１０ｃを置き、ブランク計測と同様にリング内で較正用のＲＩ線源１０ｂを回転させて計測（トランスミッション計測）し、これによってｔ−θメモリ５２に蓄積されたデータ（トランスミッションデータＴ（ｔ，θ））を獲得する。このトランスミッションデータＴ（ｔ，θ）をブランクデータＢ（ｔ，θ）で割れば、同時計数ライン上の吸収係数が求まり、エミッションデータＥ（ｔ，θ）をこの吸収係数で割ることにより、吸収補正ができる。なお、トランスミッション計測では、弱いＲＩ線源が更に測定対象物による吸収を受け、得られる光子対検出が減少する為、統計精度は著しく劣化する。その対策として、Ｔ（ｔ，θ）／Ｂ（ｔ，θ）をフィルタ等でスムージングしておく。
【００３１】
このようにして感度補正および吸収補正を行った後の真の投影データＰ（ｔ，θ）は、エミッションデータＥ（ｔ，θ）、ブランクデータＢ（ｔ，θ）およびトランスミッションデータＴ（ｔ，θ）を用いて、
Ｐ＝（Ｅ／Ｂ）／＜Ｔ／Ｂ＞ … （１）
で表される。ここで、記号＜＞は、スムージングを意味する。３Ｄ−ＰＥＴの場合も同様である。
【００３２】
本発明は、動きのある測定対象物であっても、体動補正を行い、且つ、エミッションデータを精度良く感度補正し、これによって正確な再構成画像を得るものである。
【００３３】
次に、本実施形態の要部について説明する。本実施形態は、リングが静止しているブロック型３Ｄ−ＰＥＴである。図１は、本実施形態に係るＰＥＴのシステム構成図である。
【００３４】
本実施形態に係るＰＥＴは、（１）光子を検出する多数の光子検出器からなるリング１２１、測定対象物１００の位置・方位を測定し位置方位データを出力する体動測定部１２３、および、較正用線源１１０を回転させる回転手段等を備える検出部１２０と、（２）リング１２１が検出した光子が電子・陽電子対消滅に伴って発生した光子対であるか否かを判断するとともに、電子・陽電子対消滅に伴って発生した光子対を検出した光子検出器対を識別する同時計数回路１３０と、（３）同時計数回路１３０が識別した光子検出器対を結ぶ直線を極座標で表した座標値（ｘ，ｙ，θ，ψ）に変換するｘ−ｙ−θ−ψ変換部１４０と、（４）座標値（ｘ，ｙ，θ，ψ）に対応する番地に記憶されている投影データに所定値を累積加算するｘ−ｙ−θ−ψメモリ１５０と、（５）体動測定部１２３から出力された位置方位データを入力し、測定対象物１００の極座標上の存在密度分布を生成する存在密度分布生成部１９０と、（６）ｘ−ｙ−θ−ψメモリ１５０に蓄積され感度補正等された投影データに基づいて再構成画像データを生成する画像再構成装置１６０と、（７）再構成画像データに基づいて再構成画像を表示する画像表示部１７０と、（８）以上の各部を制御するとともに、ｘ−ｙ−θ−ψメモリ１５０に蓄積された投影データを読み出して所定の処理を行うコンピュータ２００と、（９）コンピュータ２００が扱うデータを記憶する記憶装置１８０と、を備える。
【００３５】
測定空間１２２をその内部に含む検出部１２０のリング１２１には多数の光子検出器がリング状に配列されており、これらの光子検出器は、測定対象物１００が置かれる測定空間１２２の方向に受光面が向けられて、測定空間１２２で発生した光子を受光する。検出部１２０は、更に体動測定部（位置方位測定手段）１２３を有しており、これによって測定空間１２２に置かれた測定対象物１００の位置および方位を測定し、位置方位データを出力する。また、検出部１２０には、ブランク計測およびトランスミッション計測に使用される較正用線源１１０およびそれを回転させる回転機構が備えられている。この回転機構は、例えば、較正用線源１１０を支持して中心軸を中心に回転可能な支持機構、その支持機構に回転を伝達するベルト、その回転を発生させるモータからなり、コンピュータ２００によって制御される。較正用線源１１０は、使用時（ブランク計測時およびトランスミッション計測時）には、リング１２１の中心軸を中心としてリング１２１面上で回転運動し、それによって較正用線源１１０から発生する光子対のうちリング１２１面上で飛行するものが光子検出器対で検出される。一方、使用されない時（エミッション計測時）には、その較正用線源１１０から発生する光子対が何れの光子検出器にも到達しない位置に待避される。あるいは、その較正用線源１１０が本装置から取り去られてもよい。
【００３６】
リング１２１の内の何れかの光子検出器が光子を受光すると、それぞれの光子検出器と同時計数回路１３０との間の信号線を通じて、その光子エネルギに応じた信号が同時計数回路１３０に送られる。リング１２１のそれぞれの光子検出器から到達した信号を受け取った同時計数回路１３０は、リング１２１の内の２つの光子検出器が電子・陽電子対消滅に伴って発生する所定のエネルギ（５１１ｋｅＶ）を有する光子対を同時検出したことを認識して、その時のこれら２つの光子検出器それぞれを示す検出器識別信号（Ｉ，Ｊ）およびリング間差信号ＲＤを出力する。
【００３７】
これら検出器識別信号（Ｉ，Ｊ）およびリング間差信号ＲＤは、ｘ−ｙ−θ−ψ変換部１４０に入力され、検出器識別信号（Ｉ，Ｊ）が示す光子対を検出した２つの光子検出器を結ぶ直線をリング１２１内の測定空間１２２に設定された極座標で表現する場合に使用する４つの変数（ｘ，ｙ，θ，ψ）それぞれを座標軸とするｘ−ｙ−θ−ψ空間上の位置に写像される。このｘ−ｙ−θ−ψ変換部１４０は、体動測定部１２３から出力された位置方位データをも入力し、測定対象物１００の体動を補償して、（ｘ，ｙ，θ，ψ）値を出力する。
【００３８】
投影データを蓄積するｘ−ｙ−θ−ψメモリ１５０は、ｘ−ｙ−θ−ψ変換部１４０から出力された（ｘ，ｙ，θ，ψ）値を入力し、この（ｘ，ｙ，θ，ψ）値に対応する番地に記憶されている投影データに、所定値（例えば ”１” ）を累積加算する。このようにして、測定対象物１００から発生した光子対についての多数の同時計数情報が、ｘ−ｙ−θ−ψメモリ１５０に蓄積されて投影データとなる。
【００３９】
体動測定部１２３から出力された位置方位データは、存在密度分布生成部１９０にも入力され、測定対象物１００の測定空間１２２における位置および方位が、ｘ−ｙ−θ−ψ変換部１４０と同様にして、ｘ−ｙ−θ−ψ空間上の位置に写像される。一定時間内における測定対象物１００の位置方位データが、このようにｘ−ｙ−θ−ψ空間上の点に次々に写像されて、その一定時間内における測定対象物１００のｘ−ｙ−θ−ψ空間上における存在密度分布が生成される。
【００４０】
測定対象物１００におけるＲＩ線源の分布密度の画像を再構成する画像再構成装置１６０は、ｘ−ｙ−θ−ψメモリ１５０に蓄積された投影データ（エミッションデータ）或いは後述する補正処理後のエミッションデータから、リング１２１内の測定空間１２２に設定した極座標系の所定のθψ方向から見た測定対象物１００における光子対発生分布を表す再構成画像データを生成し、画像表示部１７０は、この再構成画像データを入力して再構成画像を表示する。
【００４１】
以上に述べた各部は、コンピュータ２００によってバスライン２１０を介して統合されて制御等がなされる。すなわち、コンピュータ２００は、較正用線源１１０の回転と待避を制御する。また、存在密度分布生成部１９０で生成された存在密度分布関数を取得する。また、ｘ−ｙ−θ−ψメモリ１５０の各番地の全ての内容を測定に先立ってゼロクリアし、各番地に蓄積された投影データを測定終了後に取得する。また、ｘ−ｙ−θ−ψメモリ１５０から取得した投影データ（エミッションデータ、ブランクデータ、トランスミッションデータ）に基づいて、後述する補正処理（感度補正、吸収補正）を行う。また、画像再構成装置１６０にエミッションデータを送出し、画像再構成装置１６０によって生成された再構成画像データを取得する。また、再構成画像データを画像表示部１７０に送出し、再構成画像を表示させる。
【００４２】
コンピュータ２００の各処理において扱われる各データは、記憶装置１８０に記憶され、また、必要時には記憶装置１８０から読み出される。すなわち、記憶装置１８０は、エミッション計測後にｘ−ｙ−θ−ψメモリ１５０から読み出されたエミッションデータを記憶するとともに、存在密度分布生成部１９０から読み出された存在密度分布関数を記憶する。また、記憶装置１８０は、ブランク計測で得られたブランクデータや、トランスミッション計測で得られたトランスミッションデータをも記憶する。また、感度補正や吸収補正の計算に際しては、記憶装置１８０は、エミッションデータ、ブランクデータ、トランスミッションデータなどが読み出され、そして、その計算の中間データや感度補正後のエミッションデータを記憶する。また、画像再構成処理に際しては、記憶装置１８０は、感度補正後のエミッションデータが読み出され、画像再構成装置で生成された再構成画像データを記憶する。なお、記憶装置１８０は、記憶領域が複数の領域に分割されており、必要に応じて上記各データをそれぞれ異なる領域に記憶する。記憶装置１８０は、磁気ディスク、磁気テープ等であってもよいし、或いは、コンピュータ２００内にある半導体メモリからなる主記憶装置であってもよい。
【００４３】
以上に述べた各部は、それぞれ独立して動作するのではなく、コンピュータ２００の指示に従い、互いに連係をもって動作する。すなわち、コンピュータ２００は、（１）較正用線源１１０を回転させてブランクデータ（基準感度補正データ）を取得する基準感度補正データ生成モジュール２０１と、（２）エミッションデータおよび存在密度分布関数を取得する測定モジュール２０２と、（３）測定対象物１００の体動を補償した感度補正データを生成する感度補正データ生成モジュール２０３と、（４）エミッションデータを感度補正する感度補正モジュール２０４と、（５）画像再構成処理と画像表示を行わせる画像再構成モジュール２０５と、を有し、それぞれのモジュールは、対応する処理ステップにおいて、各部の動作を統合制御等する。
【００４４】
基準感度補正データ生成モジュール２０１は、ブランク計測に先立ち、較正用線源１１０の回転を指示し、ｘ−ｙ−θ−ψメモリ１５０のデータをゼロクリアする。そして、ブランク計測が終了した後、ｘ−ｙ−θ−ψメモリ１５０に蓄積された基準感度補正データを獲得し、その基準感度補正データを記憶装置１８０に記憶させる。
【００４５】
また、エミッションデータを取得する測定モジュール２０２は、エミッション計測に先立ち、較正用線源１１０の停止と待避を指示し、ｘ−ｙ−θ−ψメモリ１５０のデータをゼロクリアする。エミッション計測が終了した後、ｘ−ｙ−θ−ψメモリ１５０に蓄積されたエミッションデータ、および、存在密度分布生成部１９０で生成された存在密度分布関数を獲得し、そのエミッションデータと存在密度分布関数とを記憶装置１８０に記憶させる。
【００４６】
また、感度補正データ生成モジュール２０３は、記憶装置１８０に記憶されている基準感度補正データと存在密度分布関数とを読み出し、これらのデータに基づいて感度補正データを生成し、その感度補正データを記憶装置１８０に記憶させる。
【００４７】
また、感度補正モジュール２０４は、記憶装置１８０に記憶されているエミッションデータと感度補正データとを読み出し、これらのデータに基づいてエミッションデータを感度補正する。さらに、吸収補正を行ってもよい。そして、補正されたエミッションデータを記憶装置１８０に記憶する。
【００４８】
また、画像再構成モジュール２０５は、記憶装置１８０に記憶されている感度補正されたエミッションデータを読み出し、そのデータを画像再構成装置１６０に送り、画像再構成装置１６０で処理されて得られる再構成画像データを取得し、そのデータを画像表示部１７０に送って画像表示させる。
【００４９】
次に、本実施形態に係るＰＥＴの検出部について詳細に説明する。図２は、本実施形態に係るＰＥＴの検出部の断面図である。検出部１２０は、多数の光子検出器Ｄからなるリング１２１と、測定対象物１００の位置および方位を測定する体動測定部１２３等を備える。
【００５０】
測定空間１２２を内部に有するリング１２１は、ブロック型の単層リングが多層（この図では５層）に配置された構造であり、それぞれの層は、それぞれ複数個の光子検出器Ｄからなるブロックが複数個リング状に配列され、それぞれのブロックのそれぞれの光子検出器Ｄは、測定対象物１００が置かれる測定空間１２２の方向に受光面が向けられている。遮蔽シールド１２４は、光子対が測定空間１２２外に洩れるのを防ぐためのものであり、単層リングのそれぞれを互いに隔てるものではない。したがって、異なる単層リング間で光子対を検出することが可能である。
【００５１】
測定対象物１００の位置および方位を測定する体動測定部１２３は、例えば、それぞれ発光素子と受光素子とからなる光学式の測距センサ１２３ａないし１２３ｃ、および、これらからの出力に基づいて位置方位データを出力する体動データ処理部１２３ｄからなる。測距センサ１２３ａないし１２３ｃそれぞれは、リング１２１との相対的な位置関係が固定されており、測定対象物１００の所定箇所（この図では、測定対象物１００である人体の頭部）に設けられたマーカ１０１ａないし１０１ｃそれぞれに光ビームを照射してその反射光を受光し、これらマーカ１０１ａないし１０１ｃそれぞれとの距離を測定する。３つの測距センサ１２３ａないし１２３ｃそれぞれから得られた距離データは、体動データ処理部１２３ｄに入力され、測定対象物１００の位置と方位とが求められる。
【００５２】
また、体動測定部１２３として撮像カメラを用い、得られた画像を解析して測定対象物１００の位置と方位とを求めてもよい。あるいは、体動測定部１２３として、測定対象物１００の所定箇所に設けられた加速度センサも用いてもよく、この場合、加速度センサからの出力に基づいて、測定対象物１００の位置と方位とを求めることができる。
【００５３】
次に、体動測定部１２３による測定対象物１００の位置および方位の測定および存在密度分布生成についてより詳細に説明する。図３は、本実施形態に係るＰＥＴにおける体動測定の説明図である。
【００５４】
測定対象物１００の位置は、リング１２１内の測定空間１２２に設定されたＸ−Ｙ−Ｚ直交座標系における測定対象物１００中の所定点（例えば、測定対象物１００の中心点）の座標値（Ｘ，Ｙ，Ｚ）で表される。また、測定対象物１００の方位は、そのＸ−Ｙ−Ｚ直交座標系における測定対象物１００の向き（Θ，Ψ）で表される。ここで、Θは、Ｘ−Ｙ平面上に投影された方位のＸ軸からの角度であり、Ψは、Ｘ−Ｙ平面からの仰角である。この測定対象物１００の位置および方位を表す位置方位データ（Ｘ，Ｙ，Ｚ，Θ，Ψ）は、体動データ処理部１２３ｄにおいて、測距センサ１２３ａないし１２３ｃから出力された距離データに基づいて求められる。
【００５５】
なお、体動測定部１２３は、測定対象物１００の位置および方位について基準値（例えば、測定開始時における測定対象物１００の位置および方位）からの変位量を測定し出力してもよい。この体動測定部１２３は、エミッション計測の測定開始時刻から測定終了時刻まで、一定時間（例えば、数ｍ秒ないし数十ｍ秒）間隔で測定対象物１００の位置および方位の変位量を測定し、各時刻ｔ_ｋにおける位置方位データ（ΔＸ_ｋ，ΔＹ_ｋ，ΔＺ_ｋ，ΔΘ_ｋ，ΔΨ_ｋ）を体動測定部１２３で得る（ｋ＝０，１，．．．，ｎ，．．．）。このように測定された測定対象物１００の位置および方位を表す位置方位データは、ｘ−ｙ−θ−ψ変換部１４０に入力されるとともに、存在密度分布生成部１９０に入力される。
【００５６】
各時刻ｔ_ｋにおける位置方位データ（ΔＸ_ｋ，ΔＹ_ｋ，ΔＺ_ｋ，ΔΘ_ｋ，ΔΨ_ｋ）を入力した存在密度分布生成部１９０は、ｘ−ｙ−θ−ψ変換部１４０と同じ要領で、各位置方位データをｘ−ｙ−θ−ψ空間上の位置（Δｘ_ｋ，Δｙ_ｋ，Δθ_ｋ，Δψ_ｋ）に写像する（ｋ＝０，１，．．．，ｎ，．．．）。図４は、本実施形態に係るＰＥＴにおける位置方位データおよびｘ−ｙ−θ−ψ空間上の位置の図表である。このようにして、一定時間内におけるｘ−ｙ−θ−ψ空間上の写像位置の分布を求めることによって、測定対象物１００のｘ−ｙ−θ−ψ空間上での存在密度分布を表す存在密度分布関数Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）が生成される。
【００５７】
この存在密度分布関数Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）は、測定対象物１００の同一地点で発生して同一方向に飛行した光子対がリング１２１で検出されて同時計数情報がｘ−ｙ−θ−ψメモリ１５０に蓄積されるに際して、測定対象物１００が基準位置・基準方位（例えば、エミッション計測開始時における位置・方位）にあるときには、同時計数情報が座標値（ｘ，ｙ，θ，ψ）に対応する番地に格納されるが、測定の最中に測定対象物１００が変位（Δｘ，Δｙ，Δθ，Δψ）なる位置・方位にあるときには、その同時計数情報が座標値（ｘ＋Δｘ，ｙ＋Δｙ，θ＋Δθ，ψ＋Δψ）に対応する番地に格納され、その割合が値Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）であることを示すものである。
【００５８】
なお、この存在密度分布関数Ｃθψは、一般にθ値およびψ値によって、すなわち、投影方向によって異なる。図５は、本実施形態に係るＰＥＴにおける測定対象物の各位置・各方位における存在密度分布の説明図である。この図は、測定対象物（人体頭部）１００の頂部から見た図であり、Ｌａ，ＬｂおよびＬｃそれぞれの投影方向における存在密度分布関数Ｃθψを、それぞれ、Ｃａ，ＣｂおよびＣｃとして表したものである。なお、存在密度分布関数は実際には４変数の関数ではあるが、この図では簡略化して表している。この図に示すように、測定対象物１００の体動が、或方向または或回転方向に偏っていると、存在密度分布関数Ｃθψは、θとψによって異なる関数となる。
【００５９】
次に、体動測定部１２３により測定された測定対象物の体動を補償して感度補正データを生成し更に感度補正する方法について説明する。尚、以下に述べる感度補正データ生成および感度補正それぞれは、コンピュータ２００の感度補正データ生成モジュール２０３および感度補正モジュール２０４それぞれによって処理される。図６は、本実施形態に係るＰＥＴにおける感度補正データ生成および感度補正のデータフロー図である。
【００６０】
エミッション計測中に存在密度分布生成部１９０で生成された存在密度分布関数Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）と、基準感度補正データ生成モジュール２０１で求められた基準感度補正データＢ（ｘ，ｙ，θ，ψ）とから、エミッション計測が終了した後に、測定対象物１００の体動を補償した感度補正データＢ１（ｘ，ｙ，θ，ψ）を、
【００６１】
【数１】

【００６２】
なる関係式で求める。ここで、４重の総和演算Σは、測定対象物１００の体動の範囲に亘って計算する。この計算は、測定対象物１００からの同一の光子対が、測定値対象物１００の体動に従って異なる光子検出器対で検出されることがあっても、その光子対を検出する可能性のある光子検出器対について存在密度分布Ｃθψで表される重みを考慮して平均化した検出感度を算出することを意味している。
【００６３】
以上より、エミッション計測で得られたエミッションデータＥ（ｘ，ｙ，θ，ψ）を、測定対象物１００の体動を補償した感度補正データＢ１（ｘ，ｙ，θ，ψ）で感度補正することによって、感度補正後のエミッションデータＥ１（ｘ，ｙ，θ，ψ）が、
Ｅ１＝Ｅ／Ｂ１ … （３）
で得られる。このエミッションデータＥ１は、測定対象物１００の体動を補正してエミッション計測するのみでなく、測定対象物１００の体動を補償したブランクデータＢ１で感度補正したものであるので、正確な感度補正がなされ、これに基づいて画像再構成すれば正確な再構成画像を得ることができる。
【００６４】
もし、さらに吸収補正を行う場合には、エミッション計測と同じ要領でトランスミッション計測を行い、体動補償のなされたトランスミッションデータＴ、および、トランスミッション計測中の存在密度分布関数を獲得し、これらに基づいて感度補正データＢ２を（２）式の要領で算出する。そして、真の線和Ｐ１は、
Ｐ１＝（Ｅ／Ｂ１）／＜Ｔ／Ｂ２＞ … （４）
なる関係式で求められる。ここで、記号＜＞は、スムージングを意味する。
【００６５】
次に、本実施例に係るＰＥＴの作用および画像再構成方法について説明する。図７は、本実施形態に係るＰＥＴにおける計測のフローチャートである。
【００６６】
先ず、ステップＳ１では、コンピュータ２００の基準感度補正データ生成モジュール２０１からの指示により、ブランク計測を行う。すなわち、リング１２１内に測定対象物１００を置かないで、較正用線源１１０を中心軸を中心にしてリング１２１面上で回転させて模擬的な平行光を生成する。その状態で、その較正用線源１１０から発生する光子対は、リング１２１を構成する多数の光子検出器により検出され、同時計数回路１３０でエネルギ弁別され、ｘ−ｙ−θ−ψ変換部１４０で、光子対を検出した２つの光子検出器を結ぶ直線に対応する（ｘ，ｙ，θ，ψ）値に変換され、その（ｘ，ｙ，θ，ψ）値に対応するｘ−ｙ−θ−ψメモリ１５０の番地に投影データとして蓄積される。このようにして投影データとして蓄積された基準感度補正データＢ（ｘ，ｙ，θ，ψ）は、基準感度補正データ生成モジュール２０１により、読み出されて記憶装置１８０に記憶される。
【００６７】
続いて、ステップＳ２では、コンピュータ２００の測定モジュール２０２からの指示により、エミッション計測を行う。すなわち、較正用線源１１０を待避させ、ＲＩ線源を注入された測定対象物１００をリング１２１内の測定空間１２２に置いて、この状態で、測定対象物１００から発生した光子対が検出されて、体動測定部１２３から出力された測定対象物１００の位置および方位を示す位置方位データに基づいて体動補正されて、ｘ−ｙ−θ−ψメモリ１５０に投影データが蓄積される。同時に、位置方位データは存在密度分布生成部１９０に入力されて、存在密度分布関数Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）が生成される。このようにして投影データとして蓄積されたエミッションデータＥ（ｘ，ｙ，θ，ψ）および存在密度分布関数Ｃθψ（Δｘ，Δｙ，Δθ，Δψ）は、測定モジュール２０２により、読み出されて記憶装置１８０に記憶される。
【００６８】
もし、さらに吸収補正をも行う場合には、トランスミッション計測を行う。すなわち、ＲＩ線源が注入されていない測定対象物１００を測定空間１２２に置き、較正用線源１１０を回転させ、この状態で、エミッション計測と同じ要領で測定し、トランスミッションデータＴと存在密度分布関数とを獲得し、記憶装置１８０に記憶する。
【００６９】
なお、ブランク計測（ステップＳ１）、エミッション計測（ステップＳ２）およびトランスミッション計測は、如何なる順序で行われても構わない。
【００７０】
以上の各計測が終了すると、続いてステップＳ３で、これらの計測で獲得されたデータに基づいて、感度補正データ生成を行う。感度補正データ生成に際しては、コンピュータ２００の感度補正データ生成モジュール２０３は、記憶装置１８０から、エミッション計測時の存在密度分布関数Ｃθψおよび基準感度補正データＢ（ｘ，ｙ，θ，ψ）を読み出して、感度補正データＢ１（ｘ，ｙ，θ，ψ）を（２）式に基づいて算出する。
【００７１】
続いて、ステップＳ４で感度補正を行う。コンピュータ２００の感度補正モジュール２０４は、記憶装置１８０から読み出したエミッションデータＥ（ｘ，ｙ，θ，ψ）を、ステップＳ３で求めた感度補正データＢ１（ｘ，ｙ，θ，ψ）で割ることにより、リング１２１を構成する複数の光子検出器間の感度ムラを補正し、感度補正後のエミッションデータＥ１（ｘ，ｙ，θ，ψ）を得る。
【００７２】
もし、さらに吸収補正をも行う場合には、トランスミッションデータＴを、トランスミッション計測時の存在密度分布関数に基づいて、ステップＳ３と同じ要領で補正し、真の投影データＰ１を（４）式に基づいて算出する。
【００７３】
そして、ステップＳ４で得られた感度補正後のエミッションデータＥ１（または、真の投影データＰ１）に基づいて、ステップＳ５で画像再構成を行う。すなわち、コンピュータ２００の画像再構成モジュール２０５は、この感度補正後のエミッションデータＥ１（または、真の投影データＰ１）を画像再構成装置１６０に送出し、画像再構成装置１６０において生成された再構成画像データを獲得し、そして、その再構成画像データを画像表示部１７０に送出して再構成画像を表示させる。このようにすることにより、より正確な電子・陽電子対消滅の発生頻度の空間分布を測定することができる。
【００７４】
本発明は、上述した３Ｄ−ＰＥＴに限定されるものではなく、他のタイプのＰＥＴにも適用可能である。
【００７５】
例えば、本発明は、多層のリングからなる２Ｄ−ＰＥＴや、１層リングからなる２Ｄ−ＰＥＴにも適用できる。多層のリングからなる２Ｄ−ＰＥＴは、実質的には、１層リングからなる２Ｄ−ＰＥＴが複数あるものと同様に考えることができるので、後者について説明する。この場合も、測定対象物が動くことは、リングが動くことと同等であり、特に、測定対象物がリング中心軸方向に動くことは、リングがその中心軸方向に動くことと同等である。したがって、測定対象物が固定されている場合には、その１断面しか測定できないのに対して、測定対象物が動く場合には、その体動の範囲内で、多数の断面の測定が可能となる。すなわち、模擬的な３Ｄ−ＰＥＴとも言える。但し、実際にはリングが１層であるので、３Ｄ−ＰＥＴに比べれば検出感度は低い。したがって、２Ｄ−ＰＥＴであっても、上記実施形態と同様にして、エミッション計測最中に測定対象物の存在密度分布関数Ｃθψを存在密度分布生成部で生成するとともに、その体動を補償したエミッションデータＥをｘ−ｙ−θ−ψメモリに蓄積し、その体動を補償した感度補正データＢ１を（２）式で算出し、感度補正後のエミッションデータＥ１を（３）式で求めて、画像再構成を行うことができる。
【００７６】
また、リングが回転運動やウォブリング運動などの運動を行う場合であっても、同様に適用可能である。この場合、リングに固定した体動測定部によって測定対象物の位置・方位を測定し、或いは、リングの回転位置および測定対象物の位置・方位の双方を測定し、これによって求められるリングと測定対象物との相対的な動きに基づいて、上記実施形態と同様にして体動補償と感度補正とを施した真の線和を求めることができる。
【００７７】
また、上記実施形態では、存在密度分布生成部は、エミッション計測中に存在密度分布を蓄積・生成したが、エミッション計測中には各時刻における座標値（Δｘ，Δｙ，Δθ，Δψ）を時系列データとして蓄積しエミッション計測終了後に存在密度分布関数Ｃθψを算出してもよい。
【００７８】
【発明の効果】
以上、詳細に説明したとおり、本発明に係るポジトロンＣＴ装置では、基準感度補正データ生成手段によって、較正用線源回転手段により回転させられている較正用線源から発生した光子対がリングによって検出され蓄積された投影データが基準感度補正データとして獲得される。位置方位測定手段によって、測定空間に置かれた測定対象物のリングに対する相対的な位置および方位が測定されて位置方位データが出力され、測定手段によって測定対象物についての投影データが測定対象物の動きを補償して測定データとして獲得されるとともに、存在密度分布生成手段によって測定対象物の存在密度分布関数が生成される。そして、感度補正データ生成手段によって存在密度分布関数に基づいて測定対象物の動きを補償した感度補正データが生成され、感度補正手段によって測定データが感度補正され、画像再構成手段によって再構成画像を得る。以上のような構成とすることにより、測定対象物の体動を補償した測定データ（エミッションデータ）が得られるだけでなく、体動補償した感度補正データが得られ、測定中に測定対象物が固定されず動く場合であっても、正確な感度補正データによって感度補正を行うことができ、したがって、アーティファクト発生のない正確な再構成画像を得ることができる。
【図面の簡単な説明】
【図１】実施形態に係るＰＥＴのシステム構成図である。
【図２】実施形態に係るＰＥＴの検出部の断面図である。
【図３】実施形態に係るＰＥＴにおける体動測定の説明図である。
【図４】実施形態に係るＰＥＴにおける位置方位データおよびｘ−ｙ−θ−ψ空間上の位置の図表である。
【図５】実施形態に係るＰＥＴにおける測定対象物の各位置・各方位における存在密度分布の説明図である。
【図６】実施形態に係るＰＥＴにおける感度補正データ生成および感度補正のデータフロー図である。
【図７】実施形態に係るＰＥＴにおける計測のフローチャートである。
【図８】２Ｄ−ＰＥＴのシステム構成図である。
【図９】２Ｄ−ＰＥＴのリングにおける光子対検出と画像再構成の説明図である。
【図１０】２Ｄ−ＰＥＴのリングにおける光子対検出の説明図である。
【図１１】ブロック型ＰＥＴのリングの構成図である。
【図１２】３Ｄ−ＰＥＴのシステム構成図である。
【図１３】３Ｄ−ＰＥＴのリングの構成図である。
【図１４】エミッション計測の説明図である。
【図１５】ブランク計測の説明図である。
【図１６】トランスミッション計測の説明図である。
【符号の説明】
１００…測定対象物、１１０…較正用線源、１２０…検出部、１２１…リング、１２２…測定空間、１２３…体動測定部、１２４…遮蔽シールド、１３０…同時計数回路、１４０…ｘ−ｙ−θ−ψ変換部、１５０…ｘ−ｙ−θ−ψメモリ、１６０…画像再構成装置、１７０…画像表示部、１８０…記憶装置、１９０…存在密度分布生成部、２００…コンピュータ。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a positron CT apparatus that measures a substance distribution in a measurement object by detecting a pair of photons emitted along with annihilation of an electron-positron pair generated by an RI source introduced to the measurement object. And its image reconstruction method.
[0002]
[Prior art]
A positron emission computed tomography (hereinafter referred to as PET) is applied to a study of a living body or a disease or a clinical test, and is used for a positron emission nuclide (hereinafter referred to as an RI source) injected into a body. It is a device for imaging distribution and seeing biological functions.
[0003]
RI ray sources, biological materials such as FDG (18 ^F- fluorodeoxyglucose) involved in glucose metabolism in the dopamine and the body involved in neurotransmission, or for example, drugs in new development, partially added Used. PET makes it possible to observe the distribution, consumption, or temporal change of such a substance in a living body. PET can also measure basic metabolism of a living body such as cerebral blood flow and oxygen consumption.
[0004]
The detection section of such a PET is composed of a large number of photon detectors (referred to as rings) arranged in a ring shape, and inside the ring, an object to be measured such as a human body into which an RI ray source has been injected or inhaled. Is placed. Positrons emitted from the RI source in the object to be measured are immediately combined with nearby electrons, and a pair of photons (gamma rays) each having an energy of 511 keV are emitted in opposite directions. By simultaneously counting the pair of photons by the ring, it is possible to specify on which straight line the electron-positron annihilation has occurred. PET accumulates such coincidence counting information and performs image reconstruction processing to create a distribution image of the RI source.
[0005]
The time required for such PET measurement depends on the half-life of the RI source and may range from tens of minutes to several hours. Conventionally, the object to be measured has a fixed position and orientation during this measurement time. Had to be. When the measurement time is long as described above, it is extremely painful for the human body and other animals as the measurement object, and the position and orientation of the measurement object may change during the measurement time. is there. Accuracy cannot be measured if the measurement object moves even a little during measurement, and the image obtained by reconstructing the image will have an artifact due to the movement of the measurement object, and obtain an accurate reconstructed image I couldn't do that.
[0006]
In addition, even if the measurement target can be fixed, the measurement target is often measured in the awake state instead of the anesthesia state. Physiological state changes. This change has an adverse effect on the measurement, and again an accurate measurement cannot be made.
[0007]
Therefore, measurement is performed without maintaining the physiological state without fixing the measurement target, that is, by accepting movement (hereinafter referred to as body movement) relating to the position and orientation of the measurement target, and during the measurement, There is known a technique of measuring body movement and correcting coincidence counting information using the body movement information (for example, JP-A-2-209133 and JP-A-4-128679).
[0008]
In order to improve the spatial measurement resolution of the ring photon detector, the ring is rotated or wobbling, but the measurement resolution is also improved by moving the object to be measured. . The body motion correction technique is also effective in this regard.
[0009]
[Problems to be solved by the invention]
By employing the above-described body motion correction technology, it has become possible to obtain a more accurate reconstructed image than before even when the measurement target moves.
[0010]
On the other hand, in general, there is uneven sensitivity among a large number of photon detectors constituting a ring. As described above, a photon pair generated from an object to be measured is detected, and based on the coincidence counting information (emission data) accumulated as it is. If the image reconstruction processing is performed, an accurate RI source distribution image cannot be obtained. Thus, sensitivity correction data (blank data) is obtained using a calibration source, and sensitivity correction is performed by dividing the emission data by the blank data.
[0011]
However, when the measurement object moves, there is a problem in performing the sensitivity correction simply as described above. In other words, even if a photon pair is generated at the same position and flies in the same direction for the measurement target, if the position and direction of the measurement target are different, the photon pair will be different photon detector pairs having different detection sensitivities. Is received at. Nevertheless, in the conventional sensitivity correction, since the sensitivity correction is performed on the emission data with the same blank data, accurate sensitivity correction cannot be performed. Therefore, in order to accurately measure a moving measurement target, a sensitivity correction technique suitable for this case has been required together with a body motion correction technique.
[0012]
The present invention has been made in order to solve the above problems, and even when the measurement object moves during the measurement, the body motion is corrected, and the accurate sensitivity correction is performed by performing the accurate sensitivity correction. An object of the present invention is to provide a positron CT apparatus capable of obtaining a reconstructed image and a method for reconstructing the image.
[0013]
[Means for Solving the Problems]
The positron CT apparatus according to claim 1 detects a pair of photons generated due to annihilation of an electron-positron pair in the measurement space by a ring including a plurality of photon detectors arranged around a predetermined axis around the measurement space. Then, a predetermined number is cumulatively added to an address corresponding to a coordinate value expressed by polar coordinates set in a measurement space for a straight line connecting two photon detectors that have detected respective photons of the photon pair, based on the accumulated projection data. A positron CT apparatus for measuring the spatial distribution of the frequency of annihilation of an electron-positron pair, comprising: (1) a calibration source for generating a photon pair associated with the annihilation of an electron-positron pair; And (2) instructing the calibration source rotating means to rotate the calibration source so that the calibration source rotates. do it Means for detecting a pair of photons generated from the calibration source and acquiring accumulated projection data as reference sensitivity correction data in the state where the calibration is performed, and (3) an object to be measured placed in the measurement space. Position and orientation measurement means for outputting position and orientation data in accordance with the relative position and orientation of the ring with respect to (1), and (4) detecting a pair of photons generated from the measurement object in a state where there is no calibration source in the measurement space. Measuring means for compensating for relative position and azimuth changes based on the position and orientation data to obtain projection data accumulated during the measurement time as measurement data; and (5) inputting the position and orientation data during the measurement time The relative position and orientation represented by the position and orientation data to coordinate values on polar coordinates, and an existence density distribution generation method for generating an existence density distribution on the polar coordinates of the measurement object during the measurement time. And (6) sensitivity correction data generating means for generating sensitivity correction data based on a result of weighted average calculation of the reference sensitivity correction data based on the existence density distribution, and (7) sensitivity based on the sensitivity correction data based on the measured data. (8) Image reconstruction for calculating the spatial distribution of the frequency of electron-positron annihilation occurring in the measurement object based on the measurement data sensitivity-corrected by the sensitivity correction means, and performing image reconstruction. Means.
[0014]
The operation of this device is as follows. The reference sensitivity correction data generating means obtains, as reference sensitivity correction data, projection data in which a pair of photons generated from the calibration source rotated by the calibration source rotating means is detected by the ring and accumulated. The position and orientation measurement means measures the relative position and orientation of the measurement object placed in the measurement space with respect to the ring, and outputs position and orientation data.The measurement means outputs projection data about the measurement object to the measurement object. Is acquired as measurement data by compensating for the movement of the object, and an existence density distribution representing the movement distribution of the measurement target object is generated by the existence density distribution generation unit. Then, sensitivity correction data for compensating for the movement of the object to be measured is generated by the sensitivity correction data generating means, the sensitivity of the measurement data is accurately corrected by the sensitivity correcting means, and an accurate reconstructed image is obtained by the image reconstructing means.
[0015]
In the positron CT apparatus according to claim 2, the existence density distribution generating means accumulates the coordinate values at each time during the measurement time, and generates the existence density distribution after the measurement time elapses based on the accumulated coordinate values. In the positron CT apparatus according to the third aspect, the existence density distribution generating means sequentially generates the existence density distribution based on the coordinate values at each time during the measurement time.
[0016]
According to a fourth aspect of the present invention, there is provided an image reconstruction method for a positron CT apparatus, wherein a ring formed by a plurality of photon detectors arranged around a predetermined axis around a measurement space is generated when an electron-positron pair disappears in the measurement space. And a predetermined number is accumulated and added to an address corresponding to a coordinate value expressed by polar coordinates set in a measurement space with respect to a straight line connecting two photon detectors that have detected each photon of the photon pair. A method for reconstructing an image of a positron CT apparatus for measuring a spatial distribution of occurrence frequency of electron-positron pair annihilation based on projected data, comprising: (1) a calibration method for generating a photon pair accompanying the annihilation of an electron-positron pair With the source rotated relative to the ring about a predetermined axis in the measurement space, a pair of photons generated from the calibration source is detected, and the accumulated projection data is used as a reference sensitivity compensation. (1) measuring the relative position and orientation of a measurement object placed in the measurement space with respect to the ring without the calibration source in the measurement space; Detecting photon pairs generated from an object, compensating for changes in the relative position and orientation, and acquiring projection data accumulated during the measurement time as measurement data, and determining the relative position and orientation during the measurement time. A second step of mapping the coordinate values on the polar coordinates to generate an existence density distribution of the measurement object on the polar coordinates during the measurement time, and (3) a result of calculating the weighted average of the reference sensitivity correction data by the existence density distribution (4) a fourth step of correcting the sensitivity of the measurement data based on the sensitivity correction data, and (5) a sensitivity correction by the fourth step. A fifth step of calculating a spatial distribution of the frequency of occurrence of electron-positron annihilation in the measurement object based on the corrected measurement data and performing image reconstruction.
[0017]
The operation of the method is as follows. In the first step, a ring detects a photon pair generated from the calibration source rotated by the calibration source rotating means, and the accumulated projection data is obtained as reference sensitivity correction data. In the second step, position and orientation data according to the relative position and orientation of the measurement object placed in the measurement space with respect to the ring are acquired, and a density distribution representing the motion distribution of the measurement object is generated. The projection data of the measurement target is obtained as measurement data by compensating for the movement of the measurement target. Then, in a third step, sensitivity correction data that compensates for the movement of the measurement object is generated, in a fourth step, the sensitivity of the measurement data is accurately corrected, and in a fifth step, an accurate reconstructed image is obtained.
[0018]
In the image reconstruction method for a positron CT apparatus according to claim 5, the second step stores the coordinate values at each time during the measurement time, and calculates the existence density distribution after the measurement time has elapsed based on the stored coordinate values. Generate. In the image reconstruction method for a positron CT apparatus according to claim 6, the second step sequentially generates an existence density distribution based on coordinate values at each time during a measurement time.
[0019]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings. In the description of the drawings, the same elements will be denoted by the same reference symbols, without redundant description.
[0020]
Prior to the description of the main parts of the embodiment, the configuration, operation, image reconstruction method, sensitivity correction method, and the like of PET to which the embodiment is applied will be described.
[0021]
First, the configuration and operation of 2D-PET will be described. FIG. 8 is a system configuration diagram of 2D-PET.
[0022]
A ring in which n photon detectors D _k (k = 1, 2, 3,..., N) are arranged in a ring around a central axis surrounding a measurement object in the measurement space is a photon. Each of the detectors D _k (k = 1, 2, 3,..., N) and the coincidence circuit 30 are connected by a signal line. When one of the photon detectors of the ring 20 detects a photon, a signal corresponding to the photon energy is sent to the coincidence counting circuit 30 through a signal line, and the coincidence counting circuit 30 Based on the obtained signal, it is recognized that two photon detectors in the ring 20 have simultaneously detected a photon pair having a predetermined photon energy (511 keV) emitted from the RI source 10, and these two photon detectors at that time are detected. A detector identification signal (I, J) indicating each of the two photon detectors is output.
[0023]
These detector identification signals (I, J) are input to the t-θ conversion unit 40, and when a straight line L1 connecting two photon detectors that have detected a photon pair is expressed by polar coordinates in a measurement space in the ring 20. It is converted into a mapping position (T, θ ′) on a t-θ plane using each of the two variables (t, θ) used as a coordinate axis. Here, T is the distance between the straight line and the coordinate origin, and θ ′ is the angle of the perpendicular to the straight line from the coordinate origin with respect to the main axis of the coordinate system. The t-θ memory 50 cumulatively adds 1 to the projection data stored at the address corresponding to (T, θ ′). In this way, coincidence counting information on one pair of photons generated by the RI source 10 is stored in the t-θ memory 50.
[0024]
In this way, the coincidence information of a large number of photon pairs generated by the RI source 10 and detected by the ring 20 is stored as projection data in the t-θ memory 50. FIG. 9 is a diagram for explaining photon pair detection and image reconstruction in a ring. An image representing the distribution of the RI source 10 viewed from the θ ′ direction (direction of the straight line L1 in FIG. 9) of the polar coordinate system set in the measurement space in the ring 20 is a polar coordinate system of the accumulated projection data. Is reconstructed from the projection data corresponding to each T point on θ = θ ′ by the image reconstruction unit 60, and the reconstructed image is displayed on the image display unit 70.
[0025]
However, not all photon pairs generated in the measurement space are detected. FIG. 10 is an explanatory diagram of photon pair detection in a ring. Photon pairs from the RI source 10 are emitted from any location where the RI source 10 is present and fly in any direction, including those occurring on and along the ring 20 plane. Only the photon pairs that have flew to the ring 20 can be detected at the ring 20. For example, as shown in FIG. 9, even if the photon pair is generated from the RI source 10 existing in the ring 20, a photon pair flying in directions L2 to L4 other than the direction along the surface of the ring 20 is not detected. . The photon pair that has flown in the direction L1 along the surface of the ring 20 is detected by one of the photon detectors I and J, is converted to (T, θ ′) by the t-θ conversion unit 40, and is stored in the t-θ memory 50. Is added to the projection data at the address corresponding to.
[0026]
In addition, the photon detector has a block type configuration because it is easy to manufacture and inexpensive. In this type of ring, as shown in FIG. 11, each of a plurality of blocks B _p has a large number of photon detectors D _pq arranged in parallel, and these blocks B _p are arranged in a ring shape (p = 1,2,3, ..., q = 1,2,3, ...). This block type PET differs only in the arrangement of the photon detectors, and there is no difference in the operation principle, accumulation of projection data, image reconstruction, and the like from the above description.
[0027]
In the above-mentioned PET, since only the photon pairs flying in the direction along the ring surface are detected from the photon pairs emitted from the RI source and flying in all directions, the photon pairs emitted from the RI source are captured. There is a problem that the probability of detection is low, the detection sensitivity is low and the statistical noise is large. In order to improve the detection sensitivity, it is conceivable to inject a large amount of an RI source into the measurement target, but there is a limit when the measurement target is a living body.
[0028]
Therefore, in such a case, a three-dimensional PET (3D-PET) is used. FIG. 12 is a 3D-PET system configuration diagram, and FIG. 13 is a 3D-PET ring configuration diagram. The 3D-PET multilayer ring 21 is a _single- layer ring R ₁ , R ₂ ,... Of a photon detector similar to that shown in FIG. . . . Is obtained by multi-layer arranged R _m, the photon pairs flying in the direction of the straight line L1 emitted by RI radiation source 10, different single-layer ring _{_{R p, R q (p ≠}} q) 2 single photons belonging to each Coincidence can also be achieved by a detector. When any two photon detectors in the multilayer ring 21 simultaneously detect a pair of photons of energy 511 keV emitted from the RI source 10, the coincidence circuit 31 generates a detector identification indicating each of the two photon detectors. A signal (I, J) and a difference signal RD (Ring Difference) between the two single-layer rings to which the two photon detectors belong are output. These detector identification signals (I, J) and the ring difference signal RD are input to the xy-θ-ψ conversion unit 41, and a straight line L1 connecting two photon detectors that have detected a photon pair is formed into a multilayer ring 21. The mapping position (x, y, θ, 上の) on the xy-θ-ψ space with each of the four variables (x, y, θ, ψ) used when expressed in polar coordinates in the measurement space within ). Here, θ and ψ represent the direction of the straight line L1, and x and y represent positions in a rectangular coordinate system on a projection plane perpendicular to the straight line L1. The xy-θ-ψ memory 51 accumulatively adds 1 to the projection data stored at the address corresponding to (x, y, θ, ψ). In this way, coincidence information on a pair of photons generated by the RI source 10 is stored as projection data in the xy-θ-ψ memory 51, and the image reconstruction unit 61 uses the projection data to store the coincidence information. The reconstructed image is displayed on the image display unit 71.
[0029]
In any of the above types of PET, the detection sensitivity of each of the photon detectors is not constant and has many irregularities. Therefore, the sensitivity correction of the photon detector is performed as follows. That is, as shown in FIG. 14, the measurement is performed with the measurement target 10a to which the RI source is administered placed in the ring, and the coincidence counting information is accumulated via the coincidence circuit 32 and the t-θ converter 42. This is called emission measurement, and the data stored in the t-θ memory 52 is called emission data, and is represented by E (t, θ). Separately, as shown in FIG. 15, in a state where there is no object to be measured in the ring, the calibration RI-ray source 10b is rotated around the center axis of the ring in the ring, and the simulated parallel light is emitted. (Blank measurement) is performed by making each incident on each photon detector. The data (blank data B (t, θ)) stored in the t-θ memory 52 in this manner represents the detection sensitivity variation of the photon detector pair. Then, sensitivity correction is performed by dividing the emission data E (t, θ) by the blank data B (t, θ).
[0030]
Further, the photon absorption of the object to be measured is corrected as follows. As shown in FIG. 16, a measurement target 10c to which the RI source is not administered is placed in the ring at the same position as that at the time of emission measurement, and the calibration RI source 10b is rotated in the ring similarly to the blank measurement. (Transmission measurement), thereby acquiring data (transmission data T (t, θ)) stored in the t-θ memory 52. If the transmission data T (t, θ) is divided by the blank data B (t, θ), the absorption coefficient on the coincidence line is obtained. The emission data E (t, θ) is divided by the absorption coefficient to obtain the absorption coefficient. Can be corrected. In the transmission measurement, the statistical accuracy is significantly degraded because the weak RI radiation source is further absorbed by the object to be measured and the number of photon pairs detected is reduced. As a countermeasure, T (t, θ) / B (t, θ) is smoothed by a filter or the like.
[0031]
The true projection data P (t, θ) after the sensitivity correction and the absorption correction are performed as described above are emission data E (t, θ), blank data B (t, θ), and transmission data T (t, θ). θ),
P = (E / B) / <T / B> (1)
It is represented by Here, the symbol <> means smoothing. The same applies to the case of 3D-PET.
[0032]
The present invention is intended to perform body motion correction even for a moving measurement object and to accurately correct sensitivity of emission data, thereby obtaining an accurate reconstructed image.
[0033]
Next, the main part of the present embodiment will be described. This embodiment is a block type 3D-PET in which the ring is stationary. FIG. 1 is a system configuration diagram of the PET according to the present embodiment.
[0034]
The PET according to the present embodiment includes (1) a ring 121 including a large number of photon detectors for detecting photons, a body movement measurement unit 123 for measuring the position and orientation of the measurement target 100, and outputting position and orientation data, and (2) determining whether or not the photon detected by the ring 121 is a photon pair generated due to the annihilation of the electron / positron pair; Polar coordinates represent a straight line connecting the coincidence counting circuit 130 for identifying the photon detector pair that has detected the photon pair generated due to the annihilation of the electron / positron pair, and (3) the photon detector pair identified by the coincidence counting circuit 130. An xy-θ-ψ conversion unit 140 for converting into coordinate values (x, y, θ, ψ); and (4) a projection stored at an address corresponding to the coordinate values (x, y, θ, ψ). X-y- for cumulatively adding a predetermined value to data an θ-１５０ memory 150, (5) an existence density distribution generation unit 190 that receives the position and orientation data output from the body motion measurement unit 123 and generates an existence density distribution on the polar coordinates of the measurement target 100, and (6) ) An image reconstruction device 160 that generates reconstructed image data based on projection data that has been stored in the xy-θ-ψ memory 150 and that has been subjected to sensitivity correction, etc., and (7) a reconstructed image based on the reconstructed image data. And (8) a computer 200 that controls each of the above units, reads projection data stored in the xy-θ-ψ memory 150, and performs a predetermined process, and (9) A storage device 180 for storing data handled by the computer 200.
[0035]
A large number of photon detectors are arranged in a ring on the ring 121 of the detection unit 120 including the measurement space 122 therein, and these photon detectors are arranged in the direction of the measurement space 122 where the measurement target 100 is placed. The light receiving surface is turned to receive photons generated in the measurement space 122. The detection unit 120 further includes a body motion measurement unit (position and orientation measurement unit) 123, which measures the position and orientation of the measurement target 100 placed in the measurement space 122, and outputs position and orientation data. . Further, the detection unit 120 includes a calibration source 110 used for blank measurement and transmission measurement, and a rotation mechanism for rotating the calibration source 110. The rotation mechanism includes, for example, a support mechanism that supports the calibration source 110 and is rotatable about a central axis, a belt that transmits rotation to the support mechanism, and a motor that generates the rotation. Is done. In use (during blank measurement and transmission measurement), the calibration source 110 rotates around the center axis of the ring 121 on the surface of the ring 121, whereby the photon pairs generated from the calibration source 110 are rotated. Among them, the one flying on the surface of the ring 121 is detected by the photon detector pair. On the other hand, when it is not used (during emission measurement), the photon pair generated from the calibration radiation source 110 is evacuated to a position where it does not reach any photon detector. Alternatively, the calibration source 110 may be removed from the device.
[0036]
When any photon detector in the ring 121 receives a photon, a signal corresponding to the photon energy is sent to the coincidence circuit 130 through a signal line between each photon detector and the coincidence circuit 130. . The coincidence circuit 130 that has received the signal arriving from each photon detector of the ring 121 has a predetermined energy (511 keV) generated when two photon detectors of the ring 121 are annihilated by electron-positron pair annihilation. Recognizing that the photon pairs have been detected at the same time, a detector identification signal (I, J) indicating each of these two photon detectors at that time and an inter-ring difference signal RD are output.
[0037]
The detector identification signal (I, J) and the ring-to-ring difference signal RD are input to the xy-θ-ψ conversion unit 140, and two of the two photon pairs detected by the detector identification signal (I, J) are detected. Xy-θ-ψ with each of four variables (x, y, θ, 座標) used as a coordinate axis when expressing a straight line connecting the photon detectors in polar coordinates set in the measurement space 122 in the ring 121. It is mapped to a position in space. The xy-θ-ψ conversion unit 140 also receives the position and orientation data output from the body movement measurement unit 123, compensates for the body movement of the measurement target 100, and (x, y, θ, ψ ) Output the value.
[0038]
The xy-θ-ψ memory 150 that stores the projection data receives the (x, y, θ, ψ) value output from the xy-θ-ψ conversion unit 140, and receives the (x, y, ，) value. The predetermined value (for example, “1”) is cumulatively added to the projection data stored at the address corresponding to the θ, ψ) value. In this way, a large number of coincidence information on the photon pairs generated from the measurement object 100 is accumulated in the xy-θ-ψ memory 150 and becomes projection data.
[0039]
The position and orientation data output from the body movement measurement unit 123 is also input to the existence density distribution generation unit 190, and the position and orientation of the measurement target 100 in the measurement space 122 are determined by the xy-θ-ψ conversion unit 140. Similarly, it is mapped to a position on the xy-θ-ψ space. The position and orientation data of the measurement target 100 within a certain time is mapped one after another to the points on the xy-θ-ψ space, and the xy-θ of the measurement target 100 within the certain time is thus obtained. A density distribution in the − における space is generated.
[0040]
The image reconstructing device 160 for reconstructing an image of the distribution density of the RI source on the measurement object 100 uses projection data (emission data) stored in the xy-θ-ψ memory 150 or after the correction processing described later. From the emission data, reconstructed image data representing the photon pair generation distribution in the measurement target 100 viewed from a predetermined θψ direction of the polar coordinate system set in the measurement space 122 in the ring 121 is generated. The reconstructed image data is input and the reconstructed image is displayed.
[0041]
The components described above are integrated and controlled by the computer 200 via the bus line 210. That is, the computer 200 controls the rotation and retraction of the calibration radiation source 110. In addition, the existence density distribution function generated by the existence density distribution generation unit 190 is obtained. Further, all the contents of each address of the xy-θ-ψ memory 150 are cleared to zero before the measurement, and the projection data accumulated at each address is acquired after the measurement is completed. Further, based on the projection data (emission data, blank data, transmission data) acquired from the xy-θ-ψ memory 150, correction processing (sensitivity correction and absorption correction) described later is performed. In addition, emission data is transmitted to the image reconstruction device 160, and reconstructed image data generated by the image reconstruction device 160 is obtained. Further, the reconstructed image data is sent to the image display unit 170, and the reconstructed image is displayed.
[0042]
Each data handled in each process of the computer 200 is stored in the storage device 180, and is read from the storage device 180 when necessary. That is, the storage device 180 stores the emission data read from the xy-θ-ψ memory 150 after the emission measurement, and stores the existence density distribution function read from the existence density distribution generation unit 190. The storage device 180 also stores blank data obtained by blank measurement and transmission data obtained by transmission measurement. In calculating the sensitivity correction and the absorption correction, the storage device 180 reads emission data, blank data, transmission data, and the like, and stores intermediate data of the calculation and emission data after the sensitivity correction. In the image reconstruction processing, the storage device 180 reads the emission data after the sensitivity correction, and stores the reconstructed image data generated by the image reconstruction device. The storage device 180 has a storage area divided into a plurality of areas, and stores the data in different areas as needed. The storage device 180 may be a magnetic disk, a magnetic tape, or the like, or may be a main storage device including a semiconductor memory in the computer 200.
[0043]
The units described above do not operate independently, but operate in cooperation with each other according to instructions from the computer 200. That is, the computer 200 (1) obtains the reference sensitivity correction data generation module 201 that obtains blank data (reference sensitivity correction data) by rotating the calibration radiation source 110, and (2) obtains the emission data and the existence density distribution function. (3) a sensitivity correction data generation module 203 that generates sensitivity correction data that compensates for the body movement of the measurement object 100; (4) a sensitivity correction module 204 that performs sensitivity correction on emission data; An image reconstructing module 205 for performing an image reconstructing process and an image display, and each module performs integrated control of the operation of each unit in a corresponding processing step.
[0044]
Prior to the blank measurement, the reference sensitivity correction data generation module 201 instructs the rotation of the calibration radiation source 110 and clears the data in the xy-θ-ψ memory 150 to zero. Then, after the blank measurement is completed, the reference sensitivity correction data accumulated in the xy-θ-ψ memory 150 is obtained, and the reference sensitivity correction data is stored in the storage device 180.
[0045]
Prior to emission measurement, the measurement module 202 for acquiring emission data instructs the calibration source 110 to stop and evacuate, and clears the data in the xy-θ-ψ memory 150 to zero. After the emission measurement is completed, the emission data stored in the xy-θ-ψ memory 150 and the existence density distribution function generated by the existence density distribution generation unit 190 are obtained, and the emission data and the existence density distribution are obtained. And the function are stored in the storage device 180.
[0046]
The sensitivity correction data generation module 203 reads out the reference sensitivity correction data and the existence density distribution function stored in the storage device 180, generates sensitivity correction data based on these data, and stores the sensitivity correction data. The information is stored in the device 180.
[0047]
Further, the sensitivity correction module 204 reads the emission data and the sensitivity correction data stored in the storage device 180, and corrects the sensitivity of the emission data based on these data. Further, absorption correction may be performed. Then, the corrected emission data is stored in the storage device 180.
[0048]
Further, the image reconstruction module 205 reads out the sensitivity-corrected emission data stored in the storage device 180, sends the data to the image reconstruction device 160, and processes the reconstruction data obtained by the image reconstruction device 160. The image data is acquired, and the data is sent to the image display unit 170 to be displayed as an image.
[0049]
Next, the PET detection unit according to the present embodiment will be described in detail. FIG. 2 is a cross-sectional view of the PET detection unit according to the present embodiment. The detection unit 120 includes a ring 121 including a large number of photon detectors D, a body motion measurement unit 123 that measures the position and orientation of the measurement target 100, and the like.
[0050]
The ring 121 having the measurement space 122 therein has a structure in which block-type single-layer rings are arranged in multiple layers (five layers in this figure), and each layer is a block including a plurality of photon detectors D. Are arranged in a ring shape, and each photon detector D of each block has its light receiving surface directed in the direction of the measurement space 122 where the measurement object 100 is placed. The shielding shield 124 is for preventing the photon pairs from leaking out of the measurement space 122, and does not separate the single-layer rings from each other. Therefore, it is possible to detect photon pairs between different single-layer rings.
[0051]
The body movement measurement unit 123 that measures the position and orientation of the measurement target 100 is, for example, an optical distance measurement sensor 123a to 123c including a light emitting element and a light receiving element, and a position and orientation based on outputs from these sensors. The body motion data processing unit 123d outputs data. Each of the distance measuring sensors 123a to 123c has a fixed relative positional relationship with the ring 121, and is provided at a predetermined position of the measurement object 100 (in this figure, the head of the human body which is the measurement object 100). Each of the markers 101a to 101c is irradiated with a light beam to receive the reflected light, and the distance to each of the markers 101a to 101c is measured. The distance data obtained from each of the three distance measurement sensors 123a to 123c is input to the body motion data processing unit 123d, and the position and orientation of the measurement target 100 are obtained.
[0052]
Further, an image capturing camera may be used as the body movement measuring unit 123, and the obtained image may be analyzed to obtain the position and the direction of the measurement target 100. Alternatively, an acceleration sensor provided at a predetermined position of the measurement target 100 may be used as the body movement measurement unit 123. In this case, the position and orientation of the measurement target 100 are determined based on the output from the acceleration sensor. You can ask.
[0053]
Next, the measurement of the position and orientation of the measurement target 100 and the generation of the existence density distribution by the body motion measurement unit 123 will be described in more detail. FIG. 3 is an explanatory diagram of body motion measurement in PET according to the present embodiment.
[0054]
The position of the measurement target 100 is a coordinate value of a predetermined point (for example, the center point of the measurement target 100) in the measurement target 100 in the XYZ rectangular coordinate system set in the measurement space 122 in the ring 121. It is represented by (X, Y, Z). The azimuth of the measuring object 100 is represented by the direction (Θ, Ψ) of the measuring object 100 in the XYZ orthogonal coordinate system. Here, Θ is the angle of the azimuth projected on the XY plane from the X axis, and Ψ is the elevation angle from the XY plane. The position and orientation data (X, Y, Z, Θ, Ψ) representing the position and orientation of the measurement target 100 are based on the distance data output from the distance measurement sensors 123a to 123c in the body motion data processing unit 123d. Desired.
[0055]
Note that the body motion measurement unit 123 may measure and output a displacement amount from a reference value (for example, the position and orientation of the measurement target 100 at the start of measurement) with respect to the position and orientation of the measurement target 100. The body movement measurement unit 123 measures the displacement of the position and orientation of the measurement target 100 at regular intervals (for example, several milliseconds to several tens of milliseconds) from the measurement start time to the measurement end time of the emission measurement. , Position and orientation data (ΔX _k , ΔY _k , ΔZ _k , ΔΘ _k , ΔΨ _k ) at each time tk are _obtained by the body motion measuring unit 123 (k = 0, 1,..., N,...). . The position and orientation data indicating the position and orientation of the measurement target 100 measured in this way are input to the xy-θ-ψ conversion unit 140 and are also input to the existence density distribution generation unit 190.
[0056]
Position orientation data at each time _{_{_{t k (ΔX k, ΔY k}}} , ΔZ k, ΔΘ k, ΔΨ k) the density distribution generating unit 190 inputs the in the same manner as x-y-θ-ψ converting section 140, Each position and orientation data is mapped to a position (Δx _k , Δy _k , Δθ _k , Δψ _k ) in the xy-θ-ψ space (k = 0, 1,..., N,...). FIG. 4 is a table showing position and orientation data and positions in the xy-θ-ψ space in PET according to the present embodiment. In this way, by determining the distribution of the mapping position in the xy-θ-ψ space within a certain time, the presence of the measurement object 100 representing the existence density distribution in the xy-θ-ψ space is obtained. A density distribution function Cθψ (Δx, Δy, Δθ, Δψ) is generated.
[0057]
The existence density distribution function Cθψ (Δx, Δy, Δθ, Δψ) is obtained by detecting, at the ring 121, a pair of photons generated at the same point on the measurement object 100 and flying in the same direction, and the coincidence count information is xy- When the measurement object 100 is at the reference position / reference azimuth (for example, the position / azimuth at the start of emission measurement) when accumulated in the θ-ψ memory 150, the coincidence information is represented by coordinate values (x, y, θ, ψ), but when the measurement object 100 is at the position or orientation where the displacement (Δx, Δy, Δθ, Δψ) occurs during the measurement, the coincidence information is represented by the coordinate value (x + Δx, It is stored at an address corresponding to y + Δy, θ + Δθ, 番 + Δψ, and indicates that the ratio is a value Cθψ (Δx, Δy, Δθ, Δψ).
[0058]
Note that the existence density distribution function Cθψ generally varies depending on the θ value and the ψ value, that is, depending on the projection direction. FIG. 5 is an explanatory diagram of the existence density distribution at each position and each direction of the measurement target in the PET according to the present embodiment. This diagram is a diagram viewed from the top of the measurement object (human head) 100, in which the existence density distribution functions Cθψ in the projection directions of La, Lb, and Lc are expressed as Ca, Cb, and Cc, respectively. It is. Although the existence density distribution function is actually a function of four variables, it is simplified in FIG. As shown in this figure, when the body movement of the measurement object 100 is biased in a certain direction or a certain rotation direction, the existence density distribution function Cθψ becomes a function different depending on θ and ψ.
[0059]
Next, a method of compensating for the body movement of the measurement object measured by the body movement measuring unit 123 to generate sensitivity correction data and further correcting the sensitivity will be described. The sensitivity correction data generation and sensitivity correction described below are processed by the sensitivity correction data generation module 203 and the sensitivity correction module 204 of the computer 200, respectively. FIG. 6 is a data flow diagram of sensitivity correction data generation and sensitivity correction in PET according to the present embodiment.
[0060]
The existence density distribution function Cθψ (Δx, Δy, Δθ, Δψ) generated by the existence density distribution generation unit 190 during emission measurement, and the reference sensitivity correction data B (x, y) obtained by the reference sensitivity correction data generation module 201 , Θ, ψ), after the emission measurement is completed, the sensitivity correction data B1 (x, y, θ, ψ) that compensates for the body movement of the measurement object 100 is
[0061]
(Equation 1)

[0062]
It is determined by the following relational expression. Here, the quadruple sum operation Σ is calculated over the range of the body movement of the measurement object 100. This calculation may detect the same photon pair from the measurement object 100 even though the same photon pair may be detected by a different photon detector pair according to the movement of the measurement object 100. This means that the averaged detection sensitivity is calculated in consideration of the weight represented by the existence density distribution Cθψ for the photon detector pair.
[0063]
As described above, the sensitivity of the emission data E (x, y, θ, ψ) obtained by the emission measurement is corrected by the sensitivity correction data B1 (x, y, θ, ψ) that compensates for the body movement of the measurement object 100. Accordingly, the emission data E1 (x, y, θ, ψ) after the sensitivity correction is
E1 = E / B1 (3)
Is obtained. The emission data E1 is obtained not only by performing the emission measurement by correcting the body movement of the measuring object 100 but also by correcting the sensitivity by using blank data B1 in which the body movement of the measuring object 100 is compensated. If an image is reconstructed based on this, an accurate reconstructed image can be obtained.
[0064]
If the absorption correction is to be further performed, the transmission measurement is performed in the same manner as the emission measurement, and the transmission data T subjected to body motion compensation and the existence density distribution function during the transmission measurement are obtained. The sensitivity correction data B2 is calculated according to the formula (2). And the true line sum P1 is
P1 = (E / B1) / <T / B2> (4)
It is obtained by the following relational expression. Here, the symbol <> means smoothing.
[0065]
Next, the operation of PET and the image reconstruction method according to the present embodiment will be described. FIG. 7 is a flowchart of measurement in PET according to the present embodiment.
[0066]
First, in step S1, blank measurement is performed according to an instruction from the reference sensitivity correction data generation module 201 of the computer 200. In other words, the calibration source 110 is rotated on the surface of the ring 121 around the central axis without placing the measurement object 100 in the ring 121 to generate simulated parallel light. In this state, the photon pairs generated from the calibration source 110 are detected by a large number of photon detectors constituting the ring 121, energy discriminated by the coincidence circuit 130, and the xy-θ-ψ conversion unit 140 Are converted to (x, y, θ, ψ) values corresponding to a straight line connecting the two photon detectors that have detected the photon pairs, and xy− corresponding to the (x, y, θ, ψ) values. The projection data is stored in the address of the θ- θ memory 150 as projection data. The reference sensitivity correction data B (x, y, θ, ψ) thus accumulated as projection data is read out by the reference sensitivity correction data generation module 201 and stored in the storage device 180.
[0067]
Subsequently, in step S2, emission measurement is performed according to an instruction from the measurement module 202 of the computer 200. That is, the calibration source 110 is evacuated, and the measurement target 100 into which the RI source is injected is placed in the measurement space 122 in the ring 121. In this state, the photon pair generated from the measurement target 100 is detected. Then, the body movement is corrected based on the position and orientation data indicating the position and orientation of the measurement object 100 output from the body movement measurement unit 123, and the projection data is accumulated in the xy-θ-ψ memory 150. At the same time, the position and orientation data is input to the existence density distribution generation unit 190, and an existence density distribution function Cθψ (Δx, Δy, Δθ, Δψ) is generated. The emission data E (x, y, θ, ψ) and the existence density distribution function Cθψ (Δx, Δy, Δθ, Δψ) thus accumulated as projection data are read out by the measurement module 202 and stored in the storage device. 180.
[0068]
If the absorption correction is also performed, transmission measurement is performed. That is, the measurement object 100 into which the RI source is not injected is placed in the measurement space 122, and the calibration source 110 is rotated. In this state, the measurement is performed in the same manner as the emission measurement, and the transmission data T and the existence density distribution are measured. The function is obtained and stored in the storage device 180.
[0069]
Note that the blank measurement (Step S1), the emission measurement (Step S2), and the transmission measurement may be performed in any order.
[0070]
Upon completion of each of the above measurements, subsequently, in step S3, sensitivity correction data is generated based on the data obtained in these measurements. When generating the sensitivity correction data, the sensitivity correction data generation module 203 of the computer 200 reads the existence density distribution function Cθψ and the reference sensitivity correction data B (x, y, θ, ψ) at the time of emission measurement from the storage device 180. , The sensitivity correction data B1 (x, y, θ, ψ) is calculated based on the equation (2).
[0071]
Subsequently, sensitivity correction is performed in step S4. The sensitivity correction module 204 of the computer 200 divides the emission data E (x, y, θ, ψ) read from the storage device 180 by the sensitivity correction data B1 (x, y, θ, ψ) obtained in step S3. Thus, the sensitivity unevenness between the plurality of photon detectors constituting the ring 121 is corrected, and the emission data E1 (x, y, θ, 感度) after the sensitivity correction is obtained.
[0072]
If the absorption correction is further performed, the transmission data T is corrected in the same manner as in step S3 based on the existence density distribution function at the time of transmission measurement, and the true projection data P1 is calculated based on equation (4). And calculate.
[0073]
Then, based on the emission data E1 after sensitivity correction obtained in step S4 (or the true projection data P1), image reconstruction is performed in step S5. That is, the image reconstruction module 205 of the computer 200 sends the emission data E1 (or the true projection data P1) after the sensitivity correction to the image reconstruction device 160, and the reconstruction generated by the image reconstruction device 160. The image data is obtained, and the reconstructed image data is sent to the image display unit 170 to display the reconstructed image. This makes it possible to more accurately measure the spatial distribution of the frequency of annihilation of electron-positron pairs.
[0074]
The present invention is not limited to the above-described 3D-PET, but can be applied to other types of PET.
[0075]
For example, the present invention can be applied to 2D-PET including a multilayer ring and 2D-PET including a single-layer ring. A 2D-PET made of a multi-layered ring can be considered substantially the same as a plurality of 2D-PETs made of a single-layered ring, so the latter will be described. Also in this case, the movement of the measurement object is equivalent to the movement of the ring. In particular, the movement of the measurement object in the center axis direction of the ring is equivalent to the movement of the ring in the center axis direction. Therefore, when the object to be measured is fixed, only one cross section can be measured. On the other hand, when the object to be measured moves, many cross sections can be measured within the range of the body movement. Become. That is, it can be said that it is a simulated 3D-PET. However, since the ring is actually one layer, the detection sensitivity is lower than that of 3D-PET. Therefore, even in the case of 2D-PET, the existence density distribution function Cθψ of the measurement object is generated by the existence density distribution generation unit during the emission measurement and the body motion is compensated for in the same manner as in the above embodiment. The data E is stored in the xy-θ-ψ memory, the sensitivity correction data B1 that compensates for the body motion is calculated by the equation (2), and the emission data E1 after the sensitivity correction is obtained by the equation (3). Image reconstruction can be performed.
[0076]
Further, the present invention can be similarly applied to a case where the ring performs a motion such as a rotation motion and a wobbling motion. In this case, the position and orientation of the object to be measured are measured by the body motion measurement unit fixed to the ring, or both the rotational position of the ring and the position and orientation of the object to be measured are measured. Based on the relative movement with respect to the object, a true line sum subjected to body motion compensation and sensitivity correction can be obtained in the same manner as in the above embodiment.
[0077]
Further, in the above embodiment, the existence density distribution generating unit accumulates / generates the existence density distribution during the emission measurement. However, during the emission measurement, the existence density distribution generation unit chronologically stores the coordinate values (Δx, Δy, Δθ, Δψ) at each time. The existence density distribution function Cθψ may be calculated after accumulating the data and completing the emission measurement.
[0078]
【The invention's effect】
As described above in detail, in the positron CT apparatus according to the present invention, the reference sensitivity correction data generation unit detects the pair of photons generated from the calibration source rotated by the calibration source rotation unit by the ring. The accumulated projection data is obtained as reference sensitivity correction data. The position and orientation measurement means measures the relative position and orientation of the measurement object placed in the measurement space with respect to the ring and outputs position and orientation data, and the measurement means outputs projection data of the measurement object to the measurement object. The motion is compensated for and obtained as measurement data, and the density distribution generating means generates the density distribution function of the object to be measured. Then, sensitivity correction data in which the movement of the measurement target is compensated for is generated by the sensitivity correction data generation means based on the existence density distribution function, the sensitivity correction means corrects the measurement data, and the image reconstruction means converts the reconstructed image. obtain. With the above configuration, not only the measurement data (emission data) in which the body movement of the measurement object is compensated but also the sensitivity correction data in which the body movement is compensated can be obtained, and the measurement object can be measured during the measurement. Even in the case of moving without being fixed, the sensitivity can be corrected using the accurate sensitivity correction data, and therefore, an accurate reconstructed image without artifacts can be obtained.
[Brief description of the drawings]
FIG. 1 is a system configuration diagram of a PET according to an embodiment.
FIG. 2 is a cross-sectional view of a PET detection unit according to the embodiment.
FIG. 3 is an explanatory diagram of body motion measurement in PET according to the embodiment.
FIG. 4 is a table showing position and orientation data and positions in xy-θ-ψ space in PET according to the embodiment.
FIG. 5 is an explanatory diagram of an existence density distribution at each position and each direction of a measurement target in PET according to the embodiment.
FIG. 6 is a data flow diagram of sensitivity correction data generation and sensitivity correction in PET according to the embodiment.
FIG. 7 is a flowchart of measurement in PET according to the embodiment.
FIG. 8 is a system configuration diagram of 2D-PET.
FIG. 9 is an explanatory diagram of photon pair detection and image reconstruction in a 2D-PET ring.
FIG. 10 is an explanatory diagram of photon pair detection in a 2D-PET ring.
FIG. 11 is a configuration diagram of a ring of a block type PET.
FIG. 12 is a system configuration diagram of 3D-PET.
FIG. 13 is a configuration diagram of a ring of 3D-PET.
FIG. 14 is an explanatory diagram of emission measurement.
FIG. 15 is an explanatory diagram of blank measurement.
FIG. 16 is an explanatory diagram of transmission measurement.
[Explanation of symbols]
Reference Signs List 100: measurement object, 110: calibration source, 120: detection unit, 121: ring, 122: measurement space, 123: body movement measurement unit, 124: shielding shield, 130: coincidence circuit, 140: xy -Θ-ψ conversion unit, 150: xy-θ-ψ memory, 160: image reconstruction device, 170: image display unit, 180: storage device, 190: existence density distribution generation unit, 200: computer.

Claims

A ring composed of a plurality of photon detectors arranged around a predetermined axis surrounding the measurement space detects a photon pair generated along with the annihilation of the electron-positron pair in the measurement space, and detects each photon of the photon pair. A predetermined number is cumulatively added to an address corresponding to a coordinate value expressed in polar coordinates set in the measurement space for a straight line connecting the two photon detectors that have detected the photon detectors, and the electron-positron pair annihilation is performed based on the accumulated projection data. Positron CT apparatus for measuring the spatial distribution of the frequency of occurrence of
Calibration source rotating means for rotating a calibration source for generating photon pairs due to annihilation of electron / positron pairs relative to the ring about the predetermined axis in the measurement space,
Instruct the calibration source rotating means to rotate the calibration source, and in a state where the calibration source is rotating, detect a photon pair generated from the calibration source. Reference sensitivity correction data generating means for acquiring the accumulated projection data as reference sensitivity correction data,
Position and orientation measurement means for outputting position and orientation data according to the relative position and orientation of the measurement object placed in the measurement space with respect to the ring,
In the state where the calibration source is not present in the measurement space, a photon pair generated from the measurement object is detected, and the change of the relative position and the direction is compensated based on the position and orientation data during the measurement time. Measurement means for acquiring the projection data accumulated in the measurement data as measurement data;
Input the position and orientation data during the measurement time and map the relative position and orientation represented by the position and orientation data to coordinate values on the polar coordinates, and on the polar coordinates of the measurement object during the measurement time. A density distribution generating means for generating a density distribution in
Based on the result of the weighted average calculation of the reference sensitivity correction data in the existence density distribution, a sensitivity correction data generating unit that generates sensitivity correction data,
Sensitivity correction means for correcting the measurement data based on the sensitivity correction data,
Image reconstruction means for calculating the spatial distribution of the frequency of occurrence of electron-positron annihilation in the measurement object based on the measurement data corrected for sensitivity by the sensitivity correction means, and performing image reconstruction.
A positron CT apparatus comprising:

The existence density distribution generating means accumulates the coordinate values at each time during the measurement time, and generates the existence density distribution after the measurement time elapses based on the accumulated coordinate values. The positron CT apparatus according to claim 1.

2. The positron CT apparatus according to claim 1, wherein the existence density distribution generation unit sequentially generates the existence density distribution based on the coordinate values at each time during the measurement time. 3.

A ring composed of a plurality of photon detectors arranged around a predetermined axis surrounding the measurement space detects a photon pair generated along with the annihilation of the electron-positron pair in the measurement space, and detects each photon of the photon pair. A predetermined number is cumulatively added to an address corresponding to a coordinate value represented by the polar coordinate set in the measurement space for a straight line connecting the two photon detectors that have detected the photon detectors, and the electron-positron pair annihilation is performed based on the accumulated projection data. An image reconstruction method of a positron CT apparatus for measuring a spatial distribution of occurrence frequency of
Generated from the calibration source in a state where the calibration source for generating a photon pair due to the annihilation of an electron / positron pair is rotated relative to the ring about the predetermined axis in the measurement space. A first step of detecting a photon pair to be obtained and acquiring the accumulated projection data as reference sensitivity correction data;
In the absence of the calibration source in the measurement space, the relative position and orientation of the measurement object placed in the measurement space with respect to the ring are measured, and a photon pair generated from the measurement object is detected. The projection data accumulated during the measurement time is obtained as measurement data by compensating for the change in the relative position and orientation, and the relative position and orientation are coordinated on the polar coordinates during the measurement time. Mapping to a value, generating a density distribution on the polar coordinates of the measurement object during the measurement time, a second step,
A third step of generating sensitivity correction data based on a result of weighted average calculation of the reference sensitivity correction data with the existence density distribution;
A fourth step of correcting the measurement data for sensitivity based on the sensitivity correction data;
A fifth step of calculating a spatial distribution of the frequency of occurrence of electron-positron pair annihilation in the measurement object based on the measurement data corrected for sensitivity in the fourth step, and performing image reconstruction;
An image reconstruction method for a positron CT apparatus, comprising:

The second step accumulates the coordinate values at each time during the measurement time, and generates the existence density distribution after the lapse of the measurement time based on the accumulated coordinate values. Item 5. An image reconstruction method for a positron CT apparatus according to Item 4.

5. The image reconstruction method for a positron CT apparatus according to claim 4, wherein in the second step, the existence density distribution is sequentially generated based on the coordinate values at each time during the measurement time.