CN111325761A - 基于三次b样条插值的血管内超声图像分割与测量方法 - Google Patents

基于三次b样条插值的血管内超声图像分割与测量方法 Download PDF

Info

Publication number
CN111325761A
CN111325761A CN202010118408.4A CN202010118408A CN111325761A CN 111325761 A CN111325761 A CN 111325761A CN 202010118408 A CN202010118408 A CN 202010118408A CN 111325761 A CN111325761 A CN 111325761A
Authority
CN
China
Prior art keywords
points
cubic
contour
spline interpolation
ultrasound image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010118408.4A
Other languages
English (en)
Inventor
朱国忠
童基均
黎凯
陈灵宇
李越
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yihui Jiaxing Technology Co ltd
Original Assignee
Yihui Jiaxing Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yihui Jiaxing Technology Co ltd filed Critical Yihui Jiaxing Technology Co ltd
Priority to CN202010118408.4A priority Critical patent/CN111325761A/zh
Publication of CN111325761A publication Critical patent/CN111325761A/zh
Priority to CN202110156416.2A priority patent/CN113160254A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4053Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开一种基于三次B样条插值的血管内超声图像分割与测量方法,包括:获取待处理的心血管内超声图像,并在图上校准像素的实际单位;通过手动在图上取少量的点,利用三次B样条插值函数生成分割轮廓;若需要进一步精确分割轮廓,可计算生成样条的位置控制点和弧度控制点,手动微调得到精确位置;得到精确轮廓后,可通过轮廓上的点计算得到过质心的最短内径和最长内径;若需要进一步得到精确轮廓,可手动提取内径;通过两次分割得到的内外膜轮廓,可计算得到封闭轮廓包含的面积和血管狭窄率。

Description

基于三次B样条插值的血管内超声图像分割与测量方法
技术领域
本发明涉及医学图像处理技术领域,尤其涉及一种基于三次B样条插值的血管内超声图像分割与测量方法。
背景技术
近年来,全球冠心病的发病率和死亡率呈上升趋势,冠状动脉疾病中常见的动脉粥样硬化(atherosclerosis,AS)已成为危害人类健康的首要疾病之一。通过对IVUS(心血管内超声)图像内外膜分割,不仅可对血管腔的面积进行精确的测量,而且可以发现早期的冠状动脉粥样硬化斑块,进而对斑块的偏心距、壁厚及成分含量进行检测。
目前主流的分割方法是手动或者半自动勾画IVUS图像的内外膜,但在严重的动脉粥样硬化病变时,典型的3层血管壁结构消失,中层变薄,且胶原纤维聚集,内膜增厚,这导致人工勾画IVUS图像非常费时,并且不同用户间存在差异。而当血管壁钙化时,会产生超声高反射性声影,遮盖住其后方的组织,这又使大多数半自动分割方法分割不精确。
发明内容
本发明提供一种基于三次B样条插值的血管内超声图像分割与测量方法,可以在内外膜上取少量的点,用三次B样条插值自动生成目标轮廓及最短、最长内径,比现有的纯手工勾画技术要节省很多时间。
为此,本发明所采用的具体技术方案如下:
一种基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于,包括以下步骤:
1)获取待处理的心血管内超声图像,并在图上校准像素的实际单位;
2)通过手动在图上取少量的点,利用三次B样条插值函数生成分割轮廓;
3)得到分割轮廓后,通过分割轮廓上的点计算得到过质心的最短内径和最长内径;
4)通过两次分割得到的内外膜轮廓,计算得到封闭轮廓包含的面积和血管狭窄率。
进一步地:在步骤1)中,校准在图像上选取两个点,由点坐标计算得到两点之间的像素个数,取定两点之间的实际长度,得到实际像素单位。
进一步地:所述的三次B样条插值函数采用下列公式:
B(t)=P0(1-t)3+3P1t(1-t)2+3P2(1-t)t2+P3t3,t∈[0,1]
其中,P0,P1,P2,P3,分别为轮廓上相邻的4个点,t的值由所***点的个数决定。
进一步地:在步骤2)中,每两个手动选取的点之间,***多个点。
进一步地:将所有手动取点和差值点放入集合s[i],其中i表示点的序号,i=1,2,3….N。
进一步地:手动取点需按顺时针或者逆时针的顺序。
进一步地:在步骤2)中,若需要进一步精确分割轮廓,可计算生成样条的位置控制点和弧度控制点,手动微调得到精确位置。
进一步地:得到样条轮廓的位置控制点和弧度控制点的方法为:位置控制点即为手动取点,弧度控制点需将手动所取相邻点之间用直线段连接,得到线段集合a[i],计算每条线段的中点,将相邻中点之间用直线段连接,得到线段集合b[i],将每个线段依次平移到其中点与手动取点重合,则该线段的两端点即为弧度控制点。
进一步地:在步骤3)和步骤4)中所述的计算得到轮廓面积和过质心的最短内径和最长内径的具体步骤分别为:
9.1、利用向量积求轮廓所包含图形的面积:
Figure BDA0002392194650000031
其中s[i].x和s[i].y分别为点i的横坐标和纵坐标;
9.2、利用上一步求得的面积计算轮廓质心:
Figure BDA0002392194650000032
9.3、采用如下方法计算得到过质心的最短和最长内径的两点:
9.3.1、外循环从i=1,2,3…N,循环计算s[i]与质心centroid两点构成的直线的参数,以及s[i]到质心centroid的距离,记为basic[i];
9.3.2、在外循环中的第i次中,内循环从j=1,2,3…N,循环计算s[j]到(第i次得到的)直线距离linedistance[i,j],s[j]到s[i]的距离pointdistance[j,i],以及s[j]到质心centroid的距离cendistance[j];在每次内循环第j次后,如果linedistance[i,j]小于basic[i]并且pointdistance[j,i]大于cendistance[j],则i,j为过轮廓质心的直线上的两个端点;
9.3.3、在外循环遍历i结束后,计算所有端点对的距离L[i],最短的端点对即为最短内径,最长的端点对即为最长内径。
进一步地:在计算得到内膜面积A1和外膜面积A2后,采用如下公式计算血管狭窄率:
Rate=(A2-A1)/A2*100;意义在于Rate为百分制单位,数值越接近100%,狭窄程度越大。
本发明与现有技术相比,具有的有益效果是:
(1)通过利用3次B样条插值技术,自动的生成血管内外膜的大致轮廓,结合手工微调,能够得到精确轮廓,且同时自动生成轮廓内的最短和最长内径,比现有的纯手工勾画技术要节省很多时间。
(2)利用分割好的轮廓点,可以自动地提供内外膜的面积、内径长度以及血管狭窄率。
附图说明
图1本发明实施例中内外膜手动取点后的自动分割及测量图;
图2本发明实施例中3次B样条的位置控制点和弧度控制点示例图;
图3本发明实施例中方法流程图。
具体实施方式
下面结合实施例和附图来详细说明本发明,并将本发明的实验中相关数据公开阐述,进而将本发明的方法详细透彻的传达给本来领域的技术人员。
本发明的具体实现步骤如下:
1)获取待处理的心血管内超声图像,并在图上校准像素的实际单位。
如图1所示,在图像上选取两个点,由点坐标计算得到两点之间的像素个数,取定两点之间的实际长度,得到实际像素单位。
2)通过手动在图上取少量的点,手动取点需按顺时针或者逆时针的顺序,利用三次B样条插值函数生成分割轮廓,如图2所示,每两个手动选取的点之间,***100个点,将所有手动取点和差值点放入集合s[i],其中i表示点的序号,i=1,2,3….N。若需要进一步精确分割轮廓,可计算生成B样条的位置控制点和弧度控制点,手动微调得到精确位置。
2.1)其中三次B样条插值函数采用下列公式:
B(t)=P0(1-t)3+3P1t(1-t)2+3P2(1-t)t2+P3t3,t∈[0,1]
这里P0,P1,P2,P3,分别为轮廓上相邻的4个点,t的值由所***点的个数决定。
2.2)其中得到样条轮廓的位置控制点和弧度控制点的方法为:位置控制点即为手动取点,弧度控制点需将手动所取相邻点之间用直线段连接,得到线段集合a[i],计算每条线段的中点,将相邻中点之间用直线段连接,得到线段集合b[i],将每个线段依次平移到其中点与手动取点重合,则该线段的两端点即为弧度控制点。
3)得到分割轮廓后,可计算得到封闭轮廓包含的面积、过质心的最短内径和最长内径;
4)通过两次分割得到的内外膜轮廓计算得到的内外膜面积,可计算血管狭窄率。
其中得到轮廓面积和过质心的最短内径和最长内径的具体步骤分别为:
1.利用向量积求轮廓所包含图形的面积:
Figure BDA0002392194650000051
其中s[i].x和s[i].y分别为点i的横坐标和纵坐标。
2.利用上一步求得的面积计算轮廓质心:
Figure BDA0002392194650000052
3.利用上一步求得的质心和轮廓上的点集s[i],由于轮廓是不规则图形,因此采用如下方法计算得到过质心的最短和最长内径的两点:
3.1)外循环从i=1,2,3…N,循环计算s[i]与质心centroid两点构成的直线的参数,以及s[i]到质心centroid的距离,记为basic[i];
3.2)在外循环中的第i次中,内循环从j=1,2,3…N,循环计算s[j]到(第i次得到的)直线距离linedistance[i,j],s[j]到s[i]的距离pointdistance[j,i],以及s[j]到质心centroid的距离cendistance[j]。在每次内循环第j次后,如果linedistance[i,j]小于basic[i]并且pointdistance[j,i]大于cendistance[j],则i,j为过轮廓质心的直线上的两个端点;
3.3)在外循环遍历i结束后,计算所有端点对的距离L[i],最短的端点对即为最短内径,最长的端点对即为最长内径;
3.4)用如下公式计算血管狭窄率:
利用步骤1的面积公式求得内膜轮廓所包含的面积A1和外膜轮廓所包含的面积A2,再计算血管狭窄率:Rate=(A2-A1)/A2*100;意义在于Rate为百分制单位,数值越接近100%,狭窄程度越大。
以上所述,仅为本发明的较佳实施举例,并不用于限制本发明,凡在本发明精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于,包括以下步骤:
1)获取待处理的心血管内超声图像,并在图上校准像素的实际单位;
2)通过手动在图上取少量的点,利用三次B样条插值函数生成分割轮廓;
3)得到分割轮廓后,通过分割轮廓上的点计算得到过质心的最短内径和最长内径;
4)通过两次分割得到的内外膜轮廓,计算得到封闭轮廓包含的面积和血管狭窄率。
2.如权利要求1所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:在步骤1)中,校准在图像上选取两个点,由点坐标计算得到两点之间的像素个数,取定两点之间的实际长度,得到实际像素单位。
3.如权利要求1所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:所述的三次B样条插值函数采用下列公式:
B(t)=P0(1-t)3+3P1t(1-t)2+3P2(1-t)t2+P3t3,t∈[0,1]
其中,P0,P1,P2,P3,分别为轮廓上相邻的4个点,t的值由所***点的个数决定。
4.如权利要求1所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:在步骤2)中,每两个手动选取的点之间,***多个点。
5.如权利要求4所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:将所有手动取点和差值点放入集合s[i],其中i表示点的序号,i=1,2,3….N。
6.如权利要求2所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:手动取点需按顺时针或者逆时针的顺序。
7.如权利要求1所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:在步骤2)中,若需要进一步精确分割轮廓,可计算生成样条的位置控制点和弧度控制点,手动微调得到精确位置。
8.如权利要求7所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:得到样条轮廓的位置控制点和弧度控制点的方法为:位置控制点即为手动取点,弧度控制点需将手动所取相邻点之间用直线段连接,得到线段集合a[i],计算每条线段的中点,将相邻中点之间用直线段连接,得到线段集合b[i],将每个线段依次平移到其中点与手动取点重合,则该线段的两端点即为弧度控制点。
9.如权利要求1所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:在步骤3)和步骤4)中所述的计算得到轮廓面积和过质心的最短内径和最长内径的具体步骤分别为:
9.1、利用向量积求轮廓所包含图形的面积:
Figure FDA0002392194640000021
其中s[i].x和s[i].y分别为点i的横坐标和纵坐标;
9.2、利用上一步求得的面积计算轮廓质心:
Figure FDA0002392194640000022
9.3、采用如下方法计算得到过质心的最短和最长内径的两点:
9.3.1、外循环从i=1,2,3…N,循环计算s[i]与质心centroid两点构成的直线的参数,以及s[i]到质心centroid的距离,记为basic[i];
9.3.2、在外循环中的第i次中,内循环从j=1,2,3…N,循环计算s[j]到(第i次得到的)直线距离linedistance[i,j],s[j]到s[i]的距离pointdistance[j,i],以及s[j]到质心centroid的距离cendistance[j];在每次内循环第j次后,如果linedistance[i,j]小于basic[i]并且pointdistance[j,i]大于cendistance[j],则i,j为过轮廓质心的直线上的两个端点;
9.3.3、在外循环遍历i结束后,计算所有端点对的距离L[i],最短的端点对即为最短内径,最长的端点对即为最长内径。
10.如权利要求9所述基于三次B样条插值的血管内超声图像分割与测量方法,其特征在于:在计算得到内膜面积A1和外膜面积A2后,采用如下公式计算血管狭窄率:
Rate=(A2-A1)/A2*100;意义在于Rate为百分制单位,数值越接近100%,狭窄程度越大。
CN202010118408.4A 2020-02-26 2020-02-26 基于三次b样条插值的血管内超声图像分割与测量方法 Pending CN111325761A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010118408.4A CN111325761A (zh) 2020-02-26 2020-02-26 基于三次b样条插值的血管内超声图像分割与测量方法
CN202110156416.2A CN113160254A (zh) 2020-02-26 2021-02-04 基于三次b样条插值的血管内超声图像分割与测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010118408.4A CN111325761A (zh) 2020-02-26 2020-02-26 基于三次b样条插值的血管内超声图像分割与测量方法

Publications (1)

Publication Number Publication Date
CN111325761A true CN111325761A (zh) 2020-06-23

Family

ID=71173022

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202010118408.4A Pending CN111325761A (zh) 2020-02-26 2020-02-26 基于三次b样条插值的血管内超声图像分割与测量方法
CN202110156416.2A Pending CN113160254A (zh) 2020-02-26 2021-02-04 基于三次b样条插值的血管内超声图像分割与测量方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202110156416.2A Pending CN113160254A (zh) 2020-02-26 2021-02-04 基于三次b样条插值的血管内超声图像分割与测量方法

Country Status (1)

Country Link
CN (2) CN111325761A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112669449A (zh) * 2020-12-31 2021-04-16 浙江理工大学 基于3d重建技术的cag和ivus精准联动分析方法及***
CN113469941A (zh) * 2021-05-27 2021-10-01 武汉楚精灵医疗科技有限公司 一种在超声胆胰管检查中胆胰管宽度的测量方法
WO2022109905A1 (zh) * 2020-11-25 2022-06-02 苏州润迈德医疗科技有限公司 调节血管轮廓及中心线的方法及存储介质
CN117982225A (zh) * 2024-04-07 2024-05-07 柏意慧心(杭州)网络科技有限公司 左心室流出道csept几何参数的测量方法、装置、设备及介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7940974B2 (en) * 2006-11-21 2011-05-10 General Electric Company Method and system for adjusting 3D CT vessel segmentation
CN101283929B (zh) * 2008-06-05 2010-06-16 华北电力大学 一种血管三维模型的重建方法
CN103455999B (zh) * 2012-06-05 2018-05-08 北京工业大学 基于内超声图像序列的血管壁边缘自动检测方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022109905A1 (zh) * 2020-11-25 2022-06-02 苏州润迈德医疗科技有限公司 调节血管轮廓及中心线的方法及存储介质
CN112669449A (zh) * 2020-12-31 2021-04-16 浙江理工大学 基于3d重建技术的cag和ivus精准联动分析方法及***
CN113469941A (zh) * 2021-05-27 2021-10-01 武汉楚精灵医疗科技有限公司 一种在超声胆胰管检查中胆胰管宽度的测量方法
CN113469941B (zh) * 2021-05-27 2022-11-08 武汉楚精灵医疗科技有限公司 一种在超声胆胰管检查中胆胰管宽度的测量方法
CN117982225A (zh) * 2024-04-07 2024-05-07 柏意慧心(杭州)网络科技有限公司 左心室流出道csept几何参数的测量方法、装置、设备及介质

Also Published As

Publication number Publication date
CN113160254A (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN111325761A (zh) 基于三次b样条插值的血管内超声图像分割与测量方法
JP6792768B2 (ja) 血流予備量比の計算方法およびシステム
US11783498B2 (en) Feature-based registration method
AU2010257203B2 (en) Fast anatomical mapping using ultrasound images
CN107115108B (zh) 快速计算血管压力差的方法及***
MX2007012541A (es) Registro de imagenes de un organo usando rasgos anatomicos fuera del organo.
CN108742587B (zh) 基于病史信息获取血流特征值的方法及装置
JP7074952B2 (ja) 血管の圧力差の修正方法、装置及び設備
CN109009061B (zh) 基于血压修正获取血流特征值的计算方法及装置
CN108742570B (zh) 基于冠脉优势类型获取血管压力差的装置
JP2021529577A (ja) 血管圧差を取得する方法及び装置
WO2012140147A2 (en) Processing ultrasound images
WO2023151280A1 (zh) 双模冠脉血管图像三维融合的方法和融合***
CN112535466A (zh) 一种基于血管影像的血流储备分数计算方法
Gurmeric et al. A new 3-D automated computational method to evaluate in-stent neointimal hyperplasia in in-vivo intravascular optical coherence tomography pullbacks
JP2010502245A (ja) 左心室のねじれを測定するシステム及び方法
CN108742667B (zh) 基于身体质量指数获取血流特征值的方法和装置
CN109009037B (zh) 基于个体性别获取血管压力差的计算方法及装置
CN116524134B (zh) 基于血管内图像的三维血管建模方法及ffr计算方法和***
CN108784676B (zh) 基于年龄信息获取压力差的方法及装置
US10332255B2 (en) Method for assessing stenosis severity in a lesion tree through stenosis mapping
CN108742586B (zh) 基于糖尿病病史信息获取血流特征值的方法及装置
WO2020210948A1 (zh) 获取冠脉血流量及血流速度的方法和装置
CN112336378B (zh) 一种用于动物超声诊断的m型超声心动图处理方法和***
US20230000461A1 (en) Ultrasound slice enhancement

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200623

WD01 Invention patent application deemed withdrawn after publication