Sequence Of Projection Images Research Articles

A Fast Image Guide Registration Supported by Single Direction Projected CBCT

Image registration is an important research topic in medical image-guided therapy, which is dedicated to registering the high-dose imaging sequences with low-dose/faster means. Registering computer tomography (CT) scanning sequences with cone beam computer tomography (CBCT) scanning sequences is a typical application and has been widely used in CBCT-guided radiotherapy. The main problem is the difference in image clarity of these two image sequences. To solve this problem, for the single projection image sequence matching tasks encountered in medical practice, a novel local quality based curved section encoding strategy is proposed in this paper, which is called the high-quality curved section (HQCS). As an optimized cross-section regularly encoded along the sequence of image, this curved section could be used in order to solve the matching problem. Referencing the independent ground truth provided by medical image physicians, with an experiment combined with the four most widely used indicators used on image registration, matching performance of HQCS on CT/CBCT datasets was tested with varying clarity. Experimental results show that the proposed HQCS can register the CT/CBCT effectively and outperforms the commonly used methods. Specifically, the proposed HQCS has low time complexity and higher scalability, which indicates that the application enhanced the task of diagnosis.

Motion reduction for quantitative brain sodium MR imaging with a navigated flexible twisted projection imaging sequence at 9.4 T.

Quantitative measurement of the tissue sodium concentration (TSC) provides a metric for tissue cell volume fraction for monitoring tumor responses to therapy and neurodegeneration in the brain as well as applications outside the central nervous system such as the fixed charge density in cartilage. Despite the low detection sensitivity of the sodium MR signal compared to the proton signal and the requirement for a long repetition time to minimize longitudinal magnetization saturation, acquisition time has been reduced to less than 10 min for a nominal isotropic voxel size of 3.3 mm with the improved acquisition efficiency of twisted projection imaging (TPI) at 9.4 T. However, patient motion can degrade the accuracy of the quantification even within these acquisition times. Our goal has been to improve the robustness of quantitative sodium MR imaging by minimizing the impact of motion that may occur even in cooperative patients. We present a method to spatially encode a lower resolution navigator echo after encoding the free induction decay signal for the quantitative image at no time penalty. Both the imaging and navigator data are sampled with flexTPI readout trajectories. Navigator images are generated at a higher temporal resolution (∼1 min) albeit at lower spatial resolution (8 mm) than the quantitative high-resolution images. The multiple volumes of navigator echo images are then aligned to extract the translational and rotational motion parameters assuming rigid-body motion. These parameters are used to align the k-space data during the acquisition of each volume of the quantitative images. Our results show significantly reduced image blurring with this method when the subject's head moved randomly by up to 7° between the navigator acquisitions.

Continuous MR bone density measurement using water- and fat-suppressed projection imaging (WASPI) for PET attenuation correction in PET-MR

Due to the lack of signal from solid bone in normal MR sequences for the purpose of MR-based attenuation correction, investigators have proposed using the ultrashort echo time (UTE) pulse sequence, which yields signal from bone. However, the UTE-based segmentation approach might not fully capture the intra- and inter-subject bone density variation, which will inevitably lead to bias in reconstructed PET images. In this work, we investigated using the water- and fat-suppressed proton projection imaging (WASPI) sequence to obtain accurate and continuous attenuation for bones. This approach is capable of accounting for intra- and inter-subject bone attenuation variations. Using data acquired from a phantom, we have found that that attenuation correction based on the WASPI sequence is more accurate and precise when compared to either conventional MR attenuation correction or UTE-based segmentation approaches.

A New Workflow for x-ray fluorescence tomography: MAPSToTomoPy.

X-ray fluorescence tomography involves the acquisition of a series of 2D x-ray fluorescence datasets between which a specimen is rotated. At the Advanced Photon Source at Argonne National Laboratory, the workflow at beamlines 2-ID-E and 21-ID-D (the Bionanoprobe, a cryogenic microscope system) has included the use of the program MAPS for obtaining elemental concentrations from 2D images, and the program TomoPy which was developed to include several tomographic reconstruction methods for x-ray transmission data. In the past, fluorescence projection images from an individual chemical element were hand-assembled into a 3D dataset for reconstruction using interactive tools such as ImageJ. We describe here the program MAPSToTomoPy, which provides a graphical user interface (GUI) to control a workflow between MAPS and TomoPy, with tools for visualizing the sinograms of projection image sequences from particular elements and to use these to help correct misalignments of the rotation axis. The program also provides an integrated output of the 3D distribution of the detected elements for subsequent 3D visualization packages.

The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis.

This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images. Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses. The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly suggestive for malignancy (BIRADS 5) indicating the required variety of shapes and margins of these models. The assessment of the BIRADS scores for all observers indicated good agreement based on Kendall's coefficient for both the 2D and the tomosynthesis evaluations. The paired analysis of the BIRADS scores between 2D and tomosynthesis for each observer revealed consistent behavior for the real and simulated masses. A database of 3D mass models, with variety of shapes and margins, was validated for the realism of their appearance for 2D digital mammography and for breast tomosynthesis. This database is suitable for use in future observer performance studies whether in virtual clinical trials or in patient images with simulated lesions.

Mapping tissue sodium concentration in the human brain: A comparison of MR sequences at 9.4 Tesla

Sodium is the second most abundant MR-active nucleus in the human body and is of fundamental importance for the function of cells. Previous studies have shown that many pathophysiological conditions induce an increase of the average tissue sodium concentration. To date, several MR sequences have been used to measure sodium. The aim of this study was to evaluate the performance and suitability of five different MR sequences for quantitative sodium imaging on a whole-body 9.4Tesla MR scanner. Numerical simulations, phantom experiments and in vivo imaging on healthy subjects were carried out. The results demonstrate that, of these five sequences, the Twisted Projection Imaging sequence is optimal for quantitative sodium imaging, as it combines a number of features which are particularly relevant in order to obtain high quality quantitative images of sodium. These include: ultra-short echo times, efficient k-space sampling, and robustness against off-resonance effects. Mapping of sodium in the human brain is a technique not yet fully explored in neuroscience. Ultra-high field sodium MRI may provide new insights into the pathogenesis of neurological disorders, and may help to develop new and disease-specific biomarkers for the early diagnosis and therapeutic intervention before irreversible brain damage has taken place.

The Sinogram Polygonizer for Reconstructing 3D Shapes

This paper proposes a novel approach, the sinogram polygonizer, for directly reconstructing 3D shapes from sinograms (i.e., the primary output from X-ray computed tomography (CT) scanners consisting of projection image sequences of an object shown from different viewing angles). To obtain a polygon mesh approximating the surface of a scanned object, a grid-based isosurface polygonizer, such as Marching Cubes, has been conventionally applied to the CT volume reconstructed from a sinogram. In contrast, the proposed method treats CT values as a continuous function and directly extracts a triangle mesh based on tetrahedral mesh deformation. This deformation involves quadratic error metric minimization and optimal Delaunay triangulation for the generation of accurate, high-quality meshes. Thanks to the analytical gradient estimation of CT values, sharp features are well approximated, even though the generated mesh is very coarse. Moreover, this approach eliminates aliasing artifacts on triangle meshes.

Tomographic reconstruction of heat release rate perturbations induced by helical modes in turbulent swirl flames

Swirl flows with vortex breakdown are widely used in industrial combustion systems for flame stabilization. This type of flow is known to sustain a hydrodynamic instability with a rotating helical structure, one common manifestation of it being the precessing vortex core. The role of this unsteady flow mode in combustion is not well understood, and its interaction with combustion instabilities and flame stabilization remains unclear. It is therefore important to assess the structure of the perturbation in the flame that is induced by this helical mode. Based on principles of tomographic reconstruction, a method is presented to determine the 3-D distribution of the heat release rate perturbation associated with the helical mode. Since this flow instability is rotating, a phase-resolved sequence of projection images of light emitted from the flame is identical to the Radon transform of the light intensity distribution in the combustor volume and thus can be used for tomographic reconstruction. This is achieved with one stationary camera only, a vast reduction in experimental and hardware requirements compared to a multi-camera setup or camera repositioning, which is typically required for tomographic reconstruction. Different approaches to extract the coherent part of the oscillation from the images are discussed. Two novel tomographic reconstruction algorithms specifically tailored to the structure of the heat release rate perturbations related to the helical mode are derived. The reconstruction techniques are first applied to an artificial field to illustrate the accuracy. High-speed imaging data acquired in a turbulent swirl-stabilized combustor setup with strong helical mode oscillations are then used to reconstruct the 3-D structure of the associated perturbation in the flame.

Respiratory Motion Compensation Using Diaphragm Tracking for Cone-Beam C-Arm CT: A Simulation and a Phantom Study

Long acquisition times lead to image artifacts in thoracic C-arm CT. Motion blur caused by respiratory motion leads to decreased image quality in many clinical applications. We introduce an image-based method to estimate and compensate respiratory motion in C-arm CT based on diaphragm motion. In order to estimate respiratory motion, we track the contour of the diaphragm in the projection image sequence. Using a motion corrected triangulation approach on the diaphragm vertex, we are able to estimate a motion signal. The estimated motion signal is used to compensate for respiratory motion in the target region, for example, heart or lungs. First, we evaluated our approach in a simulation study using XCAT. As ground truth data was available, a quantitative evaluation was performed. We observed an improvement of about 14% using the structural similarity index. In a real phantom study, using the artiCHEST phantom, we investigated the visibility of bronchial tubes in a porcine lung. Compared to an uncompensated scan, the visibility of bronchial structures is improved drastically. Preliminary results indicate that this kind of motion compensation can deliver a first step in reconstruction image quality improvement. Compared to ground truth data, image quality is still considerably reduced.

Coalescence analysis for evolving foamsviaoptical flow computation on projection image sequences

A novel image-processing procedure is proposed for the analysis of sequences of two-dimensional projection images. Sudden events like the merging of bubbles in an evolving foam can be detected and spatio-temporally located in a given projection image sequence. The procedure is based on optical flow computations extended by a forward-backward check for each time step. Compared with prior methods, efficient suppression of noise or false events is achieved owing to uniform foam motion, and the reliability of detection is thus increased. The applicability of the proposed procedure in combination with synchrotron radiography is illustrated by a series of characteristic studies of foams of different kind. First, the detection of single-bubble collapses in aqueous foams is considered. Second, a spatial distribution of coalescence events in metals foamed in casting molds is estimated. Finally, the structural stability of polymer foams containing admixed solid nanoparticles is examined.

The simulation of 3D microcalcification clusters in 2D digital mammography and breast tomosynthesis

This work proposes a new method of building 3D models of microcalcification clusters and describes the validation of their realistic appearance when simulated into 2D digital mammograms and into breast tomosynthesis images. A micro-CT unit was used to scan 23 breast biopsy specimens of microcalcification clusters with malignant and benign characteristics and their 3D reconstructed datasets were segmented to obtain 3D models of microcalcification clusters. These models were then adjusted for the x-ray spectrum used and for the system resolution and simulated into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. Six radiologists were asked to distinguish between 40 real and 40 simulated clusters of microcalcifications in two separate studies on 2D mammography and tomosynthesis datasets. Receiver operating characteristic (ROC) analysis was used to test the ability of each observer to distinguish between simulated and real microcalcification clusters. The kappa statistic was applied to assess how often the individual simulated and real microcalcification clusters had received similar scores ("agreement") on their realistic appearance in both modalities. This analysis was performed for all readers and for the real and the simulated group of microcalcification clusters separately. "Poor" agreement would reflect radiologists' confusion between simulated and real clusters, i.e., lesions not systematically evaluated in both modalities as either simulated or real, and would therefore be interpreted as a success of the present models. The area under the ROC curve, averaged over the observers, was 0.55 (95% confidence interval [0.44, 0.66]) for the 2D study, and 0.46 (95% confidence interval [0.29, 0.64]) for the tomosynthesis study, indicating no statistically significant difference between real and simulated lesions (p > 0.05). Agreement between allocated lesion scores for 2D mammography and those for the tomosynthesis series was poor. The realistic appearance of the 3D models of microcalcification clusters, whether malignant or benign clusters, was confirmed for 2D digital mammography images and the breast tomosynthesis datasets; this database of clusters is suitable for use in future observer performance studies related to the detectability of microcalcification clusters. Such studies include comparing 2D digital mammography to breast tomosynthesis and comparing different reconstruction algorithms.

Impact of gradient timing error on the tissue sodium concentration bioscale measured using flexible twisted projection imaging

The rapid biexponential transverse relaxation of the sodium MR signal from brain tissue requires efficient k-space sampling for quantitative imaging in a time that is acceptable for human subjects. The flexible twisted projection imaging (flexTPI) sequence has been shown to be suitable for quantitative sodium imaging with an ultra-short echo time to minimize signal loss. The fidelity of the k-space center location is affected by the readout gradient timing errors on the three physical axes, which is known to cause image distortion for projection-based acquisitions. This study investigated the impact of these timing errors on the voxel-wise accuracy of the tissue sodium concentration (TSC) bioscale measured with the flexTPI sequence. Our simulations show greater than 20% spatially varying quantification errors when the gradient timing errors are larger than 10 μs on all three axes. The quantification is more tolerant of gradient timing errors on the Z-axis. An existing method was used to measure the gradient timing errors with <1 μs error. The gradient timing error measurement is shown to be RF coil dependent, and timing error differences of up to ∼16 μs have been observed between different RF coils used on the same scanner. The measured timing errors can be corrected prospectively or retrospectively to obtain accurate TSC values.

Visualization and quantification of intestinal transit and motor function by real‐time tracking of19F labeled capsules in humans

A combined (19)F and (1)H MRI framework for the assessment of human intestinal transit and motor function is presented. This framework consists of silicone coated polychlorotrifluoroethylene capsules filled with perfluoro-[15]-crown-5-ether as (19)F marker, a flexible (19)F surface coil and a (19)F projection imaging sequence, allowing for real-time tracking of a single or multiple capsules. The capsules (length 11.5 mm, Ø 7.2 mm) contain 140 μL perfluoro-[15]-crown-5-ether and were tested for cytotoxicity and leakage prior to oral administration. A balanced SSFP projection sequence was implemented, yielding a temporal resolution of 133 ms. Optional multi-frequency excitation, allowing for interleaved tracking of differently labeled (19)F capsules, was incorporated. The passage of the (19)F capsules through intestinal sections was monitored in two healthy volunteers. Capsule coordinates were successfully coregistered with anatomical reference scans. Intestinal motility, residence times, lengths and forward velocities were determined. Simultaneous tracking of two capsules allowed for the assessment of peristaltic patterns with correction for respiratory motion. By providing the means for real-time multiple capsule tracking and high resolution anatomical imaging, the presented multinuclear imaging framework has the potential to provide important supplemental information for physiological and pharmaceutical research.

Quantitative sodium imaging with a flexible twisted projection pulse sequence

The quantification of sodium MR images from an arbitrary intensity scale into a bioscale fosters image interpretation in terms of the spatially resolved biochemical process of sodium ion homeostasis. A methodology for quantifying tissue sodium concentration using a flexible twisted projection imaging sequence is proposed that allows for optimization of tradeoffs between readout time, signal-to-noise ratio efficiency, and sensitivity to static field susceptibility artifacts. The gradient amplitude supported by the slew rate at each k-space radius regularizes the readout gradient waveform design to avoid slew rate violation. Static field inhomogeneity artifacts are corrected using a frequency-segmented conjugate phase reconstruction approach, with field maps obtained quickly from coregistered proton imaging. High-quality quantitative sodium images have been achieved in phantom and volunteer studies with real isotropic spatial resolution of 7.5 x 7.5 x 7.5 mm(3) for the slow T(2) component in approximately 8 min on a clinical 3-T scanner. After correcting for coil sensitivity inhomogeneity and water fraction, the tissue sodium concentration in gray matter and white matter was measured to be 36.6 +/- 0.6 micromol/g wet weight and 27.6 +/- 1.2 micromol/g wet weight, respectively.

Coalescence measurements for evolving foams monitored by real-time projection imaging

Real-time radiographic projection imaging together with novel spatio-temporal image analysis is presented to be a powerful technique for the quantitative analysis of coalescence processes accompanying the generation and temporal evolution of foams and emulsions. Coalescence events can be identified as discontinuities in a spatio-temporal image representing a sequence of projection images. Detection, identification of intensity and localization of the discontinuities exploit a violation criterion of the Fourier shift theorem and are based on recursive spatio-temporal image partitioning. The proposed method is suited for automated measurements of discontinuity rates (i.e., discontinuity intensity per unit time), so that large series of radiographs can be analyzed without user intervention. The application potential is demonstrated by the quantification of coalescence during the formation and decay of metal foams monitored by real-time x-ray radiography.

Rotational micro‐CT using a clinical C‐arm angiography gantry

Rotational angiography (RA) gantries are used routinely to acquire sequences of projection images of patients from which 3D renderings of vascular structures are generated using Feldkamp cone-beam reconstruction algorithms. However, these systems have limited resolution (<4 lp/mm). Micro-computed tomography (micro-CT) systems have better resolution (>10 lp/mm) but to date have relied either on rotating object imaging or small bore geometry for small animal imaging, and thus are not used for clinical imaging. The authors report here the development and use of a 3D rotational micro-angiography (RMA) system created by mounting a micro-angiographic fluoroscope (MAF) [35 microm pixel, resolution >10 microp/mm, field of view (FOV)=3.6 cm] on a standard clinical FPD-based RA gantry (Infinix, Model RTP12303J-G9E, Toshiba Medical Systems Corp., Tustin, CA). RA image sequences are obtained using the MAF and reconstructed. To eliminate artifacts due to image truncation, lower-dose (compared to MAF acquisition) full-FOV (FFOV) FPD RA sequences (194 microm pixel, FOV=20 cm) were also obtained to complete the missing data. The RA gantry was calibrated using a helical bead phantom. To ensure high-quality high-resolution reconstruction, the high-resolution images from the MAF were aligned spatially with the lower-dose FPD images, and the pixel values in the FPD image data were scaled to match those of the MAF. Images of a rabbit with a coronary stent placed in an artery in the Circle of Willis were obtained and reconstructed. The MAF images appear well aligned with the FPD images (average correlation coefficient before and after alignment: 0.65 and 0.97, respectively) Greater details without any visible truncation artifacts are seen in 3D RMA (MAF-FPD) images than in those of the FPD alone. The FWHM of line profiles of stent struts (100 microm diameter) are approximately 192+/-21 and 313+/-38 microm for the 3D RMA and FPD data, respectively. In addition, for the dual-acquisition 3D RMA, FFOV FPD data need not be of the highest quality, and thus may be acquired at lower dose compared to a standard FPD acquisition. These results indicate that this system could provide the basis for high resolution images of regions of interest in patients with a reduction in the integral dose compared to the standard FPD approach.

Generation of digital textured surface models from hologram recordings

Digital sensors and fast digital image processing facilitate the use of pulsed holography for 3D surface measurement of moving objects. The real image of a hologram is reconstructed optically. A sequence of high-resolution projection images of the real image with a varying distance to the hologram is recorded digitally. Focus detection in this image sequence by digital image processing yields the shape of the recorded object. The image intensity serves as a precise pixel-matching texture. An application of this concept is the generation of a textured 3D computer model of a facial surface from a portrait hologram.

3D flow estimation using perspective and parallel projection images

AbstractIn this paper, the authors discuss a new method for using perspective projection image sequences and approximated parallel projection images that were generated from multiple perspective projection images to estimate the 3D flow representing the motion of a target in space. This technique uses the light ray information of multiple perspective projection images that were obtained by capturing the target at different positions to generate the approximated parallel projection images, which are projected in multiple ray directions. It then uses the optical flows that were obtained from both the perspective projection image sequence and approximated parallel projection image sequence to estimate the 3D flow of the measurement points. A feature of this technique is that using images that were projected in the ray directions of the measurement points that were established on a perspective projection image as the approximated parallel projection image sequence corresponding to the perspective projection image sequence makes it unnecessary to map pixel values between stereo images, which is required in the conventional technique. The authors verified through experiments that this technique enabled approximated parallel projection images to be generated and the 3D flow to be estimated satisfactorily by using them. © 2006 Wiley Periodicals, Inc. Syst Comp Jpn, 37(7): 47–60, 2006; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/scj.20258

Computerized mass detection for digital breast tomosynthesis directly from the projection images

Digital breast tomosynthesis (DBT) has recently emerged as a new and promising three-dimensional modality in breast imaging. In DBT, the breast volume is reconstructed from 11 projection images, taken at source angles equally spaced over an arc of 50 degrees. Reconstruction algorithms for this modality are not fully optimized yet. Because computerized lesion detection in the reconstructed breast volume will be affected by the reconstruction technique, we are developing a novel mass detection algorithm that operates instead on the set of raw projection images. Mass detection is done in three stages. First, lesion candidates are obtained for each projection image separately, using a mass detection algorithm that was initially developed for screen-film mammography. Second, the locations of a lesion candidate are backprojected into the breast volume. In this feature volume, voxel intensities are a combined measure of detection frequency (e.g., the number of projections in which a given lesion candidate was detected), and a measure of the angular range over which a given lesion was detected. Third, features are extracted after reprojecting the three-dimensional (3-D) locations of lesion candidates into projection images. Features are combined using linear discriminant analysis. The database used to test the algorithm consisted of 21 mass cases (13 malignant, 8 benign) and 15 cases without mass lesions. Based on this database, the algorithm yielded a sensitivity of 90% at 1.5 false positives per breast volume. Algorithm performance is positively biased because this dataset was used for development, training, and testing, and because the number of algorithm parameters was approximately the same as the number.of patient cases. Our results indicate that computerized mass detection in the sequence of projection images for DBT may be effective despite the higher noise level in those images.

Dynamic small animal lung imaging via a postacquisition respiratory gating technique using micro-cone beam computed tomography 1

Micro computed tomography is an important tool for small animal imaging. On many occasions, it is desirable to image lungs in a live instead of postmortem small animal to perform a pulmonary physiology study. Because the lungs are moving, gating with respect to the ventilatory phase has to be performed to reduce motion artifacts. Precapture ventilation gating may be difficult to achieve in some situations, which motivates us to propose and implement a simple postacquisition gating method. Rats were used as the subjects in this study. A sequence of low-dose projection images were acquired at 30 frames per second for each view angle. During each capture sequence the rat undergoes multiple ventilation cycles. Using the sequence of projection images, an automated region of interest algorithm, based on integrated grayscale intensity, tracts the ventilatory phase of the lungs. In the processing of an image sequence, multiple projection images are identified at a particular phase and averaged to improve the signal-to-noise ratio. The resulting averaged projection images from different view angles are input to a Feldkamp cone-beam algorithm reconstruction algorithm to obtain isotropic image volumes. Reconstructions with reduced movement artifacts are obtained. In the gated reconstruction, registration of the bone is much better, the edge of the lung is clearly defined, and structures within the lung parenchyma are better resolved. Also, different phases of a breathing cycle can be reconstructed from one single tomographic scan by the proposed gating method. A postacquisition gating method using the phase information encoded in the 2-dimensional cone beam projections is proposed. This method is simple to implement and does not require additional experimental set-up to monitor the respiration. It may find applications in lung tumor detection, dynamic pulmonary physiology studies, and the respiratory systems modeling. Minimal motion artifact data sets improve qualitative and quantitative analysis techniques that are useful in physiologic studies of pulmonary structure and function.