Table of contents

The development of dedicated small animal PET (positron emission tomography) scanners has led to significantly higher spatial resolution and comparable sensitivity to clinical scanners. However, it is not clear whether we are approaching the fundamental limit of spatial resolution. This work aims to understand what is currently limiting spatial resolution during data formation and collection and how to apply that knowledge to obtain the best possible resolution for small animal PET without sacrificing sensitivity. Monte Carlo simulations were performed of the interactions of a 511 keV photon in a variety of detector materials to evaluate the modulation transfer function of the materials. Positron range, non-colinearity and pixel size were modelled to determine the contribution of additional components of data formation and collection on the complete modulation transfer function of a PET system. These simulations are shown to predict the intrinsic detector resolution of current high resolution systems very well. They also show that current detectors are not limited by inter-crystal scatter. An intrinsic resolution of 0.5 mm can be achieved, but would require a detector with a pixel size of around 250 µm that can be read out unambiguously. It is shown that a range of different detector materials, both scintillators and semiconductors, can be used in these high-resolution detectors. While this design relies on thin (∼3 mm) pieces of material, stacks of the material are shown to simultaneously provide spatial resolution near 0.5 mm and 60% efficiency. This work has shown that detectors with significantly better resolution and sensitivity can be developed for small animal PET applications.

Currently the only method to assess liver preservation injury is based on liver appearance and donor medical history. Previous work has shown that high-frequency ultrasound could detect ischemic cell death due to changes in cell morphology. In this study, we use high-frequency ultrasound integrated backscatter to assess liver damage in experimental models of liver ischemia. Ultimately, our goal is to predict organ suitability for transplantation using high-frequency imaging and spectral analysis techniques. To examine the effects of liver ischemia at different temperatures, livers from Wistar rats were surgically excised, immersed in phosphate buffer saline and stored at 4 and 20 °C for 24 h. To mimic organ preservation, livers were excised, flushed with University of Wisconsin (UW) solution and stored at 4 °C for 24 h. Preservation injury was simulated by either not flushing livers with UW solution or, before scanning, allowing livers to reach room temperature. Ultrasound images and corresponding radiofrequency data were collected over the ischemic period. No significant increase in integrated backscatter (∼2.5 dBr) was measured for the livers prepared using standard preservation conditions. For all other ischemia models, the integrated backscatter increased by 4–9 dBr demonstrating kinetics dependent on storage conditions. The results provide a possible framework for using high-frequency imaging to non-invasively assess liver preservation injury.

Voltage-sensitive dyes are an important tool in visualizing electrical activity in cardiac tissue. Until today, they have mainly been applied in cardiac electrophysiology to subsurface imaging. In the present study, we assess different imaging methods used in optical tomography with respect to their effectiveness in visualizing 3D cardiac activity. To achieve this goal, we simulate optical signals produced by excitation fronts initiated at different depths inside the myocardial wall and compare their properties for various imaging modes. Specifically, we consider scanning and broad-field illumination, including trans- and epi-illumination. We focus on the lateral optical resolution and signal intensity, as a function of the source depth. Optical diffusion theory is applied to derive a computationally efficient approximation of the point-spread function and to predict voltage-sensitive signals. Computations were performed both for fluorescent and absorptive voltage-sensitive dyes. Among all the above-mentioned methods, fluorescent coaxial scanning yields the best resolution (<2.5 mm) and gives the most information about the intramural cardiac activity.

A better understanding of gastric accommodation and gastric perception requires knowledge of regional gastric geometry and local gastric tension throughout the stomach. An analytic method based on medical imaging data was developed in this study to describe the three-dimensional (3D) rat stomach geometry and tension distribution. The surface principal radii of curvatures were simulated and the surface tension was calculated in the glandular and non-glandular region of the stomach at pressures from 0 Pa to 800 Pa. The radii of curvature and tension distribution in the stomach were non-homogeneous. The radii of curvature in the glandular stomach were larger than those in the non-glandular region at pressures less than 100 Pa (P < 0.001). When the pressure increased to more than 200 Pa, the radii of curvature in the non-glandular stomach was larger than in the glandular stomach (P < 0.05). The curvature and tension distribution mapping using medical imaging technology and 3D models can be used to characterize and distinguish the physical behaviour in separate regions of the stomach.

Light transmission data collected around an object show large variation with source–detector separation owing to the presence of single or multiple inhomogeneous regions in the object. This variation in the measured intensity is made use of to reconstruct regions of the inhomogeneous inclusions. In addition, it is possible to select a set of data from the above which is most likely least affected by the presence of the inhomogeneity, and estimate reasonably accurately the background optical properties from it. The reconstructed region is found to always contain the inhomogeneity and is of size approximately 140% by area of the inhomogeneity. With the regions to be reconstructed a priori known, a model-based iterative reconstruction procedure for reconstructing the optical properties of the region converged five times faster than without such information. It is also shown that whereas for the full object, a view-based propagation–backpropagation reconstruction procedure failed to converge, owing to large underdeterminacy of the problem, a smaller problem attempting to reconstruct a priori specified regions of interest converged and did so faster than a non-view-based approach for similar objects. Reconstruction results are presented from simulated transmitted intensity data from the following objects with regions of inhomogeneity in both absorption and scattering: (i) single centrally located inhomogeneity, (ii) two off-centred inhomogeneous regions of equal size and contrast (iii) two off-centred inhomogeneous regions of unequal size and equal contrast and (iv) two off-centred inhomogeneous regions of unequal size and contrast. Whereas the model-based iterative image reconstruction procedure gave good convergence in the first, second and third cases, in the fourth case the reconstructions failed to recover the exact numerical value of the optical properties in the higher contrast region.

Radiation therapy treatment planning is based on the calculation of the absorbed dose in the patient domain. For exact dose calculations, the solution of three coupled Boltzmann transport equations (BTEs) is needed to cover the transport of photons, electrons and positrons. In many situations, however, two coupled systems for photons and electrons are enough. The use of numerical methods in finding the exact solution of the unknown particle fluxes is necessary. In the stationary case, the BTE has six variables, three spatial, two directional and one energy variable. In this paper, we describe an approach in which the finite element method (FEM) is used to solve the six-dimensional problem. For the coupled photon–electron system, the variational formulation and the existence and uniqueness of the solution are derived. We simulate the solution of two coupled BTEs describing the travelling of photons and electrons in two spatial dimensions. The results are compared to Monte Carlo calculations with good agreement.

The contribution of high linear energy transfer (L) charged particles to dosimetric and microdosimetric characteristics in a clinical proton beam was experimentally studied using an ionization chamber and track etched detectors. The particles mentioned are produced by proton nuclear interactions; at the Bragg peak region slowed down protons also contribute in the L region above several keV µm⁻¹. Due to these particles the biological weighted effective dose (BWED) of the beam changes with depth. The spectra of particles with L above 7 keV µm⁻¹ were established by means of track etched detectors, which permitted us to determine their contribution to dosimetric and microdosimetric characteristics of clinical proton beams. The studies were realized in the clinical proton beam of the JINR Dubna Phasotron, with a primary energy of 205 MeV. The relative contribution to the absorbed dose of the particles with L above 7 keV µm⁻¹ increases from several per cent at the beam entrance to several tens of per cent at the Bragg peak region. The relative biological weighted efficiency (RBWE) for radiotherapy has been calculated using a biological weighting function. It increases with depth from 1.02 at the beam entrance to about 1.25 at the Bragg peak region.

The purpose of this study was to investigate the dosimetric characteristics (energy dependence, linearity, fading, reproducibility, etc) of MOSFET detectors for in vivo dosimetry in the kV x-ray range. The experience of MOSFET in vivo dosimetry in a pre-clinical study using the Alderson phantom and in clinical practice is also reported. All measurements were performed with a Gulmay D3300 kV unit and TN-502RDI MOSFET detectors. For the determination of correction factors different solid phantoms and a calibrated Farmer-type chamber were used. The MOSFET signal was linear with applied dose in the range from 0.2 to 2 Gy for all energies. Due to fading it is recommended to read the MOSFET signal during the first 15 min after irradiation. For long time intervals between irradiation and readout the fading can vary largely with the detector. The temperature dependence of the detector signal was small (0.3% °C⁻¹) in the temperature range between 22 and 40 °C. The variation of the measuring signal with beam incidence amounts to ±5% and should be considered in clinical applications. Finally, for entrance dose measurements energy-dependent calibration factors, correction factors for field size and irradiated cable length were applied. The overall accuracy, for all measurements, was dominated by reproducibility as a function of applied dose. During the pre-clinical in vivo study, the agreement between MOSFET and TLD measurements was well within 3%. The results of MOSFET measurements, to determine the dosimetric characteristics as well as clinical applications, showed that MOSFET detectors are suitable for in vivo dosimetry in the kV range. However, some energy-dependent dosimetry effects need to be considered and corrected for. Due to reproducibility effects at low dose levels accurate in vivo measurements are only possible if the applied dose is equal to or larger than 2 Gy.

A novel alpha-particle irradiator has recently been developed that provides the ability to characterize cell response. The irradiator is comprised of a collimated, planar alpha-particle source which, from below, irradiates cells cultured on a track-etch material. Cells are imaged using phase-contrast microscopy before and following irradiation to obtain geometric information and survival rates; these can be used with data from alpha-particle track images to assess cell response. A key step in this process is determining cell location within the pre-irradiation images. Although this can be done completely by a human observer, the number of images requiring analysis makes the process time-consuming and tedious. To reduce the potential human error and decrease user interaction time, a semi-automated, computer-aided method of cell detection has been developed. The method employs a two-level adaptive thresholding technique to obtain size and position information about potential cell cytoplasms and nuclei. Proximity and geometry-based thresholds are then used to mark structures as cells. False-positive detections from the automated algorithm are due mostly to imperfections in the track-etch background, camera effects and cellular residue. To correct for these, a human observer reviews all detected structures, discarding false positives. When analysing two randomly selected cell dish image databases, the semi-automated method detected 92–94% of all cells and 94–97% of cells with a well-defined cytoplasm and nucleus while reducing human workload by 32–83%.

Independent treatment verification for high dose rate (HDR) brachytherapy is needed to ensure that the treatment proceeds as prescribed. In this paper, we investigate the feasibility of a proposed real-time source position verification process. This process provides immediate confirmation of the source position during the treatment, so that the treatment can be aborted and modified if necessary. We show that an array of dosimeters placed on the patient's skin can independently verify the position in three dimensions. This verification was demonstrated by using a diamond detector placed in several locations on the surface of an anthropomorphic phantom. A mathematical algorithm was constructed to estimate the location of the source given a measured data set in the presence of tissue heterogeneity. The accuracy of the source localization was found to increase with the number of detectors used to compute the estimation of the source position. The resolution to which the 12 detectors can identify the location of the source was within 3 mm.

Several Monte Carlo simulators are currently available for positron emission tomography (PET). Because each code has been described in a different way, it is difficult to know which one is best suited to a specific application. To help clarify the capabilities and accuracy of different codes dedicated to PET simulations, we propose a uniform description of the code features. This description specifies features pertaining to the models used for simulating the physics of PET and for describing a PET acquisition, to the acceleration strategies and to the technical characteristics of the code implementation. To assess the code accuracy, we suggest validation procedures based on NEMA phantoms involving standard physical parameters and simulation of a complex activity distribution. A test characterizing the statistical properties of detected coincidences is also described. The proposed code description and validation procedures are illustrated by considering the SimSET and PET-EGS codes. These codes differ in many features, including models for randoms and dead time, and source description. Despite these differences, both codes yielded data with properties close to those of real data. Depending on the intended application, one code might be preferred however. Indeed, only PET-EGS allows for accurate modelling of count rates while SimSET is more computationally efficient. The proposed code description and validation procedures might help determine which code is most appropriate for a specific application.

Ultra-wideband (UWB) electromagnetic pulses of nanosecond duration, or nanopulses, are of considerable interest to the communications industry and are being explored for various applications in biotechnology and medicine. The propagation of a nanopulse through biological matter has been computed using the finite difference-time domain (FDTD) method. The approach required the reparametrization of existing Cole–Cole model-based descriptions of dielectric properties of biological matter in terms of the Debye model without loss of accuracy. Several tissue types have been considered. Results show that the electromagnetic field inside biological tissue depends on incident pulse rise time and width. Rise time dominates pulse behaviour inside tissue as conductivity increases. It has also been found that the amount of energy deposited by 20 kV m⁻¹ nanopulses is insufficient to change the temperature of the exposed material for pulse repetition rates of 1 MHz or less, consistent with recent experimental results.

The purpose of this study was the characterization of local magnetic susceptibility deviations by spectral analysis of their induced magnetic field inhomogeneities. Magnetic resonance spectra and related signal decay curves of local susceptibility deviations were simulated for different volume fractions and compositions of the object within the VOI. The size or composition of the object was varied at constant volume fraction, constant object size, or at constant 'magnetic strength' (defined as the product of the volume and the volume susceptibility of the object). Experimental spectra were acquired for individual metal spherical particles and a spherical air cavity. Where possible, spectra were used to characterize objects in terms of volume and composition. By simulations, a numerical relation was determined between the spectral broadening and the object's volume and composition. Comparison of spectra for various spherical objects showed the possibility of characterization with respect to size and composition. Experimental results confirmed the numerical results to a large extent, although the characterization was compromised by background signal decay, low volume fractions and limitations in signal-to-noise. In conclusion, spectral description of the field inhomogeneities related to small objects allows characterization of such objects with respect to size and composition. Practical applicability of the simulation results depends on background signal decay and volume fraction of the object.

Compliance testing of electronic article surveillance (EAS) devices requires that induced current densities in central nervous system (CNS) tissues, i.e. brain and the spinal cord, be less than the prescribed safety limits. Even though ferromagnetic cores are mostly used for activation/deactivation of embedded magnetic tags, assumed equivalent air-core coils with guessed increased number of ampere turns have always been used to calculate the magnetic fields for the proximal region to which a customer is exposed. We show that at low frequencies up to several kilohertz, duality of electric and magnetic circuits may be exploited such that the shaped high reluctance core is modelled as though it was a higher conductivity electric circuit of the corresponding shape. The proposed procedure is tested by examples of two magnetic cores typical of countertop activation/deactivation devices. The equivalent exposure magnetic fields obtained from the dual electric fields are shown to be in excellent agreement (within ±5%) with those measured for these ferromagnetic EAS devices. The previously proposed impedance method is then used to calculate the induced current densities for a 1.974 × 1.974 × 2.93 mm resolution anatomic model of a human. For the two considered EAS systems using excitation currents of 5000 A turns at 200 Hz, the maximum 1 cm² area-averaged induced current densities in the CNS tissues are calculated and found to be less than the ICNIRP safety limits.

Cryogen spray cooling (CSC) is a heat extraction process that protects the epidermis from thermal damage during dermatologic laser surgery. The objective of the present work is to investigate radial and temporal variations in the heat transferred through the surface of a skin phantom during CSC. A fast-response thermal sensor is used to measure surface temperatures every 1 mm across a 16 mm diameter of the sprayed surface of the phantom. An analytical expression based on Fourier's law and Duhamel's theorem is used to compute surface heat fluxes from temperature measurements. Results show that radial and temporal variations of the boundary conditions have a strong influence on the homogeneity of heat extraction from the skin phantom. However, there is a subregion of uniform cooling whose size is time dependent. It is also observed that the surface heat flux undergoes a marked dynamic variation, with a maximum heat flux occurring at the centre of the sprayed surface early in the spurt followed by a quick decrease. The study shows that radial and temporal variations of boundary conditions must be taken into account and ideally controlled to guarantee uniform protection during CSC of human skin.