JP7154012B2

JP7154012B2 - 画像処理プログラム、画像処理装置、および画像処理方法

Info

Publication number: JP7154012B2
Application number: JP2018018024A
Authority: JP
Inventors: 一男浜田; 和宏西川
Original assignee: Nihon Medi Physics Co Ltd
Current assignee: Nihon Medi Physics Co Ltd
Priority date: 2018-02-05
Filing date: 2018-02-05
Publication date: 2022-10-17
Anticipated expiration: 2037-09-14
Also published as: JP2019053022A

Description

本発明は、画像処理プログラム、画像処理装置、および画像処理方法に関する。

癌の骨転移を検出する手法として、骨転移診断用の放射性医薬品投与により得られた二次元の骨シンチグラムを用いる方法が、広く用いられている。骨シンチグラムを用いた診断では、平面上に投影された画像上で異常集積部位を確認することにより、癌の骨転移部位の検出が行われている。

また最近では、骨シンチグラムから骨領域を抽出した上で、骨領域上における癌の転移部位を抽出し、癌の骨転移領域の面積と骨領域全体の面積の比に基づいて算出されたBone Scan Index（以下、ＢＳＩ）という指標が提案され、臨床現場で用いられている。ＢＳＩは、癌の骨への浸潤の広がりを表す指標であり、治療効果判定や経過観察に役立つとされている（非特許文献１）。

Yusuf E. Erdi他「Quantitative Bone Metastases Analysis Based on Image Segmentation」The Journal of NUCLEAR MEDICINE 1997, 38:1401-1406.

上述したとおり、従来より用いられているＢＳＩは、骨シンチグラムにおける癌の転移部位の面積、すなわち画素数に基づいて計算された指標である。このことは、非特許文献１の１４０２頁右欄１～２８行には、骨への転移部分の関与を領域のピクセル数で計算することが記載され、同欄１８～１９行に「式３は、各骨の領域の大きさは、骨の重量に比例しているという重要な前提を述べている」と記載されていることからも分かる。つまり、従来は、２次元の領域の大きさに基づいてＢＳＩを求めていた。

そのため、例えば、シンチグラムの奥行き方向に複数の腫瘍が存在している場合等には、それらの腫瘍の占める面積が計算に用いられず、ＢＳＩが過小評価された数値となってしまう場合があった。

本発明は、上記背景に鑑み、癌の骨領域への浸潤の程度を表す精度の高い指標を提供することを目的とする。

本発明のプログラムは、三次元核医学画像に基づいて骨転移に関する指標を求めるためのプログラムであって、コンピュータに、被験者の三次元形態画像および三次元核医学画像のデータを入力するステップと、前記三次元形態画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域に対応する前記三次元核医学画像における前記骨領域を特定するステップと、前記三次元核医学画像における前記骨領域の画素値を正規化するステップと、前記正規化後の前記三次元核医学画像における前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求めるステップとを実行させる。

この構成により、三次元核医学画像における骨領域の体積と骨転移領域の体積とに基づいて、骨転移に関する指標を求めるので、奥行き方向における情報を考慮した精度の高い指標を求めることができる。

本発明のプログラムにおいて、前記骨領域に対応する前記三次元核医学画像における前記骨領域を特定するステップは、前記三次元形態画像における前記骨領域のデータに基づいて前記骨領域のマスク画像を生成し、前記骨領域のマスク画像と三次元核医学画像とが重なる部分として、前記三次元核医学画像における前記骨領域を特定しても良い。例えば、前記骨領域のマスク画像を前記三次元核医学画像に乗じて前記三次元核医学画像における前記骨領域を特定してもよい。

この構成により、三次元形態画像に基づいて求めた骨領域のマスク画像を用いて、三次元核医学画像における骨領域（の画素）を適切に特定することができる。

本発明の別の態様のプログラムは、三次元核医学画像に基づいて骨転移に関する指標を求めるためのプログラムであって、コンピュータに、被験者の三次元核医学画像のデータを入力するステップと、前記三次元核医学画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域の画素値を正規化するステップと、前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求めるステップとを実行させる。

本発明のプログラムにおいて、前記骨領域を抽出するステップは、前記三次元核医学画像における放射性医薬品の生理的集積部位を特定するステップと、前記三次元核医学画像において所定の閾値以上の画素値を有する画素を抽出するステップと、前記抽出された画素で構成される領域から、前記放射性医薬品の生理的集積部位を除外した領域を、前記骨領域として抽出するステップとを実行させてもよい。

所定の閾値以上の画素値を有する画像から、膀胱や腎臓等の放射性医薬品の生理的集積部位を除去することにより、適切に、三次元核医学画像から骨領域を抽出することができる。

本発明のプログラムにおいて、前記骨転移に関する指標を求めるステップは、被験者の全身の骨領域の体積に基づいて前記指標を求めてもよい。また、被験者の一部の骨領域の体積に基づいて前記指標を求めてもよい。

本発明の別の態様のプログラムは、三次元核医学画像の画像処理を行うためのプログラムであって、コンピュータに、被験者の三次元形態画像および三次元核医学画像のデータを入力するステップと、前記三次元形態画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域に対応する前記三次元核医学画像における前記骨領域を特定するステップと、前記三次元核医学画像の前記骨領域の画素値を正規化するステップと、前記三次元核医学画像の前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、方向の異なる２つ以上の画像を並べて表示し、表示された画像上に前記骨転移領域を表示するステップとを実行させる。

この構成により、見る方向の異なる２つ以上の三次元核医学画像が表示されるので、骨転移領域を立体的に把握することができ、診断に役立てることができる。なお、骨転移領域を表示する三次元核医学画像の方向は、任意の方向とすることができる。方向の指定は、操作者が表示したい方向を画面上で指定して行うことができる。

本発明の別の態様のプログラムは、三次元核医学画像の画像処理を行うためのプログラムであって、コンピュータに、被験者の三次元核医学画像のデータを入力するステップと、前記三次元核医学画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域の画素値を正規化するステップと、前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、方向の異なる２つ以上の画像を並べて表示し、表示された画像上に前記骨転移領域を表示するステップとを実行させる。

本発明のプログラムにおいて、前記正規化するステップは、被験者の体重、骨重量、または骨ミネラル量に関するデータを用いて正規化を行ってもよい。

本発明の画像処理装置は、被験者の三次元形態画像および三次元核医学画像のデータを入力する画像データ入力部と、前記三次元形態画像に基づいて少なくとも一部の骨領域を抽出する骨領域抽出部と、前記骨領域に対応する前記三次元核医学画像における前記骨領域を特定する骨領域特定部と、前記三次元核医学画像における前記骨領域の画素値を正規化する正規化部と、前記三次元核医学画像における前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求める骨転移領域抽出部と、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求める指標値計算部とを備える。

本発明の画像処理装置は、被験者の三次元核医学画像のデータを入力する画像データ入力部と、前記三次元核医学画像に基づいて少なくとも一部の骨領域を抽出する骨領域抽出部と、前記骨領域の画素値を正規化する正規化部と、前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求める骨転移領域抽出部と、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求める指標値計算部とを備える。

本発明の画像処理方法は、三次元核医学画像に基づいて骨転移に関する指標を求めるための方法であって、被験者の三次元形態画像および三次元核医学画像のデータを入力するステップと、前記三次元形態画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域に対応する前記三次元核医学画像における前記骨領域を特定するステップと、前記三次元核医学画像における前記骨領域の画素値を正規化するステップと、正規化後の前記三次元核医学画像における前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求めるステップとを備える。

本発明の画像処理方法は、三次元核医学画像に基づいて骨転移に関する指標を求めるための方法であって、被験者の三次元核医学画像のデータを入力するステップと、前記三次元核医学画像に基づいて少なくとも一部の骨領域を抽出するステップと、前記骨領域の画素値を正規化するステップと、前記骨領域において所定の閾値以上の画素値を有する骨転移領域を求めるステップと、前記骨領域の体積と前記骨転移領域の体積とに基づいて、骨転移に関する指標を求めるステップとを備える。

本発明によれば、奥行き方向における情報を考慮して、骨転移に関する精度の高い指標を求めることができる。

第１の実施の形態の画像処理装置の構成を示す図である。第１の実施の形態の画像処理装置の動作を示す図である。モニタに表示された出力結果の例を示す図である。骨転移領域抽出部にて求めた骨転移領域のＶＯＩを出力した例を示す図である。第２の実施の形態の画像処理装置の構成を示す図である。第２の実施の形態の画像処理装置の動作を示す図である。

以下、本発明の実施の形態の画像処理プログラム、画像処理装置、及び画像処理方法について説明する。
（第１の実施の形態）
本実施の形態では、三次元核医学画像としてＳＰＥＣＴ画像を用い、三次元形態画像としてＣＴ画像を用いる例を取り上げて説明する。現在は、ＳＰＥＣＴとＣＴとを一体化したＳＰＥＣＴ／ＣＴ装置が開発されており、ＳＰＥＣＴ検査の施行に合わせてＣＴ画像を撮影することが可能である。本実施の形態では、ＳＰＥＣＴ／ＣＴ装置を用いて撮影したＳＰＥＣＴ画像とＣＴ画像に基づいて処理を行うものとする。なお、本実施の形態はあくまでも一例であり、本発明の内容を限定する意図では無い。

図１は、画像処理装置１の構成を示す図である。画像処理装置１は、画像データ入力部１０と、演算処理部１１と、出力部１２とを備えている。この画像処理装置１は、画像処理プログラムを実行することによって演算処理部１１に各種の機能を行わせることにより実現される。演算処理部１１は、図１に示すとおり、骨領域抽出部２０、骨領域特定部２１、正規化部２２、骨転移領域抽出部２３、指標値計算部２４の機能を有する。

画像処理装置１のハードウェアの一例は、ＣＰＵ、ＲＡＭ、ＲＯＭ、ハードディスク、モニタ、キーボード、通信インターフェース、周辺機器インターフェース等を備えたコンピュータである。画像処理プログラムは、例えば、ＲＯＭに記憶されている。ＣＰＵは、ＲＯＭから画像処理プログラムを読み出して画像処理プログラムを実行する。

画像処理装置１は、画像データ入力部１０を有している。画像データ入力部１０は、被験者のＳＰＥＣＴ画像およびＣＴ画像のデータを入力する機能を有する。ＳＰＥＣＴ画像は、典型的には骨代謝が亢進した造骨部位に多く吸着する放射性医薬品（例えば、^99mTc-HMDP（日本メジフィジックス株式会社製、クリアボーン(R)注））を用いて撮影した画像である。したがって、放射能量が高い領域、すなわちＳＰＥＣＴ画像の画素値が高い領域は、癌の骨転移が疑われる部分である。

画像データ入力部１０は、例えば、ネットワークインターフェースによって実現され、ＳＰＥＣＴ／ＣＴ装置から、被験者のＳＰＥＣＴ画像およびＣＴ画像のデータを受信する。また、画像データ入力部１０は、ＵＳＢのインターフェースによって実現されてもよい。ＳＰＥＣＴ画像およびＣＴ画像のデータを記憶したＵＳＢメモリをＵＳＢインターフェースに挿入することにより、ＵＳＢメモリからＳＰＥＣＴ画像およびＣＴ画像のデータの入力を受け付ける。なお、本実施の形態で入力されるＳＰＥＣＴ画像とＣＴ画像は、ＳＰＥＣＴ／ＣＴ装置によって撮影されたデータであり、両者の画像の位置合わせ（コレジストレーション）がなされている。

次に、演算処理部１１の各機能について説明する。骨領域抽出部２０は、ＣＴ画像に基づいて、三次元的に被験者の骨領域を抽出する。例えば骨領域抽出部２０は、ＣＴ画像上で閾値以上の画素値を有する領域を骨領域として抽出する。骨領域の抽出に用いる閾値は、骨領域を示すハンスフィールド値として当業者に一般に知られている値を用いることができ、本実施の形態では、１５２を用いている。骨領域抽出部２０は、骨領域を求めたら、その骨領域の体積を算出する。

なお、閾値による抽出では、骨の内部が抜けてしまったりノイズが生じてしまったりする場合がある。そこで、骨領域抽出部２０は、閾値以上の画素値を有する画素により構築した画像データに、公知の方法によるクロージング処理を加えて残った画素を、骨領域として定めてもよい。また、ノイズ除去の方法としてボクセル数が一定閾値以下のボリュームを除去する方法やノイズをマウスやスタイラスペン等を用いて手動で指定して除去する方法を用いてもよい。

骨領域特定部２１は、骨領域抽出部２０にて求めた三次元の骨領域に対応する領域のみＳＰＥＣＴ画像の画素値を維持し、それ以外の画素値を０とする処理を行い、ＳＰＥＣＴ画像における骨領域の画素を抽出する。具体的な処理の例としては、骨領域抽出部２０にて求めた三次元の骨領域の画素値を１、それ以外の画素値を０とするマスク画像を生成し、ＳＰＥＣＴ画像に乗じる。ＳＰＥＣＴ画像とＣＴ画像は位置合わせがなされているので、マスク画像とＳＰＥＣＴ画像の対応する画素の画素値を乗じることで、ＳＰＥＣＴ画像において、骨領域の画素値のみを得ることができる。

正規化部２２は、ＳＰＥＣＴ画像の画素値を正規化する機能を有する。正規化は、撮影条件等が異なる検査結果を比較するために行う前処理である。正規化部２２は、投与した放射能（放射性医薬品の量）が全身に均一に分布し、かつ排泄されていないという状況と比較して、各画素の放射能濃度が何倍になるかを示す量に変換する。一例において、正規化部２２は、投与した放射性医薬品の放射能量と被験者の体重を用いて正規化を行う。具体的には、下記の式（１）により、ＳＰＥＣＴ画像の各画素値を正規化する。

（ＳＰＥＣＴ画素値）／（被験者に投与した放射能量／被験者の体重）・・・（１）

骨転移領域抽出部２３は、癌の転移が認められる領域を抽出する機能を有する。骨転移領域抽出部２３は、具体的には、ＳＰＥＣＴ画像の骨領域から、所定の閾値以上の画素値を有する画素を骨転移領域として抽出する。骨転移領域抽出部２３は、抽出した骨転移領域の体積を求める。

指標値計算部２４は、骨領域の体積に対する骨転移領域の体積を、骨への浸潤の程度を示す指標として求める。すなわち、下記の式（２）により、指標値を求める。

（骨転移領域の体積の合計）／（骨体積）×１００・・・（２）

なお、本実施の形態では、式（２）の分母である骨領域の体積（式（２）では「骨体積」と表示）としては、全身の骨領域の体積を用いるが、必ずしも全身の骨領域の体積を用いなくてもよい。

出力部１２は、求めた指標値を例えばモニタに出力する。このときに、出力部１２は、指標値と共に、被験者のＳＰＥＣＴ画像、ＣＴ画像およびその合成画像をモニタに出力してもよい。

図２は、第１の実施の形態の画像処理装置１の動作の例を示す図である。
画像処理装置１は、まず、ＳＰＥＣＴ／ＣＴ装置から被験者のＳＰＥＣＴ画像およびＣＴ画像のデータを読み込む（Ｓ１０）。画像処理装置１の骨領域抽出部２０は、ＣＴ画像に基づいて被験者の三次元の骨領域を抽出すると共に、抽出した骨領域の体積を求める（Ｓ１１）。上述したとおり、本実施の形態においては、骨領域の抽出には、ハンスフィールド値の閾値として１５２を用い、１５２以上のハンスフィールド値を有する画素を骨領域として抽出している。

次に、画像処理装置１の骨領域特定部２１は、ＣＴ画像に基づいて求めた骨領域に対応するＳＰＥＣＴ画像の骨領域を特定する（Ｓ１２）。具体的には、ＣＴ画像に基づいて求めた骨領域のマスク画像を生成し、マスク画像をＳＰＥＣＴ画像に乗じることにより、ＳＰＥＣＴ画像から骨領域に対応する部分の領域を特定する。膀胱や腎臓等の生理的集積部位では、異常がなくてもＳＰＥＣＴ画像の画素値が高くなる場合があるが、ＣＴ画像に基づいて求めた骨領域のマスク画像を乗じる方法により、生理的集積部位を除去して適切な診断を行うことが可能となる。

画像処理装置１の正規化部２２は、ＳＰＥＣＴ画像の骨領域における画素値を正規化する（Ｓ１３）。本実施の形態において、具体的には、正規化部２２は、被験者の体重のデータを用いて、上述した式（１）によって正規化を行う。骨転移領域抽出部２３は、正規化されたＳＰＥＣＴ画像において、所定の閾値以上の画素値を有する骨転移領域を抽出する（Ｓ１４）。ここで、抽出される骨転移領域は、ＳＰＥＣＴ画像の骨領域から抽出されるので、三次元のデータである。また、骨転移領域抽出部２３は、抽出された骨転移領域の体積を求める。

そして、画像処理装置１の指標値計算部２４は、骨領域の体積と骨転移領域の体積とに基づいて、上述した式（２）により、骨転移に関する指標値を計算する（Ｓ１５）。画像処理装置１の出力部１２は、求めた指標値と被験者のＳＰＥＣＴ画像、ＣＴ画像をモニタに出力する（Ｓ１６）。

図３は、モニタに表示された出力結果の例を示す図である。図３において、右端の「ＣＴ」と記載された画像はＣＴ画像、右から２番目の「ＮＭ」と記載された画像はＳＰＥＣＴ画像である。なお、「ＮＭ」は、核医学（Nuclear Medicine）の略である。右から３番目の「Ｆｕｓｅｄ」と記載された画像は、ＣＴ画像とＳＰＥＣＴ画像を重畳した画像であり、左端の「ＶＯＩ」と記載された画像は、骨転移領域として抽出されたＶＯＩをハイライト表示した画像である。図３における「ＶＯＩ」と記載された画像では、５箇所のＶＯＩが設定されている。

画像の下に表示されているデータにおいて、「ＭＢＶ」は体積に該当する数値、「ＰｉｘｅｌＮｕｍ」は画素数、「ＳＵＶｍａｘ」は、対象部位内でのＳＵＶの最大値を示している。ここで、「ＳＵＶ（standardized uptake value）」とは、放射性薬剤の腫瘍や臓器への集積の強さを表すための指標であり、画像で測定した腫瘍や臓器の放射能濃度÷（放射能投与量÷体重）で求められる。また、図３に示す例では、骨転移領域のＶＯＩの体積（Total=7.072）を、全身の骨部分として抽出された領域の体積（CT Masked Bone=2087.664）で割って、指標値として0.33875が求められている。

また、出力部１２は、骨転移領域抽出部２３にて骨転移領域を抽出した画像を出力してもよい。図４は、骨転移領域抽出部２３にて求めた骨転移領域のＶＯＩを出力した例を示す図である。図４には４つの画像が写っているが、左から順に、被験者の骨盤付近の長軸断画像、矢状面画像、冠状面画像、ＭＩＰ（最大値投影）画像である。図４に矢印で示すように、骨転移領域を抽出した出力画像では、所定の閾値以上の画素値を有する部位にＶＯＩが設定されている。この出力画像では、様々な方向からの画像においてＶＯＩが表示されており、三次元の骨転移領域の形状を立体的に把握することができるので、診断に役立てることができる。

なお、本例では、ＳＰＥＣＴ画像の長軸断画像、矢状面画像、冠状面画像、ＭＩＰ（最大値投影）画像を並べて表示する例を挙げたが、長軸断画像、矢状面画像、冠状面画像のうちの少なくとも２つを並べて表示することとしてもよい。この構成によっても、骨転移領域を立体的に把握することができる。また、表示するＳＰＥＣＴ画像は、上記の４つの画像に限らず、任意の方向のＳＰＥＣＴ画像を表示することとしてもよい。この場合、操作者が見たい方向を指定できるようなユーザ・インターフェースを備えてもよい。このようなユーザ・インターフェースとしては、操作者が、ＳＰＥＣＴ画像が表示されたタッチパネルをスワイプすることにより、スワイプした方向に画像が回転するようなインターフェースが考えられる。
以上、第１の実施の形態の画像処理装置１の構成および動作について説明した。

第１の実施の形態の画像処理装置１は、ＣＴ画像で求めた三次元の骨領域の体積と、ＳＰＥＣＴ画像で求めた三次元の骨転移領域の体積とに基づいて、骨転移に関する指標値を求めるので、奥行き方向における情報を考慮した精度の高い指標を求めることができる。

（第２の実施の形態）
図５は、第２の実施の形態の画像処理装置２の構成を示す図である。第２の実施の形態の画像処理装置２は、核医学画像を用いて骨転移領域に関する指標を求める装置である。本実施例では、核医学画像として、ＳＰＥＣＴ画像を用いた場合を例として説明する。第２の実施の形態の画像処理装置２の基本的な構成は、第１の実施の形態の画像処理装置１と同じである。同じ機能を有する構成要素には、同じ番号を付与している。

第２の実施の形態の画像処理装置２は、演算処理部１１が、ＳＰＥＣＴ画像に基づいて被験者の骨領域を抽出する骨領域抽出部２５を有している。骨領域抽出部２５は、ＳＰＥＣＴ画像において所定の閾値以上の画素値を有する領域を骨領域として抽出する。ただし、この処理だけだと、膀胱や腎臓等の生理的集積部位を骨領域として誤って抽出してしまうおそれがあるので、骨領域抽出部２５は、生理的集積部位を骨領域から除去する処理を行う。

本実施の形態においては、画像処理装置２はＳＰＥＣＴ画像を表示し、そのＳＰＥＣＴ画像上において、ユーザにマウスやスタイラスペン等を用いて生理的集積部位を指定させることによって、生理的集積部位を特定する。画像処理装置２は、部位指定部１３を有している。部位指定部１３は、ＳＰＥＣＴ画像上において、ユーザが指定した点の座標を認識し、骨領域抽出部２５にその座標データを渡す。骨領域抽出部２５は、ユーザが指定した点を含むひとまとまりの領域を生理的集積部位として特定する。ひとまとまりの領域とは、値の近い画素がまとまった領域である。この構成により、ＳＰＥＣＴ画像から骨領域を抽出することができる。

本実施の形態では、ユーザが生理的集積部位を指定する例を挙げたが、生理的集積部位のモデルとなる画像データを記憶しておき、パターンマッチングにより、生理的集積部位を特定する等の方法をとることもできる。また、多数の被験者の生理的集積部位のデータがある場合には、それらを教師データとして学習して求めたモデルを用いて、新しい被験者のＳＰＥＣＴ画像から生理的集積部位を求めることができる。なお、被験者のＣＴ画像がある場合には、ＣＴ画像の情報を用いて、生理的集積部位を除外することとしてもよい。演算処理部１１のその他の機能は、第１の実施の形態の画像処理装置１と同じである。

図６は、第２の実施の形態の画像処理装置２の動作の例を示す図である。
画像処理装置２は、まず、ＳＰＥＣＴ装置から被験者のＳＰＥＣＴ画像のデータを読み込む（Ｓ２０）。画像処理装置２の骨領域抽出部２５は、ＳＰＥＣＴ画像に基づいて被験者の三次元の骨領域を抽出すると共に、抽出した骨領域の体積を求める（Ｓ２１）。

次に、画像処理装置２の正規化部２２は、ＳＰＥＣＴ画像の骨領域における画素値を正規化する（Ｓ２２）。具体的には、正規化部２２は、被験者の体重のデータを用いて、上述した式（１）によって正規化を行う。骨転移領域抽出部２３は、正規化されたＳＰＥＣＴ画像において、所定の閾値以上の画素値を有する骨転移領域を抽出する（Ｓ２３）。ここで、抽出される骨転移領域は、ＳＰＥＣＴ画像の骨領域から抽出されるので、三次元のデータである。また、骨転移領域抽出部２３は、抽出された骨転移領域の体積を求める（Ｓ２３）。

そして、画像処理装置２の指標値計算部２４は、骨領域の体積と骨転移領域の体積とに基づいて、上述した式（２）により、骨転移に関する指標値を計算する（Ｓ２４）。画像処理装置２の出力部１２は、求めた指標値と被験者のＳＰＥＣＴ画像をモニタに出力する（Ｓ２５）。以上、第２の実施の形態の画像処理装置２の構成および動作について説明した。

第２の実施の形態の画像処理装置２は、ＳＰＥＣＴ画像で求めた三次元の骨領域の体積と三次元の骨転移領域の体積とに基づいて、骨転移に関する指標値を求めるので、奥行き方向における情報を考慮した精度の高い指標を求めることができる。また、第２の実施の形態の画像処理装置２は、ＳＰＥＣＴ画像を用いて処理を行えるので、ＣＴ画像のデータのない被験者についても、骨転移に関する指標値を求めることが可能である。

（変形例）
以上、本発明の画像処理装置、画像処理プログラムおよび画像処理方法について、実施の形態を挙げて詳細に説明したが、本発明は上記した実施の形態に限定されるものではない。

上記した実施の形態では、ＳＰＥＣＴ／ＣＴ装置を用いて、ＳＰＥＣＴ画像とＣＴ画像を取得する例を挙げたが、ＳＰＥＣＴ画像とＣＴ画像は、別々に撮影した画像を用いることも可能である。この場合、ＳＰＥＣＴ画像の骨領域における画素値を取得するステップ（Ｓ１２）において、ＳＰＥＣＴ画像とＣＴ画像の位置合わせを行う必要がある。

上記した実施の形態では、核医学画像としてＳＰＥＣＴ画像を用いる例を挙げて説明したが、例えばＮａＦ^１８等投与によるＰＥＴ画像を用いた指標値の計算にも適用することができる。この場合、ＰＥＴ画像とＣＴ画像は、ＰＥＴ／ＣＴ装置を用いて同時に取得してもよいし、別々に撮影した画像を用いてもよい。

上記した実施の形態では、被験者の骨領域を求めるための形態画像としてＣＴ画像を用いる例を挙げて説明したが、形態画像はＣＴ画像に限られず、例えば、ＭＲＩ画像を用いてもよい。

上記した実施の形態では、被験者の体重を用いて正規化を行う例を挙げたが、ＳＰＥＣＴ画像の正規化は、体重以外のデータを用いて行ってもよい。例えば、正規化部２２は、被験者の骨重量を用いて、ＳＰＥＣＴ画像の正規化を行うことができる。これは、投与した放射能（放射性医薬品の量）が骨組織に均一に分布し、かつ排泄されていないという状況と比較して、各画素の放射能濃度が何倍になるかを示す量に変換するものであり、具体的には、次の式（３）によって計算する。

（ＳＰＥＣＴ画素値）／（被験者に投与した放射能量／被験者の骨重量）・・・（３）

なお、骨重量の情報は、例えば、ＣＴ画像に基づいて求めた骨領域のデータに基づいて推定することができる。ＣＴ画像がない場合には、被験者の性別、年齢、身長、体重等から推定することができる。

また、正規化部２２は、骨ミネラル量（ＢＭＣ：Bone Mineral Contents）または骨ミネラル密度（ＢＭＤ：Bone Mineral Density）を用いて、ＳＰＥＣＴ画像の正規化を行ってもよい。具体的には、次の式（４）（５）によって計算することができる。

（ＳＰＥＣＴ画素値）／（被験者に投与した放射能量／骨ミネラル量）・・・（４）

（ＳＰＥＣＴ画素値）／（被験者に投与した放射能量／骨ミネラル密度）・・・（５）

なお、骨ミネラル量、骨ミネラル密度は、被験者の性別、年齢、身長、体重等から推定することができる。日本人についてはデータベースを用いた先行研究（例えば、European Journal of Clinical Nutrition (2001) 55, p.462470）がある。この文献によると、骨ミネラル量や骨ミネラル密度を、次の式から推定することができることが示されている。

骨ミネラル量（男性；kg）= -1.81 - 0.0015×年齢＋1.89×身長＋0.017×体重
骨ミネラル量（女性；kg）= -1.05 - 0.009×年齢＋1.57×身長＋0.017×体重
骨ミネラル密度（男性；g/m²）= 0.934 - 0.00081×年齢＋0.003×体重
骨ミネラル密度（女性；g/m²）= 0.824 - 0.00368×年齢＋0.137×身長＋0.0026×体重

これらの式において、年齢の単位は年、体重の単位はｋｇ、身長の単位はｍである。

本発明によって求めた骨転移に関する指標は、被験者の経過観察や重症度評価に有効に用いる事ができる。
ここで、骨転移に関する指標を計算する範囲は、全身を対象としても良いし、一部としても良い。身体の一部を対象として指標を求める場合は、経過観察等で比較に用いる各データの計算範囲を揃えておく必要がある。例えば、右上腕骨に転移の見られた患者について、当該転移領域の経過観察を詳細に行う場合は、各時間点における画像で右上腕骨のデータのみを用いて骨転移に関する指標を算出し、相互比較を行えば良い。

１，２画像処理装置
１０画像データ入力部
１１演算処理部
１２出力部
１３部位指定部
２０骨領域抽出部
２１骨領域特定部
２２正規化部
２３骨転移領域抽出部
２４指標値計算部
２５骨領域抽出部

Claims

コンピュータに、
被験者のＣＴ画像および骨代謝が亢進した造骨部位に吸着する放射性医薬品を用いて撮影したＳＰＥＣＴ画像のデータを入力するステップと、
前記ＣＴ画像に基づいて少なくとも一部の三次元の骨領域を抽出するステップと、
前記三次元の骨領域に対応する前記ＳＰＥＣＴ画像における前記三次元の骨領域を特定するステップと、
前記ＳＰＥＣＴ画像における前記三次元の骨領域の画素値を正規化するステップと、
前記正規化後の前記ＳＰＥＣＴ画像における前記三次元の骨領域において所定の閾値以上の画素値を有する三次元の骨転移領域を求めるステップと、
前記三次元の骨領域の体積と前記三次元の骨転移領域の体積とに基づいて、奥行き方向における情報を反映した骨転移に関する指標を求めるステップと、
を実行させる画像処理プログラム。
前記三次元の骨領域に対応する前記ＳＰＥＣＴ画像における前記三次元の骨領域を特定するステップは、前記ＣＴ画像における三次元の骨領域のデータに基づいて三次元の骨領域のマスク画像を生成し、
前記三次元の骨領域のマスク画像と前記ＳＰＥＣＴ画像とが重なる部分として、前記ＳＰＥＣＴ画像における前記三次元の骨領域を特定する請求項１に記載の画像処理プログラム。
前記骨転移に関する指標を求めるステップは、被験者の全身の三次元の骨領域の体積に基づいて前記指標を計算する請求項１または２に記載の画像処理プログラム。
前記骨転移に関する指標を求めるステップは、被験者の一部の三次元の骨領域の体積に基づいて前記指標を計算する請求項１または２に記載の画像処理プログラム。
２つ以上の異なる断面画像を並べて表示し、表示された断面画像上に前記三次元の骨転移領域に対応する領域を表示するステップをさらに含む、請求項１乃至４のいずれか１項に記載の画像処理プログラム。
前記正規化するステップは、被験者の体重、骨重量、または骨ミネラル量に関するデータを用いて正規化を行う請求項１乃至５のいずれか１項に記載の画像処理プログラム。
被験者のＣＴ画像および骨代謝が亢進した造骨部位に吸着する放射性医薬品を用いて撮影したＳＰＥＣＴ画像のデータを入力する画像データ入力部と、
前記ＣＴ画像に基づいて少なくとも一部の三次元の骨領域を抽出する三次元の骨領域抽出部と、
前記三次元の骨領域に対応する前記ＳＰＥＣＴ画像における前記三次元の骨領域を特定する骨領域特定部と、
前記ＳＰＥＣＴ画像における前記三次元の骨領域の画素値を正規化する正規化部と、
前記ＳＰＥＣＴ画像における前記三次元の骨領域において所定の閾値以上の画素値を有する三次元の骨転移領域を求める骨転移領域抽出部と、
前記三次元の骨領域の体積と前記三次元の骨転移領域の体積とに基づいて、奥行き方向における情報を反映した骨転移に関する指標を求める指標値計算部と、
を備える画像処理装置。
被験者のＣＴ画像および骨代謝が亢進した造骨部位に吸着する放射性医薬品を用いて撮影したＳＰＥＣＴ画像のデータを入力するステップと、
前記ＣＴ画像に基づいて少なくとも一部の三次元の骨領域を抽出するステップと、
前記三次元の骨領域に対応する前記ＳＰＥＣＴ画像における前記三次元の骨領域を特定するステップと、
前記ＳＰＥＣＴ画像における前記三次元の骨領域の画素値を正規化するステップと、
正規化後の前記ＳＰＥＣＴ画像における前記三次元の骨領域において所定の閾値以上の画素値を有する三次元の骨転移領域を求めるステップと、
前記三次元の骨領域の体積と前記三次元の骨転移領域の体積とに基づいて、奥行き方向における情報を反映した骨転移に関する指標を求めるステップと、
を備える画像処理方法。