Accurate analytical expressions for collimator resolution and sensitivity are important tools in the optimization of SPECT systems. However, presently known expressions for the sensitivity of converging collimators either diverge near the focal point or focal line(s), or are only valid on the collimator axis. As a result, these expressions are unsuitable to calculate volumetric sensitivity for e.g. short-focal length collimators that focus inside the object to enhance sensitivity. To also enable collimator optimization for these geometries, we here present non-diverging sensitivity formulas for astigmatic, cone beam and fan beam collimators that are applicable over the full collimator's field-of-view. The sensitivity was calculated by integrating previously derived collimator response functions over the full detector surface. Contrary to common approximations, the varying solid angle subtended by different detector pixels was fully taken into account which results in a closed-form non-diverging formula for the sensitivity. We validated these expressions using ray-tracing simulations of a fan beam and an astigmatic cone beam collimator and found close agreement between the simulations and the sensitivity expression. The largest differences with the simulation were found close to the collimator, where sensitivity depends on the exact placement of holes and septa, while our expression represents an average over all possible placements as is common practice for analytical sensitivity expressions. We checked that average differences between the analytical expression and simulations reduced to less than 1% of the maximum sensitivity when we averaged our simulations over different septa locations. Moreover, we found that our new expression reduced to the traditional diverging formula under certain assumptions. Therefore, the newly derived sensitivity expression may enable the optimization of converging collimators for a wide range of applications, in particular when the focus is close to, or in, the object of interest.@en

The recently developed versatile emission computed tomography (VECTor) technology enables high-energy SPECT and simultaneous SPECT and PET of small animals at sub-mm resolutions. VECTor uses dedicated clustered pinhole collimators mounted in a scanner with three stationary large-area NaI(Tl) gamma detectors. Here, we develop and validate dedicated image reconstruction methods that compensate for image degradation by incorporating accurate models for the transport of high-energy annihilation gamma photons. Ray tracing software was used to calculate photon transport through the collimator structures and into the gamma detector. Input to this code are several geometric parameters estimated from system calibration with a scanning 99mTc point source. Effects on reconstructed images of (i) modelling variable depth-of-interaction (DOI) in the detector, (ii) incorporating photon paths that go through multiple pinholes ('multiple-pinhole paths' (MPP)), and (iii) including various amounts of point spread function (PSF) tail were evaluated. Imaging 18F in resolution and uniformity phantoms showed that including large parts of PSFs is essential to obtain good contrast-noise characteristics and that DOI modelling is highly effective in removing deformations of small structures, together leading to 0.75 mm resolution PET images of a hot-rod Derenzo phantom. Moreover, MPP modelling reduced the level of background noise. These improvements were also clearly visible in mouse images. Performance of VECTor can thus be significantly improved by accurately modelling annihilation gamma photon transport.@en

In recent years, breast imaging using radiolabelled molecules has attracted significant interest. Our group has proposed a multi-pinhole molecular breast tomosynthesis (MP-MBT) scanner to obtain 3D functional molecular breast images at high resolutions. After conducting extensive optimisation studies using simulations, we here present a first prototype of MP-MBT and evaluate its performance using physical phantoms. The MP-MBT design is based on two opposing gamma cameras that can image a lightly compressed pendant breast. Each gamma camera consists of a 250 × 150 mm2 detector equipped with a collimator with multiple pinholes focusing on a line. The NaI(Tl) gamma detector is a customised design with 3.5 mm intrinsic spatial resolution and high spatial linearity near the edges due to a novel light-guide geometry and the use of square PMTs. A volume-of-interest is scanned by translating the collimator and gamma detector together in a sequence that optimises count yield from the scan region. Derenzo phantom images showed that the system can reach 3.5 mm resolution for a clinically realistic 99mTc activity concentration in an 11-minute scan, while in breast phantoms the smallest spheres visible were 6 mm in diameter for the same scan time. To conclude, the experimental results of the novel MP-MBT scanner showed that the setup had sub-centimetre breast tumour detection capability which might facilitate 3D molecular breast cancer imaging in the future.@en

Planar molecular breast imaging (MBI) is rapidly gaining in popularity in diagnostic oncology. To add 3D capabilities, we introduce a novel molecular breast tomosynthesis (MBT) scanner concept based on multi-pinhole collimation. In our design, the patient lies prone with the pendant breast lightly compressed between transparent plates. Integrated webcams view the breast through these plates and allow the operator to designate the scan volume (e.g. a whole breast or a suspected region). The breast is then scanned by translating focusing multi-pinhole plates and NaI(Tl) gamma detectors together in a sequence that optimizes count yield from the volume-of-interest. With simulations, we compared MBT with existing planar MBI. In a breast phantom containing different lesions, MBT improved tumour-to-background contrast-to-noise ratio (CNR) over planar MBI by 12% and 111% for 4.0 and 6.0 mm lesions respectively in case of whole breast scanning. For the same lesions, much larger CNR improvements of 92% and 241% over planar MBI were found in a scan that focused on a breast region containing several lesions. MBT resolved 3.0 mm rods in a Derenzo resolution phantom in the transverse plane compared to 2.5 mm rods distinguished by planar MBI. While planar MBI cannot provide depth information, MBT offered 4.0 mm depth resolution. Our simulations indicate that besides offering 3D localization of increased tracer uptake, multi-pinhole MBT can significantly increase tumour-to-background CNR compared to planar MBI. These properties could be promising for better estimating the position, extend and shape of lesions and distinguishing between single and multiple lesions.@en

Accurate gamma photon transport simulations of emission tomography systems are important to optimise system geometries and for iterative image reconstruction. Monte Carlo simulation (MCS) is widely established for this purpose but has the disadvantage of being prohibitively slow. Voxelized ray tracing (VRT) can be used as an alternative but the accuracy of VRT needs to be assessed for each simulation task at hand. The aim of this work is to propose and validate a dedicated VRT code for a novel radionuclide-based multi-pinhole molecular breast tomosynthesis (MP-MBT) scanner. The MP-MBT system images radionuclide distributions in a mildly compressed breast using two opposing gamma cameras, each equipped with a focusing multi-pinhole collimator, that slide along opposite sides of the breast. VRT simulates gamma photon transport by tracing rays efficiently through the voxelized phantom, collimator, and detector volumes using Siddon's raytracing algorithm, accelerated by dual-grid methods. To assess its accuracy, we compare point spread functions (PSFs) calculated with VRT for different voxel sizes with those generated by the established MCS toolkit GATE. Furthermore, VRT and MCS-simulated projections of realistic anthropomorphic XCAT phantoms with different compressed breast sizes are compared, as well as reconstructed images obtained from these projections. With VRT, PSFs for MP-MBT can be simulated accurately when the fine voxel size of the VRT's dual-grid is 1/8 mm. Reaching a similar deviation from noiseless PSFs takes 29 300 times longer with full MCS than with VRT. Furthermore, XCAT phantom simulations show that VRT-generated projections are very close to MCS-generated low-noise projections when these are corrected for scatter by the triple energy window method. However, we also find that primary gamma photons from the torso may in some cases reach the detector, meaning that torso activity should not be neglected in VRT. Finally, reconstructed images obtained from projections generated by VRT and MCS are visually very similar and have no significant difference in contrast and noise characteristics. We conclude that VRT can accurately and efficiently simulate MP-MBT even though it neglects scattered photons originating from the torso.@en

Accurate gamma photon transport simulations of emission tomography systems are important to optimise system geometries and for iterative image reconstruction. Monte Carlo simulation (MCS) is widely established for this purpose but has the disadvantage of being prohibitively slow. Voxelized ray tracing (VRT) can be used as an alternative but the accuracy of VRT needs to be assessed for each simulation task at hand. The aim of this work is to propose and validate a dedicated VRT code for a novel radionuclide-based multi-pinhole molecular breast tomosynthesis (MP-MBT) scanner. The MP-MBT system images radionuclide distributions in a mildly compressed breast using two opposing gamma cameras, each equipped with a focusing multi-pinhole collimator, that slide along opposite sides of the breast. VRT simulates gamma photon transport by tracing rays efficiently through the voxelized phantom, collimator, and detector volumes using Siddon's raytracing algorithm, accelerated by dual-grid methods. To assess its accuracy, we compare point spread functions (PSFs) calculated with VRT for different voxel sizes with those generated by the established MCS toolkit GATE. Furthermore, VRT and MCS-simulated projections of realistic anthropomorphic XCAT phantoms with different compressed breast sizes are compared, as well as reconstructed images obtained from these projections. With VRT, PSFs for MP-MBT can be simulated accurately when the fine voxel size of the VRT's dual-grid is 1/8 mm. Reaching a similar deviation from noiseless PSFs takes 29 300 times longer with full MCS than with VRT. Furthermore, XCAT phantom simulations show that VRT-generated projections are very close to MCS-generated low-noise projections when these are corrected for scatter by the triple energy window method. However, we also find that primary gamma photons from the torso may in some cases reach the detector, meaning that torso activity should not be neglected in VRT. Finally, reconstructed images obtained from projections generated by VRT and MCS are visually very similar and have no significant difference in contrast and noise characteristics. We conclude that VRT can accurately and efficiently simulate MP-MBT even though it neglects scattered photons originating from the torso.@en

Accurate gamma photon transport simulations of emission tomography systems are important to optimise system geometries and for iterative image reconstruction. Monte Carlo simulation (MCS) is widely established for this purpose but has the disadvantage of being prohibitively slow. Voxelized ray tracing (VRT) can be used as an alternative but the accuracy of VRT needs to be assessed for each simulation task at hand. The aim of this work is to propose and validate a dedicated VRT code for a novel radionuclide-based multi-pinhole molecular breast tomosynthesis (MP-MBT) scanner. The MP-MBT system images radionuclide distributions in a mildly compressed breast using two opposing gamma cameras, each equipped with a focusing multi-pinhole collimator, that slide along opposite sides of the breast. VRT simulates gamma photon transport by tracing rays efficiently through the voxelized phantom, collimator, and detector volumes using Siddon's raytracing algorithm, accelerated by dual-grid methods. To assess its accuracy, we compare point spread functions (PSFs) calculated with VRT for different voxel sizes with those generated by the established MCS toolkit GATE. Furthermore, VRT and MCS-simulated projections of realistic anthropomorphic XCAT phantoms with different compressed breast sizes are compared, as well as reconstructed images obtained from these projections. With VRT, PSFs for MP-MBT can be simulated accurately when the fine voxel size of the VRT's dual-grid is 1/8 mm. Reaching a similar deviation from noiseless PSFs takes 29 300 times longer with full MCS than with VRT. Furthermore, XCAT phantom simulations show that VRT-generated projections are very close to MCS-generated low-noise projections when these are corrected for scatter by the triple energy window method. However, we also find that primary gamma photons from the torso may in some cases reach the detector, meaning that torso activity should not be neglected in VRT. Finally, reconstructed images obtained from projections generated by VRT and MCS are visually very similar and have no significant difference in contrast and noise characteristics. We conclude that VRT can accurately and efficiently simulate MP-MBT even though it neglects scattered photons originating from the torso.@en

Recently, we proposed and optimized dedicated multi-pinhole molecular breast tomosynthesis (MBT) that images a lightly compressed breast. As MBT may also be performed with other types of collimators, the aim of this paper is to optimize MBT with fan beam and slit-slat collimators and to compare its performance to that of multi-pinhole MBT to arrive at a truly optimized design. Using analytical expressions, we first optimized fan beam and slit-slat collimator parameters to reach maximum sensitivity at a series of given system resolutions. Additionally, we performed full system simulations of a breast phantom containing several tumours for the optimized designs. We found that at equal system resolution the maximum achievable sensitivity increases from pinhole to slit-slat to fan beam collimation with fan beam and slit-slat MBT having on average a 48% and 20% higher sensitivity than multi-pinhole MBT. Furthermore, by inspecting simulated images and applying a tumour-to-background contrast-to-noise (TB-CNR) analysis, we found that slit-slat collimators underperform with respect to the other collimator types. The fan beam collimators obtained a similar TB-CNR as the pinhole collimators, but the optimum was reached at different system resolutions. For fan beam collimators, a 6–8 mm system resolution was optimal in terms of TB-CNR, while with pinhole collimation highest TB-CNR was reached in the 7–10 mm range.@en

In this thesis, we developed tools to model and simulate single gamma imaging modalities. These tools were subsequently used to simulate and optimize a design for Molecular Breast Tomosynthesis (MBT). The chapters in this thesis can be divided into two groups. The first group consists of Chapters 3 and 4 detailing technical topics concerning modelling and simulation of single gamma emitter imaging modalities, while the other chapters deal with MBT and its optimization. We will first discuss the chapters from the first group individually and then discuss the latter group. @en

Imaging of 99mTc-labelled tracers is gaining popularity for detecting breast tumours. Recently, we proposed a novel design for molecular breast tomosynthesis (MBT) based on two sliding focusing multi-pinhole collimators that scan a modestly compressed breast. Simulation studies indicate that MBT has the potential to improve the tumour-to-background contrast-to-noise ratio significantly over state-of-the-art planar molecular breast imaging. The aim of the present paper is to optimize the collimator-detector geometry of MBT. Using analytical models, we first optimized sensitivity at different fixed system resolutions (ranging from 5 to 12 mm) by tuning the pinhole diameters and the distance between breast and detector for a whole series of automatically generated multi-pinhole designs. We evaluated both MBT with a conventional continuous crystal detector with 3.2 mm intrinsic resolution and with a pixelated detector with 1.6 mm pixels. Subsequently, full system simulations of a breast phantom containing several lesions were performed for the optimized geometry at each system resolution for both types of detector. From these simulations, we found that tumour-to-background contrast-to-noise ratio was highest for systems in the 7 mm–10 mm system resolution range over which it hardly varied. No significant differences between the two detector types were found.@en