Development of an in vitro setup for flow studies in a stented carotid artery bifurcation

To investigate flow conditions in a double-layered carotid artery stent, a bench-top in vitro flow setup including a bifurcation phantom was designed and fabricated. The geometry of the tissue-mimicking phantom was based on healthy individuals. Two identical phantoms were created using 3D-printing techniques and molding with PVA-gel. In one of them, a clinically available CGuard double-layer stent was inserted. Measurements were performed using both continuous and pulsatile flow conditions. Blood flow studies were performed using echoPIV: a novel ultrasound-based technique combined with particle image velocimetry. A maximum deviation of 3% was visible between desired and measured flow patterns. The echoPIV measurements showed promising results on visualization and quantification of blood flow in and downstream the stent. Further research could demonstrate the effects of a double-layered stent on blood flow patterns in a carotid bifurcation in detail.

Opportunities for personalised follow-up in breast cancer: the gap between daily practice and recurrence risk

PURPOSE

Follow-up guidelines barely diverge from a one-size-fits-all approach, even though the risk of recurrence differs per patient. However, the personalization of breast cancer care improves outcomes for patients. This study explores the variation in follow-up pathways in the Netherlands using real-world data to determine guideline adherence and the gap between daily practice and risk-based surveillance, to demonstrate the benefits of personalized risk-based surveillance compared with usual care.

METHODS

Patients with stage I-III invasive breast cancer who received surgical treatment in a general hospital between 2005 and 2020 were selected from the Netherlands Cancer Registry and included all imaging activities during follow-up from hospital-based electronic health records. Process analysis techniques were used to map patients and activities to investigate the real-world utilisation of resources and identify the opportunities for improvement. The INFLUENCE 2.0 nomogram was used for risk prediction of recurrence.

RESULTS

In the period between 2005 and 2020, 3478 patients were included with a mean follow-up of 4.9 years. In the first 12 months following treatment, patients visited the hospital between 1 and 5 times (mean 1.3, IQR 1-1) and received between 1 and 9 imaging activities (mean 1.7, IQR 1-2). Mammogram was the prevailing imaging modality, accounting for 70% of imaging activities. Patients with a low predicted risk of recurrence visited the hospital more often.

CONCLUSIONS

Deviations from the guideline were not in line with the risk of recurrence and revealed a large gap, indicating that it is hard for clinicians to accurately estimate this risk and therefore objective risk predictions could bridge this gap.

Dynamic Computed Tomography Angiography for capturing vessel wall motion: A phantom study for optimal image reconstruction

BACKGROUND

Reliably capturing sub-millimeter vessel wall motion over time, using dynamic Computed Tomography Angiography (4D CTA), might provide insight in biomechanical properties of these vessels. This may improve diagnosis, prognosis, and treatment decision making in vascular pathologies.

PURPOSE

The aim of this study is to determine the most suitable image reconstruction method for 4D CTA to accurately assess harmonic diameter changes of vessels.

METHODS

An elastic tube (inner diameter 6 mm, wall thickness 2 mm) was exposed to sinusoidal pressure waves with a frequency of 70 beats-per-minute. Five flow amplitudes were set, resulting in increasing sinusoidal diameter changes of the elastic tube, measured during three simulated pulsation cycles, using ECG-gated 4D CTA on a 320-detector row CT system. Tomographic images were reconstructed using one of the following three reconstruction methods: hybrid iterative (Hybrid-IR), model-based iterative (MBIR) and deep-learning based (DLR) reconstruction. The three reconstruction methods where based on 180 degrees (half reconstruction mode) and 360 degrees (full reconstruction mode) raw data. The diameter change, captured by 4D CTA, was computed based on image registration. As a reference metric for diameter change measurement, a 9 MHz linear ultrasound transducer was used. The sum of relative absolute differences (SRAD) between the ultrasound and 4D CTA measurements was calculated for each reconstruction method. The standard deviation was computed across the three pulsation cycles.

RESULTS

MBIR and DLR resulted in a decreased SRAD and standard deviation compared to Hybrid-IR. Full reconstruction mode resulted in a decreased SRAD and standard deviations, compared to half reconstruction mode.

CONCLUSIONS

4D CTA can capture a diameter change pattern comparable to the pattern captured by US. DLR and MBIR algorithms show more accurate results than Hybrid-IR. Reconstruction with DLR is >3 times faster, compared to reconstruction with MBIR. Full reconstruction mode is more accurate than half reconstruction mode.

Geometrical Changes of the Aorta as Predictors for Thromboembolic Events After EVAR With the Anaconda Stent-Graft

PURPOSE

Thromboembolic events (TE), including limb graft occlusion (LGO) and distal limb embolization (DLE), are common complications after endovascular aneurysm repair (EVAR). The aim of this study was to find predictors for TE in patients treated with the Anaconda stent-graft for infrarenal aneurysms.

MATERIALS AND METHODS

Geometrical and anatomical variables were retrospectively analyzed in a consecutive Anaconda cohort. Pre- and postoperative CT scans were used to derive geometrical parameters length, curvature, torsion, and tortuosity index (TI) from the center lumen lines (CLLs). Limb characteristics, pre-to-post EVAR and mid-term-follow-up changes in the parameters were evaluated for their predictive value for TE.

RESULTS

Eighty-four patients (mean age 74±8.3 years, 74 men) were enrolled. The risk of TE was lowered with pre-to-post implant decreasing TI (steps of 0.05: OR: 1.30, 95% CI: 1.01-1.66, p=0.04), pre-to-post implant decreasing mean curvature (OR: 1.08, 95% CI: 1.01-1.16, p=0.03), and a larger degree of circumferential common iliac artery (CIA) calcification (OR: 0.98, 95% CI: 0.97-1.00, p=0.03). The only LGO predictor was the caudal relocation of maximal curvature after EVAR (OR: 1.01, 95% CI: 1.00-1.01, p=0.04). Preventors of DLE were CIA diameter (OR: 0.87, 95% CI: 0.76-0.99, p=0.04), circumferential CIA calcification (OR: 0.97, 95% CI: 0.95-1.00, p=0.03), mean and maximal curvature of the preoperative aortoiliac trajectory (OR: 0.86, 95% CI: 0.79-0.94, p<0.01 and OR: 0.97, 95% CI: 0.95-1.00, p=0.03, respectively) and pre-to-postoperative decrease in mean curvature (OR: 1.11, 95% CI: 1.02-1.21, p=0.02). Midterm TE predictors were length (OR: 0.95, 95% CI: 0.89-1.01, p=0.08) and torsion maximum location (OR: 1.01, 95% CI: 0.99-1.01, p=0.10).

CONCLUSION

The present study confirms that treatment of infrarenal AAA with an Anaconda stent-graft is related to a relatively high TE rate which decreases with a pre-to-postoperative reduction in curvature and TI, and a larger degree of circumferential CIA calcification. In other words, more aortoiliac straightening and more circumferential CIA calcification may prevent TE development after EVAR with this stent-graft.

Renal and Visceral Artery Configuration During the First Year of Follow-Up After Fenestrated Aortic Aneurysm Repair Using the Anaconda Stent-graft: A Prospective Longitudinal Multicenter Study With ECG-Gated CTA Scans

OBJECTIVE

The performance of fenestrated endovascular aortic aneurysm repair (FEVAR) may be compromised by complications related to the dynamic vascular environment. The aim of this study was to analyze the behavior of FEVAR bridging stent configurations during the cardiac cycle and during follow-up to improve our understanding on treatment durability.

DESIGN

Twenty-one patients presenting with complex abdominal aortic aneurysms (AAAs; 9 juxtarenal/6 pararenal/3 paravisceral/1 thoracoabdominal aortic aneurysm type IV), treated with a fenestrated Anaconda (Terumo Aortic, Inchinnan, Scotland, UK) with Advanta V12 bridging stents (Getinge, Merrimack, NH, USA), were prospectively enrolled in a multicenter observational cohort study and underwent electrocardiogram (ECG)-gated computed tomographic angiography (CTA) preoperatively, at discharge, 7-week, and 12-month follow-ups.

METHODS

Fenestrated endovascular aortic aneurysm repair stability was assessed considering the following variables: branch angle as the angle between the aorta and the target artery, end-stent angle as the angle between the end of the bridging stent and the native artery downstream from it, curvature and tortuosity index (TI) to describe the bending of the target artery. Body-bridging stent stability was assessed considering bridging stent flare lengths, the distances between the proximal sealing stent-ring and fenestrations and the distance between the fenestration and first apposition in the target artery.

RESULTS

Renal branch angles significantly increased after FEVAR toward a perpendicular position (right renal artery from median 60.9°, inter quartile range [IQR]=44.2-84.9° preoperatively to 94.4°, IQR=72.6-99.8°, p=0.001 at 12-month follow-up; left renal artery [LRA], from 63.7°, IQR=55.0-73.0° to 94.3°, IQR=68.2-105.6°, p<0.001), while visceral branch angles did not. The mean dynamic curvature only decreased for the LRA from preoperative (3.0, IQR=2.2-3.8 m-1) to 12-month follow-up (1.9, IQR=1.4-2.6 m-1, p=0.027). The remaining investigated variables did not seem to show any changes over time in this cohort.

CONCLUSIONS

Fenestrated endovascular aortic aneurysm repair for complex AAAs using the Anaconda fenestrated stent-graft and balloon-expandable Advanta V12 bridging stents demonstrated stable configurations up to 12-month follow-up, except for increasing renal branch angles toward perpendicular orientation to the aorta, yet without apparent clinical consequences in this cohort.

CLINICAL IMPACT

This study provides detailed information on the cardiac-pulsatility-induced (dynamic) and longitudinal geometry deformations of the target arteries and bridging stents after fenestrated endovascular aortic aneurysm repair (FEVAR) up to 12-month follow-up. The configuration demonstrated limited dynamic and longitudinal deformations in terms of branch angle, end-stent angle, curvature, and tortuosity index (TI), except for the increasing renal branch angles that go toward a perpendicular orientation to the aorta. Overall, the results suggest that the investigated FEVAR configurations are stable and durable, though careful consideration of increasing renal branch angles and significant geometry alterations is advised.

The influence of electrocardiogram-gated computed tomography reconstruction into 8 or 10 cardiac phases on cardiac-pulsatility-induced motion quantification of stent grafts in the aorta

Objective

The goal of this study was to determine to what extent aortic stent graft motion quantification is comparable between electrocardiogram (ECG)-gated computed tomography (CT) scans with reconstructions into 8 and 10 cardiac phases on CT scanners from two different vendors.

Methods

An experimental setup that induces motion of an aortic stent graft, according to a predefined aortic blood pressure wave, was placed in two CT scanners of different vendors. The stent graft motion was captured using an ECG-gated CT technique and quantified using dedicated analysis algorithms. The calculated motion amplitudes and total traveled path lengths of stent segmentations were compared between scans reconstructed into 8 and 10 phases and between the scanners, after validation with sensor measurements and repeated measurements.

Results

No difference in motion amplitudes in z-direction (craniocaudal direction) was observed between the reconstructions into 8 and 10 phases (0.02 mm; 95% confidence interval [CI], -0.01 to 0.05 mm; P = .358). The z-amplitudes differed by 0.04 mm (95% CI, 0.01-0.07 mm; P = .003) between the different CT scanners. Path lengths differed 0.07 mm (95% CI, 0.01-to 0.13 mm; P = .013) between the reconstructions into 8 and 10 phases and 0.13 mm (95% CI, 0.06-0.17 mm; P < .001) between the different scanners.

Conclusions

The motion amplitudes can accurately be compared between 8 and 10 phases and between the two scanners, without differences larger than the voxel size of 0.3 × 0.3 × 0.5 mm. Clinical motion analysis results of different ECG-gated CT scans and CT scanners can be compared up to the accuracy of the CT scan.

The Effect of Partial Electrical Insulation of the Tip and Active Needle Length of Monopolar Irreversible Electroporation Electrodes on the Electric Field Line Pattern and Temperature Gradient to Improve Treatment Control

Unintentional local temperature effects can occur during irreversible electroporation (IRE) treatment, especially near the electrodes, and most frequently near the tip. Partial electrical insulation of the IRE electrodes could possibly control these temperature effects. This study investigated and visualized the effect of partial electrical insulation applied to the IRE electrodes on the electric field line pattern and temperature gradient. Six designs of (partial) electrical insulation of the electrode tip and/or active needle length (ANL) of the original monopolar 19G IRE electrodes were investigated. A semolina in castor oil model was used to visualize the electric field line pattern in a high-voltage static electric field. An optical method to visualize a change in temperature gradient (color Schlieren) was used to image the temperature development in a polyacrylamide gel. Computational models were used to support the experimental findings. Around the electrode tip, the highest electric field line density and temperature gradient were present. The more insulation was applied to the electrodes, the higher the resistance. Tip and ANL insulation together reduced the active area of and around the electrodes, resulting in a visually enlarged area that showed a change in temperature gradient. Electrically insulating the electrode tip together with an adjustment in IRE parameter settings could potentially reduce the uncontrollable influence of the tip and may improve the predictability of the current pathway development.

Machine learning based model to diagnose obstructive coronary artery disease using calcium scoring, PET imaging, and clinical data

INTRODUCTION

Accurate risk stratification in patients with suspected stable coronary artery disease is essential for choosing an appropriate treatment strategy. Our aim was to develop and validate a machine learning (ML) based model to diagnose obstructive CAD (oCAD).

METHOD

We retrospectively have included 1007 patients without a prior history of CAD who underwent CT-based calcium scoring (CACS) and a Rubidium-82 PET scan. The entire dataset was split 4:1 into a training and test dataset. An ML model was developed on the training set using fivefold stratified cross-validation. The test dataset was used to compare the performance of expert readers to the model. The primary endpoint was oCAD on invasive coronary angiography (ICA).

RESULTS

ROC curve analysis showed an AUC of 0.92 (95% CI 0.90-0.94) for the training dataset and 0.89 (95% CI 0.84-0.93) for the test dataset. The ML model showed no significant differences as compared to the expert readers (p ≥ 0.03) in accuracy (89% vs. 88%), sensitivity (68% vs. 69%), and specificity (92% vs. 90%).

CONCLUSION

The ML model resulted in a similar diagnostic performance as compared to expert readers, and may be deployed as a risk stratification tool for obstructive CAD. This study showed that utilization of ML is promising in the diagnosis of obstructive CAD.

Preoperative imaging in primary hyperparathyroidism patients using 4DCT subtraction maps, a report of three cases

Four-dimensional computed tomography (4DCT) is one of the preoperative imaging modalities that can be used to localize a parathyroid adenoma in primary hyperparathyroidism patients however, sensitivity differs in literature and could be improved especially for multiglandular hyperplasia or double adenomas. The most robust feature on the 4DCT for the differentiation between parathyroid adenoma and thyroid gland tissue is arterial enhancement. To make this better visible, we have developed a subtraction map that shows arterial enhancement as a color scale to increase sensitivity for 4DCT. In this report of 3 cases, we present the usefulness of this subtraction map in a 54-year-old male, a 57-year-old female and a 51-year-old male. Subtraction maps may increase sensitivity for 4DCT, especially for multiglandular hyperplasia or double adenomas.

Development and External Validation of a PET Radiomic Model for Prognostication of Head and Neck Cancer

AIM

To build and externally validate an [¹⁸F]FDG PET radiomic model to predict overall survival in patients with head and neck squamous cell carcinoma (HNSCC).

METHODS

Two multicentre datasets of patients with operable HNSCC treated with preoperative afatinib who underwent a baseline and evaluation [¹⁸F]FDG PET/CT scan were included (EORTC: n = 20, Unicancer: n = 34). Tumours were delineated, and radiomic features were extracted. Each cohort served once as a training and once as an external validation set for the prediction of overall survival. Supervised feature selection was performed using variable hunting with variable importance, selecting the top two features. A Cox proportional hazards regression model using selected radiomic features and clinical characteristics was fitted on the training dataset and validated in the external validation set. Model performances are expressed by the concordance index (C-index).

RESULTS

In both models, the radiomic model surpassed the clinical model with validation C-indices of 0.69 and 0.79 vs. 0.60 and 0.67, respectively. The model that combined the radiomic features and clinical variables performed best, with validation C-indices of 0.71 and 0.82.

CONCLUSION

Although assessed in two small but independent cohorts, an [¹⁸F]FDG-PET radiomic signature based on the evaluation scan seems promising for the prediction of overall survival for HNSSC treated with preoperative afatinib. The robustness and clinical applicability of this radiomic signature should be assessed in a larger cohort.

The Influence of Irreversible Electroporation Parameters on the Size of the Ablation Zone and Thermal Effects: A Systematic Review

Introduction: The aim of this study was to review the effect of irreversible electroporation parameter settings on the size of the ablation zone and the occurrence of thermal effects. This insight would help to optimize treatment protocols and effectively ablate a tumor while controlling the occurrence of thermal effects. Methods: Various individual studies report the influence of variation in electroporation parameters on the ablation zone size or occurrence of thermal effects. However, no connections have yet been established between these studies. With the aim of closing the gap in the understanding of and personalizing irreversible electroporation parameter settings, a systematic review was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. A quality assessment was performed using an in-house developed grading tool based on components of commonly used grading domains. Data on the electroporation parameters voltage, number of electrodes, inter-electrode distance, active needle length, pulse length/number/protocol/frequency, and pulse interval were extracted. Ablation zone size and temperature data were grouped per parameter. Spearman correlation and linear regression were used to define the correlation with outcome measures. Results: A total of 7661 articles were screened, of which 18 preclinical studies (animal and phantom studies) met the inclusion criteria. These studies were graded as moderate (4/18) and low (14/18) quality. Only the applied voltage appeared to be a significant linear predictor of ablation zone size: length, surface, and volume. The pulse number was moderately but nonlinearly correlated with the ablation zone length. Thermal effects were more likely to occur for higher voltages (≥2000 V), higher number of electrodes, and increased active needle length. Conclusion: Firm conclusions are limited since studies that investigated and precisely reported the influence of electroporation parameters on the ablation zone size and thermal effects were scarce and mostly graded low quality. High-quality studies are needed to improve the predictability of the combined effect of variation in parameter combinations and optimize irreversible electroporation treatment protocols.

Unsupervised convolutional autoencoders for 4D transperineal ultrasound classification

Purpose

4D Transperineal ultrasound (TPUS) is used to examine female pelvic floor disorders. Muscle movement, like performing a muscle contraction or a Valsalva maneuver, can be captured on TPUS. Our work investigates the possibility for unsupervised analysis and classification of the TPUS data.

Approach

An unsupervised 3D-convolutional autoencoder is trained to compress TPUS volume frames into a latent feature vector (LFV) of 128 elements. The (co)variance of the features are analyzed and statistical tests are performed to analyze how features contribute in storing contraction and Valsalva information. Further dimensionality reduction is applied (principal component analysis or a 2D-convolutional autoencoder) to the LFVs of the frames of the TPUS movie to compress the data and analyze the interframe movement. Clustering algorithms ( K -means clustering and Gaussian mixture models) are applied to this representation of the data to investigate the possibilities of unsupervised classification.

Results

The majority of the features show a significant difference between contraction and Valsalva. The (co)variance of the features from the LFVs was investigated and features most prominent in capturing muscle movement were identified. Furthermore, the first principal component of the frames from a single TPUS movie can be used to identify movement between the frames. The best classification results were obtained after applying principal component analysis and Gaussian mixture models to the LFVs of the TPUS movies, yielding a 91.2% accuracy.

Conclusion

Unsupervised analysis and classification of TPUS data yields relevant information about the type and amount of muscle movement present.

[18F]FDG-PET/CT radiomics for the identification of genetic clusters in pheochromocytomas and paragangliomas

Objectives

Based on germline and somatic mutation profiles, pheochromocytomas and paragangliomas (PPGLs) can be classified into different clusters. We investigated the use of [¹⁸F]FDG-PET/CT radiomics, SUV_max and biochemical profile for the identification of the genetic clusters of PPGLs.

Methods

In this single-centre cohort, 40 PPGLs (13 cluster 1, 18 cluster 2, 9 sporadic) were delineated using a 41% adaptive threshold of SUV_peak ([¹⁸F]FDG-PET) and manually (low-dose CT; ldCT). Using PyRadiomics, 211 radiomic features were extracted. Stratified 5-fold cross-validation for the identification of the genetic cluster was performed using multinomial logistic regression with dimensionality reduction incorporated per fold. Classification performances of biochemistry, SUV_max and PET(/CT) radiomic models were compared and presented as mean (multiclass) test AUCs over the five folds. Results were validated using a sham experiment, randomly shuffling the outcome labels.

Results

The model with biochemistry only could identify the genetic cluster (multiclass AUC 0.60). The three-factor PET model had the best classification performance (multiclass AUC 0.88). A simplified model with only SUV_max performed almost similarly. Addition of ldCT features and biochemistry decreased the classification performances. All sham AUCs were approximately 0.50.

Conclusion

PET radiomics achieves a better identification of PPGLs compared to biochemistry, SUV_max, ldCT radiomics and combined approaches, especially for the differentiation of sporadic PPGLs. Nevertheless, a model with SUV_max alone might be preferred clinically, weighing model performances against laborious radiomic analysis. The limited added value of radiomics to the overall classification performance for PPGL should be validated in a larger external cohort.

Key Points

• Radiomics derived from [¹⁸F]FDG-PET/CT has the potential to improve the identification of the genetic clusters of pheochromocytomas and paragangliomas.

• A simplified model with SUV_maxonly might be preferred clinically, weighing model performances against the laborious radiomic analysis.

• Cluster 1 and 2 PPGLs generally present distinctive characteristics that can be captured using [¹⁸F]FDG-PET imaging. Sporadic PPGLs appear more heterogeneous, frequently resembling cluster 2 PPGLs and occasionally resembling cluster 1 PPGLs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-09034-5.

Effect of temporal sampling protocols on myocardial blood flow measurements using Rubidium-82 PET

BACKGROUND

A variety of temporal sampling protocols is used worldwide to measure myocardial blood flow (MBF). Both the length and number of time frames in these protocols may alter MBF and myocardial flow reserve (MFR) measurements. We aimed to assess the effect of different clinically used temporal sampling protocols on MBF and MFR quantification in Rubidium-82 (Rb-82) PET imaging.

METHODS

We retrospectively included 20 patients referred for myocardial perfusion imaging using Rb-82 PET. A literature search was performed to identify appropriate sampling protocols. PET data were reconstructed using 14 selected temporal sampling protocols with time frames of 5-10 seconds in the first-pass phase and 30-120 seconds in the tissue phase. Rest and stress MBF and MFR were calculated for all protocols and compared to the reference protocol with 26 time frames.

RESULTS

MBF measurements differed (P ≤ 0.003) in six (43%) protocols in comparison to the reference protocol, with mean absolute relative differences up to 16% (range 5%-31%). Statistically significant differences were most frequently found for protocols with tissue phase time frames < 90 seconds. MFR did not differ (P ≥ 0.11) for any of the protocols.

CONCLUSIONS

Various temporal sampling protocols result in different MBF values using Rb-82 PET. MFR measurements were more robust to different temporal sampling protocols.

Evaluation of electrocardiogram-gated computed tomography angiography to quantify changes in geometry and dynamic behavior of the iliac artery after placement of the Gore Excluder Iliac Branch Endoprosthesis

BACKGROUND

The GORE^® EXCLUDER^® Iliac Branch Endoprosthesis (IBE) is designed to treat iliac aneurysms with preservation of blood flow through the internal iliac artery (IIA). Little is known about the influence of IBE placement on the IIA geometry. This study aimed to provide detailed insights in the dynamic behavior and geometry of the common iliac artery (CIA) and IIA trajectory and how these are influenced after treatment with an IBE.

METHODS

Pre- and postoperative electrocardiogram-gated computed tomography angiography (ECG-gated CTA) scans were acquired in a prospective study design and analyzed with in-house written algorithms designed for aorto-iliac and endoprosthesis deformation evaluation. Cardiac pulsatility-induced motion patterns and pathlengths were computed by tracking predefined locations on the aorto-iliac tract. Centerlines through the CIA-IIA trajectory were used to investigate the static and dynamic geometry, including curvature, torsion, length and Tortuosity Index (TI).

RESULTS

Fourteen CIA-IIA trajectories were analyzed before and after IBE placement. Cardiac pulsatility-induced motion and pathlengths increased after IBE placement, especially at mid IIA and the first IIA bifurcation (P≤0.04). After IBE placement, static and dynamic curvature, length and TI decreased significantly (P<0.05). Furthermore, the average dynamic torsion increased significantly (P=0.030). The remaining geometrical outcomes were not statistically significant.

CONCLUSIONS

The placement of an IBE device stiffens and straightens the CIA-IIA trajectory. Its relation with clinical outcome is yet to be investigated, which can be done thoroughly with the ECG-gated CTA algorithms used in this study.

Predicting Aesthetic Outcome of the Nuss Procedure in Patients with Pectus Excavatum

Patients suffering from pectus excavatum often experience psychosocial distress due to perceived anomalies in their physical appearance. The ability to visually inform patients about their expected aesthetic outcome after surgical correction is still lacking. This study aims to develop an automatic, patient-specific model to predict aesthetic outcome after the Nuss procedure. Patients prospectively received preoperative and postoperative 3-dimensional optical surface scanning of their chest during the Nuss procedure. A prediction model was composed based on nonlinear least squares data-fitting, regression methods and a 2-dimensional Gaussian function with adjustable amplitude, variance, rotation, skewness, and kurtosis components. Morphological features of pectus excavatum were extracted from preoperative images using a previously developed surface analysis tool to generate a patient-specific model. Prediction accuracy was evaluated through cross-validation, utilizing the mean root squared deviation and maximum positive and negative deviations as performance measures. The prediction model was evaluated on 30 (90% male) prospectively imaged patients. The model achieved an average root mean squared deviation of 6.3 ± 2.0 mm, with average maximum positive and negative deviations of 12.7 ± 6.1 and -10.2 ± 5.7 mm, respectively, between the predicted and actual postoperative aesthetic result. Our developed 2-dimensional Gaussian model based on 3-dimensional optical surface images is a clinically promising tool to predict postsurgical aesthetic outcome in patients with pectus excavatum. Prediction of the aesthetic outcome after the Nuss procedure potentially improves information provision and expectation management among patients. Further research should assess whether increasing the sample size may reduce deviations and improve performance.

Development of a dynamic myocardial perfusion phantom model for tracer kinetic measurements

BACKGROUND

Absolute myocardial perfusion imaging (MPI) is beneficial in the diagnosis and prognosis of patients with suspected or known coronary artery disease. However, validation and standardization of perfusion estimates across centers is needed to ensure safe and adequate integration into the clinical workflow. Physical myocardial perfusion models can contribute to this clinical need as these can provide ground-truth validation of perfusion estimates in a simplified, though controlled setup. This work presents the design and realization of such a myocardial perfusion phantom and highlights initial performance testing of the overall phantom setup using dynamic single photon emission computed tomography.

RESULTS

Due to anatomical and (patho-)physiological representation in the 3D printed myocardial perfusion phantom, we were able to acquire 22 dynamic MPI datasets in which ^99mTc-labelled tracer kinetics was measured and analyzed using clinical MPI software. After phantom setup optimization, time activity curve analysis was executed for measurements with normal myocardial perfusion settings (1.5 mL/g/min) and with settings containing a regional or global perfusion deficit (0.8 mL/g/min). In these measurements, a specific amount of activated carbon was used to adsorb radiotracer in the simulated myocardial tissue. Such mimicking of myocardial tracer uptake and retention over time satisfactorily matched patient tracer kinetics. For normal perfusion levels, the absolute mean error between computed myocardial blood flow and ground-truth flow settings ranged between 0.1 and 0.4 mL/g/min.

CONCLUSION

The presented myocardial perfusion phantom is a first step toward ground-truth validation of multimodal, absolute MPI applications in the clinical setting. Its dedicated and 3D printed design enables tracer kinetic measurement, including time activity curve and potentially compartmental myocardial blood flow analysis.

In Vivo Quantification of Cardiac-Pulsatility-Induced Motion Before and After Double-Branched Endovascular Aortic Arch Repair

The Relay^®Branch stent-graft (Terumo Aortic, Sunrise, FL, USA) offers a custom-made endovascular solution for complex aortic arch pathologies. In this technical note, a modified electrocardiography (ECG)-gated computed tomography (CT)-based algorithm was applied to quantify cardiac-pulsatility-induced changes of the aortic arch geometry and motion before and after double-branched endovascular repair (bTEVAR) of an aortic arch aneurysm. This software algorithm has the potential to provide novel and clinically relevant insights in the influence of bTEVAR on aortic anatomy, arterial compliance, and stent-graft dynamics.

Systematic Review of Health Economic Evaluations Focused on Artificial Intelligence in Healthcare: The Tortoise and the Cheetah

OBJECTIVES

This study aimed to systematically review recent health economic evaluations (HEEs) of artificial intelligence (AI) applications in healthcare. The aim was to discuss pertinent methods, reporting quality and challenges for future implementation of AI in healthcare, and additionally advise future HEEs.

METHODS

A systematic literature review was conducted in 2 databases (PubMed and Scopus) for articles published in the last 5 years. Two reviewers performed independent screening, full-text inclusion, data extraction, and appraisal. The Consolidated Health Economic Evaluation Reporting Standards and Philips checklist were used for the quality assessment of included studies.

RESULTS

A total of 884 unique studies were identified; 20 were included for full-text review, covering a wide range of medical specialties and care pathway phases. The most commonly evaluated type of AI was automated medical image analysis models (n = 9, 45%). The prevailing health economic analysis was cost minimization (n = 8, 40%) with the costs saved per case as preferred outcome measure. A total of 9 studies (45%) reported model-based HEEs, 4 of which applied a time horizon >1 year. The evidence supporting the chosen analytical methods, assessment of uncertainty, and model structures was underreported. The reporting quality of the articles was moderate as on average studies reported on 66% of Consolidated Health Economic Evaluation Reporting Standards items.

CONCLUSIONS

HEEs of AI in healthcare are limited and often focus on costs rather than health impact. Surprisingly, model-based long-term evaluations are just as uncommon as model-based short-term evaluations. Consequently, insight into the actual benefits offered by AI is lagging behind current technological developments.

In Vitro Geometry Analysis of Fenestrations in Endovascular Aneurysm Repair

PURPOSE

Changes in the flared end of balloon-expandable covered stent (BECS) may precede BECS-associated complications but are not regularly assessed with computed tomographic angiography (CTA) after fenestrated endovascular aneurysm repair (FEVAR). Validation of the flare geometric analysis (FGA) and assessment of intraobserver and interobserver variability are investigated in this study.

METHODS

Two series of 3 BeGraft BECSs (Bentley InnoMed GmbH, Hechingen, Germany) and 1 series of 3 Advanta V12 BECSs (Getinge AB, Göteborg, Sweden) were deployed in 3 side branches (45°, 60°, and 90° aortic branch angles) of an aorta phantom model. A standard post-FEVAR CTA scan was acquired. Computed tomographic angiography-derived measurements consisted of centerline reconstructions and placement of 3-dimensional coordinate markers by 2 observers in a vascular workstation. Flare geometric analysis calculates 3 BECS parameters: the circumferential flare-to-fenestration distance (FFD), which is the distance from the proximal end of the flare to fenestration, and diameters at the proximal end of the flare (Dflare) and at the fenestration (Dfenestration). Computed tomographic angiography-derived measurements were validated against microscopy measurements. Bland-Altman plots were used to determine the intraobserver and interobserver variability of the BECS parameters and intraclass correlation coefficient (ICC).

RESULTS

For each BECS, the FFD at 4 equidistant quadrants of the circumference, Dflare, and Dfenestration were calculated. The mean difference and repeatability coefficient (RC) of the validation were 0.8 (2.1) mm for FFD, 0.4 (1.0) mm for Dflare, and -0.2 (1.2) mm for Dfenestration. The mean intraobserver and interobserver difference (RC) was 0.5 (1.6) mm and 0.7 (2.6) mm for FFD, 0.1 (0.6) mm and 0.1 (0.7) mm for Dflare, and -0.1 (0.8) mm and -0.8 (1.0) mm for Dfenestration. The mean ICC of intraobserver variability was 0.86 for FFD, 0.94 for Dflare, and 0.78 for Dfenestration. The mean ICC of interobserver variability was 0.77 for FFD, 0.92 for Dflare, and 0.48 for Dfenestration.

CONCLUSION

This study showed that FGA of the flared ends of BECS can be performed with high accuracy in a phantom model, with good intraobserver and interobserver variability. Flare geometric analysis can be used to determine flare geometry of the BECS on standard post-FEVAR CTA scans.

Design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomic quantification in hybrid imaging; a feasibility study

Background

Accuracy and precision assessment in radiomic features is important for the determination of their potential to characterize cancer lesions. In this regard, simulation of different imaging conditions using specialized phantoms is increasingly being investigated. In this study, the design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomics quantification in hybrid imaging is presented.

Methods

A modular multimodality imaging phantom was constructed that could simulate different patterns of heterogeneous uptake and enhancement patterns in positron emission tomography (PET), single‐photon emission computed tomography (SPECT), computed tomography (CT), and magnetic resonance (MR) imaging. The phantom was designed to be used as an insert in the standard NEMA‐NU2 IEC body phantom casing. The entire phantom insert is composed of three segments, each containing three separately fillable compartments. The fillable compartments between segments had different sizes in order to simulate heterogeneous patterns at different spatial scales. The compartments were separately filled with different ratios of ^99mTc‐pertechnetate, ¹⁸F‐fluorodeoxyglucose ([¹⁸F]FDG), iodine‐ and gadolinium‐based contrast agents for SPECT, PET, CT, and T₁‐weighted MR imaging respectively. Image acquisition was performed using standard oncological protocols on all modalities and repeated five times for repeatability assessment. A total of 93 radiomic features were calculated. Variability was assessed by determining the coefficient of quartile variation (CQV) of the features. Comparison of feature repeatability at different modalities and spatial scales was performed using Kruskal‐Wallis‐, Mann‐Whitney U‐, one‐way ANOVA‐ and independent t‐tests.

Results

Heterogeneous uptake and enhancement could be simulated on all four imaging modalities. Radiomic features in SPECT were significantly less stable than in all other modalities. Features in PET were significantly less stable than in MR and CT. A total of 20 features, particularly in the gray‐level co‐occurrence matrix (GLCM) and gray‐level run‐length matrix (GLRLM) class, were found to be relatively stable in all four modalities for all three spatial scales of heterogeneous patterns (with CQV < 10%).

Conclusion

The phantom was suitable for simulating heterogeneous uptake and enhancement patterns in [¹⁸F]FDG‐PET, ^99mTc‐SPECT, CT, and T₁‐weighted MR images. The results of this work indicate that the phantom might be useful for the further development and optimization of imaging protocols for radiomic quantification in hybrid imaging modalities.

Value of SiPM PET in myocardial perfusion imaging using Rubidium-82

BACKGROUND

PET scanners using silicon photomultipliers with digital readout (SiPM PET) have an improved temporal and spatial resolution compared to PET scanners using conventional photomultiplier tubes (PMT PET). However, the effect on image quality and visibility of perfusion defects in myocardial perfusion imaging (MPI) is unknown. Our aim was to determine the value of a SiPM PET scanner in MPI.

METHODS

We prospectively included 30 patients who underwent rest and regadenoson-induced stress Rubidium-82 (Rb-82) MPI on the D690 PMT PET (GE Healthcare) and within three weeks on the Vereos SiPM PET (Philips Healthcare). Two expert readers scored the image quality and assessed the existence of possible defects. In addition, interpreter's confidence, myocardial blood flow (MBF), and myocardial flow reserve (MFR) values were compared.

RESULTS

Image quality improved (P = 0.03) using the Vereos as compared to the D690. Image quality of the Vereos and the D690 was graded fair in 20% and 10%, good in 60% and 50%, and excellent in 20% and 40%, respectively. Defect interpretation and interpreter's confidence did not differ between the D690 and the Vereos (P > 0.50). There were no significant differences in rest MBF (P ≥ 0.29), stress MBF (P ≥ 0.11), and MFR (P ≥ 0.51).

CONCLUSION

SiPM PET provides an improved image quality in comparison with PMT PET. Defect interpretation, interpreter's confidence, and absolute blood flow measurements were comparable between both systems. SiPM PET is therefore a reliable technique for MPI using Rb-82.

TRIAL REGISTRATION

ToetsingOnline NL63853.075.17. Registered 13 November, 2017.

Development of a dedicated 3D printed myocardial perfusion phantom: proof-of-concept in dynamic SPECT

We aim to facilitate phantom-based (ground truth) evaluation of dynamic, quantitative myocardial perfusion imaging (MPI) applications. Current MPI phantoms are static representations or lack clinical hard- and software evaluation capabilities. This proof-of-concept study demonstrates the design, realisation and testing of a dedicated cardiac flow phantom. The 3D printed phantom mimics flow through a left ventricular cavity (LVC) and three myocardial segments. In the accompanying fluid circuit, tap water is pumped through the LVC and thereafter partially directed to the segments using adjustable resistances. Regulation hereof mimics perfusion deficit, whereby flow sensors serve as reference standard. Seven phantom measurements were performed while varying injected activity of ^99mTc-tetrofosmin (330–550 MBq), cardiac output (1.5–3.0 L/min) and myocardial segmental flows (50–150 mL/min). Image data from dynamic single photon emission computed tomography was analysed with clinical software. Derived time activity curves were reproducible, showing logical trends regarding selected input variables. A promising correlation was found between software computed myocardial flows and its reference (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document}ρ= − 0.98; p = 0.003). This proof-of-concept paper demonstrates we have successfully measured first-pass LV flow and myocardial perfusion in SPECT-MPI using a novel, dedicated, myocardial perfusion phantom.

Strength testing of low-cost 3D-printed transtibial prosthetic socket

Measurement and production of traditional prosthetic sockets are time-consuming, labor-intensive, and highly dependent on the personnel involved. An alternative way to make prostheses is using computer-aided design (CAD) and computer-aided manufacturing (CAM). Fused Filament Fabrication (FFF) may be an alternative to make low-cost prosthetic sockets. This study investigates the tensile properties of potential printing materials suitable for FFF according to ISO527 (Standard Test Method for Tensile Properties of Plastics). To ensure that FFF-printed sockets are safe for patient usage, the structural integrity of the 3D-printed prosthesis will be investigated according to ISO10328 (International Standard Structural Testing of Lower Limb Prostheses). Tough PLA was the most suitable print material according to ISO 527 testing. The Tough PLA printed socket completed 2.27 million cycles and a static test target value of 4025 N. Future research remains necessary to continue testing new potential materials, improve print settings, and improve the socket design for the production of FFF-printed transtibial prosthetic sockets. FFF using Tough PLA can be used to create transtibial prostheses that almost comply with the International Standard for Structural Testing of Lower Limb Prostheses.

Hemodynamic changes after intracranial aneurysm growth

OBJECTIVE

For accurate risk assessment of unruptured intracranial aneurysms, it is important to understand the underlying mechanisms that lead to rupture. It is known that hemodynamic anomalies contribute to aneurysm growth and rupture, and that growing aneurysms carry higher rupture risks. However, it is unknown how growth affects hemodynamic characteristics. In this study, the authors assessed how hemodynamic characteristics change over the course of aneurysm growth.

METHODS

The authors included patients with observed aneurysm growth on longitudinal MRA in the period between 2012 and 2016. Patient-specific vascular models were created from baseline and follow-up images. Subsequently, intraaneurysmal hemodynamic characteristics were computed using computational fluid dynamics. The authors computed the normalized wall shear stress, oscillatory shear index, and low shear area to quantify hemodynamic characteristics. Differences between baseline and follow-up measurements were analyzed using paired t-tests.

RESULTS

Twenty-five patients with a total of 31 aneurysms were included. The aneurysm volume increased by a median (IQR) of 26 (9-39) mm3 after a mean follow-up period of 4 (range 0.4-10.9) years. The median wall shear stress decreased significantly after growth. Other hemodynamic parameters did not change significantly, although large individual changes with large variability were observed.

CONCLUSIONS

Hemodynamic characteristics change considerably after aneurysm growth. On average, wall shear stress values decrease after growth, but there is a large variability in hemodynamic changes between aneurysms.

Target vessel displacement during fenestrated and branched endovascular aortic repair and its implications for the role of traditional computed tomography angiography roadmaps

Background

This retrospective study quantifies target vessel displacement during fenestrated and branched endovascular aneurysm repair due to the introduction of stiff guidewires and stent graft delivery systems. The effect that intraoperative vessel displacement has on the usability of computed tomography angiography (CTA) roadmaps is also addressed.

Methods

Patients that underwent fenestrated or branched EVAR were included in this retrospective study. Two imaging datasets were collected from each patient: (I) preoperative CTA and (II) intraoperative contrast-enhanced cone beam computed tomography (ceCBCT) acquired after the insertion of the stiff guidewire and stent graft delivery system. After image registration, the 3D coordinates of the ostium of the celiac artery, superior mesenteric artery, right renal artery and left renal artery were recorded in both the CTA and the ceCBCT dataset by two observers. The three-dimensional displacement of the ostia of the target vessels was calculated by subtracting the coordinates of CTA and ceCBCT from one another. Additionally, the tortuosity index and the maximum angulation of the aorta were calculated.

Results

In total 20 patients and 77 target vessels were included in this study. The ostium of the celiac, superior mesenteric, right renal and left renal artery underwent non-uniform three-dimensional displacement with mean absolute displacement of 8.2, 7.7, 8.2 and 6.2 mm, respectively. The average displacement of all different target vessels together was 7.8 mm. A moderate correlation between vessel displacement and the maximum angulation of the aortoiliac segment was found (Spearman's ρ=0.45, P<0.05).

Conclusions

The introduction of stiff endovascular devices during fenestrated or branched EVAR causes significant, non-uniform displacement of the ostium of the visceral and renal target vessels. Consequently, preoperative CTA roadmaps based on bone registration are suboptimal to guide target vessel catheterization during these procedures.

Assessment of Dynamic Change of Coronary Artery Geometry and Its Relationship to Coronary Artery Disease, Based on Coronary CT Angiography

To investigate the relationship between dynamic changes of coronary artery geometry and coronary artery disease (CAD) using computed tomography (CT). Seventy-one patients underwent coronary CT angiography with retrospective electrocardiographic gating. End-systolic (ES) and end-diastolic (ED) phases were automatically determined by dedicated software. Centerlines were extracted for the right and left coronary artery. Differences between ES and ED curvature and tortuosity were determined. Associations of change in geometrical parameters with plaque types and degree of stenosis were investigated using linear mixed models. The differences in number of inflection points were analyzed using Wilcoxon signed-rank tests. Tests were done on artery and segment level. One hundred thirty-seven arteries (64.3%) and 456 (71.4%) segments were included. Curvature was significantly higher in ES than in ED phase for arteries (p = 0.002) and segments (p < 0.001). The difference was significant only at segment level for tortuosity (p = 0.005). Number of inflection points was significantly higher in ES phase on both artery and segment level (p < 0.001). No significant relationships were found between degree of stenosis and plaque types and dynamic change in geometrical parameters. Non-invasive imaging by cardiac CT can quantify change in geometrical parameters of the coronary arteries during the cardiac cycle. Dynamic change of vessel geometry through the cardiac cycle was not found to be related to the presence of CAD.

Body weight-dependent Rubidium-82 activity results in constant image quality in myocardial perfusion imaging with PET

BACKGROUND

Clinical practice shows degrading image quality in heavier patients who undergo myocardial perfusion imaging (MPI) with Rubidium-82 (Rb-82) PET when using a fixed tracer activity. Our aim was to derive and validate a patient-specific activity protocol resulting in a constant image quality in PET MPI.

METHODS

We included 251 patients who underwent rest MPI with Rb-82 PET (Discovery 670, GE Healthcare). 132 patients were included retrospectively and were scanned using a fixed activity of 740 MBq. The total number of measured prompts was normalized to activity and correlated to body weight, mass per body length and body mass index to find the best predicting parameter. Next, a patient-specific activity was derived and subsequently validated in 119 additional patients. Image quality was scored by three experts on a four-point scale.

RESULTS

Both image quality and prompts decreased in heavier patients when using a fixed activity (p < .005). Body weight was used to derive a new activity formula: Activity = 8.3 MBq/kg. When applying this formula, both measured prompts and scored image quality became independent of body weight (p > .60).

CONCLUSION

Administrating a Rb-82 activity that linearly depends on body weight resulted in a constant image quality across all patients and is recommended.

Value of regional myocardial flow measurements using Rubidium-82 PET

Funding Acknowledgements

Type of funding sources: None.

Introduction

The combination of myocardial blood flow (MBF) measurements using Rubidium-82 (Rb-82) PET and visual assessment of the PET images is increasingly used due to its high diagnostic and prognostic value. Typically, flow measurements are calculated and used for the myocardium as a whole (global). However, small regional flow deficits may go unnoticed when only looking at global flow values.

Purpose

To compare the diagnostic value of regional and global myocardial flow measurements using Rb-82 PET in the detection of obstructive CAD.

Methods

 We retrospectively included 1034 patients with no history of coronary artery disease (CAD) referred for rest and regadenoson-induced stress Rb-82 PET/CT. MBFs were calculated using Lortie’s one-tissue compartment model. Myocardial flow reserve (MFR) was calculated as the ratio of MBF during stress and rest. Regional flow was determined per vessel and per segment. Vessel MFR was defined as the lowest flow reserve of LAD, LCX and RCA territories and segmental MFR as the lowest flow reserve in all 17 segments. Follow-up data were obtained from medical records. Patients were classified to have obstructive CAD if follow-up included a positive invasive coronary angiography (ICA), percutaneous coronary intervention (PCI), coronary artery bypass grafting (CABG) or all cause death. Receiver-operating characteristic (ROC) analyses were constructed to compare the diagnostic value of global and regional flow values.

Results

Follow-up was obtained in all 1034 patients and the median follow-up time was 2.1 years. Myocardial flow reserve values were significantly lower (p &lt; 0.001) in the 128 patients classified with obstructive CAD than in the 906 patients without obstructive CAD: global MFR  (median 1.9 [interquartile range 1.6-2.4] vs. 2.5 [2.1-2.9]); vessel MFR (1.6 [1.3-2.1]  vs. 2.3 [1.9-2.6]); Segmental MFR (1.3 [0.9-1.7] vs. 1.9 [1.6-2.2]). The area under the curve of vessel MFR (0.79 ± 0.02) and segmental MFR (0.81 ± 0.02) were similar but significantly (p &lt; 0.001) larger than the area of global MFR (0.75 ± 0.03), as shown in the Figure.

Conclusion

The diagnostic value improved with the use of regional MFR instead of global MFR measurements in the detection of obstructive CAD. Therefore, it seems that visual assessment of PET images can best be combined with regional flow measurements either on a per vessel or a per segment basis in Rubidium-82 PET myocardial perfusion imaging.

Prediction of obstructive coronary artery disease after Rb-82 PET myocardial perfusion imaging and coronary artery calcium scoring using machine learning

Funding Acknowledgements

Type of funding sources: None.

Background

Accurate risk stratification in patients with suspected stable coronary artery disease (CAD) is essential for choosing an appropriate treatment strategy but remains challenging in clinical practice.

Purpose

Our aim was to develop and validate a risk model to predict the presence of obstructive CAD after Rubidium-82 PET and a coronary artery calcium score (CACS) scan using a machine learning (ML) algorithm.

Methods

We retrospectively included 1007 patients without prior cardiovascular history and  a low-intermediate pre-test likelihood, referred for rest and regadenoson-induced stress Rubidium-82 PET combined with a CACS scan. Multiple features were included in the ML model; PET derived features such as summed difference score and flow values, CACS, cardiovascular risk factors (cigarette smoking, hypertension, hypercholesterolemia, diabetes, positive family history of CAD), medication; age; gender; body mass index; creatinine serum values; and visual PET interpretation. An XGBoost ML algorithm was developed using a subset of 805 patients to predict obstructive CAD by using 5-fold cross validation in combination with a grid search. Obstructive CAD during follow-up was defined as a significant stenosis during invasive coronary angiography, a percutaneous coronary intervention or a coronary artery bypass graft procedure. The ML algorithm was validated with unseen data of the remaining 202 patients.

Results

Application of the XGBoost algorithm resulted in an area under the curve (AUC) of 0.93 using the training data (n = 805) and an AUC of 0.89 using the unseen data (n = 202) in predicting obstructive CAD. The strongest predictors were the CAC-scores and quantitative PET derived features. The classical risk factors and medication hardly provided an added value in the prediction of obstructive CAD.

Conclusion

The developed ML algorithm is able to provide individualized risk stratification by predicting the probability of obstructive CAD.  Although validation with a larger dataset could result in a more well defined performance range, this model already shows potential to be implemented in the diagnostic workflow.

The Influence of the Exclusion of Central Necrosis on [18F]FDG PET Radiomic Analysis

BACKGROUND

Central necrosis can be detected on [¹⁸F]FDG PET/CT as a region with little to no tracer uptake. Currently, there is no consensus regarding the inclusion of regions of central necrosis during volume of interest (VOI) delineation for radiomic analysis. The aim of this study was to assess how central necrosis affects radiomic analysis in PET.

METHODS

Forty-three patients, either with non-small cell lung carcinomas (NSCLC, n = 12) or with pheochromocytomas or paragangliomas (PPGL, n = 31), were included retrospectively. VOIs were delineated with and without central necrosis. From all VOIs, 105 radiomic features were extracted. Differences in radiomic features between delineation methods were assessed using a paired t-test with Benjamini-Hochberg multiple testing correction. In the PPGL cohort, performances of the radiomic models to predict the noradrenergic biochemical profile were assessed by comparing the areas under the receiver operating characteristic curve (AUC) for both delineation methods.

RESULTS

At least 65% of the features showed significant differences between VOI_vital-tumour and VOI_gross-tumour (65%, 79% and 82% for the NSCLC, PPGL and combined cohort, respectively). The AUCs of the radiomic models were not significantly different between delineation methods.

CONCLUSION

In both tumour types, almost two-third of the features were affected, demonstrating that the impact of whether or not to include central necrosis in the VOI on the radiomic feature values is significant. Nevertheless, predictive performances of both delineation methods were comparable. We recommend that radiomic studies should report whether or not central necrosis was included during delineation.

Appearance of the levator ani muscle subdivisions on 3D transperineal ultrasound

BACKGROUND

The levator ani muscle (LAM) consists of different subdivisions, which play a specific role in the pelvic floor mechanics. The aim of this study is to identify and describe the appearance of these subdivisions on 3-Dimensional (3D) transperineal ultrasound (TPUS). To do so, a study designed in three phases was performed in which twenty 3D TPUS scans of vaginally nulliparous women were assessed. The first phase was aimed at getting acquainted with the anatomy of the LAM subdivisions and its appearance on TPUS: relevant literature was consulted, and the TPUS scan of one patient was analyzed to identify the puborectal, iliococcygeal, puboperineal, pubovaginal, and puboanal muscle. In the second phase, the five LAM subdivisions and the pubic bone and external sphincter, used as reference structures, were manually segmented in volume data obtained from five nulliparous women at rest. In the third phase, intra- and inter-observer reproducibility were assessed on twenty TPUS scans by measuring the Dice Similarity Index (DSI).

RESULTS

The mean inter-observer and median intra-observer DSI values (with interquartile range) were: puborectal 0.83 (0.13)/0.83 (0.10), puboanal 0.70 (0.16)/0.79 (0.09), iliococcygeal 0.73 (0.14)/0.79 (0.10), puboperineal 0.63 (0.25)/0.75 (0.22), pubovaginal muscle 0.62 (0.22)/0.71 (0.16), and the external sphincter 0.81 (0.12)/0.89 (0.03).

CONCLUSION

Our results show that the LAM subdivisions of nulliparous women can be reproducibly identified on 3D TPUS data.

Experimental validation of absolute SPECT/CT quantification for response monitoring in patients with coronary artery disease

BACKGROUND

Quantitative SPECT enables absolute quantification of uptake in perfusion defects. The aim of this experimental study is to assess quantitative accuracy and precision of a novel iterative reconstruction technique (Evolution; GE Healthcare) for the potential application of response monitoring using ^99mTc-tetrofosmin SPECT/CT in patients with coronary artery disease (CAD).

METHODS

Acquisitions of an anthropomorphic torso phantom with cardiac insert containing defects (with varying sizes), filled with ^99mTc-pertechnetate, were performed on a SPECT/CT (Discovery 670 Pro, GE Healthcare). Subsequently, volumes of interest of the defects were manually drawn on CT to assess the recovery coefficient (RC). Bull's eye plots were composed to evaluate the uptake per segment. Finally, ^99mTc-tetrofosmin SPECT/CT scans of 10 CAD patients were used to illustrate clinical application.

RESULTS

The phantom study indicated that Evolution showed convergence after 7 iterations and 10 subsets. The average repeatability deviation of all configurations was 2.91% and 3.15% (%SD mean) for filtered (Butterworth) and unfiltered data, respectively. The accuracy after post-filtering was lower compared to the unfiltered data with a mean (SD) RC of 0.63 (0.05) and 0.70 (0.07), respectively (p < 0.05). More artificial defects were found on Bull's eye plots created with the unfiltered data compared to filtered data. Eight out of ten patients showed significant changes in uptake before and after treatment (p < 0.05).

CONCLUSION

Quantification of ^99mTc-tetrofosmin SPECT/CT seems feasible for CAD patients when 7 iterations (10 subsets), Butterworth post-filtering (cut off frequency 0.52 in cycles/cm, order of 5) and manual CT-delineation are applied. However, future prospective patient studies are required for clinical application.

The Automatic Quantification of Morphological Features of Pectus Excavatum Based on Three-Dimensional Images

Visual examination and quantification of severity are essential for clinical decision making in patients with pectus excavatum. Yet, visual assessment is prone to inter- and intra-observer variability and current quantitative methods are inadequate. This study aims to develop and evaluate a novel, automatic and non-invasive method to objectively quantify pectus excavatum morphology based on three-dimensional images. Key steps of the automatic analysis are normalization of image orientation, slicing, and computation of the morphological features encompassing pectus depth, width, length, volume, position, steepness, flaring, asymmetry and mean cross-sectional area. A digital phantom mimicking a patient with pectus excavatum was used to verify the analysis method. Prospective three-dimensional imaging and subsequent surface analysis in patients with pectus excavatum was performed to assess clinical feasibility. Verification of the developed analysis tool demonstrated 100% reproducibility of all morphological feature values. Calculated parameters compared to the predetermined phantom dimensions were accurate for all but four features. The pectus width, length, volume and steepness showed an error of 4 mm (4%), 2 mm (2%), 12 mL (5%) and 1 degree (3%), respectively. Prospective imaging of 52 patients (88% males) demonstrated the feasibility of the developed tool to quantify morphological features of pectus excavatum in the clinical setting. Mean duration to calculate all features in one patient was 7.6 seconds. We have developed and presented a non-invasive pectus excavatum surface analysis tool, that is feasible to automatically quantify morphological features based on three-dimensional images with promising accuracy and reproducibility.

Cerenkov Luminescence Imaging in Prostate Cancer: Not the Only Light That Shines

Cerenkov luminescence imaging (CLI) is a novel imaging technology that might have the ability to assess surgical margins intraoperatively during prostatectomy using ⁶⁸Ga-prostate-specific membrane antigen (⁶⁸Ga-PSMA-11). This study evaluated the accuracy of CLI compared with histopathology and, as an exploratory objective, investigated the characteristics of the identified chemiluminescence signal. Methods: After intravenous injection of a mean ⁶⁸Ga-PSMA-11 activity of 69 MBq intraoperatively, all excised specimens were imaged with CLI. Areas of increased signal were marked for histopathologic comparison and scored for the likelihood of being a positive surgical margin (PSM) using a 5-point Likert scale. In addition, the chemiluminescence signal was investigated in 3 radioactive and 3 nonradioactive specimens using CLI. Results: In 15 patients, the agreement between CLI and histopathology was 60%; this improved to 83% when including close surgical margins (≤1 mm). In 6 hot spots, CLI correctly identified PSMs on histopathology, located at the apex and mid prostate. In all 15 patients, an increased signal at the prostate base was observed, without the presence of the primary tumor in this area in 8 patients. This chemiluminescence signal was also observed in nonradioactive prostate specimens, with a half-life of 48 ± 11 min. The chemiluminescence hampered the visual interpretation of 4 PSMs at the base. Conclusion: CLI was able to correctly identify margin status, including close margins, in 83% of the cases. The presence of a diathermy-induced chemiluminescent signal hampered image interpretation, especially at the base of the prostate. In the current form, CLI is most applicable to detect PSMs and close margins in the apex and mid prostate.

A review on imaging techniques and quantitative measurements for dynamic imaging of cerebral aneurysm pulsations

Measurement of intracranial aneurysm wall motion may refine the current rupture risk estimation. A golden standard for measuring aneurysm pulsation is lacking. The aim is to evaluate magnitudes of aneurysm pulsation as published in current literature. Embase and PubMed were searched for publications containing quantitative measures of cardiac-cycle related cerebral aneurysm pulsation (no date or language restrictions). Eleven studies were included, covering 197 unruptured and untreated cerebral aneurysms. Quantitative pulsation measurements were extracted from the studies. Characteristics of the study population and aneurysms were taken into account, as well as the imaging modality, scanning technique and data processing methods used. A meta-analysis was performed of studies with similar methodologies and individual IA measures and locations. The magnitude of the absolute volume pulsations varied between 14 ± 9 mm³ and 106 ± 123 mm³ and the mean relative volume change varied between 5 and 36%. The meta-analysis revealed a positive correlation between size and absolute volume change. The relative volume change in Basilar artery aneurysms seems smaller. No authors were contacted for original study data and articles only describing visual pulsations were excluded. The variation in methodologies impedes an accurate estimation of the magnitude of IA pulsations. Validation of aneurysm pulsation measurement is crucial prior to clinical studies evaluating IA pulsatility in relation to IA rupture risk. Prerequisite is a reliable and robust imaging method with high spatial and temporal resolution and standardization of the image analysis methods.

Clinical value of machine learning-based interpretation of I-123 FP-CIT scans to detect Parkinson's disease: a two-center study

PURPOSE

Our aim was to develop and validate a machine learning (ML)-based approach for interpretation of I-123 FP-CIT SPECT scans to discriminate Parkinson's disease (PD) from non-PD and to determine its generalizability and clinical value in two centers.

METHODS

We retrospectively included 210 consecutive patients who underwent I-123 FP-CIT SPECT imaging and had a clinically confirmed diagnosis. Linear support vector machine (SVM) was used to build a classification model to discriminate PD from non-PD based on I-123-FP-CIT striatal uptake ratios, age and gender of 90 patients. The model was validated on unseen data from the same center where the model was developed (n = 40) and consecutively on data from a different center (n = 80). Prediction performance was assessed and compared to the scan interpretation by expert physicians.

RESULTS

Testing the derived SVM model on the unseen dataset (n = 40) from the same center resulted in an accuracy of 95.0%, sensitivity of 96.0% and specificity of 93.3%. This was identical to the classification accuracy of nuclear medicine physicians. The model was generalizable towards the other center as prediction performance did not differ thereby obtaining an accuracy of 82.5%, sensitivity of 88.5% and specificity of 71.4% (p = NS). This was comparable to that of nuclear medicine physicians (p = NS).

CONCLUSION

ML-based interpretation of I-123-FP-CIT scans results in accurate discrimination of PD from non-PD similar to visual assessment in both centers. The derived SVM model is therefore generalizable towards centers using comparable acquisition and image processing methods and implementation as diagnostic aid in clinical practice is encouraged.

In vitro performance of echoPIV for assessment of laminar flow profiles in a carotid artery stent

Purpose: Detailed blood flow studies may contribute to improvements in carotid artery stenting. High-frame-rate contrast-enhanced ultrasound followed by particle image velocimetry (PIV), also called echoPIV, is a technique to study blood flow patterns in detail. The performance of echoPIV in presence of a stent has not yet been studied extensively. We compared the performance of echoPIV in stented and nonstented regions in an in vitro flow setup.

Approach: A carotid artery stent was deployed in a vessel-mimicking phantom. High-frame-rate contrast-enhanced ultrasound images were acquired with various settings. Signal intensities of the contrast agent, velocity values, and flow profiles were calculated.

Results: The results showed decreased signal intensities and correlation coefficients inside the stent, however, PIV analysis in the stent still resulted in plausible flow vectors.

Conclusions: Velocity values and laminar flow profiles can be measured in vitro in stented arteries using echoPIV.

Day-to-day variability of [68Ga]Ga-PSMA-11 accumulation in primary prostate cancer: effects on tracer uptake and visual interpretation

PURPOSE

Prostate-specific membrane antigen (PSMA) agents, such as [⁶⁸Ga]Ga-PSMA-11, have an unprecedented accuracy in staging prostate cancer (PCa) and detecting disease recurrence. PSMA PET/CT may also be used for response monitoring by displaying molecular changes, instead of morphological changes alone. However, there are still limited data available on the variability in biodistribution and intra-prostatic uptake of PSMA targeting radiotracers. Therefore, the aim of this study was to assess the repeatability of [⁶⁸Ga]Ga-PSMA-11 uptake in primary PCa patients in a 4-week interval.

METHODS

Twenty-four primary PCa patients were prospectively included, who already were scheduled for [⁶⁸Ga]Ga-PSMA-11 PET/CT scan on clinical indication (≥ cT3, Gleason score ≥ 7 or PSA ≥ 20 ng/mL). These patients received two [⁶⁸Ga]Ga-PSMA-11 PET/CT scans with a 4-week interval. No treatment was started in between the scans. Semiquantitative measurements (SUL_max, SUL_mean, and SUL_peak) were determined in the prostate tumor, normal tissues, and blood pool. The repeatability coefficient of every region was determined. All scans were visually analyzed by two nuclear medicine physicians.

RESULTS

Within-subject coefficient of variation of [⁶⁸Ga]Ga-PSMA-11 uptake between the two scans was on average 10% in the prostate tumor, normal tissues (liver, kidney, parotid), and blood pool. The repeatability coefficient of the prostate tumor was 18% for SUL_peak and 22% for SUL_max. Lesion uptake was visually different in 5 patients, though not clinically relevant.

CONCLUSION

Results of test-retest [⁶⁸Ga]Ga-PSMA-11 PET/CT scans in a 4-week interval show that [⁶⁸Ga]Ga-PSMA-11 uptake is repeatable, with a clinical irrelevant variation in tumor and physiological distribution. Based on the presented repeatable uptake, [⁶⁸Ga]Ga-PSMA-11 PET/CT scans can potentially be used for disease surveillance and therapy response monitoring. Changes in uptake larger than the RC are therefore likely to reflect actual biological changes in PSMA expression. Trial registration NL8263 at Trialregister.nl retrospectively registered on 03-01-2020. https://www.trialregister.nl/trial/8263.

Adding the temporal domain to PET radiomic features

Background

Radiomic features, extracted from positron emission tomography, aim to characterize tumour biology based on tracer intensity, tumour geometry and/or tracer uptake heterogeneity. Currently, radiomic features are derived from static images. However, temporal changes in tracer uptake might reveal new aspects of tumour biology. This study aims to explore additional information of these novel dynamic radiomic features compared to those derived from static or metabolic rate images.

Methods

Thirty-five patients with non-small cell lung carcinoma underwent dynamic [¹⁸F]FDG PET/CT scans. Spatial intensity, shape and texture radiomic features were derived from volumes of interest delineated on static PET and parametric metabolic rate PET. Dynamic grey level cooccurrence matrix (GLCM) and grey level run length matrix (GLRLM) features, assessing the temporal domain unidirectionally, were calculated on eight and sixteen time frames of equal length. Spearman’s rank correlations of parametric and dynamic features with static features were calculated to identify features with potential additional information. Survival analysis was performed for the non-redundant temporal features and a selection of static features using Kaplan-Meier analysis.

Results

Three out of 90 parametric features showed moderate correlations with corresponding static features (ρ≥0.61), all other features showed high correlations (ρ>0.7). Dynamic features are robust independent of frame duration. Five out of 22 dynamic GLCM features showed a negligible to moderate correlation with any static feature, suggesting additional information. All sixteen dynamic GLRLM features showed high correlations with static features, implying redundancy. Log-rank analyses of Kaplan-Meier survival curves for all features dichotomised at the median were insignificant.

Conclusion

This study suggests that, compared to static features, some dynamic GLCM radiomic features show different information, whereas parametric features provide minimal additional information. Future studies should be conducted in larger populations to assess whether there is a clinical benefit of radiomics using the temporal domain over traditional radiomics.

Intracranial aneurysm growth: consistency of morphological changes

OBJECTIVE

Previous studies have shown a relation between growth and rupture of intracranial aneurysms. Additionally, several morphological characteristics are frequently measured to estimate rupture risk. Little is known about how the rupture risk is associated with morphological characteristic changes during growth. The aim of this study was to provide insights into how morphological characteristics, associated with rupture, change during an aneurysm's growth.

METHODS

The authors retrospectively identified patients with longitudinal MRA images of unruptured growing aneurysms. The MRA images had an in-plane resolution of 0.2-0.5 mm and a slice thickness of 0.2-0.75 mm. Therefore, growth was defined as an increase of at least 0.5 mm in two directions or 1 mm in one direction. Using the MRA images, the authors semiautomatically segmented the aneurysm and the perianeurysmal vasculature. Twelve morphological characteristics were automatically measured. These characteristics were related to size (diameter, height, width, neck diameter, volume, surface area, aspect ratio, height-width ratio, and bottleneck factor) and shape (ellipticity index, nonsphericity index, and undulation index) of the aneurysm. Morphological characteristics before and after growth were compared using the Wilcoxon signed-rank test.

RESULTS

The authors included 31 patients with 38 growing aneurysms. The aneurysms' growth was detected after a mean of 218 weeks (range 23-567 weeks). A significant increase was seen in all size-related characteristics, and the bottleneck factor also significantly increased (from a median of 1.00 [IQR 0.85-1.04] to 1.03 [IQR 0.93-1.18]), while the ellipticity index decreased (from a median of 0.26 [IQR 0.25-0.28] to 0.25 [IQR 0.24-0.26]). The changes in size ratios and shape indices varied largely among patients. Larger aneurysms more often showed an increase in shape ratios.

CONCLUSIONS

Although aneurysm growth, size-related characteristics, bottleneck factor, and ellipticity index changed significantly during growth, most size ratios and shape indices showed inconsistent changes among aneurysms. This suggests that, for an accurate rupture prediction, morphological parameters need to be reassessed after growth.

Effect of irreversible electroporation parameters and the presence of a metal stent on the electric field line pattern

The final ablation zone created with irreversible electroporation (IRE) depends on the size, shape and strength of the electric field that is influenced by several parameters. A profound understanding of the effect of IRE parameter alterations on the electric field are a prerequisite for a safe and effective treatment. Here, we demonstrate a semolina in castor oil model that enables visualization of the static electric field developed by a high-voltage generator between two needle-electrodes. We intuitively visualize the variation in electric field line pattern for selected IRE parameters; active needle length, inter-needle distance, applied voltage and presence of a nearby metal stent, by cameras in three dimensions. The observations were compared to and supported by two-dimensional numerical simulations of the electric field. Our semolina model visualizes the disturbance of the electric field by a metal stent, potentially leading to an incomplete tumour ablation between the needles. The reduction in electric field strength and the area at risk for incomplete tumour ablation are confirmed by the numerical simulations. The semolina model provides insight in the fundamental physics of the electric field, the effect of alterations in IRE parameter combinations and presence of a metal stent within the ablation zone.

Improving Lives in Three Dimensions: The Feasibility of 3D Printing for Creating Personalized Medical Aids in a Rural Area of Sierra Leone

The aim of this feasibility study was to investigate how a 3D printer could be put to its best use in a resource-limited healthcare setting. We have examined whether a 3D printer can contribute to making prostheses, braces, or splints for patients who underwent major limb amputation because of complex wounds, for example, due to burns and subsequent scarring, accidents, conflicts, or congenital abnormalities. During a 3-month period, we investigated the benefits of customized, 3D-printed arm prostheses, splints, and braces in Sierra Leone. Using a handheld 3D scanner and a 3D printer, patient-specific medical aids were designed, manufactured, and tested. Questionnaires regarding patient satisfaction and the functionality of the prostheses were used for a short-term follow-up. Four esthetic prostheses were designed: two prostheses of the hand, one of the forearm, and one of the entire arm. Follow-ups were conducted after 3 to 4 weeks to investigate the quality of the prostheses and to complete a patient questionnaire. Even though the prostheses primarily fulfill esthetic needs, they also exhibit some degree of functionality. In addition, four splints for hands and arms were made to prevent scar contractures after skin transplantation. Finally, a brace for a young boy with kyphoscoliosis was manufactured. The boy has accepted the brace and will be followed up in the months to come. Long-term follow-up is required to prove the sustainability of the 3D-printed brace and prosthetic arms. Further research into how to sustain and refine this project is underway.

Geometrical changes in Anaconda endograft limbs after endovascular aneurysm repair: A potential predictor for limb occlusion

The emergence of limb occlusion after endovascular aneurysm repair may be related to the conformational changes between the endograft structure and the patient's anatomy. This study analyzed detailed geometric changes of Anaconda endograft (Terumo Aortic, Inchinnan, Scotland, UK) limbs during the cardiac cycle-based computed tomography on serial imaging after graft implantation. Fifteen patients (mean age 72.8 ± 3.7 years; 14 men) underwent postoperative electrocardiogram-gated computed tomography scans according to a prospective study design between April 2014 and May 2017. Changes in curvature, length of the limbs, and distances between successive stent rings (inter-ring distance) of the endograft limbs during a 2-year follow-up period were quantified using meticulous image processing methods involving image registration, centerline extraction, and model-based stent-ring segmentation. From discharge to 24 months, mean curvature increased significantly by 9.6 m^-1 (standard deviation [SD], 11.1 m^-1; 95% confidence interval [CI], 3.4 to 15.8 m^-1; P = .002) for the right limbs and by 6.1 m^-1 (SD 9.4 m^-1; 95% CI, 0.8 to 11.5 m^-1; P = .21) for the left limbs. The length of the right limbs decreased significantly, by 9.5 mm (SD 7.6 mm; 95% CI, 3.5 to 15.6 mm; P = .002); the length of the left limbs decreased by 10.1 mm (SD 5.1 mm; 95% CI, 5.9 to 14.2 mm; P < .001). The minimal inter-ring distance decreased by 0.36 mm (SD 0.26 mm; 95% CI, 0.17 to 0.55 mm; P < .001) for the right limbs and 0.35 mm (SD 0.19 mm; 95% CI, 0.21 to 0.49 mm; P < .001) for the left limbs. Cardiac pulsatility-induced changes in curvature, limb length, and inter-ring distance were negligible (2%, 0.3% and 0.3%, respectively). Changes in the geometry of the Anaconda endograft limbs after endovascular aortic aneurysm repair were observed during a 2-year follow-up manifest as an increase in curvature, shortening of the stent-graft limbs, and a corresponding decrease in inter-ring distance. These stent-graft conformational changes could result in inward folding of the graft fabric, which may relate to the emergence of limb occlusion. Further investigation of these metrics in a larger cohort involving patients with and without occlusions may allow determination of their predictive value.

(In)comparability of Carotid Artery Stent Characteristics: A Systematic Review on Assessment and Comparison with Manufacturer Data

Purpose

Carotid stent (CS) characteristics, such as radial force, scaffolding and flexibility, are continuously modified by stent manufacturers aiming to improve stent performance. Since manufacturers’ definitions and assessment methods are not disclosed, it is unknown how characteristics of different CSs relate to each other or to published literature. We examined in vitro methodological techniques used to measure CS characteristics and assessed comparability between published papers and outcomes as provided by the manufacturers.

Methods

A systematic review was conducted in MEDLINE, Embase, Cochrane, and Scopus databases. Studies reporting on in vitro investigations of predefined characteristics of CS used in current everyday clinical practice were included. The predefined characteristics were radial force, scaffolding, flexibility, foreshortening, side-branch preservation and visibility. Eight manufacturers of 10 currently used CS were contacted and data on the predefined device characteristics was requested.

Results

12 published articles were included and six stent manufacturers provided data on six stents (two refused to share data). Used methodologies to measure stent characteristics in published literature and manufacturer data varied greatly for all included characteristics except foreshortening. The number of different units of measurement to express outcomes ranged from two for foreshortening to six for radial force.

Conclusion

A variety of methodologies and outcome measures is used to quantify CS characteristics, which hampers comparisons between published studies and manufacturer data. Future studies are encouraged to synchronize methodologies and outcome measures. Manufacturers are encouraged up to increase transparency of applied testing methodologies and outcomes.

Electronic supplementary material

The online version of this article (10.1007/s00270-020-02499-1) contains supplementary material, which is available to authorized users.

A Systematic Review for the Design of In Vitro Flow Studies of the Carotid Artery Bifurcation

PURPOSE

In vitro blood flow studies in carotid artery bifurcation models may contribute to understanding the influence of hemodynamics on carotid artery disease. However, the design of in vitro blood flow studies involves many steps and selection of imaging techniques, model materials, model design, and flow visualization parameters. Therefore, an overview of the possibilities and guidance for the design process is beneficial for researchers with less experience in flow studies.

METHODS

A systematic search to in vitro flow studies in carotid artery bifurcation models aiming at quantification and detailed flow visualization of blood flow dynamics results in inclusion of 42 articles.

RESULTS

Four categories of imaging techniques are distinguished: MRI, optical particle image velocimetry (PIV), ultrasound and miscellaneous techniques. Parameters for flow visualization are categorized into velocity, flow, shear-related, turbulent/disordered flow and other parameters. Model materials and design characteristics vary between study type.

CONCLUSIONS

A simplified three-step design process is proposed for better fitting and adequate match with the pertinent research question at hand and as guidance for less experienced flow study researchers. The three consecutive selection steps are: flow parameters, image modality, and model materials and designs. Model materials depend on the chosen imaging technique, whereas choice of flow parameters is independent from imaging technique and is therefore only determined by the goal of the study.

Quantitative imaging: systematic review of perfusion/flow phantoms

Background

We aimed at reviewing design and realisation of perfusion/flow phantoms for validating quantitative perfusion imaging (PI) applications to encourage best practices.

Methods

A systematic search was performed on the Scopus database for “perfusion”, “flow”, and “phantom”, limited to articles written in English published between January 1999 and December 2018. Information on phantom design, used PI and phantom applications was extracted.

Results

Of 463 retrieved articles, 397 were rejected after abstract screening and 32 after full-text reading. The 37 accepted articles resulted to address PI simulation in brain (n = 11), myocardial (n = 8), liver (n = 2), tumour (n = 1), finger (n = 1), and non-specific tissue (n = 14), with diverse modalities: ultrasound (n = 11), computed tomography (n = 11), magnetic resonance imaging (n = 17), and positron emission tomography (n = 2). Three phantom designs were described: basic (n = 6), aligned capillary (n = 22), and tissue-filled (n = 12). Microvasculature and tissue perfusion were combined in one compartment (n = 23) or in two separated compartments (n = 17). With the only exception of one study, inter-compartmental fluid exchange could not be controlled. Nine studies compared phantom results with human or animal perfusion data. Only one commercially available perfusion phantom was identified.

Conclusion

We provided insights into contemporary phantom approaches to PI, which can be used for ground truth evaluation of quantitative PI applications. Investigators are recommended to verify and validate whether assumptions underlying PI phantom modelling are justified for their intended phantom application.