GPU-based RFA simulation for minimally invasive cancer treatment of liver tumours

PurposeRadiofrequency ablation (RFA) is one of the most popular and well-standardized minimally invasive cancer treatments (MICT) for liver tumours, employed where surgical resection has been contraindicated. Less-experienced interventional radiologists (IRs) require an appropriate planning tool for the treatment to help avoid incomplete treatment and so reduce the tumour recurrence risk. Although a few tools are available to predict the ablation lesion geometry, the process is computationally expensive. Also, in our implementation, a few patient-specific parameters are used to improve the accuracy of the lesion prediction.MethodsAdvanced heterogeneous computing using personal computers, incorporating the graphics processing unit (GPU) and the central processing unit (CPU), is proposed to predict the ablation lesion geometry. The most recent GPU technology is used to accelerate the finite element approximation of Penne’s bioheat equation and a three state cell model. Patient-specific input parameters are used in the bioheat model to improve accuracy of the predicted lesion.ResultsA fast GPU-based RFA solver is developed to predict the lesion by doing most of the computational tasks in the GPU, while reserving the CPU for concurrent tasks such as lesion extraction based on the heat deposition at each finite element node. The solver takes less than 3 min for a treatment duration of 26 min. When the model receives patient-specific input parameters, the deviation between real and predicted lesion is below 3 mm.ConclusionA multi-centre retrospective study indicates that the fast RFA solver is capable of providing the IR with the predicted lesion in the short time period before the intervention begins when the patient has been clinically prepared for the treatment.

[1]  Steve B. Jiang,et al.  GPU-based fast low-dose cone beam CT reconstruction via total variation. , 2010, Journal of X-Ray Science and Technology.

[2]  C.A. Cain,et al.  A cylindrical-section ultrasound phased-array applicator for hyperthermia cancer therapy , 1988, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[3]  Hui-Hong Liang,et al.  Efficacy of combining temperature- and power-controlled radiofrequency ablation for malignant liver tumors. , 2010, Chinese journal of cancer.

[4]  Karl Johan Åström,et al.  PID Controllers: Theory, Design, and Tuning , 1995 .

[5]  Andrea Schenk,et al.  Fast Estimation of the Vascular Cooling in RFA Based on Numerical Simulation , 2010 .

[6]  Christian Schumann,et al.  Interactive multi-criteria planning for radiofrequency ablation , 2015, International Journal of Computer Assisted Radiology and Surgery.

[7]  H. H. Penns Analysis of tissue and arterial blood temperatures in the resting human forearm , 1948 .

[8]  A. Jemal,et al.  Global cancer transitions according to the Human Development Index (2008-2030): a population-based study. , 2012, The Lancet. Oncology.

[9]  Chai Hong Yeong,et al.  Robot-assisted radiofrequency ablation of primary and secondary liver tumours: early experience , 2013, European Radiology.

[10]  Panagiotis Samaras,et al.  Quantitative Perfusion Analysis of Malignant Liver Tumors: Dynamic Computed Tomography and Contrast-Enhanced Ultrasound , 2012, Investigative radiology.

[11]  J. Stoer,et al.  Introduction to Numerical Analysis , 2002 .

[12]  Gerry Leversha,et al.  Introduction to numerical analysis (3rd edn), by J. Stoer and R. Bulirsch. Pp. 744. £49. 2002. ISBN 0 387 95452 X (Springer-Verlag). , 2004, The Mathematical Gazette.

[13]  H. H. Pennes Analysis of tissue and arterial blood temperatures in the resting human forearm. 1948. , 1948, Journal of applied physiology.

[14]  Luc Soler,et al.  Radiofrequency ablation of hepatic tumors: simulation, planning, and contribution of virtual reality and haptics , 2005, Computer methods in biomechanics and biomedical engineering.

[15]  S. Nahum Goldberg,et al.  Radiofrequency tumor ablation: principles and techniques , 2001 .

[16]  P. Wolf,et al.  Temperature-controlled and constant-power radio-frequency ablation: what affects lesion growth? , 1999, IEEE Transactions on Biomedical Engineering.

[17]  R Materne,et al.  Hepatic perfusion parameters in chronic liver disease: dynamic CT measurements correlated with disease severity. , 2001, AJR. American journal of roentgenology.

[18]  Steve B. Jiang,et al.  Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation , 2011, Physics in medicine and biology.

[19]  Hervé Delingette,et al.  Parameter Estimation for Personalization of Liver Tumor Radiofrequency Ablation , 2014, ABDI@MICCAI.

[20]  Hervé Delingette,et al.  Efficient Lattice Boltzmann Solver for Patient-Specific Radiofrequency Ablation of Hepatic Tumors , 2015, IEEE Transactions on Medical Imaging.

[21]  Tuomas Alhonnoro,et al.  Image-based multi-scale modelling and validation of radio-frequency ablation in liver tumours , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  ECR 2016 Book of Abstracts - C - Scientific and Educational Exhibits , 2016, Insights into Imaging.

[23]  Judith Muehl,et al.  Effects of needle placement inaccuracies in hepatic radiofrequency tumor ablation , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[24]  S. Goldberg,et al.  Radiofrequency tumor ablation: principles and techniques. , 2001, European journal of ultrasound : official journal of the European Federation of Societies for Ultrasound in Medicine and Biology.

[25]  Andrea Borsic,et al.  GPU-Accelerated real time simulation of Radio Frequency Ablation thermal dose , 2014, 2014 40th Annual Northeast Bioengineering Conference (NEBEC).

[26]  Martin A. Eisenberg,et al.  On finite element integration in natural co‐ordinates , 1973 .

[27]  J. Bruix,et al.  Management of hepatocellular carcinoma: An update , 2011, Hepatology.

[28]  Klaus H. Maier-Hein,et al.  The Medical Imaging Interaction Toolkit: challenges and advances , 2013, International Journal of Computer Assisted Radiology and Surgery.

[29]  Huang‐Wen Huang,et al.  Influence of blood vessel on the thermal lesion formation during radiofrequency ablation for liver tumors. , 2013, Medical physics.

[30]  Y Nagata,et al.  Clinical usefulness of determining the rate of thermal clearance within heated tumors. , 1996, Japanese journal of clinical oncology.

[31]  Heinz-Otto Peitgen,et al.  Numerical Simulation of Radio Frequency Ablation with State Dependent Material Parameters in Three Space Dimensions , 2006, MICCAI.

[32]  Christian Schumann,et al.  GPU-based Real-Time Approximation of the Ablation Zone for Radiofrequency Ablation , 2011, IEEE Transactions on Visualization and Computer Graphics.

[33]  Bernhard Kainz,et al.  Vessel Segmentation for Ablation Treatment Planning and Simulation , 2010, MICCAI.

[34]  M. Guevara-Valdivia,et al.  Radiofrequency ablation of multiple accessory pathways. , 2002, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[35]  Qing He,et al.  Transcatheter Arterial Chemoembolization: History for More than 30 Years , 2012, ISRN gastroenterology.

[36]  J Werner,et al.  Temperature profiles with respect to inhomogeneity and geometry of the human body. , 1988, Journal of applied physiology.

[37]  Tingying Peng,et al.  A Three-State Mathematical Model of Hyperthermic Cell Death , 2010, Annals of Biomedical Engineering.

[38]  Sheldon K Hall,et al.  Cell death, perfusion and electrical parameters are critical in models of hepatic radiofrequency ablation , 2015, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[39]  Max Grossman,et al.  Professional CUDA C Programming , 2014 .