JP6899666B2 - 造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 - Google Patents
造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 Download PDFInfo
- Publication number
- JP6899666B2 JP6899666B2 JP2017036115A JP2017036115A JP6899666B2 JP 6899666 B2 JP6899666 B2 JP 6899666B2 JP 2017036115 A JP2017036115 A JP 2017036115A JP 2017036115 A JP2017036115 A JP 2017036115A JP 6899666 B2 JP6899666 B2 JP 6899666B2
- Authority
- JP
- Japan
- Prior art keywords
- scan
- subject
- contrast
- control unit
- timing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000002591 computed tomography Methods 0.000 title claims description 157
- 238000000034 method Methods 0.000 title claims description 88
- 239000002872 contrast media Substances 0.000 claims description 80
- 210000004204 blood vessel Anatomy 0.000 claims description 49
- 210000001631 vena cava inferior Anatomy 0.000 claims description 49
- 210000000056 organ Anatomy 0.000 claims description 40
- 230000000877 morphologic effect Effects 0.000 claims description 36
- 201000010099 disease Diseases 0.000 claims description 35
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 35
- 238000003384 imaging method Methods 0.000 claims description 31
- 210000002376 aorta thoracic Anatomy 0.000 claims description 25
- 238000012937 correction Methods 0.000 claims description 16
- 210000000702 aorta abdominal Anatomy 0.000 claims description 11
- 210000000709 aorta Anatomy 0.000 claims description 9
- 210000000038 chest Anatomy 0.000 claims description 9
- 238000007619 statistical method Methods 0.000 claims description 8
- 210000004705 lumbosacral region Anatomy 0.000 claims description 6
- 238000010801 machine learning Methods 0.000 claims description 6
- 210000000115 thoracic cavity Anatomy 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 5
- 230000000004 hemodynamic effect Effects 0.000 description 21
- 210000002796 renal vein Anatomy 0.000 description 17
- 210000003240 portal vein Anatomy 0.000 description 14
- 210000001758 mesenteric vein Anatomy 0.000 description 13
- 238000012360 testing method Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 8
- 238000000611 regression analysis Methods 0.000 description 8
- 210000003462 vein Anatomy 0.000 description 6
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000000747 cardiac effect Effects 0.000 description 4
- 230000004087 circulation Effects 0.000 description 4
- 238000010968 computed tomography angiography Methods 0.000 description 4
- 210000005246 left atrium Anatomy 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 230000004660 morphological change Effects 0.000 description 3
- 230000002792 vascular Effects 0.000 description 3
- 241000702481 Chloris striate mosaic virus Species 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000001990 intravenous administration Methods 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 206010002091 Anaesthesia Diseases 0.000 description 1
- 206010007559 Cardiac failure congestive Diseases 0.000 description 1
- 208000003037 Diastolic Heart Failure Diseases 0.000 description 1
- 206010019280 Heart failures Diseases 0.000 description 1
- 206010037423 Pulmonary oedema Diseases 0.000 description 1
- 206010038669 Respiratory arrest Diseases 0.000 description 1
- 238000000692 Student's t-test Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000037005 anaesthesia Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000001269 cardiogenic effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 229940039231 contrast media Drugs 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000000326 densiometry Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012048 forced swim test Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000001631 haemodialysis Methods 0.000 description 1
- 230000002440 hepatic effect Effects 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 210000001363 mesenteric artery superior Anatomy 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000004738 parenchymal cell Anatomy 0.000 description 1
- 208000005333 pulmonary edema Diseases 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 210000001321 subclavian vein Anatomy 0.000 description 1
- 238000012353 t test Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明は、造影ＣＴスキャン装置、造影ＣＴスキャン装置の作動方法、及び造影ＣＴスキャン像を得るための方法に関する。
被験者に造影剤を投与してＣＴスキャンを行うダイナミック造影ＣＴ（Ｄｙｎａｍｉｃ Ｃｏｎｔｒａｓｔ−Ｅｎｈａｎｃｅｄ Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ：ＤＣＥ−ＣＴ）は、通常のＣＴスキャンと比較して、標的血管、標的臓器または標的疾患について、アーティファクトがより少ない、高いコントラストの造影画像を得ることができる。ＤＣＥ−ＣＴでは、造影剤を被験者に投与し、造影剤が標的血管、標的臓器または標的疾患に到達したときにＣＴスキャンを行うことにより、標的血管、標的臓器または標的疾患を精度よく観察することができる。この際、特に小さな病変を描出するために、適切な造影時相（Ｃｏｎｔｒａｓｔ Ｅｎｈａｎｃｅｍｅｎｔ Ｐｈａｓｅ：ＣＥＰ）で病変対臓器コントラスト（Ｌｅｇｉｏｎ−ｔｏ−Ｏｒｇａｎ Ｃｏｎｔｒａｓｔ：ＬＯＣ）を最大化することが重要である。
また、微細な主要栄養血管の描出や３ＤＣＴ血管造影画像を生成するために、動脈対静脈コントラスト（Ａｒｔｅｒｉａｌ−ｔｏ−Ｖｅｎｏｕｓ Ｃｏｎｔｒａｓｔ：ＡＶＣ）を最大化することが重要である。ＡＶＣは、造影剤の投与量及び投与速度が固定されている場合、スキャンタイミングに依存する。
しかし、ＤＣＥ−ＣＴにおける造影時相は、造影剤の投与について、同じ投与プロトコル及び投与後に同じスキャンタイミングでＣＴスキャンを行ったとしても、被験者の血行動態によって異なる。例えば、被験者が心不全を有する場合、造影剤の循環速度は低くなる。また、被験者の体格が大きいほど、造影剤が標的血管、標的臓器または標的疾患に到達する時間が長くなる。
様々な投与プロトコルにおいて、標的臓器または病変について、適切なＣＥＰでＣＴスキャンを行うために推奨されるスキャンタイミングがこれまでに報告されている。しかし、造影剤の投与後、標的血管、標的臓器または標的疾患のＣＥＰに合わせた所定のスキャンタイミングでＣＴスキャンを行う固定スキャンタイミング（Ｆｉｘｅｄ Ｓｃａｎ Ｔｉｍｉｎｇ：ＦＳＴ）法（非特許文献１）では、スキャンタイミングは経験的なものであり、前述のように、造影剤が投与されてから標的血管、標的臓器または標的疾患に到達する時間が被験者の血行動態によって異なるため、ＤＣＥ−ＣＴにおいて、ＣＥＰが常に最適となり、ＬＯＣまたはＡＶＣが最大となるようにＣＴスキャンを行うことができない。例えば最適なタイミングに対して数秒から数十秒程度異なるタイミングでＣＴスキャンを行っても、標的血管、標的臓器または標的疾患以外の血管等の造影効果が大きくなり、最適な造影ＣＴスキャン像を得ることができない。そこで、被験者の血行動態を考慮してＤＣＥ−ＣＴにおけるＣＥＰを最適化し、またはＡＶＣを最大化するために、ボーラストラッキング法や試験投与法などのいくつかの方法がこれまで提案されている（非特許文献２）。
ボーラストラッキング法は、造影剤を被験者に投与した後、低線量で所定の間隔でＣＴスキャンを行うことにより、投与された造影剤を追跡し、標的血管、標的臓器または標的疾患に到達したタイミングで、必要な解像度が得られる線量でＣＴスキャンを行う技術である。従って、ボーラストラッキング法によれば、ある一定以上の造影効果で標的血管、標的臓器または標的疾患のＤＣＥ−ＣＴスキャンを行うことが可能である。しかし、この方法では、放射線被曝量の増加、複雑なスキャンプロトコル及び放射線技師の技量が結果に影響するなどのいくつかの欠点を有する。さらに、動脈内の造影剤の到達タイミングは、ボーラストラッキング法によって観察可能であるが、動脈の真のピークまたは標的血管における実質臓器のコントラストの最大化を本方法では予測することができないため、本方法では、標的臓器における最大ＬＯＣを得ることができない。
試験投与法は、被験者に投与する造影剤の一部をまず投与して、標的血管、標的臓器または標的疾患における造影剤の時間濃度曲線（Ｔｉｍｅ Ｄｅｎｓｉｔｙ Ｃｕｒｖｅ：ＴＤＣ）を試験スキャンで測定し、実際のＤＣＥ−ＣＴスキャンにおける標的血管、標的臓器または標的疾患のＴＤＣの定量予測を行う技術である。しかし、試験投与された造影剤は標的血管、標的臓器または標的疾患に残留したり、拡散したりするため、ＤＣＥ−ＣＴスキャン像に悪影響を及ぼす。また、試験スキャンにより、被験者の放射線被曝量が増加する。さらに、造影剤の注入量は実際のＣＴスキャンと試験スキャンとでは全く異なるため、造影剤の血中濃度が異なり、試験スキャンで予測されたＴＤＣは常に実際のＣＴスキャンにおけるＴＤＣと合致するとは限らない。
Bae KT, Heiken JP, Brink JA., "Aortic and hepatic contrast medium enhancement at CT Part I, Prediction with a computer model", Radiology, 1998; 207(3): 647-55
Henzler T, Meyer M, Reichert M, et al., "Dual-energy CT angiography of the lungs: Comparison of test bolus and bolus tracking techniques for the determination of scan delay", European Journal of Radiology, 2012; 81: 132-138
Hallett RL1, Fleischmann D., "Tools of the trade for CTA: MDCT scanners and contrast medium injection protocols", Tech. Vasc. Interv. Radiol. 2006; 9(4): 134-42
Aurigemma GP, Gaasch WH, "Clinical practice: Diastolic heart failure", N. Engl. J. Med. 2004; 351: 1097-1105
Zhu Y, Xu H, Zhu X, et al., "Which can predict left ventricular size and systolic function: cardiothoracic ratio or transverse cardiac diameter", J. Xray Sci. Technol., 2015; 23(5): 557-65
Echt M, Duweling J, Ganer OH, et al., "Effective compliance of the total vascular bed and the intrathoracic compartment derived from changes in central venous pressure induced by volume changes in man", Circulation Research, 1974; 34: 61-68
Cheriex EC, Leunissen KM, Janssen JH, et al., "Echography of the inferior vena cava is a simple and reliable tool for estimation of ‘dry weight’ in haemodialysis patients", Nephrol. Dial. Transplant., 1989; 4: 563-568
Kent A, Bahner DP, Boulger CT, et al., "Sonographic evaluation of intravascular volume status in the surgical intensive care unit: A prospective comparison of subcla-vian vein and inferior vena cava collapsibility index", J. Surf. Res., 2013; 184: 561-566
Figueras J, Weil MH, "Blood volume prior to and following treatment of acute cardiogenic pulmonary edema", Circulation 1978; 57: 349-355
Duvekot JJ, Cheriex EC, Tan WD, et al., "Measurement of anterior-posterior diameter of inferior vena cava by ultrasonography: A new non-invasive method to assess acute changes in vascular filling state", Cardiovasc. Res., 1994; 28: 1269-1272
Iijima T1, Ueyama H, Oi Y, et al., "Determination of the standard value of circulating blood volume during anesthesia using pulse dye-densitometry: a multicenter study in Japan", J. Anesth., 2005; 19(3): 193-8
Fleischmann D, Hittmair K, "Mathematical analysis of arterial enhancement and optimization of bolus geometry for CT angiography using the discrete Fourier transform", J. Comput. Assist. Tomogr., 1999; 23(3): 474-84
Fleischmann D, Rubin GD, Bankier AA, et al., "Improved uniformity of aortic enhancement with customized contrast medium injection protocols at CT angiography", Radiology, 2000; 214(2): 363-71
Yamashita Y, Komohara Y, Takahashi M, et al., "Abdominal CT: evaluation of optimal doses of intravenous contrast material − a prospective randomized study", Radiology, 2000; 216: 718-723
Bae KT, "Peak contrast enhancement in CT and MR angiography: when does it occur and why: Pharmacokinetic study in a porcine model", Radiology, 2003; 227: 809-816
Yamada A, et al., "8110 − Quantitative prediction of contrast enhancement phases in abdominal dynamic contrast-enhanced CT using morphologic prediction factors on pre-contrast", CT. ECR, 2017
追加的な放射線の被曝及び造影剤の追加的投与なく被験者の血行動態を推測することは、適切なＣＥＰでより良好なＤＣＥ−ＣＴ画像を得るために重要な課題である。
その一方、患者の血行動態は、心機能障害の指標としての心拡大や、患者の循環血液量の指標としての下大静脈の寸法など、画像化における形態学的特徴に反映されうる。従って、ＤＣＥ−ＣＴにおけるＣＥＰは、患者の血行動態の指標として、造影前ＣＴで観察されたこれらの形態学的変化、いわゆる形態学的予測因子（Ｍｏｒｐｈｏｌｏｇｉｃ Ｐｒｅｄｉｃｔｉｏｎ Ｆａｃｔｏｒ：ＭＰＦ）を利用して定量的に予測可能であると考えられる。
本発明は、被験者の血行動態及び造影時相についての定量的指標として造影前ＣＴから得られた形態学的予測因子（ＭＰＦ）及びダイナミック時相指数（ＤＰＩ）を用いた、ＤＣＥ−ＣＴにおけるＣＥＰの定量的予測の可能性を示し、血管のＣＴスキャン像の取得において、固定スキャンタイミング（ＦＳＴ）法と比較して、修正スキャンタイミング（ＭＳＴ）法がＤＣＥ−ＣＴの最良のスキャンタイミングの定量的かつ汎用的予測が可能であることを示し、改善されたコントラストを有するＣＴスキャン像を取得することができる造影ＣＴスキャン装置、造影ＣＴスキャン装置の作動方法、及び造影ＣＴスキャン像を得るための方法を提供することである。
本発明に係るＣＴスキャン装置は、ＣＴスキャン撮像ユニットと、制御ユニットと、を含み、制御ユニットが、被験者の基本的情報及び標的血管、標的臓器または標的疾患のダイナミック時相指数を受け取るように構成され、ＣＴスキャン撮像ユニットが、被験者の第１のＣＴスキャンを実行するように構成され、制御ユニットが、ＣＴスキャン撮像ユニットから第１のＣＴスキャン像を受信するように構成され、制御ユニットが、第１のＣＴスキャン像から、被験者の少なくとも１つの形態学的予測因子を測定するように構成され、制御ユニットが、被験者の基本情報、ダイナミック時相指数、及び形態学的予測因子に基づいて、造影剤投与の開始から第２回目以降のＣＴスキャンを行うまでの時間を、修正スキャンタイミングとして決定するように構成され、ＣＴスキャン撮像ユニットが、造影剤の投与の開始から修正スキャンタイミングが経過したときに第２回目以降のＣＴスキャンを実行するように構成され、制御ユニットが、ＣＴスキャン撮像ユニットから第２回目以降のＣＴスキャン像を受信するように構成される。
本発明のＣＴスキャン装置が、造影剤投与ユニットをさらに含み、造影剤投与ユニットが、制御ユニットの命令により、被験者に造影剤を投与するように構成される。
本発明のＣＴスキャン装置において、被験者の基本的情報が、性別、年齢、身長及び体重のうち少なくとも１つを含む。
本発明のＣＴスキャン装置において、形態学的予測因子が、下行大動脈の面積及び下大静脈の短軸径のうち少なくとも１つを含む。
本発明のＣＴスキャン装置において、形態学的予測因子が、腰椎長、胸郭横断径、心胸郭比、大動脈長、上行及び下行大動脈の面積、腹部大動脈の面積、下大静脈の長軸径並びに下大静脈の短軸径のうち少なくとも１つを含む。
本発明のＣＴスキャン装置において、修正スキャンタイミングが、ダイナミック時相指数、被験者の基本情報及び形態学的予測因子の少なくとも１つを説明変数として、統計解析及びもしくは非線形機械学習により決定された方程式または非線形モデルから得られる。
本発明のＣＴスキャン装置において、修正スキャンタイミングを求めるための方程式が、ダイナミック時相指数、被験者の基本情報及び形態学的予測因子の少なくとも１つからなる線形方程式または非線形モデルである。
本発明のＣＴスキャン装置において、修正スキャンタイミングを求めるための方程式が、対象とする被験者の母集団に応じて統計解析を行う被験者のサンプルを選択して係数を最適化したものである。
本発明のＣＴスキャン装置において、修正スキャンタイミングが下記の方程式、
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
本発明に係る、ＣＴスキャン撮像ユニットと、造影剤投与ユニットと、制御ユニットと、を含むＣＴスキャン装置の作動方法は、制御ユニットが被験者の基本的情報及び標的血管、標的臓器または標的疾患のダイナミック時相指数を受け取る段階と、制御ユニットが、ＣＴスキャン撮像ユニットに、被験者の第１のＣＴスキャンを実行させる命令を送信する段階と、制御ユニットが、ＣＴスキャン撮像ユニットから第１のＣＴスキャン像を受信する段階と、制御ユニットが、第１のＣＴスキャン像から、被験者の少なくとも１つの形態学的予測因子を測定する段階と、制御ユニットが、被験者の基本情報、ダイナミック時相指数、及び形態学的予測因子に基づいて、造影剤投与の開始から第２回目以降のＣＴスキャンを行うまでの時間を、修正スキャンタイミングとして決定する段階と、制御ユニットが、造影剤投与ユニットに、造影剤の投与を開始させる命令を送信する段階と、制御ユニットが、ＣＴスキャン撮像ユニットに、造影剤投与の開始から、修正スキャンタイミングが経過したときに第２回目以降のＣＴスキャンを実行させる命令を送信する段階と、制御ユニットが、ＣＴスキャン撮像ユニットから第２回目以降のＣＴスキャン像を受信する段階と、を含む。
本発明に係る方法において、被験者の基本的情報が、性別、年齢、身長及び体重のうち少なくとも１つを含む。
本発明に係る方法において、形態学的予測因子が、下行大動脈の面積及び下大静脈の短軸径のうち少なくとも１つを含む。
本発明に係る方法において、形態学的予測因子が、腰椎長、胸郭横断径、心胸郭比、大動脈長、上行及び下行大動脈の面積、腹部大動脈の面積、下大静脈の長軸径並びに下大静脈の短軸径のうち少なくとも１つを含む。
本発明に係る方法において、修正スキャンタイミングが、ダイナミック時相指数、被験者の基本的情報及び形態学的予測因子の少なくとも１つを説明変数として、統計解析及びもしくは非線形機械学習により決定された方程式または非線形モデルから得られる。
本発明に係る方法において、修正スキャンタイミングを求めるための方程式が、ダイナミック時相指数、被験者の基本的情報及び形態学的予測因子の少なくとも１つからなる線形方程式または非線形モデルである。
本発明に係る方法において、修正スキャンタイミングを求めるための方程式が、対象とする被験者の母集団に応じて統計解析を行う被験者のサンプルを選択して係数を最適化したものである。
本発明に係る方法において、修正スキャンタイミングが下記の方程式、
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
本発明に係る造影ＣＴスキャン像を得るための方法は、造影剤投与前に、被験者に対して第１のＣＴスキャンを行う段階と、第１のＣＴスキャンにおいて得られたＣＴスキャン像から、被験者の少なくとも１つの形態学的予測因子を測定する段階と、形態学的予測因子、被験者の基本的情報及び標的血管、標的臓器または標的疾患について予め決定されたダイナミック時相指数に基づいて、造影剤投与の開始から第２回目以降のＣＴスキャンを行うまでの時間を、修正スキャンタイミングとして決定する段階と、造影剤を被験者に投与し、造影剤の投与開始から修正スキャンタイミングが経過したときに第２回目以降のＣＴスキャンを行う段階と、を含む。
本発明に係る方法において、被験者の基本的情報が、性別、年齢、身長及び体重のうち少なくとも１つを含む。
本発明に係る方法において、形態学的予測因子が、下行大動脈の面積及び下大静脈の短軸径のうち少なくとも１つを含む。
本発明に係る方法において、形態学的予測因子が、腰椎長、胸郭横断径、心胸郭比、大動脈長、上行及び下行大動脈の面積、腹部大動脈の面積、下大静脈の長軸径並びに下大静脈の短軸径のうち少なくとも１つを含む。
本発明に係る方法において、修正スキャンタイミングが、ダイナミック時相指数、被験者の基本的情報及び形態学的予測因子の少なくとも１つを説明変数として、統計解析及びもしくは非線形機械学習により決定された方程式または非線形モデルから得られる。
本発明に係る方法において、修正スキャンタイミングを求めるための方程式が、ダイナミック時相指数、被験者の基本的情報及び形態学的予測因子の少なくとも１つからなる線形方程式または非線形モデルである。
本発明に係る方法において、修正スキャンタイミングを求めるための方程式が、対象とする被験者の母集団に応じて統計解析を行う被験者のサンプルを選択して係数を最適化したものである。
修正スキャンタイミングが下記の方程式、
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）から得られる。
本発明によれば、造影前ＣＴスキャンで観察されたＭＰＦから得られた修正スキャンタイミング（ＭＴＦ）を適用して造影ＣＴスキャンを行うことにより、適切なＣＥＰで改善されたコントラストを有する良好なＤＣＥ−ＣＴ画像を取得することができる造影ＣＴスキャン装置、造影ＣＴスキャン装置の作動方法、及び造影ＣＴスキャン像を得るための方法を提供することができる。
本発明に係る、ダイナミック造影ＣＴ（ＤＣＥ−ＣＴ）による造影ＣＴスキャン像を得るための方法について説明する。
まず、造影ＣＴスキャンを行う被験者の形態学的予測因子（ＭＰＦ）を測定するために、被験者に造影剤を投与せずに、ＣＴスキャン撮像ユニットを用いて造影前ＣＴスキャンを行う。図１は、造影前ＣＴスキャンを行って得られた像を、ＭＰＦの測定例とともに示している。ＭＰＦは、図１に白色の矢印及び円で視覚的に示している。ＭＰＦの例として、図１ａは、腰椎長を、図１ｂは胸郭横断径を、図１ｃは心胸郭比（Ｃ／Ｔ）を、図１ｄは下行大動脈の面積を、図１ｅは下大静脈の長軸径を、図１ｆは下大静脈の短軸径を、それぞれ示す。しかし、ＭＰＦはこれらに限定されるものではなく、被験者の血行動態に影響を及ぼしうる体の様々な器官の寸法、例えば大静脈長、上行大動脈の面積及び腹部大動脈の面積を含みうる。
下行大動脈の面積は、例えば心最大横断面と同じ断面で測定されうる。下大静脈の長軸径及び短軸径は、例えば左右腎静脈流入部と肝部下大静脈との間で測定されうる。上行大動脈及び腹部大動脈の面積は、例えばその最大横断面において測定されうる。
これらのＭＰＦは、造影前ＣＴスキャン像をビューアに表示し、放射線技師または医師が測定箇所を決定して測定を行うことができる。また、ＣＴスキャン撮像ユニットに接続されたコンピュータなどの制御ユニットが造影前ＣＴスキャン像を受信して、画像処理技術を用いて自動的にＭＰＦを測定することも可能である。
次に、標的血管、標的臓器または標的疾患のダイナミック時相指数（ＤＰＩ）を選択する。ＤＰＩは、造影ＣＴスキャンにおいて、標的血管、標的臓器または標的疾患の造影時相についての定量的指標である。被験者に造影剤を投与すると、造影剤は、投与位置から血流によって循環される。従って、標的血管、標的臓器または標的疾患への造影剤の到達時間は、投与位置に近いほど短くなる。そのため、ＤＰＩは血管または臓器ごとに固有の値を有し、例えば右腎静脈を標的血管として造影ＣＴスキャンを行う場合、右腎静脈に対応する、あらかじめ求められたＤＰＩを選択する。
図２に、ＤＰＩ算出に用いる標的血管のＣＴスキャン像の例を示す。図２のＣＴスキャン像において、標的血管の関心領域を矢印で指した白い円形のマーカーで示した。図２ａは左心房を、図２ｂは腹部大動脈を、図２ｃは右腎静脈を、図２ｄは上腸間膜静脈を、図２ｅは下大静脈を、図２ｆは門脈を、それぞれ示している。
各スキャンタイミングにおける標的血管の造影効果を、左心房（ＣＬＡ（ｔ））、腹部大動脈（ＣＡＡＯ（ｔ））、右腎静脈（ＣＲＲＶ（ｔ））、上腸間膜静脈（ＣＳＭＶ（ｔ））、下大静脈（ＣＩＶＣ（ｔ））及び門脈（ＣＰＶ（ｔ））とすると、各スキャンタイミングｔにおけるＤＰＩは、以下の方程式を用いて、ＤＣＥ−ＣＴにおける標的血管内の造影効果から計算することができる。
ＤＰＩ（ｔ）＝｛ＣＬＡ（ｔ）＋ＣＡＡＯ（ｔ）＋ＣＲＲＶ（ｔ）＋ＣＳＭＶ（ｔ）＋ＣＩＶＣ（ｔ）｝／Ｃｍａｘ（ｔ）
標的血管、標的臓器または標的疾患におけるＤＰＩは、各標的血管等におけるコントラストが最大となる経過時間ｔにおけるＤＰＩ（ｔ）の値であり、標的血管、標的臓器、標的疾患ごとに任意に設定されうる。推奨される標的血管、標的臓器、標的疾患ごとのＤＰＩは、事前に十分な数の母集団から得た多時相造影ＣＴデータにより別途決定される。
次に、被験者のＭＰＦ、性別、年齢、身長及び体重等の基本的情報、標的血管、標的臓器または標的疾患のＤＰＩの少なくとも１つに基づいて、修正スキャンタイミング（ＭＳＴ）を決定する。ＭＳＴは、造影剤の投与開始から、造影ＣＴスキャンを行うまでの時間であり、被験者のＭＰＦ、性別、年齢、身長及び体重等の基本的情報並びに標的血管、標的臓器または標的疾患のＤＰＩ（説明変数）が考慮された値である。従って、ＭＳＴを用いることで、血行動態が被験者ごとに異なっても、常に最適なタイミングで標的血管、標的臓器または標的疾患の造影ＣＴスキャンを行うことができる。
上述の説明変数からＭＳＴを決定する方法は、種々のものが考えられるが、例えば、上述の説明変数から回帰分析などの統計解析及び人工ニューラルネットワークなどの非線形機械学習により決定された方程式とすることができ、例えば、回帰分析を行って得られた説明変数の少なくとも１つを含む線形方程式、例えば線形１次方程式とすることができる。また、ＭＳＴを決定する方程式は、人工ニューラルネットワークなどで表される非線形モデルから得ることもできる。ＭＳＴの決定に用いるＭＰＦのうち、特に下行大動脈の面積及び下大静脈の短軸径の寄与が大きいため、ＭＳＴを決定する方程式は、下行大動脈の面積及び下大静脈の短軸径のみを説明変数とした線形１次方程式としてもよい。好適には、ＭＳＴは、以下の方程式に説明変数を代入して得ることができる。
ＭＳＴ［秒］＝
１２．４×ＤＰＩ
−４５．６
＋１．７４×（性別（男＝１、女＝０））
−０．１０×（腰椎長［ｍｍ］）
＋０．１９×（胸郭横断径［ｍｍ］）
＋０．３６×（心胸郭比）
＋０．０１×（下行大動脈の面積［ｍｍ２］）
−０．００１×（下大静脈の長軸径［ｍｍ］）
＋０．２６×（下大静脈の短軸径［ｍｍ］）
ＭＳＴ［秒］＝
１２．４×ＤＰＩ
−４５．６
＋１．７４×（性別（男＝１、女＝０））
−０．１０×（腰椎長［ｍｍ］）
＋０．１９×（胸郭横断径［ｍｍ］）
＋０．３６×（心胸郭比）
＋０．０１×（下行大動脈の面積［ｍｍ２］）
−０．００１×（下大静脈の長軸径［ｍｍ］）
＋０．２６×（下大静脈の短軸径［ｍｍ］）
ＭＳＴは、例えば予め被験者の基本的情報及び標的血管、標的臓器または標的疾患のＤＰＩが入力された制御ユニットが計算して求めることができる。また、対象とする被験者の母集団（国、地域、人種）に応じて回帰分析などの統計解析を行うサンプルを選択することにより、方程式の係数を最適化することができ、より精度の高い修正スキャンタイミングを得ることが可能である。
次に、被験者に造影剤を投与する。造影剤は、ＣＴスキャンにおいて大きな造影効果を得ることができるように、Ｘ線に対する吸収係数の大きな元素、例えばヨウ素を含むことができる。造影剤は、被験者に予め取り付けられたカニューレ等を介して、例えば被験者の正中肘静脈に経静脈投与することができる。精度よく造影剤を被験者に投与することができるように、例えばポンプを有する造影剤投与ユニットを用いて、１秒間に３ｍＬの投与速度で１００ｍＬの造影剤を投与することができる。造影剤の投与開始のタイミング及び投与速度は、造影剤投与ユニットに接続された制御ユニットを用いて、適切に制御することができる。
次に、造影剤の投与開始から前述のＭＳＴが経過したときに、ＣＴスキャン撮像ユニットによって、被験者の造影ＣＴスキャンを行う。これにより、標的血管、標的臓器または標的疾患に最も適したタイミングで、造影ＣＴスキャンを行うことができる。造影ＣＴスキャンを複数回行う場合、造影ＣＴスキャンを行うたびに造影前ＣＴスキャンを事前に行ってＭＳＴを求めるのではなく、１回の造影前ＣＴスキャンで得られたＭＳＴを繰り返し適用して複数回の造影ＣＴスキャンを行うことも可能である。造影ＣＴスキャンは、作業者が開始命令を入力してもよい。しかし、より精度よく造影ＣＴスキャン像を得ようとする場合、造影剤投与ユニットによる造影剤の投与開始からＭＳＴが経過したときに、制御ユニットがＣＴスキャン撮像ユニットに造影ＣＴスキャンを行う命令を送信してもよい。
このような方法により、被験者の血行動態の違いに関わらず、最適なタイミングで標的血管、標的臓器または標的疾患の造影ＣＴスキャン像を得ることができる。
このような造影ＣＴスキャン像を得るための造影ＣＴスキャン装置は、例えば、制御ユニットと、ＣＴスキャン撮像ユニットとを含むものとすることができる。制御ユニットは、ＣＴスキャン撮像ユニットに、被験者の造影前ＣＴスキャンを実行する命令を送信し、ＣＴスキャン撮像ユニットからのＣＴスキャン像を受信する。次いで、制御ユニットは、ＣＴスキャン像からＭＰＦを測定する。ＭＰＦの測定は、作業者が制御ユニットを制御して行ってもよいし、画像処理技術を用いて、制御ユニットが自動的にＣＴスキャン像からＭＰＦを測定してもよい。次いで、制御ユニットは、あらかじめ入力された被験者の基本的情報及び標的血管、標的臓器または標的疾患のＤＰＩ、並びにＭＰＦから、例えば前述の方程式に基づいてＭＳＴを決定する。次いで、制御ユニットは、被験者への造影剤の投与開始からＭＳＴが経過したときに、ＣＴスキャン撮像ユニットに造影ＣＴスキャンを実行する命令を送信し、得られたＣＴスキャン像を受信する。
造影ＣＴスキャン装置はさらに、制御ユニットの制御により、所定の投与速度で被験者に造影剤を投与することができる造影剤投与ユニットを含むものとしてもよい。この場合、制御ユニットが造影剤投与ユニットに造影剤の投与開始命令を送信し、ＭＳＴが経過したときにＣＴスキャン撮像ユニットに造影ＣＴスキャンを実行する命令を送信することができるため、精度よく最適なタイミングで造影ＣＴスキャンを行うことができる。
次に、被験者のＭＰＦ及び基本的情報並びにＤＰＩからＭＳＴを決定した方法について説明する。
被験者のＭＰＦ及び基本的情報並びにＤＰＩとＭＳＴとの間の相関関係を決定するために、外科手術前の術前検査として、上腹部多時相造影ＣＴを行った６０人の被験者（平均年齢６７．３歳、男女比＝３１：２９）について、遡及的研究を行った。
６４列ＭＤＣＴ（米国ミシガン州ＧＥ社のＬｉｇｈｔ Ｓｐｅｅｄ ＶＣＴ）を用いて、呼吸停止下において、正中肘静脈に留置された２２Ｇ静脈カニューレを介して、造影剤を経静脈的投与（１秒間に３ｍＬの一定速度で約３３秒、全投与量１００ｍＬ）を行った。造影剤投与前、及び経静脈的投与開始から２２、２８、３４、４０、４６、５２及び５８秒経過後に、被験者についてＣＴスキャンを行った。ＣＴデータは、厚さ２．５ｍｍのスライスで、軸位断方向に画像再構成した。
造影剤は、特に制限はされないが、例えば、ヨード造影剤である、イオヘキソール、イオパミドール、イオキシラン、イオプロミド、イオトロラン、イオジキサノール、イオペントール、イオベルソール、イオメプロール、イオビトリドール、アミドトリゾエート、メトリゾエート、イオキサグレート、オヨビイオキシタラメートからなる群から選択されうる。造影剤の単位容量は、特に制限はされないが、例えば、約８０ｍＬから約１５０ｍＬの容量でありうる。造影剤の単位ヨード濃度は、特に制限はされないが、例えば、例えば、２４０ｍｇＩから３７０ｍｇＩでありうる。
ＣＴスキャン装置は、前述の装置のほかに、特に制限はされないが、例えば、東芝メディカルシステム株式会社製のＡｑｕｉｌｉｏｎ Ｏｎｅ（登録商標）、ＧＥヘルスケア・ジャパン株式会社製のＲｅｖｏｌｕｔｉｏｎ ＣＴ、シーメンスヘルスケア株式会社製のＳＯＭＡＴＯＮ Ｄｅｆｉｎｉｔｉｏｎ、株式会社フィリップス エレクトロニクス ジャパン製のＩｎｇｅｎｕｉｔｙ、株式会社日立メディコ製のＳＣＥＮＡＲＩＡ等が使用されうる。
造影剤投与装置は、特に制限されないが、例えば、株式会社根本杏林堂製のＤＵＡＬ ＳＨＯＴ ＧＸ ７、バイエル薬品株式会社製のＳｔｅｌｌａｎｔ等が使用されうる。
各スキャンタイミングにおける標的血管の造影効果（左心房（ＣＬＡ（ｔ））、腹部大動脈（ＣＡＡＯ（ｔ））、右腎静脈（ＣＲＲＶ（ｔ））、上腸間膜静脈（ＣＳＭＶ（ｔ））、下大静脈（ＣＩＶＣ（ｔ））及び門脈（ＣＰＶ（ｔ）））を、ＤＣＥ−ＣＴ画像から、ＤＩＣＯＭビューア上で、２人の放射線医師によって設定された関心領域から得た。図２のＣＴスキャン像において、標的血管の関心領域を白い円形のマーカーで示した。図２ａは左心房を、図２ｂは腹部大動脈を、図２ｃは右腎静脈を、図２ｄは上腸間膜静脈を、図２ｅは下大静脈を、図２ｆは門脈を、それぞれ示している。
各スキャンタイミングｔにおけるＤＰＩを、以下の方程式を用いて、ＤＣＥ−ＣＴにおける標的血管内の造影効果から計算した。
ＤＰＩ（ｔ）＝｛ＣＬＡ（ｔ）＋ＣＡＡＯ（ｔ）＋ＣＲＲＶ（ｔ）＋ＣＳＭＶ（ｔ）＋ＣＩＶＣ（ｔ）｝／Ｃｍａｘ（ｔ）
図３Ａは、このようにして得られたＤＰＩとスキャンタイミングとの間の関係をプロットした散布図である。ＤＰＩが大きいほど、造影剤の投与位置から遠く、造影剤の到達に時間がかかるため、ＤＰＩは、スキャンタイミングと線形的な相関（ｒ＝０．８６）を有することが分かる。被験者の血行動態が全く同じであれば、スキャンタイミングとＤＰＩは一対一の関係を示すため、同一のスキャンタイミングではＤＰＩはほぼ同じとなり、固定スキャンタイミング（ＦＳＴ）法であっても十分精度よく最適なスキャンタイミングで標的血管、標的臓器または標的疾患について造影ＣＴスキャンを行うことができる。しかし、実際には図３Ａに示すように、同じスキャンタイミングであっても、ＤＰＩは大きく異なる。これは、被験者の血行動態がそれぞれ異なり、造影剤の循環に影響を与えたためである。
次に、各被験者について、腰椎長、胸郭横断径、心胸郭比、下行大動脈の面積、下大静脈の長軸径及び短軸径をＭＰＦとして、各被験者について造影前ＣＴスキャンで測定した。
次に、スキャンタイミング、性別、年齢、身長及び体重のような被験者の基本的情報、並びにＭＰＦを含む様々な説明変数を用いて、ＤＰＩについて段階的回帰分析を行った。回帰分析においては、Ｐ値が０．０５未満である場合には、統計的に有意であるとして考慮した。全ての統計的分析は、ＭＡＴＬＡＢソフトウェア（米国 Ｍａｔｈｗｏｒｋｓ社、２０１６）を使用して行った。
回帰分析の結果、スキャンタイミング（ｒ２＝０．７３３、Ｐ＜０．００１）、下行大動脈の面積（Ｐ＜０．００１）、胸郭横断径（Ｐ＜０．００１）、下大静脈の短軸径（Ｐ＜０．００１）、腰椎長（Ｐ＜０．００１）、下大動脈の長軸径（Ｐ＝０．０３）及び性別（Ｐ＝０．０３）のパラメータが、ＤＰＩについて統計的に有意な説明因子であることが明らかになった。
回帰分析から得られたＤＰＩについての最終的な回帰方程式は以下の通りであった。
ＤＰＩ＝
０．０８１×（スキャンタイミング［秒］）
−０．１４１×（性別（男＝１、女＝０））
＋０．００８×（腰椎長［ｍｍ］）
−０．０１５×（胸郭横断径［ｍｍ］）
−０．０２９×（心胸郭比）
−０．００１×（下行大動脈の面積［ｍｍ２］）
＋０．０１２×（下大静脈の長軸径［ｍｍ］）
−０．０２１×（下大静脈の短軸径［ｍｍ］）
＋３．６９
ＤＰＩ＝
０．０８１×（スキャンタイミング［秒］）
−０．１４１×（性別（男＝１、女＝０））
＋０．００８×（腰椎長［ｍｍ］）
−０．０１５×（胸郭横断径［ｍｍ］）
−０．０２９×（心胸郭比）
−０．００１×（下行大動脈の面積［ｍｍ２］）
＋０．０１２×（下大静脈の長軸径［ｍｍ］）
−０．０２１×（下大静脈の短軸径［ｍｍ］）
＋３．６９
最終的な回帰方程式の決定係数（Ｒ２）は０．８７３であった。図３Ｂに、上記回帰方程式からシミュレートされたＤＰＩを横軸に、実際に観察されたＤＰＩを縦軸にとってプロットした散布図を示す。観察されたＤＰＩは、スキャンタイミング及び様々な形態学的予測因子（ｒ＝０．９３）を含む回帰モデルによってシミュレートされたＤＰＩと線形的な相関関係を示した。図３Ａと比較すると明らかなように、図３Ｂの結果は、ＤＰＩとスキャンタイミングのみとの相関関係よりも明らかに高いことを示している。
この結果から、スキャンタイミングだけでなく、被験者の基本的情報及び造影前ＣＴスキャン像から得られたＭＰＦも、ＤＣＥ−ＣＴにおけるＤＰＩについての顕著な予測因子となりうることが明らかになった。従って、上記の方程式によって、ＤＰＩから被験者ごとに異なる血行動態の影響を排除し、ＤＣＥ−ＣＴにおいて、ＤＰＩを、最適なＣＥＰを得るための定量的指標とすることができることが分かった。
回帰分析の結果から、様々なＭＰＦのうち、胸郭横断径、心胸郭比、下行大動脈の面積、下大静脈の短軸径、及び性別がＤＰＩに対して負の相関を示した。この結果は、これらのＭＰＦが、患者内の血液循環を遅くする因子として相関を有する可能性があることを意味する。心駆出率の低下または循環血液量の増加は、ＣＴスキャンにおけるこれらの形態学的変化に寄与しうる。一方、腰椎長及び下大静脈の長軸径は、ＤＰＩに対して弱い正の相関を示した。この結果は、これらのＭＰＦが患者の血液循環を速くする因子として相関を有する可能性があることを意味する。
一方、ＤＰＩは、前述のように標的血管、標的臓器または標的疾患について、既知の固有の値を有する。従って、上述のＤＰＩについての回帰方程式を変形してスキャンタイミングを左辺に移項し、標的血管、標的臓器または標的疾患についてＤＰＩを得ることができるスキャンタイミングを修正スキャンタイミング（ＭＳＴ）と定義すると、次の方程式が得られる。
ＭＳＴ［秒］＝
１２．４×ＤＰＩ
−４５．６
＋１．７４×（性別（男＝１、女＝０））
−０．１０×（腰椎長［ｍｍ］）
＋０．１９×（胸郭横断径［ｍｍ］）
＋０．３６×（心胸郭比）
＋０．０１×（下行大動脈の面積［ｍｍ２］）
−０．００１×（下大静脈の長軸径［ｍｍ］）
＋０．２６×（下大静脈の短軸径［ｍｍ］）
ＭＳＴ［秒］＝
１２．４×ＤＰＩ
−４５．６
＋１．７４×（性別（男＝１、女＝０））
−０．１０×（腰椎長［ｍｍ］）
＋０．１９×（胸郭横断径［ｍｍ］）
＋０．３６×（心胸郭比）
＋０．０１×（下行大動脈の面積［ｍｍ２］）
−０．００１×（下大静脈の長軸径［ｍｍ］）
＋０．２６×（下大静脈の短軸径［ｍｍ］）
この方程式に、被験者のＭＰＦ及び基本的情報並びに標的血管、標的臓器または標的疾患についての所定のＤＰＩを代入することで、ＭＳＴを決定することができる。
次に、全ての被験者について、右腎静脈、上腸間膜静脈、及び門脈などの各標的血管における動脈対静脈コントラスト（ＡＶＣ）を最大化するスキャンタイミング（ＳＴ）及びＤＰＩの平均値を、ＳＴｍａｘ及びＤＰＩｍａｘとして計算した。これらの平均値（ＳＴｍａｘ及びＤＰＩｍａｘ）を使用して、標的血管内のＡＶＣを２つの方法でシミュレートした。１つは、特定のＤＰＩ（ＤＰＩｍａｘ）及び各被験者のＭＰＦからＭＳＴを決定する上述の方程式を利用してスキャンタイミングを決定する、修正スキャンタイミング（ＭＳＴ）法である。もう１つは、ＳＴｍａｘを用いた固定スキャンタイミング（ＦＳＴ）法である。観察されたデータ間を線形補完した時間濃度曲線を、シミュレーションにおいて使用した。
ＦＳＴ法とＭＳＴ法との比較において、シミュレートされたＡＶＣの平均値の違いを、ｔ検定を用いて統計的に分析した。その際、Ｐ値が０．０５未満である場合、統計的に有意であるとして考慮した。全ての数学的計算及び統計分析は、ＭＡＴＬＡＢソフトウェア（米国 Ｍａｔｈｗｏｒｋｓ社、２０１６）を使用して行った。
造影剤の投与開始後、被験者の標的血管におけるＡＶＣが最大となるスキャンタイミングの平均値ＳＴｍａｘは、右腎静脈では２５．６秒であり、門脈では２７．８秒であり、上腸間膜静脈では２９．３秒であった。被験者の標的血管のＡＶＣが最大となるＤＰＩの平均値ＤＰＩｍａｘは、右腎静脈では１．４であり、門脈では１．６であり、上腸間膜静脈では１．７であった。
図４Ａは、右腎静脈を標的血管とした場合の時間−ＤＰＩ−ＡＶＣ曲線である。図４Ｂは、門脈を標的血管とした場合の時間−ＤＰＩ−ＡＶＣ曲線である。図４Ｃは、上腸間膜静脈を標的血管とした場合の時間−ＤＰＩ−ＡＶＣ曲線である。ＣＴスキャンを行ったスキャンタイミングごとに異なるマーカーで、ＤＰＩとＡＶＣとの関係を散布図にプロットした。スキャンタイミング及び標的血管が同じであっても、被験者の血行動態によりＡＶＣが変動するため、時間−ＡＶＣデータのピークは不明瞭であった。ＤＰＩ−ＡＶＣの鋭いピークは、それぞれ右腎静脈（図４Ａ）、門脈（図４Ｂ）、及び上腸間膜静脈（図４Ｃ）のＤＰＩｍａｘの周囲で見られた。
ＦＳＴ法及びＭＳＴ法による、シミュレートされた標的血管内の平均値及び標準偏差はそれぞれ、右腎静脈において１８１±５９．１ＨＵ及び２０１±７７．１ＨＵ（Ｐ＝０．００２）、門脈において２２０±６５．５ＨＵ及び２３２±７１．９ＨＵ（Ｐ＝０．０２９）、上腸間膜静脈において２４４±７０．４ＨＵ及び２５４±７１．５ＨＵ（Ｐ＝０．００９）であった。ＨＵはＨｏｕｎｓｆｉｅｌｄ Ｕｎｉｔの略であり、Ｘ線吸収係数である。
シミュレートされた、標的血管におけるＡＶＣは、ＭＳＴ法ではＦＳＴ法と比較して顕著に高かった。ＦＳＴ法とＭＳＴ法とを比較した、右腎静脈におけるＡＶＣのブランド−アルトマンプロットを図５Ａに、門脈におけるＡＶＣのブランド−アルトマンプロットを図５Ｂに、上腸間膜静脈におけるＡＶＣのブランド−アルトマンプロットを図５Ｃに、それぞれ示す。ＨＵはＨｏｕｎｓｆｉｅｌｄ Ｕｎｉｔの略であり、Ｘ線吸収係数である。
図５Ａから５Ｃは、ＭＳＴ法によるＡＶＣの値であるＡＶＣｍｓｔと、ＦＳＴ法によるＡＶＣの値であるＡＶＣｆｓｔとの差を縦軸にプロットし、ＡＶＣｍｓｔとＡＶＣｆｓｔの平均値を横軸にプロットした。右腎静脈（図５Ａ）、門脈（図５Ｂ）及び上腸間膜静脈（図５Ｃ）で、ＦＳＴ法によるＡＶＣよりもＭＳＴ法によるＡＶＣの方が高いという固定バイアスが見られた。
以上の結果から、ＭＳＴ法はＦＳＴ法と比較して、ＤＣＥ−ＣＴにおける標的血管のＡＶＣを向上させることができることが明らかになった。様々な標的血管、標的臓器または標的疾患において、様々な投与プロトコルの下で、より良好なＡＶＣ及びＬＯＣまたはより良好なコントラストを得るために推奨されるいくつかのＦＳＴが報告されている。しかし、これらのスキャンタイミングは、被験者の集団内の最良のスキャンタイミングとなりうるが、個々の被験者の血行動態はＦＳＴ法では考慮されていないため、個々の被験者については、必ずしもＦＳＴ法によるスキャンタイミングが最良になるとは限らない。ＦＳＴ法で用いられるＦＳＴが、ある集団内で決定されたとしても、ＭＳＴ法は、ＤＣＥ−ＣＴにおける最良のスキャンタイミングの予測において、ＦＳＴ法よりも有利である。従って、造影前ＣＴにおけるＭＰＦを用いた被験者の血行動態の推定が、ＭＳＴ法による高いＡＶＣを有する、より良好な血管画像を得るための最良のスキャンタイミングの予測に有効であり、有用であると結論付けることができる。
実施形態において、被験者の血行動態の指標として、造影前ＣＴによって得られるＭＰＦを用いることができることを示した。超音波検査法や核医学検査においても、被験者の血行動態の正確な評価が必要になる可能性がある。しかし、これらの検査法では、血行動態についての従来の評価は、利用可能性が低く、追加的な時間や高額な検査医療費などのいくつかの欠点を有していた。この点において、ＭＰＦを用いたＭＳＴ法は、被験者の血行動態を予測するためにさらなる検査を必要とせず、最適な解決手段となりうる。さらに、造影前ＣＴによるＭＰＦの測定及びＭＳＴ法における修正スキャンタイミングの計算は複雑なものではなく、医療現場において幅広く利用可能である。
以上の点から、本発明は、ＭＰＦ及びＤＰＩを使用するという２つの点で特に新規性を有する。本発明は、被験者の血行動態が、画像上の形態学的変化に反映されうることを示した。従来は、ＤＰＩなどのＤＣＥ−ＣＴにおける造影時相の定量的方法を利用することはできなかったため、画像上のＭＰＦをＤＣＥ−ＣＴにおける最良のスキャンタイミングの定量的予測に利用することは、本発明により初めて可能になった。標的血管内のＤＰＩｍａｘは、造影剤の到達時間に従って次第に増加するため、ＤＰＩは、ＤＣＥ−ＣＴの造影時相に関する定量的かつ汎用的指標となりうる。そのため、ＭＳＴ法は、本発明の実施形態で示した血管画像の取得以外の、様々な標的器官や病変における最良のスキャンタイミングの予測に適用可能である。
さらに、ＭＳＴ法から得られた修正スキャンタイミングは、その他の方法にも使用可能となる。例えば、ＭＳＴ法とボーラストラッキング法とを組み合わせれば、ＭＳＴ法に依って最適なスキャンタイミングを予測することができるため、より少ないサンプリング画像またはより少ない放射線被曝量で、高いコントラストのＤＣＥ−ＣＴ画像を得ることが可能になりうる。
従って、造影時相の定量的指標であるＤＰＩ及び造影前ＣＴから得られるＭＰＦを用いて、ＤＣＥ−ＣＴにおけるＣＥＰ及び最良のスキャンタイミングの定量的かつ汎用的な予測が可能となる。
Claims (17)
- ＣＴスキャン撮像ユニットと、制御ユニットと、を含むＣＴスキャン装置であって、
前記制御ユニットが、被験者の基本的情報及び標的血管、標的臓器または標的疾患のダイナミック時相指数を受け取るように構成され、
前記ＣＴスキャン撮像ユニットが、前記被験者の第１のＣＴスキャンを実行するように構成され、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットから第１のＣＴスキャン像を受信するように構成され、
前記制御ユニットが、前記第１のＣＴスキャン像から、前記被験者の少なくとも１つの形態学的予測因子を測定するように構成され、
前記制御ユニットが、前記被験者の基本情報、前記ダイナミック時相指数、及び前記形態学的予測因子に基づいて、造影剤投与の開始から第２回目以降のＣＴスキャンを行うまでの時間を、修正スキャンタイミングとして決定するように構成され、
前記ＣＴスキャン撮像ユニットが、前記造影剤の投与の開始から前記修正スキャンタイミングが経過したときに前記第２回目以降のＣＴスキャンを実行するように構成され、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットから第２回目以降のＣＴスキャン像を受信するように構成された、
ＣＴスキャン装置。 - 前記ＣＴスキャン装置が、造影剤投与ユニットをさらに含み、
前記造影剤投与ユニットが、前記制御ユニットの命令により、前記被験者に造影剤を投与するように構成された、請求項１に記載のＣＴスキャン装置。 - 前記被験者の基本的情報が、性別、年齢、身長及び体重のうち少なくとも１つを含む、請求項１に記載のＣＴスキャン装置。
- 前記形態学的予測因子が、下行大動脈の面積及び下大静脈の短軸径のうち少なくとも１つを含む、請求項１に記載のＣＴスキャン装置。
- 前記形態学的予測因子が、腰椎長、胸郭横断径、心胸郭比、大動脈長、上行及び下行大動脈の面積、腹部大動脈の面積、下大静脈の長軸径並びに下大静脈の短軸径のうち少なくとも１つを含む、請求項１に記載のＣＴスキャン装置。
- 前記修正スキャンタイミングが、前記ダイナミック時相指数、前記被験者の基本情報及び前記形態学的予測因子の少なくとも１つを説明変数として、統計解析もしくは非線形機械学習により決定された方程式または非線形モデルから得られる、請求項１に記載のＣＴスキャン装置。
- 前記修正スキャンタイミングを求めるための方程式が、前記ダイナミック時相指数、前記被験者の基本情報及び前記形態学的予測因子の少なくとも１つを含む線形方程式または非線形モデルである、請求項６に記載のＣＴスキャン装置。
- 前記修正スキャンタイミングを求めるための方程式が、対象とする被験者の母集団に応じて統計解析を行う被験者のサンプルを選択して係数を最適化した、請求項６または７に記載のＣＴスキャン装置。
- 前記修正スキャンタイミングが下記の方程式から得られる、請求項６に記載のＣＴスキャン装置：
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）。 - ＣＴスキャン撮像ユニットと、造影剤投与ユニットと、制御ユニットと、を含むＣＴスキャン装置の作動方法であって、
前記制御ユニットが被験者の基本的情報及び標的血管、標的臓器または標的疾患のダイナミック時相指数を受け取る段階と、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットに、前記被験者の第１のＣＴスキャンを実行させる命令を送信する段階と、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットから第１のＣＴスキャン像を受信する段階と、
前記制御ユニットが、前記第１のＣＴスキャン像から、前記被験者の少なくとも１つの形態学的予測因子を測定する段階と、
前記制御ユニットが、前記被験者の基本情報、前記ダイナミック時相指数、及び前記形態学的予測因子に基づいて、造影剤投与の開始から第２回目以降のＣＴスキャンを行うまでの時間を、修正スキャンタイミングとして決定する段階と、
前記制御ユニットが、前記造影剤投与ユニットに、造影剤の投与を開始させる命令を送信する段階と、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットに、造影剤投与の開始から、前記修正スキャンタイミングが経過したときに前記第２回目以降のＣＴスキャンを実行させる命令を送信する段階と、
前記制御ユニットが、前記ＣＴスキャン撮像ユニットから第２回目以降のＣＴスキャン像を受信する段階と、
を含む、方法。 - 前記被験者の基本的情報が、性別、年齢、身長及び体重のうち少なくとも１つを含む、請求項１０に記載の方法。
- 前記形態学的予測因子が、下行大動脈の面積及び下大静脈の短軸径のうち少なくとも１つを含む、請求項１０に記載の方法。
- 前記形態学的予測因子が、腰椎長、胸郭横断径、心胸郭比、大動脈長、上行及び下行大動脈の面積、腹部大動脈の面積、下大静脈の長軸径並びに下大静脈の短軸径のうち少なくとも１つを含む、請求項１０に記載の方法。
- 前記修正スキャンタイミングが、前記ダイナミック時相指数、前記被験者の基本的情報及び前記形態学的予測因子の少なくとも１つを説明変数として、統計解析もしくは非線形機械学習により決定された方程式または非線形モデルから得られる、請求項１０に記載の方法。
- 前記修正スキャンタイミングを求めるための方程式が、前記ダイナミック時相指数、前記被験者の基本的情報及び前記形態学的予測因子の少なくとも１つを含む線形方程式または非線形モデルである、請求項１４に記載の方法。
- 前記修正スキャンタイミングを求めるための方程式が、対象とする被験者の母集団に応じて統計解析を行う被験者のサンプルを選択して係数を最適化した、請求項１４または１５に記載の方法。
- 前記修正スキャンタイミングが下記の方程式から得られる、請求項１４に記載の方法：
修正スキャンタイミング［秒］＝１２．４×（ダイナミック時相指数）−４５．６＋１．７４×（性別（男＝１、女＝０））−０．１０×（腰椎長［ｍｍ］）＋０．１９×（胸郭横断径［ｍｍ］）＋０．３６×（心胸郭比）＋０．０１×（下行大動脈の面積［ｍｍ２］）−０．００１×（下大静脈の長軸径［ｍｍ］）＋０．２６×（下大静脈の短軸径［ｍｍ］）。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017036115A JP6899666B2 (ja) | 2017-02-28 | 2017-02-28 | 造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017036115A JP6899666B2 (ja) | 2017-02-28 | 2017-02-28 | 造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018139901A JP2018139901A (ja) | 2018-09-13 |
| JP6899666B2 true JP6899666B2 (ja) | 2021-07-07 |
Family
ID=63527217
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017036115A Active JP6899666B2 (ja) | 2017-02-28 | 2017-02-28 | 造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6899666B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110084429A (zh) * | 2019-04-29 | 2019-08-02 | 东软医疗系统股份有限公司 | 扫描时间的预测方法、装置、存储介质及电子设备 |
| JP7312611B2 (ja) * | 2019-05-30 | 2023-07-21 | キヤノンメディカルシステムズ株式会社 | Ｘ線診断装置 |
| CN110584699B (zh) * | 2019-09-30 | 2023-09-26 | 东软医疗系统股份有限公司 | 一种ct扫描的控制方法、装置及终端设备 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5687208A (en) * | 1995-10-06 | 1997-11-11 | Bhb General Partnership | Method of and apparatus for predicting computed tomography contrast enhancement with feedback |
| JP4537681B2 (ja) * | 2003-09-24 | 2010-09-01 | 株式会社東芝 | 血流解析装置 |
| US20140086461A1 (en) * | 2012-09-25 | 2014-03-27 | The Johns Hopkins University | Method and system for determining time-based index for blood circulation from angiographic imaging data |
| JP6479346B2 (ja) * | 2014-06-05 | 2019-03-06 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
| CN104287764B (zh) * | 2014-09-11 | 2017-05-31 | 沈阳东软医疗系统有限公司 | 一种ct灌注成像方法和设备 |
| CN104688259B (zh) * | 2015-02-16 | 2018-01-09 | 袁小东 | Ct扫描灌注方法及装置 |
-
2017
- 2017-02-28 JP JP2017036115A patent/JP6899666B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2018139901A (ja) | 2018-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9949650B2 (en) | Fractional flow reserve estimation | |
| US8285014B2 (en) | Measuring blood volume with C-arm computed tomography | |
| Greupner et al. | Head-to-head comparison of left ventricular function assessment with 64-row computed tomography, biplane left cineventriculography, and both 2-and 3-dimensional transthoracic echocardiography: comparison with magnetic resonance imaging as the reference standard | |
| Davis et al. | 4D digital subtraction angiography: implementation and demonstration of feasibility | |
| US9633454B2 (en) | Fluid-dynamic analysis of a vascular tree using angiography | |
| Lawler et al. | MDCT evaluation of the coronary arteries, 2004: how we do it—data acquisition, postprocessing, display, and interpretation | |
| Fursevich et al. | Bariatric CT imaging: challenges and solutions | |
| US11406339B2 (en) | System and method for determining vascular velocity using medical imaging | |
| EP2838430B1 (en) | Perfusion imaging | |
| JP6899666B2 (ja) | 造影ｃｔスキャン装置、造影ｃｔスキャン装置の作動方法、及び造影ｃｔスキャン像を得るための方法 | |
| Mao et al. | Dual-standard reference values of left ventricular volumetric parameters by multidetector CT angiography | |
| Entezari et al. | Analysis of the thoracic aorta using a semi-automated post processing tool | |
| Tanaka et al. | Detection of pulmonary embolism based on reduced changes in radiographic lung density during cardiac beating using dynamic flat-panel detector: an animal-based study | |
| Kidoh et al. | Novel contrast-injection protocol for coronary computed tomographic angiography: contrast-injection protocol customized according to the patient’s time-attenuation response | |
| Seifarth et al. | Introduction of an individually optimized protocol for the injection of contrast medium for coronary CT angiography | |
| Zellerhoff et al. | Measurement of cerebral blood volume using angiographic C-arm systems | |
| Jiang et al. | Preoperative evaluation of anomalous pulmonary venous connection using dual-source computed tomography: Comparison with echocardiography | |
| Iezzi et al. | Tailoring protocols for chest CT applications: when and how? | |
| Farrelly et al. | Low dose dual-source CT angiography of the thoracic aorta | |
| Kondo et al. | Prospective ECG-gated coronary 320-MDCT angiography with absolute acquisition delay strategy for patients with persistent atrial fibrillation | |
| Nasis et al. | Assessment of left ventricular volumes, ejection fraction and regional wall motion with retrospective electrocardiogram triggered 320-detector computed tomography: a comparison with 2D-echocardiography | |
| Stauder et al. | Coronary artery bypass grafts: assessment of graft patency and native coronary artery lesions using 16-slice MDCT | |
| US20150282779A1 (en) | Treating an Ischemic Stroke | |
| Bruners et al. | Assessment of global left and right ventricular function using dual-source computed tomography (DSCT) in comparison to MRI: an experimental study in a porcine model | |
| Reynolds et al. | Toward improved 3D carotid artery imaging with Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200225 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210126 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210301 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210413 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210517 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210615 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6899666 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |