Antj Alexis Kotte

Calculation of change in brain temperatures due to exposure to a mobile phone

In this study we evaluated for a realistic head model the 3D temperature rise induced by a mobile phone. This was done numerically with the consecutive use of an FDTD model to predict the absorbed electromagnetic power distribution, and a thermal model describing bioheat transfer both by conduction and by blood flow. We calculated a maximum rise in brain temperature of 0.11 degrees C for an antenna with an average emitted power of 0.25 W, the maximum value in common mobile phones, and indefinite exposure. Maximum temperature rise is at the skin. The power distributions were characterized by a maximum averaged SAR over an arbitrarily shaped 10 g volume of approximately 1.6 W kg(-1). Although these power distributions are not in compliance with all proposed safety standards, temperature rises are far too small to have lasting effects. We verified our simulations by measuring the skin temperature rise experimentally. Our simulation method can be instrumental in further development of safety standards.

Label Fusion in Atlas-Based Segmentation Using a Selective and Iterative Method for Performance Level Estimation (SIMPLE)

In a multi-atlas based segmentation procedure, propagated atlas segmentations must be combined in a label fusion process. Some current methods deal with this problem by using atlas selection to construct an atlas set either prior to or after registration. Other methods estimate the performance of propagated segmentations and use this performance as a weight in the label fusion process. This paper proposes a selective and iterative method for performance level estimation (SIMPLE), which combines both strategies in an iterative procedure. In subsequent iterations the method refines both the estimated performance and the set of selected atlases. For a dataset of 100 MR images of prostate cancer patients, we show that the results of SIMPLE are significantly better than those of several existing methods, including the STAPLE method and variants of weighted majority voting.

Integrating a MRI scanner with a 6 MV radiotherapy accelerator: dose deposition in a transverse magnetic field

Integrating magnetic resonance imaging (MRI) functionality with a radiotherapy accelerator can facilitate on-line, soft-tissue based, position verification. A technical feasibility study, in collaboration with Elekta Oncology Systems and Philips Medical Systems, led to the preliminary design specifications of a MRI accelerator. Basically the design is a 6 MV accelerator rotating around a 1.5 T MRI system. Several technical issues and the clinical rational are currently under investigation. The aim of this paper is to determine the impact of the transverse 1.5 T magnetic field on the dose deposition. Monte Carlo simulations were used to calculate the dose deposition kernel in the presence of 1.5 T. This kernel in turn was used to determine the dose deposition for larger fields. Also simulations and measurements were done in the presence of 1.1 T. The pencil beam dose deposition is asymmetric. For larger fields the asymmetry persists but decreases. For the latter the distance to dose maximum is reduced by approximately 5 mm, the penumbra is increased by approximately 1 mm, and the 50% isodose line is shifted approximately 1 mm. The dose deposition in the presence of 1.5 T is affected, but the effect can be taken into account in a conventional treatment planning procedure. The impact of the altered dose deposition for clinical IMRT treatments is the topic of further research.

A flexible algorithm for construction of 3-D vessel networks for use in thermal modeling

A new algorithm for the construction of artificial blood vessel networks is presented. The algorithm produces three-dimensional (3-D) geometrical representations of both arterial and venous networks. The key ingredient of the algorithm is a 3-D potential function defined in the tissue volume. This potential function controls the paths by which points are connected to existing vessels, thereby producing new vessel segments. The potential function has no physiological interpretation, but, by adjustment of parameters governing the potential, it is possible to produce networks that have physiologically meaningful geometrical properties. If desired, the veins can be generated counter current to the arteries. Furthermore, the potential function allows fashioning of the networks to the presence of bone or air cavities. The resulting networks can be used for thermal simulations of hyperthermia treatment.

SEGMENTATION OF THE PROSTATE IN MR IMAGES BY ATLAS MATCHING

Prostate cancer treatment by radiation therapy requires an accurate localisation of the prostate. For the treatment planning, primarily computed tomography (CT) images are used, but increasingly magnetic resonance (MR) images are added, because of their soft-tissue contrast. In current practice at our hospital, a manual delineation of the prostate is made, based on the CT and MR scans, which is a labour-intensive task. We propose an automatic segmentation method, based on non-rigid registration of a set of prelabelled MR atlas images. The algorithm consists of three stages. Firstly, the target image is nonrigidly registered with each atlas image, using mutual information as the similarity measure. After that, the best registered atlas images are selected by comparing the mutual information values after registration. Finally, the segmentation is obtained by averaging the selected deformed segmentations and thresholding the result. The method is evaluated on 22 images by calculating the overlap of automatic and manual segmentations. This results in a median Dice similarity coefficient of 0.82

Temperature simulations in tissue with a realistic computer generated vessel network

The practical use of a discrete vessel thermal model for hyperthermia treatment planning requires a number of choices with respect to the unknown part of the patients vasculature. This work presents a study of the thermal effects of blood flow in a simple tissue geometry with a detailed artificial vessel network. The simulations presented here demonstrate that an incomplete discrete description of the detailed network results in a better prediction of the temperature distribution than is obtained using the conventional bio-heatsink equation. Therefore, efforts to obtain information on the positions of the large vessels in an individual hyperthermia patient will be rewarded with a more accurate prediction of the temperature distribution.

Tests of the geometrical description of blood vessels in a thermal model using counter-current geometries

We have developed a thermal model, for use in hyperthermia treatment planning, in which blood vessels are described as geometrical objects; 3D curves with associated diameters. For the calculation of the heat exchange with the tissue an analytic result is used. To arrive at this result some assumptions were made. One of these assumptions is a cylindrically symmetric temperature distribution. In this paper the behaviour of the model is examined for counter-current vessel geometries for which this assumption is not valid. Counter-current vessel pairs intersecting a circular tissue slice are tested. For these 2D geometries vessel spacing, tissue radius and resolution are varied, as well as the position of the vessel pair with respect to the discretized tissue grid. The simulation results are evaluated by comparison of the different heat flow rates with analytical predictions. The tests show that for a fixed vessel configuration the accuracy is not a simple decreasing function of the voxel dimensions, but is also sensitive to the position of the configuration with respect to the discretized tissue grid.

Accuracy of geometrical modelling of heat transfer from tissue to blood vessels

We have developed a thermal model in which blood vessels are described as geometrical objects, 3D curves with associated diameters. Here the behaviour of the model is examined for low resolutions compared with the vessel diameter and for strongly curved vessels. The tests include a single straight vessel and vessels describing the path of a helix embedded in square tissue blocks. The tests show the excellent behaviour of our discrete vessel implementation.

Dose uniformity of ferromagnetic seed implants in tissue with discrete vasculature: a numerical study on the impact of seed characteristics and implantation techniques.

The results from simulations with a new three-dimensional treatment planning system for interstitial hyperthermia with ferromagnetic seeds are presented in this study. The thermal model incorporates discrete vessel structures as well as a heat sink and enhanced thermal conductivity. Both the discrete vessels and the ferroseeds are described parametrically in separate calculation spaces. This parametric description has the advantage of an arbitrary orientation of the structures within the tissue grid, easy manipulation of the structures and independence from the resolution of the tissue voxels (tissue calculation space). The power absorption of the self-regulating seeds is according to empirical data. The thermal effects of an unlimited number of thin layers surrounding the seed (coatings, catheters) can be modelled. The initial calculations have been performed for an array of 12 identical ferromagnetic seeds in a tissue volume with a computer generated artificial vessel network spanning four vessel generations in both the arterial and venous tree. The heterogeneously distributed large isolated vessels impair the temperature distribution significantly, indicating the limited accuracy of continuum models. Simulations with different types of ferromagnetic seeds have confirmed that the efforts of previous studies to optimize the self-regulating temperature control and the implantation techniques of the ferroseeds will improve the homogeneity of the temperature distribution in the target volume. Multifilament seeds implanted in brachytherapy needles and tubular seeds appear to be the most favourable configurations. The division of long seeds into shorter segments with the appropriate Curie temperature will further improve the homogeneity of the temperature distribution without increasing the average temperature in the volume of interest. Given the proper thermal tissue data, the model presented in this study will prove to be a useful tool in making choices for the implant geometry, seed spacing and Curie temperature.

Registration of structurally dissimilar images in MRI-based brachytherapy.

A serious challenge in image registration is the accurate alignment of two images in which a certain structure is present in only one of the two. Such topological changes are problematic for conventional non-rigid registration algorithms. We propose to incorporate in a conventional free-form registration framework a geometrical penalty term that minimizes the volume of the missing structure in one image. We demonstrate our method on cervical MR images for brachytherapy. The intrapatient registration problem involves one image in which a therapy applicator is present and one in which it is not. By including the penalty term, a substantial improvement in the surface distance to the gold standard anatomical position and the residual volume of the applicator void are obtained. Registration of neighboring structures, i.e. the rectum and the bladder is generally improved as well, albeit to a lesser degree.