Tomographic 3D reconstruction through linear inversion for breast microwave imaging

In the present thesis we address the problem of detecting and localizing a small spherical target with characteristic electrical properties inside a volume of cylindrical shape, representing female breast, with MWI. One of the main works of this project is to properly extend the existing linear inversion algorithm from planar slice to volume reconstruction; results obtained, under the same conditions and experimental setup are reported for the two different approaches. Preliminar comparison and performance analysis of the reconstruction algorithms is performed via numerical simulations in a software-created environment: a single dipole antenna is used for illuminating the virtual breast phantom from different positions and, for each position, the corresponding scattered field value is registered. Collected data are then exploited in order to reconstruct the investigation domain, along with the scatterer position, in the form of image called pseudospectrum. During this process the tumor is modeled as a dielectric sphere of small radius and, for electromagnetic scattering purposes, it's treated as a point-like source. To improve the performance of reconstruction technique, we repeat the acquisition for a number of frequencies in a given range: the different pseudospectra, reconstructed from single frequency data, are incoherently combined with MUltiple SIgnal Classification (MUSIC) method which returns an overall enhanced image. We exploit multi-frequency approach to test the performance of 3D linear inversion reconstruction algorithm while varying the source position inside the phantom and the height of antenna plane. Analysis results and reconstructed images are then reported. Finally, we perform 3D reconstruction from experimental data gathered with the acquisition system in the microwave laboratory at DIFA, University of Bologna for a recently developed breast-phantom prototype; obtained pseudospectrum and performance analysis for the real model are reported.

Determination of Pelvic Orientation from Ultrasound Images Using Patch-SSMs and a Hierarchical Speed of Sound Compensation Strategy

In the field of computer assisted orthopedic surgery (CAOS) the anterior pelvic plane (APP) is a common concept to determine the pelvic orientation by digitizing distinct pelvic landmarks. As percutaneous palpation is - especially for obese patients - known to be error-prone, B-mode ultrasound (US) imaging could provide an alternative means. Several concepts of using ultrasound imaging to determine the APP landmarks have been introduced. In this paper we present a novel technique, which uses local patch statistical shape models (SSMs) and a hierarchical speed of sound compensation strategy for an accurate determination of the APP. These patches are independently matched and instantiated with respect to associated point clouds derived from the acquired ultrasound images. Potential inaccuracies due to the assumption of a constant speed of sound are compensated by an extended reconstruction scheme. We validated our method with in-vitro studies using a plastic bone covered with a soft-tissue simulation phantom and with a preliminary cadaver trial.

Comparison of three optical tracking systems in a complex navigation scenario

Three-dimensional rotational X-ray imaging with the SIREMOBIL Iso-C3D (Siemens AG, Medical Solutions, Erlangen, Germany) has become a well-established intra-operative imaging modality. In combination with a tracking system, the Iso-C3D provides inherently registered image volumes ready for direct navigation. This is achieved by means of a pre-calibration procedure. The aim of this study was to investigate the influence of the tracking system used on the overall navigation accuracy of direct Iso-C3D navigation. Three models of tracking system were used in the study: Two Optotrak 3020s, a Polaris P4 and a Polaris Spectra system, with both Polaris systems being in the passive operation mode. The evaluation was carried out at two different sites using two Iso-C3D devices. To measure the navigation accuracy, a number of phantom experiments were conducted using an acrylic phantom equipped with titanium spheres. After scanning, a special pointer was used to pinpoint these markers. The difference between the digitized and navigated positions served as the accuracy measure. Up to 20 phantom scans were performed for each tracking system. The average accuracy measured was 0.86 mm and 0.96 mm for the two Optotrak 3020 systems, 1.15 mm for the Polaris P4, and 1.04 mm for the Polaris Spectra system. For the Polaris systems a higher maximal error was found, but all three systems yielded similar minimal errors. On average, all tracking systems used in this study could deliver similar navigation accuracy. The passive Polaris system showed ? as expected ? higher maximal errors; however, depending on the application constraints, this might be negligible.

Effect of beam hardening on arterial enhancement in thoracoabdominal CT angiography with increasing patient size: an in vitro and in vivo study

To assess the effect of beam hardening on arterial enhancement in thoracoabdominal computed tomographic (CT) angiography in various body sizes in a phantom and in a clinical study.

Copper filtration in pediatric digital X-ray imaging: its impact on image quality and dose

The effect of copper (Cu) filtration on image quality and dose in different digital X-ray systems was investigated. Two computed radiography systems and one digital radiography detector were used. Three different polymethylmethacrylate blocks simulated the pediatric body. The effect of Cu filters of 0.1, 0.2, and 0.3 mm thickness on the entrance surface dose (ESD) and the corresponding effective doses (EDs) were measured at tube voltages of 60, 66, and 73 kV. Image quality was evaluated in a contrast-detail phantom with an automated analyzer software. Cu filters of 0.1, 0.2, and 0.3 mm thickness decreased the ESD by 25-32%, 32-39%, and 40-44%, respectively, the ranges depending on the respective tube voltages. There was no consistent decline in image quality due to increasing Cu filtration. The estimated ED of anterior-posterior (AP) chest projections was reduced by up to 23%. No relevant reduction in the ED was noted in AP radiographs of the abdomen and pelvis or in posterior-anterior radiographs of the chest. Cu filtration reduces the ESD, but generally does not reduce the effective dose. Cu filters can help protect radiosensitive superficial organs, such as the mammary glands in AP chest projections.

Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.

Validation of the Swiss Monte Carlo Plan for a static and dynamic 6 MV photon beam

Monte Carlo (MC) based dose calculations can compute dose distributions with an accuracy surpassing that of conventional algorithms used in radiotherapy, especially in regions of tissue inhomogeneities and surface discontinuities. The Swiss Monte Carlo Plan (SMCP) is a GUI-based framework for photon MC treatment planning (MCTP) interfaced to the Eclipse treatment planning system (TPS). As for any dose calculation algorithm, also the MCTP needs to be commissioned and validated before using the algorithm for clinical cases. Aim of this study is the investigation of a 6 MV beam for clinical situations within the framework of the SMCP. In this respect, all parts i.e. open fields and all the clinically available beam modifiers have to be configured so that the calculated dose distributions match the corresponding measurements. Dose distributions for the 6 MV beam were simulated in a water phantom using a phase space source above the beam modifiers. The VMC++ code was used for the radiation transport through the beam modifiers (jaws, wedges, block and multileaf collimator (MLC)) as well as for the calculation of the dose distributions within the phantom. The voxel size of the dose distributions was 2mm in all directions. The statistical uncertainty of the calculated dose distributions was below 0.4%. Simulated depth dose curves and dose profiles in terms of [Gy/MU] for static and dynamic fields were compared with the corresponding measurements using dose difference and γ analysis. For the dose difference criterion of ±1% of D(max) and the distance to agreement criterion of ±1 mm, the γ analysis showed an excellent agreement between measurements and simulations for all static open and MLC fields. The tuning of the density and the thickness for all hard wedges lead to an agreement with the corresponding measurements within 1% or 1mm. Similar results have been achieved for the block. For the validation of the tuned hard wedges, a very good agreement between calculated and measured dose distributions was achieved using a 1%/1mm criteria for the γ analysis. The calculated dose distributions of the enhanced dynamic wedges (10°, 15°, 20°, 25°, 30°, 45° and 60°) met the criteria of 1%/1mm when compared with the measurements for all situations considered. For the IMRT fields all compared measured dose values agreed with the calculated dose values within a 2% dose difference or within 1 mm distance. The SMCP has been successfully validated for a static and dynamic 6 MV photon beam, thus resulting in accurate dose calculations suitable for applications in clinical cases.

How reliable are Hounsfield-unit measurements in forensic radiology?

OBJECTIVE: To assess the reliability of computed tomography (CT) numbers, also known as Hounsfield-units (HU) in the differentiation and identification of forensically relevant materials and to provide instructions to improve the reproducibility of HU measurements in daily forensic practice. MATERIALS AND METHODS: We scanned a phantom containing non-organic materials (glass, rocks and metals) on three different CT scanners with standardized parameters. The t-test was used to assess the influence of the scanner, the size and shape of different types of regions-of-interest (ROI), the composition and shape of the object, and the reader performance on HU measurements. Intra-class correlation coefficient was used to assess intra- and inter-reader reliability. RESULTS: HU values did not change significantly as a function of ROI-shape or -size (p>0.05). Intra-reader reliability reached ICC values >0.929 (p<0.001). Inter-reader reliability was also excellent with an ICC of 0.994 (p<0.001). Four of seven objects yielded significantly different CT numbers at different levels within the object (p<0.05). In 6/7 objects the HU changed significantly from CT scanner to CT scanner (p<0.05). CONCLUSION: Reproducible CT number measurements can be achieved through correct ROI-placement and repeat measurements within the object of interest. However, HU may differ from CT-scanner to CT-scanner. In order to obtain comparable CT numbers we suggest that a dedicated Forensic Reference Phantom be developed.

A multi-center study: Intra-scan and inter-scan variability of diffusion spectrum imaging

The objective of this study was to investigate whether it is possible to pool together diffusion spectrum imaging data from four different scanners, located at three different sites. Two of the scanners had identical configuration whereas two did not. To measure the variability, we extracted three scalar maps (ADC, FA and GFA) from the DSI and utilized a region and a tract-based analysis. Additionally, a phantom study was performed to rule out some potential factors arising from the scanner performance in case some systematic bias occurred in the subject study. This work was split into three experiments: intra-scanner reproducibility, reproducibility with twin-scanner settings and reproducibility with other configurations. Overall for the intra-scanner and twin-scanner experiments, the region-based analysis coefficient of variation (CV) was in a range of 1%-4.2% and below 3% for almost every bundle for the tract-based analysis. The uncinate fasciculus showed the worst reproducibility, especially for FA and GFA values (CV 3.7-6%). For the GFA and FA maps, an ICC value of 0.7 and above is observed in almost all the regions/tracts. Looking at the last experiment, it was found that there is a very high similarity of the outcomes from the two scanners with identical setting. However, this was not the case for the two other imagers. Given the fact that the overall variation in our study is low for the imagers with identical settings, our findings support the feasibility of cross-site pooling of DSI data from identical scanners.

A 20-channel receive-only mouse array coil for a 3 T clinical MRI system

A 20-channel phased-array coil for MRI of mice has been designed, constructed, and validated with bench measurements and high-resolution accelerated imaging. The technical challenges of designing a small, high density array have been overcome using individual small-diameter coil elements arranged on a cylinder in a hexagonal overlapping design with adjacent low impedance preamplifiers to further decouple the array elements. Signal-to-noise ratio (SNR) and noise amplification in accelerated imaging were simulated and quantitatively evaluated in phantoms and in vivo mouse images. Comparison between the 20-channel mouse array and a length-matched quadrature driven small animal birdcage coil showed an SNR increase at the periphery and in the center of the phantom of 3- and 1.3-fold, respectively. Comparison with a shorter but SNR-optimized birdcage coil (aspect ratio 1:1 and only half mouse coverage) showed an SNR gain of twofold at the edge of the phantom and similar SNR in the center. G-factor measurements indicate that the coil is well suited to acquire highly accelerated images.

The low-energy β− and electron emitter 161Tb as an alternative to 177Lu for targeted radionuclide therapy

The low-energy β− emitter 161Tb is very similar to 177Lu with respect to half-life, beta energy and chemical properties. However, 161Tb also emits a significant amount of conversion and Auger electrons. Greater therapeutic effect can therefore be expected in comparison to 177Lu. It also emits low-energy photons that are useful for gamma camera imaging. The 160Gd(n,γ)161Gd→161Tb production route was used to produce 161Tb by neutron irradiation of massive 160Gd targets (up to 40 mg) in nuclear reactors. A semiautomated procedure based on cation exchange chromatography was developed and applied to isolate no carrier added (n.c.a.) 161Tb from the bulk of the 160Gd target and from its stable decay product 161Dy. 161Tb was used for radiolabeling DOTA-Tyr3-octreotate; the radiolabeling profile was compared to the commercially available n.c.a. 177Lu. A 161Tb Derenzo phantom was imaged using a small-animal single-photon emission computed tomography camera. Up to 15 GBq of 161Tb was produced by long-term irradiation of Gd targets. Using a cation exchange resin, we obtained 80%–90% of the available 161Tb with high specific activity, radionuclide and chemical purity and in quantities sufficient for therapeutic applications. The 161Tb obtained was of the quality required to prepare 161Tb–DOTA-Tyr3-octreotate. We were able to produce 161Tb in n.c.a. form by irradiating highly enriched 160Gd targets; it can be obtained in the quantity and quality required for the preparation of 161Tb-labeled therapeutic agents.

Topographic and age-dependent distribution of subchondral bone density in the elbow joints of clinically normal dogs

OBJECTIVE: To investigate topographic and age-dependent adaptation of subchondral bone density in the elbow joints of healthy dogs by means of computed tomographic osteoabsorptiometry (CTOAM). Animals-42 elbow joints of 29 clinically normal dogs of various breeds and ages. PROCEDURES: Subchondral bone densities of the humeral, radial, and ulnar joint surfaces of the elbow relative to a water-hydroxyapatite phantom were assessed by means of CTOAM. Distribution patterns in juvenile, adult, and geriatric dogs (age, < 1 year, 1 to 8 years, and > 8 years, respectively) were determined and compared within and among groups. RESULTS: An area of increased subchondral bone density was detected in the humerus distomedially and cranially on the trochlea and in the olecranon fossa. The ulna had maximum bone densities on the anconeal and medial coronoid processes. Increased bone density was detected in the craniomedial region of the joint surface of the radius. A significant age-dependent increase in subchondral bone density was revealed in elbow joint surfaces of the radius, ulna, and humerus. Mean subchondral bone density of the radius was significantly less than that of the ulna in paired comparisons for all dogs combined and in adult and geriatric, but not juvenile, dog groups. CONCLUSIONS AND CLINICAL RELEVANCE: An age-dependent increase in subchondral bone density at the elbow joint was revealed. Maximal relative subchondral bone densities were detected consistently at the medial coronoid process and central aspect of the humeral trochlea, regions that are commonly affected in dogs with elbow dysplasia.

Endoscopic Image Overlay for the Targeting of Hidden Anatomy in Laparoscopic Visceral Surgery

Limitations associated with the visual information provided to surgeons during laparoscopic surgery increases the difficulty of procedures and thus, reduces clinical indications and increases training time. This work presents a novel augmented reality visualization approach that aims to improve visual data supplied for the targeting of non visible anatomical structures in laparoscopic visceral surgery. The approach aims to facilitate the localisation of hidden structures with minimal damage to surrounding structures and with minimal training requirements. The proposed augmented reality visualization approach incorporates endoscopic images overlaid with virtual 3D models of underlying critical structures in addition to targeting and depth information pertaining to targeted structures. Image overlay was achieved through the implementation of camera calibration techniques and integration of the optically tracked endoscope into an existing image guidance system for liver surgery. The approach was validated in accuracy, clinical integration and targeting experiments. Accuracy of the overlay was found to have a mean value of 3.5 mm ± 1.9 mm and 92.7% of targets within a liver phantom were successfully located laparoscopically by non trained subjects using the approach.

Assessment of spline-based 2D-3D registration for image-guided spine surgery

2D-3D registration of pre-operative 3D volumetric data with a series of calibrated and undistorted intra-operative 2D projection images has shown great potential in CT-based surgical navigation because it obviates the invasive procedure of the conventional registration methods. In this study, a recently introduced spline-based multi-resolution 2D-3D image registration algorithm has been adapted together with a novel least-squares normalized pattern intensity (LSNPI) similarity measure for image guided minimally invasive spine surgery. A phantom and a cadaver together with their respective ground truths were specially designed to experimentally assess possible factors that may affect the robustness, accuracy, or efficiency of the registration. Our experiments have shown that it is feasible for the assessed 2D-3D registration algorithm to achieve sub-millimeter accuracy in a realistic setup in less than one minute.

Computer-assisted, fluoroscopy-based ventral spondylodesis of thoracolumbar fractures

OBJECTIVE: To design and evaluate a novel computer-assisted, fluoroscopy-based planning and navigation system for minimally invasive ventral spondylodesis of thoracolumbar fractures. MATERIALS AND METHODS: Instruments and an image intensifier are tracked with the SurgiGATE navigation system (Praxim-Medivision). Two fluoroscopic images, one acquired from anterior-posterior (AP) direction and the other from lateral-medial (LM) direction, are used for the complete procedure of planning and navigation. Both of them are calibrated with a custom-made software to recover their projection geometry and to co-register them to a common patient reference coordinate system, which is established by attaching an opto-electronically trackable dynamic reference base (DRB) on the operated vertebra. A bi-planar landmark reconstruction method is used to acquire deep-seated anatomical landmarks such that an intraoperative planning of graft bed can be interactively done. Finally, surgical actions such as the placement of the stabilization devices and the formation of the graft bed using a custom-made chisel are visualized to the surgeon by superimposing virtual instrument representations onto the acquired images. The distance between the instrument tip and each wall of the planned graft bed are calculated on the fly and presented to the surgeon so that the surgeon could formalize the graft bed exactly according to his/her plan. RESULTS: Laboratory studies on phantom and on 27 plastic vertebras demonstrate the high precision of the proposed navigation system. Compared with CT-based measurement, a mean error of 1.0 mm with a standard deviation of 0.1 mm was found. CONCLUSIONS: The proposed computer assisted, fluoroscopy-based planning and navigation system promises to increase the accuracy and reliability of minimally invasive ventral spondylodesis of thoracolumbar fractures.