Interventional, Multimodal & Auxiliary Engineering
Engineering/Interventional/Safety Thursday, 20 May 2021

Oral Session - Interventional, Multimodal & Auxiliary Engineering
Engineering/Interventional/Safety
Thursday, 20 May 2021 18:00 - 20:00
  • An MR Safe Steerable Catheter for MR-guided Endovascular Interventions
    Mohamed E. M. K. Abdelaziz1, Libaihe Tian1, Thomas Lottner2, Simon Reiss2, Klaus Düring3, Guang-Zhong Yang4, Michael Bock2, and Burak Temelkuran1
    1Hamlyn Centre for Robotic Surgery, Imperial College London, London, United Kingdom, 2Dept. of Radiology, Medical Physics, Medical Center, University of Freiburg, Freiburg, Germany, 3MaRVis Interventional GmBH, Krün, Germany, 4Shanghai Jiao Tong University, Shanghai, China
    Real-time MR image guidance of a bespoke 7F steerable catheter was demonstrated in a vessel phantom in a clinical 3T MR scanner to evaluate its mechanical efficacy and MR visibility. The catheter could be probed efficiently and was visible over its entirety thanks to its MR visible pull-wires.
    (a) Laser cut catheter tip with the positions of the additional tip markers labeled. (b) Rendering of the handle showing the dual-dials, mechanical locking mechanism for shape locking the catheter's distal end and a proximal luer lock for contrast agent injection and/or guidewire insertion through the central lumen of the catheter. (c) Catheter tip manipulation in 4-directions using the dials demonstrated in the handle design. (d) Handle functional prototype.
    Sequence of steps showing the right renal artery probing (a-c) and subsequent contrast injection (d)
  • Real-Time Slice Steering for MR-Guided Interventions Using Endovascular Devices Equipped with Passive MRI Markers
    Daniel Christopher Hoinkiss1, Han Nijsink2, Paul Borm3, Sabrina Haase1, Jan Strehlow1, Jurgen Futterer2, and Torben Pätz1
    1Fraunhofer Institute for Digital Medicine MEVIS, Bremen, Germany, 2Radboud University Medical Centre (Radboudumc), Nijmegen, Netherlands, 3Nano4imaging GmbH, Aachen, Germany
    The presented workflow for automatic interventional slice steering combines real-time tracking information based on passive MRI markers with pre-calculated therapy planning for smooth and accurate device monitoring that increases precision during MR-guided interventions.
    Figure 5: a) Automatic slice steering based on reconstructed wire geometry (top) and vascular information (bottom) with a non-moving wire. A slight jittering of the real-time slices can be seen in the top row where orientation of the slices changes between timepoints (see differences in blue circle). b) and c) Time series of real-time images, acquired with slice steering and while performing guidewire motion. The experiments were performed at different MRI sites using two different MRI phantoms.
    Figure 3: An additional option allows to use a-priori information of pre-processed vessel masks and the envisioned path for intervention (a). A smoothed path representation is synchronized to the reconstructed wire geometry (b) and sample points are calculated based on the interventional path (c) with which the imaging plane information can be calculated (d). This allows to better align the real-time MRI slices to the vascular characteristics.
  • Minimal Artifact Actively Shimmed Metallic Needles for Interventional MRI
    Saikat Sengupta1, Xinqiang Yan1, Tamarya Hoyt2, Anthony Gunderman3, and Yue Chen3
    1Department of Radiology, Vanderbilt University Institute of Imaging Science, Nashville, TN, United States, 2Department of Radiology, Vanderbilt University Medical Center, Nashville, TN, United States, 3Department of Mechanical Engineering, University of Arkansas, Fayetteville, AR, United States
    Signal voids caused by metallic probes are a challenge in many interventional MRI applications. Here, we demonstrate effective recovery of metallic probe induced signal loss using a novel, actively shimmed metallic probe at 3 Tesla.

    Figure 4

    Results of active shimming, 1 x 1 x 1 mm3 3D GRE images and fieldmaps. Excellent recovery of lost signal is achieved in all orientations using pre-estimated shim currents. The width of the signal void approaches the needle width in all cases. Fieldmaps show correction of the underlying ΔB0. Note that the regions closest to the rod with field information in the ‘With Shim’ case have no corresponding field information in the ‘No Shim’ case due to signal loss. The percentage of signal recovered (not accounting for the needle itself) is indicated for all orientations.

    Figure 2

    Active shimming hardware (a) Shim insert design with bevel and orthogonal slots for shim wires (b) Shim insert inside the titanium needle (c) Needle placed in holder with compartment for electronics (d) Phantoms with guide holes at different angles used for experiments.

  • In Vivo Susceptibility-based Positive Contrast Imaging of MR Compatible Metallic Devices Based on Modified Slab-Selective 3D SPACE Sequence
    Caiyun Shi1,2, Dong Liang1,2,3, Zhilang Qiu1, Xin Liu1,2, Yanjie Zhu1,2, and Haifeng Wang1,2
    1Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China, 2Shenzhen College of Advanced Technology, University of Chinese Academy of Sciences, Shenzhen, China, 3Research Centre for Medical AI, Shenzhen Institutes of Advanced Technology, Chinese Academy of Science, Shenzhen, China
    Susceptibility-based positive contrast MR imaging exhibits excellent efficacy for visualizing the MR compatible metallic devices. In this work, a novel method is developed to realize the 3D susceptibility-based positive contrast MR imaging on in vivo human patient.
    Figure 1.(a) Diagram of the modified 3D SPACE sequence. Compared to the conventional SPACE sequence, each readout in the modified SPACE sequence is shifted by a short time Tshift. The nature of the SPACE sequence helps to reduce the rapid signal dephasing and phase wrapping due to the high susceptibility of the metallic devices. By processing two sets of imaging data (with and without Tshift), the local field map (phase change) induced by the susceptibility of interventional devices can be calculated. (b) Variable-flip-angle of the 3D SPACE with refocusing pulse.
    Figure 5.Representative results of the in vivo experiments of a human patient with one tumour of the scapula. (a) Magnitude images generated by the proposed method; (b~c) The positive contrast images generated by the proposed 3D SPACE-based method, 2D FSE based sequence; (d) Magnitude images generated by the SUMO, susceptibility gradient mapping using the original resolution; (e~f) The positive contrast images generated by the 3D SUMO, 2D GRASP (gradient echo acquisition for super-paramagnetic particles).
  • Improvements in in vivo imaging and temperature mapping using passive “propeller-beanie” antenna in transcranial MR-guided focused ultrasound
    Xinqiang Yan1,2, Steven P. Allen3, Craig H. Meyer3, and William A. Grissom1,2,4
    1Vanderbilt University Institute of Imaging Science, Nashville, TN, United States, 2Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States, 3Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, United States, 4Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States
    The “propeller-beanie” antenna method can significantly alleviate the dark band artifacts in in vivo imaging in a Insightec tcMRgFUS system, which improves temperature precision and may enable the use of diffusion imaging to monitor treatment. 
    Figure 3 In vivo transmit RF field (B1+) maps in coronal, sagittal and axial planes. Top: without passive wires. Bottom: with 11-cm-long passive wires. The full-FOV maps including the water bath are shown in the first row of each case, and maps masked to the brain are shown in the second row of each case. By using 11-cm wires, the B1+ uniformity was considerably improved. However, it is also noted that the B1+ is still not so uniform as that of in simulation. This could be attributed to the fact that wires deviated from the central position in the axial plane.

    Figure 2 Optimization setup to evaluate the best length for the passive wires (left) and the simulated RF transmit magnitude (|B1+|) maps. The optimal wire length was 11 cm (highest B1+ efficiency near the middle of the brain) which was used for the in vivo experiments.

  • An Anthropomorphic Pelvis Phantom for Prostate Brachytherapy and Biopsy
    Dominik F. Bauer1,2, Eva Oelschlegel1,2, Alena-Kathrin Golla1,2, Anne Adlung1,2, Tom Russ1,2, Ingo Hermann1,2, Irène Brumer1,2, Julian Rosenkranz3, Fabian Tollens4, Sven Clausen5, Philipp Aumüller5, Lothar R. Schad1,2, Dominik Nörenberg4, and Frank G. Zöllner1,2
    1Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany, 2Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany, 3Fraunhofer Institute for Manufacturing Engineering and Automation, Project Group for Automation in Medicine and Biotechnology, Mannheim, Germany, 4Department of Radiology and Nuclear Medicine, University Medical Center Mannheim, Mannheim, Germany, 5Department of Radiation Oncology, University Medical Center Mannheim, Heidelberg University, Mannheim, Germany
    We present an anthropomorphic pelvis phantom with lesions for transperineal and transrectal prostate needle interventions. The phantom is puncturable and suitable for multiparametric MRI (mpMRI) and CT imaging.
    Figure 1: (a) 3D model of the organ structures. (b) Real organ structures. (c) Final phantom: Organ structures embedded in ballistic gelatin.

    Figure 2: (a) Body casting mold. (b) Prostate casting mold including prostate lesions and urethra. The lesions and urethra were placed inside the mold before molding. (c) Bladder and Bladder casting mold.

  • A generic framework for real-time 3D motion estimation from highly undersampled k-space using deep learning
    Maarten Terpstra1,2, Matteo Maspero1,2, Tom Bruijnen1,2, Joost Verhoeff1, Jan Lagendijk1, and Cornelis A.T. van den Berg1,2
    1Department of Radiotherapy, University Medical Center Utrecht, Utrecht, Netherlands, 2Computational Imaging Group for MR diagnostics & therapy, University Medical Center Utrecht, Utrecht, Netherlands
    We propose a deep learning approach for real-time motion estimation from highly accelerated NUFFT-reconstructed acquisitions. At R=30, the model produces accurate motion fields within 200 ms, including k-space acquisition. The proposed model even generalizes to 4D-CT without retraining.
    Figure 3: Typical motion reconstruction at R~18. The computed DVFs were computed using NUFFTreconstructed respiratory-resolved images with 70 respiratory phases (top row, R~18), while the bottom row shows optical flow computed on compressed sense reconstructions. Compared to these DVFs, our model achieves a mean EPE of \(2.48\pm0.42\) mm compared to compressed sense optical flow. However, some motion artifacts remain at the anterior side due to undersampling artifacts.
    Figure 1: Schematic of 4D reconstruction and multi-level motion estimation. Long radial acquisitions (7.3 min.) with motion mixing are respiratory-binned and reconstructed to three spatial resolutions by cropping k-space around k0. These images serve as training data for our multi-resolution motion model -- together with optical flow to inhale, exhale, and midvent positions -- to learn a three-dimensional deformation vector field. The lowest resolution CNN has 3x3x3 3D convolutional filters, the middle CNN has 5x5x5 filters, and the final CNN again uses 3x3x3 filters.
  • Real-time deep artifact suppression using recurrent U-nets for interactive Cardiac Magnetic Resonance imaging.
    Olivier Jaubert1,2, Javier Montalt-Tordera2, Dan Knight2,3, Gerry J. Coghlan2,3, Simon Arridge1, Jennifer Steeden2, and Vivek Muthurangu2
    1Department of Computer Science, UCL, London, United Kingdom, 2Centre for Cardiovascular Imaging, UCL, London, United Kingdom, 3Department of Cardiology, Royal Free London NHS Foundation Trust, London, United Kingdom
    A deep learning based framework using a 2D recurrent residual U-Net trained on multiple orientations is proposed to reconstruct an interactively acquired bSSFP tiny golden angle radial sequence for catheter guidance in patients.
    Figure 5. Animation. Proof of concept interactive acquisition in a healthy subject. Changes of orientations are performed interactively both abruptly and through continuous motion between RVOT, PA and SAX. Short transition periods can be observed before convergence to good image quality.
    Figure 4. Animation. Pulmonary artery (top) and right ventricular outflow tract (bottom) views in two different catheterized patients where the balloon can be seen. The videos are shown for the conventional scan, the gridded images (input to the network), retrospective compressed sensing reconstruction (CS with temporal TV regularization) and proposed real time images.
  • Deep Learning-Driven Automatic Scan Plane Alignment for Needle Tracking in MRI-Guided Interventions
    Xinzhou Li1,2, Yu-Hsiu Lee3, David S. Lu1, Tsu-Chin Tsao3, and Holden H. Wu1,2
    1Radiological Sciences, University of California, Los Angeles, Los Angeles, CA, United States, 2Bioengineering, University of California, Los Angeles, Los Angeles, CA, United States, 3Mechanical and Aerospace Engineering, University of California, Los Angeles, Los Angeles, CA, United States
    This work used a deep learning-based needle localization algorithm in a new automatic workflow to realign the MRI scan plane with the needle. In one degree-of-freedom needle insertion experiments, the MRI scan plane was rapidly and accurately aligned with the needle.
    Figure 4: (a) The needle feature at the entry point and after insertion were displayed for initial plane. The incomplete needle feature after insertion was caused by misalignment between the initial plane with the needle trajectory. Orientation difference (dθ) and Hausdorff distance (HD) are 19.2° and 13.9 mm. The reference needle was extracted by segmenting needle feature on a high-resolution 3D confirmation scan. (b) After executing the proposed automated workflow, the final plane was aligned with the needle and the complete needle feature is visible. dθ and HD were 1.8° and 1.8 mm.
    Figure 2: The workflow for automatic MRI scan plane alignment with the needle using Mask R-CNN and the scan plane control (SPC) module. (I) An initial plane was manually selected based on the needle feature at the entry point on 3-plane localizer images. (II-IV) Scans 1-3 automatically started using the plane selection from the previous step. Mask R-CNN needle localization results were used to automatically select a new plane that aligned with the needle.
  • Ultra-quality 4D-MRI synthesis using deep learning-based deformable image registration
    Haonan Xiao1, Tian Li1, Jiang Zhang1, Ruiyan Ni1, Ge Ren1, Yibao Zhang2, Weiwei Liu2, Weihu Wang2, Hao Wu2, Victor Lee3, Andy Cheung3, Hing-Chiu Chang3, and Jing Cai1
    1The Hong Kong Polytechnic University, Hong Kong, China, 2Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Radiation Oncology, Peking University Cancer Hospital & Institute, Beijing, China, 3The University of Hong Kong, Hong Kong, China
    The synthetic ultra-quality 4D-MRIs showed accurate tumor motion trajectory with significantly improved image quality than the traditional 4D-MRI.
    Figure 2: (a) Visual comparison between (1) original 4D MR, (2) synthetic 4D T1w, (3) synthetic 4D T2w, (4) synthetic 4D DWI (b=50), and (5) synthetic 4D DWI (b=800) images. Tumors were indicated by the yellow arrows. (b) Tumor motion during 2 respiratory cycles in superior-inferior (SI), anterior-posterior (AP), and mid-lateral (ML) direction measured from different 4D-MRIs. Synthetic 4D-MRIs overall showed a well-matched tumor motion with the original T1w 4D-MRI.
    Figure 1: (a) Training: reference DVFs were calculated between the starting frame and following frames and predicted DVFs were generated from the DL model with the two frames as inputs. The training loss was defined as the difference between the two sets of DVFs. (b) Application: appropriate starting frames were registered to other frames using the DL model, and the predicted DVFs warped the 3D MRIs to generate MRI frames at corresponding respiratory phases.
Back to Top
Digital Poster Session - MR-Guided Interventions: Methods
Engineering/Interventional/Safety
Thursday, 20 May 2021 19:00 - 20:00
  • Exploration of the Surgical Placement of the Local Pituitary Coil for Microadenomas
    Jiahao Lin1,2, Siyuan Liu2, Rock Hadley3, Marvin Bergsneider4, Giyarpuram N Prashant4, Sophie Peeters4, Robert Candler2,5, and Kyunghyun Sung1
    1Department of Radiological Sciences, David Geffen School of Medicine, University of California, Los Angeles, LOS ANGELES, CA, United States, 2Department of Electrical and Computer Engineering, University of California, Los Angeles, LOS ANGELES, CA, United States, 3Department of Radiology, University of Utah, Salt Lake City, UT, United States, 4Department of Neurosurgery, University of California, Los Angeles, LOS ANGELES, CA, United States, 5California NanoSystems Institute, Los Angeles, CA, United States
    We construct an agar phantom to examine the SNR variation as the coil plane rotates from 0° to 90°, respect to B0. Our local pituitary coil has improved SNR with a factor, ranged from 3.3 to 4.2, compared to the commercial Siemens head coil as long as the positioning of the coil to be within 0° and 70°.
    Figure 5, Local pituitary coil and Siemens head coil combined SNR, versus simulation SNR. Siemens commercial head coil SNR is not affected by the local coil angle.
    Figure 1, local pituitary coil placement against pituitary gland, and the illustration coil angle with respect to the scanner bed (B0 field).
  • Toward automatic lesion transmurality assessment using machine learning: a proof of concept in preclinical EP studies under MRI-guidance
    Valéry Ozenne1,2,3,4, Pierre Bour2,3,4, Marylène Delcey2,3,4, Nicolas Cedilnik5, Maxime Sermesant5, and Bruno Quesson2,3,4
    1Centre de Résonance Magnétique des Systèmes Biologiques, UMR 5536, CNRS, Bordeaux, France, 2IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Bordeaux, France, 3Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, Bordeaux, France, 4INSERM, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, Bordeaux, France, 5Université Côte d’Azur, Inria, Epione, Sophia Antipolis, France
    The purpose of this work is to evaluate the feasibility of automatic in-line segmentation in routine preclinical EP studies with applications for roadmap segmentation and lesion transmurality assessment
    Figure 1: Schematic view of the EP workflow under MRI guidance. A 3D bSSFP acquisition is used as a reference for the roadmap & EP mapping. The cavities are segmented and converted to meshes. Then catheter navigation is performed using active catheter tracking followed by EP mapping. RFA ablation is performed under MRI temperature monitoring followed by post-ablation imaging without contrast agent. Evaluation of the lesion size and transmurality is done on the post ablation images after epi/endo segmentation. The segmentation steps could benefit from machine learning.
    Figure 2: Specific features of the 3D images. Top: 3D bSSFP acquisition or 3D roadmap acquisition after catheter introduction into the cavities. The artefacts of the catheter (red arrows) are larger than the actual size of the catheter, in particular at the vicinity of the tip. The epicardial fat in hyper signal( blue arrows). Bottom: 3D long-TI acquisition performed after RFA. The created lesion can be visualized in hyper signal (yellow arrows) compared to the surrounding tissue attributed to edema (green arrows). The epicardial fat is highly visible in anterior and posterior view.
  • Towards Catheter-based Intra-Arterial Spin Labeling for Perfusion Measurements
    Kevin Waescher1, Simon Reiss1, Ali Caglar Özen1,2, Thomas Lottner1, and Michael Bock1
    1Dept. of Radiology, Medical Physics, Medical Center University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany, 2German Consortium for Translational Cancer Research Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
    In this study, we present in vitro results to demonstrate that arterial spin labelling using a local catheter based labelling coil provides an alternative technique for perfusion measurements during cardiovascular interventions without the use of exogenous contrast agents.
    Figure 2: Frequency dependence of the labelling signal fitted with a Lorentzian (a) along with the power dependence (b). Image of the catheter (c) and map of the relative signal difference created by the labelling (d). The difference image indicates the transmit B1 profile of the labelling coil.
    Figure 3: Image of the perfusion phantom setup (a) and maps of the relative signal difference between the labelled and unlabeled case for two different flow velocities (c and d). A plot of the normalized signal difference shows the exponential decay within the perfusion phantom.
  • Prototype Platform for Real-time MR-guided Brain Clot Evacuation
    Robert Moskwa1, Azam Ahmed2, and Walter Block1
    1Medical Physics, University of Wisconsin-Madison, Madison, WI, United States, 2Neurological Surgery, University of Wisconsin-Madison, Madison, WI, United States
    We present an MR-guidance prototype designed to provide neurosurgeons with an interface having similar viewpoints to those commonly used in stereotactic operating room settings. Flexible and interactive interface that leverages real-time MRI sequences. 
    Three displays, all individually acquired. Images are acquired in real time using a real-time Spiral GRE sequence. Each display can show a different view plane in real time by applying interactive tools. For example, a sample catheter trajectory is shown in red. The top-left axial display was rotated. The top-right view shows a perpendicular, but collinear display. The bottom image is perpendicular to the trajectory and shows the probe’s eye view. Parameters also adjustable.
    Three reference planes and a main display. Images are acquired in real time using a real-time Spiral GRE sequence. The main display can be changed to show a different view plane in real time by manipulating a graphic slice Rx on any of the three reference planes. For example, a sample plane holding the catheter trajectory is shown in red on the sagittal display in red. By selecting a perpendicular plane, shown in green, the probe’s eye view can be observed in real time. Parameters adjustable.
  • PETRA subtraction-based MRA to assess middle-cerebral-artery stenosis before and after treatment with angioplasty
    Feifei Zhang1, Yuncai Ran1, Shujian Li1, Jinxia Zhu2, Xuemei Gao1, Jingliang Cheng1, and Chengcheng Zhu3
    1The first Affiliated Hospital of Zhengzhou University, Zhengzhou, China, Zhengzhou, China, 2MR Collaboration, Siemens Healthcare Ltd., Beijing, China, Beijing, China, 3Department of Radiology, University of Washington, Seattle, Seattle, WA, United States
    PETRA-MRA is a promising non-invasive technique to evaluate MCA stenoses and could be an alternative technique for patient follow-up assessments after stent angioplasty.
    Fig1. Example results in a 63-year-old female patient who had left stent angioplasty in the left middle-cerebral-artery M1 segment. Preoperative (top row) and post-operative (bottom row) images are shown for time-of-flight magnetic resonance angiography (TOF-MRA, left), pointwise encoding time reduction with a radial acquisition (middle, PETRA-MRA), and digital subtraction angiography (right). Preoperative images indicated a similar degrees of stenosis (blue arrow). Obvious signal loss was observed in post-operative scans, with PETRA-MRA outperforming TOF-MRA.
    Table 1. Comparison of PETRA-MRA/TOF-MRA and DSA in measurements of stenosis of MCA
  • Deep Learning-Based Needle Tracking Trained on Bloch-Simulated Data and Evaluated on Clinical Real-Time bSSFP Images
    Ralf Vogel1,2, Dieter Ritter2, Jonathan Weine2,3, Jonas Faust2,4, Elodie Breton5, Julien Garnon5,6, Afshin Gangi5,6, Andreas Maier1, and Florian Maier2
    1Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany, 2Siemens Healthcare, Erlangen, Germany, 3TU Dortmund, Dortmund, Germany, 4Universität Heidelberg, Heidelberg, Germany, 5ICube UMR7357, University of Strasbourg, CNRS, FMTS, Strasbourg, France, 6Imagerie Interventionnelle, Hôpitaux Universitaires de Strasbourg, Strasbourg, France

    Synthetic training data was successfully used to train a U-Net for image-based needle tracking, utilizing automatically generated ground-truth labels. The evaluation on clinical patient data indicates that synthetic MR images can replace patient data and animal experiments.

    Figure 1: Needle predictions on clinical interventional images (red overlays), generated by U-net trained with synthetic data only.
    Figure 5: (a) random selection of simulated images generated using the AustinMan model. All images were simulated in the transversal plane, tilted by up to ±45° to the coronal plane. The needle was either placed in liver tissue from a reasonable angle or completely random (11, 12, 14, 15). (b) Image including the automatically generated needle label (green).
  • Deep Learning-based MR-only Radiation Therapy Planning for Head&Neck and Pelvis
    Florian Wiesinger1, Sandeep Kaushik1, Mathias Engström2, Mika Vogel1, Graeme McKinnon3, Maelene Lohezic1, Vanda Czipczer4, Bernadett Kolozsvári4, Borbála Deák-Karancsi4, Renáta Czabány4, Bence Gyalai4, Dorottya Hajnal4, Zsófia Karancsi4, Steven F. Petit5, Juan A. Hernandez Tamames5, Marta E. Capala5, Gerda M. Verduijn5, Jean-Paul Kleijnen5, Hazel Mccallum6, Ross Maxwell6, Jonathan J. Wyatt6, Rachel Pearson6, Katalin Hideghéty7, Emőke Borzasi7, Zsófia Együd7, Renáta Kószó7, Viktor Paczona7, Zoltán Végváry7, Suryanarayanan Kaushik3, Xinzeng Wang3, Cristina Cozzini1, and László Ruskó4
    1GE Healthcare, Munich, Germany, 2GE Healthcare, Stockholm, Sweden, 3GE Healthcare, Waukesha, WI, United States, 4GE Healthcare, Budapest, Hungary, 5Erasmus MC, Rotterdam, Netherlands, 6Newcastle University, Newcastle, United Kingdom, 7University of Szeged, Szeged, Hungary
    Deep Learning provides powerful tools to address unsolved problems and unmet needs of MR-only Radiation Therapy Planning (RTP) in terms of synthetic CT conversion (required for acquired dose calculation) and time-consuming organ-at-risk (OAR) delineation.
    Figure 1: In-phase ZTE without (top) and with 2xFOV extension plus DL image reconstruction (2nd row) for head&neck, together with corresponding DL derived synthetic CT (3rd row) and true CT (bottom row).
    Figure 3: 2D T2 PROPELLER based automated OAR segmentation in the head & neck (i.e. brain, brainstem, eyes, lens, optic nerves, chiasm, pituitary gland, cochlea, parotid glands, mandible, oral cavity, submandibular glands, larynx, spinal cord, and body contour).
  • Neural network based denoising of high temporal resolution cine images for tumor tracking in MR-guided radiotherapy
    Florian Friedrich1,2, Juliane Hörner-Rieber3, Peter Bachert1,2, Mark E. Ladd1,2,4, and Benjamin R. Knowles1
    1Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany, 2Department of Physics and Astronomy, Heidelberg University, Heidelberg, Germany, 3Department of Radiation Oncology, University Hospital of Heidelberg, Heidelberg, Germany, 4Faculty of Medicine, Heidelberg University, Heidelberg, Germany
    Image denoising and artifact suppression using U-net increases tumor tracking stability in undersampled Cartesian and radial cine acquisitions on MR-linac systems.
    Figure 2: Denoising of undersampled Cartesian images of a liver tumor (arrow).
    Figure 3: Denoising of undersampled radial images of a liver tumor (arrow).
  • Towards higher accuracy mapping of MRI to electron density using a 3D deep CNN for MRI-only radiotherapy treatment planning
    Jessica E Scholey1, Abhejit Rajagopal2, Elena Grace Vasquez3, Atchar Sudhyadhom4, and Peder Eric Zufall Larson2
    1Department of Radiation Oncology, University of California, San Francisco, San Francisco, CA, United States, 2Department of Radiology, University of California, San Francisco, San Francisco, CA, United States, 3Department of Physics, University of California, Berkeley, Berkeley, CA, United States, 4Department of Radiation Oncology, Harvard Medical School, Boston, MA, United States
    A 3D deep convolutional neural network was used to synthetize CTs (acquired at MV energies) from MRIs for more mapping of accurate electron density (versus those acquired at kV energies) for radiotherapy treatment planning
    Figure 2: Transverse images for four representative test datasets with their corresponding a) MRI, b) MVCT, c) sMVCT, and d) HU difference (sMVCT-MVCT) maps.
    Figure 3: Transverse images of a) dose distribution calculated on the sMVCT, b) the identical plan calculated on the actual MVCT, c) the difference in dose distributions, and d) the dose volume histograms for the target and several organs-at-risk. In panel (d), DVHs are shown as dotted and solid lines for the sMVCT plan and MVCT plan, respectively.
  • Multi-Task MR Simulation for Abdominal Radiation Treatment Planning: Technical Development
    Junzhou Chen1,2, Pei Han1,2, Fei Han3, Zhehao Hu1,2, Nan Wang1,2, Wensha Yang4, Anthony G Christodoulou1,2, Debiao Li1,2, and Zhaoyang Fan1,2,5
    1Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, CA, United States, 2Department of Bioengineering, University of California, Los Angeles, Los Angeles, CA, United States, 3Siemens Medical Solutions USA, Inc., Los Angeles, CA, United States, 4Department of Radiation Oncology, University of Southern California, Los Angeles, CA, United States, 5Department of Radiology, University of Southern California, Los Angeles, CA, United States
    We present a MR platform for radiation therapy planning in the abdomen using MR Multitasking, which is able to generate volumetric, multi-contrast, respiratory motion-resolved images under free-breathing with less than a 10 min scan time. 
    Fig.1 A) The pulse sequence diagram of the MR technique, with the red arrows being the navigator line locations. B) Illustration of the Cartesian spiral-in trajectory with Cartesian readout in Kx and spirals in Ky-Kz (PE-Partition) plane. C) Illustration of three-way tensor factorization to generate a multi-dimensional image 𝒜.
    Figure 3[GIF]. In-vivo MT multi-dimensional image from two volunteers. A user can view any contrasts, any motion states and any slice available in this image.
  • Image quality comparisons of novel and commercial coil setups in MRI for head and neck radiotherapy simulation
    Evangelia Kaza1, Jeffrey P Guenette2, Christian V Guthier1, Steven Hatch3, Alexander Marques3, Lisa Singer1, and Jonathan D Schoenfeld1
    1Radiation Oncology, Brigham and Women’s Hospital, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, United States, 2Division of Neuroradiology, Brigham and Women’s Hospital, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, United States, 3Radiation Oncology, Brigham and Women’s Hospital, Dana-Farber Cancer Institute, Boston, MA, United States
    A new UltraFlexLarge18 setup was more spacious and yielded higher SNR relative to a diagnostic head/neck coil than a commercially recommended FlexLarge4 coil setup for head images acquired in treatment position. Its clinical application could benefit head/neck radiotherapy MR simulations.
    Figure 1. a) A healthy volunteer in a radiotherapy immobilization mask fixed on a Qfix Portrait board placed on a modified Qfix Insight board. Two UFL18 coils attached by Velcro straps over and under the Insight board encompassed the volunteer’s head. b) The same volunteer and immobilization with the head surrounded by two attached FL4 coils placed in a Qfix Insight holder. c) The same volunteer in standard head/neck diagnostic position using a dedicated HN20 coil, where treatment masks do not fit.
    Figure 3. Graphs (top) and boxplots (bottom) of ratios of image quality parameters obtained using the flexible coil arrangements relative to the parameters obtained using a diagnostic HN20 coil, over all subjects and imaging sequences. From left to right: ratio of SNR, artifact size (AS) and CNR from the UFL18 (blue) and FL4 (red) coil images over SNR, AS and CNR from the HN20 coil images. ***: statistically significant difference (p < 0.001) between SNR_UFL18/SNR_HN20 and SNR_FL4/SNR_HN20.
  • Real-Time B0 Correction in a MRI Guided Radiotherapy System
    Austen Curcuru1, Deshan Yang1,2, and Michael Gach1,2,3
    1Biomedical Engineering, Washinton University in Saint Louis, Saint Louis, MO, United States, 2Radiation Oncology, Washinton University in Saint Louis, Saint Louis, MO, United States, 3Radiology, Washinton University in Saint Louis, Saint Louis, MO, United States
    Radiotherapy gantry rotation can cause artifacts and B0 fluctuations that impact tumor tracking and RT/MRI isocenter coincidence. Sequence parameters were adjusted at run time to reduce these issues by integrating a navigator into a bSSFP sequence and measuring B0 variations in real-time.
    Figure 2: Resonant frequency offset during gantry rotation starting at 30° and rotating around counterclockwise back to 33°. Offsets were calculated using a navigator inserted into a bSSFP sequence (red) and by measuring the peak offset in a FID sequence (blue).
    Figure 4: Illustration of the ViewRay MRI-Linac gantry showing the six mu-metal buckets and rotating passive shim tray. The 60° radial spacing between the buckets results in the periodic B0 offsets observed at different gantry angles.
  • Respiratory Motion Detection and Reconstruction Using CAPTURE and Deep Learning for a 0.35T MRI-LINAC System: An Initial Study
    Sihao Chen1, Cihat Eldeniz1, Weijie Gan1, Ulugbek Kamilov1, Deshan Yang1, Michael Gach1, and Hongyu An1
    1Washington University in St. Louis, Saint Louis, MO, United States
    A self-navigated respiratory motion detection (CAPTURE) and a deep learning 4D reconstruction method were used to derive the 3D deformable motion field for a 0.35 T MRI-LINAC system. Promising results were obtained despite the low SNR at 0.35 T. 
    Figure 3. Images reconstructed using non-MoCo, MCNUFFT, CS and P2P from a healthy subject. The top row shows phase 1, corresponding to the end-of-expiration, and the bottom row shows phase 5, corresponding to the end-of-inspiration. Blue arrows point to a region of comparison for image quality among the three methods.
    Figure 5. Images reconstructed using P2P at 0.35T and at 3T for different subjects for a rough qualitative assessment.
  • Joint Radial Trajectory Correction for Fast T2* Mapping on an MR-Linac
    Wajiha Bano1,2, Will Holmes1,2, Mohammad Golbabaee3, Alison Tree1,2, Uwe Oelfke1,2, and Andreas Wetscherek1,2
    1Joint Department of Physics, The Institute of Cancer Research, London, United Kingdom, 2The Royal Marsden NHS Foundation Trust, London, United Kingdom, 3Computer Science Department, The University of Bath, Bath, United Kingdom
    The joint gradient delay correction and T2* mapping approach outperforms the existing trajectory auto correction method. This will facilitate integration of T2* mapping for hypoxia imaging in an MR-linac treatment work flow.
    Figure 2: T2* maps obtained from fully sampled numerical phantom with the noise level of 0.1 for gradient delays [1 -1] (top) and [1 2] (bottom).
    Figure 4: T2* maps (Top) and T2*-weighted images (TE=5 ms) reconstructed from the in-vivo dataset without correction and with the sequential, respectively joint approach. Arrows highlight a banding artifact, which disappeared on one side after trajectory correction but partially remained on the other side.
  • Diffusion changes in prostate cancer patients undergoing radiation treatment on an MR-Linac system: Preliminary findings
    Colleen Bailey1,2, Rachel W Chan1, Jay Detsky3,4, Danny Vesprini3,4, and Angus Z Lau1,2
    1Odette Cancer Centre, Sunnybrook Research Institute, Toronto, ON, Canada, 2Medical Biophysics, University of Toronto, Toronto, ON, Canada, 3Department of Radiation Oncology, Sunnybrook Health Sciences Centre, Toronto, ON, Canada, 4Department of Radiation Oncology, University of Toronto, Toronto, ON, Canada
    Prostate cancer patients were scanned at five radiation treatment time points on an MR-Linac. High b-value imaging revealed visible lesions in 6/8 patients. Significant ADC increases were observed in 3/8 patients starting at the third treatment time point.
    Figure 3 Selected parameter maps with heterogeneous low ADC region identified as tumour (red arrow). The ADC and IVIM Dt maps have similar features and become less hypointense at later time points. The VERDICT intracellular fraction fI is higher in this region and decreases over time. The data are not always sufficient to fit the R parameter, particularly outside of the tumour ROI, resulting in dark sections in these maps.
    Figure 4 (a) Summary of parameter values in the tumour ROI over the course of treatment for one patient. The ADC and D­t are similar and increase at later treatment time points (* indicates statistically significant increase over Scan 1). The fp from IVIM and VERDICT parameters show more variability and do not demonstrate consistent or statistically significant changes, although the fI trend is inversely related to the ADC trend. (b) Changes in the ADC throughout treatment for all eight patients.
  • MRI-based tumor localisation after tantalum clip placement for proton beam therapy planning of uveal melanoma
    Myriam Jaarsma-Coes1,2, Teresa Ferreira1, Marina Marinkovic2, Khanh Vu2, Gre Luyten2, Coen Rasch3, Berit Verbist1, and Jan-Willem Beenakker1,2
    1Radiology, Leiden University Medical Center, Leiden, Netherlands, 2Ophthalmology, Leiden University Medical Center, Leiden, Netherlands, 3Radiotherapy, Leiden University Medical Center, Leiden, Netherlands
    A dedicated MRI protocol helps to improve proton beam radiation planning for uveal melanoma patients, as it provides a 3D visualisation of the radiographic marker-tumor relation compared to the conventional peroperative optical measurement.

    A,B). Coronal and sagittal view of the tumor are used to plan the scan to locate the clips.

    C) 2DT1 scan trough the base clearly show the clips and is used to plan the subsequent scans perpendicular to the tumor base and through centre of the tumour and centre of clip C/D.

    D,E). T2 spin echo sequence and T1 gradient echo showing the tumour (dagger), retinal detachment (double dagger) and clip (arrow). The clip location is better visualized on the gradient echo sequence, while a more accurate measurement can be made on the spin echo sequence.

    F) T2 image of clip D.

    Summary of the differences in clip-tumor distances between the optical peroperative measurement and the MR-measurement. For 51% of the clips no clinically significant difference was found.
  • Deep learning based synthetic CT skull for transcranial MRgFUS interventions using 3D V-net–Transfer learning implications
    Pan Su1,2, Sijia Guo2,3, Steven Roys2,3, Florian Maier4, Thomas Benkert4, Himanshu Bhat1, Elias R. Melhem2, Dheeraj Gandhi2, Rao P. Gullapalli2,3, and Jiachen Zhuo2,3
    1Siemens Medical Solutions USA, Inc., Malvern, PA, United States, 2Department of Diagnostic Radiology and Nuclear Medicine, University of Maryland, School of Medicine, Baltimore, MD, United States, 3Center for Metabolic Imaging and Therapeutics (CMIT), University of Maryland Medical Center, Baltimore, MD, United States, 4Siemens Healthcare GmbH, Erlangen, Germany
    Deep learning can be used to generate synthetic CT skull, thereby simplifying workflow of tcMRgFUS. 3D V-Net can take advantage of contextual information in volumetric images. Furthermore, the pre-trained model can be applied in dataset acquired using different sequence/protocols.
    Figure 4: 3D cinematic rendering of the reference CT skull and synthetic CT skull from a representative subject.
    Figure 5: Transfer learning for 0.8mm3 isotropic synthetic CT skull: results of applying the trained model above to dual echo spiral UTE data (0.8mm3 isotropic) acquired on different 3T scanner. ab) dual echo spiral UTE; c) reference CT skull; d) results of directly applying the previous model trained on radial UTE without any retraining of the spiral UTE data; e) results of applying the model after retraining with single spiral UTE data. Bottom are zoom-in windows showing the details in posterior skulls. Note that all CTs here are 0.8mm3 isotropic, as CT were coregistered to UTE space.
  • FMRI of Post High-Frequency Focused Ultrasound Ablation of ViM Shows Reduced Ipsilateral Thalamic and Cortical Motor Activation
    Anna Crawford1, Mark Lowe1, Sean Nagel2, Daniel Lockwood1, Emmanuel Obusez1, Andre Machado2, and Stephen Jones1
    1Imaging Institute, Cleveland Clinic Foundation, Cleveland, OH, United States, 2Neurological Institute, Cleveland Clinic Foundation, Cleveland, OH, United States
    HIFU treatment targets the thalamus using measurements and landmarks. We explore an alternative method using 7T task-related fMRI and present preliminary data of BOLD activation in patients before and after HIFU treatment.
    Figure 2: Mean statistical difference maps of pre-HIFU minus post-HIFU, using a finger tapping fMRI paradigm. Arrows show regions with decreased activation after HIFU: green, motor cortex in precentral gyrus; red, motor regions of thalamus; blue, subthalamic nucleus; magenta, globus pallidus; orange, putamen. Threshold for top panel is p=0.07; for bottom panel is p=0.02.
    Figure 1: Single subject sample of finger tapping time series data from a voxel in the right motor cortex. The maximum and minimum values as indicated by red dots are 1396 and 1252 respectively.
  • Comparison of Automated Thalamic Segmentation Techniques: Applications in MRgFUS Planning
    Kain Kyle1, Jerome Maller2, Yael Barnett3, Stephen Tisch3, Benjamin Jonker3, Michael Barnett1, Arkiev D'Souza1, and Chenyu Wang1
    1University of Sydney, Sydney, Australia, 2GE Healthcare, Sydney, Australia, 3St Vincent's Hospital Sydney, Sydney, Australia
    In this investigation the clinical utility of thalamic segmentation tools FreeSurfer and THOMAS for use in MRgFUS planning.  We demonstrate that the ventral lateral nuclei segmented by FreeSurfer and THOMAS are consistent with the MRgFUS ablation target chosen by the treating clinician. 
    Figure 2. Comparison of typical Freesurfer and THOMAS thalamic segmentations. TopFreesurfer thalamic segmentation VLp (red) and VLa(green). Bottom THOMAS WMnMPRAGE thalamic segmentation, VLp (yellow) and VLa (blue). Ablation segmentation is overlayed in white.
    Figure 1. Volume of ablation in the FreeSurfer VLa, VLp and whole Thalamus, as a percentage of the total ablation volume.
  • Reduced-FOV 3D MR-ARFI with a joint reconstruction for localizing the focused ultrasound beam in neuromodulation
    Huiwen Luo1,2, Michelle K. Sigona1,2, Li Min Chen2,3, Charles F. Caskey2,3, and William A. Grissom1,2,3
    1Biomedical Engineering, Vanderbilt University, Nashville, TN, United States, 2Vanderbilt University Institute of Imaging Science, Nashville, TN, United States, 3Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States
    MR-ARFI can encode the ultrasound-induced displacements into the phase of an MR image. A two-minute reduced-FOV 3D MR-ARFI scan with a joint image reconstruction method at 3 Tesla was demonstrated to image and localize the entire focus in FUS neuromodulation with a low FUS duty-cycle of 0.85%.
    Figure 2. (a) 3D MR-ARFI sequence overview; (b) The proposed k-space undersampling scheme.

    Figure 4. Representative reconstructed magnitude images (FUS-ON, positive MEGs) (a) and corresponding ARF displacement maps (b).

Back to Top
Digital Poster Session - Multimodal Imaging & Auxiliary Devices
Engineering/Interventional/Safety
Thursday, 20 May 2021 19:00 - 20:00
  • A compact and clonable ultrasound-based sensor system to monitor physiological motion
    Bruno Madore1, Cheng-Chieh Cheng2, and Frank Preiswerk3
    1Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, United States, 2Computer Science and Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan, 3Amazon Robotics, Westborough, MA, United States
    An ultrasound-based sensor system was built that is compact, clonable and GUI-controlled. A key strength of these sensors is that they attach to the skin and can follow a patient through different procedures, possibly linking them in ways that take physiological motion into account.
    Fig. 3: The sensor system is controlled through a GUI. Its initial state is shown in (a), and a snapshot of a real-time acquisition in shown in (b). Instead of raw data (b), one can select velocity (c) or displacement (d) for the GUI’s display window, and switch between display types during the real-time acquisition if so desired. Raw data are saved to file and can be retrospectively re-reconstructed and displayed, which happens to be how (c) and (d) were generated; note the gold color of the selected ‘retro’ button in (c) in (d), as opposed to its pale-blue non-selected color of (a) and (b).
    Fig 1: a) A compact, clonable sensor system was developed for ultrasound-based sensors, to monitor physiological motion during diagnostic imaging. Up to four sensors can be used at a time (red arrow), and an RF antenna can be used for synchronization with MRI (green arrow). While this system is still one-of-a-kind, it was built to be clonable, to enable possible collaborations across institutions and geographical locations.
  • Acoustically transparent and low-profile head coil for high precision magnetic resonance guided focused ultrasound at 3 T
    Isabelle Saniour1, Fraser Robb2, Victor Taracila2, Jana Vincent2, Henning U. Voss1, Michael G. Kaplitt3, J. Levi Chazen1, and Simone Angela Winkler1
    1Weill Cornell Medicine, New York, NY, United States, 2MR Engineering, GE Healthcare, Aurora, OH, United States, 3Department of Neurological Surgery, Weill Cornell Medicine, New York, NY, United States
     In vivo MR results with a very thin, ultra-flexible and acoustic transparent head coil showed increase of the SNR over the body coil by a factor of 7.3 and 7.6 in a brain image with and without the presence of the transducer, respectively for MRgFUS procedure.
    Figure 3. a) Normalized radial acoustic intensity magnitude for 1MHz (left); 650kHz (center); 220kHz (right) along r-coordinate passing the focal point. b) 2D map of the intensity magnitude showing the effects of the RF coil on the magnitude and the focal point (1MHz (left); 650kHz (center); 220kHz (right). For each frequency, a magnified view of the focal point is displayed in the dotted rectangle. For the frequency of 220kHz, a magnified color bar is displayed. First row : With FUS-Flex coil present at 20 mm ; Second row : at 80 mm ; Third row : reference without coil .
    Figure 1. a) Illustration of a patient positioned inside a focused-ultrasound transducer. b) Photograph of the FUS-Flex coil.
  • A dedicated 3-ch breast coil for Microwave Tomography at 3T
    Xiaoyu Yang1, Thomas Eastlake1, Tsinghua Zheng1, Shireen Geimer2, Timothy Raynolds2, Paul Meaney2, and Keith Paulsen2
    1Innovation, Quality Electrodynamics, LLC, Mayfield Village, OH, United States, 2Thayer School of Engineering, Dartmouth College, Hanover, NH, United States
    Combining MR and Microwave Tomography (MT) for breast imaging may provide high sensitivity with an improve specificity. There is no devoted MR coil for the MT application.  A dedicated 3-ch MR coil is proposed for superior MRI imaging quality that enables the future MT MR volunteer evaluation.
    Figure 2. the 3-ch coil diagrams and photo
    Figure 3. The matching/feeding schematic
  • A Dedicated Multichannel Head Coil Array for PET Insert on 3 T MRI
    Jo Lee1,2, Ziru Sang1,2, Yongfeng Yang1,2, Xiaoliang Zhang3, and Ye Li1,2
    1Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China, 2Shenzhen Key Laboratory for MRI, Shenzhen, China, 3Department of Biomedical Engineering, State University of New York, Buffalo, NY, United States
    A quadrature birdcage/47Rx head coil array designed and built for a dedicated PET insert was presented and tested on a 3T MRI system. The coil performances have been verified through flip-angle maps, SNR, and anatomical images.
    Fig.5. 0.6 mm isotropic modulated flip angle technique in refocused imaging with extended echo train (MATRIX) images were acquired using the quadrature birdcage/47Rx head coil array at R = 7-fold acceleration. The acceleration direction was on phase encoding and slice phase encoding. H = head; P = posterior; R = right; A = anterior.
    Fig.1. The (a.) application scenario and (b.) coil wiring of the quadrature birdcage/47Rx head coil array. The birdcage's size was displayed in red, the coverage of the surface coil was displayed in blue, and the size of a single surface loop was displayed in yellow.
  • Vector Modulator Based Automated Active Compensation of Direct Feedthrough in Magnetic Particle Imaging
    Bilal Tasdelen1,2, Mustafa Utkur1,2, Ergin Atalar1,2, and Emine Ulku Saritas1,2
    1Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey, 2National Magnetic Resonance Research Center (UMRAM), Ankara, Turkey
    We propose a vector modulator based active compensation method and a lookup-table-based algorithm for reducing the direct feedthrough interference by more than 40 dB within two seconds. We demonstrate a significant increase in detection sensitivity.
    Figure 3: (A) Boxplot of the simulation results showing the performance of the algorithm with and without measurement noise added. There is a high correlation (ρ=0.45) between the final cancellation and the initial feedthrough. (B) Monitoring the feedthrough level over 20 seconds, while heating the coil. Green dashed boxes mark the time intervals where the algorithm was running. The threshold for running the algorithm was -50 dBV. Different initial isolations are set up by adjusting the gradiometer knob.
    Figure 1: (A) Diagram of the AWR (bottom) and vector modulator (up). (B) AWR coil (bottom) and the small pickup coil positioned on top. (C) Vector modulator. An Arduino Uno Rev 3 was used to establish the communication between the rheostats and the computer where the experiments were orchestrated. (D) Photo of the overall system.
  • Interventional Magnetic Particle Imaging Open-bore Scanner Design
    Michael Seeg1, Martin Rueckert1, Stefan Herz2, Thomas Kampf1,3, Thorsten Bley2, Volker Behr1, and Patrick Vogel1
    1Experimental Physics V, Julius-Maximilians-University, Wuerzburg, Germany, 2Department of Diagnostic and Interventional Radiology, University Hospital, Wuerzburg, Germany, 3Department of Diagnostic and Interventional Neuroradiology, University Hospital, Wuerzburg, Germany
    An open bore magnetic particle imaging scanner concept for interventional treatments is presented. As a radiation-free method, based on magnetic particle imaging with a field free line, it may potentially help reducing the health risk for medical staff and patients.
    Current setup of the scanner device. Additionally a receive coil is inserted with a diameter of 136 mm.
    Scanner design with only the drive coils shown. In z-direction an FFL can be seen, which can be guided through the scanner along the symmetric axis. Simulated with software for particle imaging.19
  • MR safe electromagnetic direct current motor
    Lorne W Hofstetter1, J Rock Hadley1, Robb Merrill1, and Dennis L Parker1
    1Department of Radiology and Imaging Sciences, University of Utah, Salt Lake City, UT, United States
    A novel direct current electromagnetic motor design that can safely operate inside the MR scanner bore was developed, tested, and characterized. This new actuator technology may help enable robotic assisted interventions under MR-guidance.
    Figure 1. Diagram of a standard permanent magnet DC motor in (a) and an MR-compatible electromagnetic motor in (b). The magnetic components in (a) are labeled in red. The corresponding functional elements for the MR-compatible motor are labeled as blue in (b).
    Figure 2. Images of MR safe DC motor and assembly are shown. The major motor components are shown in (a). An exploded view of the rotor, commutators, and end cap assembly is shown in (b). The assembled motor is shown in (c) with the twisted and shielded power supply cables shown in (d).
  • Real-time electric field estimation in transcranial magnetic stimulation using deep learning and magnetic resonance imaging
    Guoping Xu1,2, Yogesh Rathi2,3, Joan A Camprodon3,4, and Lipeng Ning2,3
    1Wuhan Institute of Technology, Wuhan, China, 2Brigham and Women's Hospital, Boston, MA, United States, 3Harvard Medical School, Boston, MA, United States, 4Massachusetts General Hospital, Boston, MA, United States
    We proposed a deep-neural-network based approach for real-time prediction of TMS-evoked E-field using subject specific MRI. The predicted E-field is similar to the result estimated using finite element methods.
    Fig. 2: The magnitude [V/m] of E-field from VCM (the first column) and the proposed method (the middle column), and the absolute difference between VCM and the proposed method (the last column).
    Fig. 3: Evaluation results on brain surface: (a) target overlapping coefficient, (b) E-field peak distance, (c) E-filed correlation coefficient, (d) maximum absolute error.
  • A wearable “RF-EEG Cap” for full head coverage concurrent TMS/EEG/fMRI experiments at 3T: a feasibility study
    Lucia Navarro de Lara1,2, Padmavathi Sundaram 1,2, Lena Nohava3, Elmar Laistler3, Mohammad Daneshzand1,2, Lawrence Wald1,2, Jason Stockmann1,2, and Aapo Nummenmaa1,2
    1Martinos Center - MGH, Charlestown, MA, United States, 2Harvard Medical School, Boston, MA, United States, 3High Field MR Center, Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, Vienna, Austria
    A 2-channel “RF-EEG Cap” prototype for 3T was successfully built using flexible coaxial technology to allow concurrent TMS/fMRI/EEG applications. The effects on the SNR, B0, EEG data quality and EPI image were found minimal.
    Fig2.
    Fig5.
  • Quantification of HD-EEG net attenuation on 18F-FDG static PET: a phantom and human brain study
    Erica Silvestri1,2,3, Marco Castellaro1,4, Alessandra Zorz5, Andrea Bettinelli 1,5, Cristina Campi6, Marta Paiusco5, Diego Cecchin2,7, and Alessandra Bertoldo1,2
    1Department of Information Engineering, University of Padova, Padova, Italy, 2Padova Neuroscience Center, University of Padova, Padova, Italy, 3Department of Neuroscience, University of Padova, Padova, Italy, 4Department of Neurosciences, Biomedicine and Movement Sciences, University of Verona, Verona, Italy, 5Department of Medical Physics, Veneto Institute of Oncology IOV - IRCCS, Padova, Italy, 6Department of Math, University of Padova, Padova, Italy, 7Department of Medicine, Unit of Nuclear Medicine, University of Padova, Padova, Italy
    HD-EEG net attenuates PET signal in NAC images in both phantom and human brain, respectively of 3.01% (±3.69%) and of 0.53% (±5.89%). AC images show an opposite behaviour: attenuation correction factor tends to be overestimated in CT- and UTE-derived UMAP obtained with HD-EEG net.
    Figure 3: Spatial pattern of relative % differences between PET images acquired with and without HD-EEG net (blue-red scale) over-imposed on phantom CT image or human brain T1w MPRAGE image (grey scale).
    Figure 1: Schematic representation of the two acquisition protocols.
  • Triggering pulse generation by non-contact monitoring of heart rate using an antenna with 600 MHz resonance
    Kazuya Okamoto1, Takafumi Ohishi2, Sojuro Kato1, Ryuichi Nanaumi3, Kenji Oyama3, and Kazuhiko Fukutani3
    1Advanced MRI Development PJ Team, Canon Medical Systems Corp., Kawasaki, Japan, 2Advanced Technology Research Department, Research and Development Center, Canon Medical Systems Corp., Kawasaki, Japan, 3Medical Products Technology Development Center, R&D Headquarters, Canon Inc., Tokyo, Japan
       New non-contact cardiac and respiratory gating method have been developed, which includes an antenna monitoring heart and respiration movement, and independent triggering pulse generation unit with signal processing software.
    Figure 2:Triggering pulse generation unit.
    Figure 3: Waveforms of ECG (blue), the antenna signals before IIR filter (orange) and the IIR filtered antenna signals (green). The inflection points of the IIR waveform and the R-wave in ECG are shown in green and blue circles, respectively
  • MRI data transmission via fifth-generation (5G) cellular networks
    Nicolai Spicher1, Ramon Barakat1, and Thomas Martin Deserno1
    1Peter L. Reichertz Institute for Medical Informatics, TU Braunschweig and Hannover Medical School, Braunschweig, Germany
    We made use of fifth-generation (5G) cellular networks in the current 3GPP Release 15 for transmitting DICOM data and reached more than 50 Mbit/s upload speed using consumer-grade hard- and software. This holds great potential for cellular transmission of data from portable MRI scanners.
    Figure 1: Experimental setup. The laptop is connected via ethernet cable to the 5G router and displays its web-interface.
    Table 1: Results of experiments. Transmission times are given as mean ± standard deviation. The reference measurement was perfomed via ethernet connection between the DICOM client and server, therefore quality quantities of cellular networks are not available.
  • Test Bed for High Speed Serial Wireless MRI Data Studies
    Audrey Chan1, Fraser Robb2, John Pauly3, Shreyas Vasanawala4, and Greig Cameron Scott5
    1Victoria University of Wellington, Wellington, New Zealand, 2GE Healthcare, Aurora, OH, United States, 3Stanford University, Stanford, CA, United States, 4Radiology, Stanford University, Stanford, CA, United States, 5Electrical Engineering, Stanford University, Stanford, CA, United States
    A test system is prototyped to perform pseudo-random binary sequence (PRBS) generation and clock/data recovery to assess high speed serial link performance of wireless and optical signal integrity for MRI data interfacing. 
    Figure 1: Clock and Data Recovery board prototypes. a) Micrel SY87701 can perform clock & data recovery from 28Mbps to 1.3Gbps, but requires a reference clock. b) Analog Devices ADN2817 operates from 10Mbps to 2.7Gbps, integrates a front end limiting amplifier for small signals, and has a programmable pseudo-random binary sequence (PRBS) generator.
    Figure 2: a) PRBS generator uses an Arduino module to configure the ADN2817 CDR device as a PRBS generator. The clock is provided by an SiLabs Si571 programmable clock generator eval board. b) Demonstration of PRBS generation from 400-800 Mbps to a 3.2GHz amplitude shifted keyed modulator and direct detection by an ADL6012 envelope detector. The detector output is yellow, and clock/data (pink/blue) is recovered by the SY87701. The Si571 clock reference is the green trace.
  • Wireless Powered Frequency Encoding of Locally Detected MR Images for Remote Signal Transmission
    Wei Qian1 and Chunqi Qian1
    1Michigan State University, East Lansing, MI, United States
    Wireless MR signal transmission by a compact and wirelessly powered parametric oscillator.
    (a) A double frequency parametric resonator included a rectangular conductor with a center conductor. Its top and bottom gaps were filled by variable capacitors. (b) A single frequency resonator with the same dimension was tuned by two capacitors. (c) The two resonators were overlapped through a substrate. (d) The modulator was placed on top of a phantom containing 1% agarose, before being inserted into the center of a 7T magnet. During MR signal acquisition, pumping power was applied on the antenna to activate the modulator whose oscillation signal was received by the volume coil.
    When the relative SNR of MR images was measured at different distance separations, the amplified resonator (gray) could retain constant sensitivity for transmission attenuation up to 21 dB (at 4.5 cm separation). In comparison, the frequency modulator (yellow) could retain constant sensitivity for transmission attenuation up to 34 dB (at 11 cm separation).
  • Single board computer as a satellite-linked, deep learning capable pocket MR workstation: a feasibility study
    Keerthi Sravan Ravi1,2, John Thomas Vaughan Jr.2, and Sairam Geethanath2
    1Biomedical Engineering, Columbia University, New York, NY, United States, 2Columbia Magnetic Resonance Research Center, New York, NY, United States
    We implemented an optimized Tensorflow-Lite segmentation model on Raspberry Pi (RPi) reducing 75% disk space. Data was uploaded across three internet access methods including a satellite link. We also ran a DICOM viewer on RPi.
    Figure 5. Comparing mobile and fixed workstations across six dimensions. The Raspberry Pi based mobile workstation (this work, first column) is the most affordable and most portable. *An example of a commercially available portable MRI workstation is from Hyperfine (https://hyperfine.io/).
    Figure 2. Illustration of the three hardware setups discussed in this work. The Raspberry Pi was used in combination with keyboard/mouse and display monitor input/output peripherals. The three Internet connectivity modalities studied in this work are fixed broadband (WiFi), cellular broadband (3G/LTE) and satellite-link. Picture at the bottom is an example of the hardware setup highlighted in yellow.
  • STAR MRI With Non-Magnetic, Integrated Circulator based on Switched Transmission Lines
    Aravind Nagulu1, Ahmed Kord1, Gehua Tong1, Michael Garwood2, Lance DelaBarre2, Djaudat Idiyatullin2, SungMin Sohn3, J. Thomas Vaughan1, and Harish Krishnaswamy1
    1Columbia University, New York, NY, United States, 2University of Minnesota, Minneapolis, MN, United States, 3Arizona State Univeristy, Tempe, AZ, United States
    Our circulator which is designed to operate with a 1.5 Tesla eight-channel TEM head coil and achieves a transmission loss of 6.1dB and 5.8dB in the TX and the RX paths respectively, while achieving a TX-to-RX isolation of 75dB at 64MHz.
    Figure 1: Concept diagram of the (a) switched transmission line gyrator and (b) 3-port circulator.
    Figure 5: Measured TX-to-RX leakage when the circulator is terminated with 3 of the 8 channels of a 1.5 T MR coil. The dotted red line depicts the isolation presented by the circulator when the MR coil is loaded with a human head.
  • A 5.7mW 1.5T Silicon Germanium Pre-Amp, and its Effects on Image Distortion
    Christopher Vassos1, Fraser Robb2, Shreyas Vasanawala3, John Pauly1, and Greig Scott1
    1Electrical Engineering, Stanford University, Stanford, CA, United States, 2GE Healthcare, Aurora, OH, United States, 3Radiology, Stanford University, Stanford, CA, United States
    A low power (5.7mW) Silicon Germanium Amplifier was constructed and its effects on image distortion shown via benchtop measurement. Image distortions from device non-linearity began at peak input powers of-25dBm. Reducing barriers to wireless receive arrays by reducing power consumption. 
    Amplifier Schematic. The BFP840 acts as the primary amplification device and the BFP183 acts as a cascode device.
    Magnitude and Difference Images (x10). Peak input powers are (left to right): -40dBm, -30dBm, -25dBm, -20dBm. Difference images taken with regard to the reference dataset. A significant change in image contrast is evident upon reaching input powers of –25dBm, further evident at –20dBm. .
  • Hybrid Pair Ratio Adjustable Power Splitter using off-shelf components and easy-to-replace microstrip phase shifter
    Yue Zhu1,2, William A Grissom1,2,3,4, John C Gore1,2,3,4, and Xinqiang Yan1,2
    1Vanderbilt University Institute of Imaging Science, Nashville, TN, United States, 2Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN, United States, 3Biomedical Engineering, Vanderbilt University, Nashville, TN, United States, 4Electrical and Computer Engineering, Vanderbilt University, Nashville, TN, United States
    A low loss ratio adjustable power splitter (RAPS) was designed with a Wilkinson Splitter. It has potential problems of challenging and laborious building process and safety hazard. In this abstract, a new RAPS designed with hybrid couplers and printed microstrip lines mitigated the problems. 
    Figure 1. a) Diagram of the power splitter model. The final ratio between Vout1 and Vout2 are adjusted by varying the microstrip line. b) Transmission matrices for the hybrid coupler (Thc), the transmission line (Ttm), and the RAPS (TRAPS). The scattering matrix of RAPS (SRAPS) is derived from TRAPS. We can see from the S-parameter matrix, for any ratio, HP-RAPS is a lossless circuit with all ports well matched and output ports isolated. The output ratio has a tangent dependence on the length difference between the two microstrip lines (S13/S14 = tan[(φ1φ2)/2]).
    Figure 5. Bench tests of the boards with designed magnitude ratios of a) 1:2, b) 1:4, and c) 1:8. The matching and isolation are similar to that of the 1:1 ratio board. The tested magnitude ratios are 1:2.4, 1:4.8, and 1:6.9, respectively.
  • High Speed Serial ASK Signal Integrity for Wireless MRI
    Greig Cameron Scott1, Audrey Chan2, Wonje Lee3, Fraser Robb4, John Pauly5, and Shreyas Vasanawala6
    1Electrical Engineering, Stanford University, Stanford, CA, United States, 2Victoria University of Wellington, Wellington, New Zealand, 3Pediatric Radiology, Stanford University, Stanford, CA, United States, 4GE Healthcare, Aurora, OH, United States, 5Stanford University, Stanford, CA, United States, 6Radiology, Stanford University, Stanford, CA, United States
    Antenna components and amplitude shift keyed (ASK) modulation is assessed for very short range high speed data links.  Prototypes are developed using serializer/deserializer components to achieve rates of 600 Mbps and signal integrity is evaluated.
    Figure 4: Mock bore tests using a pair of 5 GHz Biquad antennas at 200, 500 and 600 Mbps show progressive signal degradation of the demodulated bit pattern and eye diagrams, though signal quality was sufficient for clock recovery.
    Figure 5: Mock Bore multi-channel deployment and spectrum. Four ASK modulators at 200 Mbps with 3.2, 3.8, 4.4 and 5.0 GHz carriers and Biquad antennas are deployed. A single Taoglas UWB is used for signal reception. The spectrum shows the SINC-like nulls at 200 MHz offsets, and how further signal lobes would over-run to neighboring bands as adjacent channel interference.
  • A Sip in the Bore: How to Make Coffee with an MRI System
    Michael Bock1, Thomas Lottner1, and Ali Caglar Özen1,2
    1Radiology - Medical Physics, University Medical Center, Medical Faculty, University Freiburg, Freiburg, Germany, 2German Consortium for Translational Cancer Research (DKTK), Freiburg, Germany

    Energy harvesting from MRI gradients can provide power levels on the order of a kilowatt.

    Figure 1: Left: The induction coil (yellow) is being connected to a modified drip coffee machine on the patient table of a 3T MRI system. (Right) Coffee was ready at the rear end of the magnet about 5 min after the start of a pulse sequence with an oscillating z-gradient field.