Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The First Application of a Gd3Al2Ga3O12:Ce Single-Crystal Scintillator to Neutron Radiography
Next Article in Special Issue
Recovering the Magnetic Image of Mars from Satellite Observations
Previous Article in Journal
Quantitative Evaluation of Soil Structure and Strain in Three Dimensions under Shear Using X-ray Computed Tomography Image Analysis
Previous Article in Special Issue
On a Variational and Convex Model of the Blake–Zisserman Type for Segmentation of Low-Contrast and Piecewise Smooth Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimization-Based Meta-Learning Model for MRI Reconstruction with Diverse Dataset

1
Department of Mathematics, University of Florida, Gainesville, FL 32611, USA
2
Department of Mathematics and Statistics, Georgia State University, Atlanta, GA 30303, USA
*
Author to whom correspondence should be addressed.
Author list in alphabetical order.
J. Imaging 2021, 7(11), 231; https://doi.org/10.3390/jimaging7110231
Submission received: 23 September 2021 / Revised: 26 October 2021 / Accepted: 28 October 2021 / Published: 31 October 2021
(This article belongs to the Special Issue Inverse Problems and Imaging)

Abstract

:
This work aims at developing a generalizable Magnetic Resonance Imaging (MRI) reconstruction method in the meta-learning framework. Specifically, we develop a deep reconstruction network induced by a learnable optimization algorithm (LOA) to solve the nonconvex nonsmooth variational model of MRI image reconstruction. In this model, the nonconvex nonsmooth regularization term is parameterized as a structured deep network where the network parameters can be learned from data. We partition these network parameters into two parts: a task-invariant part for the common feature encoder component of the regularization, and a task-specific part to account for the variations in the heterogeneous training and testing data. We train the regularization parameters in a bilevel optimization framework which significantly improves the robustness of the training process and the generalization ability of the network. We conduct a series of numerical experiments using heterogeneous MRI data sets with various undersampling patterns, ratios, and acquisition settings. The experimental results show that our network yields greatly improved reconstruction quality over existing methods and can generalize well to new reconstruction problems whose undersampling patterns/trajectories are not present during training.

1. Introduction

Deep learning methods have demonstrated promising performance in a variety of image reconstruction problems. However, deep learning models are often trained for specific tasks and require the training samples to follow the corresponding distribution. In particular, the source-domain/training samples and target-domain/testing samples need to be drawn from the same distribution [1,2,3,4]. In practice, these data sets are often collected at different sources and exhibit substantial heterogeneity, and thus the samples may follow related but different distributions in real-world applications [5,6]. Therefore, the robust and efficient training of deep neural networks using such data sets is theoretically important and practically relevant in the application of deep learning-based methods.
Meta-learning provides a unique paradigm to achieve robust and efficient neural network training [1,4,7,8,9,10]. Meta-learning is known as learning-to-learn and aims to quickly learn unseen tasks from the experience of learning episodes that cover the distribution of relevant tasks. In a multiple-task scenario, given a family of tasks, meta-learning has been proven to be a useful tool for extracting task-agnostic knowledge and improving the learning performance of new tasks from that family [11,12]. We leverage this feature of meta-learning for network training where the MRI training data are acquired by using different under-sampling patterns (e.g., Cartesian mask, Radial mask, Poisson mask), under-sampling ratios, and different settings of the scanning parameters, which result in different levels of contrast (e.g., T1-weighted, T2-weighted, proton-density (PD), and Flair). These data are vastly heterogeneous and can be considered as being from various tasks. Thus, our goal is to develop a robust and generalizable image reconstruction method in the meta-learning framework to leverage such large-scale heterogeneous data for MRI reconstruction.
Our approach can be outlined as follows. First, we introduce a variational model rendering a nonconvex nonsmooth optimization problem for image reconstruction. In our variational model, the regularization term is parameterized as a structured deep neural network where the network parameters can be learned during the training process. We then propose a learnable optimization algorithm (LOA) with rigorous convergence guarantees to solve this optimization problem. Then, we construct a deep reconstruction network by following this LOA exactly; namely, each phase of this LOA-induced network is exactly one iteration of the LOA. This approach is inspired by [13], but the LOA developed in the present work is computationally more efficient than that in [13]: the safeguard variable in this work is updated only if necessary, which can significantly reduce computational cost while retaining the convergence guarantee.
Second, to improve network robustness and mitigate the overfitting issue, we explicitly partition the network parameters of the regularization into two parts: a task-invariant part and a task-specific part. The former extracts common prior information of images from different tasks in the training data and learns the task-invariant component of the regularization. The latter, on the other hand, is obtained from another network which exploits proper task-specific components (also called meta-knowledge) of the regularization for different tasks. The hyperparameters of this network are also learned during training. Furthermore, we split the available data into two sets: the training set and the validation set. Then, we introduce a bilevel optimization model for learning network parameters. Specifically, the lower-level problem (also known as inner problem) finds the task-invariant part of the regularization with the fixed task-specific part on the training dataset, whereas the upper-level (outer) problem seeks for the optimal task-specific part of the regularization parameter using the validation dataset. This approach greatly increases the robustness of the learned regularization, meaning that the trained LOA-induced deep reconstruction network can generalize well to unseen tasks.
As demonstrated by our numerical experiments in Section 5, our proposed framework yields much improved image qualities using diverse data sets of various undersampling trajectories and ratios for MRI image reconstruction. The reason is that effective regularization can integrate common features and prior information from a variety of training samples from diverse data sets, but they need to be properly weighed against the data fidelity term obtained in specific tasks (i.e., undersampling trajectory and ratios). Our contributions can be summarized as follows:
  • An LOA inspired network architecture—our network architecture exactly follows a proposed LOA with guaranteed convergence. Thus, the network is more interpretable, parameter-efficient, and stable than existing unrolling networks.
  • Adaptive design of regularization—our adaptive regularizer consists of a task-invariant part and a task-specific part, both of which can be appropriately trained from data.
  • Improved network robustness and generalization ability—we improve the robustness of the network parameter training process by posing it as a bilevel optimization using training data in the lower-level and validation data in the upper-level. This approach also improves the generalization ability of the trained network so that it can be quickly adapted to image reconstruction with new unseen sampling trajectories and produces high-quality reconstructions.
The remainder of the paper is organized as follows. In Section 2, we discuss related work for both optimization-based meta-learning and deep unrolled networks for MRI reconstructions. We propose our meta-learning model and the neural network in Section 3 and describe the implementation details in Section 4. Section 5 provides the numerical results of the proposed method. Section 6 concludes the paper.

2. Related Work

In recent years, meta-learning methods have demonstrated promising results in various fields with different techniques [12]. Meta-learning techniques can be categorized into three groups [14,15,16]: metric-based methods [17,18,19], model-based methods [20,21,22,23], and optimization-based methods [8,24,25]. Optimization-based methods are often cast as a bilevel optimization problem and exhibit relatively better generalizability for wider task distributions. We mainly focus on optimization-based meta-learning in this paper. For more comprehensive literature reviews and developments of meta-learning, we refer the readers to the recent surveys [12,16].
Optimization-based meta-learning methods have been widely used in a variety of deep learning applications [8,24,25,26,27,28,29,30,31]. The network training problems in these meta-learning methods are often cast as the bilevel optimization of a leader variable and a follower variable. The constraint of the bilevel optimization is that the follower variable is optimal for the lower-level problem for each fixed leader variable, and the ultimate goal of bilevel optimization is to find the optimal leader variable (often, the corresponding optimal follower variable as well) that minimizes the upper-level objective function under the constraint. The lower-level problem is approximated by one or a few gradient descent steps in many existing optimization-based meta learning applications, such as Model-Agnostic Meta-Learning (MAML) [8], and a large number of followup works of MAML proposed to improve generalization using similar strategy [9,15,27,29,30,32,33,34]. Deep bilevel learning [35] seeks to obtain better generalization than when trained on one task and generalize well to another task. The model is used to optimize a regularized loss function to find network parameters from the training set and identify hyperparameters so that the network performs well on the validation dataset.
When the unseen tasks lie in inconsistent domains with the meta-training tasks, as revealed in [36], the generalization behavior of the meta-learner will be compromised. This phenomenon partially arises from the meta-overfitting on the already seen meta-training tasks, which is identified as a memorization problem in [34]. A meta-regularizer forked with information theory is proposed in [34] to handle the memorization problem by regulating the information dependency during the task adaption. MetaReg [4] decouples the entire network into the feature network and task network, where the meta-regularization term is only applied to the task network. They first update the parameters of the task network with a meta-train set to obtain the domain-aligned task network and then update the parameters of the meta-regularization term on the meta-test set to learn the cross-domain generalization. In contrast to MetaReg, Feature-Critic Networks [37] exploit the meta-regularization term to pursue a domain-invariant feature extraction network. The meta-regularization is designed as a feature-critic network that takes the extracted feature as an input. The parameters of the feature extraction network are updated by minimizing the new meta-regularized loss. The auxiliary parameters in the feature-critic network are learned by maximizing the performance gain over the non-meta case. To effectively evaluate the performance of the meta-learner, several new benchmarks [38,39,40] were developed under more realistic settings that operate well on diverse visual domains. As mentioned in [39], the generalization to unseen tasks within multimodal or heterogeneous datasets remains a challenge to the existing meta-learning methods.
The aforementioned methods pursue domain generalization for the classification networks that learned a regularization function to learn cross-domain generalization. Our proposed method was developed to solve the inverse problem, and we construct an adaptive regularization that not only learns the universal parameters among tasks but also the task-aware parameters. The designated adaptive regularizer assists the generalization ability of the deep model so that the well-trained model can perform well on heterogeneous datasets of both seen and unseen tasks.

3. Proposed Method

3.1. Preliminaries

We first provide the background of compressed sensing MRI (CS-MRI), the image reconstruction problem, and the learned optimization algorithm to solve the image reconstruction problem. CS-MRI accelerates MRI data acquisition by under-sampling the k-space (Fourier space) measurements. The under-sampled k-space measurement are related to the image by the following formula [41]:
y = P F x + n ,
where y C p represents the measurements in k-space with a total of p sampled data points, x C N × 1 is the MR image to be reconstructed with N pixels, F C N × N is the 2D discrete Fourier transform (DFT) matrix, and P R p × N ( p < N ) is the binary matrix representing the sampling trajectory in k-space. n is the acquisition noise in k-space.
Solving x from (noisy) under-sampled data y according to (1) is an ill-posed problem. An effective strategy to elevate the ill-posedness issue is to incorporate prior information to the reconstruction. The variational method is one of the most effective ways to achieve this. The general framework of this method is to minimize an objective function that consists of a data fidelity term and a regularization term as follows:
x ¯ = arg min x 1 2 P F x y 2 + R ( x ) ,
where the first term is data fidelity, which ensures consistency between the reconstructed x and the measured data y , and the second term R ( x ) is the regularization term, which introduces prior knowledge to the image to be reconstructed. In traditional variational methods, R ( x ) is a hand-crafted function such as Total Variation (TV) [42]. The advances of the optimization techniques allowed more effective algorithms to solve the variational models with theoretical justifications. However, hand-crafted regularizers may be too simple to capture subtle details and satisfy clinic diagnostic quality.
In recent years, we have witnessed the tremendous success of deep learning in solving a variety of inverse problems, but the interpretation and generalization of these deep-learning-based methods still remain the main concerns. As an improvement over generic black-box-type deep neural networks (DNNs), several classes of learnable optimization algorithms (LOAs) inspired neural networks, known as unrolling networks, which unfold iterative algorithms to multi-phase networks and have demonstrated promising solution accuracy and efficiency empirically [43,44,45,46,47,48,49,50]. However, many of them are only specious imitations of the iterative algorithms and hence lack the backbone of the variational model and any convergence guarantee.
In light of the substantial success of deep learning and the massive amount of training data now available, we can parameterize the regularization term as a deep convolutional neural network (CNN) that learns from training samples. LOA-induced reconstruction methods have been successfully applied to CS-MRI to solve inverse problems with a learnable regularizer:
arg min x 1 2 P F x y 2 + R ( x ; Θ ) .
where R ( x ; Θ ) is the regularization parameterized as a deep network with parameter Θ . Depending on the specific parametric form of R ( x ; Θ ) and the optimization scheme used for unrolling, several unrolling networks have been proposed in recent years. For example, the variational network (VN) [51] was introduced to unroll the gradient descent algorithm and parametrize the regularization as a combination of linear filters and nonlinear CNNs. MoDL [52] proposed a weight sharing strategy in a recursive network to learn the regularization parameters by unrolling the conjugate gradient method. ADMM-Net [53] mimics the celebrated alternating direction method of multipliers; the regularizer is designed to be L 1 -norm replaced by a piecewise linear function. ISTA-Net [54] considers the regularizer as the L 1 -norm of a convolutional network. The network unrolls several phases iteratively, and each phase mimics one iteration of iterative shrinkage thresholding algorithm (ISTA) [55,56]. However, these networks only superficially mimic the corresponding optimization schemes, but they lack direct relations to the original optimization method or variational model and do not retain any convergence guarantee. In this work, we first develop a learnable optimization algorithm (LOA) for (3) with comprehensive convergence analysis and obtain an LOA-induced network by following the iterative scheme of the LOA exactly.

3.2. LOA-Induced Reconstruction Network

In this section, we first introduce a learned optimization algorithm (LOA) to solve (3) where the regularization network parameter Θ is fixed. As Θ is fixed in (3), we temporarily omit this in the derivation of the LOA below and write R ( x ; Θ ) as R ( x ) for notation simplicity.
In this work, to incorporate sparsity along with the learned features, we parameterize the function R ( x ) = κ · r ( x ) , where κ > 0 is a weight parameter that needs be chosen properly depending on the specific task (e.g., noise level, undersampling ratio, etc.), and r is a regularizer parameterized as a composition of neural networks and can be adapted to a broad range of imaging applications. Specifically, we parameterize r as the composition of the l 2 , 1 norm and a learnable feature extraction operator g ( x ) . That is, we set r in (11) to be
r ( x ) : = g ( x ) 2 , 1 = j = 1 m g j ( x ) .
Here, “:=” stands for “defined as”. g ( · ) = ( g 1 ( · ) , , g m ( · ) ) , g j ( · ) = g j ( · ; θ ) is parametrized as a convolutional neural network (CNN) for j = 1 , , m , and θ is the learned and fixed network parameter in r ( · ; θ ) , as mentioned above. We also consider κ to be learned and fixed as θ for now, and we discuss how to learn both of them in the next subsection. We use a smooth activation function in g as formulated in (20), which renders g a smooth but nonconvex function. Due to the nonsmooth · 2 , 1 , r is therefore a nonsmooth nonconvex function.
Since the minimization problem in (11) is nonconvex and nonsmooth, we need to derive an efficient LOA to solve it. Here, we first consider smoothing the l 2 , 1 norm that for any fixed g ( x ) is
r ε ( x ) = j = 1 m g j ( x ) 2 + ε 2 ε .
We denote R ε = κ · r ε . The LOA derived here is inspired by the proximal gradient descent algorithm and iterates the following steps to solve the smoothed problem:
z t + 1 = x t α t f ( x t )
x t + 1 = prox α t R ε t ( z t + 1 ) ,
where ε t denotes the smoothing parameter ε at the specific iteration t, and the proximal operator is defined as prox α g ( b ) : = arg min x x b + α g ( x ) in (6b). A quick observation from (5) is that R ε R as ε diminishes, so later we intentionally push ε t 0 at Line 16 in Algorithm 1. Then, one can readily show that R ε ( x ) R ( x ) R ε ( x ) + ε for all x and ε > 0 . From Algorithm 1, line 16 automatically reduces ε , and the iterates will converge to the solution of the original nonsmooth nonconvex problem (11)—this is clarified precisely in the convergence analysis in Appendix A.
Since R ε t is a complex function involving a deep neural network, its proximal operator does not have a closed form and cannot be computed easily in the subproblem in (6b). To overcome this difficulty, we consider to approximate R ε t by
R ^ ε t ( z t + 1 ) = R ε t ( z t + 1 ) + R ε t ( z t + 1 ) , x z t + 1 + 1 2 β t x z t + 1 2 .
Then, we update u t + 1 = prox α t R ^ ε t ( z t + 1 ) to replace (6b); therefore, we obtain
u t + 1 = z t + 1 τ t R ε t ( z t + 1 ) , where τ t = α t β t α t + β t .
In order to guarantee the convergence of the algorithm, we introduce the standard gradient descent of ϕ ε t (where ϕ ε t : = f + R ε t ) at x :
v t + 1 = arg min x f ( x t ) , x x t + R ε ( x t ) , x x t + 1 2 α t x x t 2 ,
which yields
v t + 1 = x t α t ϕ ε t ( x t ) ,
to serve as a safeguard for the convergence. Specifically, we set x t + 1 = u t + 1 if ϕ ε t ( u t + 1 ) ϕ ε t ( v t + 1 ) ; otherwise, we set x t + 1 = v t + 1 . Then, we repeat this process.
Our algorithm is summarized in Algorithm 1. The prior term with unknown parameters has the exact residual update itself which improves the learning and training process [57]. The condition checking in Line 5 is introduced to make sure that it is in the energy descending direction. Once the condition in Line 5 fails, the process moves to v t + 1 , and the line search in Line 12 guarantees that the appropriate step size can be achieved within finite steps to make the function value decrease. From Line 3 to Line 14, we consider that it solves a problem of minimizing ϕ ε t with ε t fixed. Line 15 is used to update the value of ε t depending on a reduction criterion. The detailed analysis of this mechanism and in-depth convergence justification is shown in Appendix A. The corresponding unrolling network exactly follows Algorithm 1 and thus shares the same convergence property. Compared to LDA [13], which computes both candidates u t + 1 , v t + 1 at every iteration and then chooses the one that achieves a smaller function value, we propose the criteria above in Line 5 for updating x t + 1 , which potentially saves extra computational time for calculating the candidate v t + 1 and potentially mitigates the frequent alternations between the two candidates. Besides, the smoothing method proposed in this work is more straightforward than smoothing in dual space [13] while still keeping provable convergence, as shown in Theorem A5.
The proposed LOA-induced network is a multi-phase network whose architecture exactly follows the proposed LOA (Algoirthm 1) in the way that each phase corresponds to one iteration of the algorithm. Specifically, we construct a deep network, denoted by F Θ , that follows Algorithm 1 exactly for a user-specified number of T iterations. Here, Θ denotes the set of learnable parameters in F Θ , which includes the regularization network parameter θ , weight κ , and other algorithmic parameters of Algorithm 1. Therefore, for any input under-sampled k-space measurement y , F Θ ( y ) executes the LOA (Algorithm 1) for T iterations and generates an approximate solution of the minimization problem (11):
F Θ ( y ) arg min x ϕ Θ ( x , y ) : = f ( x , y ) + R ( x ; Θ ) .
where we use “≈” since F Θ follows only finitely many steps of the optimization algorithm to approximate the solution. It is worth emphasizing that this approach can be readily applied to a much broader class of image reconstruction problems as long as f is (possibly nonconvex and) continuously differentiable with the Lipschitz continuous gradient. In the next subsection, we develop a meta-learning based approach for the robust training of the network parameter Θ .

3.3. Bilevel Optimization Algorithm for Network Training

In this section, we consider the parameter training problem of the LOA-induced network F Θ . Specifically, we develop a bilevel optimization algorithm to train our network parameters Θ from diverse data sets to improve network robustness and generalization ability.
Algorithm 1: Algorithmic Unrolling Method with Provable Convergence
1:
Input: Initial x 0 , 0 < ρ , γ < 1 , and ε 0 , σ > 0 . Max total phases T or tolerance ϵ tol > 0 .
2:
for t = 0 , 1 , 2 , , T 1 do
3:
z t + 1 = x t α t f ( x t )
4:
u t + 1 = z t + 1 τ t R ε t ( z t + 1 ) ,
5:
if ϕ ε t ( x t ) a u t + 1 x t and ϕ ε t ( u t + 1 ) ϕ ε t ( x t ) 1 a u t + 1 x t 2 then
6:
  set x t + 1 = u t + 1 ,
7:
else
8:
   v t + 1 = x t α t ϕ ε t ( x t ) ,
9:
  if ϕ ε t ( v t + 1 ) ϕ ε t ( x t ) 1 a v t + 1 x t 2 holds then
10:
   set x t + 1 = v t + 1 ,
11:
  else
12:
   update α t ρ α t , then go to 8,
13:
  end if
14:
end if
15:
if   ϕ ε t ( x t + 1 ) < σ γ ε t , set ε t + 1 = γ ε t ; otherwise, set ε t + 1 = ε t .
16:
if  σ ε t < ϵ tol , terminate.
17:
end for
18:
Output:  x t .
Recall that the LOA-induced network F Θ exactly follows Algorithm 1, which is designed to solve the variational model (11) containing learnable regularization R ( x ; Θ ) . As shown in Section 3.2, we design R ( x ; Θ ) = κ · r ( x ; Θ ) , where r is learned to capture the intrinsic property of the underlying common features across all different tasks. To account for the large variations in the diverse training/validation data sets, we introduce a task-specific parameter ω i to approximate the proper κ for the ith task. Specifically, for the ith task, the weight κ is set to σ ( ω i ) ( 0 , 1 ) , where σ ( · ) is the sigmoid function. Therefore, κ = σ ( ω i ) finds the proper weight of r for the i-th task according to its specific sampling ratio or pattern. The parameters ω i are to be optimized in conjunction with Θ through the hyperparameter tuning process below.
Suppose that we are given M data pairs { ( y m , x m * ) } m = 1 M for the use of training and validation, where y m is the observation, which is the partial k-space data in our setting, and x m * is the corresponding ground truth image. The data pairs are then sampled into N tasks { D τ i } i = 1 N , where each D τ i represents the collection of data pairs in the specific task τ i . In each task τ i , we further divide the data into the task-specific training set D τ i t r and validation set D τ i v a l . The architecture of our base network exactly follows the LOA (Algorithm 1) developed in the previous section with learnable parameters θ and a task-specific parameter ω i for the ith task. More precisely, for one data sample denoted by ( y j ( i ) , x j ( i ) ) in task τ i with index j, we propose the algorithmic unrolling network for task τ i as
F θ , ω i ( y j ( i ) ) arg min x f ( x , y j ( i ) ) + σ ( ω i ) r ( x ; θ ) ,
where θ denotes the learnable common parameters across different tasks with task-invariant representation, whereas ω i is a task-specific parameter for task τ i . The weight σ ( ω i ) represents the weight of r associated with the specific task τ i . In our proposed network, Θ is the collection of ( θ , ω i ) for task i = 1 , N . We denote ω to be the set { ω i } i = 1 N . The detailed architecture of this network is illustrated in Section 3.2. We define the task-specific loss
τ i ( θ , ω i ; D τ i ) : = j = 1 | D τ i | F θ , ω i ( y j ( i ) ) , x j ( i ) ,
where | D τ i | represents the cardinality of D τ i and
F θ , ω i ( y j ( i ) ) , x j ( i ) : = 1 2 F θ , ω i ( y j ( i ) ) x j ( i ) 2 .
For the sake of preventing the proposed model from overfitting the training data, we introduce a novel learning framework by formulating the network training as a bilevel optimization problem to learn ω and θ in (12) as
min θ , ω = { ω i : i [ N ] } i = 1 N τ i ( θ ( ω ) , ω i ; D τ i v a l )
s . t . θ ( ω ) = arg min θ i = 1 N τ i ( θ , ω i ; D τ i t r ) .
In (15), the lower-level optimization learns the task-invariant parameters θ of the feature encoder with the fixed task-specific parameter ω i on the training dataset, and the upper-level adjusts the task-specific parameters { ω i } so that the task-invariant parameters θ can perform robustly on the validation dataset as well. For simplicity, we omit the summation and redefine L ( θ , ω ; D ) : = i = 1 N τ i ( θ , ω ; D ) and then briefly rewrite (15) as
min θ , ω L ( θ ( ω ) , ω ; D v a l ) s . t . θ ( ω ) = arg min θ L ( θ , ω ; D t r ) .
Then, we relax (16) into a single-level constrained optimization where the lower-level problem is replaced with its first-order necessary condition following [58]
min θ , ω L ( θ , ω ; D v a l ) s . t . θ L ( θ , ω ; D t r ) = 0 .
which can be further approximated by an unconstrained problem by a penalty term as
min θ , ω { L ˜ ( θ , ω ; D t r , D v a l ) : = L ( θ , ω ; D v a l ) + λ 2 θ L ( θ , ω ; D t r ) 2 } .
We adopt the stochastic gradients of the loss functions on mini-batch data sets in each iteration. In our model, we need to include the data pairs of multiple tasks in one batch; therefore, we propose the cross-task mini-batches when training. At each training iteration, we randomly sample the training data pairs B τ i t r = { ( y j ( i ) , x j ( i ) ) D τ i t r } j = 1 J t r and the validation pairs B τ i v a l = { ( y j ( i ) , x j ( i ) ) D τ i v a l } j = 1 J v a l on each task τ i . Then, the overall training and validation mini-batches B t r and B v a l used in every training iteration are composed of the sampled data pairs from the entire set of tasks; i.e., B t r = i = 1 N { B τ i t r } and B v a l = i = 1 N { B τ i v a l } . Thus in each iteration, we have N · J t r and N · J v a l data pairs used for training and validation, respectively. To solve the minimization problem (17), we utilize the stochastic mini-batch alternating direction method summarized in Algorithm 2, which is modified from [58].
As analyzed in [58], this penalty-type method has linear time complexity without computing the Hessian of the low level and only requires a constant space since we only need to store the intermediate θ , ω at each training iteration, which is suitable for solving the large-scale bilevel optimization problem. Algorithm 2 relaxes the bi-level optimization problem to a single-level constrained optimization problem by using the first-order necessary condition, which is not equivalent to the original problem but is much easier and efficient to solve. In the inner-loop (Line 5–9) of Algorithm 2, we continue minimizing the converted single-level optimization function (18) with respect to θ for K steps and then ω once alternatively until the condition with tolerance δ in Line 5 fails. The basic idea behind the condition in Line 5 arises from the first-order necessary condition as we would like to push the gradient of L ˜ toward 0. Furthermore, at Line 11 of the outer loop (Line 2–11), we decrease the tolerance δ . Combining Line 5 and 11 guarantees that each time the inner loop terminates, the gradients of L ˜ with respect to θ and ω become increasingly close to 0. The parameter δ t o l is used to control the accuracy of the entire algorithm, and the outer-loop will terminate when δ is sufficiently small (i.e., δ δ t o l ). In addition, λ is the weight for the second constraint term of (18); in the beginning, we set λ to be small to achieve a quick starting convergence, then gradually increase its value to emphasize the constraint.
Algorithm 2: Stochastic mini-batch alternating direction penalty method to solve problem (17)
1:
Input D τ i t r , D τ i v a l , δ t o l > 0 .
2:
Initialize θ , ω , δ , λ > 0 and ν δ ( 0 , 1 ) , ν λ > 1 .
3:
while δ > δ t o l do
4:
    Sample cross-task training batch B t r = i = 1 N { ( y j ( i ) , x j ( i ) ) D τ i t r } j = 1 : J t r
5:
    Sample cross-task validation batch B v a l = i = 1 N { ( y j ( i ) , x j ( i ) ) D τ i v a l } j = 1 : J v a l
6:
    while  θ L ˜ ( θ , ω ; B t r , B v a l ) 2 + ω L ˜ ( θ , ω ; B t r , B v a l ) 2 > δ  do
7:
        for  k = 1 , 2 , , K (inner loop) do
8:
            θ θ ρ θ k θ L ˜ ( θ , ω ; B t r , B v a l )
9:
        end for
10:
         ω ω ρ ω ω L ˜ ( θ , ω ; B t r , B v a l )
11:
    end while
12:
    update  δ ν δ δ , λ ν λ λ
13:
end while
14:
output: θ , ω .

4. Implementation

4.1. Feature Extraction Operator

We set the feature extraction operator g to be a vanilla l-layer CNN with the component-wise nonlinear activation function φ and no bias, as follows:
g ( x ) = w l φ φ ( w 3 φ ( w 2 φ ( w 1 x ) ) ) ,
where { w q } q = 1 l denote the convolution weights consisting of d kernels with identical spatial kernel size, and ∗ denotes the convolution operation. Here, φ is constructed to be the smoothed rectified linear unit as defined below:
φ ( x ) = 0 , i f x δ , 1 4 δ x 2 + 1 2 x + δ 4 , i f δ < x < δ , x , i f x δ ,
where the prefixed parameter δ is set to be 0.001 in our experiment. The default configuration of the feature extraction operator is set as follows: the feature extraction operator g consists of l = 3 convolution layers and all convolutions are with 4 kernels of a spatial size of 3 × 3 .

4.2. Setups

As our method introduces an algorithmic unrolling network, there exists a one-to-one correspondence between the algorithm iterations and the neural network phases (or blocks). Each phase of the forward propagation can be viewed as one algorithm iteration, which motivates us to imitate the iterating of the optimization algorithm and use a stair training strategy [13]. At the first stage, we start training the network parameters using one phase, then after the the loss converges, we add more phases (one phase each time) then continue the training process. We repeat this procedure and stop it when the loss does not decrease any further when we add more blocks. We minimize the loss for 100 epochs/iterations each time using the SGD-based optimizer Adam [59] with β 1 = 0.9 , β 2 = 0.999 , and the initial learning rate set to 10 3 as well as a mini-batch size of 8. The Xavier Initializer [60] is used to initialize the weights of all convolutions. The initial smoothing parameter ε 0 is set to be 0.001 and then learned together with other network parameters. The input x 0 of the unrolling network is obtained by the zero-filling strategy [61]. The deep unrolling network was implemented using the Tensorflow toolbox [62] in the Python programming language.

5. Numerical Experiments

5.1. Dataset

To validate the performance of the proposed method, the data we used were from Multimodal Brain Tumor Segmentation Challenge 2018 [63], in which the training dataset contains four modalities (T1, T1 c , T2 and FLAIR )scanned from 285 patients and the validation dataset contains images from 66 patients, each with a volume size of 240 × 240 × 155 . Each modality consists of two types of gliomas: 75 volumes of low-grade gliomas (LGG) and 210 volumes of high-grade gliomas (HGG). Our implementation involved HGG MRI in two modalities—T1 and T2 images—and we chose 30 patients from each modality in the training dataset to train our network. In the validation dataset, we randomly picked 15 patients as our validation data and 6 patients in the training dataset as testing data, which were distinct from our training set and validation set. We cropped the 2D image size to be 160 × 180 in the center region and picked 10 adjacent slices in the center of each volume, resulting in a total of 300 images as our training data, 150 images as our validation data, and a total of 60 images as testing data. The amount of data mentioned here is for a single task, but since we emploedy multi-task training, the total number of images in each dataset should be multiplied by the number of tasks. For each 2D slice, we normalized the spatial intensity by dividing the maximum pixel value.

5.2. Experiment Settings

All the experiments were implemented on a Windows workstation with an Intel Core i9 CPU at 3.3GHz and an Nvidia GTX-1080Ti GPU with 11 GB of graphics card memory via TensorFlow [64]. The parameters in the proposed network were initialized by using Xavier initialization [65]. We trained the meta-learning network with four tasks synergistically associated with four different CS ratios—10%, 20%, 30%, and 40%—and tested the well-trained model on the testing dataset with the same masks of these four ratios. We used 300 training data for each CS ratio, amounting to a total of 1200 images in the training dataset. The results for T1 and T2 MR reconstructions are shown in Table 1 and Table 2, respectively. The associated reconstructed images are displayed in Figure 1 and Figure 2. We also tested the well-trained meta-learning model on unseen tasks with radial masks for unseen ratios of 15%, 25%, and 35% and random Cartesian masks with ratios of 10%, 20%, 30%, and 40%. The task-specific parameters for the unseen tasks were retrained for different masks with different sampling ratios individually with fixed task-invariant parameters θ . In this experiments, we only needed to learn ω i for three unseen CS ratios with radial mask and four regular CS ratios with Cartesian masks. The experimental training proceeded with fewer data and iterations, where we used 100 MR images with 50 epochs. For example, to reconstruct MR images with a CS ratio of 15% from the radial mask, we fixed the parameter θ and retrained the task-specific parameter ω on 100 raw data points with 50 epochs, then tested with renewed ω on our testing data set with raw measurements sampled from the radial mask with a CS radial of 15%. The results associated with radial masks are shown in Table 3 and Table 4, Figure 3 and Figure 4 for T1 and T2 images, respectively. The results associated with Cartesian masks are listed in Table 5 and reconstructed images are displayed in Figure 5.
We compared our proposed meta-learning method with conventional supervised learning, which was trained with one task at each time and only learned the task-invariant parameter θ without the task-specific parameter ω i . The forward network of conventional learning unrolled Algorithm 1 with 11 phases, which was the same as meta-learning. We merged the training set and validation set, resulting in 450 images for the training of the conventional supervised learning. The training batch size was set as 25 and we applied a total of 2000 epochs, while in meta-learning, we applied 100 epochs with a batch size of 8. The same testing set was used in both meta-learning and conventional learning to evaluate the performance of these two methods.
We made comparisons between meta-learning and the conventional network on the seven different CS ratios (10%, 20%, 30%, 40%, 15%, 25%, and 35%) in terms of two types of random under-sampling patterns: radial sampling mask and Cartesian sampling mask. The parameters for both meta-learning and conventional learning networks were trained via the Adam optimizer [59], and they both learned the forward unrolled task-invariant parameter θ . The network training of the conventional method used the same network configuration as the meta-learning network in terms of the number of convolutions, depth and size of CNN kernels, phase numbers and parameter initializer, etc. The major difference in the training process between these two methods is that meta-learning is performed for multi-tasks by leveraging the task-specific parameter ω i learned from Algorithm 2, and the common features among tasks are learned from the feed-forward network that unrolls Algorithm 1, while conventional learning solves the task-specific problem by simply unrolling the forward network via Algorithm 1, where both training and testing are implemented on the same task. To investigate the generalizability of meta-learning, we tested the well-trained meta-learning model on MR images in different distributions in terms of two types of sampling masks with various trajectories. The training and testing of conventional learning were applied with the same CS ratios; that is, if the conventional method was trained with a CS ratio 10%, then it was also tested on a dataset with a CS ratio of 10%, etc.
Because MR images are represented as complex values, we applied complex convolutions [66] for each CNN; that is, every kernel consisted of a real part and imaginary part. Three convolutions were used in g , where each convolution contained four filters with a spatial kernel size of 3 × 3 . In Algorithm 1, a total of 11 phases can be achieved if we set the termination condition ϵ tol = 1 × 10 3 , and the parameters of each phase are shared except for the step sizes. For the hyperparameters in Algorithm 1, we chose an initial learnable step size α 0 = 0.01 , τ 0 = 0.01 , ε 0 = 0.001 , and we set prefixed values of a = 10 5 , σ = 10 3 , ρ = 0.9 , and γ = 0.9 . The principle behind the choices of those parameters is based on the convergence of the algorithm and effectiveness of the computation. The parameter 0 < ρ < 1 is the reduction rate of the step size during the line search used to guarantee the convergence. The parameter 0 < γ < 1 at step 15 is the reduction rate for ε . In Algorithm 1, from step 2 to step 14, the smoothing level ε is fixed. When the gradient of the smoothed function is small enough, we reduce ε by a fraction factor γ to find an approximate accumulation point of the original nonsmooth nonconvex problem. We chose a larger a in order to have more iterations k for which u k + 1 satisfies the conditions in step 5, so that there would be fewer iterations requiring the computation of v k + 1 . Moreover, the scheme for computing u k + 1 is in accordance with the residual learning architecture that has been proven effective for reducing training error.
In Algorithm 2, we set ν δ = 0.95 and the parameter δ was initialized as δ 0 = 1 × 10 3 and stopped at value δ t o l = 4.35 × 10 6 , and a total of 100 epochs were performed. To train the conventional method, we set 2000 epochs with the same number of phases, convolutions, and kernel sizes as used to train the meta-learning approach. The initial λ was set as 1 × 10 5 and ν λ = 1.001 .
We evaluated our reconstruction results on the testing data sets using three metrics: peak signal-to-noise ratio (PSNR) [67], structural similarity (SSIM) [68], and normalized mean squared error (NMSE) [69]. The following formulations compute the PSNR, SSIM, and NMSE between the reconstructed image x and ground truth x * . PSNR can be induced by the mean square error (MSE) where
P S N R ( x , x * ) = 20 log 10 max ( | x * | ) M S E ( x , x * ) ,
where N is the total number of pixels of the ground truth and MSE is defined by M S E ( x , x * ) = 1 N x * x 2 .
S S I M ( x , x * ) = ( 2 μ x μ x * + C 1 ) ( 2 σ x x * + C 2 ) ( μ x 2 + μ x * 2 + C 1 ) ( σ x 2 + σ x * 2 + C 2 ) ,
where μ x , μ x * represent local means, σ x , σ x * denote standard deviations, σ x x * represents the covariance between x and x * , C 1 = ( k 1 L ) 2 , C 2 = ( k 2 L ) 2 are two constants which avoid the zero denominator, and k 1 = 0.01 , k 2 = 0.03 . L is the largest pixel value of MR image.
N M S E ( x , x * ) = x x * 2 2 x 2 2 ,
where NMSE is used to measure the mean relative error. For detailed information of these three metrics mentioned above, please refer to [67,68,69].

5.3. Experimental Results with Different CS Ratios in Radial Mask

In this section, we evaluate the performance of well-trained meta-learning and conventional learning approaches. Table 1, Table 2 and Table 5 report the quantitative results of averaged numerical performance with standard deviations and associated descaled task-specific meta-knowledge σ ( ω i ) . From the experiments implemented with radial masks, we observe that the average PSNR value of meta-learning improved by 1.54 dB in the T1 brain image for all four CS ratios compared with the conventional method, and for the T2 brain image, the average PSNR of meta-learning improved by 1.46 dB. Since the general setting of meta-learning aims to take advantage of the information provided from each individual task, with each task associated with an individual sampling mask that may have complemented sampled points, the performance of the reconstruction from each task benefits from other tasks. Smaller CS ratios will inhibit the reconstruction accuracy, due to the sparse undersampled trajectory in raw measurement, while meta-learning exhibits a favorable potential ability to solve this issue even in the situation of insufficient amounts of training data.
In general supervised learning, training data need to be in the same or a similar distribution; heterogeneous data exhibit different structural variations of features, which hinder CNNs from extracting features efficiently. In our experiments, raw measurements sampled from different ratios of compressed sensing display different levels of incompleteness; these undersampled measurements do not fall in the same distribution but they are related. Different sampling masks are shown at the bottom of Figure 1 and Figure 3, and these may have complemented sampled points, in the sense that some of the points which a 40 % sampling ratio mask did not sample were captured by other masks. In our experiment, different sampling masks provided their own information from their sampled points, meaning that four reconstruction tasks helped each other to achieve an efficient performance. Therefore, this explains why meta-learning is still superior to conventional learning when the sampling ratio is large.
Meta-learning expands a new paradigm for supervised learning—the purpose is to quickly learn multiple tasks. Meta-learning only learns task-invariant parameters once for a common feature that can be shared with four different tasks, and each σ ( ω i ) provides task-specific weighting parameters according to the principle of “learning to learn”. In conventional learning, the network parameter needs to be trained four times with four different masks since the task-invariant parameter cannot be generalized to other tasks, which is time-intensive. From Table 1 and Table 2, we observe that a small CS ratio needs a higher value of σ ( ω i ) . In fact, in our model (11), the task-specific parameters behave as weighted constraints for task-specific regularizers, and the tables indicate that lower CS ratios require larger weights to be applied for the regularization.
A qualitative comparison between conventional and meta-learning methods is shown in Figure 1 and Figure 2, displaying the reconstructed MR images of the same slice for T1 and T2, respectively. We label the zoomed-in details of HGG in the red boxes. We observe evidence that conventional learning is more blurry and loses sharp edges, especially with lower CS ratios. From the point-wise error map, we find that meta-learning has the ability to reduce noises, especially in some detailed and complicated regions, compared to conventional learning.
We also tested the performance of meta-learning with two-thirds of the training and validation data for the T1-weighted image used in the previous experiment, denoted as “meta-learning”. For conventional learning, the network was also trained by using two-thirds of the training samples in the previous experiment. The testing dataset remained the same as before. These results are displayed in Table 1, where we denote the reduced data experiments as “meta-learning * ” and “conventional * ”. These experiments reveal that the accuracy of test data decreases when we reduce the training data size, but it is not a surprise that meta-learnining * still outperforms conventional learning * , and even conventional learning.
To verify the reconstruction performance of the proposed LOA 1, we compared the proposed conventional learning with ISTA-Net + [54], which is a state-of-the-art deep unfolded network for MRI reconstruction. We retrained ISTA-Net + with the same training dataset and testing dataset as conventional learning on the T1-weighted image. For a fair comparison, we used the same number of convolution kernels, the same dimension of kernels for each convolution during training, and the same phase numbers as conventional learning. The testing numerical results are listed in Table 1 and the MRI reconstructions are displayed in Figure 1. We can observe that the conventional learning which unrolls Algorithm 1 outperforms ISTA-Net + in any of the CS ratios. From the corresponding point-wise absolute error, the conventional learning attains a much lower error and much better reconstruction quality.

5.4. Experimental Results with Different Unseen CS Ratios in Different Sampling Patterns

In this section, we test the generalizability of the proposed model for unseen tasks. We fixed the well-trained task-invariant parameter θ and only trained ω i for sampling ratios of 15%, 25%, and 35% with radial masks and sampling ratios of 10%, 20%, 30%, and 40% with Cartesian masks. In this experiment, we only used 100 training data points for each CS ratio and applied a total of 50 epochs. The averaged evaluation values and standard deviations are listed in Table 3 and Table 4 for reconstructed T1 and T2 brain images, respectively, with radial masks, and Table 5 shows the qualitative performance for the reconstructed T2 brain image with random Cartesian sampling masks applied. In the T1 image reconstruction results, meta-learning showed an improvement of 1.6921 dB in PSNR for the 15% CS ratio, 1.6608 dB for the 25% CS ratio, and 0.5764 dB for the 35% ratio compared to the conventional method, showing the tendency that the level of reconstruction quality for lower CS ratios improved more than higher CS ratios. A similar trend was found for T2 reconstruction results with different sampling masks. The qualitative comparisons are illustrated in Figure 3, Figure 4 and Figure 5 for T1 and T2 images tested with unseen CS ratios in radial masks and T2 images tested with Cartesian masks with regular CS ratios, respectively. In the experiments conducted with radial masks, meta-learning was superior to conventional learning, especially at a CS ratio of 15%—one can observe that the detailed regions in red boxes maintained their edges and were closer to the true image, while the conventional method reconstructions are hazier and lost details in some complicated tissues. The point-wise error map also indicates that meta-learning has the ability to suppress noises.
Training with Cartesian masks is more difficult than radial masks, especially for conventional learning, where the network is not very deep since the network only applies three convolutions each with four kernels. Table 5 indicates that the average performance of Meta-learning improved about 1.87 dB compared to conventional methods with T2 brain images. These results further demonstrate that meta-learning has the benefit of parameter efficiency, and the performance is much better than conventional learning even if we apply a shallow network with a small amount of training data.
The numerical experimental results discussed above show that meta-learning is capable of fast adaption to new tasks and has more robust generalizability for a broad range of tasks with heterogeneous, diverse data. Meta-learning can be considered as an efficient technique for solving difficult tasks by leveraging the features extracted from easier tasks.
Next, we empirically demonstrate the convergence of Algorithm 1 in Figure 6. This shows that the objective function value ϕ decreases and the PSNR value for testing data increases steadily as the number of phases increases, which indicates that the learned algorithm is indeed minimizing the learned function as we desired.

5.5. Future Work and Open Challenges

Deep optimization-based meta-learning techniques have shown great generalizability, but there are several open challenges that can be discussed and can potentially be addressed in future work. A major issue is the memorization problem, since the base learner needs to be optimized for a large number of phases and the training algorithm contains multiple gradient steps; furthermore, the computation is very expensive in terms of time and memory costs. In addition to reconstructing MRI through different trajectories, another potential application for medical imaging could be multi-modality reconstruction and synthesis. Capturing images of anatomy with multi-modality acquisitions enhances the diagnostic information and could be cast as a multi-task problem that could benefit from meta-learning.

6. Conclusions

In this paper, we put forward a novel deep model for MRI reconstructions via meta-learning. The proposed method has the ability to solve multi-tasks synergistically, and the well-trained model could generalize well to new tasks. Our baseline network is constructed by unfolding an LOA, which inherits the convergence property, improves the interpretability, and promotes the parameter efficiency of the designed network structure. The designated adaptive regularizer consists of a task-invariant learner and task-specific meta-knowledge. Network training follows a bilevel optimization algorithm that minimizes task-specific parameter ω in the upper level for the validation data and minimizes task-invariant parameters θ in the lower level for the training data with fixed ω . The proposed approach is the first model designed to solve the inverse problem by applying meta-training on the adaptive regularization in the variational model. We consider the recovery of undersampled raw data across different sampling trajectories with various sampling patterns as different tasks. Extensive numerical experiments on various MRI datasets demonstrate that the proposed method generalizes well at various sampling trajectories and is capable of fast adaption to unseen trajectories and sampling patterns. The reconstructed images achieve higher quality compared to conventional supervised learning for both seen and unseen k-space trajectory cases.

Author Contributions

All authors are contributed equally. Conceptualization, W.B., Y.C., X.Y. and Q.Z.; methodology, W.B., Y.C., X.Y. and Q.Z.; software, W.B. and Q.Z.; validation, Y.C. and X.Y.; formal analysis, W.B., Y.C., X.Y. and Q.Z.; investigation, W.B., Y.C., X.Y. and Q.Z.; resources, Y.C.; data curation, W.B. and Q.Z.; writing—original draft preparation, W.B. and Q.Z.; writing—review and editing, X.Y. and Y.C.; visualization, X.Y. and Y.C.; supervision, Y.C.; project administration, Y.C.; funding acquisition, X.Y. and Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by NSF grants DMS-1319050, DMS-1719932, DMS-1818886, DMS-1925263 and University of Florida AI Catalyst Grants.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on Multimodal Brain Tumor Segmentation Challenge (BraTs) 2018 [63] at https://www.med.upenn.edu/sbia/brats2018/data.html (accessed on 29 October 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Convergence Analysis

We make the following assumptions regarding f and g throughout this work:
  • (a1): f is differentiable and (possibly) nonconvex, and f is L f -Lipschitz continuous.
  • (a2): Every component of g is differentiable and (possibly) nonconvex, and g is L g -Lipschitz continuous.
  • (a3): sup x X g ( x ) M for some constant M > 0 .
  • (a4): ϕ is coercive, and ϕ * = min x X ϕ ( x ) > .
First, we state the Clark subdifferential [13] of r ( x ) in Lemma A1, and we show that the gradient of r ε is Lipschitz continuous in Lemma A2.
Lemma A1.
Let r ( x ) be defined in (4); then, the Clarke subdifferential of r at x is
r ( x ) = { j I 0 g i ( x ) w j + j I 1 g j ( x ) g j ( x ) g j ( x ) | w j R d , Π ( w j ; C ( g i ( x ) ) ) 1 , j I 0 } ,
where I 0 = { j [ m ] | g j ( x ) = 0 } , I 1 = [ m ] I 0 , and Π ( w ; C ( A ) ) is the projection of w onto C ( A ) which stands for the column space of A .
Lemma A2.
The gradient of r ε is Lipschitz continuous with constant m ( L g + 2 M 2 ε ) .
Proof. 
From r ε ( x ) = j = 1 m ( g j ( x ) 2 + ε 2 ) 1 2 ε , it follows that
r ε ( x ) = j = 1 m g j ( x ) g j ( x ) ( g j ( x ) 2 + ε 2 ) 1 2 .
For any x 1 , x 2 X , we first define h ( x ) = g j ( x ) ( g j ( x ) 2 + ε 2 ) 1 2 , so h ( x ) < 1 .
h ( x 1 ) h ( x 2 )
= g j ( x 1 ) g j ( x 1 ) 2 + ε 2 g j ( x 2 ) g j ( x 2 ) 2 + ε 2
= g j ( x 1 ) g j ( x 1 ) 2 + ε 2 g j ( x 1 ) g j ( x 2 ) 2 + ε 2 + g j ( x 1 ) g j ( x 2 ) 2 + ε 2 g j ( x 2 ) g j ( x 2 ) 2 + ε 2
g j ( x 1 ) g j ( x 2 ) 2 + ε 2 g j ( x 1 ) 2 + ε 2 g j ( x 1 ) 2 + ε 2 g j ( x 2 ) 2 + ε 2 + g j ( x 1 ) g j ( x 2 ) g j ( x 2 ) 2 + ε 2
g j ( x 1 ) g j ( x 1 ) 2 + ε 2 g j ( x 2 ) 2 + ε 2 g j ( x 1 ) 2 + ε 2 g j ( x 2 ) 2 + ε 2 + 1 ε g j ( x 1 ) g j ( x 2 )
1 ε g j ( x 2 ) 2 + ε 2 g j ( x 1 ) 2 + ε 2 + 1 ε g j ( x 1 ) g j ( x 2 )
1 ε g j ( x 2 ) 2 g j ( x 1 ) 2 g j ( x 2 ) 2 + ε 2 + g j ( x 1 ) 2 + ε 2 + 1 ε g j ( x 1 ) g j ( x 2 )
1 ε g j ( x 2 ) + g j ( x 1 ) g j ( x 2 ) 2 + ε 2 + g j ( x 1 ) 2 + ε 2 < 1 g j ( x 2 ) g j ( x 1 ) + 1 ε g j ( x 1 ) g j ( x 2 )
1 ε g j ( x 2 ) g j ( x 1 ) + 1 ε g j ( x 1 ) g j ( x 2 )
= 2 ε g j ( x 1 ) g j ( x 2 ) .
where to obtain (A3f) we used g j ( x 1 ) g j ( x 1 ) 2 + ε 2 < 1 and 1 g j ( x 1 ) 2 + ε 2 < 1 ε .
Therefore, we have
r ε ( x 1 ) r ε ( x 2 )
= j = 1 m g j ( x 1 ) h ( x 1 ) g j ( x 2 ) h ( x 2 )
= j = 1 m g j ( x 1 ) h ( x 1 ) g j ( x 2 ) h ( x 1 ) + g j ( x 2 ) h ( x 1 ) g j ( x 2 ) h ( x 2 )
j = 1 m g j ( x 1 ) g j ( x 2 ) h ( x 1 ) + g j ( x 2 ) h ( x 1 ) h ( x 2 ) j = 1 m g j ( x 1 ) g j ( x 2 ) h ( x 1 ) + g j ( x 2 ) h ( x 1 ) h ( x 2 )
j = 1 m g j ( x 1 ) g j ( x 2 ) + g j ( x 2 ) 2 ε g j ( x 1 ) g j ( x 2 ) by ( A 3 ) and
m ( L g x 1 x 2 + M 2 ε · M x 1 x 2 ) ,
where the first term of the last inequality is due to the L g -Lipschitz continuity of g j . The second term is because of g j ( x 1 ) g j ( x 2 ) = g j ( x ˜ ) ( x 1 x 2 ) for some x ˜ X due to the mean value theorem and g j ( x ˜ ) sup x X g j ( x ) M . Therefore, we obtain
r ε ( x 1 ) r ε ( x 2 ) m ( L g + 2 M 2 ε ) x 1 x 2 .
Lemma A3.
Let ε , η , τ t , a > 0 , 0 < ρ < 1 , and choose the initial x 0 X . Suppose the sequence { x t } is generated by executing Lines 3–14 of Algorithm 1 with fixed ε t = ε and that 0 < δ < L ε 2 / a + L ε < 1 exists such that α t δ L ε > 0 , where L ε = L f + L h + m L g + 2 M 2 ε and ϕ * : = min x X ϕ ( x ) . Then, the following statements hold:
  • ϕ ε ( x t ) 0 as t .
  • max { t N | ϕ ε ( x t + 1 ) η } max { a L ε 2 δ 2 , a 3 } { ϕ ε ( x 0 ) ϕ * + ε } .
Proof. 
1. In each iteration, we compute u t + 1 = z t + 1 τ t σ ( ω i ) r ε t ( z t + 1 ) .
1.1.
In the case the condition
ϕ ε ( x t ) a u t + 1 x t a n d ϕ ε ( u t + 1 ) ϕ ε ( x t ) 1 a u t + 1 x t 2
holds with a > 0 , we put x t + 1 = u t + 1 , and we have ϕ ε ( u t + 1 ) ϕ ε ( x t ) .
1.2.
Otherwise, we compute v t + 1 = x t α t ϕ ε ( x t ) , where α t is found through the line search until the criteria
ϕ ε ( v t + 1 ) ϕ ε ( x t ) 1 a v t + 1 x t 2
holds, and then put x t + 1 = v t + 1 . From Lemma A2, we know that the gradient r ε ( x ) is Lipschitz continuous with constant m ( L g + 2 M 2 ε ) . Furthermore, we assumed in (a1) that f is L f -Lipschitz continuous. Hence, putting L ε = L f + m ( L g + 2 M 2 ε ) , we find that ϕ ε is L ε -Lipschitz continuous, which implies
ϕ ε ( v t + 1 ) ϕ ε ( x t ) + ϕ ε ( x t ) , v t + 1 x t + L ε 2 v t + 1 x t 2 .
Furthermore, by the optimality condition of
v t + 1 = arg min x f ( x t ) , x x t + σ ( ω i ) r ε ( x t ) , x x t + 1 2 α t x x t 2 ,
we have
ϕ ε ( x t ) , v t + 1 x t + 1 2 α t v t + 1 x t 2 0 .
Combining (A8) and (A9) and v t + 1 = x t α t ϕ ε ( x t ) in line 8 of Algorithm 1 yields
ϕ ε ( v t + 1 ) ϕ ε ( x t ) 1 2 α t L ε 2 v t + 1 x t 2 .
Therefore, it is sufficient for α t 1 2 / a + L ε for the criteria (A7) to be satisfied. This process only take finitely many iterations since we can find a finite t such that ρ t α t 1 2 / a + L ε , and through the line search, we can obtain ϕ ε ( v t + 1 ) ϕ ε ( x t ) .
Therefore, in either case of 1.1. or 1.2. where we take x t + 1 = u t + 1 or v t + 1 , we can obtain
ϕ ε ( x t + 1 ) ϕ ε ( x t ) , for all t 0 .
Now, from case 1.1., (A6) gives
ϕ ε ( x t ) 2 a 2 u t + 1 x t 2 a 3 ϕ ε ( x t ) ϕ ε ( u t + 1 ) ,
therefore if x t + 1 = u t + 1 we get ϕ ε ( x t ) 2 a 3 ϕ ε ( x t ) ϕ ε ( x t + 1 ) .
From case 1.2. and v t + 1 = x t α t ϕ ε ( x t ) , we have
ϕ ε ( v t + 1 ) ϕ ε ( x t ) 1 a v t + 1 x t 2 = 1 a α t 2 ϕ ε ( x t ) 2
ϕ ε ( x t ) 2 a α t 2 ϕ ε ( x t ) ϕ ε ( v t + 1 ) ,
then if x t + 1 = v t + 1 , we have ϕ ε ( x t ) 2 a α t 2 ϕ ε ( x t ) ϕ ε ( x t + 1 ) .
Since δ L ε α t 1 2 / a + L ε , we have
ϕ ε ( v t + 1 ) 2 a L ε 2 δ 2 ϕ ε ( v t + 1 ) ϕ ε ( x t + 1 ) .
Combining (A12b) and (A14) and selecting C = max { a L ε 2 δ 2 , a 3 } , we obtain
ϕ ε ( x t ) 2 C ϕ ε ( x t ) ϕ ε ( x t + 1 ) .
Summing up (A15) for t = 0 , , T , we have
t = 0 T ϕ ε ( x t ) 2 C ϕ ε ( x 0 ) ϕ ε ( x T + 1 ) .
Combined with the fact that ϕ ε ( x ) ϕ ( x ) ε ϕ * ε for every x X , we have
t = 0 T ϕ ε ( x t ) 2 C ϕ ε ( x 0 ) ϕ * + ε .
The right-hand side is a finite constant, and hence by letting t , we know that ϕ ε ( x t ) 0 , which proves the first statement.
2. Denote κ : = max { t N | ϕ ε ( x t + 1 ) η } ; then, we know that ϕ ε ( x t + 1 ) η for all t κ 1 . Hence, we have
κ η 2 t = 0 κ 1 ϕ ε ( x t + 1 ) 2 = t = 1 κ ϕ ε ( x t ) 2 C ϕ ε ( x 0 ) ϕ * + ε .
which implies the second statement. □
Lemma A4.
Suppose that the sequence { x t } is generated by Algorithm 1 with an initial guess x 0 . Then, for any t 0 , we have ϕ ε t + 1 ( x t + 1 ) + ε t + 1 ϕ ε t ( x t ) + ε t .
Proof. 
To prove this statement, we can prove
ϕ ε t + 1 ( x t + 1 ) + ε t + 1 ϕ ε t ( x t + 1 ) + ε t ϕ ε t ( x t ) + ε t .
The second inequality is immediately obtained from (A11). Now, we prove the first inequality.
For any ε > 0 , denote
r ε , j ( x ) = g j ( x ) 2 2 + ε 2 ε .
Since ϕ ε ( x ) = f ( x ) + σ ( ω i ) j = 1 m r ε , j ( x ) , it suffices to show that
r ε t + 1 , j ( x t + 1 ) + ε t + 1 r ε t , j ( x t + 1 ) + ε t
If ε t + 1 = ε t , then the two quantities above are identical and the first inequality holds. Now, suppose ε t + 1 = γ ε t ε t ; then,
r ε t + 1 , j ( x t + 1 ) + ε t + 1 = g j ( x ) 2 2 + ε t + 1 2 g j ( x ) 2 2 + ε t 2 = r ε t , j ( x t + 1 ) + ε t ,
which implies the first inequality of (A19). □
Theorem A5.
Suppose that { x t } is the sequence generated by Algorithm 1 with any initial x 0 , ϵ tol = 0 and T = . Let { x t l + 1 } be the subsequence that satisfies the reduction criterion in step 15 of Algorithm 1, i.e., ϕ ε t l ( x t l + 1 ) σ ε t l γ for t = t l and l = 1 , 2 , . Then { x t l + 1 } has at least one accumulation point, and every accumulation point of { x t l + 1 } is a clarke stationary point of min x ϕ ( x ) : = f ( x ) + σ ( ω i ) r ( x ) .
Proof. 
By Lemma A4 and ϕ ( x ) ϕ ε ( x ) + ε for all ε > 0 and x X , we know that
ϕ ( x t ) ϕ ε t ( x t ) + ε t ϕ ε 0 ( x 0 ) + ε 0 < .
Since ϕ is coercive, we know that { x t } is bounded, and the selected subsequence { x t l + 1 } is also bounded and has at least one accumulation point.
Note that ϕ ε t l ( x t l + 1 ) σ ε t l γ = σ ε 0 γ l + 1 0 as l . Let { x p + 1 } be any convergent subsequence of { x t l + 1 } and denote ε p as the corresponding ε t used in Algorithm 1 that generates x p + 1 . Then, there exists x * X such that x p + 1 x * as ε p 0 , and ϕ ε p ( x p + 1 ) 0 as p .
Note that the Clarke subdifferential of ϕ at x * is given by ϕ ( x * ) = f ( x * ) + σ ( ω i ) r ( x * ) :
ϕ ( x ^ ) = { f ( x ^ ) + σ ( ω i ) j I 0 g j ( x ^ ) w j + σ ( ω i ) j I 1 g j ( x ^ ) g j ( x ^ ) g j ( x ^ ) | Π ( w j ; C ( g j ( x ^ ) ) ) 1 , j I 0 } ,
where I 0 = { j [ m ] | g i ( x ^ ) = 0 } and I 1 = [ m ] I 0 .
If j I 0 , we have g j ( x ) = 0 g j ( x ) = 0 : then,
r ε ( x ) = j I 0 g j ( x ) g j ( x ) g j ( x ) 2 + ε 2 1 2 + j I 1 g j ( x ) g j ( x ) g j ( x ) 2 + ε 2 1 2
= 0 + j I 1 g j ( x ) g j ( x ) g j ( x ) 2 + ε 2 1 2
Therefore, we obtain
ϕ ε p ( x p + 1 ) = f ( x p + 1 ) + σ ( ω i ) j I 1 g j ( x ) g j ( x ) g j ( x ) 2 + ε p 2 1 2 .
Comparing (A24) and (A26), we can see that the first term on the right-hand side of (A26) converges to that of (A24), due to the fact that x p + 1 x ^ and the continuity of f . Together with the continuity of g i and g i , the last term of (A24) converges to the last term of (A26) as ε p 0 and g j ( x ) > 0 . Furthermore, apparently 0 is a special case of the second term in (A26). Hence, we know that
dist ( ϕ ε p ( x p + 1 ) , ϕ ( x ^ ) ) 0 ,
as p . Since ϕ ε p ( x p + 1 ) 0 and ϕ ( x ^ ) is closed, we conclude that 0 ϕ ( x ^ ) . □

References

  1. Li, D.; Yang, Y.; Song, Y.Z.; Hospedales, T.M. Learning to generalize: Meta-learning for domain generalization. In Proceedings of the Thirty-Second AAAI Conference on Artificial Intelligence, New Orleans, LA, USA, 2–7 February 2018. [Google Scholar]
  2. Quiñonero-Candela, J.; Sugiyama, M.; Lawrence, N.D.; Schwaighofer, A. Dataset Shift in Machine Learning; Mit Press: Cambridge, MA, USA, 2009. [Google Scholar]
  3. Ben-David, S.; Blitzer, J.; Crammer, K.; Pereira, F. Analysis of representations for domain adaptation. Adv. Neural Inf. Process. Syst. 2007, 19, 137. [Google Scholar]
  4. Balaji, Y.; Sankaranarayanan, S.; Chellappa, R. Metareg: Towards domain generalization using meta-regularization. Adv. Neural Inf. Process. Syst. 2018, 31, 998–1008. [Google Scholar]
  5. Fan, J.; Han, F.; Liu, H. Challenges of big data analysis. Natl. Sci. Rev. 2014, 1, 293–314. [Google Scholar] [CrossRef] [Green Version]
  6. Day, O.; Khoshgoftaar, T.M. A survey on heterogeneous transfer learning. J. Big Data 2017, 4, 1–42. [Google Scholar] [CrossRef]
  7. Munkhdalai, T.; Yu, H. Meta networks. In Proceedings of the International Conference on Machine Learning, Sydney, Australia, 6–11 August 2017; pp. 2554–2563. [Google Scholar]
  8. Finn, C.; Abbeel, P.; Levine, S. Model-agnostic meta-learning for fast adaptation of deep networks. In Proceedings of the International Conference on Machine Learning, Sydney, NSW, Australia, 6–11 August 2017; pp. 1126–1135. [Google Scholar]
  9. Rusu, A.A.; Rao, D.; Sygnowski, J.; Vinyals, O.; Pascanu, R.; Osindero, S.; Hadsell, R. Meta-Learning with Latent Embedding Optimization. In Proceedings of the International Conference on Learning Representations, Vancouver, BC, Canada, 30 April–3 May 2018. [Google Scholar]
  10. Yao, H.; Huang, L.K.; Zhang, L.; Wei, Y.; Tian, L.; Zou, J.; Huang, J. Improving generalization in meta-learning via task augmentation. In Proceedings of the International Conference on Machine Learning, Virtual Event, 18–24 July 2021; pp. 11887–11897.
  11. Thrun, S.; Pratt, L. Learning to learn: Introduction and overview. In Learning to Learn; Springer: Berlin/Heidelberg, Germany, 1998; pp. 3–17. [Google Scholar]
  12. Hospedales, T.M.; Antoniou, A.; Micaelli, P.; Storkey, A.J. Meta-Learning in Neural Networks: A Survey. arXiv 2020, arXiv:2004.05439. [Google Scholar]
  13. Chen, Y.; Liu, H.; Ye, X.; Zhang, Q. Learnable Descent Algorithm for Nonsmooth Nonconvex Image Reconstruction. SIAM J. Imaging Sci. 2021, 14, 1532–1564. [Google Scholar] [CrossRef]
  14. Yao, H.; Wu, X.; Tao, Z.; Li, Y.; Ding, B.; Li, R.; Li, Z. Automated relational meta-learning. arXiv 2020, arXiv:2001.00745. [Google Scholar]
  15. Lee, Y.; Choi, S. Gradient-based meta-learning with learned layerwise metric and subspace. In Proceedings of the International Conference on Machine Learning, Stockholm, Sweden, 10–15 July 2018; pp. 2927–2936. [Google Scholar]
  16. Huisman, M.; van Rijn, J.N.; Plaat, A. A survey of deep meta-learning. Artif. Intell. Rev. 2021, 54, 4483–4541. [Google Scholar] [CrossRef]
  17. Koch, G.; Zemel, R.; Salakhutdinov, R. Siamese neural networks for one-shot image recognition. In Proceedings of the ICML Deep Learning Workshop, Lille, France, 6– 11 July 2015; Volume 2. [Google Scholar]
  18. Vinyals, O.; Blundell, C.; Lillicrap, T.; Wierstra, D. Matching networks for one shot learning. Adv. Neural Inf. Process. Syst. 2016, 29, 3630–3638. [Google Scholar]
  19. Snell, J.; Swersky, K.; Zemel, R.S. Prototypical networks for few-shot learning. arXiv 2017, arXiv:1703.05175. [Google Scholar]
  20. Mishra, N.; Rohaninejad, M.; Chen, X.; Abbeel, P. A simple neural attentive meta-learner. arXiv 2017, arXiv:1707.03141. [Google Scholar]
  21. Ravi, S.; Larochelle, H. Optimization as a Model for Few-Shot Learning. 2016. Available online: https://openreview.net/forum?id=rJY0-Kcll (accessed on 29 October 2021).
  22. Qiao, S.; Liu, C.; Shen, W.; Yuille, A.L. Few-shot image recognition by predicting parameters from activations. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–23 June 2018; pp. 7229–7238. [Google Scholar]
  23. Graves, A.; Wayne, G.; Danihelka, I. Neural turing machines. arXiv 2014, arXiv:1410.5401. [Google Scholar]
  24. Rajeswaran, A.; Finn, C.; Kakade, S.; Levine, S. Meta-Learning with Implicit Gradients. 2019. Available online: https://papers.nips.cc/paper/2019/hash/072b030ba126b2f4b2374f342be9ed44-Abstract.html (accessed on 29 October 2021).
  25. Li, Z.; Zhou, F.; Chen, F.; Li, H. Meta-sgd: Learning to learn quickly for few-shot learning. arXiv 2017, arXiv:1707.09835. [Google Scholar]
  26. Antoniou, A.; Edwards, H.; Storkey, A. How to train your MAML. arXiv 2018, arXiv:1810.09502. [Google Scholar]
  27. Nichol, A.; Achiam, J.; Schulman, J. On first-order meta-learning algorithms. arXiv 2018, arXiv:1803.02999. [Google Scholar]
  28. Finn, C.; Rajeswaran, A.; Kakade, S.; Levine, S. Online meta-learning. In Proceedings of the International Conference on Machine Learning, Long Beach, CA, USA, 9–15 June 2019; pp. 1920–1930. [Google Scholar]
  29. Grant, E.; Finn, C.; Levine, S.; Darrell, T.; Griffiths, T. Recasting gradient-based meta-learning as hierarchical bayes. arXiv 2018, arXiv:1801.08930. [Google Scholar]
  30. Finn, C.; Xu, K.; Levine, S. Probabilistic model-agnostic meta-learning. arXiv 2018, arXiv:1806.02817. [Google Scholar]
  31. Yoon, J.; Kim, T.; Dia, O.; Kim, S.; Bengio, Y.; Ahn, S. Bayesian model-agnostic meta-learning. In Proceedings of the 32nd International Conference on Neural Information Processing Systems, Montréal, QC, Canada, 3–8 December 2018; pp. 7343–7353. [Google Scholar]
  32. Vuorio, R.; Sun, S.H.; Hu, H.; Lim, J.J. Multimodal model-agnostic meta-learning via task-aware modulation. arXiv 2019, arXiv:1910.13616. [Google Scholar]
  33. Yao, H.; Wei, Y.; Huang, J.; Li, Z. Hierarchically structured meta-learning. In Proceedings of the International Conference on Machine Learning, Virtual, 18–24 July 2021; pp. 7045–7054. [Google Scholar]
  34. Yin, M.; Tucker, G.; Zhou, M.; Levine, S.; Finn, C. Meta-Learning without Memorization. arXiv 2019, arXiv:1912.03820. [Google Scholar]
  35. Jenni, S.; Favaro, P. Deep bilevel learning. In Proceedings of the European Conference on Computer Vision (ECCV), Munich, Germany, 8–14 September 2018; pp. 618–633. [Google Scholar]
  36. Chen, W.Y.; Liu, Y.C.; Kira, Z.; Wang, Y.C.; Huang, J.B. A Closer Look at Few-shot Classification. In Proceedings of the International Conference on Learning Representations, New Orleans, LA, USA, 6–9 May 2019. [Google Scholar]
  37. Li, Y.; Yang, Y.; Zhou, W.; Hospedales, T. Feature-critic networks for heterogeneous domain generalization. In Proceedings of the International Conference on Machine Learning, Long Beach, CA, USA, 9–15 June 2019; pp. 3915–3924. [Google Scholar]
  38. Rebuffi, S.A.; Bilen, H.; Vedaldi, A. Learning multiple visual domains with residual adapters. In Proceedings of the Advances in Neural Information Processing Systems, Long Beach, CA, USA, 4–9 December 2017. [Google Scholar]
  39. Triantafillou, E.; Zhu, T.; Dumoulin, V.; Lamblin, P.; Evci, U.; Xu, K.; Goroshin, R.; Gelada, C.; Swersky, K.J.; Manzagol, P.A.; et al. Meta-Dataset: A Dataset of Datasets for Learning to Learn from Few Examples. In Proceedings of the International Conference on Learning Representations, Addis Ababa, Ethiopia, 26–30 April 2020. [Google Scholar]
  40. Yu, T.; Quillen, D.; He, Z.; Julian, R.; Hausman, K.; Finn, C.; Levine, S. Meta-world: A benchmark and evaluation for multi-task and meta reinforcement learning. In Proceedings of the Conference on Robot Learning, Virtual, 16–18 November 2020; pp. 1094–1100. [Google Scholar]
  41. Haldar, J.P.; Hernando, D.; Liang, Z.P. Compressed-sensing MRI with random encoding. IEEE Trans. Med Imaging 2010, 30, 893–903. [Google Scholar] [CrossRef] [Green Version]
  42. Blomgren, P.; Chan, T. Color TV: Total variation methods for restoration of vector-valued images. IEEE Trans. Image Process. 1998, 7, 304–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Lundervold, A.S.; Lundervold, A. An overview of deep learning in medical imaging focusing on MRI. Z. Für Med. Phys. 2019, 29, 102–127. [Google Scholar] [CrossRef]
  44. Liang, D.; Cheng, J.; Ke, Z.; Ying, L. Deep magnetic resonance image reconstruction: Inverse problems meet neural networks. IEEE Signal Process. Mag. 2020, 37, 141–151. [Google Scholar] [CrossRef]
  45. Sandino, C.M.; Cheng, J.Y.; Chen, F.; Mardani, M.; Pauly, J.M.; Vasanawala, S.S. Compressed sensing: From research to clinical practice with deep neural networks: Shortening scan times for magnetic resonance imaging. IEEE Signal Process. Mag. 2020, 37, 117–127. [Google Scholar] [CrossRef] [PubMed]
  46. McCann, M.T.; Jin, K.H.; Unser, M. Convolutional neural networks for inverse problems in imaging: A review. IEEE Signal Process. Mag. 2017, 34, 85–95. [Google Scholar] [CrossRef] [Green Version]
  47. Zhou, S.K.; Greenspan, H.; Davatzikos, C.; Duncan, J.S.; van Ginneken, B.; Madabhushi, A.; Prince, J.L.; Rueckert, D.; Summers, R.M. A Review of Deep Learning in Medical Imaging: Image Traits, Technology Trends, Case Studies with Progress Highlights, and Future Promises. 2020. [CrossRef]
  48. Singha, A.; Thakur, R.S.; Patel, T. Deep Learning Applications in Medical Image Analysis. Biomed. Data Min. Inf. Retr. Methodol. Tech. Appl. 2021, 6, 293–350. [Google Scholar]
  49. Chandra, S.S.; Bran Lorenzana, M.; Liu, X.; Liu, S.; Bollmann, S.; Crozier, S. Deep learning in magnetic resonance image reconstruction. J. Med Imaging Radiat. Oncol. 2021, 5, 564–577. [Google Scholar] [CrossRef]
  50. Ahishakiye, E.; Van Gijzen, M.B.; Tumwiine, J.; Wario, R.; Obungoloch, J. A survey on deep learning in medical image reconstruction. Intell. Med. 2021. Available online: https://www.sciencedirect.com/science/article/pii/S2667102621000061 (accessed on 29 October 2021). [CrossRef]
  51. Hammernik, K. Learning a variational network for reconstruction of accelerated MRI data. Magn. Reson. Med. 2018, 79, 3055–3071. [Google Scholar] [CrossRef]
  52. Aggarwal, H.K.; Mani, M.P.; Jacob, M. MoDL: Model-Based Deep Learning Architecture for Inverse Problems. IEEE Trans. Med Imaging 2019, 38, 394–405. [Google Scholar] [CrossRef] [PubMed]
  53. Sun, J.; Li, H.; Xu, Z. Deep ADMM-Net for compressive sensing MRI. Adv. Neural Inf. Process. Syst. 2016, 29, 10–18. [Google Scholar]
  54. Zhang, J.; Ghanem, B. ISTA-Net: Interpretable optimization-inspired deep network for image compressive sensing. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 23 June 2018; pp. 1828–1837. [Google Scholar]
  55. Lions, P.L.; Mercier, B. Splitting algorithms for the sum of two nonlinear operators. SIAM J. Numer. Anal. 1979, 16, 964–979. [Google Scholar] [CrossRef]
  56. Beck, A.; Teboulle, M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems. SIAM J. Imaging Sci. 2009, 2, 183–202. [Google Scholar] [CrossRef] [Green Version]
  57. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 30 June 2016; pp. 770–778. [Google Scholar]
  58. Mehra, A.; Hamm, J. Penalty method for inversion-free deep bilevel optimization. arXiv 2019, arXiv:1911.03432. [Google Scholar]
  59. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. arXiv 2015, arXiv:1412.6980. [Google Scholar]
  60. Glorot, X.; Bengio, Y. Understanding the difficulty of training deep feedforward neural networks. In Proceedings of the International Conference on Artificial Intelligence and Statistics. Society for Artificial Intelligence and Statistics, Virtual Event, 13–15 April 2021. [Google Scholar]
  61. Bernstein, M.A.; Fain, S.B.; Riederer, S.J. Effect of windowing and zero-filled reconstruction of MRI data on spatial resolution and acquisition strategy. J. Magn. Reson. Imaging Off. J. Int. Soc. Magn. Reson. Med. 2001, 14, 270–280. [Google Scholar] [CrossRef] [Green Version]
  62. Abadi, M.; Agarwal, A.; Barham, P.; Brevdo, E.; Chen, Z.; Citro, C.; Corrado, G.S.; Davis, A.; Dean, J.; Devin, M.; et al. TensorFlow: Large-Scale Machine Learning on Heterogeneous Systems. 2015. Software. Available online: tensorflow.org (accessed on 29 October 2021).
  63. Menze, B.H.; Jakab, A.; Bauer, S.; Kalpathy-Cramer, J.; Farahani, K.; Kirby, J.; Burren, Y.; Porz, N.; Slotboom, J.; Wiest, R.; et al. The multimodal brain tumor image segmentation benchmark (BRATS). IEEE Trans. Med Imaging 2014, 34, 1993–2024. [Google Scholar] [CrossRef] [PubMed]
  64. Abadi, M.; Barham, P.; Chen, J.; Chen, Z.; Davis, A.; Dean, J.; Devin, M.; Ghemawat, S.; Irving, G.; Isard, M.; et al. Tensorflow: A system for large-scale machine learning. In Proceedings of the 12th {USENIX} Symposium on Operating Systems Design and Implementation ({OSDI} 16), Savannah, GA, USA, 2–4 November 2016; pp. 265–283. [Google Scholar]
  65. 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, Sardinia, Italy, 13–15 May 2010; pp. 249–256. [Google Scholar]
  66. Wang, S.; Cheng, H.; Ying, L.; Xiao, T.; Ke, Z.; Zheng, H.; Liang, D. DeepcomplexMRI: Exploiting deep residual network for fast parallel MR imaging with complex convolution. Magn. Reson. Imaging 2020, 68, 136–147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Hore, A.; Ziou, D. Image quality metrics: PSNR vs. SSIM. In Proceedings of the 2010 20th International Conference on Pattern Recognition, Istanbul, Turkey, 23–26 August 2010; pp. 2366–2369. [Google Scholar]
  68. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [Green Version]
  69. Poli, A.; Cirillo, M. On the use of the normalized mean square error in evaluating dispersion model performance. Atmos. Environ. Part A Gen. Top. 1993, 27, 2427–2434. [Google Scholar] [CrossRef]
Figure 1. The pictures (from top to bottom) display the reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning, conventional learning, and ISTA-Net + with four different CS ratios of 10%, 20%, 30%, 40% (from left to right). The top-right image is the ground truth fully-sampled image.
Figure 1. The pictures (from top to bottom) display the reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning, conventional learning, and ISTA-Net + with four different CS ratios of 10%, 20%, 30%, 40% (from left to right). The top-right image is the ground truth fully-sampled image.
Jimaging 07 00231 g001
Figure 2. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learningwith four different CS ratios of 10%, 20%, 30%, 40% (from left to right). The top-right iage is the ground truth fully-sampled image.
Figure 2. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learningwith four different CS ratios of 10%, 20%, 30%, 40% (from left to right). The top-right iage is the ground truth fully-sampled image.
Jimaging 07 00231 g002
Figure 3. The pictures (from top to bottom) display the T1 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learning. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with three different unseen CS ratios of 15%, 25%, and 35% (from left to right). Conventional learning was trained and tested with the same CS ratios of 15%, 25%, and 35%. The top-right image is the ground truth fully-sampled image.
Figure 3. The pictures (from top to bottom) display the T1 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learning. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with three different unseen CS ratios of 15%, 25%, and 35% (from left to right). Conventional learning was trained and tested with the same CS ratios of 15%, 25%, and 35%. The top-right image is the ground truth fully-sampled image.
Jimaging 07 00231 g003
Figure 4. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learning. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and test edwith three different unseen CS ratios of 15%, 25%, and 35% (from left to right). Conventional learning was trained and tested with the same CS ratios of 15%, 25%, and 35%. The top-right image is the ground truth fully-sampled image.
Figure 4. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated radial masks for meta-learning and conventional learning. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and test edwith three different unseen CS ratios of 15%, 25%, and 35% (from left to right). Conventional learning was trained and tested with the same CS ratios of 15%, 25%, and 35%. The top-right image is the ground truth fully-sampled image.
Jimaging 07 00231 g004
Figure 5. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated Cartesian masks for meta-learning and conventional learningwith four different CS ratios of 10%, 20%, 30%, and 40% (from left to right). The top-right image is the ground truth fully-sampled image.
Figure 5. The pictures (from top to bottom) display the T2 brain image reconstruction results, zoomed-in details, point-wise errors with a color bar, and associated Cartesian masks for meta-learning and conventional learningwith four different CS ratios of 10%, 20%, 30%, and 40% (from left to right). The top-right image is the ground truth fully-sampled image.
Jimaging 07 00231 g005
Figure 6. Convergence behavior of Algorithm 1 on T1 weighted MRI image reconstruction with four different CS ratios using radial mask. Left: Objective function value ϕ versus phase number. Right: PSNR value versus phase number.
Figure 6. Convergence behavior of Algorithm 1 on T1 weighted MRI image reconstruction with four different CS ratios using radial mask. Left: Objective function value ϕ versus phase number. Right: PSNR value versus phase number.
Jimaging 07 00231 g006
Table 1. Quantitative evaluations of the reconstructions of T1 brain image associated with various sampling ratios of radial masks. Conventional * and meta-learning * are trained with two-thirds of the dataset used in training conventional and meta-learning approaches, respectively.
Table 1. Quantitative evaluations of the reconstructions of T1 brain image associated with various sampling ratios of radial masks. Conventional * and meta-learning * are trained with two-thirds of the dataset used in training conventional and meta-learning approaches, respectively.
CS RatioMethodsPSNRSSIMNMSE σ ( ω i )
ISTA-Net + [54]21.2633 ± 1.03170.5487 ± 0.04400.1676 ± 0.0253
Conventional * 21.6947 ± 1.02640.5689 ± 0.04040.1595 ± 0.0240
10%Conventional21.7570 ± 1.06770.5650 ± 0.04120.0259 ± 0.0082
Meta-learning * 22.9633 ± 1.09690.5962 ± 0.04150.0194 ± 0.00650.9339
Meta-learning23.2672 ± 1.12290.6101 ± 0.04360.0184 ± 0.00670.9218
ISTA-Net + [54]26.2734 ± 1.01150.7068 ± 0.03640.0944 ± 0.0155
Conventional * 26.4639 ± 1.02330.7107 ± 0.03570.0924 ± 0.0154
20%Conventional26.6202 ± 1.16620.7121 ± 0.03970.0910 ± 0.0169
Meta-learning * 27.9381 ± 1.11210.7541 ± 0.03600.0063 ± 0.00230.8150
Meta-learning28.2944 ± 1.21190.7640 ± 0.03770.0058 ± 0.00220.7756
ISTA-Net + [54]28.8309 ± 1.31370.7492 ± 0.04070.0708 ± 0.0142
Conventional * 29.2923 ± 1.31940.7522 ± 0.03990.0671 ± 0.0136
30%Conventional29.5034 ± 1.44460.7557 ± 0.04080.0657 ± 0.0143
Meta-learning * 30.8691 ± 1.58970.8310 ± 0.03940.0033 ± 0.00150.6359
Meta-learning31.1417 ± 1.58660.8363 ± 0.03850.0031 ± 0.00140.6501
ISTA-Net + [54]30.7282 ± 1.54820.8008 ± 0.04280.0572 ± 0.0127
Conventional * 31.3761 ± 1.58920.8035 ± 0.04200.0532 ± 0.0121
40%Conventional31.4672 ± 1.63900.8111 ± 0.04220.0029 ± 0.0014
Meta-learning * 32.7330 ± 1.63860.8623 ± 0.03580.0022 ± 0.00100.6639
Meta-learning32.8238 ± 1.70390.8659 ± 0.03700.0022 ± 0.00100.6447
Table 2. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of radial masks.
Table 2. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of radial masks.
CS RatioMethodsPSNRSSIMNMSE σ ( ω i )
10%Conventional23.0706 ± 1.24690.5963 ± 0.03490.2158 ± 0.0347
Meta-learning24.0842 ± 1.38630.6187 ± 0.03800.0112 ± 0.01170.9013
20%Conventional27.0437 ± 1.06130.6867 ± 0.02610.1364 ± 0.0213
Meta-learning28.9118 ± 1.07170.7843 ± 0.02400.0122 ± 0.00300.8742
30%Conventional29.5533 ± 1.09270.7565 ± 0.02650.1023 ± 0.0166
Meta-learning31.4096 ± 0.98140.8488 ± 0.02170.0069 ± 0.00190.8029
40%Conventional32.0153 ± 0.94020.8139 ± 0.02380.0770 ± 0.0128
Meta-learning33.1114 ± 1.01890.8802 ± 0.02100.0047 ± 0.00150.7151
Table 3. Quantitative evaluations of the reconstructions of T1 brain image associated with various sampling ratios of radial masks. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with unseen ratios of 15%, 25%, and 35%. The conventional method was subjected to regular training and testing with the same CS ratios of 15%, 25%, and 35%.
Table 3. Quantitative evaluations of the reconstructions of T1 brain image associated with various sampling ratios of radial masks. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with unseen ratios of 15%, 25%, and 35%. The conventional method was subjected to regular training and testing with the same CS ratios of 15%, 25%, and 35%.
CS RatioMethodsPSNRSSIMNMSE σ ( ω i )
15%Conventional24.6573 ± 1.02440.6339 ± 0.03820.1136 ± 0.0186
Meta-learning26.3494 ± 1.01020.7088 ± 0.03520.0090 ± 0.00300.9429
25%Conventional28.4156 ± 1.23610.7533 ± 0.03680.0741 ± 0.0141
Meta-learning30.0764 ± 1.46450.8135 ± 0.03800.0040 ± 0.00170.8482
35%Conventional31.5320 ± 1.52420.7923 ± 0.04200.0521 ± 0.0119
Meta-learning32.1084 ± 1.64810.8553 ± 0.03790.0025 ± 0.00110.6552
Table 4. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of radial masks. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with unseen ratios of 15%, 25%, and 35%. Conventional method was subjected to regular training and testing with the same CS ratios of 15%, 25%, and 35%.
Table 4. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of radial masks. Meta-learning was trained with CS ratios of 10%, 20%, 30%, and 40% and tested with unseen ratios of 15%, 25%, and 35%. Conventional method was subjected to regular training and testing with the same CS ratios of 15%, 25%, and 35%.
CS RatioMethodsPSNRSSIMNMSE σ ( ω i )
15%Conventional24.8921 ± 1.23560.6259 ± 0.02850.1749 ± 0.0280
Meta-learning26.7031 ± 1.25530.7104 ± 0.03180.0205 ± 0.00520.9532
25%Conventional29.0545 ± 1.19800.7945 ± 0.02920.1083 ± 0.0173
Meta-learning30.0698 ± 0.99690.8164 ± 0.02350.0093 ± 0.00220.8595
35%Conventional31.5201 ± 1.00210.7978 ± 0.02360.0815 ± 0.0129
Meta-learning32.0683 ± 0.92040.8615 ± 0.02090.0059 ± 0.00140.7388
Table 5. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of random Cartesian masks.
Table 5. Quantitative evaluations of the reconstructions of T2 brain image associated with various sampling ratios of random Cartesian masks.
CS RatioMethodsPSNRSSIMNMSE σ ( ω i )
10%Conventional20.8867 ± 1.29990.5082 ± 0.04750.0796 ± 0.0242
Meta-learning22.0434 ± 1.35550.6279 ± 0.04440.0611 ± 0.01880.9361
20%Conventional22.7954 ± 1.28190.6057 ± 0.04120.0513 ± 0.0157
Meta-learning24.7162 ± 1.39190.6971 ± 0.03800.0329 ± 0.01010.8320
30%Conventional24.2170 ± 1.23960.6537 ± 0.03600.0371 ± 0.0117
Meta-learning26.4537 ± 1.34710.7353 ± 0.03400.0221 ± 0.00680.6771
40%Conventional25.3668 ± 1.32790.6991 ± 0.02880.1657 ± 0.0265
Meta-learning27.5367 ± 1.41070.7726 ± 0.02970.0171 ± 0.00500.6498
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bian, W.; Chen, Y.; Ye, X.; Zhang, Q. An Optimization-Based Meta-Learning Model for MRI Reconstruction with Diverse Dataset. J. Imaging 2021, 7, 231. https://doi.org/10.3390/jimaging7110231

AMA Style

Bian W, Chen Y, Ye X, Zhang Q. An Optimization-Based Meta-Learning Model for MRI Reconstruction with Diverse Dataset. Journal of Imaging. 2021; 7(11):231. https://doi.org/10.3390/jimaging7110231

Chicago/Turabian Style

Bian, Wanyu, Yunmei Chen, Xiaojing Ye, and Qingchao Zhang. 2021. "An Optimization-Based Meta-Learning Model for MRI Reconstruction with Diverse Dataset" Journal of Imaging 7, no. 11: 231. https://doi.org/10.3390/jimaging7110231

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop