• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 14
    • 下载费用:30 金币  

    重庆时时彩万能直选: 由CT影像确定主动脉缩窄处压力差的装置与方法.pdf

    关 键 词:
    CT 影像 确定 主动脉 缩窄处 压力 装置 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710021582.5

    申请日:

    2017.01.12

    公开号:

    CN106803251A

    公开日:

    2017.06.06

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 专利申请权的转移IPC(主分类):G06T 7/00登记生效日:20171218变更事项:申请人变更前权利人:西安电子科技大学变更后权利人:西安电子科技大学变更事项:地址变更前权利人:710071 陕西省西安市雁塔区太白南路2号变更后权利人:710071 陕西省西安市太白南路2号变更事项:申请人变更后权利人:广东省心血管病研究所|||实质审查的生效IPC(主分类):G06T 7/00申请日:20170112|||公开
    IPC分类号: G06T7/00(2017.01)I; G06T7/12(2017.01)I 主分类号: G06T7/00
    申请人: 西安电子科技大学
    发明人: 黄力宇; 杨茂青; 黄美萍; 李军; 庄建
    地址: 710071 陕西省西安市雁塔区太白南路2号
    优先权:
    专利代理机构: 陕西电子工业专利中心 61205 代理人: 田文英;王品华
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710021582.5

    授权公告号:

    ||||||

    法律状态公告日:

    2018.01.05|||2017.06.30|||2017.06.06

    法律状态类型:

    专利申请权、专利权的转移|||实质审查的生效|||公开

    摘要

    本发明公开了一种由CT影像确定主动脉缩窄处压力差的装置和方法。本发明的装置包括:数据读取???、二维切片???、分割???、三维重建???、主动脉缩窄判别???、特征提取???、分类???、结果显示???。本发明的步骤包括:1、读取CT数据,2、对CT数据进行二维切片处理,3、分割主动脉图像,4、测量主动脉直径比,5、判别是否存在缩窄,6、从三维主动脉模型中提取特征,7、对主动脉缩窄程度进行分类,8、显示主动脉缩窄处压差值。本发明利用机器学习算法对CT影像形态学特征进行分类,可由CT影像直接得到较为准确的主动脉缩窄处压差值,提高了CT影像在主动脉缩窄诊断中的准确性。

    权利要求书

    1.一种由CT影像确定主动脉缩窄处压力差的装置,包括数据读取???、二维切片???、
    分割???、三维重建???、主动脉缩窄判别???、特征提取???、分类???、结果显示???,
    其中:
    所述的数据读取???,用于读入格式为.dcm或.raw的原始胸部CT数据;
    所述的二维切片???,用于将读入的原始胸部CT数据分别映射为不同灰度值像素点的
    二维切片图像矩阵;
    所述的分割???,用于从二维切片图像矩阵中选择包含降主动脉末端的切片图像矩
    阵,选择包含降主动脉末端的切片图像矩阵,以所选的切片图像矩阵的中心点为中心,生成
    n*n厘米的正方形框,n的大小介于图像矩阵长度的五分之一到四分之一之间;利用S[i,j]
    =I[p,q]×G[0.25]公式,计算平滑后图像各像素点坐标处的灰度值,去除二维切片图像噪
    声,利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁边缘图像;
    将血管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管内部区域图像,利
    用公式,计算血管内部区域图像的质心;载入已标记血管区域切
    片图像矩阵之外的最下层切片图像矩阵,在切片图像矩阵上确定一个种子点元素,该种子
    点元素的横坐标与上层切片图像矩阵血管内部区域的质心横坐标相等,纵坐标与上层切片
    图像矩阵血管内部区域的质心纵坐标相等,以血管壁图像的边缘为界限,以种子点元素坐
    标确定种子点进行8邻域自适应区域生长,得到切片图像矩阵上的主动脉管内区域图像;其
    中,S[i,j]表示平滑后图像中位于[i,j]处的像素点的灰度值,I[p,q]表示二维切片图像中
    位于[p,q]处的像素点的灰度值,G[0.25]表示标准差为0.25的高斯函数,M表示血管内部区
    域图像的质心点,∫表示积分操作,f(x,y)表示血管内部区域图像像素点(x,y)处的灰度值;
    所述的三维重建???,用于将切片图像矩阵上的主动脉管内区域图像导入到三维重建
    软件中进行三维体绘制,得到绘制后的三维主动脉模型;
    所述的主动脉缩窄判别???,用于在三维主动脉模型中,每间隔1mm测量一次主动脉直
    径,分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表中,若直径
    比统计表中的值均大于0.8,则存在缩窄,否则,不存在缩窄;
    所述的特征提取???,用于将三维主动脉模型导入医学图像处理软件,将软件计算得
    到的主动脉血管的最大梯度dmax值作为特征1的值;将主动脉缩窄位于升主动脉位置时的
    标识赋值为0,将主动脉缩窄位于降主动脉位置时的标识赋值为-1,将主动脉缩窄位于主动
    脉弓位置时的标识赋值为1,将赋值后的标识作为特征2的值;将三维主动脉模型导入医学
    图像处理软件,将软件测量得到的主动脉血管的最缩窄处直径作为特征3的值;将三维主动
    脉模型导入医学图像处理软件,将软件测量并计算得到的最缩窄处直径与降主动脉直径的
    比值作为特征4的值;利用R1=πR2/(SQRT((H×W)/3600)公式,计算主动脉缩窄面积比,将
    计算结果作为特征5的值;利用R2=D/SQRT(SQRT((H×W)/3600)公式,计算主动脉缩窄比
    率,将计算结果作为特征6的值;其中,R1表示主动脉图像的缩窄面积比,π表示圆周率,R表
    示主动脉图像最缩窄处半径,/表示除法操作,SQRT表示开方操作,H表示患者身高,W表示患
    者体重,R2表示主动脉图像的缩窄比率,D表示主动脉图像中最缩窄处直径;
    所述的分类???,用于将6个特征值输入到主缩压差模型中,主缩压差模型输出与6个
    特征值相对应的主动脉缩窄处的压差值;
    所述的结果显示???,用于显示主缩压差模型得到的主动脉缩窄处压差值。
    2.一种由CT影像确定主动脉缩窄处压力差的方法,包括如下步骤:
    (1)读取CT数据:
    数据读取??槎寥敫袷轿?dcm或.raw的原始胸部CT数据;
    (2)对CT数据进行二维切片处理:
    二维切片??榻寥氲脑夹夭緾T数据,分别映射为不同灰度值像素点的二维切片图
    像矩阵;
    (3)分割主动脉图像:
    (3a)分割??榇佣衅枷窬卣笾?,选择包含降主动脉末端的切片图像矩阵,以所
    选的切片图像矩阵的中心点为中心,生成n*n厘米的正方形框,n的大小介于图像矩阵长度
    的五分之一到四分之一之间;
    (3b)利用高斯平滑公式,计算平滑后图像各像素点坐标处的灰度值,去除二维切片图
    像噪声,利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁边缘图
    像;
    (3c)将血管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管内部区
    域图像,利用Image Moments公式,计算血管内部区域图像的质心;
    (3d)载入已标记血管区域切片图像矩阵之外的最下层切片图像矩阵,利用主动脉切片
    图像分割方法,得到切片图像矩阵上的主动脉管内区域图像;
    (3e)判断是否已载入所有切片图像矩阵,若是,则执行步骤(4),否则,执行步骤(3d);
    (4)测量主动脉直径比:
    (4a)三维重建??榻衅枷窬卣笊系闹鞫龉苣谇蛲枷?,导入到三维重建软件中
    进行三维体绘制,得到绘制后的三维主动脉模型;
    (4b)在三维主动脉模型中,每间隔1mm测量一次主动脉直径;
    (4c)分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表中;
    (5)主动脉缩窄判别??榕卸现本侗韧臣票碇惺欠翊嬖谛∮?.8的值,若是,则执行步
    骤(6),否则,执行步骤(7);
    (6)从三维主动脉模型中提取特征:
    (6a)特征提取??榻鞫瞿P?,导入医学图像处理软件,将软件计算得到的主
    动脉血管的最大梯度dmax值作为特征1的值;
    (6b)将主动脉缩窄位于升主动脉位置时的标识赋值为0,将主动脉缩窄位于降主动脉
    位置时的标识赋值为-1,将主动脉缩窄位于主动脉弓位置时的标识赋值为1,将赋值后的标
    识作为特征2的值;
    (6c)将三维主动脉模型导入医学图像处理软件,将软件测量得到的主动脉血管的最缩
    窄处直径作为特征3的值;
    (6d)将三维主动脉模型导入医学图像处理软件,将软件测量并计算得到的最缩窄处直
    径与降主动脉直径的比值作为特征4的值;
    (6e)利用缩窄面积比公式,计算主动脉图像主动脉缩窄面积比,将计算结果作为特征5
    的值;
    (6f)利用缩窄比率公式,计算主动脉图像缩窄比率,将计算结果作为特征6的值;
    (7)对主动脉缩窄程度进行分类:
    (7a)分类??榻?个特征值输入到主缩压差模型中;
    (7b)主缩压差模型输出与6个特征值相对应的主动脉缩窄处的压差值;
    (8)结果显示??橄允局魉跹共钅P偷玫降闹鞫鏊跽ρ共钪?。
    3.根据权利要求2所述的由CT影像确定主动脉缩窄处压力差的方法,其特征在于,步骤
    (3b)中所述的高斯平滑公式如下:
    S[i,j]=I[p,q]×G[0.25]
    其中,S[i,j]表示平滑后图像中位于[i,j]处的像素点的灰度值,I[p,q]表示二维切片
    图像中位于[p,q]处的像素点的灰度值,G[0.25]表示标准差为0.25的高斯函数。
    4.根据权利要求2所述的由CT影像确定主动脉缩窄处压力差的方法,其特征在于,步骤
    (3c)中所述的Image Moments公式如下:
    <mrow> <mi>M</mi> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mi>&infin;</mi> </msubsup> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mi>&infin;</mi> </msubsup> <mi>x</mi> <mi>y</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow>
    其中,M表示血管内部区域图像的质心点,∫表示积分操作,f(x,y)表示血管内部区域图
    像像素点(x,y)处的灰度值。
    5.根据权利要求2所述的由CT影像确定主动脉缩窄处压力差的方法,其特征在于,步骤
    (3d)中所述的主动脉切片图像分割方法的步骤如下:
    第一步,在切片图像矩阵上确定一个种子点元素,该种子点元素的横坐标与上层切片
    图像矩阵血管内部区域的质心横坐标相等,纵坐标与上层切片图像矩阵血管内部区域的质
    心纵坐标相等;
    第二步,以血管壁图像的边缘为界限,以种子点元素坐标确定种子点进行8邻域自适应
    区域生长。
    6.根据权利要求2所述的由CT影像确定主动脉缩窄处压力差的方法,其特征在于,步骤
    (6e)中所述的缩窄面积比公式如下:
    R1=πR2/(SQRT((H×W)/3600)
    其中,R1表示主动脉图像的缩窄面积比,π表示圆周率,R表示主动脉图像最缩窄处半
    径,/表示除法操作,SQRT表示开方操作,H表示患者身高,W表示患者体重。
    7.根据权利要求2所述的由CT影像确定主动脉缩窄处压力差的方法,其特征在于,步骤
    (6f)中所述的缩窄比率公式如下:
    R2=D/SQRT(SQRT((H×W)/3600)
    其中,R2表示主动脉图像的缩窄比率,D表示主动脉图像中最缩窄处直径。

    说明书

    由CT影像确定主动脉缩窄处压力差的装置与方法

    技术领域

    本发明属于电子技术领域,更进一步涉及医学影像处理、人工智能和机器学习技
    术领域中的一种由计算机断层扫描CT(Computed Tomography)影像确定主动脉缩窄处压力
    差的装置及方法。本发明针对医学临床中诊断主动脉缩窄时,需要对压力差进行量化的情
    形,利用确定主动脉缩窄处压力差的装置,实现了由CT影像定量得到主动脉缩窄处两端的
    压力差。

    背景技术

    主动脉缩窄是先天性心脏病中一种常见的畸形,目前临床上由CT诊断主动脉缩窄
    主要通过从形态学特征对主动脉缩窄严重程度进行定性估测,要想实现准确地评估主动脉
    缩窄严重程度,需要使用压力导丝有创地伸入病灶内部,测量前后两端的压差值,这种测量
    方法不但对患者有较大的损伤,临床操作难度较大,同时费用也比较昂贵。

    中国科学院半导体研究所在其申请的专利文献“一种光纤心内压导丝”(公开号:
    CN105054916A,申请号:CN201510609653.4,申请日:2015年11月18日)中公开了一种光纤心
    内压导丝。该装置包括光纤和固定在光纤端面上的测压腔,实现了将导丝伸入人体内部测
    量压力差。该装置存在的不足是,必须将装置伸入至人体内部,有创侧压的方法对患者有一
    定的损伤,并且费用较为昂贵。

    Matthew J.Budoff,A Shittu,S Roy,在其发表的论文“Use of cardiovascular
    computed tomography in the diagnosis and management of coarctation of the
    aorta”(Journal of Thoracic&Cardiovascular Surgery,2013,146(1):229-32)中描述了
    由CT影像评估主动脉缩窄的方法。该方法首先对病人的病灶部位进行CT扫描,并将CT影像
    导入至计算机的影像系统,其次,从计算机上读取CT影像上的特征,如:缩窄的比率、缩窄所
    在的位置等,实现对主动脉缩窄的定性评估。该方法存在的不足之处在于,无法定量地评估
    主动脉缩窄的程度,不能由CT影像得到主动脉缩窄处的压力差,这一判别主动脉缩窄的金
    标准。

    苏州润心医疗科技有限公司在其申请的专利文献“基于心脏CT图像的冠状动脉血
    流储备分数计算方法”(公开号:CN106023202A,申请号:CN201610339892.7,申请日:2016年
    05月20日)中公开了一种由CT图像定量确定冠状动脉狭窄程度的方法。该方法首先提取心
    肌图像,对冠状动脉精确分割,然后对冠状动脉体数据做边缘检测,生成冠状动脉三角网格
    模型,最后实现对冠状动脉狭窄程度的定量评估。该方法存在的不足之处是,所针对的血管
    只是冠状动脉,冠状动脉生理结构与主动脉差距大,且冠状动脉狭窄往往由于血栓等堵塞
    而形成,而主动脉缩窄是由于血管生理结构的改变,上述方法不适合用来评估主动脉缩窄。

    综上所述,现有的方法只能由CT影像对主动脉缩窄情况进行粗略的估计,若要得
    到主动脉缩窄处两端的压力差只能采取有创测量的方法,对患者的损伤较大,并且费用较
    为昂贵。

    发明内容

    本发明的目的是针对医学临床中诊断主动脉缩窄时,需要对压力差进行量化的情
    形,提供了一种由CT影像直接得到缩窄处压力差的装置和方法,可无创、定量、快速地根据
    CT影像评估主动脉缩窄的严重程度,为医生的诊断和治疗提供了有力的辅助,并可降低病
    人的治疗成本。

    实现本发明目的的思路是,分割出先心病患者术前CT图像中的主动脉血管,并将
    其进行三维可视化重建,对重建得到的模型进行主动脉缩窄判定,对判定为缩窄的数据进
    一步提取特征,并使用机器学习的方法将特征映射到具体的压差值范围内,得到主动脉缩
    窄处的压差值,为临床医生提供参考。

    为实现上述目的,本发明的装置包括数据读取???、二维切片???、分割???、三
    维重建???、主动脉缩窄判别???、特征提取???、分类???、结果显示???,其中:

    所述的数据读取???,用于读入格式为.dcm或.raw的原始胸部CT数据;

    所述的二维切片???,用于将读入的原始胸部CT数据分别映射为不同灰度值像素
    点的二维切片图像矩阵;

    所述的分割???,用于从二维切片图像矩阵中选择包含降主动脉末端的切片图像
    矩阵,选择包含降主动脉末端的切片图像矩阵,以所选的切片图像矩阵的中心点为中心,生
    成n*n厘米的正方形框,n的大小介于图像矩阵长度的五分之一到四分之一之间;利用S[i,
    j]=I[p,q]×G[0.25]公式,计算平滑后图像各像素点坐标处的灰度值,去除二维切片图像
    噪声,利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁边缘图
    像;将血管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管内部区域图
    像,利用公式,计算血管内部区域图像的质心;载入已标记血管区
    域切片图像矩阵之外的最下层切片图像矩阵,在切片图像矩阵上确定一个种子点元素,该
    种子点元素的横坐标与上层切片图像矩阵血管内部区域的质心横坐标相等,纵坐标与上层
    切片图像矩阵血管内部区域的质心纵坐标相等,以血管壁图像的边缘为界限,以种子点元
    素坐标确定种子点进行8邻域自适应区域生长,得到切片图像矩阵上的主动脉管内区域图
    像;其中,S[i,j]表示平滑后图像中位于[i,j]处的像素点的灰度值,I[p,q]表示二维切片
    图像中位于[p,q]处的像素点的灰度值,G[0.25]表示标准差为0.25的高斯函数,M表示血管
    内部区域图像的质心点,∫表示积分操作,f(x,y)表示血管内部区域图像像素点(x,y)处的
    灰度值;

    所述的三维重建???,用于将切片图像矩阵上的主动脉管内区域图像导入到三维
    重建软件中进行三维体绘制,得到绘制后的三维主动脉模型;

    所述的主动脉缩窄判别???,用于在三维主动脉模型中,每间隔1mm测量一次主动
    脉直径,分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表中,若
    直径比统计表中的值均大于0.8,则存在缩窄,否则,不存在缩窄;

    所述的特征提取???,用于将三维主动脉模型导入医学图像处理软件,将软件计
    算得到的主动脉血管的最大梯度dmax值作为特征1的值;将主动脉缩窄位于升主动脉位置
    时的标识赋值为0,将主动脉缩窄位于降主动脉位置时的标识赋值为‐1,将主动脉缩窄位于
    主动脉弓位置时的标识赋值为1,将赋值后的标识作为特征2的值;将三维主动脉模型导入
    医学图像处理软件,将软件测量得到的主动脉血管的最缩窄处直径作为特征3的值;将三维
    主动脉模型导入医学图像处理软件,将软件测量并计算得到的最缩窄处直径与降主动脉直
    径的比值作为特征4的值;利用R1=πR2/(SQRT((H×W)/3600)公式,计算主动脉缩窄面积
    比,将计算结果作为特征5的值;利用R2=D/SQRT(SQRT((H×W)/3600)公式,计算主动脉缩
    窄比率,将计算结果作为特征6的值;其中,R1表示主动脉图像的缩窄面积比,π表示圆周率,
    R表示主动脉图像最缩窄处半径,/表示除法操作,SQRT表示开方操作,H表示患者身高,W表
    示患者体重,R2表示主动脉图像的缩窄比率,D表示主动脉图像中最缩窄处直径;

    所述的分类???,用于将6个特征值输入到主缩压差模型中,主缩压差模型输出与
    6个特征值相对应的主动脉缩窄处的压差值;

    所述的结果显示???,用于显示主缩压差模型得到的主动脉缩窄处压差值。

    本发明的方法包括如下步骤:

    (1)读取CT数据:

    数据读取??槎寥敫袷轿?dcm或.raw的原始胸部CT数据;

    (2)对CT数据进行二维切片处理:

    二维切片??榻寥氲脑夹夭緾T数据,分别映射为不同灰度值像素点的二维切
    片图像矩阵;

    (3)分割主动脉图像:

    (3a)分割??榇佣衅枷窬卣笾?,选择包含降主动脉末端的切片图像矩阵,
    以所选的切片图像矩阵的中心点为中心,生成n*n厘米的正方形框,n的大小介于图像矩阵
    长度的五分之一到四分之一之间;

    (3b)利用高斯平滑公式,计算平滑后图像各像素点坐标处的灰度值,去除二维切
    片图像噪声,利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁边
    缘图像;

    (3c)将血管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管内
    部区域图像,利用Image Moments公式,计算血管内部区域图像的质心;

    (3d)载入已标记血管区域切片图像矩阵之外的最下层切片图像矩阵,利用主动脉
    切片图像分割方法,得到切片图像矩阵上的主动脉管内区域图像;

    (3e)判断是否已载入所有切片图像矩阵,若是,则执行步骤(4),否则,执行步骤
    (3d);

    (4)测量主动脉直径比:

    (4a)三维重建??榻衅枷窬卣笊系闹鞫龉苣谇蛲枷?,导入到三维重建软
    件中进行三维体绘制,得到绘制后的三维主动脉模型;

    (4b)在三维主动脉模型中,每间隔1mm测量一次主动脉直径;

    (4c)分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表
    中;

    (5)主动脉缩窄判别??榕卸现本侗韧臣票碇惺欠翊嬖谛∮?.8的值,若是,则执
    行步骤(6),否则,执行步骤(7);

    (6)从三维主动脉模型中提取特征:

    (6a)特征提取??榻鞫瞿P?,导入医学图像处理软件,将软件计算得到
    的主动脉血管的最大梯度dmax值作为特征1的值;

    (6b)将主动脉缩窄位于升主动脉位置时的标识赋值为0,将主动脉缩窄位于降主
    动脉位置时的标识赋值为-1,将主动脉缩窄位于主动脉弓位置时的标识赋值为1,将赋值后
    的标识作为特征2的值;

    (6c)将三维主动脉模型导入医学图像处理软件,将软件测量得到的主动脉血管的
    最缩窄处直径作为特征3的值;

    (6d)将三维主动脉模型导入医学图像处理软件,将软件测量并计算得到的最缩窄
    处直径与降主动脉直径的比值作为特征4的值;

    (6e)利用缩窄面积比公式,计算主动脉图像主动脉缩窄面积比,将计算结果作为
    特征5的值;

    (6f)利用缩窄比率公式,计算主动脉图像缩窄比率,将计算结果作为特征6的值;

    (7)对主动脉缩窄程度进行分类:

    (7a)分类??榻?个特征值输入到主缩压差模型中;

    (7b)主缩压差模型输出与6个特征值相对应的主动脉缩窄处的压差值;

    (8)结果显示??橄允局魉跹共钅P偷玫降闹鞫鏊跽ρ共钪?。

    本发明与现有的技术相比具有以下优点:

    第一,由于本发明的装置采用CT图像分割???、CT图像特征提取???、CT图像分类
    ???,实现CT影像直接得到较为准确的主动脉缩窄处压差值,克服了现有技术导丝测压法
    对患者损伤大,临床操作难度大,费用比较昂贵的不足,使得本发明可以无创、快速地测量
    压差,降低了临床操作难度,节省了诊断成本。

    第二,由于本发明的方法采用分割原始胸部CT数据中的主动脉、从三维主动脉模
    型中提取特征、对主动脉缩窄程度进行分类,实现提取CT影像特征,得到较为准确的主动脉
    缩窄处压差值,克服了现有技术仅能从CT形态学特征对主动脉缩窄严重程度进行定性估
    测,无法定量评估的不足,使得本发明大大提高了CT影像对于诊断主动脉缩窄的准确性。

    附图说明

    图1是本发明装置的方框图;

    图2是本发明方法的流程图;

    图3是本发明方法中分割主动脉图像的流程图;

    图4是本发明方法中分割切片图像矩阵结果图,其中黑线标识区域图像为分割后
    得到的主动脉管内区域图像;

    图5是本发明方法中分割后切片的三维重建,得到的绘制后三维主动脉模型。

    具体实施方式

    下面结合附图对本发明做详细描述。

    参照附图1,对本发明的装置进行清楚、完整地描述。

    本发明的装置包括数据读取???、二维切片???、分割???、三维重建???、主动
    脉缩窄判别???、特征提取???、分类???、结果显示???。

    数据读取???,用于读入格式为.dcm或.raw的原始胸部CT数据。

    二维切片???,用于将读入的原始胸部CT数据分别映射为不同灰度值像素点的二
    维切片图像矩阵。

    分割???,用于从二维切片图像矩阵中选择包含降主动脉末端的切片图像矩阵,
    选择包含降主动脉末端的切片图像矩阵,以所选的切片图像矩阵的中心点为中心,生成n*n
    厘米的正方形框,n的大小介于图像矩阵长度的五分之一到四分之一之间;利用S[i,j]=I
    [p,q]×G[0.25]公式,计算平滑后图像各像素点坐标处的灰度值,去除二维切片图像噪声,
    利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁边缘图像;将血
    管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管内部区域图像,利用
    公式,计算血管内部区域图像的质心;载入已标记血管区域切片
    图像矩阵之外的最下层切片图像矩阵,在切片图像矩阵上确定一个种子点元素,该种子点
    元素的横坐标与上层切片图像矩阵血管内部区域的质心横坐标相等,纵坐标与上层切片图
    像矩阵血管内部区域的质心纵坐标相等,以血管壁图像的边缘为界限,以种子点元素坐标
    确定种子点进行8邻域自适应区域生长,得到切片图像矩阵上的主动脉管内区域图像;其
    中,S[i,j]表示平滑后图像中位于[i,j]处的像素点的灰度值,I[p,q]表示二维切片图像中
    位于[p,q]处的像素点的灰度值,G[0.25]表示标准差为0.25的高斯函数,M表示血管内部区
    域图像的质心点,∫表示积分操作,f(x,y)表示血管内部区域图像像素点(x,y)处的灰度值。

    三维重建???,用于将切片图像矩阵上的主动脉管内区域图像导入到三维重建软
    件中进行三维体绘制,得到绘制后的三维主动脉模型。

    主动脉缩窄判别???,用于在三维主动脉模型中,每间隔1mm测量一次主动脉直
    径,分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表中,若直径
    比统计表中的值均大于0.8,则存在缩窄,否则,不存在缩窄。

    特征提取???,用于将三维主动脉模型导入医学图像处理软件,将软件计算得到
    的主动脉血管的最大梯度dmax值作为特征1的值;将主动脉缩窄位于升主动脉位置时的标
    识赋值为0,将主动脉缩窄位于降主动脉位置时的标识赋值为‐1,将主动脉缩窄位于主动脉
    弓位置时的标识赋值为1,将赋值后的标识作为特征2的值;将三维主动脉模型导入医学图
    像处理软件,将软件测量得到的主动脉血管的最缩窄处直径作为特征3的值;将三维主动脉
    模型导入医学图像处理软件,将软件测量并计算得到的最缩窄处直径与降主动脉直径的比
    值作为特征4的值;利用R1=πR2/(SQRT((H×W)/3600)公式,计算主动脉缩窄面积比,将计
    算结果作为特征5的值;利用R2=D/SQRT(SQRT((H×W)/3600)公式,计算主动脉缩窄比率,
    将计算结果作为特征6的值;其中,R1表示主动脉图像的缩窄面积比,π表示圆周率,R表示主
    动脉图像最缩窄处半径,/表示除法操作,SQRT表示开方操作,H表示患者身高,W表示患者体
    重,R2表示主动脉图像的缩窄比率,D表示主动脉图像中最缩窄处直径。

    分类???,用于将6个特征值输入到主缩压差模型中,主缩压差模型输出与6个特
    征值相对应的主动脉缩窄处的压差值。

    结果显示???,用于显示主缩压差模型得到的主动脉缩窄处压差值。

    参照附图2对本发明的方法做进一步的详细描述。

    步骤1,读取CT影像。

    数据读取??槎寥敫袷轿?dcm或.raw的原始胸部CT数据。

    步骤2,对CT数据进行二维切片处理。

    二维切片??榻寥氲脑夹夭緾T数据分别映射为不同灰度值像素点的二维切
    片图像矩阵,采用3DMed软件将原始胸部CT数据根据灰度值映射为二维切片图像矩阵。

    步骤3,分割主动脉图像。

    第一步,分割??榇佣衅枷窬卣笾?,选择包含降主动脉末端的切片图像矩
    阵,以所选的切片图像矩阵的中心点为中心,生成n*n厘米的正方形框,n的大小介于图像矩
    阵长度的五分之一到四分之一之间。

    第二步,利用高斯平滑公式,计算平滑后图像各像素点坐标处的灰度值,去除二维
    切片图像噪声,利用医学图像软件,提取正方形框内的降主动脉末端血管边缘,得到血管壁
    边缘图像。

    所述的高斯平滑公式如下:

    S[i,j]=I[p,q]×G[0.25]

    其中,S[i,j]表示平滑后图像中位于[i,j]处的像素点的灰度值,I[p,q]表示二维
    切片图像中位于[p,q]处的像素点的灰度值,G[0.25]表示标准差为0.25的高斯函数。

    第三步,将血管壁边缘图像所包围的所有内部点的灰度值置为255,组成一个血管
    内部区域图像,利用Image Moments公式,计算血管内部区域图像的质心。

    所述的Image Moments公式如下:


    其中,M表示血管内部区域图像的质心点,∫表示积分操作,f(x,y)表示血管内部区
    域图像像素点(x,y)处的灰度值。

    第四步,载入已标记血管区域切片图像矩阵之外的最下层切片图像矩阵,在载入
    的切片图像矩阵上确定一个种子点元素,其横坐标与上层切片图像矩阵血管内部区域的质
    心横坐标相等,纵坐标与上层切片图像矩阵血管内部区域的质心纵坐标相等,以第三步得
    到的血管壁图像的边缘为界限,以种子点元素坐标确定种子点进行8邻域自适应区域生长,
    得到切片图像矩阵上的主动脉管内区域图像。

    第五步,判断是否已载入所有切片图像矩阵,若是,则执行步骤4,否则,执行第四
    步。

    步骤4,测量主动脉直径比。

    三维重建??榻衅枷窬卣笊系闹鞫龉苣谇蛲枷?,导入到三维重建软件中
    进行三维体绘制,得到绘制后的三维主动脉模型。

    在三维主动脉模型中,每间隔1mm测量一次主动脉直径。

    分别将相邻的两次测量的主动脉直径做比,并将比值结果存入直径比统计表中。

    步骤5,主动脉缩窄判别???,判断直径比统计表中是否存在小于0.8的值,若是,
    则执行步骤(6),否则,执行步骤(7)。

    步骤6,从三维主动脉模型中提取特征。

    特征提取??榻鞫瞿P偷既胍窖枷翊砣砑?,将软件计算得到的主动
    脉血管的最大梯度dmax值作为特征1的值。

    将主动脉缩窄位于升主动脉位置时的标识赋值为0,将主动脉缩窄位于降主动脉
    位置时的标识赋值为‐1,将主动脉缩窄位于主动脉弓位置时的标识赋值为1,将赋值后的标
    识作为特征2的值。

    将三维主动脉模型导入医学图像处理软件,将软件测量得到的主动脉血管的最缩
    窄处直径作为特征3的值。

    将三维主动脉模型导入医学图像处理软件,将软件测量并计算得到的最缩窄处直
    径与降主动脉直径的比值作为特征4的值。

    利用缩窄面积比公式,计算主动脉图像主动脉缩窄面积比,将计算结果作为特征5
    的值。

    所述的缩窄面积比公式如下:

    R1=πR2/(SQRT((H×W)/3600)

    其中,R1表示主动脉图像的缩窄面积比,π表示圆周率,R表示主动脉图像最缩窄处
    半径,/表示除法操作,SQRT表示开方操作,H表示患者身高,W表示患者体重。

    利用缩窄比率公式,计算主动脉图像缩窄比率,将计算结果作为特征6的值。

    所述的缩窄比率公式如下:

    R2=D/SQRT(SQRT((H×W)/3600)

    其中,R2表示主动脉图像的缩窄比率,D表示主动脉图像中最缩窄处直径。

    步骤7,对主动脉缩窄程度进行分类:

    从医院数据库中采集500例主动脉缩窄患者的病例数据组成训练集,每例病例数
    据包含患者的CT影像,身长,体重,主动脉缩窄处压差四种数据。

    使用步骤6中的方法,定量计算患者的六个特征。

    根据患者主动脉缩窄处的压差确定该训练数据的标签label值,具体为,压差从0-
    5mmHg、6-10mmHg、10-15mmHg…95-100mmHg,标签label值分别为1、2、3…20。

    将所有病例的特征及相应的label值作为训练集,输入到分类器中,所有500例数
    据导入完毕后,运行分类器得到分类模型model。

    选取100例病例作为测试集,提取每一例测试病人的六个特征,输入到分类器的
    model模型进行分类,将得到的结果与该病例实际的压差进行比较,统计分类的准确率,准
    确率在80%以上则表明该模型具有很好的可用性。

    分类??榻?个特征值输入到主缩压差模型中。

    主缩压差模型输出与6个特征值相对应的主动脉缩窄处的压差值,具体是根据主
    缩压差模型的输出分类结果确定具体的压差值,输出结果为0表示该患者主动脉缩窄压差
    为0‐5mmHg,输出为1则表示该患者主动脉缩窄压差为5–10mmHg,以此类推,输出为20表示患
    者的主动脉缩窄压差为95–100mmHg。

    步骤8,结果显示??橄允局魉跹共钅P偷玫降闹鞫鏊跽ρ共钪?。

    关于本文
    本文标题:由CT影像确定主动脉缩窄处压力差的装置与方法.pdf
    链接地址://www.4mum.com.cn/p-6001327.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 黑龙江时时官网lm0 快三稳赚 重庆时时现场开奖结果软件 秒速时时彩计划稳赚软件 时时彩定胆规律技巧 七星彩今天预测号 重庆时时三星彩走势图360 澳门骰子玩法图 北京pk赛车3码追号计划 赛车大小单双怎么看 pk10最牛稳赚技巧 聚富视界软件 11选3稳赚投注技巧 赛车6码滚雪球盈利表 大乐透中奖计算器 新时时历史开奖