• 中国科技论文统计源期刊
  • 中国科技核心期刊
  • 中国高校优秀期刊
  • 安徽省优秀科技期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

新型肾结石临床-影像组学模型在预测经皮肾镜结石清除率的应用价值

吴宇 付永 谷明利 刘建民

引用本文:
Citation:

新型肾结石临床-影像组学模型在预测经皮肾镜结石清除率的应用价值

    作者简介: 吴宇(1989-),男,硕士,主治医师
    通讯作者: 刘建民, liu_john_jm@sina.com
  • 基金项目:

    蚌埠医学院自然科学研究重点项目 BYKY2019164ZD

  • 中图分类号: R692.4

Application value of a new clinical-radiomics model of kidney stone in predicting the stone-free rate of percutaneous nephrolithotomy

    Corresponding author: LIU Jian-min, liu_john_jm@sina.com
  • CLC number: R692.4

  • 摘要: 目的构建预测经皮肾镜手术的结石清除率(SFR)的新型临床-影像组学模型并进行验证。方法回顾性收集113例行经皮肾镜取石术(PCNL)病人的相关资料, 根据术后复查的泌尿系CT或尿路平片, 将其分为结石清除组和结石残留。收集病人临床和影像组学资料, 利用图像分析软件及计算机程序设计语言工具, 将病人的CT图像划取感兴趣区域并提取出120个影像组学特征。对训练组进行变量选择, 得到最佳的特征选集, 采用多因素logistic回归分析构建新型临床-影像组学预测模型, 使用曲线下面积(AUC)评估模型的预测效果。结果113例病人术后1个月复查泌尿系影像, 结石清除者68例, 结石残留者45例, 总体SFR为60.2%。2组病人性别、术后血白细胞(WBC)、住院时间、结石长度、结石宽度及Guy's分级(GSS)差异均有统计学意义(P < 0.05~P < 0.01)。单因素logistic分析显示, 性别、GSS、术后血WBC、结石长度、结石宽度差异均具有统计学意义(P < 0.05~P < 0.01)。多因素logistic分析显示, 性别、GSS、术后血WBC为PCNL术后SFR的独立预测因子(P < 0.05~P < 0.01)。将Lasso回归筛选出的有意义的14个影像组学特征行单因素和多因素logistic分析, 结果显示, 最大三维直径、球度差异均具有统计学意义(P < 0.05和P < 0.01)。将最大三维直径、球度纳入构建临床-影像组学预测模型, 在训练集中曲线下面积为0.923, 在验证集中曲线下面积为0.876, 均优于GSS。结论构建的新型临床-影像组学模型结合临床特征, 有助于预测PCNL病人的SFR, 可为泌尿科医生术前沟通及手术方式的选择提供参考。
  • 图 1  VOI的划取(白圈为划取工具)

    图 2  Lasso回归筛选影像组学特征

    表 1  病人一般资料比较(x±s)

    特征 结石清除组(n=68) 结石残留组(n=45) t P
    年龄/岁 53.84±12.82 49.53±14.52 1.66 >0.05
    51 19 12.34 < 0.01
    体质量/kg 70.85±12.46 66.39±13.12 1.82 >0.05
    发热 4 5 1.01 >0.05
    结石长度/cm 2.30±1.12 3.87±2.46 4.60 < 0.01
    结石宽度/cm 1.54±0.72 2.48±1.41 4.66 < 0.01
    GSS
      Ⅰ 9 0 34.45 < 0.01
      Ⅱ 45 12
      Ⅲ 12 21
      Ⅳ 2 12
    肾穿刺
      双侧 19 13 5.16 >0.05
      左侧 31 12
      右侧 18 20
    术前WBC(×109/L) 6.09±1.35 6.81±2.62 1.92 >0.05
    术前中性粒细胞(×109/L) 0.15±0.50 0.20±0.55 0.50 >0.05
    术前胱抑素C/(mg/L) 1.04±0.31 1.11±0.43 1.01 >0.05
    术前尿酸/(μmol/L) 307.79±103.41 286.91±88.94 1.11 >0.05
    手术方式
      双侧PCNL 2 0 3.41 >0.05
      左侧PCNL 37 19
      右侧PCNL 29 26
    手术时间/min 137.27±44.91 150.02±42.95 1.50 >0.05
    术后WBC(×109/L) 8.04±2.99 9.43±3.73 2.19 < 0.05
    术后胱抑素C/(mg/L) 0.98±0.29 1.11±0.55 1.64 >0.05
      住院时间/d 12.00±3.45 16.40±9.06 3.63 < 0.01
    △示χ2
    下载: 导出CSV

    表 2  肾结石病人临床特征的单因素logistic回归分析(训练集)

    变量 B SE Waldχ2 P OR(95%CI)
    性别(对照=女) -1.43 0.41 12.36 < 0.01 0.24(0.11~0.54)
    GSS 1.78 0.36 24.79 < 0.01 5.93(3.11~12.63)
    术后血WBC 0.12 0.06 4.09 < 0.05 1.13(1.01~1.28)
    结石长度 0.51 0.14 13.30 < 0.01 1.66(1.30~2.24)
    结石宽度 0.83 0.22 14.34 < 0.01 2.30(1.55~3.67)
    术前血WBC 0.19 0.11 3.04 >0.05 1.21(0.99~1.52)
    下载: 导出CSV

    表 3  肾结石病人临床特征的多因素logistic回归分析(训练集)

    变量 B SE Waldχ2 P OR(95%CI)
    性别(对照=女) -1.96 0.59 10.97 < 0.01 0.14(0.04~0.41)
    GSS 1.88 0.54 12.34 < 0.01 6.56(2.48~20.24)
    术后血WBC 0.20 0.08 5.50 < 0.05 1.22(1.04~1.45)
    结石长度 -0.25 0.34 0.53 >0.05 0.78(0.41~1.56)
    结石宽度 0.3 0.56 0.29 >0.05 1.35(0.44~3.96)
    下载: 导出CSV

    表 4  病人影像组学特征的单因素回归分析(训练集)

    变量 B SE Waldχ2 P OR(95%CI)
    3D形状特征—伸长率 0.17 0.19 0.83 >0.05 1.19(0.82~1.73)
    3D形状特征—最大二维直径 1.85 0.34 29.37 < 0.01 6.36(3.46~13.19)
    3D形状特征—最大三维直径 0.66 0.22 9.42 < 0.01 1.94(1.30~3.03)
    3D形状特征—球度 -2.04 0.37 29.75 < 0.01 0.13(0.06~0.26)
    一阶统计量特征—最小值 -0.49 0.24 4.17 < 0.05 0.61(0.36~0.93)
    一阶统计量特征—健壮的平均绝对偏差 -0.04 0.20 0.04 >0.05 0.96(0.65~1.40)
    一阶统计量特征—总能量 1.19 0.39 9.33 < 0.01 3.29(1.74~8.02)
    灰度共生矩阵—集群阴影 0.25 0.19 1.58 >0.05 1.28(0.88~1.90)
    灰度共生矩阵—相关度 0.29 0.20 2.07 >0.05 1.34(0.91~2.02)
    灰度共生矩阵—信息相关度量2 -1.05 0.25 17.60 < 0.01 0.35(0.21~0.56)
    灰度共生矩阵—最大概率 -1.17 0.43 7.38 < 0.01 0.31(0.12~0.65)
    灰度形状矩阵—大区域高灰度强调 0.11 0.21 0.30 >0.05 1.12(0.75~1.69)
    灰度依赖矩阵—大依赖低灰度强调 -1.11 0.24 0.20 >0.05 0.90(0.53~1.34)
    邻域灰度差矩阵—粗糙度 -3.22 0.71 20.71 < 0.01 0.04(0.01~0.16)
    下载: 导出CSV

    表 5  病人影像组学特征的多因素回归分析(训练集)

    变量 B SE Waldχ2 P OR(95%CI)
    3D形状特征-最大二维直径 1.12 0.70 2.59 >0.05 3.07(0.86~13.20)
    3D形状特征-最大三维直径 -1.05 0.46 5.23 < 0.05 0.35(0.13~0.78)
    3D形状特征-球度 -2.04 0.66 9.53 < 0.01 0.13(0.03~0.40)
    一阶统计量特征-最小值 -0.80 0.45 3.19 >0.05 0.45(0.17~0.98)
    一阶统计量特征-总能量 0.40 0.61 0.42 >0.05 1.49(0.53~5.84)
    灰度共生矩阵-信息相关度量2 -0.65 0.52 1.61 >0.05 0.52(0.18~1.36)
    灰度共生矩阵-最大概率 -1.51 0.79 3.71 >0.05 0.22(0.04~0.87)
    邻域灰度差矩阵-粗糙度 0.43 0.91 0.22 >0.05 1.53(0.18~6.32)
    下载: 导出CSV
  • [1] TURK C, NEISIUS A, PETRIK A, et al. EAU Guideline on Urolithiasis. EAU Annual Congress Barcelona 2019: EAU Guidelines Office, Arnhem, The Netherlands (2019). Available online at: http://uroweb.org/guidelines/compilations-of-all-guidelines/ (accessed December 14, 2019).
    [2] 黄健. 中国泌尿外科和男科疾病诊断治疗指南2019版[M]. 北京: 科学出版社, 2020: 245.
    [3] OKHUNOV Z, FRIEDLANDER JI, GEORGE AK, et al. S.T.O.N.E. nephrolithometry: novel surgical classification system for kidney calculi[J]. Urology, 2013, 81(6): 1154. doi: 10.1016/j.urology.2012.10.083
    [4] THOMAS K, SMITH NC, HEGARTY N, et al. The Guy's stone score--grading the complexity of percutaneous nephrolithotomy procedures[J]. Urology, 2011, 78(2): 277. doi: 10.1016/j.urology.2010.12.026
    [5] SMITH A, AVERCH TD, SHAHROUR K, et al. A nephrolithometric nomogram to predict treatment success of percutaneous nephrolithotomy[J]. J Urol, 2013, 190(1): 149. doi: 10.1016/j.juro.2013.01.047
    [6] 陈路, 陈艾琪, 刘浩, 等. 磁共振影像组学在鉴别中低危和高危前列腺癌中的应用[J]. 蚌埠医学院学报, 2022, 47(1): 90.
    [7] LAMBIN P, LEIJENAAR R, DEIST TM, et al. Radiomics: the bridge between medical imaging and personalized medicine[J]. Nat Rev Clin Oncol, 2017, 14(12): 749. doi: 10.1038/nrclinonc.2017.141
    [8] LIMKIN EJ, SUN R, DERCLE L, et al. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology[J]. Ann Oncol, 2017, 28(6): 1191. doi: 10.1093/annonc/mdx034
    [9] 索杰, 段万里, 邓骞, 等. 肾下盏解剖对软式输尿管镜治疗肾下盏结石术后清石率的影响及预测模型构建[J]. 现代泌尿外科杂志, 2021, 26(11): 914.
    [10] CHOO MS, UHMN S, KIM JK, et al. A prediction model using machine learning algorithm for assessing stone-free status after single session shock wave lithotripsy to treat ureteral stones[J]. J Urol, 2018, 200(6): 1371. doi: 10.1016/j.juro.2018.06.077
    [11] VAN GRIETHUYSEN J, FEDOROV A, PARMAR C, et al. Computational radiomics system to decode the radiographic phenotype[J]. Cancer Res, 2017, 77(21): e104. doi: 10.1158/0008-5472.CAN-17-0339
    [12] 董建民, 刘建新. 基于PyRadiomics的医疗影像纹理获取原型系统集成[J]. 软件导刊, 2020, 19(7): 223.
    [13] 金鹏飞, 杨丽勤, 包婕, 等. 影像组学联合临床指标预测PI-RADS V2.13分病变中有临床意义前列腺癌[J]. 临床放射学杂志, 2022, 41(12): 2274.
    [14] HANLEY JA, MCNEIL BJ. The meaning and use of the area under a receiver operating characteristic (ROC) curve[J]. Radiology, 1982, 143(1): 29. doi: 10.1148/radiology.143.1.7063747
    [15] BOZZINI G, VERZE P, ARCANIOLO D, et al. A prospective randomized comparison among SWL, PCNL and RIRS for lower calyceal stones less than 2 cm: a multicenter experience: a better understanding on the treatment options for lower pole stones[J]. World J Urol, 2017, 35(12): 1967. doi: 10.1007/s00345-017-2084-7
    [16] 项立波, 何炜, 江敦勤, 等. 输尿管软镜取石术与经皮肾镜取石术治疗2~3 cm肾结石的有效性和安全性分析[J]. 蚌埠医学院学报, 2020, 45(7): 895.
    [17] KARAKOYUNLU N, GOKTUG G, ENER NC, et al. A comparison of standard PCNL and staged retrograde FURS in pelvis stones over 2 cm in diameter: a prospective randomized study[J]. Urolithiasis, 2015, 43(3): 283. doi: 10.1007/s00240-015-0768-2
    [18] WU WJ, OKEKE Z. Current clinical scoring systems of percutaneous nephrolithotomy outcomes[J]. Nat Rev Urol, 2017, 14(8): 459. doi: 10.1038/nrurol.2017.71
    [19] 崔振宇, 高彦君, 杨文增, 等. 可视标准通道联合F4.8可视穿刺超微经皮肾镜在肾多发结石中的应用(附46例报告)[J]. 中国内镜杂志, 2017, 23(9): 99.
    [20] 郭凡, 苏醒, 马帅军, 等. 微创经皮肾镜碎石取石术治疗肾结石影响一期清石的多因素分析[J]. 中国微创外科杂志, 2021, 21(6): 486.
    [21] RIVERA M, VIERS B, COCKERILL P, et al. Pre- and postoperative predictors of infection-related complications in patients undergoing percutaneous nephrolithotomy[J]. J Endourol, 2016, 30(9): 982.
    [22] YE Z, ZENG G, YANG H, et al. The status and characteristics of urinary stone composition in China[J]. BJU Int, 2020, 125(6): 801.
    [23] 杨杰, 王坤杰, 王佳. 经皮肾镜碎石取石术清石率影响因素分析[J]. 四川大学学报(医学版), 2019, 50(3): 445.
    [24] XUN Y, CHEN M, LIANG P, et al. A novel clinical-radiomics model pre-operatively predicted the stone-free rate of flexible ureteroscopy strategy in kidney stone patients[J]. Front Med (Lausanne), 2020, 7: 576925.
    [25] 刘欢, 邓远忠. Guy's结石评分和改良Clavien分级方法预测PCNL术后结石清除率和并发症的研究[J]. 医学信息, 2019, 32(9): 101.
    [26] 黄忠江, 姜增誉, 李健丁, 等. 基于增强CT影像组学联合机器学习鉴别均质性肾透明细胞癌与肾乏脂肪血管平滑肌脂肪瘤[J]. 实用医学杂志, 2021, 37(17): 2266.
  • [1] 李奔张华杨庆峰朱铁军杨涛杨辉刘坤 . 经皮肾镜取石术治疗无积水肾结石32例临床分析. 蚌埠医学院学报, 2014, 39(8): 1093-1094.
    [2] 徐久平邵明峰华和园余子强邹建安陈久发 . 微通道与标准通道经皮肾镜取石术治疗复杂性肾结石的比较. 蚌埠医学院学报, 2013, 37(8): 956-958,961.
    [3] 项立波何炜江敦勤陈弓王华方勇袁存和陆兆祥娄彦亭黄玉良刘军石栋 . 输尿管软镜取石术与经皮肾镜取石术治疗2~3 cm肾结石的有效性和安全性分析. 蚌埠医学院学报, 2020, 45(7): 895-897, 901. doi: 10.13898/j.cnki.issn.1000-2020.07.014
    [4] 梁建奇林灼怡 . 肥胖对经皮肾镜取石术临床疗效的影响. 蚌埠医学院学报, 2013, 37(7): 836-838.
    [5] 陈德钢刘成益徐鹏程许全超范欣欢张桃福李增稳郭亮 . 经肾上盏径路和经肾中下盏径路行经皮肾镜取石术的疗效比较. 蚌埠医学院学报, 2013, 37(9): 1099-1102.
    [6] 葛亮徐航庄惠强邓超雄 . 彩超引导下肾盏穹窿部穿刺在经皮肾结石取石术中的应用. 蚌埠医学院学报, 2016, 41(12): 1602-1604. doi: 10.13898/j.cnki.issn.1000-2200.2016.12.020
    [7] 郭一俊易绮尤佳周任远 . 经皮肾镜和输尿管软镜治疗肾结石的疗效比较. 蚌埠医学院学报, 2016, 41(9): 1181-1184. doi: 10.13898/j.cnki.issn.1000-2200.2016.09.020
    [8] 娄小华徐九云吴立新 . 联合评估策略下萧氏双C护理模式在经皮肾镜碎石术围手术期病人中的应用. 蚌埠医学院学报, 2023, 48(8): 1163-1167. doi: 10.13898/j.cnki.issn.1000-2200.2023.08.032
    [9] 戚敏俊曹志刚徐郑许露伟梁凯刘军苏江浩 . 经皮肾镜钬激光碎石术联合体外冲击波碎石术治疗复杂性肾结石46例. 蚌埠医学院学报, 2012, 36(11): 1310-1311,1314.
    [10] 顾恒王昌兵袁宇峰李博叶大文 . 超声引导经皮肾镜气压弹道联合超声碎石术治疗肾结石128例. 蚌埠医学院学报, 2013, 37(6): 724-726.
    [11] 陈芳芳 . 肾结石微创术后复发风险模型构建及护理干预研究. 蚌埠医学院学报, 2021, 46(8): 1133-1137. doi: 10.13898/j.cnki.issn.1000-2200.2021.08.037
    [12] 张建柱李庆文 . 复杂性肾结石46例的开放性手术治疗. 蚌埠医学院学报, 2012, 36(8): 964-966.
    [13] 戴廷山刘伟曾建华王二朋张长虹 . 肾窦内肾盂加肾后下肾实质切开治疗巨大鹿角形肾结石. 蚌埠医学院学报, 2007, 32(1): 40-41.
    [14] 张华吴斌李奔朱铁军杨庆峰杨涛 . B超引导下局麻经皮肾镜取石术治疗结石性脓肾. 蚌埠医学院学报, 2015, 40(9): 1189-1191. doi: 10.13898/j.cnki.issn.1000-2200.2015.09.015
    [15] 刘恒超陶新全申勇胡永全99mTc-DTPA肾动态显像评价上尿路结石患肾功能. 蚌埠医学院学报, 2009, 34(10): 888-890.
    [16] 刘定益俞家顺王健唐崎周燕峰李文敏胡桑 . L型无萎缩性肾切开结合输尿管镜治疗全鹿角形结石15例. 蚌埠医学院学报, 2019, 44(2): 225-227. doi: 10.13898/j.cnki.issn.1000-2200.2019.02.026
    [17] 谢海龙关超谷明利方文革许海斌 . 电子输尿管软镜联合钬激光治疗肾、输尿管上段结石的效果观察. 蚌埠医学院学报, 2016, 41(5): 619-622. doi: 10.13898/j.cnki.issn.1000-2200.2016.05.021
    [18] 王君梁宇郑金一李海波王茂富董志强李松 . ESWL和F-URS在不同体质量儿童肾结石病人中的应用. 蚌埠医学院学报, 2022, 47(7): 905-908. doi: 10.13898/j.cnki.issn.1000-2200.2022.07.016
    [19] 苏清华林浩群刘建梁善玲杨少伟王健 . 普通超声引导经皮肾镜取石术的效果分析. 蚌埠医学院学报, 2017, 42(1): 78-80. doi: 10.13898/j.cnki.issn.1000-2200.2017.01.022
    [20] 马艳陈冬梅常菲菲 . 快速康复外科理念在经皮肾镜取石术围手术期护理中的应用. 蚌埠医学院学报, 2018, 43(4): 547-549. doi: 10.13898/j.cnki.issn.1000-2200.2018.04.036
  • 加载中
图(2)表(5)
计量
  • 文章访问数:  299
  • HTML全文浏览量:  121
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-06
  • 录用日期:  2022-10-12
  • 刊出日期:  2023-08-15

新型肾结石临床-影像组学模型在预测经皮肾镜结石清除率的应用价值

    通讯作者: 刘建民, liu_john_jm@sina.com
    作者简介: 吴宇(1989-),男,硕士,主治医师
  • 1. 蚌埠医学院第二附属医院 泌尿外科, 安徽 蚌埠 233040
  • 2. 蚌埠医学院第一附属医院 泌尿外科, 安徽 蚌埠 233004
基金项目:  蚌埠医学院自然科学研究重点项目 BYKY2019164ZD

摘要: 目的构建预测经皮肾镜手术的结石清除率(SFR)的新型临床-影像组学模型并进行验证。方法回顾性收集113例行经皮肾镜取石术(PCNL)病人的相关资料, 根据术后复查的泌尿系CT或尿路平片, 将其分为结石清除组和结石残留。收集病人临床和影像组学资料, 利用图像分析软件及计算机程序设计语言工具, 将病人的CT图像划取感兴趣区域并提取出120个影像组学特征。对训练组进行变量选择, 得到最佳的特征选集, 采用多因素logistic回归分析构建新型临床-影像组学预测模型, 使用曲线下面积(AUC)评估模型的预测效果。结果113例病人术后1个月复查泌尿系影像, 结石清除者68例, 结石残留者45例, 总体SFR为60.2%。2组病人性别、术后血白细胞(WBC)、住院时间、结石长度、结石宽度及Guy's分级(GSS)差异均有统计学意义(P < 0.05~P < 0.01)。单因素logistic分析显示, 性别、GSS、术后血WBC、结石长度、结石宽度差异均具有统计学意义(P < 0.05~P < 0.01)。多因素logistic分析显示, 性别、GSS、术后血WBC为PCNL术后SFR的独立预测因子(P < 0.05~P < 0.01)。将Lasso回归筛选出的有意义的14个影像组学特征行单因素和多因素logistic分析, 结果显示, 最大三维直径、球度差异均具有统计学意义(P < 0.05和P < 0.01)。将最大三维直径、球度纳入构建临床-影像组学预测模型, 在训练集中曲线下面积为0.923, 在验证集中曲线下面积为0.876, 均优于GSS。结论构建的新型临床-影像组学模型结合临床特征, 有助于预测PCNL病人的SFR, 可为泌尿科医生术前沟通及手术方式的选择提供参考。

English Abstract

  • 肾结石作为泌尿外科常见疾病,部分国家其患病率可达10%以上[1],每年相关的治疗费用高达数十亿美元。经皮肾镜碎石取石术(percutaneous nephrolithotomy,PCNL)治疗复杂性结石、肾盂输尿管连接部大结石及鹿角形结石成为首选[2]。由于结石的负荷和分布、集合系统的解剖结构及术者的操作熟练度是因人而异的,所以术后结石残留在所难免。为满足临床医生需求,有国外学者提出S.T.O.N.E评分[3]、Guy′s分级(Guy′s score, GSS)[4]与CROES图表计数系统[5]用于术前评估肾结石的复杂程度及预测术后残石的发生率,有助于术前医生对于病人病情的评估及交流,然而这3种评分系统对结石清除率(stone-free rate,SFR)的预测精准度仍有争议。近年来,伴随影像技术的提升及影像组学的研究深入,其临床应用越发广泛。该影像组学的新方法已被证明能够影响和改变肿瘤领域的诊断和治疗策略[6-8]。此外,多项研究[9-10]表明,基于影像组学或机器学习的预测模型对于某些手术治疗,如输尿管软镜碎石和体外冲击波碎石的术后结果预测更加精准。基于此,本研究开发一种结合影像组学特征和临床指标的新型PCNL预测模型,特别是对于鹿角形肾结石,用于PCNL的术前评估和预测SFR。

    • 回顾性分析2018年1月至2019年6月我科行PCNL治疗的113例肾结石病人临床及影像学资料,所有病人均采用标准单通道,按照术后是否有结石残留分为结石清除组和结石残留组。排除标准:严重的心肺功能、凝血功能异常或脊柱畸形不能耐受手术者;重度尿路感染或肾功能不全者;合并同侧泌尿系统肿瘤的病人。病人术后一个月复查尿路平片,确定病人肾结石清除的评估标准为:未发现肾结石表现或残石碎片 < 4 mm。收集所有病人临床特征和术前影像学资料,依据病人CT图像获取相应的影像组学数据。采用随机种子法将病人数据以7∶ 3比例随机分为训练集79例和验证集34例。GSS依据结石部位、情况以及肾脏解剖情况分为Ⅰ~Ⅳ级。Ⅰ级为肾盂或中/下盏单发结石,具有简单的集合系统解剖结构;Ⅱ级为集合系统解剖正常且上盏单发结石、集合系统解剖正常且多发结石、集合系统解剖异常且单发结石;Ⅲ级为集合系统解剖异常且多发结石、部分鹿角形结石、肾盏憩室结石;Ⅳ级为完全鹿角形结石或伴有任何的脊柱损伤或畸形。

    • (1) 一般情况,包括年龄、性别、术前1周内是否发热;(2)结石位置,包括上盏、中盏、下盏、肾盂、输尿管;(3)结石大小,根据本院CT报告单描述的结石长度与宽度;(4)术前尿常规及尿培养;(5)术前、术后血常规及生化常规;(6)病人住院时间。

    • 收集病人泌尿系CT影像资料,CT影像扫描参数:管电压100~120 kV,自动管电流200~350 mA,旋转时间0.5 s,层厚5 mm。按以下流程分析:(1)通过开源软件3D-slicer (version 4.9.0;www.slicer.org)导入DICOM格式CT影像数据,构建Radiomics数据库。(2)手动逐层划取结石感兴趣区域(volume of interest, VOI), 在勾画结石轮廓时,设置阈值区分肾脏及周围组织,提高VOI的准确性。影像组学特征数据源于VOI,为降低人为误差,每位病人的结石轮廓线由2名泌尿外科住院医师共同勾画,汇总后由1名泌尿外科主任医师确定病人的全部结石被包含在VOI中(见图 1),将生成的VOI保存后得到包含病人影像组学数据特征的nrrd文件及其标签文件。(3)打开计算机程序设计语言工具Python(version3.8),利用PyRadiomics语言分析包提取并计算所有特征值。PyRadiomics[11]是一个开源的分析平台,包含目前所有的影像纹理分析算法,利用相对应的Dicom格式开源模块,对数据文件进行解码和数据获取,通过定义及批量计算将数据处理,最终获得量化后的影像组学特征[12]。本研究共提取出120个影像组学特征,包括以下几大类:(1)一阶统计量特征19个;(2)2D形状特征10个;(3)3D形状特征16个;(4)灰度共生矩阵24个;(5)灰度长度矩阵16个;(6)灰度形状矩阵16个;(7)邻域灰度差矩阵5个;(8)灰度依赖矩阵14个。

      图  1  VOI的划取(白圈为划取工具)

    • Lasso回归是一种用于多共线数据或高维数据的高级变量选择算法。既往研究[13]表明,将Lasso回归引入到多共线数据或高维数据分析中,在简化模型复杂性的同时,提高模型的预测精准度。使用10倍交叉验证法选择λ值,以误差最小的λ值为标准,所对应的因素即为筛选出来影响清石率的因素。

    • 在训练集中,使用logistic回归方程构建预测模型,该模型结合了用于预测无结石状态的影像组学特征和临床因素。通过受试者工作特征(ROC)曲线计算的曲线下面积(AUC)用于评估已建立模型的辨别性能[14]。描绘校准曲线以评估临床-影像组学模型的预测准确性,然后进行Hosmer-Lemeshow拟合优度检验。采用十折交叉验证,计算AUC值并评估模型的预测效果, 将构建好的最佳预测模型应用在验证集中,来进一步验证模型的预测效果。

    • 麻醉取截石体位,常规消毒、铺巾。输尿管镜下将输尿管导管插入患侧,将导尿管与输尿管导管相固定,建立人工肾积水。更换俯卧位,垫高腰部,重新消毒、铺巾。在B超引导下,根据结石大小位置和肾积水程度选择穿刺点,穿刺成功后置入金属导丝,两步法将合适大小的筋膜扩张器扩张至目标盏内,保留外鞘。采用瑞士EMS超声气压弹道碎石机碎石后吸出或冲出,探查视野无明显结石残留及出血,留置输尿管支架管后退镜。留置肾造瘘管一根并清洁包扎固定造瘘口。

    • 采用t检验、秩和检验、χ2检验、Lasso回归分析、logistic回归分析和ROC曲线分析。

    • 113例病人术后一个月复查泌尿系影像,结石清除者68例,结石残留者45例,总体SFR为60.2%。2组病人性别、术后血白细胞(WBC)、住院时间、结石长度、结石宽度及GSS差异均有统计学意义(P < 0.05~P < 0.01)(见表 1)。

      特征 结石清除组(n=68) 结石残留组(n=45) t P
      年龄/岁 53.84±12.82 49.53±14.52 1.66 >0.05
      51 19 12.34 < 0.01
      体质量/kg 70.85±12.46 66.39±13.12 1.82 >0.05
      发热 4 5 1.01 >0.05
      结石长度/cm 2.30±1.12 3.87±2.46 4.60 < 0.01
      结石宽度/cm 1.54±0.72 2.48±1.41 4.66 < 0.01
      GSS
        Ⅰ 9 0 34.45 < 0.01
        Ⅱ 45 12
        Ⅲ 12 21
        Ⅳ 2 12
      肾穿刺
        双侧 19 13 5.16 >0.05
        左侧 31 12
        右侧 18 20
      术前WBC(×109/L) 6.09±1.35 6.81±2.62 1.92 >0.05
      术前中性粒细胞(×109/L) 0.15±0.50 0.20±0.55 0.50 >0.05
      术前胱抑素C/(mg/L) 1.04±0.31 1.11±0.43 1.01 >0.05
      术前尿酸/(μmol/L) 307.79±103.41 286.91±88.94 1.11 >0.05
      手术方式
        双侧PCNL 2 0 3.41 >0.05
        左侧PCNL 37 19
        右侧PCNL 29 26
      手术时间/min 137.27±44.91 150.02±42.95 1.50 >0.05
      术后WBC(×109/L) 8.04±2.99 9.43±3.73 2.19 < 0.05
      术后胱抑素C/(mg/L) 0.98±0.29 1.11±0.55 1.64 >0.05
        住院时间/d 12.00±3.45 16.40±9.06 3.63 < 0.01
      △示χ2

      表 1  病人一般资料比较(x±s)

    • 表 1中差异有统计学意义的变量进一步行单因素和多因素logistic分析。单因素logistic分析显示,性别、GSS、术后血WBC、结石长度、结石宽度差异均具有统计学意义(P < 0.05~P < 0.01)(见表 2)。多因素logistic分析显示,性别、GSS、术后血WBC为PCNL术后SFR的独立预测因子(P < 0.05~P < 0.01)(见表 3)。

      变量 B SE Waldχ2 P OR(95%CI)
      性别(对照=女) -1.43 0.41 12.36 < 0.01 0.24(0.11~0.54)
      GSS 1.78 0.36 24.79 < 0.01 5.93(3.11~12.63)
      术后血WBC 0.12 0.06 4.09 < 0.05 1.13(1.01~1.28)
      结石长度 0.51 0.14 13.30 < 0.01 1.66(1.30~2.24)
      结石宽度 0.83 0.22 14.34 < 0.01 2.30(1.55~3.67)
      术前血WBC 0.19 0.11 3.04 >0.05 1.21(0.99~1.52)

      表 2  肾结石病人临床特征的单因素logistic回归分析(训练集)

      变量 B SE Waldχ2 P OR(95%CI)
      性别(对照=女) -1.96 0.59 10.97 < 0.01 0.14(0.04~0.41)
      GSS 1.88 0.54 12.34 < 0.01 6.56(2.48~20.24)
      术后血WBC 0.20 0.08 5.50 < 0.05 1.22(1.04~1.45)
      结石长度 -0.25 0.34 0.53 >0.05 0.78(0.41~1.56)
      结石宽度 0.3 0.56 0.29 >0.05 1.35(0.44~3.96)

      表 3  肾结石病人临床特征的多因素logistic回归分析(训练集)

    • 在训练集中,通过Lasso回归分析提取的120个影像组学特征,采用最小标准和10倍交叉验证法,通过引入惩罚调整参数(λ)对模型系数进行压缩,使无关系数趋向于零,从而实现变量的自动筛选,筛选出14个影像组学特征有意义(见图 2)。

      图  2  Lasso回归筛选影像组学特征

    • 将Lasso回归筛选出的有意义的变量进一步行单因素和多因素logistic分析,结果显示,最大三维直径、球度差异均具有统计学意义(P < 0.05和P < 0.01)(见表 45)。

      变量 B SE Waldχ2 P OR(95%CI)
      3D形状特征—伸长率 0.17 0.19 0.83 >0.05 1.19(0.82~1.73)
      3D形状特征—最大二维直径 1.85 0.34 29.37 < 0.01 6.36(3.46~13.19)
      3D形状特征—最大三维直径 0.66 0.22 9.42 < 0.01 1.94(1.30~3.03)
      3D形状特征—球度 -2.04 0.37 29.75 < 0.01 0.13(0.06~0.26)
      一阶统计量特征—最小值 -0.49 0.24 4.17 < 0.05 0.61(0.36~0.93)
      一阶统计量特征—健壮的平均绝对偏差 -0.04 0.20 0.04 >0.05 0.96(0.65~1.40)
      一阶统计量特征—总能量 1.19 0.39 9.33 < 0.01 3.29(1.74~8.02)
      灰度共生矩阵—集群阴影 0.25 0.19 1.58 >0.05 1.28(0.88~1.90)
      灰度共生矩阵—相关度 0.29 0.20 2.07 >0.05 1.34(0.91~2.02)
      灰度共生矩阵—信息相关度量2 -1.05 0.25 17.60 < 0.01 0.35(0.21~0.56)
      灰度共生矩阵—最大概率 -1.17 0.43 7.38 < 0.01 0.31(0.12~0.65)
      灰度形状矩阵—大区域高灰度强调 0.11 0.21 0.30 >0.05 1.12(0.75~1.69)
      灰度依赖矩阵—大依赖低灰度强调 -1.11 0.24 0.20 >0.05 0.90(0.53~1.34)
      邻域灰度差矩阵—粗糙度 -3.22 0.71 20.71 < 0.01 0.04(0.01~0.16)

      表 4  病人影像组学特征的单因素回归分析(训练集)

      变量 B SE Waldχ2 P OR(95%CI)
      3D形状特征-最大二维直径 1.12 0.70 2.59 >0.05 3.07(0.86~13.20)
      3D形状特征-最大三维直径 -1.05 0.46 5.23 < 0.05 0.35(0.13~0.78)
      3D形状特征-球度 -2.04 0.66 9.53 < 0.01 0.13(0.03~0.40)
      一阶统计量特征-最小值 -0.80 0.45 3.19 >0.05 0.45(0.17~0.98)
      一阶统计量特征-总能量 0.40 0.61 0.42 >0.05 1.49(0.53~5.84)
      灰度共生矩阵-信息相关度量2 -0.65 0.52 1.61 >0.05 0.52(0.18~1.36)
      灰度共生矩阵-最大概率 -1.51 0.79 3.71 >0.05 0.22(0.04~0.87)
      邻域灰度差矩阵-粗糙度 0.43 0.91 0.22 >0.05 1.53(0.18~6.32)

      表 5  病人影像组学特征的多因素回归分析(训练集)

    • 将多因素logistic回归分析有意义的变量依次纳入回归方程,得出预测PCNL术后SFR的概率模型。

      临床概率模型:P=1/1+EXP(-4.896+1.78×GSS)

      影像组学概率模型:P=1/1+EXP(-3.392-0.4 156×最大三维直径-1.999 2×球度+1.03×GSS)

    • 纳入GSS的模型训练集AUC值为0.843 (0.83~0.95), 诊断阈值为0.46,灵敏度为76.3%,特异度为83.5%;验证集AUC值为0.697(0.64~0.99), 诊断阈值为0.45,灵敏度为67.2%,特异度为73.3%;纳入影像组学两个指标(最大三维直径,球度)综合预测效果达0.923(0.85~0.97),诊断阈值为0.87,灵敏度为99.5%,特异度为91.3%;在验证队列中验证建立的模型,AUC值为0.876(0.61~0.98),诊断阈值为0.45,灵敏度为79.7%,特异度为83.2%;优于GSS。

    • 肾结石发病率居泌尿系结石首位。对于>2 cm的肾结石,PCNL的SFR相对ESWL及输尿管软镜更有优势,并且其手术费用更低[15-16]。但PCNL病人术后结石残留及并发症的发生不可避免,导致病人手术无效,需要再次治疗[17]。这带来巨大的安全风险和医疗资源浪费,对医患双方均造成很大困扰,所以术前评估PCNL手术效果,选择合适的病人尤为重要。伴随多种结石评分系统的相继提出,较常用的有GSS、CROES评分、S.T.O.N.E评分以及S-ReSC分级评分系统[18]。而一个评分指标易获取、评分标准易学习、评价结果稳定且一致性良好的结石评分系统具有最佳的临床应用前景。本研究选择GSS,考虑其影像学数据获取及分析简单且稳定,在应用过程中完善相关定义,如完全与不完全鹿角型结石划分、肾脏解剖异常的概念,可减少误差、避免分歧,提高SFR预测的准确性。

      本研究旨在开发一种新型临床-影像组学预测模型,该模型是临床因素和影像组学特征的结合体。通过收集导致结石残留的变量,预测PCNL术后的SFR,从而针对不同肾结石病人提供个体化的治疗。本研究结果显示,结石残留组GSS为Ⅲ级46.7%、Ⅳ级26.7%的构成比明显高于结石清除组,表明GSS级别越高术后越容易出现结石残留。GSS常被用来评估肾脏解剖结构、结石负荷和分布等因素对PCNL术后结石残留的评分指标。对于多发性或鹿角形结石, 因结石负荷大、分布散或角度刁钻, 如平行盏、盏颈狭长等, 从而导致取石困难、手术时间延长及术后结石残留多等特点[19-20]

      本研究中,将临床特征有意义的变量进行单因素及多因素logistic分析,结果显示,性别、GSS、术后血WBC升高是PCNL术后SFR的独立危险因素。一般多发性结石、鹿角型结石或感染性结石的病人手术时间较长,碎石过程中因肾盂黏膜损坏、血管破裂及感染性物质吸收增加,导致术后血WBC升高[21]。YE等[22]发现草酸钙结石和尿酸结石多见于男性,碳酸磷结石和磷酸镁铵结石在女性病人中更常见。而草酸钙及尿酸结石相对更坚硬, 导致击碎难度增加, 若碎石不彻底,则直接引起术后残石率增高。而术后感染或结石残留的病人,其住院天数也会相应的增加。

      有文献[23]报道,随着结石体积及表面积的增加,结石负荷相应变大导致清石率下降。本研究将Lasso回归筛选出有意义的影像组学变量进行单因素和多因素logistic分析,结果显示,最大三维直径和球度是PCNL术后SFR的独立危险因素,均纳入预测模型中。

      将影像组学特征数值标准化处理后,得出两个模型的公式系数,将GSS和影像组学两个指标(最大三维直径、球度)依次纳入。纳入影像组学两个指标的综合预测效果达0.923,验证队列中验证建立模型的AUC值为0.876,均优于GSS,与XUN等[24-25]研究相接近, 也间接反映该模型的准确性。

      随着精准医疗的快速发展,影像组学的应用越发广泛[26]。计算机辅助成像技术的快速发展使从数字医学图像中高通量提取定量特征成为可能,可以直观地理解为将视觉影像信息转化为深层次的特征来进行量化研究。随着影像学技术和设备的不断提升,针对肾结石PCNL清石率的研究也越来越多。尽管影像组学在评估结石方面具有一定的优势,但目前仍无法广泛应用。本研究既要追求模型的准确性,也要重视模型的实用性,故保留部分临床因素的分析。PCNL术后结局受到术前、术中、术后等多方因素的影响。故本研究存在以下不足:首先,这是一项回顾性研究,可能会导致选择偏倚;其次,样本量相对不足及在勾画VOI时会存在一定误差;此外,该模型是基于单中心数据建立的,需要前瞻性的多中心研究来进一步验证我们的结果;最后,由于实际限制,本研究无法对多发性结石和下盏结石进行全面回顾,这可能是未来前瞻性研究的目标。

      综上所述,本研究开发的新型临床-影像组学模型对于经皮肾镜病人的术前SFR具有良好的预测能力,可协助泌尿科医生与病人沟通,提供术前咨询并预测手术结果,最终可以选择合适的手术治疗方法来清除肾结石。

参考文献 (26)

目录

    /

    返回文章
    返回