Biomedical subjects

Ruola Ning

Publications and source records attributed to Ruola Ning.

13 recordsLinked to original sources

PubMed2006

Volume fusion for two-circular-orbit cone-beam tomography.

By using the Feldkamp-Davis-Kress (FDK) algorithm, we can efficiently produce a digital volume, called the FDK volume, from cone-beam data acquired along a circular scan orbit. Due to the insufficiency of the cone-beam data set, the FDK volume suffers from nonuniform reproduction exactness. Specifically, the midplane (on the scan-orbit plane) can be exactly reproduced, and the reproduction exactness of off-midplanes decreases as the distance from the midplane increases. We describe the longitudinal falling-off degradation by a hatlike function and the spatial distribution over the object domain by an exactness volume. With two orthogonal circular scan orbits, we can reconstruct two FDK volumes and generate two exactness volumes. We propose a volume fusion scheme to combine the two FDK volumes into a single volume. Let Va and Vb denote the two FDK volumes, let Ea and Eb denote the exactness volumes for orbits Gamma(a) and Gamma(b), respectively, then the volume fusion is defined by Vab=VaWa+VbWb, with Wa=Ea/(Ea+Eb) and Wb=1-Wa. In the result, the overall reproduction exactness of Vab is expected to outperform that of Va, or Vb, or (Va+Vb)/2. In principle, this volume-fusion scheme is applicable for general cone-beam tomography with multiple nonorthogonal and noncircular orbits.

Algorithms↗

PubMed2006

Early detection and measurement of urothelial tumors in mice.

OBJECTIVES: To establish reliable noninvasive in vivo methods to detect, measure, and monitor experimentally induced urothelial tumors in mice. METHODS: UPII-SV40T transgenic mice reliably develop bladder tumors by expression of simian virus 40 large T antigen specifically in bladder urothelium through the use of the uroplakin II promoter. Two wild-type and 10 UPII-SV40T transgenic mice were monitored for microhematuria two to three times weekly using dipstick analysis. A unique flat panel detector-based cone beam computed tomography (FPD-CBCT) system imaged the urinary tracts of anesthetized mice after tail vein injection of an iodinated contrast agent (Omnipaque) that is excreted in urine. Within 10 seconds, the FPD-CBCT system acquired 290 two-dimensional images, which produced three-dimensional volumes with true isotropic resolution (180 microm)3 using a filtered back projection-based modified Feldkamp reconstruction algorithm. Amira, version 3.1.1-1, for MacOSX was used for data analysis and advanced visualization of the three-dimensional reconstructed FPD-CBCT images. RESULTS: Hematuria was present in UPII-SV40T transgenic mice at 32 days of age; the wild-type animals exhibited no hematuria. Filling defects, associated with histologically confirmed tumors, in the bladders of the UPII-SV40T transgenic mice were visualized in the reconstructed FPD-CBCT images 1 to 45 minutes after contrast agent injection. Longitudinal FPD-CBCT imaging sessions showed the tumor position, volume, and growth. CONCLUSIONS: The combination of early detection of hematuria and high-resolution in vivo FPD-CBCT imaging of murine bladder tumors enabled accurate longitudinal assessment of tumor growth and progression in individual animals. This approach could provide an important alternative to serial sacrifice experimental designs, while refining statistical power and reducing animal use.

Animals↗

PubMed2005

Supergridded cone-beam reconstruction and its application to point-spread function calculation.

In cone-beam computed tomography (CBCT), the volumetric reconstruction may in principle assume an arbitrarily fine grid. The supergridded cone-beam reconstruction refers to reconstructing the object domain or a subvolume thereof with a grid that is finer than the proper computed tomography sampling grid (as determined by gantry geometry and detector discreteness). This technique can naturally reduce the voxelization effect, thereby retaining more details for object reproduction. The grid refinement is usually limited to two or three refinement levels because the detail pursuit is eventually limited by the detector discreteness. The volume reconstruction is usually targeted to a local volume of interest due to the cubic growth in a three-dimensional (3D) array size. As an application, we used this technique for 3D point-spread function (PSF) measurement of a CBCT system by reconstructing edge spread profiles in a refined grid. Through an experiment with a Teflon ball on a CBCT system, we demonstrated the supergridded volume reconstruction (based on a Feldcamp algorithm) and the CBCT PSF measurement (based on an edge-blurring technique). In comparison with a postreconstruction image refinement technique (upsampling and interpolation), the supergridded reconstruction could produce better PSFs (in terms of a smaller FWHM and PSF fitting error).

Algorithms↗

PubMed2005

Forest representation of vessels in cone-beam computed tomographic angiography.

Cone-beam computed tomographic angiography (CBCTA) provides a fast three-dimensional (3D) vascular imaging modality, aiming at digitally representing the spatial vascular structure in an angiographic volume. Due to the finite coverage of cone-beam scan, as well as the volume cropping in volumetric image processing, an angiographic volume may fail to contain a whole vascular tree, but rather consist of a multitude of vessel segments or subtrees. As such, it is convenient to represent multitudinal components by a forest. The vessel tracking issue then becomes component characterization/identification in the forest. The forest representation brings several conveniences for vessel tracking: (1) to sort and count the vessels in an angiographic volume, for example, according to spatial occupancy and skeleton pathlength; (2) to single out a vessel and perform in situ 3D measurement and 3D visualization in the support space; (3) to delineate individual vessels from the original angiographic volume; and (4) to cull the forest by getting rid of non-vessels and small vessels. A 3D skeletonization is used to generate component skeletons. For tree construction from skeletons, we suggest a pathlength-based procedure, which lifts the restrictions of unit-width skeleton and root determination. We experimentally demonstrate the forest representation of a dog's carotid arteries in a CBCTA system. In principle, the forest representation is useful for managing vessels in both 2D angiographic images and 3D angiographic volumes.

Angiography↗

PubMed2005

Image denoising based on wavelets and multifractals for singularity detection.

This paper presents a very efficient algorithm for image denoising based on wavelets and multifractals for singularity detection. A challenge of image denoising is how to preserve the edges of an image when reducing noise. By modeling the intensity surface of a noisy image as statistically self-similar multifractal processes and taking advantage of the multiresolution analysis with wavelet transform to exploit the local statistical self-similarity at different scales, the pointwise singularity strength value characterizing the local singularity at each scale was calculated. By thresholding the singularity strength, wavelet coefficients at each scale were classified into two categories: the edge-related and regular wavelet coefficients and the irregular coefficients. The irregular coefficients were denoised using an approximate minimum mean-squared error (MMSE) estimation method, while the edge-related and regular wavelet coefficients were smoothed using the fuzzy weighted mean (FWM) filter aiming at preserving the edges and details when reducing noise. Furthermore, to make the FWM-based filtering more efficient for noise reduction at the lowest decomposition level, the MMSE-based filtering was performed as the first pass of denoising followed by performing the FWM-based filtering. Experimental results demonstrated that this algorithm could achieve both good visual quality and high PSNR for the denoised images.

Algorithms↗

PubMed2004

Three-dimensional point spread function measurement of cone-beam computed tomography system by iterative edge-blurring algorithm.

With separability assumed, we decompose a three-dimensional point spread function (3D PSF) into two-dimensional (2D) PSFs and further into one-dimensional (ID) PSFs. Based on the observation of the location invariance of a step edge under convolution, we propose a rectification procedure to automatically establish the step-edge function from a blurred edge profile. The ID PSF is modelled as a single-parameter Gaussian function, which is determined by iteratively blurring a step-edge function into a spread edge profile. A plastic solid ball (diameter approximately 6 mm) is used to provide double-edged rectangular functions along scanlines passing through the ball centre, and correspondingly, the reconstructed digital volume provides the blurred rectangular profiles. Experimenting with a cone-beam computed tomography system, we demonstrate the iterative edge-blurring algorithm for PSF measurement. By repositioning the ball phantom in the object support space, we measure the system's spatial variance in terms of full-width-at-half-maximum (FWHM) of the local PSFs. Specifically, we obtained the FWHMs for three specific locations at (0, 0, -40 mm), (0, 0, 0) and (0, 0, 40 mm), which are given by 0.92 +/- 0.10 mm, 0.65 +/- 0.08 mm and 0.93 +/- 0.10 mm, respectively.

Algorithms↗

PubMed2004

Breast volume denoising and noise characterization by 3D wavelet transform.

Breast imaging through cone-beam computed tomography provides a digital breast volume, with which the three-dimensional (3D) breast tissues can be analyzed. Data denoising, as a preprocessing step for subsequent volumetric breast segmentation is always needed. In this paper, we report a volumetric denoising technique by a separable 3D wavelet transform (WT), i.e. a '2D WT plus 1D WT' scheme. Specifically, the scheme performs two-dimensional (2D) wavelet denoising on a stack of slice images of the breast volume, followed by one-dimensional (1D) wavelet denoising along the stacking direction. The denoising is achieved by wavelet decomposition, high-pass subband attenuation, and wavelet synthesis. A one-level 3D WT produces eight subbands occupying the octants of the 3D wavelet space. Multilevel WT also provides a multiresolution representation of breast volume, i.e. a sequence of low-pass subbands. In general, most noise and irregularity features are imparted into the high-pass subbands, which are removed or reduced for denoising purpose. Meanwhile, the information in a subband can be characterized in terms of energy, variance, and entropy. Through 3D visualization, the spatial structure in a subband can also be visually perceived. Experimental demonstration with the breast volume reconstructed from a specimen is provided.

Mammography↗

PubMed2004

Image denoising based on multiscale singularity detection for cone beam CT breast imaging.

It was recently reported that the real-time flat panel detector-based cone-beam computed tomography (CBCT) breast imaging can help improve the detectability of small breast tumors with an X-ray dose comparable to that of the conventional mammography. In this paper, an efficient denoising algorithm is proposed to further reduce the X-ray exposure level required by a CBCT scan to acquire acceptable image quality. The proposed wavelet-based denoising algorithm possesses three significant characteristics: 1) wavelet coefficients at each scale are classified into two categories: irregular coefficients, and edge-related and regular coefficients; 2) noise in irregular coefficients is reduced as much as possible without producing artifacts to the denoised images; and 3) for the edge-related and regular coefficients, if they are at the first decomposition level, they are further denoised, otherwise, no modifications are made to them so as to obtain good visual quality for diagnosis. By applying the proposed denoising algorithm to the filtered projection images, the X-ray exposure level necessary for the CBCT scan can be reduced by up to 60% while obtaining clinically acceptable image quality. This denoising result indicates that in the clinical application of CBCT breast imaging, the patient radiation dose can be significantly reduced.

Algorithms↗

PubMed2004

X-ray scatter correction algorithm for cone beam CT imaging.

Developing and optimizing an x-ray scatter control and reduction technique is one of the major challenges for cone beam computed tomography (CBCT) because CBCT will be much less immune to scatter than fan-beam CT. X-ray scatter reduces image contrast, increases image noise and introduces reconstruction error into CBCT. To reduce scatter interference, a practical algorithm that is based upon the beam stop array technique and image sequence processing has been developed on a flat panel detector-based CBCT prototype scanner. This paper presents a beam stop array-based scatter correction algorithm and the evaluation results through phantom studies. The results indicate that the beam stop array-based scatter correction algorithm is practical and effective to reduce and correct x-ray scatter for a CBCT imaging task.

Algorithms↗

PubMed2003

Filling the Radon domain in computed tomography by local convex combination.

Radon data interpolation is a necessary procedure in computed tomography (CT), especially for reconstruction from divergent beam scanning. In a polar-grid representation, the Radon data of a fanbeam projection are populated on an arc, rather on a radial line. Collectively, the Radon data generated from a fanbeam CT system are unevenly populated: The population becomes sparser as the polar distance increases. In CT reconstruction, the Fourier central slice theorem requires a radial scanline full of Radon data. Therefore the vacant entries of a scanline must be filled by interpolation. In addition, interpolation is also required in polar-to-Cartesian conversion. In this paper we propose a practical interpolation technique for filling the vacant entries by local convex combination. It is a linear interpolant that generates a value for a grid point from the available data lying in its neighborhood, by a weighted average, with the weights corresponding to the inverse distances. In fact, the linear convex combination serves as a general flat-smoothing operation in filling a vacancy. Specifically, this technique realizes a variety of linear interpolations, including nearest-neighbor replication, two-point collinear, three-point triangulation, and four-point quadrilateral, and local extrapolation, in a unified framework. Algorithms and a simulation demonstration are provided.

Algorithms↗

PubMed2003

Why should breast tumour detection go three dimensional?

Although x-ray mammography is widely developed for breast tumour detection, it suffers from spatial superposition in its two-dimensional (2D) representation of a three-dimensional (3D) breast structure. Accordingly, 3D breast imaging, such as cone-beam computed tomography (CT), arises at the historic moment. In this paper, we theoretically elucidate the spatial superposition effect associated with x-ray mammography on breast tumour detection. This explanation is based on the line integral of x-ray traversing a composite breast model. As a result, we can characterize the difficulty of detecting small tumours in terms of local intensity contrast in x-ray images. In comparison, we also introduce cone-beam CT breast imaging for 3D breast volume representation, which offers advantages for breast mass segmentation and measurement. The discussion is demonstrated with an experiment with a breast surgical specimen. In conclusion, we strongly believe that 3D volumetric representation allows for more accurate breast tumour detection.

Algorithms↗

PubMed2003

Flat panel detector-based cone beam computed tomography with a circle-plus-two-arcs data acquisition orbit: preliminary phantom study.

Cone beam computed tomography (CBCT) has been investigated in the past two decades due to its potential advantages over a fan beam CT. These advantages include (a) great improvement in data acquisition efficiency, spatial resolution, and spatial resolution uniformity, (b) substantially better utilization of x-ray photons generated by the x-ray tube compared to a fan beam CT, and (c) significant advancement in clinical three-dimensional (3D) CT applications. However, most studies of CBCT in the past are focused on cone beam data acquisition theories and reconstruction algorithms. The recent development of x-ray flat panel detectors (FPD) has made CBCT imaging feasible and practical. This paper reports a newly built flat panel detector-based CBCT prototype scanner and presents the results of the preliminary evaluation of the prototype through a phantom study. The prototype consisted of an x-ray tube, a flat panel detector, a GE 8800 CT gantry, a patient table and a computer system. The prototype was constructed by modifying a GE 8800 CT gantry such that both a single-circle cone beam acquisition orbit and a circle-plus-two-arcs orbit can be achieved. With a circle-plus-two-arcs orbit, a complete set of cone beam projection data can be obtained, consisting of a set of circle projections and a set of arc projections. Using the prototype scanner, the set of circle projections were acquired by rotating the x-ray tube and the FPD together on the gantry, and the set of arc projections were obtained by tilting the gantry while the x-ray tube and detector were at the 12 and 6 o'clock positions, respectively. A filtered backprojection exact cone beam reconstruction algorithm based on a circle-plus-two-arcs orbit was used for cone beam reconstruction from both the circle and arc projections. The system was first characterized in terms of the linearity and dynamic range of the detector. Then the uniformity, spatial resolution and low contrast resolution were assessed using different phantoms mainly in the central plane of the cone beam reconstruction. Finally, the reconstruction accuracy of using the circle-plus-two-arcs orbit and its related filtered backprojection cone beam volume CT reconstruction algorithm was evaluated with a specially designed disk phantom. The results obtained using the new cone beam acquisition orbit and the related reconstruction algorithm were compared to those obtained using a single-circle cone beam geometry and Feldkamp's algorithm in terms of reconstruction accuracy. The results of the study demonstrate that the circle-plus-two-arcs cone beam orbit is achievable in practice. Also, the reconstruction accuracy of cone beam reconstruction is significantly improved with the circle-plus-two-arcs orbit and its related exact CB-FPB algorithm, as compared to using a single circle cone beam orbit and Feldkamp's algorithm.

Algorithms↗

PubMed2002

Cone-beam volume CT breast imaging: feasibility study.

X-ray projection mammography, using a film/screen combination, or digital techniques, has proven to be the most effective imaging modality currently available for early detection of breast cancer. However, the inherent superimposition of structures makes a small carcinoma (a few millimeters in size) difficult to detect when it is occult or in dense breasts, leading to a high false-positive biopsy rate. Cone-beam x-ray-projection-based volume imaging using flat panel detectors (FPDs) may allow obtaining three-dimensional breast images, resulting in more accurate diagnosis of structures and patterns of lesions while eliminating the hard compression of breasts. This article presents a novel cone-beam volume computed tomographic breast imaging (CBVCTBI) technique based on the above techniques. Through a variety of computer simulations, the key issues of the system and imaging techniques were addressed, including the x-ray imaging geometry and corresponding reconstruction algorithms, x-ray characteristics of breast tissue and lesions, x-ray setting techniques, the absorbed dose estimation, and the quantitative effect of x-ray scattering on image quality. The preliminary simulation results support the proposed CVBCTBI modality for breast imaging in respect to its feasibility and practicability. The absorbed dose level is comparable to that of current mammography and will not be a prominent problem for this imaging technique. Compared to conventional mammography, the proposed imaging technique with isotropic spatial resolution will potentially provide significantly better low-contrast detectability of breast tumors and more accurate location of breast lesions.

Algorithms↗