(57)【特許請求の範囲】
【請求項１】
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）薬剤を投与した後の第１の時間帯に前記被験者の検査部位を撮影したＲＩ画像に係る第１の撮影画像データと、前記第１の時間帯から所定時間後の第２の時間帯に前記被験者の前記検査部位を撮影したＲＩ画像に係る第２の撮影画像データとを記憶する記憶手段と、
前記第１の撮影画像データから複数の第１の代表値を定め、各第１の代表値を前記複数の第１の代表値の中の最大値で正規化して第１正規化データを得るとともに、前記第２の撮影画像データから複数の第２の代表値を定め、各第２の代表値を前記複数の第２の代表値の中の最大値で正規化して第２正規化データを得る正規化手段と、
前記第１正規化データと前記第２正規化データとに基づいて、前記検査対象部位からの前記ＲＩ薬剤の洗い出し量を評価するための指標を算出する算出手段と、
前記算出手段によって算出された指標を出力する出力手段と、を備える画像解析装置。
【請求項２】
前記ＲＩ薬剤は、ＲＩ心筋血流製剤であることを特徴とする請求項１記載の画像解析装置。
【請求項３】
前記ＲＩ薬剤は、９９ｍＴｃ−ＳＥＳＴＡＭＩＢＩであることを特徴とする請求項１記載の画像解析装置。
【請求項４】
前記検査部位は心臓であることを特徴とする請求項１〜３のいずれかに記載の画像解析装置。
【請求項５】
前記正規化手段は、前記第１の撮影画像データ及び前記第２の撮影画像データのそれぞれについて、前記検査対象部位を複数のセグメントに分割し、各セグメントの代表値を前記複数の第１及び第２の代表値とし、セグメントごとに正規化して前記第１及び第２正規化データを得、
前記算出手段は、前記セグメントごとに正規化された前記第１及び第２正規化データに基づいて、各セグメントについて前記指標を算出することを特徴とする請求項１〜４のいずれかに記載の画像解析装置。
【請求項６】
前記指標は、前記第１正規化データから前記第２正規化データへの変化率を示す指標であることを特徴とする請求項１〜５のいずれかに記載の画像解析装置。
【請求項７】
前記指標は、以下の式により算出されるＷＯＩＮＤｘ（ｎ）である、請求項５に記載の画像解析装置。
【数６】
0005710166.tif 000007
Ｅ＿ＳＥＧ（ｎ）：早期像の各セグメントの代表値
Ｌ＿ＳＥＧ（ｎ）：後期像の各セグメントの代表値
ｎ：セグメントの識別子
Ｅ＿ＳＥＧ＿ＭＡＸ：早期像の各セグメントの代表値の最大値
Ｌ＿ＳＥＧ＿ＭＡＸ：後期像の各セグメントの代表値の最大値
【請求項８】
前記第１の撮影画像データは、前記第１の時間帯に撮影したＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｎｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像を変換したアレイデータであり、
前記第２の撮影画像データは、前記第２の時間帯に撮影したＳＰＥＣＴ画像を変換したアレイデータであり、
前記画像解析装置は、前記第１及び第２の撮影画像データに基づいてポーラーマップを表示させる表示手段を、さらに備え、
前記第１及び第２の正規化データは、前記ポーラーマップのセグメントごとに正規化されたデータであり、
前記算出手段は、前記ポーラーマップのセグメント単位で前記指標を算出することを特徴とする請求項１〜７のいずれかに記載の画像解析装置。
【請求項９】
画像解析装置が行う方法であって、
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）薬剤を投与した後の第１の時間帯に前記被験者の検査部位を撮影したＲＩ画像に係る第１の撮影画像データと、前記第１の時間帯から所定時間後の第２の時間帯に前記被験者の前記検査部位を撮影したＲＩ画像に係る第２の撮影画像データとを記憶するステップと、
前記第１の撮影画像データから複数の第１の代表値を定め、各第１の代表値を前記複数の第１の代表値の中の最大値で正規化して第１正規化データを得るとともに、前記第２の撮影画像データから複数の第２の代表値を定め、各第２の代表値を前記複数の第２の代表値の中の最大値で正規化して第２正規化データを得るステップと、
前記第１正規化データと前記第２正規化データとに基づいて、前記検査対象部位からの前記ＲＩ薬剤の洗い出し量を評価するための指標を算出するステップと、
前記算出された指標を出力するステップと、を行う画像解析方法。
【請求項１０】
画像解析装置のためのコンピュータプログラムであって、
前記画像解析装置に、
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）薬剤を投与した後の第１の時間帯に前記被験者の検査部位を撮影したＲＩ画像に係る第１の撮影画像データと、前記第１の時間帯から所定時間後の第２の時間帯に前記被験者の前記検査部位を撮影したＲＩ画像に係る第２の撮影画像データとを記憶するステップと、
前記第１の撮影画像データから複数の第１の代表値を定め、各第１の代表値を前記複数の第１の代表値の中の最大値で正規化して第１正規化データを得るとともに、前記第２の撮影画像データから複数の第２の代表値を定め、各第２の代表値を前記複数の第２の代表値の中の最大値で正規化して第２正規化データを得る正規化ステップと、
前記第１正規化データと前記第２正規化データとに基づいて、前記検査対象部位からの前記ＲＩ薬剤の洗い出し量を評価するための指標を算出する指標算出ステップと、
前記算出された指標を出力するステップと、を実行させるコンピュータプログラム。
【請求項１１】
前記正規化ステップは、前記第１の撮影画像データ及び前記第２の撮影画像データのそれぞれについて、前記検査対象部位を複数のセグメントに分割し、各セグメントの代表値を前記複数の第１及び第２の代表値とし、セグメントごとに正規化して前記第１及び第２正規化データを得、
前記指標算出ステップは、前記セグメントごとに正規化された前記第１及び第２正規化データに基づいて、各セグメントについて前記指標を算出することを特徴とする請求項１０に記載のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）画像の解析技術に関し、特に同一被験者の同一検査部位を異なる時刻に撮影したＲＩ画像を対比して、その検査対象部位から流出したＲＩ量を評価する技術に関する。
【背景技術】
【０００２】
ＲＩ薬剤としてＲＩ心筋血流製剤（例えば、９９ｍＴｃ−ＳＥＳＴＡＭＩＢＩ、あるいは９９ｍＴｃ−Ｔｅｔｒｏｆｏｓｍｉｎ）を被験者に投与すると、急性心筋梗塞治療後や重篤な狭心症患者であれば、心筋からＲＩ薬剤が流出する洗い出し（Ｗａｓｈｏｕｔ）の亢進が観測される。そこで、ＲＩ薬剤を被験者に投与して、その後の異なる時間帯に被験者の心臓をＳＰＥＣＴ撮影した複数の画像を用いて、心疾患などの診断に役立てる検査方法がある。例えば、ＲＩ薬剤投与後の第１の時間帯に撮影した早期像のカウント値と、第１の時間帯より後の第２の時間帯に撮影した後期像のカウント値とから求まる、心臓から洗い出されたＲＩ薬剤量の割合を示す洗い出し率（Ｗａｓｈｏｕｔ  ｒａｒｅ，ＷＲ）が知られている（例えば非特許文献１）。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】A.M.Peters. A Unified Approach to Quantification by Kinetic Analysis in NuclearMedicine. J Nucl Med. 1993;34:706-13.
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、ＷＲは定量的な指標であるが、算出された指標値は早期像及び後期像の画像を目視で対比したときに行う定性的な評価とは大きく乖離することが多かった。例えば、ＷＲの値が同じであっても、目視評価では被験者によって正常と判断されたり、異常と判断されたりする場合があった。これは、撮像時間間隔が異なることや、被験者ごとに洗い出されるＲＩ薬剤量に個人差があるからである。そのため、臨床現場においては、早期像及び後期像の画像を目視で対比した定性的な評価と良好に一致する定量的な評価指標が求められていた。
【０００５】
そこで、本発明の目的は、早期像及び後期像を目視で対比した定性的な評価と一致する定量的な評価指標を提供することである。
【課題を解決するための手段】
【０００６】
本発明の一つの実施態様に従う画像解析装置は、被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）薬剤を投与した後の第１の時間帯に前記被験者の検査部位を撮影したＲＩ画像に係る第１の撮影画像データと、前記第１の時間帯から所定時間後の第２の時間帯に前記被験者の前記検査部位を撮影したＲＩ画像に係る第２の撮影画像データとを記憶する記憶手段と、前記第１の撮影画像データを正規化して第１正規化データを得るとともに、前記第２の撮影画像データを正規化して第２正規化データを得る正規化手段と、前記第１正規化データと前記第２正規化データとに基づいて、前記検査対象部位からの前記ＲＩ薬剤の洗い出し量を評価するための指標を算出する算出手段と、前記算出手段によって算出された指標を出力する出力手段と、を備える。
【０００７】
好適な実施態様では、前記ＲＩ薬剤は、ＲＩ心筋血流製剤であってもよい。
【０００８】
好適な実施態様では、前記ＲＩ薬剤は、９９ｍＴｃ−ＳＥＳＴＡＭＩＢＩであってもよい。
【０００９】
好適な実施態様では、前記検査部位は心臓であってもよい。
【００１０】
好適な実施態様では、前記正規化手段は、前記第１の撮影画像データ及び前記第２の撮影画像データのそれぞれについて、前記検査対象部位を複数のセグメントに分割し、各セグメントの代表値を用いてセグメントごとに正規化して、前記第１及び第２正規化データを得る。そして、前記算出手段は、前記セグメントごとに正規化された前記第１及び第２正規化データに基づいて、各セグメントについて前記指標を算出するようにしてもよい。
【００１１】
好適な実施態様では、前記指標は、前記第１正規化データから前記第２正規化データへの変化率を示す指標であってもよい。
【００１２】
好適な実施態様では、前記指標は、以下の式により算出されるＷＯＩＮＤｘ（ｎ）であってもよい。
【数１】
0005710166.tif 000002
Ｅ＿ＳＥＧ（ｎ）：早期像の各セグメントの代表値
Ｌ＿ＳＥＧ（ｎ）：後期像の各セグメントの代表値
ｎ：セグメントの識別子
Ｅ＿ＳＥＧ＿ＭＡＸ：早期像の各セグメントの代表値の最大値
Ｌ＿ＳＥＧ＿ＭＡＸ：後期像の各セグメントの代表値の最大値
【００１３】
好適な実施態様では、前記第１の撮影画像データは、前記第１の時間帯に撮影したＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｎｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像を変換したアレイデータであり、前記第２の撮影画像データは、前記第２の時間帯に撮影したＳＰＥＣＴ画像を変換したアレイデータであり、前記画像解析装置は、前記第１及び第２の撮影画像データに基づいてポーラーマップを表示させる表示手段を、さらに備えてもよい。そして、前記第１及び第２の正規化データは、前記ポーラーマップのセグメントごとに正規化されたデータであり、前記算出手段は、前記ポーラーマップのセグメント単位で前記指標を算出するようにしてもよい。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施形態に係る画像解析装置１の構成図である。
【図２】アレイデータの生成手順を示すフローチャートである。
【図３】心臓の断層画像を模式的に示す図である。
【図４】ＷＯＩＮＤｘ算出の全体処理手順を示すフローチャートである。
【図５】正規化処理の具体的な手順を示すフローチャートである。
【図６】アレイデータから正規化セグメントデータへの変換の説明図である。
【図７】ポーラーマップ及びＷＯＩＮＤｘの表示例である。
【発明を実施するための形態】
【００１５】
以下、本発明の一実施形態に係る画像解析装置について、図面を参照して説明する。
【００１６】
図１は、本実施形態に係る画像解析装置１の構成図を示す。本実施形態に係る画像解析装置１は、被験者にＲＩ薬剤を投与した後の異なる時間帯に撮影された、被験者の検査部位の複数のＲＩ画像を解析して、検査対象部位から流出したＲＩ量を評価する指標を算出する。本実施形態では、特に、ＲＩ薬剤として９９ｍＴｃ−ＳＥＳＴＡＭＩＢＩを用いて、図示しないＳＰＥＣＴ撮像装置で撮影した心臓のＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｎｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像を用いて解析を行う。
【００１７】
この画像解析装置１は、例えば汎用的なコンピュータシステムにより構成され、以下に説明する画像解析装置１の個々の構成要素または機能は、例えば、コンピュータプログラムを実行することにより実現される。このコンピュータプログラムは、コンピュータ読み取り可能な記録媒体に格納可能である。
【００１８】
画像解析装置１は、同図に示すように、画像データ記憶部１１と、アレイデータ生成部１３と、アレイデータ記憶部１５と、正規化処理部１６と、指標算出部１７と、指標記憶部１８と、表示処理部１９とを備える。
【００１９】
画像データ記憶部１１には、図示しないＳＰＥＣＴ撮像装置で撮影した心臓のＳＰＥＣＴ画像のデータが記憶されている。ここでは、ＳＰＥＣＴ画像のデータとしては、ＲＩ薬剤を投与した後の第１の時間帯に被験者の心臓を撮影した早期像データ１１１と、第１の時間帯から所定時間後の第２の時間帯に同じ被験者の心臓を撮影した後期像データ１１３とが記憶されている。早期像は、例えば、ＲＩ薬剤投与１時間後に撮影したＳＰＥＣＴ画像であり、後期像は、例えば、ＲＩ薬剤投与６時間後に撮影したＳＰＥＣＴ画像である。
【００２０】
早期像データ１１１及び後期像データ１１３は、いずれも、３次元のボクセルデータである。各ボクセルの値は、ＳＰＥＣＴ撮像装置でカウントされたカウント値である。
【００２１】
アレイデータ生成部１３は、画像データ記憶部１１に記憶されている早期像データ１１１及び後期像データ１１３に基づいて、ポーラーマップ（ブルズアイ）を生成するための早期像アレイデータ１５１及び後期像アレイデータ１５３をそれぞれ生成する。
【００２２】
ここで、図２はアレイデータ生成部１３が行うアレイデータの生成手順を示すフローチャートである。同図に従って、アレイデータ生成部１３の処理手順を詳細に説明する。アレイデータ生成部１３は、早期像データ１１１及び後期像データ１１３に対してそれぞれ以下の処理を適用して早期像アレイデータ１５１及び後期像アレイデータ１５３を生成する。
【００２３】
まず、アレイデータ生成部１３は、処理対象となる画像（早期像データ１１１または後期像データ１１３）の画像データを読み込む（Ｓ１１）。
【００２４】
アレイデータ生成部１３は、対象の画像における心臓領域の範囲を選択し、その範囲を処理範囲として設定する（Ｓ１２）。この心臓領域の選択は、オペレータの指示に従って行っても良いし、アレイデータ生成部１３が自動的に心臓領域を検出して行っても良い。
【００２５】
図３は、心臓の断層画像を模式的に示す。同図Ａは垂直長軸画像の模式図であり、同図Ｂは短軸画像の模式図である。オペレータは、例えば、表示処理部１９によって同図Ａ，Ｂのような画像が表示装置２２に表示されているときに、入力装置２１を用いてそれぞれの心臓領域を選択する。例えば、オペレータは、同図Ａの垂直長軸画像を用いて、心基部ｂｅ及び心尖部ａｅを指定する。また、同図Ｂの短軸画像を用いて心臓の中心ｃｅ及び半径ｒｅ（外周）を設定する。
【００２６】
図２に戻ると、アレイデータ生成部１３は、短軸画像の画像データに基づいて、心基部ｂｅと心尖部ａｅの間の長さに応じて、短軸画像を所定枚数（ここでは１６枚）に正規化する（Ｓ１３）。
【００２７】
アレイデータ生成部１３は、短軸画像において、図３Ｂに示すように、心臓の中心ｃｅから３時方向の半径ｒｅを起点として、７．５度おき、４８方向の半径ｒｅについて、それぞれの半径ｒｅ上で最も大きい画素値を抽出する（Ｓ１４）。アレイデータ生成部１３は、この処理を、正規化によって得られた１６枚の短軸画像のそれぞれについて行う（Ｓ１５）。
【００２８】
これにより、４８×１６のアレイデータが抽出される。このアレイデータに基づいて、後述するポーラーマップが表示される。
【００２９】
アレイデータ生成部１３により生成された早期像アレイデータ１５１（第１の撮影画像データ）及び後期像アレイデータ１５３（第２の撮影画像データ）は、それぞれアレイデータ記憶部１５に保存される。
【００３０】
正規化処理部１６は、早期像アレイデータ１５１及び後期像アレイデータ１５３をそれぞれ正規化する。すなわち、正規化処理部１６は早期像アレイデータ１５１（第１回撮影画像データ）を正規化して早期像の正規化セグメントデータ（第１正規化データ）を得るとともに、後期像アレイデータ１５３（第２の撮影画像データ）を正規化して後期像の正規化セグメントデータ（第２正規化データ）を得る。正規化処理部１６は、早期像アレイデータ１５１（第１の撮影画像データ）及び後期像アレイデータ１５３（第２の撮影画像データ）のそれぞれについて、心臓を複数のセグメントに分割し、各セグメントの代表値を用いてセグメントごとに正規化して、早期像及び後期像の正規化セグメントデータ（第１及び第２正規化データ）を得る。早期像及び後期像の正規化セグメントデータは、心臓を複数のセグメントに分割して得られたセグメントごとに正規化されたデータである。
【００３１】
指標算出部１７は、早期像の正規化セグメントデータ（第１正規化データ）と後期像の正規化セグメントデータ（第２正規化データ）とに基づいて、被験者の心臓からのＲＩ薬剤の洗い出し量を評価するための指標を算出する。指標算出部１７は、セグメントごとに、各セグメントについてＲＩ薬剤の洗い出し量を評価するための指標を算出する。ここで算出する指標をＷａｓｈｏｕｔＩＮＤＥＸ（ＷＯＩＮＤｘ）とよぶ。
【００３２】
図４は、ＷＯＩＮＤｘ算出の全体処理手順を示すフローチャートである。同図を用いて、ＷＯＩＮＤｘ算出の全体処理手順を説明する。
【００３３】
まず、正規化処理部１６は、早期像アレイデータ１５１に対して、後述する正規化処理を適用して正規化を行い、早期像の正規化セグメントデータを生成する（Ｓ２１）。
【００３４】
同様に正規化処理部１６は、後期像アレイデータ１５３に対して、後述する正規化処理を適用して正規化を行い、後期像の正規化セグメントデータを生成する（Ｓ２３）。なお、早期像の正規化と後期像の正規化は、いずれを先に行っても良い。
【００３５】
そして、指標算出部１７は、ステップＳ２１において正規化された早期像の正規化セグメントデータと、ステップＳ２３において正規化された後期像の正規化セグメントデータとを用いて、ＷＯＩＮＤｘを算出する（Ｓ２５）。以下、ステップＳ２１〜Ｓ２５の各処理を詳細に説明する。
【００３６】
まず、ステップＳ２１及びＳ２３の正規化処理の詳細を説明する。ステップＳ２１及びＳ２３は、対象となるデータが異なるのみで、処理内容は共通である。図５は、この正規化処理の具体的な手順を示すフローチャートである。図６は、アレイデータから正規化セグメントデータへの変換の説明図である。以下、図５及び図６を用いて正規化処理について説明する。
【００３７】
正規化処理部１６は、対象のアレイデータａｒ［ｉ，ｊ］（早期像アレイデータ１５１または後期像アレイデータ１５３）において、最大値が１００になるように、アレイデータ内のそれぞれの値を正規化する（％Ｕｐｔａｋｅ）（Ｓ３１）。
【００３８】
次に、ステップＳ３１で正規化された対象のアレイデータｎ＿ａｒ［ｉ，ｊ］の各値を、ポーラーマップを作成する際の１７セグメントに分け、各セグメントの代表値ＳＥＧ（ｎ）（ｎはセグメントの識別子）からなるセグメントデータ（早期像セグメントデータまたは後期像セグメントデータ）を求める（Ｓ３３）。ここで、対象のアレイデータの各値がどのセグメントに属する値であるかは予め定められている。なお、セグメント数は必ずしも１７セグメントでなくても良く、任意のセグメント数でよい。また、各セグメントの代表値ＳＥＧ［ｎ］は、例えば、正規化されたアレイデータｎ＿ａｒ［ｉ，ｊ］の平均値、最大値、最小値、中間値あるいは最頻値などの各セグメントを代表する値でよい。
【００３９】
正規化処理部１６は、ステップＳ３３で定めたセグメントデータＳＥＧ（ｎ）をさらに正規化して、正規化セグメントデータＳＶ［ｎ］を生成する（Ｓ３５）。例えば、正規化処理部１６は、以下の式（１）により、各セグメントの代表値ＳＥＧ（ｎ）のなかの最大値（ＳＥＧ＿ＭＡＸ）が１００になるように、１７セグメントの代表値ＳＥＧ（ｎ）を正規化して、各セグメントの正規化された代表値ＳＶ（ｎ）を生成する（セグメントごとの％Ｕｐｔａｋｅ）。
【数２】
0005710166.tif 000003
【００４０】
本実施形態では、ステップＳ３１で％Ｕｐｔａｋｅに正規化したものをさらにステップＳ３５で正規化している。これにより、比較対象となるデータのベースが揃うので、被験者ごとの個人差によるばらつきを抑制する効果ある。なお、このように正規化を繰り返すことは必須ではない。例えば、２回の正規化処理のうちステップＳ３１を省略し、ステップＳ３５のみを行うようにしても良い。
【００４１】
ここで、早期像のセグメントデータをＥ＿ＳＥＧ（ｎ）、後期像のセグメントデータをＬ＿ＳＥＧ（ｎ）で表す。さらに、早期像のセグメントデータＥ＿ＳＥＧ（ｎ）の中の最大値をＥ＿ＳＥＧ＿ＭＡＸとし、後期像のセグメントデータＬ＿ＳＥＧ（ｎ）の中の最大値をＬ＿ＳＥＧ＿ＭＡＸとする。このとき、早期像の正規化セグメントデータＥ＿ＳＶ（ｎ）、及び後期像の正規化セグメントデータＬ＿ＳＶ（ｎ）は、それぞれ、以下の式（２）、及び式（３）で定義される。
【数３】
0005710166.tif 000004
【数４】
0005710166.tif 000005
【００４２】
指標算出部１７は、上記の早期像の正規化セグメントデータＥ＿ＳＶ（ｎ）及び後期像の正規化セグメントデータＬ＿ＳＶ（ｎ）を用いて、以下の式（４）により、セグメント別に洗い出し指標ＷＯＩＮＤｘ（ｎ）を算出する。
【数５】
0005710166.tif 000006
【００４３】
ＷＯＩＮＤｘ（ｎ）は、セグメント別の早期像の正規化セグメントデータＥ＿ＳＶ（ｎ）（第１正規化データ）から後期像の正規化セグメントデータＬ＿ＳＶ（ｎ）（第２正規化データ）への変化率を示す指標である。すなわち、ＷＯＩＮＤｘ（ｎ）は、セグメントごとに、早期像を基準とし、早期像から後期像へカウント値が変化した変化率を表したものである。その際に、早期像及び後期像のカウント値は、それぞれの画像内で正規化した値を用いて比較しているので、ＲＩが放出する放射能の減衰及び正常領域の正常な洗い出し（生理的な洗い出し）の影響を除去することができる。
【００４４】
指標算出部１７が算出したＷＯＩＮＤｘは、指標記憶部１８に保存される。
【００４５】
表示処理部１９は、早期像アレイデータ１５１及び後期像アレイデータ１５３（前記第１及び第２の撮影画像データ）に基づいてポーラーマップを生成して、表示装置２２に表示させる。さらに表示処理部１９は、指標記憶部１８に記憶されているＷＯＩＮＤｘを表示装置２２に表示させる。それらの表示例を図７に示す。
【００４６】
同図Ａは早期像アレイデータ１５１に基づく早期像のポーラーマップであり、同図Ｂは後期像アレイデータ１５３に基づく後期像のポーラーマップであり、同図ＣはＷＯＩＮＤｘをセグメント別に表示させた表示例である。ＷＯＩＮＤｘは、数値が大きいほど、心臓から洗い出されたＲＩの量が多いことを示す。
【００４７】
この結果、同図Ａ、Ｂを見た医師の評価によると、洗い出し量の多いセグメントが図中にＡで示す領域である。そして、同図Ｃによれば、この領域Ａ内のセグメントのＷＯＩＮＤｘが、他のセグメントよりも大きな値となっていることがわかる。すなわち、ＷＯＩＮＤｘは、早期像及び後期像を目視で対比した定性的な評価と良好に一致する。
【００４８】
一般に、検査を行う施設が異なると、早期像と後期像との時間間隔も異なることが多い。そのため、ＲＩが放出する放射能の減衰の程度及び正常領域から洗い出される正常なＲＩ量が異なるため、異なる施設での検査結果の比較が困難であった。しかしながら、本実施形態に係るＷＯＩＮＤｘは、時間経過によるＲＩが放出する放射能の減衰及び正常領域の正常な洗い出しの影響を除去しているので、異なる施設での検査結果であっても比較が可能である。つまり、異なる施設での検査結果を同じ閾値を用いて評価することができる。
【００４９】
また、正規化において％Ｕｐｔａｋｅを用いることにより、ＲＩの核種の違いがあっても正常領域と異常領域との差を客観的、かつ定量的に評価することが可能となる。従って、％Ｕｐｔａｋｅによって正規化したデータを用いて算出したＷＯＩＮＤｘの臨床的用途は広いと考えられる。
【００５０】
さらに、本実施形態では、ＷＯＩＮＤｘを算出する際に、早期像のセグメントデータの中の最大値Ｅ＿ＳＥＧ＿ＭＡＸ、及び後期像のセグメントデータの最大値Ｌ＿ＳＥＧ＿ＭＡＸがいずれも１００になるように正規化した正規化セグメントデータＥ＿ＳＶ（ｎ）、Ｌ＿ＳＶ（ｎ）を用いている。このため、早期像と後期像とを比較する際のベースが揃うので、後期像において、早期像と比較したときに洗い出しが亢進しているセグメントを定量的に評価することができる。
【００５１】
また、式（４）では、分母及び分子ともにＥ＿ＳＥＧ＿ＭＡＸまたはＬ＿ＳＥＧ＿ＭＡＸで正規化した値を用いているので（式（２）、（３）参照）、洗い出し量が比較的少ない正常領域では、ＷＯＩＮＤｘは０に近い値となる。そのため、ＷＯＩＮＤｘを用いれば、正常領域と洗い出しが亢進している領域とは明瞭に区別することができる。これは、正常領域と異常領域との値の差がはっきりしていないＷＲとは対照的である。
【００５２】
さらに、後期像は早期像と同じ条件で収集するとカウント値が少なく、異常部位の過大評価が問題になる。そこで、後期像の画質を維持するためには後期像の収集時間を長くする必要があった。本実施形態のＷＯＩＮＤｘでは、カウント値を正規化しているので、このような収集時間の長短の影響を受けない。
【００５３】
上述した本発明の実施形態は、本発明の説明のための例示であり、本発明の範囲をそれらの実施形態にのみ限定する趣旨ではない。当業者は、本発明の要旨を逸脱することなしに、他の様々な態様で本発明を実施することができる。
【００５４】
例えば、検査対象部位は心臓以外の臓器でも良い。さらに、本実施形態ではＳＰＥＣＴ画像に基づいた画像解析について説明したが、本発明はＰＥＴ（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ）画像に対して適用することもできる。
【符号の説明】
【００５５】
１  画像解析装置
１１  画像データ記憶部
１３  アレイデータ生成部
１５  アレイデータ記憶部
１６  正規化処理部
１７  指標算出部
１８  指標記憶部
１９  表示処理部
１１１  早期像データ
１１３  後期像データ
１５１  早期像アレイデータ
１５３  後期像アレイデータ
【図１】
0005710166.tif 000008
【図２】
0005710166.tif 000009
【図３】
0005710166.tif 000010
【図４】
0005710166.tif 000011
【図５】
0005710166.tif 000012
【図６】
0005710166.tif 000013
【図７】
0005710166.tif 000014