Automated size-specific dose estimates framework in thoracic CT using convolutional neural network based on U-Net model.

PURPOSE: This study aimed to develop an automated method that uses a convolutional neural network (CNN) for calculating size-specific dose estimates (SSDEs) based on the corrected effective diameter (Deff corr ) in thoracic computed tomography (CT). METHODS: Transaxial images obtained from 108 adult patients who underwent non-contrast thoracic CT scans were analyzed. To calculate the Deff corr according to Mihailidis et al., the average relative electron densities for lung, bone, and other tissues were used to correct the lateral and anterior-posterior dimensions. The CNN architecture based on the U-Net algorithm was used for automated segmentation of three classes of tissues and the background region to calculate dimensions and Deff corr values. Then, 108 thoracic CT images and generated segmentation masks were used for network training. The water-equivalent diameter (Dw ) was determined according to the American Association of Physicists in Medicine Task Group 220. Linear regression and Bland-Altman analysis were performed to determine the correlations between SSDEDeff corr(automated) , SSDEDeff corr(manual) , and SSDEDw . RESULTS: High agreement was obtained between the manual and automated methods for calculating the Deff corr SSDE. The mean values for the SSDEDeff corr(manual) , SSDEDw , and SSDEDeff corr(automated) were 14.3 ± 2.1 mGy, 14.6 ± 2.2 mGy, and 14.5 ± 2.4 mGy, respectively. The U-Net model was successfully trained and used to accurately predict SSDEs, with results comparable to manual-labeling results. CONCLUSION: The proposed automated framework using a CNN offers a reliable and efficient solution for determining the Deff corr SSDE in thoracic CT.

Evaluation of patient radiation dose and risk of cancer from CT examinations.

Computed tomography (CT) examinations have been increasingly requested and become the major sources of patient exposure. The cancer risk from CT scans is contingent upon the amount of absorbed dose of organs. This study aims to determine the organ doses and risk of cancer incidence and mortality from CT examinations at high dose (cumulative effective dose, CED ≥ 100 mSv) in a single day to low dose (CED < 100 mSv) from common CT procedures. Data were gathered from two academic centers of patients aged 15 to 75 years old performed CT examinations during the period of 5 years. CED and organ dose were calculated using Monte Carlo simulation software. Lifetime attributable risk (LAR) was determined following Biological Effects of Ionizing Radiation (BEIR) VII report based on life table and baseline cancer rates of Thai population. At high dose, the highest LAR for breast cancer incidence in young female was 82 per 100,000 exposed patients with breast dose of 148 mGy (CT whole abdomen). The highest LAR for liver cancer incidence in male patient was 72 per 100,000 with liver dose of 133 mGy (multiple CT scans). At low dose, the highest average LAR for breast cancer incidence in young female was 23 per 100,000 while for liver cancer incidence in male patients was 22 per 100,000 (CTA whole aorta). Even though the LAR of cancer incidence and mortality was less than 100 per 100,000, they should not be neglected. The risk of cancer incidence may be increased in later life, particularly in young patients.

Comparison of absorbed doses to the tumoral and non-tumoral liver in HCC patients undergoing 99mTc-MAA and 90Y-microspheres radioembolization.

PURPOSE: This study aimed to determine the absorbed doses in the tumoral-liver and non-tumoral liver of hepatocellular carcinoma (HCC) patients undergoing radioembolization with Yttrium-90 (90Y) resin microspheres, and compared with those derived from 99mTc-MAA using the partition model. METHODS: A total of 42 HCC patients (28 males and 14 females, mean age 65 ± 11.51 years) who received 45 treatment sessions with 90Y-microspheres between 2016 and 2021 were included. Pre-treatment 99mTc-MAA and post-treatment 90Y-bremsstrahlung SPECT/CT were acquired for each patient. Semi-automated segmentation of regions of interest (ROIs) was performed using MIM Encore software to determine the tumor-liver ratio (TLR) encompassing the liver volume, tumoral-liver, and lungs, and verified by both nuclear medicine physician and interventional radiologist. A partition dosimetry model was used to estimate the administered activity of 90Y-microspheres and the absorbed doses to the tumoral-liver and non-tumoral liver. The student's paired t test and Bland-Altman plot were used for the statistical analysis. RESULTS: The mean TLR values obtained from 99mTc-MAA SPECT/CT and 90Y-bremsstrahlung SPECT/CT were 4.78 ± 3.51 and 2.73 ± 1.18, respectively. The mean planning administered activity of 90Y-microspheres based on 99mTc-MAA SPECT/CT was 1.56 ± 0.80 GBq, while the implanted administered activity was 2.53 ± 1.23 GBq (p value < 0.001). The mean absorbed doses in the tumoral-liver estimated from 99mTc-MAA and 90Y-bremsstrahlung SPECT/CT were 127.44 ± 4.36 Gy and 135.98 ± 6.30 Gy, respectively. The corresponding mean absorbed doses in the non-tumoral liver were 34.61 ± 13.93 Gy and 55.04 ± 16.36 Gy. CONCLUSION: This study provides evidence that the administered activity of 90Y-microspheres, as estimated from 90Y-bremsstrahlung SPECT/CT, was significantly higher than that estimated from 99mTc-MAA SPECT/CT resulted in increased absorbed doses in both the tumoral-liver and non-tumoral liver. However, 99mTc-MAA SPECT/CT remains a valuable planning tool for predicting the distribution of 90Y-microspheres in liver cancer treatment.

Comparative post-therapeutic dosimetry between 2D planar-based and hybrid-based methods for personalized Lu-177 treatment.

PURPOSE: This study aims to compare the calculated absorbed dose in target organs and tumors obtained using the different imaging protocols and the calculation methodologies implemented by HERMES HybridViewer dosimetry software for 177Lu-PSMA I&T and 177Lu-DOTATATE therapy. METHODS: Multiple time-point whole-body planar images and one SPECT/CT image were acquired from 18 patients including 177Lu-PSMA I&T (13 patients) and 177Lu-DOTATATE treatment (5 patients) after administration of 3.80-8.58 GBq injected activity. The regions of interest were drawn in the whole body, kidneys, liver, urinary bladder, salivary glands, and tumors to determine the time-integrated activity (TIA) in source organs. Absorbed doses in target organs were calculated according to the Medical Internal Radiation Dose (MIRD) scheme using the HERMES HybridViewer dosimetry integrated with OLINDA/EXM V.2.1 that utilizes the non-uniform rational B-splines (NURBS) for computational digital phantom. RESULTS: The planar-based dosimetry showed a higher dose per injected activity compared to the hybrid-based dosimetry, primarily due to organ overlap. The highest difference in absorbed dose between the imaging scenarios was observed in the spleen with a variation of up to 51.6%, while the difference for other target organs and tumors was less than 40%. CONCLUSION: The dosimetry calculation derived from the 2D planar-based method consistently demonstrates a significantly higher absorbed dose in organs and tumors compared with the hybrid-based method. However, the hybrid method outperforms the planar method in terms of tumor visualization and overlap-free organ delineation.

Dual-time-point dynamic 68Ga-PSMA-11 PET/CT for parametric imaging generation in prostate cancer.

PURPOSE: To investigate the optimal dual-time-point (DTP) approaches using dynamic 68Ga-PSMA-11 PET/CT imaging to generate parametric images for prostate cancer patients. METHODS: Fifteen patients with prostate cancer were intravenously administered 68Ga-PSMA-11 of 181.9 ± 47.2 MBq, followed by an immediate 60 min dynamic PET/CT scan. List-mode data were reconstructed into 25 timeframes (6 × 10 s, 8 × 30 s, and 11 × 300 s) and corrected for motion and partial volume effect. DTP parametric images were generated using different interval time points of 5 min and 10 min, with a minimum of 30 min time interval. Net influx rates (Ki) were calculated through the fitting of a single irreversible two-tissue compartmental model. Intraclass correlation coefficient (ICC) values between DTP protocols and 60 min Ki were obtained. Lesion-to-background ratios (LBRs) of Ki and standardized uptake value (SUV) images in each DTP protocol were determined. RESULTS: The DTP protocol of 5-10 min with a 40-45 min interval showed the highest ICC of 0.988 compared with the 60 min Ki, whereas the ICC values for the intervals of 0-5 min with 55-60 min and 0-10 min with 50-60 min were 0.941. The LBRs of the 60 min Ki, 5-10 min with 40-45 min Ki, 0-5 min with 55-60 min Ki, 0-10 min with 50-60 min Ki, SUVmean, and SUVmax images were 29.53 ± 27.33, 13.05 ± 15.28, 45.15 ± 53.11, 45.52 ± 70.31, 19.77 ± 23.43, and 25.06 ± 30.07, respectively. CONCLUSION: The 0-5 min with 55-60 min DTP parametric imaging exhibits a comparable Ki to 60 min parametric imaging and remarkable image quality and contrast than SUV imaging, enhancing prostate cancer diagnosis while maintaining time efficiency.

Bone turnover prediction in patients with chronic kidney disease (CKD) undergoing hemodialysis using shortened dynamic 18F-NaF PET/CT Ki-Patlak.

This study investigated whether Ki-Patlak derived from a shortened scan time for dynamic 18F-NaF PET/CT in chronic kidney disease (CKD) patients undergoing hemodialysis can provide predictive accuracy comparable to that obtained from a longer scan. Twenty-seven patients on chronic hemodialysis, involving a total of 42 scans between December 2021 and August 2023 were recruited. Dynamic 18F-NaF PET/CT scans, lasting 60-90 min, were immediately acquired post-injection, covering the mid-twelfth thoracic vertebra to the pelvis region. Ki-Patlak analysis was performed on bone time-activity curves at 15, 30, 45, 60, and 90 min in the lumbar spine (L1-L4) and both anterior iliac crests. Spearman's rank correlation (rs) and interclass correlation coefficient were used to assess the correlation and agreement of Ki-Patlak between shortened and standard scan times. Bone-specific alkaline phosphatase (BsAP) and tartrate-resistant acid phosphatase isoform 5b (TRAP5b) were tested for their correlation with individual Ki-Patlak. Strong correlations and good agreement were observed between Ki-Patlak values from shortened 30-min scans and longer 60-90-min scans in both lumbar spine (rs = 0.858, p < 0.001) and anterior iliac crest regions (rs = 0.850, p < 0.001). The correlation between BsAP and Ki-Patlak in the anterior iliac crests was weak and statistically insignificant. This finding suggests that a proposed shortened dynamic 18F-NaF PET/CT scan is effective in assessing bone metabolic flux in CKD patients undergoing hemodialysis, offering a non-invasive alternative approach for bone turnover prediction.

A Novel method of seizure onset zone localization by serial Tc-99 m ECD brain perfusion SPECT clearance patterns.

In this prospective study, we postulate that there is a difference between clearance of [99mTc]Tc- ethyl cysteinate dimer (ECD) in the seizure onset zone (SOZ) and other brain areas and thus SOZ localization by clearance patterns might become a potential novel method for SOZ localization in epilepsy. The parametric images of brain ECD clearance were generated by linear regression model analysis from serial brain SPECT scans from 30 to 240 min after ECD injection (7-times point) in 7 patients with drug-resistant epilepsy and 3 normal volunteers. Clearance patterns of the SOZ confirmed by good surgical outcome or consensus with other investigations were analyzed quantitatively and semi-quantitatively by visual grading (slower or faster washout than contralateral brain regions). The average [99mTc]Tc-ECD clearance rates of SOZs were + 1.08% ± 2.57%/hr (wash in), -7.02% ± 2.56%/hr (washout), and -5.37% ± 1.71%/hr (washout) in ictal, aura and interictal states, respectively. Paired t-tests between the SOZ and contralateral regions showed statistically significant difference (p = 0.039 in interictal state). Clearance patterns that can define the SOZs were 1) wash in and slow washout on ictal slope, 2) fast washout on aura slope and interictal slope with 100% (6/6), 100% (2/2) and 75% (6/8) localization using ictal, aura, and interictal slope maps, respectively. Our study provided the evidence that clearance pattern methods are potential additive diagnostic tools for SOZ localization when routine one-time point SPECT are unable to define the SOZ.

Patient dosimetry of 177Lu-PSMA I&T in metastatic prostate cancer treatment: the experience in Thailand.

PURPOSE: This study aimed to determine the radiation dosimetry for 177Lu-PSMA imaging and therapy (I&T) in Thai patients who were treated for metastatic prostate cancer. METHODS: Whole-body planar images acquired at immediately, 4 and 24 h after 177Lu-PSMA I&T injection (range 4.44-8.51 GBq) were collected from 12 treatment cycles of 8 prostate cancer patients. Region of interests (ROIs) were manually contoured on the whole-body, liver, spleen, urinary bladder, lacrimal glands, parotid, and submandibular glands to determine time-integrated activity (TIA) in source organs and fitted time-activity curves using mono-exponential extrapolation. The S values calculated utilizing non-uniform rational B-splines (NURBS) computational phantoms were extracted from the OLINDA/EXM v. 2.0 to calculate the absorbed dose coefficient in target organs according to the Medical Internal Radiation Dose (MIRD) scheme. The absorbed doses to bone marrow were estimated using the planar two-compartment image-based method by separating the high-uptake and low-uptake compartment. The spherical model was used to calculate the lacrimal gland absorbed doses. RESULTS: Mean absorbed dose coefficients to the kidneys, bone marrow, liver, urinary bladder, spleen, lacrimal glands, parotid, and submandibular glands were 0.81 ± 0.24, 0.02 ± 0.01, 0.13 ± 0.10, 0.27 ± 0.25, 0.16 ± 0.07, 3.62 ± 1.78, 0.21 ± 0.14, and 0.09 ± 0.07 Gy/GBq, respectively. Dose constraints for the kidneys (23 Gy) and bone marrow (2 Gy) were not reached in any patients. The absorbed dose in lacrimal glands calculated by the NURBS computational phantoms was slightly lower than the calculation based on the Cristy-Eckerman computational phantoms using OLINDA/EXM v. 1.0 by 6.37 ± 0.14%. CONCLUSION: Dosimetry results in this study suggested that 177Lu-PSMA I&T treatment with higher activities and more cycles is possible without the risk of damaging normal organs in prostate cancer patients.

Pharmacokinetic Modeling of 18F-FDOPA PET in the Human Brain for Early Parkinson's Disease

Objectives: Early detection is essential for the treatment approaches of Parkinson's disease (PD). Clinical criteria alone may be insufficient to distinguish early PD from other conditions. This study aimed to investigate the transfer rate constants of 6-18F-fluoro-L-dopa (18F-FDOPA) in positron emission tomography (PET) brain images as a sensitive parameter to detect early PD. Methods: Retrospective 18F-FDOPA PET data of five patients with early PD were collected. PET data were acquired for 90 min after intravenous injection of 306-379 MBq 18F-FDOPA, and reconstructed into a series of 18 five-minute frames. Reoriented PET images were coregistered and normalized with the PET brain template on the statistical parametric mapping. The 18F-FDOPA activity concentrations were measured in the striatum, caudate, and putamen on both sides: Contralateral (as PD) and ipsilateral (as control) to the main motor symptoms. The pharmacokinetic model was generated using the SAAM II simulation software. The transfer rate constants across the blood-brain barrier (forward, K1 and reverse, k2) and decarboxylation rate constants (k3) were estimated in these regions. Results: The activity uptakes in the contralateral striatum (0.0323%±0.0091%) and putamen (0.0169%±0.0054%) were significantly lower than the control (0.0353%±0.0086%, 0.0199%±0.0054%, respectively). The K1 and k3 were significantly lower in the contralateral striatum and putamen (p<0.05). There were no significant differences in any transfer rate constants in the caudate. Conclusion: The transfer rate constants (K1 and k3) of 18F-FDOPA on the contralateral striatum and putamen were significantly lower than the control. These biokinetic data could be potential indicators for quantitative detection of early PD diagnosis.

Renal 99mTc-DMSA pharmacokinetics in pediatric patients.

99mTc-DMSA is one of the most commonly used pediatric nuclear medicine imaging agents. Nevertheless, there are no pharmacokinetic (PK) models for 99mTc-DMSA in children, and currently available pediatric dose estimates for 99mTc-DMSA use pediatric S values with PK data derived from adults. Furthermore, the adult PK data were collected in the mid-70's using quantification techniques and instrumentation available at the time. Using pediatric imaging data for DMSA, we have obtained kinetic parameters for DMSA that differ from those applicable to adults. METHODS: We obtained patient data from a retrospective re-evaluation of clinically collected pediatric SPECT images of 99mTc-DMSA in 54 pediatric patients from Boston's Children Hospital (BCH), ranging in age from 1 to 16 years old. These were supplemented by prospective data from twenty-three pediatric patients (age range: 4 months to 6 years old). RESULTS: In pediatric patients, the plateau phase in fractional kidney uptake occurs at a fractional uptake value closer to 0.3 than the value of 0.5 reported by the International Commission on Radiological Protection (ICRP) for adult patients. This leads to a 27% lower time-integrated activity coefficient in pediatric patients than in adults. Over the age range examined, no age dependency in uptake fraction at the clinical imaging time was observed. Female pediatric patients had a 17% higher fractional kidney uptake at the clinical imaging time than males (P < 0.001). CONCLUSIONS: Pediatric 99mTc-DMSA kinetics differ from those reported for adults and should be considered in pediatric patient dosimetry. Alternatively, the differences obtained in this study could reflect improved quantification methods and the need to re-examine DMSA kinetics in adults.

Records and redescription of a mygalomorph spider genus ignored for over 100 years with a new species: the genus Atmetochilus Simon, 1887 (Araneae, Nemesiidae) in Thailand.

The wishbone spider of genus Atmetochilus Simon, 1887 (Nemesiidae: Bemmerinae) containing six species and is known from Asia. We describe a new Thai species, Atmetochilus songsangchotei sp. nov. Since the holotype of the type species of Atmetochilus, A. fossor, is presumed lost, we describe fresh specimens (ARA-2018-132) and designate one as as neotype; similarly, because a holotype (only syntypes) of Atmetochilus atriceps in NHMUK was not designated, we chose the most intact specimen of the syntypes (NHMUK 1895.9.21.16) and designate it as lectotype. From a photo of the type specimen of Atmetochilus bifidus (sternum) from ZSI, we transfer Atmetochilus bifidus back to Damarchus.

Determination of radiation dose and low-dose protocol for digital chest tomosynthesis using radiophotoluminescent (RPL) glass dosimeters.

PURPOSE: This study aimed to determine a low-dose protocol for digital chest tomosynthesis (DTS). METHODS: Five simulated nodules with a CT number of approximately 100 HU with size diameter of 3, 5, 8, 10, and 12 mm were inserted into an anthropomorphic chest phantom (N1 Lungman model), and then scanned by DTS system (Definium 8000) with varying tube voltage, copper filter thickness, and dose ratio. Three radiophotoluminescent (RPL) glass dosimeters, type GD-352 M with a dimension of 1.5 × 12 mm, were used to measure the entrance surface air kerma (ESAK) in each protocol. The effective dose (ED) was calculated using the recorded total dose-area-product (DAP). The signal-to-noise ratio (SNR) was determined for qualitative image quality evaluation. The image criteria and nodule detection capability were scored by two experienced radiologists. The selected low-dose protocol was further applied in a clinical study with 30 pulmonary nodule follow-up patients. RESULTS: The average ESAK obtained from the standard default protocol was 1.68 ± 0.15 mGy, while an ESAK of 0.47 ± 0.02 mGy was found for a low-dose protocol. The EDs for the default and low-dose protocols were 313.98 ± 0.72 µSv and 100.55 ± 0.28 µSv, respectively. There were small non-significant differences in the image criteria and nodule detection scoring between the low-dose and default protocols interpreted by two radiologists. The effective dose of 98.87 ± 0.08 µSv was obtained in clinical study after applying the low-dose protocol. CONCLUSIONS: The low-dose protocol obtained in this study can substantially reduce radiation dose while preserving an acceptable image quality compared to the standard protocol.

Re-evaluation of pediatric 18F-FDG dosimetry: Cristy-Eckerman versus UF/NCI hybrid computational phantoms.

Because of the concerns associated with radiation exposure at a young age, there is an increased interest in pediatric absorbed dose estimates for imaging agents. Almost all reported pediatric absorbed dose estimates, however, have been determined using adult pharmacokinetic data with radionuclide S values that take into account the anatomical differences between adults and children based upon the older Cristy-Eckerman (C-E) stylized phantoms. In this work, we use pediatric model-derived pharmacokinetics to compare absorbed dose and effective dose estimates for 18F-FDG in pediatric patients using S values generated from two different geometries of computational phantoms. Time-integrated activity coefficients of 18F-FDG in brain, lungs, heart wall, kidneys and liver, retrospectively, calculated from 35 pediatric patients at the Boston's Children Hospital were used. The absorbed dose calculation was performed in accordance with the Medical Internal Radiation Dose method using S values generated from the University of Florida/National Cancer Institute (UF/NCI) hybrid phantoms, as well as those from C-E stylized computational phantoms. The effective dose was computed using tissue-weighting factors from ICRP Publication 60 and ICRP Publication 103 for the C-E and UF/NCI, respectively. Substantial differences in the absorbed dose estimates between UF/NCI hybrid pediatric phantoms and the C-E stylized phantoms were found for the lungs, ovaries, red bone marrow and urinary bladder wall. Large discrepancies in the calculated dose values were observed in the bone marrow; ranging between -26% to +199%. The effective doses computed by the UF/NCI hybrid phantom S values were slightly different than those seen using the C-E stylized phantoms with percent differences of -0.7%, 2.9% and 2.5% for a newborn, 1 year old and 5 year old, respectively. Differences in anatomical modeling features among computational phantoms used to perform Monte Carlo-based photon and electron transport simulations for 18F, and very likely for other radionuclides, impact internal organ dosimetry computations for pediatric nuclear medicine studies.

Dose Estimation in Pediatric Nuclear Medicine.

The practice of nuclear medicine in children is well established for imaging practically all physiologic systems but particularly in the fields of oncology, neurology, urology, and orthopedics. Pediatric nuclear medicine yields images of physiologic and molecular processes that can provide essential diagnostic information to the clinician. However, nuclear medicine involves the administration of radiopharmaceuticals that expose the patient to ionizing radiation and children are thought to be at a higher risk for adverse effects from radiation exposure than adults. Therefore it may be considered prudent to take extra care to optimize the radiation dose associated with pediatric nuclear medicine. This requires a solid understanding of the dosimetry associated with the administration of radiopharmaceuticals in children. Models for estimating the internal radiation dose from radiopharmaceuticals have been developed by the Medical Internal Radiation Dosimetry Committee of the Society of Nuclear Medicine and Molecular Imaging and other groups. But to use these models accurately in children, better pharmacokinetic data for the radiopharmaceuticals and anatomical models specifically for children need to be developed. The use of CT in the context of hybrid imaging has also increased significantly in the past 15 years, and thus CT dosimetry as it applies to children needs to be better understood. The concept of effective dose has been used to compare different practices involving radiation on a dosimetric level, but this approach may not be appropriate when applied to a population of children of different ages as the radiosensitivity weights utilized in the calculation of effective dose are not specific to children and may vary as a function of age on an organ-by-organ bias. As these gaps in knowledge of dosimetry and radiation risk as they apply to children are filled, more accurate models can be developed that allow for better approaches to dose optimization. In turn, this will lead to an overall improvement in the practice of pediatric nuclear medicine by providing excellent diagnostic image quality at the lowest radiation dose possible.

Pharmacokinetic modeling of [(18)F]fluorodeoxyglucose (FDG) for premature infants, and newborns through 5-year-olds.

BACKGROUND: Absorbed dose estimates for pediatric patients require pharmacokinetics that are, to the extent possible, age-specific. Such age-specific pharmacokinetic data are lacking for many of the diagnostic agents typically used in pediatric imaging. We have developed a pharmacokinetic model of [(18)F]fluorodeoxyglucose (FDG) applicable to premature infants and to 0- (newborns) to 5-year-old patients, which may be used to generate model-derived time-integrated activity coefficients and absorbed dose calculations for these patients. METHODS: The FDG compartmental model developed by Hays and Segall for adults was fitted to published data from infants and also to a retrospective data set collected at the Boston Children's Hospital (BCH). The BCH data set was also used to examine the relationship between uptake of FDG in different organs and patient weight or age. RESULTS: Substantial changes in the structure of the FDG model were required to fit the pediatric data. Fitted rate constants and fractional blood volumes were reduced relative to the adult values. CONCLUSIONS: The pharmacokinetic models developed differ substantially from adult pharmacokinetic (PK) models which can have considerable impact on the dosimetric models for pediatric patients. This approach may be used as a model for estimating dosimetry in children from other radiopharmaceuticals.

Automated tumour boundary delineation on (18)F-FDG PET images using active contour coupled with shifted-optimal thresholding method.

This study presents an automatic method to trace the boundary of the tumour in positron emission tomography (PET) images. It has been discovered that Otsu's threshold value is biased when the within-class variances between the object and the background are significantly different. To solve the problem, a double-stage threshold search that minimizes the energy between the first Otsu's threshold and the maximum intensity value is introduced. Such shifted-optimal thresholding is embedded into a region-based active contour so that both algorithms are performed consecutively. The efficiency of the method is validated using six sphere inserts (0.52-26.53 cc volume) of the IEC/2001 torso phantom. Both spheres and phantom were filled with (18)F solution with four source-to-background ratio (SBR) measurements of PET images. The results illustrate that the tumour volumes segmented by combined algorithm are of higher accuracy than the traditional active contour. The method had been clinically implemented in ten oesophageal cancer patients. The results are evaluated and compared with the manual tracing by an experienced radiation oncologist. The advantage of the algorithm is the reduced erroneous delineation that improves the precision and accuracy of PET tumour contouring. Moreover, the combined method is robust, independent of the SBR threshold-volume curves, and it does not require prior lesion size measurement.