Geometrical response of multihole collimators

A complete theory of camera multihole collimators is presented. The geometrical system response is determined in closed form in frequency space. This closed form accounts for the known efficiency and resolution formulae for parallel beam, fan beam, cone beam and astigmatic collimators as well as for the most frequent hole array patterns and hole shapes. The point spread function in the space domain for a certain collimator and source position can be calculated via a discrete fast Fourier transform. Beside the complete theoretical definition of the response of multihole collimators, this theory allows the definition of accurate models of the geometrical response for SPECT reconstruction and it is suitable for designing new collimators.     

Quinine improves the results of intensive chemotherapy in myelodysplastic syndromes expressing P glycoprotein: results of a randomized study

A method is developed which makes it possible to scan and reconstruct an object with cone beam x-rays in a spiral scan path with area detectors much shorter than the length of the object. The method is mathematically exact. If only a region of interest of the object is to be imaged, a top circle scan at the top level of the region of interest and a bottom circle scan at the bottom level of the region of interest are added. The height of the detector is required to cover only the distance between adjacent turns in the spiral projected at the detector. To reconstruct the object, the Radon transform for each plane intersecting the object is computed from the totality of the cone beam data. This is achieved by suitably combining the cone beam data taken at different source positions on the scan path; the angular range of the cone beam data required at each source position can be determined easily with a mask which is the spiral scan path projected on the detector from the current source position. The spiral scan algorithm has been successfully validated with simulated cone beam data.     

Ionization chamber shift correction and surface dose measurements in electron beams

Cylindrical ionization chambers produce perturbations (gradient and fluence) in the medium, and hence the point of measurement is not accurately defined in electron beam dosimetry. The gradient perturbation is often corrected by a shift method depending on the type of ion chamber. The shift is in the range of 0.33-0.85 times the inner radius (r) of the ion chamber, upstream from the centre of the chamber, depending upon the dosimetry protocol. This variation in shift causes the surface dose to be uncertain due to the high dose gradient. An investigation was conducted to estimate the effective point of measurement of cylindrical ion chambers in electron beams. Ionization measurements were taken with the ion chamber in air and in a phantom at source to chamber distances of <100 cm and >100 cm respectively. The data in air and in the phantom were fitted with the inverse square and electron depth dose functions, respectively. The intersection of the two functions provides an accurate estimate of the ion chamber shift and the surface dose. Our results show that the shift correction for an ion chamber is energy dependent. The measured shifts vary from 0.9r to 0.5r between 6 MeV and 20 MeV beams respectively. The surface dose measured with the ion chambers and mathematically determined values are in agreement to within 3%. The method presented in this report is unambiguous, fast and reliable for the estimation of surface dose and the shift needed in electron beam dosimetry.     

An orthotopic model of human pancreatic cancer in severe combined immunodeficient mice: potential application for preclinical studies

The Compton camera can collect SPECT data with high efficiency due to electronic collimation. The data acquired from a Compton camera are projections of source activity along cones and are approximated in this paper by cone-surface integrals. This paper proposes the use of an orthogonal spherical expansion to convert the cone-surface integrals into plane integrals. The conversion technique is efficient. Once the plane integrals are obtained, a 3D image can be reconstructed by the 3D Radon inversion formula. The algorithm is implemented and computer simulations are used to demonstrate the efficiency and accuracy of the proposed reconstruction algorithm.     

Pixel driven implementation of filtered backprojection for reconstruction of fan beam SPECT data using a position dependent effective projection bin length

Filtered backprojection is commonly implemented as a pixel driven algorithm in which the density is reconstructed on an array of grid points that are usually associated with the centres of square image pixels. In fan beam geometry, this conventional pixel driven approach using two-bin linear interpolation leads to an inefficient use of the projection data due to magnification at the projection line. For typical count limited SPECT data, this results in increased reconstructed image noise. We propose an alternative type of pixel driven algorithm that makes efficient use of projection data by averaging over all bins within a position dependent "effective' projection bin interval. The effective bin interval is equivalent to the image pixel length magnified on the projection line. This heuristic method leads to more complete use of the projection data and results in reconstructed images with superior noise properties compared with the conventional method while presenting similar spatial resolution characteristics.     

Automatic three-dimensional matching of CT-SPECT and CT-CT to localize lung damage after radiotherapy

The aim of this study was to develop a fast and clinically robust automatic method to register SPECT and CT scans of the lungs. METHODS: CT and SPECT scans were acquired in the supine position from 20 patients with healthy lungs. After partial irradiation of the lungs by radiotherapy, the scans were repeated. Two matching methods were compared: a conventional method with external skin markers and a new method using chamfer matching of the lung contours. In the latter method, a unique value for the SPECT threshold, needed for segmentation of the SPECT lungs, was determined by iteratively applying the chamfer matching algorithm. RESULTS: The new technique for CT-SPECT matching could be implemented in a fully automatic manner and required less than 2 min. No large systematic shifts or rotations were present between the matches obtained with the marker method and the lung contour method for healthy or partially irradiated lungs. For healthy lungs, the number of ventilation SPECT counts outside the CT-defined lung was taken as a measure for a good match. This number of outside counts was slightly lower for the new method than for the conventional method, which indicates that the accuracy of the new method is at least comparable to the conventional method. For ventilation, a systematic difference between the results of the matching methods, a small translation in the anterior --> posterior direction, could be attributed to an inconsistency of the marker positions (2 mm). For perfusion, a somewhat larger anterior --> posterior shift was found, which was attributed to the gravity force. CT-CT correlation on the lung contours using chamfer matching was tested with the same dataset. For accurate matching, the CT slices encompassing the diaphragm had to be deleted. CONCLUSION: The new method based on lung contour matching is a fast, automatic procedure and allows accurate clinical follow-up.     

Quantitative measurement of cerebral blood flow by the microsphere model with super-early 123I-IMP brain SPECT

Cerebral blood flow (CBF) was quantitatively measured in 6 healthy young volunteers based on "super-early" acquisition of N-isopropyl-p-[123I]iodoamphetamine (IMP) brain SPECT obtained 4-6 min after IMP injection with a three-head rotating gamma camera and the microsphere (MS) model. The ratio of radioactivity (count/pixel/min) in the conventional early SPECT image (taken 25-55 min after IMP injection) to that in the "super-early" image for each brain region negatively correlated with regional CBF value obtained with the "super-early" MS method. This indicates that wash-out of IMP from the regions with higher CBF is faster than that from the regions with lower CBF and that CBF values are underestimated with the conventional MS method in regions with higher CBF. Regional CBF was quantitatively measured with the "super-early" MS method and the ARG method, a recently developed method based on two-compartment model. The mean cortical CBF was 52.5 +/- 7.0 (ml/100 g/min, mean +/- SD) with the "super-early" MS method and 47.5 +/- 3.3 with the ARG method. The CBF values obtained with the "super-early" MS method agreed with those previously reported with positron emission tomography. Since the MS method is theoretically the simplest model, the "super-early" MS method can be applied various disorders of the central nervous system where the behavior of IMP is not fully understood.     

An analytic model of pinhole aperture penetration for 3D pinhole SPECT image reconstruction

Photons penetrate the attenuating material close to the aperture of pinhole collimators in nuclear medicine, broadening the tails of point spread functions (PSFs) and degrading the resolution of planar and SPECT images. An analytic approximation has been developed that models this penetration contribution to the PSF for knife-edge point pinhole apertures. The approximation has the form exp(-gamma r), where r is the distance on the detector surface from the projection of the point source through the pinhole. The rolloff coefficient gamma is a function of the photon energy, point source location and the design parameters of the collimator. There was excellent agreement between measured values of gamma from photon transport simulations of I-131 point sources (364 keV emission only) and theoretical predictions from the analytic formula. Predicted gamma values from the analytic formula averaged 25% greater than measured values from experimental I-131 point source acquisitions. Photon transport simulations were performed that modelled the 364 keV and less abundant 637 and 723 keV emissions and scatter within the scintillation crystal. Measured gamma values from these simulations averaged 12% greater than the experimental values, indicating that about half of the error between the analytic formula and the experimental measurements was due to unmodelled 637 and 723 keV emissions. The remaining error may be due in part to scatter in the pinhole region and backscatter from gamma camera components behind the scintillation crystal. The analytic penetration model was used in designing Metz filters to compensate for penetration blur and these filters were applied to the projection data as part of 3D SPECT image reconstruction. Image resolution and contrast were improved in simulated and experimental I-131 tumour phantom studies. This analytic model of pinhole aperture penetration can be readily incorporated into iterative 3D SPECT pinhole reconstruction algorithms.     

Second-order correction to the Bigeleisen-Mayer equation due to the nuclear field shift

The nuclear field shift affects the electronic, rotational, and vibrational energies of polyatomic molecules. The theory of the shifts in molecular spectra has been studied by Schlembach and Tiemann [Schlembach, J. & Tiemann, E. (1982) Chem. Phys. 68, 21]; measurements of the electronic and rotational shifts of the diatomic halides of Pb and Tl have been made by Tiemann et al. [Tiemann, E., Kn?ckel, H. & Schlembach, J. (1982) Ber. Bunsenges. Phys. Chem. 86, 821]. These authors have estimated the relative shifts in the harmonic frequencies of these compounds due to the nuclear field shift to be of the order of 10(-6). I have used this estimate of the relative shift in vibrational frequency to calculate the correction to the harmonic oscillator approximation to the isotopic reduced partition-function ratio (208)Pb(32)S/(207)Pb(32)S. The correction is 0.3% of the harmonic oscillator value at 300 K. In the absence of compelling evidence to the contrary, it suffices to calculate the nuclear field effect on the total isotopic partition-function ratio from its shift of the electronic zero point energy and the unperturbed molecular vibration.     

Red, green, and red-green hybrid pigments in the human retina: correlations between deduced protein sequences and psychophysically measured spectral sensitivities

To analyze the human red, green, and red-green hybrid cone pigments in vivo, we studied 41 male dichromats, each of whose X chromosome carries only a single visual pigment gene (single-gene dichromats). This simplified arrangement avoids the difficulties of complex opsin gene arrays and overlapping cone spectral sensitivities present in trichromats and of multiple genes encoding identical or nearly identical cone pigments in many dichromats. It thus allows for a straightforward correlation between each observer's spectral sensitivity measured at the cornea and the amino acid sequence of his visual pigment. For each of the 41 single-gene dichromats we determined the amino acid sequences of the X-linked cone pigment as deduced from its gene sequence. To correlate these sequences with spectral sensitivities in vivo, we determined the Rayleigh matches to different red/green ratios for 29 single-gene dichromats and measured psychophysically the spectral sensitivity of the remaining green (middle wavelength) or red (long wavelength) cones in 37 single-gene dichromats. Cone spectral sensitivity maxima obtained from subjects with identical visual pigment amino acid sequences show up to a approximately 3 nm variation from subject to subject, presumably because of a combination of inexact (or no) corrections for variation in preretinal absorption, variation in photopigment optical density, optical effects within the photoreceptor, and measurement error. This variation implies that spectral sensitivities must be averaged over multiple subjects with the same genotype to obtain representative values for a given pigment. The principal results of this study are that (1) approximately 54% of the single-gene protanopes (and approximately 19% of all protanopes) possess any one of several 5'red-3'green hybrid genes that encode anomalous pigments and that would be predicted to produce protanomaly if present in anomalous trichromats; (2) the alanine/serine polymorphism at position 180 in the red pigment gene produces a spectral shift of approximately 2.7 nm; (3) for each exon the set of amino acids normally associated with the red pigment produces spectral shifts to longer wavelengths, and the set of amino acids normally associated with the green pigment produces spectral shifts to shorter wavelengths; and (4) changes in exons 2, 3, 4, and 5 from green to red are associated with average spectral shifts to long wavelengths of approximately 1 nm (range, -0.5 to 2.5 nm), approximately 3.3 nm (range, -0.5 to 7 nm), approximately 2.8 nm (range, -0.5 to 6 nm), and approximately 24.9 nm (range, 22.2-27.6 nm).     

Comparison of four motion correction techniques in SPECT imaging of the heart: a cardiac phantom study

The aim of this study was to evaluate the accuracy of four different motion correction techniques in SPECT imaging of the heart. METHODS: We evaluated three automated techniques: the cross-correlation (CC) method, diverging squares (DS) method and two-dimensional fit method and one manual shift technique (MS) using a cardiac phantom. The phantom was filled with organ concentrations of 99mTc closely matching those seen in patient studies. The phantom was placed on a small sliding platform connected to a computer-controlled stepping motor. Linear, random, sinusoidal and bounce motions of magnitude up to 2 cm in the axial direction were simulated. Both single- and dual-detector 90 degrees acquisitions were acquired using a dual 90 degrees detector system. Data were acquired over 180 degrees with 30 or 15 frames/detector (single-/dual-head) at 30 sec/frame in a 64x64 matrix. RESULTS: The simulated single-detector system, CC method, failed to accurately correct for any of the simulated motions. The DS technique overestimated the magnitude of phantom motion, particularly for images acquired between 45 degrees left anterior oblique and 45 degrees left posterior oblique. The two-dimensional and MS techniques accurately corrected for motion. The simulated dual 90 degrees detector system, CC method, only partially tracked random or bounce cardiac motion and failed to detect sinusoidal motion. The DS technique overestimated motion in the latter half of the study. Both the two-dimensional and MS techniques provided superior tracking, although no technique was able to accurately track the rapid changes in cardiac location simulated in the random motion study. Average absolute differences between true and calculated position of the heart on single- and dual 90 degrees -detectors were 1.7 mm and 1.5 mm for the two-dimensional and MS techniques, respectively. The corresponding values for the DS and CC techniques were 5.7 and 8.9 mm, respectively. CONCLUSION: Of the four techniques evaluated, manual correction by an experienced technologist proved to be the most accurate, although results were not significantly different from those observed with the two-dimensional method. Both techniques accurately determined cardiac location and permitted artifact-free reconstruction of the simulated cardiac studies.     

A simple algorithm for planar image registration in radiation therapy

A simple algorithm is presented for planar image registration and the method is applied to the simulator and portal image registration for patient setup verification in radiation therapy. Basically, the algorithm follows the concept proposed by Balter et al. [Med. Phys. 19, 329-334 (1992)], which converts the problem of open curve registration into matching a series of points along the curves. Balter's algorithm consists of three steps: (1) to determine a common starting point for each curve pair, (2) acquire two corresponding point sets along each curve, and (3) obtain a global transform matrix by matching two point sets. We integrate all three steps into one simple procedure which fits the sampled points along the intended curve pair by taking the relative path length shift as an independent fitting parameter. After being modified, the algorithm is able to take the different magnification factors of images into account, and it avoids curvature calculations. Numerical simulation as well as clinical and phantom images have been utilized to test the accuracy of the algorithm. The typical errors are less than 1 mm in translation and 1 degree in rotation. We also made a comparison study with the chamfer method. The results of the two methods agree to within 0.5 mm in translation and 0.5 degree in rotation.     

A method of correlating and merging cerebral morphology and function by a special head holder

A method is introduced which records and correlates topographically identical morphological and functional image data by means of a special head holder system which is adaptable to different modalities (MRI, PET, SPECT). It is based on a commercially available thermoplastic headmask. A thermoplastic form is used to mold an individual headmask (time required: 7 min) for all modalities providing a fixation of the patient during acquisition. This method guarantees an exact repositioning and therefore the same slice orientation of the images. The in-plane correlation is performed by adapting standard offset parameters determined with a homogeneous head phantom. By fusion (merging) of the brain outline contours or images themselves the accuracy was tested in MRI vs. PET resp. SPECT (transaxial accuracy: < or = 2.0 mm, axial: < or = 3.0 mm in patients with MRI-slice thickness of 6.0 mm).     

Interactive Failure of Two Impacting Beams

A complete analysis of an inelastic beam-to-beam impact is presented. Both beams are of solid, rectangular cross-sections. The problem is formulated based on the momentum conservation and the kinematic and dynamic continuity conditions at the moving wavefront. Closed-form solutions are obtained for transient transverse velocities, deflection profiles, and tensile strains based on the moderately large deflection theory of the beam. Three different regimes of the solution are distinguished, depending on relative values of mass ratios and wave speed ratios. Critical impact velocities to break either of the beams are determined for both the striking and struck beams by assuming both beams fail by tensile necking. However, location of the fracture point depends on relative values of various parameters. It can be right in the contact zone or away from it. It is also shown that after one beam breaks, the other beam will deform further without breaking.     

Oxygen inhalation can differentiate gaseous from nongaseous microemboli detected by transcranial Doppler ultrasound

BACKGROUND AND PURPOSE: Clinically silent circulating microemboli can be detected by transcranial Doppler sonography. The composition of these emboli in different clinical conditions is unclear. METHODS: We performed 1-hour transcranial Doppler sonographic recordings from the middle cerebral or posterior cerebral artery in 20 patients with mechanical prosthetic heart valves, in 78 patients with an arterial embolic source, and in 20 control subjects. During 30 minutes of this recording, the patients inspired room air and 6 L of oxygen per minute via a loosely fitting facial mask; during the remaining 30 minutes, they breathed room air only. RESULTS: There was a significant decline of embolic signals (ES) under oxygen in the patients with mechanical prosthetic cardiac valves (144 ES without oxygen versus 63 ES with oxygen; P = .002) but not in the patients with arterial embolic sources (145 ES without oxygen versus 135 ES with oxygen; P = NS). In the control subjects, no ES were found. CONCLUSIONS: ES in patients with mechanical prosthetic cardiac valves correspond mainly to gas bubbles. Oxygen inhibits the cavitation process of mechanical prosthetic heart valves or speeds up redissolution of gas bubbles generated by cavitation. In contrast, solid microemboli originating from thrombus or atheroma cannot be suppressed by oxygen inhalation. This simple method of oxygen inhalation should help to clarify the composition of microemboli in various clinical and experimental settings.     

Oculo-facio-cardio-dental (OFCD) syndrome

OBJECTIVE: This study examined the effects on SPECT quantitation caused by erroneous size and position of the attenuation map and inaccurate pixel size used in the Chang algorithm. METHODS: Projection data of a three-dimensional head phantom were simulated with a uniform attenuation coefficient of 0.15/cm for the inside of the phantom. Images were reconstructed using the filtered backprojection algorithm without attenuation compensation and the Chang algorithm with different attenuation maps. Quantitative comparison then was performed between the reconstructed images and the phantom. RESULTS: The pixel values obtained for noisy data by using the first-order Chang algorithm with an accurate attenuation map were less than 10% different from the true values and the left-right asymmetry was under 5%. Small errors in the geometric parameters of the attenuation map, however, caused considerable quantitative inaccuracy in the reconstructed image. For example, a 0.64-cm error in the size of the map caused 10% deviation from the true value and a 0.64-cm shift of the position of the map towards the left produced 10% left-right pixel value asymmetry. CONCLUSION: The accuracy of the Chang algorithm critically depends on the geometric parameters. For a uniform attenuator with symmetric geometry, such as the human brain, a true left-right symmetry in the pixel value can be altered significantly by a small error in the geometric parameters, while symmetry can be maintained with no attenuation compensation.     

Accurate portal dose measurement with a fluoroscopic electronic portal imaging device (EPID) for open and wedged beams and dynamic multileaf collimation

Measuring portal dose with an electronic portal imaging device (EPID) in external beam radiotherapy can be used to perform routine dosimetric quality control checks on linear accelerators and to verify treatments (in vivo dosimetry). An accurate method to measure portal dose images (PDIs) with a commercially available fluoroscopic EPID has been developed. The method accounts for (i) the optical 'cross talk' within the EPID structure, (ii) the spatially nonuniform EPID response and (iii) the nonlinearity of the EPID response. The method is based on a deconvolution algorithm. Measurement of the required input data is straightforward. The observed nonlinearity of the EPID response was largely due to the somewhat outdated EPID electronics. Nonlinearity corrections for more modern systems are expected to be smaller. The accuracy of the method was assessed by comparing PDIs measured with the EPID with PDIs measured with a scanning ionization chamber in a miniphantom, located at the same position as the fluorescent screen. For irradiations in open, wedged and intensity modulated 25 MV photon beams (produced with dynamic multileaf collimation) EPID and ionization chamber measurements agreed to within 1% (1 SD).     

Quantitative microwave imaging with a 2.45-GHz planar microwave camera

Efforts are being made to develop a new type of CT system that can scan volumes over a large range within a short time with thin slice images. One of the most promising approaches is the combination of helical scanning with multi-slice CT, which involves several detector arrays stacked in the z direction. However, the algorithm for image reconstruction remains one of the biggest problems in multi-slice CT. Two helical interpolation methods for single-slice CT, 360LI and 180LI, were used a starting points and extended to multi-slice CT. The extended methods, however, had a serious image quality problem due to the following three reasons: (1) excessively close slice positions of the complementary and direct data, resulting in a larger sampling interval; (2) the existence of several discontinuous changeovers in pairs of data samples for interpolation; and (3) the existence of cone angles. Therefore we have proposed a new algorithm to overcome the problem. It consists of the following three parts: (1) optimized sampling scan; (2) filter interpolation; and (3) fan-beam reconstruction. Optimized sampling scan refers to a special type of multi-slice helical scan developed to shift the slice position of complementary data and to acquire data with a much smaller sampling interval in the z direction. Filter interpolation refers to a filtering process performed in the z direction using several data. The normal fan-beam reconstruction technique is used. The section sensitivity profile (SSP) and image quality for four-array multi-slice CT were investigated by computer simulations. Combinations of three types of optimized sampling scan and various filter widths were used. The algorithm enables us to achieve acceptable image quality and spatial resolution at a scanning speed that is about three times faster than that for single-slice CT. The noise characteristics show that the proposed algorithm efficiently utilizes the data collected with optimized sampling scan. The new algorithm allows suitable combinations of scan and filter parameters to be selected to meet the purpose of each examination.     

Use of the selective linogram in cardiac tomography quality control

BACKGROUND: Quality control for detection of patient motion is essential in tomographic myocardial imaging. Despite significant limitations, the summation image or conventional "linogram" has long been advocated as a useful image in the detection of vertical motion. In this study a new quality control image entitled the "selective linogram" is proposed to replace the summation image in routine cardiac single-photon emission cardiac tomography (SPECT) quality control. The selective linogram is constructed in a manner somewhat analogous to the sinogram. In the sinogram, each row represents a different projection angle; in the selective linogram each column represents a different projection angle. METHODS AND RESULTS: After selection of eight motion-free studies from acquisitions at our clinical center, vertical motion of various types (bounces, shifts, and creep) were added to the projection frames. Summation image and selective linogram quality control images from these motion-containing studies and the original motion-free studies were presented in a blinded manner to two observers for scoring of patient motion. The selective linogram was significantly more accurate in allowing detection of vertical motion than was the summation image (accuracy 89% vs 47%). CONCLUSIONS: The selective linogram image is markedly superior to the summation image for the detection of vertical patient motion during cardiac SPECT. This new technique can be a valuable aid in SPECT quality control.