An objective comparison of 3-D image interpolation methods

Specialised CT scan protocols for 3-D pre-operative planning of total hip replacement

X-ray computer tomography (CT) provides an accurate source of information in orthopaedics. Many computer aided orthopaedic surgery systems are based on CT images; thus, obtaining high resolution images is important. However, this may result in an excessive radiation dose for the patient. This study is aimed at developing a special CT scanning protocol for the hip region that can be adopted in clinical practise for 3-D pre-operative planning of total hip replacement surgery. Optimisation of CT acquisition parameters is investigated for both axial and spiral CT and the two resulting protocols are compared in terms of effective radiation dose to the patient. Results show that spiral CT with D=3 and P=1.5 in regions with higher morphological and density gradients and D=5 and P=1.5 in regions where the morphology is more regular degrades the image quality slightly but allows acquisition of a higher number of images at comparable costs, increasing the longitudinal resolution of the acquired data set. The effective dose is comparable to that of a standard pelvic CT exam. Adjusting the axial CT scan parameters the effective dose can be reduced, however lowering the accuracy of 3-D bone geometry reconstruction.

Application of an adaptive control grid interpolation technique to morphological vascular reconstruction

The problem of interslice magnetic resonance (MR) image reconstruction arises in a broad range of medical applications. In such cases, there is a need to approximate information present in the original subject that is not reflected in contiguously acquired MR images because of hardware sampling limitations. In the context of vascular morphology reconstruction, this information is required in order for subsequent visualization and computational analysis of blood vessels to be most effective. Toward that end we have developed a method of vascular morphology reconstruction based on adaptive control grid interpolation (ACGI) to function as a precursor to visualization and computational analysis. ACGI has previously been implemented in addressing various problems including video coding and tracking. This paper focuses on the novel application of the technique to medical image processing. ACGI combines features of optical flow-based and block-based motion estimation algorithms to enhance insufficiently dense MR data sets accurately with a minimal degree of computational complexity. The resulting enhanced data sets describe vascular geometries. These reconstructions can then be used as visualization tools and in conjunction with computational fluid dynamics (CFD) simulations to offer the pressure and velocity information necessary to quantify power loss. The proposed ACGI methodology is envisioned ultimately to play a role in surgical planning aimed at producing optimal vascular configurations for successful surgical outcomes.

Binary morphological shape-based interpolation applied to 3-D tooth reconstruction

In this paper, we propose an interpolation algorithm using a mathematical morphology morphing approach. The aim of this algorithm is to reconstruct the n-dimensional object from a group of (n - 1)-dimensional sets representing sections of that object. The morphing transformation modifies pairs of consecutive sets such that they approach in shape and size. The interpolated set is achieved when the two consecutive sets are made idempotent by the morphing transformation. We prove the convergence of the morphological morphing. The entire object is modeled by successively interpolating a certain number of intermediary sets between each two consecutive given sets. We apply the interpolation algorithm for three-dimensional tooth reconstruction.

Quantitative evaluation of convolution-based methods for medical image interpolation

Interpolation is required in a variety of medical image processing applications. Although many interpolation techniques are known from the literature, evaluations of these techniques for the specific task of applying geometrical transformations to medical images are still lacking. In this paper we present such an evaluation. We consider convolution-based interpolation methods and rigid transformations (rotations and translations). A large number of sinc-approximating kernels are evaluated, including piecewise polynomial kernels and a large number of windowed sinc kernels, with spatial supports ranging from two to ten grid intervals. In the evaluation we use images from a wide variety of medical image modalities. The results show that spline interpolation is to be preferred over all other methods, both for its accuracy and its relatively low computational cost.

Serial automated three-dimensional intravascular ultrasound analysis of the self-expanding Radius stent

Automated 3-dimensional intravascular ultrasound (IVUS) analysis was used to assess status of the treated coronary artery immediately and 6 months after placement of a self-expanding Radius stent in 15 patients. Serial 3-dimensional IVUS analysis demonstrated gradual stent expansion that countered neointimal proliferation and preserved the lumen.

Survey: interpolation methods in medical image processing

Image interpolation techniques often are required in medical imaging for image generation (e.g., discrete back projection for inverse Radon transform) and processing such as compression or resampling. Since the ideal interpolation function spatially is unlimited, several interpolation kernels of finite size have been introduced. This paper compares 1) truncated and windowed sinc; 2) nearest neighbor; 3) linear; 4) quadratic; 5) cubic B-spline; 6) cubic; g) Lagrange; and 7) Gaussian interpolation and approximation techniques with kernel sizes from 1 x 1 up to 8 x 8. The comparison is done by: 1) spatial and Fourier analyses; 2) computational complexity as well as runtime evaluations; and 3) qualitative and quantitative interpolation error determinations for particular interpolation tasks which were taken from common situations in medical image processing. For local and Fourier analyses, a standardized notation is introduced and fundamental properties of interpolators are derived. Successful methods should be direct current (DC)-constant and interpolators rather than DC-inconstant or approximators. Each method's parameters are tuned with respect to those properties. This results in three novel kernels, which are introduced in this paper and proven to be within the best choices for medical image interpolation: the 6 x 6 Blackman-Harris windowed sinc interpolator, and the C2-continuous cubic kernels with N = 6 and N = 8 supporting points. For quantitative error evaluations, a set of 50 direct digital X rays was used. They have been selected arbitrarily from clinical routine. In general, large kernel sizes were found to be superior to small interpolation masks. Except for truncated sinc interpolators, all kernels with N = 6 or larger sizes perform significantly better than N = 2 or N = 3 point methods (p << 0.005). However, the differences within the group of large-sized kernels were not significant. Summarizing the results, the cubic 6 x 6 interpolator with continuous second derivatives, as defined in (24), can be recommended for most common interpolation tasks. It appears to be the fastest six-point kernel to implement computationally. It provides eminent local and Fourier properties, is easy to implement, and has only small errors. The same characteristics apply to B-spline interpolation, but the 6 x 6 cubic avoids the intrinsic border effects produced by the B-spline technique. However, the goal of this study was not to determine an overall best method, but to present a comprehensive catalogue of methods in a uniform terminology, to define general properties and requirements of local techniques, and to enable the reader to select that method which is optimal for his specific application in medical imaging.

Subunity coordinate translation with Fourier transform to achieve efficient and quality three-dimensional medical image interpolation

A new approach to the interpolation of three-dimensional (3D) medical images is presented. Instead of going through the conventional interpolation scheme where the continuous function is first reconstructed from the discrete data set and then resampled, the interpolation is achieved with a subunity coordinate translation technique. The original image is first transformed into the spatial-frequency domain. The phase of the transform is then modified with n-1 linear phase terms in the axial direction to achieve n-1 subunity coordinate translations with a distance 1/n, where n is an interpolation ratio, following the phase shift theorem of Fourier transformation. All the translated images after inverse Fourier transformation are then interspersed in turn into the original image. Since windowing plays an important role in the process, different window functions have been studied and a proper recommendation is provided. The interpolation quality produced with the present method is as good as that with the sampling (sinc) function, while the efficiency, thanks to the fast Fourier transformation, is very much improved. The approach has been validated with both computed tomography (CT) and magnetic resonance (MR) images. The interpolations of 3D CT and MR images are demonstrated.

A task-specific evaluation of three-dimensional image interpolation techniques

Image interpolation is an important operation that is widely used in medical imaging, image processing, and computer graphics. A variety of interpolation methods are available in the literature. However, their systematic evaluation is lacking. In a previous paper, we presented a framework for the task-independent comparison of interpolation methods based on certain image-derived figures of merit using a variety of medical image data pertaining to different parts of the human body taken from different modalities. In this work, we present an objective task-specific framework for evaluating interpolation techniques. The task considered is how the interpolation methods influence the accuracy of quantification of the total volume of lesions in the brain of multiple sclerosis (MS) patients. Sixty lesion-detection experiments coming from ten patient studies, two subsampling techniques and the original data, and three interpolation methods are carried out, along with a statistical analysis of the results.