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.     

Experimental study on the location of energy windows for scatter correction by the TEW method in 201Tl imaging]

To investigate validity of scatter correction by the TEW method in 201Tl imaging, we performed an experimental study using the gamma camera with the capability to perform the TEW method and a plate source with a defect. Images were acquired with the triple energy window which is recommended by the gamma camera manufacturer. The result of the energy spectrum showed that backscattered photons were included within the lower sub-energy window and main energy window, and the spectral shapes in the upper half region of the photopeak (70 keV) were not changed greatly by the source shape and the thickness of scattering materials. The scatter fraction calculated using energy spectra and, visual observation and the contrast values measured at the defect using planar images also showed that substantial primary photons were included in the upper sub-energy window. In TEW method (for scatter correction), two sub-energy windows are expected to be defined on the part of energy region in which total counts mainly consist of scattered photons. Therefore, it is necessary to investigate the use of the upper sub-energy window on scatter correction by the TEW method in 201Tl imaging.     

Transmission-based scatter correction of 180 degrees myocardial single-photon emission tomographic studies

Meaningful comparison of single-photon emission tomographic (SPET) reconstructions for data acquired over 180 degrees or 360 degrees can only be performed if both attenuation and scatter correction are applied. Convolution subtraction has appeal as a practical method for scatter correction; however, it is limited to data acquired over 360 degrees. A new algorithm is proposed which can be applied equally well to data acquired over 180 degrees or 360 degrees. The method involves estimating scatter based on knowledge of reconstructed transmission data in combination with a reconstructed estimate of the activity distribution, obtained using attenuation correction with broad beam attenuation coefficients. Processing is implemented for planes of activity parallel to the projection images for which a simplified model for the scatter distribution may be applied, based on the measured attenuation. The appropriate broad beam (effective) attenuation coefficients were determined by considering the scatter buildup equation. It was demonstrated that narrow beam attenuation coefficients should be scaled by 0.75 and 0.65 to provide broad beam attenuation coefficients for technetium-99m and thallium-201 respectively. Using a thorax phantom, quantitative accuracy of the new algorithm was compared with conventional transmission-based convolution subtraction (TDCS) for 360 degrees data. Similar heart to lung contrasts were achieved and correction of 180 degrees data yielded a 10.4% error for cardiac activity compared to 5.2% for TDCS. Contrast for myocardium to ventricular cavity was similarly good for scatter-corrected 180 degrees and 360 degrees data, in contrast to attenuation-corrected data, where contrast was significantly reduced. The new algorithm provides a practical method for correction of scatter applicable to 180 degrees myocardial SPET.     

Combined scatter and attenuation correction for 201Tl myocardial perfusion SPECT using OS-EM algorithm

There are two possible ways to obtain scatter-corrected images with the ML-EM (maximum likelihood expectation maximization) algorithm: one is the subtraction of scatter estimate si from projection data pi, and then (pi-si) is used for scatter-corrected projection data (denoted as SC(T)); the other method is the addition of scatter estimate si to the projections calculated from the reconstructed image without performing data subtraction (SC(E)). This paper investigated these two ML-EM algorithms of combined scatter and attenuation correction on 201Tl myocardial perfusion SPECT imaging. Scatter windows were placed one full width at half maximum (FWHM) below and above the photopeak centerline. The scatter fraction in the primary peak was estimated using trapezoidal approximation by the triple energy window method. Phantom and clinical images were reconstructed using 6 iterations of ordered subsets EM algorithm (OS-EM). A cylindrical phantom with a cold-rod insert and a heart/thorax phantom with liver insert were used to evaluate scatter and the attenuation compensation technique. A cylindrical phantom filled with uniform 201Tl solution was used to evaluate statistical noise. The percent root-mean-square uncertainty (%RMSU) was used as a quantitative measure of noise amplification. %RMSU showed that the SC(E) method amplified noise less in comparison with the SC(T) method, however, no significant difference in image quality was observed between the two methods. In conclusion, both the SC(T) and SC(E) methods provided significant and similar improvement in the removal of scatter in 201Tl myocardial perfusion SPECT imaging.     

A CT assisted method for absolute quantitation of internal radioactivity

A method is described for the determination of radioactivity (microCi or MBq) at an organ site within an object or patient. Using both anatomic image data (CT or MRI scans) and planar gamma camera images, activity at depth is determined using a matrix inversion method based on least squares. The result of the inversion analysis was the unknown set of n linear (uniform) activity densities representative of each organ within the phantom or patient. The problem was overdetermined since the number of unknown activity densities (microCi/cm) was much less than the number of analysis points (N) within the nuclear image. This method, defined as the CT assisted matrix inversion (CAMI) technique, was accurate to within 15% for a three "organ" plastic phantom, wherein the organs were right circular cylinders having activities of 74 to 508 microCi (or 2.74 MBq to 18.8 MBq). This accuracy included image quantitation effects, particularly assumptions concerning attenuation correction. The average absolute percent error of the estimated activity in four distinct radioactive volumes in the phantom was 9.8%. It was found that the background activity within the phantom was estimated to be too high if sampling regions near strong sources were used in the analysis (scatter effect). This was minimized by going at least 2 cm away from such sources. By applying the method to a monoclonal antibody clinical study, activities within the patient's major organs such as liver, spleen, and kidney could be estimated, even in cases where the organ could not be visualized. Here, the CAMI algorithm gave internally consistent results for the patient's left and right lung linear activity concentrations. The CAMI technique resolves the problem of tissue superimposition using depth information from 3-D CT and is applicable in cases where a number of organs overlap in the gamma camera image. Thus, the method should be generally useful to nuclear image quantitation and the estimation of absorbed radiation doses in patients. One particular application is the estimation of radiation doses in radioimmunotherapy (RIT).     

Slot technique--an alternative method of scatter reduction in radiography

The most common method of scatter reduction is the use of an antiscatter grid. Its disadvantage is the absorption of a certain percentage of primary radiation in the lead strips of the grid and the fact that due to the limited thickness of the lead strips their scatter absorption is also limited. A possibility for avoiding this disadvantage is offered by the so-called slot technique, ie, the successive exposure of the subject with a narrow fan beam provided by slots in rather thick lead plates. In this paper the results of a comparison between grid and slot technique regarding dose to the patient, scatter reduction, image quality and the effect of automatic exposure control are reported.     

Two-dimensional scatter integration method for brachytherapy dose calculations in 3D geometry

In brachytherapy clinical practice, applicator shielding and tissue heterogeneities are usually not explicitly taken into account. None of the existing dose computational methods are able to reconcile accurate dose calculation in complex three-dimensional (3D) geometries with high efficiency and simplicity. We propose a new model that performs two-dimensional integration of the scattered dose component. The model calculates the effective primary dose at the point of interest and estimates the scatter dose as a superposition of the scatter contributions from pyramid-shaped minibeams. The approach generalizes a previous scatter subtraction model designed to calculate the dose for axial points in simple cylindrically symmetric geometry by dividing the scattering volume into spatial regions coaxial with the source-to-measurement point direction. To allow for azimuthal variation of the primary dose, these minibeams were divided into equally spaced azimuthally distributed pyramidal volumes. The model uses precalculated scatter-to-primary ratios (SPRs) for collimated isotropic sources. Effective primary dose, which includes the radiation scattered in the source capsule, is used to achieve independence from the source structure. For realistic models of the 192Ir HDR and PDR sources, the algorithm agrees with Monte Carlo within 2.5% and for the 125I type 6702 seed within 6%. The 2D scatter integration (2DSI) model has the potential to estimate the dose behind high-density heterogeneities both accurately and efficiently. The algorithm is much faster than Monte Carlo methods and predicts the dose around sources with different gamma-ray energies and differently shaped capsules with high accuracy.     

Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction

Accurate scatter compensation in SPECT can be performed by modelling the scatter response function during the reconstruction process. This method is called reconstruction-based scatter compensation (RBSC). It has been shown that RBSC has a number of advantages over other methods of compensating for scatter, but using RBSC for fully 3D compensation has resulted in prohibitively long reconstruction times. In this work we propose two new methods that can be used in conjunction with existing methods to achieve marked reductions in RBSC reconstruction times. The first method, coarse-grid scatter modelling, significantly accelerates the scatter model by exploiting the fact that scatter is dominated by low-frequency information. The second method, intermittent RBSC, further accelerates the reconstruction process by limiting the number of iterations during which scatter is modelled. The fast implementations were evaluated using a Monte Carlo simulated experiment of the 3D MCAT phantom with 99mTc tracer, and also using experimentally acquired data with 201Tl tracer. Results indicated that these fast methods can reconstruct, with fully 3D compensation, images very similar to those obtained using standard RBSC methods, and in reconstruction times that are an order of magnitude shorter. Using these methods, fully 3D iterative reconstruction with RBSC can be performed well within the realm of clinically realistic times (under 10 minutes for 64 x 64 x 24 image reconstruction).     

Timedelay based all-frequency correction method for discrete spectrum

Many spectrum correction methods have been developed, but their performance degrades signif- icantly when they are applied to the correction of low frequency component ( LFC ). It owns to that the model underlying the conventional approaches neglects the interference of the negative frequency in the real signal. A new approach for the correction of the LFC is proposed, which suits all kinds of symmetrical windows. It divides a signal into three sections and makes use of the first spectrum line of each section. Then this approach is modified so that it is also applicable to the correction of the high frequency component. Thus a timedelay-based all-frequency correction method is proposed. The simulation results show that this method is simple and feasible. By this method, the accurate frequency, amplitude and phase of the spectral line can be obtained whether it is close to or far from OHz.     

Role of the radionuclide brain image in the diagnosis of brainstem gliomas

The noninvasive detection of brainstem gliomas remains difficult. Eleven of our patients with proven brainstem gliomas had radionuclide brain imaging prior to the initiation of therapy or confirmation of the diagnosis; six studies were positive. Pneumoencephalography remains the most reliable diagnostic test for brainstem glioma, and is invariably required for confirmation. Although angiography is useful in the evaluation of vasocularity it may not detect small infiltrating lesions. Radionuclide brain imaging is useful in the initial workup of patients with suspected brainstem gliomas.     

Digital image ratio: a new radiographic method for quantifying changes in alveolar bone. Part II: Clinical application

BACKGROUND: Men with benign prostatic hyperplasia can be treated with alpha 1-adrenergic-antagonist drugs that relax prostatic smooth muscle or with drugs that inhibit 5 alpha-reductase and therefore reduce tissue androgen concentrations. However, the effects of the two types of drugs have not been compared. METHODS: We compared the safety and efficacy of placebo, terazosin (10 mg daily), finasteride (5 mg daily), and the combination of both drugs in 1229 men with benign prostatic hyperplasia. American Urological Association symptom scores and peak urinary-flow rates were determined at base line and periodically for one year. RESULTS: The mean changes from base line in the symptom scores in the placebo, finasteride, terazosin, and combination-therapy groups at one year were decreases of 2.6, 3.2, 6.1, and 6.2 points, respectively (P<0.001 for the comparisons of both terazosin and combination therapy with finasteride and with placebo). The mean changes at one year in the peak urinary-flow rates were increases of 1.4, 1.6, 2.7, and 3.2 ml per second, respectively (P<0.001 for the comparisons of both terazosin and combination therapy with finasteride and with placebo). Finasteride had no more effect on either measure than placebo. In the placebo group, 1.6 percent of the men discontinued the study because of adverse effects, as did 4.8 to 7.8 percent of the men in the other three groups. CONCLUSIONS: In men with benign prostatic hyperplasia, terazosin was effective therapy, whereas finasteride was not, and the combination of terazosin and finasteride was no more effective than terazosin alone.     

计算动力触探修正系数的多项式拟合法

给出了用于计算动力触探修正系数的多项式拟合法.首先通过比对确定了拟合函数,指出用二次多项式拟合是考虑精确性和实用性的最佳契合点,通过数学分析介绍了双因素函数拟合的嵌套方法,并给出了重型动力触探修正系数在考虑杆长和实测锤击数条件下的二元二次多项式拟合函数及其行列式形式.通过拟合出的修正系数与规范所给的修正系数之间的误差分析判定拟合函数精度较好,适于工程应用.并且此公式适用于计算机计算.     

Effects of scatter and attenuation correction on quantitative assessment of regional cerebral blood flow with SPECT

The pyruvate dehydrogenase complex (PDC) plays a key role in the anaerobic mitochondrial metabolism of the parasitic nematode Ascaris suum. A cDNA coding for an A. suum pyruvate dehydrogenase kinase (APDK) has been cloned and sequenced from poly(A)+ RNA isolated from adult A. suum muscle.2 APDK exhibited significant sequence identity to mammalian PDKs. Nucleotide sequence analysis of the APDK cDNA revealed a 22-nucleotide spliced leader, characteristic of many nematode mRNAs, a 5'-UTR of 6 nucleotides, an open reading frame of 1197 nucleotides, and a 3'-UTR of 101 nucleotides that included a putative polyadenylation signal. The open reading frame predicted a protein of 399 amino acids with a molecular weight of 45,402 that included a putative 18-aminoacid leader peptide. Recombinant APDK (rAPDK) was functionally expressed in Escherichia coli with a his tag at its N-terminus and purified to apparent homogeneity on Ni-NTA-agarose. Recombinant APDK was a dimer and was not autophosphorylated and its activity was stimulated in the presence of APDK-deficient adult A. suum muscle PDC presumably by the binding of APDK to the dihydrolipoyl transacetylase (E2) core of the complex. After binding to the core, rAPDK activity was stimulated by elevated NADH/NAD+ and acetyl CoA/CoA ratios within the same ranges as observed for the native APDK. Immunoblotting suggested that native APDK focused as a series of 43-kDa spots (pI 6.1-6.8) on two-dimensional gels of the purified adult A. suum muscle PDC.     

Simple correction method for k-space trajectory deviations in MRI

STUDY OBJECTIVES: Recently, pulmonologists have performed thoracoscopy under local anesthesia using rigid thoracoscopes or flexible bronchoscopes. The latter allow greater access within the pleural cavity but are difficult to manipulate. The Olympus LTF semiflexible fiberoptic thoracoscope combines features of both instruments, having a solid body and a flexible terminal section. In the first study with this instrument, we evaluated ease of use and compared diagnostic yield with closed needle biopsy. PATIENTS: Twenty-four patients with pleural effusion were investigated. SETTING: Scottish University Hospital. DESIGN: Thoracoscopy was performed in the bronchoscopy suite after premedication with atropine and papaveretum. Following a standard Abram's needle biopsy, the LTF thoracoscope was inserted through a flexible introducer (Olympus Optical Co Ltd; Tokyo, Japan). The pleura was inspected and biopsy specimens were taken of suspicious areas. RESULTS: The final diagnosis was malignant pleural effusion in 16 of 24 patients. Ten of 16 were positive by Abram's biopsy, giving a sensitivity of 62%. Thirteen of 16 were positive by fiberoptic thoracoscopy, giving an improved sensitivity of 81%. The LTF thoracoscope was easy to use for pulmonologists experienced in rigid thoracoscopy and flexible bronchoscopy. Excellent views of the pleura were obtained from a single entry point. The procedure was well tolerated and no complications were encountered. CONCLUSION: The LTF thoracoscope allows excellent pleural access but a larger biopsy channel (currently 2 mm) might increase the accuracy of diagnosis.     

The derivation of tissue-maximum ratio from percentage depth dose requires peak scatter factor to be considered a function of source-to-surface distance

A formula for the calculation of tissue-maximum ratio (TMR) from percentage depth dose (PDD) and peak scatter factor (PSF) is derived from first principles using a simple geometric model for the case when the field size for PDD and PSF is defined at the surface. The derivation is carried out in two ways: (a) taking field size for PDD and PSF as defined at the depth of maximum dose and then applying a conversion factor, and (b) by a direct derivation. The first of these methods yields a formula which agrees with BJR Supplement 25, but the latter yields a result which differs from it. Numerically, this difference is insignificant, but it has implications for the theoretical basis of the conversion formulae. The difference arises due to the translation of field size from one depth to another when calculating PSF: two different values of source-to-surface distance (SSD) yield two apparently different PSFs for the same size of field at the depth of maximum dose. Disagreements of this type are prevalent throughout the standard conversion formulae given in BJR Supplement 25 when field size for PDD and PSF is defined at the surface rather than at the depth of dose maximum. These disagreements are illustrated here using the conversion of PDD from one SSD to another as an example. The difficulty is overcome by considering PSF to be a function of SSD as well as field size.     

Electron-beam CT: the effect of using a correction function on coronary artery calcium quantitation

The relationships between 4 bacterial and 3 viral antibody titers and morbidity (undifferentiated fever (UF)) and mortality were investigated in recently weaned beef calves. Blood samples from 100 animals that required treatment for UF (Cases) and 100 healthy control animals (Controls) were obtained: upon arrival at the feedlot (Arrival), at the time of selection as a Case or Control (Selection), and at approximately 33 d of the feeding period (Convalescent). Seroconversion to Pasteurella haemolytica antileukotoxin was associated with an increased risk of UF (OR = 2.83); however, seroconversion to bovine herpesvirus-1 G-IV glycoprotein was associated with a decreased risk of UF (OR = 0.43). Higher Arrival bovine viral diarrhea virus antibody titer was associated with a decreased risk of UF (OR = 0.83). Increases in Mycoplasma alkalescens antibody titer after Arrival were associated with an increased risk of UF (OR = 1.10). Higher Arrival Haemophilus somnus antibody titer and increases in Haemophilus somnus antibody titer after Arrival were both associated with a decreased risk of UF (OR = 0.76 and OR = 0.78). The odds of overall mortality (OR = 5.09) and hemophilosis mortality (OR = 11.31) in Cases were significantly (P < 0.05) higher than in the Controls. Higher Arrival bovine herpesvirus-1 antibody titer was associated with an increased risk of mortality (OR = 1.30). Protective immunity to Pasteurella haemolytica antileukotoxin, Haemophilus somnus, bovine herpesvirus-1 G-IV glycoprotein, bovine viral diarrhea virus, and Mycoplasma spp. may be necessary to reduce the occurrence of UF. Animals with UF are at an increased risk of overall and hemophilosis mortality.