Institute for Artificial Intelligence in Medicine (IKIM), University Hospital Essen, Essen, Germany.
Department of Diagnostic and Interventional Radiology and Neuroradiology, University Hospital Essen, Essen, Germany.
Eur J Nucl Med Mol Imaging. 2023 Jun;50(7):2196-2209. doi: 10.1007/s00259-023-06163-x. Epub 2023 Mar 2.
The aim of this study was to systematically evaluate the effect of thresholding algorithms used in computer vision for the quantification of prostate-specific membrane antigen positron emission tomography (PET) derived tumor volume (PSMA-TV) in patients with advanced prostate cancer. The results were validated with respect to the prognostication of overall survival in patients with advanced-stage prostate cancer.
A total of 78 patients who underwent [Lu]Lu-PSMA-617 radionuclide therapy from January 2018 to December 2020 were retrospectively included in this study. [Ga]Ga-PSMA-11 PET images, acquired prior to radionuclide therapy, were used for the analysis of thresholding algorithms. All PET images were first analyzed semi-automatically using a pre-evaluated, proprietary software solution as the baseline method. Subsequently, five histogram-based thresholding methods and two local adaptive thresholding methods that are well established in computer vision were applied to quantify molecular tumor volume. The resulting whole-body molecular tumor volumes were validated with respect to the prognostication of overall patient survival as well as their statistical correlation to the baseline methods and their performance on standardized phantom scans.
The whole-body PSMA-TVs, quantified using different thresholding methods, demonstrate a high positive correlation with the baseline methods. We observed the highest correlation with generalized histogram thresholding (GHT) (Pearson r (r), p value (p): r = 0.977, p < 0.001) and Sauvola thresholding (r = 0.974, p < 0.001) and the lowest correlation with Multiotsu (r = 0.877, p < 0.001) and Yen thresholding methods (r = 0.878, p < 0.001). The median survival time of all patients was 9.87 months (95% CI [9.3 to 10.13]). Stratification by median whole-body PSMA-TV resulted in a median survival time from 11.8 to 13.5 months for the patient group with lower tumor burden and 6.5 to 6.6 months for the patient group with higher tumor burden. The patient group with lower tumor burden had significantly higher probability of survival (p < 0.00625) in eight out of nine thresholding methods (Fig. 2); those methods were SUVmax50 (p = 0.0038), SUV ≥3 (p = 0.0034), Multiotsu (p = 0.0015), Yen (p = 0.0015), Niblack (p = 0.001), Sauvola (p = 0.0001), Otsu (p = 0.0053), and Li thresholding (p = 0.0053).
Thresholding methods commonly used in computer vision are promising tools for the semiautomatic quantification of whole-body PSMA-TV in [Ga]Ga-PSMA-11-PET. The proposed algorithm-driven thresholding strategy is less arbitrary and less prone to biases than thresholding with predefined values, potentially improving the application of whole-body PSMA-TV as an imaging biomarker.
本研究旨在系统评估计算机视觉中用于量化前列腺特异性膜抗原正电子发射断层扫描(PSMA-PET)衍生肿瘤体积(PSMA-TV)的阈值算法的效果,该效果与晚期前列腺癌患者的总生存预后相关。
回顾性纳入 2018 年 1 月至 2020 年 12 月期间接受 [Lu]Lu-PSMA-617 放射性核素治疗的 78 例患者。在放射性核素治疗前获取 [Ga]Ga-PSMA-11 PET 图像,用于分析阈值算法。所有 PET 图像均首先使用预评估的专有软件解决方案进行半自动分析,作为基线方法。随后,应用五种基于直方图的阈值方法和两种在计算机视觉中广泛应用的局部自适应阈值方法,以量化分子肿瘤体积。使用不同的阈值方法定量得到的全身 PSMA-TVs 与基线方法具有高度正相关性。我们观察到与广义直方图阈值法(GHT)的相关性最高(Pearson r(r),p 值(p):r = 0.977,p < 0.001)和 Sauvola 阈值法(r = 0.974,p < 0.001),与 Multiotsu 法(r = 0.877,p < 0.001)和 Yen 阈值法(r = 0.878,p < 0.001)的相关性最低。所有患者的中位总生存期为 9.87 个月(95%CI [9.3 至 10.13])。按中位全身 PSMA-TV 分层,肿瘤负荷较低的患者组中位生存时间为 11.8 至 13.5 个月,肿瘤负荷较高的患者组中位生存时间为 6.5 至 6.6 个月。在八种阈值方法中的八种方法中,肿瘤负荷较低的患者组具有更高的生存概率(p < 0.00625);这些方法是 SUVmax50(p = 0.0038)、SUV≥3(p = 0.0034)、Multiotsu(p = 0.0015)、Yen(p = 0.0015)、Niblack(p = 0.001)、Sauvola(p = 0.0001)、Otsu(p = 0.0053)和 Li 阈值法(p = 0.0053)(图 2)。
计算机视觉中常用的阈值方法是半自动化量化 [Ga]Ga-PSMA-11-PET 全身 PSMA-TV 的有前途的工具。与使用预定义值的阈值法相比,该算法驱动的阈值策略不那么任意,也不易产生偏差,可能会提高全身 PSMA-TV 作为成像生物标志物的应用。