The PETRRA positron camera: design, characterization and results of a physical evaluation

The PETRRA positron camera is a large-area (600 mm x 400 mm sensitive area) prototype system that has been developed through a collaboration between the Rutherford Appleton Laboratory and the Institute of Cancer Research/Royal Marsden Hospital. The camera uses novel technology involving the coupling of 10 mm thick barium fluoride scintillating crystals to multi-wire proportional chambers filled with a photosensitive gas. The performance of the camera is reported here and shows that the present system has a 3D spatial resolution of approximately 7.5 mm full-width-half-maximum (FWHM), a timing resolution of approximately 3.5 ns (FWHM), a total coincidence count-rate performance of at least 80-90 kcps and a randoms-corrected sensitivity of approximately 8-10 kcps kBq(-1) ml. For an average concentration of 3 kBq ml(-1) as expected in a patient it is shown that, for the present prototype, approximately 20% of the data would be true events. The count-rate performance is presently limited by the obsolete off-camera read-out electronics and computer system and the sensitivity by the use of thin (10 mm thick) crystals. The prototype camera has limited scatter rejection and no intrinsic shielding and is, therefore, susceptible to high levels of scatter and out-of-field activity when imaging patients. All these factors are being addressed to improve the performance of the camera. The large axial field-of-view of 400 mm makes the camera ideally suited to whole-body PET imaging. We present examples of preliminary clinical images taken with the prototype camera. Overall, the results show the potential for this alternative technology justifying further development.

Statistical pixelwise inference models for planar data analysis: an application to gamma-camera uniformity monitoring

In this paper two tests based on statistical models are presented and used to assess, quantify and provide positional information of the existence of bias and/or variations between planar images acquired at different times but under similar conditions. In the first test a linear regression model is fitted to the data in a pixelwise fashion, using three mathematical operators. In the second test a comparison using z-scoring is used based on the assumption that Poisson statistics are valid. For both tests the underlying assumptions are as simple and few as possible. The results are presented as parametric maps of either the three operators or the z-score. The z-score maps can then be thresholded to show the parts of the images which demonstrate change. Three different thresholding methods (naive, adaptive and multiple) are presented: together they cover almost all the needs for separating the signal from the background in the z-score maps. Where the expected size of the signal is known or can be estimated, a spatial correction technique (referred to as the reef correction) can be applied. These tests were applied to flood images used for the quality control of gamma camera uniformity. Simulated data were used to check the validity of the methods. Real data were acquired from four different cameras from two different institutions using a variety of acquisition parameters. The regression model found the bias in all five simulated cases and it also found patterns of unstable regions in real data where visual inspection of the flood images did not show any problems. In comparison the z-map revealed the differences in the simulated images from as low as 1.8 standard deviations from the mean, corresponding to a differential uniformity of 2.2% over the central field of view. In all cases studied, the reef correction increased significantly the sensitivity of the method and in most cases the specificity as well. The two proposed tests can be used either separately or in combination and are capable of showing trends and/or the magnitude of difference between images acquired under similar conditions with high positional and statistical precision. In addition to gamma camera quality control, they could be applied to any pair (or set) of registered planar images to detect subtle changes, e.g. a set of scintigrams or conventional radiographs of a patient before, during and after treatment.

Scaling images using their background ratio. An application in statistical comparisons of images

Comparison of two medical images often requires image scaling as a pre-processing step. This is usually done with the scaling-to-the-mean or scaling-to-the-maximum techniques which, under certain circumstances, in quantitative applications may contribute a significant amount of bias. In this paper, we present a simple scaling method which assumes only that the most predominant values in the corresponding images belong to their background structure. The ratio of the two images to be compared is calculated and its frequency histogram is plotted. The scaling factor is given by the position of the peak in this histogram which belongs to the background structure. The method was tested against the traditional scaling-to-the-mean technique on simulated planar gamma-camera images which were compared using pixelwise statistical parametric tests. Both sensitivity and specificity for each condition were measured over a range of different contrasts and sizes of inhomogeneity for the two scaling techniques. The new method was found to preserve sensitivity in all cases while the traditional technique resulted in significant degradation of sensitivity in certain cases.

Automated CT marker segmentation for image registration in radionuclide therapy

In this paper a novel, automated CT marker segmentation technique for image registration is described. The technique, which is based on analysing each CT slice contour individually, treats the cross sections of the external markers as protrusions of the slice contour. Knowledge-based criteria, using the shape and dimensions of the markers, are defined to enable marker identification and segmentation. Following segmentation, the three-dimensional (3D) markers' centroids are localized using an intensity-weighted algorithm. Finally, image registration is performed using a least-squares fit algorithm. The technique was applied to both simulated and patient studies. The patients were undergoing 131I-mIBG radionuclide therapy with each study comprising several 99mTc single photon emission computed tomography (SPECT) scans and one CT marker scan. The mean residual 3D registration errors (+/- 1 SD) computed for the simulated and patient studies were 1.8 +/- 0.3 mm and 4.3 +/- 0.5 mm respectively.

Continuous angiotensin II infusion increases tumour: normal blood flow ratio in colo-rectal liver metastases

Insufficient blood flow within colo-rectal hepatic metastases is a factor which may limit drug delivery to, and thus the response of, these tumours to regional chemotherapy. Loco-regional flow may be manipulated pharmacologically to enhance the tumour blood flow relative to that of the normal liver. However, as yet, only transient effects have been studied. Patients receiving regional chemotherapy for unresectable hepatic disease were given a 45 min regional infusion of the vasoconstrictor Angiotensin II. Intrahepatic blood flow distribution was assessed serially by Positron Emission Tomography (PET) imaging together with the trapping tracer copper(II) pyruvaldehyde bis(N-4-methylthiosemicarbazone) (Cu-PTSM) labelled using copper-62. Eleven lesions in nine patients were studied, with no adverse effects. Prior to Angiotensin II administration tumour blood flow was generally found to be greater than that of liver (10/11 lesions; 8/9 patients; median TNR 1.3, iqr 0.9-2.5). A significant increase in relative flow to tumour was seen in response to 10 min Angiotensin II infusion in most cases (7/11 lesions; 7/9 patients; median TNR 2.1, iqr 1.4-4.1; P = 0.008), which appeared to be sustained throughout the 45 min infusion period (median TNR 1.85, iqr 1.3-3.8; P = 0.03). These effects were accompanied by transient elevation of mean arterial pressure, but no change in pulse rate. These observations suggest that prolonged regional vasoconstrictor administration could prove useful in the management of unresectable colo-rectal hepatic metastases, and that further development of vascular manipulation to enhance tumour targeting and drug delivery is warranted.

An automated technique for SPECT marker-based image registration in radionuclide therapy

An automated technique for marker-based image registration in radionuclide therapy is described. This technique is based on localization of the centroids of external markers and on establishing correspondence between the individual markers of the two studies to be registered. Localization of the centroids of markers relies on segmenting the markers using iterative thresholding. Thresholding is locally adaptive in order to account for study-dependent conditions (e.g. crossover between adjacent markers and markers with varying radioactive concentrations). Following marker segmentation, the centroids of the markers are computed based on an intensity-weighted method. Finally, prior to the least-squares fit registration, the markers of the two sets are matched to achieve one-to-one correspondence. The technique was applied to both simulated and patient studies resulting in mean residual three-dimensional registration errors (+/- 1SD) of 1.7 +/- 0.1 mm and 3.5 +/- 0.3 mm respectively. The technique was compared with a semi-automated approach and no significant difference was found between the mean residual three-dimensional registration errors (t = 0.281. p = 0.8). This automated marker-based image registration technique provides robust and accurate registration. Although it was developed as part of a programme to generate three-dimensional dose distributions for radionuclide therapy, it could be useful for other clinical applications.

62Cu-PTSM and PET used for the assessment of angiotensin II-induced blood flow changes in patients with colorectal liver metastases

The aim of this study was to establish a quantitative positron emission tomography (PET) method for investigating angiotensin II (AII)-induced changes in blood flow distribution in the liver. This was in order to evaluate the role of vascular manipulation applied to locoregional chemotherapy treatment in patients with colorectal liver metastases. The tracer selected was copper-62 (II) pyruvaldehyde bis-(N4-methyl)thiosemicarbazone (62Cu-PTSM), which exhibits high first-pass extraction and tissue retention following intra-arterial administration. The short half-life of the tracer and its availability from a 62Zn/62Cu generator enabled short-interval repeat PET scans on patients in a single imaging session. Distribution of tracer within the liver was imaged in a single view using a PET camera with rotating large-area detectors. By optimisation of the acquisition protocol, it was possible to acquire sufficient data to produce good-quality images and to quantify tracer uptake with an accuracy of <10%. Reproducibility of the imaging method was assessed in a single patient in whom three consecutive 62Cu-PTSM PET scans were obtained, and in whom no vascular manipulation was performed. Sets of scans (before, during and immediately after a 45-min AII infusion) were obtained in nine patients to assess blood flow changes associated with prolonged vascular manipulation. Significant individual responses, varying in both the magnitude and the duration of flow change, were observed in the majority of cases (7/11 lesions; 7/9 patients). These findings illustrate the potential of 62Cu-PTSM and PET for pharmacological studies. The wide range of individual patient responses to AII infusion suggests that PET blood flow assessment would be of value for selecting patients in whom this procedure may be effective.

Effects of cellular repair and proliferation on targeted radionuclide therapy: a modelling study

A previous targeted radionuclide therapy modelling study has been extended to include the radiobiological effects of cellular repair and proliferation. Dose distributions have been converted into biologically effective dose (BED) distributions using a previously published formulation. With suitable estimated parameters, corrected tumour control probability (TCP) values were derived. The dependence of BED on the physical half-life of the radionuclide was also modelled. Results indicate that the TCP is greater when a shorter physical half-life is employed.

Intercomparison of four reconstruction techniques for positron volume imaging with rotating planar detectors

Four reconstruction techniques for positron volume imaging have been evaluated for scanners based on rotating planar detectors using measured and simulated data. The four techniques compared are backproject then filter (BPF), the 3D reprojection (3D RP) method for 3D filtered backprojection (FBP), Fourier rebinning (FORE) in conjunction with 2D FBP (FORE + 2D FBP) and 3D ordered subsets expectation maximization (3D OSEM). The comparison was based on image resolution and on the trade-off between contrast and noise. In general FORE + 2D FBP offered a better contrast-noise trade-off than 3D RP, whilst 3D RP offered a better trade-off than BPF. Unlike 3D RP, FORE + 2D FBP did not suffer any contrast degradation effect at the edges of the axial field of view, but was unable to take as much advantage from high-accuracy data as the other methods. 3D OSEM gave the best contrast at the expense of greater image noise. BPF, which demonstrated generally inferior contrast-noise behaviour due to use of only a subset of the data, gave more consistent spatial resolution over the field of view than the projection-data based methods, and was best at taking full advantage of high-accuracy data.

Fast accurate iterative reconstruction for low-statistics positron volume imaging

A fast accurate iterative reconstruction (FAIR) method suitable for low-statistics positron volume imaging has been developed. The method, based on the expectation maximization-maximum likelihood (EM-ML) technique, operates on list-mode data rather than histogrammed projection data and can, in just one pass through the data, generate images with the same characteristics as several ML iterations. Use of list-mode data preserves maximum sampling accuracy and implicitly ignores lines of response (LORs) in which no counts were recorded. The method is particularly suited to systems where sampling accuracy can be lost by histogramming events into coarse LOR bins, and also to sparse data situations such as fast whole-body and dynamic imaging where sampling accuracy may be compromised by storage requirements and where reconstruction time can be wasted by including LORs with no counts. The technique can be accelerated by operating on subsets of list-mode data which also allows scope for simultaneous data acquisition and iterative reconstruction. The method is compared with a standard implementation of the EM-ML technique and is shown to offer improved resolution, contrast and noise properties as a direct result of using improved spatial sampling, limited only by hardware specifications.

Thyroid volume measurement in thyrotoxic patients: comparison between ultrasonography and iodine-124 positron emission tomography

The aim of this paper was to compare ultrasound (US) assessment of thyroid volume with that obtained by positron emission tomography (PET), in patients scheduled for adaptive radioiodine therapy, in which 50 Gy was prescribed to the functional PET volume. Firstly a pilot study was performed to ascertain the optimum method for US assessment of thyroid volume. Then 17 comparative measurements of thyroid volume by US and PET were made on 15 patients (two male and thirteen female, ages 28-73 years) with suspected Graves' disease. This comparison showed that in normal sized and enlarged thyroid glands (n=13), the ratio of functional PET to anatomical US volume was approximately 2:3. However, using the same ellipsoid model, PET and US assessment of anatomical volume agreed within the measurement errors. Owing to the presence of nodules and non-uniform distribution of radioiodine, the functional PET volume and anatomical US volume are often not equivalent. If high-resolution emission tomography (e.g. PET) is unavailable, the comparative data presented in this paper could be used to derive the functional volume from the US volume for calculating functional thyroid dose in hyperthyroid patients undergoing radioiodine therapy.

Radiation dosimetry for 131I-mIBG therapy of neuroblastoma

This paper describes the methodology which can be used to determine whole-body, red marrow, blood, bladder, liver, and tumour doses delivered during 131I-mIBG therapy of neuroblastoma. The methodology is based on the Physics Protocol used in a multi-centre study undertaken by the United Kingdom Children's Cancer Study Group (UKCCSG). In this study, the estimates of the doses delivered, using 2.4-12.1 GBq 131I-mIBG, were in the following ranges: whole body, 0.14-0.65 mGy MBq-1; red marrow, 0.17-0.63 mGy MBq-1; blood, 0.04-0.17 mGy MBq-1; bladder, 2.2-5.3 mGy MBq-1; liver, 0.3-1.9 mGy MBq-1; and tumour, 0.2-16.6 mGy MBq-1.

A prototype rotating slat collimator for single photon emission computed tomography

A collimator consisting of a series of highly attenuating parallel slats has been constructed and used in conjunction with a gamma-camera to approximately measure planar projections of a given radionuclide distribution. The enlarged solid angle of acceptance afforded by the slat collimator gave rise to an increased geometric efficiency of between 12 and 28 times that observed with a low-energy high-resolution (LEHR) parallel-hole collimator. When the slats rotated over the face of the detector and the camera gantry turned about the object, sufficient projections were acquired to reconstruct a three-dimensional (3-D) image using the inversion of the 3-D radon transform. The noise behavior of an algorithm for implementing this inversion was studied analytically and the resulting relationship has been verified by computer simulation. The substantially improved geometric efficiency of the slat collimator translated to improvements in reconstructed signal-to-noise ratio (SNR) by, at best, up to a factor of 2.0 with respect to standard parallel-hole collimation. The spatial resolution achieved with the slat collimator was comparable to that obtained with a LEHR collimator and no significant differences were observed in terms of scatter response. Accurate image quantification was hindered by the spatially variant response of the slat collimator.

Absolute 99Tcm-DMSA renal uptake in children: optimum time to scan

Absolute renal uptake was measured at 2, 4 and 6 h in 27 patients in order to determine a more convenient time for uptake compared with the 6-h measurement proposed by other authors. Measurements made at 2 and 4 h would need to be increased by 20 and 6%, respectively, to convert to the value at 6 h. Measurements at 4 h are a reasonable compromise to achieve a high-count/low-background image in a reasonable time and to obtain a good estimate of the predicted 6 h uptake at a single scanning session. The percent renal uptake at 6 h, U(6), can be derived from the percent uptake measured at t hours after injection, U(t), using the following formula: U(6) = U(t) x CF(t), where CF(t) = 1.0 + 0.03 (6 - t) for 4 < or = t < or = 6.

Imaging metastatic testicular germ cell tumours with 18FDG positron emission tomography: prospects for detection and management

The aim of this study was to investigate the role of positron emission tomography (PET) with [18F]fluoro-2-deoxyglucose (18FDG) in metastatic testicular germ cell tumours. Twenty-one patients with stage II-IV testicular germ cell tumours were imaged by PET with a multiwire proportional chamber PET system and 18FDG. Avid 18FDG uptake was seen in metastatic disease from primary seminoma and malignant teratoma. Normal tissue uptake was seen in differentiated teratoma or necrotic, fibrotic tissue. 18FDG standard uptake values and tumour to normal tissue ratios were 6.0 +/- 1.4 and 1.7 +/- 0.4 (mean +/- 1SD), respectively, for malignant tissue. Reduction of 18FDG tumour to normal tissue ratios from pre-treatment to on-treatment scans was predictive of response (n = 3). No significant reduction in 18FDG uptake was seen in patients not responding to therapy (n = 2). These results suggest a role for 18FDG PET in the detection and management of metastatic testicular germ cell tumours.

The experimental evaluation of a prototype rotating slat collimator for planar gamma camera imaging

A collimator consisting of a series of parallel slats has been constructed and used in conjunction with a conventional gamma camera to collect one-dimensional projections of the radioisotope distribution being imaged. With the camera remaining stationary, the collimator was made to rotate continuously over the face of the detector and the projections acquired were used to reconstruct a planar image by the theory of computed tomography. The propagation of noise on image reconstruction was largely offset by the increased geometric efficiency that resulted from the enlarged solid angle of acceptance afforded by the slat collimator. For a uniform disc of activity the signal to noise ratio (SNR) at a point in an image reconstructed by convolution and backprojection is shown to be given by [formula:see text] and Q1(xi) is the one-dimensional filter function in Fourier space. Improved noise behaviour was observed for images acquired with the slat collimator compared to those acquired with a low-energy high-resolution (LEHR) collimator for small distributions of activity. Spatial resolution with the slat collimator was approximately equal to that obtained with an LEHR collimator and improved contrast was observed in images of small hot regions.

Blood flow measurement in extremity soft tissue sarcoma with technetium-99m hexamethyl-propyleneamineoxime and single photon emission computed tomography

Blood flow measurements were made in 28 patients with soft tissue sarcoma of the extremities to investigate the prognostic significance of tumour vascularity. Four patients with benign tumours also underwent blood flow measurement. Mean and maximum tumour blood flow was calculated from technetium-99m hexamethyl-propyleneamineoxime uptake measured using single photon emission computed tomography (SPECT), tumour volume measured from SPECT transaxial image reconstructions and cardiac output assessed with Doppler ultrasonography. Twenty-seven malignant lesions and one benign tumour showed increased uptake of isotope relative to surrounding tissues. Mean sarcoma blood flow varied between 1 and 33 ml per 100 ml tumour per min, and maximum flow between 5 and 57 ml per 100 ml per min. Fourteen patients developed progressive disease during the first year of follow-up. Eight of 11 patients with a high isotope uptake ratio, eight of 12 with a high mean blood flow and eight of 14 with a high maximum blood flow relative to the respective medians for the series showed disease progression.

Dose-response study on thyrotoxic patients undergoing positron emission tomography and radioiodine therapy

With the acknowledged problems associated with assessment of functioning thyroid mass and hence radiation dose, our policy had been to give 75 MBq iodine-131 at 6-monthly intervals to patients with Graves' disease until they became euthyroid. Since positron emission tomography (PET) has been available at this hospital, the radiation dose to the thyroid has been calculated with an accuracy of approximately 20%, the thyroid mass being determined from an iodine-124 PET scan. A dose-response study has been carried out on 65 patients who have received single or cumulative radiation doses of < 80 Gy. The results show that patients who receive a low radiation dose (< 20 Gy) at their first treatment have a high probability of remaining toxic at 12 months. In contrast, patients who receive higher radiation doses (> 40 Gy) at their first treatment have a high probability of control. The probability of becoming euthyroid increases more rapidly with increasing radiation dose than the probability of becoming hypothyroid. Following this dose-response study, a new treatment protocol has been introduced. A 124I PET tracer study prior to 131I therapy will be performed to enable a prescribed thyroid dose of 50 Gy to be delivered to patients with Graves' disease. Further 131I therapy will only be considered if patients are still toxic at 12 months.

Quantitative single-photon emission tomography for tumour blood flow measurement in bronchial carcinoma

A single-photon emission tomography (SPET) technique for the absolute measurement of tumour perfusion is described. Phantom studies have shown that source-background ratios are dependent upon source size and radial position within the phantom. A means of correcting source-background count ratios for these variables has been developed and used to correct tumour-lung ratios obtained in 28 patients with bronchial carcinomas who underwent technetium-99m hexamethylpropyleneamine oxime (99mTc-HMPAO) SPET. On SPET images, the normal lung appears as a relatively homogeneous background. The relationship between 99mTc background concentration (kBq/ml) and counts/pixel was determined from phantom studies and the tumour 99mTc concentration from the background 99mTc concentration and corrected tumour-lung ratio. The total activity of the lipophilic 99mTc-HMPAO species injected was measured. The activity reaching the systemic circulation (Asys) was obtained by subtracting the activity trapped in the pulmonary circulation (obtained from background 99mTc concentration and lung volume). Tumour blood flow may then be calculated from fraction of Asys contained in the tumour provided cardiac output and extraction fraction are known. Blood flow through the central region of tumours ranged from zero to 59.0 (mean 14.1) ml min-1 100 g-1 and through the whole tumour from 0.6 to 68.0 (mean 20.6) ml min-1 100 g-1.

Radiation dose assessment in radioiodine therapy. Dose-response relationships in differentiated thyroid carcinoma using quantitative scanning and PET

Dose-response charts have been constructed to determine the tumouricidal dose for differentiated thyroid carcinoma metastases and thus enable precise activities of radioiodine to be prescribed in order to maximise tumour kill and minimise morbidity. Tumour and normal residual thyroid absorbed doses from radioiodine-131 have been determined with increased precision using a dual-headed whole-body rectilinear scanner with special high-resolution low-sensitivity collimators. Improved accuracy in the estimation of functioning tumour mass has been achieved using positron emission tomography (PET) with a low-cost large area PET camera. Dose-response data have been obtained for 33 patients. Following near-total thyroidectomy and 3.0 GBq 131I, a mean absorbed dose of 410 Gy achieved complete ablation of thyroid remnants in 75% of patients. Patients who had persistent uptake in the thyroid region on subsequent radioiodine scanning had received a mean dose of only 83 Gy. Cumulative absorbed doses in excess of 100 Gy were found to eradicate cervical node metastases. Patients with bone metastases, who generally have a poor prognosis, were found to have received doses of the order of only 20 Gy to the tumour deposits. The dose-response data explain the spectrum of clinical responses to fixed activities of radioiodine. In future, they will enable precise prescription of radioiodine to achieve tumouricidal doses whilst avoiding the morbidity and expense of ineffective therapy.

Phase I/II study of iodine 131 metaiodobenzylguanidine in chemoresistant neuroblastoma: a United Kingdom Children's Cancer Study Group investigation

PURPOSE

The goal of this study was to evaluate the toxicity of iodine 131 metaiodobenzylguanidine (mIBG) in metastatic neuroblastoma.

PATIENTS AND METHODS

A multicenter phase I study of 131I mIBG has been undertaken by the United Kingdom Children's Cancer Study Group (UKCCSG) in children with advanced chemoresistant neuroblastoma. Activity prescription was based on a prescribed whole-body radiation dose, which was established for individual patients by performing an initial tracer investigation with 75 MBq of 131I mIBG. An activity was derived from this pharmacokinetic study that would deliver an initial whole-body-absorbed radiation dose of 1 Gy. Subsequent dose escalations were based on observed toxicity.

RESULTS

Twenty-five patients, aged 1 to 10 years, were treated with prescribed whole-body dose levels of 1.0 Gy (n = 2), 2.0 Gy (n = 13), and 2.5 Gy (n = 10). This necessitated administration of 2.4 to 12.1 GBq of activity. Hematologic, hepatic, kidney, and adrenal toxicity were observed, with bone marrow suppression being the principal dose-limiting toxicity. Bone marrow toxicity increased with prescribed whole-body-absorbed radiation dose, with 80% of patients developing grade 3 or 4 thrombocytopenia at a prescribed whole-body radiation dose of 2.5 Gy. Objective evidence of tumor response was seen in soft tissue (primary or nodal disease), bone, and bone marrow, with an overall response rate of 33% (partial response, n = 8; static disease, n = 9; progressive disease, n = 7).

CONCLUSIONS

This study has established an effective method of activity prescription that predicts subsequent toxicity, with the maximally tolerated dose being sufficient activity to deliver a whole-body-absorbed radiation dose of 2.5 Gy. The objective response rate is comparable to other single agents in chemoresistant neuroblastoma and suggests that 131I mIBG may be a useful method for targeting radiotherapy in metastatic neuroblastoma.

Treatment planning for 131I-mIBG radiotherapy of neural crest tumours using 124I-mIBG positron emission tomography

Patients designated to receive 131I-meta-iodobenzylguanadine (mIBG) for the treatment of neural crest tumours have been scanned with 124I-mIBG using the MUP-PET positron camera. Uptake was detected in tumour sites in lung, liver and abdomen. The tomographic images produced have allowed estimates to be made of the concentration of mIBG in both tumour and normal tissue. From these data it is possible to predict the radiation doses that would be achieved using therapy levels (up to 11 GBq) of 131I-mIBG. The levels of tumour uptake are between 0.5 and 2.0 kBq/g indicating that the radiation doses to tumour would be in the range 3 Gy to 7.5 Gy.

The treatment of resistant neuroblastoma with 131I-mIBG: alternative methods of dose prescription

A UK multi-centre study has been carried out to collect medical and dosimetry data from treatments with 131I-metaiodobenzylguanidine (mIBG) for patients suffering from resistant neuroblastoma. All data have been acquired in a standardised way, following strict physics and clinical protocols. The accuracy of three different methods of dose prescription was studied. The results show that the most accurate method involved the use of an initial tracer study to determine the therapeutic activity required to deliver a predetermined absorbed whole-body (WB) dose.

Monte Carlo modelling of the performance of a rotating slit-collimator for improved planar gamma-camera imaging

Planar imaging with a gamma camera is currently limited by the performance of the collimator. Spatial resolution and sensitivity trade off against each other; it is not possible with conventional parallel-hole collimation to have high geometric sensitivity and at the same time excellent spatial resolution unless field-of-view is sacrificed by using fan- or cone-beam collimators. We propose a rotating slit-collimator which collects one-dimensional projections from which the planar image may be reconstructed by the theory of computed tomography. The performance of such a collimator is modelled by Monte Carlo methods and images are reconstructed by a convolution and backprojection technique. The performance is compared with that of a conventional parallel-hole collimator and it is shown that higher spatial resolution with increased sensitivity is possible with the slit-collimator. For a point source a spatial resolution of some 6 mm at a distance of 100 mm from the collimator with a x7 sensitivity compared with a parallel-hole collimator was achieved. Applications to bone scintigraphy are modelled and an improved performance in hot-spot imaging is demonstrated. The expected performance in cold-spot imaging is analytically investigated. The slit-collimator is not expected to improve cold-spot imaging. Practical design considerations are discussed.

Anatomically derived attenuation coefficients for use in quantitative single photon emission tomography studies of the thorax

Elimination of errors due to poor attenuation correction is an essential part of any quantitative single photon emission tomography (SPET) technique. Attenuation coefficients (mu Tc) for use in attenuation correction of SPET data were determined using technetium 99m and cobalt 57 flood sources and using topographical information obtained from computed tomography (CT) scans and magnetic resonance (MR) images. In patients with carcinoma of the bronchus, the mean attenuation coefficient for 99mTc was 0.096 cm-1 when determined across a transverse section of the thorax at the level of the tumour by means of a 57Co flood source (13 patients) and 0.093 and 0.074 cm-1 as determined from CT scans for points in the centre of the tumour and contralateral normal lung, respectively (21 patients). In 18 patients with breast tumours, the mean attenuation coefficient for 99mTc was 0.110 and 0.076 cm-1 when determined from MRI cross-sections for points in the centre of the tumour and normal contralateral lung, respectively. This indicates significant overcorrection for attenuation when the conventional value of 0.12 cm-1 is used. A value in the range 0.08-0.09 cm-1 would be more appropriate for SPET studies of the thorax. An alternative approach to quantitative region of interest (ROI) analysis is to perform attenuation correction appropriate to the centre of each ROI (using topographical information derived from CT or MRI) on non-attenuation-corrected reconstructions.

Measurements of blood-brain barrier permeability in patients undergoing radiotherapy and chemotherapy for primary cerebral lymphoma

Positron emission tomography (PET) has been used to measure changes in regional blood-brain barrier (BBB) permeability in patients with primary cerebral lymphoma undergoing radiotherapy and chemotherapy. The method employed is to measure the rate of wash-out of a radioactive tracer (68Ga-EDTA) from blood into brain tissue using time-sequence PET imaging. Preliminary studies carried out on patients with more common primary cerebral tumours show that time-activity data are reproducible to approximately 10%. Measurements made in 2 patients with primary cerebral lymphoma treated with initial chemotherapy showed significant changes in permeability in the region of the tumour. Within 5 weeks of the start of treatment, permeability values reached the levels of normal brain. No changes in BBB permeability in normal brain were seen immediately after radiotherapy.

Dosimetry of iodine 131 metaiodobenzylguanidine for treatment of resistant neuroblastoma: results of a UK study

In 1987, the United Kingdom Children's Cancer Study Group (UKCCSG) set up a multi-centre study to investigate the toxicity of iodine 131 metaiodobenzylguanidine (mIBG) in the treatment of resistant neuroblastoma. Since December 1987, 25 children suffering from neuroblastoma have been treated with 131I-mIBG at six UK centres. All centres followed standardised physics and clinical protocols to provide consistent toxicity and dosimetry data. These protocols describe the methods employed for both the tracer study using 131I-mIBG and the subsequent therapy. Whole-body dosimetry calculations were performed on data from the tracer study. The activity administered for therapy was the amount predicted to deliver a predefined whole-body dose. Estimates of doses delivered to various organs during treatment are given in Table 1.

Improved quantification of radionuclide uptake using deconvolution and windowed subtraction techniques for scatter compensation in single photon emission computed tomography

A comparison of two methods of scatter compensation in single photon emission computed tomography (SPECT) imaging is made on the basis of improvements in quantification. The methods, scatter-window subtraction and constrained deconvolution of an average point source response function (PSRF), are described; the theoretical basis of each method is also briefly assessed. Improvements in relative quantification offered by each method are measured by examining ratios of counts in hot cylinders to counts in a slightly radioactive background. The cylinder/background concentration ratio was varied by a factor of five; the sizes of the cylinders remained constant. Keeping the diameter of the cylinders constant allowed for an assessment of the effect of concentration differences, without the conflicting effect of variation in the size of the hot source. Results showed that while both scatter-window subtraction and constrained deconvolution offer quantification improvements in the final image, the method of prereconstruction deconvolution of a planar PSRF from each planar projection is substantially more successful than either scatter-window subtraction or other methods of implementing the deconvolution procedure.

The effects of single dose oral hydralazine on blood flow through human lung tumours

Hydralazine has been shown to reduce tumour blood flow and to potentiate the cytotoxicity of melphalan and bioreductive agents in mice. In order to determine whether such a strategy might have clinical potential, a study was undertaken to investigate the effects of hydralazine on blood flow through human tumours. Twenty-two patients with carcinoma of the bronchus received a single oral dose of hydralazine in the range 25 to 150 mg (0.37-2.86 mg/kg) according to age and acetylator status. Tumour blood flow was assessed by single photon emission computed tomography (SPECT) performed 10 min following intravenous 99Tcm-HMPAO on two occasions 2-8 days apart, the second being performed 60 min after hydralazine administration. In 20 evaluable patients, hydralazine caused a 38% increase in blood flow through the whole tumour (p = 0.007) and a 28% increase in flow through the tumour centre (p = 0.03) with greater increases occurring in patients sustaining greater falls in peripheral resistance. Tumour vascular resistance fell indicating active vasodilation in arterioles supplying tumours. Side-effects due to hydralazine were reported by eight patients.

Thyroid imaging using positron emission tomography--a comparison with ultrasound imaging and conventional scintigraphy in thyrotoxicosis

Forty-six comparative studies were performed on 41 patients with hyperthyroidism. Clinically these comprised two groups: those with Graves' disease, and those with multinodular goitre. All patients underwent an ultrasound examination and positron emission tomography (PET) using 124I, then gamma camera pinhole imaging following their 131I therapy administration. Although the 131I pinhole imaging was not performed for diagnostic purposes, there was good correlation (78% agreement) between it and 124I PET in determining relative lobe size. Hence either imaging modality could be used as an indicator of the relative radiation dose delivered to each thyroid lobe at a macroscopic level. In terms of gland morphology the PET images corresponded well to the high resolution ultrasound images (78% agreement), unlike the pinhole images which correlated poorly (only 28% agreement). The results showed that PET imaging gives better anatomical and physiological detail than 131I pinhole imaging. In 77% of cases where the pinhole image showed a uniform distribution of radioisotope, the improved spatial resolution of the PET images revealed non-homogeneous distribution indicating a non-uniform distribution of radiation dose. Since all dosimetry calculations are based on the assumption of uniform distribution of radioiodine, this non-uniformity could possibly have important consequences in the outcome of radioiodine therapy in thyrotoxicosis.

An evaluation of 99mTc-HMPAO uptake in cerebral gliomas--a comparison with X-ray CT

Nineteen patients with biopsy-proven cerebral gliomas were studied with 99mTc-HMPAO single photon emission tomography (SPECT) imaging and X-ray computed tomography (CT). The uptake of 99mTc-HMPAO was correlated with tumour size and morphology as shown by X-ray CT, and overall patient survival. It appears that uptake of 99mTc-HMPAO is associated with larger, ill-defined tumours and was an adverse factor in patient survival. In those tumours with normal or increased uptake, 99mTc-HMPAO imaging is useful in distinguishing the tumour margin from surrounding oedema.

The performance of a multiwire proportional chamber positron camera for clinical use

The Rutherford Appleton Laboratory clinical positron camera consists of two opposed multiwire proportional chambers (MWPCS) mounted on a rotating gantry capable of performing routine nuclear medicine studies. The system has operated since the end of 1986 with complete reliability. It has a sensitivity of 37 kcps MBq-1 cm3 per axial cm for a 20 cm diameter uniformly filled cylinder of activity. The best spatial resolution obtainable is about 6 mm, although in practice images are smoothed in order to reduce statistical noise with a corresponding decrease in resolution. Cross-plane rays are utilised during image reconstruction and the resulting three-dimensional images exhibit the same spatial resolution in three orthogonal directions over a large cylindrical field of view 15 cm high by 30 cm in diameter. The maximum data-taking rate is limited to 1.5 to 2 kcps at present due to deadtime in the read-out system. The performance of the system is described with particular reference to the problems of imaging with the large fractions of random and scattered events which are a consequence of using large-area detectors. Images of phantoms and patients are presented and proposed modifications to the camera are described.

Radiation dose assessment in radioiodine therapy. 2. Practical implementation using quantitative scanning and PET, with initial results on thyroid carcinoma

We have designed special high-resolution, low-sensitivity collimators for a dual-headed whole-body scanner for imaging and quantifying therapy levels of iodine-131. In addition, we have used positron emission tomography (PET) with a low-cost large-area PET camera to achieve improved accuracy in the estimate to tumor mass. The physical performance of these two imaging systems is described. In order to illustrate the practical implementation of these systems for the assessment of radiation dose to normal and tumour tissue during radioiodine therapy, three clinical examples are reported, and a summary of the initial clinical results obtained from 16 patients with carcinoma of the thyroid is presented. The dose to normal thyroid remnants for patients undergoing ablation ranged from 16 to 400 Gy, while the dose to involved neck nodes ranged from 2.5 to 33 Gy for patients undergoing post-ablation radioiodine therapy. In one patient with distant metastasis in the spine, a dose of 100 Gy was achieved. The techniques described in this paper can be used to determine if sufficient activity can be accumulated in tumours to provide a therapeutic effect while minimising irradiation of normal tissues by avoiding administrations which do not provide tumouricidal radiation doses.

Technetium-99m HMPAO and SPECT in the assessment of blood flow in human lung tumours

In order to assess the blood flow patterns through human lung tumours, 20 patients received 400-750 MBq 99TcmHMPAO intravenously 10 min before single photon emission computed tomography (SPECT). Ratios of uptake in the whole tumour relative to normal lung ranged from 0.35 to 1.53 (mean 1.01) with eight tumours showing less uptake than normal lung and ten showing greater uptake. In one patient the tumour was not distinguishable from surrounding lung and in another a large pleural effusion prevented evaluation. Tumour: lung ratios for central tumour regions ranged from 0 to 1.83 (mean 0.80) with 13 showing lower uptake than normal lung and five showing greater uptake. Duplicate scans were performed in eight patients demonstrating satisfactory reproducibility. This technique provides a simple and reproducible method for the assessment of tumour blood flow.

Radiation dose assessments in radioiodine (131I) therapy. 1. The necessity for in vivo quantitation and dosimetry in the treatment of carcinoma of the thyroid

In order to destroy thyroid cancer metastases by radioiodine an average tissue dose of 80-300 Gy is needed. Such high doses can be expected, following the administration of the conventional 5.5 GBq of 131I, only if both the percentage uptake per gram in the target tissue and the effective half life of the radioiodine in it are higher than well-defined threshold values, and if every dimension of the tissue exceeds several millimeters. The fulfillment of such favourable conditions in actual clinical cases can only be confirmed by in vivo quantitation of the absorbed dose achieved as a result of the administration of radioiodine.

Image quantification with a large area multiwire proportional chamber positron camera (MUP-PET)

A large area multi wire proportional chamber positron camera system is under evaluation for clinical Nuclear Medicine investigation using isotopes both from the house generators (68Ga) and from remote cyclotrons (18F, 124I, 64Cu and 81Rb). Images are reconstructed using a fully three dimensional (3D) algorithm and exhibit equal resolution in all three orthogonal directions. The axial field of view of 15 cm is large enough to allow coverage of the whole brain. This paper discusses the performance of the MUP-PET system with particular emphasis on quantification and the development of an attenuation correction scheme for a large area detector system. The distribution of scattered and accidental events is investigated and observed to contribute a fairly uniform background to the image. Dead time correction factors are calculated from an analytical expression obtained by considering the various sources of dead time present in the system. Following correction for attenuation, accidentals, scattered events and dead time, reconstructed phantom images show a strong linear correlation (r = 0.998) between count density and regional isotope concentration. The extension of the methods to clinical studies is discussed.

A simple method of producing parenchymal renograms using parametric imaging

Standard dynamic renography procedures are used to produce time-activity curves for defined regions of interest around the kidneys. Further analysis of the acquired data using Fitted Retention and Excretion Equations (FREE) to determine single kidney glomerular filtration rate (GFR) and kidney transit time is performed at the Royal Marsden Hospital, Sutton. The disadvantage of these methods is that the renal function is assessed for the whole kidney and not on a regional basis. The aim of this study was to use parametric imaging to extend the standard analysis so that the parenchyma could be identified and analysed separately from the collecting system. The parenchyma pixel time-activity curves have shorter time-to-peak (TTP) values than the collecting system. TTP parametric images were used to set limits on the TTP values and thus define the parenchyma pixels. The time-activity curves from the parenchyma were further analysed using FREE to produce parametric images of both GFR and effective blood volume. Finally, the original dynamic data were re-analysed so that time-activity and FREE curves were generated for the parenchyma alone. A study of 15 patients (30 kidneys) showed that the extended analysis indicated parenchymal isotope clearance. By comparing the time-activity and FREE curves of the whole kidney with those from the parenchyma alone, a more useful analysis of parenchymal function was obtained. This proved to be particularly useful for patients with pelvi-ureteric junction (PUJ) obstruction.

Initial experience with Tc-99m-HM-PAO in the study of brain tumors

A preliminary study of the distribution of the 99mTc complex of hexamethylpropylene amine oxime (HM-PAO) in 12 patients with brain neoplasms before, during, and after radiotherapy has been performed. Untreated brain tumors were found to exhibit a range of 99mTc-HM-PAO uptake, varying from areas of markedly increased isotope activity to photopenic areas, when compared to normal brain tissue. A ratio of 99mTc-HM-PAO tumor uptake to contralateral normal tissue uptake was calculated prior to and during radiotherapy. This ratio tended to return towards unity in lesions responding to therapy. A predictable alteration in whole brain 99mTc-HM-PAO uptake during radiotherapy was not demonstrated. Unlike the radiolabeled amines, 99mTc-HM-PAO localizes in primary tumors, probably indicating that its uptake mechanism is independent of non specific amine receptors. 99mTc-HM-PAO may be useful in the study of brain tumor physiology and response to therapy.

Deconvolution of emission tomographic data: a clinical evaluation

A method of improving the quality of images in single photon emission computed tomography (SPECT) is demonstrated using transaxial images of the liver and brain. Deconvolution of the nuclear medicine data by a point source response function (PSRF) acquired previously in a scattering medium attempts to compensate for scattered radiation within the patient. The average geometric response of the collimator of the gamma camera is also compensated for with this technique. Three patients with known metastatic lesions in the liver and three with primary lesions in the brain were imaged. Clinical assessment of reconstructed slices both before and after deconvolution demonstrates that compensating for the effects of scatter and of collimator blurring leads to enhanced detail of pathological lesions. In all cases, cold lesions seen prior to deconvolution were enhanced in detail and, in addition, new lesions were seen with this technique.

A comparison of deconvolution and windowed subtraction techniques for scatter compensation in SPECT

Three procedures for the removal of Compton-scattered data in SPECT by constrained deconvolution are presented. The first is a deconvolution of a 2-D measured PSRF containing scatter from a single reconstructed transaxial image; the second is a deconvolution of a 2-D measured point-source response function (PSRF) from each frame of projection data prior to reconstruction; the third involves deconvolution of a 3-D measured PSRF from a stack of reconstructed slices. Results of applying these procedures to data obtained from a phantom containing cold cylinders and to data from a cold spot-resolution phantom are presented and are shown to be superior to the results of correcting for scatter by scatter-window substraction. Both 3-D deconvolution from reconstructed images and 2-D deconvolution from projection data show major improvements in image contrast, resolution, and quantitation. Improvements are especially marked for small (1.0-3.0 cm) cold sources.

Three-dimensional display of data obtained by single photon emission computed tomography

Three-dimensional display of radioactive distributions has been achieved using data generated by single photon emission computed tomography (SPECT). Computer techniques for constructing images are presented and the special problems associated with manipulating SPECT data are reviewed. The potential role of three-dimensional imaging in nuclear medicine is discussed and SPECT studies from several illustrative cases are presented.

Measurement of radiation dose to the thyroid using positron emission tomography

Measurements of the functioning volume of thyroid tissue have been made in 22 patients undergoing radioiodine therapy for thyrotoxicosis, using a prototype multiwire proportional counter positron camera. Tomographic images were produced of the distribution of 124I in the thyroid. Functioning volumes were found to be in the range 21-79 cm3 with volume errors of the order of +/- 4% to +/- 14%. Radioiodine uptake varied from 28% to 98%. Using a value of 6 days for the effective half-life of radioiodine in hyperactive thyroids, radiation doses from a standard therapy administration of 75 MBq of 131I varied from 11 to 48 Gy (compared with a recommended 50-70 Gy). In five cases PET imaging showed a non-uniform distribution of radioiodine in thyroids thought to have uniform uptake from conventional pinhole scintigraphy.