Search

In this paper, we propose the use of Recurrent Inference Machines (RIMs) to perform T1 and T2 mapping. The RIM is a neural network framework that learns an iterative inference process based on the signal model, similar to conventional statistical methods for quantitative MRI (QMRI), such as the Maximum Likelihood Estimator (MLE). This framework combines the advantages of both data-driven and model-based methods, and, we hypothesize, is a promising tool for QMRI. Previously, RIMs were used to solve linear inverse reconstruction problems. Here, we show that they can also be used to optimize non-linear problems and estimate relaxometry maps with high precision and accuracy. The developed RIM framework is evaluated in terms of accuracy and precision and compared to an MLE method and an implementation of the Residual Neural Network (ResNet). The results show that the RIM improves the quality of estimates compared to the other techniques in Monte Carlo experiments with simulated data, test-retest analysis of a system phantom, and in-vivo scans. Additionally, inference with the RIM is 150 times faster than the MLE, and robustness to (slight) variations of scanning parameters is demonstrated. Hence, the RIM is a promising and flexible method for QMRI. Coupled with an open-source training data generation tool, it presents a compelling alternative to previous methods.

Background: T2 mapping is increasingly used to quantify cartilage degeneration in knee osteoarthritis (OA), yet reproducibility studies in a multicenter setting are limited. The purpose of this study was to determine the longitudinal reproducibility and multicenter variation of cartilage T2 mapping, using various MRI equipment and acquisition protocols. Methods: In this prospective multicenter study, four traveling, healthy human subjects underwent T2 mapping twice at five different centers with a 6-month-interval. Centers had various MRI scanners, field strengths, and T2 mapping acquisition protocols. Mean T2 values were calculated in six cartilage regions of interest (ROIs) as well as an average value per patient. A phantom was scanned once at each center. To evaluate longitudinal reproducibility, intraclass correlation coefficients (ICC), root-mean-square coefficient of variation (RMS-CV), and a Bland-Altman plot were used. To assess the variation of in vivo and phantom T2 values across centers, ANOVA was performed. Results: ICCs of the T2 mapping measurements per ROI and the ROI's combined ranged from 0.73 to 0.91, indicating good to excellent longitudinal reproducibility. RMS-CVs ranged from 1.1% to 1.5% (per ROI) and 0.6% to 1.6% (ROIs combined) across the centers. A Bland-Altman plot did not reveal a systematic error. Evident, but consistent, discrepancies in T2values were observed across centers, both in vivo and in the phantom. Conclusions: The results of this study suggest that T2mapping can be used to longitudinal assess cartilage degeneration in multicenter studies. Given the differences in absolute cartilag

Purpose: To improve image quality of multi-contrast imaging with the proposed Autocalibrated Parallel Imaging Reconstruction for Extended Multi-Contrast Imaging (APIR4EMC). Methods: APIR4EMC reconstructs multi-contrast images in an autocalibrated parallel imaging reconstruction framework by adding contrasts as virtual coils. Compensation of signal evolution along the echo train of different contrasts is performed to improve signal prediction for missing samples. As a proof of concept, we performed prospectively accelerated phantom and in-vivo brain acquisitions with T1, T1-fat saturated (Fatsat), T2, PD, and FLAIR contrasts. The k-space sampling patterns of these acquisitions were jointly optimized. Images were jointly reconstructed with the proposed APIR4EMC method as well as individually with GRAPPA. Root mean square error (RMSE) to fully sampled reference images and g-factor maps were computed for both methods in the phantom experiment. Visual evaluation was performed in the in-vivo experiment. Results: Compared to GRAPPA, APIR4EMC reduced artifacts and improved SNR of the reconstructed images in the phantom acquisitions. Quantitatively, APIR4EMC substantially reduced noise amplification (g-factor) as well as RMSE compared to GRAPPA. Signal evolution compensation reduced artifacts. In the in-vivo experiments, 1 mm3 isotropic 3D images with contrasts of T1, T1-Fatsat, T2, PD, and FLAIR were acquired in as little as 7.5 min with the acceleration factor of 9. Reconstruction quality was consistent with the phantom results. Conclusion: Compared to single contrast reconstruction with GRAPPA, APIR4EMC reduces artifacts and noise amplification in accelerated multi-contrast imaging.

Background: Patellofemoral pain (PFP) is a common knee condition and possible precursor of knee osteoarthritis (OA). Inflammation, leading to an increased perfusion, or increased volume of the infrapatellar fat pad (IPFP) may induce knee pain. The aim of the study was to compare quantitative dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) parameters, as imaging biomarkers of inflammation, and volume of the IPFP between patients with PFP and controls and between patients with and without IPFP edema or joint effusion. Methods: Patients with PFP and healthy controls were included and underwent non-fat suppressed 3D fast-spoiled gradient-echo (FSPGR) and DCE-MRI. Image registration was applied to correct for motion. The IPFP was delineated on FSPGR using Horos software. Volume was calculated and quantitative perfusion parameters were extracted by fitting extended Tofts' pharmacokinetic model. Differences in volume and DCE-MRI parameters between patients and controls were tested by linear regression analyses. IPFP edema and effusion were analyzed identically. Results: Forty-three controls and 35 PFP patients were included. Mean IPFP volume was 26.04 (4.18) mL in control subjects and 27.52 (5.37) mL in patients. Median Ktrans was 0.017 (0.016) min-1 in control subjects and 0.016 (0.020) min-1 in patients. None of the differences in volume and perfusion parameters were statistically significant. Knees with effusion showed a higher perfusion of the IPFP compared to knees without effusion in patients only. Conclusions: The IPFP has been implicated as source of knee pain, but higher DCE-MR blood perfusion, an imaging biomarker of inflammation, and larger volume are not associated with PFP. Patient's knees with effusion showed a higher perfusion, pointing towards inflammation.

Purpose: Quantification of the T2 ∗ relaxation time constant is relevant in various magnetic resonance imaging applications. Mono- or bi-exponential models are typically used to determine these parameters. However, in case of complex, heterogeneous tissues these models could lead to inaccurate results. We compared a model, provided by the fractional-order extension of the Bloch equation with the conventional models. Methods: Axial 3D ultra-short echo time (UTE) scans were acquired using a 3.0 T MRI and a 16-channel surface coil. After image registration, voxel-wise T2 ∗ was quantified with mono-exponential, bi-exponential and fractional-order fitting. We evaluated all three models repeatability and the bias of their derived parameters by fitting at various noise levels. To investigate the effect of the SNR for the different models, a Monte-Carlo experiment with 1000 repeats was performed for different noise levels for one subject. For a cross-sectional investigation, we used the mean fitted values of the ROIs in five volunteers. Results: Comparing the mono-exponential and the fractional order T2 ∗ maps, the fractional order fitting method yielded enhanced contrast and an improved delineation of the different tissues. In the case of the bi-exponential method, the long T2 ∗ component map demonstrated the anatomy clearly with high contrast. Simulations showed a nonzero bias of the parameters for all three mathematical models. ROI based fitting showed that the T2 ∗ values were different depending on the applied method, and they differed most for the patellar tendon in all subjects. Conclusions: In high SNR cases, the fractional order and bi-exponential models are both performing well with low bias. However, in all observed cases, one of the bi-exponential components has high standard deviation in T2 ∗. The bi-exponential model is suitable for T2 ∗ mapping, but we recommend using the fractional order model for cases of low SNR.

Purpose: To optimize the diffusion-weighting b values and postprocessing pipeline for hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region. Methods: Optimized diffusion-weighting b value sets ranging between 5 and 30 b values were constructed by optimizing the Cramér-Rao lower bound of the hybrid intravoxel incoherent motion diffusion kurtosis imaging model. With this model, the perfusion fraction, pseudodiffusion coefficient, diffusion coefficient, and kurtosis were estimated. Sixteen volunteers were scanned with a reference b value set and 3 repeats of the optimized sets, of which 1 with volunteers swallowing on purpose. The effects of (1) b value optimi

Background: Quantitative MRI of patellar tendinopathy (PT) can be challenging due to spatial variation of T2* relaxation times. Purpose: 1) To compare T2* quantification using a standard approach with analysis in specific tissue compartments of the patellar tendon. 2) To evaluate test–retest reliability of different methods for fitting ultrashort echo time (UTE)-relaxometry data. Study Type: Prospective. Subjects: Sixty-five athletes with PT. Field Strength/Sequence: 3D UTE scans covering the patellar tendon were acquired using a 3.0T scanner and a 16-channel surface coil. Assessment: Voxelwise median T2* was quantified with monoexponential, fractional-order, and biexponential fitting. We applied two methods for T2* analysis: first, a standard approach by analyzing all voxels covering the proximal patellar tendon. Second, within subregions of the patellar tendon, by using thresholds on biexponential fitting parameter percentage short T2* (0–30% for mostly long T2*, 30–60% for mixed T2*, and 60–100% for mostly short T2*). Statistical Tests: Average test–retest reliability was assessed in three athletes using coefficients-of-variation (CV) and coefficients-of-repeatability (CR). Results: With standard image analysis, we found a median [interquartile range, IQR] monoexponential T2* of 6.43 msec [4.32–8.55] and fractional order T2* 4.39 msec [3.06–5.78]. The percentage of short T2* components was 52.9% [35.5–69.6]. Subregional monoexponential T2* was 13.78 msec [12.11–16.46], 7.65 msec [6.49–8.61], and 3.05 msec [2.52–3.60] and fractional order T2* 11.82 msec [10.09–14.44], 5.14 msec [4.25–5.96], and 2.19 msec [1.82–2.64] for 0–30%, 30–60%, and 60–100% short T2*, respectively. Biexponential component short T2* was 1.693 msec [1.417–2.003] for tissue with mostly short T2* and long T2* of 15.79 msec [13.47–18.61] for mostly long T2*. The average CR (CV) was 2 msec (15%), 2 msec (19%) and 10% (22%) for monoexponential, fractional order and percentage short T2*, respectivel

Objective: Infrapatellar fat pad (IPFP) fat-suppressed T2 (T2FS) hyperintense regions on MRI are an important imaging feature of knee osteoarthritis (OA) and are thought to represent inflammation. These regions are also common in non-OA subjects, and may not always be linked to inflammation. Our aim was to evaluate quantitative blood perfusion parameters, as surrogate measure of inflammation, within T2FS-hyperintense regions in patients with OA, with patellofemoral pain (PFP) (supposed OA precursor), and control subjects. Methods: Twenty-two knee OA patients, 35 PFP patients and 43 healthy controls were included and underwent MRI, comprising T2 and DCE-MRI sequences. T2FS-hyperintense IPFP regions were delineated and a reference region was drawn in adjacent IPFP tissue with normal signal intensity. After fitting the extended Tofts pharmacokinetic model, quantitative DCE-MRI perfusion parameters were compared between the two regions within subjects in each subgroup, using a paired Wilcoxon signed-rank test. Results: T2FS-hyperintense IPFP regions were present in 16 of 22 (73%) OA patients, 13 of 35 (37%) PFP patients, and 14 of 43 (33%) controls. DCE-MRI perfusion parameters were significantly different between regions with and without a T2FS-hyperintense signal in OA patients, demonstrating higher Ktrans compared to normal IFPF tissue (0.039 min−1 versus 0.025 min−1, p = 0.017) and higher Ve (0.157 versus 0.119, p = 0.010). For PFP patients and controls no significant differences were found. Conclusions: IPFP T2FS-hyperintense regions are associated with higher perfusion in knee OA patients in contrast to identically appearing regions in PFP patients and controls, pointing towards an inflammatory pathogenesis in OA only. Key Points: • Morphologically identical appearing T2FS-hyperintense infrapatellar fat pad regions show different perfusion in healthy subjects, subjects with patellofemoral pain, and subjects with knee osteoarthritis. • Elevated DCE-MRI perfusion pa

Objective: Subchondral bone changes, characterized by increased bone turnover and vascularity, are believed to stimulate progression and pain in knee osteoarthritis (OA). The objective of this study was to evaluate the bone perfusion in knee OA using quantitative dynamic contrast enhanced MRI (DCE-MRI). Design: Unicompartmental knee OA patients were included and underwent 3 Tesla DCE-MRI and T2-weighted MRI. Quantitative DCE-MRI analysis of Ktrans and Kep, representing perfusion parameters, was performed to evaluate differences between the most and least affected knee compartment. First, DCE-MRI parameter differences between epimetaphyseal and subchondral bone in both femur and tibia were assessed. Second, DCE-MRI parameters in subchondral bone marrow lesions (BMLs) were compared to surrounding subchondral bone without BMLs. Results: Twenty-three patients were analyzed. Median Ktrans and Kep in epimetaphyseal bone were significantly higher (p < 0.05) in the most affected (Ktrans: 0.014; Kep: 0.054 min−1) compared to least affected (Ktrans: 0.010; Kep: 0.016 min−1) compartment. For subchondral bone, DCE-MRI parameters were significantly higher (p < 0.05) in the most affected (Ktrans: 0.019; Kep: 0.091 min−1) compared to least affected (Ktrans: 0.014; Kep: 0.058 min−1) compartment as well. Subchondral BMLs detected on fat-saturated T2-weighted images were present in all patients. Median Ktrans (0.091 vs 0.000 min−1) and Kep (0.258 vs 0.000 min−1) were significantly higher within subchondral BMLs compared to surrounding subchondral bone without BMLs (p < 0.001). Conclusions: Increased perfusion parameters in epimetaphyseal bone, subchondral bone and BMLs are observed in unicompartmental knee OA. BMLs likely account for most of the effect of the higher bone perfusion in knee OA.

Purpose: To reduce artifacts and scan time of GRASE imaging by selecting an optimal sampling pattern and jointly reconstructing gradient echo and spin echo images. Methods: We jointly reconstruct images for the different echo types by considering these as additional virtual coil channels in the novel Autocalibrated Parallel Imaging Reconstruction with Sampling Pattern Optimization for GRASE (APIR4GRASE) method. Besides image reconstruction, we identify optimal sampling patterns for the acquisition. The selected optimal patterns were validated on phantom and in-vivo acquisitions. Comparison to the conventional GRASE without acceleration, and to the GRAPPA reconstruction with a single echo type was also performed. Results: Using identified optimal sampling patterns, APIR4GRASE eliminated modulation artifacts in both phantom and in-vivo experiments; mean square error (MSE) was reduced by 78% and 94%, respectively, compared to the conventional GRASE with similar scan time. Both artifacts and g-factor were reduced compared to the GRAPPA reconstruction with a single echo type. Conclusion: APIR4GRASE substantially improves the speed and quality of GRASE imaging over the state-of-the-art, and is able to reconstruct both spin echo and gradient echo images.

Wall shear stress (WSS), the frictional force exerted on endothelial cells by blood flow, is hypothesised to influence atherosclerotic plaque growth and composition. We developed a methodology for image registration of MR and histology images of advanced human carotid plaques and corresponding WSS data, obtained by MRI and computational fluid dynamics. The image registration method requires four types of input images, in vivo MRI, ex vivo MRI, photographs of transversally sectioned plaque tissue and histology images. These images are transformed to a shared 3D image domain by applying a combination of rigid and non-rigid registration algorithms. Transformation matrices obtained from registration of these images are used to transform subject-specific WSS data to the shared 3D image domain as well. WSS values originating from the 3D WSS map are visualised in 2D on the corresponding lumen locations in the histological sections and divided into eight radial segments. In each radial segment, the correlation between WSS values and plaque composition based on histological parameters can be assessed. The registration method was successfully applied to two carotid endarterectomy specimens. The resulting matched contours from the imaging modalities had Hausdorff distances between 0.57 and 0.70 mm, which is in the order of magnitude of the in vivo MRI resolution. We simulated the effect of a mismatch in the rigid registration of imaging modalities on WSS results by relocating the WSS data with respect to the stack of histology images. A 0.6 mm relocation altered the mean WSS values projected on radial bins on average by 0.59 Pa, compared to the output of original registration. This mismatch of one image slice did not change the correlation between WSS and plaque thickness. In conclusion, we created a method to invest

Purpose: High-resolution three-dimensional (3D) structural MRI is useful for delineating complex or small structures of the body. However, it requires long acquisition times and high SAR, limiting its clinical use. The purpose of this work is to accelerate the acquisition of high-resolution images by combining compressed sensing and parallel imaging (CSPI) on a 3D-GRASE sequence and to compare it with a (CS)PI 3D-FSE sequence. Several sampling patterns were investigated to assess their influence on image quality. Methods: The proposed k-space sampling patterns are based on two undersampled k-space grids, variable density (VD) Poisson-disc, and VD pseudo-random Gaussian, and five different trajectories described in the literature. Bloch simulations are performed to obtain the transform point spread function and evaluate the coherence of each sampling pattern. Image resolution was assessed by the full-width at half-maximum (FWHM). Prospective CSPI 3D-GRASE phantom and in vivo experiments in knee and brain are carried out to assess image quality, SNR, SAR, and acquisition time compared to PI 3D-GRASE, PI 3D-FSE, and CSPI 3D-FSE acquisitions. Results: Sampling patterns with VD Poisson-disc obtain the lowest coherence for both PD-weighted and T2 -weighted acquisitions. VD pseudo-random Gaussian obtains lower FWHM, but higher sidelobes than VD Poisson-disc. CSPI 3D-GRASE reduces acquisition time (43% for PD-weighted and 40% for T2 -weighted) and SAR (∼45% for PD-weighted and T2 -weighted) compared to CSPI 3D-FSE. Conclusions: CSPI 3D-GRASE reduces acquisition time compared to a CSPI 3DFSE acquisition, preserving image quality. The design of the sampling pattern is crucial for image quality in CSPI 3D-GRASE image acquisitions.

Purpose: To identify the optimal combination of pharmacokinetic model and arterial input function (AIF) for quantitative analysis of blood perfusion in the patellar bone using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Materials and Methods: This method design study used a random subset of five control subjects from an Institutional Review Board (IRB)-approved case–control study into patellofemoral pain, scanned on a 3T MR system with a contrast-enhanced time-resolved imaging of contrast kinetics (TRICKS) sequence. We systematically investigated the reproducibility of pharmacokinetic parameters for all combinations of Orton and Parker AIF models with Tofts, Extended Tofts (ETofts), and Brix pharmacokinetic models. Furthermore, we evaluated if the AIF should use literature parameters, be subject-specific, or group-specific. Model selection was based on the goodness-of-fit and the coefficient of variation of the pharmacokinetic parameters inside the patella. This extends previous studies that were not focused on the patella and did not evaluate as many combinations of arterial and pharmacokinetic models. Results: The vascular component in the ETofts model could not reliably be recovered (coefficient of variation [CV] of vp >50%) and the Brix model parameters showed high variability of up to 20% for kel across good AIF models. Compared to group-specific AIF, the subject-specific AIF's mostly had higher residual. The best reproducibility and goodness-of-fit were obtained by combining Tofts' pharmacokinetic model with the group-specific Parker AIF. Conclusion: We identified several good combinations of pharmacokinetic models and AIF for quantitative analysis of perfusion in the patellar bone. The recommended combination is Tofts pharmacokinetic model combined with a group-specific Parker AIF model. Level of Evidence: 2. Technical Efficacy: Stage 1. J. Magn. Reson. Imaging 2018;47:848–858.

The most widespread technique used to register sets of medical images consists of selecting one image as fixed reference, to which all remaining images are successively registered. This pairwise scheme requires one optimization procedure per pair of images to register. Pairwise mutual information is a common dissimilarity measure applied to a large variety of datasets. Alternative methods, called groupwise registrations, have been presented to register two or more images in a single optimization procedure, without the need of a reference image. Given the success of mutual information in pairwise registration, we adapt one of its multivariate versions, called total correlation, in a groupwise context. We justify the choice of total correlation among other multivariate versions of mutual information, and provide full implementation details. The resulting total correlation measure is remarkably close to measures previously proposed by Huizinga et al. based on principal component analysis. Our experiments, performed on five quantitative imaging datasets and on a dynamic CT imaging dataset, show that total correlation yields registration results that are comparable to Huizinga’s methods. Total correlation has the advantage of being theoretically justified, while the measures of Huizinga et al. were designed empirically. Additionally, total correlation offers an alternative to pairwise mutual information on quantitative imaging datasets.

Background: Altered perfusion might play an important role in the pathophysiology of patellofemoral pain (PFP), a common knee complaint with unclear pathophysiology. Purpose: To investigate differences in dynamic contrast-enhanced (DCE)-MRI perfusion parameters between patients with PFP and healthy control subjects. Population/Subjects/Phantom/Specimen

A large number of mathematical models have been proposed to describe the measured signal in diffusion-weighted (DW) magnetic resonance imaging (MRI). However, model comparison to date focuses only on specific subclasses, e.g. compartment models or signal models, and little or no information is available in the literature on how performance varies among the different types of models. To address this deficiency, we organized the ‘White Matter Modeling Challenge’ during the International Symposium on Biomedical Imaging (ISBI) 2015 conference. This competition aimed to compare a range of different kinds of models in their ability to explain a large range of measurable in vivo DW human brain data. Specifically, we assessed the ability of models to predict the DW signal accurately for new diffusion gradients and b values. We did not evaluate the accuracy of estimated model parameters, as a ground truth is hard to obtain. We used the Connectome scanner at the Massachusetts General Hospital, using gradient strengths of up to 300 mT/m and a broad set of diffusion times. We focused on assessing the DW signal prediction in two regions: the genu in the corpus callosum, where the fibres are relatively straight and parallel, and the fornix, where the configuration of fibres is more complex. The challenge participants had access to three-quarters of the dataset and their models were ranked on their ability to predict the remaining unseen quarter of the data. The challenge provided a unique opportunity for a quantitative comparison of diverse methods from multiple groups worldwide. The comparison of the challenge entries reveals interesting trends that could potentially influence the next generation of diffusion-based quantitative MRI techniques. The first is that signal models do not necessarily outperform tissue models; in fact, of those tested, tissue models rank highest on average. The second is that assuming a non-Gaussian (rather than purely Gaussian) noise model provides lit

The goal of this paper is to develop a method for assessment of microstructural properties of the fornix in conventional (low resolution, single non-zero b-value) diffusion-weighted magnetic resonance imaging (DW-MRI) data. For this purpose, a bi-tensor model, comprising of an isotropic and an anisotropic diffusion compartment, was fitted to the diffusion-weighted images. Two subject-specific constraints were studied to so

Purpose This paper presents and studies a framework for reliable modeling of diffusion MRI using a data-acquisition adaptive prior. Methods Automated relevance determination estimates the mean of the posterior distribution of a rank-2 dual tensor model exploiting Jeffreys prior (JARD). This data-acquisition prior is based on the Fisher information matrix and enables the assessment whether two tensors are mandatory to describe the data. The method is compared to Maximum Likelihood Estimation (MLE) of the dual tensor model and to FSL's ball-and-stick approach. Results Monte Carlo experiments demonstrated that JARD's volume fractions correlated well with the ground truth for single and crossing fiber configurations. In single fiber configurations JARD automatically reduced the volume fraction of one compartment to (almost) zero. The variance in fractional anisotropy (FA) of the main tensor component was thereby reduced compared to MLE. JARD and MLE gave a comparable outcome in data simulating crossing fibers. On brain data, JARD yielded a smaller spread in FA along the corpus callosum compared to MLE. Tract-based spatial statistics demonstrated a higher sensitivity in detecting age-related white matter atrophy using JARD compared to both MLE and the ball-and-stick approach. Conclusions The proposed framework offers accurate and precise estimation of diffusion properties in single and dual fiber regions.

Purpose To evaluate the reproducibility and sensitivity of the modified CINE inversion recovery (mCINE-IR) acquisition on rats for measuring the myocardial T1 at 7 Tesla. Materials and Methods The recently published mCINE-IR acquisition on humans was applied on rats for the first time, enabling the possibility of translational studies with an identical sequence. Simulations were used to study signal evolution and heart rate dependency. Gadolinium phantoms, a heart specimen and a healthy rat were used to study reproducibility. Two cryo-infarcted rats were scanned to measure late gadolinium enhancement (LGE). Results In the phantom reproducibility studies the T1 measurements had a maximum coefficient of variation (COV) of 1.3%. For the in vivo reproducibility the COV was below 5% in the anterior cardiac segments. In simulations with phantoms and specimens, a heart rate dependency of approximately 0.5 ms/bpm was present. The T1 maps of the cryo-infarcted rats showed a clear lowering of T1 in de LGE region. Conclusion The results show that mCINE-IR is highly reproducible and that the sensitivity allows detecting T1 changes in the rat myocardium.

In small animal imaging studies, when the locations of the micro-structures of interest are unknown a priori, there is a simultaneous need for full-body coverage and high resolution. In MRI, additional requirements to image contrast and acquisition time will often make it impossible to acquire such images directly. Recently, a resolution enhancing postprocessing technique called super-resolution reconstruction (SRR) has been demonstrated to improve visualization and localization of micro-structures in small animal MRI by combining multiple low-resolution acquisitions. However, when the field-of-view is large relative to the desired voxel size, solving the SRR problem becomes very expensive, in terms of both memory requirements and computation time. In this paper we introduce a novel local approach to SRR that aims to overcome the computational problems and allow researchers to efficiently explore both global and local characteristics in whole-body small animal MRI. The method integrates state-of-The-Art image processing techniques from the areas of articulated atlas-based segmentation, planar reformation, and SRR. A proof-of-concept is provided with two case studies involving CT, BLI, and MRI data of bone and kidney tumors in a mouse model. We show that local SRR-MRI is a computationally efficient complementary imaging modality for the precise characterization of tumor metastases, and that the method provides a feasible high-resolution alternative to conventional MRI.