In this study, a Generative Adversarial Network is used for detection of abnormalities in MR brain images such as lesions, artifacts etc. Given a query image, a generative model that is trained to create normal appearing brain images is used to find a best match. Since abnormalities cannot be reproduced accurately by the generative model, pathologies and artifacts become apparent.
Neural networks have shown great potential for automatic detection of abnormal structures in medical images - both pathologies1 and artifacts2. The main requirement for such convincing performance, however, is the availability of large amounts of high-quality annotations for each abnormality to be detected. For many pathologies (in particular rare diseases) and artifacts, such annotations are usually not available.
Here we demonstrate a general abnormality detection technique (e.g. lesions, artifacts etc.) based on a Generative Adversarial Network3 (GAN) that does not require annotations of abnormal structures. It relies on an unsupervised learning approach with the aim to train a generative model to cover the distribution of all normal images, thereby allowing for the detection of all possible abnormalities.
A large dataset of normal MP-RAGE brain MR images was generated based on the ADNI dataset4. Images with motion artifacts were discarded. All remaining images were normalized on a volume level and cropped to a matrix size of 160 x 256. Only the central 30 slices were selected to reduce the variability of the data to be modeled, resulting in training and testing datasets of 43,320 and 750 images, respectively. Data augmentation was realized using random translation, rotation, horizontal flipping, and deformation. A symmetric GAN was then trained to learn the manifold of these MR images.
The Generator consisted of 6 transposed convolutional layers (kernel size k=4, stride s=2), each followed by batch normalization and a Recitified Linear Unit (ReLU). A vector of 200 random numbers provided the input to the Generator. The Discriminator consisted for 6 convolutional layers (k=4, s=2), where batch normalization and leaky ReLUs were used. The GAN was trained for 90 epochs using the Adam optimizer. An additional regularization term was implemented that penalized large gradients of the discriminator on real data, i.e. on the true data distribution5.
During evaluation, all parameters of the networks were kept fixed, and backpropagation was used to determine the optimal input vector z that produced the GAN-generated image most similar to the query image. The loss function for this search was given by
$$L(z)=(1-\alpha)\sum_{x,y}|I_q(x,y)-I_g(x,y;z)|^2 + \alpha|D(I_g(x,y;z))-1|^2$$
where $$$\alpha$$$ is a weighting factor (chosen as 0.02), $$$I_q$$$ and $$$I_g$$$ are the query and generated images, respectively, and $$$D$$$ is the Discriminator. While the first term enforces similar image content in query and generated image, the second term ensures that generated images lie on the manifold of normal images.
To simulate abnormal tissue
structures, local hyper- and hypointensities were added to a subset of the test
images. To simulate typical MR artifacts, fold-over artifacts due to
incorrectly placed FOV as well as susceptibility-induced signal voids were
added to another subset of the test images.
1. Pereira S, et al. Brain tumor segmentation using convolutional neural networks in MRI images. IEEE Trans. Med. Imaging 2016;35:1240–51.
2. Demesmaeker R, et al. Deep Learning for Automated Medical Image Quality Assessment: Proof of Concept in Whole-Heart Magnetic Resonance Imaging. Proc ISMRM 2018; no. 4784.
3. Goodfellow I, et al. Generative adversarial nets. Advances in neural information processing systems 2014:2672-2680.
4. Jack Jr CR, et al. The Alzheimer's disease neuroimaging initiative (ADNI): MRI methods. J Magn Reson Imaging 2008;27(4):685-91.
5. Mescheder L. et al. Which Training Methods for GANs do actually Converge? International Conference on Machine Learning 2018.