(57)【要約】【課題】頭部内の位置によらずなめらかな輪郭を設定することができる核医学診断装置を提供する。【解決手段】本発明の一実施形態に係る核医学診断装置は、頭部の核医学投影データにもとづいて、スライスごとにサイノグラムの左右端を抽出する輪郭抽出部３４と、スライスごとに、輪郭抽出部により抽出されたサイノグラムの左右端の曲線の粗さに応じた精度で曲線をカーブフィッティングしてスライスごとに輪郭を設定する輪郭設定部３５と、輪郭設定部により設定された輪郭を有する減弱マップを生成する減弱マップ生成部３６と、を備えたものである。【選択図】  図２
2015152494.tif 000002
【特許請求の範囲】
【請求項１】
頭部の核医学投影データにもとづいて、スライスごとにサイノグラムの左右端を抽出する輪郭抽出部と、
スライスごとに、前記輪郭抽出部により抽出されたサイノグラムの左右端の曲線の粗さに応じた精度で前記曲線をカーブフィッティングしてスライスごとに輪郭を設定する輪郭設定部と、
前記輪郭設定部により設定された輪郭を有する減弱マップを生成する減弱マップ生成部と、
を備えた核医学診断装置。
【請求項２】
前記輪郭設定部は、
スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を取得し、前記大脳部領域に属するスライスに対する前記カーブフィッティングの精度より粗い精度で前記頭頂部領域および前記小脳部領域に属するスライスに対する前記カーブフィッティングを行う、
請求項１記載の核医学診断装置。
【請求項３】
前記輪郭設定部は、前記カーブフィッティングとしてフーリエ近似を用い、
前記カーブフィッティングの精度の精粗は、前記フーリエ近似の次数の高低である、
請求項１または２記載の核医学診断装置。
【請求項４】
前記輪郭設定部は、
前記頭頂部領域に属するスライスのうち、サイノグラムの左右端間の平均距離が第１距離閾値以下であり、かつサイノグラムの左右端間の距離の標準偏差が第１標準偏差閾値以上であるスライスに対しては前記大脳部領域に属するスライスに対する前記カーブフィッティングの精度より粗い精度で前記カーブフィッティングを行う一方、前記平均距離が前記第１距離閾値より大きいまたは前記標準偏差が前記第１標準偏差閾値より小さいスライスに対しては前記大脳部領域に属するスライスと同一の精度で前記カーブフィッティングを行う、
請求項２または３に記載の核医学診断装置。
【請求項５】
前記輪郭設定部は、
頭頂部から小脳部に向かうスライス順にサイノグラムを抽出し、抽出したサイノグラムごとに左右端間の平均距離を求め、スライス間の前記平均距離の変化率が変化率閾値以下となるスライスから前記頭頂部までのスライスを前記頭頂部領域に属するスライスとする、
請求項２ないし４のいずれか１項に記載の核医学診断装置。
【請求項６】
前記輪郭設定部は、
小脳部の下端から頭頂部に向かうスライス順にサイノグラムを抽出し、抽出したサイノグラムごとに左右端間の平均距離および左右端間の距離の標準偏差を求め、前記平均距離が第２距離閾値以下であるか、または前記標準偏差が第２標準偏差閾値以上であるスライスを前記小脳部領域に属するスライスとする、
請求項２ないし５のいずれか１項に記載の核医学診断装置。
【請求項７】
前記輪郭設定部は、
ユーザによる入力部を介した指示にもとづいて、前記スライスのそれぞれが前記頭頂部領域、前記小脳部領域および前記大脳部領域のいずれに属するかの情報を取得する、
請求項２ないし４のいずれか１項に記載の核医学診断装置。
【請求項８】
前記輪郭設定部は、
ユーザによる入力部を介した指示にもとづいて、前記輪郭抽出部により抽出されたサイノグラムの左右端の前記曲線の粗さに応じた前記カーブフィッティングの精度を設定する、
請求項１ないし７のいずれか１項に記載の核医学診断装置。
【請求項９】
前記輪郭抽出部は、
ユーザによる入力部を介した所定の割合の指示を受け、 各スライスのサイノグラムについて最大計数値に対して前記所定の割合以下の計数値が割り当てられた画素をカットする画像処理により、スライスごとにサイノグラムの左右端を抽出する、
請求項１ないし８のいずれか１項に記載の核医学診断装置。
【請求項１０】
前記頭部の核医学投影データにもとづいて、前記減弱マップを用いて減弱補正処理を施した再構成画像を生成し表示部に表示させる再構成部、
をさらに備えた請求項１ないし９のいずれか１項に記載の核医学診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、核医学診断装置に関する。
【背景技術】
【０００２】
核医学診断装置は、放射性同位元素（Radio Isotope、以下ＲＩという）を含む薬品（血流マーカ、トレーサ）が生体内の特定組織や臓器に選択的に取り込まれる性質を利用して、生体内に分布したＲＩから放射されるガンマ線を、生体外に配設されたガンマ線の検出器で検出するようになっている。ガンマ線の検出結果は、ガンマ線の線量分布を画像化することによる核医学画像の生成や、体内臓器等の機能の診断などに利用される。
【０００３】
ガンマ線は、生体内で減弱される。このため、ガンマ線の検出結果には、この生体内での減弱の影響が含まれている。この種の生体内での減弱の影響を補正する方法としてよく用いられる方法に、使用する核種のガンマ線エネルギーの減弱係数の分布を示す減弱計数マップ（以下、減弱マップという）を生成し、この減弱マップにもとづいてガンマ線の検出結果を補正する方法がある。この方法によれば、ガンマ線の生体内での減弱の影響を補正する（以下、減弱補正するという）ことができる。このため、減弱の補正を行わない場合にくらべ、核医学画像をより高精度に生成することができる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−８１６４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、減弱マップの輪郭形状と補正対象となる核医学画像の輪郭形状とが適切な対応関係にないと、減弱補正後の核医学画像にもとづく診断を正確に行うことが難しくなってしまう。たとえば、頭部の核医学画像を減弱補正する場合では、減弱マップの輪郭形状に応じて、減弱補正後の核医学画像にもとづいて求められる脳血流（ＣＢＦ、Cerebral Blood Flow）の値が変化してしまうことが知られている。具体的には、減弱マップの輪郭が核医学画像内の頭部の外縁よりも大きすぎると、吸収を過補正してしまい脳血流値の過大評価につながるおそれがある一方、減弱マップの輪郭が核医学画像内の頭部の外縁よりも小さすぎると、吸収が十分に補正できず脳血流の過小評価につながるとともにアーティファクトの原因となってしまう。
【０００６】
このため、減弱マップの輪郭形状は、減弱補正後の核医学画像にもとづく病理診断における定量評価の精度や誤診の可能性に影響する。したがって、減弱マップは、適切かつなめらかな輪郭形状を有するように生成することが重要である。
【０００７】
スライスごとに減弱マップの輪郭を設定する方法として、たとえば投影データからスライスごとにサイノグラムの左右端を仮に自動抽出し、この抽出した左右端の曲線を所定の精度でカーブフィッティングすることによってスライスごとになめらかな輪郭を設定する方法がある。しかし、各スライスのカウント数は、頭部内の位置に応じて異なる。このため、この方法では、所定の精度に適合したスライスの輪郭はなめらかにフィッティング可能である一方、他の位置のスライスの輪郭は粗くガタついてしまう。
【課題を解決するための手段】
【０００８】
本発明の一実施形態に係る核医学診断装置は、上述した課題を解決するために、頭部の核医学投影データにもとづいて、スライスごとにサイノグラムの左右端を抽出する輪郭抽出部と、スライスごとに、輪郭抽出部により抽出されたサイノグラムの左右端の曲線の粗さに応じた精度で曲線をカーブフィッティングしてスライスごとに輪郭を設定する輪郭設定部と、輪郭設定部により設定された輪郭を有する減弱マップを生成する減弱マップ生成部と、を備えたものである。
【図面の簡単な説明】
【０００９】
【図１】本発明の一実施形態に係る核医学診断装置の一例を示すブロック図。
【図２】制御部のＣＰＵ等の演算装置による機能実現部の構成例を示す概略的なブロック図。
【図３】（ａ）はサイノグラムの生データの一例を示す説明図、（ｂ）は輪郭抽出部により自動抽出されたサイノグラムの左右端（自動抽出左右端）の一例を示す説明図、（ｃ）は輪郭設定部により得られるサイノグラムの左右端のフィッティングカーブ（近似左右端）の一例を示す説明図、（ｄ）は輪郭設定部により近似左右端にもとづいて輪郭を設定する方法を説明するための図、（ｅ）は輪郭設定部により近似左右端にもとづいて設定された輪郭の一例を示す説明図。
【図４】頭部のサイノグラムの一例を示す説明図。
【図５】（ａ）は頭頂部付近スライスのサイノグラムの自動抽出左右端の一例を示す図、（ｂ）は（ａ）に示す自動抽出左右端に対して高次のフーリエ近似により求められた近似左右端の一例を示す説明図、（ｃ）は（ａ）に示す自動抽出左右端に対して低次のフーリエ近似により求められた近似左右端の一例を示す説明図。
【図６】全てのスライスについて高次のフーリエ近似により求められた近似左右端にもとづいて設定された輪郭と減弱補正前の核医学画像（再構成画像）との重畳画像の一例を示す説明図。
【図７】大脳部領域に属するスライスに対しては高次のフーリエ近似を行う一方、頭頂部領域に属するスライスおよび小脳部領域に属するスライスに対しては低次のフーリエ近似を行うことにより求められた近似左右端にもとづいて設定された輪郭と、減弱補正前の核医学画像（再構成画像）との重畳画像の一例を示す説明図。
【図８】図１に示す制御部により、頭部内の位置によらずなめらかな輪郭を設定する際の手順の一例を示すフローチャート。
【図９】図８のステップＳ４で輪郭設定部により実行される輪郭設定処理の手順を示すサブルーチンフローチャート。
【図１０】図９のステップＳ４１で輪郭設定部により実行される、頭頂部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャート。
【図１１】図９のステップＳ４２で輪郭設定部により実行される、小脳部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャート。
【図１２】図９のステップＳ４３で輪郭設定部により実行される、大脳部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャート。
【図１３】スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を、ユーザに設定させるための画像の一例を示す説明図。
【図１４】頭部内の位置によらずなめらかな輪郭を設定する際の他の手順の一例を示すフローチャート。
【図１５】図１４のステップＳ１０３で輪郭設定部により実行される輪郭設定処理の手順を示すサブルーチンフローチャート。
【発明を実施するための形態】
【００１０】
本発明に係る核医学診断装置の実施の形態について、添付図面を参照して説明する。本実施形態に係る核医学診断装置は、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）やＰＥＴ（Positron Emission Tomography）などのガンマ線検出器を備えた各種装置に適用することが可能である。以下の説明では、本発明に係る核医学診断装置として２検出器型のガンマ線検出器回転型のＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置を用いる場合の一例について示す。なお、ガンマ線検出器回転型ＳＰＥＣＴ装置としては、検出器が１つまたは３以上のものであってもよい。また、本実施形態に係る核医学診断装置の画像処理は被検体の様々な部位の核医学投影データに対して適用可能である。以下の説明では被検体の頭部の核医学投影データに対して適用する場合の例について示す。
【００１１】
図１は、本発明の一実施形態に係る核医学診断装置１の一例を示すブロック図である。
【００１２】
核医学診断装置１は、スキャナ装置２および画像処理装置３を有する。なお、画像処理装置３はスキャナ装置２とデータ送受信可能に接続されていればよく、同一の部屋や建屋に設けられずともよい。
【００１３】
スキャナ装置２は、ガンマ線検出部１１および１２、回転部１３を有するガントリ１４、回転駆動装置１５、データ収集部１６を有する。
【００１４】
ガンマ線検出部１１は、被検体Ｐに対して薬品に含まれて投与されたテクネシウムなどのＲＩ（放射性同位元素）から放射されるガンマ線を検出する検出部である。なお、ガンマ線検出部１２はガンマ線検出部１１と同様の構成および作用を有するため、説明を省略する。
【００１５】
ガンマ線検出部１１としては、シンチレータ型検出器を用いてもよいし、半導体型検出器を用いてもよい。
【００１６】
シンチレータ型検出器を用いてガンマ線検出部１１を構成する場合は、ガンマ線検出部１１は、ガンマ線の入射角度を規定するためのコリメータ、コリメートされたガンマ線が入射すると瞬間的な閃光を発するシンチレータ、ライトガイド、シンチレータから射出された光を検出するための２次元に配列された複数の光電子増倍管、およびシンチレータ用電子回路などを有する。シンチレータは、たとえばタリウム活性化ヨウ化ナトリウムＮａＩ（Ｔｌ）により構成される。
【００１７】
シンチレータ用電子回路は、ガンマ線が入射する事象（イベント）が発生するごとに、複数の光電子増倍管の出力にもとづいて複数の光電子増倍管により構成される検出面内におけるガンマ線の入射位置情報（位置情報）、入射強度情報および入射時刻情報を生成しデータ収集部１６に出力する。この位置情報は、検出面内の２次元座標の情報であってもよいし、あらかじめ検出面を複数の分割領域（１次セル）に仮想的に分割しておき（たとえば１０２４×１０２４個に分割しておき）、どの１次セルに入射があったかを示す情報であってもよい。
【００１８】
一方、半導体型検出器を用いてガンマ線検出部１１を構成する場合は、ガンマ線検出部１１は、コリメータ、コリメートされたガンマ線を検出するための２次元に配列された複数のガンマ線検出用半導体素子（以下、半導体素子という）および半導体用電子回路などを有する。半導体素子は、たとえばＣｄＴｅやＣｄＺｎＴｅ（ＣＺＴ）により構成される。
【００１９】
半導体用電子回路は、ガンマ線が入射する事象（イベント）が発生するごとに、半導体素子の出力にもとづいて入射位置情報、入射強度情報および入射時刻情報を生成しデータ収集部１６に出力する。この位置情報は、複数の半導体素子（たとえば１０２４×１０２４個）のうちのどの半導体素子に入射したかを示す情報である。
【００２０】
すなわち、ガンマ線検出部１１は、イベントごとに入射位置情報、入射強度情報および入射時刻情報を出力する。また、位置情報は、１次セルのどの位置にガンマ線が入射したかを示す情報および検出面内の２次元座標の情報の少なくとも一方である。以下の説明では、ガンマ線検出部１１が位置情報として検出面内のどの位置にガンマ線が入射したかを示す情報を出力する場合の例について示す。
【００２１】
ガンマ線検出部１１および１２は、データ収集部１６により撮像タイミングを制御される。
【００２２】
回転部１３はガントリ１４に支持され、ガンマ線検出部１１および１２を保持する。回転部１３が回転駆動装置１５を介してデータ収集部１６に制御されて所定の回転軸ｒ周り（ｚ軸（体軸）周り）に回転することにより、ガンマ線検出部１１および１２は回転軸ｒの周りを回転する。
【００２３】
回転駆動装置１５は、データ収集部１６に制御されて、回転部１３を所定の回転軸ｒ（ｚ軸）の周りに回転させる。
【００２４】
被検体Ｐは、天板１７に載置される。天板駆動装置１８は、データ収集部１６に制御されて、天板１７をｙ軸方向に昇降動させる。また、天板駆動装置１８は、データ収集部１６に制御されて、回転部１３の中央部分の開口部の撮影領域へ天板１７をｚ軸方向に沿って移送する。
【００２５】
データ収集部１６は、画像処理装置３により制御されて、ガンマ線検出部１１および１２、回転駆動装置１５および天板駆動装置１８を制御することにより、被検体Ｐのスキャンを実行する。
【００２６】
データ収集部１６は、ガンマ線検出部１１および１２のそれぞれの出力をたとえばリストモードで収集する。リストモードでは、ガンマ線の検出位置情報、強度情報、ガンマ線検出部１１および１２と被検体Ｐとの相対位置を示す情報（ガンマ線検出部１１および１２の位置や角度など）、およびガンマ線の検出時刻がガンマ線の入射イベントごとに収集される。
【００２７】
一方、画像処理装置３は、図１に示すように、制御部２１、表示部２２、入力部２３および記憶部２４を有する。
【００２８】
制御部２１は、ＣＰＵ等の演算装置、ＲＡＭおよびＲＯＭをはじめとする記憶媒体などにより構成され、この記憶媒体に記憶されたプログラムに従って画像処理装置３の動作を制御する。
【００２９】
制御部２１のＣＰＵ等の演算装置は、ＲＯＭをはじめとする記憶媒体に記憶された画像処理プログラムおよびこのプログラムの実行のために必要なデータをＲＡＭへロードし、このプログラムに従って、頭部内の位置によらずなめらかな輪郭を設定するための処理を実行する。
【００３０】
制御部２１のＲＡＭは、ＣＰＵ等の演算装置が実行するプログラムおよびデータを一時的に格納するワークエリアを提供する。制御部２１のＲＯＭをはじめとする記憶媒体は、画像処理装置３の起動プログラム、画像処理プログラムや、これらのプログラムを実行するために必要な各種データを記憶する。
【００３１】
なお、ＲＯＭをはじめとする記憶媒体は、磁気的もしくは光学的記録媒体または半導体メモリなどの、ＣＰＵ等の演算装置により読み取り可能な記録媒体を含んだ構成を有し、これら記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介してダウンロードされるように構成してもよい。
【００３２】
表示部２２は、たとえば液晶ディスプレイやＯＬＥＤ（Organic Light Emitting Diode）ディスプレイなどの一般的な表示出力装置により構成され、制御部２１の制御に従って核医学診断画像などの各種情報を表示する。
【００３３】
入力部２３は、たとえばキーボード、タッチパネル、テンキーなどの一般的な入力装置により構成され、ユーザの操作に対応した操作入力信号を制御部２１に出力する。
【００３４】
記憶部２４は、磁気的もしくは光学的記録媒体または半導体メモリなどの、ＣＰＵ等の演算装置により読み取り可能な記録媒体を含んだ構成を有し、これら記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介してダウンロードされるように構成してもよい。記憶部２４は、制御部２１により制御されて表示画素ごとの計数値や核医学画像、変化率閾値ｔｈΔをはじめとする各種閾値などを記憶する。
【００３５】
図２は、制御部２１のＣＰＵ等の演算装置による機能実現部の構成例を示す概略的なブロック図である。なお、この機能実現部は、ＣＰＵ等の演算装置を用いることなく回路などのハードウエアロジックによって構成してもよい。
【００３６】
図２に示すように、制御部２１のＣＰＵ等の演算装置は、画像処理プログラムによって、少なくともスキャン制御部３１、ファンパラ変換部３２、前処理部３３、輪郭抽出部３４、輪郭設定部３５、減弱マップ生成部３６および再構成部３７として機能する。この各部３１〜３８は、ＲＡＭの所要のワークエリアを、データの一時的な格納場所として利用する。なお、ファンパラ変換部３２は、パラレルビームコリメータを用いた場合は処理不要である。
【００３７】
スキャン制御部３１は、ユーザから入力部２３を介してスキャン計画の実行指示を受けて、スキャン計画にもとづいてスキャナ装置２を制御することにより、被検体Ｐの頭部のスキャンを実行する。この結果、被検体Ｐの頭部から放出されたガンマ線の情報（核医学投影データ）がスキャナ装置２からデータ収集部１６を介してファンパラ変換部３２に与えられる。
【００３８】
ガンマ線検出部１１および１２がファンビームコリメータを有しファンビーム投影データを収集する場合、ファンパラ変換部３２は、データ収集部１６から受けたファンビーム投影データをパラレルビーム投影データに変換する。一方、ガンマ線検出部１１および１２がパラレルビーム投影データを収集するよう構成されている場合は、ファンパラ変換部３２は、データ収集部１６から受けたパラレルビーム投影データをそのまま前処理部３３に与える。
【００３９】
前処理部３３は、入力部２３を介してユーザにより設定された、またはあらかじめ記憶部２４に記憶された前処理条件情報を取得する。前処理条件情報には、前処理の対象となる画像データ（パラレルデータ）の情報、処理範囲の設定情報、前処理フィルターの情報などが含まれる。前処理部３３は、取得した前処理条件情報に従って、均一性補正、回転中心補正、前処理フィルター処理などの前処理を実行する。
【００４０】
図３（ａ）はサイノグラムの生データの一例を示す説明図であり、（ｂ）は輪郭抽出部３４により自動抽出されたサイノグラムの左右端（自動抽出左右端）４１の一例を示す説明図であり、（ｃ）は輪郭設定部３５により得られるサイノグラムの左右端のフィッティングカーブ（近似左右端）４２の一例を示す説明図である。さらに、図３（ｄ）は輪郭設定部３５により近似左右端４２にもとづいて輪郭４３を設定する方法を説明するための図であり、（ｅ）は輪郭設定部３５により近似左右端４２にもとづいて設定された輪郭４３の一例を示す説明図である。
【００４１】
輪郭抽出部３４は、頭部の核医学投影データにもとづいて、スライスごとにサイノグラムの左右端を自動抽出する（図３（ａ）、（ｂ）参照）。
【００４２】
この自動抽出左右端４１を抽出する方法としては、たとえば頭部の核医学投影データの最大計数値に対して所定の割合（たとえば３％、５％、１０％、１５％など）以下の計数値が割り当てられた画素をカットする画像処理により、サイノグラムの自動抽出左右端４１を抽出する方法を用いることができる。
【００４３】
また、再構成部３７は、前処理部３３による前処理後の画像データに含まれる入射位置情報にもとづいて、表示画素のそれぞれについて入射ガンマ線を光子数として計数した計数値を割り当てる（分配する）ことにより、被検体Ｐの頭部の減弱補正前の核医学画像を生成する。輪郭抽出部３４は、この核医学画像に対して上記計数値にもとづく画像処理を行うことにより頭部の輪郭を抽出し、抽出した輪郭に対応するサイノグラムの左右端を自動抽出左右端４１としてもよい。
【００４４】
また、輪郭抽出部３４は、ユーザによって入力部２３を介した指示に応じて上記計数値にもとづく画像処理における所定の割合（左右端の自動抽出の閾値）を設定し、あるいは変更してもよい。また、所定の割合はスライスごとに異なって設定されてもよい。
【００４５】
次に、輪郭設定部３５は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線の粗さに応じた精度で、左右端の曲線をカーブフィッティングすることにより、スライスごとに近似左右端（フィッティングカーブ）４２を求める（図３（ｃ）参照）。
【００４６】
たとえば、輪郭設定部３５は、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を取得し、大脳部領域に属するスライスに対するカーブフィッティングの精度より粗い精度で頭頂部領域に属するスライスの一部または全部および小脳部領域に属するスライスの一部または全部に対するカーブフィッティングを行う。
【００４７】
カーブフィッティング法としては、従来種々のものが知られており、これらのうち任意のものを使用することが可能である。
【００４８】
以下の説明では、カーブフィッティング法としてフーリエ近似法を用いる場合の例について示す。フーリエ近似法を用いる場合、カーブフィッティングの精度の精粗は、フーリエ近似の次数の高低である。すなわち、大脳部領域に属するスライスに対しては高次（たとえば５〜１１次）のフーリエ近似を行う一方、頭頂部領域に属するスライスの一部または全部および小脳部領域に属するスライスの一部または全部に対しては低次（たとえば２〜５次）のフーリエ近似を行う。
【００４９】
なお、高次、低次の次数はそれぞれあらかじめ設定された次数を用いてもよいし、ユーザにより入力部２３を介して各次数を設定されてもよい。また、頭頂部領域と小脳部領域でそれぞれ異なる次数を用いてもよい。また、頭頂部領域を複数に分割し、頭頂部に近づくほど次数が下がるようにしてもよい。同様に、小脳部領域を複数に分割し、小脳部下端に近づくほど次数が下がるようにしてもよい。
【００５０】
また、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報は、輪郭設定部３５が自らサイノグラムデータにもとづいて判定して取得してもよいし、ユーザによる入力部２３を介した指示にもとづいて取得してもよい。
【００５１】
輪郭設定部３５が自らサイノグラムデータにもとづいて各スライスがどの領域に属するかを判定する方法としては、たとえばサイノグラムの左右端の間の距離Ｌ（θ）の平均ＭＬや距離Ｌ（θ）の標準偏差ＳＤＬについて閾値判定する方法などが挙げられる。この自動判定方法の詳細については図１０−１２を用いて後述する。
【００５２】
そして、輪郭設定部３５は、サイノグラムの近似左右端４２のデータのみをバックプロジェクションするように投影方向に向かう直線を考え、次の角度θの直線との交点を輪郭の構成点を求めていく（図３（ｄ）参照）。１８０度分のデータに対応する輪郭の構成点を求めることにより、輪郭設定部３５はスライスごとに輪郭４３を設定することができる（図３（ｅ）参照）。
【００５３】
輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接、この曲線の粗さに応じた精度で左右端の曲線をカーブフィッティングすることにより、スライスごとに近似左右端４２を求めてもよい。
【００５４】
減弱マップ生成部３６は、たとえば水などの均一吸収体を仮定して、輪郭設定部３５により設定された輪郭４３を有する減弱マップ（使用する核種のガンマ線エネルギーの減弱係数の分布を示す減弱係数マップ）を生成する。
【００５５】
再構成部３７は、前処理部３３による前処理後の画像データに含まれる入射位置情報にもとづいて、表示画素のそれぞれについて入射ガンマ線を光子数として計数した計数値を割り当てる（分配する）ことにより、被検体Ｐの頭部の減弱補正前の核医学画像を生成する。再構成部３７は、この減弱補正前の核医学画像（再構成画像）に対して輪郭設定部３５により設定された輪郭４３を重畳した重畳画像を表示部２２に表示させる。重畳画像には、輪郭４３とともに、あるいは輪郭４３に変えて、減弱マップを重畳してもよい。また、減弱マップを重畳する場合には、核医学画像についてはファンパラ変換前にＴＥＷ（Triple Energy Window）処理により散乱線補正を行ったものを用いる一方、減弱マップについてはファンパラ変換前にＴＥＷ処理を行わないものを用いるようにしてもよい。
【００５６】
ユーザは、この重畳画像を視認することにより、輪郭設定部３５により設定された輪郭４３が適切か否かを判断することができる。また、ユーザは、輪郭設定部３５により設定された輪郭４３が不適切であると考える場合は、入力部２３を介して輪郭設定部３５に対して輪郭４３の修正を指示することができる。このときユーザはさらに、輪郭設定部３５によるカーブフィッティングの精度をスライスごとに新たに設定することができる。輪郭４３の修正が指示された場合は、輪郭設定部３５は、たとえばユーザにより新たに設定されたカーブフィッティングの精度（たとえばフーリエ近似であれば新たな次数）を用いて近似左右端４２を再度求めて輪郭４３を再設定する。
【００５７】
また、輪郭設定部３５は、特に頭頂部領域や小脳部領域に属するスライスについて、１つのスライスに対して複数の次数でフーリエ近似を行って複数の輪郭４３を生成してもよい。この場合、再構成部３７は、再構成画像に対して複数の輪郭４３を重畳して表示部２２に表示させる。この場合、ユーザは複数の輪郭４３から所望の輪郭４３を選択することができる。また、減弱マップ生成部３６は、ユーザから複数の輪郭４３から選択された輪郭４３を有する減弱マップを生成する。
【００５８】
なお、ユーザによる輪郭４３の適切性の確認を省略してもよい場合は、再構成部３７は重畳画像を生成せずともよい。
【００５９】
そして、再構成部３７は、減弱マップにもとづいてガンマ線の検出結果を補正する（減弱補正処理を施す）ことによりガンマ線の生体内での減弱の影響を補正した画像（減弱補正した核医学画像（再構成画像））を生成し、表示部２２に表示させる。
【００６０】
続いて、輪郭設定部３５による近似左右端４２の算出処理について説明する。
【００６１】
図４は、頭部のサイノグラムの一例を示す説明図である。また、図５（ａ）は頭頂部付近スライスのサイノグラムの自動抽出左右端４１の一例を示す図であり、（ｂ）は（ａ）に示す自動抽出左右端４１に対して高次のフーリエ近似により求められた近似左右端４２の一例を示す説明図であり、（ｃ）は（ａ）に示す自動抽出左右端４１に対して低次のフーリエ近似により求められた近似左右端４２の一例を示す説明図である。
【００６２】
また、図６は、全てのスライスについて高次のフーリエ近似により求められた近似左右端４２にもとづいて設定された輪郭４３と減弱補正前の核医学画像（再構成画像）との重畳画像の一例を示す説明図である。また、図７は、大脳部領域に属するスライスに対しては高次のフーリエ近似を行う一方、頭頂部領域に属するスライスおよび小脳部領域に属するスライスに対しては低次のフーリエ近似を行うことにより求められた近似左右端４２にもとづいて設定された輪郭４３と、減弱補正前の核医学画像（再構成画像）との重畳画像の一例を示す説明図である。図６および図７は、頭部のアキシャル断面のスライスを頭頂部上端から小脳部下端まで一覧表示する場合の表示部２２の表示画面の一例である。
【００６３】
図４に示すように、大脳部領域に属するスライスの計数値は、頭頂部領域や小脳部領域に属するスライスの計数値よりも大きい。このため、大脳部領域に属するスライスの自動抽出左右端４１の曲線の粗さ（がたつき）は少なくなめらかとなると考えられる。一方、頭頂部領域や小脳部領域に属するスライスの計数値は、頭頂部上端および小脳部下端に近づくほど少なくなっていく。このため、頭頂部上端および小脳部下端に近づくほど、自動抽出左右端４１の曲線の粗さ（がたつき）は大きくなっていくと考えられる。
【００６４】
たとえば、図５（ａ）に示すように、頭頂部上端付近のスライスのサイノグラムの自動抽出左右端４１は、がたつきの大きい曲線となる。このがたつきの大きい自動抽出左右端４１に対して高次のフーリエ近似を行ってしまうと、精度よく曲線あてはめがおこなわれるために、がたつきが再現された近似左右端４２が得られてしまう（図５（ｂ）、図６参照）。
【００６５】
そこで、本実施形態に係る核医学診断装置１の輪郭設定部３５は、計数値が少なく自動抽出左右端４１の曲線が粗い場合は、自動抽出左右端４１に対するカーブフィッティングの精度をあえて粗くする（フーリエ近似の次数を下げる）ことにより、曲線あてはめの追従性を劣化させてなめらかな近似左右端４２を得ることでなめらかな輪郭４３を設定する（図５（ｃ）、図７参照）。
【００６６】
このため、計数値が少なく自動抽出左右端４１の曲線が粗いスライスであってもなめらかな輪郭４３を設定することができる。したがって、単に自動抽出左右端４１にもとづく輪郭を用いる場合に比べ、頭部内の位置によらずなめらかな輪郭４３を設定することができ、計数値が少なく自動抽出左右端４１の曲線が粗いスライスにおける減弱マップの精度を大幅に向上させることができる。
【００６７】
次に、本実施形態に係る核医学診断装置１の動作の一例について説明する。
【００６８】
図８は、図１に示す制御部２１により、頭部内の位置によらずなめらかな輪郭４３を設定する際の手順の一例を示すフローチャートである。図８において、Ｓに数字を付した符号はフローチャートの各ステップを示す。
【００６９】
まず、ステップＳ１において、スキャン制御部３１は、ユーザから入力部２３を介してスキャン計画の実行指示を受けてスキャンを実行する。この結果、被検体Ｐから放出されたガンマ線の情報がスキャナ装置２からデータ収集部１６を介してファンパラ変換部３２に与えられる。そして、ガンマ線検出部１１および１２がファンビームコリメータを有しファンビーム投影データを収集する場合、ファンパラ変換部３２は、データ収集部１６から受けたファンビーム投影データをパラレルビーム投影データに変換して前処理部３３に与える。なお、ガンマ線検出部１１および１２がパラレルビーム投影データを収集するよう構成されている場合は、ファンパラ変換部３２は、データ収集部１６から受けたパラレルビーム投影データをそのまま前処理部３３に与える。
【００７０】
次に、ステップＳ２において、前処理部３３は、入力部２３を介してユーザにより設定された、またはあらかじめ記憶部２４に記憶された前処理条件情報を取得し、前処理条件情報に従って、均一性補正、回転中心補正、前処理フィルター処理などの前処理を実行する。
【００７１】
次に、ステップＳ３において、輪郭抽出部３４は、頭部の核医学投影データにもとづいて、スライスごとに自動抽出左右端４１を抽出する（図３（ａ）、（ｂ）参照）。
【００７２】
次に、ステップＳ４において、輪郭設定部３５は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線の粗さに応じた精度で、左右端の曲線をカーブフィッティングすることにより、スライスごとに近似左右端４２を求めて輪郭４３を設定する（図３（ｅ）、図７参照）。
【００７３】
次に、ステップＳ５において、減弱マップ生成部３６は、たとえば水などの均一吸収体を仮定して、輪郭設定部３５により設定された輪郭４３を有する減弱マップを生成する
【００７４】
次に、ステップＳ６において、再構成部３７は、減弱補正前の再構成画像に対して輪郭設定部３５により設定された輪郭４３を重畳した重畳画像を表示部２２に表示させる。なお、この重畳画像には、輪郭４３とともに、あるいは輪郭４３に変えて、減弱マップを重畳してもよい。
【００７５】
次に、ステップＳ７において、輪郭設定部３５は、輪郭４３を修正すべきか否かを判定する。輪郭４３を修正すべき場合はステップＳ３に戻る。一方、輪郭４３の修正が不要な場合はステップＳ８に進む。なお、このステップＳ７は省略されてもよい。
【００７６】
輪郭４３を修正すべき場合としては、たとえば重畳画像を確認したユーザが、輪郭設定部３５により設定された輪郭４３が不適切であると考える場合であって、入力部２３を介して輪郭設定部３５に対して輪郭４３の修正を指示した場合が挙げられる。このときユーザはさらに、輪郭設定部３５によるフーリエ近似の次数をスライスごとに新たに設定することができる。ステップＳ７でＹＥＳと判定されてステップＳ３に戻る場合、ステップＳ４において輪郭設定部３５は、ステップＳ７であらたに設定された次数で各スライスのフーリエ近似を行う。
【００７７】
また、輪郭設定部３５は、重畳画像における輪郭４３のがたつきの度合いをたとえば近似楕円からの乖離度合いなどにもとづいて画像解析により自動判定してもよい。この場合、がたつき度合いが閾値以上のスライスがある場合には輪郭４３を修正すべき場合と判定してステップＳ３に戻り、ステップＳ４において輪郭設定部３５は、がたつき度合いが閾値以上のスライスについて、前回実行時よりも次数を下げてフーリエ近似を行うとよい。
【００７８】
そして、ステップＳ８において、再構成部３７は、減弱マップにもとづいてガンマ線の検出結果を補正する（減弱補正処理を施す）ことによりガンマ線の生体内での減弱の影響を補正した画像（減弱補正した核医学画像（再構成画像））を生成し、表示部２２に表示させる。
【００７９】
以上の手順により、頭部内の位置によらず、なめらかな輪郭４３を設定することができる。
【００８０】
図９は、図８のステップＳ４で輪郭設定部３５により実行される輪郭設定処理の手順を示すサブルーチンフローチャートである。
【００８１】
ステップＳ４１において、輪郭設定部３５は、頭頂部領域に属するスライスの輪郭４３を設定する。
次に、輪郭設定部３５は、小脳部領域に属するスライスの輪郭４３の設定（ステップＳ４２）、および大脳部領域に属するスライスの輪郭４３の設定（ステップＳ４３）を行い、図８のステップＳ５に進む。
【００８２】
続いて、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報について、輪郭設定部３５が自らサイノグラムデータにもとづいて判定して取得する場合の例について図１０−１２を用いて説明する。
【００８３】
図１０は、図９のステップＳ４１で輪郭設定部３５により実行される、頭頂部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャートである。
【００８４】
ステップＳ４１１において、輪郭設定部３５は、頭頂部上端のスライスから順にサイノグラムを抽出する。
【００８５】
次に、ステップＳ４１２において、輪郭設定部３５は、スライスごとに、サイノグラムの左右端の間の距離Ｌ（θ）を角度ごとに算出する（図３（ｄ）参照）。
【００８６】
次に、ステップＳ４１３において、輪郭設定部３５は、スライスごとに、距離Ｌ（θ）の平均ＭＬを算出する。
【００８７】
次に、ステップＳ４１４において、輪郭設定部３５は、平均距離ＭＬのスライス間の変化率ΔＭＬ／Δｚが変化率閾値ｔｈΔ以上か否かを判定する。ΔＭＬ／Δｚ≧変化率閾値ｔｈΔの場合は頭頂部領域に属するスライスと判断し、ステップＳ４１５に進む。一方、ΔＭＬ／Δｚ＜変化率閾値ｔｈΔの場合はステップＳ４１９に進む。
【００８８】
次に、ステップＳ４１５において、輪郭設定部３５は、頭頂部領域に属すると判断したスライスについて、距離Ｌ（θ）の標準偏差ＳＤＬを算出する。
【００８９】
次に、ステップＳ４１６において、輪郭設定部３５は、平均距離ＭＬが第１距離閾値ｔｈＭＬ１以下であり、かつ標準偏差ＳＤＬが第１標準偏差閾値ｔｈＳＤＬ１以上であるか否かを判定する。ＭＬ≦ｔｈＭＬ１かつＳＤＬ≧ｔｈＳＤＬ１を満たすスライスは、サイノグラムの左右端間の距離が小さく、かつ左右端の曲線のがたつきが大きいスライスであり、頭頂部領域に属するスライスの中でも頭頂部上端に近いスライスであると考えられる。ＭＬ≦ｔｈＭＬ１かつＳＤＬ≧ｔｈＳＤＬ１である場合は、ステップＳ４１７に進む。一方、ＭＬ＞ｔｈＭＬ１またはＳＤＬ＜ｔｈＳＤＬ１である場合は、ステップＳ４１８に進む。
【００９０】
次に、ステップＳ４１７において、輪郭設定部３５は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定してステップＳ４１１に戻る。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してステップＳ４１１に戻ってもよい。
【００９１】
一方、ステップＳ４１８において、ＭＬ＞ｔｈＭＬ１またはＳＤＬ＜ｔｈＳＤＬ１を満たすスライスについては、ＭＬ≦ｔｈＭＬ１かつＳＤＬ≧ｔｈＳＤＬ１をスライスに対する次数よりも高い次数でフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してステップＳ４１１に戻る。
【００９２】
なお、ここでいう高い次数とは、大脳部領域に属するスライスに対して適用される次数と同じであってもよいし、大脳部領域に属するスライスに対して適用される次数よりも低い次数であってもよい。
【００９３】
他方、ステップＳ４１４でΔＭＬ／Δｚ＜変化率閾値ｔｈΔと判定されると、ステップＳ４１９において、輪郭設定部３５は、このスライスについては大脳部領域に属すると判断して頭頂部からのスライス順のサイノグラムの抽出を終了し、図９のステップＳ４２に進む。
【００９４】
なお、頭頂部領域に属するスライスの全てに対して同一次数の低次のフーリエ近似を用いてもよい。この場合、ステップＳ４１６およびＳ４１８は省略される。
【００９５】
以上の手順により、頭頂部領域に属するスライスに対して低次のフーリエ近似を用いて輪郭４３を設定することができる。
【００９６】
図１１は、図９のステップＳ４２で輪郭設定部３５により実行される、小脳部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャートである。
【００９７】
ステップＳ４２１において、輪郭設定部３５は、小脳部下端のスライスから逆順にサイノグラムを抽出する。
【００９８】
次に、ステップＳ４２２において、輪郭設定部３５は、スライスごとに、サイノグラムの左右端の間の距離Ｌ（θ）を角度ごとに算出する（図３（ｄ）参照）。
【００９９】
次に、ステップＳ４２３において、輪郭設定部３５は、スライスごとに、距離Ｌ（θ）の平均ＭＬおよび距離Ｌ（θ）の標準偏差ＳＤＬを算出する。
【０１００】
次に、ステップＳ４２４において、輪郭設定部３５は、平均距離ＭＬが第２距離閾値ｔｈＭＬ２以下であるか、または標準偏差ＳＤＬが第２標準偏差閾値ｔｈＳＤＬ２以上であるか否かを判定する。ＭＬ≦ｔｈＭＬ２またはＳＤＬ≧ｔｈＳＤＬ２を満たすスライスは、サイノグラムの左右端間の距離が小さいか、または左右端の曲線のがたつきが大きいスライスであり、小脳部領域に属するスライスの中でも小脳部下端に近いスライスであると考えられる。
【０１０１】
なお、頭頂部領域判定に用いる第１距離閾値ｔｈＭＬ１と小脳部領域判定に用いる第２距離閾値ｔｈＭＬ２は異なることが好ましい。これは、頭頂部領域に属するスライスよりも小脳部領域に属するスライスのほうが計数値が多く、サイノグラムの左右端が大きいためである（図４、６、７参照）。第１距離閾値ｔｈＭＬ１および第２距離閾値ｔｈＭＬ２としては、あらかじめ設定した値を用いてもよいし、たとえば全スライスの中で最も大きい平均距離ＭＬに対してそれぞれ異なる定数を乗じて定めてもよい。
【０１０２】
ＭＬ≦ｔｈＭＬ２またはＳＤＬ≧ｔｈＳＤＬ２である場合は、小脳部領域に属するスライスと判断し、ステップＳ４２５に進む。一方、ＭＬ＞ｔｈＭＬ１かつＳＤＬ＜ｔｈＳＤＬ１である場合は、ステップＳ４２６に進む。
【０１０３】
次に、ステップＳ４２５において、輪郭設定部３５は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定してステップＳ４１１に戻る。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してステップＳ４２１に戻ってもよい。なお、ここでいう低次の次数としては、頭頂部領域に属するスライスに対して図１０のステップＳ４１７で適用される次数と同一でもよいし、異なっていてもよい。
【０１０４】
他方、ステップＳ４２４でＭＬ＞ｔｈＭＬ１かつＳＤＬ＜ｔｈＳＤＬ１と判定されると、ステップＳ４２６において、輪郭設定部３５は、このスライスについては大脳部領域に属すると判断して小脳部下端からの逆順のサイノグラムの抽出を終了し、図９のステップＳ４３に進む。
【０１０５】
以上の手順により、小脳部領域に属するスライスに対して低次のフーリエ近似を用いて輪郭４３を設定することができる。
【０１０６】
図１２は、図９のステップＳ４３で輪郭設定部３５により実行される、大脳部領域に属するスライスの輪郭設定処理の手順を示すサブルーチンフローチャートである。
【０１０７】
図１０および図１１に示す手順において抽出されなかったスライスは、大脳部領域に属するスライスであると考えられる。そこで、ステップＳ４３１において、輪郭設定部３５は、この大脳領域に属するスライスについて、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して高次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定してステップＳ４１１に戻る。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接高次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してステップＳ４２１に戻ってもよい。なお、ここでいう高次の次数は、頭頂部領域および小脳部領域に属するスライスに対して適用される次数よりも高いものとする。
【０１０８】
以上の手順により、大脳部領域に属するスライスに対して低次のフーリエ近似を用いて輪郭４３を設定することができる。
【０１０９】
続いて、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報について、ユーザによる入力部２３を介した指示にもとづいて取得する場合の例について図１３−１５を用いて説明する。
【０１１０】
図１３は、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を、ユーザに設定させるための画像の一例を示す説明図である。
【０１１１】
スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報について、ユーザによる入力部２３を介した指示にもとづいて取得する場合には、領域境界の設定を受け付けるための画像（以下、領域境界設定用の画像という）を表示部２２に表示させる。
【０１１２】
領域境界設定用の画像として、たとえば図１３に示すような頭部のサジタル断面画像を表示部２２に表示する場合、ユーザは、入力部２３を介してカーソル５１を操作し、各領域の境界を設定する。また、領域境界設定用の画像として、減弱補正前の再構成画像に対して輪郭抽出部３４により抽出された自動抽出左右端４１に応じた自動抽出輪郭を重畳した重畳画像を表示部２２に表示させてもよい。この場合でも、ユーザはカーソル５１を操作して各領域の境界となるスライス間にカーソル５１を置けばよい。輪郭設定部３５は、このカーソル５１の位置の情報にもとづいて、各スライスが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を取得する。
【０１１３】
図１４は、頭部内の位置によらずなめらかな輪郭４３を設定する際の他の手順の一例を示すフローチャートである。図１４において、Ｓに数字を付した符号はフローチャートの各ステップを示す。図１４は、図８に示した手順において、スライスの所属領域情報をユーザによる入力部２３を介した指示にもとづいて取得する場合の手順の一例である。図８と同等のステップには同一符号を付し、重複する説明を省略する。
【０１１４】
ステップＳ３で自動抽出左右端４１を抽出されると、ステップＳ１０１において、再構成部３７は、領域境界設定用の画像を表示部２２に表示させる（図１３参照）。
【０１１５】
次に、ステップＳ１０２において、輪郭設定部３５は、領域境界設定用の画像に対してユーザにより設置されたカーソル５１の位置の情報にもとづいて、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報を取得する。
【０１１６】
次に、ステップＳ１０３において、輪郭設定部３５は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線の粗さに応じた精度で、左右端の曲線をカーブフィッティングすることにより、スライスごとに近似左右端４２を求めて輪郭４３を設定する。
【０１１７】
図１５は、図１４のステップＳ１０３で輪郭設定部３５により実行される輪郭設定処理の手順を示すサブルーチンフローチャートである。
【０１１８】
ステップＳ１０３１において、輪郭設定部３５は、ユーザ指示に応じて設定した頭頂部領域に属するスライスについて、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定する。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してもよい。
【０１１９】
次に、ステップＳ１０３２において、輪郭設定部３５は、ユーザ指示に応じて設定した小脳部領域に属するスライスについて、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定する。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接低次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してもよい。
【０１２０】
次に、ステップＳ１０３３において、輪郭設定部３５は、ユーザ指示に応じて設定した大脳部領域に属するスライスについて、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線に対して、高次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求め（図５（ｃ）、図７参照）、輪郭４３を設定する。なお、輪郭設定部３５は、自動抽出左右端４１の左右端の曲線に対して直接高次のフーリエ近似を行うことにより、スライスごとに近似左右端４２を求めて輪郭４３を設定してもよい。
【０１２１】
なお、各領域を自動設定する場合と同様に、頭頂部領域および小脳部領域に属するスライスに適用される次数は大脳部領域に属するスライスに適用される次数より低ければよく、頭頂部領域および小脳部領域に属するスライスに適用される次数は互いに異なってもよい。また、ユーザ指示に応じて各領域を設定する場合にも、ステップＳ１０３１において図１０のステップＳ４１５−Ｓ４１８の処理を行ってもよい。すなわち、ユーザ指示に応じて各領域を設定する場合にも、頭頂部領域の特に頭頂部上端に位置するスライスに対してのみ低次のフーリエ近似を適用してもよいし、頭頂部上端付近のスライスと頭頂部領域の他の位置のスライスとで異なる次数のフーリエ近似を適用してもよい。
【０１２２】
図１４、図１５に示すように、スライスのそれぞれが頭頂部領域、小脳部領域および大脳部領域のいずれに属するかの情報について、ユーザによる入力部２３を介した指示にもとづいて取得する場合であっても、頭部内の位置によらずなめらかな輪郭４３を設定することができる。
【０１２３】
本実施形態に係る核医学診断装置１は、自動抽出左右端４１にもとづいて生成された３６０度分のサイノグラムの左右端の曲線の粗さに応じた精度で、左右端の曲線をカーブフィッティングすることにより、スライスごとに近似左右端４２を求めることができる。このため、計数値が少なく自動抽出左右端４１の曲線が粗くなってしまうスライスであっても、なめらかな輪郭４３を設定することができる。したがって、単に自動抽出左右端４１にもとづく輪郭を用いる場合に比べ、頭部内の位置によらずなめらかな輪郭４３を設定することができ、計数値が少なく自動抽出左右端４１の曲線が粗いスライスにおける減弱マップの精度を大幅に向上させることができる。
【０１２４】
よって、本実施形態に係る核医学診断装置１によれば、Ｘ線ＣＴ装置などの他のモダリティのデータを用いることなく、容易に適切な形状を有する減弱マップを生成することができる。このため、減弱補正後の核医学画像にもとづいて求められる脳血流（ＣＢＦ）の値の精度を向上させることができるとともに、誤診の可能性を大幅に低下させることができる。
【０１２５】
なお、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【０１２６】
たとえば、輪郭設定部３５による輪郭４３の設定は、自動抽出左右端４１にもとづく輪郭の画像を確認したユーザによる指示に応じて実行されてもよい。この場合、図８のステップＳ３の後ステップＳ４の前（図１４のステップＳ３の後ステップＳ１０１の前）に、再構成部３７は、減弱補正前の再構成画像に対して、輪郭抽出部３４により抽出された自動抽出左右端４１に応じた自動抽出輪郭、または輪郭設定部３５が全てのスライスに対して同一の次数でフーリエ近似することにより設定した輪郭４３を重畳した重畳画像を表示部２２に表示させる。ユーザは、この重畳画像を確認して、低次のフーリエ近似が必要なスライスが有るか否かを判定し、必要な場合にのみステップＳ４以降（図１４のステップＳ１０１以降）を実行するよう入力部２３を介して指示する。そして、輪郭設定部３５は、この指示を受けた場合のみ、ステップＳ４以降（図１４のステップＳ１０１以降）を実行すればよい。
【０１２７】
また、本発明の実施形態では、フローチャートの各ステップは、記載された順序に沿って時系列的に行われる処理の例を示したが、必ずしも時系列的に処理されなくとも、並列的あるいは個別実行される処理をも含むものである。
【符号の説明】
【０１２８】
１  核医学診断装置
２２  表示部
２３  入力部
３４  輪郭抽出部
３５  輪郭設定部
３６  減弱マップ生成部
３７  再構成部
４３  輪郭
【図１】
2015152494.tif 000003
【図２】
2015152494.tif 000004
【図３】
2015152494.tif 000005
【図８】
2015152494.tif 000006
【図９】
2015152494.tif 000007
【図１０】
2015152494.tif 000008
【図１１】
2015152494.tif 000009
【図１２】
2015152494.tif 000010
【図１４】
2015152494.tif 000011
【図１５】
2015152494.tif 000012
【図４】
2015152494.tif 000013
【図５】
2015152494.tif 000014
【図６】
2015152494.tif 000015
【図７】
2015152494.tif 000016
【図１３】
2015152494.tif 000017