(57)【要約】【課題】骨シンチグラフィー画像の自動正規化の安定性を向上させる。【解決手段】好適な実施形態は、骨シンチグラフィー画像を表す画像データの画素値ヒストグラムを作成することと；前記画素値ヒストグラムに基づいて、画素値に関する複数の閾値を設定することと；前記設定した複数の閾値の各々に対して画素値平均値を計算することと；前記計算した画素値平均値を、値の大きな順序に並べることと；前記順序に並べられた前記画素値平均値の集合の少なくとも一部に基づいて、前記画像データを正規化するための基準値を決定することと；を含み、前記基準値を決定することが、前記順序に並べられた前記画素値平均値の集合のうち、前記画素値平均値が小さな領域を近似する直線を一つ決定することと；前記直線に基づいて前記基準値を計算することと；を含む。
2017010119.tif 000002
【特許請求の範囲】
【請求項１】
（ａ）核医学の手法により得られる骨シンチグラフィー画像を正規化するための方法であって、
（ｂ）前記骨シンチグラフィー画像を表す画像データを読み込むことと；
（ｃ）前記画像データの画素値ヒストグラムを作成することと；
（ｄ）前記画素値ヒストグラムに基づいて、画素値に関する複数の閾値を設定することと；
（ｅ）前記設定した複数の閾値の各々に対して画素値平均値を計算することと；
（ｆ）前記計算した画素値平均値を、値の大きな順序に並べることと；
（ｇ）前記順序に並べられた前記画素値平均値の集合の少なくとも一部に基づいて、前記画像データを正規化するための基準値を決定することと；
を含み、前記基準値を決定することが、
（ｈ）前記順序に並べられた前記画素値平均値の集合のうち、前記画素値平均値が小さな領域を近似する直線を一つ決定することと；
（ｉ）前記直線に基づいて前記基準値を計算することと；
を含むことを特徴とする、方法。
【請求項２】
前記直線を一つ決定することが、前記画素値平均値の集合から、値の大きな画素値平均値をいくつか除去すると共に、残った部分を近似する近似直線を計算し、さらに、該近似直線と、該近似直線の計算に用いた画素値平均値の集合との相関係数を計算することを含む、請求項１に記載の方法であって、該方法は更に、
（ｊ）前記除去する画素値平均値の個数を変えて、それぞれの場合に前記近似直線及び前記相関係数を計算し、該計算した複数の相関係数の値に基づいて、前記直線を一つ決定することを含む、方法。
【請求項３】
前記画素値平均値が小さな領域は、前記順序に並べられた前記画素値平均値の集合の値の変化率に基づいて決定される、請求項１に記載の方法。
【請求項４】
前記基準値は、前記直線の切片に基づいて定められる、請求項１に記載の方法。
【請求項５】
前記基準値に基づいて、前記画像データを正規化すると共に、該正規化した画像データを表示することを更に含む、請求項１から４のいずれかに記載の方法。
【請求項６】
前記画素値ヒストグラムに基づいて、画素値に関する複数の閾値を設定することを、多重閾値法を使って行う、請求項１から５のいずれかに記載の方法。
【請求項７】
処理手段及び記憶手段を備えるシステムであって、前記記憶手段がプログラム命令を格納し、前記プログラム命令は、前記処理手段に実行されることにより、前記システムに、上記請求項のいずれかに記載の方法を遂行させるように構成される、システム。
【請求項８】
システムの処理手段に実行されることにより、前記システムに、上記請求項のいずれかに記載の方法を遂行させるように構成されるプログラム命令を備える、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、核医学の手法により得られるシンチグラフィー画像を正規化することに関する。
【背景技術】
【０００２】
核医学検査で用いられるシンチグラフィーは、核医学画像化技術の一種であり、放射性薬剤を体内に投与し、その核種の崩壊に起因して放出されるガンマ線をＳＰＥＣＴ装置等で捕らえて画像化する技術である。シンチグラフィーでは、経過観察や画像の相互比較等を行う目的で、画素値の正規化、すなわち、他の画像と比較可能な値への画素値の変換が行われる。正規化法には、術者が目視にてレベル調整を行うといった手動による方法と、自動的に正規化を行う方法とが存在する。【０００３】
手動による方法は、特別なプログラムを用いる必要が無いといったメリットを有するものの、結果が術者の技術レベルに左右されるといった欠点を有している。例えば、目視にてレベル調整する方法では、術者の感覚に結果が左右されやすく、その結果、ばらつきが生じやすい。また、胸椎にＲＯＩをとってその最大値を利用して正規化する方法もあるが、この方法は、作成したＲＯＩに異常集積が存在している場合には、利用する事ができない。【０００４】
下掲特許文献１及び非特許文献１には、シンチグラフィーの応用分野の一つである骨シンチグラフィー画像について、正規化の基準値を半自動的に決定する手法が開示されている。これらの文献には、下掲非特許文献２に開示される多重閾値法を用いて、骨シンチグラフィー画像のヒストグラムを解析し、放射性薬剤の異常集積に対応する画素値と正常集積に対応する画素値との境界値を特定し、その境界値を基準として骨シンチグラフィー画像を正規化することが、記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】ＷＯ２００７／０６２１３５
【非特許文献】
【０００６】
【非特許文献１】Shiraishi J et al. Development of a computer-aided diagnostic scheme for detection of interval changes in successive whole-body bone scans.  Med Phys 34:25-36, 2007.
【非特許文献２】M. L. Gger, K. Doi, and H. MacMahon, "lmagc featurc analysis and compurcr-aided diagnosis in digital radiography. 3. Automated detection of nodules in peripberal lung ficlds," Mcd. Phys. 15, 158-166 (1988).
【発明の概要】
【０００７】
しかし、本願発明者が検討したところ、特許文献１及び非特許文献１の手法は、正規化の結果に時折安定性を欠くことがあった。そこで本願発明者は、これらの従来技術よりも安定的に正規化を行う手法を開発することを志し、本願発明を為すに至った。
【０００８】
特許文献１及び非特許文献１の技術は、画素値が高い領域の範囲や値に依存して、正規化のための基準値を定めることが特徴である。すなわち、放射性医薬品の異常集積が生じている領域に基づいて、基準値を定めることが特徴である。本願発明者は、この特徴が、安定性を欠く原因になっていると考えた。
【０００９】
そして本願発明者は、正常集積領域に基づいて、正規化のための基準値を定めることにより、正規化の安定性を高めることができることを見出した。
【００１０】
本願発明の好適な実施形態の一例は、次のような方法である。この方法は、核医学の手法により得られる骨シンチグラフィー画像を正規化するための方法であって、
・  前記骨シンチグラフィー画像を表す画像データを読み込むことと；
・  前記画像データの画素値ヒストグラムを作成することと；
・  前記画素値ヒストグラムに基づいて、画素値に関する複数の閾値を設定することと；
・  前記設定した複数の閾値の各々に対して画素値平均値を計算することと；
・  前記計算した画素値平均値を、値の大きな順序に並べることと；
・  前記順序に並べられた前記画素値平均値の集合の少なくとも一部に基づいて、前記画像データを正規化するための基準値を決定することと；
を含み、前記基準値を決定することが、
・  前記順序に並べられた前記画素値平均値の集合のうち、前記画素値平均値が小さな領域を近似する直線を一つ決定することと；
・  前記直線に基づいて前記基準値を計算することと；
を含むことを特徴とする。
【００１１】
実施形態によっては、前記直線を一つ決定することが、前記画素値平均値の集合から、値の大きな画素値平均値をいくつか除去すると共に、残った部分を近似する近似直線を計算し、さらに、該近似直線と、該近似直線の計算に用いた画素値平均値の集合との相関係数を計算することを含んでもよい。
【００１２】
実施形態によっては、前記方法は、前記除去する画素値平均値の個数を変えて、それぞれの場合に前記近似直線及び前記相関係数を計算し、該計算した複数の相関係数の値に基づいて、前記直線を一つ決定することを含んでもよい。
【００１３】
実施形態によっては、前記画素値平均値が小さな領域は、前記順序に並べられた前記画素値平均値の集合の値の変化率に基づいて決定されてもよい。
【００１４】
実施形態によっては、前記基準値は、前記直線の切片に基づいて定められてもよい。
【００１５】
実施形態によっては、前記方法は、前記基準値に基づいて、前記画像データを正規化すると共に、該正規化した画像データを表示することを更に含んでもよい。
【００１６】
実施形態によっては、前記画素値ヒストグラムに基づいて画素値に関する複数の閾値を設定することを、多重閾値法を使って行ってもよい。
【００１７】
本願発明の好適な実施形態には、処理手段及び記憶手段を備えるシステムであって、前記記憶手段がプログラム命令を格納し、前記プログラム命令が、前記処理手段に実行されることにより、前記システムに、前記いずれかの方法を遂行させるように構成される、システムが含まれる。
【００１８】
本願発明の好適な実施形態には、システムの処理手段に実行されることにより、前記システムに、前記いずれかの方法を遂行させるように構成されるプログラム命令を備える、コンピュータプログラムが含まれる。
【００１９】
本願発明の好適な実施形態には、次のようなシステムが含まれる。このシステムは、核医学の手法により得られる骨シンチグラフィー画像を正規化するためのシステムであって、
・  前記骨シンチグラフィー画像を表す画像データを読み込む手段と；
・  前記画像データの画素値ヒストグラムを作成する手段と；
・  前記画素値ヒストグラムに基づいて、画素値に関する複数の閾値を設定する手段と；
・  前記設定した複数の閾値の各々に対して画素値平均値を計算する手段と；
・  前記計算した画素値平均値を、値の大きな順序に並べる手段と；
・  前記順序に並べられた前記画素値平均値の集合の少なくとも一部に基づいて、前記画像データを正規化するための基準値を決定する手段と；
を備え、前記基準値を決定する手段が、
・  前記順序に並べられた前記画素値平均値の集合のうち、前記画素値平均値が小さな領域を近似する直線を一つ決定する手段と；
・  前記直線に基づいて前記基準値を計算する手段と；
を備えることを特徴とする。
【００２０】
本発明の好適な具現化形態のいくつかを、特許請求の範囲に含まれる請求項に特定している。しかしこれらの請求項に特定される構成が、本願明細書及び図面に開示される新規な技術思想の全てを含むとは限らない。出願人は、現在の請求項に記載されているか否かに関わらず、本願明細書及び図面に開示される新規な技術思想の全てについて、特許を受ける権利を有することを主張するものであることを記しておく。
【図面の簡単な説明】
【００２１】
【図１】本明細書で開示される様々な処理を実行しうるハードウェアの例である、システム１００の主な構成を説明するための図である。
【図２】本明細書で開示される画像正規化処理を説明するためのフローチャートである。
【図３】骨シンチグラフィー画像の画素値ヒストグラムの例である。
【図４】ソートした画素値平均値をグラフ表示した例である。
【図５】ソートされた画素値平均値のデータを近似する直線を決定する処理の例を説明するためのフローチャートである。
【図６】図６Ａは、正規化前の画像データの表示例であり、図６Ｂは、正規化後の画像データの表示例である。
【図７】図７Ａは、ある骨シンチグラフィー画像についての、核医学専門医による手動正規化値と従来技術による自動正規化値との相関をみるための図であり、図７Ｂは、同じ骨シンチグラフィー画像についての、核医学専門医による手動正規化値と本発明による自動正規化値との相関をみるための図である。
【詳細説明】
【００２２】
以下、添付図面を参照しつつ、好適な実施例を用いて、本発明をより詳しく説明する。
【００２３】
図１は、本明細書で開示される様々な処理を実行しうるハードウェアの例である、システム１００の主な構成を説明するための図である。図１に描かれるように、システム１００は、ハードウェア的には一般的なコンピュータと同様であり、ＣＰＵ１０２，主記憶装置１０４，大容量記憶装置１０６，ディスプレイ・インターフェース１０７，周辺機器インタフェース１０８，ネットワーク・インターフェース１０９などを備えることができる。一般的なコンピュータと同様に、主記憶装置１０４としては高速なＲＡＭ（ランダムアクセスメモリ）を使用することができ、大容量記憶装置１０６としては、安価で大容量のハードディスクやＳＳＤなどを用いることができる。システム１００には、情報表示のためのディスプレイを接続することができ、これはディスプレイ・インターフェース１０７を介して接続される。またシステム１００には、キーボードやマウス、タッチパネルのようなユーザインタフェースを接続することができ、これは周辺機器インタフェース１０８を介して接続される。ネットワーク・インターフェース１０９は、ネットワークを介して他のコンピュータやインターネットに接続するために用いられることができる。
【００２４】
大容量記憶装置１０６には、オペレーティングシステム（ＯＳ）１１０や、本明細書で開示される特徴的な処理を提供するための、画像正規化プログラム１２０などが格納されていることができる。システム１００の最も基本的な機能は、ＯＳ１１０がＣＰＵ１０２に実行されることにより提供される。また、本明細書で開示される特徴的な処理は、画像正規化プログラム１２０に含まれるプログラム命令群の少なくとも一部がＣＰＵ１０２に実行されることにより提供される。よく知られているように、プログラムの実装形態には様々なものがあり、それらのバリエーションは全て、本願で開示される発明の範囲に含まれるものである。
【００２５】
大容量記憶装置１０６にはさらに、画像正規化プログラム１２０による正規化処理の対象となる骨シンチグラフィー画像データ１３０や、画像正規化プログラム１２０による処理の結果として生成される正規化基準値データ１３２、画像データ１３０を正規化した画像データ１３４なども格納されていることができる。
【００２６】
システム１００は、図１に描かれた要素の他にも、電源や冷却装置など通常のコンピュータシステムが備える装置と同様の構成を備えることができる。コンピュータシステムの実装形態には、記憶装置の分散・冗長化や仮想化、複数ＣＰＵの利用、ＣＰＵ仮想化、ＤＳＰなど特定処理に特化したプロセッサの使用、特定の処理をハードウェア化してＣＰＵに組み合わせることなど、様々な技術を利用した様々な形態のものが知られている。本願で開示される発明は、どのような形態のコンピュータシステム上に搭載されてもよく、コンピュータシステムの形態によってその範囲が限定されることはない。本明細書に開示される技術思想は、一般的に、（１）処理手段に実行されることにより、当該処理手段を備える装置またはシステムに、本明細書で説明される各種の処理を遂行させるように構成される命令を備えるプログラム、（２）当該処理手段が当該プログラムを実行することにより実現される装置またはシステムの動作方法、（３）当該プログラム及び当該プログラムを実行するように構成される処理手段を備える装置またはシステムなどとして具現化されることができる。前述のように、ソフトウェア処理の一部はハードウェア化される場合もある。
【００２７】
また、システム１００の製造販売時や起動時には、データ１３０〜１３４は、大容量記憶装置１０６の中に記憶されていない場合が多いことに注意されたい。骨シンチグラフィー画像データ１３０は、例えば周辺機器インタフェース１０８やネットワーク・インターフェース１０９を介して外部の装置からシステム１００に転送されるデータであってもよい。実施形態によっては、データ１３２，１３４は、プログラム１２０がＣＰＵ１０２に実行されることを通じて形成されたものであってもよい。また、プログラム１２０やＯＳ１１０の実装形態によっては、データ１３２及び／又は１３４は大容量記憶装置１０６に格納されず、主記憶装置１０４にしか格納されない場合もある。本願で開示される発明の範囲は、データ１３０〜１３４の存在の有無によって限定されるものではないことを、念のために記しておく。
【００２８】
次に、本明細書で開示される実施例において処理の対象となる骨シンチグラフィー画像データ１３０について簡単に説明しておく。骨シンチグラフィー画像データ１３０は、例えば、骨無機質の基本組成であるハイドロキシアパタイト結晶に吸着するような放射性医薬品を用いて施行されるＳＰＥＣＴ検査によって得られる、２次元画像データであることができる。より具体的には、例えば、放射性医薬品として99mTc-HMDPを被験者に静脈内投与し、体内から放射される放射線をＳＰＥＣＴ装置により検出して、その放射線カウント値に基づいて形成した２次元画像データであることができる。通常、これらの画像を構成する画素の各画素値は放射能カウント値に対応する値を有しており、すなわち各画素値は放射能の強度を表している。場合によっては、画像データ１３０は、時系列の２次元データや３次元画像データ、時系列の３次元データであることもある。
【００２９】
次に、図２以降を用いて、本明細書で開示される特徴的な処理の例をいくつか説明していく。
【００３０】
図２は、本明細書で開示される画像正規化処理２００を説明するためのフローチャートである。このフローチャートで紹介される処理は、例えば、画像正規化プログラム１２０のプログラム命令が、ＣＰＵ１０２に実行されることにより、システム１００が遂行する処理であることができる。
【００３１】
ステップ２０２は処理の開始を示す。ステップ２０４では、本実施例で処理対象となる骨シンチグラフィー画像データ１３０が読み込まれる。例えば、画像正規化プログラム１２０に含まれるプログラム命令群の指示に従って、ＣＰＵ１０２が、画像データ１３０の少なくとも一部を、大容量記憶装置１０６から主記憶装置１０４へとコピーする。
【００３２】
なお、実施形態によっては、骨シンチグラフィー画像データ１３０は、記憶媒体に格納されたものを専用の読み取り装置により読み取り、周辺機器インタフェース１０８を介してシステム１００に取り込まれたものであることができる。また別の実施形態においては、骨シンチグラフィー画像データ１３０は、搬送波に重畳されたデータ信号として、ネットワーク・インターフェース１０９を介して取り込まれたものであってよい。実施例によっては、取り込まれた骨シンチグラフィー画像データ１３０は、一旦大容量記憶装置１０６に格納し、主記憶装置１０４へコピーしても良いし、直接主記憶装置１０４に格納してそのまま後続の処理に供しても良い。
【００３３】
ステップ２０６では、画像データ１３０の画素値のヒストグラムを作成する。作成したヒストグラムのデータは、以降のステップの処理に用いたり、表示に用いたりするなどのために、主記憶装置１０４又は大容量記憶装置１０６に格納されてもよい。参考までに、図３に、ある実際の骨シンチグラフィー画像データから作成した画素値ヒストグラムを表示した例を載せた。ヒストグラムの横軸は階級値であり、画素値に関連する値である。縦軸は度数である。
【００３４】
ステップ２０８では、ステップ２０６で作成した画素値ヒストグラムを用いて、画素値に関する複数の閾値を設定する。好適な実施形態では、複数の閾値を設定するために多重閾値法を利用する。具体的には例えば次のように行う。まず、階級値ゼロからある階級値までのヒストグラム下の面積が、全ヒストグラム下面積の所定割合となるような階級値を複数求める。例えば、全ヒストグラム下面積を例えば１％に分割する様な階級値Ａ１を求める。次に、この分割割合を、例えば１％〜１００％まで、順次例えば１％刻みで増加させていき、それぞれに対応する階級値Ａ１，Ａ２，・・・Ａ１００を求める。そして、求めた階級値を閾値とする。このようにして複数の閾値を定め、特定の画素を抽出する手法を、多重閾値法という。
【００３５】
この例において、ヒストグラムのピーク付近では、階級値の僅かの違いで度数が大きく変動するであろうから、階級値Ａｉと階級値Ａｉ＋１の値の差は非常に小さくなるだろう。一方、階級値の大きなところでは、階級値が変わっても度数があまり変化せず、度数自体も小さいことから、ＡｉとＡｉ＋１の値の差は、比較的大きくなるだろう。
【００３６】
符号２１０−２１８で表されるループでは、ステップ２０８で定めた複数の閾値の各々に対して画素値平均値を計算する。
【００３７】
ステップ２１２では、画像データ１３０から、現在のループ位置における閾値（例えばＡｉ）以上の画素値を有する画素群を抽出する。ステップ２１４はオプションの処理なので、後で説明する。ステップ２１６では、ステップ２１２で抽出した画素群の画素値の平均値を計算する。ステップ２１２及びステップ２１６を、ステップ２０８で定めた複数の閾値の全てに対して行い、それぞれ対応する画素値平均値を計算する。計算結果は、以降のステップの処理に用いたり、表示に用いたりするなどのために、主記憶装置１０４又は大容量記憶装置１０６に格納されてもよい。
【００３８】
実施形態によっては、次に説明するステップ２１４を、ステップ２１２とステップ２１６の間に実行してもよい。
【００３９】
ステップ２１４では、ステップ２１２で抽出した画素群につき、隣り合う画素をつないでクラスターを形成する。そして、形成した各クラスターの画素数を調べ、当該画素数が所定の条件を満たさないクラスターに含まれる画素を特定する。そして、ステップ２１４を実行する場合、ステップ２１６における画素値平均値の計算において、所定の条件を満たさないクラスターに含まれる画素については、計算から除外することとしてもよい。一例として具体例を載せると、例えば、画素数が１０以下のクラスターに含まれる画素については、ステップ２１６における画素値平均値計算から除外することとしてもよい。他にも、例えば、画素数が所定の範囲にないクラスターに含まれる画素については、画素値平均値の計算から除外するといった処理を行ってもよい。
【００４０】
ステップ２２０では、符号２１０−２１８で表されるループで計算した複数の画素値平均値を、値の大きな順序に並べる。すなわち降順にソート処理を行う。ここで、ソートした画素値平均値をグラフ表示してもよい。グラフ表示の例を図４に載せた。このグラフの横軸の値は、大きさの順番を表す数値であり、例えば、グラフ中で横軸の値が１である点は、その画素値平均値が全てのデータ点の中で最も大きな値であることを示す。同様に、例えば、グラフ中で横軸の値が１０である点は、その画素値平均値が全てのデータ点の中で１０番目に大きな値であることを示している。なお、図４中に引かれた直線４１０については、図５のステップ５１４に関連して後述する。
【００４１】
ステップ２２２では、図４のようにグラフで表すことが可能な、ソートされた画素値平均値のデータを近似する直線を一つ決定する。ここで決定される直線は、以降のステップで、画像データ１３０を正規化するための基準値を計算する基礎となる。ステップ２２２で重要なことは、画素値平均値の値が小さな領域のデータを用いて近似直線を決定することである。本願発明の基本的な技術思想は、シンチグラフィー画像中で正常集積を表す領域のデータに基づいて、正規化のための基準値を定めることであるが、図４のグラフで正常集積を表す領域とは、画素値平均値の値が小さな領域であるからである。
【００４２】
複数のデータ点の近似直線を求める手法は既知であるが、近似直線を求める領域として、画素値平均値の値が小さな領域を自動的に使うようにするには、工夫が必要である。本発明の好適な実施形態の一つは、この課題を、図５に例示したようなアルゴリズムを使って解決している。
【００４３】
図５は、図４のようにグラフで表すことが可能な、ソートされた画素値平均値のデータを近似する直線を一つ決定する処理の例であって、画素値平均値の値が小さな領域を自動的に使って近似直線を決定する処理の例５００を説明するためのフローチャートである。ステップ５０２は処理の開始を示す。符号５０４−５１２はループ処理となり、処理が５０６に戻る度に、ステップ５０８で近似直線を計算するためのデータ範囲が変更される。ステップ５０８では、ステップ５０６で設定された範囲のデータを用いて近似直線が計算され、ステップ５１０では、ステップ５０８で求めた近似直線と、当該近似直線の計算に用いたデータ（すなわちステップ５０６で設定された範囲のデータ）との相関係数が計算される。
【００４４】
本実施例では、ステップ５０６における近似範囲の設定は次のように行われる。まず、処理５００が最初に５０６に進んだ時には、全てのデータ範囲を、近似直線を計算するための範囲として設定する。処理５００が次に５０６に戻ってきたときは、画素値平均値の値が最も大きなデータ点を除いて、近似直線を計算するための範囲を再設定する。更に次に処理５００が次に５０６に戻ってきたときは、前回の範囲の中で画素値平均値の値が最も大きなデータ点を除いて、近似直線を計算するための範囲を再設定する。
【００４５】
つまり、例えば、データ点の数が１００個あり、これらを図４のようにグラフ表示した場合、処理５００が最初に５０６に進んだ時に設定されるデータ範囲は、横軸の値で表すと、１−１００となる。処理５００が次に５０６に戻ってきたときに設定されるデータ範囲は、横軸の値で２−１００となる。処理５００が次に５０６に戻ってきたときに設定されるデータ範囲は、横軸の値で３−１００となる。ステップ５０８，５１０ではそれぞれ、設定された範囲毎に、近似直線と相関係数が計算される。
【００４６】
なお、実施形態に応じて、処理５００が最初に５０６に進んだ時に設定されるデータ範囲は、画素値平均値の値が最も大きなデータ点を必ずしも含む必要はなく、また、各データ範囲は、画素値平均値の値が最も小さなデータ点を含む必要は必ずしもないことを、念のために述べておく。更に、データ範囲に含まれるデータ点の個数には、下限値が定められていてもよい。例えば、上の例においても、処理５００が最初に５０６に進んだ時に設定されるデータ範囲は、横軸の値で例えば３−９０などであってもよく、符号５０４−５１２で表されるループを抜ける直前に設定されるデータ範囲は、横軸の値で例えば８１−９０などであってもよい。
【００４７】
従って、符号５０４−５１２で表されるループを抜けると、設定した範囲の情報と、対応する直線の情報と、対応する相関係数の情報とが、複数計算されている。これらの情報は、以降のステップの処理に用いたり、表示に用いたりするなどのために、主記憶装置１０４又は大容量記憶装置１０６に格納されてもよい。
【００４８】
ステップ５１４では、ループ５０４−５１２で計算された情報に基づいて、近似直線を一つ決定する。一つの例では、ループ５０４−５１２において、相関係数が最も大きくなった直線と決定してもよい。図５の直線４１０は、この手法によって決定した直線を図示したものである。
【００４９】
別の例では、ループ５０４−５１２において、相関係数が最初にピークを呈した直線を、ステップ５１４で決定する直線としてもよい。相関係数のピークを決定するにあたっては、適切なスムージング処理を適用した後に、ピークを決定してもよい。実施形態によっては、最初のピークではなく、２番目や３番目等、別の順番のピークに対応する直線を、ステップ５１４で決定する直線としてもよい。
【００５０】
さらに別の例では、ループ５０４−５１２において、相関係数が所定の値を上まわった時点で、そのときの直線をステップ５１４で決定する直線とするようにされてもよい。
【００５１】
いずれの手法を用いても、決定された直線の計算の基となったデータは、画素値平均値の値が小さな領域のデータになっている。つまり、図４のグラフでいうと、決定された直線の計算の基となったデータは、横軸の値が大きな領域のデータになっている。
【００５２】
ステップ５１６は、処理の終了を表す。
【００５３】
なお、ステップ２２２で近似直線を決定する手法は上記に限られない。例えば、実施形態によっては、単純に、図４のようにグラフ化したデータ点のうち、予め定めた範囲を、近似直線を計算する範囲と決めてしまってもよい。例えば、データ点が存在する横軸最大値の、例えば３０％−８０％の範囲を、近似直線を計算する範囲と決めてしまってもよい。
【００５４】
別の例として、データを図４のようにグラフ化した上で、隣接するデータ点の間の変化率を求め、その変化率の絶対値が所定の閾値を下回った点を、近似直線を計算する範囲の下限と定めるという手法もある。図４から見て取れるように、横軸の値が小さい領域では変化率の絶対値は大きく、横軸の値が大きくなるにつれて徐々に小さくなってゆく。そこで、変化率の絶対値の閾値を十分小さくとっておけば、画素値平均値の値が小さな領域のみが含まれるように、近似範囲を設定することが可能になる。
【００５５】
ステップ２２４では、ステップ２２２で決定した近似直線に基づいて、データ１３０を正規化するための基準値を決定する。
【００５６】
ある実施形態では、この基準値は、ステップ２２２で決定した近似直線のＹ切片の値である。従って、図４に例示した近似直線４１０の場合、この基準値は、およそ６５である。
【００５７】
別の実施形態では、この基準値は、当該近似直線の基になったデータ範囲における、縦軸最大値である。
【００５８】
基準値を決定する手法は、これらの例に限られない。しかし、どの場合であっても、ステップ２２２で決定した近似直線に基づいて、基準値を定めることが必要である。
【００５９】
決定した基準値は、以降のステップの処理に用いたり、表示に用いたりするなどのために、主記憶装置１０４又は大容量記憶装置１０６に格納されてもよい。
【００６０】
ステップ２２６では、ステップ２２４で決定した基準値を用いて、骨シンチグラフィー画像データ１３０を正規化する。本発明の実施形態において、基準値をどのように用いるかには特に制限はなく、様々なバリエーションが考えられる。一つの例であるが、骨シンチグラフィー画像データ１３０の各画素の画素値を例えば１−１０２４に正規化することとし、上記基準値が、そのスケールのなかで、特定の値（例えば３００等）になるように、全体の値を調節することとしてもよい。正規化した画像データは、更なる処理に用いたり、表示に用いたりするなどのために、主記憶装置１０４又は大容量記憶装置１０６に格納されてもよい。
【００６１】
ステップ２２８では、正規化した画像データの表示が行われる。表示例を図６Ａと６Ｂに載せた。図６Ａが正規化処理前の画像データ１３０を表示したものであり、図６Ｂが、画像データ１３０を上記の処理により正規化したものを表示したものである。これらの画像に細かな構造が現れていないのは、特許出願に用いることのできる画像データは二値しか使えないという制限のためである。本来は、グレースケールやカラーを用いた多値形式で、高精細な画像を表示することが可能である。
【００６２】
図７は、同じ画像データの同じ画素の画素値に対して、ある熟達した核医学専門医が手動で正規化を行った値を横軸に、自動で正規化した値を縦軸にとって、グラフにプロットしたものである。図７Ａで手動正規化値と比較されているものは、特許文献１及び非特許文献１の手法による自動正規化値であり、図７Ｂで手動正規化値と比較されているものは、本願発明の手法による自動正規化値である。データを直線近似して相関をみると、図７Ａの場合はR2 = 0.485であるに対して、図７Ｂの場合はR2 = 0.7063となり、本願発明による自動正規化値の方が、熟達した核医学専門医による正規化値との相関が高いことが示された。すなわち、本願発明による自動正規化値の方が、従来技術よりもばらつきが少ないことが示された。
【００６３】
以上、好適な実施例を用いて本発明を詳しく説明してきたが、上記の説明や添付図面は、本発明の範囲を限定する意図で提示されたものではなく、むしろ、法の要請を満たすために提示されたものである。本発明の実施形態には、ここで紹介されたもの以外にも、様々なバリエーションが存在する。例えば、明細書又は図面に示される各種の数値もいずれも例示であり、これらの数値は発明の範囲を限定する意図で提示されたものではない。
【００６４】
明細書又は図面に紹介した各種の実施例に含まれている個々の特徴は、その特徴が含まれることが直接記載されている実施例と共にしか使用できないものではなく、ここで説明された他の実施例や説明されていない各種の具現化例においても、組み合わせて使用可能である。特にフローチャートで紹介された処理の順番は、必ず紹介された順番で実行しなければならないわけではなく、実施者の好みや必要性に応じて、順序を入れ替えたり並列的に同時実行したり、さらに複数のブロックを一体不可分に実装したり、適当なループとして実行したりするように実装してもよい。これらのバリエーションは、全て、本願で開示される発明の範囲に含まれるものであり、処理の実装形態によって発明の範囲が限定されることはない。請求項に特定される処理の記載順も、処理の必須の順番を特定しているわけではなく、例えば処理の順番が異なる実施形態や、ループを含んで処理が実行されるような実施形態なども、請求項に係る発明の範囲に含まれるものである。現在の特許請求の範囲で特許請求がなされているか否かに関わらず、出願人は、発明の思想を逸脱しない全ての形態について、特許を受ける権利を有することを主張するものであることを記しておく。
【符号の説明】
【００６５】
１００    システム
１０４    主記憶装置
１０６    大容量記憶装置
１０７    ディスプレイ・インターフェース
１０８    周辺機器インタフェース
１０９    ネットワーク・インターフェース
１２０    画像正規化プログラム
１３０    骨シンチグラフィー画像データ
【図１】
2017010119.tif 000003
【図２】
2017010119.tif 000004
【図３】
2017010119.tif 000005
【図４】
2017010119.tif 000006
【図５】
2017010119.tif 000007
【図６】
2017010119.tif 000008
【図７】
2017010119.tif 000009
【国際調査報告】
2017010119.tif 000010
2017010119.tif 000011