CN111281410B - 一种基于造影技术的脑血管疾病成像方法及应用的造影剂 - Google Patents
一种基于造影技术的脑血管疾病成像方法及应用的造影剂 Download PDFInfo
- Publication number
- CN111281410B CN111281410B CN202010174701.2A CN202010174701A CN111281410B CN 111281410 B CN111281410 B CN 111281410B CN 202010174701 A CN202010174701 A CN 202010174701A CN 111281410 B CN111281410 B CN 111281410B
- Authority
- CN
- China
- Prior art keywords
- contrast agent
- polylactic acid
- cerebrovascular
- image
- dimensional
- 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.)
- Expired - Fee Related
Links
- 239000002872 contrast media Substances 0.000 title claims abstract description 72
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 208000026106 cerebrovascular disease Diseases 0.000 title claims abstract description 18
- 238000002601 radiography Methods 0.000 title description 5
- 229920000747 poly(lactic acid) Polymers 0.000 claims abstract description 52
- 239000004626 polylactic acid Substances 0.000 claims abstract description 52
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 47
- 230000002490 cerebral effect Effects 0.000 claims abstract description 31
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims abstract description 24
- 239000010931 gold Substances 0.000 claims abstract description 24
- 229910052737 gold Inorganic materials 0.000 claims abstract description 24
- 239000002105 nanoparticle Substances 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 20
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 claims abstract description 17
- 229910052740 iodine Inorganic materials 0.000 claims abstract description 17
- 239000011630 iodine Substances 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 17
- 229920003171 Poly (ethylene oxide) Polymers 0.000 claims abstract description 16
- 125000000956 methoxy group Chemical group [H]C([H])([H])O* 0.000 claims abstract description 16
- -1 polyoxyethylene Polymers 0.000 claims abstract description 16
- 125000001436 propyl group Chemical group [H]C([*])([H])C([H])([H])C([H])([H])[H] 0.000 claims abstract description 16
- 239000007864 aqueous solution Substances 0.000 claims abstract description 12
- 230000011218 segmentation Effects 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 10
- KPUWHANPEXNPJT-UHFFFAOYSA-N disiloxane Chemical group [SiH3]O[SiH3] KPUWHANPEXNPJT-UHFFFAOYSA-N 0.000 claims abstract description 8
- 239000003607 modifier Substances 0.000 claims abstract description 5
- 239000012071 phase Substances 0.000 claims description 31
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 26
- 239000000243 solution Substances 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 6
- 238000002360 preparation method Methods 0.000 claims description 5
- 238000003325 tomography Methods 0.000 claims description 5
- HEDRZPFGACZZDS-UHFFFAOYSA-N Chloroform Chemical compound ClC(Cl)Cl HEDRZPFGACZZDS-UHFFFAOYSA-N 0.000 claims description 4
- 229960001701 chloroform Drugs 0.000 claims description 4
- 238000009210 therapy by ultrasound Methods 0.000 claims description 4
- 210000003462 vein Anatomy 0.000 claims description 4
- 239000008215 water for injection Substances 0.000 claims description 4
- 239000008346 aqueous phase Substances 0.000 claims description 3
- 238000004132 cross linking Methods 0.000 claims description 2
- 238000001723 curing Methods 0.000 claims 1
- 210000004872 soft tissue Anatomy 0.000 abstract description 17
- 238000003708 edge detection Methods 0.000 abstract description 14
- 230000000694 effects Effects 0.000 abstract description 9
- 238000001914 filtration Methods 0.000 abstract description 8
- 238000010253 intravenous injection Methods 0.000 abstract description 2
- 238000002591 computed tomography Methods 0.000 description 16
- 230000006872 improvement Effects 0.000 description 15
- 238000001514 detection method Methods 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 239000002086 nanomaterial Substances 0.000 description 3
- 230000002035 prolonged effect Effects 0.000 description 3
- 201000008450 Intracranial aneurysm Diseases 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- NTHXOOBQLCIOLC-UHFFFAOYSA-N iohexol Chemical compound OCC(O)CN(C(=O)C)C1=C(I)C(C(=O)NCC(O)CO)=C(I)C(C(=O)NCC(O)CO)=C1I NTHXOOBQLCIOLC-UHFFFAOYSA-N 0.000 description 2
- 229960001025 iohexol Drugs 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 208000019553 vascular disease Diseases 0.000 description 2
- 208000022211 Arteriovenous Malformations Diseases 0.000 description 1
- 206010003226 Arteriovenous fistula Diseases 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 206010065559 Cerebral arteriosclerosis Diseases 0.000 description 1
- 208000034710 Cerebral arteriovenous malformation Diseases 0.000 description 1
- 208000002263 Intracranial Arteriovenous Malformations Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 208000009443 Vascular Malformations Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000005744 arteriovenous malformation Effects 0.000 description 1
- 201000000034 arteriovenous malformations of the brain Diseases 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 206010008087 cerebral arteritis Diseases 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 201000002676 cerebral atherosclerosis Diseases 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000003517 fume Substances 0.000 description 1
- 230000002008 hemorrhagic effect Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000007917 intracranial administration Methods 0.000 description 1
- 201000005851 intracranial arteriosclerosis Diseases 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000002198 surface plasmon resonance spectroscopy Methods 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/501—Apparatus 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61K—PREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
- A61K49/00—Preparations for testing in vivo
- A61K49/04—X-ray contrast preparations
- A61K49/0433—X-ray contrast preparations containing an organic halogenated X-ray contrast-enhancing agent
- A61K49/0438—Organic X-ray contrast-enhancing agent comprising an iodinated group or an iodine atom, e.g. iopamidol
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood 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和血管多尺度因子σ，设计边缘检测算子，对血管进行多尺度边缘检测，滤除没有边缘结构的含有造影剂的软组织噪声。
作为本发明的进一步改进，所述边缘检测算子的公式如下：
作为本发明的进一步改进，在步骤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)所示的分段线性函数对脑血管横截面二维图像进行分段线性变换：
如对原始输入的脑血管横截面二维图像I(x,y)进行变换时，a、b代表需要调整的灰度范围区间的端点值，c、d代表调整好的灰度范围区间的端点，决定对脑血管横截面二维图像中某个灰度区间的像素进行拉伸或者是压缩，是灰度变换的重点。采用此种变换，能够对脑血管横截面二维图像进行灰度改善，提升血管像素的灰度值，压缩非血管像素的灰度值，从而有效的抑制了非血管像素的灰度值，保留了过渡区间的灰度值，并一定程度上增强了血管目标的灰度值。
当造影剂在血管中流动时，造影剂不可避免的会渗出血管到血管周围的软组织，如造影剂渗出达到一定程度时，血管周围的组织也会在CT造影过程中被保留下来，成为影响血管影像的主要干扰因素，尤其在相比较来说比较细的血管区域这种干扰更加明显。渗透到软组织中的造影剂会影响视觉效果，但是，含有造影剂的软组织是扩散式分布的，没有固定的区域边缘，血管却是具有明确的边缘结构。根据这个特性，本发明设计一个边缘检测算子专门针对血管边缘结构进行检测，滤除没有边缘结构的含有造影剂的软组织区域。
因此，作为本发明的进一步改进，在步骤S2中，所述图像增强处理是利用梯度算子|▽t|、曲率检测因子λ2和血管多尺度因子σ，设计边缘检测算子，对血管进行多尺度边缘检测，滤除没有边缘结构的含有造影剂的软组织噪声。所述边缘检测算子的公式如式(2)所示：
采用式(2)进行边缘检测增强时，由于血管在边缘区域会产生灰度的陡变，陡变度可以利用梯度算子进行检测，从而得到边缘明确的血管，很大程度的平滑血管边缘，凸显血管结构。当检测到血管内的像素时，曲率检测因子λ2的值相对较大，导致梯度T较小，从而不影响血管的增强；当检测到血管外的像素点时，由于软组织没有固定的边缘和走向，所以不存在曲率检测因子λ2，从而通过梯度T过滤这些像素点，实现了软组织噪声的去除。
作为本发明的进一步改进，在步骤S3中，所述血管分割采用基于二阶微分的拉普拉斯算子，如式(3)所示，对脑血管横截面预处理图像进行分割：
优选地，在采用基于二阶微分的拉普拉斯算子进行血管分割时，先采用高斯滤波器对脑血管横截面预处理图像进行平滑处理，如式(4)所示：
式中，σ为图像的标准差，在实际平滑过程中，σ过小，则偏离中心的所有像素权重会非常的小，当σ小到一定程度时，滤波操作则会退化为对图像进行点运算，无法起到平滑噪声的作；但如果σ过大，邻域较小时，在邻域内的高斯模板将会退化为均值模板，会造成边缘模糊的效果，所以本发明σ值优选为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所述每个层面的脑血管横截面二维图像进行分段线性变换：
式中，a、b代表需要调整的灰度范围区间的端点值，c、d代表调整好的灰度范围区间的端点；
采用采用下式所示的边缘检测算子对进行边缘检测增强处理，得到每个层面的脑血管横截面预处理图像：
式中，|▽t|为梯度算子、λ2为曲率检测因子，σ为血管多尺度因子；
S3.采用下式所示的高斯滤波函数对步骤S2所述每个层面的脑血管横截面预处理图像进行平滑处理：
式中，σ为图像的标准差，取值为0.7；
然后采用下式所示的基于二阶微分的拉普拉斯算子对平滑后的每个层面的脑血管横截面预处理图像进行血管分割，得到每个层面的脑血管图像：
S4.采用光线投射法将步骤S3所述每个层面的脑血管图像进行三维重建，得到脑血管三维图像。
综上所述，本发明采用自制的造影剂，将金纳米粒子聚乳酸溶液作为油相，非离子型碘造影剂水溶液作为水相，得到金纳米粒子和聚乳酸包裹非离子型碘造影剂的聚乳酸微泡，再以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷修饰聚乳酸微泡表面，能够提高造影剂在血液中的生物相容性和稳定性，并能够增强造影，延长体内血液循环时间。以此造影剂进行CT造影成像，对采集的二维图像进行分阶段线性灰度变换和边缘检测算子滤除软组织噪声处理，再采用基于二阶微分的拉普拉斯算子进行血管分割，最后将二维血管图像进行三维重建得到脑血管三维图像。能够有效的抑制非血管像素的灰度值，提升血管像素的灰度值，并有效滤除血管周围的软组织区域，从而得到高精度的脑血管二维图像。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
Claims (3)
1.一种基于造影技术的脑血管疾病成像方法中应用的造影剂，其特征在于，所述造影剂为以含有金纳米粒子的聚乳酸为油相，以含有非离子型碘造影剂的水溶液为水相，以2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷为改性剂，制备成的硅氧烷修饰的聚乳酸微泡；所述造影剂的制备方法包括以下步骤：
(1)将金纳米粒子加入到含有聚乳酸的三氯甲烷溶液中，作为油相；
(2)将非离子型碘造影剂加入到注射用水，得到质量分数为70％-90％的造影剂水溶液，作为水相；
(3)将步骤(2)所述水相加入到步骤(1)所述油相中，超声处理得到聚乳酸微泡；
(4)将步骤(3)所述聚乳酸微泡加入到浓度为4％-8％(w/v)的2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷的乙醇溶液中，使溶液中的乙醇完全挥发掉，此过程中，2-[甲氧基(聚氧乙烯)丙基]-三甲氧基硅烷在所述聚乳酸微泡表面固化交联，得到硅氧烷修饰的聚乳酸微泡；
所述造影剂的使用方法为：
S1.从静脉注射所述造影剂，待所述造影剂到达脑血管区域后，对脑血管区域进行X射线断层扫描，得到不同层面的脑血管横截面二维图像；
S2.对步骤S1中每个所述层面的脑血管横截面二维图像依次进行灰度变换和图像增强处理，得到每个层面的脑血管横截面预处理图像；
S3.对步骤S2所述每个层面的脑血管横截面预处理图像进行血管分割，得到每个层面的脑血管图像；
S4.将步骤S3所述每个层面的脑血管图像进行三维重建，得到脑血管三维图像。
2.根据权利要求1所述的脑血管疾病成像方法中应用的造影剂，其特征在于，在步骤(1)中，所述油相中金纳米粒子的浓度为2-6mg/ml，聚乳酸的浓度为30-60mg/ml。
3.根据权利要求1所述的脑血管疾病成像方法中应用的造影剂，其特征在于，在步骤(3)中，所述聚乳酸微泡和乙醇溶液的体积比为1:2.5-1:4。
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 CN111281410A (zh) | 2020-06-16 |
| CN111281410B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117649371B (zh) * | 2024-01-30 | 2024-04-09 | 西安交通大学医学院第一附属医院 | 一种脑血管介入手术模拟器图像处理方法及装置 |
Citations (5)
| 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 | 博动医学影像科技（上海）有限公司 | 血管图像的处理方法、装置、计算机存储介质及成像设备 |
-
2020
- 2020-03-13 CN CN202010174701.2A patent/CN111281410B/zh not_active Expired - Fee Related
Patent Citations (5)
| 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 | 博动医学影像科技（上海）有限公司 | 血管图像的处理方法、装置、计算机存储介质及成像设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111281410A (zh) | 2020-06-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bayat et al. | Background removal and vessel filtering of noncontrast ultrasound images of microvasculature | |
| US20110150309A1 (en) | Method and system for managing imaging data, and associated devices and compounds | |
| US8452067B2 (en) | Method for generating bone mask | |
| US20060025680A1 (en) | Method and apparatus for contrast enhanced medical imaging | |
| JP4363833B2 (ja) | 局所血流動態に関するインデックスを演算する方法及び装置 | |
| Chen et al. | Radiation dose and image quality of 70 kVp cerebral CT angiography with optimized sinogram-affirmed iterative reconstruction: comparison with 120 kVp cerebral CT angiography | |
| CN107527341A (zh) | 血管造影图像的处理方法和系统 | |
| JP2019512295A (ja) | 撮像方法およびシステム | |
| Longo et al. | Assessment of hessian-based Frangi vesselness filter in optoacoustic imaging | |
| CN109350089B (zh) | 一种基于ct图像的甲状腺区域自动分割方法 | |
| CN111281410B (zh) | 一种基于造影技术的脑血管疾病成像方法及应用的造影剂 | |
| JP5289966B2 (ja) | インターベンショナル手順中、シルエットレンダリング及び画像を表示するための画像処理システム及び方法 | |
| JP4302180B2 (ja) | 局所血流動態に関するインデックスを演算する方法及び装置 | |
| JP2008100121A (ja) | 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及び記憶媒体 | |
| 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 | |
| Rollano-Hijarrubia et al. | Improving the imaging of calcifications in CT by histogram-based selective deblurring | |
| Huang et al. | Non-local Hessian-based weighted filter for 2D noisy angiogram image enhancement with cardiovascular tree structure preservation | |
| Oezdemir et al. | Three-dimensional super-resolution ultrasound imaging of chicken embryos-A validation framework for analysis of microvascular morphology | |
| ÖZKAN | 3D Lung vessel segmentation in computed tomography angiography images | |
| Kierski et al. | Acoustic Molecular Imaging Beyond the Diffraction Limit In Vivo | |
| CN113052805A (zh) | X射线血管造影序列图像的静态背景分离方法 | |
| Lin et al. | Clinical photoacoustic computed tomography of breast cancer treated with neoadjuvant chemotherapy | |
| Hong | Care of Maintenance Hemodialysis Patients through Intelligent Algorithm-Based Low-Dose Digital Subtraction Angiography |
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 |
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 |
|
| CB02 | Change of applicant information | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wu Weiquan Inventor after: Li Yuange Inventor after: Yang Zhenquan Inventor before: Yang Zhenquan |
|
| 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 |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210413 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |