Monte Carlo dose calculation improvements for low energy electron beams using eMC

The electron Monte Carlo (eMC) dose calculation algorithm in Eclipse (Varian Medical Systems) is based on the macro MC method and is able to predict dose distributions for high energy electron beams with high accuracy. However, there are limitations for low energy electron beams. This work aims to improve the accuracy of the dose calculation using eMC for 4 and 6 MeV electron beams of Varian linear accelerators. Improvements implemented into the eMC include (1) improved determination of the initial electron energy spectrum by increased resolution of mono-energetic depth dose curves used during beam configuration; (2) inclusion of all the scrapers of the applicator in the beam model; (3) reduction of the maximum size of the sphere to be selected within the macro MC transport when the energy of the incident electron is below certain thresholds. The impact of these changes in eMC is investigated by comparing calculated dose distributions for 4 and 6 MeV electron beams at source to surface distance (SSD) of 100 and 110 cm with applicators ranging from 6 x 6 to 25 x 25 cm(2) of a Varian Clinac 2300C/D with the corresponding measurements. Dose differences between calculated and measured absolute depth dose curves are reduced from 6% to less than 1.5% for both energies and all applicators considered at SSD of 100 cm. Using the original eMC implementation, absolute dose profiles at depths of 1 cm, d(max) and R50 in water lead to dose differences of up to 8% for applicators larger than 15 x 15 cm(2) at SSD 100 cm. Those differences are now reduced to less than 2% for all dose profiles investigated when the improved version of eMC is used. At SSD of 110 cm the dose difference for the original eMC version is even more pronounced and can be larger than 10%. Those differences are reduced to within 2% or 2 mm with the improved version of eMC. In this work several enhancements were made in the eMC algorithm leading to significant improvements in the accuracy of the dose calculation for 4 and 6 MeV electron beams of Varian linear accelerators.

The Swiss IMRT dosimetry intercomparison using a thorax phantom

In 2008, a national intensity modulated radiation therapy (IMRT) dosimetry intercomparison was carried out for all 23 radiation oncology institutions in Switzerland. It was the aim to check the treatment chain focused on the planning, dose calculation, and irradiation process.

Automated detection of fiducial screws from CT/DVT volume data for image-guided ENT surgery

This paper presents an automated solution for precise detection of fiducial screws from three-dimensional (3D) Computerized Tomography (CT)/Digital Volume Tomography (DVT) data for image-guided ENT surgery. Unlike previously published solutions, we regard the detection of the fiducial screws from the CT/DVT volume data as a pose estimation problem. We thus developed a model-based solution. Starting from a user-supplied initialization, our solution detects the fiducial screws by iteratively matching a computer aided design (CAD) model of the fiducial screw to features extracted from the CT/DVT data. We validated our solution on one conventional CT dataset and on five DVT volume datasets, resulting in a total detection of 24 fiducial screws. Our experimental results indicate that the proposed solution achieves much higher reproducibility and precision than the manual detection. Further comparison shows that the proposed solution produces better results on the DVT dataset than on the conventional CT dataset.

Integrating optical fiber force sensors into microforceps for ORL microsurgery

The delicate anatomy of the ear require surgeons to use great care when operating on its internal structures. One example for such an intervention is the stapedectomy, where a small crook shaped piston is placed in the oval window of the cochlea and connected to the incus through crimping thus bypassing the diseased stapes. Performing the crimp process with the correct force is necessary since loose crimps poorly transmit sound whereas tight crimps will eventually result in necrosis of the incus. Clinically, demand is high to reproducibly conduct the crimp process through a precise force measurement. For this reason, we have developed a fiber Bragg grating (FBG) integrated microforceps for use in such interventions. This device was calibrated, and tested in cadaver preparations. With this instrument we were able to measure for the first time forces involved in crimping a stapes prosthesis to the incus. We also discuss a method of attaching and actuating such forceps in conjunction with a robot currently under development in our group. Each component of this system can be used separately or combined to improve surgical accuracy, confidence and outcome.

Low-dose pulmonary CT angiography: reduced radiation exposure and iodine load at low tube kilovoltage

CT pulmonary angiography is the currently accepted standard in ruling out acute pulmonary embolism. Issues of radiation dose received by patients via CT have been extensively disputed by radiologists and reported by the media. In recent years there has been considerable research performed to find ways for reducing radiation exposure from CT. Herein, we will discuss specific measures that have been shown to be valuable for CT pulmonary angiography. The limitations and the potential benefits of reduced CT peak tube kilovoltage will be detailed as this method is capable of reducing both radiation exposure and iodine load to the patient simultaneously. We discuss some of the emerging tools, which will hopefully play a significant role in wider acceptance of low-dose CT pulmonary angiography protocols.

Multimodal target correction by local bone registration: a PET/CT evaluation

PET/CT guidance for percutaneous interventions allows biopsy of suspicious metabolically active bone lesions even when no morphological correlation is delineable in the CT images. Clinical use of PET/CT guidance with conventional step-by-step technique is time consuming and complicated especially in cases in which the target lesion is not shown in the CT image. Our recently developed multimodal instrument guidance system (IGS) for PET/CT improved this situation. Nevertheless, bone biopsies even with IGS have a trade-off between precision and intervention duration which is proportional to patient and personnel exposure to radiation. As image acquisition and reconstruction of PET may take up to 10 minutes, preferably only one time consuming combined PET/CT acquisition should be needed during an intervention. In case of required additional control images in order to check for possible patient movements/deformations, or to verify the final needle position in the target, only fast CT acquisitions should be performed. However, for precise instrument guidance accounting for patient movement and/or deformation without having a control PET image, it is essential to be able to transfer the position of the target as identified in the original PET/CT to a changed situation as shown in the control CT.

ECM versus ICP for point registration

Iterative Closest Point (ICP) is a widely exploited method for point registration that is based on binary point-to-point assignments, whereas the Expectation Conditional Maximization (ECM) algorithm tries to solve the problem of point registration within the framework of maximum likelihood with point-to-cluster matching. In this paper, by fulfilling the implementation of both algorithms as well as conducting experiments in a scenario where dozens of model points must be registered with thousands of observation points on a pelvis model, we investigated and compared the performance (e.g. accuracy and robustness) of both ICP and ECM for point registration in cases without noise and with Gaussian white noise. The experiment results reveal that the ECM method is much less sensitive to initialization and is able to achieve more consistent estimations of the transformation parameters than the ICP algorithm, since the latter easily sinks into local minima and leads to quite different registration results with respect to different initializations. Both algorithms can reach the high registration accuracy at the same level, however, the ICP method usually requires an appropriate initialization to converge globally. In the presence of Gaussian white noise, it is observed in experiments that ECM is less efficient but more robust than ICP.

Monte Carlo implementation, validation, and characterization of a 120 leaf MLC

Recently, the new high definition multileaf collimator (HD120 MLC) was commercialized by Varian Medical Systems providing high resolution in the center section of the treatment field. The aim of this work is to investigate the characteristics of the HD120 MLC using Monte Carlo (MC) methods.

A virtual source model for kilo-voltage cone beam CT: source characteristics and model validation

The purpose of this investigation was to study the source characteristics of a clinical kilo-voltage cone beam CT unit and to develop and validate a virtual source model that could be used for treatment planning purposes.

Skeletal muscle (1) H MRSI before and after prolonged exercise. I. muscle specific depletion of intramyocellular lipids

Aim of the study was to determine distribution and depletion patterns of intramyocellular lipids (IMCL) in leg muscles before and after two types of standardized endurance exercise. ¹H-magnetic resonance spectroscopic imaging was performed (1) in the thigh of eight-trained cyclists after exercising on an ergometer for 3 h at 52 ± 8% of maximal speed and (2) in the lower leg of eight-trained runners after exercising on a treadmill for 3 h at 49 ± 3% of maximal workload. Pre-exercise IMCL contents were reduced postexercise in 11 out of 13 investigated upper and lower leg muscles (P < 0.015 for all). A strong linear correlation with a slope of ∼0.5 between pre-exercise IMCL content and IMCL depletion was found. IMCL depletion differed strongly between muscles. Absolute and also relative IMCL reduction was significantly higher in muscles with predominantly slow fibers compared to those with fast fibers. Creatine levels and fiber orientation were stable and unchanged after exercise, while trimethyl-ammonium groups increased. This is presented in the accompanying paper. In conclusion, a systematic comparison of metabolic changes in cross sections of the upper and lower leg was performed. The results imply that pre-exercise IMCL levels determine the degree of IMCL depletion after exercise.

Skeletal muscle (1) H MRSI before and after prolonged exercise. II. visibility of free carnitine

Carnitine (Car) buffers excess acetyl-CoA through the formation of acetylCar (AcCar). AcCar's acetyl group (AG-AcCar) gives rise to a peak at 2.13 ppm in ¹H MR spectra of skeletal muscle, whereas the trimethylammonium (TMA) groups of both, AcCar and Car, are thought to contribute to the TMA peak at 3.23 ppm. Surprisingly, in previous studies both resonances, AG-AcCar and TMA, increased after exercise. The aim of this study was to assess if the exercise-related TMA increase correlated with AcCar production. Magnetic resonance spectroscopic imaging (pulse repetition time/echo time = 1200/35 ms) was performed before and after prolonged exercise in the lower leg and thigh of eight runners and eight cyclists, respectively. TMA and AG-AcCar increased after exercise (P < 0.001). TMA's increase correlated with the AG-AcCar increase (R² = 0.73, P < 0.001, lower leg; R² = 0.28, P < 0.001, thigh). The correlation of ΔTMA with ΔAG-AcCar suggests that the TMA increase is due to AcCar formation. As total Car (Car + AcCar) remains unchanged with exercise, these findings suggest that the contribution of free Car to the TMA peak is limited and, therefore, is partly invisible in muscle ¹H MR spectra. This indicates that the biochemically relevant cytosolic content of free Car is considerably lower than the overall concentration determined by radioisotopic assays, a potentially important result with respect to regulation of substrate oxidation.

A simple, economical, and effective portable paediatric mock circulatory system

Ventricular assist devices (VADs) and total artificial hearts have been in development for the last 50 years. Since their inception, simulators of the circulation with different degrees of complexity have been produced to test these devices in vitro. Currently, a new path has been taken with the extensive efforts to develop paediatric VADs, which require totally different design constraints. This paper presents the manufacturing details of an economical simulator of the systemic paediatric circulation. This simulator allows the insertion of a paediatric VAD, includes a pumping ventricle, and is adjustable within the paediatric range. Rather than focusing on complexity and physiological simulation, this simulator is designed to be simple and practical for rapid device testing. The simulator was instrumented with medical sensors and data were acquired under different conditions with and without the new PediaFlowTM paediatric VAD. The VAD was run at different impeller speeds while simulator settings such as vascular resistance and stroke volume were varied. The hydraulic performance of the VAD under pulsatile conditions could be characterized and the magnetic suspension could be tested via manipulations such as cannula clamping. This compact mock loop has proven to be valuable throughout the PediaFlow development process and has the advantage that it is uncomplicated and can be manufactured cheaply. It can be produced by several research groups and the results of different VADs can then be compared easily.

Quantitative mapping of T2 using partial spoiling

Fast quantitative MRI has become an important tool for biochemical characterization of tissue beyond conventional T1, T2, and T2*-weighted imaging. As a result, steady-state free precession (SSFP) techniques have attracted increased interest, and several methods have been developed for rapid quantification of relaxation times using steady-state free precession. In this work, a new and fast approach for T2 mapping is introduced based on partial RF spoiling of nonbalanced steady-state free precession. The new T2 mapping technique is evaluated and optimized from simulations, and in vivo results are presented for human brain at 1.5 T and for human articular cartilage at 3.0 T. The range of T2 for gray and white matter was from 60 msec (for the corpus callosum) to 100 msec (for cortical gray matter). For cartilage, spatial variation in T2 was observed between deep (34 msec) and superficial (48 msec) layers, as well as between tibial (33 msec), femoral, (54 msec) and patellar (43 msec) cartilage. Excellent correspondence between T2 values derived from partially spoiled SSFP scans and the ones found with a reference multicontrast spin-echo technique is observed, corroborating the accuracy of the new method for proper T2 mapping. Finally, the feasibility of a fast high-resolution quantitative partially spoiled SSFP T2 scan is demonstrated at 7.0 T for human patellar cartilage.