CN109447967B - 一种颅内动脉瘤图像的分割方法及系统 - Google Patents

一种颅内动脉瘤图像的分割方法及系统 Download PDF

Info

Publication number
CN109447967B
CN109447967B CN201811260220.2A CN201811260220A CN109447967B CN 109447967 B CN109447967 B CN 109447967B CN 201811260220 A CN201811260220 A CN 201811260220A CN 109447967 B CN109447967 B CN 109447967B
Authority
CN
China
Prior art keywords
image
intracranial
blood vessel
tumor
intracranial aneurysm
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.)
Active
Application number
CN201811260220.2A
Other languages
English (en)
Other versions
CN109447967A (zh
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.)
Union Strong Beijing Technology Co ltd
Original Assignee
Union Strong Beijing 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 Union Strong Beijing Technology Co ltd filed Critical Union Strong Beijing Technology Co ltd
Priority to CN201811260220.2A priority Critical patent/CN109447967B/zh
Publication of CN109447967A publication Critical patent/CN109447967A/zh
Application granted granted Critical
Publication of CN109447967B publication Critical patent/CN109447967B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • 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/10116X-ray 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/30016Brain
    • 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/30096Tumor; Lesion
    • 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)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本说明书实施例提供一种颅内动脉瘤图像的分割方法及系统。该方法包括:从待分割的颅内载瘤血管图像上,截取局部三维图像;获取所述局部三维图像中载瘤血管图像的最大内接圆,计算颅内载瘤血管图像的中心线和半径;基于所述颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割。该方法提供了一种分割精度较好、分割效率较高的颅内动脉瘤图像的分割方法,实现了颅内动脉瘤图像的自动分割。

Description

一种颅内动脉瘤图像的分割方法及系统
技术领域
本说明书涉及医学影像领域,尤其涉及一种颅内动脉瘤图像的分割方法及系统。
背景技术
颅内动脉瘤是由颅内动脉内腔的局部异常扩张所致动脉壁的一种瘤状突起,是一种常见的血管性疾病。据报道,颅内未破裂动脉瘤在我国成人中患病率高达7%,破裂后造成蛛网膜下腔出血,可导致严重残疾或死亡。2014年国家统计局数据显示,急性脑血管病是我国人口第二大死亡原因。动脉瘤性蛛网膜下腔出血是继缺血性脑卒中和高血压脑出血之后最常见的急性脑血管疾病,死残率高达64%,约15%的病人院前死亡,不同经济发展水平地区的救治水平差异很大,已经成为引起我国居民死亡的最常见原因之一。由此可见,颅内动脉瘤的临床教学研究,对于颅内动脉瘤的治疗,具有促进意义。
现有技术中,对颅内动脉瘤图像进行分割,去除载瘤血管图像,得到纯粹的颅内动脉瘤图像,常基于操作者的经验,采用手动方法进行分割,分割精度较低,而且分割速度慢。
因此,需要一种颅内动脉瘤图像的分割方法,能够提高颅内动脉瘤图像的分割精度及分割效率。
发明内容
本说明书实施例提供一种颅内动脉瘤图像的分割方法及系统,用于解决以下技术问题:颅内动脉瘤分割精度低、分割速度慢。
为解决上述技术问题,本说明书实施例是这样实现的:
本说明书实施例提供的一种颅内动脉瘤图像的分割方法,包括:
从待分割的颅内载瘤血管图像上,截取局部三维图像;
获取所述局部三维图像中载瘤血管图像的最大内接圆,计算颅内载瘤血管图像的中心线和半径;
基于所述颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割。
进一步地,所述从待分割的颅内载瘤血管图像上,截取局部三维图像,具体包括:
从待分割的颅内载瘤血管图像上,确定种子点和两个定位点,获得种子点坐标和定位点坐标;
根据种子点坐标和定位点坐标,将包含部分或全部颅内动脉瘤图像和部分颅内载瘤血管图像的局部三维图像区域截取出来,截取时以种子点坐标和定位点坐标确定的最小长方体为截取范围。
进一步地,所述获取所述局部三维图像中载瘤血管图像的最大内接圆,计算颅内载瘤血管图像的中心线和半径,具体包括:
在所述局部三维图像所确定的颅内载瘤血管图像腔体内,沿着颅内载瘤血管图像,逐点计算颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
进一步地,所述基于所述颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割,具体包括
以中心线和种子点为初始条件,对所述截取的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像;
利用颅内载瘤血管图像的中心线和半径,进行颅内载瘤血管图像表面重建;
利用所述待分割的颅内载瘤血管图像对膨胀后的颅内动脉瘤图像进行分割,然后利用所述重建的颅内载瘤血管表面图像对膨胀后的颅内动脉瘤图像进行分割,以去除载瘤血管图像;
对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割。
进一步地,所述对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割,具体包括:
利用种子点坐标,采用区域生长方法,对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割。
进一步地,所述对切割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割,具体包括:
以种子点坐标和剥离掉载瘤血管的颅内动脉瘤图像进行区域生长,完成颅内动脉瘤图像与颅内载瘤血管图像的分割。
本说明书实施例提供一种颅内动脉瘤图像的分割系统,包括:
输入接口,用于待分割的颅内载瘤血管图像的输入;
处理工作站,基于颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割;
输出单元,将分割的颅内动脉瘤图像进行输出。
进一步地,所述基于颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割,包括:
从待分割的颅内载瘤血管图像上,截取局部三维图像;
获取所述局部三维图像中载瘤血管图像的最大内接圆,计算颅内载瘤血管图像的中心线和半径;
基于所述颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割。
进一步地,所述从待分割的颅内载瘤血管图像上,截取局部三维图像,具体包括:
从待分割的颅内载瘤血管图像上,确定种子点和两个定位点,获得种子点坐标和定位点坐标;
根据种子点坐标和定位点坐标,将包含部分或全部颅内动脉瘤图像和部分颅内载瘤血管图像的局部三维图像区域截取出来,截取时以种子点坐标和定位点坐标确定的最小长方体为截取范围。
进一步地,所述获取所述局部三维图像中载瘤血管的最大内接圆,计算颅内载瘤血管图像的中心线和半径,具体包括:
在所述局部三维图像所确定的颅内载瘤血管图像腔体内,沿着颅内载瘤血管图像,逐点计算颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
进一步地,所述基于所述颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割,具体包括
以中心线和种子点为初始条件,对所述截取的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像,该膨胀后的颅内动脉瘤图像包含完整颅内动脉瘤的图像;
利用颅内载瘤血管图像的中心线和半径,进行颅内载瘤血管图像表面重建;
利用所述待分割的颅内载瘤血管图像对膨胀后的颅内动脉瘤图像进行分割,然后利用所述重建的颅内载瘤血管表面图像对膨胀后的颅内动脉瘤图像进行分割,以去除载瘤血管图像;
对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割。
进一步地,所述对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割,具体包括:
利用种子点坐标,采用区域生长方法,对分割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割。
进一步地,所述对切割后的颅内动脉瘤图像进行重建,完成颅内动脉瘤图像的分割,具体包括:
以种子点坐标和剥离掉载瘤血管的颅内动脉瘤图像进行区域生长,完成颅内动脉瘤图像与颅内载瘤血管图像的分割。
本说明书实施例采用的上述至少一个技术方案能够达到以下有效效果:本说明书实施例采用基于中心线的分割算法,实现颅内动脉瘤图像的自动分割,能够快速进行颅内动脉瘤图像的自动分割,提高颅内动脉瘤图像分割精度及分割效率。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本说明书提供的一种颅内动脉瘤图像的分割方法的示意图;
图2为本说明书提供的二维空间中两点确定最小长方形示意图;
图3为本说明书提供的二维空间中三点确定最小长方形示意图;
图4为本说明书提供的一种颅内动脉瘤图像的分割方法的效果示意图;
图5为本说明书提供的一种颅内动脉瘤图像的分割方法的流程图;
图6为本说明书提供的一种颅内动脉瘤图像的分割系统的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本说明书实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1为本说明书提供的一种颅内动脉瘤图像的分割方法的示意图。该方法包括:
步骤S101:从待分割的颅内载瘤血管图像上,截取局部三维图像。
DSA和MRA图像是常见的颅内动脉瘤图像,但由于该图像中含有背景等干扰因素,不利于临床教学中进行研究,因此,需要进行颅内动脉瘤图像的处理。
DSA(Digital Subtraction Angiography,数字减影血管造影技术)是一种使X射线序列图片中的血管可视化的技术。其基本原理是将注入造影剂前后拍摄的两帧X射线图像经数字化输入图像计算机,通过减影、增强和再成像过程来获得清晰的纯血管影像,同时实时地呈现血管影。DSA技术在图像质量、判断血流方向和优势供血等方面是其他检查手段所不能比拟的,因此被称为血管性疾病诊断的“金标准”。
MRA(Magnetic Resonance Angiography,磁共振血管造影),指磁共振血管成像,磁共振可以行血管造影,即显示血管,可发现血管狭窄和闭塞的部位。其基本原理是基于饱和效应、流入增强效应、流动去相位效应。MRA是将预饱和带置于3D层块的头端以饱和静脉血流,反向流动的动脉血液进入3D层块,因未被饱和从而产生MR信号。扫描时将一个较厚容积分割成多个薄层激发,减少激发容积厚度以减少流入饱和效应,且能保证扫描容积范围,获得数层相邻层面的薄层图像,使图像清晰,血管的细微结构显示好,空间分辨力提高。
基于DSA的三维图像,通过区域生长的办法,可以实现颅内载瘤血管图像的分割,该方法可以有效减少噪声干扰,提升预算效率。而基于MRA的三维图像,通过区间二值的方法,可以实现颅内载瘤血管图像的分割,该方法可以快速准确的分割颅内载瘤血管图像,使血管边界更加清晰。
利用上述方法得到的分割的颅内载瘤血管图像,需要进一步进行处理,实现颅内载瘤血管上动脉瘤图像的分割,进而进行颅内动脉瘤的教学研究。
利用上述得到的颅内载瘤血管图像,确定种子点和定位点。种子点和定位点均为空间坐标,为了便于区分,将生长的起始点定义为种子点,将选取在载瘤血管图像上的点定义为定位点。种子点的选取,可以选在动脉瘤体图像表面,也可以选取在动脉瘤体图像之内。而定位点选在与动脉瘤图像相交的载瘤血管图像上面。由于颅内动脉瘤包括常规侧边瘤和分叉血管瘤,因此定位点的确定,要根据颅内动脉瘤的类型,采用不同的定位点确定方法。对于常规侧边瘤而言,需要在颅内载瘤血管的上游和下游给出两点即可,一般选取离颅内动脉瘤图像5-10mm范围内选取两点;对于分叉血管瘤,则需要在颅内载瘤血管图像的上游给出一个定位点,下游在各分支上分别给出一个定位点,共三个定位点即可。其中,上游定位点为定位点1,下游定位点为定位点2,对于分叉血管,下游定位点包括两个定位点。定位点可以放置在颅内载瘤血管图像的表面或载瘤血管图像以内都可以,两者没有区别。
局部三维图像的截取,是根据种子点坐标和定位点坐标,确定的最小长方体,进行横向及纵向的像素增量延伸,使延伸后能够包括全部颅内动脉瘤图像,用延伸后确定的长方体区域截取局部三维图像。
图2为本说明书提供的二维空间中两点确定最小长方形示意图。在三维空间中按照类似方法,根据两个定位点及种子点,确定最小长方体。
图3为本说明书提供的二维空间中三点确定最小长方形示意图。在三维空间中,按照类似方法,根据三个定位点及种子点,确定最小长方体。
步骤S102:获取所述局部三维图像中载瘤血管图像的最大内接圆,计算颅内载瘤血管图像的中心线和半径。
颅内载瘤血管图像的中心线的提取,对颅内动脉瘤图像的分割至关重要。本说明书实施例利用种子点坐标和定位点坐标,计算颅内载瘤血管图像的中心线和半径。首先根据种子点坐标和定位点坐标,将包含部分或全部颅内动脉瘤图像和部分载瘤血管图像的局部三维图像区域截取出来。截取时以种子点和定位点确定的最小长方体为截取范围。然后在上述局部三维图像区域所确定的颅内载瘤血管图像腔体内,沿着颅内载瘤血管图像,逐点计算颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为两个定位点间颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
步骤S103:基于颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割。
本说明书实施例是利用颅内动脉血管图像切割颅内动脉瘤图像,来实现颅内动脉瘤图像的分割。因此首先需要进行动脉血管图像的重建和动脉瘤图像的生成。具体过程如下:利用前述获取的种子点坐标和定位点坐标,截取局部三维图像,以中心线和种子点为初始条件,对步骤S102截取的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像。兼顾计算效率和动脉瘤大小,上述预设值可以选择为16次,上述局部三维图像膨胀16次后,得到包含完整颅内动脉瘤图像的膨胀后的颅内动脉瘤图像;以前述计算出来的颅内动脉瘤图像的中心线和半径,以中心线上的点为球心,以该点处的半径为半径,沿着中心线逐点生成颅内载瘤血管图像表面。利用上述得到的颅内载瘤血管表面图像和原始待分割动脉瘤图像,切割颅内动脉瘤图像,实现颅内载瘤血管图像与颅内动脉瘤图像的剥离。以种子点坐标和剥离掉载瘤血管的动脉瘤图像进行区域生长,实现颅内动脉瘤图像与颅内载瘤血管图像的分割,得到完整、干净的颅内动脉瘤图像。
图4为本说明书提供的颅内动脉瘤图像分割的效果示意图。从该示意图中,可以看到,利用本说明书提供的方法,可以实现动脉瘤图像的准确分割。
图5为本说明书提供的一种颅内动脉瘤图像的分割方法的流程图。
步骤S501:从待分割的颅内载瘤血管图像上,确定种子点坐标和定位点坐标。
DSA和MRA图像作为颅内动脉瘤的常用图像形式,由于存在背景等干扰因素,因此在临床教学研究等应用中,需要首先进行颅内载瘤血管图像的分割。基于DSA的三维图像,通过区域生长的办法,可以实现颅内载瘤血管图像的分割;而基于MRA的三维图像,通过区间二值的方法,可以实现颅内载瘤血管图像的分割。
在本说明书的实施例中,种子点坐标和定位点坐标的确定,可以实现载瘤血管图像中心线和半径的计算,生成颅内动脉瘤图像。根据前述步骤S102的方法确定种子点坐标和定位点坐标,从颅内动脉瘤图像上,截取局部三维图像。
步骤S502:形态学膨胀,生成膨胀的颅内动脉瘤图像。
形态学膨胀的目的是为了获得包含全部颅内动脉瘤图像的图像。具体实现时,采用如下方法:以中心线和种子点为初始条件,利用前述获得的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像。兼顾计算效率和动脉瘤大小,上述预设值可以选择为16次,上述局部三维图像膨胀16次后,得到包含完整颅内动脉瘤图像的膨胀后的颅内动脉瘤图像。
步骤S503:利用种子点坐标和定位点坐标,计算颅内载瘤血管图像的中心线和半径。
从上述截取的局部图像包含的颅内载瘤血管图像腔体内,沿着颅内载瘤血管图像,逐点计算颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为两个定位点间颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
步骤S504:基于上述颅内载瘤血管图像的中心线和半径,进行颅内载瘤血管图像的重建。
为了有效进行动脉瘤图像的分割,需要进行颅内载瘤血管图像的重建。基于上述颅内载瘤血管图像的中心线和半径,以中心线上的点为球心,以该点处的半径为半径,沿中心线逐点生成颅内载瘤血管图像表面,进行颅内载瘤血管图像表面重建。
步骤S505:用重建的颅内载瘤血管图像表面,切割前述得到的完整的颅内动脉瘤图像,得到分割的颅内动脉瘤图像。
前述利用中心线,进行了颅内载瘤血管图像的重建,用重建的颅内载瘤血管图像表面,切割前述得到的完整的颅内动脉瘤图像,去除载瘤血管图像,实现颅内载瘤血管图像与颅内动脉瘤图像的分割,得到分割的颅内动脉瘤图像。
步骤S506:区域生长,得到分割的颅内动脉瘤图像。
上述分割的颅内动脉瘤图像,由于存在边界不清等干扰,因此需要进一步对分割的颅内动脉瘤图像进行处理,去除干扰因素,以得到干净、完整的颅内动脉瘤图像。上述得到的分割的颅内动脉瘤图像为二值图像,因此采用区域生长的方法,以种子点坐标为生长点,获得连通的动脉瘤图像,最终得到分割完整、干净的颅内动脉瘤图像。
图6为本说明书提供的一种颅内动脉瘤图像的分割系统的示意图。该分割系统具体包括:
输入接口,用于待分割的颅内载瘤血管图像的输入;
处理工作站,基于颅内载瘤血管图像的中心线和半径,进行颅内动脉瘤图像的分割;
输出单元,将分割的颅内动脉瘤图像进行输出。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置、电子设备、非易失性计算机存储介质实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书实施例提供的装置、电子设备、非易失性计算机存储介质与方法是对应的,因此,装置、电子设备、非易失性计算机存储介质也具有与对应方法类似的有益技术效果,由于上面已经对方法的有益技术效果进行了详细说明,因此,这里不再赘述对应装置、电子设备、非易失性计算机存储介质的有益技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本说明书一个或多个实施例时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本领域内的技术人员应明白,本说明书实施例可提供为方法、系统、或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书是参照根据本说明书实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本说明书实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (8)

1.一种颅内动脉瘤图像的分割方法,其特征在于,包括以下步骤:
从待分割的颅内载瘤血管图像上,截取局部三维图像,具体包括:从所述待分割的颅内载瘤血管图像上,确定种子点和两个定位点,获得种子点坐标和定位点坐标;根据所述种子点坐标和所述定位点坐标,将包含部分或全部颅内动脉瘤图像和部分颅内载瘤血管图像的局部三维图像区域截取出来,截取时以所述种子点坐标和所述定位点坐标确定的最小长方体为截取范围;
获取所述局部三维图像中载瘤血管图像的最大内接圆,计算所述颅内载瘤血管图像的中心线和半径;
以所述中心线和种子点为初始条件,对所述截取的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像;
利用所述颅内载瘤血管图像的中心线和半径,进行所述颅内载瘤血管图像表面重建;
利用所述待分割的颅内载瘤血管图像对所述膨胀后的颅内动脉瘤图像进行分割,然后利用所述重建的颅内载瘤血管表面图像对膨胀后的颅内动脉瘤图像进行分割,以去除载瘤血管图像;
对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割。
2.如权利要求1所述的方法,其特征在于,所述获取所述局部三维图像中载瘤血管图像的最大内接圆,计算所述颅内载瘤血管图像的中心线和半径,具体包括:
在所述局部三维图像所确定的颅内载瘤血管图像腔体内,沿着所述颅内载瘤血管图像,逐点计算所述颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为所述颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
3.如权利要求1所述的方法,其特征在于,所述对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割,具体包括:
利用种子点坐标,采用区域生长方法,对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割。
4.如权利要求1所述的方法,其特征在于,所述对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割,具体包括:
以种子点坐标和剥离掉载瘤血管的颅内动脉瘤图像进行区域生长,完成所述颅内动脉瘤图像与所述颅内载瘤血管图像的分割。
5.一种颅内动脉瘤图像的分割系统,其特征在于,包括:
输入接口,用于待分割的颅内载瘤血管图像的输入;
处理工作站,从所述待分割的颅内载瘤血管图像上,截取局部三维图像,具体包括:从所述待分割的颅内载瘤血管图像上,确定种子点和两个定位点,获得种子点坐标和定位点坐标;根据所述种子点坐标和所述定位点坐标,将包含部分或全部颅内动脉瘤图像和部分颅内载瘤血管图像的局部三维图像区域截取出来,截取时以所述种子点坐标和所述定位点坐标确定的最小长方体为截取范围;
获取所述局部三维图像中载瘤血管图像的最大内接圆,计算所述颅内载瘤血管图像的中心线和半径;
以所述中心线和种子点为初始条件,对所述截取的局部三维图像,进行形态学膨胀,得到膨胀后的颅内动脉瘤图像;
利用所述颅内载瘤血管图像的中心线和半径,进行所述颅内载瘤血管图像表面重建;
利用所述待分割的颅内载瘤血管图像对膨胀后的颅内动脉瘤图像进行分割,然后利用所述重建的颅内载瘤血管表面图像对所述膨胀后的颅内动脉瘤图像进行分割,以去除载瘤血管图像;
对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割;
输出单元,将分割后的颅内动脉瘤图像进行输出。
6.如权利要求5所述的系统,其特征在于,所述获取所述局部三维图像中载瘤血管图像的最大内接圆,计算所述颅内载瘤血管图像的中心线和半径,具体包括:
在所述局部三维图像所确定的颅内载瘤血管图像腔体内,沿着所述颅内载瘤血管图像,逐点计算所述颅内载瘤血管图像腔体内的内接圆,以最大内接圆的圆心连线作为所述颅内载瘤血管图像的中心线,以最大内接圆的半径作为该点处的颅内载瘤血管图像的半径。
7.如权利要求5所述的系统,其特征在于,所述对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割,具体包括:
利用种子点坐标,采用区域生长方法,对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割。
8.如权利要求5所述的系统,其特征在于,所述对分割后的颅内动脉瘤图像进行重建,完成所述颅内动脉瘤图像的分割,具体包括:
以种子点坐标和剥离掉载瘤血管的颅内动脉瘤图像进行区域生长,完成所述颅内动脉瘤图像与所述颅内载瘤血管图像的分割。
CN201811260220.2A 2018-10-26 2018-10-26 一种颅内动脉瘤图像的分割方法及系统 Active CN109447967B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811260220.2A CN109447967B (zh) 2018-10-26 2018-10-26 一种颅内动脉瘤图像的分割方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811260220.2A CN109447967B (zh) 2018-10-26 2018-10-26 一种颅内动脉瘤图像的分割方法及系统

Publications (2)

Publication Number Publication Date
CN109447967A CN109447967A (zh) 2019-03-08
CN109447967B true CN109447967B (zh) 2022-03-01

Family

ID=65548665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811260220.2A Active CN109447967B (zh) 2018-10-26 2018-10-26 一种颅内动脉瘤图像的分割方法及系统

Country Status (1)

Country Link
CN (1) CN109447967B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109345585B (zh) * 2018-10-26 2021-11-30 强联智创(北京)科技有限公司 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN111785381B (zh) * 2020-07-27 2024-03-29 北京市神经外科研究所 一种支架模拟方法、装置以及设备
CN111815622B (zh) * 2020-07-27 2024-04-02 北京市神经外科研究所 一种支架模拟中心线的优化方法、装置以及设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1973301A (zh) * 2004-06-23 2007-05-30 皇家飞利浦电子股份有限公司 用于显示与3d管状对象的参数有关的信息的图像处理系统
CN101290685A (zh) * 2007-04-20 2008-10-22 美国西门子医疗解决公司 冠状动脉的三维建模
CN102239495A (zh) * 2008-12-04 2011-11-09 皇家飞利浦电子股份有限公司 用于采集医学图像数据的方法、装置和计算机程序产品

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101152104A (zh) * 2006-04-21 2008-04-02 美国西门子医疗解决公司 半自动主动脉瘤分析用的系统和方法
US9830427B2 (en) * 2011-06-20 2017-11-28 Siemens Healthcare Gmbh Method for intracranial aneurysm analysis and endovascular intervention planning
CN103198202B (zh) * 2012-12-19 2017-05-17 首都医科大学 颅内动脉瘤介入治疗支架植入图像仿真方法
US9367667B2 (en) * 2014-04-30 2016-06-14 Siemens Aktiengesellschaft Method and system for advanced aneurysm analysis
JP5890055B1 (ja) * 2015-07-09 2016-03-22 株式会社アルム 血管画像処理装置、血管画像処理プログラム、および血管画像処理方法
CN107392872B (zh) * 2017-07-27 2020-04-28 强联智创(北京)科技有限公司 一种微导管塑形器的生成方法及生成系统
CN108030550B (zh) * 2017-12-26 2020-05-01 成都真实维度科技有限公司 一种基于虚拟成像的动脉瘤的瘤颈角度计算方法
CN108335304B (zh) * 2018-02-07 2021-08-03 华侨大学 一种腹部ct扫描序列图像的主动脉瘤分割方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1973301A (zh) * 2004-06-23 2007-05-30 皇家飞利浦电子股份有限公司 用于显示与3d管状对象的参数有关的信息的图像处理系统
CN101290685A (zh) * 2007-04-20 2008-10-22 美国西门子医疗解决公司 冠状动脉的三维建模
CN102239495A (zh) * 2008-12-04 2011-11-09 皇家飞利浦电子股份有限公司 用于采集医学图像数据的方法、装置和计算机程序产品

Also Published As

Publication number Publication date
CN109447967A (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN109493348B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN109345585B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN109472780B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN109472823B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN109448004B (zh) 一种基于中心线的颅内血管图像的截取方法及系统
CN109448003B (zh) 一种颅内动脉血管图像分割方法及系统
CN109389637B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN109447967B (zh) 一种颅内动脉瘤图像的分割方法及系统
CN109584997B (zh) 一种颅内动脉瘤图像的形态学参数的测量方法及系统
CN111081378B (zh) 一种动脉瘤破裂风险评估方法及系统
CN109584169A (zh) 一种基于中心线的颅内血管图像的截取方法及系统
CN111127428A (zh) 一种基于颅脑影像数据的目标区域的提取方法及系统
CN111863263B (zh) 一种模拟方法、装置以及设备
CN111091563A (zh) 一种基于颅脑影像数据的目标区域的提取方法及系统
CN111223089B (zh) 一种动脉瘤的检测方法、装置及计算机可读存储介质
CN111105404B (zh) 一种基于颅脑影像数据的目标位置的提取方法及系统
CN109584261B (zh) 一种颅内动脉瘤图像的分割方法及系统
CN110503642B (zh) 一种基于dsa图像的定位方法及系统
CN109472803B (zh) 一种颅内动脉血管分割方法及系统
CN109377504B (zh) 一种颅内动脉血管图像分割方法及系统
CN110517244B (zh) 一种基于dsa图像的定位方法及系统
CN112734726B (zh) 一种血管造影的分型方法、装置以及设备
CN111863262B (zh) 一种模拟方法、装置以及设备
CN111862062B (zh) 一种中心线优化的方法、装置以及设备
CN111785381B (zh) 一种支架模拟方法、装置以及设备

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
GR01 Patent grant
GR01 Patent grant