Reconstruction of multi-animal PET acquisitions with anisotropically variant PSF.

Among other factors such as random, attenuation and scatter corrections, uniform spatial resolution is key to performing accurate quantitative studies in Positron emission tomography (PET). Particularly in preclinical PET studies involving simultaneous acquisition of multiple animals, the degradation of image resolution due to the depth of interaction (DOI) effect far from the center of the Field of View (FOV) becomes a significant concern. In this work, we incorporated a spatially-variant resolution model into a real time iterative reconstruction code to obtain accurate images of multi-animal acquisition. We estimated the spatially variant point spread function (SV-PSF) across the FOV using measurements and Monte Carlo (MC) simulations. The SV-PSF obtained was implemented in a GPU-based Ordered subset expectation maximization (OSEM) reconstruction code, which includes scatter, attenuation and random corrections. The method was evaluated with acquisitions from two preclinical PET/CT scanners of the SEDECAL Argus family: a Derenzo phantom placed 2 cm off center in the 4R-SuperArgus, and a multi-animal study with 4 mice in the 6R-SuperArgus. The SV-PSF reconstructions showed uniform spatial resolution without significant increase in reconstruction time, with superior image quality compared to the uniform PSF model.

Fluorescence diffuse optical tomography using the split Bregman method.

PURPOSE: Standard image reconstruction methods for fluorescence Diffuse Optical Tomography (fDOT) generally make use of L2-regularization. A better choice is to replace the L2 by a total variation functional that effectively removes noise while preserving edges. Among the wide range of approaches available, the recently appeared Split Bregman method has been shown to be optimal and efficient. Furthermore, additional constraints can be easily included. We propose the use of the Split Bregman method to solve the image reconstruction problem for fDOT with a nonnegativity constraint that imposes the reconstructed concentration of fluorophore to be positive. METHODS: The proposed method is tested with simulated and experimental data, and results are compared with those yielded by an equivalent unconstrained optimization approach based on Gauss-Newton (GN) method, in which the negative part of the solution is projected to zero after each iteration. In addition, the method dependence on the parameters that weigh data fidelity and nonnegativity constraints is analyzed. RESULTS: Split Bregman yielded a reduction of the solution error norm and a better full width at tenth maximum for simulated data, and higher signal-to-noise ratio for experimental data. It is also shown that it led to an optimum solution independently of the data fidelity parameter, as long as the number of iterations is properly selected, and that there is a linear relation between the number of iterations and the inverse of the data fidelity parameter. CONCLUSIONS: Split Bregman allows the addition of a nonnegativity constraint leading to improve image quality.

Detection of visual activation in the rat brain using 2-deoxy-2-[(18)F]fluoro-D: -glucose and statistical parametric mapping (SPM).

PURPOSE: This study was designed to assess changes in brain glucose metabolism in rats after visual stimulation. MATERIALS AND METHODS: We sought to determine whether visual activation in the rat brain could be detected using a small-animal positron emission tomography (PET) scanner and 2-deoxy-2-[(18)F]fluoro-D: -glucose (FDG). Eleven rats were divided into two groups: (a) five animals exposed to ambient light and (b) six animals stimulated by stroboscopic light (10 Hz) with one eye covered. Rats were injected with FDG and, after 45 min of visual stimulation, were sacrificed and scanned for 90 min in a dedicated PET tomograph. Images were reconstructed by a three-dimensional ordered subset expectation maximization algorithm (1.8 mm full width at half maximum). A region-of-interest (ROI) analysis was performed on 14 brain structures drawn on coronal sections. Statistical parametric mapping (SPM) adapted for small animals was also carried out. Additionally, the brains of three rats were sliced into 20-microm sections for autoradiography. RESULTS: Analysis of ROI data revealed significant differences between groups in the right superior colliculus, right thalamus, and brainstem (p < or = 0.05). SPM detected the same areas as the ROI approach. Autoradiographs confirmed the existence of hyperactivation in the left superior colliculus and auditory cortex. CONCLUSIONS: To our knowledge, this is the first report that uses FDG-PET and SPM analysis to show changes in rat brain glucose metabolism after a visual stimulus.

Sinogram bow-tie filtering in FBP PET reconstruction.

Low-pass filtering of sinograms in the radial direction is the most common practice to limit noise amplification in filtered back projection (FBP) reconstruction of positron emission tomography studies. Other filtering strategies have been proposed to prevent the loss in resolution due to low-pass radial filters, although results have been diverse. Using the well-known properties of the Fourier transform of a sinogram, the authors defined a binary mask that matches the expected shape of the support region in the Fourier domain of the sinogram ("bow tie"). This mask was smoothed by a convolution with a ten-point Gaussian kernel which not only avoids ringing but also introduces a pre-emphasis at low frequencies. A new filtering scheme for FBP is proposed, comprising this smoothed bow-tie filter combined with a standard radial filter and an axial filter. The authors compared the performance of the bow-tie filtering scheme with that of other previously reported methods: Standard radial filtering, angular filtering, and stackgram-domain filtering. All the quantitative data in the comparisons refer to a baseline reconstruction using a ramp filter only. When using the smallest size of the Gaussian kernel in the stackgram domain, the authors achieved a noise reduction of 33% at the cost of degrading radial and tangential resolutions (14.5% and 16%, respectively, for cubic interpolation). To reduce the noise by 30%, the angular filter produced a larger degradation of contrast (3%) and tangential resolution (46% at 10 mm from the center of the field of view) and showed noticeable artifacts in the form of circular blurring dependent on the distance to the center of the field of view. For a similar noise reduction (33%), the proposed bow-tie filtering scheme yielded optimum results in resolution (gain in radial resolution of 10%) and contrast (1% increase) when compared with any of the other filters alone. Experiments with rodent images showed noticeable image quality enhancement when using the proposed bow-tie filtering scheme.

PeneloPET, a Monte Carlo PET simulation tool based on PENELOPE: features and validation.

Monte Carlo simulations play an important role in positron emission tomography (PET) imaging, as an essential tool for the research and development of new scanners and for advanced image reconstruction. PeneloPET, a PET-dedicated Monte Carlo tool, is presented and validated in this work. PeneloPET is based on PENELOPE, a Monte Carlo code for the simulation of the transport in matter of electrons, positrons and photons, with energies from a few hundred eV to 1 GeV. PENELOPE is robust, fast and very accurate, but it may be unfriendly to people not acquainted with the FORTRAN programming language. PeneloPET is an easy-to-use application which allows comprehensive simulations of PET systems within PENELOPE. Complex and realistic simulations can be set by modifying a few simple input text files. Different levels of output data are available for analysis, from sinogram and lines-of-response (LORs) histogramming to fully detailed list mode. These data can be further exploited with the preferred programming language, including ROOT. PeneloPET simulates PET systems based on crystal array blocks coupled to photodetectors and allows the user to define radioactive sources, detectors, shielding and other parts of the scanner. The acquisition chain is simulated in high level detail; for instance, the electronic processing can include pile-up rejection mechanisms and time stamping of events, if desired. This paper describes PeneloPET and shows the results of extensive validations and comparisons of simulations against real measurements from commercial acquisition systems. PeneloPET is being extensively employed to improve the image quality of commercial PET systems and for the development of new ones.

Design and performance evaluation of a coplanar multimodality scanner for rodent imaging.

This work reports on the development and performance evaluation of the VrPET/CT, a new multimodality scanner with coplanar geometry for in vivo rodent imaging. The scanner design is based on a partial-ring PET system and a small-animal CT assembled on a rotatory gantry without axial displacement between the geometric centers of both fields of view (FOV). We report on the PET system performance based on the NEMA NU-4 protocol; the performance characteristics of the CT component are not included herein. The accuracy of inter-modality alignment and the imaging capability of the whole system are also evaluated on phantom and animal studies. Tangential spatial resolution of PET images ranged between 1.56 mm at the center of the FOV and 2.46 at a radial offset of 3.5 cm. The radial resolution varies from 1.48 mm to 1.88 mm, and the axial resolution from 2.34 mm to 3.38 mm for the same positions. The energy resolution was 16.5% on average for the entire system. The absolute coincidence sensitivity is 2.2% for a 100-700 keV energy window with a 3.8 ns coincident window. The scatter fraction values for the same settings were 11.45% for a mouse-sized phantom and 23.26% for a rat-sized phantom. The peak noise equivalent count rates were also evaluated for those phantoms obtaining 70.8 kcps at 0.66 MBq/cc and 31.5 kcps at 0.11 MBq/cc, respectively. The accuracy of inter-modality alignment is below half the PET resolution, and the image quality of biological specimens agrees with measured performance parameters. The assessment presented in this study shows that the VrPET/CT system is a good performance small-animal imager, while the cost derived from a partial ring detection system is substantially reduced as compared with a full-ring PET tomograph.

Computed Tomography-Based Biomarker for Longitudinal Assessment of Disease Burden in Pulmonary Tuberculosis.

PURPOSE: Computed tomography (CT) images enable capturing specific manifestations of tuberculosis (TB) that are undetectable using common diagnostic tests, which suffer from limited specificity. In this study, we aimed to automatically quantify the burden of Mycobacterium tuberculosis (Mtb) using biomarkers extracted from x-ray CT images. PROCEDURES: Nine macaques were aerosol-infected with Mtb and treated with various antibiotic cocktails. Chest CT scans were acquired in all animals at specific times independently of disease progression. First, a fully automatic segmentation of the healthy lungs from the acquired chest CT volumes was performed and air-like structures were extracted. Next, unsegmented pulmonary regions corresponding to damaged parenchymal tissue and TB lesions were included. CT biomarkers were extracted by classification of the probability distribution of the intensity of the segmented images into three tissue types: (1) Healthy tissue, parenchyma free from infection; (2) soft diseased tissue, and (3) hard diseased tissue. The probability distribution of tissue intensities was assumed to follow a Gaussian mixture model. The thresholds identifying each region were automatically computed using an expectation-maximization algorithm. RESULTS: The estimated longitudinal course of TB infection shows that subjects that have followed the same antibiotic treatment present a similar response (relative change in the diseased volume) with respect to baseline. More interestingly, the correlation between the diseased volume (soft tissue + hard tissue), which was manually delineated by an expert, and the automatically extracted volume with the proposed method was very strong (R2 ≈ 0.8). CONCLUSIONS: We present a methodology that is suitable for automatic extraction of a radiological biomarker from CT images for TB disease burden. The method could be used to describe the longitudinal evolution of Mtb infection in a clinical trial devoted to the design of new drugs.

Extraction of the respiratory signal from small-animal CT projections for a retrospective gating method.

We propose a retrospective respiratory gating algorithm to generate dynamic CT studies. To this end, we compared three different methods of extracting the respiratory signal from the projections of small-animal cone-beam computed tomography (CBCT) scanners. Given a set of frames acquired from a certain axial angle, subtraction of their average image from each individual frame produces a set of difference images. Pixels in these images have positive or negative values (according to the respiratory phase) in those areas where there is lung movement. The respiratory signals were extracted by analysing the shape of the histogram of these difference images: we calculated the first four central and non-central moments. However, only odd-order moments produced the desired breathing signal, as the even-order moments lacked information about the phase. Each of these curves was compared to a reference signal recorded by means of a pneumatic pillow. Given the similar correlation coefficients yielded by all of them, we selected the mean to implement our retrospective protocol. Respiratory phase bins were separated, reconstructed independently and included in a dynamic sequence, suitable for cine playback. We validated our method in five adult rat studies by comparing profiles drawn across the diaphragm dome, with and without retrospective respiratory gating. Results showed a sharper transition in the gated reconstruction, with an average slope improvement of 60.7%.

FIRST: Fast Iterative Reconstruction Software for (PET) tomography.

Small animal PET scanners require high spatial resolution and good sensitivity. To reconstruct high-resolution images in 3D-PET, iterative methods, such as OSEM, are superior to analytical reconstruction algorithms, although their high computational cost is still a serious drawback. The higher performance of modern computers could make iterative image reconstruction fast enough to be viable, provided we are able to deal with the large number of probability coefficients for the system response matrix in high-resolution PET scanners, which is a difficult task that prevents the algorithms from reaching peak computing performance. Considering all possible axial and in-plane symmetries, as well as certain quasi-symmetries, we have been able to reduce the memory requirements to store the system response matrix (SRM) well below 1 GB, which allows us to keep the whole response matrix of the system inside RAM of ordinary industry-standard computers, so that the reconstruction algorithm can achieve near peak performance. The elements of the SRM are stored as cubic spline profiles and matched to voxel size during reconstruction. In this way, the advantages of 'on-the-fly' calculation and of fully stored SRM are combined. The on-the-fly part of the calculation (matching the profile functions to voxel size) of the SRM accounts for 10-30% of the reconstruction time, depending on the number of voxels chosen. We tested our approach with real data from a commercial small animal PET scanner. The results (image quality and reconstruction time) show that the proposed technique is a feasible solution.

[Technical limitations of the positron emission tomography (PET) for small laboratory animals]. / Limitaciones tecnológicas de la tomografía por emisión de positrones (PET) para pequeños animales de laboratorio.

Visualization and quantification of organ function by PET in small laboratory animals is becoming an outstanding tool for characterization of phenotype of transgenic and knock-out animals, for the study of animal models of human diseases, and for the development of new therapeutic drugs and diagnostic biochemical probes. To be able to make use of the PET with small laboratory animals in the same way as it is operated with humans it is necessary to account for the volumetric scale factor as well as for the requirement of maintaining the counting statistics. This work sketches the problems that these requirements represent for the technical design of the scanners and for the execution of the experiments. Finally, some characteristics of commercially available scanners (microPET/ FOCUS, HiDAC, eXplore VISTA, MOSAIC, YAP-(S)PET and rPET) are briefly discussed.

Improved quantification for local regions of interest in preclinical PET imaging.

In Positron Emission Tomography, there are several causes of quantitative inaccuracy, such as partial volume or spillover effects. The impact of these effects is greater when using radionuclides that have a large positron range, e.g. (68)Ga and (124)I, which have been increasingly used in the clinic. We have implemented and evaluated a local projection algorithm (LPA), originally evaluated for SPECT, to compensate for both partial-volume and spillover effects in PET. This method is based on the use of a high-resolution CT or MR image, co-registered with a PET image, which permits a high-resolution segmentation of a few tissues within a volume of interest (VOI) centered on a region within which tissue-activity values need to be estimated. The additional boundary information is used to obtain improved activity estimates for each tissue within the VOI, by solving a simple inversion problem. We implemented this algorithm for the preclinical Argus PET/CT scanner and assessed its performance using the radionuclides (18)F, (68)Ga and (124)I. We also evaluated and compared the results obtained when it was applied during the iterative reconstruction, as well as after the reconstruction as a postprocessing procedure. In addition, we studied how LPA can help to reduce the 'spillover contamination', which causes inaccurate quantification of lesions in the immediate neighborhood of large, 'hot' sources. Quantification was significantly improved by using LPA, which provided more accurate ratios of lesion-to-background activity concentration for hot and cold regions. For (18)F, the contrast was improved from 3.0 to 4.0 in hot lesions (when the true ratio was 4.0) and from 0.16 to 0.06 in cold lesions (true ratio = 0.0), when using the LPA postprocessing. Furthermore, activity values estimated within the VOI using LPA during reconstruction were slightly more accurate than those obtained by post-processing, while also visually improving the image contrast and uniformity within the VOI.

Development and validation of an open source quantification tool for DSC-MRI studies.

MOTIVATION: This work presents the development of an open source tool for the quantification of dynamic susceptibility-weighted contrast-enhanced (DSC) perfusion studies. The development of this tool is motivated by the lack of open source tools implemented on open platforms to allow external developers to implement their own quantification methods easily and without the need of paying for a development license. MATERIALS AND METHODS: This quantification tool was developed as a plugin for the ImageJ image analysis platform using the Java programming language. A modular approach was used in the implementation of the components, in such a way that the addition of new methods can be done without breaking any of the existing functionalities. For the validation process, images from seven patients with brain tumors were acquired and quantified with the presented tool and with a widely used clinical software package. The resulting perfusion parameters were then compared. RESULTS: Perfusion parameters and the corresponding parametric images were obtained. When no gamma-fitting is used, an excellent agreement with the tool used as a gold-standard was obtained (R(2)>0.8 and values are within 95% CI limits in Bland-Altman plots). CONCLUSION: An open source tool that performs quantification of perfusion studies using magnetic resonance imaging has been developed and validated using a clinical software package. It works as an ImageJ plugin and the source code has been published with an open source license.

Synthesis, structure, and pharmacological evaluation of the stereoisomers of furnidipine.

The synthesis and pharmacological activities of the four stereoisomers of methyl tetrahydrofuran-2-ylmethyl 2,6-dimethyl-4-(2'-nitrophenyl)-1,4-dihydropyridine-3,5- dicarboxylate(furnidipine) are reported. The four isomers were synthesized by a modified Hantzsch synthesis by reaction of (-)- or (+)-tetrahydrofuran-2-ylmethyl 3-aminocrotonate and methyl 2-[(2'-nitrophenyl)methylene]acetoacetate or, alternatively, by reaction of (-)- or (+)-tetrahydrofuran-2-ylmethyl 2-[(2'-nitrophenyl)methylene]acetoacetate and methyl 3-aminocrotonate. The 1:1 diastereomeric mixtures thus obtained were separated by chromatography, using poly(D-phenylglycine) as the chiral stationary phase. The enantiomeric purity of the stereoisomers was determined by a high-performance liquid chromatography-chiral stationary phase technique (HPLC-CSP). Attempts to obtain crystals of a single stereoisomer failed in different solvents, while methanol crystallization of the product obtained from (+/-)-tetrahydrofuran-2-ylmethyl 2-[(2'-nitrophenyl)methylene]acetoacetate and methyl 3-aminocrotonate yielded good-quality crystals of the most insoluble racemate which proved to be a mixture of the (SS)/(RR) enantiomers by X-ray crystallography. Conformational analysis of the stereoisomers, assuming rotation of the aryl substituent and ester groups, shows small energy differences (about 4 kcal.mol-1) between the most and the least favorable conformations. Binding studies were performed using [3H]isradipine as a reference ligand. The results showed stereospecificity of the furnidipine isomers in brain, ileum, and cardiac tissues, the (SS)- and (SR)-isomers clearly being more potent than their (RR)- and (RS)-enantiomers. The (SS)- and (SR)-isomers were also more selective on cerebral tissue when compared with ileal and cardiac preparations.

Reaction of 2-bromomethylazoles and TosMIC: A domino process to azolopyrimidines. Synthesis of core tricycle of the variolins alkaloids.

A new reaction of N-protected 2-bromomethylazoles and tosylmethyl isocyanide (TosMIC) leading to the preparation of azolopyrimidines is described. This domino sequence was used to synthesize the pyrido[3',2':4,5]pyrrolo[1,2-c]pyrimidine core of alkaloids variolins from 4-methoxy-2-methylpyrrolo[2,3-b]pyrimidine in two steps.

Automatic scoring of sister chromatid exchanges by image analysis in a dose response experiment.

A system which automatically selects second division metaphases and then, automatically scores the number of SCEs of each cell is described. In an initial set of experiments, the performance of the components of the system was measured using a data set in which metaphases had been visually classified as either 2nd division or other; and in 2nd division metaphases, every SCE had been marked on a hard copy. SCE scoring had a true positive rate of about 75% and a false positive rate of about 1.5 false SCEs per metaphase analyzed. Second division detection had a true positive rate of 80% and a false positive rate of about 10% of the non-2nd division cells. Next, the overall system was compared to human visual scoring in a dose-response experiment by analyzing the effect of mitomycin C on human chromosomes scored visually by two observers and by the fully automatic scoring. Human visual scoring and machine analysis showed similar dose responses, but the variability between them was considerable.

Experiment assessment of mass effects in the rat: implications for small animal PET imaging.

In vivo imaging using positron emission tomography (PET) is important in the development of new radiopharmaceuticals in rodent animal models for use as biochemical probes, diagnostic agents, or in drug development. We have shown mathematically that, if small animal imaging studies in rodents are to have the same "quality" as human PET studies, the same number of coincidence events must be detected from a typical rodent imaging "voxel" as from the human imaging voxel. To achieve this using the same specific activity preparation, we show that roughly the same total amount of radiopharmaceutical must be given to a rodent as to a human subject. At high specific activities, the mass associated with human doses, when administered to a rodent, may not decrease the uptake of radioactivity at non saturable sites or sites where an enzyme has a high capacity for a substrate. However, in the case of binding sites of low density such as receptors, the increased mass injected could saturate the receptor and lead to physiologic effects and non-linear kinetics. Because of the importance of the mass injected for small animal PET imaging, we experimentally compared high and low mass preparations using ex vivo biodistribution and phosphorimaging of three compounds: 2-fluoro-2-deoxyglucose (FDG), 6-fluoro-L-metatyrosine (FMT) and one receptor-directed compound, the serotonin 5HT1A receptor ligand, trans-4-fluoro-N-[2-[4-(2-methoxylphenyl) piperazino]ethyl]-N-(2-pyridyl) cyclohexane- carboxamide (FCWAY). Changes in the mass injected per rat did not affect the distribution of FDG, FMT, and FCWAY in the range of 0.6-1.9 nmol per rat. Changes in the target to nontarget ratio were observed for injected masses of FCWAY in the range of approximately 5-50 nmol per rat. If the specific activity of such compounds and/or the sensitivity of small animal scanners are not increased relative to human studies, small animal PET imaging will not correctly portray the "true" tracer distribution. These difficulties will only be exacerbated in animals smaller than the rat, e.g., mice.

Monitoring the correction of glycogen storage disease type 1a in a mouse model using [(18)F]FDG and a dedicated animal scanner.

Monitoring gene therapy of glycogen storage disease type 1a in a mouse model was achieved using [(18)F]FDG and a dedicated animal scanner. The G6Pase knockout (KO) mice were compared to the same mice after infusion with a recombinant adenovirus containing the murine G6Pase gene (Ad-mG6Pase). Serial images of the same mouse before and after therapy were obtained and compared with wild-type (WT) mice of the same strain to determine the uptake and retention of [(18)F]FDG in the liver. Image data were acquired from heart, blood pool and liver for twenty minutes after injection of [(18)F]FDG. The retention of [(18)F]FDG was lower for the WT mice compared to the KO mice. The mice treated with adenovirus-mediated gene therapy had retention similar to that found in age-matched WT mice. These studies show that FDG can be used to monitor the G6Pase concentration in liver of WT mice as compared to G6Pase KO mice. In these mice, gene therapy returned the liver function to that found in age matched WT controls as measured by the FDG kinetics in the liver compared to that found in age matched wild type controls.

Performance characteristics of a compact position-sensitive LSO detector module.

We assembled a compact detector module comprised of an array of small, individual crystals of lutetium oxyorthosilicate:Ce (LSO) coupled directly to a miniature, metal-can, position-sensitive photomultiplier tube (PSPMT). We exposed this module to sources of 511-keV annihilation radiation and beams of 30- and 140-keV photons and measured spatial linearity; spatial variations in module gain, energy resolution, and event positioning; coincidence timing; the accuracy and sensitivity of identifying the crystal-of-first-interaction at 511 keV; and the effects of intercrystal scatter and LSO background radioactivity. The results suggest that this scintillator/phototube combination should be highly effective in the coincidence mode and can be used, with some limitations, to image relatively low-energy single photon emitters. Photons that are completely absorbed on their first interaction at 511 keV are positioned by the module at the center of a crystal. Intercrystal scatter events, even those that lead to total absorption of the incident photon, are placed by the module in a regular "connect-the-dot" pattern that joins crystal centers. As a result, the accuracy of event positioning can be made to exceed 90%, though at significantly reduced sensitivity, by retaining only events that occur within small regions-of-interest around each crystal center and rejecting events that occur outside these regions in the connect-the-dot pattern.

Probe efficiency improvement with remote and transmission line tuning and matching.

In certain nuclear magnetic resonance (NMR) applications in which accessibility to the gantry is limited, performing optimal tuning and matching represents a major problem. Here, we discuss a method of tuning NMR probe circuits and matching their impedances which uses cables with different impedance values. This simple but efficient method may be advantageous compared with much more difficult perfect tuning and matching.

MRI visualization of small structures using improved surface coils.

In this paper we present the spatial resolution enhancement and noise reduction level achieved with an optimized inductively coupled surface coil specifically designed for our experiments. The technique of designing and implementing customized coils for magnetic resonance imaging of very small structures is described. We have designed a low cost prototype of an inductively coupled circular surface coil, tuned for 1H magnetic resonance imaging at 200 MHz. The coil is mounted on a customized teflon support. The inductive coupling used in this coil improves the signal-to-noise ratio by reducing various loss mechanisms (specially the dielectric losses). Test images have been acquired to determine the evolution of induced articular lesions in a rabbit animal model, as well as brain tumors in rats. The images show high spatial resolution, excellent B1 field homogeneity and no "hot spots". Comparing these images with those acquired with conventional coils, one finds better spatial resolution and signal-to-noise ratio, as well as larger field of view with less intense illumination artifact. The methodology can be used in any application that requires high quality imaging of small structures.