买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：西北工业大学

摘要：本发明公开了一种面向CTA影像的智能冠状动脉狭窄评估方法，首先，以语义分割网络U‑Net为基础，提出了一种基于注意力残差网络的AttentionRes‑UNet冠脉狭窄分割网络模型用于冠脉血管的自动分割。其次，利用冠脉CPR原图与上一步得到的分割图相结合进行分类预测，在Res‑Net50基础模型上进行改进，提出了一种多尺度卷积神经网络Res2Net冠脉狭窄分类网络模型，用于判断冠脉血管的狭窄程度。本发明按照美国CASS冠脉病变血管诊断标准将冠状动脉狭窄类型分为轻度狭窄、中度狭窄和重度狭窄，以此实现对冠状动脉狭窄的自动评估。

主权项：1.一种面向CTA影像的智能冠状动脉狭窄评估方法，其特征在于，包括如下步骤：步骤1：基于最短路径分支删除的冠状动脉中心线提取；步骤1-1：采用体素十二方向的删除点，十二个方向代表立方体的十二条边，将它们分为四组，每组包含三个方向；步骤1-2：将三维骨架线转化为具有有向边和节点的有向图，其中，只有单个相邻点的体素点被称为端点，两个相邻点之间的体素点被称为普通链接点，若有多于两个相邻点，则被视为节点；普通链接点之间的边构成一段有向边；根据有向图中节点数量N构建n×n的邻接矩阵Z，其表达式如式1所示： 在冠状动脉中心线的搜索过程中，第i个节点到第j个节点的路径长度被用表示为zij；若节点A与节点B之间的路径长度比经过节点C的路径长度更短，则表明节点A和节点B是相邻的；当搜索完所有节点后，将它们依次连接起来，得到冠状动脉中心线；步骤2：多角度三线性插值曲面重建算法；步骤2-1：基于Frenet框架的曲面重建算法；将冠状动脉的中心线视为中心点在整个三维空间中的运动轨迹；Frenet-Serret公式如式2所示 其中，t表示单位切向量；N表示单位法向量，是通过将弧长参数进行微分而得出的矢量；B作为副法向量，是由切向量T和法向量N的外积构成；dds是用来计算弧长的积分，κ和τ分别为所描述曲线的曲率和挠率；在冠脉血管进行曲面重建时，为了计算每个中心点的Frenet坐标，将中心线视为离散的中心点，使其离散化；假设冠脉血管中心线存在三个连续点M1、M2、M3，将三个点的坐标分别记作[x1，y1，z1]、[x2，y2，z2]、[x3，y3，z3]，并通过点与点之间的关系来确定每个点的切向量，从而得到每个点处的平均切向量，从而实现对冠脉血管的拉直曲面重建，其表达式如式3所示： 通过计算切向量和法向量之间的正交关系，可以将它们的乘积记作点的法向量，其表达式如式4所示： 将切向量T和法向量N的平均值记为该点的副法向量，其表达式如式5所示： 通过对三个向量进行标准正交化，计算出中心点M2处的Frenet坐标，根据该坐标即可获得中心点M2的灰度值；通过遍历中心线内的所有中心点，并使用Frenet公式计算出它们在TNB坐标系中的位置和灰度值，最后将这些灰度值按顺序拼接起来，就可以得到冠脉血管曲面重建后的图像；步骤2-2：三线性插值算法；三线性插值需要与其相邻的八个点进行插值；点Px,y,z在一个三维空间中，其中点P的左边分量点x分别向上和向下取整后得到xL与xH；同理，分别对分量y和分量z取整，得到分量坐标后，对y方向进行插值，先对V1进行计算，计算过程如式5所示：V1＝yH-yIxH，yH,zH5+y-yLIxH，yH,zH同理计算出V2、V3和V4，根据计算结果对x方向进行线性插值，计算过程如式6所示：V5＝xH-xV4+x-xLV3V6＝xH-xV1+x-xLV26最后再对z方向进行插值计算即可得到点P处的灰度值，计算过程如式7所示：Vo＝zH-zV6+z-zLV57在提取到所有点的灰度值后，再对中心线数据进行均匀采样，删除重复的中心点；确保采样后的中心点数据的最小间隔为单个体素块；步骤2-3：多角度曲面重建算法；改变重建体数据的感兴趣矢量对冠状动脉进行曲面重建；沿着X轴方向从0°到180°中均匀选取五个方向作为映射方向对冠脉血管进行曲面重建；步骤3：构建AttentionRes-UNet冠脉狭窄分割网络；AttentionRes-UNet冠脉狭窄分割网络是在U-Net网络框架上添入残差模块和空间注意力机制；首先将残差模块替换基础网络中编码和解码阶段中的卷积模块；其次，在解码阶段，将注意力机制添加到跳跃连接之前；最后通过一个1×1的卷积进行通道调整，输出模型最终的预测结果；步骤3-1：AttentionRes-UNet中的残差模块；将U-Net网络原始编码阶段和解码阶段的卷积模块改为残差模块；AttentionRes-UNet中的残差模块中包括三个卷积层，第一个卷积层的作用是对特征通道进行降维，然后在第二个11卷积层中恢复特征通道，其次利用残差块恒等映射的方法来计算，最后将其叠加，以此作为残差模块的输出；步骤3-2：AttentionRes-UNet中的注意力机制；将AttentionGate模块通过跳跃连接集成到U-Net网络的上采样模块中，实现空间上的注意力机制，突出CPR图像中需要分割的血管区域，抑制背景中的噪声、运动伪影以及其他可能存在的血管区域；1空间注意力使用空间注意力机制SAM，将注意力分配至冠脉血管区域，从而实现增强目标区域并且抑制背景区域的目的；首先对一个H×W×C的特征F进行平均池化操作与最大池化操作，在这个过程中得到两个新的H×W×1的特征，将二者分别记为同时将这两个新特征拼接在一块；再对其进行一个7×7的卷积操作和Sigmod激活函数，即得到权重系数Ms；最后将Ms与输入的特征F相乘，得到的结果即为比输入特征F拥有更多冠脉血管位置的有效特征F′；SAM权重系数Ms计算公式如式8，其中，f7*7代表经过7×7的卷积层； 2AttentionRes-UNet的注意力机制；AttentionRes-UNet算法在U-Net网络的解码阶段引入注意力模块，在跳跃连接之前添加注意力机制；设xl为编码过程中所提取到的特征，gl为解码过程中所提取到的特征，首先将xl和gl相加，使得二者重复区域拥有更大的激活值；其次利用空间注意力SAM计算出注意力权重系数σ1，从而获得冠脉CPR图中要分割的血管特征信息，抑制背景中的噪音以及运动伪影等无用信息；注意力机制的输出特征为编码阶段xl与ɑl按元素相乘的结果最后将该特征与解码阶段获得的特征进行拼接融合；注意力机制输出特征的计算公式如式9所示： 采用加性注意力来获得ɑl，加性注意力公式如式10与式11所示：Ql＝ψσ1Wxxl+Wggl+bg+bψ10al＝σ2Ql11其中，σ1和σ2分别对应ReLu激活函数和Sigmoid激活函数，Wx与Wg为注意力权重，bg与bψ为注意力中的偏执项；步骤4：基于多尺度卷积神经网络Res2Net冠脉狭窄分类网络模型；在Res-Net50网络框架里构建多通道的残差连接，并采用前期融合方法将多个与目标相关的特征进行融合；首先，将冠脉CPR原图与分割后的图像进行融合，使网络可以获得更多特征从而提高分类预测的准确率；其次，对Res-Net50网络中的残差模块进行重构，增大每一层网络的感受野；步骤4-1：多尺度卷积神经网络Res2Net中的特征融合；将冠脉血管原始CPR图像与分割后的冠脉血管CPR图像进行特征融合，增加输入到网络中的特征维度，从而使网络获得更多冠脉血管的特征信息，提高图像分类的准确率；采用前期融合方法对冠脉血管狭窄程度分类网络进行特征融合：在输入层上做融合，先融合多层特征，再输入到深度学习的网络当中；步骤4-2：多尺度卷积神经网络Res2Net中的Res2Net模块；采用残差模块Res2Net，在Res-Net网络中构建出具有等级制的残差连接；Res2Net模块中，将3×3的卷积层替换成分组小卷积，并且通过类似残差连接的方式，将不同的卷积组逐层串连起来，形成一种新的神经网络模块；在Res2Net模块中，特征图x首先经过1×1的卷积后，被分为s个子集，每个子集用xi表示，其中i∈{1,2,3,…,s}；拆分后每个子集的通道数是原来的1s，空间大小与原有的特征图保持一致；除了第一个特征子集外，每个特征子集xi都有与之对应的3×3卷积层，用Ki表示，输出定义为yi；在处理过程中，为了减少参数并增加s的数量，x1不经过3×3的卷积处理，即y1＝x1，其余特征子集xi和Ki-1相加后一同被送入Ki进行处理；因此yi由公式12表示：

全文数据：

权利要求：

百度查询： 西北工业大学 一种面向CTA影像的智能冠状动脉狭窄评估方法