Impact of beam-hardening corrections on proton relative stopping power estimates from single- and dual-energy CT

A major contributing factor to proton range uncertainty is the conversion of computed tomography (CT) Hounsfield units (HU) to proton relative stopping power (RSP). This uncertainty is heightened in the presence of X-ray beam-hardening artifact (BHA), which has two manifestations: cupping and streaking, especially in and near bone tissue. This uncertainty can affect the accuracy of proton RSP calculation for treatment planning in proton radiotherapy. Dual-energy CT (DECT) and iterative beam-hardening correction (iBHC) both show promise in mitigating CT BHA. This present work attempts to analyze the relative robustness of iBHC and DECT techniques on both manifestations of BHA. The stoichiometric method for HU to RSP conversion was used for single-energy CT (SECT) and DECT-based monochromatic techniques using a tissue substitute phantom. Cupping BHA was simulated by measuring the HU of a bone substitute plug in wax/3D-printed phantoms of increasing size. Streaking BHA was simulated by placing a solid water plug between two bone plugs in a wax phantom. Finally, the effect of varying calibration phantom size on RSP was calculated in an anthropomorphic head phantom. The RSP decreased -0.002 cm^-1 as phantom size increased for SECT but remained largely constant when iBHC applied or with DECT techniques. The RSP varied a maximum of 2.60% in the presence of streaking BHA in SECT but was reduced to 1.40% with iBHC. For DECT techniques, the maximum difference was 2.40%, reduced to 0.6% with iBHC. Comparing calibration phantoms of 20- and 33-cm diameter, maximum voxel differences of 5 mm in the water-equivalent thickness were observed in the skull but reduced to 1.3 mm with iBHC. The DECT techniques excelled in mitigating cupping BHA, but streaking BHA still could be observed. The use of iBHC reduced RSP variation with BHA in both SECT and DECT techniques.

Accuracy of proton stopping power estimation of silicone breast implants with single and dual-energy CT calibration techniques

A major contributing factor to proton range uncertainty is the conversion of computed tomography (CT) Hounsfield Units (HU) to proton relative stopping power (RSP). This uncertainty is elevated with implanted devices, such as silicone breast implants when computed with single energy CT (SECT). In recent years, manufacturers have introduced implants with variations in gel cohesivity. Deriving the RSP for these implants from dual-energy CT (DECT) can result in a marked reduction of the error associated with SECT. In this study, we investigate the validity of DECT calibration of HU to RSP on silicone breast implants of varying cohesivity levels. A DECT capable scanner was calibrated using the stoichiometric method of Bourque et al for SECT and DECT using a tissue substitute phantom. Three silicone breast implants of increasing gel cohesivity were measured in a proton beam of clinical energy to determine ground-truth RSP and water equivalent thickness (WET). These were compared to SECT-derived RSP at three CT spectrum energies and DECT with two energy pairs (80/140 kVp and 100/140 kVp) as obtained from scans with and without an anthropomorphic phantom. The RSP derived from parameters estimates from CT vendor-specific software (syngo.via) was compared. The WET estimates from SECT deviated from MLIC ground truth approximately +11%-19%, which would result in overpenetration if used clinically. Both the Bourque calibration and syngo.via WET estimates from DECT yielded error ≤0.5% from ground truth; no significant difference was found between models of varying gel cohesivity levels. WET estimates without the anthropomorphic phantom were significantly different than ground truth for the Bourque calibration. From these results, gel cohesivity had no effect on proton RSP. User-generated DECT calibration can yield comparably accurate RSP estimates for silicone breast implants to vendor software methods. However, care must be taken to account for beam hardening effects.

Solution of the Michaelis-Menten equation using the decomposition method

We present a low-order recursive solution to the Michaelis-Menten equation using the decomposition method. This solution is algebraic in nature and provides a simpler alternative to numerical approaches such as differential equation evaluation and root-solving techniques that are currently used to compute substrate concentration in the Michaelis-Menten equation. A detailed characterization of the errors in substrate concentrations computed from decomposition, Runge-Kutta, and bisection methods over a wide range of s(0) : K(m) values was made by comparing them with highly accurate solutions obtained using the Lambert W function. Our results indicated that solutions obtained from the decomposition method were usually more accurate than those from the corresponding classical Runge-Kutta methods. Moreover, these solutions required significantly fewer computations than the root-solving method. Specifically, when the stepsize was 0.1% of the total time interval, the computed substrate concentrations using the decomposition method were characterized by accuracies on the order of 10(-8) or better. The algebraic nature of the decomposition solution and its relatively high accuracy make this approach an attractive candidate for computing substrate concentration in the Michaelis-Menten equation.

Establishing a threshold for 131I bioassay in nuclear medicine personnel

Since the late 1970's, manufacturers in nuclear medicine have reformulated the I solution to reduce the volatility of the iodine. There has also been an increase in use of the iodide in encapsulated form. Per the requirement of the current U.S. Nuclear Regulatory Commission (U.S. NRC) regulation, with the available results on the volatility of the reformulated radioiodine, we review the I bioassay program for nuclear medicine workers. Our analysis shows the threshold quantity for bioassay monitoring for the routine use of I in nuclear medicine is much higher than the criteria set in U.S. NRC Regulatory Guide 8.20. The latter is a broad bioassay guideline for the general usage of radioactive iodine. For treatment of thyroid carcinoma and hyperthyroidism, a single therapeutic I dose large enough to yield a detectable thyroid burden is very unlikely to occur in a nuclear medicine clinic. Accidental ingestion or inhalation would be an exception to our conclusion. Based on this analysis, we propose a new bioassay policy for the routine use of I in nuclear medicine clinics.

The effects of dialysis on 131I kinetics and dosimetry in thyroid cancer patients--a pharmacokinetic model

Currently, an accepted post-surgical treatment of patients with thyroid carcinoma is administration of an ablative dose of I. This treatment is well established based on extensive experience and modeling. However, for patients with renal disease, reduced iodine removal rates result in controversial thyroid doses and potentially excessive red bone marrow doses. There are differences of opinion regarding I dose recommendations ranging from a reduction in dose to an increase in dose compared with conventional amounts. Determination of suitable doses must take into account varying dialysis protocols and absorbed dose considerations to the thyroid and sensitive tissues such as red bone marrow. The specific aim of this study was to develop a simple yet comprehensive compartmental model for I kinetics in patients with thyroid carcinoma and end stage renal disease, which accounts for dialysis and provides absorbed dose estimates for the thyroid as well as the red bone marrow. STELLA, a compartmental modeling software program, was used to develop a kinetic model that includes the blood pool, thyroid, gastrointestinal tract, kidneys, bladder, and a conventional dialysis machine. Benchmarking was performed to demonstrate the validity of the model with data obtained from ICRP 30 and MIRD Dose Estimate Report No. 5. Iodine kinetics were simulated for normal patients, thyroid cancer patients, and patients with thyroid cancer and renal failure undergoing two standard types of dialysis, hemodialysis and continuous ambulatory peritoneal dialysis (CAPD). Results in this work show that thyroid doses to patients with thyroid cancer and renal failure on hemodialysis or CAPD are slightly higher than doses to patients with thyroid cancer and normal renal function. These results further indicate that red bone marrow doses to patients with thyroid cancer and renal failure on dialysis can be significantly higher than red bone marrow doses to patients with thyroid cancer and normal renal function, and thus these patients could benefit from a reduction in administered activity. Thyroid doses and red bone marrow doses to patients on standard hemodialysis depend on both dialysis frequency and the time interval between administration and first dialysis. The results in this study provide guidelines on how much activity a patient on dialysis should receive based on thyroid and red bone marrow absorbed dose (Gy MBq) considerations. This study should help to clarify some of the contradictory recommendations regarding I dose for thyroid carcinoma patients with renal failure.

Request form history, clinical indication, and yield of brain magnetic resonance studies

PURPOSE

To investigate whether improved clinical history allows the radiologist to better predict the pretest probability of obtaining a positive or negative result from a magnetic resonance (MR) examination.

MATERIALS AND METHODS

Six neuroradiologists prospectively reviewed 100 consecutive requests for brain MR examinations and sequentially assessed 1) quality of written history, 2) degree of indication for requested study, and 3) any pertinent new information found during chart review that may have altered the degree of indication. MR yield was correlated with the degree of indication assessed before and after chart review.

RESULTS

Most request form histories were judged as poor (63%), and chart review reduced the overall indications for MR examinations, as there was a tendency for high-indication requests to migrate to the low-indication category. Based on request form history alone, the yields for low- and high-indication studies were 13% and 37%, respectively. Correlations between MR yield and indication after chart review improved significantly (P < 0.05) with 2% and 61% for low and high indications, respectively. Sensitivity and specificity for a positive MR yield were 71% and 62%, respectively, for the indication judged by the request history alone, and 96% and 80%, respectively, after chart review. Positive and negative prediction rates were 37% and 87%, respectively, for the indication judged by the request history alone, and 61% and 98%, respectively, for the indication judged after chart review.

CONCLUSION

Based on our limited data, most request form histories were inadequate, and essential information available in the chart before MR examinations was frequently missing from the request forms. When adequate information was provided, the indication for the studies as judged by the radiologists predicted the MR yield more accurately, particularly for those requests with low indication. Therefore, our study suggests that MR imaging (MRI) may be used more effectively when pertinent clinical history is available. However, our study is limited and further studies are needed to confirm our results.

Qualitative assessment of cervical spinal stenosis: observer variability on CT and MR images

BACKGROUND AND PURPOSE

Several studies have been undertaken to validate quantitative methods of evaluating cervical spinal stenosis. This study was performed to assess the degree of interobserver and intraobserver agreement in the qualitative evaluation of cervical spinal stenosis on CT myelograms and MR images.

METHODS

Cervical MR images and CT myelograms of 38 patients were evaluated retrospectively. Six neuroradiologists with various backgrounds and training independently assessed the level, degree, and cause of stenosis on either MR images or CT myelograms. Unknown to the evaluators, 16 of the patients were evaluated twice to determine intraobserver variability.

RESULTS

Interobserver agreement among the radiologists with regard to level, degree, and cause of stenosis on CT myelograms showed kappa values of 0.50, 0.26, and 0.32, respectively, and on MR images showed kappa values of 0.60, 0.31, and 0.22, respectively. Intraobserver agreement with regard to level, degree, and cause of stenosis on CT myelograms showed mean kappa values of 0.69, 0.41, and 0.55, respectively, and on MR images showed mean kappa values of 0.80, 0.37, and 0.40, respectively.

CONCLUSION

MR imaging and CT myelographic evaluation of cervical spinal stenosis by using current qualitative methods results in significant variation in image interpretation.

Pitfalls of radiation survey after brachytherapy implant removal preceding Tl-201 study

An unexpected elevated postimplant radiation survey is described in an elderly patient with an interstitial low-dose-rate iridium-192 (Ir-192) needle implant for endometrial cancer. The elevated activity was related to prolonged clearance of Tl-201 from a cardiac study that had been performed 7 days earlier. The Tl-201 accumulated in the soft tissue, particularly the colon, resulting in increased survey readings over the abdomen and raising concern that an Ir-192 source remained within the patient. This case shows that delayed excretion of a diagnostic radionuclide agent can cause elevated activity high enough to confound postradiotherapy implant survey readings. The estimated surface exposure from a single iridium source left in the pelvis was determined using a phantom study. Possible factors causing decreased excretion of Tl-201 in a patient with heart disease, arteriosclerotic vascular disease, previous pelvic radiation therapy, and a brachytherapy procedure are discussed. A preloading radiation survey is recommended in patients who have had previous nuclear medicine studies involving radionuclides with long half-lives.

Parameter estimation using a direct solution of the integrated Michaelis-Menten equation

A novel method of estimating enzyme kinetic parameters is presented using the Lambert omega function coupled with nonlinear regression. Explicit expressions for the substrate and product concentrations in the integrated Michaelis-Menten equation were obtained using the omega function which simplified kinetic parameter estimation as root-solving and numerical integration of the Michaelis-Menten equation were avoided. The omega function was highly accurate in describing the substrate and product concentrations in the integrated Michaelis-Menten equation with an accuracy of the order of 10(-16) when double precision arithmetic was used. Progress curve data from five different experimental systems were used to demonstrate the suitability of the omega function for kinetic parameter estimation. In all cases, the kinetic parameters obtained using the omega function were almost identical to those obtained using the conventional root-solving technique. The availability of highly efficient algorithms makes the computation of omega simpler than root-solving or numerical integration. The accuracy and simplicity of the omega function approach make it an attractive alternative for parameter estimation in enzyme kinetics.