(57)【要約】【課題】患者に将来イベントが生じる可能性を示す指標を精度よく算出する。【解決手段】情報処理装置１００は、心機能診断薬を被験者に投与して得られるＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）を記憶した被験者データ記憶部１１３と、被験者データ記憶部１１３から被験者のＨ／Ｍ比を取得し、複数の被験者に対する５年後の生死の実績に応じて定められた関数に被験者データ記憶部１１３から取得したＨ／Ｍ比を代入することにより、被験者に所定のイベントが生じる可能性を示す指標を算出する生存率算出部１１９とを備える。【選択図】図７
2014215222.tif 000002
【特許請求の範囲】
【請求項１】
第１の被験者に心機能診断薬を投与して得たＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）を含む被験者データを記憶した記憶部を有する情報処理装置に実行させるためのコンピュータプログラムであって、
前記記憶部から前記被験者データの前記Ｈ／Ｍ比を取得するステップと、
複数の第２の被験者に対する所定のイベントの成否の実績に応じて定められた関数に、前記記憶部から取得した前記Ｈ／Ｍ比を代入することにより、前記第１の被験者に前記所定のイベントが生じる可能性を示す指標を算出するステップと、を前記情報処理装置に実行させるためのコンピュータプログラム。
【請求項２】
前記被験者データは、前記Ｈ／Ｍ比に加えて、被験者の年齢、性別、心不全の重症度を示す指標、及び左室機能を示す指標、のうちの少なくとも一つ以上をさらに含み、
前記算出ステップでは、前記関数に、前記Ｈ／Ｍ比に加えて、前記年齢、前記性別、前記心不全の重症度を示す指標、及び前記左室機能を示す指標、のうちの少なくとも一つ以上とを代入することにより、前記第１の被験者に前記所定のイベントが生じる可能性を示す指標を算出する請求項１に記載のコンピュータプログラム。
【請求項３】
前記関数には、前記第１の被験者の、前記Ｈ／Ｍ比、前記年齢、前記性別、前記心不全の重症度を示す指標、及び前記左室機能を示す指標のみが代入される請求項１又は請求項２に記載のコンピュータプログラム。
【請求項４】
前記所定のイベントが生じる可能性を示す指標は、５年生存率である請求項１から請求項３の何れか一項に記載のコンピュータプログラム。
【請求項５】
前記心機能診断薬は、ＭＩＢＧ（１２３Ｉ−３−ｉｏｄｏｂｅｎｚｙｌｇｕａｎｉｄｉｎｅ）である、請求項１から請求項４の何れか一項に記載のコンピュータプログラム。
【請求項６】
前記第１の被験者の前記被験者データに含まれるＨ／Ｍ比は、第１の撮像環境の所定の撮像条件におけるＨ／Ｍ比を、撮像環境に依存しない前記撮像条件から定まるＨ／Ｍ比へ変換する変換関数により変換されたＨ／Ｍ比であって、
前記変換関数は、
前記第１の撮像環境でファントムを前記所定の撮像条件で撮像したファントム画像のデータであるファントムデータと、心臓ＲＯＩ及び縦隔ＲＯＩを有するデジタルファントムのデータであるデジタルファントムデータとに基づいて、前記ファントム画像に前記デジタルファントムの位置合わせを行って、前記心臓ＲＯＩ及び前記縦隔ＲＯＩを設定し、
前記心臓ＲＯＩ及び前記縦隔ＲＯＩが設定されている前記ファントム画像のファントムデータに基づいて、前記第１の撮像環境における前記ファントムのＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）である第１のファントムＨ／Ｍ比を算出し、
前記第１のファントムＨ／Ｍ比及び前記撮像条件から定まるＨ／Ｍ比から導出された、前記変換関数である、請求項１から請求項５の何れか一項に記載のコンピュータプログラム。
【請求項７】
前記第１の被験者の前記被験者データに含まれるＨ／Ｍ比は、第１の撮像環境におけるＨ／Ｍ比を、第２の撮像環境におけるＨ／Ｍ比へ変換する変換関数によって変換されたＨ／Ｍ比であって、
前記変換関数は、
前記第１の撮像環境でファントムを撮像した第１のファントム画像及び前記第２の撮像環境で前記ファントムを撮像した第２のファントム画像のデータであるファントムデータと、心臓ＲＯＩ及び縦隔ＲＯＩを有するデジタルファントムのデータであるデジタルファントムデータとに基づいて、前記第１のファントム画像に前記デジタルファントムの位置合わせを行って、前記心臓ＲＯＩ及び前記縦隔ＲＯＩを設定し、かつ、前記第２のファントム画像に前記デジタルファントムの位置合わせを行って、前記心臓ＲＯＩ及び前記縦隔ＲＯＩを設定し、
前記心臓ＲＯＩ及び前記縦隔ＲＯＩが設定されている前記第１のファントム画像のファントムデータに基づいて、前記第１の撮像環境における前記ファントムのＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）である第１のファントムＨ／Ｍ比を算出し、かつ、前記心臓ＲＯＩ及び前記縦隔ＲＯＩが設定されている前記第２のファントム画像のファントムデータに基づいて、前記第２の撮像環境における前記ファントムのＨ／Ｍ比である第２のファントムＨ／Ｍ比を算出し、
前記第１のファントムＨ／Ｍ比及び前記第２のファントムＨ／Ｍ比に基づいて導出された、前記変換関数である請求項１から請求項５の何れか一項に記載のコンピュータプログラム。
【請求項８】
前記コンピュータプログラムは、
前記変換関数を導出するステップ、を前記情報処理装置にさらに実行させる、請求項６又は請求項７に記載のコンピュータプログラム。
【請求項９】
第１の被験者に心機能診断薬を投与して得たＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）を含む被験者データを記憶する記憶部と、
前記記憶部から前記被験者データの前記Ｈ／Ｍ比を取得する取得手段と、
複数の第２の被験者に対する所定のイベントの成否の実績に応じて定められた関数に、前記記憶部から取得した前記Ｈ／Ｍ比を代入することにより、前記第１の被験者に前記所定のイベントが生じる可能性を示す指標を算出する算出手段と、を備える情報処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心不全の患者に対し所定のイベントが生じる可能性を示す指標を算出する技術に関する。
【背景技術】
【０００２】
従来より、心不全の患者に対する検査結果に基づき、将来起こりうるイベントに関する指標を算出するための関数を導出し、当該関数により患者に将来イベントが生じる可能性を示す指標を算出する処理が提案されている。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかし、従来提案されている関数では、将来イベントが生じる可能性を示す指標を精度よく算出することができなかった。
【０００４】
本発明の目的は、患者に将来イベントが生じる可能性を示す指標を精度よく算出することである。
【課題を解決するための手段】
【０００５】
本発明の一つの実施態様に従うコンピュータプログラムは、第１の被験者に心機能診断薬を投与して得たＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）を含む被験者データを記憶した記憶部を有する情報処理装置に実行させるためのコンピュータプログラムであって、前記記憶部から前記被験者データの前記Ｈ／Ｍ比を取得するステップと、複数の第２の被験者に対する所定のイベントの成否の実績に応じて定められた関数に、前記記憶部から取得した前記Ｈ／Ｍ比を代入することにより、前記第１の被験者に前記所定のイベントが生じる可能性を示す指標を算出するステップと、を前記情報処理装置に実行させる。
【０００６】
また、本発明の一つの実施態様に従うコンピュータプログラムは、複数の第２の被験者の被験者データと、前記複数の第２の被験者に対する所定のイベントの成否の実績を示すイベントデータとを記憶した記憶部を有する情報処理装置に実行させるためのコンピュータプログラムであって、前記被験者データは、各被験者に心機能診断薬を投与して得たＨｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比（Ｈ／Ｍ比）を含み、前記記憶部から前記複数の第２の被験者の被験者データ、及びイベントデータを取得するステップと、前記記憶部から取得した、前記複数の第２の被験者の被験者データ及びイベントデータに基づいて導出され、第１の被験者に心機能診断薬を投与して得たＨ／Ｍ比を代入することにより、前記第１の被験者に将来前記所定のイベントが生じる可能性を示す指標を算出するための関数を導出するステップと、を前記情報処理装置に実行させる。
【図面の簡単な説明】
【０００７】
【図１】本発明の一実施形態に係る情報処理装置１の全体構成を示す。
【図２】ファントム５の一例を示す。
【図３】ファントム５の４層の各アクリル板のそれぞれの平面図を示す。
【図４】変換関数を示す。
【図５】変換関数を算出する手順を示すフローチャートを示す。
【図６】Ｈ／Ｍ比を変換する手順を示すフローチャートを示す。
【図７】本発明の他の実施形態に係る情報処理装置１００の全体構成を示す。
【図８】５年生存率関数を導出する手順を示すフローチャートを示す。
【図９】５年生存率を算出する手順を示すフローチャートを示す。
【発明を実施するための形態】
【０００８】
以下、本発明の第１の実施形態に係る情報処理装置について、図面を参照して説明する。
【０００９】
本実施形態に係る情報処理装置は、異なる撮像環境で撮像された画像から算出された診断指標であるＨ／Ｍ比を互いに比較可能にするための処理を行う。
【００１０】
本実施形態では、心機能診断薬が投与された被験者を、所定の撮像装置を含むある撮像環境で撮像したプラナ画像、及びファントムに同じ薬剤を保持させて撮像したプラナ画像が用いられる。以下の実施形態では、心機能診断薬としてＭＩＢＧ（１２３Ｉ−３−ｉｏｄｏｂｅｎｚｙｌｇｕａｎｉｄｉｎｅ）を用いた場合について説明する。
【００１１】
図１は、本実施形態に係る情報処理装置１の全体構成を示す。同図に示すように、情報処理装置１は、情報処理装置本体１０に入力装置２、出力装置３及び入出力装置４の一つ以上が接続されている。情報処理装置本体１０は、例えば、プロセッサ及びメモリを備えた汎用的なコンピュータシステムにより構成され、以下に説明する情報処理装置本体１０内の個々の構成要素または機能は、例えば、コンピュータプログラムを実行することにより実現される。このコンピュータプログラムは、コンピュータが読み取り可能な記録媒体に格納することができる。入力装置２は、例えば、キーボード、ポインティングデバイスなどでよい。出力装置３は、例えば、ディスプレイ装置、プリンタなどでよい。入出力装置４は、例えば、大容量記憶装置あるいはネットワークインタフェース装置などで良い。
【００１２】
情報処理装置本体１０は、ファントム画像データ記憶部１１と、デジタルファントムデータ記憶部１３と、Ｈ／Ｍ比データ記憶部１５と、変換関数記憶部１７と、被験者画像データ記憶部１９と、位置合わせ処理部２１と、ＲＯＩ設定部２３と、Ｈ／Ｍ比算出部２５と、変換関数算出部２７と、Ｈ／Ｍ比変換処理部２９とを備える。
【００１３】
ファントム画像データ記憶部１１は、異なる撮像環境で所定の撮像条件に調整されたファントムを撮像したファントム画像のデータであるファントム画像データを記憶する。ファントム画像データは、例えば、入力装置２または入出力装置４を介して取得する。ファントム画像データ記憶部１１には、例えば、第１の撮像環境、第２の撮像環境などの撮像環境別、及び撮像条件別にファントム画像データが記憶されている。撮像環境とは、例えば、撮像装置が設置されている施設、撮像装置の機種、撮像時に用いられるコリメータの種類等により特定される。撮像条件とは、ファントムの種類（後述するファントムの構造の違い）及び撮像方向（後述する平板状のファントムの平面視で表または裏）などにより定まる。
【００１４】
本実施形態では、所定の撮像条件でファントムを撮像したファントム画像が用いられる。すなわち、ファントム画像は、各撮像環境で複数の共通の撮像条件で撮像されたファントムの画像である。本実施形態では、例えば、各撮像環境で４つの異なる撮像条件で撮像されたファントム画像が用いられる。
【００１５】
以下、本実施形態で共通に用いられるファントムの構造について説明する。ファントムは、所定の厚みを有し、人体の胸部の断面を模した平板形状を有する。ファントムは、その内部に液体が収容可能な空間を有する。ファントム内部の空間には、被験者に投与する薬剤と同じ薬剤を含む溶液を収容することができる。ファントム内部の空間は、人体（バックグラウンド）領域と、その人体領域内にさらに設けられた心臓領域、甲状腺領域、肝臓領域、肺領域などの各臓器に対応する領域とがある。それぞれの領域は互いに連通している。以下、具体例を用いてさらに詳細に説明する。なお、ここで必要とされる臓器の領域は、用いる薬剤の種類に応じて異なるようにしても良い。
【００１６】
図２は、ファントム５の一例を示す。同図Ａは、ファントム５の平面図であり、同図Ｂは、ファントム５の側面図である。
【００１７】
同図に示すファントム５は、４枚の透明のアクリル板５１〜５４を階層状に重ねた４層構造で構成される。各アクリル板５１〜５４の構成は後述する。４枚のアクリル板５１〜５４は、さらに、同じ材質で構成される透明な板材である外層板５０，５０で挟まれている。外層板５０の厚みは例えば１０ｍｍである。ファントム５は、同図Ａに示すような平面視において、人体領域６１と、心臓領域６２と、肝臓領域６３と、甲状腺領域６４と、肺領域６５（６５ａ，６５ｂ）とを有する。ただし、後述するように、これらの各領域６１〜６５はすべて連通している。以下、この連通している領域を内部空間と呼ぶ。ファントム５は、さらに注入／排出口６０を有する。ファントムの内部空間は、注入／排出口６０を介して外部と連通している。つまり、注入／排出口６０を介して、ファントム５の内部空間に、液体の注入及び排出を行うことができる。注入／排出口６０はネジ式の栓６０ａによって封止され、内部空間が密閉される。
【００１８】
図３は、４層の各アクリル板のそれぞれの平面図である。
【００１９】
同図Ａは第１層のアクリル板５１を示す。第１層のアクリル板５１は、人体領域６１を形成するための開口５１ａを有する。第１層のアクリル板５１の厚みは例えば１５ｍｍである。
【００２０】
同図Ｂは第２層のアクリル板５２を示す。第２層のアクリル板５２は、心臓領域６２、肝臓領域６３、甲状腺領域６４、及び肺領域６５をそれぞれ形成するための開口５２ｂ，５２ｃ，５２ｄ，５２ｅを有する。第２層のアクリル板５２の厚みは例えば５ｍｍである。
【００２１】
同図Ｃは第３層のアクリル板５３を示す。第３層のアクリル板５３は、心臓領域６２及び肝臓領域６３を形成するための開口５３ｂ，５３ｃを有する。第３層のアクリル板５３の厚みは例えば５ｍｍである。
【００２２】
同図Ｄは第４層のアクリル板５４を示す。第４層のアクリル板５４は、肝臓領域６３を形成するための開口５４ｃを有する。第４層のアクリル板５４の厚みは例えば５ｍｍである。
【００２３】
アクリル板５１〜５４を重ねると、平面視において同じ場所に位置する各開口が厚み方向でつながり、内部空間が形成される。その結果、人体領域６１、心臓領域６２、肝臓領域６３、甲状腺領域６４、及び肺領域６５は、それぞれ厚み方向の幅が異なるので、それに応じて各領域６１〜６５に存在する溶液の厚み方向の深さが異なる。例えば、心臓領域６２は、開口５１ａ，５２ｂ，５３ｂで形成され、甲状腺領域６４は、開口５１ａ、５２ｄで形成される。その結果、均一な薬剤濃度の溶液で内部空間が充填されると、各領域６１〜６５はそれぞれの厚みに応じて単位面積あたり異なる放射能量が蓄積されることになる。
【００２４】
第１〜第４層のアクリル板５１〜５４及び外層板５０，５０の厚みは、上述したものに限られない。従って、これと異なる厚みを有する第１〜第４層のアクリル板５１〜５４及び外層板５０，５０でファントムを構成しても良い。本実施形態では、第１〜第４層の厚みの異なる複数のタイプのファントムを用いて撮像したファントム画像を用いる。すなわち、ファントム画像データ記憶部１１には、少なくとも、各撮像環境で、２種類のファントムを表及び裏から撮像して得られる、撮像条件が異なる４枚のファントム画像のデータが記憶される。
【００２５】
図１に戻ると、デジタルファントムデータ記憶部１３は、デジタルファントムデータを記憶する。デジタルファントムデータは、例えば、入力装置２または入出力装置４を介して取得する。デジタルファントムデータは、上述した実物のファントム５の平面視における形状を示すデータである。すなわち、デジタルファントムデータは、図２Ａに示すファントム５内の各領域６１〜６５の相対的な位置及び形状を示す。デジタルファントムデータには、さらに、所定領域にＲＯＩ（Ｒｅａｇｉｏｎ  ｏｆ  Ｉｎｔｅｒｅｓｔ）を設定するためのデータを有する。例えば、本実施形態であれば、デジタルファントムデータには、心臓ＲＯＩ及び縦隔ＲＯＩを設定するための領域を示すデータが含まれる。
【００２６】
Ｈ／Ｍ比データ記憶部１５は、Ｈ／Ｍ比算出部２５によって算出されたＨ／Ｍ比のデータを記憶する。Ｈ／Ｍ比データ記憶部１５には、入力装置２または入出力装置４から取得したＨ／Ｍ比のデータを記憶しても良い。
【００２７】
変換関数記憶部１７は、変換関数算出部２７によって算出された変換関数を示すデータが記憶される。例えば、変換関数記憶部１７には、変換関数と、変換前後の撮像環境の識別情報とが対応づけて記憶される。撮像環境の識別情報としては、例えば、施設名、撮像装置の機種及びコリメータの種類等の組み合わせでも良い。ある撮像環境から、後述する理論上のＨ／Ｍ比への変換関数の場合、変換前の撮像環境と理論値への変換であることを示す情報が対応づけられている。変換関数記憶部１７には、入力装置２または入出力装置４から取得した変換関数のデータを記憶しても良い。
【００２８】
被験者画像データ記憶部１９は、ＭＩＢＧが投与された被験者を撮像したプラナ画像である被験者画像の画像データを記憶する。被験者画像データは、例えば、入力装置２または入出力装置４を介して取得する。被験者画像の各画素値はＲＩのカウント値を示す。各被験者画像は、それぞれ撮像環境を示す情報と対応づけられている。
【００２９】
位置合わせ処理部２１は、ファントム画像データ記憶部１１に保存されているファントム画像データ及びデジタルファントムデータ記憶部１３に保存されているデジタルファントムデータを取得し、ファントム画像に対してデジタルファントムデータの位置を合わせる位置合わせを行う。この位置合わせ処理は、例えば公知の位置合わせアルゴリズムを用いることができる。ここで、デジタルファントムデータは、ファントム画像の心臓領域及び縦隔領域に対応する心臓ＲＯＩ及び縦隔ＲＯＩを有する。従って、ファントム画像データにデジタルファントムが位置合わせされると、自動的に心臓ＲＯＩ及び縦隔ＲＯＩも設定される。すなわち、位置合わせ処理部２１は、デジタルファントムを位置合わせすることによりＲＯＩを設定する。
【００３０】
これにより、ファントム画像に対するＲＯＩ設定が標準化される。そのため、ＲＯＩ設定に関しては術者間のバラツキがなくなり、Ｈ／Ｍ比の精度が向上する。
【００３１】
位置合わせ処理部２１は、被験者画像データ記憶部１９に保存されている被験者画像に対しても同様に、デジタルファントムの位置合わせを行ってもよい。
【００３２】
ＲＯＩ設定部２３は、被験者画像データ記憶部１９から被験者画像データを取得し、被験者画像の所定の領域にＲＯＩを設定する。例えば、ＲＯＩ設定部２３は、被験者画像の心臓領域及び縦隔領域に、それぞれ心臓ＲＯＩ及び縦隔ＲＯＩを設定する。このＲＯＩ設定は、例えば、特許文献１に記載の手法を用いても良い。なお、オペレターが手作業でＲＯＩを設定する場合には、ＲＯＩ設定部２３による処理を省略することができる。
【００３３】
位置合わせ処理部２１で被験者画像に対してデジタルファントムの位置合わせを行ったときは、ＲＯＩ設定部２３によるＲＯＩ設定を省略しても良い。
【００３４】
Ｈ／Ｍ比算出部２５は、ファントム画像データまたは被験者画像データに基づいて、Ｈ／Ｍ比を算出する。Ｈ／Ｍ比算出部２５は、例えば、ファントム画像または被験者画像に設定された心臓ＲＯＩ及び縦隔ＲＯＩ内の画素値（カウント値）に基づいて、Ｈ／Ｍ比を算出する。例えば、Ｈ／Ｍ比算出部２５は、心臓ＲＯＩのカウント値の平均値、及び縦隔ＲＯＩ内のカウント値の平均値を算出し、縦隔ＲＯＩ内のカウント値の平均値を分母として両者の比を算出する。Ｈ／Ｍ算出部２５は、対象画像がファントム画像及び被験者画像のいずれであっても、対象画像に基づいてＨ／Ｍ比を算出する。Ｈ／Ｍ比算出部２５が算出したＨ／Ｍ比は、出力装置３に表示させるようにしてもよい。
【００３５】
変換関数算出部２７は、ある撮像環境から求まるＨ／Ｍ比を、別の撮像環境から求まるＨ／Ｍ比または理論上のＨ／Ｍ比へ変換する変換関数を算出する。算出された変換関数のデータは、変換関数記憶部１７に保存される。
【００３６】
例えば、図４に示すように、変換関数算出部２７は、ある撮像環境（第１の撮像環境）で撮像した画像から算出したＨ／Ｍ比を、別の撮像環境（第２の撮像環境）で撮像した画像から算出されるべきＨ／Ｍ比へ変換する変換関数を算出する。変換関数算出部２７は、Ｈ／Ｍ比データ記憶部１５から第１の撮像環境で複数の撮像条件で撮像したファントム画像から算出した複数のＨ／Ｍ比、及び第２の撮像環境で上記の撮像条件と共通の撮像条件で撮像したファントム画像から算出された複数のＨ／Ｍ比のデータを読み出す。そして、変換関数算出部２７は、共通の撮像条件のＨ／Ｍ比を対応づけて、図４に示すグラフ上にプロットし、回帰分析を行って両者を変換するための関数を導出する。
【００３７】
この変換関数によれば、例えば、第１の撮像環境を有する施設と、第２の撮像環境を有する施設など、異なる施設で撮像された画像から算出されたＨ／Ｍ比が比較可能となる。あるいは、同一施設であっても、撮像装置を交換した場合など撮像環境を変更したときに、変更前の撮像環境で撮像された画像から算出したＨ／Ｍ比と、変更後の撮像装置で撮像された画像から算出したＨ／Ｍ比とを比較可能にするための変換関数を、変換関数算出部２７を用いて算出できる。
【００３８】
また、変換関数算出部２７は、ある撮像環境で撮像した画像から算出したＨ／Ｍ比を、撮像環境に依存しない理論値に標準化する関数を算出する。変換関数算出部２７は、Ｈ／Ｍ比データ記憶部１５から、ある撮像環境で複数の撮像条件で撮像したファントム画像から算出したＨ／Ｍ比を読み出す。このとき、撮像条件として、ファントムの構成（例えば材質、厚さ等）、及び内部空間に収容されている薬剤の種類、濃度等が予め定められた状態に調整されている。この撮像条件に基づいて、ファントムによる放射線の散乱及び吸収なども考慮した上で、そのファントム画像から求められるＨ／Ｍ比が理論的に定まる。そこで、変換関数算出部２７は、撮像されたファントム画像から求めたＨ／Ｍ比と理論上のＨ／Ｍ比とに基づいて、回帰分析を行って両者を変換するための変換関数を導出する。
【００３９】
上述のように、デジタルファントムを用いてファントム画像にＲＯＩを自動設定することによりＨ／Ｍ比の精度が向上している。本実施形態では、このＨ／Ｍ比を用いて変換関数を算出しているので、変換関数も精度よく算出することができる。
【００４０】
Ｈ／Ｍ比変換処理部２９は、変換関数記憶部１７から所望の変換関数のデータを取得し、Ｈ／Ｍデータ記憶部１５に保存されているＨ／Ｍ比の変換を行う。例えば、Ｈ／Ｍ比変換処理部２９は、第１の撮像環境で撮像した被験者画像から算出したＨ／Ｍ比を、第２の撮像環境で撮像した被験者画像から算出されるべきＨ／Ｍ比へ変換するためには、それに対応する変換関数のデータを変換関数記憶部１７から取得し、その変換関数を用いてＨ／Ｍデータ記憶部１５に保存されているＨ／Ｍ比の変換処理を行う。あるいは、Ｈ／Ｍ比変換処理部２９は、同様にして、ある撮像環境で撮像した被験者画像から算出したＨ／Ｍ比を、理論上のＨ／Ｍ比へ変換する。
【００４１】
上述のように、変換関数の精度が高いので、被験者画像から求めたＨ／Ｍ比の変換も高精度で行うことができる。
【００４２】
上記のような構成を備える本実施形態の情報処理装置は、以下に説明する手順で変換関数の算出及びＨ／Ｍ比の変換処理を行う。
【００４３】
図５は、変換関数を算出する手順を示すフローチャートである。ここでは、第１の撮像環境のＨ／Ｍ比を第２の撮像環境のＨ／Ｍ比へ変換する関数、及び第１の撮像環境のＨ／Ｍ比を理論上のＨ／Ｍ比へ変換する関数を算出する処理について説明する。
【００４４】
第１及び第２の撮像環境において、所定の撮像条件で撮像された複数のファントム画像の画像データをファントム画像データ記憶部１１に保存する（Ｓ１１）。この画像データは、各ファントム画像を撮像した撮像環境及び撮像条件を示す識別情報と対応づけて保存される。デジタルファントムデータは、デジタルファントムデータ記憶部１３に保存される（Ｓ１３）。画像データ及びデジタルファントムデータは、例えば、入出力装置４を介して外部から取得しても良い。ファントム画像データとデジタルファントムデータとが準備されると、これ以降の処理が実行可能となる。ステップＳ１１とＳ１３はいずれを先に行っても良いし、既にこれらのデータが準備されているときには、それぞれ対応するステップを省略することができる。また、本実施形態では、複数のファントム画像のデータをファントム画像データ記憶部１１に蓄積してからステップＳ１５以降を行っているが、一つのファントム画像の画像データをファントム画像データ記憶部１１に蓄積するごとにステップＳ１５以降を行うようにしても良い。
【００４５】
位置合わせ処理部２１は、一つのファントム画像を処理対象とし、ファントム画像データ記憶部１１から処理対象のファントム画像の画像データを読み出す。位置合わせ処理部２１はさらに、デジタルファントムデータ記憶部１３からデジタルファントムデータを読み出す。そして、位置合わせ処理部２１は、処理対象のファントム画像に対してデジタルファントムの位置合わせを行う。これにより、同時にファントム画像の心臓領域及び縦隔領域にＲＯＩが設定される（Ｓ１５）。
【００４６】
Ｈ／Ｍ比算出部２５が、ステップＳ１５でファントム画像に設定されたＲＯＩ内の画像のカウント値に基づいて、Ｈ／Ｍ比を算出する（Ｓ１９）。ここで算出されたＨ／Ｍ比のデータはＨ／Ｍ比データ記憶部１５に保存される。
【００４７】
ここで、処理対象のファントム画像のすべてについてＨ／Ｍ比の算出が終わるまで、ステップＳ１５へ戻り、Ｈ／Ｍ比の算出処理を繰り返す（Ｓ２１：Ｎｏ）。
【００４８】
処理対象のファントム画像のすべてのＨ／Ｍ比の算出が終わると（Ｓ２１：Ｙｅｓ）、変換関数算出部２７がＨ／Ｍ比データ記憶部１５から第１及び第２の撮像環境のファントム画像から算出したＨ／Ｍ比のデータを読み出し、２つの環境間の変換関数を算出する（Ｓ２３）。また、第１の撮像環境のＨ／Ｍ比を理論上のＨ／Ｍ比へ変換する変換関数を算出する場合は、変換関数算出部２７が予め有している、撮像条件毎の理論上のＨ／Ｍ比を用いて変換関数を算出する。変換関数算出部２７が算出した変換関数を示すデータが変換関数記憶部１７に保存される。
【００４９】
ここで算出された変換関数に関する情報は、出力装置３によって出力されても良い。
【００５０】
次に、図６は、上記の変換関数を用いて、第１の撮像環境で被験者を撮像した画像から求めたＨ／Ｍ比を変換する手順を示すフローチャートである。
【００５１】
第１の撮像環境において、ＭＩＢＧが投与された被験者を撮像した被験者画像の画像データを被験者画像データ記憶部１９に保存する（Ｓ３１）。画像データは、例えば、入出力装置４を介して外部から取得しても良い。被験者画像データが準備されると、これ以降の処理が実行可能となる。既に被験者画像データが準備されているときには、ステップＳ３１を省略することができる。
【００５２】
ＲＯＩ設定部２３は、被験者画像データ記憶部１９から処理対象の被験者画像の画像データを読み出す。ＲＯＩ設定部２３は、さらに、対象の被験者画像に対して、心臓領域及び縦隔領域にＲＯＩを設定する（Ｓ３７）。
【００５３】
Ｈ／Ｍ比算出部２５が、ステップＳ３７で設定されたＲＯＩ内の画像のカウント値に基づいて、Ｈ／Ｍ比を算出する（Ｓ３９）。ここで算出されたＨ／Ｍ比のデータはＨ／Ｍ比データ記憶部１５に保存される。
【００５４】
Ｈ／Ｍ比変換処理部２９は、Ｈ／Ｍ比データ記憶部１５からステップＳ３９で算出されたＨ／Ｍ比のデータを読み出し、変換関数記憶部１７から所望の変換関数のデータを読み出す。そして、Ｈ／Ｍ比変換処理部２９は、その変換関数を用いて、Ｈ／Ｍ比を変換する（Ｓ４１）。例えば、Ｈ／Ｍ比変換処理部２９は、変換関数記憶部１７から第１の撮像環境から第２の撮像環境へ変換するための変換関数のデータを読み出し、第１の撮像環境のＨ／Ｍ比を第２の撮像環境のＨ／Ｍ比へ変換する。あるいは、Ｈ／Ｍ比変換処理部２９は、変換関数記憶部１７から第１の撮像環境のＨ／Ｍ比を理論上のＨ／Ｍ比へ変換するための変換関数のデータを読み出し、第１の撮像環境のＨ／Ｍ比を理論値へ変換する。
【００５５】
この変換前後のＨ／Ｍ比は、出力装置３によって出力されても良い。
【００５６】
本実施形態によれば、ＭＩＢＧなどの心機能診断薬が投与された被験者を異なる撮像環境で撮像した画像から算出したＨ／Ｍ比を互いに精度良く比較することが可能となる。特に、Ｈ／Ｍ比を理論値へ変換することにより、撮像環境に依存しないＨ／Ｍ比に精度良く標準化することができる。
【００５７】
以下、本発明の第２の実施形態に係る情報処理装置について、図面を参照して説明する。なお、第１の実施形態と同一の構成については同一の参照番号を付して説明を省略する。
【００５８】
本実施形態に係る情報処理装置は、所定のイベントが生じる可能性を示す指標を算出するための関数を導出する処理を行う。ここで、本実施形態における所定のイベントとは、例えば、心不全患者の所定期間後（例えば、５年後）の生死である。指標は、所定のイベントが将来起こるか否かを示す指標であり、本実施の形態では患者が５年後に生存している可能性を示す指標（５年後生存率）である。
【００５９】
図７は、本実施形態に係る情報処理装置１０１の全体構成を示す。同図に示すように、情報処理装置１０１は、情報処理装置本体１１０に入力装置２、出力装置３及び入出力装置４の一つ以上が接続されている。情報処理装置本体１１０は、例えば、プロセッサ及びメモリを備えた汎用的なコンピュータシステムにより構成され、以下に説明する情報処理装置本体１１０内の個々の構成要素または機能は、例えば、コンピュータプログラムを実行することにより実現される。このコンピュータプログラムは、コンピュータが読み取り可能な記録媒体に格納することができる。
【００６０】
情報処理装置本体１１０は、分析用データ記憶部１１１と、被験者データ記憶部１１３と、生存率関数記憶部１１５と、ロジスティック回帰分析実行部１１７と、生存率算出部１１９とを備える。
【００６１】
分析用データ記憶部１１１は、単一または複数の施設において多数の心不全の被験者（第２の被験者に相当する。）に対し、定期的又は不定期に行った検査に関する分析用データを記憶している。分析用データは、年齢と、性別と、Ｈ／Ｍ比（Ｈｅａｒｔ／Ｍｅｄｉａｓｔｉｎｕｍ比）と、心不全の重症度を示す指標（例えば、ＮＹＨＡ  ｃｌａｓｓ（Ｎｅｗ  Ｙｏｒｋ  Ｈｅａｒｔ  Ａｓｓｏｃｉａｔｉｏｎ））と、左室機能を示す指標（例えば、ＬＶＥＦ（Ｌｅｆｔ  Ｖｅｎｔｒｉｃｕｌａｒ  Ｅｊｅｃｔｉｏｎ  Ｆｒａｃｔｉｏｎ））と、５年後の生死とから構成される。分析用データは、例えば、入力装置２または入出力装置４を介して取得する。なお、年齢、性別、Ｈ／Ｍ比、心不全の重症度を示す指標、及び左室機能を示す指標は、第２の被験者の被験者データに相当し、５年後の生死は、イベントデータに相当する。
【００６２】
年齢は、各被験者の検査時における年齢である。性別は、各被験者の性別であり、男もしくは女である。
【００６３】
Ｈ／Ｍ比は、第１の実施形態に記載のように、心機能診断薬としてＭＩＢＧ（１２３Ｉ−３−ｉｏｄｏｂｅｎｚｙｌｇｕａｎｉｄｉｎｅ）を被験者に投与した後、その被験者から放出された放射能を検出して得られるＲＩプラナ画像から取得した心／縦隔比で良い。なお、分析用データ記憶部１１１に記憶されるＨ／Ｍ比は、第１の実施形態におけるＨ／Ｍ比算出部２５によって算出されたＨ／Ｍ比のデータであっても良いし、Ｈ／Ｍ比変換処理部２９により変換された撮像環境に依存しない標準化されたＨ／Ｍ比であっても良い。更に、Ｈ／Ｍ比変換処理部２９により変換された第２の撮像環境のＨ／Ｍ比であっても良い。
【００６４】
ＮＹＨＡ  ｃｌａｓｓは、ニューヨーク心臓協会（Ｎｅｗ  Ｙｏｒｋ  Ｈｅａｒｔ  Ａｓｓｏｃｉａｔｉｏｎ：ＮＹＨＡ）が定めた心不全の症状の程度の分類であり、画像及び問診等を総合的に判断することにより、以下のように心不全の重症度を４種類に分類するものである。
ＮＹＨＡ  Ｉ度：心疾患があるが症状はなく、通常の日常生活は制限されないもの。
ＮＹＨＡ  ＩＩ度：心疾患患者で日常生活が軽度から中等度に制限されるもの。安静時には無症状だが、普通の行動で疲労・動悸・呼吸困難・狭心痛を生じる。
ＮＹＨＡ  ＩＩＩ度：心疾患患者で日常生活が高度に制限されるもの。安静時は無症状だが、平地の歩行や日常生活以下の労作によっても症状が生じる。
ＮＹＨＡ  ＩＶ度：心疾患患者で非常に軽度の活動でも何らかの症状を生ずる。安静時においても心不全・狭心症症状を生ずることもある。
【００６５】
ＬＶＥＦは、左室駆出率であり、左心室の収縮力を表す指標である。例えば、ＬＶＥＦは、超音波エコー画像の解析により取得され、一回の駆出によって左心室の心室容積の何％が駆出されるかを表す。
【００６６】
５年後の生死は、心機能検査を行った被験者が、検査から５年後に生存しているか否か（生存又は死亡）を示すデータである。即ち、５年後の生死のデータは、所定のイベントの成否の実績のデータである。
【００６７】
被験者データ記憶部１１３は、後述の５年生存率算出の対象となる被験者（第１の被験者に相当する。）の被験者データを記憶する。被験者データは、例えば、被験者の年齢と、性別と、Ｈ／Ｍ比と、ＮＹＨＡ  ｃｌａｓｓと、ＬＶＥＦとを含む。また、被験者データは、一人の被験者の被験者データであっても良いし、複数の被験者の被験者データであっても良い。被験者データは、例えば、入力装置２または入出力装置４を介して取得する。
【００６８】
なお、被験者データ記憶部１１３に記憶されるＨ／Ｍ比は、第１の実施形態におけるＨ／Ｍ比算出部２５によって算出されたＨ／Ｍ比のデータであっても良いし、Ｈ／Ｍ比変換処理部２９により変換された撮像環境に依存しない標準化されたＨ／Ｍ比であっても良い。更に、Ｈ／Ｍ比変換処理部２９により変換された第２の撮像環境のＨ／Ｍ比であっても良い。また、被験者データ記憶部１１３及び分析用データ記憶部１１１に記憶されるＨ／Ｍ比は、上記のどのＨ／Ｍ比であっても良いが、撮像環境に依存しない標準化されたＨ／Ｍ比又は第２の撮像環境のＨ／Ｍ比を両方の記憶部に記憶して、後述の５年生存率関数の導出及び５年生存率算出に用いることにより、５年生存率に関してさらに精度の高い関数を導出することができ、さらに精度の高い５年生存率を算出することができる。
【００６９】
生存率関数記憶部１１５は、ロジスティック回帰分析実行部１１７により算出された５年生存率関数を記憶する。
【００７０】
ロジスティック回帰分析実行部１１７は、分析用データ記憶部１１１に記憶されている分析用データを読み出し、公知のロジスティック回帰分析を実行して、５年生存率に関する関数（５年生存率＝Ｆ（年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、ＬＶＥＦ））を導出する。よって、５年生存率に関する関数は、所定のイベントの成否（５年後の生死）の実績に応じて定められた関数であり、将来イベントが生じる可能性を示す指標を算出するための関数である。即ち、分析用データの内、年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、及びＬＶＥＦはロジスティック回帰分析の独立変数（説明変数）であり、５年生存率はロジスティック回帰分析の従属変数（目的変数）である。
【００７１】
生存率算出部１１９は、被験者データ記憶部１１３に記憶されている被験者データと、生存率関数記憶部１１５に記憶されている５年生存率関数とを読み出して、５年生存率を算出する。具体的には、５年生存率関数に、被験者データの年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、及びＬＶＥＦを代入することにより、５年生存率を算出する。５年生存率は、何％のような表現であっても良いし、生存又は死亡のような表現であっても良い。また、算出された５年生存率は、出力装置３に表示されても良いし、出力装置３により印刷され出力されても良い。更に、算出された５年生存率は、図示せぬ記憶部に記憶されても良い。
【００７２】
上記のような構成を備える本実施形態の情報処理装置１０１は、以下に説明する手順で５年生存率関数の導出及び５年生存率の算出（所定のイベントが生じる可能性を示す指標の算出）を行う。
【００７３】
図８は、５年生存率関数を導出する手順を示すフローチャートである。まず、多数の心不全の被験者に対し行った検査に関する分析用データを分析用データ記憶部１１１に保存する（Ｓ４１）。分析用データは、入出力装置４を介して外部から取得しても良い。既に分析用データが準備されているときには、当該ステップを省略しても良い。
【００７４】
次に、ロジスティック回帰分析実行部１１７が、分析用データ記憶部１１１に記憶された分析用データを読み出し、ロジスティック回帰分析を実行し、５年生存率に関する関数（５年生存率＝Ｆ（年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、ＬＶＥＦ））を導出する（Ｓ４３）。そして、導出された５年生存率に関する関数は、５年生存率関数記憶部１１５に保存される（Ｓ４５）。
【００７５】
図９は、５年生存率を算出する手順を示すフローチャートである。まず、被験者データを被験者データ記憶部１１３に保存する（Ｓ５１）。被験者データは、入出力装置４を介して外部から取得しても良い。既に被験者データが準備されているときには、当該ステップを省略しても良い。
【００７６】
次に、生存率算出部１１９は、被験者データ記憶部１１３に保存された被験者データを取得し（Ｓ５３）、生存率関数記憶部１１５に保存された５年生存率関数を読み出す（Ｓ５５）。そして、読み出した５年生存率関数に被験者データ（年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、及びＬＶＥＦ）を代入することにより、５年生存率を算出する（Ｓ５７）。算出された５年生存率は、出力装置３に表示されても良いし、出力装置３により印刷され出力されても良い。
【００７７】
本実施の形態によれば、独立変数の候補として挙げられる多数の変数の内、年齢、性別、Ｈ／Ｍ比、ＮＹＨＡ  ｃｌａｓｓ、及びＬＶＥＦを独立変数として選択し、ロジスティック回帰分析を行っているので、５年生存率に関して精度の高い関数を導出することができる。よって、精度の高い５年生存率を算出することができる。また、Ｈ／Ｍ比に、撮像環境に依存しない標準化されたＨ／Ｍ比、又は第２の撮像環境のＨ／Ｍ比を用いることにより、５年生存率に関してさらに精度の高い関数を導出することができ、さらに精度の高い５年生存率を算出することができる。
【００７８】
上述した本発明の実施形態は、本発明の説明のための例示であり、本発明の範囲をそれらの実施形態にのみ限定する趣旨ではない。当業者は、本発明の要旨を逸脱することなしに、他の様々な態様で本発明を実施することができる。
【００７９】
また、上記の第２の実施形態では、所定のイベントが生じる可能性を示す指標として、５年生存率に関する指標を算出したが、生存率は何年であっても良いし、生存率ではなく致死性不整脈、又は、心臓突然死であっても良い。
【００８０】
また、独立変数としてのＮＹＨＡ  ｃｌａｓｓは、Ｉ度〜ＩＶ度の４つに分類していたが、Ｉ度及びＩＩ度を１カテゴリー、ＩＩＩ度及びＩＶ度を１カテゴリーとして、ＮＹＨＡ  ｃｌａｓｓを２つに分類して回帰分析を行っても良い。更に、年齢を１０歳毎に１カテゴリーとして回帰分析を行っても良い。
【００８１】
また、上記の実施形態では、Ｈ／Ｍ比を算出する処理と、５年生存率を算出する処理を分けて行ったが、５年生存率関数を予め算出し、Ｈ／Ｍ比を算出する処理に５年生存関数を組み込んで、Ｈ／Ｍ比が算出されると５年生存率が算出されるようにしても良い。
【００８２】
また、被験者データには、Ｈ／Ｍ比と、被験者の年齢、性別、心不全の重症度を示す指標、及び左室機能を示す指標の少なくともいずれか一つと、が含まれていれば良い。また、上記の実施形態では、５年生存率に関する関数の独立変数は、Ｈ／Ｍ比、被験者の年齢、性別、心不全の重症度を示す指標、及び左室機能を示す指標の５つであったが、独立変数は、Ｈ／Ｍ比の１つでも良く、Ｈ／Ｍ比に加え、被験者の年齢、性別、心不全の重症度を示す指標、及び左室機能を示す指標、のうちの少なくとも一つ以上であれば良い。
【符号の説明】
【００８３】
１、１０１  情報処理装置、５  ファントム、１０、１１０  情報処理装置本体、１１  ファントム画像データ記憶部、１３  デジタルファントムデータ記憶部、１５  Ｈ／Ｍ比データ記憶部、１７  変換関数記憶部、１９  被験者画像データ記憶部、２１  位置合わせ処理部、２３  ＲＯＩ設定部、２５  Ｈ／Ｍ比算出部、２７  変換関数算出部、２９  Ｈ／Ｍ比変換処理部、１１１  分析用データ記憶部、  １１３  被験者データ記憶部、  １１５  生存率関数記憶部、  １１７  ロジスティック回帰分析実行部、  １１９  生存率算出部
【図１】
2014215222.tif 000003
【図２】
2014215222.tif 000004
【図３】
2014215222.tif 000005
【図４】
2014215222.tif 000006
【図５】
2014215222.tif 000007
【図６】
2014215222.tif 000008
【図７】
2014215222.tif 000009
【図８】
2014215222.tif 000010
【図９】
2014215222.tif 000011