CN111281410A - 一种基于造影技术的脑血管疾病成像方法及应用的造影剂 - Google Patents

一种基于造影技术的脑血管疾病成像方法及应用的造影剂 Download PDF

Info

Publication number
CN111281410A
CN111281410A CN202010174701.2A CN202010174701A CN111281410A CN 111281410 A CN111281410 A CN 111281410A CN 202010174701 A CN202010174701 A CN 202010174701A CN 111281410 A CN111281410 A CN 111281410A
Authority
CN
China
Prior art keywords
contrast agent
cerebrovascular
polylactic acid
imaging method
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.)
Granted
Application number
CN202010174701.2A
Other languages
English (en)
Other versions
CN111281410B (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.)
Affiliated Hospital of Guangdong Medical University
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN202010174701.2A priority Critical patent/CN111281410B/zh
Publication of CN111281410A publication Critical patent/CN111281410A/zh
Application granted granted Critical
Publication of CN111281410B publication Critical patent/CN111281410B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/501Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61KPREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
    • A61K49/00Preparations for testing in vivo
    • A61K49/04X-ray contrast preparations
    • A61K49/0433X-ray contrast preparations containing an organic halogenated X-ray contrast-enhancing agent
    • A61K49/0438Organic X-ray contrast-enhancing agent comprising an iodinated group or an iodine atom, e.g. iopamidol
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge 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/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • 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

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Computer Graphics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Neurosurgery (AREA)
  • Neurology (AREA)
  • Software Systems (AREA)
  • Pulmonology (AREA)
  • Epidemiology (AREA)
  • Physiology (AREA)
  • Geometry (AREA)

Abstract

本发明提供了一种基于造影技术的脑血管疾病成像方法及应用的造影剂。所述造影剂为以含有金纳米粒子的聚乳酸为油相,以含有非离子型碘造影剂的水溶液为水相,以2‑[甲氧基(聚氧乙烯)丙基]‑三甲氧基硅烷为改性剂,制备成的硅氧烷修饰的聚乳酸微泡。采用上述造影剂静脉注射后,采集不同层面的脑血管横截面二维图像,然后依次对图像进行分阶段线性灰度变换和边缘检测算子滤除软组织噪声处理,再采用基于二阶微分的拉普拉斯算子进行血管分割,最后将二维血管图像进行三维重建得到脑血管三维图像。本发明的造影剂造影增强效果明显,且副作用小,再结合本发明的成像方法,能够得到清晰的脑血管三维图像。

Description

一种基于造影技术的脑血管疾病成像方法及应用的造影剂
技术领域
本发明属于医学图像处理及医学制剂制备技术领域,尤其涉及一种基于造影技术的脑血管疾病成像方法及应用的造影剂。
背景技术
脑血管疾病泛指脑部血管的各种疾病,包括脑动脉粥样硬化、血栓形成、狭窄、闭塞、脑动脉炎、脑动脉损伤、脑动脉瘤、颅内血管畸形、脑动静脉瘘等,其共同特点是引起脑组织的缺血或出血性意外,导致患者的残废或死亡。因此,需要通过提前检测和诊断脑血管中的危险因素,来预防心脑血管疾病的发生或者复发,从而达到降低死亡率的目的。
脑血管造影是观察脑部血管最直接的方法,能了解血管的形态、分布、粗细、移位、闭塞、狭窄等,尚可观察血管本身的病变(如脑动脉瘤和脑动静脉畸形(AVM)等)。电子计算机断层扫描(Computed Tomography,CT)是利用X射线束对人体组织器官进行一定厚度的层面进行扫描,具有扫描时间快,图像清晰等优点,广泛用于医学诊断中。由于CT的成像特点,其对实体组织如骨骼、肝脏等具有较好的成像功能,但对血管等组织的成像能力较差,为了增强CT对血管等组织的成像效果,CT造影剂开始被开发和应用。传统CT造影剂主要利用了碘溶液吸收X射线的特性,一定程度增强了CT造影效果,但大多存在生物毒性较大、造影效果不佳、制作工艺复杂等缺点。近年来,随着纳米材料技术的飞速发展,各种性能优越(如吸收能力强、散射截面大、表面等离子体共振效应明显等)的纳米材料不断被研制出来。因此,以新型的纳米材料作为造影剂进行辅助增强成像的相关研究受到了科学家们的广泛关注。
申请号为200910113084.9的中国发明专利公开了一种CT成像造影剂及其制备方法,该造影剂由纳米金颗粒为核心,外围包裹碘海醇,虽然增强了对X射线的吸收系数,但仍然存在注射剂量大、血液半衰期短的问题。而且,采用CT采集注射造影剂后的脑血管图像时,随着采集时间的延长,造影剂不可避免的会渗出至血管周围的软组织中,导致血管与周围软组织的界限变模糊,从而影响血管分割的精度。
发明内容
针对上述现有技术存在的缺陷,本发明的目的在于提供一种基于造影技术的脑血管疾病成像方法。采用X射线断层扫描采集注射造影剂后的脑血管横截面二维图像,然后依次对图像进行灰度变换、图像增强、血管分割和三维重建得到脑血管三维图像。
本发明的目的还在于提供一种上述脑血管疾病成像方法应用的造影剂,所述造影剂为以含有金纳米粒子的聚乳酸为油相,以含有非离子型碘造影剂的水溶液为水相,以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷为改性剂,制备成的硅氧烷修饰的聚乳酸微泡。
为实现上述目的,本发明采用如下的技术方案:
一种基于造影技术的脑血管疾病成像方法,包括以下步骤:
S1.从静脉注射造影剂,待造影剂到达脑血管区域后,对脑血管区域进行X射线断层扫描,得到不同层面的脑血管横截面二维图像;
S2.对步骤S1所述每个层面的脑血管横截面二维图像依次进行灰度变换和图像增强处理,得到每个层面的脑血管横截面预处理图像;
S3.对步骤S2所述每个层面的脑血管横截面预处理图像进行血管分割,得到每个层面的脑血管图像;
S4.将步骤S3所述每个层面的脑血管图像进行三维重建,得到脑血管三维图像。
作为本发明的进一步改进,在步骤S2中,所述灰度变换采用分段线性变换。
作为本发明的进一步改进,在步骤S2中,所述图像增强处理是利用梯度算子|▽t|、曲率检测因子λ2和血管多尺度因子σ,设计边缘检测算子,对血管进行多尺度边缘检测,滤除没有边缘结构的含有造影剂的软组织噪声。
作为本发明的进一步改进,所述边缘检测算子的公式如下:
Figure BDA0002410398250000031
作为本发明的进一步改进,在步骤S3中,所述血管分割采用基于二阶微分的拉普拉斯算子对脑血管横截面预处理图像进行分割。
作为本发明的进一步改进,在步骤S4中,所述三维重建采用光线投射法。
一种以上所述的基于造影技术的脑血管疾病成像方法中应用的造影剂,所述造影剂为以含有金纳米粒子的聚乳酸为油相,以含有非离子型碘造影剂的水溶液为水相,以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷为改性剂,制备成的硅氧烷修饰的聚乳酸微泡。
作为本发明的进一步改进,所述造影剂的制备方法包括以下步骤:
(1)将金纳米粒子加入到含有聚乳酸的三氯甲烷溶液中,作为油相;
(2)将非离子型碘造影剂加入到注射用水,得到质量分数为70%-90%的造影剂水溶液,作为水相;
(3)将步骤(2)所述水相加入到步骤(1)所述油相中,超声处理得到聚乳酸微泡;
(4)将步骤(3)所述聚乳酸微泡加入到浓度为4%-8%(w/v)的2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷的乙醇溶液中,使溶液中的乙醇完全挥发掉,此过程中,2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷在所述聚乳酸微泡表面固化交联,得到硅氧烷修饰的聚乳酸微泡。
作为本发明的进一步改进,在步骤(1)中,所述油相中金纳米粒子的浓度为2-6mg/ml,聚乳酸的浓度为30-60mg/ml。
作为本发明的进一步改进,在步骤(3)中,所述聚乳酸微泡和乙醇溶液的体积比为1:2.5-1:4。
有益效果
与现有技术相比,本发明提供的基于造影技术的脑血管疾病成像方法及应用的造影剂具有如下有益效果:
(1)本发明提供的基于造影技术的脑血管疾病成像方法,采用CT造影成像的方法,首先采用自制的造影剂静脉注射后,采集不同层面的脑血管横截面二维图像,然后依次对图像进行分阶段线性灰度变换和边缘检测算子滤除软组织噪声处理,再采用基于二阶微分的拉普拉斯算子进行血管分割,最后将二维血管图像进行三维重建得到脑血管三维图像。如此操作,采用分阶段线性灰度变换,能够对脑血管横截面二维图像进行灰度改善,提升血管像素的灰度值,压缩非血管像素的灰度值,从而有效的抑制了非血管像素的灰度值,保留了过渡区间的灰度值,并一定程度上增强了血管目标的灰度值。采用边缘检测算子滤除软组织区域,由于血管在边缘区域会产生灰度的陡变,陡变度可以利用梯度算子进行检测,从而得到边缘明确的血管,很大程度的平滑血管边缘,凸显血管结构。当检测到血管内的像素时,曲率检测因子λ2的值相对较大,导致梯度T较小,从而不影响血管的增强;当检测到血管外的像素点时,由于软组织没有固定的边缘和走向,所以不存在曲率检测因子λ2,从而通过梯度T过滤这些像素点,实现了软组织噪声的去除。采用二阶微分的拉普拉斯算子对于血管图像的细节可以产生更强的响应,能够分割得到高精度的血管图像。
(2)本发明提供的脑血管疾病成像方法应用的造影剂,将金纳米粒子聚乳酸溶液作为油相,非离子型碘造影剂水溶液作为水相,得到金纳米粒子和聚乳酸包裹非离子型碘造影剂的聚乳酸微泡,再以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷修饰聚乳酸微泡表面,能够提高造影剂在血液中的生物相容性和稳定性。金纳米粒子作为聚乳酸微泡的外层,不仅能够增强造影,而且能够延长体内血液循环时间从而克服了传统碘造影剂排除快的缺陷。另外,金纳米粒子很容易与各种功能生物分子作用,获得组织靶向或多模成像功能,从而可以为血管疾病的诊断提供准确、丰富的信息。
附图说明
图1为本发明提供的基于造影技术的脑血管疾病成像方法的流程框图。
具体实施方式
以下将对本发明各实施例的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例;基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例,都属于本发明所保护的范围。
CT造影是选用X射线束对造影部位以垂直照射,并利用高灵敏度的射线探测仪对造影部位以某一厚度进行环绕平扫,得到穿透人体部位后衰减过的射线束,再利用光电转换理论将接受到的射线束转换为电信号,再利用模/数变换,将电信号进行转换,最后以数字图像形式输出,得到CT图像。一般地,一次扫描可以得到连续的厚度相同的多张横切图像,这些横切图像可以进行三维重建等操作。在生成CT图像过程中,模拟信号转换的数字信号中,每一个数字信号都可以对应灰度范围内的某一灰度值,单个数字信号在灰度图像中映射为像素,像素时CT二维横切图像的处理单元。在一列横切图像中,选取连续几个大小相同、位置对应的像素构成块,将序列图像分为完全相同的长方体,这些长方体称为体素,图像进行三维操作时,是以体素为单元进行处理的。
基于此,本发明提供的一种基于造影技术的脑血管疾病成像方法如图1所示,包括以下步骤:
S1.从静脉注射造影剂,待造影剂到达脑血管区域后,对脑血管区域进行X射线断层扫描,得到不同层面的脑血管横截面二维图像;
S2.对步骤S1所述每个层面的脑血管横截面二维图像依次进行灰度变换和图像增强处理,得到每个层面的脑血管横截面预处理图像;
S3.对步骤S2所述每个层面的脑血管横截面预处理图像进行血管分割,得到每个层面的脑血管图像;
S4.将步骤S3所述每个层面的脑血管图像进行三维重建,得到脑血管三维图像。
作为本发明的进一步改进,在步骤S2中,所述灰度变换采用如式(1)所示的分段线性函数对脑血管横截面二维图像进行分段线性变换:
Figure BDA0002410398250000061
如对原始输入的脑血管横截面二维图像I(x,y)进行变换时,a、b代表需要调整的灰度范围区间的端点值,c、d代表调整好的灰度范围区间的端点,决定对脑血管横截面二维图像中某个灰度区间的像素进行拉伸或者是压缩,是灰度变换的重点。采用此种变换,能够对脑血管横截面二维图像进行灰度改善,提升血管像素的灰度值,压缩非血管像素的灰度值,从而有效的抑制了非血管像素的灰度值,保留了过渡区间的灰度值,并一定程度上增强了血管目标的灰度值。
当造影剂在血管中流动时,造影剂不可避免的会渗出血管到血管周围的软组织,如造影剂渗出达到一定程度时,血管周围的组织也会在CT造影过程中被保留下来,成为影响血管影像的主要干扰因素,尤其在相比较来说比较细的血管区域这种干扰更加明显。渗透到软组织中的造影剂会影响视觉效果,但是,含有造影剂的软组织是扩散式分布的,没有固定的区域边缘,血管却是具有明确的边缘结构。根据这个特性,本发明设计一个边缘检测算子专门针对血管边缘结构进行检测,滤除没有边缘结构的含有造影剂的软组织区域。
因此,作为本发明的进一步改进,在步骤S2中,所述图像增强处理是利用梯度算子|▽t|、曲率检测因子λ2和血管多尺度因子σ,设计边缘检测算子,对血管进行多尺度边缘检测,滤除没有边缘结构的含有造影剂的软组织噪声。所述边缘检测算子的公式如式(2)所示:
Figure BDA0002410398250000071
采用式(2)进行边缘检测增强时,由于血管在边缘区域会产生灰度的陡变,陡变度可以利用梯度算子进行检测,从而得到边缘明确的血管,很大程度的平滑血管边缘,凸显血管结构。当检测到血管内的像素时,曲率检测因子λ2的值相对较大,导致梯度T较小,从而不影响血管的增强;当检测到血管外的像素点时,由于软组织没有固定的边缘和走向,所以不存在曲率检测因子λ2,从而通过梯度T过滤这些像素点,实现了软组织噪声的去除。
作为本发明的进一步改进,在步骤S3中,所述血管分割采用基于二阶微分的拉普拉斯算子,如式(3)所示,对脑血管横截面预处理图像进行分割:
Figure BDA0002410398250000072
优选地,在采用基于二阶微分的拉普拉斯算子进行血管分割时,先采用高斯滤波器对脑血管横截面预处理图像进行平滑处理,如式(4)所示:
Figure BDA0002410398250000073
式中,σ为图像的标准差,在实际平滑过程中,σ过小,则偏离中心的所有像素权重会非常的小,当σ小到一定程度时,滤波操作则会退化为对图像进行点运算,无法起到平滑噪声的作;但如果σ过大,邻域较小时,在邻域内的高斯模板将会退化为均值模板,会造成边缘模糊的效果,所以本发明σ值优选为0.6-0.7,有助于达到较好的平滑效果。
作为本发明的进一步改进,在步骤S4中,所述三维重建采用光线投射法。该方法的思想是在投影平面对每个像素以投影角度分别发出射线,穿过体数据,对每条射线上数据值进行比较,选一个最大的值作为响应像素显示值。最大密度投影法最主要的优点是计算量小,计算速度快。
一种以上所述的基于造影技术的脑血管疾病成像方法中应用的造影剂,所述造影剂为以含有金纳米粒子的聚乳酸为油相,以含有非离子型碘造影剂的水溶液为水相,以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷为改性剂,制备成的硅氧烷修饰的聚乳酸微泡。
作为本发明的进一步改进,所述造影剂的制备方法包括以下步骤:
(1)将金纳米粒子加入到含有聚乳酸的三氯甲烷溶液中,作为油相;
(2)将非离子型碘造影剂加入到注射用水,得到质量分数为70%-90%的造影剂水溶液,作为水相;
(3)将步骤(2)所述水相加入到步骤(1)所述油相中,超声处理得到聚乳酸微泡;
(4)将步骤(3)所述聚乳酸微泡加入到浓度为4%-8%(w/v)的2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷的乙醇溶液中,使溶液中的乙醇完全挥发掉,此过程中,2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷在所述聚乳酸微泡表面固化交联,得到硅氧烷修饰的聚乳酸微泡。
作为本发明的进一步改进,在步骤(1)中,所述油相中金纳米粒子的浓度为2-6mg/ml,聚乳酸的浓度为30-60mg/ml。
作为本发明的进一步改进,在步骤(3)中,所述聚乳酸微泡和乙醇溶液的体积比为1:2.5-1:4。
将金纳米粒子聚乳酸溶液作为油相,非离子型碘造影剂水溶液作为水相,得到金纳米粒子和聚乳酸包裹非离子型碘造影剂的聚乳酸微泡。再以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷修饰聚乳酸微泡表面,能够提高造影剂在血液中的生物相容性和稳定性。金纳米粒子不仅具有增强造影的能力,而且它能够延长体内血液循环时间从而克服了传统碘造影剂排除快的缺陷。另外,金纳米粒子很容易与各种功能生物分子作用,获得组织靶向或多模成像功能,从而可以为血管疾病的诊断提供准确、丰富的信息。
实施例1
请参阅图1所示,一种基于造影技术的脑血管疾病成像方法,包括以下步骤:
S1.从静脉注射造影剂,待造影剂到达脑血管区域后,对脑血管区域进行X射线断层扫描,得到不同层面的脑血管横截面二维图像;所述造影剂通过以下步骤制备:
(1)将金纳米粒子加入到含有聚乳酸的三氯甲烷溶液中,作为油相;所述油相中金纳米粒子的浓度为4mg/ml,聚乳酸的浓度50mg/ml;
(2)将碘海醇加入到注射用水,得到质量分数为80%的造影剂水溶液,作为水相;
(3)将步骤(2)所述水相滴加到步骤(1)所述油相中,滴加完毕后立即用探头超声对样品进行超声处理,得到聚乳酸微泡;
(4)将步骤(3)所述聚乳酸微泡加入到浓度为6%(w/v)的2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷的乙醇溶液中(所述聚乳酸微泡和乙醇溶液的体积比为1:3),置于通风橱中搅拌,使溶液中的乙醇完全挥发掉,在此过程中,2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷在所述聚乳酸微泡表面固化交联,得到硅氧烷修饰的聚乳酸微泡。
S2.采用下式所示的分段线性函数对对步骤S1所述每个层面的脑血管横截面二维图像进行分段线性变换:
Figure BDA0002410398250000101
式中,a、b代表需要调整的灰度范围区间的端点值,c、d代表调整好的灰度范围区间的端点;
采用采用下式所示的边缘检测算子对进行边缘检测增强处理,得到每个层面的脑血管横截面预处理图像:
Figure BDA0002410398250000102
式中,|▽t|为梯度算子、λ2为曲率检测因子,σ为血管多尺度因子;
S3.采用下式所示的高斯滤波函数对步骤S2所述每个层面的脑血管横截面预处理图像进行平滑处理:
Figure BDA0002410398250000103
式中,σ为图像的标准差,取值为0.7;
然后采用下式所示的基于二阶微分的拉普拉斯算子对平滑后的每个层面的脑血管横截面预处理图像进行血管分割,得到每个层面的脑血管图像:
Figure BDA0002410398250000104
S4.采用光线投射法将步骤S3所述每个层面的脑血管图像进行三维重建,得到脑血管三维图像。
综上所述,本发明采用自制的造影剂,将金纳米粒子聚乳酸溶液作为油相,非离子型碘造影剂水溶液作为水相,得到金纳米粒子和聚乳酸包裹非离子型碘造影剂的聚乳酸微泡,再以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷修饰聚乳酸微泡表面,能够提高造影剂在血液中的生物相容性和稳定性,并能够增强造影,延长体内血液循环时间。以此造影剂进行CT造影成像,对采集的二维图像进行分阶段线性灰度变换和边缘检测算子滤除软组织噪声处理,再采用基于二阶微分的拉普拉斯算子进行血管分割,最后将二维血管图像进行三维重建得到脑血管三维图像。能够有效的抑制非血管像素的灰度值,提升血管像素的灰度值,并有效滤除血管周围的软组织区域,从而得到高精度的脑血管二维图像。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种基于造影技术的脑血管疾病成像方法,其特征在于,包括以下步骤:
S1.从静脉注射造影剂,待造影剂到达脑血管区域后,对脑血管区域进行X射线断层扫描,得到不同层面的脑血管横截面二维图像;
S2.对步骤S1所述每个层面的脑血管横截面二维图像依次进行灰度变换和图像增强处理,得到每个层面的脑血管横截面预处理图像;
S3.对步骤S2所述每个层面的脑血管横截面预处理图像进行血管分割,得到每个层面的脑血管图像;
S4.将步骤S3所述每个层面的脑血管图像进行三维重建,得到脑血管三维图像。
2.根据权利要求1所述的基于造影技术的脑血管疾病成像方法,其特征在于,在步骤S2中,所述灰度变换采用分段线性变换。
3.根据权利要求1所述的基于造影技术的脑血管疾病成像方法,其特征在于,在步骤S2中,所述图像增强处理是利用梯度算子|▽t|、曲率检测因子λ2和血管多尺度因子σ,设计边缘检测算子,对血管进行多尺度边缘检测,滤除没有边缘结构的含有造影剂的软组织噪声。
4.根据权利要求3所述的基于造影技术的脑血管疾病成像方法,其特征在于,所述边缘检测算子的公式如下:
Figure FDA0002410398240000011
5.根据权利要求1所述的基于造影技术的脑血管疾病成像方法,其特征在于,在步骤S3中,所述血管分割采用基于二阶微分的拉普拉斯算子对脑血管横截面预处理图像进行分割。
6.根据权利要求1所述的基于造影技术的脑血管疾病成像方法,其特征在于,在步骤S4中,所述三维重建采用光线投射法。
7.一种权利要求1至6中任一项权利要求所述的基于造影技术的脑血管疾病成像方法中应用的造影剂,其特征在于,所述造影剂为以含有金纳米粒子的聚乳酸为油相,以含有非离子型碘造影剂的水溶液为水相,以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷为改性剂,制备成的硅氧烷修饰的聚乳酸微泡。
8.根据权利要求7所述的脑血管疾病成像方法中应用的造影剂,其特征在于,所述造影剂的制备方法包括以下步骤:
(1)将金纳米粒子加入到含有聚乳酸的三氯甲烷溶液中,作为油相;
(2)将非离子型碘造影剂加入到注射用水,得到质量分数为70%-90%的造影剂水溶液,作为水相;
(3)将步骤(2)所述水相加入到步骤(1)所述油相中,超声处理得到聚乳酸微泡;
(4)将步骤(3)所述聚乳酸微泡加入到浓度为4%-8%(w/v)的2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷的乙醇溶液中,使溶液中的乙醇完全挥发掉,此过程中,2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷在所述聚乳酸微泡表面固化交联,得到硅氧烷修饰的聚乳酸微泡。
9.根据权利要求8所述的脑血管疾病成像方法中应用的造影剂,其特征在于,在步骤(1)中,所述油相中金纳米粒子的浓度为2-6mg/ml,聚乳酸的浓度为30-60mg/ml。
10.根据权利要求8所述的脑血管疾病成像方法中应用的造影剂,其特征在于,在步骤(3)中,所述聚乳酸微泡和乙醇溶液的体积比为1:2.5-1:4。
CN202010174701.2A 2020-03-13 2020-03-13 一种基于造影技术的脑血管疾病成像方法及应用的造影剂 Expired - Fee Related CN111281410B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010174701.2A CN111281410B (zh) 2020-03-13 2020-03-13 一种基于造影技术的脑血管疾病成像方法及应用的造影剂

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010174701.2A CN111281410B (zh) 2020-03-13 2020-03-13 一种基于造影技术的脑血管疾病成像方法及应用的造影剂

Publications (2)

Publication Number Publication Date
CN111281410A true CN111281410A (zh) 2020-06-16
CN111281410B CN111281410B (zh) 2021-04-13

Family

ID=71024963

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010174701.2A Expired - Fee Related CN111281410B (zh) 2020-03-13 2020-03-13 一种基于造影技术的脑血管疾病成像方法及应用的造影剂

Country Status (1)

Country Link
CN (1) CN111281410B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117649371A (zh) * 2024-01-30 2024-03-05 西安交通大学医学院第一附属医院 一种脑血管介入手术模拟器图像处理方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101732733A (zh) * 2009-12-24 2010-06-16 厦门大学 一种ct成像造影剂及其制备方法
CN102905696A (zh) * 2010-02-23 2013-01-30 首尔大学校产学协力团 表面改性的氧化钽纳米颗粒及其制备方法、利用该氧化钽纳米颗粒的x射线计算机断层摄影用造影剂及高介电薄膜
WO2014162273A1 (en) * 2013-04-03 2014-10-09 Koninklijke Philips N.V. Vessel segmentation
CN106355584A (zh) * 2016-08-30 2017-01-25 上海交通大学 基于局部熵确定阈值的眼底图像微动脉瘤自动检测方法
CN109166101A (zh) * 2018-07-24 2019-01-08 博动医学影像科技(上海)有限公司 血管图像的处理方法、装置、计算机存储介质及成像设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101732733A (zh) * 2009-12-24 2010-06-16 厦门大学 一种ct成像造影剂及其制备方法
CN102905696A (zh) * 2010-02-23 2013-01-30 首尔大学校产学协力团 表面改性的氧化钽纳米颗粒及其制备方法、利用该氧化钽纳米颗粒的x射线计算机断层摄影用造影剂及高介电薄膜
WO2014162273A1 (en) * 2013-04-03 2014-10-09 Koninklijke Philips N.V. Vessel segmentation
CN106355584A (zh) * 2016-08-30 2017-01-25 上海交通大学 基于局部熵确定阈值的眼底图像微动脉瘤自动检测方法
CN109166101A (zh) * 2018-07-24 2019-01-08 博动医学影像科技(上海)有限公司 血管图像的处理方法、装置、计算机存储介质及成像设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117649371A (zh) * 2024-01-30 2024-03-05 西安交通大学医学院第一附属医院 一种脑血管介入手术模拟器图像处理方法及装置
CN117649371B (zh) * 2024-01-30 2024-04-09 西安交通大学医学院第一附属医院 一种脑血管介入手术模拟器图像处理方法及装置

Also Published As

Publication number Publication date
CN111281410B (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
Bayat et al. Background removal and vessel filtering of noncontrast ultrasound images of microvasculature
US11399779B2 (en) System-independent quantitative perfusion imaging
US20110150309A1 (en) Method and system for managing imaging data, and associated devices and compounds
CN107527341B (zh) 血管造影图像的处理方法和系统
US8452067B2 (en) Method for generating bone mask
JP4363833B2 (ja) 局所血流動態に関するインデックスを演算する方法及び装置
Longo et al. Assessment of hessian-based Frangi vesselness filter in optoacoustic imaging
WO2022105623A1 (zh) 一种基于迁移学习的颅内血管病灶识别方法
JP2019512295A (ja) 撮像方法およびシステム
Huang et al. Improved ultrafast power Doppler imaging by using spatiotemporal non-local means filtering
CN104688269B (zh) 一种时间强度曲线的呼吸运动补偿及双峰拟合的方法
Mredhula et al. An extensive review of significant researches on medical image denoising techniques
JP5289966B2 (ja) インターベンショナル手順中、シルエットレンダリング及び画像を表示するための画像処理システム及び方法
CN111281410B (zh) 一种基于造影技术的脑血管疾病成像方法及应用的造影剂
CN110101401B (zh) 一种肝脏对比剂数字减影造影方法
JP4302180B2 (ja) 局所血流動態に関するインデックスを演算する方法及び装置
JP2008100121A (ja) 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及び記憶媒体
Son et al. Image quality of CT angiography with model-based iterative reconstruction in young children with congenital heart disease: comparison with filtered back projection and adaptive statistical iterative reconstruction
CN115984270A (zh) 一种脑灌注数据分析方法、装置及存储介质
van Andel et al. Interobserver variability in the detection of cerebral venous thrombosis using CT venography with matched mask bone elimination
Hafizah et al. New technique towards operator independent kidney ultrasound scanning
Yoshida et al. Semi-automatic detection of calcified plaque in coronary CT angiograms with 320-MSCT
Huang et al. Non-local Hessian-based weighted filter for 2D noisy angiogram image enhancement with cardiovascular tree structure preservation
Rollano-Hijarrubia et al. Improving the imaging of calcifications in CT by histogram-based selective deblurring
Oezdemir et al. Three-dimensional super-resolution ultrasound imaging of chicken embryos-A validation framework for analysis of microvascular morphology

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 430000 No. 1280-1285 Jiefang Avenue, Wuhan, Hubei Province

Applicant after: Yang Zhenquan

Address before: 250000 Shandong image center, 324 jingwuweiqi Road, Shizhong District, Jinan City, Shandong Province

Applicant before: Yang Zhenquan

CB03 Change of inventor or designer information

Inventor after: Wu Weiquan

Inventor after: Li Yuange

Inventor after: Yang Zhenquan

Inventor before: Yang Zhenquan

CB03 Change of inventor or designer information
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20210326

Address after: 524001 South 57, people's Road South, Xiashan District, Zhanjiang, Guangdong

Applicant after: AFFILIATED HOSPITAL OF GUANGDONG MEDICAL University

Address before: No. 1280-1285, Jiefang Avenue, Wuhan, Hubei 430000

Applicant before: Yang Zhenquan

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210413