(57)【要約】  診断画像化システム（１０）は、イベントを見るように配置されたセンサ（２２）のマトリクスを有している。各センサ（２２）は、関連するセンサ（２２）の出力アナログ値をデジタル数に変換するためのアナログ／デジタル変換器（２４）に接続されている。上記マトリクスにおける、イベントに応答して他のセンサ（２２）に対し最も高い出力値を有するようなセンサ（５０）が識別される。該高センサ（５０）に対して最も近い近隣者であるような外側センサ（５２）も識別される。外側センサ（５２）の出力は、対応する出力のビット数を低減するために種々の非線形平方根関数の使用により圧縮される。最も多い情報を伝達するセンサの出力は最も少なく圧縮される一方、最も少ない情報を伝達するものは最も多く圧縮される。各イベントのアドレスが発生され、ルックアップテーブル（４４）に記憶される。該ルックアップテーブル（４４）は、リアルタイムな位置決め反復アルゴリズムをオフラインで実行するために使用される。
2007518096.tif 000002
【特許請求の範囲】
【請求項１】
イベントを見るように配置されたセンサのマトリクスであって、これらセンサが前記イベントに応答する各出力を有するようなセンサのマトリクスと、
複数のアナログ／デジタル変換器であって、個々のアナログ／デジタル変換器に前記センサの各々が接続されて、関連するセンサの出力アナログ値をデジタル数に変換するようなアナログ／デジタル変換器と、
前記マトリクスにおける、前記イベントに応答して他のセンサに対し最も高い出力値を有するような高センサを識別する手段と、
前記マトリクスにおける、前記高センサに最も近い隣接者であるような複数の外側センサを識別する手段と、
前記外側センサの出力を、これら各出力のビット数を低減するように圧縮する手段と、
前記圧縮された出力によりアドレス指定されて、対応するイベント位置を取り出すルックアップテーブルと、
を有する診断画像化システム。
【請求項２】
請求項１に記載のシステムであって、前記取り出されたイベント位置が前記高センサの中心に対するものであるようなシステムにおいて、
前記イベント位置を前記高センサの中心の位置と組み合わせる位置調整手段、
を更に有することを特徴とするシステム。
【請求項３】
請求項１に記載のシステムにおいて、前記センサが光電子増倍チューブ（ＰＭＴ）であることを特徴とするシステム。
【請求項４】
請求項１に記載のシステムにおいて、前記圧縮する手段が、
前記センサの各出力に対して複数のビットを出力の可能な範囲に基づいて動的に割り当てる符号化方法選択手段、
を含んでいることを特徴するシステム。
【請求項５】
請求項１に記載のシステムにおいて、
前記外側センサの分数を、関連する出力を前記最も高い出力値によって除算することにより決定する分数決定手段と、
前記センサの分数の降順を最も高い分数から最も低い分数まで降順に決定する分数順序決定手段と、
を更に含むことを特徴とするシステム。
【請求項６】
請求項５に記載のシステムにであって、水平軸及び垂直軸の一方に沿う対称性が用いられるようなシステムにおいて、
前記最も高い分数が位置する象限に従って前記イベントの位置を対称軸に対して反転する反転手段、
を更に含むことを特徴とするシステム。
【請求項７】
請求項５に記載のシステムにおいて、前記圧縮する手段が、
符号化方法を、
前記ルックアップテーブルで利用可能なビット数、
前記センサの数、及び
前記最も高い分数の値、
のうちの少なくとも１つに基づいて動的に選択する符号化方法選択手段、
を含んでいることを特徴とするシステム。
【請求項８】
請求項７に記載のシステムにおいて、前記符号化方法選択手段がどの圧縮方法が使用されるかを示す少なくとも１つのビットを設定することを特徴とするシステム。
【請求項９】
請求項７に記載のシステムにおいて、前記高センサ及び前記外側センサが７ＰＭＴ群を規定することを特徴とするシステム。
【請求項１０】
請求項９に記載のシステムにおいて、前記符号化方法選択手段は、
前記最も高い分数及び２番目に高い分数が互いに隣接していない、又は
低い分数のうちの１つの値が第１の所定の値より大きい、
の一方の場合に前記イベントを破棄することを特徴とするシステム。
【請求項１１】
請求項９に記載のシステムにおいて、前記最も高い分数（F1）の値が第１の所定の値（V1）より大きく、且つ、平方根関数（B）を用いて６ビット分数（F1comp）に圧縮され、
F1comp=(int)(sqrt(k*(F1-V1))+0.5)
ここで、F1compは前記最も高い分数（F1）に適用された加重平方根圧縮であり、
k=(Codemax)2/F1maxであり、
Codemaxは前記最も高い分数（F1）に対する最大の選択されたコードであり、
F1maxは前記最も高い分数（F1）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１２】
請求項１１に記載のシステムにおいて、前記２番目に高い分数（F2）の値が第２の所定の値（V2）より大きく、且つ、平方根関数（C）を用いて５ビット分数（F2comp）に圧縮され、
F2comp=(int)(sqrt(k*(F2-V2+1))+0.5)
ここで、F2compは前記２番目に高い分数（F2）に適用された加重平方根圧縮であり、
k=(Codemax)2/F2maxであり、
Codemaxは前記２番目に高い分数（F2）に対する最大の選択されたコードであり、
F2maxは前記２番目に高い分数（F2）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１３】
請求項１１に記載のシステムにおいて、前記２番目に高い分数（F2）の値が第２の所定の値（V2）以下であり、且つ、平方根関数（D）を用いて５ビット分数（F2comp）に圧縮され、
F2comp=(int)(sqrt(k*F2)+0.5)
ここで、F2compは前記２番目に高い分数（F2）に適用された加重平方根圧縮であり、
k=(Codemax)2/F2maxであり、
Codemaxは前記２番目に高い分数（F2）に対する最大の選択されたコードであり、
F2maxは前記２番目に高い分数（F2）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１４】
請求項１１に記載のシステムにおいて、前記低い分数（F3〜F6）の値が、各々、平方根関数（E）を用いて３ビット分数（Fcomp）に圧縮され、
Fcomp=(int)(sqrt(k*F)+0.5)
ここで、Fcompは対応する低い分数（F3〜F6）に適用された加重平方根圧縮であり、
k=(Codemax)2/F2maxであり、
Codemaxは対応する低い分数（F3〜F6）に対する最大の選択されたコードであり、
Fmaxは対応する低い分数（F3〜F6）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１５】
請求項１１に記載のシステムにおいて、前記圧縮する手段が、前記２番目に高い分数（F2）の前記最も高い分数（F1）に対する方向的関係を示すために少なくとも１つのビットを設定することを特徴とするシステム。
【請求項１６】
請求項９に記載のシステムにおいて、前記最も高い分数（F1）の値が第１の所定値（V1）以下であると共に第２の所定値（V2）より大きく、前記分数（F1〜F6）の値が、各々、平方根関数（G）を使用することにより４ビットの分数（Fcomp）に圧縮され、
Fcomp=(int)(sqrt(k*F)+0.5)
ここで、Fcompは対応する分数（F1〜F6）に適用された加重平方根圧縮であり、
k=(Codemax)2/Fmaxであり、
Codemaxは対応する分数（F1〜F6）に対する最大の選択されたコードであり、
Fmaxは対応する分数（F1〜F6）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１７】
請求項９に記載のシステムにおいて、前記最も高い分数（F1）の値が第２の所定値（V2）以下であると共に第３の所定値（V3）より大きく、前記分数（F1〜F6）の値が、各々、平方根関数（H）を使用することにより４ビットの分数（Fcomp）に圧縮され、
Fcomp=(int)(sqrt(k*F)+0.5)
ここで、Fcompは対応するチューブ分数（F1〜F6）に適用された加重平方根圧縮であり、
k=(Codemax)2/Fmaxであり、
Codemaxは対応する分数（F1〜F6）に対する最大の選択されたコードであり、
Fmaxは対応する分数（F1〜F6）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１８】
請求項９に記載のシステムにおいて、前記最も高い分数（F1）の値が第３の所定値（V3）以下であり、前記分数（F1〜F6）の値が、各々、平方根関数（J）を使用することにより４ビットの分数（Fcomp）に圧縮され、
Fcomp=(int)(sqrt(k*F)+0.5)
ここで、Fcompは対応する分数（F1〜F6）に適用された加重平方根圧縮であり、
k=(Codemax)2/Fmaxであり、
Codemaxは対応する分数（F1〜F6）に対する最大の選択されたコードであり、
Fmaxは対応する分数（F1〜F6）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項１９】
請求項８に記載のシステムにおいて、前記高センサ及び前記外側センサが、４−ＰＭＴ群を規定することを特徴とするシステム。
【請求項２０】
請求項１９に記載のシステムにおいて、前記分数（F1〜F3）の値が、各々、平方根関数（A）を使用することにより８ビットの分数（Fcomp）に圧縮され、
Fcomp=(int)(sqrt(k*F)+0.5)
ここで、Fcompは対応するチューブ分数（F1〜F3）に適用された加重平方根圧縮であり、
k=(Codemax)2/Fmaxであり、
Codemaxは対応する分数（F1〜F3）に対する最大の選択されたコードであり、
Fmaxは対応する分数（F1〜F3）に対する最大の選択された値である、
ことを特徴とするシステム。
【請求項２１】
ガンマカメラのセンサのマトリクスを用いて放射を検出するステップであって、前記センサが、イベントを見るように配置されると共に前記イベントに応答する各出力を有するようステップと、
関連するセンサの出力アナログ値をアナログ／デジタル変換器によりデジタル数に変換するステップであって、前記センサの各々が個々のアナログ／デジタル変換器に接続されているようなステップと、
前記マトリクスにおける、前記イベントに応答して他のセンサに対し最も高い出力値を有するような高センサを識別するステップと、
前記マトリクスにおける、前記高センサに最も近い隣接者であるような複数の外側センサを識別するステップと、
前記外側センサの出力を、これら各出力のビット数を低減するように圧縮するステップと、
前記圧縮された出力によりルックアップテーブルをアドレス指定して、対応するイベント位置を取り出すステップと、
を有する診断画像化方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、診断画像化システム及び方法に関する。本発明は、ポジトロン放出トモグラフィ（ＰＥＴ）及び単光子放出トモグラフィ（ＳＰＥＣＴ）に関連して特別な用途を有し、特に、これらを参照して説明されるであろう。本発明が他の同様の用途及び診断画像化モードにも適用可能であることが分かるであろう。
【背景技術】
【０００２】
核画像化は、患者の解剖学的構造を画像化するために、放射能源を使用する。典型的には、放射性医薬品が患者に注入される。放射性医薬化合物は、予測可能なレートのガンマ線崩壊及び特性エネルギを持つ放射性アイソトープを含んでいる。放出された放射を監視及び記録するために、患者に隣接して放射検出器が配置される。しばしば、該検出器は、複数の方向からの放出放射を監視するために患者の周りで回転され又は割り出し回転（index）される。検出された位置及びエネルギ等の情報に基づいて、人体内の放射性医薬品の分布が決定され、循環系、及び選択された器官又は組織における放射性医薬の取り込み等を調べるために該分布の画像が再構築される。
【０００３】
伝統的なシンチレーション検出器においては、検出器は大きなシンチレーション結晶又は小さなシンチレーション結晶のマトリクスからなるシンチレータを有している。何れの場合においても、上記シンチレータはセンサのマトリクスにより見られる。普通に使用されるセンサは、光電子増倍チューブ（ＰＭＴ）である。放射吸収材料の格子状又はハニカム状アレイを含むコリメータを、上記シンチレータと被検体との間に配置して、当該シンチレータに衝突する放射の受入角を制限する。上記シンチレータに衝突する各放射イベントは、前記ＰＭＴにより見られるような対応する閃光（シンチレーション）を発生する。ＰＭＴの稠密六方アレイにおいては、斯かるイベントは主に最も近いＰＭＴ及び該最も近いＰＭＴを囲む６つのＰＭＴにより見られる。当該閃光元に対する個々のＰＭＴの近さは、当該ＰＭＴにより光りが見られる程度に影響を与える。イベントを見る各ＰＭＴは、対応する電気パルスを発生する。これら電気パルスの各振幅は、当該閃光からの各ＰＭＴの距離に概ね比例する。斯かるＰＭＴからの出力に基づいて、ガンマカメラが放射イベントをマッピングする。即ち、該ガンマカメラは上記シンチレータに衝突する放射線のエネルギ及び位置を決定する。
【０００４】
典型的には、イベント位置はイベント位置決めのための従来のアンガ（Anger）法を使用して決定され、該方法はイベントの発生後にＰＭＴにより出力される信号を加算及び加重する。イベント位置決めのためのアンガ法は、簡単な初期時点計算（first moment calculation）に基づくものである。エネルギは典型的には全てのＰＭＴ信号の和として測定され、位置は典型的には全てのＰＭＴ信号の“質量の中心”即ち重心として測定される。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、現在のリアルタイム位置決めアルゴリズムは、既知の反復方法のようには正確ではない。それでも、斯かる反復方法はリアルタイムな画像発生及び処理をサポートするには計算的に厳しい。かくして、イベントの一層正確な位置を決定するために反復方法を使用するのを可能にするような方法及び装置に対する需要がある。本発明は、上述した課題及びその他を克服するような新たな改善された画像化装置及び方法を提供するものである。
【課題を解決するための手段】
【０００６】
本発明の一態様によれば、診断画像化システムが開示される。イベントを見るためにセンサのマトリクスが配置され、これらセンサは当該イベントに応答するような出力を各々有する。これらセンサの各々は、対応するセンサの出力アナログ値をデジタル数に変換する個々のアナログ／デジタル変換器に接続される。或る手段が、上記マトリクスにおける、当該イベントに応答して他のセンサに対し最高の出力値を持つような高いセンサを識別する。或る手段が、当該マトリクスにおける、上記高いセンサの最も近い隣接者であるような複数の外側センサを識別する。或る手段が、上記外側センサの出力を圧縮して、各出力のビット数を低減する。斯かる圧縮された出力によりルックアップテーブルがアドレス指定され、対応するイベント位置を取り出す。
【０００７】
本発明の他の態様によれば、診断画像化のための方法が開示される。放射はガンマカメラのセンサのマトリクスにより検出され、これらセンサはイベントを見るように配置されると共に、斯かるイベントに応答するような対応する出力を有する。関連するセンサの出力アナログ値はアナログ／デジタル変換器によりデジタル数に変換され、これらセンサの各々は個々のアナログ／デジタル変換器に接続される。上記マトリクスにおける、当該イベントに応答して他のセンサに対し最高の出力を有するような高いセンサが識別される。当該マトリクスにおける、該高いセンサに最も近い隣接者であるような複数の外側センサが識別される。これら外側センサの出力は、対応する出力のビット数を低減するために圧縮される。斯かる圧縮された出力によりルックアップテーブルがアドレス指定され、対応するイベント位置を取り出す。
【０００８】
本発明の一つの利点は、反復アルゴリズムの使用による一層正確なイベント位置の推定にある。
【０００９】
他の利点は、適切に位置決めすることができないようなイベントのためにイベントの区別／推定を実行することにある。
【００１０】
本発明の更に他の利点及び利益は、当業者にとり、好ましい実施例の下記の詳細な説明を精読及び理解することにより明らかとなるであろう。
【００１１】
本発明は、種々の構成要素及び構成要素の配置、並びに種々のステップ及びステップの配列の形態をとり得る。また、図面は好ましい実施例を図示する目的のためのみのものであって、本発明を限定するものと見なしてはならない。
【発明を実施するための最良の形態】
【００１２】
図１を参照すると、ＰＥＴスキャナ１０は複数の検出器ヘッド、即ちガンマカメラ又は検出器１２を含み、これら検出器は検査領域１６内に配置された対象１４（好ましくは、放射性核種分布を含む）に対面すると共に、好ましくは該対象の周りで移動するように取り付けられる。検出器１２の各々はシンチレータ２０を含み、該シンチレータは放射イベント（例えば、シンチレータ２０に衝突する上記放射性核種の分布からの放射線）を閃光即ちシンチレーションに変換する。光電子増倍チューブ（ＰＭＴ）２２等のセンサのマトリクスが、シンチレータ２０からの上記閃光を見る又は受信するように配置されている。好ましい実施例においては、上記センサ２２のマトリクスはＰＭＴの稠密六方配列である。しかしながら、他のセンサ及び詰め込み配列も考えられる。一般的に、放射は、放射性アイソトープの崩壊において生じるようなガンマ量子を生成する。斯かる崩壊量子は、好ましくはドーピングされたヨウ化ナトリウム（ＮａＩ）を含むシンチレータ２０を打撃し、シンチレーションを生じる。該シンチレーションからの光は複数のセンサ２２上に分散される。所与のセンサ２２が見る又は受信する特定のシンチレーションからの光の量は、当該イベントからのセンサの距離に伴って累進的に減少する傾向にある。センサ２２の各々は、受信された閃光に応答して対応する出力信号（例えば、アナログ電気パルス）を発生するが、該出力信号は受信された閃光のものと比例する。好ましくは、センサ２２の各々はアナログ／デジタル（Ａ／Ｄ）変換器２４に電気的に接続され、該変換器は各アナログ出力をデジタル信号に変換する。当該情報はプロセッサ２８に伝達され、該プロセッサは各検出器ヘッドに対して発生する各シンチレーションイベントの位置（ｘ，ｙ）及び／又はエネルギ（ｚ）を測定又は決定する。所与のイベントに関して、該イベントに応答したセンサの出力信号の積分（即ち、センサ２２の出力パルスに関し時間に対して振幅又は強度をプロットする曲線の下の面積又はその幾らかの部分を見つける処理）、該イベントに応答するセンサ出力信号のピーク振幅、又は観測される光の量、従って該イベントのエネルギに比例する若しくは適切に関係する何らかの他の尺度により、各センサ２２の出力値がオプションとして決定される。イベントのシンチレータ２０上の位置は、二次元（２Ｄ）デカルト座標系において公称のｘ及びｙ座標により分解及び／又は決定される。しかしながら、他の座標系も考えられる。
【００１３】
図１を続けて参照すると、スキャナ１０は所望に応じてＳＰＥＣＴモード又はＰＥＴモードの何れかにおいて選択的に動作することができる。ＳＰＥＣＴモードでは、カメラ１２にコリメータが取り付けられ（図示略）、該コリメータは放射の受入を特定の方向に、即ち既知の光線に沿って制限する。このように、放射が検出される当該シンチレータ２０上の決定された位置及びカメラ１２の角度位置が、それに沿って各放射イベントが発生したような光線を定める。これらの光線の軌道及びカメラ角度位置（例えば、角度位置レゾルバ３０により得られる）は、再生プロセッサ３２に伝達され、該プロセッサは上記光線を体積的画像表現に逆投影又は再生し、該表現が画像メモリ３４に記憶される。
【００１４】
ＰＥＴモードでは、前記コリメータは除去される。このように、単一のシンチレーションイベントの位置は光線を規定しない。しかしながら、ＰＥＴ走査に使用される放射性核種は消滅イベントの対象となり、該イベントにおいて放射の２つの光子は直径方向的に反対方向に、即ち１８０度離れて同時に放出される。一致検出器３６は、２つのカメラ１２上のシンチレーションが同時に発生する場合を検出する。２つの同時的シンチレーションの位置は、消滅イベントによる光線の端点を定める。光線又は軌道計算器３８は、同時に受信されたシンチレーションイベントの各対から、対象１４を経る対応する光線を計算する。光線計算器３８からの上記光線軌道は、体積的画像表現に再生するために再生プロセッサ３２に伝達される。結果としての画像表現は体積画像メモリ３４に記憶される。
【００１５】
ＳＰＥＣＴモード及びＰＥＴモードの両モードにおいて、ビデオプロセッサ４０はモニタ４２上での表示のために、上記画像表現データを処理及び／又はフォーマットする。
【００１６】
図１を続けて参照すると共に、図２を更に参照すると、プロセッサ２８は、反復的アルゴリズムを使用してリアルタイムなイベント位置決めを実行するために発生されたルックアップテーブル４４を含んでいる。典型的には、画像化に先行して、ＰＭＴ２２は既知のアイソトープの１つにより校正される。全てのＰＭＴと同様に、高いチューブ５０の中心の位置は、先験的情報（例えば、物理的測定）から正確に分かる。各チューブ２２に対する校正情報はルックアップテーブル４４に記憶される。ルックアップテーブル４４は、高いチューブの中心に対する或るイベントの位置を識別するために使用される。更に詳細には、プロセッサ２８はＰＭＴ値決定手段４６を含み、該手段はＰＭＴ２２の初期値を決定する。最高出力決定手段又はアルゴリズム４８は、全ての他のＰＭＴ２２に対して最高の出力値を持つような高いチューブ５０と、該高いチューブ５０に隣接する複数の最も近い近隣チューブ又は外側チューブ５２とを（例えば、図５に見られるようなＰＭＴ部分集合又はＰＭＴ群５４）識別する。好ましくは、最も近い近隣チューブ５２の最大数は６とし、中心のセンサが最高の出力値を持つものとなるような７ＰＭＴの稠密六方群（hexagonal close-packed cluster）を画定する。しかしながら、上記群５４は、例えば３又は５等の異なる数のチューブを含んでも良い。分数決定手段５６は、群５４のチューブ５０、５２の出力を高いチューブ５０の値により除算し、従って斯かる変数の１つ、例えば高いチューブ５０の値を除去する。該演算の結果、群５４内の外側チューブ５２の数に依存して、３〜６個の１０ビット分数Ｆ１〜Ｆ６が得られる。分数順序決定手段５８は、チューブ５２の分数Ｆ１〜Ｆ６の降順（例えば最高の分数Ｆ１、２番目に高い分数、低い分数Ｆ３〜Ｆ６）を識別する。分数Ｆ１〜Ｆ６はＰＭＴ値メモリ６０に記憶される。分数Ｆ１〜Ｆ６は、上記高いチューブに対して当該イベントの位置を識別すると共に、ルックアップテーブル４４をアドレス指定するために使用することができる。
【００１７】
典型的には、各アナログ／デジタル変換器２４は１０ビットの強度値を出力する。結果として、圧縮なしでは、ルックアップテーブル４４は１０２１のメモリアドレスを有するであろう。圧縮手段７０は圧縮アルゴリズムを動的に適用して、ルックアップテーブル４４をアドレス指定するためのビット数を低減する。例えば、１０ビット分数Ｆ１〜Ｆ６は、以下に詳述する種々の規準に基づいて、ルックアップテーブル４４内の利用可能な空間に圧縮される。即ち、好ましくは、低い分数を高い分数よりも一層大きな程度に圧縮するような動的圧縮が適用される。更に、最高のチューブ５０に対して次に高いチューブ（又は複数のチューブ）のイベント空間位置を符号化する必要性を除去するために対称性が使用される。勿論、一層大きな（例えば１０２１の）ルックアップテーブルの作成が可能になった場合は、当該圧縮ステップを削除することができるようになることが考えられる。
【００１８】
図２を続けて参照すると、符号化選択手段７２は、各イベントに対して、ルックアップテーブル４４に関して利用可能なアドレスビット数、群５４におけるチューブ５０、５２の数、最高の分数Ｆ１の大きさ及び他の評価規準に基づいて符号化方法を動的に選択する。最も多い情報を伝達するチューブ５２は最も少なく圧縮される一方、最も少ない情報を伝達するものは最も多く圧縮される。好ましくは、各分数を符号化するのに要するビット数を低減するために、チューブ分数Ｆ１〜Ｆ６の非線形圧縮（例えば、加重平方根関数）が用いられる。量子化エラーを最小化するために、チューブ分数範囲に基づいて、異なる非線形圧縮アルゴリズムが使用される。好ましくは、量子化エラーを最小化するために値が圧縮されねばならないダイナミックレンジを減少させるために、信号のバイアス処理が実行された場合（例えば、一定数を最高の分数から減算することができる）を示すためにアドレスビットが使用される。
【００１９】
図２を続けて参照すると、アドレス発生手段７４はイベントのアドレスを発生する。結果としてのイベントアドレスは、ルックアップテーブル４４にアクセスして、前記最高のチューブ中心に対するイベント位置を取り出す。位置決定手段８０は当該イベントのｘ，ｙ位置を含む１６ビットアドレスを取り出し、これを別個の８ビットのｘ及びｙ値に分解する。該高いチューブ位置の中心のｘ，ｙ位置は上記ルックアップテーブルの出力を用いて位置調整手段８２により調整され、当該イベントの最終位置を導出する。
【００２０】
図５の実施例において、高いチューブ５０は４チューブ群内に位置する。即ち、該高いチューブ５０は３つの近隣チューブ、即ち外側チューブ５２を有している。水平軸に沿う水平対称性が、事例（ケース）の数を２なる係数により減少させるために使用される。４チューブ群５４は、外側チューブ５２が当該検出器の上に配置されたかのように解釈される。好ましくは、最高の分数Ｆ１が常に上の象限に存在することを示すためにフラグが設定される。
【００２１】
再び図２を参照して、反転手段８４は、最高の分数Ｆ１が位置する象限に従って、上記ルックアップテーブルからのイベント位置を対称軸に対して反転する。例えば、ｘ位置の調整の符号はｙ軸に対する最高分数Ｆ１の位置に基づいて選択される一方、ｙ補正の符号はｘ軸に対する位置に基づいて調整される。図示した４群の例では、位置決定手段８０は、符号化されたアドレスにおいてルックアップを実行し、高いチューブ５０の中心に対するｘ，ｙ位置を取り出す。反転が実行されない場合、位置調整手段８２は最終ｘ，ｙ位置を、高いチューブ５０の中心位置に相対ｘ，ｙ位置を加えたものに等しいとして計算する。反転が実行される場合、位置調整手段８２は最終ｘ，ｙ位置を、高いチューブ５０の中心位置から相対ｘ，ｙ位置を引いたものに等しいとして計算する。勿論、垂直対称性が適宜使用されることも考えられる。
【００２２】
図２及び３を参照すると共に、更に図４を参照すると、ステップ９０において、符号化方法選択手段７２は符号化方法６を選択してルックアップテーブル４４を動的に選択する。好ましくは、例えばビット２５〜２７等の３つのアドレスビットが、チューブ分数F1〜F3がどの様に符号化及び圧縮されているかを示すために使用される。符号化方法６において、各分数F1、F2、F3は、
【数１】
2007518096.tif 000003
なる平方根関数（Ａ）を用いて８ビット（２５６コード）に変換される、ここで、Ｆcompは対応するチューブ分数F1〜F6に適用された加重平方根圧縮であり、ｋはチューブ関数を圧縮するために使用される一定スカラ、ｋ＝（Codemax）２／Fmaxであり、Codemaxは２bits−１＝２５５に等しい、各分数F1〜F6に対する最大コードであり、Fmaxは各分数F1〜F6に対する最大値＝１０２３であり、ｋ＝(Codemax)２／Fmax＝２５５２／１０２３＝63.56304985である。
【００２３】
符号化方法６により符号化された各分数F1〜F6の復号即ち伸張は、
【数２】
2007518096.tif 000004
なる関数を用いることにより達成される。
【００２４】
図７の実施例において、高いチューブ５０は７つのチューブの群５４内に配置されている。即ち、高いチューブ５０は６つの隣接する、即ち外側のチューブ５２を有している。図５の実施例と同様に、水平対称が、３つの最高分数位置のみが必要とされるように使用される。反転手段８４は、最高分数F1が下側チューブのうちの１つであるなら、チューブ信号を反転する。
【００２５】
図３を再び参照すると、ステップ９２において、郡５４が７つのチューブを有し、最高の分数F1が２５５より大きい場合、例えば当該イベントが高いチューブ５０の中心から離れて発生し、隣接するチューブ５２の１以上に近づき、ＰＭＴの二重点ｄ又は三重点ｔを示す場合、符号化方法選択手段７２は符号化方法０、１又は２を選択する。ステップ９４において、最高の分数Ｆ１及び２番目に大きい分数Ｆ２が互いに隣接していないと判定された場合、符号化方法選択手段７２は当該イベントを破棄する。ステップ９６において、低い分数Ｆ３〜Ｆ６のうちの１つが２５５より大きいと判定された場合（これは、当該イベントが汚染されていることを示す）、符号化方法選択手段７２は該イベントを破棄する。
【００２６】
図６を参照すると、ビット２５〜２７は、どの様にチューブ分数が符号化及び圧縮されているかを示すために使用される。該符号化ＩＤは、最高の分数Ｆ１を有するチューブの位置を示すような図７のチューブ図を用いて設定される。殆どの情報は最高の及び２番目に高い分数Ｆ１、Ｆ２に含まれているが、好ましくは、他のアドレスビット（例えばビット２４）も最高の分数Ｆ１に対する２番目に高い分数Ｆ２の方向的関係（時計方向又は反時計方向）を示すために使用される。これは、最高の及び２番目に高い分数Ｆ１、Ｆ２に対して一層低い分数Ｆ３〜Ｆ６に対してよりも多くのアドレスビットが割り当てられるのを可能にする。好ましくは、ビット２４は、通常、ＣＣＷ（例えば、２番目に高い分数Ｆ２は最高の分数Ｆ１に対して反時計方向である）に設定される。ビット２４は、２番目に高い分数Ｆ２が最高の分数Ｆ１に対して時計方向である場合にはクリアされる。
【００２７】
続けて図６を参照すると、ビット１７〜２２は最高の分数Ｆ１に対して指定される。圧縮は、最高の分数Ｆ１から２５５（０ｘＦＦ）を減算すると共に該結果を下記平方根関数（Ｂ）（入力範囲は０〜７６８［０ｘ３００］）を用いて６ビットに変換することにより実行される。
【数３】
2007518096.tif 000005
ここで、Ｆ１compは対応するチューブ分数F1に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆ１ｍａｘ、であり、Codemaxは２bits−１＝６３に等しいような分数F1に対する最大コードであり、Ｆ１maxは７６８に等しいような分数Ｆ１に対する最大値であり、ｋ＝６３２／７６８＝5.16796875である。
【００２８】
符号化関数（Ｂ）により符号化された分数Ｆ１の復号即ち伸張は、
【数４】
2007518096.tif 000006
なる関数を用いることにより達成される。
【００２９】
図６を続けて参照すると、ビット１２〜１６は２番目に高い分数Ｆ２に対して指定されている。符号化方法選択手段７２は、上記２番目に高い分数Ｆ２の値に基づいて符号化方法を選択する。上記２番目に高い分数Ｆ２が２００（０ｘＣ８）より大きい場合、圧縮は、該２番目に高い分数Ｆ２の値から２０１（０ｘＣ９）を減算すると共に、結果の値を下記平方根関数（Ｃ）（入力範囲は、０〜８２２［０ｘ３３６］）を用いて５ビット（３２コード）に圧縮することにより実行される。ビット２３は、減算が実行されたことを示すために１に設定される。
【数５】
2007518096.tif 000007
ここで、Ｆ２compは上記２番目に高い分数Ｆ２に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆ２ｍａｘ、であり、Codemaxは２bits−１＝３１に等しいような分数Ｆ２に対する最大コードであり、Ｆ２maxは８２２に等しいような分数Ｆ２に対する最大値であり、ｋ＝３１２／８２２＝1.169099757である。
【００３０】
符号化関数（Ｃ）により符号化された分数Ｆ２の復号即ち伸張は、
【数６】
2007518096.tif 000008
なる関数を用いることにより達成される。
【００３１】
上記２番目に高い分数Ｆ２が２００（０ｘＣ８）以下である場合、当該圧縮は、該２番目に高い分数の値を下記平方根関数（Ｄ）（入力範囲は、０〜２００［０ｘＣ８］）を用いて５ビット（３２コード）に圧縮すると共にビット２３をクリアすることにより実行される。
【数７】
2007518096.tif 000009
ここで、Ｆ２compはチューブ分数Ｆ２に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆ２ｍａｘ、であり、Codemaxは２bits−１＝３１に等しいような分数Ｆ２に対する最大コードであり、Ｆ２maxは２００に等しいような分数Ｆ２に対する最大値であり、ｋ＝３１２／２００＝4.805である。
【００３２】
符号化関数（Ｄ）により符号化された分数Ｆ２の復号即ち伸張は、
【数８】
2007518096.tif 000010
なる関数を用いることにより達成される。
【００３３】
図６を続いて参照すると、ビット０〜１１は低い分数Ｆ３〜Ｆ６、例えば３番目ないし６番目の分数に対して指定されている。低い分数Ｆ３〜Ｆ６は、各々、平方根関数Ｅを用いて３ビット（８つのコード）に変換される。分数Ｆ３〜Ｆ６の最大値は２５５（０ｘＦＦ）に等しいと仮定される。即ち、
【数９】
2007518096.tif 000011
ここで、Ｆcompは対応するチューブ分数Ｆ３〜Ｆ６に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆｍａｘ、であり、Codemaxは２bits−１＝７に等しいような低い分数Ｆ３〜Ｆ６に対する最大コードであり、Ｆmaxは２５５に等しいような低い分数Ｆ３〜Ｆ６に対する最大値であり、ｋ＝７２／２５５＝0.192156863である。
【００３４】
符号化関数（Ｅ）により符号化された分数Ｆ３〜Ｆ６の復号即ち伸張は、
【数１０】
2007518096.tif 000012
なる関数を用いることにより達成される。
【００３５】
再び図３を参照すると共に、図８〜９を更に参照すると、ステップ９８において、符号化方法選択手段７２は、群５４が７つのチューブ５０、５２を有し、最高の分数Ｆ１が２００より大きく２５５以下である場合、符号化方法３を選択する。分数Ｆ１〜Ｆ６は、各々、下記平方根関数Ｇを用いて４ビット（１６コード）に変換される。尚、最大分数Ｆ１は２５５（０ｘＦＦ）に等しいと仮定する。
【数１１】
2007518096.tif 000013
ここで、Ｆcompは対応するチューブ分数Ｆ１〜Ｆ６に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆｍａｘ、であり、Codemaxは２bits−１＝１５に等しいような分数Ｆ１〜Ｆ６の各々に対する最大コードであり、Ｆmaxは２５５に等しいような各分数Ｆ１〜Ｆ６に対する最大値であり、ｋ＝１５２／２５５＝0.882352941である。
【００３６】
符号化関数（Ｇ）により符号化された分数Ｆ１〜Ｆ６の復号即ち伸張は、
【数１２】
2007518096.tif 000014
なる関数を用いることにより達成される。
【００３７】
図３及び図８〜９を続いて参照すると、ステップ１００において、符号化方法選択手段７２は、群５４が７つのチューブ５０、５２を有し、最高の分数Ｆ１が１６０より大きく２００以下である（例えば、イベントが高いチューブ５０の中心近傍で発生した）場合、符号化方法４を選択する。分数Ｆ１〜Ｆ６は、下記平方根関数Ｈを用いて４ビット（１６コード）に変換される。最大の分数Ｆ１は２００（０ｘＣ８）に等しいと仮定する。
【数１３】
2007518096.tif 000015
ここで、Ｆcompは対応するチューブ分数Ｆ１〜Ｆ６に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆｍａｘ、であり、Codemaxは２bits−１＝１５に等しいような各分数Ｆ１〜Ｆ６に対する最大コードであり、Ｆmaxは２００に等しいような各分数Ｆ１〜Ｆ６に対する最大値であり、ｋ＝１５２／２００＝1.125である。
【００３８】
符号化関数（Ｈ）により符号化された分数Ｆ１〜Ｆ６の復号即ち伸張は、
【数１４】
2007518096.tif 000016
なる関数を用いることにより達成される。
【００３９】
図３及び図８〜９を続いて参照すると、ステップ１０２において、符号化方法選択手段７２は、群５４が７つのチューブ５０、５２を有し、最高の分数Ｆ１が１６０以下である場合、符号化方法５を選択する。分数Ｆ１〜Ｆ６は、各々、下記平方根関数Ｊを用いて４ビット（１６コード）に変換される。分数Ｆ１〜Ｆ６の最大値は１６０（０ｘＡ０）に等しいと仮定する。
【数１５】
2007518096.tif 000017
ここで、Ｆcompは対応するチューブ分数Ｆ１〜Ｆ６に適用される加重平方根圧縮であり、ｋは当該チューブ分数を圧縮するために使用される一定スカラ、ｋ＝(Codemax)２／Ｆｍａｘ、であり、Codemaxは２bits−１＝１５に等しいような各分数Ｆ１〜Ｆ６に対する最大コードであり、Ｆmaxは１６０に等しいような各分数Ｆ１〜Ｆ６に対する最大値であり、ｋ＝１５２／１６０＝1.40625である。
【００４０】
符号化関数（Ｊ）により符号化された分数Ｆ１〜Ｆ６の復号即ち伸張は、
【数１６】
2007518096.tif 000018
なる関数を用いることにより達成される。
【００４１】
以上、本発明を好ましい実施例を参照して説明した。明らかなことに、上述した詳細な説明を精読及び理解すれば、当業者は変更例及び変形例を思いつくであろう。本発明は、添付請求項又はそれらの均等物の範囲内に入る限りにおいて、全ての斯かる変形例及び変更例をも含むと見なされることを意図するものである。
【図面の簡単な説明】
【００４２】
【図１】図１は、診断画像化システムの概念図である。
【図２】図２は、診断画像化システムの一部の概略図である。
【図３】図３は、本発明によるソフトウェアプログラムの一部である。
【図４】図４は、本発明による一符号化方法の概念図である。
【図５】図５は、水平対称性を使用する４チューブ群の概念図である。
【図６】図６は、本発明による他の符号化方法の概念図である。
【図７】図７は、イベントが外側チューブの近くで発生した場合の７チューブ群の概念図である。
【図８】図８は、本発明による他の符号化方法の概念図である。
【図９】図９は、イベントが高いチューブの中心の近くで発生した場合の７チューブ群の概念図である。
【図１】
2007518096.tif 000019
【図２】
2007518096.tif 000020
【図３】
2007518096.tif 000021
【図４】
2007518096.tif 000022
【図５】
2007518096.tif 000023
【図６】
2007518096.tif 000024
【図７】
2007518096.tif 000025
【図８】
2007518096.tif 000026
【図９】
2007518096.tif 000027
【国際調査報告】
2007518096.tif 000028
2007518096.tif 000029