Gadgetron implementation for PET/MR motion correction

Transcription

1 implementation for motion correction Thomas Küstner School of Biomedical Engineering & Imaging Sciences, King s College London, St. Thomas Hospital, London, UK Dept. of Diagnostic and Interventional Radiology, University Hospital Tübingen, Tübingen, Germany Institute of Signal Processing and System Theory, University of Stuttgart, Stuttgart, Germany CCP PET-MR Software Meeting

2 Motivation Aim: capture and correction motion to improve quality in clinical feasible setting MR uncorrected corrected PET

3 Motivation Aim: capture and correction motion to improve quality in clinical feasible setting Küstner et al. Med Image Anal 2017 Küstner et al. ISMRM 2015, 2016, 2017 MR PET estimate motion motion model correct motion uncorrected corrected

4 data accumul. k-space splitting image accumul. image slicer MR PET Motion Correction Pipeline PET MR Sensors Motion Model 90 s further diagnostic sequences sensor fusion listmode ECG, Belt,... 4D k-space vendor pipeline Sensor fusion raw data emitter Reconstruction and Motion Correction raw data injector Sensor vendor pipeline image emitter image injector CASToR listmode-based PET correction and recon Host PC TCP/IP

5 [1] Hansen and Sorensen MRM 2013 Integration into scanner system 1 workstation MRIR client-server architecture

6 Communication connection via TCP/IP VB20P: TCP tunnel VE11P: network switch style-sheet forwarded as config file fully controlled by WIP parameters in sequence

7 ICE integration emitter: send out raw data injector: receive reconstructed image

8 ICE integration emitter: Scan injector: send out raw data receive raw data

9 CSMoCoLab Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition µmap processing PET reconstruction image slicer docker pull kspaceastronauts/cs_moco_lab

10 CSMoCoLab MR side Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition buffer data into image and navigator to process full k-space wait for ACQ_END flag since VE11P: not always properly set and detected by workaround: overwrite close() method to perform handover

11 CSMoCoLab MR side Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition PCA-based Power-based sensors (camera, belt, ) sensor fusion respiratory cardiac respiration extraction from N-D self-navigation signal cardiac motion from N-D self-navigation signal or ECG PCA-based extraction Power-based extraction

12 CSMoCoLab MR side Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition respiratory gating: equidistant or equicount cardiac gating: linear or diastolic view-sharing amongst gates Gaussian or linear soft-weighting

13 [1] Uecker ISMRM 2015 [2] Küstner ICASSP 2015 CSMoCoLab MR side Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition FOCUSS BART 1 image motion model Compressed Sensing reconstruction of subsampled and gated k-space include motion correction 2

14 CSMoCoLab Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition Local All-Pass 1,2 elastix 3 motion model non-rigid image registration elastix: cubic B-Spline Local All-Pass: optical flow (global non-rigid deformation modelled as local rigid displacement) motion model: send back to scanner as 3 measurement DICOMs [1] Gilliam ISBI 2016 [2] Küstner ISMRM MoCo 2017 [3] Klein TMI 2010

15 CSMoCoLab MR side Scanner for raw data Reconstruction pipeline - MRIR emitter functor injector functor for postprocessing Host Accumulator Navigation Gating & Population Reconstruction Registration MRIImage- ToAcquisition convert back to (now fully sampled) acquired k-space data header information of one acquisition stored in a Singleton class (global access)

16 CSMoCoLab PET side Scanner for raw data scanner µmap µmap processing listmode Reconstruction pipeline - MRIR emitter functor PET reconstruction injector functor for postprocessing image slicer Host scanner µmap recon. motionresolved MR motion model deform 1-Point DIXON µmap motion-state matching µmaps deform scanner generated µmap (2-Point DIXON scan) with derived MR motion model extract from motion-resolved MR image (1-Point DIXON)

17 CSMoCoLab PET side Scanner for raw data scanner µmap µmap processing listmode Reconstruction pipeline - MRIR emitter functor PET reconstruction injector functor for postprocessing image slicer Host listmode motion model e7 tools CASToR motion-corrected PET motion-corrected PET reconstruction e7 tools (Siemens): image-based motion correction CASToR: listmode-based motion correction

18 Control manage your configs and Gadgets set up config parametrize it run config store output

19 liver 2D MIP spatial 3D thorax temporal Results 66y melanoma patient with liver metastasis PET corrected uncorrected MR resp.+ card. resp. only cardiac only MR: TA=90s, TE=1.23ms, TR=2.6ms, 1.9mm 3 iso. PET: OPL-EM, 7 iter., 7 subsets, 4mm Gaussian

20 PET MR Conclusions capture and track motion under free movement clinical feasible setup minimal blocking of MR for MC with 100% scan efficiency sensor fusion for full PET scan time coverage streamlined processing in improve PET quanitification and detection mean improvements in a cohort of 37 patients: SUVmax: 25% SNR: 10% contrast: 26% FWHM: 28% Slope: 99% uncorrected corrected

21 Acknowledgements sites.google.com/site/ kspaceastronauts/ Küstner: implementation for motion correction