Deep learning for MR imaging and analysis
Shanshan Wang Paul C. Lauterbur Research Center for Biomedical Imaging Shenzhen Institutes of Advanced Technology (SIAT)
2020-04-15
Deep learning for MR imaging and analysis Shanshan Wang Paul C. - - PowerPoint PPT Presentation
Deep learning for MR imaging and analysis Shanshan Wang Paul C. Lauterbur Research Center for Biomedical Imaging Shenzhen Institutes of Advanced Technology (SIAT) 2020-04-15 Learning Reconstruction and Analysis Image reconstruction a.
2020-04-15
a. Stroke lesion segmentation b. Breast tumor classification and segmentation c. Cervical cancer classification
a. Summary b. Linear reconstruction c. Non-linear iterative reconstruction d. Deep learning-based reconstruction a. Background
Correlated &interacting
Short imaging time may cause issues like low resolution and SNR, while long time may cause issues like claustrophobia, motion artifact and signal distortion
MR Image Sampling Data (K-Space) Recontruction
Scan time acquisition lines repetition
➢
MR Physics (1970’s)
➢ Hardware (2000’s)
➢
Image reconstruction from incomplete K-space data (past decade)
Encoding Reconstruction Radio Frequency Pulse Phased Array Coil K-Space Data Image for diagnosis
Sanity (relation to measurements) Prior or regularization y : Given measurements x : Unknown to be recovered E : Encoding matrix Pr(x)
Incomplete K-space data Reconstructed image
Sparsity domain
Big dataset collection Deep prior learning
Big datasets 1 st phase linear reconstruction (IFFT, SENSE, SMASH, GRAPPA, ….) 2 nd phase Nonlinear iterative reconstruction (CS, low-rank, dictionary learning, …)
3 rd phase Deep learning MR reconstruction (CNN, ADMM-NET,VN-net, Automap, MoDL, U-NET, …)
Ky Kx
Image domain
[1] Pruessmann, Klaas P., et al. Magnetic Resonance in Medicine, 1999. [2] Sodickson, Daniel K., and Warren J.
[3] Griswold, Mark A., et al. Magnetic Resonance in Medicine, 2002. [4] Lustig et al. Magnetic Resonance in Medicine, 2007. [5] Jianhua Luo, Shanshan Wang, et al. Journal of Magnetic Resonance 224 (2012): 82-93 [1] Justin H, et al. Magnetic Resonance in Medicine, 2016 [2] Zhi-Pei Liang, et al. IEEE Transactions in Medical Image, 2003. [3] Zhou, Yihang, et al. IEEE International Symposium on Biomedical Imaging, 2015. [4] Shanshan Wang,et al, IEEE Transactions Medical Imaging, 37(1):251- 261, 2018 [1] Shanshan Wang, et al. . IEEE International Symposium on Biomedical Imaging, 2016. [2] Bo Zhu, Nature, 2018. [3] Florian Knoll, et al. Magnetic Resonance Imaging, 2018. [4] Yang, Guang, et al. IEEE Transactions in Medical Image, 2018. [5] J. Schlemper, et al. IEEE Transactions in Medical Image, 2018.
Partial K-space
Wavelet transform / TV Inverse transform
➢ Layer deconvolution spectral (LDS) analysis method
Jianhua Luo, Shanshan Wang, et al. Journal of Magnetic Resonance 224 (2012): 82-93.
from the image containing truncation artifact
parameters from truncated k- space data
recover Missing k-space data
image through Inverse transform
➢ Main steps Convolution/deconvolution is always a very powerful tool !
Initial ZF image containing truncation artefacts (with truncation frequency c = 64). (b–d) are respectively the images after removing the artefacts in (a) using the Hamming window, the TV and the proposed LDS methods.
MR images. (a and c) Represent respectively a stomach MR image and a brain MR image having truncation artefacts. (b and d) Show the images after removing artefacts in (a) and (c), respectively, using the proposed LDS method.
Jianhua Luo, Shanshan Wang, et al. Journal of Magnetic Resonance 224 (2012): 82-93.
➢ Pros:
➢ Cons:
𝑛 × 1 measurements
𝑐 = 𝐵 𝑌 𝑜 × 1 vector 𝑙 # of non-zeros
Wakin, Michael, et al. "Compressive imaging for video representation and coding." Picture Coding Symposium. Vol. 1. No. 13. 2006. Takhar, Dharmpal, et al. "A new compressive imaging camera architecture using optical-domain compression." Computational Imaging IV. Vol. 6065. International Society for Optics and Photonics, 2006.
Sparse unknown vector
Dynamic imaging time course
Example 1 Example 2 Prior knowledge can be roughly categorized as non-adaptive and adaptive ones. Non-adaptive: Fixed transform, statistical modelling, model fitting, low-rank. Adaptive: Dictionary learning, data-driven tight frame
Shanshan Wang, et al. IEEE Transactions on Image Processing 22.12 (2013): 5214-5225. Shanshan Wang, et al. Signal Processing 93.9 (2013): 2696-2708. Dong Pei, Shanshan Wang*, et al. IEEE Transactions on Image Processing 25.11 (2016): 5035-5049 Qiegen Liu, Shanshan Wang, et al. IEEE Transactions on Image Processing 22.12 (2013): 4652-4663
Shanshan Wang, Dong Liang,et al, IEEE Transactions on Medical Imaging,37(1):251-261, 2018 Shanshan Wang, Dong Liang, et al, BioMed Research International(SCI), 2860643, 2016 Qiegen Liu, Shanshan Wang, et al, IEEE Transactions on Medical Imaging, 2013
adaptive joint sparse coding
frame learning for MR imaging
reconstruction and accelerate the convergence.
smallest reconstruction error. DL-PI SparseSENSE Sparse BLIP CaLM MRI Proposed Label
IFR-CS
Shanshan Wang, Dong Liang, et al. Physics in Medicine and Biology, 2016,61, 3291-3316
𝑣
2 + 𝜇𝑀1 𝑣 𝑀1
𝐽
𝑞𝐽 − 𝑔 2 2 + 𝜈 𝐽 − 𝐽𝑢
The technology has applied for a patent : 201410452350.1 Physics in Medicine and Biology(SCI)
Detected structure Feature descriptor Residual image
Article and 2016 Research Highlight by Physics in Medicine and Biology.
Reference IRM-TV DLMRI Proposed Reconstruction Error Enlarged region
➢ Pros:
➢ Cons:
Linear
CS-MRI
DL-MRI Recon
S Wang, L. Ying, D Liang etal,. IEEE International Symposium on Biomedical Imaging 2016: 514-517.
Offline Training Online Reconstruction Input Reconstruction Training data
Θ
S Wang, L. Ying, D Liang, “Accelerating Magnetic Resonance Imaging via Deep Learning,” IEEE-ISBI 2016: 514-517.
Initialize CS Reconstruction model
Network Prediction Reconstructed MR image
Reference Sampling mask Initialization Network Output Final reconstruction Error
➢ 3T scanner ➢ 32-channel coil ➢ T1-weighted
(spoiled GRE)
➢ TE min full ➢TR = 7.5ms ➢ 256×256, ➢ thickness = 17mm ➢ R = 3
S Wang, L. Ying, D Liang, “Accelerating Magnetic Resonance Imaging via Deep Learning,” IEEE-ISBI 2016: 514-517.
Sun, Jian, Huibin Li, and Zongben Xu. "Deep ADMM-Net for compressive sensing MRI." Advances in Neural Information Processing Systems. 2016.
(MLP, U-NET, GAN, DAGAN, GANCS…)
(Deep ALOHA ,DeepSPIRiT, FCRNN, k-space learning…)
(DIMENSION, AUTOMAP, D5C5, SMS…)
➢ Single channel validation ➢ Multi-coil MRI ➢ Multi-contrast MRI ➢ Dynamic-MR Imaging
Model-based CSMRI learning reconstruction
End-to-end data driven learning reconstruction
Y X 𝑇 = 𝑌 + 𝑗𝑍 X 𝑇 = 𝑌2 + 𝑍2
∠𝑇 = tan−1 𝑍 𝑌
𝑇(𝑙𝑧, 𝑙𝑦) =
𝑧
𝑦
𝐽 𝑦, 𝑧 𝑓
−𝑗2𝜌(𝑙𝑧𝑧+𝑙𝑦𝑦) 𝑂
𝑆𝑓𝑏𝑚(𝐽) Imaginary(𝐽) 𝑁𝑏𝑜𝑗𝑢𝑣𝑒𝑓(𝐽) 𝑄ℎ𝑏𝑡𝑓(𝐽)
Shanshan Wang, Dong Liang, et al, DeepcomplexMRI, Magnetic Resonance Imaging 68, 136-147 Shanshan Wang, Huitao Cheng, et al. ISMRM, 2018;
Shanshan Wang, Hairong Zheng, Dong Liang, et al, DeepcomplexMRI, Magnetic Resonance Imaging 68, 136-147 Shanshan Wang, Dong Liang, et al. ISMRM, Paris, France 2018;
Implemented in tensorflow as real-valued networks representing real and imaginary components with complex-valued operations and initializations.
This figure shows : The comparison of SPIRiT, L1-SPIRiT, VN and the proposed method with real convolution (rc) and with complex convolution half parameter (cc/0.5) and with the same parameter (cc/1) under 1D random (top two lines) and uniform (bottom two lines) sampling at a sampling rate
33%.
PSNR RMSE SSIM CS-MRI 29.8495 0.1369 0.8483 DCCNN 31.8295 0.1090 0.8747 Complex-DCCNN 32.9362 0.0960 0.8759
This figure shows : (a)the reference image and the under-sampled image; reconstructed images by (b) CS-MRI, (c) DCCNN, (d) proposed Complex-DCCNN and their corresponding error maps. The sampling pattern is 2D Poisson disc and the acceleration factor is 5x.
(a) (b) (c) (d)
Shanshan Wang et al. International Society for Magnetic Resonance in Medicine, 2017.
Shanshan Wang et al. International Society for Magnetic Resonance in Medicine, 2017.
Net acceleration factor = 4
Ground truth Deep Learning GRAPPA SPIRiT L1-SPIRiT 7.69 95.80 48.27 40.48
correlation explored with learned filters
arg min
𝒀 1 2
𝐁𝒀 − 𝒁 2
2 + λ𝑡𝓆 𝚾𝑡𝒀 ) + λ𝑑𝑝𝑗𝑚𝑡𝓆(𝚾𝑑𝑝𝑗𝑚𝑡𝒀)
Yanxia Cheng, et al. Shanshan Wang*, MICCAI 2019. Oral (Top 3%)
Code: https://github.com/yanxiachen/ConvDe-AliasingNet. Calibration free No explicit Sensitivity Calculation
PSNR/SSIM
Yanxia Chen Shanshan Wang*, MICCAI 2019.
35/50
k
X
1
9
=
2
15
= + noise
k
X
Trained EDAEP model
Mean
LS solver +1 k
X
Real/imaginary part prior gradient prior gradient
prior gradient
Final prior gradient
f F T
u
1
2
m loops/iteration
MEDAEP multi-contrast MRI reconstruction
Xiangshun Liu, Shanshan Wang*, et al, ISBI 2020
proposed 42.51/0.977 proposed 42.96/0.977 proposed 45.54/0.984 zero-filled PD 24.97/0.658 zero-filled T1 24.48/0.618 zero-filled T2 28.89/0.697 DLMRI 36.90/0.927 DLMRI 37.65/0.913 DLMRI 38.75,0.939
Ground truth PD random 6.7x Ground truth T1 random 6.7x Ground truth T2 random 6.7x
Shanshan Wang, Ziwen Ke, et al. International Society for Magnetic Resonance in Medicine, 2018.
MSE PSNR SSIM Input 0.004213 23.7540 0.7746 D2C2 0.000406 33.9164 0.9708 KI 0.000223 36.5118 0.9818 Reference Undersampled image D2C2 Recon D2C2 Error KI Error KI Recon
Shanshan Wang, Ziwen Ke, et al. International Society for Magnetic Resonance in Medicine, 2018.
Shanshan Wang, Ziwen Ke, Hairong Zheng, Dong Liang, et al. NMR in Biomedicine, DOI:10.1002/nbm.4131 .
Code:https://github.com/Keziwen/DIMENSION
Shanshan Wang, Hairong Zheng, Dong Liang, et al. “ NMR in Biomedicine, DOI:10.1002/nbm.4131 ..
➢ Pros:
➢ Cons:
➢ Handling hard samples ➢ Exploring multi-modality ➢ Designing light network
Deep Neural Network T1C T2W Deep Neural Network with Massive Number
Light Network Deep Neural Network Samples with Large range of stroke lesion scales
Motivation and focus Issue 1 Issue 2 Issue 3
features and context information by the famous encoder-decoder structure
0.5x 0.5x 0.5x 0.5x 2x 2x 2x 2x
Encoder-Decoder T1 GroundTruth
➢ Contribution:
Full use of different levels of features
Address the challenges with big variety of lesion scales
More fine structures captured
Hao Yang& Shanshan Wang*, MICCAI 2019.
c
32 32 256 64 32× 2 64 512 128 128 128× 4 1024 256 256 256× 5 2048 512 512 1 1 32 32 64 64 128 256 128 64 32 1024 256 256 128 256* 4 256 128 64 32 Conv LSTM Conv LSTM Conv LSTM Conv LSTM 256*9 C1 C2 C3 C4 64× 3 256 128 64 32
(a) Cross-Level fusion and Context Inference Network (CLCI-Net)
Conv 1× 1 Stride 1 BN ReLU c Concatenation Conv 3× 3 Stride 1 BN ReLU Up Conv 2× 2 Max Pooling Conv 3× 3 Stride 2 BN ReLU Conv 3× 3 Stride 4 BN ReLU Conv 3× 3 Stride 8 BN ReLU Conv 3× 3 Stride 16 BN ReLU
n
... ...
Cn
(b) Cross-Level fusion (CLF)
Conv 1× 1
(c) Extension of Atrous Spatial Pyramid Pooling (ASPP)
Conv 3× 3 rate 2 Conv 3× 3 rate 4 Conv 3× 3 rate 6 Image Poolin g
...
ASPP Conv 1× 1 Stride 1 Sigmoid
GT Ours Baseline DenseUnet DeepLabv3+ PSPNet FCN-8s T1
a
Lesions with different scales
Hao Yang &Shanshan Wang*, MICCAI 2019.
➢ Contribution:
the network size.
Module (FSM). Capture long-range dependencies
Kehan Qi&Shanshan Wang*, MICCAI 2019.
FSM Encoder Decoder
X-Block Maxpooling 2×2 Upsampling 2×2 and Conv 3×3 Feature Map Convolution Layer Skip and Concate
64 128 256 512 1024 512 256 128 64
compromising ➢ the performances
Kehan Qi&Shanshan Wang*, MICCAI 2019. U-net SegNet PSPNet DenseU Deeplab V3+ ResUNet Ours Groundtruth MR
➢ Limitations of pure 2D or 3D network: Spatial context or Resource occupation
Yongjin Zhou&Shanshan Wang*, IEEE/ACM Transactions on Computational Biology and Bioinformatics 2019.
based networks require too much computing resources, we propose D-UNet.
Up- Sampling 2D Max- Pooling 2d Concatenate BN Conv2d Conv2d BN Conv3d BN Conv3d BN Max- pooling 3d 192×192×1
(a) D-Unet architecture
96×96×2 ×32 96×96×2 ×64 96×96×2 ×64 96×96 ×64 48×48×1 ×64 192×192×4 ×32 48×48×1 ×128 48×48 ×128
(b) Dimension Transform Block
HxWxC
Global pooling FC ReLu FC Sigmoid Mul
1×1×C 1×1×C/r 1×1×C/r 1×1×C Squeeze-and-Excitation Block H×W×C
Global pooling FC ReLu FC Sigmoid Mul
1×1×C 1×1×C/r 1×1×C/r 1×1×C
Conv3d (1×1×1) Squeeze Conv2d (3×3) H×W×D×C H×W×D×1 H×W×D H×W×C H×W×C Add
H×W×C Dimension Transform Block
➢ Contribution:
Extracted 3D information from MRI data and reduced resource consumption.
Loss). Reduce the time required to train the network.
Yongjin Zhou&Shanshan Wang*, IEEE/ACM Transactions on Computational Biology and Bioinformatics 2019.
− − − = − − =
= =
b f
N i N i
p p g if p p g p FL
1 1
), 1 log( ) 1 ( 1 ), log( ) 1 ( ) , (
= =
+ + + − =
N i i i N i i i
g p g p g p DL
1 2 2 1
2 1 ) , (
➢ Traditional Focal Loss(FL) and Dice Loss(DL): ➢ Proposed Enhance Mixing Loss(EML):
Ground Truth Ours UNet Deeplab v3+ PSPNet Original SegNet (1) (2) (3) (4) (5) (6) (7)
Yongjin Zhou&Shanshan Wang*, IEEE/ACM Transactions on Computational Biology and Bioinformatics 2019.
Method DSC Precision Recall SegNet 0.329±0.251 0.385±0.288 0.332±0.265 PSP 0.446±0.263 0.500±0.291 0.470±0.278 Deeplab v3+ 0.453±0.291 0.563±0.325 0.446±0.303 UNet 0.497±0.291 0.551±0.330 0.504±0.304 D-UNet 0.535±0.276 0.633±0.296 0.524±0.291
➢ Multi-modal breast MRI: High-sensitive contrast-enhanced MRI (T1C) + High-specific T2-weighted MRI (T2W) → Accurate diagnosis. ➢ Challenges in multi-modal image segmentation: Effective multi-modal information fusion.
Organs Dense glandular tissues
Cheng Li &Shanshan Wang*, MICCAI 2019.
Overall architecture Master-assistant cross- modal learning framework Cross-modal supervision learning module
Cheng Li &Shanshan Wang*, MICCAI 2019.
performance achieved.
is important.
respective modalities.
Cheng Li &Shanshan Wang*, MICCAI 2019.
Challenges Existing Methods
Limitations
[1] Afshar P, Mohammadi A, Plataniotis K N, et al. From Hand-Crafted to Deep Learning-based Cancer Radiomics: Challenges and Opportunities[J]. arXiv preprint arXiv:1808.07954, 2018.
Two types of tumor
.
Overall workflow
Xiran Jiang, Shanshan Wang*, et al, IEEE/ACM Transactions on Computational Biology and Bioinformatics 2019 Wenqing Hua, Shanshan Wang*, et al, Biomedical signal processing and control, 2020
Siemens 3T United Imaging 3T
Research center Named after the Nobel Laureate
Research bases
➢ Shenzhen Biomedical Imaging Technology Engineering Laboratory ➢ Shenzhen Key Laboratory of Magnetic Resonance Imaging ➢ Guangdong Province Biological Medical Equipment Technology Innovation Platform
Research projects
➢ The National 973 Program ➢ National Key Technology R&D Program ➢ The National Natural Science Foundation of China ➢ The Chinese Academy of Sciences ➢ …
Cooperation company
➢ Cheng Li ➢ Taohui Xiao ➢ Yanxia Chen ➢ Hao Yang ➢ Kehan Qi ➢ Haoyun Liang ➢ Chuyu Rong ➢ Weijian Huang ➢ Yu Gong ➢ Zhenkun Peng ➢ Chen Hu ➢ Rui Yang ➢ Juan Zou
➢ Hairong Zheng(Chinese Academy of Sciences) ➢ David Dagan Feng (The University of Sydney) ➢ Dong Liang (Chinese Academy of Sciences) ➢ Lesilie Ying (University at Buffalo, The State University of New York) ➢ Qiegen Liu (NanChang University) ➢ Yong Xia (Northwestern University of Technology) ➢ Meiyun Wang (Henan Provincial Peoples Hospital) ➢ Zaiyi Liu (Guangdong Provincial Peoples Hospital) ➢ Xiran Jiang (China Medical University) ➢ Yuemin Zhu (University of Lyon 1) ➢ Wanqing Li (University of Wollongong)
➢ National Natural Science Foundation of China ➢ Natural Science Foundation of Guangdong Province ➢ Shenzhen Science and Technology Innovation Committee ➢ …
Code: http://www.escience.cn/people/ShanshanWang/project73837.html