PubMed Health⌕ Search

Biomedical subjects

Keh-Shih Chuang

Publications and source records attributed to Keh-Shih Chuang.

12 recordsLinked to original sources

Polygon interpolation for serial cross sections.

In this paper, a new technique for contour interpolation between slices is presented. We assumed that contour interpolation is equivalent to the interpolation of a polygon that approximates the object shape. The location of each polygon vertex is characterized by a set of parameters. Polygon interpolation can be performed on these parameters. These interpolated parameters are then used to reconstruct the vertices of the new polygon. Finally, the contour is approximated from this polygon using a cubic spline interpolation. This new technique takes into account the shape, the translation, the size, and the orientation of the object's contours. A comparison with regular shape-based interpolation is made on several object contours. The preliminary results show that this new method yields a better contour and is computationally more efficient than shape-based interpolation. This technique can be applied to gray-level images too. The interpolation result of an MR image does not show artifact of intermediate substance commonly seen in a typical linear gray-level interpolation.

Algorithms↗

A blurring index for medical images.

This study was undertaken to investigate a useful image blurring index. This work is based on our previously developed method, the Moran peak ratio. Medical images are often deteriorated by noise or blurring. Image processing techniques are used to eliminate these two factors. The denoising process may improve image visibility with a trade-off of edge blurring and may introduce undesirable effects in an image. These effects also exist in images reconstructed using the lossy image compression technique. Blurring and degradation in image quality increases with an increase in the lossy image compression ratio. Objective image quality metrics [e.g., normalized mean square error (NMSE)] currently do not provide spatial information about image blurring. In this article, the Moran peak ratio is proposed for quantitative measurement of blurring in medical images. We show that the quantity of image blurring is dependent upon the ratio between the processed peak of Moran's Z histogram and the original image. The peak ratio of Moran's Z histogram can be used to quantify the degree of image blurring. This method produces better results than the standard gray level distribution deviation. The proposed method can also be used to discern blurriness in an image using different image compression algorithms.

Algorithms↗

Reevaluation of the emergency planning zone for nuclear power plants in Taiwan using MACCS2 code.

According to government regulations, the emergency planning zone (EPZ) of a nuclear power plant (NPP) must be designated before operation and reevaluated every 5 years. Corresponding emergency response planning (ERP) has to be made in advance to guarantee that all necessary resources are available under accidental releases of radioisotope. In this study, the EPZ for each of the three operating NPPs, Chinshan, Kuosheng, and Maanshan, in Taiwan was reevaluated using the MELCOR Accident Consequence Code System 2 (MACCS2) developed by Sandia National Laboratory. Meteorological data around the nuclear power plant were collected during 2003. The source term data including inventory, sensible heat content, and timing duration, were based on previous PRA information of each plant. The effective dose equivalent and thyroid dose together with the related individual risk and societal risk were calculated. By comparing the results to the protective action guide and related safety criteria, 1.5, 1.5, and 4.5km were estimated for Chinshan, Kuosheng, and Maanshan NPPs, respectively. We suggest that a radius of 5.0km is a reasonably conservative value of EPZ for each of the three operating NPPs in Taiwan.

Journal Article↗

Fuzzy c-means clustering with spatial information for image segmentation.

A conventional FCM algorithm does not fully utilize the spatial information in the image. In this paper, we present a fuzzy c-means (FCM) algorithm that incorporates spatial information into the membership function for clustering. The spatial function is the summation of the membership function in the neighborhood of each pixel under consideration. The advantages of the new method are the following: (1) it yields regions more homogeneous than those of other methods, (2) it reduces the spurious blobs, (3) it removes noisy spots, and (4) it is less sensitive to noise than other techniques. This technique is a powerful method for noisy image segmentation and works for both single and multiple-feature data with spatial information.

Algorithms↗

Scatter correction for 3D PET using beam stoppers combined with dual-energy window acquisition: a feasibility study.

Fully three-dimensional (3D) positron emission tomography (PET) can achieve high sensitivity of coincidence events, but the absence of inter-slice septa inevitably leads to increased scattered events. The scattered events can represent as much as 50% of the total detected events. In this research, we proposed a scatter correction method for 3D PET based on beam stoppers and dual-energy window acquisition. The beam stoppers were placed surrounding the object to attenuate primary beams. The scatter fractions were directly estimated at those blocked lines of response and then the entire scatter fraction distribution was recovered using the dual-energy window ratio as reference. The performance was evaluated by using Monte Carlo simulations of various digital phantoms. For the Utah phantom study, the proposed method accurately estimated the scatter fraction distribution, and improved image contrast and quantification based on four different quality indices as performance measures. For the non-homogeneous Zubal phantom, the simulated results also demonstrated that the proposed method achieved a better restoration of image contrast than the dual-energy window method. We conclude that the proposed scatter correction method could effectively suppress various kinds of scattered events, including multiple scatter and scatter from outside the field of view.

Computer Simulation↗

A maximum likelihood expectation maximization algorithm with thresholding.

The maximum likelihood expectation maximization (MLEM) algorithm has several advantages over the conventional filtered back-projection (FBP) for image reconstruction. However, the slow convergence and the high computational cost for its practical implementation have limited its clinical applications. This study proposes the incorporation of a thresholding technique in both the MLEM and ordered subsets EM (OSEM) algorithm to accelerate convergence. The threshold is set to c*m, where m is the mean pixel value of the whole image. The reconstruction time is proportional to the total number of pixels, so a thresholding technique that nullifies the value of a pixel if it falls below a threshold, can effectively remove the non-active pixels and substantially accelerate reconstruction. Preliminary tests on simulated PET data reveal that the thresholding technique accelerates the convergence rate and reduce error in the reconstructed image. The reconstruction performance improves with the increase of the threshold level and the MSE reaches minimum for c value equals to about 1.

Algorithms↗

A three-dimensional registration method for automated fusion of micro PET-CT-SPECT whole-body images.

Micro positron emission tomography (PET) and micro single-photon emission computed tomography (SPECT), used for imaging small animals, have become essential tools in developing new pharmaceuticals and can be used, among other things, to test new therapeutic approaches in animal models of human disease, as well as to image gene expression. These imaging techniques can be used noninvasively in both detection and quantification. However, functional images provide little information on the structure of tissues and organs, which makes the localization of lesions difficult. Image fusion techniques can be exploited to map the functional images to structural images, such as X-ray computed tomography (CT), to support target identification and to facilitate the interpretation of PET or SPECT studies. Furthermore, the mapping of two functional images of SPECT and PET on a structural CT image can be beneficial for those in vivo studies that require two biological processes to be monitored simultaneously. This paper proposes an automated method for registering PET, CT, and SPECT images for small animals. A calibration phantom and a holder were used to determine the relationship among three-dimensional fields of view of various modalities. The holder was arranged in fixed positions on the couches of the scanners, and the spatial transformation matrix between the modalities was held unchanged. As long as objects were scanned together with the holder, the predetermined matrix could register the acquired tomograms from different modalities, independently of the imaged objects. In this work, the PET scan was performed by Concorde's microPET R4 scanner, and the SPECT and CT data were obtained using the Gamma Medica's X-SPECT/CT system. Fusion studies on phantoms and animals have been successfully performed using this method. For microPET-CT fusion, the maximum registration errors were 0.21 mm +/- 0.14 mm, 0.26 mm +/- 0.14 mm, and 0.45 mm +/- 0.34 mm in the X (right-left), Y (upper lower), and Z (rostral-caudal) directions, respectively; for the microPET-SPECT fusion, they were 0.24 mm +/- 0.14 mm, 0.28 mm +/- 0.15 mm, and 0.54 mm +/- 0.35 mm in the X, Y, and Z directions, respectively. The results indicate that this simple method can be used in routine fusion studies.

Algorithms↗

A statistical method for evaluation quality of medical images: a case study in bit discarding and image compression.

Many studies have been performed on quality evaluation for subtle differences in medical images. However, only limited success has been achieved. In this paper, medical images were prior manipulated by denoising, lossy compression and filtering. The Moran statistics is then applied to extract spatial information of images and using Kolmogorov-Smirnov (KS) test to determine whether the manipulated and original images differ significantly. Results show that on average discarding 1-2 bits in T1 and CR images or 2-3 bits in T2 and body CT images are indistinguishable. This method is also applied to a reconstructed MR, body CT image and an electronic SMPTE (Society of Motion Picture and Television Engineer) phantom from lossy image compression software. Compression ratios of 16:1 for a MR image, 8-9:1 for a cropped body CT image, 7:1 and 5:1 for high- and low-resolution regions in electronic phantom is proved undifferentiated from original. The proposed method is useful for complementing the human visual system, to optimize the performance of image compression technique.

Data Interpretation, Statistical↗

Quality degradation in lossy wavelet image compression.

The objective of this study was to develop a method for measuring quality degradation in lossy wavelet image compression. Quality degradation is due to denoising and edge blurring effects that cause smoothness in the compressed image. The peak Moran z histogram ratio between the reconstructed and original images is used as an index for degradation after image compression. The Moran test is applied to images randomly selected from each medical modality, computerized tomography, magnetic resonance imaging, and computed radiography and compressed using the wavelet compression at various levels. The relationship between the quality degradation and compression ratio for each image modality agrees with previous reports that showed a preference for mildly compressed images. Preliminary results show that the peak Moran z histogram ratio can be used to quantify the quality degradation in lossy image compression. The potential for this method is applications for determining the optimal compression ratio (the maximized compression without seriously degrading image quality) of an image for teleradiology.

Diagnostic Imaging↗

A novel image quality index using Moran I statistics.

Measurement of image quality is very important for various applications such as image compression, restoration and enhancement. Conventional methods (e.g., mean squared error; MSE) use error summation to measure quality change pixel by pixel and do not correlate well with subjective quality measurement. This is due to the fact that human eyes extract structural information from the viewing field. In this study a new quality index using a Moran I statistics is proposed. The Moran statistic that measures the sharpness from a local area is a good index of quality as most image processing techniques alter the smoothness of the image. Preliminary results show that the new quality index outperforms the MSE significantly under various types of image distortions.

Algorithms↗

Determination of beam intensity in a single step for IMRT inverse planning.

In intensity modulated radiotherapy (IMRT), targets are treated by multiple beams at different orientations each with spatially-modulated beam intensities. This approach spreads the normal tissue dose to a greater volume and produces a higher dose conformation to the target. In general, inverse planning is used for IMRT treatment planning. The inverse planning requires iterative calculation of dose distribution in order to optimize the intensity profile for each beam and is very computation intensive. In this paper, we propose a single-step method utilizing a figure of merit (FoM) to estimate the beam intensities for IMRT treatment planning. The FoM of a ray is defined as the ratio between the delivered tumour dose and normal tissue dose and is a good index for the dose efficacy of the ray. To maximize the beam utility, it is natural to irradiate the tumour with intensity of each ray proportional to the value of the FoM. The nonuniform beam intensity profiles are then fixed and the weights of the beam are determined iteratively in order to yield a uniform tumour dose. In this study, beams are employed at equispaced angles around the patient. Each beam with its field size that just covers the tumour is divided into a fixed number of beamlets. The FoM is calculated for each beamlet and this value is assigned to be the beam intensity. Various weighting factors are incorporated in the FoM computation to accommodate different clinical considerations. Two clinical datasets are used to test the feasibility of the algorithm. The resultant dose-volume histograms of this method are presented and compared to that of conformal therapy. Preliminary results indicate that this method reduces the critical organ doses at a small expense of uniformity in tumour dose distribution. This method estimates the beam intensity in one single step and the computation time is extremely fast and can be finished in less than one minute using a regular PC.

Algorithms↗

Comprehensive measurement for clinical dosimetry of tertiary multileaf collimator.

To initiate the use of a tertiary multileaf collimator (MLC) in the clinic, a set of dosimetry data for clinical use of the MLC, the secondary field size jaw and the MLC tracked by the jaw were measured. The dose calculation technique from the commissioned jaw field data was established. The dosimetry characteristics included absolute output (Dw), collimator scatter factor (Sc), total scatter factor (Scp), phantom scatter factor (Sp), percentage depth dose (PDD), tissue-maximum ratio (TMR), and peak scatter factor (PSF). The absolute output of the MLC field was +5% to +2% greater than that of the same jaw size field from 4 x 4 to 24 x 24 cm2 fields. The variation of Sc and Scp ranged from 4 x 4 to 24 x 24 cm2 fields and were less than that of the jaw fields, while the Sp, PDD and TMR values remained the same. Importantly, when the MLC-only field was performed without the collimator jaws tracking close to the field segments, greater output was delivered, and PSFs should be used to calculate the MLC field output.

Humans↗