Evaluation of dose recalculation vs dose deformation in a commercial platform for deformable image registration with a computational phantom

Deformable image registration (DIR) is an important tool for mapping both dose and contours of a new set of planning images when recurrent, that is, adaptive radiotherapy, or further treatment, that is, re-treatment, is required. The aim of this study was to evaluate the need for plan recalculation in deformed anatomies and to develop a reliable workflow for validating the DIR algorithm to be applied to dose warping such as dose accumulation (DA) for adaptive radiation therapy, and dose summation (DS) in the case of re-treatment. A set of 3 computational phantoms was developed to validate the application of B-Spline-based registrations for dose mapping among the various computed tomography image data sets. Two different versions were defined for each phantom to simulate clinical needs: adaptive radiotherapy and re-treatment; a DIR was performed to obtain a displacement vector field (DVF) for dose applications. Comparison of calculated and deformed doses was carried out by means of known markers inside the virtual phantoms. The differences were evaluated using a 3% dose index as acceptance criteria. A paired Wilcoxon signed-rank test was carried out to test the statistical significance of differences within the markers. Significant differences were only observed for the deformed dose in the DA test; no significant differences were observed for recalculated dose values in the DA and DS tests. The dose index is in accordance with these results. The warping dose process obtained by registering various sets of images was validated; a recalculation approach seems to be more accurate for DA purposes and image-guided adaptive radiation therapy (IGART) applications.     

Boundary-based warping of brain MR images

The goal of this work was to develop a warping technique for mapping a brain image to another image or atlas data, with minimum user interaction and independent of gray level information. We have developed and tested three different methods for warping magnetic resonance (MR) brain images. We utilize a deformable contour to extract and warp the boundaries of the two images. A mesh-grid coordinate system is constructed for each brain, by applying a distance transformation to the resulting contours, and scaling. In the first method (MGC), the first image is mapped to the second image based on a one-to-one mapping between different layers defined by the mesh-grid. In the second method (IDW), the corresponding pixels in the two images are found using the above mesh-grid system and a local inverse-distance weights interpolation. In the third proposed method (TSB), a subset of grid points is used for finding the parameters of a spline transformation, which defines the global warping. The warping methods were applied to clinical MR consisting of diffusion-weighted and T2-weighted images of the human brain. The IDW and TSB methods were superior in ranking of diagnostic quality of the warped MR images to the MGC (P < 0.01) as defined by a neuroradiologist. The deformable contour warping produced excellent diagnostic quality for the diffusion-weighted images coregistered and warped to T2 weighted images. J. Magn. Reson. Imaging 2000;12:417-429.     

Research on time-course differential imaging by bone scintigraphy

The temporal image subtraction technique was applied to bone scintigraphy, using Photoshop (commercially available image processing software) and Morpher (public domain warping software). For the temporal subtraction images, 81 subtraction images (19 cases) were prepared by a method used to subtract the previous images from the current ones. Registration of the current and previous images was performed by manual operation using Photoshop, and warping was done using the warping function of Morpher. In addition, difference images prepared after correcting the distributions of radioactive isotopes of the current and previous images using the count of the pelvic region were also examined. Compared with manual operation, alignment of images by warping improved registration and reduced the generation of pseudo-images of subtraction images. The rate of identification of abnormal accumulation-enhanced regions and subjective evaluation by doctors was improved for warping more than for manual operation. Furthermore, abnormal hot regions, which are difficult to find in film images, could be found in three subtraction images. In addition, it was confirmed that abnormal hot regions become more visible in many cases by preparing subtraction images after correcting the count between images using the count of the pelvic region. Thus, it is suggested that the temporal image subtraction technique in bone scintigraphy enables more accurate observation of enhancement of or changes in abnormal hot regions, which will support diagnostic reading. It is considered that enhancement of or changes in abnormal hot regions will be more accurately understood through further detailed discussion in the future.     

Three-dimensional reconstruction of fine vascularity in ultrasound breast imaging using contrast-enhanced spatial compounding: in vitro analyses

RATIONALE AND OBJECTIVES: Ultrasound image quality can be improved by imaging an object (here: the female breast) from different viewing angles in one image plane. With this technique, which is commonly referred to as spatial compounding, a more isotropic resolution is achieved while speckle noise and further artifacts are reduced. We present results obtained from a combination of spatial compounding with contrast-enhanced ultrasound imaging in three dimensions to reduce contrast specific artifacts (depth dependency, shadowing, speckle) and reconstruct vascular structures. MATERIALS AND METHODS: We used a conventional ultrasound scanner and a custom made mechanical system to rotate an ultrasound curved array probe around an object (360 degrees , 36 transducer positions). For 10 parallel image planes, ultrasound compound images were generated of a flow-mimicking phantom consecutively supplied with water and contrast agent. These compound images were combined to form a volume dataset and postprocessed to obtain a sonographic subtraction angiography. RESULTS: Image quality was significantly improved by spatial compounding for the native (ie, without contrast agent), and, in particular, for the contrast-enhanced case. After subtracting the native images from the contrast-enhanced ones, only structures supplied with contrast agent remain. This technique yields much better results for compound images than for conventional ultrasound images because speckle noise and an anisotropic resolution affect the latter. CONCLUSIONS: With the presented approach contrast specific artifacts can be eliminated efficiently, and a subtraction angiography can be computed. A speckle reduced three-dimensional reconstruction of submillimeter vessel structures was achieved for the first time. In the future, this technique can be applied in vivo to image the vascularity of cancer in the female breast.     

Geometric distortion correction of high-resolution 3 T diffusion tensor brain images.

Diffusion-weighted images based on echo planar sequences suffer from distortions due to field inhomogeneities from susceptibility differences as well as from eddy currents arising from diffusion gradients. In this paper, a novel approach using nonlinear warping based on optic flow to correct distortions of baseline and diffusion weighted echo planar images (EPI) acquired at 3 T is presented. The distortion correction was estimated by warping the echo planar images to the anatomically correct T2-weighted fast spin echo images (T2-FSE). A global histogram intensity matching of the T2-FSE precedes the base line EPI image distortion correction. A local intensity-matching algorithm was used to transform labeled T2-FSE regions to match intensities of diffusion-weighted EPI images prior to distortion correction of these images. Evaluation was performed using three methods: (i) visual comparison of overlaid contours, (ii) a global mutual information index, and (iii) a local distance measure between homologous points. Visual assessment and the global index demonstrated a decrease in geometrical distortion and the distance measure showed that distortions are reduced to a subvoxel level. In conclusion, the warping algorithm is effective in reducing geometric distortions, enabling generation of anatomically correct diffusion tensor images at 3 T.     

New method of evaluating edge-preserving adaptive filters for computed tomography (CT): digital phantom method

To evaluate the characteristics of edge-preserving adaptive filters for selectively eliminating noise without affecting resolution in low-dose scanning, we have developed a digital phantom image and evaluated noise statistical values, noise characteristics, and resolution characteristics. The results confirmed that edge-preserving adaptive filters function as smoothing filters in low-contrast regions containing noise, permitting the density resolution to be improved, while the strength of the smoothing filter is reduced to maintain spatial resolution in high-contrast regions containing small structures. It has therefore been confirmed that edge-preserving adaptive filters function as filters for selectively eliminating only the noise elements that are increased when the exposure dose is reduced and that such filters are effective for improving image quality. Using such digital phantom images, images acquired using conditions that are difficult to set in actual CT scanning can be obtained and images specifically for the evaluation target can easily be generated. In addition, the noise level, frequency distribution of the noise, and resolution characteristics of the objects present in the input image can be freely set. It is concluded that evaluation of processing using a digital phantom image is effective for evaluating image processing.     

Semiautomatic nonrigid registration for the prostate and pelvic MR volumes

RATIONALE AND OBJECTIVES: Three-dimensional (3D) nonrigid image registration for potential applications in prostate cancer treatment and interventional magnetic resonance (iMRI) imaging-guided therapies were investigated. MATERIALS AND METHODS: An almost fully automated 3D nonrigid registration algorithm using mutual information and a thin plate spline (TPS) transformation for MR images of the prostate and pelvis were created and evaluated. In the first step, an automatic rigid body registration with special features was used to capture the global transformation. In the second step, local feature points (FPs) were registered using mutual information. An operator entered only five FPs located at the prostate center, left and right hip joints, and left and right distal femurs. The program automatically determined and optimized other FPs at the external pelvic skin surface and along the femurs. More than 600 control points were used to establish a TPS transformation for deformation of the pelvic region and prostate. Ten volume pairs were acquired from three volunteers in the diagnostic (supine) and treatment positions (supine with legs raised). RESULTS: Various visualization techniques showed that warping rectified the significant pelvic misalignment by the rigid-body method. Gray-value measures of registration quality, including mutual information, correlation coefficient, and intensity difference, all improved with warping. The distance between prostate 3D centroids was 0.7 +/- 0.2 mm after warping compared with 4.9 +/- 3.4 mm with rigid-body registration. CONCLUSION: Semiautomatic nonrigid registration works better than rigid-body registration when patient position is changed greatly between acquisitions. It could be a useful tool for many applications in the management of prostate.     

Retrospective geometric correlation of MR, CT, and PET images

Magnetic resonance imaging, computed tomographic, and positron emission tomographic studies of the brain provide complementary information, and many patients undergo more than one of these studies during the course of their diagnostic workup and treatment. A new technique for quantitative geometric correlation of such studies makes it possible to create integrated multimodality images by mapping features from one image onto an image obtained with another modality. The coordinate transformation between any pair of images is found by a semiautomatic algorithm for matching models of the patient's external surface as depicted in the two data sets. The resultant hybrid images, which combine complementary features of different studies, are often more useful for diagnosis and treatment planning than are the original single-modality images. The algorithm can also be used for spatial registration of baseline studies with follow-up images created with the same modality, which allows tracking of a lesion to detect subtle interval changes in size and shape. This technique can be applied to images acquired in routine clinical practice, since it is completely retrospective and does not necessitate special positioning or landmarking of the patient.     

Three-dimensional method for comparing in vivo interventional MR images of thermally ablated tissue with tissue response

PURPOSE: To investigate the ability of magnetic resonance (MR) to monitor radio-frequency (RF) ablation treatments by comparing MR images of thermal lesions to histologically assayed cellular damage. We developed a new methodology using three-dimensional registration for making spatial correlations. MATERIALS AND METHODS: A low-field, open MRI system was used to guide an ablation probe into rabbit thigh muscle and acquire MR volumes after ablation. After fixation, we sliced and photographed the tissue at 3-mm intervals, using a specially designed apparatus, to obtain a volume of tissue images. Histologic samples were digitized using a video microscopy system. For our three-dimensional registration method, we used the tissue images as the reference, and registered histology and MR images to them using two different computer alignment steps. First, the MR volume was aligned to the volume of tissue images by registering needle fiducials placed near the tissue of interest. Second, we registered the histology images with the tissue images using a two-dimensional warping technique that aligned internal features and the outside boundary of histology and tissue images. RESULTS: The MR and histology images were very well aligned, and registration accuracy, determined from displacement of needle fiducials, was 1.32 +/- 0.39 mm (mean +/- SD), which compared favorably to the MR voxel dimensions (0.70 mm in-plane and 3.0 mm thick). A preliminary comparison of MR and tissue response showed that the region inside the elliptical hyperintense rim in MR closely corresponds to the region of necrosis as established by histology, with a mean absolute distance between MR and histology boundaries of 1.17 mm, slightly smaller than the mean registration error. The MR region slightly overestimated the region of necrosis, with a mean signed distance between boundaries of 0.85 mm. CONCLUSION: Our results suggest that our methodology can be used to achieve three-dimensional registration of histology and in vivo MR images. In MR lesion images, the inner border of the hyperintense region corresponds to the border of irreversible cell damage. This is good evidence that during RF ablation treatments, iMRI lesion images can be used for real-time feedback.     

Reconstruction from truncated projections in CT using adaptive detruncation

If the object exceeds the field of measurement (FOM) of a given CT scanner, severe artifacts may result. In this work, we propose an adaptive detruncation (ADT) method to reconstruct images from medical CT projections which are truncated in the transaxial direction. The truncated projections are extrapolated by estimating the convex hull of the patient. The ADT method allows us not only to achieve artifact-free images in the FOM but also to extend the images beyond the FOM, and can therefore be very attractive, for example, in PET/CT scanners for attenuation correction.     

Toward atlas-assisted automatic interpretation of MRI morphological brain scans in the presence of tumor

RATIONALE AND OBJECTIVES: Determination of distorted brain anatomy surrounding a tumor causing a mass effect is much more difficult than interpretation of normal brain scans, particularly because this distortion is not easily predictable a tumor may be located in any place and vary substantially in size, shape, and radiological appearance. The objective of our work is to provide a qualitative means for rapid estimation of brain anatomy distorted by tumor. MATERIALS AND METHODS: Toward achieving this objective, we use an electronic and deformable brain atlas of gross anatomy along with a fast atlas-to-data warping technique. The deformed atlas determines the distorted anatomy surrounding a tumor and can be used for structure labeling (naming). The warping algorithm uses the Talairach transformation followed by three-dimensional nonlinear tumor deformation based on a geometric assumption that the tumor, delineated on radiological images, compresses its surrounding tissues radially. RESULTS: The approach is implemented and a dedicated application is developed. It processes the data automatically in five steps: (1) load data, (2) set the Talairach landmarks and perform the Talairach transformation, (3) segment the tumor, (4) warp the scan nonlinearly in three dimensions, and (5) explore the scan. The approach is very fast, and a magnetic resonance imaging scan is processed in 10-15 seconds on a standard personal computer. It is fully automatic and gives the user control over the entire process. CONCLUSION: Despite its limitations in modeling and validation, this practical solution provides a rapid and potentially useful qualitative assessment of anatomy deformed by a mass effect tumor.     

Multishot diffusion-weighted FSE using PROPELLER MRI.

A method for obtaining diffusion-weighted images that are free from the artifacts associated with echo-planar acquisitions, such as signal pile-up and geometric warping, is introduced. It uses an ungated, multishot fast spin-echo (FSE) acquisition that is self-navigated. The phase of the refocusing pulses is alternated to minimize non-Carr-Purcell-Meiboom-Gill (CPMG) artifacts. Several reconstruction methods are combined to make this method robust against motion artifacts. Examples are shown of clinical diffusion-weighted imaging and high-resolution diffusion tensor imaging.     

Adaptive keyhole methods for dynamic magnetic resonance image reconstruction.

Dynamic magnetic resonance imaging (MRI) acquires a sequence of images for the visualization of the temporal variation of tissue or organs. Keyhole methods such as Fourier keyhole (FK) and keyhole SVD (KSVD) are the most popular methods for image reconstruction in dynamic MRI. This paper provides a class of adaptive keyhole methods, called adaptive FK (AFK) and adaptive KSVD (AKSVD), for dynamic MRI reconstruction. The proposed methods are based on the conventional Fourier encoding and SVD encoding schemes. Instead of the conventional keyhole methods' duplication of un-acquired data from the reference images, the proposed methods use a temporal model to depict the inter-frame dynamic changes and to estimate the un-acquired data in each successive frame. Because the model is online identified from the acquired data, the proposed methods do not require the pre-imaging process, the navigator signals, and any prior knowledge of the imaged objects. Furthermore, the new methods use the conventional keyhole encoding schemes without the bias to any particular object characters, and the temporal model for updating information is in the general form of AR process without the preference to any particular motion types. Hence, the proposed methods are designed as a generic approach to dynamic MRI, other than for any specific class of objects. Studies on dynamic MRI data set show that the new methods can produce images with much lower reconstruction error than the conventional FK and KSVD.     

Correcting spatial distortion in histological images.

We described an interactive method for correcting spatial distortion in histology samples, applied them to a large set of image data, and quantitatively evaluated the quality of the corrections. We demonstrated registration of histology samples to photographs of macroscopic tissue samples and to MR images. We first described methods for obtaining corresponding fiducial and anatomical points, including a new technique for determining boundary correspondence points. We then describe experimental methods for tissue preparation, including a technique for adding color-coded internal and boundary ink marks that are used to validate the method by measuring the registration error. We applied four different transformations with internal and boundary correspondence points, and measured the distance error between other internal ink fiducials. A large number of boundary points, typically 20-30, and at least two internal points were required for accurate warping registration. Interior errors with the transformation methods were ordered: thin plate spline (TPS) approximately non-warping相似文献   

"Motion-frozen" display and quantification of myocardial perfusion.

Gated myocadial perfusion SPECT (MPS) incorporates functional and perfusion information of the left ventricle (LV). To improve the image quality and accuracy of gated MPS we propose to eliminate the influence of cardiac LV motion in the display and quantification by a novel "motion-frozen" (MF) technique. METHODS: Three-dimensional LV contours were identified on images of the individual time phases. Three-dimensional phase-to-phase motion vectors were derived by sampling of the epi- and endocardial surfaces. A nonlinear image warping (thin-plate spline) was applied to warp all image phases to fit the end-diastolic (ED) phase. Warped images were created to provide the LV image in the ED phase but containing counts from an arbitrary number of time intervals. MF quantification has been performed using the same phase-to-phase motion vectors. MF normal perfusion limits were created from (99m)Tc sestamibi gated MPS studies of 40 females and 40 males with low likelihood (<5%) of coronary artery disease. All MF processing was completely automated. In the initial evaluation, we assessed the display quality and quantification of stress images using MF processing in 51 consecutive patients with 16-frame electrocardiographic gating and available coronary angiography. RESULTS: The display quality was significantly better for MF images as assessed visually. The MF images had the appearance of ED frames but were less noisy and of higher resolution than the summed images. MF images had higher maximum count values in the LV (116% +/- 6%) and higher contrast (12.5 +/- 7.7 vs. 9.5 +/- 3.2) than the corresponding summed images. The area under the receiver operator characteristic curve for prediction of stenoses > or = 70% by the MF method was 0.92 +/- 0.04 versus 0.89 +/- 0.04 by standard quantification (P = not significant). The computation time for automated MF quantification and warping was <25 s for each case. CONCLUSION: We have developed a novel technique for display and quantification of gated myocardial perfusion images, which retrospectively eliminates blur due to cardiac motion. Such processing of gated MPS appears to improve the effective resolution of images. Initial evaluation indicates that it may improve the accuracy of gated MPS in detection of coronary artery disease.     

Ultrasound motion estimation using a hierarchical feature weighting algorithm.

The quality of ultrasonic images is usually influenced by speckle noises and the temporal decorrelation of the speckle patterns. Most traditional motion estimation algorithms are not suitable for speckle tracking in medical ultrasonic images which usually have a low signal-to-noise ratio (SNR). This paper proposes a new motion estimation algorithm that is designed for assessing the dense velocity fields of soft tissue motion in a sequence of ultrasonic B-mode images. We design a hierarchical maximum a posteriori estimator together with an adaptive feature weighted mechanism to estimate the motion field from an ultrasonic image sequence. The proposed method was compared with several existing motion estimation methods via a series of experiments with synthetic speckle image sequences. Performance was also tested on in vivo ultrasonic images. The experimental results show that motion can be assessed with better accuracy than other methods for synthetic speckle images and a good correspondence with clinicians' observations has also been achieved for clinical ultrasonic images.     

Correction for geometric distortion and N/2 ghosting in EPI by phase labeling for additional coordinate encoding (PLACE).

Echo-planar imaging (EPI) is vulnerable to geometric distortion and N/2 ghosting. These artifacts can be analyzed with an intuitive k-t space tool, and here we propose a simple method for their correction. In a slightly modified additional EPI acquisition, we sample the k-t space with a shift in k(y) by adding a small area to the phase-encoding (PE) gradient. Physically, the added gradient area creates a relative phase ramp across the object and directly encodes the undistorted original y-coordinate of each voxel into a phase difference between two distorted complex images, in a method called "phase labeling for additional coordinate encoding" (PLACE). The phase information is then used to map the mismapped signals back to their original locations for geometric and intensity correction. Smoothing of expanded complex data matrix effectively reduces noise in the differential phase map and allows subpixel warping. The two acquired images can also be averaged to effectively suppress the N/2 ghost. Efficient correction for both artifacts can be achieved with three acquisitions. These acquisitions can also serve as reference scans to correct for geometric distortion and/or N/2 ghost artifacts on all images in a time series. The technique was successfully demonstrated in phantom and animal studies.     

Magnetic field changes in the human brain due to swallowing or speaking

Variations in the magnetic field in the human brain caused by the processes of swallowing or speaking are measured. In both processes, motion of the pharyngeal muscles, especially the tongue and jaw, alter the susceptibility-induced magnetic field distribution at the brain slice being imaged. This leads to image warping, compromising the analysis of a time series of images, such as in functional magnetic resonance imaging (fMRI). These dynamic changes are assessed by acquiring a time series of images using a gradient-echo asymmetric-spin-echo sequence (GREASE), a technique in which two images are acquired for each excitation–one during the gradient echo, and one during the latter part of the spin echo. The NMR phase difference between the two images is a measure of the magnetic field distribution. A series of brain images, acquired with this sequence while the subject either swallows or speaks, indicated negative magnetic field changes of up to O.087 ppm in the inferior region of the brain for both speaking and swallowing, and in some speech, additional positive changes of up to O.056 ppm in the frontal region of the brain were indicated.     

MR imaging of periosteal and cortical changes of bone

The changes seen in the periosteum and cortical bone are fundamental radiographic features of bone disease. The basic radiographic findings used for diagnosis of bone lesions (patterns of cortical destruction and of periosteal new bone formation) can be well identified with magnetic resonance (MR) imaging. The authors used comparative radiographic, computed tomographic, and MR images to illustrate patterns of periosteal reaction (simple, laminated, spiculated, Codman triangle), geographic and permeative cortical destruction, cortical erosion, cortical expansion and continuity, and intraosseous and extraosseous calcification. The only feature not well demonstrated by MR imaging is pattern or extent of soft-tissue calcification. Although MR images are not required for diagnosis of most peripheral bone lesions, when they are obtained, these fundamental diagnostic features should not be ignored.     

Modern gray-scale sonography of the breast

Recently, the diagnosis of breast diseases by ultrasound has changed radically. It is no longer a complementary modality to mammography but a separate method to investigate breast disease. Innovative high-resolution ultrasound allows more specific diagnosis of breast tumors. Tissue-harmonic imaging not only uses the transmitted, fundamental frequency to obtain an image but also the harmonic frequency. The harmonic signal is processed by the ultrasound system with the result of better delineation of tissue structures and spatial compounding assembles an image from multiple images taken from different angles of echo waves. The effect is the reduction of artifacts with optimized contrast. Finally the advanced speckle reduction technique is used to smooth and homogenize the image. Additionally continuous advancement of new high-resolution linear transducers is responsible for the essential improvement of image quality. In conclusion, it is recommended to integrate all of the described modalities in order to obtain diagnostically conclusive image quality. This article demonstrates the new techniques and applications exemplified using images.