(57)【要約】【課題】効果的なタイミング校正を可能とする。【解決手段】実施形態のイメージングシステムは、ポジトロン放出放射性同位元素線源と、対消滅ターゲットと、検出器要素と、計算部と、校正部とを含む。ポジトロン放出放射性同位元素線源は、イメージングシステムのイメージング領域の中または近傍に配置される。対消滅ターゲットは、ポジトロン放出放射性同位元素線源から所定の位置に配置される。検出器要素は、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出する。計算部は、検出された同時計数事象の対に基づき検出器要素用の校正時間オフセットを計算する。校正部は、計算された校正時間オフセットで検出器要素を校正する。ポジトロン放出放射性同位元素線源は円筒形状の対消滅ターゲットの内側に設けられる。【選択図】図１１
2020034568.tif 000002
【特許請求の範囲】
【請求項１】
イメージングシステムであって、
前記イメージングシステムのイメージング領域の中または近傍に配置されたポジトロン放出放射性同位元素線源と、
前記ポジトロン放出放射性同位元素線源から所定の位置に配置された円筒形状の対消滅ターゲットと、
前記対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出する検出器要素と、
前記検出された同時計数事象の対に基づき前記検出器要素用の校正時間オフセットを計算する計算部と、
前記計算された校正時間オフセットで前記検出器要素を校正する校正部とを含み、
前記ポジトロン放出放射性同位元素線源は前記円筒形状の対消滅ターゲットの内側に設けられる、
イメージングシステム。
【請求項２】
前記円筒形状の直径は少なくとも１０ｃｍ以上である、請求項１に記載のイメージングシステム。
【請求項３】
前記ポジトロン放出放射性同位元素線源と前記対消滅ターゲットを含むファントムは中空形状である、請求項１又は２に記載のイメージングシステム。
【請求項４】
前記ポジトロン放出放射性同位元素線源は、前記円筒形状の前記対消滅ターゲットの内面上に実質的に均一に堆積される、請求項１乃至３のいずれか一項に記載のイメージングシステム。
【請求項５】
前記ポジトロン放出放射性同位元素線源は、前記円筒形状の対消滅ターゲットの内面に塗布されている、請求項１乃至３のいずれか一項に記載のイメージングシステム。
【請求項６】
液状のポジトロン放出放射性同位元素線源が、ゼラチンまたは粘着タイプの材料と混合される、あるいは前記円筒形状の前記対消滅ターゲットの内面は吸収材料から作られる、請求項１乃至３のいずれか一項に記載のイメージングシステム。
【請求項７】
前記計算部は、前記ポジトロン放出放射性同位元素線源の大きさに基づいて前記イメージング領域の大きさを決定する、請求項１乃至６のいずれか一項に記載のイメージングシステム。
【請求項８】
前記計算部は、前記検出された同時計数事象の対から少なくとも１つの対タイミング差ヒストグラムを生成するよう、かつ前記少なくとも１つの対タイミング差ヒストグラムを用いて、前記検出器要素用の校正時間オフセットを計算する、請求項１に記載のイメージングシステム。
【請求項９】
前記校正時間オフセットは、前記対タイミング差ヒストグラムの重心に対応する、請求項８に記載のイメージングシステム。
【請求項１０】
前記計算部は、差分最小化技術を適用することにより前記校正時間オフセットを計算する、請求項８に記載のイメージングシステム。
【請求項１１】
前記計算部は、前記対タイミング差ヒストグラムの測定された重心と各検出器要素に対する予測される対のタイミング重心との差が予め設定された閾値未満になるまで、または各検出器要素に対する前記対タイミング差ヒストグラムの重心の絶対値が予め設定された閾値未満になるまで、検出器要素毎に前記校正時間オフセットを繰り返し更新する、請求項９に記載のイメージングシステム。
【請求項１２】
前記計算部は、同時計数線毎に前記対タイミング差ヒストグラムを複数の領域に分割し、各領域は、前記対消滅ターゲットと前記同時計数線の１つの交点を表す、請求項８に記載のイメージングシステム。
【請求項１３】
前記計算部は、交点から前記同時計数線の端点を表す各検出器要素までの距離に基づき計算される異なる時間シフトを、各領域に適用する、請求項１２に記載のイメージングシステム。
【請求項１４】
前記計算部は、前記対消滅ターゲットと任意の同時計数線の複数の交点が存在する場合、後続の交点との間隔が、前記対タイミング差ヒストグラムの検査から互いに異なる領域を識別することができるように十分に広くなるよう、同時計数対を形成するための前記イメージング領域を制限する、請求項１２に記載のイメージングシステム。
【請求項１５】
前記ポジトロン放出放射性同位元素線源はＧｅ−６８／Ｇａ−６８、Ｎａ−２２、Ｆ−１８、ＦＤＧおよびＲｂ−８２のうちの１つである、請求項１乃至１４のいずれか一項に記載のイメージングシステム。
【請求項１６】
前記対消滅ターゲットの位置と向きを決定するために、前記対消滅ターゲットを画像化する処理部を更に備える、請求項１乃至１５のいずれか一項に記載のイメージングシステム。
【請求項１７】
前記計算部は、前記画像から決定された前記位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算する、請求項１６に記載のイメージングシステム。
【請求項１８】
前記計算部は、前記検出された同時計数事象の対を用いてウォーク補正係数を計算し、かつ前記校正部は、前記ウォーク補正係数を用いて前記検出器要素を校正する、請求項１乃至１７のいずれか一項に記載のイメージングシステム。
【請求項１９】
前記ポジトロン放出放射性同位元素線源は、実質的に前記イメージングシステムのアイソセンタに配置され、
前記計算部は、前記ポジトロン放出放射性同位元素線源における、または線源ホルダーにおける対消滅から生じる前記同時計数事象の対を用いてウォーク補正係数を計算し、かつ
前記校正部は、前記ウォーク補正係数を用いて前記検出器要素を校正する、請求項１又は２に記載のイメージングシステム。
【請求項２０】
前記計算部は、反対側の検出器対で検出された事象が個別に処理されるように、前記ウォーク補正係数の計算を並行処理にする、請求項１９に記載のイメージングシステム。
【請求項２１】
イメージングシステムのイメージング領域の中または近傍に配置されたポジトロン放出放射性同位元素線源と、
前記ポジトロン放出放射性同位元素線源から所定の位置に配置された円筒形状の対消滅ターゲットと、
を含むファントムであって、
前記ポジトロン放出放射性同位元素線源は前記円筒形状の対消滅ターゲットの内側に設けられ、
前記イメージングシステムは、前記対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出する検出器要素と、前記検出された同時計数事象の対に基づき前記検出器要素用の校正時間オフセットを計算する計算部と、前記計算された校正時間オフセットで前記検出器要素を校正する校正部と、を含む、
ファントム。
【請求項２２】
イメージングシステムを校正する校正方法であって、
前記イメージングシステムのイメージング領域の中または近傍にポジトロン放出放射性同位元素線源を配置するステップと、
前記ポジトロン放出放射性同位元素線源から所定の位置に円筒形状の対消滅ターゲットを配置するステップと、
前記対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出するステップと、
前記検出された同時計数事象の対に基づき前記イメージングシステムにおける検出器要素用の校正時間オフセットを計算するステップと、
前記計算された校正時間オフセットで前記検出器要素を校正するステップとを含み、
前記ポジトロン放出放射性同位元素線源を配置するステップは、前記ポジトロン放出放射性同位元素線源を前記円筒形状の対消滅ターゲットの内側に設けることを含む、
校正方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、イメージングシステム、ファントムおよび校正方法に関する。
【背景技術】
【０００２】
ＰＥＴ画像の品質は、飛行時間（Time-Of-Flight：ＴＯＦ）ポジトロン放射断層撮影（Positron  Emission  Tomography：ＰＥＴ）技術を用いることにより改善されてきた。しかし、検出器の正確なタイミング校正は、ＴＯＦ  ＰＥＴシステムにとって極めて重要である。正確な、ロバストな、そしてあまり時間を必要としない結晶ベースのタイミング校正法が、ＴＯＦ  ＰＥＴのための最良のタイミング情報を得るために必要である。
【０００３】
現在のタイミング校正法は、（１）光パルスを光電子増倍管（Photomultiplier  Tube：ＰＭＴ）に導入する方法または電気パルスをプリアンプに導入する方法と、（２）ＰＭＴに結合されたプラスチックシンチレータに埋め込まれた放射性線源を使用する方法と、（３）回転する線状放射性線源を使用する方法と、（４）散乱ファントム内の点状放射性線源を使用する方法とがある。
【０００４】
図１は、放射性同位元素から放出されているポジトロンの例を示す。図１に示すように、従来のタイミング校正は、ポジトロン３と電子４が消滅（対消滅（annihilation））した後の同時発生の５１１ｋｅＶのガンマ線５を用いて実行される。図１に示すように、ポジトロン３は、放射性同位元素２からかなり大きい運動エネルギーで放出され（通常、〜ＭｅＶ）、対消滅前にいくらかの距離を移動する。従来技術では、一般的に、線源１および線源周りの材料は、放射性同位元素において、または通常ポジトロン放出同位元素２の数ミリメートル内にある周囲の材料において、対消滅が発生するように配置される。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−２８１８１６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明が解決しようとする課題は、効果的なタイミング校正を可能とするイメージングシステムおよび校正方法を提供することである。
【課題を解決するための手段】
【０００７】
実施形態のイメージングシステムは、ポジトロン放出放射性同位元素線源と、対消滅ターゲットと、検出器要素と、計算部と、校正部とを含む。ポジトロン放出放射性同位元素線源は、イメージングシステムのイメージング領域の中または近傍に配置される。対消滅ターゲットは、円筒形状であり、ポジトロン放出放射性同位元素線源から所定の位置に配置される。検出器要素は、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出する。計算部は、検出された同時計数事象の対に基づき検出器要素用の校正時間オフセットを計算する。校正部は、計算された校正時間オフセットで検出器要素を校正する。ポジトロン放出放射性同位元素線源は円筒形状の対消滅ターゲットの内側に設けられる。
【図面の簡単な説明】
【０００８】
【図１】図１は、放射性同位元素から放出されているポジトロンの例を示す図である。
【図２】図２は、対消滅位置の分布に対応する１／１０値幅プロットを図示した図である。
【図３】図３は、一実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す図である。
【図４】図４は、別の実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す図である。
【図５】図５は、別の実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す図である。
【図６】図６は、一実施形態による複数のポジトロン放出線源およびシールドを有するイメージングシステムを示す図である。
【図７Ａ】図７Ａは、複数の追加実施形態によるシールドおよびスキャナカバーを有していないイメージングシステムを示す図である。
【図７Ｂ】図７Ｂは、複数の追加実施形態によるシールドおよびスキャナカバーを有していないイメージングシステムを示す図である。
【図８Ａ】図８Ａは、ポジトロン放出線源が、イメージング装置の画像再構成領域のアイソセンタにある例を示す図である。
【図８Ｂ】図８Ｂは、異なる検出器対に対する飛行時間ヒストグラムを示す図である。
【図９Ａ】図９Ａは、同時計数線および対応する未校正タイミングスペクトル（対タイミング差ヒストグラム）の例を図示した図である。
【図９Ｂ】図９Ｂは、対消滅ターゲットの半径および同時計数対の画像再構成領域の半径を図示した図である。
【図１０Ａ】図１０Ａは、２つの起こり得る同時計数線を図示した図である。
【図１０Ｂ】図１０Ｂは、同時計数線に対応するタイミングスペクトル（対タイミング差ヒストグラム）を図示した図である。
【図１１】図１１は、一実施形態による方法を示したフローチャートを示す図である。
【図１２】図１２は、結晶毎のタイミングオフセットを計算するステップをより詳細に見て示したフローチャートを示す図である。
【図１３】図１３は、結晶毎のウォーク補正係数を計算するステップをより詳細に見て示したフローチャートを示す図である。
【図１４】図１４は、一実施形態によるイメージングシステムに対応するブロック図を示した図である。
【図１５】図１５は、一実施形態による計算装置を示すブロック図を示した図である。
【発明を実施するための形態】
【０００９】
本明細書に記載の実施形態およびこれらに付随する多くの利点をより完全に理解することは、以下の詳細な説明を参照し添付図面と関連付けて考えれば、容易にできる。
【００１０】
本実施形態は、概してポジトロン放出同位元素およびポジトロン対消滅ターゲット（「ターゲット」は「目標」とも表記する場合がある）を用いた飛行時間ポジトロン放射断層撮影（Positron  Emission  Tomography：ＰＥＴ）用のタイミング校正を行う装置および方法に関する。
【００１１】
本実施形態は、イメージングシステムを校正する装置および方法を説明する。
【００１２】
本明細書に開示する実施形態は、イメージングシステムを校正する方法を提供する。方法は、ポジトロン放出放射性同位元素線源をイメージングシステムのイメージング領域の中または近傍に配置するステップと、ポジトロン放出放射性同位元素線源から所定の距離だけ離れた位置に対消滅ターゲットを配置するステップと、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象（event）の対を検出するステップと、検出された同時計数事象の対に基づきイメージングシステムにおける検出器要素用の校正時間オフセットを計算するステップと、完成した校正時間オフセットで検出器要素を校正するステップとを含む。
【００１３】
別の実施形態の方法によると、方法は、検出された同時計数事象の対から少なくとも１つの対タイミング差ヒストグラムを生成するステップをさらに含む。
【００１４】
別の実施形態の方法によると、計算ステップは、少なくとも１つの対タイミング差ヒストグラムを用いて、検出器要素用の校正時間オフセットを計算するステップをさらに含む。
【００１５】
別の実施形態の方法によると、校正時間オフセットは、対タイミング差ヒストグラムの重心（質量中心（centroid））に対応する。
【００１６】
別の実施形態の方法によると、校正時間オフセットを計算するステップは、差分最小化技術を適用するステップをさらに含む。
【００１７】
別の実施形態の方法によると、方法は、各検出器要素に対する対タイミング差ヒストグラムの質量中心の絶対値が予め設定された閾値未満になるまで、各検出器要素に対する校正時間オフセットを繰り返し更新するステップをさらに含む。
【００１８】
別の実施形態の方法によると、ポジトロン放出放射性同位元素線源は、点状線源および線状線源のうちの１つである。
【００１９】
別の実施形態の方法によると、ポジトロン放出放射性同位元素線源は、イメージングシステムの画像再構成領域の外側にある。
【００２０】
別の実施形態の方法によると、対消滅ターゲットは、イメージングシステムの内側のボアカバーであるか、あるいはイメージングシステム内の検出器の表面である。
【００２１】
別の実施形態の方法によると、方法は、対消滅ターゲットをスキャナの画像再構成領域内に吊り下げるステップを含む。
【００２２】
別の実施形態の方法によると、対消滅ターゲットは、円筒形シェルの形状をしている。
【００２３】
別の実施形態の方法によると、方法は、ポジトロン放出放射性同位元素線源を円筒形シェル内に配置するステップをさらに含む。
【００２４】
別の実施形態によると、方法は、ポジトロン放出放射性同位元素線源と対消滅ターゲットとの間の領域を部分的に真空状態にするステップをさらに含む。このステップは、ポジトロン飛程を延ばす効果があり、これにより対消滅ターゲットにおける対消滅事象の数が増加する。
【００２５】
別の実施形態によると、方法は、対消滅ターゲットをポジトロン放出放射性同位元素線源から少なくとも５ｃｍの位置に配置するステップをさらに含む。
【００２６】
別の実施形態によると、方法は、同時計数線毎に対タイミング差ヒストグラムを少なくとも１つの領域に分割するステップをさらに含み、各領域は、対消滅ターゲットと同時計数線の１つの交点を表す。
【００２７】
別の実施形態によると、方法は、交点から同時計数線の端点を表す各検出器要素までの距離に基づき計算される異なる時間シフトを、各領域に適用するステップをさらに含む。このステップにより、対消滅の実際の位置に対し補正することができる。
【００２８】
別の実施形態の方法によると、方法は、ポジトロン放出放射性同位元素線源を円筒形シェルの内面上に実質的に均一に堆積させるステップをさらに含む。
【００２９】
別の実施形態の方法によると、方法は、ポジトロン放出放射性同位元素線源を樹脂内に封入するステップをさらに含む。
【００３０】
別の実施形態の方法によると、方法は、ポジトロン放出放射性同位元素線源を液体の状態で内面に塗布するステップをさらに含む。
【００３１】
別の実施形態の方法によると、方法は、液状のポジトロン放出放射性同位元素線源をゼラチンと混合するステップ、粘着タイプの材料を液状のポジトロン放出放射性同位元素線源と混合するステップ、および内面を吸収材料から作るステップの少なくとも１つをさらに含む。
【００３２】
別の実施形態の方法によると、方法は、対消滅ターゲットと任意の同時計数線の複数の交点が存在する場合、後続の交点との間隔が、対タイミング差ヒストグラムの検査から互いに異なる領域を識別することができるように十分に広くなるよう、同時計数対を形成するための画像再構成領域を配置するステップをさらに含む。対消滅ターゲットとの各交点のヒストグラムにおいて、１つの独立した「隆起部」があることになる。
【００３３】
別の実施形態の方法によると、ポジトロン放出線源はＧｅ−６８／Ｇａ−６８である。この線源は、運動エネルギーが高く半減期が比較的長いため使用される。
【００３４】
別の実施形態の方法によると、ポジトロン放出線源は、Ｎａ−２２、Ｆ−１８、ＦＤＧ、またはＲｂ−８２のうちの１つである。
【００３５】
別の実施形態の方法によると、方法は、対消滅ターゲットの位置と向きを決定するために、対消滅ターゲットをイメージングシステムで画像化するステップをさらに含む。画像化は、イメージングシステムがＰＥＴ、ＰＥＴ／ＣＴ、またはＰＥＴ／ＭＲシステムであるかに応じて、ＰＥＴ、ＣＴイメージング、またはＭＲであってもよい。
【００３６】
別の実施形態の方法によると、方法は、３次元レーザースキャナを用いて、内側のボアカバーの位置と向きおよび検出器表面の位置と向きの少なくとも１つを決定するステップをさらに含む。
【００３７】
別の実施形態の方法によると、方法は、画像から決定された位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算するステップをさらに含む。
【００３８】
別の実施形態の方法によると、方法は、３次元レーザースキャナから決定された決定位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算するステップをさらに含む。
【００３９】
別の実施形態の方法によると、方法は、検出された同時計数事象の対を用いてウォーク補正係数を計算するステップと、ウォーク補正係数を用いて検出器要素を校正するステップとをさらに含む。
【００４０】
別の実施形態の方法によると、方法は、ポジトロン放出放射性同位元素線源を実質的にイメージングシステムのアイソセンタに配置するステップと、ポジトロン放出放射性同位元素線源における、または線源ホルダーにおける対消滅から生じる同時計数事象の対を用いてウォーク補正係数を計算するステップと、ウォーク補正係数を用いて検出器要素を校正するステップとをさらに含む。
【００４１】
別の実施形態の方法によると、方法は、反対側の検出器対で検出された事象が個別に処理されるよう、ウォーク補正係数の計算を並列処理するステップをさらに含む。この考え方は、アイソセンタにおいて線源を使用することにより、反対側の検出器対からのデータだけを用いてウォーク補正係数を並列処理することが容易になるというものである。
【００４２】
別の実施形態の方法によると、円筒形シェルの直径は少なくとも１０ｃｍである。この値は、後述の決定に、かつＰＥＴシステム用の最先端のタイミング分解能、例えば、目標の２００ピコ秒に基づく。
【００４３】
本明細書に開示する実施形態は、イメージングシステムを提供する。イメージングシステムは、イメージングシステムのイメージング領域の中または近傍に配置されたポジトロン放出放射性同位元素線源と、ポジトロン放出放射性同位元素線源から所定の距離だけ離れた位置に配置された対消滅ターゲットと、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出するよう構成されたイメージング装置と、検出された同時計数事象の対に基づきイメージング装置の検出器要素用の校正時間オフセットを計算するよう構成された計算回路と、完成した校正時間オフセットで検出器要素を校正するよう構成された校正回路とを含む。
【００４４】
別の実施形態のシステムによると、計算回路は、検出された同時計数事象の対から少なくとも１つの対タイミング差ヒストグラムを生成するよう、かつ少なくとも１つの対タイミング差ヒストグラムを用いて、検出器要素用の校正時間オフセットを計算するようさらに構成される。
【００４５】
別の実施形態のシステムによると、校正時間オフセットは、対タイミング差ヒストグラムの質量中心に対応する。
【００４６】
別の実施形態のシステムによると、計算回路は、差分最小化技術を適用することにより校正時間オフセットを計算するようさらに構成される。
【００４７】
別の実施形態のシステムによると、計算回路は、対タイミング差ヒストグラムの測定された質量中心と各検出器要素に対する予測される対のタイミング質量中心との差が予め設定された閾値未満になるまで、または各検出器要素に対する対タイミング差ヒストグラムの質量中心の絶対値が予め設定された閾値未満になるまで、検出器要素毎に校正時間オフセットを繰り返し更新するようさらに構成される。
【００４８】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源は、点状線源および線状線源のうちの１つである。
【００４９】
別の実施形態のシステムによると、ポジトロン放出放射性線源は、イメージングシステムの画像再構成領域の外側にある。
【００５０】
別の実施形態のシステムによると、対消滅ターゲットは、イメージングシステムの内側のボアカバーであるか、あるいはイメージングシステム内の検出器の表面である。
【００５１】
別の実施形態のシステムによると、対消滅ターゲットは、スキャナの画像再構成領域内に吊り下げられる。
【００５２】
別の実施形態のシステムによると、対消滅ターゲットは、円筒形シェルの形状をしている。
【００５３】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源は、円筒形シェル内に配置される。
【００５４】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源と対消滅との間の領域は、部分的に真空状態にされる。
【００５５】
別の実施形態のシステムによると、対消滅ターゲットは、ポジトロン放出放射性同位元素線源から少なくとも５ｃｍの位置に置かれる。
【００５６】
別の実施形態のシステムによると、計算回路は、同時計数線毎に対タイミング差ヒストグラムを少なくとも１つの領域に分割するようさらに構成され、各領域は、対消滅ターゲットと同時計数線の１つの交点を表す。
【００５７】
別の実施形態のシステムによると、計算回路は、交点から同時計数線の端点を表す各検出器要素までの距離に基づき計算される異なる時間シフトを、各領域に適用するようさらに構成される。
【００５８】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源は、円筒形シェルの内面上に実質的に均一に堆積させられる。
【００５９】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源は、樹脂内に封入される。
【００６０】
別の実施形態のシステムによると、計算回路は、ポジトロン放出放射性同位元素線源を液体の状態で内面に塗布するようさらに構成される。
【００６１】
別の実施形態のシステムによると、液状のポジトロン放出放射性同位元素線源が、ゼラチンまたは粘着タイプの材料と混合される、あるいは内面は吸収材料から作られる。
【００６２】
別の実施形態のシステムによると、計算回路は、対消滅ターゲットと任意の同時計数線の複数の交点が存在する場合、後続の交点との間隔が、対タイミング差ヒストグラムの検査から互いに異なる領域を識別することができるように十分に広くなるよう、同時計数対を形成するための画像再構成領域を制限するようさらに構成される。
【００６３】
別の実施形態のシステムによると、ポジトロン放出線源はＧｅ−６８／Ｇａ−６８である。
【００６４】
別の実施形態のシステムによると、ポジトロン放出線源は、Ｎａ−２２、Ｆ−１８、ＦＤＧ、またはＲｂ−８２のうちの１つである。
【００６５】
別の実施形態のシステムによると、イメージング装置は、対消滅ターゲットの位置と向きを決定するために、対消滅ターゲットをイメージングシステムで画像化するようさらに構成される。
【００６６】
別の実施形態のシステムによると、計算回路は、３次元レーザースキャナを用いて、内側のボアカバーの位置と向きおよび検出器表面の位置と向きの少なくとも１つを決定するようさらに構成される。
【００６７】
別の実施形態のシステムによると、計算回路は、画像から決定された位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算するようさらに構成される。
【００６８】
別の実施形態のシステムによると、計算回路は、３次元レーザースキャナから決定された決定位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算するようさらに構成される。
【００６９】
別の実施形態のシステムによると、計算回路は、検出された同時計数事象の対を用いてウォーク補正係数を計算するようさらに構成され、かつ校正回路は、ウォーク補正係数を用いて検出器要素を校正するようさらに構成される。
【００７０】
別の実施形態のシステムによると、ポジトロン放出放射性同位元素線源は、実質的にイメージングシステムのアイソセンタに配置され、計算回路は、ポジトロン放出放射性同位元素線源における、または線源ホルダーにおける対消滅から生じる同時計数事象の対を用いてウォーク補正係数を計算するようさらに構成され、かつ校正回路は、ウォーク補正係数を用いて検出器要素を校正するようさらに構成される。
【００７１】
別の実施形態のシステムによると、計算回路は、反対側の検出器対で検出された事象が個別に処理されるように、ウォーク補正係数の計算を並行処理にするようさらに構成される。
【００７２】
別の実施形態のシステムによると、円筒形シェルの直径は少なくとも１０ｃｍである。
【００７３】
一実施形態によると、イメージングシステムは、イメージングシステムのイメージング領域の中または近傍に配置されたポジトロン放出放射性同位元素線源と、ポジトロン放出放射性同位元素線源から所定の位置に配置された対消滅ターゲットと、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出する検出器要素と、検出された同時計数事象の対に基づき前記検出器要素用の校正時間オフセットを計算する計算部と、計算した校正時間オフセットで前記検出器要素を校正する校正部とを含む。
【００７４】
別の実施形態によると、計算部は、検出された同時計数事象の対から少なくとも１つの対タイミング差ヒストグラムを生成するよう、かつ少なくとも１つの対タイミング差ヒストグラムを用いて、検出器要素用の校正時間オフセットを計算する。
【００７５】
別の実施形態によると、校正時間オフセットは、対タイミング差ヒストグラムの重心に対応する。
【００７６】
別の実施形態によると、計算部は、差分最小化技術を適用することにより校正時間オフセットを計算する。
【００７７】
別の実施形態によると、計算部は、対タイミング差ヒストグラムの測定された重心と各検出器要素に対する予測される対のタイミング重心との差が予め設定された閾値未満になるまで、または各検出器要素に対する対タイミング差ヒストグラムの重心の絶対値が予め設定された閾値未満になるまで、検出器要素毎に校正時間オフセットを繰り返し更新する。
【００７８】
別の実施形態によると、ポジトロン放出放射性同位元素線源は、点状線源および線状線源のうちの１つである。
【００７９】
別の実施形態によると、ポジトロン放出放射性線源は、イメージングシステムの画像再構成領域の外側にある。
【００８０】
別の実施形態によると、対消滅ターゲットは、イメージングシステムの内側のボアカバーであるか、あるいはイメージングシステム内の前記検出器の表面である。
【００８１】
別の実施形態によると、対消滅ターゲットは、円筒形状をしている。
【００８２】
別の実施形態によると、ポジトロン放出放射性同位元素線源は、円筒形状内に配置される。
【００８３】
別の実施形態によると、ポジトロン放出放射性同位元素線源と対消滅ターゲットとの間の領域は、部分的に真空状態にされる。
【００８４】
別の実施形態によると、対消滅ターゲットは、ポジトロン放出放射性同位元素線源から少なくとも５ｃｍ以上離れた位置に置かれる。
【００８５】
別の実施形態によると、計算部は、同時計数線毎に対タイミング差ヒストグラムを少なくとも１つの領域に分割し、各領域は、対消滅ターゲットと同時計数線の１つの交点を表す。
【００８６】
別の実施形態によると、計算部は、交点から同時計数線の端点を表す各検出器要素までの距離に基づき計算される異なる時間シフトを、各領域に適用する。
【００８７】
別の実施形態によると、ポジトロン放出放射性同位元素線源は、円筒形状の内面上に実質的に均一に堆積させられる。
【００８８】
別の実施形態によると、液状のポジトロン放出放射性同位元素線源が、ゼラチンまたは粘着タイプの材料と混合される、あるいは円筒形状の内面は吸収材料から作られる。
【００８９】
別の実施形態によると、計算部は、対消滅ターゲットと任意の同時計数線の複数の交点が存在する場合、後続の交点との間隔が、対タイミング差ヒストグラムの検査から互いに異なる領域を識別することができるように十分に広くなるよう、同時計数対を形成するための画像再構成領域を制限する。
【００９０】
別の実施形態によると、ポジトロン放出放射性同位元素線源はＧｅ−６８／Ｇａ−６８、Ｎａ−２２、Ｆ−１８、ＦＤＧ、またはＲｂ−８２のうちの１つである。
【００９１】
別の実施形態によると、対消滅ターゲットの位置と向きを決定するために、対消滅ターゲットを画像化する処理部を更に備える。
【００９２】
別の実施形態によると、計算部は、画像から決定された位置と向きを用いて、同時計数線毎に予測される対のタイミング差を計算する。
【００９３】
別の実施形態によると、計算部は、検出された同時計数事象の対を用いてウォーク補正係数を計算し、かつ校正部は、ウォーク補正係数を用いて検出器要素を校正する。
【００９４】
別の実施形態によると、ポジトロン放出放射性同位元素線源は、実質的にイメージングシステムのアイソセンタに配置され、計算部は、ポジトロン放出放射性同位元素線源における、または線源ホルダーにおける対消滅から生じる同時計数事象の対を用いてウォーク補正係数を計算し、かつ校正部は、ウォーク補正係数を用いて検出器要素を校正する。
【００９５】
別の実施形態によると、計算部は、反対側の検出器対で検出された事象が個別に処理されるように、ウォーク補正係数の計算を並行処理にする。
【００９６】
別の実施形態によると、円筒形状の直径は少なくとも１０ｃｍ以上である。
【００９７】
別の実施形態によると、イメージングシステムを校正する校正方法であって、ポジトロン放出放射性同位元素線源をイメージングシステムのイメージング領域の中または近傍に配置するステップと、ポジトロン放出放射性同位元素線源から所定の位置に対消滅ターゲットを配置するステップと、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対を検出するステップと、検出された同時計数事象の対に基づき前記イメージングシステムにおける検出器要素用の校正時間オフセットを計算するステップと、計算した校正時間オフセットで前記検出器要素を校正するステップとを含む。
【００９８】
本実施形態は、ＴＯＦ  ＰＥＴスキャナによる再構成画像内の統計雑音の効果的な低減を可能とする優れたタイミング分解能を提供し、画像品質を向上させる。同時計数タイミングオフセットは、ＴＯＦ  ＰＥＴスキャナのための優れたタイミング分解能を実現するために、個々の結晶毎に決定される。さらに、同時計数タイミングのエネルギー依存の特徴を表すウォーク補正係数は、ＴＯＦ  ＰＥＴスキャナのためのタイミング分解能を最適化するために、個々の結晶毎に校正される。
【００９９】
タイミングオフセットおよびウォーク補正係数は、個々の検出されたガンマ線の測定タイムスタンプを校正タイムスタンプに変換するために用いられる校正パラメータである。例えば、以下で表される。
ｔ校正＝ｔ測定＋ｔオフセット＋Ｗ×（Ｅ測定−Ｅ５１１ｋｅＶ）    ［式１］
ここで、
ｔ校正は校正されたタイムスタンプ、
ｔ測定は測定されたタイムスタンプ、
ｔオフセットはタイミングオフセット、
Ｗはウォーク補正係数、
Ｅ測定はガンマ線の測定された積算信号、
Ｅ５１１ｋｅＶは５１１ｋｅＶのガンマ線に対する積算信号である。
【０１００】
式１において、ウォーク補正が、測定エネルギーとタイムウォークとの間の線形依存性を示す単一の係数を用いて行われる。線形補正は一般的に良い結果が得られ、かつ関心エネルギー範囲（すなわち、ＰＥＴにおける５１１ｋｅＶ近傍）にわたって十分である一方で、より高次の多項式または他の関数形式を用いてウォーク補正を行うことができる。
【０１０１】
タイミングオフセット校正は、１つの結晶が同時計数のガンマ光子により多数の他の結晶と対になる際に取得してもよい。この観点から、同時計数のガンマ光子の広く分布する線源により、各結晶が多数の他の結晶と確実に対になる。計数要求値を下げるために、同時計数ガンマ光子の線源は、関心の同時計数線（Line  Of  Response：ＬＯＲ）に沿った方向に短い空間広がりを有してもよい。
【０１０２】
本実施形態において、同時計数ガンマ光子のそのように分布する線源は、ポジトロン放出線源および空間的に離された対消滅ターゲットを用いて作成される。線源とターゲットとの間の所望の離隔距離は、例えば、数１０ｃｍのオーダーである。対消滅ターゲットのいくつかの例は、１）放射線源の外側に置かれ、かつスキャナの画像再構成領域（Field  Of  View：ＦＯＶ）内に吊り下げられたシールド、２）ＰＥＴスキャナカバー、および３）ＰＥＴ検出器である。
【０１０３】
以下に説明するように、本実施形態を用いて同時計数ガンマ光子の分布線源を作成することで、先の実施形態より一定の利点がある。本実施形態ではポジトロンが、空間的に離れたターゲットにおいて対消滅する前に、相当な距離を一般的に空気または部分真空を通って移動するという事実により、本実施形態は先の実施形態と区別される。この文脈において空間的に離れたとは、例えば、ポジトロン放出線源から数センチメートル離れていることに相当する。本実施形態におけるポジトロン放出線源と対消滅ターゲットとの間の距離は、一般的に約５ｃｍ以上である。以前は、ポジトロン放出体は、水、アクリル（プラスチック）、アルミニウム、鋼鉄、または銅等の材料によって埋め込まれ、あるいは含まれ、ポジトロン対消滅のほぼ全てが、ポジトロン放出線源を直接取り囲む材料の中で発生していた。ポジトロンが対消滅の前に移動する固有の距離を、ポジトロン飛程と呼ぶ。ポジトロン飛程は、（放出されたポジトロンの平均運動エネルギーの違いにより）同位元素とポジトロンがその中を移動する（例えば、実効原子密度、実効原子量、および実効原子番号により示される）材料との両方に依存する。図２は、対消滅位置の分布に対応する１／１０値幅プロットを図示したものである。例として図２に、いくつかのポジトロン放出放射性核種に対する水中でのポジトロン対消滅箇所の分布を示す。これらの分布から計算される１／１０値幅（Full-Width-at-Tenth-Maximum：ＦＷＴＭ）が、ポジトロン飛程を説明する１つの方法である。また、例えば二乗平均距離を計算する等、他の技術を用いてポジトロン飛程を説明してもよい。以下の表１に、様々な材料についての、多数の異なる同位元素に対するＦＷＴＭにおけるポジトロン飛程の例を示す。
【０１０４】
【表１】
2020034568.tif 000003
【０１０５】
上述のように、本実施形態では、ポジトロン放出線源と対消滅ターゲットの距離は、約５ｃｍ以上である。これが可能である理由は、本実施形態では、ポジトロンが、放出される非常に短い距離の中で空気または部分真空に到達することができるように、ポジトロン放出線源が実装されているからである。一旦ポジトロンが空気に到達すると、ポジトロン飛程は（表１に示されているように）メートルのオーダーであり、ほとんどが空間的に離された対消滅ターゲットに到達する。ポジトロン放出線源が固体状または液体状に埋め込まれる先の実施形態では、飛程はミリメートルのオーダーである。表１に挙げた最高のエネルギーの放出体（Ｇｅ−６８／Ｇａ−６８）および最長の飛程を有する固体状または液体状の材料（水）を取り上げると、先の実施形態における放出と対消滅との間の最大距離は約５ｍｍ以下であり、これは本実施形態で通常用いられるポジトロン放出線源と対消滅ターゲットとの間の最小離隔距離の１０分の１より小さい。
【０１０６】
十分な数のポジトロンが対消滅ターゲットに確実に到達するようにするため、ポジトロンを単に少し減衰させるよう線源ホルダーを設計してもよい。そのようなシステムの一例は、極薄肉のステンレス鋼の皮下注射針（１６−ＲＷゲージまたはより薄肉のもの等）である。さらに、放射性同位元素は、比較的長いポジトロン飛程（すなわち、Ｇｅ−６８／Ｇａ−６８等の高い運動エネルギー）を有してもよい。放射性同位元素は、液体の状態で生成され、かつ（ポジトロンに対して）低減衰の材料（吸収紙等）内または別の材料の表面上に含まれてもよい。
【０１０７】
データ収集の後、対消滅位置補正をデータに適用することにより、タイミングオフセットを校正することができる。次に、対消滅位置補正およびタイミングオフセット補正をデータに適用することにより、タイミングウォークを校正することができる。
【０１０８】
タイミング校正は、放射性線源の外側のシールドにおいて発生するポジトロン対消滅より得られるデータを用いて、実行することができる。例えば、放射性線源がＦＯＶの中に置かれる場合、放射性線源とシールドの距離は、線源（ホルダー）とシールドにおけるポジトロン対消滅が区別されるように十分に大きくとる必要がある。放射性線源の外側のシールドは、密封して真空状態にしてもよく、こうすることで空気中のポジトロン対消滅を減らし、シールドに到達するポジトロンの数を増やすことができる。
【０１０９】
図３は、一実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す。ここで図面を参照すると、同じ参照番号はいくつかの図面にわたり同一または対応する部品を示しており、図３はＰＥＴスキャナにおける検出器モジュール４０の配置の例を示す。開示する実施形態の範囲内で、検出器モジュールをいくつでも使用することができる。図３に示す例では、シールド４１が検出器４０のリング内に配置される。シールド４１と検出器４０の間にスキャナカバー４２がある。図３に示す例では、シールド４１は対消滅ターゲットである。シールド内に、ポジトロン放出線源４３およびホルダー４４が置かれる。この例では、ポジトロンはホルダー４４を離脱し、シールド４１において対消滅４５する。
【０１１０】
この図では、ポジトロン対消滅は、放射性線源の外側のシールド４１の中で発生する。スキャナ中心に置いた線状線源を用いた。図３は、左にスキャナの正面図、右にスキャナの側面図を示している。上述のように、ポジトロン放出線源は、例えば、１６ＲＷゲージ（またはそれより薄肉の）ステンレス鋼皮下注射管の中に封入されたＧｅ−６８／Ｇａ−６８であってもよい。対消滅ターゲットを、（対消滅領域の小さい空間広がりのための）短い距離の後でポジトロンを停止させるために十分高い密度となるよう設計してもよいが、対消滅事象により放出される５１１ｋｅＶガンマ線の減衰を最小限にするように設計してもよい。また対消滅ターゲットは、プラスチック（数ｍｍ）または薄い（〜ｍｍ）金属シリンダ（球形や直角プリズム等の他の形状も可能）であってもよい。
【０１１１】
図４は、別の実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す。図４に、ＰＥＴにおける検出器モジュール４０およびスキャナカバー４２の配置の別の例示実施形態を示す。図４は、左にスキャナの正面図、右にスキャナの側面図を示している。ポジトロン対消滅は、放射性線源５３の外側のシールド４１において発生する。この例では、ポジトロン放出同位元素５３が、対消滅目標４１の内面に塗布されている。線源は、内面に均一に堆積させた、または表面の硬化樹脂／エポキシに含まれるＧｅ−６８／Ｇａ−６８であってもよく、あるいは表面に塗布されたフルオロデオキシグルコース（１８Ｆ）（Fluorodeoxyglucose：ＦＤＧ）であってもよい。ＦＤＧの接着力を増すために、ＦＤＧをゼラチンに混合してもよく、またはターゲットの内面を吸収材料で作ってもよい。また、その他の液体線源も可能である。
【０１１２】
図５は、別の実施形態によるポジトロン放出線源およびシールドを有するイメージングシステムを示す。図５に、ＰＥＴスキャナにおける検出器モジュール４０およびスキャナカバー４２の配置の別の例示実施形態を示す。図５は、左にスキャナの正面図、右にスキャナの側面図を示している。ポジトロン対消滅４５は、放射性線源６３の外側のシールド４２において発生する。この例では、ポジトロン放出同位元素６３が、スキャナの画像再構成領域の外側に置かれる。さらに、対消滅目標４１を密閉し、密閉されたシールド（例えば、シリンダ）から空気を（少なくとも部分的に）真空引きするための手段を設けることもできる。これにより、ポジトロン飛程が延びることになり、（空気中ではなく）目標４１において対消滅するポジトロンの数が増加する。シールドの密閉および真空引きは、他の実施形態にも同様に適用できる。
【０１１３】
図６は、一実施形態による複数のポジトロン放出線源およびシールドを有するイメージングシステムを示す。図６に、ＰＥＴスキャナにおける検出器モジュール４０およびスキャナカバー４２の配置の別の例示実施形態を示す。図６は、左にスキャナの正面図、右にスキャナの側面図を示している。ポジトロン対消滅４５は、放射性線源７３の外側のシールド４１において発生する。例えば、長いポジトロン飛程を有するあらゆるポジトロン放出線源を用いることができる。
【０１１４】
またタイミング校正は、スキャナカバーまたはＰＥＴ検出器において発生するポジトロン対消滅からなるデータを用いて実行することもできる。図７Ａおよび図７Ｂは、複数の追加実施形態によるシールドおよびスキャナカバーを有していないイメージングシステムを示す。図７Ａおよび図７Ｂに、ポジトロン対消滅８１／８２がスキャナカバー８１、またはＰＥＴ検出器８２のどれかにおいて発生する、この実施形態の例を示す。十分に長いポジトロン飛程を有するポジトロン放出線源を、この例で用いることができる。図７Ａはスキャナカバーが含まれる例を示し、図７Ｂはスキャナカバーが取り除かれている例を示す。
【０１１５】
図８Ａは、ポジトロン放出線源が、イメージング装置の画像再構成領域のアイソセンタにある例を示す。図８Ｂは、異なる検出器対に対する飛行時間ヒストグラムを示す。図８Ａおよび図８Ｂに、シールドが所定の位置になく、対消滅がスキャナカバー９１において発生する別の例を示す。この例では、点状線源９２がスキャナＦＯＶのアイソセンタに置かれる。スキャナカバーにおけるポジトロン対消滅９３を、例えば、タイミングオフセット校正のために用いてもよい。線源／線源ホルダーにおけるポジトロン対消滅をタイミングウォーク校正のために用いてもよい。濃い線９４で示すように、いくつかのポジトロンが線源（ホルダー）において対消滅し、一方薄い線９５で示すように、いくつかのポジトロンがスキャナカバーにおいて対消滅する。
【０１１６】
図８Ｂは、様々な検出器対に対する飛行時間（Time  Of  Flight：ＴＯＦ）ヒストグラム（対タイミング差ヒストグラム）を示している。線源（ホルダー）におけるポジトロン対消滅は、ＴＯＦヒストグラム（対タイミング差ヒストグラム）の中央のピークとして現れる。スキャナカバーにおけるポジトロン対消滅は、異なる距離（例えば、検出器差）を有する検出器対によって検出することができ、ＴＯＦヒストグラム（対タイミング差ヒストグラム）におけるサイドローブとして現れる。例えば、検出器差が長くなることで、検出器間を結ぶ直線がスキャナカバー９１を通るようになると、そのスキャナカバー９１で対消滅が生じることから、スキャナカバー９１に対応するサイドロープがＴＯＦヒストグラムに現れる。さらに、検出器間を結ぶ直線が点状線源９２近傍を通るようになると、点状線源９２近傍で生じる対消滅に対応したピークのあるＴＯＦヒストグラムとなる。
【０１１７】
図９Ａは、同時計数線および対応する未校正タイミングスペクトル（対タイミング差ヒストグラム）の例を図示したものである。図９Ｂは、対消滅ターゲットの半径および同時計数対の画像再構成領域の半径を図示したものである。図９Ａに、同時計数線（Line  Of  Response：ＬＯＲ）１０２および対応する（未校正の）タイミングスペクトル（対タイミング差ヒストグラム）の例を示す。図９Ｂは、３つのピークが中央のＬＯＲ１０５に沿って分解されるように対消滅ターゲット（ｒ目標）の半径が選択される場合に、適当なＦＯＶ１０７が、ｒ目標に等しい交点間の最短距離を有することに基づいて規定できることを示している。例えば、対消滅ターゲットを有する交点間の距離１０６がｒ目標に等しいとすると、ｒＦＯＶ＝（√３）／２×ｒ目標である。任意の所与のＬＯＲに沿って、対消滅ターゲット（または線源／線源ホルダー）との複数の交点が、未校正のタイミング分解能により決定される距離だけ離されることが好ましい。例えば、図９Ａに示す２つのＬＯＲ１０２を考えると、この例において、３つの交点が線源／ホルダーを通過するＬＯＲに対して分解されるように対消滅ターゲットの半径が選択される場合に、図９Ｂは収集および分析のための適当なＦＯＶをいかにして規定することができるかを示している（特により狭いその他のＦＯＶも可能である）。
【０１１８】
結晶毎にタイミングオフセットを決定するためのデータの処理は、一般的に反復プロセスである。各繰り返しの間、タイミング差のヒストグラム（対タイミング差ヒストグラム）が、各結晶において検出された同時計数事象の全てに対して形成される。これは、測定された時間差と、形状が与えられた際に予測される時間差との間の差のヒストグラムである。例えば、スキャナの中央における１つの点状線源に対して、予測される時間差は、線源を通過する全てのＬＯＲに対して零である。他の線源位置に対して、予測される時間差を計算することができる（例えば、線源／ターゲットの位置／向きは周知で先験的なものにすることができ、あるいは位置／向きはイメージングにより、または３次元レーザ走査法により決定することができる）。（一般的に質量中心から）時間差のヒストグラム（対の時間差のヒストグラム）に基づき、各結晶のタイミングオフセットは各繰り返しにおいて更新され、例えば、ａ）一定の繰り返し数だけ、ｂ）各繰り返しにおいて変化が前もって設定された閾値未満になるまで、またはｃ）全てのヒストグラムの質量中心の絶対値が前もって設定された閾値未満になるまで、繰り返しが継続される。
【０１１９】
図１０Ａは、２つの起こり得る同時計数線を図示したものである。図１０Ｂは、同時計数線に対応するタイミングスペクトル（対タイミング差ヒストグラム）を図示したものである。図１０Ａに示すように、あるＬＯＲが対消滅ターゲット（または線源／線源ホルダー）との複数の交点を有している場合、異なる「予測された時間差」をタイミングスペクトルの異なる領域用に規定することができる。この例では、測定された時間差と「予測される時間差」の差を、ＬＯＲ／ターゲットの形状およびタイミングスペクトルの領域（すなわち、「補正された時間差」）に基づき事象毎に計算することができる。次に、時間差スペクトルをＬＯＲ毎に生成して合計し、任意の結晶に対する最終の時間差スペクトルを形成することができ、または、任意の結晶における事象に対する「補正された時間差」の全てをヒストグラム化することにより、最終の時間差スペクトルを直接生成することができる。
【０１２０】
図１０Ａは、２つの起こり得るＬＯＲ１１１および１１２を例示している。ＬＯＲ１１１はターゲットおよび／または線源／線源ホルダー１１５との３つの交点１１６〜１１８を有する。ＬＯＲ１１２は、ターゲットおよび／または線源／線源ホルダー１１５との２つの交点１１９および１２０を有する。図１０Ｂは対応するタイミングスペクトル（対タイミング差ヒストグラム）ＬＯＲ１１３およびＬＯＲ１１４を示し、それぞれＬＯＲ１１１および１１２に対応する。さらに、各未校正タイミングスペクトル（対タイミング差ヒストグラム）１１３および１１４は、交点にそれぞれ対応する領域１２１〜１２５を有する。例えば、タイミングスペクトル（対タイミング差ヒストグラム）１１３における領域１は、交点１１８に対応する。さらに、図１０Ｂにおける各領域は、異なる「予測されたタイミング差」を有する。
【０１２１】
図１１は、一実施形態による方法を示したフローチャートを示す。図１１に、一実施形態による結晶毎のタイミングオフセットおよびウォーク補正係数を計算する方法を示すフローチャートを示す。ステップＳ１００１において、線源および対消滅ターゲットが準備される。線源および対消滅ターゲットの準備は、図３〜７に記載の実施形態のいずれに対応していてもよい。
【０１２２】
ステップＳ１００２において、対消滅ターゲットがスキャナのＦＯＶの中に置かれる。特に、ポジトロン放出線源を、スキャナのＦＯＶの中または外側に置いてもよい。
【０１２３】
ステップＳ１００３において、対消滅ターゲットの位置が測定される。例えば、線源（ホルダー）においてポジトロン対消滅を発生させる場合は、線源の正確な位置が、既知であるまたは決定されることが好ましい。例えば、位置をイメージングから決定することができる。シールドにおいてポジトロン対消滅を発生させる場合は、シールドを線源の外側に配置することが好ましい。シールドの形状および位置は、既知であるまたは決定されることが好ましい。例えば、形状および位置は、イメージングまたは３次元レーザースキャンから決定することができる。スキャナカバーにおいてポジトロン対消滅を発生させる場合、スキャナカバーの形状および位置は既知であることが好ましい。検出器においてポジトロン対消滅を発生させる場合、スキャナカバーは取り除かれ、検出器の位置は既知であることが好ましい。このように発生させることのそれぞれは、代替の実施形態では、対消滅ターゲットの位置の推定に基づき決定してもよい。この場合、推定の精度がタイミング校正の質に影響を及ぼすことになる。
【０１２４】
複数の対消滅位置が単一のＬＯＲに沿ってありうるなら、異なる位置で発生する事象をタイミングスペクトル（対タイミング差ヒストグラム）において区別できるように、対消滅位置となりうる異なる位置の間、例えば、線源とシールドの間、または線源とスキャナカバーの間に十分な空間があることが好ましい。
【０１２５】
ステップＳ１００４において、データ収集が実行される。特に、十分な数の同時計数事象が、適切なＦＯＶで取得される。データ収集ＦＯＶは、対消滅ターゲットとほぼ同じサイズかそれよりもわずかに大きくする必要がある。しかし、対消滅ターゲットの形状により単一のＬＯＲに沿って複数の対消滅が起こり得るなら、空間的に分離可能な事象のみを取得するようにＦＯＶを制限することができる（例えば、図９Ａおよび図９Ｂに関する説明を参照）。さらに、ポジトロン対消滅事象の数は、結晶のそれぞれに対してタイミング差のヒストグラム（対タイミング差ヒストグラム）からピーク位置を校正するのに十分多いことが好ましい。
【０１２６】
ステップＳ１００１、Ｓ１００２、およびＳ１００４の結果が、生の同時計数リストモードデータ５００である。このデータは、ステップＳ１００５で使用され、事象毎に対消滅ターゲットと同時計数線との交点が計算される。ステップＳ１００６において、適切な交点が、事象毎にタイミングに基づき対消滅位置として選択される。ステップＳ１００７において、事象毎に対消滅位置補正がタイミングに適用される。これらのステップの結果が、対消滅位置に対して補正された補正同時計数リストモードデータ５０１である。
【０１２７】
このデータ５０１は、ステップＳ１００８で使用され、結晶毎のタイミングオフセット５０３が計算される。図１２は、結晶毎のタイミングオフセットを計算するステップをより詳細に見て示したフローチャートを示す。図１２に、結晶毎のタイミングオフセットの計算のより詳細な図を示す。特に、図１２は、結晶毎のタイミングオフセットを計算するための反復法を説明している。
【０１２８】
結晶毎のタイミングオフセットは、１つの結晶が反対側の多数の結晶と対になる際に、あらゆる位置で発生するポジトロン対消滅からなるデータから計算することができる。対消滅位置補正は、全ての事象のＴＯＦの差に適用される（ＬＯＲが複数の場所でターゲット／線源と交差するなら、時間差について異なる領域を規定する必要がある場合がある）。次に、結晶毎のタイミングオフセットが、結晶毎のＴＯＦヒストグラム（対タイミング差ヒストグラム）のピーク位置を見つけることにより計算される。結晶毎のタイミングオフセットに対するＴＯＦの差が計算され、シーケンスが収束するまでプロセスが繰り返される。結晶毎の最終のタイミングオフセットは、全ての繰り返しにおいて校正された結晶毎のタイミングオフセットにわたる合計である。
【０１２９】
このプロセスを図１２に示す。特に、ステップＳ２００１において、初期のタイミングヒストグラム（対タイミング差ヒストグラム）が結晶毎に生成される。ステップＳ２００２において、結晶毎のタイミングオフセットが推定される。ステップＳ２００３において、リストモードデータが、結晶毎のタイミングオフセットに基づき更新される。ステップＳ２００４において、結晶毎のタイミングヒストグラム（対タイミング差ヒストグラム）が生成される。ステップＳ２００５において、結晶毎の新たなタイミングオフセットが推定される。最後に、ステップＳ２００６において、現在の繰り返しと前の繰り返しとの間のタイミングオフセット差が閾値よりも小さいか否か、または繰り返しの最大数に達したかどうかが判断される。これらの判断結果のいずれかが肯定的な結果であれば、結晶毎のタイミングオフセット５０３が生成される。これらの判断結果のいずれも肯定的な結果でなければ、フローはステップＳ２００３に戻り、プロセスを繰り返す。
【０１３０】
一旦、結晶毎のタイミングオフセットが生成されると、フローは図１１のステップＳ１００９に進む。このステップにおいて、リストモードデータが、結晶毎の新たなタイミングオフセットに基づき更新される。これにより、対消滅位置およびタイミングオフセットに対して補正された同時計数リストモードデータ５０２が生成される。
【０１３１】
次に、このデータ５０２は、結晶毎のウォーク補正係数を計算するためにステップＳ１０１０で使用される。図１３は、結晶毎のウォーク補正係数を計算するステップをより詳細に見て示したフローチャートを示す。図１３に、結晶毎のウォーク補正係数の計算のより詳細な図を示す。特に、この図は、結晶毎のウォーク補正係数を計算するための多項式フィッティングについて説明している。
【０１３２】
結晶毎のタイミングウォーク補正係数は、対消滅位置補正およびタイミングオフセット補正を適用した後に、任意の位置で発生するポジトロン対消滅からなるデータから校正してもよい。データがスキャナのアイソセンタにおいて発生するポジトロン対消滅からなる場合、並列処理を用いてデータ解析プロセスをスピードアップしてもよい。例えば、データを多数の区分（例えば、４０の検出器モジュールからなるシステムに対して２０の検出器対）に分割してもよい。データの異なる区分に対するウォーク補正係数を、同時に校正することができる。データの区分毎に、対消滅位置補正およびタイミングオフセット補正をデータに適用してもよい。特定の結晶毎に、ＴＯＦ対エネルギー曲線を生成してもよく、この特定の結晶と反対側の任意の結晶とを結ぶＬＯＲが考慮される。特定の結晶毎に、その結晶用のＴＯＦ対エネルギー曲線に適した多項式を適用することにより、結晶毎のウォーク補正係数が計算される。
【０１３３】
このプロセスは図１３に示されている。特に、ステップＳ３００１において、結晶毎のタイミング対エネルギー曲線が計算される。ステップＳ３００２において、ウォーク補正係数が、各結晶に対してタイミング対エネルギー曲線に適した多項式を適用することにより計算される。これにより、結晶毎のウォーク補正係数５０４が生成される。
【０１３４】
図１４は、一実施形態によるイメージングシステムに対応するブロック図を示したものである。図１４に、イメージングシステムの例示実施形態を示す。この図では、イメージングシステムは、例えば、ＰＥＴスキャナ内に検出器モジュール１８０を含む。必要に応じて、スキャナカバー１８４およびシールド１８１が含まれる。これらの要素（１８０、１８１、および１８４）のそれぞれが、対消滅目標１８３として機能してもよい。一旦、対消滅ターゲットにおけるポジトロンの対消滅の結果生じる同時計数事象の対が検出されると、この情報が計算回路１８５に入力される。計算回路１８５はこの情報を用いて、検出された同時計数事象の対に基づき、イメージングシステム内の検出器要素毎に校正時間オフセットを計算する。次に、完成した校正時間オフセットが校正回路１８６に入力され、校正回路１８６は、完成した校正時間オフセットで式１に基づいて検出器要素を校正する。例えば、結晶毎のオフセット値（およびウォーク係数）が、検出器内の回路内のメモリに読み込まれる。次に、検出された事象毎に、式１を用いて、測定されたタイムスタンプが「校正された」タイムスタンプに変換される。計算回路１８５および校正回路１８６は、マイクロプロセッサを有するプログラム式コンピュータにより、またはハードウェア回路により実装されてもよい。
【０１３５】
本実施形態により、これまで見られなかった大きな利点が得られる。例えば、光パルスをＰＭＴに導入する場合と比べ、本実施形態では、追加のハードウェアが不要であり、検出器および電子システム全体が校正される。
【０１３６】
ＰＭＴに結合されたプラスチックシンチレータに埋め込まれた放射性線源を使用する場合と比べ、本実施形態では、追加の収集電子機器が不要であり、より正確なタイミング校正を実現することができる。
【０１３７】
回転する線状放射性線源を有するタイミング校正アルゴリズムを使用する場合と比べ、本実施形態では、静止線源を使用してもよいため、実装がより簡単である。つまり、電動式の線源回転装置は不要である。
【０１３８】
散乱ブロック内に点状線源を有するタイミング校正アルゴリズムを使用する場合と比べ、本実施形態により多数の利点が得られる。例えば、本実施形態では、臨床での利用と同じエネルギーウインドウをタイミング校正に使用することができ、これによりタイミング校正の精度を向上させることができる。さらに、同じデータをタイミングオフセット校正とウォーク補正係数校正の両方に使用してもよい。さらに、１つの結晶が、散乱ブロックを用いる場合よりも（高い統計データで）多くの結晶と対になるため、アルゴリズムがより早く収束する。例えば、点線源から所定の位置に配置された対消滅ターゲットにより生じる対消滅を検出し、タイミング校正を行うことができることから、点線源のみの対消滅を検出してタイミング校正を行う場合と比較して、タイミング校正に関する統計データの収集を効率的に行うことができる。これは、対消滅ターゲットにより生じる対消滅の事象を統計データとして活用でき、点線源のみの対消滅の事象を用いる場合よりも、高い統計データを収集しているとも言える。
【０１３９】
図１４に示す計算回路および校正回路等の上で説明した処理と、図１１のＳ１００３においてＰＥＴイメージングにより対消滅ターゲットの位置を測定するために用いられる画像再構成との少なくとも特定の部分を、少なくとも１つのマイクロプロセッサを有するある形態のコンピュータを用いることにより、あるいはプロセッサを用いることにより、実装または支援することができる。当業者には自明であるが、コンピュータプロセッサは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ）、または他の複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）等の個別論理ゲートとして実装することができる。ＦＰＧＡまたはＣＰＬＤの実装は、ＶＨＤＬ、Ｖｅｒｉｌｏｇまたはその他のハードウェア記述言語でコード化してもよく、またそのコードを直接ＦＰＧＡまたはＣＰＬＤ内の電子メモリの中に、または別の電子メモリとして、格納してもよい。さらに電子メモリは、ＲＯＭ、消去可能ＰＲＯＭ（Erasable  PROM：ＥＰＲＯＭ）、電気的消去可能ＰＲＯＭ（Electrically  Erasable  PROM：ＥＥＰＲＯＭ）、フラッシュメモリ等の不揮発性のものでもよい。また電子メモリは、スタティックまたはダイナミックＲＡＭなどの揮発性のものでもよい。また、ＦＰＧＡまたはＣＰＬＤと電子メモリの間の対話のみならず電子メモリを管理するためにマイクロコントローラやマイクロプロセッサ等のプロセッサを備えてもよい。
【０１４０】
あるいは、コンピュータプロセッサは、本明細書に記載の機能を実行する１組のコンピュータ可読命令を含むコンピュータプログラムであって、上記の非一時的電子メモリおよび／またはハードディスクドライブ、ＣＤ、ＤＶＤ、フラッシュドライブ、もしくは他の任意の公知の記憶媒体のうちのいずれかに格納されているプログラムを実行してもよい。さらに、上記のコンピュータ可読命令は、米国インテル社によるＸｅｎｏｎ（登録商標）プロセッサまたは米国ＡＭＤ社によるＯｐｔｅｒｏｎ（登録商標）プロセッサ等のプロセッサ、ならびにＭｉｃｒｏｓｏｆｔ  ＶＩＳＴＡ（登録商標）、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ（登録商標）、ＬＩＮＵＸ（登録商標）、Ａｐｐｌｅ  ＭＡＣ−ＯＳＸ（登録商標）、および当業者に公知の他のオペレーティングシステム等のオペレーティングシステムと一緒に実行する、ユーティリティアプリケーション、バックグラウンドデーモン、もしくはオペレーティングシステムの構成要素、またはそれらの組み合わせとして提供することができる。
【０１４１】
さらに、実施形態の特定の機能は、コンピュータを利用したシステム（図１５）を用いて実現することができる。図１５は、一実施形態による計算装置を示すブロック図を示したものである。コンピュータ１０００は、情報を通信するためのバスＢまたは他の通信機構と、情報を処理するためのバスＢに結合されたプロセッサ／ＣＰＵ１００４とを含む。またコンピュータ１０００は、プロセッサ／ＣＰＵ１００４により実行される情報および命令を格納するためのバスＢに結合された、ランダムアクセスメモリ（Random  Access  Memory：ＲＡＭ）または他の動的記憶装置（例えば、ダイナミックＲＡＭ（Dynamic  RAM：ＤＲＡＭ）、スタティックＲＡＭ（Static  RAM：ＳＲＡＭ）、およびシンクロナスＤＲＡＭ（Synchronous  DRAM：ＳＤＲＡＭ）等のメインメモリ／記憶装置１００３を含む。さらに、記憶装置１００３は、ＣＰＵ１００４による命令の実行の間に一時的変数または他の中間情報を格納するために用いてもよい。また、コンピュータ１０００は、読み出し専用メモリ（Read  Only  Memory：ＲＯＭ）または他の静的記憶装置（例えば、ＣＰＵ１００４のため静的情報および命令を格納するためのバスＢに結合された、プログラマブルＲＯＭ（Programmable  ROM：ＰＲＯＭ）、消去可能ＰＲＯＭ（Erasable  PROM：ＥＰＲＯＭ）、および電気的消去可能ＰＲＯＭ（Electrically  Erasable  PROM：ＥＥＰＲＯＭ））をさらに含んでもよい。
【０１４２】
またコンピュータ１０００は、大容量記憶装置１００２等の情報および命令を格納する１つ以上の記憶装置を制御するためにバスＢに結合されたディスク制御部と、ドライブ装置１００６（例えば、フロッピー（登録商標）ディスクドライブ、読み出し専用コンパクトディスクドライブ、読み出し／書き込みコンパクトディスクドライブ、コンパクトディスクジュークボックス、テープドライブ、およびリムーバブル光磁気ドライブ）とを含む。記憶装置は、適切なデバイスインタフェース（例えば、小型コンピュータシステムインタフェース（Small  Computer  System  Interface：ＳＣＳＩ）、統合デバイスエレクトロニクス（Integrated  Device  Electronics:ＩＤＥ）、拡張ＩＤＥ（Enhanced-IDE：Ｅ−ＩＤＥ）、直接メモリアクセス（Direct  Memory  Access：ＤＭＡ）、またはウルトラＤＭＡ）を使用して、コンピュータシステム１０００に追加されてもよい。
【０１４３】
またコンピュータシステム１０００は、専用論理デバイス（例えば、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ））、または構成可能論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複雑プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、およびフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））を含んでもよい。
【０１４４】
またコンピュータシステム１０００は、コンピュータユーザに情報を表示するためのブラウン管（Cathode  Ray  Tube：ＣＲＴ）等のディスプレイを制御する、バスＢに結合されたディスプレイ制御部を含んでもよい。コンピュータシステムは、キーボードやポインティングデバイス等の、コンピュータユーザと対話するとともにプロセッサに情報を提供するための入力装置を含む。ポインティングデバイスは、例えば、プロセッサに指示情報やコマンド選択を伝達しかつディスプレイ上のカーソルの動きを制御するためのマウス、トラックボール、またはポインティングスティックであってよい。さらに、プリンタにより、コンピュータシステムによって格納および／または生成されたデータのリスト出力を印刷してもよい。
【０１４５】
コンピュータ１０００は、記憶装置１００３等のメモリに含まれる１つ以上の命令からなる１つ以上のシーケンスを実行するＣＰＵ１００４に応じて、本発明の処理ステップの少なくとも一部を実行する。そのような命令は、大容量記憶装置１００２やリムーバルブメディア１００１等の別のコンピュータ可読媒体から記憶装置に読み込まれてもよい。記憶装置１００３に含まれる命令シーケンスを実行するために、マルチプロセッシング構成の１つ以上のプロセッサを使用してもよい。他の実施形態では、ソフトウェア命令の代わりまたはソフトウェア命令との組み合わせで、ハードワイヤード回路を使用してもよい。したがって、実施形態は、ハードウェアの回路構成およびソフトウェアのいかなる特定の組み合わせにも限定されない。
【０１４６】
上述したように、コンピュータ１０００は、本発明の教示に従ってプログラムされた命令を保持し、データ構造、表、記録、または本明細書に記載の他のデータを収容するための少なくとも１つのコンピュータ可読媒体１００１またはメモリを含む。コンピュータ可読媒体の例として、コンパクトディスク、ハードディスク、フロッピーディスク、テープ、光磁気ディスク、ＰＲＯＭ（ＥＰＲＯＭ、ＥＥＰＲＯＭ、フラッシュＥＰＲＯＭ）、ＤＲＡＭ、ＳＲＡＭ、ＳＤＲＡＭ、または任意の他の磁気媒体、コンパクトディスク（例えば、ＣＤ−ＲＯＭ）、またはコンピュータが読み出すことができる任意の他の媒体がある。
【０１４７】
コンピュータ可読媒体のいずれか１またはその組み合わせに格納されるものとして、本発明は、主処理装置１００４を制御し、本発明を実現するための１または複数のデバイスを作動させ、主処理装置１００４が人間ユーザと対話できるようにするための、ソフトウェアを含む。そのようなソフトウェアは、デバイスドライバ、オペレーティングシステム、開発ツール、およびアプリケーションソフトウェアを含んでもよいが、これらに限定されるものではない。そのようなコンピュータ可読媒体は、本発明を実現する際に実行される処理の全部または一部（処理が分散されている場合）を実行するための、本発明のコンピュータプログラム製品をさらに含む。
【０１４８】
本発明のコンピュータ可読媒体におけるコンピュータコード要素は、スクリプ卜、解釈可能なプログラム、ダイナミックリンクライブラリ（Dynamic  Link  Library：ＤＬＬ)、Ｊａｖａ（登録商標）クラス、および完全な実行可能プログラムを含むが、これらに限定されるものではなく、任意の解釈可能なまたは実行可能なコード機構であってもよい。さらに、性能、信頼性、および／またはコストを改善するために、本発明の処理の部分を分散させてもよい。
【０１４９】
本明細書で使用する用語「コンピュータ可読媒体」は、実行用のＣＰＵ１００４に命令を与えることに関与する任意の媒体を指す。コンピュータ可読媒体は、不揮発性媒体または揮発性媒体を含むが、これらに限定されるものではなく、多くの形態を取ることができる。例えば、不揮発性媒体は、大容量記憶装置１００２またはリムーバブルメディア１００１等の光ディスク、磁気ディスク、および光磁気ディスクを含む。揮発性媒体は、記憶装置１００３等のダイナミックメモリを含む。
【０１５０】
実行用のＣＰＵ１００４に対する１つ以上の命令からなる１つ以上のシーケンスを実行する際に、様々な形態のコンピュータ可読媒体を用いてもよい。例えば、命令を、最初に、リモートコンピュータの磁気ディスクに収容してもよい。バスＢに結合された入力は、データを受信し、データをバスＢに乗せることができる。バスＢは、データを記憶装置１００３へ伝達し、そこからＣＰＵ１００４は命令を取り出して実行する。記憶装置１００３が受信した命令は、必要に応じて、ＣＰＵ１００４が実行する前後どちらにでも大容量記憶装置１００２に格納してもよい。
【０１５１】
またコンピュータ１０００は、バスＢに結合された通信インタフェース１００５も含む。通信インタフェース１００４は、例えば、ローカルエリアネットワーク（Local  Area  Network：ＬＡＮ）、またはインターネット等の別の通信ネットワークに接続されるネットワークにつながる双方向データ通信を可能にする。例えば、通信インタフェース１００５は、任意のパケット交換ＬＡＮに取り付けるネットワークインタフェースカードであってもよい。他の例として、通信インタフェース１００５は、非対称デジタル加入回線（Asymmetrical  Digital  Subscriber  Line：ＡＤＳＬ）カード、総合サービスデジタル網（Integrated  Services  Digital  Network：ＩＳＤＮ）カード、または対応するタイプの通信回線へのデータ通信接続を可能にするモデムであってもよい。また無線リンクが実装されてもよい。そのような実施態様のいずれにおいても、通信インタフェース１００５は、様々な種類の情報を表すデジタルデータストリームを搬送する電気信号、電磁信号、または光信号を送受信する。
【０１５２】
ネットワークは、一般的に、１つ以上のネットワークを介して他のデータデバイスへのデータ通信を可能にする。例えば、ネットワークは、ローカルネットワーク（例えば、ＬＡＮ）を通じて、または通信ネットワークを介して通信サービスを提供するサービスプロバイダによって操作される機器を通じて、別のコンピュータに接続できるようにしてもよい。ローカルネットワークおよび通信ネットワークは、例えば、デジタルデータストリームを搬送する電気信号、電磁信号、または光信号、および関連する物理層(例えば、カテゴリ５ケーブル、同軸ケーブル、光ファイバ等)を使用する。さらに、ネットワークは、パーソナル携帯情報機器（Personal  Digital  Assistant：ＰＤＡ)、ラップトップコンピュータ、または携帯電話等の携帯機器に接続できるようにしてもよい。
【０１５３】
上記説明において、フローチャートの中のいずれのプロセス、記述またはブロックも、プロセスにおける特定の論理機能またはステップを実行するための１つ以上の実行可能な命令を含むコードのモジュール、セグメントまたは一部を表すものとして理解されるべきである。当業者には自明であるが、代替の実施態様が本発明の例示実施形態の範囲に含まれ、関係する機能によっては実質的に同時にまたは逆の順序とすることも含め、先に示したまたは述べたものから順序を変えて機能を実行してもよい。
【０１５４】
以上、説明したとおり本実施形態によれば、効果的なタイミング校正を行うことができる。
【０１５５】
ある特定の実施形態を説明したが、これらの実施形態は、単なる事例として提示したものであり、発明の範囲を限定するものではない。実際、本明細書で説明した新規の方法、装置、およびシステムは、様々な別の態様で具体化が可能である。さらには、本発明の要旨を逸脱することなく、本明細書で説明した方法、装置、およびシステムの態様において様々な省略、置換、および変更が可能である。添付の特許請求の範囲およびそれらの均等物は、本発明の範囲や要旨に収まると考えられる態様または変更を有効範囲に含むためのものである。
【図１】
2020034568.tif 000004
【図２】
2020034568.tif 000005
【図３】
2020034568.tif 000006
【図４】
2020034568.tif 000007
【図５】
2020034568.tif 000008
【図６】
2020034568.tif 000009
【図７Ａ】
2020034568.tif 000010
【図７Ｂ】
2020034568.tif 000011
【図８Ａ】
2020034568.tif 000012
【図８Ｂ】
2020034568.tif 000013
【図９Ａ】
2020034568.tif 000014
【図９Ｂ】
2020034568.tif 000015
【図１０Ａ】
2020034568.tif 000016
【図１０Ｂ】
2020034568.tif 000017
【図１１】
2020034568.tif 000018
【図１２】
2020034568.tif 000019
【図１３】
2020034568.tif 000020
【図１４】
2020034568.tif 000021
【図１５】
2020034568.tif 000022