Denoised and texture enhanced MVCT to improve soft tissue conspicuity

PURPOSE

MVCT images have been used in TomoTherapy treatment to align patients based on bony anatomies but its usefulness for soft tissue registration, delineation, and adaptive radiation therapy is limited due to insignificant photoelectric interaction components and the presence of noise resulting from low detector quantum efficiency of megavoltage x-rays. Algebraic reconstruction with sparsity regularizers as well as local denoising methods has not significantly improved the soft tissue conspicuity. The authors aim to utilize a nonlocal means denoising method and texture enhancement to recover the soft tissue information in MVCT (DeTECT).

METHODS

A block matching 3D (BM3D) algorithm was adapted to reduce the noise while keeping the texture information of the MVCT images. Following imaging denoising, a saliency map was created to further enhance visual conspicuity of low contrast structures. In this study, BM3D and saliency maps were applied to MVCT images of a CT imaging quality phantom, a head and neck, and four prostate patients. Following these steps, the contrast-to-noise ratios (CNRs) were quantified.

RESULTS

By applying BM3D denoising and saliency map, postprocessed MVCT images show remarkable improvements in imaging contrast without compromising resolution. For the head and neck patient, the difficult-to-see lymph nodes and vein in the carotid space in the original MVCT image became conspicuous in DeTECT. For the prostate patients, the ambiguous boundary between the bladder and the prostate in the original MVCT was clarified. The CNRs of phantom low contrast inserts were improved from 1.48 and 3.8 to 13.67 and 16.17, respectively. The CNRs of two regions-of-interest were improved from 1.5 and 3.17 to 3.14 and 15.76, respectively, for the head and neck patient. DeTECT also increased the CNR of prostate from 0.13 to 1.46 for the four prostate patients. The results are substantially better than a local denoising method using anisotropic diffusion.

CONCLUSIONS

The authors showed that it is feasible to extract more soft tissue contrast information from the noisy MVCT images using a nonlocal means 3D block matching method in combination with saliency maps, revealing information that was originally unperceivable to human observers.

Application of an independent dose calculation software for estimating the impact of inter-fractional setup shifts in Helical Tomotherapy treatments

The purpose of this study is to validate the capability of in-house independent point dose calculation software to be used as a second check for Helical Tomotherapy treatment plans. The software performed its calculations in homogenous conditions (using the Cheese phantom, which is a cylindrical phantom with radius 15 cm and length 18 cm) using a factor-based algorithm. Fifty patients, who were treated for pelvic (10), prostate (14), lung (10), head & neck (12) and brain (4) cancers, were used. Based on the individual patient kVCT images and the pretreatment MVCT images for each treatment fraction, the corresponding daily patient setup shifts in the IEC-X, IEC-Y, and IEC-Z directions were registered. For each patient, the registered fractional setup shifts were grouped into systematic and random shifts. The average systematic dosimetric variations showed small dose deviation for the different cancer types (1.0%-3.0%) compared to the planned dose. Of the fifty patients, only three had percent differences larger than 5%. The average random dosimetric variations showed relatively small dose deviations (0.2%-1.1%) compared to the planned dose. None of the patients had percent differences larger than 5%. By examining the individual fractions of each patient, it is observed that only in 31 out of 1358 fractions the percent differences exceeded the border of 5%. These results indicate that the overall dosimetric impact from systematic and random variations is small and that the software is a capable platform for independent point dose validation for the Helical Tomotherapy modality.

Diagnosis of liver cancer and cirrhosis by the fluorescence spectra of blood and urine

Liver cancer or hepato cellular carcinoma (HCC) is a serious malady with only 10% survival rate and is fatal next only to pancreatic cancer. This disease is conventionally detected and diagnosed by ultra sound, CT or MRI scans which are quite expensive. Also the discrimination between cirrhosis and HCC, by these imaging techniques, is poor. The conventional tissue biopsy is quite invasive and painful. In the new diagnostic procedure presented in this paper we have obtained fluorescence emission spectra with excitation at 400 nm and the synchronous emission spectra (Δλ = 10 nm) for a set of blood and urine samples (Normal control N = 25, Liver Malfunction N = 58). Based on the ratio fluorometric parameters, all the three liver maladies (minor inflammation like Hepatitis C, serious diseases like Cirrhosis and hepatoma) could be detected and discriminated with an accuracy of about 80%. Hence this inexpensive, non invasive, optical technique may have significant impact in screening, diagnosis and also prognosis of HCC in large segment of people in the populous Asian countries.

Quantitative characterization of tomotherapy MVCT dosimetry

Megavoltage computed tomography (MVCT) is used as image guidance for patient setup in almost every tomotherapy treatment. Frequent use of ionizing radiation for image guidance has raised concern of imaging dose. The purpose of this work is to quantify and characterize tomotherapy MVCT dosimetry. Our dose calculation was based on a commissioned dose engine, and the calculation result was compared with film measurement. We studied dose profiles, center dose, maximal dose, surface dose, and mean dose on homogeneous cylindrical water phantoms of various diameters for various scanning parameters, including 3 different jaw openings (of nominal value J4, J1, and J0.1) and couch speeds (fine, normal, and coarse). The comparison between calculation and film measurement showed good agreement. In particular, the thread pattern on the film of the helical delivery matched very well with calculation. For the J1 jaw and coarse imaging mode, the maximum difference between calculation and measurement was about 6% of the center dose. Calculation on various sizes of synthesized phantoms showed that the center dose decreases almost linearly as the phantom diameter increases, and that the fine mode (couch speed of 4mm/rotation) received twice the dose of the normal mode (couch speed of 8mm/rotation) and 3 times that of the coarse mode (couch speed of 12mm/rotation) as expected. The maximal dose ranged from 100% to ∼200% of the center dose, with increasing ratios for larger phantoms, smaller jaws, and faster couch speed. For all jaw settings and couch speeds, the mean dose and average surface dose vary from 95% to 125% of the center dose with increasing ratios for larger phantoms. We present a quantitative dosimetric characterization of the tomotherapy MVCT in terms of scanning parameters, phantom size, center dose, maximal dose, surface dose, and mean dose. The results can provide an overall picture of dose distribution and a reference data set that enables estimation of CT dose index for the tomotherapy MVCT.

B-Spline Registration Based on New Concept of an Intelligent Masking Procedure and GPU Computations for the Head and Neck Adaptive Tomotherapy

The deformable image registration (DIR) procedure has been optimized for helical tomotherapy. The data on registration shifts obtained on matching planning image with pre-treatment megavoltage CT are used in our software for acceleration of the first step (rigid registration) of the DIR procedure and for implementation of the B-Spline algorithm with intelligent masking. Priorities of the masks were automatically calculated based on disagreement detected during rigid registration. Evaluation tasks included: (a) comparison of accuracy and rate for schemes of pre-registered and non-registered images; (b) qualification of the effectiveness of the intelligent masking process, and (c) determination of acceleration of achievable with GPU computing. A specially designed head and neck phantom used for evaluation included structures with controlled changes of position, volume, density, and shape. Re-contouring procedures were performed with an Adaptive Planning software (Tomotherapy Inc.). No statistical difference was observed in accuracy of DIR based on structure position match on the tomotherapy unit and non pre-registered images (p > 0.7). Using pre-registered data reduces the total time required for execution of the elastic registration procedure by 5%. These data are also necessary for intelligent masking procedure during B-Spine registration. Intelligent masking procedure increases accuracy of the registration for a masked structure (p < 0.04) without decreasing the accuracy in non-masked tissues and additionally reduces the total time by 13%. GPU computations speed up procedure 30 times. GPU computing of the DIR in current status of our investigation could be realized in a relatively short time after pre-treatment imaging. The proposed approach can be used in the routine assessment of anatomic changes occurring in healthy tissue during the course of radiotherapy. Further developments will be concentrated on the full integration of DIR computations in the imaging and treatment process of helical tomotherapy.