Bayesian CT Perfusion Imaging
in Ischemic Stroke

Thomas Pickersgill, Global Clinical Marketing Manager, Healthcare IT, Canon Medical Systems
Jeroen Tijhaar, Application Development Specialist, CT, Canon Medical Systems
Chloe Steveson, Senior Application Development Specialist, CT, Canon Medical Systems
Erin Angel, Senior Manager, Clinical Collaborations, Canon Medical Systems USA

Introduction

Stroke accounts for more than one in every 6 deaths in the United States1. In the treatment of cerebral ischemia, it is often said that “time is brain.” The faster the clinician can detect areas of decreased blood flow and determine the optimal treatment pathway, the better the patient’s chances for survival and recovery.
Advanced imaging now appears to have established its role in the standard workup of acute stroke. Recently published trials have demonstrated the benefit of advanced imaging for selecting patients who may benefit from endovascular treatment after the traditional time window of 6 hours after ictus2,3,4. Moreover, recent updates of the stroke definitions by the American Heart Association and American Stroke Association have large shifted from a clinical diagnosis to an advanced imagingbased diagnosis of stroke and stroke subtypes5.
In one clinical trial the researchers demonstrated that 45% of patients treated with the clot removal procedure within the 16-hour time frame achieved functional independence 90 days after treatment, compared with 17% in the control group. Thrombectomy was also associated with improved survival: 14% of the treated group deceased within 90 days of the study, compared with 26% in the control group4.
Similar to the aforementioned trial, another multicenter stroke treatment study investigated the response to therapy for stroke patients. This trial monitored the impact of mechanical thrombectomy patients who were last seen well 6-24 hour prior to treatment. Patients in the thrombectomy group were less impaired at 90 days posttreatment compared to the standard care group. Rates of functional independence at 90 days were 49% in the thrombectomy group, compared to 3% in the control group3.
For perfusion analysis to be effective and accurate, the data needs to be taken with an imaging modality capable of acquiring temporally uniform dynamic images of the entire anatomy, have sufficient temporal sampling, and the perfusion algorithm can represent flow characteristics independent of the image acquisition and contrast injection processes.
CT perfusion was traditionally confined to imaging only a portion of the brain. With the introduction of the Aquilion ONE area detector CT in 2007, the scope of cerebral perfusion analysis was changed by enabling dynamic volume imaging of the entire brain with isophasic and physiological uniformity.
Dynamic volume, whole brain imaging on Aquilion ONE was recently paired with the Bayesian perfusion algorithm to produce advanced CT perfusion imaging. Bayesian is the latest advancement in CT perfusion algorithm science. This algorithm was originally developed for MRI DWI and PWI analysis6. The Bayesian algorithm is delay insensitive and adapts for contrast injection timing and flow characteristics. A key feature of the Bayesian algorithm is a lower noise CTP map7, thus improving the visualization and quantification of CTP maps used for clinical diagnosis.
This paper reviews the principles of brain perfusion, discusses methods of measuring perfusion, and describes the unique benefits of whole brain CT perfusion imaging used in conjunction with the Bayesian perfusion algorithm.
The primary application of brain perfusion imaging in CT is to determine if the patient has had a stroke and to assess the viability of the brain tissue. Stroke diagnosis and management is a significant focus in today’s hospitals because stroke is the fifth leading cause of death in the United States and a major cause of serious long-term disability8,9.
When a patient presents with stroke symptoms, a non-contrast CT is typically used to visualize bleeding and to rule out haemorrhagic stroke. However, in 85% of strokes there is no intra-cranial bleed10. Once a haemorrhagic stroke has been ruled out, a CT perfusion examination will be performed. Perfusion measurements are used to visualize effects of ischemic stroke. Perfusion measurements help distinguish which areas of the brain show signs of irreversible repair and which areas of the brain may be saved through intervention. In this way, perfusion analysis can help clinicians estimate treatment response and develop therapeutic pathways designed specifically for individual patients. Routinely after a haemorrhagic stroke has been ruled out an angiogram of the carotid arteries will also be performed to assess possible causes of the patient symptoms10.
Cerebral CT perfusion is discussed primarily in the context of characterizing stroke, but perfusion measurements are also valuable for a range of clinical applications such as the evaluation of vasospasm11, vasculitis12, assessment of perfusion after head trauma13 and determining microvascular permeability in brain tumours14.

Pathophysiology

Healthy brain tissue relies on continuous flow of oxygenated blood and requires a precise pressure balance to remain viable. Blood is supplied to the brain via four major arteries: the left and right carotid and vertebral arteries. When the brain is working properly, all blood flows into the brain from arterial inputs, remains in the cerebral vasculature (arteries, capillaries or vessels) and exits the organ through the veins (Figure 1).
Since the brain is surrounded by the skull, there is little room to expand to accommodate pressure changes. To maintain the delicate pressure balance, the cerebral vasculature has a unique autoregulation mechanism to automatically adjust flow. Blood flow is regulated primarily through vasodilation, vasoconstriction and collateralization. Vasodilation and vasoconstriction are the automatic dilation and constriction of vessels to regulate blood flow and maintain blood pressure. If a region of the brain is receiving inadequate blood supply, the vessels will automatically dilate to restore blood flow to the region. Likewise, if there is too much blood pressure in an area, the vessels will constrict to reduce flow. Collateralization is the redirection of blood to a region using small “detour” vessels called collateral arteries. Collateral arteries are typically small (often closed) arteries that can open, expand, or extend to redirect blood around a blockage. Autoregulation can compensate for small or transient changes in blood pressure, but once the vessels have reached their expansion limit, and if collateral flow is insufficient, the brain tissue will become ischemic (deprived of blood flow). When an area of the brain becomes ischemic, there are two possible consequences for the cells in that region. If the cells have been deprived of blood for an extended time, they will eventually become irreversibly damaged and die. These cells are unsalvageable and form a region called the ischemic core. Surrounding the infarct core are cells that are starved for blood but have not yet died and if blood flow can be restored, the cells will recover. This salvageable region is called the ischemic penumbra (Figure 2).
For stroke patients, perfusion imaging is used to visualize ischemic tissue, and to quantify the infarct core and penumbra. The infarct core is used to diagnose acute stroke or confirm suspected diagnosis of stroke. The ischemic penumbra is measured to determine whether the patient is a good candidate for revascularization therapy. Revascularization can be performed interventionally using a clot retrieval device or by using pharmaceutical agents such as thrombolytic reperfusion agents. Thrombolytic agents can dissolve blood clots and return flow to ischemic regions. However, in areas of ischemic core, these agents can cause increased risk of hemorrhage. In recent trials, thrombolytics have proven effective if administered up to 24h after stroke onset3,4.
Canon Medical Systems’ Neuro ONE protocol allows full stroke workup in a single examination that includes physiological and anatomical information about the entire brain15.
This uniquely comprehensive exam combines whole brain dynamic perfusion maps to analyse blood flow and characterize brain tissue viability, as well as a 4D CT Digital Subtraction Angiogram (4D CT DSA) to help visualize obstructed vessels for treatment planning12,15-18. The entire exam is performed in less than five minutes with low contrast dose (50 mL of intravenous contrast) and low radiation dose (typically less than 5 mSv for the entire exam).

CT or MR

In MR, perfusion weighted imaging (PWI) can be combined with diffusion- weighted imaging (DWI) to characterize the penumbra. Perfusion CT can be used to quickly visualize stroke effects by monitoring the flow of blood through the cerebral vasculature. Studies have demonstrated that in terms of patient selection for reperfusion therapies, CT and MR are equally proficient in characterizing regions of infarct and penumbra19,20.
Traditionally CT perfusion was limited by partial brain coverage (typically 4 cm or less), partial volume effect of vessels, and artefact in the posterior fossa (Figure 3). Today’s CT scanners are now capable of high resolution, three-dimensional imaging, essentially eliminating the previous concerns. The partial brain coverage limitations of CT perfusion were overcome with the introduction of Aquilion ONE. Dynamic volume perfusion imaging can be performed for the entire brain with 16 cm of anatomical coverage and with temporal uniformity (Figure 4). The Aquilion ONE technology has overcome the previous limitations of CT perfusion, making it an attractive alternative to MR for stroke evaluation. Another clear advantage of CT for acute stroke workup is that there is no contraindication for patients with metallic or electrical implants such as aneurysm clips or pacemakers (which can be common in stroke patients). Also, patient monitors or ventilators containing metal can be used during the CT exam. In the case of the Neuro ONE protocol on Aquilion ONE, the entire stroke workup can be combined into a single exam, thus minimizing contrast and radiation dose, and minimizing exam time.
Although advanced imaging now seems to have established its role in the standard workup of acute stroke, there is still debate whether CT perfusion imaging or MR imaging should be used to visualize the infarct core and penumbra21,22. On the one hand, CT perfusion provides a more practical approach that does not delay treatment decisions since patients already undergo non-contrast CT and CT angiography imaging. On the other hand, MR diffusion weighted imaging (DWI) is still the reference standard in most centers for defining the ischemic core and there may be a mismatch between infarct measurements using CT perfusion and DWI23. It has been postulated that these differences may be a result of high image noise and low signal to noise ratios in CT perfusion imaging compared to DWI24.

CT Perfusion – Basic Principles

To measure cerebral perfusion using CT, intravenous contrast agent is administered to the blood stream and a series of CT images are acquired over time to observe arterial input, tissue uptake, and venous outflow of the contrast agent. Since blood remains in the vasculature, the tissue uptake is actually a measure of the blood in small vasculature and capillaries. As the contrast-labelled blood enters the anatomy, the contrast density increases to a peak enhancement and then decreases as that blood washes out of the region. Contrast density is determined by measuring the temporal change in CT number of an input artery, the brain tissue and an output vein. For each measurement, a graph is generated which plots the CT number (Hounsfield unit) versus time. These graphs, called time density curves, represent the uptake and washout of contrast-labelled blood in an artery, in brain tissue and in a vein (Figure 5).
In order to generate cerebral perfusion maps, it is necessary to generate a separate tissue TDC for every voxel in the brain. For simplicity, a single artery input and venous output are tracked throughout the dynamic scan. This results in one arterial TDC, one venous TDC and several tissue TDCs — one for each voxel in the brain. To generate the arterial and venous TDCs, regions of interest (ROIs) are automatically placed in optimally chosen large vessels that are perpendicular to the scan acquisition plane. Typically, the arterial ROI is chosen in an anterior cerebral artery or a middle cerebral artery. The venous ROI is often placed over the superior sagittal sinus, transverse sinus or torcular herophili.

Temporal Sampling

To perfectly sample the time density curves, imaging would have to be performed continuously throughout the exam (like filming with a video camera). But to minimize dose, images are only acquired frequently enough to reliably reconstruct the curves. To accomplish this, images are acquired intermittently as the contrast agent moves through the brain. CT numbers are recorded at each time point and time density curves are generated (Figure 5). The Aquilion ONE has the advantage that the entire brain can be imaged without table movement and, therefore, the sampling rate is not limited by the time it would take for the table to move the length of the brain. This dynamic volume imaging allows the sampling rate to be adjusted, depending on the application, to ensure an accurate measurement of the TDC while maintaining low dose.
sou
Figure 5   In order to create brain perfusion maps, the Time Density Curve (TDC) is measured in an artery, in a vein, and at each voxel in the brain tissue to visualize the uptake of contrast-labeled blood.

Temporal Uniformity

To obtain accurate and complete perfusion maps, the scanner must image the entire brain at a single instance in time so that the flowing intravenous contrast does not have time to change during the acquisition. With 16 cm of coverage, the Aquilion ONE CT scanner is capable of dynamic volume imaging of the entire brain with temporal uniformity so that all contrast-labelled blood is visualized at the same point in time.

Quantitative Maps

In perfusion studies, the TDCs are used to calculate several parameters at each location in the brain thus generating quantitative perfusion-related maps, which help clinicians characterize cerebral pathophysiology. (Figure 6). Each of the perfusion maps can indicate pathologies or different flow irregularities. Longer TTP or Delay may indicate delayed flow from collateral autoregulation. MTT, TTP, and Delay visualize vasospasm or vessel stenosis that cannot otherwise be seen.
Consider, for example, a patient with an ischemic stroke caused by an occlusion in the middle cerebral artery. The occlusion will reduce the perfusion to the affected region of the brain as demonstrated by decreased CBF. This reduced flow will trigger an autoregulatory response and feeding arteries will dilate to restore flow to the region. This vasodilation is depicted as increased MTT (analogous to the length of time it takes for a log to float down the wider parts of a river as opposed to the rapids). If blood is restored via collateral circulation, the more tortuous path of collateral vessels means it will take more time for blood to reach the region. This delayed time of arrival is visible in the TTP and Delay maps. Autoregulation maintains blood supply to the region as demonstrated by an increase or maintained CBV until the vessels can no longer dilate and collateralization cannot maintain cerebral reserve, at which point CBV will decrease demonstrating infarct. In this example, the mismatch between CBF and CBV represents the ischemic penumbra.

Perfusion algorithms

There are several ways to calculate perfusion parameters from the TDC’s, including the max. slope method, gamma-variate, moments-method, deconvolution and Bayesian methods.

SVD+

There are several types of deconvolution that can be used for perfusion analyses. Singular Value Decomposition (SVD+) is a delay insensitive SVD algorithm that uses an innovative technique to account for delayed blood flow and perform calculations with fast computation times. SVD algorithms are often viewed as the most accurate method because they rely on less assumptions than other methods, however these algorithms are sensitive to noise and in some occasions, can result in underestimation of high MTT and CBF values and overestimation of low CBF values in certain conditions.

Bayesian

Canon Medical Systems has now introduced the Bayesian post processing algorithm to the CT perfusion application. The Bayesian method is a probabilistic algorithm, which is based on the Bayes theorem25. The Bayesian algorithm aids clinicians to evaluate CT perfusion images with accurate perfusion mapping which enables a more confident diagnosis.
Bayes theory of probability was first published in the 18th century. The theory allows the combination of experimental data and priori information about the parameters of a model (such as CBF and MTT values), to generate a robust probability distribution for these parameters26. Thanks to this probabilistic approach, the parameter estimations are more robust25.
The Bayesian method uses a probabilistic approach to CT perfusion, which is delay insensitive. The algorithm calculates a delay in CT perfusion independent of MTT values, giving highly accurate estimates of delay values. This has been validated in a number of studies6,7,27-29. Sometimes when performing CT perfusion high levels of noise often prevent the accurate estimation of penumbra volumes when using deconvolution based methods. In comparison, the Bayesian method is more robust against noise30.
In CT perfusion studies, Bayesian has been reported to outperform SVD based algorithms, in particular for CBF and MTT computation7. On VitreaTM, the Bayesian algorithm produces similar CT perfusion maps as SVD+ (CBV, CBF, MTT, TTP, Delay).
The Bayesian and SVD+ algorithms are available for all types of acquisition protocols from all CT systems.

Radiation and contrast dose

Being a dynamic scan, it is imperative that CT perfusion protocols are designed to maintain low radiation exposure and low doses of iodinated contrast, while ensuring diagnostic quality of the perfusion maps. The Neuro ONE protocol is designed with patient safety in mind. The protocol combines multiple exam types (CT perfusion and full brain 4D CT angiography) into a single acquisition while maintaining a low radiation dose of 5 mSv or less. These low doses are achievable because of a combination of scanner characteristics such as fast rotation time and low dose intermittent scan acquisition. By combining exam types into a single acquisition, the Neuro ONE protocol enables a single low dose injection of iodinated contrast (around 50 mL) administered with an injection rate of 5 mL/sec.

Conclusion

The Aquilion ONE provides a uniquely comprehensive exam to reduce diagnosis time for patients experiencing serious cerebrovascular conditions, such as stroke. By pairing low dose whole brain imaging with the Bayesian perfusion algorithm, the Aquilion ONE produces advanced and accurate CT perfusion maps for evaluating cerebral blood flow and brain tissue viability.
Bayesian methods are less sensitive to noise and could result in better signal-to-noise ratios and image quality, and therefore may improve diagnostic performance for infarct detection with the potential to reduce radiation dose.
Recent clinical trials have demonstrated the value of CT perfusion imaging to define infarct and regions of penumbra in stroke patients, to enable clinical decisions and treatment planning. The same studies have shown that patients who received treatment up to 24h after the onset of a stroke benefited significantly from improved functional outcomes in comparison to control groups3,4.
As a result, the Aquilion ONE and Vitrea solution produces maps providing more information physicians can use to accurately and quickly diagnose cerebrovascular disorders.

Clinical Example

sou
Images courtesy of Prof. Prokop, Radboud University Medical Center, Nijmegen, the Netherlands

Figure 7  False positive perfusion deficit was seen on the original perfusion maps in right hemisphere. The study was then processed using Bayesian based perfusion maps, which were read as true negative. Which correlates with the DWI images.

References

  1. Benjamin, EJ, et al. 2018. Correction to: Heart Disease and Stroke Statistics—2018 Update: AReport from the American Heart Association. Circulation, 137(12).
  2. Powers WJ, Rabinstein AA, Ackerson T, et al. 2018 Guidelines for the Early Management of Patient with Acute Ischemic Stroke: A Guideline for Healthcare Professionals from the American Heart Association/American Stroke Association. Stroke. 2018;49(3):e46-e11
  3. Nogueira RG, et al. "Thrombectomy 6 to 24 Hours after Stroke with a Mismatch between Deficit and Infarct". The New England Journal of Medicine. 2018; 378: 11-21.3.
  4. Albers GW et al, Thrombectomy for stroke at six to 16 hours with selection by perfusion imaging. New England Journal of Medicine, February 22, 2018
  5. Sacco RL, Kasner SE, Broderick JP, et al. AN updated definition of stroke for the 21st century: a statement for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. 2018;49(3)E46-E110
  6. Boutelier T, Kudo K, Pautot F, Sasaki M. Bayesian hemodynamic parameter estimation by bolus tracking perfusion weighted imaging. IEEE Trans Med Imaging.2012,31:1381-95.
  7. Sasaki M, Kudo K, Boutelier T, Pautot F, Christensen S, Uwano I, Goodwin J, Higuchi S, Ito K, Yamashita F. Assessment of the accuracy of a Bayesian estimation algorithm for perfusion CT by using a digital phantom. Neuroradiology 2013. DOI 10.1007/ sOO234-013-1237-7.
  8. Kochanek KD, Xu JQ, Murphy SL, Arias E. Mortality in the United States, 2013. NCHS Data Brief, No. 178. Hyattsville, MD: National Center for Health Statistics, Centers for Disease Control and Prevention, Department of Health and Human Services; 2014.
  9. Mozzafarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, et al., on behalf of the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Heart disease and stroke statistics—2016 update: a report from the American Heart Association. Circulation 2016;133(4):e38–360.
  10. Wintermark, M. et al., 2013. Imaging Recommendations for Acute Stroke and Transient Ischemic Attack Patients: A Joint Statement by the American Society of Neuroradiology, the American College of Radiology, and the Society of Neuro Interventional Surgery. American Journal of Neuroradiology, 34(11).
  11. Shankar, J.S. & Lum, C., 2011. Whole brain CT perfusion on a 320-slice CT scanner. Indian Journal of Radiology and Imaging, 21(3), p.209.
  12. Yahyavi-Firouz-Abadi N, Wynn BL, Rybicki FJ, et al. Steroidresponsive large vessel vasculitis: application of whole-brain 320-detector row dynamic volume CT angiography and perfusion. AJNR Am J Neuroradiol. 2009;30(7)1409-1411.
  13. Wintermark M, van Mlle G, Schnyder P, et al. Admission perfusion CT: prognostic value in patients with severe head trauma. Radiology. 2006;232 (1):211-220
  14. Chen, T. et al., 2014. Preliminary study of whole-brain CT perfusion imaging in patients with intracranial tumors adjacent to large blood vessels. Clinical Radiology, 69(1).
  15. Angel E. SVD+ Dynamic Volume CT: Delay Insensitive Brain Perfusion Analysis. White paper. Toshiba America Medical Systems, Inc. 2010
  16. San Millän Ruiz D, Murphy K, Gailloud P. 320- Multidetector row whole-head dynamic subtracted CT angiography and wholebrain CT perfusion before and after carotid artery stenting: Technical note EurJ Radiol. 2009. Available at: http://www. ncbi.n/m.nih.gov/pubmed/19410406 [Accessed January 29, 2010].
  17. Klingebiel R, Siebert E, Diekmann S, et al. 4- D Imaging in cerebrovascular disorders by using 320-slice CT: feasibility and preliminary clinical experience. Acad Radiol.2009;16(2)123-129
  18. Salomon EJ Barrett J Willems PWA et al Dynamic CT angiography and CT perfusion employing a 320 - detector row CT: protocol and current clinical applications Klin Neuroradiol. 2000;19(3):187-196
  19. Liu HL, Pu Y, Liu Y, et al. Cerebral blood flow measurement by dynamic contrast MRI using singular value decomposition with an adaptive threshold. Magn Reson Med. 1999;42(1):167–172
  20. Wirestam R, Andersson L, Ostergaard L, et al. Assessment of regional cerebral blood flow by dynamic susceptibility contrast MRI using different deconvolution techniques. Magn Reson Med. 2000;43(5):691–700.18. Konstas AA, Lev MH. CT perfusion imaging of acute stroke: the need for arrival time delay insensitive, and standardized postprocessing algorithms? Radiology. 2010;254(1):22-25
  21. García-Bermejo P, Castaño C, Dávalos A. Multimodal CT versus MRI in Selecting Acute Stroke Patients for Endovascular Treatment. Interv Neurol. 2013;1(2):65–76.
  22. Warach SJ, Luby M, Albers GW, et al. Acute Stroke Imaging Research Roadmap III Imaging Selection and Outcomes in Acute Stroke Reperfusion Clinical Trials: Consensus Recommendations and Further Research Priorities. Stroke. 2016;47(5):1389–1398.
  23. Copen WA, Morais LT, Wu O, et al. In Acute Stroke, Can CT Perfusion-Derived Cerebral Blood Volume Maps Substitute for Diffusion-Weighted Imaging in Identifying the Ischemic Core? PLoS One. 2015;10(7):e0133566.
  24. González RG. Low signal, high noise and large uncertainty make CT perfusion unsuitable for acute ischemic stroke patient selection for endovascular therapy. J Neurointerv Surg. jnis.bmj. com; 2012;4(4):242–245.oen
  25. Boutelier, T et al. Olea Imagein. Issue 1. March 2016. Olea Medical
  26. Mouridsen K, Friston K, Hjort N, Gyldensted L, Ostergaard L, Kiebel S. Bayesian estimation of cerebral perfusion using a physiological model of microvasculature. NeuroImage 2006,33:570– 579.
  27. Kudo K, Boutelier T, Pautot F, et al. Bayesian analysis of perfusion-weighted imaging to predict infarct volume: comparison with singular value decomposition. Magn Reson Med Sci. 2014;13(1):45–50.
  28. Smit EJ, Prokop WM. Noise Reduction in Image Data. World Patent. 2017.
  29. Kudo K, Christensen S, Sasaki M, Ostergaard L, Shirato H, Ogasawara K, Wintermark M, Warach S. Accuracy and reliability assessment of CT and MR perfusion analysis software using a digital phantom. Radiology 2013,267:201-211.Sasaki M,
  30. Fabrice Pautot, K.Kudo, Timothé Boutelier, Makoto Sasaki. SNR and Acquisition Duration Explain Erroneous Mean Treansit Time in Acute Stroke CT Perfusion Imaging. in ISC. 2012. New Orleans.
  – Kudo, M. Sasaki, K. Yamada, S. Momoshima, H. Utsunomiya, H. Shirato, and K. Ogasawara. Differences in CT perfusion maps generated by different commercial software: Quantitative analysis by using identical source data of acute stroke patients, Radiology, vol. 254, no. 1, pp. 200–209, 2010

Download this PDF file

Contact Us