1 Introduction
Image registration is the process of aligning two or more images to achieve pointwise spatial correspondence. This is a fundamental step for many medical image analysis tasks and has been an active field of research for decades. Since recently, deep learning based approaches have been successfully employed for image registration [2, 4, 6, 7, 13, 16]. They have shown promising results in a wide range of application. However, capturing large motion and deformation with deep learning based registration is still an open challenge. In common iterative image registration approaches, this is typically addressed with a multilevel coarsetofine registration strategy [1, 11, 15]. Starting on a coarse grid with smoothed and downsampled versions of the input images a deformation field is computed which is subsequently prolongated on the next finer level as a initial guess. Hereby, a coarse level alignment is obtained first that typically captures the large motion components and which is later improved on finer levels for the alignment of more local details. Most of the recently presented deep learning based approaches also make use of a multilevel strategy as they are based on the UNet architecture [2, 6, 7, 13]. Thereby, the first half of the ”U” is used to generate features on different scales starting at the highest resolution and reducing the resolution through pooling operations. In this procedure, however, only feature maps on different levels are calculated but neither different image resolutions are used nor deformation fields are computed. Only a few approaches implement a multiresolution or hierarchical strategy in the sense of multilevel strategies associated with conventional methods. In [8] the authors proposed an architecture which is divided into a global and a local network, which are optimized together. In [4] a multilevel strategy is incorporated into the training of a Unet. Here, a CNN is grown and trained progressively levelbylevel. In [16] a patch based approach is presented, where multiple CNNs (ConvNets) are combined additive into a larger architecture for performing coarsetofine image registration of patches. The results from the patches are then combined into a deformation field warping the whole image. In this work, we address this challenge and present a multilevel strategy for deep learning based image registration to advance stateoftheart approaches. The contribution of this paper includes:

We present deep learning based multilevel registration that is able to compensate and handle large deformations by computing deformation fields on different scales and functionally compose them.

Our method is a theoretically sound and a direct transference of coarsetofine registration from conventional, iterative registration schemes to the deep learning based methods.

We do not rely on patches. We take the whole image information into account and always consider the full field of view on all levels.

A robust and fast registration method for the complex task of inhaletoexhale registration validated on a large dataset of 270 thoracic CT scan pairs of the multicenter COPDGene study and on the publicly available DIRLAB dataset [3].
2 Method
Our deep learning based framework for deformable image registration consists of two main building blocks. The first one is the specific design of the convolutional neural network and the loss function. In general, several architectures together with different distance measures, regularizer and penalty terms can be used. However, we focus on a UNet based architecture, combined with a loss function that has shown good results for the task of pulmonary registration
[14]. The second main building block is the embedding into a multilevel approach from coarse to fine. In the following, we give a brief outline of the variational setup, then we describe our particular architecture and loss function and, finally, we present its embedding into a multilevel approach.Variational Registration Approach:
Following [10], let denote the fixed image and moving image, respectively, and let be a domain modeling the field of view of . We aim to compute a deformation that aligns the fixed image and the moving image on the field of view such that and are similar for . The deformation is defined as a mimimizer of a suitable cost function that typically takes the form
(1) 
with socalled distance measure that quantifies the similarity of fixed image and deformed moving image and socalled a regularizer that forces smoothness of the deformation typically by penalizing of spatial derivatives. Typical examples for the distance measure are, e.g., the squared norm of the difference image (SSD), cross correlation (CC) or mutual information (MI). In our experiments, we follow the approach of [14] using the edge based normalized gradient fields distance measure (NGF) and second order curvature regularization.
Deep Learning based Image Registration
In contrast to conventional registration [10], we do not employ iterative optimization during inference of new unseen images but use a convolutional neural network (CNN) that takes images and as input and yields the deformation as output. Thus, in the context of CNNs we can consider
as a function of a trainable CNN model parameter vector
and input images , , i.e.. In an unsupervised learning approach, we set up a loss function
that depends on , and , and then is learned by training, i.e., minimizing the expected value of among a set of representative input images w.r.t. . A natural choice would . However, in our particular application, we have additional information available during training and we perform a weakly supervised approach. To this end, we define our loss function as suggested in [7](2) 
with binary segmentation masks and of the fixed and warped moving image, respectively. Note that these segmentations are only used to evaluate the loss function for training and their are not used as network input.
Single Level Architecture
Our CNN
is based on a UNet which takes the concatenated 3D moving and fixed image as input and predicts a 3D dense displacement field. The network consists of three resolution levels starting with 16 filters in the first layer, which are doubled after each downsampling step. We apply 3D convolutions in both encoder and decoder stage with a kernel size of 3 followed by a batch normalization and a ReLU layer. For downsampling the feature maps during the encoder path, an
average pooling operation with a stride of 2 is used. Transposed convolutions upsample and halve the feature maps in the decoder path. At the final layer, a
convolution is used to map each 16 component feature vector to a three dimensional displacement vector.Multilevel Deep Learning based Registration
Multilevel continuation and scale space techniques have been proven very efficient in conventional variational registration approaches to avoid local minima, to reduce topological changes or foldings and to speed up runtimes [1, 9, 11, 15]. However, beside carrying over these properties, our major motivation here is, to overcome the limitation of deep learning based registration to small and local deformations.
We follow the ideas of standard multilevel registration and compute coarse grid solutions that are prolongated and refined on the subsequent finer level. To this end, first we create image pyramids for with coarsest level . We start on finest level and subsequently halve image size and resolution from level to level. Registration starts on coarsest level and we compute deformation from images and as network input. On all finer levels , we incorporate the deformations from all preceding coarse levels as initial guess. Therefore, we combine them by functional composition and warp the moving image at current level. Let denote the cellcentered image grid on level , we compute the warped moving with
and use it together with fixed image as network input, yielding the deformation field on the current level. The final output deformation is then given by composition of the whole sequence of coarsetofine solutions, i.e.,
. To evaluate deformations and images at nongrid grid points, we use trilinear interpolation. Our scheme is summarized in Algorithm
1.In our experiments we use in particular a three level scheme (
). and we create image pyramid with three reduced resolution images generated from the original 3D images by applying a lowpass filter with a stride of two, four and eight. During training, the three networks are learned progressively. First, the network on the coarsest level is trained for a fixed amount of epochs. Afterwards, the parameters of the middle network are learned while the coarsest network stays fixed and is only used to produce the initial deformation field. The same procedure is repeated on the finest level. The same architecture is used on all levels. The convolution parameters on the coarsest level are initialized with Xavier uniform
[5]. Whereas, all other networks are using the learned parameters of the previous network as initialization. Note that the receptive field in voxel is the same for all used networks, however, due to to the decreased resolution on the coarse levels, the receptive field in mm is much higher.3 Experiments and Results
We demonstrate our deep learning based registration method by registration of inhaletoexhale lung CT scans. We use data from 500 patients for training and a disjoint set of 50 patients for validation from the COPDGene study, a large multicenter clinical trial with over 10,000 subjects with chronic obstructive pulmonary disease (COPD) [12]. The dataset was acquired across 21 imaging centers using a variety of scanner makes and models. Each patient had received two breathhold 3D CT scans, one on full inspiration (200mAs) and one at the end of normal expiration (50mAs). For all scans segmentations of the lobes are available, which were computed automatically and manually corrected and verified by trained human analysts. The original images have sizes in the range of voxels. Due to memory and time limitations, we create lowresolution images by resampling to a fixed size of voxels. The lowresolution images are then used during training for the computation of the deformation field and for evaluating of the loss function. Note that, our method is generally not limited to any fixed input size. Although we use images with voxels, the computed deformation field are defined on full fieldofview and can be evaluated on grids with arbitrary resolution by using trilinear interpolation. Consequently, we use original fullresolution images for the evaluation of our method.
Multilevel vs. Single Level
First, we evaluate our multilevel approach on a disjoint subset of 270 patients from the COPDGene study. We compare our proposed method against a single level approach with only one UNet using the images on finest level as inputs. We train both approaches for 75 epochs with the same hyperparameters. For the multilevel approach, the epochs are split equally at each level. We also evaluate the effect on how the network parameter are initialized. Therefore, we compare a Xavier initialization for all convolution parameters of all three networks against our proposed progressive learning strategy. Therefore, only the convolution parameters on the coarsest level are initialized with Xavier initialization and the training of subsequent network is started with the learned parameters form the network of the previous level.
We evaluate our method by measuring the overlap of the lobe masks. The underlying assumption is, that if a deformation yields accurate spatial correspondences, then lobe segmentations of the fixed and the warped lobe segmentation of the moving image should overlap well. Fig. 1 shows the Dice scores for each label and the average over all labels as a boxplot. Our proposed multilevel approach increase the Dice Score from to . In contrast, the single level method archive a Dice Score of . Furthermore, the multilevel approach produced less foldings ( to 2.1% ). Figure 2 shows representative qualitative results for of two scan pairs before registration and after our single level and multilevel registration. In both cases the respiratory motion was successfully recovered. Although the single level registration produces reasonable Dice scores, it does not well align the inner structures. This is also reflected by the landmark errors in the following section. Comparing the results of the pretrained initialization to the random initialization, an improvement of about in terms of the Dice Score could be reached.
(a) before  (b) single level  (c) multilevel  (d) before  (e) single level  (f) multilevel 
Comparison with stateoftheart
Additionally, we evaluate our method and compare it to others on the public available DIRLAB dataset [3]
. It is a collection of ten inspirationexpiration cases with 300 expertannotated landmarks in the lung. The landmarks are used for evaluating our deformable registration method. The mean (and standard deviation) for all ten scans for the deep learning based multiresolution approaches of Eppenhof
[4] and de Vos (DLIR) [16], the single VIRNET and our proposed method are listed in table 1. The overall average landmark error is with a standard deviation of. In contrast to the other methods, our mlVIRNET is more robust against outliers and can better handle large initial landmark distances without training on this specific dataset.
Scan  Initial  Eppenhof [4]  DLIR [16]  single VIRNET  mlVIRNET 

Case 1  3.89(2.78)  2.18(1.05)  1.27(1.16)  1.73(0.83)  1.33(0.73) 
Case 2  4.34(3.90)  2.06(0.96)  1.20(1.12)  2.38(1.11)  1.33(0.69) 
Case 3  6.94(4.05)  2.11(1.04)  1.48(1.26)  3.01(1.86)  1.48(0.94) 
Case 4  9.83(4.85)  3.13(1.60)  2.09(1.93)  4.28(2.37)  1.85(1.37) 
Case 5  7.48(5.50)  2.92(1.70)  1.95(2.10)  3.17(2.2)  1.84(1.39) 
Case 6  10.89(6.96)  4.20(2.00)  5.16(7.09)  4.85(3.04)  3.57(2.15) 
Case 7  11.03(7.42)  4.12(2.97)  3.05(3.01)  3.67(1.82)  2.61(1.63) 
Case 8  14.99(9.00)  9.43(6.28)  6.48(5.37)  5.75(3.93)  2.62(1.52) 
Case 9  7.92(3.97)  3.82(1.69)  2.10(1.66)  4.90(2.25)  2.70(1.46) 
Case 10  7.30(6.34)  2.87(1.96)  2.09(2.24)  3.49(2.21)  2.63(1.93) 
Total  8.46(6.58)  3.68(3.32)  2.64(4.32)  3.72(2.45)  2.19(1.62) 
4 Discussion and Conclusion
We presented an endtoend multilevel framework for deep learning based image registration which is able to compensate and handle large deformations by computing deformation fields on different scales. Our method takes the whole image information into account and predicts a dense 3D deformation field. We validated our framework on the challenging task of large motion inhaletoexhale registration using large image data of the multicenter COPDGene study. We have shown that our proposed method archives better results than the comparable single level variant. In particular with regard to the alignment of inner lung structures and the presence of foldings. Only less than voxel positions of the images showed a folding. Additionally, we demonstrated that using the network parameter of the previous level as initialization, yields to better registration results. Moreover, we demonstrated the transferability of our approach to new datasets by evaluating our learned method on the publicly available DIRLAB dataset and showing a lower landmark error than other deep learning based registration methods.
Acknowledgements:
We gratefully acknowledge the COPDGene Study for providing the data used. COPDGene is funded by Award
Number U01 HL089897 and Award Number U01 HL089856 from the National Heart, Lung, and Blood Institute. The COPDGene project is also supported by the COPD Foundation through contributions made to an Industry Advisory Board comprised of AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Novartis, Pfizer, Siemens and Sunovion.
References

[1]
Bajcsy, R., Kovačič, S.: Multiresolution elastic matching. Computer vision, graphics, and image processing
46(1), 1–21 (1989)  [2] Balakrishnan, G., Zhao, A., Sabuncu, M.R., Guttag, J., Dalca, A.V.: Voxelmorph: a learning framework for deformable medical image registration. IEEE TMI (2019)
 [3] Castillo, R., Castillo, E., Fuentes, D., Ahmad, M., Wood, A.M., Ludwig, M.S., Guerrero, T.: A reference dataset for deformable image registration spatial accuracy evaluation using the copdgene study archive. Physics in Medicine & Biology 58(9), 2861 (2013)
 [4] Eppenhof, K.A., Lafarge, M.W., Pluim, J.P.: Progressively growing convolutional networks for endtoend deformable image registration. In: Medical Imaging 2019: Image Processing. vol. 10949, p. 109491C. International Society for Optics and Photonics (2019)

[5]
Glorot, X., Bengio, Y.: Understanding the difficulty of training deep feedforward neural networks. In: Proceedings of the thirteenth international conference on artificial intelligence and statistics. pp. 249–256 (2010)
 [6] Hering, A., Heldmann, S.: Unsupervised learning for large motion thoracic ct followup registration. In: SPIE Medical Imaging: Image Processing. vol. 10949, p. 109491B (2019)
 [7] Hering, A., Kuckerts, S., Heldmann, S., Heinrich, M.P.: Enhancing labeldriven deep deformable image registration with local distance metrics for stateoftheart cardiac motion tracking. In: BVM 2019. pp. 309–314. Springer (2019)

[8]
Hu, Y., Modat, M., Gibson, E., Ghavami, N., Bonmati, E., Moore, C.M., Emberton, M., Noble, J.A., Barratt, D.C., Vercauteren, T.: Labeldriven weaklysupervised learning for multimodal deformable image registration. In: Proc. of ISBI 2018. pp. 1070–1074. IEEE (2018)
 [9] Kabus, S., Lorenz, C.: Fast elastic image registration. Medical Image Analysis for the Clinic: A Grand Challenge pp. 81–89 (2010)
 [10] Modersitzki, J.: FAIR: Flexible Algorithms for Image Registration. SIAM (2009)
 [11] Modersitzki, J., Haber, E.: COFIR: Coarse and Fine Image Registration, chap. 14, pp. 277–288. Computational Science & Engineering: RealTime PDEConstrained Optimization, SIAM (2007)
 [12] Regan, E.A., Hokanson, J.E., Murphy, J.R., Make, B., Lynch, D.A., Beaty, T.H., CurranEverett, D., Silverman, E.K., Crapo, J.D.: Genetic epidemiology of copd (copdgene) study design. COPD: Journal of Chronic Obstructive Pulmonary Disease 7(1), 32–43 (2011)
 [13] Rohé, M.M., Datar, M., Heimann, T., Sermesant, M., Pennec, X.: Svfnet: Learning deformable image registration using shape matching. In: MICCAI 2017. pp. 266–274. Springer (2017)

[14]
Rühaak, J., Polzin, T., Heldmann, S., Simpson, I.J., Handels, H., Modersitzki, J., Heinrich, M.P.: Estimation of large motion in lung ct by integrating regularized keypoint correspondences into dense deformable registration. IEEE TMI
36(8), 1746–1757 (2017)  [15] Schnabel, J.A., Rueckert, D., Quist, M., Blackall, J.M., CastellanoSmith, A.D., Hartkens, T., Penney, G.P., Hall, W.A., Liu, H., Truwit, C.L., et al.: A generic framework for nonrigid registration based on nonuniform multilevel freeform deformations. In: MICCAI 2001. pp. 573–581. Springer (2001)
 [16] de Vos, B.D., Berendsen, F.F., Viergever, M.A., Sokooti, H., Staring, M., Išgum, I.: A deep learning framework for unsupervised affine and deformable image registration. Medical image analysis 52, 128–143 (2019)