Computation and visualization of three-dimensional soft tissue motion in the orbit

This work presents a method to measure the soft tissue motion in three dimensions in the orbit during gaze. It has been shown that two-dimensional (2-D) quantification of soft tissue motion in the orbit is effective in the study of orbital anatomy and motion disorders. However, soft tissue motion is a three-dimensional (3-D) phenomenon and part of the kinematics is lost in any 2-D measurement. Therefore, T1-weighted magnetic resonance (MR) imaging volume sequences are acquired during gaze and soft tissue motion is quantified using a generalization of the Lucas and Kanade optical flow algorithm to three dimensions. New techniques have been developed for visualizing the 3-D flow field as a series of color-texture mapped 2-D slices or as a combination of volume rendering for display of the anatomy and scintillation rendering for the display of the motion field. We have studied the performance of the algorithm on four-dimensional volume sequences of synthetic motion, simulated motion of a static object imaged by MR, an MR-imaged rotating object and MR-imaged motion in the human orbit during gaze. The accuracy of the analysis is sufficient to characterize motion in the orbit and scintillation rendering is an effective visualization technique for 3-D motion in the orbit.     

Automatic scoliosis detection based on local centroids evaluationon moire topographic images of human backs

This paper presents a technique for automating human scoliosis detection by computer based on moiré topographic images of human backs. Scoliosis is a serious disease often suffered by teenagers. For prevention, screening is performed at schools in Japan employing a moiré method in which doctors inspect moiré images of subjects' backs visually. The inspection of a large number of moiré images collected by the school screening causes exhaustion of doctors and leads to misjudgment. Computer-aided diagnosis of scoliosis has, therefore, been requested eagerly by orthopedists. To automate the inspection process, unlike existent three-dimensional techniques, displacement of local centroids is evaluated two-dimensionally between the left-hand side and the right-hand side of the moiré images in the present technique. The technique was applied to real moiré images to draw a distinction between normal and abnormal cases. According to the leave-out method, the entire 120 image data (60 normal and 60 abnormal) were separated into three data sets. The linear discriminant function based on Mahalanobis distance was defined on the two-dimensional feature space employing one of the data sets containing 40 moiré images and classified 80 images in the remaining two sets. The technique finally achieved the average classification rate of 88.3%.     

Retrospective shading correction based on entropy minimization

Shading is a prominent phenomenon in microscopy, manifesting itself via spurious intensity variations not present in the original scene. The elimination of shading effects is frequently necessary for subsequent image processing tasks, especially if quantitative analysis is the final goal. While most of the shading effects may be minimized by setting up the image acquisition conditions carefully and capturing additional calibration images, object-dependent shading calls for retrospective correction. In this paper a novel method for retrospective shading correction is proposed. Firstly, the image formation process and the corresponding shading effects are described by a linear image formation model, which consists of an additive and a multiplicative parametric component. Secondly, shading correction is performed by the inverse of the image formation model, whose shading components are estimated retrospectively by minimizing the entropy of the acquired images. A number of tests, performed on artificial and real microscopical images, show that this approach is efficient for a variety of differently structured images and as such may have applications in and beyond the field of microscopical imaging.     

Ridge-based vessel segmentation in color images of the retina

A method is presented for automated segmentation of vessels in two-dimensional color images of the retina. This method can be used in computer analyses of retinal images, e.g., in automated screening for diabetic retinopathy. The system is based on extraction of image ridges, which coincide approximately with vessel centerlines. The ridges are used to compose primitives in the form of line elements. With the line elements an image is partitioned into patches by assigning each image pixel to the closest line element. Every line element constitutes a local coordinate frame for its corresponding patch. For every pixel, feature vectors are computed that make use of properties of the patches and the line elements. The feature vectors are classified using a kappaNN-classifier and sequential forward feature selection. The algorithm was tested on a database consisting of 40 manually labeled images. The method achieves an area under the receiver operating characteristic curve of 0.952. The method is compared with two recently published rule-based methods of Hoover et al. and Jiang et al. The results show that our method is significantly better than the two rule-based methods (p < 0.01). The accuracy of our method is 0.944 versus 0.947 for a second observer.     

Three-dimensional guide-wire reconstruction from biplane image sequences for integrated display in 3-D vasculature

Using three-dimensional rotational X-ray angiography (3DRA), three-dimensional (3-D) information of the vasculature can be obtained prior to endovascular interventions. However, during interventions, the radiologist has to rely on fluoroscopy images to manipulate the guide wire. In order to take full advantage of the 3-D information from 3DRA data during endovascular interventions, a method is presented that yields an integrated display of the position of the guide wire and vasculature in 3-D. The method relies on an automated method that tracks the guide wire simultaneously in biplane fluoroscopy images. Based on the calibrated geometry of the C-arm, the 3-D guide-wire position is determined and visualized in the 3-D coordinate system of the vasculature. The method is evaluated in an intracranial anthropomorphic vascular phantom. The influence of the angle between projections, distortion correction of the projection images, and accuracy of geometry knowledge on the accuracy of 3-D guide-wire reconstruction from biplane images is determined. If the calibrated geometry information is used and the images are corrected for distortion, a mean distance to the reference standard of 0.42 mm and a tip distance of 0.65 mm is found, which means that accurate guide-wire reconstruction from biplane images can be performed.     

A General Framework for Geometry-Driven Evolution Equations

This paper presents a general framework to generate multi-scale representations of image data. The process is considered as an initial value problem with an acquired image as initial condition and a geometrical invariant as driving force of an evolutionary process. The geometrical invariants are extracted using the family of Gaussian derivative operators. These operators naturally deal with scale as a free parameter and solve the ill-posedness problem of differentiation. Stability requirements for numerical approximation of evolution schemes using Gaussian derivative operators are derived and establish an intuitive connection between the allowed time-step and scale. This approach has been used to generalize and implement a variety of nonlinear diffusion schemes. Results on test images and medical images are shown.     

Hierarchical interpolation

A state-of-the-art report on reversible (or lossless, or error-free, or information preserving) compression of medical images is presented. Reversible compression consists of two steps, decorrelation and coding. Methods for intraframe decorrelation of 2-D images can be divided into three classes, viz., transform, predictive, and multiresolution decorrelation. Transform decorrelation methods are found to be unsuited for reversible compression. DPCM (differential pulse code modulation) the best predictive method, decorrelates well, but has a great disadvantage in its parameter dependency     

Retrospective motion correction in digital subtraction angiography: a review

Digital subtraction angiography (DSA) is a well-established modality for the visualization of blood vessels in the human body. A serious disadvantage of this technique, inherent to the subtraction operation, is its sensitivity to patient motion. The resulting artifacts frequently reduce the diagnostic value of the images. Over the past two decades, many solutions to this problem have been put forward. In this paper, the authors give an overview of the possible types of motion artifacts and the techniques that have been proposed to avoid them. The main purpose of this paper is to provide a detailed review and discussion of retrospective motion correction techniques that have been described in the literature, to summarize the conclusions that can be drawn from these studies, and to provide suggestions for future research     

MR venography of the human brain using susceptibility weighted imaging at very high field strength

Objective We investigate the implications of high magnetic field strength on MR venography based on susceptibility-weighted imaging (SWI) and estimate the optimum echo time to obtain maximum contrast between blood and brain tissue. Materials and methods We measured tissue contrast and relaxation times at 7 T of gray matter, white matter, and venous blood in vivo. Results relaxation times of gray matter, white matter, and venous blood in vivo yielded 32.9 ± 2.3, 27.7 ± 4.3, and 7.4 ± 1.4 ms, respectively. Optimum TE was found to be 15 ms which is supported by theoretical considerations. Using this optimum TE, we acquired 3D high resolution datasets with a large volume coverage in a short measurement time that show very detailed microanatomical structures of the human brain such as intracortical veins and laminar cortical substructures. Conclusions By applying optimised vessel filters (vesselness filter and vessel enhancing diffusion) whole brain MR venograms can be obtained at 7 T with a significantly reduced measurement time compared to 3 T.     

Automatic detection of abnormalities in chest radiographs using local texture analysis

A fully automatic method is presented to detect abnormalities in frontal chest radiographs which are aggregated into an overall abnormality score. The method is aimed at finding abnormal signs of a diffuse textural nature, such as they are encountered in mass chest screening against tuberculosis (TB). The scheme starts with automatic segmentation of the lung fields, using active shape models. The segmentation is used to subdivide the lung fields into overlapping regions of various sizes. Texture features are extracted from each region, using the moments of responses to a multiscale filter bank. Additional "difference features" are obtained by subtracting feature vectors from corresponding regions in the left and right lung fields. A separate training set is constructed for each region. All regions are classified by voting among the k nearest neighbors, with leave-one-out. Next, the classification results of each region are combined, using a weighted multiplier in which regions with higher classification reliability weigh more heavily. This produces an abnormality score for each image. The method is evaluated on two databases. The first database was collected from a TB mass chest screening program, from which 147 images with textural abnormalities and 241 normal images were selected. Although this database contains many subtle abnormalities, the classification has a sensitivity of 0.86 at a specificity of 0.50 and an area under the receiver operating characteristic (ROC) curve of 0.820. The second database consist of 100 normal images and 100 abnormal images with interstitial disease. For this database, the results were a sensitivity of 0.97 at a specificity of 0.90 and an area under the ROC curve of 0.986.