(57)【要約】【課題】経験的なフィッティング処理に依ることなく一相注入のための適正な注入プロトコルを得る。【解決手段】テスト注入プロトコルに従って被検体に造影剤をテスト注入するとともに、被検体がプレスキャンされる（テスト検査：Ｓ１〜Ｓ３）。このテスト検査に基づき、テスト時間濃度曲線が作成される（Ｓ４，Ｓ５）。テスト時間濃度曲線に基づき、注入速度はテスト注入速度と等しい一定値とし、注入速度は種々の値に定めた場合の複数の推定時間濃度曲線が作成される（Ｓ７）。各推定時間濃度曲線において、所要の濃度維持時間を確保できる最大濃度値が求められる（Ｓ８）。この最大濃度値と所要濃度値との比に基づき、注入速度が定められる（Ｓ９）。さらに、注入速度および注入時間を乗じて造影剤注入量が求められ、造影剤注入量および注入速度よって規定される一相注入を表す複数の注入プロトコルが作成される（Ｓ１０）。【選択図】図１８
2007143880.tif 000002
【特許請求の範囲】
【請求項１】
撮影対象に対して、一定のテスト注入速度で所定のテスト注入時間に渡って造影剤をテスト注入するステップと、
前記撮影対象の関心領域を撮像し、造影剤のテスト注入に伴う画像濃度の時間変化を表すテスト時間濃度曲線を作成するステップと、
前記テスト時間濃度曲線に基づいて、注入速度を前記テスト注入速度と等しい一定値とし、前記テスト注入時間とは異なる所定の注入時間を設定した場合についての時間濃度曲線を推定して、推定時間濃度曲線を作成するステップと、
前記推定時間濃度曲線において所要の濃度維持時間を確保できる最大濃度値を求めるステップと、
前記最大濃度値および所要濃度値に基づいて、前記濃度維持時間にわたって前記所要濃度値を維持できるように、前記所定の注入時間に対応した注入速度を定めるステップと、
前記所定の注入時間およびこれに対応する注入速度に基づいて造影剤注入量を求め、この造影剤注入量およびこれに対応する前記注入速度の対を含む注入プロトコルを作成するステップとを含む、造影剤注入プロトコル決定方法。
【請求項２】
前記所定の注入時間は前記テスト注入時間の自然数倍であり、
前記推定時間濃度曲線を作成するステップは、前記テスト時間濃度曲線と、前記テスト時間濃度曲線を時間軸上で前記テスト注入時間の自然数倍だけずらした一つ以上のシフト時間濃度曲線とを重ね合わせることにより、前記所定の注入時間に対応した推定時間濃度曲線を作成するステップを含む、請求項１記載の造影剤注入プロトコル決定方法。
【請求項３】
前記推定時間濃度曲線を作成するステップは、
前記テスト時間濃度曲線に基づいて、所定の微小単位時間の造影剤注入に対応する時間濃度曲線であるインパルス応答曲線を求めるステップと、
前記所定の注入時間に応じて、前記インパルス応答曲線を時間軸上でずらして複数個重ね合わせて前記推定時間濃度曲線を求めるステップとを含む、請求項１記載の造影剤注入プロトコル決定方法。
【請求項４】
前記インパルス応答曲線を求めるステップは、
テスト注入時の注入速度の時間変化を表すテスト注入関数をフーリエ変換するステップと、
前記テスト時間濃度曲線をフーリエ変換するステップと、
前記テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることにより、注入速度の時間変化を表す注入関数と時間濃度曲線との関係を表す伝達関数を求めるステップと、
前記伝達関数を逆フーリエ変換して前記インパルス応答曲線を求めるステップとを含む、請求項３記載の造影剤注入プロトコル決定方法。
【請求項５】
前記推定時間濃度曲線を作成するステップは、
テスト注入時の注入速度の時間変化を表すテスト注入関数をフーリエ変換するステップと、
前記テスト時間濃度曲線をフーリエ変換するステップと、
前記テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることにより、注入速度の時間変化を表す注入関数と時間濃度曲線との関係を表す伝達関数を求めるステップと、
注入速度を前記テスト注入速度と等しい一定値とし、前記注入時間を前記所定の注入時間に設定した場合の注入関数である設定注入関数をフーリエ変換するステップと、
前記設定注入関数のフーリエ変換および前記伝達関数に基づいて、前記設定注入関数に対応する推定時間濃度曲線を作成するステップとを含む、請求項１記載の造影剤注入プロトコル決定方法。
【請求項６】
前記テスト時間濃度曲線のフーリエ変換に対して高周波成分を除去するフィルタ処理を行うステップをさらに含み、
前記伝達関数を求めるステップは、前記フィルタ処理されたテスト時間濃度曲線のフーリエ変換を用いて行われる、請求項４または５記載の造影剤注入プロトコル決定方法。
【請求項７】
前記所定の注入時間は、互いに異なる複数の注入時間を含み、
前記推定時間濃度曲線を作成するステップは、前記複数の注入時間にそれぞれ対応する複数の推定時間濃度曲線を作成するステップであり、
前記濃度維持時間を確保できる最大濃度値を求めるステップは、前記複数の推定時間濃度曲線のそれぞれについて、前記濃度維持時間を確保できる最大濃度値を求めるステップであり、
前記注入速度を求めるステップは、前記複数の推定時間濃度曲線のそれぞれについて求められた複数の最大濃度値と前記所要濃度値とに基づいて、前記複数の注入時間にそれぞれ対応する複数の注入速度を求めるステップであり、
前記注入プロトコルを作成するステップは、前記複数の注入速度およびそれらにそれぞれ対応する複数の造影剤注入量の複数の対をそれぞれ含む複数の注入プロトコルを作成するステップである、請求項１〜６のいずれかに記載の造影剤注入プロトコル決定方法。
【請求項８】
前記複数の注入プロトコルから、注入速度が所定範囲内の値である注入プロトコルを抽出するステップをさらに含む、請求項７記載の造影剤注入プロトコル決定方法。
【請求項９】
前記複数の注入プロトコルから、注入速度が最低の注入プロトコルを抽出するステップをさらに含む、請求項７または８記載の造影剤注入プロトコル決定方法。
【請求項１０】
前記複数の注入プロトコルから、造影剤注入量が最低の注入プロトコルを抽出するステップとをさらに含む、請求項７〜９のいずれかに記載の造影剤注入プロトコル決定方法。
【請求項１１】
撮影対象に対して一定のテスト注入速度で所定のテスト注入時間に渡って造影剤をテスト注入したときに関心領域を撮像して得られるテスト時間濃度曲線に基づいて本番検査のための造影剤注入プロトコルを演算する装置であって、
前記テスト時間濃度曲線に基づいて、注入速度を前記テスト注入速度と等しい一定値とし、前記テスト注入時間とは異なる所定の注入時間を設定した場合についての時間濃度曲線を推定して、推定時間濃度曲線を作成する推定時間濃度曲線作成手段と、
この推定時間濃度曲線作成手段によって作成された前記推定時間濃度曲線において所要の濃度維持時間を確保できる最大濃度値を求める最大濃度値演算手段と、
この最大濃度値演算手段によって演算された前記最大濃度値および所要濃度値に基づいて、前記濃度維持時間にわたって前記所要濃度値を維持できるように、前記所定の注入時間に対応した注入速度を定める注入速度演算手段と、
前記所定の注入時間およびこれに対応して前記注入速度演算手段によって演算された注入速度に基づいて造影剤注入量を求め、この造影剤注入量およびこれに対応する前記注入速度の対を含む注入プロトコルを作成する注入プロトコル作成手段とを含む、造影剤注入プロトコル演算装置。
【請求項１２】
撮影対象に造影剤を注入する造影剤注入手段と、
前記テスト注入速度およびテスト注入時間に対応したテスト注入プロトコルに従って前記撮影対象に造影剤が注入されるように前記造影剤注入手段を制御するテスト注入制御手段と、
請求項１１記載の造影剤注入プロトコル演算装置と、
この造影剤注入プロトコル演算装置によって演算された注入プロトコルに従って前記撮影対象に造影剤が注入されるように前記造影剤注入手段を制御する本番注入制御手段とを含む、造影剤注入装置。
【請求項１３】
請求項１２記載の造影剤注入装置と、
撮影対象を撮影して画像濃度に対応する信号を生成する撮影手段と、
テスト注入時に前記撮影対象の関心領域を撮影するように前記撮影手段を制御するテスト撮影制御手段と、
前記本番注入制御手段によって前記造影剤注入手段が制御されるときに、これに同期して前記撮影対象の関心領域を撮影するように前記撮影手段を制御する本番撮影制御手段とを含む、撮影システム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、被検体（典型例は人体）に造影剤を注入して当該被検体の関心領域の濃度を増強して行う撮影法に適用可能な造影剤注入プロトコル決定方法および造影剤注入プロトコル演算装置に関する。さらに、この発明は、そのような造影剤注入プロトコル演算装置を備えた造影剤注入装置および撮影システムに関する。
【背景技術】
【０００２】
動脈瘤や血管閉塞といった、血管やその周辺臓器の病変を発見する手段としてＣＴＡ（Computed Tomography Angiography；コンピュータ断層撮影血管造影法）がある。通常のＣＴ（Computed Tomography；コンピュータ断層撮影）によって得られたＣＴ画像は、臓器、血管、筋肉などの軟部組織におけるＸ線吸収値の差が小さく、低コントラストの画像となってしまう。したがって、これらの部位に関して満足な診断を行うことが困難である。それに対し、ＣＴＡは、Ｘ線吸収値が高い造影剤を静脈に投与して、注目する血管部位のＸ線吸収値が増加している状態で、連続してＣＴを行う手法であり、これによりコントラストを高めて、明瞭な臓器の像を得ることができる。
【０００３】
ＣＴＡにおける造影剤の注入は、一般的には、ヨード系造影物質を用い、その総注入量および注入速度を決定して行われる。そして、注入開始から撮影のためのスキャン開始までの遅延時間を変化させることによって、造影剤による画像濃度の増強が調節される（非特許文献１）。
しかし、造影剤の注入速度および注入量に基づいて設定される注入パラメータが適切でない場合、画像の濃度値の差が小さいために診断の役に立たなかったり、逆に濃度値が高すぎることによりコントラストがつかない画像となったりする可能性がある。ところが、造影剤注入による時間濃度変化（画像濃度の時間変化）は、患者の体重や抱えている疾患などによって大きく異なる。そのため、濃度値や造影剤の注入開始から注目領域の血管に造影剤が現れるまでの遅延時間のばらつきが大きく、各々の患者に対してＣＴＡにおける動脈の濃度変化を推定して造影剤注入パラメータを最適化することは、適切なＣＴ血管造影像を得る上で非常に重要な事項であるにも拘わらず、必ずしも容易ではない。
【０００４】
さらに、現在の医療の現場では、造影剤の投与は患者にとって負担となるにもかかわらず、安定した造影を行うために、患者に対して適正量よりやや多量の造影剤を投与しているという問題がある。
このような問題に対し、計算機による支援により時間濃度曲線を推定できれば、使用する造影剤量を適切に設定でき、過剰に投与することがなくなると期待できる。これにより患者の身体的、経済的負担も軽減できると考えられる。
【０００５】
時間濃度曲線を解析する従来研究としては、Freischmannらの手法（非特許文献１）が挙げられる。この手法は、ある患者に対して、テスト注入パラメータで造影剤を注入したときの時間濃度曲線（ＣＴ値の時間変化）が既知である場合に適用することができる。すなわち、この手法では、まず、テスト注入パラメータとそれによる時間濃度曲線とに対してフーリエ変換を適用し、その患者の注入に対する時間濃度変化の関係を伝達関数として導出する。その伝達関数を用いることで、任意の注入パラメータに対する時間濃度曲線を推定でき、また、逆に、ＣＴにおいて要求される時間濃度曲線を設定することで、それを得るための最適な注入パラメータを推定することができる。
【０００６】
しかし、この手法では、最適注入パラメータを求める際のフィッティング処理が経験的に行われている。すなわち、最適注入パラメータは、注入速度が時間に伴って非線形に変化する複雑な形状の関数によって表されるが、自動注入器を用いて行われる実際の注入では、時間経過に伴って注入速度を複雑に変化させていくことは非現実的である。そこで、注入速度を２段階に切り換える、いわゆる二相注入を行うこととして、最適注入パラメータが求まった後は、操作者の経験によって当該最適注入パラメータに適合する二相注入のパラメータを定める、フィッティング処理が行われる。
【０００７】
前記手法において二相注入が行われるようになっているのは、最適注入パラメータが二相注入にフィッティングしやすい形状となっていて、注入速度を終始一定に保持する一相注入に対するフィッティングを行うことができないからである。しかし、臨床上は、一相注入の方が現実的であるという側面もある。
【非特許文献１】Dominik Freischmann and Karl Hittmair: "Mathematical Analysis of Arterial Enhancement and Optimization of Bolus Geometry for CT Angiography Using the Discrete Fourier Transform", Journal of Computer Assisted Tomography 23(3):474-484, 1999
【非特許文献２】立入弘・稲邑清也監修、山下一也・速水昭宗編集: “診療放射線技術上巻改訂第10 版”, 南江堂,pp.83,pp.86,pp.175-177,pp.339,2001
【非特許文献３】日本医用画像工学会監修, 同・医用画像工学ハンドブック編集委員会編集, “医用画像工学ハンドブック”, 篠原出版,pp.220-221,pp.364-365,1994
【発明の開示】
【発明が解決しようとする課題】
【０００８】
そこで、この発明の目的は、経験的なフィッティング処理に依ることなく一相注入のための適正な注入プロトコルを得ることができる造影剤注入プロトコル決定方法および造影剤注入プロトコル演算装置を提供することである。
また、この発明の他の目的は、そのような造影剤注入プロトコル演算装置を用いた造影剤注入装置およびそれを用いた撮影装置を提供することである。
【課題を解決するための手段】
【０００９】
上記の目的を達成するための請求項１記載の発明は、撮影対象に対して、一定のテスト注入速度で所定のテスト注入時間に渡って造影剤をテスト注入するステップ（Ｓ２）と、前記撮影対象の関心領域を撮像し、造影剤のテスト注入に伴う画像濃度の時間変化を表すテスト時間濃度曲線を作成するステップ（Ｓ３〜Ｓ５）と、前記テスト時間濃度曲線に基づいて、注入速度を前記テスト注入速度と等しい一定値とし、前記テスト注入時間とは異なる所定の注入時間を設定した場合についての時間濃度曲線を推定して、推定時間濃度曲線を作成するステップ（Ｓ７）と、前記推定時間濃度曲線において所要の濃度維持時間を確保できる最大濃度値を求めるステップ（Ｓ８）と、前記最大濃度値および所要濃度値（濃度増強値）に基づいて、前記濃度維持時間にわたって前記所要濃度値を維持できるように、前記所定の注入時間に対応した注入速度を定めるステップ（Ｓ９）と、前記所定の注入時間およびこれに対応する注入速度に基づいて造影剤注入量を求め、この造影剤注入量およびこれに対応する前記注入速度の対を含む注入プロトコルを作成するステップ（Ｓ１０）とを含む、造影剤注入プロトコル決定方法である。なお、括弧内の英数字は後述の実施形態における対応構成要素等を表す。以下、この項において同じ。
【００１０】
この方法では、テスト注入は、一定のテスト注入速度で所定のテスト注入時間にわたって行われる。すなわち、このテスト注入は、いわゆる一相注入である。このテスト注入に対応するテスト時間濃度曲線が作成され、このテスト時間濃度曲線に基づき、注入速度をテスト注入速度と等しい一定値とし、かつ、注入時間はテスト注入時間とは異ならせた場合についての推定時間濃度曲線が作成される。さらに、推定時間濃度曲線において、所要の濃度維持時間を確保できる最大濃度値が求められる。そして、この最大濃度値と所要濃度値（本番検査の撮影時において必要とされる濃度値）とに基づいて、前記所定の注入時間に対応した注入速度が求められる。この注入速度は、前記濃度維持時間にわたって所要濃度値を維持できる注入速度であり、具体的には、前記最大濃度値に対する所要濃度値の比をテスト注入速度に乗じることによって求めることができる。さらに、前記所定の注入時間およびこれに対応する注入速度に基づいて造影剤注入量（注入時間および注入速度の積）が求められる。こうして、造影剤注入量およびこれに対応する注入速度によって規定される一相注入のプロトコルを得ることができる。
【００１１】
この方法においては、注入プロトコルを得る過程において、経験的なフィッティングを必要としていない。そのため、所要の濃度維持時間にわたって所要濃度値を維持することができる適正な注入プロトコルを確実に作成することができる。しかも、医療の現場において現実的に適用することができる一相注入プロトコルを作成できる。
請求項２記載の発明は、前記所定の注入時間は前記テスト注入時間の自然数倍であり、前記推定時間濃度曲線を作成するステップは、前記テスト時間濃度曲線と、前記テスト時間濃度曲線を時間軸上で前記テスト注入時間の自然数倍だけずらした一つ以上のシフト時間濃度曲線とを重ね合わせることにより、前記所定の注入時間に対応した推定時間濃度曲線を作成するステップ（Ｓ７）を含む、請求項１記載の造影剤注入プロトコル決定方法である。
【００１２】
この方法では、前記所定の注入時間をテスト注入時間の自然数倍に限定することにより、テスト時間濃度曲線を時間軸上でずらして重ね合わせる簡単な処理によって推定時間濃度曲線を作成することができる。
請求項３記載の発明は、前記推定時間濃度曲線を作成するステップは、前記テスト時間濃度曲線に基づいて、所定の微小単位時間の造影剤注入に対応する時間濃度曲線であるインパルス応答曲線（関数）を求めるステップ（Ｓ２１〜Ｓ２５）と、前記所定の注入時間に応じて、前記インパルス応答曲線を時間軸上でずらして複数個重ね合わせて前記推定時間濃度曲線を求めるステップ（Ｓ２６）とを含む、請求項１記載の造影剤注入プロトコル決定方法である。
【００１３】
この方法では、所定の微小単位時間の造影剤注入に対する撮影対象の応答を表すインパルス応答曲線が求められる。このインパルス応答曲線を時間軸上でずらして複数個重ね合わせることにより、任意の注入時間に対応する推定時間濃度曲線を作成することができる。
請求項４記載の発明は、前記インパルス応答曲線を求めるステップは、テスト注入時の注入速度の時間変化を表すテスト注入関数をフーリエ変換するステップ（Ｓ２１）と、前記テスト時間濃度曲線をフーリエ変換するステップ（Ｓ２２）と、前記テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることにより、注入速度の時間変化を表す注入関数と時間濃度曲線との関係を表す伝達関数を求めるステップ（Ｓ２４）と、前記伝達関数を逆フーリエ変換して前記インパルス応答曲線を求めるステップ（Ｓ２５）とを含む、請求項３記載の造影剤注入プロトコル決定方法である。
【００１４】
この方法では、テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることによって、注入関数と時間濃度曲線との関係を表す伝達関数が求められる。この伝達関数を逆フーリエ変換すれば、前記インパルス応答曲線が得られる。
請求項５記載の発明は、前記推定時間濃度曲線を作成するステップは、テスト注入時の注入速度の時間変化を表すテスト注入関数をフーリエ変換するステップ（Ｓ２１）と、前記テスト時間濃度曲線をフーリエ変換するステップ（Ｓ２２）と、前記テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることにより、注入速度の時間変化を表す注入関数と時間濃度曲線との関係を表す伝達関数（Ｓ２４）を求めるステップと、注入速度を前記テスト注入速度と等しい一定値とし、前記注入時間を前記所定の注入時間に設定した場合の注入関数である設定注入関数をフーリエ変換するステップ（Ｓ３１，Ｓ３２）と、前記設定注入関数のフーリエ変換および前記伝達関数に基づいて、前記設定注入関数に対応する推定時間濃度曲線を作成するステップ（Ｓ３３，Ｓ３４）とを含む、請求項１記載の造影剤注入プロトコル決定方法である。
【００１５】
この方法では、テスト注入関数およびテスト時間濃度曲線のフーリエ変換を求めるとともに、注入時間をテスト注入時間とは異なる所定の注入時間に設定した場合の注入関数である設定注入関数のフーリエ変換が求められる。そして、テスト注入関数およびテスト時間濃度曲線の各フーリエ変換の比をとることによって伝達関数が求められる。さらに、この伝達関数と前記設定注入関数のフーリエ変換とに基づき、前記設定注入関数に対応する推定時間濃度曲線が作成される。伝達関数は、任意の注入関数に対する時間濃度曲線の関係を表すものであるため、この伝達関数を用いることにより、時間濃度曲線を推定することができる。より具体的には、伝達関数および設定注入関数のフーリエ変換に基づいて、周波数領域の推定時間濃度曲線が作成され（Ｓ３３）、これをさらに逆フーリエ変換することにより、時間領域の推定時間濃度曲線を得ることができる（Ｓ３４）。
【００１６】
請求項６記載の発明は、前記テスト時間濃度曲線のフーリエ変換に対して高周波成分を除去するフィルタ処理を行うステップ（Ｓ２３）をさらに含み、前記伝達関数を求めるステップは、前記フィルタ処理されたテスト時間濃度曲線のフーリエ変換を用いて行われる、請求項４または５記載の造影剤注入プロトコル決定方法である。
この方法では、テスト時間濃度曲線のフーリエ変換に対して、高周波成分を除去するフィルタ処理が行われる。これにより、請求項４記載の発明と組み合わせる場合には、高周波成分を含まない良好なインパルス応答曲線を作成することができる。また、請求項５記載の発明と組み合わせる場合には、高周波成分を含まない伝達関数を得ることができるので、良好な推定時間濃度曲線を得ることができる。
【００１７】
請求項７記載の発明は、前記所定の注入時間は、互いに異なる複数の注入時間を含み、
前記推定時間濃度曲線を作成するステップは、前記複数の注入時間にそれぞれ対応する複数の推定時間濃度曲線を作成するステップであり、前記濃度維持時間を確保できる最大濃度値を求めるステップは、前記複数の推定時間濃度曲線のそれぞれについて、前記濃度維持時間を確保できる最大濃度値を求めるステップであり、前記注入速度を求めるステップは、前記複数の推定時間濃度曲線のそれぞれについて求められた複数の最大濃度値と前記所要濃度値とに基づいて、前記複数の注入時間にそれぞれ対応する複数の注入速度を求めるステップであり、前記注入プロトコルを作成するステップは、前記複数の注入速度およびそれらにそれぞれ対応する複数の造影剤注入量の複数の対をそれぞれ含む複数の注入プロトコルを作成するステップである、請求項１〜６のいずれかに記載の造影剤注入プロトコル決定方法である。
【００１８】
この方法により、複数の注入時間に対応する注入速度が求められ、さらにそれらに対応する造影剤注入量が求められる。その結果、造影剤注入量および注入速度の組をそれぞれ含む複数の注入プロトコルを得ることができる。この複数の注入プロトコルの中から、注入速度や造影剤注入量を考慮して、適切な注入プロトコルを選択すればよい。
請求項８記載の発明は、前記複数の注入プロトコルから、注入速度が所定範囲内の値である注入プロトコルを抽出するステップ（Ｓ１２）をさらに含む、請求項７記載の造影剤注入プロトコル決定方法である。
【００１９】
この方法では、注入速度が所定範囲内の値である注入プロトコルを抽出することができるので、前記複数の注入プロトコルを、さらに、実際に適用可能な注入プロトコルに絞り込むことができる。
請求項９記載の発明は、前記複数の注入プロトコルから、注入速度が最低の注入プロトコルを抽出するステップ（Ｓ１２）をさらに含む、請求項７または８記載の造影剤注入プロトコル決定方法である。
【００２０】
この方法によれば、注入速度が最低の注入プロトコルが抽出されるので、撮影対象（患者）に対する負担の少ない注入プロトコルを選択することができる。
請求項１０記載の発明は、前記複数の注入プロトコルから、造影剤注入量が最低の注入プロトコルを抽出するステップ（Ｓ１２）とをさらに含む、請求項７〜９のいずれかに記載の造影剤注入プロトコル決定方法である。
【００２１】
この方法によれば、造影剤注入量が最低の注入プロトコルが抽出されるので、撮影対象（患者）に対する負担の少ない注入プロトコルを選択することができる。
請求項１１記載の発明は、撮影対象に対して一定のテスト注入速度で所定のテスト注入時間に渡って造影剤をテスト注入したときに関心領域を撮像して得られるテスト時間濃度曲線に基づいて本番検査のための造影剤注入プロトコルを演算する装置（２）であって、前記テスト時間濃度曲線に基づいて、注入速度を前記テスト注入速度と等しい一定値とし、前記テスト注入時間とは異なる所定の注入時間を設定した場合についての時間濃度曲線を推定して、推定時間濃度曲線を作成する推定時間濃度曲線作成手段（３７）と、この推定時間濃度曲線作成手段によって作成された前記推定時間濃度曲線において所要の濃度維持時間を確保できる最大濃度値を求める最大濃度値演算手段（３８）と、この最大濃度値演算手段によって演算された前記最大濃度値および所要濃度値に基づいて、前記濃度維持時間にわたって前記所要濃度値を維持できるように、前記所定の注入時間に対応した注入速度を定める注入速度演算手段（３９）と、前記所定の注入時間およびこれに対応して前記注入速度演算手段によって演算された注入速度に基づいて造影剤注入量を求め、この造影剤注入量およびこれに対応する前記注入速度の対を含む注入プロトコルを作成する注入プロトコル作成手段（４０）とを含む、造影剤注入プロトコル演算装置である。
【００２２】
この構成によれば、テスト注入の結果に基づいて、所定の濃度維持時間にわたって所要濃度値を確保できる適正な一相注入プロトコルを自動的に得ることができる。
むろん、この造影剤注入プロトコル演算装置の発明についても、前述の造影剤注入プロトコル決定方法の発明に関連して説明したのと同様な変形を施すことができる。
前述のような造影剤注入プロトコル演算装置は、所定のコンピュータプログラムをコンピュータに実行させることによって実現することができる。この場合、コンピュータプログラムは、コンピュータ上で実行されることにより、当該コンピュータを前記の各機能実現手段として機能させるものである。このようなコンピュータプログラムは、コンピュータによる読み取りが可能な記録媒体の形態で提供されてもよいし、通信回線を伝搬する信号の形態で提供されてもよい。
【００２３】
請求項１２記載の発明は、撮影対象に造影剤を注入する造影剤注入手段（３）と、前記テスト注入速度およびテスト注入時間に対応したテスト注入プロトコルに従って前記撮影対象に造影剤が注入されるように前記造影剤注入手段を制御するテスト注入制御手段（３２，Ｓ１，Ｓ２）と、請求項１１記載の造影剤注入プロトコル演算装置と、この造影剤注入プロトコル演算装置によって演算された注入プロトコルに従って前記撮影対象に造影剤が注入されるように前記造影剤注入手段を制御する本番注入制御手段（３２，Ｓ１３）とを含む、造影剤注入装置である。
【００２４】
この構成により、テスト注入の結果に基づき、適正な一相注入プロトコルが求められ、さらに、求められた注入プロトコルに従って、撮影対象に対して造影剤を注入することができる。
請求項１３記載の発明は、請求項１２記載の造影剤注入装置と、撮影対象を撮影して画像濃度に対応する信号を生成する撮影手段（１）と、テスト注入時に前記撮影対象の関心領域を撮影するように前記撮影手段を制御するテスト撮影制御手段（３１，Ｓ３）と、前記本番注入制御手段によって前記造影剤注入手段が制御されるときに、これに同期して前記撮影対象の関心領域を撮影するように前記撮影手段を制御する本番撮影制御手段（３１，Ｓ１３）とを含む、撮影システムである。
【００２５】
この構成により、テスト検査および本番検査を、造影剤注入装置および撮影装置を制御することにより、適切な一相注入プロトコルを適用して、自動的に行うことができる。
【発明を実施するための最良の形態】
【００２６】
以下では、この発明の実施の形態を、添付図面を参照して詳細に説明する。
１．撮影システムの構成
図１は、この発明の一実施形態に係る撮影システムの構成例を示すブロック図である。この撮影システムは、ＣＴスキャナ１と、操作コンソール２と、自動注入器３とを備えている。ＣＴスキャナ１は、被検体（医療現場における患者）の断面画像を撮影するための装置である。自動注入器３は、被検体に血管造影剤を自動注入するための装置である。操作コンソール２は、ＣＴスキャナ１および自動注入器３との間でデータを授受し、これらを制御する。より具体的には、操作コンソール２は、ＣＴスキャナ１に対して、スキャン範囲、スキャン開始時間その他の制御データを与えるとともに、ＣＴスキャナ１から撮影結果を表すデータを受信する。また、操作コンソール２は、自動注入器３に対して注入プロトコル（この実施形態では、造影剤注入量および注入速度を含む。）を設定するとともに、自動注入器３からの動作情報（注入開始、注入終了等）を受信する。
【００２７】
操作コンソール２は、コンピュータとしての基本構成を有している。より具体的には、二次元画像を表示する表示装置５と、キーボード、マウスその他の操作入力部６と、ＣＰＵを含むデータ処理部７と、各種のデータやプログラムを記録する記憶部８と、ＣＴスキャナ１および自動注入器３に対するデータ授受のためのインタフェース部９とを備えている。記憶部８は、ＲＡＭ、ＲＯＭ、ハードディスクその他の適切な記憶媒体で構成されており、少なくともその記憶領域の一部はデータの書き込みが可能なものである。
【００２８】
データ処理部７は、記憶部８に格納されたプログラムに従って動作し、これにより、データ処理部７および記憶部８などによって、この実施形態に係る造影剤注入プロトコル演算装置が構成されている。また、この造影剤注入プロトコル演算装置と自動注入器３とによって、この実施形態に係る造影剤注入装置が構成されている。
自動注入器３は、被検体の血管に造影剤を注入するための注入ヘッド１０と、この注入ヘッド１０の動作を制御するための制御ユニット１１と、操作コンソール２との間でデータを授受するためのインタフェース部１２とを備えている。制御ユニット１１は、操作コンソール２から受け渡される注入プロトコルに従って注入ヘッド１０を制御する。また、制御ユニット１１は、注入ヘッド１０の作動状態を監視し、造影剤の注入開始およびその完了を操作コンソール２へと通知する。
２．ＣＴ(Computed Tomography; コンピュータ断層撮影)
ＣＴ画像とは、被検体の断面画像であり、Ｘ線を用いて得られたスキャンデータに対して画像再構成処理を行うことによって得ることができる。
２−１．ＣＴの原理
図２に、ＣＴスキャナ１による走査方法の一例を示す。図２に示すのは、第３世代rotate/rotate方式による走査方式であり、被検体１５を包含する広がり角度３０〜４０°の扇状Ｘ線１６(fan beam方式)をＸ線管１７から発生させ、Ｘ線管１７に対向する位置に、検出器１８として、円弧上に５００〜９００個の検出チャンネルを持つＸｅ電離箱検出器が配置されている。Ｘ線管１７と検出器１８とは一体となって回転し、被検体１５を走査する（非特許文献２参照）。
【００２９】
Ｘ線管１７から出たＸ線１６は、被検体１５を透過し、検出器１８で検出されて、その強度が電気信号として計測される。Ｘ線１６は、被検体１５を透過している間に吸収され、その強度が弱くなる。Ｘ線１６の吸収の度合い(Ｘ線減弱係数) は、物質によって異なる。したがって、ＣＴにおいて、被検体１５は、Ｘ線減弱係数の空間的な分布ｆ(x,y,z)によって表されることになる。
【００３０】
Ｘ線管１７と検出器１８とは、ｚ＝一定の平面（被検体１５の体軸に直交する平面）上を回転し、この平面に含まれるあらゆる直線に沿って透過Ｘ線の強度を測定する。これにより、ｆ(x,y,z)｛ｚ＝一定｝という二次元の分布に関連するデータが収集される。これをスキャンといい、得られたデータをスキャンデータ（透過Ｘ線データ）という。
スキャンデータは、被検体１５の一次元投影分布の集まりである。このデータがＣＴスキャナ１から操作コンソール２に受け渡され、この操作コンソール２において画像再構成が行われる。これによって、Ｘ線減弱係数の分布を表す二次元画像が得られる。スライス面ｚを少しずつ変えて同じことをくり返せば、最終的に三次元の分布ｆ(x, y, z) が、多数の二次元画像の集合として表現されることになる。
【００３１】
三次元データを得る方式としては、ヘリカルＣＴが挙げられる。この方式は、Ｘ線１６を照射しながらＸ線管１７および検出器１８を連続回転すると同時に、被検体１５の乗る寝台（図示せず）を体軸方向に移動させ、２つの運動の合成によって被検体１５を螺旋状に走査した投影データを得るものである。
人体臓器の多くは水に近いＸ線減弱係数を持つので、水を基準にした尺度を使うと便利である。そこで、Ｘ線減弱係数の値に変換を施し、水が０、空気が−１０００になるようにする。この尺度をＣＴ値（単位：Hounsfield Unit＝H.U.）という（非特許文献３）。
２−２．ＣＴにおける画像再構成
ＣＴでの画像再構成とは、多方向から得られた被検体の一次元投影分布から、線源に関係する物理量の分布を二次元の画像として再生することである。Radonの原理に従って投影データから画像を再構成する数学的手法には、いくつかの方法があり、逆投影法、逐次近似法および解析的方法に大別される。
【００３２】
図３のように、透過Ｘ線データ（投影データ）から減弱係数の分布を推定する直観的な方法が逆投影法である。被検体を透過して得られた透過Ｘ線データの値を逆に投影し、すべての方向の逆投影を重ね合わせると、元の位置の減弱係数を推定できる。実際には、透過Ｘ線データに“ボケ” が含まれているため、そのまま使うことはできないが、基本的には、透過Ｘ線データの修正とこれを逆投影する処理とを数学的に行っている。
【００３３】
Ｘ線ＣＴでは、一般的に、フィルタ補正逆投影法またはコンボリューション逆投影法が用いられている。フィルタ補正逆投影法、コンボリューション逆投影法の具体的な手順を以下に示す。
図４に示すように、静止した直交座標系（ｘ，ｙ）に物体が置かれていて、この物体のＸ線減弱係数の（ｘ，ｙ）平面上における分布がｆ(x,y)により表されるものとする。（ｘ，ｙ）座標系に対して角度θだけ傾斜した新しい座標系（ｓ，ｔ）を定義し、そのｔ軸に沿って物体の外部からＸ線を照射して、物体に入射させるものとする。このとき、投影データｐ(s,θ)（透過Ｘ線データ）は、直線Ｌs,θ（ｔ軸に平行で座標系原点からｓ軸方向にｓだけ変位した直線）に沿った線積分を用いて、次式(1)で表される。
【００３４】
【数１】
2007143880.tif 000003
【００３５】
θを少しずつ変化させて得られる全方向の投影データｐ(s,θ) から分布ｆ(x,y)を復元することが画像再構成となる。ここで、画像ｆ(x,y)の二次元フーリエ変換は次式(2)のとおりであり、投影データｐ(s,θ) のｓに関する一次元フーリエ変換は次式(3)のとおりである。
【００３６】
【数２】
2007143880.tif 000004
【００３７】
ここで、（ｘ，ｙ）を角度θだけ回転させた座標が（ｓ，ｔ）なので、式(1)は、次式(4)のように書き換えられる。
【００３８】
【数３】
2007143880.tif 000005
【００３９】
式(3)に式(4)を代入し、（ｓ，ｔ）を（ｘ，ｙ）に変数変換すると、次式(5)が得られる。
【００４０】
【数４】
2007143880.tif 000006
【００４１】
式(2)と式(5)とを比較すると、次式(6)の関係を得る。
【００４２】
【数５】
2007143880.tif 000007
【００４３】
この式は、ｆ(x,y)のｓ軸への投影データｐ(s,θ) の一次元フーリエ変換Ｐ(w,θ)がｆ(x,y)の二次元フーリエ変換Ｆ(u,v) のｗ軸上の断面Ｆ(wcosθ,wsinθ) と一致することを示している。これを中央断面定理と呼び、この性質を用いて画像再構成を行う。
Ｆ(u,v) からｆ(x,y)を求めるフーリエ逆変換は、次式(7)のとおりである。
【００４４】
【数６】
2007143880.tif 000008
【００４５】
この式(7)において、変数変換u＝wcosθ、v＝wsinθを行うと、次式(8)が得られる。
【００４６】
【数７】
2007143880.tif 000009
【００４７】
この式(8)において、ｗに関する積分は、Ｐ(w,θ)と｜ｗ｜の積のフーリエ逆変換を表しているが、これはＰ(w,θ)のフーリエ逆変換ｐ(s,θ)と｜ｗ｜のフーリエ逆変換の畳み込みとして計算することができる。したがって、｜ｗ｜のフーリエ逆変換をｑ(s)とすると、次式(9)が得られる。q(s)は補正関数と呼ばれ、下記式(10)のとおりである。
【００４８】
【数８】
2007143880.tif 000010
【００４９】
式(9)は投影データｐ(s,θ)に補正関数ｑ(s)を畳み込みし、ｓ＝ｘcosθ+ｙsinθを満たすｓに対する畳み込みの値を角度方向について加算すればｆ(x,y)が再構成できることを示している。フィルタ補正処理を周波数領域の乗算で行う式(8)をフィルタ補正逆投影法、空間領域の畳み込みで行う式(9)をコンボリューション逆投影法と呼ぶ。
２−３．ＣＴＡngiography (ＣＴＡ；コンピュータ断層撮影血管造影法)
ＣＴＡとは、血管に造影剤を投与して行うＣＴであり、血行動態を調べたり、血管病変を診断する目的で行われる撮影法である。
【００５０】
通常のＸ線撮影では、臓器、血管、筋肉などの軟部組織は、Ｘ線吸収値の差が小さく、満足な診断を行うことが困難である。そこで、Ｘ線造影剤を血管内に直接注入することで、造影剤と人体組織のＸ線吸収値の差によって、臓器の明瞭な像を得ることができる。したがって、ＣＴＡ画像は、通常のＣＴ画像と比較して、臓器や血管のコントラストがより強調される。腹部血管造影（たとえば腹部大動脈の血管造影）では、造影剤として、軟部組織よりもＸ線吸収値の高い水溶性ヨード造影剤が用いられる（非特許文献２）。
【００５１】
一方、造影剤は時間経過に伴って血管内を運ばれていくので、造影剤を注入して得られる濃淡画像では、時間経過に伴って各部位の濃度値が変化する。ＣＴ造影像中の関心領域（注目領域）における濃度値の時間変化を表す曲線は「時間濃度曲線」と呼ばれる。
３．フーリエ変換を用いた時間濃度曲線の解析
３−１．対象データ
臨床において腹部大動脈のＣＴ造影を行う際は、まず、本番検査に先だって、テスト検査として、１０〜２０mlの造影剤を用いてテスト注入を行い、腹部大動脈のテスト時間濃度曲線(ＣＴ値の経時変化を表す時間濃度曲線) をデータとして得る。そのデータを元に、注目領域への造影剤到達時間を調べ、適切なＣＴスキャンの撮像開始時間を決定した上で、本番検査の撮像を行っている。
【００５２】
従来では、テスト時間濃度曲線は、本番検査のＣＴスキャンの開始時間を決定する目的のみで計測されている。これに対して、この実施形態では、テスト時間濃度曲線とテスト注入パラメータとの２つのデータを用いることで、各々の患者（被検体）ごとの時間濃度曲線の特徴を解析し、本番検査における時間濃度曲線のデータを推定する。
よって、この実施形態において、解析、推定を行う対象は、図５(a)(b)に示すような、ＣＴＡにおける造影剤注入設定の時系列データ（注入パラメータ）と、その注入を行って得られた濃度変化の時系列データ（時間濃度曲線）である。
【００５３】
図５(a)は自動注入器３（図１参照）に設定される造影剤注入速度の時系列データであり、これを以後、「造影剤注入関数」または単に「注入関数」（造影剤注入速度の時間変化を表す関数）と呼ぶ。テスト検査時の造影剤注入関数を区別するときには「テスト注入関数」と呼ぶ。
図５(b)は２−３節で述べたＣＴＡにより獲得した造影像の腹部大動脈部分の濃度値の時系列データ例であり、このような時系列データが「時間濃度曲線」を表すことになる。テスト検査時において、ＣＴＡは、腹部の特定断面に対して一定時間連続して行われ、その時間内で一定間隔ごとに造影像を得る。テスト検査によって得られる時間濃度曲線を区別するときには「テスト時間濃度曲線」と呼ぶことにする。
３−２．解析を行う上での基礎理論
従来手法（非特許文献１）では、図６に示すように、人体という複雑系をブラックボックスとみなして、入力と出力のみの分析、つまり、機器の造影剤注入関数injとそれによって得られる時間濃度曲線enhとの２つのデータのみを用いて分析している。
【００５４】
造影剤注入関数injとは、自動注入器３（造影剤注入機器）の設定値である「注入速度」と「総注入量」との２つの値から求められる時系列データを指し、これを系への入力と考える。入力から出力までの中間処理は一つの伝達関数で表し、以後これを「伝達関数Ｇ」と呼ぶ。
伝達関数Ｇは、あらかじめテストを行って得られたある被検体の造影剤注入関数（テスト注入関数）と、それによる時間濃度曲線（テスト時間濃度曲線）の値とを用いて導出される。この伝達関数Ｇを用いることで、当該被検体（患者）について、任意の造影剤注入関数に対応する時間濃度曲線を推定することができる。また、逆に、理想的な時間濃度曲線を得るための造影剤注入関数を推定することが可能である。
【００５５】
ここで、ブラックボックスとおいた部分を線形時不変システムであると仮定すると、次式(11)のように、動脈の時間濃度曲線enhは、系への入力となる造影剤注入関数injと、当該系を表すシステム関数ｇを畳み込み積分して求めることができ、フーリエ変換を用いて解析することが可能となる。
【００５６】
【数９】
2007143880.tif 000011
【００５７】
３−３．フーリエ変換を用いた解析と推定
式(11)は、時間領域での分析を行うための数学的関係を表している。ここで、動脈の時間濃度曲線enhは、造影剤注入関数injとシステム関数ｇとの２つの関数の畳み込みであるので、フーリエ変換を式(11) に適用すると次式(12)のように表すことができる。
【００５８】
【数１０】
2007143880.tif 000012
【００５９】
これは周波数領域での計算であり、Ｅnh、Ｉnj、Ｇは、式(11) の各関数をフーリエ変換して得られる周波数領域の関数である。
時間領域の計算も含めた伝達関数Ｇの導出過程を図７に示す。
伝達関数Ｇは、造影剤のテスト注入関数Ｉnjtestと、それによって得られるテスト時間濃度曲線Ｅnhtestとをあらかじめ測定しておくことで、次式(13)を用いて求められる。
【００６０】
【数１１】
2007143880.tif 000013
【００６１】
伝達関数Ｇが求まれば、造影剤注入関数Injから時間濃度曲線Enhを推定することは容易である。
時間濃度曲線の推定過程を図８に示す。
ある本番注入関数Ｉnjstdを入力として与えたときに得られる本番時間濃度曲線Ｅnhstdは、次式(14)で推定できる。
【００６２】
【数１２】
2007143880.tif 000014
【００６３】
式(14)によって得られたＥnhstdに逆フーリエ変換を適用すれば、時間領域の時間濃度曲線enhstdが得られる。
ここで、フーリエ変換の特性から、次式(15)によって、理想的な時間濃度曲線Enhidealを得るための理想的な注入関数Ｉnjidealを計算することが可能である。理想的な注入関数Ｉnjidealの推定過程を図９に示す。
【００６４】
【数１３】
2007143880.tif 000015
【００６５】
伝達関数Ｇの周波数スペクトルにおいて、０でない周波数成分に対しては、式(15)は有効である。しかし、一般に、伝達関数Ｇの高周波成分には０の値が含まれているため、対応するＥnhidealの高周波成分の計算を行わないようにすべきことに注意を要する。
腹部大動脈の理想的な時間濃度曲線を得るための理想的な造影剤注入関数の理論値injidealは、式(15)によって得られる注入関数Ｉnjidealを時間領域に逆フーリエ変換することによって計算することができる。
【００６６】
臨床で用いる場合は、理想的な造影剤注入関数の理論値に対してフィッティングを行うとともに、人体に注入を行っても問題のない注入パラメータを決定する必要がある。
非特許文献１に示されているとおり、理想的な造影剤注入関数に対して、現実的な造影剤注入関数をフィッティングするには、二相注入が必要になるが、実際の医療の現場では、二相注入は現実的でない場合が多い。したがって、実際には、式(15)によって理想的な造影剤注入関数を求めても、これを実現する現実的な手段がない。
４．造影剤注入パラメータの推定
図１０は、操作コンソール２に備えられたデータ処理部７（図１参照）の機能的な構成を説明するためのブロック図である。データ処理部７は、所定のプログラムを実行することによって、図１０に示された複数の機能処理部としての機能を実現するようになっている。これらの機能処理部には、ＣＴスキャナ１を制御するスキャン制御部３１と、自動注入器３を制御する注入制御部３２と、ＣＴスキャナ１によって収集された透過Ｘ線データに基づいて、被検体の断面画像を再構成する画像再構成部３３と、テスト注入によって得られたデータに基づいて本番注入時の造影剤注入プロトコルを演算する注入プロトコル演算部３４とを含む。
【００６７】
注入制御部３２は、記憶部８からテスト注入プロトコルを読み出し、これをインタフェース部９を介して自動注入器３の制御ユニット１１に受け渡すとともに、制御ユニット１１から自動注入器３の動作情報（造影剤の注入開始および注入終了に関する情報を含む。）を、インタフェース部９を介して受け取る。テスト注入プロトコルは、操作者が操作入力部６の操作によって予め入力することとして、この入力されたテスト注入プロトコルが記憶部８に格納されるようになっていてもよい。このテスト注入プロトコルは、この実施形態では、一相注入に対応する注入プロトコルであり、テスト造影剤注入量およびテスト注入速度の組によって表される。この場合、テスト注入関数は、一定のテスト注入速度がテスト注入時間（＝テスト造影剤注入量÷テスト注入速度）だけ継続する矩形の関数である。
【００６８】
スキャン制御部３１は、注入制御部３２と連携してＣＴスキャナ１を制御する。すなわち、注入制御部３２から、造影剤の注入開始および注入終了に関する情報を得て、ＣＴスキャナ１による撮影動作の開始および終了を制御する。この期間中、ＣＴスキャナ１は、被検体の関心領域を含む断面の撮影を継続し、透過Ｘ線データを出力する。この透過Ｘ線データがインタフェース部９を介して受信され、記憶部８へと蓄積される。
【００６９】
画像再構成部３３は、記憶部８に蓄積された透過Ｘ線データに基づいて、二次元画像を再構成する。これにより、時系列に従う複数枚の断面画像が再構成され、この複数枚の断面画像を表す画像データが記憶部８に書き込まれる。
注入プロトコル演算部３４は、テスト時間濃度曲線作成部３６と、推定時間濃度曲線作成部３７と、最大濃度値演算部３８と、注入速度演算部３９と、注入プロトコル作成部４０と、注入プロトコル選択部４１とを含み、本番検査の際に使用する注入プロトコルを演算する。
【００７０】
テスト時間濃度曲線作成部３６は、テスト注入時に画像再構成部３３によって作成された複数枚の断面画像中の関心領域に注目し、当該関心領域における画像濃度の時間変化を表すテスト時間濃度曲線（画像濃度の時系列データ）を作成する。このテスト時間濃度曲線は、記憶部８に格納される。
推定時間濃度曲線作成部３７は、記憶部８からテスト時間濃度曲線を読み出し、このテスト時間濃度曲線に基づいて、本番検査の際に想定される複数の造影剤注入時間（注入時間候補）に対応する推定時間濃度曲線を作成する。この推定時間濃度曲線は、造影剤の注入速度をテスト注入時と等しい値として、様々な注入時間で造影剤を被検体に注入したときにおいて、関心領域において観測されると予測される時間濃度曲線である。こうして作成された複数の推定時間濃度曲線は、記憶部８に格納される。
【００７１】
最大濃度値演算部３８は、記憶部８から、前記複数の推定時間濃度曲線を読み出し、各推定時間濃度曲線において、所定の濃度維持時間Ｔを満たすことができる最大濃度値ｈを求める。濃度維持時間Ｔは、操作入力部６からの設定操作によって、操作者によって設定することができるようになっていてもよい。この濃度維持時間Ｔは、ＣＴＡ撮影時において、関心領域の画像濃度が所定の所要濃度値Ｈ（濃度増強値）以上に保持されるべき時間である。各推定時間濃度曲線について求められた最大濃度値ｈは、記憶部８に格納されるようになっている。
【００７２】
注入速度演算部３９は、各推定時間濃度曲線について求められた最大濃度値ｈを記憶部８から読み出し、各最大濃度値ｈを前記所要濃度値Ｈと比較する。この比較結果に基づき、注入速度演算部３９は、各推定時間濃度曲線に対して（すなわち、各注入時間に対して）、前記所要濃度値Ｈを達成することができる造影剤注入速度ｖを演算する。この演算された複数の造影剤注入速度ｖは、それぞれの注入時間と対応付けられて、記憶部８に格納される。
【００７３】
注入プロトコル作成部４０は、記憶部８から前記演算された複数の注入速度ｖを読み出し、それらに対応する注入速度を乗じることによって、造影剤注入量（＝注入速度×注入時間）を求める。こうして、造影剤注入量および注入速度ｖの組からなる複数の造影剤注入プロトコルが作成され、記憶部８に格納される。
注入プロトコル選択部４１は、操作入力部６からの入力によって設定される造影剤注入条件に基づき、記憶部８に格納された前記複数の造影剤注入プロトコルの中から、適切な造影剤注入プロトコルを選択し、これを本番注入プロトコル（本番検査に使用すべき注入プロトコル）として記憶部８に格納する。この本番注入プロトコルが、注入制御部３２により、自動注入器３の制御ユニット１１へと受け渡される。
４−１．重ね合わせ、および、フーリエ変換を用いた濃度曲線の生成
人体という複雑系を線形時不変システムと仮定することにより、注入速度をテスト注入速度(ｖtest)と等しい値で一定値とし、注入時間をテスト注入時間ｔtestの自然数倍としたときの時間濃度曲線は、図１１に示すように、テスト時間濃度曲線の重ね合わせによって、時間領域で直接計算できる。より具体的には、テスト時間濃度曲線２１を時間軸方向にテスト注入時間ｔtestだけずらしたシフト時間濃度曲線２２を作成し、これらの曲線２１，２２を重ね合わせることによって、テスト注入時間の２倍の注入時間に対応した推定時間濃度曲線２３を作成することができる。同様にして、テスト注入時間の任意の自然数倍の注入時間に対応した推定時間濃度曲線の作成が可能である。
【００７４】
図１２は、推定時間濃度曲線作成部３７の構成例を説明するためのブロック図である。この推定時間濃度曲線作成部３７は、記憶部８から読み出したテスト時間濃度曲線をテスト注入時間ｔtestの自然数倍だけ時間軸方向にずらしたシフト時間濃度曲線を作成するシフト処理部４５と、重ね合わせ処理部４６とを備えている。重ね合わせ処理部４６は、シフト処理部４５によって作成されたシフト時間濃度曲線と、前記記憶部８から読み出したテスト時間濃度曲線とを重ね合わせる重ね合わせ処理を実行する。
【００７５】
シフト処理部４５は、テスト注入時間ｔtestの２倍の注入時間に対応する推定時間濃度曲線を作成するときには、テスト時間濃度曲線をテスト注入時間ｔtestだけ時間軸方向にずらした１つのシフト時間濃度曲線を作成する。また、テスト注入時間の３倍の注入時間に対応する推定時間濃度曲線を作成すべきときには、テスト時間濃度曲線をテスト注入時間ｔtestだけ時間軸方向にずらした第１のシフト時間濃度曲線と、テスト時間濃度曲線をテスト注入時間ｔtestの２倍の時間だけ時間軸方向にずらした第２のシフト時間濃度曲線とを作成する。テスト注入時間の４倍以上の自然数倍の注入時間に対応する推定時間濃度曲線を作成するときには、それに応じて、テスト注入時間ｔtestずつ時間軸上でずれた３個以上のシフト時間濃度曲線がシフト処理部４５によって作成される。こうしてシフト処理部４５で作成された１つまたは複数のシフト時間濃度曲線が、重ね合わせ処理部４６によってテスト時間濃度曲線に重ね合わせられることにより、テスト注入時間ｔtestの任意の自然数倍の注入時間に対応する推定時間濃度曲線を作成することができる。
【００７６】
時間領域の計算では、生成できる濃度曲線が、テスト注入時間の自然数倍の注入時間に関するものに制限されるという制約がある。
そこで、フーリエ変換を利用することにより、図７および図８の関係から、任意の注入関数に対する推定時間濃度曲線が得られるので、注入速度一定の条件下で注入時間を任意の値に変化させたときの推定時間濃度曲線を生成できる。
【００７７】
図１３は、推定時間濃度曲線作成部３７の他の構成例を説明するためのブロック図である。この推定時間濃度曲線作成部３７は、テスト注入関数injtestをフーリエ変換する第１フーリエ変換部５１と、テスト時間濃度曲線enhtestをフーリエ変換する第２フーリエ変換部５２と、第２フーリエ変換部５２の演算結果（周波数領域のテスト時間濃度曲線Ｅnhinj）に対して高周波成分を除去するためのフィルタ処理を行うフィルタ処理部５４と、第１フーリエ変換部５１およびフィルタ処理部５４の演算結果（周波数領域のテスト注入関数Ｉnjtestおよびテスト時間濃度曲線Ｅnhtest（フィルタ処理後のもの））に基づいて、前記式(13)に従って伝達関数Ｇを演算する伝達関数演算部５３と、この伝達関数演算部５３よって求められた伝達関数Ｇに対して逆フーリエ変換を施し、インパルス応答を求める逆フーリエ変換部５５と、このインパルス応答を時間軸上でシフトさせるシフト処理部５７と、シフト処理部５７によってシフトされたインパルス応答を逆フーリエ変換部５５によって求められたインパルス応答に重ね合わせることによって、任意の注入時間の推定時間濃度曲線を演算する重ね合わせ処理部５６と含んでいる。
【００７８】
伝達関数Ｇの逆フーリエ変換はシステム関数gとなり、このシステム関数gは被検体に対応したブラックボックス（図６参照）のインパルス応答（微小単位時間（たとえば０．２５秒）の注入に対応する被検体の応答）を表す。そこで、このインパルス応答を、推定時間濃度曲線を求めたい注入時間に応じて時間軸上でシフトさせて重ね合わせることにより、任意の注入時間に対応した推定時間濃度曲線を求めることができる。
【００７９】
フィルタ処理部５４は、たとえば、テスト時間濃度曲線のスペクトルにおいて、最大振幅に対して所定比率（たとえば１％）に満たない周波数成分を除去することによって、高周波成分を除去する。これにより、伝達関数Ｇの高周波成分に０値が含まれることを回避できる。
図１４(a)は周波数領域のテスト時間濃度曲線から高周波成分を除去せずに伝達関数Ｇを求めた場合のインパルス応答（システム関数ｇ）を示し、図１４(b)は前述の高周波除去フィルタ処理を行ったテスト時間濃度曲線に基づいて伝達関数Ｇを求めた場合のインパルス応答（システム関数ｇ）を示す。これらの比較から、フィルタ処理によって、インパルス応答の高周波成分を除去できることが判る。
【００８０】
図１５(a)〜(h)に、図１４(b)のインパルス応答を重ね合わせて種々の注入時間（１秒、２秒、３秒、４秒、５秒、１０秒、１５秒、２０秒）の推定時間濃度曲線を生成した例を示す。このように、テスト注入時間の自然数倍の注入時間に限らず、任意の注入時間の推定時間濃度曲線を生成できる。
図１６に、前述の図１２の構成によって作成した推定時間濃度曲線の一例（実線）と、図１３の構成によって作成した推定時間濃度曲線の一例（破線）とを併せて示す。これらは、同じテスト時間濃度曲線から生成したものである。
【００８１】
図１６の２つの曲線の比較から、フィルタ処理によって一部の高周波成分を除外した場合でも、重ね合わせで生成した時間濃度曲線とほぼ同じ推定時間濃度曲線を得ることができ、フィルタ処理の影響は非常に小さいことが分かる。したがって、フーリエ変換を利用することにより、注入時間に関する制約を避けながら、任意の注入時間についての時間濃度曲線を良好に推定できる。逆に、テスト注入時間の自然数倍の注入時間についてのプロトコル演算で十分な場合には、フーリエ変換を含まない簡単な処理で時間濃度曲線を良好に推定できる。
【００８２】
図１７は、推定時間濃度曲線作成部３７のさらに他の構成例を説明するためのブロック図である。この図１７において、前述の図１３示された各部と同様の機能部分には図１３の場合と同一の参照符号を付して示す。この推定時間濃度曲線作成部３７は、テスト注入関数injtestをフーリエ変換する第１フーリエ変換部５１と、テスト時間濃度曲線enhtestをフーリエ変換する第２フーリエ変換部５２と、第２フーリエ変換部５２の演算結果（周波数領域のテスト時間濃度曲線Ｅnhtest）に対して高周波成分を除去するためのフィルタ処理を行うフィルタ処理部５４と、第１フーリエ変換部５１およびフィルタ処理部５４の演算結果（周波数領域のテスト注入関数Ｉnjtestおよびテスト時間濃度曲線Ｅnhtest（フィルタ処理後のもの））に基づいて、前記式(13)に従って伝達関数Ｇを演算する伝達関数演算部５３と、テスト注入時間とは異なる複数の注入時間に対応した注入関数（設定注入関数：注入速度の時系列データ）のフーリエ変換を作成する第３フーリエ変換部５８と、この第３フーリエ変換部５８の演算結果（周波数領域の注入関数）とフィルタ処理後の伝達関数Ｇとを用い、前記式(14)に従って、周波数領域の推定時間濃度曲線Ｅnhを求める推定時間濃度曲線演算部５９と、この推定時間濃度曲線演算部５９の演算結果に対して逆フーリエ変換を施すことによって、時間領域の推定時間濃度曲線enhを求める逆フーリエ変換部６０とを備えている。
【００８３】
この構成によって、任意の注入時間に対応した注入関数によって得られるべき時間濃度曲線（推定時間濃度曲線）を演算することができる。ただし、この実施形態では、造影剤の注入を一相注入によって行うこととしているので、第３フーリエ変換部５８に与えられる注入関数は、注入速度がテスト注入速度と等しく、注入時間のみがテスト注入時間と異なる矩形の関数である。したがって、一般には、図１３の構成の場合の方が、処理が簡単になる。
【００８４】
フィルタ処理部５４の働きによって、テスト時間濃度曲線Ｅnhtestの高周波成分を予め除去しているため、伝達関数Ｇの高周波成分に０値を含むことを回避できる。
４−２．造影剤注入パラメータの推定
図１８は、推定時間濃度曲線作成部３７が前述の図１２の構成を有する場合に、本番検査のための注入プロトコルを求める処理手順を説明するためのフローチャートである。操作者は、予め、操作入力部６を操作することによって、テスト注入プロトコルを入力する（ステップＳ１）。このテスト注入プロトコルは、例えば、テスト注入速度ｖtest、テスト注入時間ｔtestおよびテスト注入量ａtestを含み、ａ＝ｖtestｔ（ａは注入量、ｖは注入速度、ｔは時間）の関係となるように、すなわち、テスト注入時間に渡って一定のテスト注入速度ｖtestで造影剤が注入されるように定められる。この場合、テスト注入関数injtestは、一定のテスト注入速度ｖtestがテスト注入時間ｔtestにわたって継続する矩形の関数で表される。
【００８５】
注入制御部３２は、テスト注入プロトコルを自動注入器３の制御ユニット１１に受け渡す。これにより、自動注入器３は、テスト注入プロトコルに従って、造影剤を被検体に対してテスト注入する（ステップＳ２）。このときの自動注入器３の動作を表す動作情報は、制御ユニット１１から注入制御部３２に通知され、この動作情報に基づき、スキャン制御部３１が、ＣＴスキャナ１を制御する。こうして、造影剤のテスト注入に並行して、被検体がプレスキャン（テスト検査のためのスキャン）される（ステップＳ３）。
【００８６】
ＣＴスキャナ１によって収集された透過Ｘ線データは、操作コンソール２の記憶部８に格納され、この透過Ｘ線データに基づき、画像再構成部３３によって、時系列に従う複数枚のＣＴＡ造影像を表す画像データが作成される（ステップＳ４）。さらに、この複数枚の画像データに基づき、テスト時間濃度曲線作成部３６により、テスト時間濃度曲線が作成される（ステップＳ５）。
【００８７】
次に、操作者は、本番検査の際に獲得したい濃度増強値ＨＵ［ΔＨＵ］と維持時間Ｔ（秒）とを操作入力部６から入力する（ステップＳ６）。
推定時間濃度曲線作成部３７は、前述のように、テスト時間濃度曲線の重ね合わせ処理によって、注入速度をテスト注入速度と等しい値で一定とし、注入時間をテスト注入時間の自然数倍とした場合に得られるべき時間濃度曲線の推定値（推定時間濃度曲線）を作成して、記憶部８に格納する（ステップＳ７）。推定時間濃度曲線作成部３７は、異なる複数の注入時間に関して同様の処理を行い、複数の推定時間濃度曲線を作成する。
【００８８】
作成される推定時間濃度曲線の例は、図１９(a)，１９(b)，１９(c)に示されている。これらの各図の上段には、注入関数inj（注入速度ｖの時間変化を表す関数）が示されており、各図の下段に推定時間濃度曲線enhが示されている。
図１９(a)は、テスト注入関数inj1と、それに対応するテスト時間濃度曲線enh1を示している。この例では、テスト注入ｖtest速度を３ml/秒として、テスト注入時間ｔtestを５秒とした例が示されている。本番検査において濃度増強値Ｈを維持すべき維持時間Ｔを確保できる最大濃度値はｈ1である。
【００８９】
図１９(b)は、注入速度ｖをテスト注入速度ｖtestと等しい一定値３ml/秒とし、注入時間ｔ2をテスト注入時間ｔtestの２倍の１０秒とした注入関数inj2の場合を示す。図１９(a)のテスト時間濃度曲線ehn1に対して、これをテスト注入時間ｔtestである５秒だけ時間軸方向にシフトさせて得られるシフト時間濃度曲線を重ね合わせることにより、図１９(b)の推定時間濃度曲線enh2が得られる。この推定時間濃度曲線enh2において維持時間Ｔを確保できる最大濃度値はｈ2である。
【００９０】
図１９(c)は、注入速度ｖをテスト注入速度ｖtestと等しい一定値３ml/秒とし、注入時間ｔ6をテスト注入時間ｔtestの６倍の３０秒とした注入関数inj6の場合を示す。図１９(a)のテスト時間濃度曲線ehn1に対して、これを５秒（テスト注入時間ｔtestの１倍）、１０秒（テスト注入時間ｔtestの２倍）、１５秒（テスト注入時間ｔtestの３倍）、２０秒（テスト注入時間ｔtestの４倍）、および２５秒（テスト注入時間ｔtestの５倍）だけ時間軸方向にそれぞれシフトさせて得られる五つのシフト時間濃度曲線を重ね合わせることにより、図１９(c)の推定時間濃度曲線enh6が得られる。この推定時間濃度曲線enh6において維持時間Ｔを確保できる最大濃度値はｈ6である。
【００９１】
推定時間濃度曲線作成部３７は、同様にして、テスト注入時間の３倍〜５倍の注入時間ｔ3〜ｔ5に対応した推定時間濃度曲線enh3〜enh5を作成する。実際には、推定時間濃度曲線ehni（ｉは自然数であり、注入時間がテスト注入時間ｔtestの何倍であるかを表すインデックスである。）は、推定時間濃度曲線ehni-1に対して、テスト時間濃度曲線enhiを時間軸方向に（ｉ−１）・ｔtestだけずらして得られる一つのシフト時間濃度曲線を重ね合わせることによって生成することができる。
【００９２】
再び図１８を参照する。
最大濃度値演算部３８は、作成された複数の（この実施形態では６個の）推定時間濃度曲線enh1〜enh6（ただし、enh1はテスト時間濃度曲線）において、前記維持時間Ｔを確保できる最大濃度値ｈ1，ｈ2，…，ｈ6を求める（ステップＳ８）。維持時間Ｔは、ＣＴスキャナ１における本番検査時の撮像に要する時間である。
【００９３】
次に、注入速度演算部３９は、各推定時間濃度曲線enh1〜enh6に対して得られた最大濃度値ｈ1〜ｈ6に対する濃度増強値Ｈの比Ｈ／ｈi（ただし、ｉ＝１，２，３，……）を求め、これにテスト注入速度ｖtestを乗じることによって、各注入時間ｔiに対応する注入速度ｖi（＝ｖtest・Ｈ／ｈi）を求める（ステップＳ９）。このようにして注入速度ｖiが求まるのは、注入関数inj(t)と濃度曲線enh(t)とが線形関係にあると仮定しているからである。
【００９４】
さらに、注入プロトコル作成部４０は、注入速度ｖiと、これに対応する注入時間ｔiとを乗じることによって、造影剤注入量ａi（ただし、ａ1＝ａtest（テスト注入時の造影剤注入量）を演算する。こうして、造影剤注入量ａiおよび注入速度ｖiの組で表される６通りの注入プロトコルが得られ（ステップＳ１０）、この６通りの注入プロトコルを表すデータが記憶部８に格納される。この注入プロトコルは、注入時間ｔiに渡って一定の注入速度ｖiが保持される一相注入に対応した矩形の注入関数を表す。
【００９５】
図２０は、注入速度演算部３９による処理の内容を説明するための図であり、図２０(a)(b)(c)は、それぞれ、図１９(a)(b)(c)に対応している。注入速度ｖiを、テスト注入速度ｖtestのＨ／ｈi倍に設定することにより、各注入時間ｔiに渡る一相注入によって、維持時間Ｔにわたって所要の濃度増強値Ｈを確保できる。
なお、図２０の各図の上段には、注入速度ｖiをテスト注入速度ｖtestのＨ／ｈi倍した場合における推定時間濃度曲線が表されており、各図の下段には、それぞれ対応する注入関数が表されている。各注入関数は、造影剤注入量およびおよび注入速度の組によって表される注入プロトコルとして、前述のとおり、記憶部８に格納されることになる。図２０の各図の上段に示す推定時間濃度曲線は、これを作成して記憶部８に格納しておくことにより、本番検査時において、ＣＴスキャナ１による撮像タイミングの制御に用いることができる。各推定時間濃度曲線において、維持時間Ｔの開始時間が撮影を開始する時間となり、維持時間Ｔの終了時間以前に撮影を終了することになる。
【００９６】
こうして６通りの注入プロトコルが求まると、図１８に示すように、注入プロトコル選択部４１は、操作者が操作入力部６からの入力によって指定する注入条件（ステップＳ１１）に適合する１つの注入プロトコルを本番検査のための注入プロトコル（本番注入プロトコル）として選択する（ステップＳ１２）。注入条件としては、たとえば造影剤注入量が最少であること、造影剤の注入速度が所定値以下であること、造影剤の注入速度が最低であること、その他、注目したいパラメータに応じて、適切な条件を設定することができる。
【００９７】
このようにして本番注入プロトコルが定まると、この本番注入プロトコルが、注入制御部３２から自動注入器３の制御ユニット１１に受け渡される。それとともに、スキャン制御部３１により、ＣＴスキャナ１の制御が行われ、自動注入器３による造影剤の注入と同期して、本番検査のための撮像がＣＴスキャナ１によって実行される（ステップＳ１３）。このときの撮像開始および終了のタイミングは、本番注入プロトコルに対応する推定時間濃度曲線に基づいて定められる。
【００９８】
図２１は、推定時間濃度曲線作成部３７が前述の図１３の構成を有する場合における処理手順を説明するためのフローチャートである。この図２１において、前述の図１８と同様の処理または操作が行われる各ステップには、図１８の場合と同一の参照符号を付して示す。この図２１に示されている処理では、図１８におけるステップＳ７の処理の代わりに、次のような処理が行われる。
【００９９】
すなわち、第１フーリエ変換部５１によって、テスト注入関数injtestのフーリエ変換Ｉnjtestが求められる（ステップＳ２１）。また、第２フーリエ変換部５２によって、テスト時間濃度曲線ehhtestのフーリエ変換Ｅhhtestが求められ（ステップＳ２２）、これに対してフィルタ処理部５４によるフィルタ処理（高周波成分除去処理）が行われる（ステップＳ２３）。こうして求められた周波数領域のテスト注入関数Ｉnjtestおよびテスト時間濃度曲線Ｅhhtestに基づき、伝達関数演算部５３により、被検体の伝達関数Ｇが、前記式(13)に従って求められる（ステップＳ２４）。さらに、この伝達関数Ｇに対して、逆フーリエ変換部５５による逆フーリエ変換が行われることによって、被検体のシステム関数ｇが求められる（ステップＳ２５）。このシステム関数ｇは、前述のとおりインパルス応答を表すものである。
【０１００】
次いで、重ね合わせ処理部５６およびシフト処理部５７の働きによって、前記インパルス応答を複数の異なる注入時間ｔiに渡って時間領域で重ね合わせることにより、それぞれの注入時間ｔiに対応する複数の推定時間濃度曲線enhiが求められ、記憶部８に格納される（ステップＳ２６）。この場合の注入時間ｔiは、テスト注入時間の自然数倍である必要はない。
【０１０１】
この後の処理は、図１８を参照して説明したステップＳ８〜Ｓ１３の処理と同様である。
図２２は、推定時間濃度曲線作成部３７が前述の図１７の構成を有する場合における処理手順を説明するためのフローチャートである。この図２２において、前述の図２１と同様の処理または操作が行われる各ステップには、図２１の場合と同一の参照符号を付して示す。この図２２に示されている処理では、図２１におけるステップＳ２５，Ｓ２６の処理の代わりに、次のような処理が行われる。
【０１０２】
すなわち、注入速度ｖiがテスト注入速度ｖtestに等しく、注入時間ｔiがテスト注入時間ｔtestとは異なる複数の注入関数inji（一相注入を表すもの）が作成される（ステップＳ３１）。そして、個々の注入関数injiに対して、第３フーリエ変換部５８により、注入関数injiのフーリエ変換Ｉnjiが求められる（ステップＳ３２）。さらに、伝達関数Ｇと周波数領域の注入関数Ｉnjiとを用いて、前記式(14)に従って、推定時間濃度曲線演算部５９により、周波数領域の推定時間濃度曲線Ｅnhiが求められる（ステップＳ３３）。これが逆フーリエ変換部５５によって逆フーリエ変換されることにより、時間領域の推定時間濃度曲線enhiが得られ、記憶部８に格納される（ステップＳ３４）。ステップＳ３３およびＳ３４の処理を複数の注入関数injiに対して繰り返し実行することによって、複数の注入時間ｔiにそれぞれ対応した複数の推定時間濃度曲線enhiが得られる。注入時間ｔiは、テスト注入時間ｔtestの自然数倍である必要はない。
【０１０３】
この後の処理は、図１８を参照して説明したステップＳ８〜Ｓ１３の処理と同様である。
５．実験
５−１．ファントムの概要
ヒトの腹部大動脈の血流を摸倣するファントムを用いて実験を行った。このファントムは、その時間濃度曲線が、ヒトの時間濃度曲線と類似していることが実験的に確認されているものである。
【０１０４】
まず、このファントムにおいて想定されている、患者内に注入された造影剤の挙動を以下に示す。
ステップ１．自動注入器を用いて患者の左の肘静脈に造影剤を注入する。造影剤が注入完了した時点で肘静脈に溜まった造影剤を後押しするため、引き続き自動注入器によって生理食塩水の注入を行う。
ステップ２．注入された造影剤は、左肘静脈の血管を通り、心臓の右心房へ流入する。
ステップ３．造影剤は右心房から右心室、肺動脈、肺、肺静脈、そして心臓の左心房、左心室へと移動する。
ステップ４．造影剤は心臓の左心室から腹部大動脈へと流入する。
ステップ５．造影剤が腹部大動脈を通過する。ＣＴＡが行われるのはこの部位（スキャン領域）であり、定点で一定時間の造影が行われて、濃度値の変動を得る。
ステップ６．造影剤が腹部大動脈を通過した後、全身をめぐって再び右心房へと流入する。
ステップ７．ステップ３へ戻る。
【０１０５】
以上の造影剤の循環を摸倣したファントムの概略図を図２３に示す。ここで、
ｖ0：自動注入器に設定されている造影剤注入速度(ml/秒)
Ｃv：ファントム内の循環速度(ml/秒)
ＴC：チューブの容積(ml)
である。
【０１０６】
このファントムは、腹部大動脈の時間濃度変化に関係する人体の各部分をモデル化したものである。ステップ１、２の注入箇所から心臓までの肘静脈はチューブに、ステップ３の右心房から左心室までの心肺系はポンプと水槽に、ステップ４、５の腹部大動脈はスキャン領域Ｐ1に、ステップ６における「全身」はスキャン領域Ｐ1から水槽までの間に、それぞれ対応している。
【０１０７】
ファントムは水で満たされており、ポンプにより、循環速度Ｃv[ml/秒]で水を循環させている。また、造影剤注入後、同一速度で２０mlの整理食塩水で後押しすることで、チューブ内の造影剤をファントムに注入している。
５−２．ファントム実験データ
ファントム実験では、ＣＴ装置としてPhilips Medical Systems（米国、Ohio 州、Cleveland 市）製Brilliance40（４０列マルチスライスＣＴ）を用い、自動注入器（Dual Shot、根本杏林堂、東京）から造影剤（Iomeron-350、ヨード含量300mgI/ml、エーザイ、東京）を注入した。注入開始後１秒ごとに６０秒まで撮像し、腹部大動脈部（図２３のスキャン領域Ｐ1)の濃度値を計測した。
【０１０８】
以下の計算では、１秒刻みの実測値のサンプリング間隔を線形補間により０．２５秒としたデータを用いて推定した。
下記表１のテスト注入パラメータに対するテスト時間濃度曲線を図２４に示す。
【０１０９】
【表１】
2007143880.tif 000016
【０１１０】
下記表２に示す種々の注入パラメータで造影剤を注入して画像濃度を実測した観測値を図２５(a)〜(e)において曲線Ｌ１で示す。各図の曲線Ｌ２，Ｌ３は、前述の図１３に示す手法によって作成した推定時間濃度曲線を示す。ただし、曲線Ｌ２はテスト時間濃度曲線（図２４参照）の６０秒までの全観測点のデータを用いた推定結果を示し、曲線Ｌ３はテスト時間濃度曲線（図２４参照）の２０秒までのデータを用い、２０秒以降の濃度値を「０」とした場合の推定結果を示す。
【０１１１】
曲線Ｌ１と曲線Ｌ２，Ｌ３との比較から、本実施形態の手法により、時間濃度曲線の実測値に適合する推定時間濃度曲線が得られることがわかる。
【０１１２】
【表２】
2007143880.tif 000017
【０１１３】
下記表３には、図２４のテスト時間濃度曲線を用いて、注入時間１０秒、２０秒および３０秒について、２００ＨＵの濃度増強を１０秒間維持するのに必要な注入速度等を推定した結果を示す。この結果から、テスト時間濃度曲線の全観測点のデータを用いた場合と、テスト時間濃度曲線の０〜２０秒までのデータを用いた場合とでほぼ同じ結果が得られることがわかる。したがって、テスト時間濃度曲線において濃度変化が有意な変化を示す期間のデータを抽出して用いることによって、演算負荷を軽減しつつ、妥当な注入パラメータを決定できることがわかる。また、テスト時間濃度曲線において濃度変化が有意な変化を示す期間のデータが収集できれば十分であるので、テスト検査時におけるスキャン時間を短く設定することが可能であり、これにより、被曝の影響を軽減しつつ、妥当な注入パラメータを決定できる。
【０１１４】
【表３】
2007143880.tif 000018
【０１１５】
６．まとめ
以上のようにこの実施形態によれば、注入速度をテスト注入速度と等しい一定値とし、注入時間をテスト注入時間とは異ならせた複数の注入関数のそれぞれについて、テスト時間濃度曲線に基づいて時間濃度曲線が推定される。そして、その推定時間濃度曲線において所要の濃度維持時間Ｔを確保できる最大濃度値ｈを求め、この最大濃度値ｈと所要濃度増強値Ｈとの比に基づいて注入速度ｖが求められる。このようにして、フィッティングを経験的に行ったりすることなく、濃度増強値Ｈを濃度維持時間Ｔに渡って維持できる複数の適正な注入プロトコルを得ることができる。しかも、この注入プロトコルは一相注入に対応するものであるので、実際の医療の現場において、現実に活用することができる。
７．変形例
前述の実施形態では、注入プロトコルに対する条件設定を行うことによって、本番検査に使用する注入プロトコルが自動的に選択されるようにしているが、操作者が、複数の注入プロトコルから本番検査に使用する注入プロトコルを手動で選択するようにしてもよい。
【０１１６】
また、前述の実施形態では、Ｘ線を利用した断層撮影に本発明を適用した例を示したが、この発明は、ＭＲＩ（磁気共鳴イメージング）や超音波画像診断などの他の撮影法のシステムに対しても適用することができる。
さらに、前述の実施形態では、ＣＴスキャナ１、操作コンソール２および自動注入器３が連携するシステムの例について説明したが、ＣＴスキャナ１や自動注入器３から切り離されたコンピュータによって、テスト注入プロトコルおよびテスト時間濃度曲線に基づいて本番注入プロトコルを演算する造影剤注入プロトコル演算装置を構成するようにしてもよいし、このようなコンピュータと自動注入器とを接続して、注入プロトコル演算機能付きの造影剤自動注入装置を構成してもよい。
【０１１７】
その他、特許請求の範囲に記載された事項の範囲で種々の設計変更を施すことが可能である。
【図面の簡単な説明】
【０１１８】
【図１】この発明の一実施形態に係る撮影システムの構成例を示すブロック図である。
【図２】ＣＴスキャナによる走査方法の一例を示す概念図である。
【図３】透過Ｘ線データ（投影データ）から減弱係数の分布を推定する方法の一つである逆投影法を説明するための概念図である。
【図４】フィルタ補正逆投影法、コンボリューション逆投影法を説明するための概念図である。
【図５】造影剤注入関数および時間濃度曲線の例を示す図である。
【図６】人体を表すシステム関数ｇの考え方を示す概念図である。
【図７】人体を表す伝達関数Ｇの導出過程を説明するための図である。
【図８】時間濃度曲線の推定過程を説明するための図である。
【図９】理想的な注入関数の推定過程を説明するための図である。
【図１０】操作コンソールに備えられたデータ処理部の機能的な構成を説明するためのブロック図である。
【図１１】テスト時間濃度曲線の重ね合わせによって推定時間濃度曲線を求める処理を説明するための図である。
【図１２】推定時間濃度曲線作成部の構成例を説明するためのブロック図である。
【図１３】推定時間濃度曲線作成部の他の構成例を説明するためのブロック図である。
【図１４】システム関数（インパルス応答）を示す図である。
【図１５】インパルス応答から作成した推定時間濃度曲線の例を示す図である。
【図１６】図１２および図１３の構成によって作成した推定時間濃度曲線の一例を示す。
【図１７】推定時間濃度曲線作成部のさらに他の構成例を説明するためのブロック図である。
【図１８】推定時間濃度曲線作成部が図１２の構成を有する場合に、本番検査のための注入プロトコルを求める処理手順を説明するためのフローチャートである。
【図１９】推定時間濃度曲線から所要の維持時間を確保できる最大濃度値を求める処理を説明するための図である。
【図２０】所要の濃度増強値と前記最大濃度値とから注入速度を求める処理を説明するための図である。
【図２１】推定時間濃度曲線作成部が図１３の構成を有する場合における処理手順を説明するためのフローチャートである。
【図２２】推定時間濃度曲線作成部が図１７の構成を有する場合における処理手順を説明するためのフローチャートである。
【図２３】人体における造影剤の循環を模倣したファントムの構成を示す概略図である。
【図２４】表１のテスト注入パラメータに対するテスト時間濃度曲線を示す
【図２５】実測により求めた時間濃度曲線および推定時間濃度曲線の例を示す図である。
【符号の説明】
【０１１９】
１    ＣＴスキャナ
２    操作コンソール
３    自動注入器
５    表示装置
６    操作入力部
７    データ処理部
８    記憶部
９    インタフェース部
１０    注入ヘッド
１１    制御ユニット
１２    インタフェース部
１５    被検体
１６    Ｘ線
１７    Ｘ線管
１８    検出器
２１    テスト時間濃度曲線
２２    シフト時間濃度曲線
２３    推定時間濃度曲線
３１    スキャン制御部
３２    注入制御部
３３    画像再構成部
３４    注入プロトコル演算部
３６    テスト時間濃度曲線作成部
３７    推定時間濃度曲線作成部
３８    最大濃度値演算部
３９    注入速度演算部
４０    注入プロトコル作成部
４１    注入プロトコル選択部
４５    シフト処理部
４６    重ね合わせ処理部
５１    第１フーリエ変換部
５２    第２フーリエ変換部
５３    伝達関数演算部
５４    フィルタ処理部
５５    逆フーリエ変換部
５６    重ね合わせ処理部
５７    シフト処理部
５８    第３フーリエ変換部
５９    推定時間濃度曲線演算部
６０    逆フーリエ変換部
【図１】
2007143880.tif 000019
【図２】
2007143880.tif 000020
【図３】
2007143880.tif 000021
【図４】
2007143880.tif 000022
【図５】
2007143880.tif 000023
【図６】
2007143880.tif 000024
【図７】
2007143880.tif 000025
【図８】
2007143880.tif 000026
【図９】
2007143880.tif 000027
【図１０】
2007143880.tif 000028
【図１１】
2007143880.tif 000029
【図１２】
2007143880.tif 000030
【図１３】
2007143880.tif 000031
【図１４】
2007143880.tif 000032
【図１５】
2007143880.tif 000033
【図１６】
2007143880.tif 000034
【図１７】
2007143880.tif 000035
【図１８】
2007143880.tif 000036
【図１９】
2007143880.tif 000037
【図２０】
2007143880.tif 000038
【図２１】
2007143880.tif 000039
【図２２】
2007143880.tif 000040
【図２３】
2007143880.tif 000041
【図２４】
2007143880.tif 000042
【図２５】
2007143880.tif 000043