特許 7034218 シミュレータ、該シミュレータを備える注入装置又は撮像システム、及びシミュレーションプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-03-03

(45)【発行日】2022-03-11

(54)【発明の名称】シミュレータ、該シミュレータを備える注入装置又は撮像システム、及びシミュレーションプログラム

(51)【国際特許分類】

   A61B 6/00 20060101AFI20220304BHJP

   A61B 6/03 20060101ALI20220304BHJP

   A61B 5/055 20060101ALI20220304BHJP

   A61M 5/142 20060101ALI20220304BHJP

【ＦＩ】

A61B6/00 331E

A61B6/03 375

A61B6/03 360D

A61B5/055 383

A61M5/142 530

【請求項の数】 11

(21)【出願番号】P 2020125181

(22)【出願日】2020-07-22

(62)【分割の表示】P 2016561247の分割

【原出願日】2015-11-26

(65)【公開番号】P2020179229

(43)【公開日】2020-11-05

【審査請求日】2020-08-17

(31)【優先権主張番号】P 2014240006

(32)【優先日】2014-11-27

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】504136568

【氏名又は名称】国立大学法人広島大学

(73)【特許権者】

【識別番号】391039313

【氏名又は名称】株式会社根本杏林堂

(74)【代理人】

【識別番号】100094112

【弁理士】

【氏名又は名称】岡部 讓

(74)【代理人】

【識別番号】100101498

【弁理士】

【氏名又は名称】越智 隆夫

(74)【代理人】

【識別番号】100107401

【弁理士】

【氏名又は名称】高橋 誠一郎

(74)【代理人】

【識別番号】100120064

【弁理士】

【氏名又は名称】松井 孝夫

(72)【発明者】

【氏名】檜垣 徹

(72)【発明者】

【氏名】粟井 和夫

(72)【発明者】

【氏名】中村 優子

(72)【発明者】

【氏名】増田 和正

(72)【発明者】

【氏名】弓場 孝治

(72)【発明者】

【氏名】根本 茂

【審査官】亀澤 智博

(56)【参考文献】

【文献】特表２０００－５０６３９８（ＪＰ，Ａ）

【文献】特開２０１３－１０６９９０（ＪＰ，Ａ）

【文献】特開２０１２－１４７９３４（ＪＰ，Ａ）

【文献】特表平０７－５０８２０２（ＪＰ，Ａ）

【文献】特表２０１３－５３３０２８（ＪＰ，Ａ）

【文献】特開２０１４－１００５５５（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１３／００２８４９４（ＵＳ，Ａ１）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ６／００ － ６／１４

Ａ６１Ｂ ５／０５５

Ａ６１Ｍ ５／００ － ５／５２

Ｇ０６Ｔ １／００

(57)【特許請求の範囲】

【請求項1】

被写体の組織における画素値の経時変化を予測するシミュレータであって、
前記被写体の情報を取得する被写体情報取得部と、
造影剤の注入プロトコルを取得するプロトコル取得部と、
前記組織の情報を取得する組織情報取得部と、
前記被写体の情報と、前記注入プロトコルと、前記組織の情報とに基づいて、前記組織の画素値の経時変化を予測する予測部とを備え、
前記組織の情報は、前記組織における第１のパラメータである前記造影剤の染み出し速度と、前記第１のパラメータとは異なる第２のパラメータである染み戻り速度と、を含む、シミュレータ。

【請求項2】

前記染み出し速度及び前記染み戻り速度は、各組織の質量に応じた毛細血管面積を各組織に割り当て、各組織の毛細血管面積と各組織の透過性の積によって算出される、請求項１に記載のシミュレータ。

【請求項3】

表示部と、前記表示部を制御する表示制御部とを備え、
前記表示制御部は、所定の組織の予測画像を前記画素値に応じた濃度の色で表示するように、前記表示部を制御する、請求項１又は２に記載のシミュレータ。

【請求項4】

前記表示制御部は、前記画素値の経時変化に応じて前記予測画像の濃淡を変更するように、前記表示部を制御する、請求項３に記載のシミュレータ。

【請求項5】

前記表示制御部は、オペレータが選択した組織の時間濃度曲線を、前記予測画像と同時に表示するように、前記表示部を制御する、請求項３又は４に記載のシミュレータ。

【請求項6】

前記表示制御部は、所定の時間に遅延時間を加算した時間における画素値に応じた濃度の色で前記予測画像を表示するように、前記表示部を制御する、請求項３から５のいずれか１項に記載のシミュレータ。

【請求項7】

被写体の組織における画素値の経時変化を予測するシミュレータであって、
前記被写体の情報を取得する被写体情報取得部と、
造影剤の注入プロトコルを取得するプロトコル取得部と、
前記組織の情報を取得する組織情報取得部と、
前記被写体の情報と、前記注入プロトコルと、前記組織の情報とに基づいて、前記組織の画素値の経時変化を予測する予測部と、
表示部と、
前記表示部を制御して所定の組織の予測画像を画素値に応じた表示様式で表示させる表示制御部とを備え、
前記組織の情報は、前記組織における前記造影剤の染み出し速度と染み戻り速度とを含み、
前記表示制御部は、オペレータが選択した組織の時間濃度曲線を、前記予測画像と同時に前記表示部に表示させる、シミュレータ。

【請求項8】

被写体の組織における画素値の経時変化を予測するシミュレータであって、
前記被写体の情報を取得する被写体情報取得部と、
造影剤の注入プロトコルを取得するプロトコル取得部と、
前記組織の情報を取得する組織情報取得部と、
前記被写体の情報と、前記注入プロトコルと、前記組織の情報とに基づき、前記造影剤の注入開始からの経過時間として規定する所定の時間に、ヘリカルスキャンによる遅延時間を加算した時間における前記組織の画素値の経時変化を予測する予測部とを備え、
前記組織の情報は、前記組織における前記造影剤の染み出し速度と染み戻り速度とを含む、シミュレータ。

【請求項9】

表示部と、前記表示部を制御して所定の組織の予測画像を画素値に応じた表示様式で表示させる表示制御部とを更に備え、
前記表示制御部は、前記造影剤の注入開始からの経過時間として規定する所定の時間に、ヘリカルスキャンによる遅延時間を加算した時間における画素値に応じた表示様式で前記予測画像を前記表示部に表示させる、請求項８に記載のシミュレータ。

【請求項10】

注入プロトコルに従って造影剤を注入する注入ヘッドと、
請求項１から９のいずれか１項に記載のシミュレータとを備える、注入装置。

【請求項11】

被写体を撮像する医療用の撮像装置と、
請求項１から９のいずれか１項に記載のシミュレータとを備える、撮像システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、被写体の組織における撮像画像の画素値（ＣＴ値）の経時変化を予測するシミュレータ、該シミュレータを備える注入装置又は撮像システム、及びシミュレーションプログラムに関する。

【背景技術】

【0002】

従来、被写体としての患者の組織を、血管内に注入される造影剤で増強してＣＴ（Computed Tomography）装置を用いて撮像する方法があった。また、被写体の体質的特徴と、造影剤の注入プロトコルとにより、造影剤によるコントラスト強度の増強レベル（画素値）を予測する予測方法が知られていた。そして、この予測方法は、被写体の組織におけるコントラスト強度の増強程度を、造影剤の注入開始時刻からの時間経過の関数として予測している（特許文献１）。

【先行技術文献】

【特許文献】

【0003】

【文献】特許第３５５３９６８号明細書

【発明の概要】

【発明が解決しようとする課題】

【0004】

特許文献１に記載の方法では、心臓及び血管が１つのコンパートメントであると仮定されている。また、その他の臓器は、血管内空間及び細胞外空間を有する１つのコンパートメントであると仮定されている。そして、各コンパートメントにおいて、造影剤の到達と同時にコンパートメント全体に造影剤が拡散すると仮定して予測が行われる。しかし、実際には組織毎に体積が異なるため、予測結果が実際の画素値の変化と大きく異なることがあった。特に、造影剤の注入量が少ない場合、造影剤の注入時間が短い場合、又は造影剤の濃度が低い場合には、実際の組織内での拡散速度が遅い。そのため、予測結果は、実際の画素値の変化と大きく異なりやすい。

【課題を解決するための手段】

【0005】

上記課題を解決するため、本発明の一例としてのシミュレータは、被写体の組織における画素値の経時変化を予測するシミュレータであって、前記被写体の情報を取得する被写体情報取得部と、造影剤の注入プロトコルを取得するプロトコル取得部と、前記組織の情報を取得する組織情報取得部と、前記被写体の情報と、前記注入プロトコルと、前記組織の情報とに基づいて、前記組織の画素値の経時変化を予測する予測部とを備え、前記組織の情報は、前記組織の情報における第１のパラメータである前記造影剤の染み出し速度と、前記第１のパラメータとは異なる第２のパラメータである染み戻り速度と、を含む。

【0006】

また、本発明の他の例としての注入装置は、注入プロトコルに従って造影剤を注入する注入ヘッドと、上記シミュレータとを備える。

【0007】

また、本発明の他の例としての撮像システムは、被写体を撮像する医療用の撮像装置と、上記シミュレータとを備える。

【0008】

また、本発明の他の例としてのシミュレーションプログラムは、被写体の組織における画素値の経時変化をコンピューターに予測させるシミュレーションプログラムであって、前記シミュレーションプログラムは前記コンピューターを、前記被写体の情報を取得する被写体情報取得部、造影剤の注入プロトコルを取得するプロトコル取得部、前記組織の情報を取得する組織情報取得部、及び前記被写体の情報と、前記注入プロトコルと、前記組織の情報とに基づいて、前記組織を血流方向に沿って分割した複数のコンパートメントのそれぞれの画素値の経時変化を予測する予測部として機能させる。

【0009】

これにより、実際の組織における画素値の経時変化に近似したより高精度の予測を行うことができる。特に、造影剤の注入量が少ない場合、造影剤の注入時間が短い場合、又は造影剤の濃度が低い場合であっても、高精度の予測を行うことができる。また、各組織内における造影剤の位置を予測することができる。そのため、最適な注入条件又は撮像条件を実際に注入する前に予測することができるので、組織の撮像に失敗することを防止できる。

【0010】

本発明のさらなる特徴は、添付図面を参照して例示的に示した以下の実施例の説明から明らかになる。

【図面の簡単な説明】

【0011】

【図1】シミュレータの概略ブロック図である。

【図2】第１実施形態に係る複数のコンパートメントを説明する図面である。

【図3】血流モデルを説明する図面である。

【図4】胃、脾臓、膵臓及び腸管のパラメータを示す表である。

【図5】シミュレータの表示部に表示される予測画像である。

【図6】表示部に表示される一例としての操作画面である。

【図7】ヘリカルスキャン選択時にシミュレータの表示部に表示される予測画像である。

【図8】注入装置及び撮像システムの概略図である。

【図9】第２実施形態に係る複数のコンパートメントを説明する図面である。

【図10】加算処理を説明するフローチャートである。

【図11】変形形態に関する血流モデルを説明する図面である。

【図12】変形形態に関するタイムデンシティカーブを示すグラフである。

【発明を実施するための形態】

【0012】

以下、本発明を実施するための例示的な実施形態を、図面を参照して詳細に説明する。ただし、以下の実施形態で説明する寸法、材料、形状、構成要素の相対的な位置等は任意であり、本発明が適用される装置の構成又は様々な条件に応じて変更できる。また、特別な記載がない限り、本発明の範囲は、以下に具体的に記載された実施形態に限定されるものではない。なお、本明細書において、上下とは重力方向における上方向と下方向とにそれぞれ対応する。

【0013】

［第１実施形態］
図１に示すように、被写体の組織における画素値の経時変化を予測するシミュレータ（灌流シミュレータ）２０は、予測部１６を備えている。そして、予測部１６は、被写体の情報と、注入プロトコルと、組織の情報とに基づいて、被写体の組織を血流方向に沿って分割した複数のコンパートメントのそれぞれの画素値の少なくとも造影剤に起因する経時変化を予測する。なお、画素値は、造影剤の他に生理食塩水及び血液等の影響を受けて変化する。

【0014】

また、シミュレータ２０は、被写体の情報を取得する被写体情報取得部１１を備えている。そして、予測部１６は、被写体情報取得部１１から、ヘモグロビン量（g/dL）及び被写体の体重（kg）等の被写体である被験者に関する被写体情報を受け取る。ここで被写体情報取得部１１は、シミュレータ２０の入力部２７を介して、オペレーターが入力した被写体情報を取得する。

【0015】

また、被写体情報取得部１１は、シミュレータ２０の記憶部２４又は外部記憶装置（サーバー）から被写体情報を取得してもよい。このようなサーバーとしては、例えば、ＲＩＳ（Radiology Information System）、ＰＡＣＳ（Picture Archiving and Communication Systems）、ＨＩＳ（Hospital Information System）、検像システム及び画像作成用ワークステーション等がある。さらに、被写体情報取得部１１は、撮像装置３（図８）又は注入装置２（図８）から被写体情報を取得してもよい。なお、被写体情報は、除脂肪体重、循環血液量、被験者番号（被験者ＩＤ）、被験者氏名、性別、生年月日、年齢、身長、血液量、血流速度、体表面積、被験者の疾病、副作用の履歴、クレアチニン値、心拍数、及び心拍出量等を含んでいてもよい。

【0016】

また、シミュレータ２０は、造影剤の注入プロトコルを取得するプロトコル取得部１２を備える。そして、予測部１６は、プロトコル取得部１２から、造影剤注入速度（mL/sec）及び造影剤注入時間（sec）等の注入プロトコルを取得する。ここでプロトコル取得部１２は、入力部２７を介してオペレーターが入力した注入プロトコルを取得する。なお、注入プロトコルは、注入方法、造影剤注入箇所、注入量、注入タイミング、造影剤濃度及び注入圧力等の注入条件に関する情報を含んでいてもよい。また、プロトコル取得部１２は、記憶部２４、外部記憶装置又は注入装置２から注入プロトコルを取得してもよい。

【0017】

特に、造影剤注入箇所は、シミュレータ２０の注入設定画面から入力することができ、標準の場合は上肢静脈が選択される。その他の注入箇所としては、肝動脈（肝動脈造影法：ＣＴＨＡ）、上腸間膜動脈（門脈造影法：ＣＴＡＰ）、右心室、上行大動脈等を選択することができる。また、注入プロトコルは、一定の造影剤注入速度の他、造影剤の後押し注入の有無、生理食塩水の注入速度、生理食塩水の注入時間、注入速度の増減、及び注入用チューブの体積等の情報を含んでいてもよい。

【0018】

また、シミュレータ２０は、被写体の組織の情報を取得する組織情報取得部１３を備えている。そして、予測部１６は、組織情報取得部１３から、組織におけるコンパートメント数（血管及び臓器の分割コンパートメント数）、組織の体積（血管腔の体積）、毛細血管の体積、細胞外液腔の体積、単位組織あたりの血流量（血流速度）、組織における造影剤の染み出し速度（毛細血管透過性表面積）、組織における造影剤の染み戻り速度（毛細血管透過性表面積）、及び組織生来の画素値等の組織情報を取得する。

【0019】

ここで組織情報取得部１３は、入力部２７を介してオペレーターが入力した被写体情報を取得する。なお、組織には、心臓（右心室及び左心室）、血管、その他の臓器及び筋肉等が含まれる。そして、予測部１６が組織生来の画素値を取得した場合、各組織の生来の画素値に基づいて造影剤による増強の程度が予測される。また、組織情報取得部１３は、記憶部２４、外部記憶装置又は注入装置２から組織情報を取得してもよい。

【0020】

また、シミュレータ２０は、薬液に関する薬液情報を取得する薬液情報取得部１４を備えている。そして、予測部１６は、薬液情報取得部１４から、造影剤濃度（mgI/mL）、造影剤量（mL）、総ヨード量（mgI）及び造影剤の半減期（造影剤排出レート）等の薬液情報を取得する。さらに、予測部１６は、総ヨード量と被写体の体重から、体重１ｋｇあたりのヨード量（mgI/kg）を算出することができる。ここで薬液情報取得部１４は、入力部２７を介してオペレーターが入力した薬液情報を取得する。なお、薬液情報は、製品名称、製品ＩＤ、化学分類、含有成分、濃度、粘度、消費期限、シリンジ容量、シリンジ耐圧、シリンダ内径、ピストンストローク及びロット番号等を含んでいてもよい。

【0021】

また、薬液情報取得部１４は、記憶部２４、外部記憶装置又は注入装置２から薬液情報を取得してもよい。さらに、薬液情報取得部１４は、注入装置２に内蔵された読取部から薬液情報を取得してもよい。そして、読取部は、注入ヘッドに搭載されるシリンジに取り付けられたデータキャリアの読み取りを行う。このデータキャリアには、ＲＦＩＤチップ、ＩＣタグ、バーコード等があり、薬液に関連する薬液情報が記憶されている。

【0022】

また、予測部１６は、入力部２７を介して管球電圧（kVp）等の検査情報を取得することができる。この検査情報は、検査番号（検査ＩＤ）、検査部位、検査日時、薬液種類、薬液名称、及び撮像条件（撮像部位等）等を含んでいてもよい。

【0023】

さらに、予測部１６は、局所灌流（ボーラス伝送）の考慮の有無及び解析時間（sec）等の付加情報を入力部２７を介して取得することができる。ここで、解析時間は、予測の対象となる時間の長さであり、時間濃度曲線（ＴＤＣカーブ）４３を示すグラフ（図６）のＸ軸の長さに対応する。また、オペレーターが局所灌流を考慮することを選択した場合、予測部１６は、組織内で毛細血管から細胞外液腔への造影剤の染み出しを考慮する。

【0024】

そして、予測部１６は、被写体の情報と、注入プロトコルと、組織の情報とに基づき、組織を血流方向に沿って分割した複数のコンパートメントのそれぞれについて、画素値の経時変化を予測する。その後、予測部１６は、時間ごとの各コンパートメントの画素値を、各組織と関連付けてシミュレータ２０の記憶部２４に記憶させる。

【0025】

また、シミュレータ２０はＣＰＵ等の制御部２５を備え、予測部１６による予測結果を記憶する記憶部２４は制御プログラム等を記憶している。そして、制御部２５は、記憶部２４に記憶された制御プログラムに従ってシミュレータ２０を制御している。また、制御部２５は、被写体情報取得部１１、プロトコル取得部１２、組織情報取得部１３、薬液情報取得部１４、及び表示部２６を制御する表示制御部１５を有している。そして、記憶部２４に実装された制御プログラムに対応して制御部２５が各種処理を実行することにより、各部が各種機能として論理的に実現される。

【0026】

また、記憶部２４には、被写体の組織における画素値の経時変化をコンピューター（制御部）に予測させるシミュレーションプログラムが記憶されている。このシミュレーションプログラムは、コンピューターを、被写体の情報を取得する被写体情報取得部１１、造影剤の注入プロトコルを取得するプロトコル取得部１２、組織の情報を取得する組織情報取得部１３、及び被写体の情報と、注入プロトコルと、組織の情報とに基づいて、組織を血流方向に沿って分割した複数のコンパートメントのそれぞれの画素値の経時変化を予測する予測部１６として機能させる。なお、このシミュレーションプログラムは、コンピューター読み取り可能な記録媒体に記憶させることができる。

【0027】

また、記憶部２４は、制御部２５が動作するためのシステムワークメモリであるＲＡＭ（Random Access Memory）、プログラム若しくはシステムソフトウェア等を格納するＲＯＭ（Read Only Memory）、又はハードディスクドライブ等を有する。なお、制御部２５は、ＣＤ（Compact Disc）及びＤＶＤ（Digital Versatile Disc）、ＣＦ（Compact Flash）カード等の可搬記録媒体、又はインターネット上のサーバー等の外部記憶媒体に記憶されたプログラムに従って各種処理を制御することもできる。

【0028】

さらに、シミュレータ２０は、各組織のコンパートメントを画素値に応じた濃度の色で表示する表示部２６を備えている。そして、表示制御部１５は、画素値の経時変化に応じて、表示部２６に表示された各組織のコンパートメントの濃淡を変更する。そのために、表示制御部１５は、記憶部２４から所定の時間におけるコンパートメントの画素値を読み出して、コンパートメントの濃淡を変更する。なお、表示部２６には、入力画面等の操作画面が表示される。また、注入プロトコル、装置の入力状態、設定状態、及び注入結果等の各種情報が表示されてもよい。

【0029】

また、シミュレータ２０の入力部２７は、被写体情報取得部１１、プロトコル取得部１２、組織情報取得部１３、及び薬液情報取得部１４に接続されている。なお、入力部２７としてはキーボード等を用いることができるが、タッチパネルを用いることにより入力部２７と表示部２６とを兼用させることもできる。

【0030】

上述したシミュレータ２０は、後述する図８に記載された、医療用の撮像装置３を備える撮像システム１００、又は造影剤を注入するための注入装置２に搭載されることができる。また、シミュレータ２０は、撮像装置３又は注入装置２に有線又は無線接続される外部コンピューターに搭載することもできる。なお、撮像装置３としては、例えば、ＭＲＩ（Magnetic Resonance Imaging）装置、ＣＴ（Computed Tomography）装置、アンギオ撮像装置、ＰＥＴ（Positron Emission Tomography）装置、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置、ＣＴアンギオ装置、ＭＲアンギオ装置、超音波診断装置、及び血管撮像装置等の各種医療用の撮像装置があるが、本明細書ではＣＴ装置について説明している。

【0031】

次に、図２を参照して予測部１６による画素値の経時変化の予測について説明する。この図２の上側には、血流方向Ｄにおいて直列に連通する血管Ａ１、臓器Ａ２、血管Ａ３及び臓器Ａ４が、それぞれ１つのコンパートメントに対応する場合の模式図が示されている。この場合、予測部１６は、造影剤が各組織に到達した直後に、それぞれの組織全体の画素値が変化するものとして予測を行う。そして、実際の拡散速度及び位置に関わらず、組織全体の画素値が変化したものと予測されてしまう。

【0032】

すなわち、図２の臓器Ａ４では、造影剤が到達した直後であり、実際には臓器Ａ４のほとんどの部分で画素値の変化が生じていないが、臓器Ａ４全体の画素値が増加したもの（白）として予測されてしまう。また、血管Ａ３では、実際には造影剤が臓器Ａ４にほとんど移動してないが、血管Ａ３全体の画素値が減少したもの（灰色）として予測されてしまう。さらに、臓器Ａ２では、実際には造影剤が残っているが、臓器Ａ２全体の画素値が減少したもの（黒）として予測されてしまう。そのため、特に実際の組織内での拡散速度が遅い場合には、組織の画素値の変化を正確に予測することができなくなってしまう。

【0033】

一方、図２の下側には第１実施形態に係る模式図が示されている。そして、血流方向Ｄにおいて直列に連通する血管Ｂ１、臓器Ｂ２、血管Ｂ３及び臓器Ｂ４が、それぞれ血流方向Ｄに沿って１５個の複数のコンパートメントに分割されている。すなわち、各組織が、組織情報取得部１３から取得された組織の分割コンパートメント数に応じて、血流方向Ｄに沿って複数のコンパートメントに分割されている。

【0034】

この場合、予測部１６は、予測対象のコンパートメントを含む組織の体積と、当該組織の毛細血管体積と、当該組織の細胞外液腔体積とを、分割コンパートメント数で除算して、各コンパートメントについて予測を行う。例えば、分割コンパートメント数が１５である場合、予測部１６は、組織の体積と、毛細血管体積と、細胞外液腔体積とのそれぞれを、１５で除算した値に基づいて画素値の経時変化を予測する。

【0035】

そのため、実際にはほとんどの部分で画素値の変化が生じていない臓器Ｂ４では、血管Ｂ３側に位置するコンパートメントのみの画素値が増加したもの（白）として予測される。同様に、血管Ｂ３では、実際には造影剤が臓器Ｂ４にほとんど移動してないので、臓器Ｂ２側に位置するコンパートメントのみの画素値が減少したもの（灰色）として予測される。さらに、臓器Ｂ２では、実際には造影剤が残っているので、血管Ｂ１側に位置するコンパートメントのみの画素値が減少したもの（黒）として予測される。これにより、実際の組織内での拡散速度が遅い場合であっても、組織の画素値の変化を正確に予測することができる。

【0036】

次に、図３を参照して具体的な画素値の変化の予測について説明する。図３に示すように、第１実施形態における被写体の組織は、右心室、大動脈、静脈、動脈、脳（頭）、上肢、右冠動脈（右冠動脈が支配的な心筋）、前下行枝（前下行枝が支配的な心筋）、回旋枝（回旋枝が支配的な心筋）、肺、肝臓、胃、脾臓、膵臓、腸管、腎臓、下肢、左心室、上行大動脈、下行大動脈、及び腹部大動脈を含む。そして、上肢静脈から注入された造影剤は、右心室、肺、左心室、大動脈（上行大動脈、下行大動脈）を介して各臓器に移動し、その後、静脈を介して右心室に到達する。そして、体内に注入された造影剤は、腎臓を介して体外に排出される。

【0037】

予測部１６は、血流方向における上流及び下流のそれぞれに向かって右心室から順に各組織の画素値の経時変化を予測する。すなわち、予測部１６は、最初に１－右心室の予測を行い、次いで右心室に対して血流方向の上流側に位置する２－大静脈及び２－静脈と、右心室に対して血流方向の下流側に位置する２－動脈とを含む第２組織群の予測を行う。その後、予測部１６は、静脈、脳、上肢、右冠動脈、前下行枝、回旋枝及び肺を含む第３組織群、動脈、静脈及び肝臓を含む第４組織群、左心室、動脈、上行大動脈、下行大動脈、腎臓及び下肢を含む第５組織群、腹部大動脈、胃、脾臓、膵臓、腸管及び動脈を含む第６組織群、動脈を含む第７組織群の順に予測を行う。なお、図３においては、各組織の名称の前にハイフンを介して付した数字によって予測の順番を示している。

【0038】

なお、予測部１６は、血流方向における上流及び下流のそれぞれに向かって造影剤の注入箇所に近い組織から順に各組織の画素値の経時変化を予測してもよい。例えば、造影剤を肝動脈に注入する場合（ＣＴＨＡ）、予測部１６は、最初に肝臓の予測を行ってもよい。その後、予測部１６は、肝臓に対して血流方向の上流側に位置する動脈と、肝臓に対して血流方向の下流側に位置する静脈とを含む組織群の予測を行う。

【0039】

また、予測部１６は、各組織（血管及び臓器）における画素値の変化を時間関数として求めるために、例えば下記式１のような微分方程式を用いている。なお、コンパートメントに流入する造影剤の濃度をＣ_１とし、コンパートメントから流出する造影剤の濃度をＣ_２とし、コンパートメントの体積をＶとし、コンパートメントにおける単位組織あたりの血流量（血流速度）をＱとしている。

【0040】

【数1】

【0041】

さらに、予測部１６は、右心室、左心室、血管以外の組織における画素値の変化を求めるために、造影剤が毛細血管から細胞外液腔へと透過する際の染み出し速度と、細胞外液腔から毛細血管へと透過する際の染み戻り速度とを考慮する。そのため、予測部１６は、例えば下記式２及び式３のような微分方程式を用いている。なお、細胞外液腔の体積をＶｅｃとし、細胞外液腔の造影剤の濃度をＣｅｃとし、毛細血管の体積をＶｉｖとし、毛細血管の造影剤の濃度をＣｉｖとし、染み出し速度をＰＳ_１とし、染み戻り速度をＰＳ_２としている。

【0042】

【数2】

【0043】

【数3】

【0044】

そして、上記微分方程式を解くことにより、注入開始からの経過時間と画素値（造影剤濃度）の変化が時間関数として求められる。この予測に用いられるパラメータの例として、図４には、胃、脾臓、膵臓及び腸管の各組織についての値を表に示している。第１実施形態において、予測部１６は、胃と、脾臓と、膵臓と、腸管とをそれぞれ別の組織として画素値の経時変化を予測している。

【0045】

図４に示すように、胃については、１２０mL以上１６０mL以下の組織体積、２mL以上５mL以下の毛細血管体積、１２mL以上１８mL以下の細胞外液腔体積、１２０mL/min以上１８０mL/min以下の単位組織あたり血流量(動脈血流速度)、１５以上２５以下の染み出し速度、及び１５以上２５以下の染み戻り速度が用いられる。

【0046】

また、脾臓については、１２０mL以上１６０mL以下の組織体積、１０mL以上１５mL以下の毛細血管体積、４５mL以上６５mL以下の細胞外液腔体積、１５０mL/min以上２５０mL/min以下の単位組織あたり血流量、１５以上２５以下の染み出し速度１５－２５、及び１５以上２５以下の染み戻り速度が用いられる。

【0047】

また、膵臓については、１２０mL以上１５０mL以下の組織体積、３mL以上６mL以下の毛細血管体積、３０mL以上５０mL以下の細胞外液腔体積、１２０mL/min以上１８０mL/min以下の単位組織あたり血流量、１５以上２５以下の染み出し速度、及び１５以上２５以下の染み戻り速度が用いられる。

【0048】

また、腸管については、１８００mL以上２０００mL以下の組織体積、３０mL以上４０mL以下の毛細血管体積、５００mL以上６００mL以下の細胞外液腔体積、０．４mL/min以上０．５mL/min以下の単位組織あたり血流量、１５０以上２５０以下の染み出し速度、及び１５０以上２５０以下の染み戻り速度が用いられる。

【0049】

ここで、染み出し速度及び染み戻り速度は、毛細血管面積と透過性（permeability）の積により算出することができる。例えば、人体内の毛細血管の総面積を８００m²と仮定して、各臓器ごとにその重量に応じた毛細血管面積を割り当てる。そして、全ての臓器の透過性を１ml/min/gと仮定して、染み出し速度及び染み戻り速度を求めることができる。

【0050】

予測部１６は、予測結果を記憶部２４（図１）に順次記憶させる。この予測結果は、組織に関連付けられた時間ごとの画素値の情報が含まれている。そして、表示部２６（図１）は、複数のコンパートメントを含む各組織の予測画像を模式的に表示する。さらに、表示制御部１５は、記憶部２４から画素値を読み出して、画素値の経時変化に応じて各コンパートメントの濃淡を変更するように表示部２６を制御する。

【0051】

図５は、一例としての予測画像における濃淡の変化を示しており、体の頭尾方向の水平断面に対応する。ただし、実際の断面とは異なり、各組織を一覧できるように全ての組織を示している。また、ウィンドウ幅３５０及びウィンドウレベル４０が設定されている。なお、ウィンドウ幅は、画素値のコントラストの範囲に対応し、ウィンドウレベルは画面の明るさに対応する。そして、ウィンドウレベルの値からウィンドウ幅の値の半分を減算して得られた値よりも画素値が小さい場合、表示部２６は黒で表示する。また、ウィンドウレベルの値にウィンドウ幅の値の半分を加算して得た値よりも画素値が大きい場合、表示部２６は白で表示する。

【0052】

図５の上側には、上肢静脈に造影剤を注入した直後の各組織の画像Ｎ１を示しており、上肢静脈は生来の画素値（濃い灰色）で示されている。そして、腹部大動脈及び腹腔動脈を含む全ての血管にはまだ造影剤が到達しておらず、全ての血管は生来の画素値（濃い灰色）で示されている。次に、図５の中央には、注入開始から約２５秒経過したときの各組織の画像Ｎ２が示されており、腹部大動脈及び腹腔動脈、内頸静脈などが特に白く染まっている。一方、上肢静脈は、造影剤が早い段階で流出したために画素値が減少して薄い灰色で示されている。

【0053】

そして、図５の下側には、注入開始から約１２０秒経過したときの各組織の画像Ｎ３が示されている。全身の血管及び臓器に造影剤が拡散し一様に分布したことで、画像Ｎ２と比較すると画素値が低下して全体が薄い灰色で示されている。

【0054】

続いて図６を参照して、表示部に表示される一例としての操作画面について説明する。図６の右側に示すように、予測部１６による予測が完了すると、表示制御部１５は、所定の時間、例えばオペレーターが選択した時間における各コンパートメントの画素値を記憶部２４から読み出す。そして、表示制御部１５は、読み出した画素値を予測画像４１に反映させ、表示部２６に表示させる。例えば、図６においては、９．９０秒の時点が選択され、当該時点の予測画像４１が表示されている。なお、初期設定では、注入時、つまり０秒の時点における予測画像４１が表示される。

【0055】

図６の予測画像４１においては、ウィンドウ幅３５０及びウィンドウレベル４０が設定され、予測画像４１の右上隅にそれぞれの値が表示されている。また、ウィンドウ幅及びウィンドウレベルの下には、画素値として－１０００HUが表示されている。これは、予測画像４１において、矢印で示されたポインターが指し示す部分に対応する位置の画素値を表示している。なお、オペレーターは、ウィンドウ幅（WW）、及びウィンドウレベル（ＷＬ）を、入力部２７を介して入力することができる。

【0056】

予測画像４１の下方には、操作ボタン４２が表示されている。この操作ボタン４２には、図６の左から順に停止ボタン、再生ボタン、２倍速再生ボタン、３倍速再生ボタン及び１０倍速再生ボタンが含まれている。そして、オペレーターが再生ボタンを選択すると、経過時間に沿って予測画像４１が連続的に動画として再生される。これにより、オペレーターは、所望の時間における、各組織内の造影剤の位置を視覚により認識することができる。

【0057】

予測画像４１の左側には、時間濃度曲線４３が表示されている。この時間濃度曲線４３は、Ｘ軸（横軸）が注入開始からの経過時間に対応し、Ｙ軸（縦軸）が画素値に対応する。また、Ｙ軸には第１軸と第２軸とがあり、複数の組織を１つのグラフ上で表示させることができる。例えば、図６では、肝臓、門脈及び肝動脈の時間濃度曲線が表示されている。

【0058】

オペレーターは、時間濃度曲線４３に表示する組織を、時間濃度曲線４３の下側にある表示オプション４４において選択できる。例えば、図６では、脳、肝臓、門脈、肝動脈及び右心室の中から、肝臓、門脈及び肝動脈が選択されている。また、体内から排出された造影剤の量を表示させることもできる。なお、他の選択可能な組織としては、上肢、肺、左心室、心筋、右冠動脈、前下行枝、回旋枝、気管支、脾臓、腸管、腎臓、下肢、膵臓、肺動脈、肺静脈、上行大動脈、下行大動脈、腹部大動脈、腹腔動脈、上腸間膜動脈、下腹部大動脈、肝静脈、腎動脈、腎静脈、脳動脈、脳静脈、上肢動脈、上肢静脈、下肢動脈、下肢静脈、上部大静脈及び下部大静脈等が挙げられる。

【0059】

また、表示オプション４４の右側には、Ｙ軸入力ボックス４５が表示されている。オペレーターは、Ｙ軸入力ボックス４５において自動設定を選択するか、又は所望の値を入力することができる。例えば、図６では、第１軸について自動設定が選択されており、第１軸の最大値と最小値が自動で設定されている。また、第２軸については最大値として１００が入力され、最小値として５０が入力されている。

【0060】

また、Ｙ軸入力ボックス４５の下側には、現在時点ボックス４６が表示されている。そして、オペレーターが現在時点ボックス４６を選択すると、時間濃度曲線４３に現在時点バー４７が表示される。この現在時点バー４７は、予測画像４１に対応する時間（図６では９．９０秒）を示している。そして、現在時点バー４７は、予測画像４１が連続的に再生されると、経過時間に対応してＸ軸に沿って移動する。

【0061】

さらに、現在時点ボックス４６の下側には、ヘリカルスキャンボックス４８が表示されている。そして、オペレーターは、ヘリカルスキャンボックス４８を選択して、寝台移動速度（cm/sec）を入力することができる。なお、図６では、８．０cm/secの寝台移動速度が入力されている。

【0062】

ヘリカルスキャンボックス４８が選択されると、表示制御部１５は、ヘリカルスキャンによる遅延時間を取得する。ここで、遅延時間は、頭を撮像してから各組織を撮像するまでの経過時間（寝台の移動時間）に対応し、予測画像４１の上端から各組織までの長さに基づいて求められる。そして、表示制御部１５は、所定の時間に遅延時間を加算した時間における画素値を記憶部２４から読み出す。すなわち、表示制御部１５は、取得した遅延時間を所定の時間（現在時点）に加算して得られた時間における各組織の画素値を読み出す。例えば、図６の予測画像４１においては、頭（脳）は現在時点である９．９０秒経過した時点の画素値を示し、右心室は１４．９０秒経過した時点の画素値を示している。

【0063】

これにより、ヘリカルスキャンを行った際の予測画像を得ることができる。例えば、現在時点を注入開始直後の時点（０秒）とした場合に、脳については注入開始直後の画素値を示し、右心室については注入開始から５秒経過した時点の画素値を示し、肝臓については注入開始から７．５秒経過した時点の画素値を示すことができる。なお、表示制御部１５は、算出することによって遅延時間を取得してもよい。また、寝台移動速度と関連付けた遅延時間を予め記憶部２４に記憶しておき、表示制御部１５が記憶部２４から遅延時間を取得してもよい。

【0064】

このようにして得られた予測画像について、図７を参照して説明する。図７は、ヘリカルスキャン時の一例としての予測画像の濃淡の変化を示している。なお、図７の画像は、体の頭尾方向の水平断面に対応するが、実際の断面とは異なり全ての組織を示している。また、ウィンドウ幅３５０及びウィンドウレベル４０が設定されている。

【0065】

図７の上側には、上肢静脈に造影剤を注入した直後の各組織の画像Ｈ１を示している。しかし、図５の画像Ｎ１とは異なり、上肢静脈は遅延時間が加算された時点の画素値を示している。そのため、上肢静脈が造影剤により白く染まっている。そして、その他の血管にはまだ造影剤が到達しておらず、その他の血管は生来の画素値（濃い灰色）で示されている。

【0066】

次に、図７の中央には、注入開始から約２５秒経過したときの各組織の画像Ｈ２が示されており、内頸静脈は造影剤が到達したために白く染まっている。しかし、図５の画像Ｎ２とは異なり、遅延時間が加算された時点の画素値を示しているため、造影剤が既に流出した腹部大動脈及び腹腔動脈は、画素値が減少して薄い灰色で示されている。また、造影剤が循環して門脈に到達したことから、門脈が白く染まっている。

【0067】

そして、図７の下側には、注入開始から約１２０秒経過したときの各組織の画像Ｈ３が示されており、全身の血管及び臓器に造影剤が拡散し一様に分布したことで、画像Ｈ２と比較すると画素値が低下して全体が薄い灰色で示されている。

【0068】

図６の説明に戻り、ヘリカルスキャンボックス４８の下側には、解析ボタン４９が表示されている。そして、オペレーターが解析ボタン４９を選択すると、予測部１６が画素値の予測を開始する。なお、予測部１６は、被写体情報と、注入プロトコルと、組織情報とを取得した時に画素値の予測を開始してもよい。

【0069】

続いて図８を参照して、シミュレータ２０を備える撮像システム１００について説明する。撮像システム１００は、造影剤を注入するための注入装置２と、注入装置２に有線又は無線で接続され且つ被写体を撮像する医療用の撮像装置３とを備えている。そして、注入装置２又は撮像装置３は、上述のシミュレータ２０を備えている。

【0070】

撮像装置３は、撮像プランに従って被験者の撮像を行う撮像部３１と、撮像装置３の全体を制御する制御装置３２と、表示部２６としてのディスプレイ３３とを有している。なお、制御装置３２とディスプレイ３３とは、一体的に構成することもできる。また、撮像装置３は、注入装置２と有線又は無線で接続され、例えば、ゲートウェイ装置（不図示）を介して接続される。

【0071】

撮像装置３の撮像プランには、撮像部位、実効管電圧、機種名、メーカー名、撮像時間、管電圧、撮像範囲、回転速度、ヘリカルピッチ、曝射時間、線量及び撮像方法等の情報を含めることができる。そして、制御装置３２は、撮像プランに従うように撮像部３１を制御して被験者を撮像する。また、制御装置３２は、ディスプレイ３３に接続されており、ディスプレイ３３には、装置の入力状態、設定状態、撮像結果、及び各種情報等が表示される。

【0072】

撮像部３１は、寝台と、被写体である被験者にＸ線を照射するＸ線源と、被験者を透過したＸ線を検出するＸ線検出器等を有している。そして、撮像部３１は、被験者にＸ線を曝射し、被験者を透過したＸ線に基づいて被験者の体内を逆投影することで、被験者の透視画像を撮像する。なお、撮像部３１は、Ｘ線に代えて、ラジオ波又は超音波を用いて撮像してもよい。また、制御装置３２は、撮像部３１及び注入装置２等と有線又は無線で通信を行うことができる。

【0073】

造影剤を注入するための注入装置２は、シリンジに充填された薬液、例えば、各種造影剤及び生理食塩水等を被写体としての被験者の体内に注入する。また、注入装置２は注入プロトコルに従って造影剤を注入する注入ヘッド２１を備えている。さらに、注入装置２は、注入ヘッド２１を保持するスタンド２２と、注入ヘッド２１に有線又は無線で接続されたコンソール２３とを有する。

【0074】

コンソール２３は、注入ヘッド２１を制御する制御装置として機能すると共に、シミュレータ２０としても機能する。また、コンソール２３は、入力部２７及び表示部２６として機能するタッチパネルを備え、注入ヘッド２１及び撮像装置３等と有線又は無線で通信を行うことができる。なお、注入装置２は、タッチパネルに代えて、表示部２６としてのディスプレイと、入力部２７としてのキーボード等のユーザインタフェースとを備えていてもよい。

【0075】

また、シミュレータ２０、入力部２７及び表示部２６は、それぞれ別体で構成することができる。例えば、注入装置２は、コンソール２３に代えて、注入ヘッド２１に接続された制御装置と、該制御装置に接続され且つ薬液の注入状況等が表示される表示部２６（タッチパネルディスプレイ等）とを有していてもよい。このような制御装置も、シミュレータ２０として機能する。また、注入ヘッド２１及び制御装置は、スタンド２２と一体的に構成することもできる。また、スタンド２２に代えて天吊部材を設け、該天吊部材を介して天井から注入ヘッド２１を天吊することもできる。

【0076】

また、注入装置２は、電源又はバッテリー、コンソール２３に接続されるハンドスイッチ、及び注入ヘッド２１を遠隔操作する遠隔操作装置等を有していてもよい。この遠隔操作装置は、注入ヘッド２１を遠隔操作して注入を開始又は停止することができる。また、電源又はバッテリーは、注入ヘッド２１又は制御装置（コンソール２３）のいずれかに設けることができ、これらとは別に設けることもできる。

【0077】

注入ヘッド２１は、造影剤が充填されたシリンジが搭載される第１保持部２１４と、造影剤を後押しするための薬液としての生理食塩水が充填されたシリンジが搭載される第２保持部２１５とを有する。また、注入ヘッド２１は、第１保持部２１４に搭載されたシリンジ内の薬液を注入プロトコルに従って押し出す駆動機構（不図示）と、第２保持部２１５に搭載されたシリンジ内の薬液を注入プロトコルに従って押し出す駆動機構（不図示）とを有している。

【0078】

また、注入ヘッド２１は、注入条件、注入状況、装置の入力状態、設定状態、及び各種注入結果等が表示されるヘッドディスプレイ２１１と、駆動機構の動作を入力するための操作部２１２とを有している。なお、ヘッドディスプレイ２１１は省略することができる。また、ヘッドディスプレイ２１１は、タッチパネル等から構成することによって、操作部２１２として用いることもできる。

【0079】

操作部２１２には、駆動機構の前進ボタン、駆動機構の後進ボタン、又は最終確認ボタン等が設けられている。そして、薬液を注入する際には、注入ヘッド２１に搭載されたシリンジの先端部にミキシングチューブ等が接続される。オペレーターは、ミキシングチューブの接続等の注入準備が完了すると最終確認ボタンを押す。これにより、注入ヘッド２１は、注入を開始できる状態で待機する。

【0080】

その後、シリンジから押し出された薬液は、ミキシングチューブ等を介して、被験者の体内へ注入される。このミキシングチューブは、造影剤と希釈薬液との混合器として機能する。なお、このような混合器としては、株式会社根本杏林堂製の「ＳＰＩＲＡＬ ＦＬＯＷ（登録商標）」等がある。

【0081】

注入ヘッド２１には、例えば、ＲＦＩＤチップ、ＩＣタグ、バーコード等のデータキャリアを有するプレフィルドシリンジ等の、種々のシリンジを搭載することができる。そして、注入ヘッド２１は、シリンジに取り付けられたデータキャリアの読み取りを行う読取部（不図示）を内蔵している。このデータキャリアには、薬液に関連する薬液情報が記憶されている。

【0082】

注入装置２は、内部又は外部のゲートウェイ装置を介して、不図示のサーバー（外部記憶装置）から情報を受信することができ、サーバーへ情報を送信することもできる。また、撮像装置３は、サーバーから情報を受信することができ、サーバーへ情報を送信することもできる。このサーバーには、予め検査オーダーが記憶されている。そして、検査オーダーは、被験者に関する被写体情報と、検査内容に関する検査情報とを含んでいる。また、サーバーは、撮像装置３から送信された画像のデータ等の撮像結果に関する情報と、注入装置２から送信された注入結果に関する情報を記憶することができる。

【0083】

このようなシミュレータ２０を備えた撮像装置３によれば、オペレーターは、ディスプレイ３３で予測画像４１を確認しながら撮像装置３を操作することができる。また、撮像装置３は、予測部１６による予測結果に応じて撮像プランを変更してもよい。例えば、オペレーターが撮像装置３に所望の画素値を入力し、当該所望の画素値が、撮像対象となる組織について予測部１６が予測した画素値と異なる場合に、撮像装置３は予測結果が所望の画素値と一致するように、管電圧又は管電流等を変更してもよい。

【0084】

また、シミュレータ２０を備えた注入装置２によれば、オペレーターは、コンソール２３で予測画像４１を確認しながら注入装置２を操作することができる。また、注入装置２は、予測部１６による予測結果に応じて注入プロトコルを変更してもよい。例えば、オペレーターが注入装置２に所望の画素値を入力し、当該所望の画素値が、撮像対象となる組織について予測部１６が予測した画素値と異なる場合に、注入装置２は予測結果が所望の画素値と一致するように、注入速度又は注入時間等を変更してもよい。

【0085】

以上説明した第１実施形態に係る発明によれば、実際の組織における画素値の経時変化に近似したより高精度の予測を行うことができる。特に、造影剤の注入量が少ない場合、造影剤の注入時間が短い場合、又は造影剤の濃度が低い場合であっても、高精度の予測を行うことができる。また、各組織内における造影剤の位置を予測することができる。さらに、オペレーターは、所望の時間における、各組織内の造影剤の位置を視覚により認識することができる。

【0086】

また、第１実施形態に係る発明によれば、注入プロトコルが、造影剤の後押し注入又は注入速度の増減等を含んでいても各組織の画素値を予測することができる。例えば、造影剤の注入速度を徐々に減少させると同時に生理食塩水の注入速度を増加させる、いわゆるクロス注入法の場合であっても、実際の撮像画像を予測することができる。さらに、新しい注入プロトコル又は撮像プランを作成した場合に、実際の撮像画像を予測することができる。

【0087】

［第２実施形態］
第２実施形態について、図９を参照して説明する。第１実施形態では、各組織がそれぞれ同じ数のコンパートメントに分割されていた。一方、第２実施形態では、大きな体積を有する組織が、より多くの数のコンパートメントに分割される。すなわち、第２実施形態に係る予測部１６は、小さな体積を有する組織（第１組織）を分割した複数のコンパートメントのそれぞれにおける画素値の経時変化を予測すると共に、大きな体積を有する組織（第２組織）を、小さな体積を有する組織よりも多くの数に分割した複数のコンパートメントのそれぞれにおける画素値の経時変化を予測する。

【0088】

なお、第２実施形態の説明においては、第１実施形態との相違点について説明し、第１実施形態で説明した構成要素については同じ参照番号を付し、その説明を省略する。特に説明した場合を除き、同じ参照符号を付した構成要素は略同一の動作及び機能を奏し、その作用効果も略同一である。

【0089】

実際の臓器は、それぞれ異なる体積を有している。例えば図９の臓器Ａ８の体積は、臓器Ａ６の体積よりも数倍大きい。そのため、臓器Ａ８の各コンパートメントの体積は、臓器Ａ６の各コンパートメントの体積よりも数倍大きくなる。そして、第１実施形態と同様に予測を行った場合、予測画像４１においては、１つのコンパートメントの画素値の変化であっても、臓器Ａ６においては小さな変化（造影剤の移動距離が短い）として表示される。一方、臓器Ａ８においては、臓器Ａ６と比較して大きな変化（造影剤の移動距離が長い）として表示される。これにより、オペレーターは、視覚により各組織内の造影剤の位置を正確に認識することが困難となってしまう。

【0090】

そこで、第２実施形態においては、大きな体積を有する組織を、小さな体積を有する組織よりも多くの数のコンパートメントに分割している。具体的に、図９においては、臓器Ｂ６を３個のコンパートメントに分割し、臓器Ｂ８を１５個のコンパートメントに分割している。これにより、臓器Ｂ６及び臓器Ｂ８のそれぞれのコンパートメントの体積が近似するため、オペレーターは、視覚により各組織内の造影剤の位置を正確に認識することができる。

【0091】

このような第２実施形態においては、各組織の体積に応じた最適な分割コンパートメント数が予め求められ、記憶部２４に記憶されている。そして、組織情報取得部１３は、記憶部２４から分割コンパートメント数を取得し、予測部１６は取得された分割コンパートメント数に基づいて画素値の変化を予測する。

【0092】

以上説明した第２実施形態に係る発明によっても、実際の組織における画素値の経時変化に近似したより高精度の予測を行うことができる。特に、造影剤の注入量が少ない場合、造影剤の注入時間が短い場合、又は造影剤の濃度が低い場合であっても、高精度の予測を行うことができる。また、各組織内における造影剤の位置を予測することができる。さらに、オペレーターは、所望の時間における、各組織内の造影剤の位置を視覚により認識することができる。

【0093】

また、第２実施形態に係る発明によれば、大きな体積を有する組織と小さな体積を有する組織とが予測画像４１に含まれる場合であっても、オペレーターは、視覚により各組織内の造影剤の位置を正確に認識することができる。

【0094】

なお、図９においては、血管Ａ５、血管Ａ７、血管Ｂ５及び血管Ｂ７がそれぞれ同じ数、すなわち１５個のコンパートメントに分割されている。しかし、大きな体積を有する血管を、小さな体積を有する血管と比較してより多くの数のコンパートメントに分割してもよい。また、大きな体積を有する組織Ｂ８と小さな体積を有する組織Ｂ６の各コンパートメントの体積が略一致するように、分割コンパートメント数を設定してもよい。

【0095】

以上、各実施形態を参照して本発明について説明したが、本発明は上記実施形態に限定されるものではない。本発明に反しない範囲で変更された発明、及び本発明と均等な発明も本発明に含まれる。また、上述の実施形態及び各変形形態は、本発明に反しない範囲で適宜組み合わせることができる。

【0096】

例えば、表示部２６は、体の水平断面のみでなく、冠状断面の予測画像４１を表示させることもできる。また、組織の分割コンパートメント数は１５には限られず、２以上の任意の数が選択される。

【0097】

また、記憶部２４に予め画素値のノイズ情報を記憶させ、表示制御部１５が記憶部２４からノイズ情報を読み出して、各コンパートメントの予測画像４１にノイズ情報を加えてもよい。このノイズ情報としては、例えば、造影剤によって白く染まった組織同士の間に発生する放射状のノイズを示した画像等がある。ノイズを示した画像を予測画像４１に重ねわせて加えることにより、より実際の撮像画像に近似した予測画像４１を得ることができる。

【0098】

なお、上記各実施形態において、表示部２６は、複数の組織が血流方向において連続する模式図を構成するようにコンパートメントを配置して、各コンパートメントを画素値に応じた濃度の色で表示していた。しかし、表示部２６は、各組織がそれぞれ単独で表示されるようにコンパートメントを配置して、各コンパートメントを画素値に応じた濃度の色で表示してもよい。また、表示制御部１５は、各組織におけるコンパートメント数が異なるように、表示部２６を制御してもよい。この場合、表示制御部１５は、オペレーターが設定したコンパートメント数、又は予め記憶部２４に記憶されたコンパートメント数を含むように各組織を表示させる。なお、表示部２６は、各コンパートメントを白黒以外の色で表示してもよい。

【0099】

また、注入ヘッド２１は、２つのシリンジを保持するタイプには限定されず、３つ以上のシリンジ保持部を有するタイプ、又は１つのみのシリンジ保持部を有するタイプであってもよい。さらに、表示制御部１５は、図６に示す予測画像４１において、所定の組織が最大画素値に到達した時に、所定の組織の画像上又は所定の組織の画像近傍に読み出した最大画素値を表示させてもよい。例えば、表示オプション４４において肝臓、門脈及び肝動脈が選択されている場合、表示制御部１５は、予測画像４１における肝臓、門脈及び肝動脈の近傍に、選択されたそれぞれの組織の最大画素値を表示させてもよい。なお、最大画素値は、白黒以外の色、例えば、青色、緑色、赤色又は黄色で表示させることができる。

【0100】

また、表示制御部１５は、図６に示す予測画像４１において、各組織（コンパートメント）を白黒の濃淡（グレースケール）以外の色、例えば、青色、緑色、赤色又は黄色の濃淡で表示してもよい。さらに、表示制御部１５は、所定の組織を、白黒以外の色の濃淡で表示してもよい。例えば、表示オプション４４において肝臓、門脈及び肝動脈が選択されている場合、表示制御部１５は、予測画像４１における肝臓、門脈及び肝動脈をそれぞれ、赤色、青色及び緑色の濃淡で表示させることができる。

【0101】

［変形形態］
予測部１６は、薬液の注入による単位組織あたりの血流量（血流速度）の変化を考慮してもよい。すなわち、薬液を注入すると、注入箇所よりも血流方向おいて下流にある組織（コンパートメント）では、注入された薬液に押されて血流速度が変化する。特に、通常の血流速度よりも速い速度で薬液を注入した場合、下流にある組織では血流速度が上昇する。そこで、予測部１６は、薬液の注入速度が通常の血流速度よりも速い場合には、注入速度から血流速度を差し引いて得られた差分を血流速度に加算し、血流速度の上昇を考慮することができる。

【0102】

そのため、予測部１６は、加算して得られた血流速度に基づいて画素値の経時変化を予測する。すなわち、予測部１６は、上述した数式を用いて画素値の経時変化を予測する際に、コンパートメントにおける単位組織あたりの血流速度Ｑに得られた差分を加算する。具体的には、図１０の加算処理のフローチャートに示すように、予測部１６は被写体情報取得部１１から被写体情報として注入箇所に対応する組織、すなわち薬液が注入される組織の通常の血流速度を取得する（Ｓ１０１）。また、予測部１６は、プロトコル取得部１２から受け取った注入プロトコルに含まれる薬液の注入速度を取得する（Ｓ１０２）。なお、予測部１６は、注入速度の後に血流速度を取得してもよい。また、予測部１６は、被写体情報取得部１１からの取得に代えて、被験者の体重から算出することによって血流速度を取得しもよい。

【0103】

次に、予測部１６は、血流速度と注入速度を比較して、注入速度が血流速度よりも速いか否かを判断する（Ｓ１０３）。そして、注入速度が血流速度以下である場合には（Ｓ１０３でＮＯ）、予測部１６が加算を行わずに処理を終了する。一方、注入速度が血流速度よりも速い場合には（Ｓ１０３でＹＥＳ）、予測部１６は、注入速度から注入箇所に対応する組織の血流速度を差し引いて差分を算出する（Ｓ１０４）。そして、予測部１６は、算出した差分を血流速度に加算し（Ｓ１０５）、処理を終了する。例えば、注入速度が２．０mL/secであり、血流速度が１．５mL/secである場合、予測部１６は、以下のように０．５mL/secを血流速度に加算する。そして、速度上昇後の血流速度として２．０mL/secを適用する。
１．５mL/sec＋（注入速度２．０mL/sec－１．５mL/sec）＝２．０mL/sec

【0104】

ここで、注入速度は単位時間当たりの薬液の注入量である。そのため、造影剤と生理食塩水とを同時に注入する場合、予測部１６は、両者の合計注入速度に基づいて加算を行う。また、通常の血流速度よりも速い速度の生理食塩水で後押し注入する場合、予測部１６は、生理食塩水の注入速度に基づいて加算を行う。加算処理の終了後、予測部１６は、加算して得られた血流速度を含む被写体の情報と、注入プロトコルと、組織の情報とに基づき、組織を血流方向に沿って分割した複数のコンパートメントのそれぞれについて、画素値の経時変化を予測する。これにより、血流速度をより実際の速度に近づけて、画素値の予測精度を向上させることができる。

【0105】

また、加算処理において、予測部１６は、複数の組織の全てを加算対象として血流速度に差分を加算できる。しかし、予測部１６は、複数の組織の一部のみを加算対象としてもよい。具体的には、注入箇所に対応する組織に対して血流方向のすぐ下流にある組織から、血流方向において注入箇所に対応する組織に至るまでの組織のみについて血流速度に差分を加算してもよい。例えば、図１１を参照して、上肢静脈から薬液を注入する場合に複数の組織の一部のみの血流速度に差分を加算する例を説明する。図１１において、右心室、動脈、肺、静脈、左心室、上行大動脈、動脈、心筋（右冠動脈が支配的な心筋、前下行枝が支配的な心筋及び回旋枝が支配的な心筋）、及び静脈は、右心室を起点とする閉回路を構成している。すなわち、これらの組織は、注入箇所に対応する組織に対して、血流方向のすぐ下流にある右心室を起点とする閉回路を構成している。そこで、予測部１６は、当該閉回路に含まれる右心室から静脈までの組織のみを加算対象とし、例えば下肢の組織の血流速度には差分を加算しない。これにより、加算の影響を閉回路内に留めることができるので、シミュレーションを簡略化して算出時間を削減できる。

【0106】

また、予測部１６は、薬液の注入開始と同時に加算対象となる全組織の血流速度に差分を加算して画素値の経時変化を予測する。しかし、予測部１６は、注入開始から所定の時間を経過した時点で血流速度に差分が加算されるものとして画素値の経時変化を予測してもよい。すなわち、注入箇所に対応する組織の下流に位置する組織においては、注入開始から所定の時間が経過してから、血流速度が変化する場合がある。そこで、予測部１６は、注入箇所に対応する組織に対しては、薬液の注入開始と同時に血流速度に差分が加算されるものとして画素値の経時変化を予測する。一方、当該組織から血流方向において距離が遠い組織に対して、予測部１６は、距離に応じて長くなる時間を経過した時点で差分が加算されるものとして画素値の経時変化を予測する。

【0107】

図１２は、変形形態に係る加算処理を行って得られた予測結果であるタイムデンシティカーブを示すグラフである。なお、図１２は肝動脈のタイムデンシティカーブを示しており、横軸が注入開始からの経過時間に対応し、縦軸が画素値に対応している。また、被験者の身長が１７０cmに設定され、被験者の体重が７０kgに設定され、造影剤濃度が３００mgI/mLに設定されている。また、図１２では、１００mLの造影剤を２．０mL/secの注入速度で５０秒間注入し、続けて３０mLの生理食塩水を２．０mL/secの注入速度で１５秒間後押し注入した場合のタイムデンシティカーブを表示している。そして、加算処理を行った場合のタイムデンシティカーブを実線で示し、加算処理を行わない場合のタイムデンシティカーブを破線で示している。

【0108】

加算処理を行わない場合、造影剤が組織に到達する時間が遅くなる。そのため、加算処理を行った場合と比較すると、遅れたタイミングで造影剤が到達したかのようにシミュレートされてしまう。そのため、破線で示されたタイムデンシティカーブにおいては、図１２の矢印Ａに示すように画素値の第２のピークが生じている。一方、加算処理を行った場合には、造影剤の循環がシミュレートされる。そのため、実線で示されたタイムデンシティカーブに示されるように、第２のピークが生じない。

【0109】

このように、加算処理を行うことによって、より実際に近い画素値を予測できる。さらに、生理食塩水による後押し注入をシミュレートすることもできる。具体的には、後押し注入によれば、造影剤の注入終了後にうっ滞していた造影剤が、より早いタイミングで押し出される。そして、図１２の実線で示すタイムデンシティカーブでは、後押し注入がシミュレートされた結果、ピークにおける画素値の予測精度が向上し、より画素値が高くなっている。

【0110】

さらに、予測部１６は、隣り合うコンパートメント間での造影剤の拡散を考慮してもよい。すなわち、隣り合うコンパートメント間において造影剤の濃度差が存在する場合、濃度の高いコンパートメントから濃度の低いコンパートメントへと造影剤の拡散が生じる。そこで、予測部１６は、濃度の高いコンパートメントの造影剤濃度を低下させて、濃度の低いコンパートメントの造影剤濃度を増加させるように画素値の予測を行い、造影剤の拡散を考慮してもよい。

【0111】

ここで、予測部１６は、濃度差が大きい場合には、造影剤濃度の低下量及び増加量を大きくする。また、予測部１６は、薬液情報取得部１４から造影剤の浸透圧を取得し、浸透圧が大きいときには、造影剤濃度の低下量及び増加量を大きくする。さらに、予測部１６は、コンパートメント同士の接触面積が大きい場合には、造影剤濃度の低下量及び増加量を大きくする。例えば、予測部１６は、異なる組織のコンパートメント同士が隣り合う場合には、接触面積が小さくなるため、造影剤濃度の低下量及び増加量を小さくする。そして、同じ組織のコンパートメント同士が隣り合う場合には、接触面積が大きくなるため、予測部１６は、造影剤濃度の低下量及び増加量を大きくする。

【0112】

なお、予測部１６は、接触面積を算出して、算出した接触面積が大きい場合には、造影剤濃度の低下量及び増加量を大きくしてもよい。また、上述した各変形形態は、本発明に反しない範囲で他の実施形態又は変形形態と適宜組み合わせることができる。

【0113】

この出願は２０１４年１１月２７日に出願された日本国特許出願第２０１４－２４０００６号からの優先権を主張し、その全内容を引用してこの出願の一部とする。

【符号の説明】

【0114】

２：注入装置、３：撮像装置、２０：シミュレータ、１１：被写体情報取得部、１２：プロトコル取得部、１３：組織情報取得部、１５：表示制御部、１６：予測部、２１：注入ヘッド、２６：表示部、１００：撮像システム、Ｄ：血流方向

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】