Neuroimage Analysis for Automated Brain Disease Diagnosis
University of North Carolina at Chapel Hill mxliu@med.unc.edu http://mingxia.web.unc.edu/
07-17-2019
Neuroimage Analysis for Automated Brain Disease Diagnosis Mingxia - - PowerPoint PPT Presentation
Neuroimage Analysis for Automated Brain Disease Diagnosis Mingxia Liu University of North Carolina at Chapel Hill mxliu@med.unc.edu http://mingxia.web.unc.edu/ 07-17-2019 Background Healthy Brain vs. Alzheimers Sulcus Sulcus Gyrus
University of North Carolina at Chapel Hill mxliu@med.unc.edu http://mingxia.web.unc.edu/
07-17-2019
Healthy Brain Alzheimer’s
Ventricle Sulcus Gyrus Sulcus Gyrus
Healthy Brain Alzheimer’s
Normal Control Mild Cognitive Impairment (MCI) Alzheimer’s Disease
Stable MCI (sMCI) Progressive MCI (pMCI)
MRI PET CSF CSF Brain Dura
fMRI PET sMRI
Neuroimaging Data
Image Preprocessing Feature Extraction/Selection Classifier Learning Machine Learning & Deep Learning
CSF PET sMRI
sMRI
Training Data
...
Patch Extraction
ata Training Data
...
Landmark Discovery
ata
Test MRI
ata
Landmark Detection
Training Data
...
Training MRIs
Pre-processed MR Images …
… CNN
Deep Multi-channel Convolutional Neural Network (CNN) Disease Classification Clinical Score Prediction * Featured Article of IEEE Journal of Biomedical and Health Informatics, 2018
0.01 p-values
Concatenated
sMRI
FC8-32 FC8-32 FC8-32
…
FC9-32L Conv1- 32@3×3×3 Conv2- 32@3×3×3 Max-pooling Conv3- 64@2×2×2 Conv4- 64@2×2×2 Max-pooling Conv5- 128@2×2×2 Conv6- 128@2×2×2 Max-pooling Conv1- 32@3×3×3 Conv2- 32@3×3×3 Max-pooling Conv3- 64@2×2×2 Conv4- 64@2×2×2 Max-pooling Conv5- 128@2×2×2 Conv6- 128@2×2×2 Max-pooling FC7-128
…
FC7-128 Conv1- 32@3×3×3 Conv2- 32@3×3×3 Max-pooling Conv3- 64@2×2×2 Conv4- 64@2×2×2 Max-pooling Conv5- 128@2×2×2 Conv6- 128@2×2×2 Max-pooling
Patch 1 Patch 2
…
Patch L
FC7-128 FC10-8L FC11-2
Class Label
Soft-max
Global Image-level Representation Local Patch-level Representation
0.431
Overall Accuracy
0.404 0.46 7 0.486 0.487
0.518
0.325
Correlation Coefficient
0.289 0.468 0.492 0.538
0.567
Classification Results for AD vs. sMCI vs. pMCI vs. NC Regression Results for MMSE
15
– Automatically and identify disease-related ROIs in the whole sMR image – Jointly learn multi-scale features and construct a classification model
Input: sMRI
Input: sMRI
1) Location proposals
1) Location proposals
Input: sMRI
2) Patch-level sub- networks (PSN) (shared weights)
PSN PSN
Class_P 64 64 64 128 128 32
PSN
1) Location proposals
… …
PSN PSN PSN PSN PSN PSN PSN Input: sMRI
2) Patch-level sub- networks (PSN) (shared weights)
PSN
Class_P 64 64 64 128 128 32
…
64 Conv_R Class_R
…
PSN PSN PSN PSN PSN PSN PSN PSN
2) Patch-level sub- networks (PSN) (shared weights) 3) Region-level sub-networks 1) Location proposals
Input: sMRI
PSN
Class_P 64 64 64 128 128 32
…
64 Conv_R Class_R
…
64 Conv_R Class_R
…
PSN PSN PSN PSN PSN PSN PSN PSN
2) Patch-level sub- networks (PSN) (shared weights) 3) Region-level sub-networks 1) Location proposals
Input: sMRI
PSN
Class_P 64 64 64 128 128 32
…
64 Conv_S Class_S 64 Conv_R Class_R 64 Conv_R Class_R
… …
Output: Class label
PSN PSN PSN PSN PSN PSN PSN PSN
3) Region-level sub-networks 4) Subject-level sub-network 1) Location proposals
Input: sMRI
2) Patch-level sub- networks (PSN) (shared weights)
Classification (1 × 1 × 1 Conv) 4 × 4 × 4 Conv 1 × 1 × 1 Conv Channel concatenation Region-level Conv 2 × 2 × 2 max pooling 3 × 3 × 3 Conv Spatial concatenation Subject-level Conv Potentially pruned sub-networks Skipped connection Conv: Convolution PSN
Class_P 64 64 64 128 128 32
Network Pruning: To remove less informative regions
1 2 3 1 2 3 1 2 3 1 2 3 4 3 4 3 4 1 2 1 2 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 1 2 1 2 1 2 3 1 2 3 1 2 3 1 2 3
(a) AD Subject #1 (b) AD Subject #2 (c) AD Subject #3 (d) AD Subject #4 (e) AD Subject #5 (f) AD Subject #6
Sagittal View Axial View Coronal View 3D View
Sagittal View Axial View Coronal View 3D View
0.5 0.6 0.7 0.8 0.9 1
Patch Region Subject ACC
0.5 0.6 0.7 0.8 0.9 1
Patch Region Subject AUC
H-FCN before network pruning H-FCN after network pruning Results of AD vs. NC classification obtained by patch-, region-, and subject-level sub- networks in H-FCN without /with network pruning
1) Global subject-level representation is more useful 2) Network pruning promotes the classification performance
CSF PET sMRI
sMRI
30
Multi-Modality Data CSF PET MRI
PET CSF MRI Missing Data
* MICCAI Young Scientist Award Nomination, 2016 * MICCAI Travel Award, 2016
31
Multi-Modality Data CSF PET MRI Multi-View Data Grouping
View 1 View 2 View 3 View 4 View 5 View 6 PET CSF MRI Missing Data
* MICCAI Young Scientist Award Nomination, 2016 * MICCAI Travel Award, 2016
Multi-Modality Data CSF PET MRI Multi-View Data Grouping
View 1 View 2 View 3 View 4 View 5 View 6 Sparse Representation Sparse Representation
Sparse Representation based Hypergraph Construction
Sparse Representation Sparse Representation Sparse Representation Sparse Representation PET CSF MRI Missing Data
View-Aligned Hypergraph Classification
Multi-Modality Data CSF PET MRI Multi-View Data Grouping
View 1 View 2 View 3 View 4 View 5 View 6 Sparse Representation Sparse Representation
Sparse Representation based Hypergraph Construction
Sparse Representation Sparse Representation Sparse Representation Sparse Representation PET CSF MRI Missing Data
Multi-View Label Fusion
Multi-Modality Data CSF PET MRI Multi-View Data Grouping
View 1 View 2 View 3 View 4 View 5 View 6 Sparse Representation Sparse Representation
Sparse Representation based Hypergraph Construction
Sparse Representation Sparse Representation Sparse Representation Sparse Representation
View-Aligned Hypergraph Classification
PET CSF MRI Missing Data
Label Space
𝐲𝑜1
MRI
𝐲𝑜1
PET
𝐲𝑜1
CSF
𝑔
𝑜1 MRI
𝑔
𝑜1 PET
𝑔
𝑜1 CSF
𝑔
𝑜2 MRI
𝑔
𝑜2 PET
𝐲𝑜2
MRI
𝐲𝑜2
PET
Missing CSF MRI PET Subject #2 CSF MRI PET Subject #1
36
𝐆, 𝛃, {𝐗𝑛}𝑛=1
𝑁
𝑁
2
𝑁
𝐺 2
𝑂
𝑁
𝑁
𝑛 𝛻𝑜,𝑜 𝑞
𝑜 𝑛 − 𝑔 𝑜 𝑞 2
𝑁
𝑁
𝑂𝑓
𝑛
𝑗,𝑗 𝑛 = 1, ∀𝑋 𝑗,𝑗 𝑛 ≥ 0.
View-aligned Regularizer Hypergraph Laplacian matrix , and .
AD vs. NC MCI vs. NC pMCI vs. sMCI pMCI vs. NC 0.0 0.1 0.2 0.3 0.4 0.5 0.6
Weights
MRI PET CSF PET+MRI MRI+CSF PET+MRI+CSF
ACC SEN SPE BAC PPV NPV AUC 50 60 70 80 90 100
Results (%)
(c) pMCI vs. sMCI classification
Zero KNN EM SVD VAHL
Pre-processed MR and PET Images Stage 2: Brain Disease Classification Stage 1: PET Image Synthesis based on sMRI Synthesizing Missing PET based on MRI for Brain Disease Diagnosis
Hybrid cycle-consistent generative adversarial network (HGAN) to impute missing PET scans based on sMRI
Synthetic PET Synthetic MRI
32 32
Residual Net Block (RNB)
Real PET
1 128 64 32 16
𝑬𝑸
1 128
64
32 16
𝑬𝑵
1 16 32 32 16 8
𝑯𝑵
⋯
RNB RNB
8 16 32
𝑯𝑸
⋯
RNB RNB
32 16 1 3×3×3 Convolution 7×7×7 Convolution 3×3×3 Deconvolution 4×4×4 Convolution Addition
Real MRI
6 6 1/0 1/0
𝐘𝑁 𝐘𝑄 𝐻1 𝐘𝑁 𝐻2 𝐘𝑄
PET (RID: 5016)
(a) GAN (b) CGAN (c) VGAN (d) HGAN (Ours) (e) Ground Truth
PET (RID: 4352)
MRI (RID: 5016)
(a) GAN (b) CGAN (c) VGAN (d) HGAN (Ours) (e) Ground Truth
MRI (RID: 4352)
PET (Real + Synthetic) MRI (Real)
24 24 24 24 24 24 24 24 24 24 24 24 … 32 32 32 64 64 64
Sub-network 1
DCM DCM DCM 16 16 16
Sub-network 2
32 32 32 64 64 64 DCM DCM DCM 16 16 16 Concatenation 32 8K …
Clinical scores at four time-points
8 32 8 128 32 BL M06 M24 M12
Fully-connected Down-sampling Copy 3×3×3 Convolution 2×2×2 Max-pooling Channel concatenation
Landmark-based Deep Network for Brain Disease Classification using MRI and PET (Real+Synthetic)
PET (Real + Synthetic)
MCI conversion prediction with complete MRI and complete (after imputation) PET
1 128 64 32 16
𝐸
3×3×3 Conv 7×7×7 Conv 3×3×3 Deconv 4×4×4 Conv Addition Difference
32 32 RNB 1 16 32 32 16 8
𝐻
RNB RNB ⋯ 6
𝐻𝑁 𝐻𝑄 𝐸𝑄 𝐸𝑁
𝔐 𝔐
𝐺
𝑄
𝐺
𝑁
𝔐c 𝔐c FG AN
Synthetic MRI Real PET Synthetic PET Real MRI Synthetic PET Real PET
1 2 K
1 2 K
⋮
1/0 𝑚2 64 64 32 16 64 64 64 32 16
𝔐c
64
Feature-consistent Component 𝐺
𝑄
Cosine Kernel
Feature-consistent Component 𝐺
𝑄
30 50 70 90
AUC ACC SEN SEP ROI LMF LDMIL DSNN (Ours)
MCI conversion prediction with complete MRI and complete (after imputation) PET
Scanners Scanning Parameters Population Noise Level Image Contrast
*Best Poster Award, MICS, 2019
50
P2 P1 PS P New Representation
Latent Representation Space
Linear Representation
Site T Site S Site 1 Site 2
*Best Poster Award, MICS, 2019
51
𝐿
∗ + 𝛽 𝐅𝑇𝑗 1 + 𝛾 𝐅𝑄𝑗 1
*Best Poster Award, MICS, 2019
52
Performance of different methods in the task of Autism vs. NC classification, with NYU as the target domain and the other four sites as the source domains
*Best Poster Award, MICS, 2019
CSF PET sMRI
sMRI
mxliu@med.unc.edu http://mingxia.web.unc.edu/