基于冠状动脉CTA图像粥样硬化斑块综合特征的患者危险度分级系统

技术领域

本发明涉及冠状动脉粥样硬化性心脏病患者危险度分级及辅助临床个体化治疗领域，具体涉及一种根据提取自冠状动脉CTA图像上粥样硬化斑块的特征进行综合权重，并利用权重积分结果对患者进行精准危险度分级的系统。

背景技术

无创的冠状动脉CTA检查能提供多样的冠状动脉解剖学测量信息，对狭窄性冠状动脉疾病(coronary artery disease,CAD)的检出率及排除率均高于侵入性的冠状动脉造影。已有研究证实冠状动脉CTA的检查结果与患者恶性心血管事件的发生率相关，CTA图像上存在斑块的CAD患者恶性心血管事件的发生率显著高于CTA图像上不存在斑块的CAD患者，CTA图像上无阳性表现的患者恶性心血管事件的发生率最低。作为临床广泛应用的无创检查手段，最大限度的对CTA图像进行解读有利于精准医疗的实施，从而对有限的医疗资源合理再分配。

目前，CTA图像的解读主要使用标准化冠状动脉疾病报告与数据系统(CoronaryArtery Disease Reporting and Data System，CAD-RADS)，依据冠状动脉最大狭窄程度将患者由低到高分为0-5级共六个危险度等级，该系统在临床已得到普及使用。然而，CAD-RADS并未综合的体现出CTA图像上观测到的斑块特征。研究证实，与临床实际的随访结果相对照，单纯依靠CAD-RADS进行危险度分级存在过高或过低的预判，临床上迫切的需要更准确的个体化分级方法来满足精准医疗的需求。

大型队列研究证实，CTA图像上的粥样硬化斑块信息，包括斑块位置、斑块数目、以及斑块成分(非钙化斑、混合斑或钙化斑)等都与CAD患者的临床结局密切相关，但目前尚未见到将这些复杂信息融合成危险度分级体系的报道。

发明内容

本发明的目的在于提供一种基于冠状动脉CTA图像粥样硬化斑块综合特征的患者危险度分级系统，该系统根据提取自冠状动脉CTA图像上粥样硬化斑块的综合特征进行权重积分，从而对CAD患者进行更精准的危险度分级，有助于辅助临床个体化诊治。

为达到上述目的，本发明采用了以下技术方案：

一种基于冠状动脉CTA图像的患者危险度分级系统，该危险度分级系统包括斑块特征积分模块及患者危险度分级模块；所述危险度分级系统通过将冠状动脉CTA图像上粥样硬化斑块依据综合特征(通常包括多个斑块的特征信息)进行权重积分，从而对相应CAD患者的危险度进行分级，以更针对性的指导治疗及调整随访间隔时间，为CAD患者个体化诊治(例如，防治恶性心血管事件)提供依据；

所述斑块特征积分模块用于对CAD患者冠状动脉CTA图像中一个或多个(例如，冠状动脉CTA图像上可见的所有单个粥样硬化斑块)粥样硬化斑块按照斑块特征的综合权重系数进行权重积分；所述斑块特征包括斑块位置、斑块成分、管腔狭窄程度及斑块高危征象，所述综合权重系数是斑块位置、斑块成分、管腔狭窄程度及斑块高危征象的权重系数的乘积；

所述患者危险度分级模块用于依据上述CAD患者冠状动脉CTA图像的相应粥样硬化斑块的权重积分结果(CTA Based CAD Score，即总积分)，确定该患者的风险等级，所述风险等级包括低风险组(总积分＜5)、中风险组(总积分为5-10)及高风险组(总积分＞10)。

优选的，所述斑块位置的权重系数是按照Leaman score对粥样硬化斑块所处冠状动脉不同分段位置进行权重而确定的，该权重系数的取值范围为0～6。

优选的，所述斑块成分的权重系数是通过对粥样硬化斑块成分CT值测量结果的定性分类(例如，非钙化斑、钙化斑、混合斑)进行权重而确定的，该权重系数的取值范围为1～2。

优选的，所述管腔狭窄程度的权重系数是通过对粥样硬化斑块阻塞冠状动脉管腔的程度进行权重而确定的，该权重系数的取值范围为1～2。

优选的，所述斑块高危征象的权重系数是通过对无斑块高危征象或出现的特定斑块高危征象(例如，低密度斑、点状钙化斑、“餐巾环”征)进行权重而确定的，该权重系数的取值范围为1～2。

优选的，所述斑块特征积分模块采用以下公式进行权重积分：

其中，l表示斑块位置的权重系数，c表示斑块成分的权重系数，s表示管腔狭窄程度的权重系数，h表示斑块高危征象的权重系数，(l

优选的，所述危险度分级系统还包括斑块特征识别模块，该模块用于分析或输入冠状动脉CTA图像粥样硬化斑块的斑块位置、斑块成分、管腔狭窄程度及斑块高危征象的权重系数。

一种冠状动脉CTA图像的斑块特征融合方法，该特征融合方法包括以下步骤：

提取冠状动脉CTA图像上粥样硬化斑块的特征，然后计算斑块特征的综合权重系数，所述斑块特征包括斑块位置、斑块成分、管腔狭窄程度及斑块高危征象，所述综合权重系数是斑块位置、斑块成分、管腔狭窄程度及斑块高危征象的权重系数的乘积。

优选的，所述特征融合方法还包括以下步骤：对一个或多个粥样硬化斑块按照斑块特征的综合权重系数进行权重积分。

本发明的有益效果体现在：

本发明提出的危险度分级系统依据CTA图像上粥样硬化斑块的综合特征进行权重积分，利用综合权重系数实现了对复杂斑块特征的有效融合，充分利用了CTA图像中的斑块特征信息，提高了CAD患者危险度分级的准确度。与CAD-RADS相比，本发明得到的危险度分级可获得更好的患者区分度及最佳重分类改善指数(Net Reclassification Index，NRI)，提高了CTA结果的临床应用价值。本发明的分级流程快捷、简便，有利于临床广泛普及使用。

附图说明

图1为斑块位置示意图；其中：1～17为冠状动脉分段编号。

图2为斑块成分示意图。

图3为斑块狭窄程度示意图；其中：a为狭窄段的最窄直径，b为狭窄远端的正常直径。

图4为斑块高危征象示意图。

图5为使用本发明的危险度分级系统得到的患者无事件生存率(A)及使用CAD-RADS分级得到的患者无事件生存率(B)。

具体实施方式

下面结合附图和实施例对本发明做进一步详细说明，所述实施例仅用于解释本发明，而非对本发明保护范围的限制。

本发明基于研究结果和文献回顾，将CTA源的多种粥样硬化斑块特征(斑块位置、斑块成分、管腔狭窄程度及斑块高危征象)通过权重积分进行融合，并依据所得总积分对CAD患者进行更加精准、精细化的危险度分级，可以为临床CAD患者个体化诊治提供依据。

(一)冠状动脉CTA图像上粥样硬化斑块信息的综合权重积分

斑块位置权重系数的确定：依据美国心脏病协会17段分段法对斑块位置进行描述(图1)，同时使用Leaman score对斑块位置进行权重，详见表1。

表1.斑块位置权重系数

斑块成分权重系数的确定：依据测得的斑块成分CT值进行定性(图2)，依据定性结果将斑块分类并进行权重；其中，CT值＞130HU定义为钙化斑，权重系数为1.1；CT值≤130HU定义为非钙化斑，权重系数为1.2；存在钙化斑及非钙化斑两种成分定义为混合斑，权重系数为1.3，详见表2。

表2.斑块成分权重系数

管腔狭窄程度权重系数的确定：依据斑块在长轴图像上对管腔的阻塞程度进行判断(图3)，并按照阻塞程度不同进行权重；其中，轻度狭窄定义为狭窄程度1％-49％，权重系数为1.0；中度狭窄定义为狭窄程度50％-69％，权重系数为1.4；重度狭窄定义为狭窄程度70％-100％或者左冠状动脉主干狭窄程度≥50％，权重系数为1.5，详见表3。

表3.管腔狭窄程度权重系数

高危斑块征象权重系数的确定：依据高危斑块征象出现与否以及出现特定的高危斑块征象进行权重(图4)；其中，无高危斑块特征出现，权重系数为1.0；低密度斑块表现为非钙化斑块且斑块成分CT值≤30HU，权重系数为1.4；点状钙化斑表现为钙化斑最大径(长径)≤3.0mm，权重系数为1.5；“餐巾环”征表现为环状纤维帽(CT值＜130HU)环绕低密度脂质核心区(CT值≤30HU)，权重系数为1.6，详见表4。

表4.斑块高危征象权重系数

患者个体总积分(CTA Based CAD Score)的确定：

将CAD患者冠状动脉CTA图像上各粥样硬化斑块的斑块位置权重系数、斑块成分权重系数、管腔狭窄程度权重系数及高危斑块征象权重系数输入计算机，由计算机按照以下公式完成综合权重系数的权重积分，得到患者个体总积分：

其中，i代表CAD患者冠状动脉CTA图像上第i个单个粥样硬化斑块，l代表斑块位置的权重系数，c代表斑块成分的权重系数，s代表管腔狭窄程度的权重系数，h代表斑块高危征象的权重系数。

计算机依据患者个体总积分，将患者分为低风险组(＜5分)、中风险组(5-10分)或高风险组(＞10分)，各组患者5年内恶性心血管事件发病率按风险等级依次升高。

(二)CAD患者冠状动脉CTA图像积分实例

实例1

通过分析某患者CTA图像，发现该患者为右冠状动脉优势型，共有4处粥样硬化斑块。第一处斑块位于右冠状动脉近段，测得该斑块CT值为110HU，管腔狭窄程度为30％，无高危征象；第二处斑块位于前降支动脉近段，测得该斑块CT值为110HU，管腔狭窄程度为70％，呈“餐巾环”征高危征象；第三处斑块位于第一对角支，测得斑块CT值为210HU，管腔狭窄程度为50％，呈“点状钙化斑”高危征象；第四处斑块位于回旋支动脉近段，测得该斑块存在CT值为90HU及210HU两种成分，管腔狭窄程度为10％，无高危征象。

上述患者第一处斑块的位置权重系数为1.0，成分权重系数为1.2，狭窄程度权重系数为1.0，高危斑块权重系数为1.0，第一处斑块综合权重系数计1.2；第二处斑块的位置权重系数3.5，成分权重系数为1.2，狭窄程度权重系数为1.5，高危斑块权重系数为1.6，第二处斑块综合权重系数计10.08；第三处斑块的位置权重系数为1.0，成分权重系数为1.1分，狭窄程度权重系数为1.4，高危斑块权重系数为1.5，第三处斑块综合权重系数计2.31分；第四处斑块的位置权重系数为1.5，成分权重系数为1.3，狭窄程度权重系数为1.0，高危斑块权重系数为1.0，第四处斑块综合权重系数计1.95。通过权重积分得到该患者个体总积分为15.54分，为一例高风险组患者，提示未来恶性心血管事件发生率较高，需要积极的治疗和密切随访。

实例2

CTA图像显示，某患者为右冠状动脉优势型，共二处粥样硬化斑块，第一处斑块为位于右冠状动脉远段的钙化斑，管腔狭窄程度为70％，无高危征象，第二处斑块为位于第一对角支的非钙化斑，管腔狭窄程度为30％，无高危征象。

依据CAD-RADS分级，该患者最大狭窄程度为70％，分级为4级高危级。依据本发明的分级系统，该患者总积分＝1.0×1.1×1.5×1.0+1.0×1.2×1.0×1.0＝2.85，分级为低风险组。而患者五年内随访并无恶性心血管事件发生，可见本发明的危险度分级更准确。

实例3

为验证本发明的危险度分级系统的优越性，回顾性纳入2014年1月至2015年9月之间3897例行CTA检查的患者的五年随访记录进行实验。

参见图5，本组患者随访事件中位数为1442天。五年的随访结果显示，CAD-RADS 0-1级(低危组)患者的无事件生存率为97.6％，2-3级(中危组)患者的无事件生存率为70.8％，4-5分组(高危组)患者的无事件生存率为37.0％；而本发明分级系统所得低风险组(患者个体总积分＜5分)五年无事件存活率为93.4％，中风险组(患者个体总积分为5-10分)五年无事件存活率为76.5％，高风险组(患者个体总积分＞10分)五年无事件存活率为33.4％。

两种分级方式均显示，随着分组升级(危险度升高)，患者发生主要心血管不良事件(Major Adverse Cardiac Event,MACE)的可能性随时间增加。

表1.Kaplan-Meier累计概率法所得患者在不同时间点不同分级系统所得的MACE发生率

如表1所示，使用CAD-RADS分级所得低危组患者在3个月、6个月、1年、2年、3年主要不良心血管累计发生率均为0，本发明分级能更精准的为低风险组预测MACE发生率。此外，与CAD-RADS相比，本发明分级得到的中风险组患者5年累计事件发生率较低，且高风险组患者5年累计事件发生率较高，体现了本发明分级系统可更加精确的区分患者危险。

表2.单纯使用危险因素、CAD-RADS与本发明分级的预测评价

如表2所示，三种分级方式中，本发明的分级系统所得C指数最高，患者危险度分级的区分度效能最高。

表3.重分类改善指数表

如表3所示，依据是否出现MACE事件，计算本发明系统相对于CAD-RADS的重分类改善指数(net reclassification improvement,NRI)为0.210，显示出本发明分级系统更为优越的患者危险度分级能力。

总之，CAD-RADS预测患者五年无事件生存率的C指数为0.785，而本发明分级系统的C指数为0.814，同时，与CAD-RADS相比，本发明分级系统的NRI值为0.210。以上结果均证实与CAD-RADS相比，本发明的分级系统更具优越性。