Monte Carlo PET simulations: effect of photon polarization on scatter estimates

High accuracy multiple scatter modelling for 3D whole body PET

A new technique for modelling multiple-order Compton scatter which uses the absolute probabilities relating the image space to the projection space in 3D whole body PET is presented. The details considered in this work give a valuable insight into the scatter problem, particularly for multiple scatter. Such modelling is advantageous for large attenuating media where scatter is a dominant component of the measured data, and where multiple scatter may dominate the total scatter depending on the energy threshold and object size. The model offers distinct features setting it apart from previous research: (1) specification of the scatter distribution for each voxel based on the transmission data, the physics of Compton scattering and the specification of a given PET system; (2) independence from the true activity distribution; (3) in principle no scaling or iterative process is required to find the distribution; (4) explicit multiple scatter modelling; (5) no scatter subtraction or addition to the forward model when included in the system matrix used with statistical image reconstruction methods; (6) adaptability to many different scatter compensation methods from simple and fast to more sophisticated and therefore slower methods; (7) accuracy equivalent to that of a Monte Carlo model. The scatter model has been validated using Monte Carlo simulation (SimSET).

Model-based scatter correction for fully 3D PET

A method is presented that directly calculates the mean number of scattered coincidences in data acquired with fully 3D positron emission tomography (PET). This method uses a transmission scan, an emission scan, the physics of Compton scatter, and a mathematical model of the scanner in a forward calculation of the number of events for which one photon has undergone a single Compton interaction. The distribution of events for which multiple Compton interactions have occurred is modelled as a linear transformation of the single-scatter distribution. Computational efficiency is achieved by sampling at rates no higher than those required by the scatter distribution and by implementing the algorithm using look-up tables. Evaluation studies in phantoms with large scatter fractions show that the method yields images with quantitative accuracy equivalent to that of slice-collimated PET in clinically useful times.

Optimization of a polarized source for in vivo x-ray fluorescence analysis of platinum and other heavy metals

The Monte Carlo method was used to optimize a polarized photon source for the x-ray fluorescence analysis of platinum and other heavy metals in vivo. The source consisted of a 140 kVp, 25 mA x-ray tube with the photons plane-polarized by 90 degrees scattering. The use of plane-polarized photons results in a significant reduction in background when the fluorescent radiation is measured along the direction of polarization. A Monte Carlo computer programme was written to simulate the production and interaction of polarized photons in order to determine the optimal polarizing material and dimensions, together with beam width and geometrical arrangement of source, polarizer and beam collimators. Calculated photon energy distributions are compared with experimental data to test the validity of the model. The best configuration of the polarization system for the in vivo analysis of platinum consisted of a 20 mm Cu polarizing block with a secondary collimator subtending a 0.1 radian angle at the polarizer.