社団法人日本放射線技術学会
核医学分科会SPECT データ処理評価用
デジタルファントム
（物理特性評価）取 扱 説 明 書

ご使用前に、必ずこの「取扱説明書」をお読みください。お読みになった後も、大切に保管してください。

発行

２００７，４

はじめに日本放射線技術学会 学術調査研究班「核医学分野におけるデジタルファントム作成研究」により作成されたデジタルファントム（物理性能評価用）をご利用頂き、ありがとうございます。この取扱説明書は、
デジタルファントム（物理性能評価用）の取り扱いについて記載してあります。このデジタルデータを安全に使用していただくために、ご使用の前にこの取扱説明書をお読みになり、十分理解してください。

核医学分野におけるデジタルファントム作成班メンバー
本村信篤（東芝メディカルシステムズ） 班長
大西英雄（県立広島大学大学院）
白川誠士（藤田保健衛生大学）
高橋正昭（中村記念病院）
柳沢正道（日本メジフィジックス）

作成協力者
尾川浩一（法政大学）： シミュレーションによるデータ作成
志田原美保（放射線医学総合研究所）： シミュレーションによるデータ作成
栗原英之（GE 横河メディカル）： DICOM フォーマット変換
金田明義（東芝医用システムエンジニアリング）： フォーマット変換など
千葉県核医学技術研究会： データ読込み（インポート）確認など
日本核医学会、日本核医学技術学会(1)

本書の一部、または全部を無断で転載することは禁止されています。

(2)

本書の内容については、将来予告なしに変更することがあります。

(3)

本書の内容については万全を期しておりますが、万一ご不審な点や、誤り、
記載漏れなど、お気づきの点がありましたらご連絡ください。

© （社）日本放射線技術学会
ALL RIGHTS RESERVED
C - 1

＊

安全上の注意
１．シグナルワードの定義
取扱説明書では安全および重要事項について、危険、警告、注意のシグナルワードを使用しています。それぞれのシグナルワードとその意味は次のようになっています。内容を良く理解してから、本文をお読みください。
シグナルワード

危険
警告
注意

意

味

この警告に従わないと死亡または重傷を招く差し迫った危険な状況を示しています。
この警告に従わないと死亡または重傷を招く可能性がある潜在的に危険な状況を示しています。
この警告に従わないと軽傷または中程度の傷害を招くことがある潜在的に危険な状態を示しています。
この警告に従わないと財産の損害を引き起こすおそれがある状況を示し

注意

ています。

２．安 全 上 の 注 意
本品をご使用いただく上で患者および操作者の安全を守るために、このデジタルデータをご利用になるデータ処理装置、システムの取扱説明書に記載の注意事項もよくお読みになるようお願いします。

注意本デジタルデータのデジタルデータの使用責任は
使用責任は、使用者にあります
使用者にあります。
にあります。本デジタルデータを
デジタルデータをデータ処理装置
データ処理装置、処理装置、シンチ
レーションカメラ装置
レーションカメラ装置などに装置などに据付などに据付け据付けることで生ることで生じた不具合
じた不具合、
不具合、故障、故障、データ消失
データ消失、消失、検査の検査の不実施など不実施などの責任は責任は全て使用者が使用者が負うことになります。
うことになります。日本放射線技術学会、
日本放射線技術学会、また機器
また機器メーカ
機器メーカはメーカは一切の一切の責任を負いません。
いません。この点
この点を御了解の御了解の上、使用をお使用をお願をお願い致します。
します。

S - 1

＊

使用・管理に関する重要情報(1)

本品は、使用者の責任のもとで維持・管理してください。

(2)

使用者の責任で、本品を変更／改造することはできます。但し、変更／改造の責任は全て使用者にあります。

(3)

この取扱説明書には、当作成班が予見し得る範囲で危険性について警告文で記載してあります。
それ以外の危険性についても十分注意してください。

(4)

本品は、マルウエア（コンピュータウイルスやワームなど、感染したコンピュータに被害を与える悪意あるソフトウエア）感染防止のセキュリティ管理がされているネットワークに接続してください。

W - 1

免責事項について(1)

本品は使用者の責任で行ってください。日本放射線技術学会、装置メーカなどは一切の責任を負
いません。

(2)

本品が変更、改造された場合は、一切の動作を保障しません。

(3)

本品により核医学データ処理装置、パーソナルコンピュータが故障、データ消失などを起こした場
合に、日本放射線技術学会、装置メーカなどは一切の責任を負いません。

(4)

マルウエア（コンピュータウイルスやワームなど、感染したコンピュータに被害を与える悪意あるソ
フトウエア）感染を原因とする次の事象について日本放射線技術学会、装置メーカなどは一切責任を負いません。
・製品に記憶されたデータ（臨床データ等）の消去・改ざん・漏洩・製品の誤動作による事故
・この製品を媒体とした他の製品への感染、および感染による被害
・その他、マルウエア感染を原因とするすべての事象

W - 2

＊

デジタルデータの許諾範囲
本品に使用しているデジタルデータは、日本放射線技術学会が著作権などすべての権利を所有していています。
したがって、使用者が本品のデジタルデータを使用される際には、下記に示す事項を守って頂く必要があります。

１．使 用 権 の 許 諾(1)

このデジタルデータに関する知的所有権は、使用者に譲渡されません。

(2)

使用者の責任で、本品を変更／改造することはできます。但し、変更／改造の責任は全て使用者にあります。

L - 1＊

このマニュアルについて
１．表 記 上 の 約 束
この取扱説明書では、安全上の注意事項を表示するシグナルワード（「安全上の注意」参照）の他に、以下のような表示を使用しています。

ご参考

操作に当たっての有効な知識・情報などの内容を示してあります。

２．本 書 に つ い て
本書は、デジタルファントム（物理性能評価用）の取り扱いについて記載してあります。

３．サポートについて
本デジタルファントムの配布に関する問合せは、下記まで電子メールでお願い致します。
現時点で技術的なサポートは実施致しません。本取扱説明書をお読みの上、使用者の責任で対応をお願い致します。
なお、将来的には核医学分科会のホームページにて技術サポートを予定しています。

問合せ
問合せ先
本村信篤（もとむらのぶとく）
E-mail nobutoku.motomura@toshiba.co.jp

注意
本デジタルデータのデジタルデータの使用責任は
使用責任は、使用者にあります
使用者にあります。
にあります。本デジタルデータのデジタルデータのデータ処理装置
データ処理装置、処理装置、シンチレーションカメラ装置
チレーションカメラ装置などへの装置などへの据付などへの据付けなどは据付けなどは使用者
けなどは使用者により使用者により行
により行って頂って頂きます。
きます。日本放射線技術学会、また機器
また機器メーカ
機器メーカはメーカはサポートのサポートの義務がありません
義務がありません。
がありません。この点
この点を御了解の御了解の上、使用をお使用をお願をお願い致します。

注意
本 DICOM 形式のデータは、すべてのメーカでの動作確認は行っていません。DICOM 対応機種であっても、動作しない場合があります。DICOM データの送信に対応していても、受信は別ソフトになっている場合なども考えられます。
また装置メーカに DICOM 対応を依頼する場合は、有償となる場合があります。

U - 1

目

次
ページ

安 全 上 の 注 意 ······························································································································································· S - 1
使用・管理に関する重要情報 ························································································································ W - 1
免責事項について ······················································································································································· W - 3
ソフトウエアの許諾範囲 ········································································································································· L - 1
このマニュアルについて ······································································································································· U - 1
第１章

作成の経緯と目的 ···································································································································· 1

第２章

デジタルファントムについて ··········································································································· 2

２．１

ファントム形状 ······················································································································································· 2

２．２

シミュレーション ····················································································································································· 4

２．３

データの種類 ························································································································································· 5

２．４

データのフォーマットおよびデータ型 ··········································································································· 6

第３章

データ読込み（インポート）方法 ································································································· 7

３．１

データのファイル名 ············································································································································· 7

３．２

DICOM 対応データ処理装置 ······················································································································· 13

３．２

Prominence Processor ···································································································································· 13

３．３

GMS-5500A ····························································································································································14

３．４

他のデータ処理装置 ······································································································································· 14

第４章

実施例 ·································································································································································· 15４．１

濃度（コントラスト） ·········································································································································· 15４．２

部分容積効果······················································································································································ 16

４．３

空間分解能（コリメータの違い） ·················································································································· 17

４．４

空間分解能（回転半径の違い） ·················································································································· 18

４．５

前処理フィルタ（Butterworth フィルタ） ··································································································· 19

４．６

再構成法（FBP と OSEM の違い） ············································································································· 21

４．７

再構成法（OSEM パラメータの違い） ······································································································· 23４．８

散乱線補正··························································································································································· 24４．９

散乱線（コリメータによる違い） ······································································································· 26

４．１０

減弱補正（均一性） ······································································································································ 28

４．１１

減弱補正（コントラスト）······················································································································· 29

- a -

＊

第１章

作成の経緯と目的

現在、核医学診断装置のメーカは機種によって異なる独自の画質の画像を作り出しているのが現状である。また、核医学検査に携わっている放射線技師の中にはそれぞれのメーカが推奨するパラメータを用いて画像処理を行って、なぜ、そのような既定値を用いているか、さまざまな状況でどのように値を変えていくべきかを理解していない場合が多く見られる。収集条件でのマトリックスサイズの決定方法、
ピクセル分解能の理解、ナイキスト周波数に代表されるようなピクセルサイズと周波数との関係、散乱・減弱の正確な意味付けとそれに伴う現象の補正効果など核医学検査、特に SPECT 画像処理に必要不
可欠な部分が、かなりなおざりにされて検査が行なわれている。その現状を把握するために第３２回秋
季学術大会の核医学分科会にてデジタル画像収集に関する理解度をアンケート方式で実施した結果、上記に述べた内容があまり十分理解されずに使用しているのが判明した。一方、データ収集や処理法を比較検討するためのデジタルファントムとして、人体の構造を模擬した MCAT（Myocardial Torso）ファント
ムなどが存在する。しかし、このファントムは人体構造を模擬しているために異なる均一性、空間分解能、
コントラスト分解能、ピクセル分解能などの基本性能を評価することが出来ない構造となっている。
このような現状から、核医学分科会は性能評価等のデジタルファントム作成は急務であると感じ、2005
年度日本放射線技術学会の班研究の一環として、我々が提案したデジタルファントム作成研究班と、おりしも同時期に構想を練っておられた日本核医学会の尾川先生らのデジタルファントム研究班との共同企画でこの SPECT 評価用デジタルファントムの立案と作成が始まった。このデジタルファントムは、基礎的な部分を理解する基礎ファントム（SPECT データ評価用デジタルファントム＜物理特性評価＞）とより臨床に近い構造を持つ脳デジタルファントムの２種類から構成されている。前者の物理的評価ファントムの立案及び作成を放射線技術学会の研究班が担当することになった。それぞれのファントムは基本的に、散乱・減弱もない理想的な条件のデータと散乱・減弱等を考慮しモンテカルロシュミレーション（EGS4）法で作成したデータのファントムを用意した。データ形式に関しては、基本的に既存の核医学専
用データ処理装置上で処理が行なえることを念頭に置き、DICOM 形式の投影データを作成した。このDICOM 形式を採用することで各処理装置にインポート（入力）出来るようにして汎用性を高めている。また、DICOM 形式以外の各種のデータ形式も併せて作成した。
今回、SPECT データ評価用デジタルファントム＜物理特性評価＞をリリースすることで、我々は核医学従事者の基本的教育や、日常業務における基本画像処理をより深く理解することができると考える。また、このファントムを利用することで画像処理の内容を十分理解し、国内の各病院における核医学画像の標準化に役立てば幸いです。尚、これらのデジタルファントム等はインターネットに広く公開し核医学
デジタル画像処理の教育・啓蒙に役立てるものと考える。
平成 19 年 3 月吉日
日本放射線技術学会核医学分科会会長
大西 英雄

- 1 -

第２章
２．１

デジタルファントムについて
ファントム形状
ファントムの形状は、容易に SPECT の基本性能を確認、評価できることを目的とした。評価対象とした基本性能は、濃度（コントラスト、部分容積効果）、位置分解能、均一性である。
Fig.1 左はファントムの全体像を示したものであり、ファントム上部は濃度（コントラスト、部分
容積効果）、位置分解能の検討を行うために、それぞれ濃度の異なるロッド領域（バックグ
ラウンド領域の濃度を１として０、２、４）を設けてある。また、ファントム下部は均一性評価用に、均一濃度領域を設けてある。Fig.1 右は各ロッドの拡大図である。60、40、20、10、4mm
径、30mm 高の円柱が連続して配置されている。
Fig.2 は各ロッド径と均一部の断面像を示す。

60mm
10mm
（均一性領域）
150mm（ロッド領域）

30mm

40mm
30mm

40mm
（均一性領域）
200mm

20mm
30mm

濃度：4

150mm
60mm

10mm
30mm

濃度：0
濃度：2
コールド
バックグラウンド
濃度：１

4mm
30mm

Fig.1 ファントム形状（ロッドの位置と大きさ）- 2 -60mm径部

40mm径部

20mm径部

10mm径部

4mm径部

均一部

Fig.2 各ロッド径と均一部の断面像

表１ 各部位のスライス位置部位

2mm/pixel (128Matrix)

4mm/pixel (64Matrix)

6mm/pixel (64Matrix)

ブランク（データ０）

１～１４

１～７

１～１４

均一部（濃度１）

１５～１９

８～９

60mm 径部

２０～３４

１１～１７

40mm 径部

３５～４９

１８～２４

20mm 径部

５０～６４

２６～３２

10mm 径部

６５～７９

３３～３９

4mm 径部

８０～９４

４１～４７

均一部（濃度１）

９５～１１４

４８～５７

３７～４７

ブランク（データ０）

１１５～１２８

５８～６４

４８～６４

- 3 -

１０

２５
４０

１５～１６
１７～２１
２２～２６
２７～３１
３２～３６
―――――

２．２

シミュレーション
２．１章で示したファントム形状にて EGS4 によるモンテカルロシミュレーションで SPECT 投影データを作成した。収集パラメータは、コリメータ（LEHR：低エネルギー高分解能、LEGP：
低エネルギー汎用）、検出器回転半径・軌道（150mm 円軌道、200mm 円軌道）、ピクセルサ
イズ（2mm、4mm、6mm）、投影ステップ角（3 度）、最大カウント（50、100）とした。RI は全て
Tc-99m（140keV の単一光電ピーク）とした。
LEHR コリメータの入射許容角度は 2.5 度、LEGP コリメータの入射許容角度は 3.6 度とした。
またエネルギー分解能は 10%とした。検出器の固有分解能は考慮していない。
それぞれにおいて全光子、一次線のみ、散乱線のみ、および Triple Energy Window（TEW）法に基づく低エネルギーサブウインドウのデータが存在する。すなわち・20%幅光電ピークの全光子データ（一次線＋散乱線）
・20%幅光電ピークの一次線データ（一次線のみ）
・20%幅光電ピークの散乱線データ（散乱線のみ）
・7%幅低エネルギー側サブウインドウの全光子データ（一次線＋散乱線）

20%幅光電ピークウインドウ

Counts7%幅低エネルギー側サブウインドウ

一次線

散乱線

50

100
Energy (keV)

- 4 -

140

２．３

データの種類ａ）モンテカルロシミュレーションで作成した SPECT 投影データ
収集パラメータを組み合わせで作成したデータの種類は下表の通り。
No.

コリメータ

回転半径

ピクセル

ステップ角

最大値

１

LEHR

200mm

4mm

3度

100

２

LEGP

200mm

4mm

3度

100

３

LEHR

150mm

4mm

3度

100

４

LEHR

150mm

2mm

3度

100

５

LEHR

200mm

6mm

3度

100

６

LEGP

200mm

6mm

3度

100

７

LEHR

200mm

4mm

3度

50

それぞれに対して、以下のエネルギーウインドウ、一次線／散乱線のデータがある。
・20%幅光電ピークウインドウの全光子データ
・20%幅光電ピークウインドウの一次線データ
・20%幅光電ピークウインドウの散乱線データ
・7%幅低エネルギー側サブウインドウの全光子データｂ）理想の SPECT 投影データ比較のため、ピクセルサイズごとに理想的な SPECT 投影データ（コリメータ開口による位置分
解能劣化なし／散乱、吸収なし／統計ノイズなし）がある。
・ピクセルサイズ 2mm での理想投影データ
・ピクセルサイズ 4mm での理想投影データ
・ピクセルサイズ 6mm での理想投影データ
ｃ）理想の SPECT 断層像データ比較のため、ピクセルサイズごとに理想的な SPECT 断層像データ（コリメータ開口による位置分解能劣化なし／散乱、吸収なし／統計ノイズなし）がある。
・ピクセルサイズ 2mm での理想断層像データ
・ピクセルサイズ 4mm での理想断層像データ
・ピクセルサイズ 6mm での理想断層像データｄ）減弱係数マップ
減弱係数マップ画像が、減弱係数とピクセルサイズの組合せである。
・ピクセルサイズ 2mm
減弱係数 0.05cm-1、0.10cm-1、0.15cm-1、0.20cm-1・ピクセルサイズ 4mm
減弱係数 0.05cm-1、0.10cm-1、0.15cm-1、0.20cm-1・ピクセルサイズ 6mm
減弱係数 0.05cm-1、0.10cm-1、0.15cm-1、0.20cm-1
- 5 -

２．４

データのフォーマットおよびデータ型

データのフォーマットは、DICOM、Prominence Processor（開発元 核医学画像処理技術カンファレンス）、GMS-5500（東芝メディカルシステムズ）、Raw-data の４種類である。
フォーマットごとにデータが作成されている。
Prominence Processor の減弱係数マップを除き、データ型は全て２バイト符号付整数である。
Prominence Processor の減弱係数マップは４バイトの実数である。
各ファイルの詳細は第３章「データ読込み（インポート）方法」を参照のこと。
＊注意＊
注意＊
DICOM フォーマットのフォーマットのデータでもデータでも、
でも、DICOM 準拠の準拠のデータ処理装置
データ処理装置に処理装置に対応できない
対応できない場合できない場合があります
場合があります。
があります。
第３章にデータ読込
データ読込み読込み（インポート）
インポート）方法を
方法を説明しますが説明しますが、しますが、日本放射線技術学会、
日本放射線技術学会、装置メーカ
装置メーカはメーカは一切のサポートをサポートを行いません。
いません。
＊注意＊
注意＊
GMSGMS-5500 フォーマットはフォーマットは、HG シリーズにシリーズには対応していません
対応していません。
していません。

- 6 -

第３章
３．１

データ読込み（インポート）方法

データのファイル名

デジタルファントムは当作成班が作成／配布した CD-ROM 内に入っている。フォーマットごとにフォルダが分かれている。以下にそれぞれのファイル名について解説する。
ａ） DICOM フォーマット
CD-ROM 内の「DICOM」フォルダにファイルが格納されている。
Ⅰ）モンテカルロシミュレーションで作成した SPECT 投影データ
収集パラメータを組み合わせで作成したデータは、下記のルールで名前が付けられたフォルダ内に格納されている。
フォルダ名
フォルダ名＝コリメータ_
コリメータ_回転半径_
回転半径_ピクセルサイズ_
ピクセルサイズ_ステップ角
ステップ角_最大カウント
最大カウント
フォルダ名

コリメータ

回転半径

ピクセル

ステップ角

最大値

HR_200_4mm_3d_100

LEHR

200mm

4mm

3度

100

GP_200_4mm_3d_100

LEGP

200mm

4mm

3度

100

HR_150_4mm_3d_100

LEHR

150mm

4mm

3度

100

HR_150_2mm_3d_100

LEHR

150mm

2mm

3度

100

HR_200_6mm_3d_100

LEHR

200mm

6mm

3度

100

GP_200_6mm_3d_100

LEGP

200mm

6mm

3度

100

HR_200_4mm_3d_50

LEHR

200mm

4mm

3度

50

各フォルダには、以下のエネルギーウインドウ、一次線／散乱線ごとのファイルがある。
Prj_M_Total.dcm： 20%幅光電ピークウインドウの全光子データ
Prj_M_True.dcm： 20%幅光電ピークウインドウの一次線データ
Prj_M_Scat.dcm： 20%幅光電ピークウインドウの散乱線データ
Prj_LoSub.dcm： 7%幅低エネルギー側サブウインドウの全光子データ
Ⅱ）理想の SPECT 投影データ
理想の SPECT 投影データは、フォルダ「Ideal_Phantom_Proj」内に格納されている。ピクセルサ
イズごとに以下のファイルがある。
Proj128_2mm_3d.dcm： ピクセルサイズ 2mm での理想投影データ
Proj64_4mm_3d.dcm： ピクセルサイズ 4mm での理想投影データ
Proj64_6mm_3d.dcm： ピクセルサイズ 6mm での理想投影データ
Ⅲ）理想の SPECT 断層像データ
理想の SPECT 断層像データは、フォルダ「Ideal_Phantom_Slice」内に格納されている。ピクセ
ルサイズごとに以下のファイルがある。
Slice128_2mm.dcm： ピクセルサイズ 2mm での理想断層像データ- 7 -Slice64_4mm.dcm： ピクセルサイズ 4mm での理想断層像データ
Slice64_6mm.dcm： ピクセルサイズ 6mm での理想断層像データ
Ⅳ）減弱係数マップ
減弱係数マップデータは、フォルダ「Attenuation_map」内に格納されている。ピクセルサイズごと、減弱係数ごとに以下のファイルがある。
・ピクセルサイズ 2mm
Mu005_128_2mm.dcm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_128_2mm.dcm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_128_2mm.dcm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu020_128_2mm.dcm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 4mm
Mu005_64_4mm.dcm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_64_4mm.dcm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_64_4mm.dcm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu0205_64_4mm.dcm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 6mm
Mu005_64_6mm.dcm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_64_6mm.dcm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_64_6mm.dcm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu020_64_6mm.dcm： 減弱係数 0.20cm-1 の減弱係数マップｂ） Prominence Processor フォーマット
CD-ROM 内の「Prominence」フォルダにファイルが格納されている。
Ⅰ）モンテカルロシミュレーションで作成した SPECT 投影データ
収集パラメータを組み合わせで作成したデータは、下記のルールで名前が付けられたフォルダ内に格納されている。
フォルダ名
フォルダ名＝コリメータ_
コリメータ_回転半径_
回転半径_ピクセルサイズ_
ピクセルサイズ_ステップ角
ステップ角_最大カウント
最大カウント
フォルダ名

コリメータ

回転半径

ピクセル

ステップ角

最大値

HR_200_4mm_3d_100

LEHR

200mm

4mm

3度

100

GP_200_4mm_3d_100

LEGP

200mm

4mm

3度

100

HR_150_4mm_3d_100

LEHR

150mm

4mm

3度

100

HR_150_2mm_3d_100

LEHR

150mm

2mm

3度

100

HR_200_6mm_3d_100

LEHR

200mm

6mm

3度

100

GP_200_6mm_3d_100

LEGP

200mm

6mm

3度

100

HR_200_4mm_3d_50

LEHR

200mm

4mm

3度

50

各フォルダには、以下のエネルギーウインドウ、一次線／散乱線ごとのファイルがある。
- 8 -Tc_M_total__M.sps： 20%幅光電ピークウインドウの全光子データ
Tc_M_true_M.sps： 20%幅光電ピークウインドウの一次線データ
Tc_M_scat__M.sps： 20%幅光電ピークウインドウの散乱線データ
Tc_LowSub_L.sps： 7%幅低エネルギー側サブウインドウの全光子データ
Ⅱ）理想の SPECT 投影データ
理想の SPECT 投影データは、フォルダ「Ideal_Phantom_01_Prj」内に格納されている。ピクセル
サイズごとに以下のファイルがある。
Proj2mm_3deg__i__M.sps： ピクセルサイズ 2mm での理想投影データ
Proj4mm_3deg__i__M.sps： ピクセルサイズ 4mm での理想投影データ
Proj6mm_3deg__i__M.sps： ピクセルサイズ 6mm での理想投影データ
Ⅲ）理想の SPECT 断層像データ
理想の SPECT 断層像データは、フォルダ「Ideal_Phantom_01_SPECT」内に格納されている。ピクセルサイズごとに以下のファイルがある。
IdealSl2mm_i.sts： ピクセルサイズ 2mm での理想断層像データ
IdealSl4mm_i.sts： ピクセルサイズ 4mm での理想断層像データ
IdealSl6mm_i.sts： ピクセルサイズ 6mm での理想断層像データ
Ⅳ）減弱係数マップ（データタイプは実数（４バイト））
減弱係数マップデータは、フォルダ「Attenuation_map」内に格納されている。ピクセルサイズごと、減弱係数ごとに以下のファイルがある。
・ピクセルサイズ 2mm
Map_Mu005_2mm.im.sts： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_2mm.im.sts： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_2mm.im.sts： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_2mm.im.sts： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 4mm
Map_Mu005_4mm.im.sts： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_4mm.im.sts： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_4mm.im.sts： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_4mm.im.sts： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 6mm
Map_Mu005_6mm.im.sts： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_6mm.im.sts： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_6mm.im.sts： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_6mm.im.sts： 減弱係数 0.20cm-1 の減弱係数マップ- 9 -ｃ） GMSGMS-5500 フォーマット （＊HG（＊HG シリーズにシリーズには適応しない
適応しない＊）
しない＊）
CD-ROM 内の「GMS-5500」フォルダにファイルが格納されている。tar フォーマットファイルのDigital_phantom_03.tar を GMS-5500 上で解凍する。手順の詳細は３．３項を参照。
Ⅰ）モンテカルロシミュレーションで作成した SPECT 投影データ
収集パラメータを組み合わせで作成したデータは、患者 Name は「Digital_phantom_01」で、患者 ID は下記のルールで付けられる。
フォルダ名
フォルダ名＝コリメータ_
コリメータ_回転半径_
回転半径_ピクセルサイズ_
ピクセルサイズ_ステップ角
ステップ角_最大カウント
最大カウント
フォルダ名

コリメータ

回転半径

ピクセル

ステップ角

最大値

HR_200_4mm_3d_100

LEHR

200mm

4mm

3度

100

GP_200_4mm_3d_100

LEGP

200mm

4mm

3度

100

HR_150_4mm_3d_100

LEHR

150mm

4mm

3度

100

HR_150_2mm_3d_100

LEHR

150mm

2mm

3度

100

HR_200_6mm_3d_100

LEHR

200mm

6mm

3度

100

GP_200_6mm_3d_100

LEGP

200mm

6mm

3度

100

HR_200_4mm_3d_50

LEHR

200mm

4mm

3度

50

エネルギーウインドウ、一次線／散乱線ごとのファイルは Title で識別される。
Tc_M_total： 20%幅光電ピークウインドウの全光子データ
00Tc_M_true： 20%幅光電ピークウインドウの一次線データ
Tc_M_scat： 20%幅光電ピークウインドウの散乱線データ
Tc_LowSub： 7%幅低エネルギー側サブウインドウの全光子データ
Ⅱ）理想の SPECT 投影データ
理想の SPECT 投影データは、患者 Name は「Digital_phantom_01」で、患者 ID は
「Ideal_Phantom_01」である。ファイルはピクセルサイズごとに Title で識別される。
2mm_3d__-i： ピクセルサイズ 2mm での理想投影データ
4mm_3d__-i： ピクセルサイズ 4mm での理想投影データ
6mm_3d__-i： ピクセルサイズ 6mm での理想投影データ
Ⅲ）理想の SPECT 断層像データ
理想の SPECT 断層像データは、患者 Name は「Digital_phantom_01」で、患者 ID は
「Ideal_Phantom_01」である。ファイルはピクセルサイズごとに Title で識別される。
IdealSl2mm： ピクセルサイズ 2mm での理想断層像データ
IdealSl4mm： ピクセルサイズ 4mm での理想断層像データ
IdealSl6mm： ピクセルサイズ 6mm での理想断層像データ
Ⅳ）減弱係数マップ
減弱係数マップデータは、患者 Name は「Digital_phantom_01」で、患者 ID は「Attenuation_map」
である。ファイルはピクセルサイズごと、減弱係数ごとに Title で識別される。
- 10 -

・ピクセルサイズ 2mm
Mu005_2mm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_2mm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_2mm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu020_2mm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 4mm
Mu005_4mm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_4mm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_4mm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu020_4mm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 6mm
Mu005_6mm： 減弱係数 0.05cm-1 の減弱係数マップ
Mu010_6mm： 減弱係数 0.10cm-1 の減弱係数マップ
Mu015_6mm： 減弱係数 0.15cm-1 の減弱係数マップ
Mu020_6mm： 減弱係数 0.20cm-1 の減弱係数マップｄ） Raw
Rawaw-data フォーマット
Raw-data はファイルのヘッダー情報がない画像データのみからなるファイルである。CD-ROM 内の「Raw-data」フォルダにファイルが格納されている。
Ⅰ）モンテカルロシミュレーションで作成した SPECT 投影データ
収集パラメータを組み合わせで作成したデータは、下記のルールで名前が付けられたフォルダ内に格納されている。
フォルダ名
フォルダ名＝コリメータ_
コリメータ_回転半径_
回転半径_ピクセルサイズ_
ピクセルサイズ_ステップ角
ステップ角_最大カウント
最大カウント
フォルダ名

コリメータ

回転半径

ピクセル

ステップ角

最大値

HR_200_4mm_3d_100

LEHR

200mm

4mm

3度

100

GP_200_4mm_3d_100

LEGP

200mm

4mm

3度

100

HR_150_4mm_3d_100

LEHR

150mm

4mm

3度

100

HR_150_2mm_3d_100

LEHR

150mm

2mm

3度

100

HR_200_6mm_3d_100

LEHR

200mm

6mm

3度

100

GP_200_6mm_3d_100

LEGP

200mm

6mm

3度

100

HR_200_4mm_3d_50

LEHR

200mm

4mm

3度

50

各フォルダには、以下のエネルギーウインドウ、一次線／散乱線ごとのファイルがある。
Prj_M_Total： 20%幅光電ピークウインドウの全光子データ
Prj_M_Ture： 20%幅光電ピークウインドウの一次線データ
Prj_M_Scat： 20%幅光電ピークウインドウの散乱線データ
Prj_LowSub： 7%幅低エネルギー側サブウインドウの全光子データ- 11 -Ⅱ）理想の SPECT 投影データ
理想の SPECT 投影データは、フォルダ「Ideal_Phantom_01_Prj」内に格納されている。ピクセル
サイズごとに以下のファイルがある。
Proj128_2mm_3d： ピクセルサイズ 2mm での理想投影データ
Proj64_4mm_3d： ピクセルサイズ 4mm での理想投影データ
Proj64_6mm_3d： ピクセルサイズ 6mm での理想投影データ
Ⅲ）理想の SPECT 断層像データ
理想の SPECT 断層像データは、フォルダ「Ideal_Phantom_01_SPECT」内に格納されている。ピクセルサイズごとに以下のファイルがある。
Ideal_Slice128_2mm： ピクセルサイズ 2mm での理想断層像データ
Ideal_Slice64_4mm： ピクセルサイズ 4mm での理想断層像データ
Ideal_Slice64_6mm： ピクセルサイズ 6mm での理想断層像データ
Ⅳ）減弱係数マップ
減弱係数マップデータは、フォルダ「Attenuation_map」内に格納されている。ピクセルサイズごと、減弱係数ごとに以下のファイルがある。
・ピクセルサイズ 2mm
Map_Mu005_128_2mm： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_128_2mm： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_128_2mm： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_128_2mm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 4mm
Map_Mu005_64_4mm： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_64_4mm： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_64_4mm： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_64_4mm： 減弱係数 0.20cm-1 の減弱係数マップ
・ピクセルサイズ 6mm
Map_Mu005_64_6mm： 減弱係数 0.05cm-1 の減弱係数マップ
Map_Mu010_64_6mm： 減弱係数 0.10cm-1 の減弱係数マップ
Map_Mu015_64_6mm： 減弱係数 0.15cm-1 の減弱係数マップ
Map_Mu020_64_6mm： 減弱係数 0.20cm-1 の減弱係数マップ- 12 -

３．２

DICOM 対応データ処理装置
デジタルファントム保存 CD-ROM の「DICOM Format」フォルダの中にあるデータを使用する。
データ読込みに関する詳細な手順は、各装置の取扱説明書を参照のこと。

注意
本 DICOM 形式のデータは、すべてのメーカでの動作確認は行っていない。DICOM 対応機種であっても、動作しない場合がある。DICOM データの送信に対応していても、受信は別ソフトになっている場合なども考えられる。
また装置メーカに DICOM 対応を依頼する場合は、有償となる場合がある。
ａ） e.soft ワークステーション
東芝メディカルシステムズ社製、シーメンス旭メディック社製のＤＩＣＯＭ対応データ処理装置である。東芝メディカルシステムズ社製の GMS-7700 も同様である。
1) デジタルファントム保存 CD-ROM を e.soft ワークステーションの CD-ROM ドライブにセットする。
2) Windows エクスプラーラを使って、CD-ROM 内の「DICOM Format」フォルダにあるファ
イルから据付けたいファイルを D:\Temp フォルダの下にコピーする。必要に応じて、
D:\Temp フォルダの下にフォルダを作ってもよい。
3) e.soft データベースのメニューバーから「被験者」－＞「ブラウザ」を選択する。患者ブ
ラウザのウインドウが表示される。
4) 患者ブラウザのメニューバーから「転送」－＞「オフラインからのインポート」を選択する。
5) ファイル選択ウインドウにて、デジタルファントムのファイルが保存されているディレク
トリィ（D:\Temp、もしくは作成したフォルダ）を選択し、対象ファイルを選択する。
6) “OK”ボタンをクリックし、データベースへの据付を開始する。
7) e.soft のファイルリスト表示画面にて「リストの更新」ボタンをクリックし、据付けたファ
イルがデータベースにあることを確認する。
8) CD-ROM を取り出す。
b） eNTEGRA （Ver. 2.5302）
2.5302）ワークステーション
GE 横河メディカル社製のＤＩＣＯＭ対応データ処理装置である。
1) eNTEGRA アプリケーションを最小化して、D/Einstein/Einstein_user/DicomP10/Import
フォルダに DICOM データをコピーする。
2) eNTEGRA の Conversion タブを選択し、さらに Import タブを選択、Change Configuration
ボタンをクリックする。
3) 表示されたウィンドウで、・ Import from DICOM Part10…
・ Local Station
・ ../../Einstein_user/DicomP10…
を選択して、Apply&Quit をクリックする。
- 13 -

4) リストにデータ名が表示されます。
それぞれをクリックして選択、または Select All ですべて選択して Import する。
5) Import されたプラナーデータを選択した状態で、BrainSPECT をスタートアイコンに飛ばし、Customize..
Customize..ボタンをクリックする。
Customize..
6) 開いたカスタマイズのウィンドウで Emission を選択し、Add/Modify の Names のプル
ダウンリストから該当データの名前を選択し、Group プルダウンから EARLY を選択、Add すると Names リストにリストアップされる。Apply
Apply&Quit
Apply&Quit で画面を閉じる。
ｃ） Odyssey ワークステーション
島津社製のＤＩＣＯＭ対応データ処理装置である。本装置は、本デジタルファントムデータの適用対象外である。適用に必要な DICOM タグの書換えなどは、使用者の自己責任の下に実施する必要がある。
ｄ） Pegasys ワークステーション
日立メディコ社製、Philips 社製のＤＩＣＯＭ対応データ処理装置である。本装置に関して、
データ読込みの検証は行われていない。

３．２

Prominence Processor
Prominence Processor は核医学画像処理技術カンファレンスが開発した Windows PC 上で動作するソフトウエアである。
核医学画像処理技術カンファレンスデジタルファントム保存 CD-ROM の「Prominence
Format」フォルダの中にあるデータを使用する。以下に操作手順を記す。
1) デジタルファントム保存 CD-ROM を Windows PC の CD-ROM ドライブにセットする。
2) Windows エクスプラーラを使って、CD-ROM 内の「Prominence Format」フォルダにあるファイルから据付けたいファイルを任意のフォルダの下にコピーする。
3) Prominence Processor を起動し、据付けたファイルが Import できることを確認する。
4) CD-ROM を取り出す。

３．３

GMS-5500A
対象は東芝メディカルシステムズ社製 GMS-5500A/DI、UI、PI である。HG シリーズはサポ
ートしていない。デジタルファントム保存 CD-ROM の「GMS Format」フォルダの中にあるデータを使用する。Windows PC と GMS-5500 がネットワーク接続されていることを前提に操作手順を記す。
1) デジタルファントム保存 CD-ROM を Windows PC の CD-ROM ドライブにセットする。
2) Windows エクスプラーラを使って、CD-ROM 内の「GMS-5500 Format」フォルダにある"
Digital_phantom_03.tar "ファイルを FTP などで GMS-5500A の/GMS_HD の下にコピーする。
3) GMS-5500A にて tar ファイルを解凍する。

- 14 -

cd /GMS_HD
tar xvf Digital_phantom_03.tar
4) GMS-5500A にて" Digital_phantom_03.tar "ファイルを削除する。
cd /GMS_HD
rm DigiPhantom_GMS5500.tar
5) CD-ROM を取り出す。

３．４

他のデータ処理装置
ヘッダーなしフォーマットは Raw Data として汎用画像処理ソフトで扱うことが可能である。デ
ジタルファントム保存 CD-ROM の「Raw Data Format」フォルダの中にあるデータを使用する。
この場合、マトリクスサイズ（ピクセルサイズが 2mm の場合は 128x128、4mm と 6mm は64x64）と投影数（=360 度/投影角度）および２バイト符号付整数データであることを考慮す
る必要がある。
再構成処理を行うために必要なソフトは各自で準備する。

- 15 -

第４章

実施例

４．１ 濃度（コントラスト）
目的：
目的： 収集条件、処理、あるいは補正などの評価法の一つとして、コントラストを算出する。
使用データ
使用データ：
データ： GP_200_4mm_3d_100_Total_M
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件：
再構成条件 前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

オリジナルデータ（ROI サイズ 100%）

オリジナルデータ（ROI サイズ 65%）

GP_200_4mm_3d_100_Total_M（ROI サイズ 100%）

GP_200_4mm_3d_100_Total_M（ROI サイズ 65%）

評価方法：
評価方法
上記に示す処理条件にて、画像再構成を行った。
オリジナルファントムデータのロッド径 40mm のスライスにおいて、関心領域- 16 -（ROI）をロッド径の 100%、65%のサイズで設定した。ROI 内の平均カウント
を次式に代入してコントラストを算出する。

C=

a−b
a+b

仮に b=0 ならば、コントラストは 1.0 となり、a と b の差がなければコントラス
トは 0.0 となる。
また、ROI のサイズの設定にも注意が必要である。ロッド径の 100%の設定では、部分容積効果の影響によりコントラストを過小評価してしまうため、70%程度に設定する。

ROI サイズ 100%
ROI サイズ 65%

コントラスト
0.53
0.66

４．２ 部分容積効果
目的：
目的： 核医学画像は空間分解能、統計ノイズの両方共に劣っているために標本化間隔が粗く、画素サ
イズを大きく取らざるを得ない（空間分解能は悪くなる）。そこで同じ量の放射能が入っているにも関わらず、標的物の大きさによっては同じ放射能には成らない。これを部分容積効果と呼び特に核医学画像ではこの現象により定量性劣化の原因となる。この効果を実験的に求める。
使用データ
使用データ：
データ： HR_200_4mm_3d_100_Total_M
GP_200_4mm_3d_100_Total_M
HR_150_4mm_3d_100_Total_M
使用データ
使用データ処理装置
データ処理装置：
処理装置 GMS-5500A
再構成スライス
再構成スライス：スライス： 7～55
再構成条件：
再構成条件： 前処理フィルタ BW order 8, 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

評価方法：
評価方法：
視覚評価：
視覚評価 上記に示すように、処理条件にて画像再構成を行った。本来，同じカウント値が入っているはずだが，ロッド径が小さくなるとロッド内のカウント値が減少している。

- 17 -

定量的評価：
定量的評価： 各ロッド部に関心領域（ROI）を設定し、ROI 内の最高カウント若しくは平均カウントを算出する。そこである基準カウントで正規化し、ロッドの径別に計数値を求める。求めた結果を横軸にロッド径を縦軸に ROI の計数値をプロットしたグラフを作成する。（今回は、データのばらつきを排除するために、同ロッド径の画像を加算し、60ｍｍ径より大きい矩形 ROI を作成し、全てのロッドに適応して、ROI 内の最高カウント算出し、60mm 径の場所での矩形 ROI 内最高カウントで正規化を行った。
本来、ロッド径が異なっても同じ計数値が入っている限りは、ロッド径の大きさにより計数値が変化しないはず。しかし、結果はロッド径が小さくなるにつれて相対値が減少している。これはロッド径自体の空間分解能と収集時のピクセル分解能との差により生じた部分容積効果である。これを補正する
ために回復係数(recovery coefficient :RC)などが定義されている。コリメータ（LEHR, LEGP）において空間分解能の悪い LEGP が部分容積効果大きく、また回転半径が大きい時も同様に分解能劣化の影響で部分容積効果が増大している。臨床では脳血流 SPECT 画像等で如実にその現象が現れている。

４．３ 空間分解能（コリメータの違い）
目的：
目的： コリメータは特定の方向から入射したγ線のみ検出できるように制御を行う。コリメータの孔径や隔壁の厚さ等、寸法の違いにより空間分解能が異なり、核医学画像自体の空間分解能に密接に関係する。コリメータの違いによる空間分解能の変化を視覚評価およびプロファイルカーブにて
確認する。
使用データ
使用データ：
データ： GP_200_4mm_3d_100_Total_M、
HR_200_4mm_3d_100_Total_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件：
再構成条件： 前処理フィルタ BW order 8, 0.5 [cycles/cm] (0.2 [cycles/pixel])- 18 -

再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

HR_200_4mm_3d_100
HR_200_4mm_3d_100_Total_M_Total_M

GP_200_4mm_3d_100_Total_M
GP_200_4mm_3d_100_Total_M

評価方法
視覚評価：
視覚評価 理想データの処理画像を基準に、それぞれのコリメータの処理画像を比較する。（注：理想データとそれぞれのコリメータの再構成範囲を統一させておく。）視覚的に分解できる内部の最小ロッド径の違いは？
定量評価：
定量評価 各処理画像における内部ロッド径画像のプロファイルカーブを作成する。プロファイルカーブより FWHM（半値幅）が算出可能な場合は、数値を比較検討する。処理画像におけるロッドの確認が困難な場合は理想データにてロッドの位置を確認する。
プロファイルカーブにどのような違いがあるか？FWHM はどのような差があるか？

HR_200_4mm_3d_100
HR_200_4mm_3d_100_Total_M_Total_M

GP_200_4mm_3d_100_Total_M
GP_200_4mm_3d_100_Total_M

４．４ 空間分解能（回転半径の違い）
目的：
目的： SPECT 収集時における回転半径が大きくなると SPECT 画像の空間分解能は大きく劣化する。
回転半径の違いによる空間分解能の変化を視覚評価およびプロファイルカーブにて確認する。
使用データ
使用データ：
データ： HR_150_4mm_3d_100_Total_M、
HR_200_4mm_3d_100_Total_M- 19 -

Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件：
再構成条件： 前処理フィルタ BW order 8, 0.5 [cycles/cm] ([0.2 cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

HR_15
HR_150_4mm_3d_10
150_4mm_3d_100
0_4mm_3d_100_Total_M_Total_M

HR_200_4mm_3d_100_Total_M
HR_200_4mm_3d_100_Total_M

評価方法
視覚評価：
視覚評価 理想データの処理画像を基準に、それぞれのコリメータの処理画像を比較する。（注：理想データとそれぞれのコリメータの再構成範囲を統一させておく。）視覚的に分解できる内部の最小ロッド径の違いは？
定量評価：
定量評価 各処理画像における内部ロッド径画像のプロファイルカーブを作成する。プロファイルカーブより FWHM（半値幅）が算出可能な場合は、数値を比較検討する。処理画像におけるロッドの確認が困難な場合は理想データにてロッドの位置を確認する。
プロファイルカーブにどのような違いがあるか？FWHM はどのような差があるか？

４．５ 前処理フィルタ（Butterworth フィルタ）目的：
目的： 通常の SPECT 収集における投影データは、計数地が少ないため多くの統計雑音を含んでおり、再構成画像の画質を劣化させる。統計雑音を低減させるために前処理フィルタが使用される。その代表的なものである Butterworth フィルタの効果を評価する。Butterworth フィルタのパラメータのうち、遮断周波数を変化させると特に画質は大きく変化する。遮断周波数の違いによる画像の変化を視覚的に確認するとともに、最適遮断周波数について検討する。
使用データ
使用データ：
データ： HR_150_4mm_3d_100_Total_M
Ideal_4mm
使用データ
使用データ処理装置
処理装置：
データ処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件
再構成条件：
条件： 前処理フィルタ BW order 8, cutoff freq. 0.3～0.7 [cycles/cm]まで変化させる（変化の幅は任意）- 20 -

再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

0.3

0.4

0.5

0.6

0.7

[cycles/cm]

評価方法：
評価方法：
視覚評価：
視覚評価 遮断周波数を変化させた場合の画像の違いを視覚的に確認する。遮断周波数を高くしていった（低くしていった）場合に画像はどのように変化するか？また、視覚的に確認できる最小ロッド径に変化はあるか？
定量評価：
定量評価
・各サイズのロッドと同じ大きさの ROI を設定（理想データを参照する）し、ROI 内のカウントを比較検討する。遮断周波数を変化させた場合、カウントはどのように変化するか？
・均一性領域にファントム直径の 80％程度の大きさの ROI を設定し、ROI 内のカウントの標準偏差（SD）または変動係数（CV）を算出する。遮断周波数を変化させた場合、標準偏差または変動係数はどのように変化するか？
・処理画像の各ロッドに ROI 各処理画像における内部ロッド径画像のプロファイルカーブを作成する。
プロファイルカーブより FWHM（半値幅）が算出可能な場合は、数値を比較検討する。処理画像におけるロッドの確認が困難な場合は理想データにてロッドの位置を確認する。プロファイルカーブにどのような違いがあるか？FWHM はどのような差があるか？

0.3

0.4

0.6

0.7

0.5

[cycles/cm]- 21 -

0.3

0.4

0.6

0.7

0.3

0.4

0.6

0.7

0.5

[cycles/cm]

0.5[cycles/cm])

４．６ 再構成法（FBP と OSEM の違い）
目的：
目的： SPECT の再構成アルゴリズムは、解析的手法と逐次近似的手法にわけられる。現在、臨床で一般的に用いられているのは、解析的手法として filtered back projection(FBP)、逐次近似的手法と- 22 -

して ordered subsets expectation maximization(OSEM)であるが、再構成アルゴリズムの違いにより画質は大きく変化する。FBP と OSEM の違いによる画像の変化を視覚的に確認するとともに、定量的に比較検討する。
使用データ
使用データ：
データ： HR_150_4mm_3d_100_Total_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件：
再構成条件： 前処理フィルタ BW order 8, 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
OSEM イタレーション 3, サブセット 10（組み合わせは任意）
散乱、減弱、分解能補正なし

OSEM

FBP

評価方法：
評価方法：
視覚評価：
視覚評価： OSEM と FBP の画像の違いを視覚的に確認する。どのような違いがあるか？
定量評価：
定量評価
・各サイズのロッドと同じ大きさの ROI を設定（理想データを参照する）し、ROI 内のカウントを比較検討する。FBP と OSEM のカウントにはどのような違いがあるか？
・均一性領域にファントム直径の 80％程度の大きさの ROI を設定し、ROI 内のカウントの標準偏差（SD）または変動係数（CV）を算出する。FBP と OSEM のカウントの標準偏差または変動係数はどのような違いがあるか？

OSEM

FBP

- 23 -

OSEM

FBP

４．７ 再構成法（OSEM パラメータの違い）
目的：
目的： OSEM において設定の必要なパラメータにはイタレーションとサブセットがあり、この組み合わせの違いにより画質は大きく変化する。イタレーションとサブセットの違いによる画像の変化を視覚的に確認するとともに、最適な組み合わせについて検討する。
使用データ
使用データ：
データ： HR_150_4mm_3d_100_Total_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成：
再構成： 前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
OSEM イタレーション（It）1，2，3，5，10，20
サブセット（S）1，3，5，10，20，30，60，120 の組み合わせ
散乱、減弱、分解能補正なし

It=3 and S=1

It=3 and S=3- 24 -

It=3 and S=5

It=3 and S=10

It=3 and S=20

It=3 and S=60

It=3 and S=120

It=3 and S=30

評価方法
視覚評価：
視覚評価： パラメータの違いによる画像の変化を視覚的に確認する。イタレーションを固定してサブセットを大きくしていった場合、画像はどのように変化するか？ また、サブセットを固定してイタレー
ションを増やしていった場合に画像はどのように変化するか？
定量評価：

・各サイズのロッドと同じ大きさの ROI を設定（理想データを参照する）し、ROI 内のカウントを比較検討する。パラメータの違いによりカウントはどのように変化するか？
・均一性領域にファントム直径の 80％程度の大きさの ROI を設定し、ROI 内のカウントの標準偏差（SD）または変動係数（CV）を算出する。パラメータの違いにより標準偏差または変動係数はどのように変化するか？
また、最適な組み合わせを見つけることは可能か？

４．８ 散乱補正目的：
目的： 散乱補正により再構成画像のコントラストは上昇する。TEW 法による散乱補正の効果を視覚的に確認するとともに、sub window にかけるバターワースフィルタのカットオフ周波数を変化させ、適正な補正が行われているかを定量的に検討する。
使用データ
使用データ：
データ： HR_150_4mm_3d_100_Total_M
HR_150_4mm_3d_100_True _M
HR_150_4mm_3d_100_ LowSub _M
Ideal_4mm- 25 -

使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件：
再構成条件： 前処理フィルタ BW order 8, 0.5 [cycles/cm] (0.2 cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱補正法：TEW 法
TEW 法ではサブウインドウデータに対し、前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 [cycles/pixel])
減弱、分解能補正なし

true

補正なし

sub window cutoff freq. 0

sub window cutoff freq.0.03

sub window cutoff freq. 005- 26 -

評価方法：
評価方法： 散乱補正の効果を視覚的に確認
ROI テンプレートを用いて ROI 内カウント値を測定し True と比較

４．９ 散乱線（コリメータによる違い）
目的：
目的： コリメータの孔径や隔壁の厚さ等、寸法の違いにより散乱線量が異なり、核医学画像自体の画質に密接に関係する。使用データは，散乱線領域の“_scat_M”を用いて処理することで，コリメー
タの違いによる散乱線量の変化を視覚評価およびプロファイルカーブにて確認する。
使用データ
使用データ：
データ GP_200_4mm_3d_100_scat_M
HR_200_4mm_3d_100_scat_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件
構成条件：
条件： 前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱、減弱、分解能補正なし

HR_200_4mm_3d_100_
HR_200_4mm_3d_100_scat_200_4mm_3d_100_scat_M
scat_M

HR_200_4mm_3d_100_
HR_200_4mm_3d_100_scat_200_4mm_3d_100_scat_M
scat_M

GP_200_4mm_3d_100_
GP_200_4mm_3d_100_scat_200_4mm_3d_100_scat_M
scat_M

GP_200_4mm
GP_200_4mm_3d_100__200_4mm_3d_100_scat_3d_100_scat_M
scat_M- 27 -

評価方法：
評価方法：
視覚評価：
視覚評価：コリメータ別の画像を良く見比べて，カウント差及び空間分解能等による画像劣化を評価する。コリメータが散乱線をより多く含んでいるかどうかを確認する。
定量評価：
定量評価：均一性領域にて減弱補正を行わないで、理想データと比較参照する。本来，散乱線等は場所依存性が無いとされているが，プロファイルカーブによる評価ではどうなのかを確認する。

４．１０ 減弱補正（均一性）
目的：
目的： SPECT 収集時には，散乱（カウントが増加）と減弱（カウントが減少）の２つの相反する現象がある。それらの現象の補正の組み合わせでの画像変化とその性質を知る。また，均一性を維持する為の散乱線補正と減弱補正係数との関係を確認する。
使用データ
使用データ：
データ GP_200_4mm_3d_100_Total_M、
GP_200_4mm_3d_100_Low_S
GP_200_4mm_3d_100_True_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 GMS-5500A
再構成スライス
再構成スライス：スライス 7～55
再構成条件
再構成条件：
条件： 前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱補正：TEW 法
減弱補正： Chang（μ=0.10[cm-1]（散乱補正無）、μ=0.15[cm-1]（散乱補正有）分解能補正なし

散乱線補正無（μ＝0.10[cm-1]）

散乱線補正有（μ＝0.15[cm-1]）（左上段：total 画像、右上段：減弱係数マップ、左下段：減弱補正後画像）- 28 -

profile curve
300

450
400
350
300
250
200
150
100
50
0

250

μ=0.05μ=0.10μ=0.15average countaverage countprofile curve200
μ=0.05μ=0.10μ=0.15

150
100
50
0

0

20

40

60

0

pixel

10

20

30

40

50

60

70

pixel

GP_200_4mm_3d_100_True_M(散乱線無)GP_200_4mm_3d_100_Total_M(散乱線有)profile curve
350

average count

300
250

μ=0.05μ=0.10μ=0.13μ=0.15

200
150
100
50
0
0

20

40
pixel

60

80

GP_200_4mm_3d_100_Total_M
GP_200_4mm_3d_100_Low_S (散乱線補正有)
評価方法：
評価方法：
視覚評価：
視覚評価：散乱線の有無による画像と減弱係数を変化させた場合の画像を比較する。特に，画像の中心部及び辺縁部の変化を確認する。散乱線補正と減弱係数の組み合わせからプロファイル曲
線が均一となる画像を見つける。
定量的評価：
定量的評価：理想データの再構成画像のプロファイルカーブを基本に，それぞれの組み合わせでのプロファイル曲線を比較検討して，最良の組み合わせ（散乱線補正の有無＋減弱係数値）を算出してみる。 最良の組み合わせは，プロファイル曲線がほぼ平坦となることを確認する。

４．１１ 減弱補正（コントラスト）
目的：
目的： 均一性での評価は前項で行ったが，実際の SPECT 画像での散乱線補正と減弱補正による画像の影響について調べる。画像の評価指標であるコントラストの変化について検討を加える。また，
減弱及び散乱線補正の長所・短所を学ぶ。
使用データ
使用データ：
データ GP_200_4mm_3d_100_Total_M、
GP_200_4mm_3d_100_True_M
Ideal_4mm
使用データ
使用データ処理装置
データ処理装置：
処理装置 Windows PC に Prominence Processor ソフトを搭載
再構成スライス
再構成スライス：スライス 7～55
再構成条件
再構成条件：
条件： 前処理フィルタ BW order 8, cutoff freq. 0.5 [cycles/cm] (0.2 [cycles/pixel])
再構成法 FBP （再構成フィルタ ramp）
散乱補正：TEW 法- 29 -

減弱補正： Chang（μ=0.10[cm-1]（散乱補正無）、μ=0.15[cm-1]（散乱補正有）分解能補正なし
1000
900
800

カウント数

700
600
散乱線補正有散乱線補正無

500
400
300
200
100
0
1

2
ROI

3

散乱線補正無散乱線補正有

評価方法：
評価方法：
視覚評価：
視覚評価：特に，散乱線補正及び減弱係数の変化によるカウントの変化を見る。また，前項②での評価も併せて行う。
定量的評価：
定量的評価：各ロッドにターゲットの大きさの 80％ROI を設定し、その ROI 中のカウントを比較する。
コントラストの式にてコントラストを算出し，散乱・減弱補正によりどのようにコントラストが変化するかを評価する。（例：縦軸にカウント，横軸に減弱係数をとり，そのグラフの傾き及びｙ切片等で評価を行う。）減弱係数及びコリメータを変化させた場合の変化も確認すると良い。

- 30 -