CN111528885B - 一种基于能谱增强ct的图像处理方法 - Google Patents
一种基于能谱增强ct的图像处理方法 Download PDFInfo
- Publication number
- CN111528885B CN111528885B CN202010295220.7A CN202010295220A CN111528885B CN 111528885 B CN111528885 B CN 111528885B CN 202010295220 A CN202010295220 A CN 202010295220A CN 111528885 B CN111528885 B CN 111528885B
- Authority
- CN
- China
- Prior art keywords
- energy spectrum
- low
- region
- patient
- energy
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 35
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- 206010061216 Infarction Diseases 0.000 claims abstract description 18
- 230000007574 infarction Effects 0.000 claims abstract description 18
- 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 16
- 239000011630 iodine Substances 0.000 claims abstract description 15
- 239000002872 contrast media Substances 0.000 claims abstract description 13
- 208000006011 Stroke Diseases 0.000 claims abstract description 9
- 238000002591 computed tomography Methods 0.000 claims description 15
- 206010053648 Vascular occlusion Diseases 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 8
- 208000021331 vascular occlusion disease Diseases 0.000 claims description 8
- 210000000988 bone and bone Anatomy 0.000 claims description 3
- 210000001175 cerebrospinal fluid Anatomy 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 208000032382 Ischaemic stroke Diseases 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 5
- 210000005013 brain tissue Anatomy 0.000 description 4
- 230000003727 cerebral blood flow Effects 0.000 description 4
- 230000010412 perfusion Effects 0.000 description 3
- 230000001154 acute effect Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004091 panning Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000002355 dual-layer Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
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/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/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/504—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 blood vessels, e.g. by angiography
-
- 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
- A61B6/5241—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于能谱增强CT的图像处理方法,利用能谱CT技术采集注射碘对比剂后患者头部的增强能谱像,并对能谱像进行计算机图像处理,从而可以更敏感的得到碘对比剂分布较低的区域,即患者缺血性卒中血管堵塞后引起的梗死区域,实现对血管闭塞后中风患者的梗死区域的精准定位。
Description
技术领域
本发明属于计算机图像技术领域,涉及一种图像处理方法,具体是一种基于能谱增强CT的图像处理方法。
背景技术
脑动脉的急性或慢性闭塞会使部分脑组织得不到充足的血液供应而梗死,从而形成缺血型卒中。缺血型卒中的急性期影像学检查中,常使用计算机断层成像(ComputerTomography)平扫(non-contrast enhancedCT)以及注射碘对比剂后的计算机断层增强成像(CTA)来进行脑组织梗死位置的诊断,其中,梗死区域一般会以低密度影呈现,即在影像上相对健康脑组织更暗。
然而,无论是CT平扫还是CTA,对梗死区域的诊断敏感性都不足够,很难界定低密度影的范围并量化其体积。因此,使用计算机断层扫描动态灌注成像(ComputerTomography Perfusion Imaging,CTP)来计算CBF(Cerebral Blood Flow)、CBV(CerebralBlood Volume)、TMAX(Time to maxenhancement)、TTP(Time to Peak)等参数图,再进一步的对这些参数图进行阈值处理,经过实验发现,可以认为CBF小于30%中位数的区域为核心梗死区,该种扫描需要128排(256层)并配有摇篮床的CT扫描仪,该种扫描仪价格昂贵,没有得到很好的普及,大部分医疗机构在急诊时,CT扫描还局限于传统的平扫像及增强像两种,无法开展CTP灌注成像。此外,这种对于梗死核心区的CBF阈值假定是基于循证医学多中心随机对照试验的统计结果,并没有直接反映脑组织的生理特性。
发明内容
本发明的目的在于提供一种基于能谱增强CT的图像处理方法,利用能谱计算机断层成像在头部的增强能谱像的处理,以定位血管闭塞后中风患者的梗死区域。
本发明的目的可以通过以下技术方案实现:
一种基于能谱增强CT的图像处理方法,具体包括以下步骤:
步骤S1,对患者注射碘对比剂进行动脉期CTA造影后,利用能谱CT技术,延迟3-7秒采集延迟期的能谱像,所述能谱像的每一个体素点均为CT值(HU)-x线能量(KeV)曲线;
步骤S2,截取能谱像任一体素点的CT值-x线能量曲线中25keV至70keV的部分,并求截取部分曲线的一阶导数曲线;
步骤S3,对该一阶导数曲线进行傅里叶变换,将频率等于0时的取值的倒数作为该体素点的低灌注参数;
步骤S4,遍历所有体素点,利用所有体素点的低灌注参数形成低灌注参数图,并与原始影像进行融合叠加;
步骤S5,去除骨骼脑脊液区域,阈值处理该低灌注参数图,得到低灌注参数高于预设阈值的区域,量化该区域的体积,得到患者的梗死区域,完成对血管闭塞后中风患者梗死区域的定位。
本发明的有益效果:本发明提供了一种基于能谱增强CT的图像处理方法,利用能谱CT技术采集注射碘对比剂后患者头部的增强能谱像,并对能谱像进行计算机图像处理,从而可以更敏感的得到碘对比剂分布较低的区域,即患者缺血性卒中血管堵塞后引起的梗死区域,实现对血管闭塞后中风患者的梗死区域的精准定位。
附图说明
下面结合附图和具体实施例对本发明作进一步详细描述。
图1是本发明的方法流程图。
图2是碘原子随光子能量增强的线衰减系数曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“开孔”、“上”、“下”、“厚度”、“顶”、“中”、“长度”、“内”、“四周”等指示方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的组件或元件必须具有特定的方位,以特定的方位构造和操作,因此不能理解为对本发明的限制。
如图1所示,本发明提供了一种基于能谱增强CT的图像处理方法,具体包括以下步骤:
步骤S1,对血管闭塞后中风患者注射碘对比剂进行动脉期CTA造影后,利用能谱CT技术,延迟3-7秒采集延迟期的能谱像。
其中,能谱CT技术是指利用双层探测器得到任意体素点的各能级的CT值,在注射碘对比剂后,使用能谱或双能量CT可以得到所谓碘图成像(iodine mapping),即碘对比剂在组织中的分布情况。
其中,能谱像的每一个体素点均为CT值(HU)-x线能量(KeV)曲线。
研究表明,碘原子随光子能量增强的线衰减系数在40keV光子能级时,有一个阶跃上升,即所谓的k-edge,如图2所示,从而,可以利用双能量CT在40KeV能量附近进行成像,捕捉到该衰减系数的上升以得到更敏感的碘对比剂的分布情况。
步骤S2,截取能谱像任一体素点的CT值-x线能量曲线中25keV至70keV的部分,并求截取部分曲线的一阶导数曲线。
步骤S3,对该一阶导数曲线进行傅里叶变换,将频率等于0时的取值的倒数作为该体素点的低灌注参数。
步骤S4,遍历所有体素点,利用所有体素点的低灌注参数形成低灌注参数图,并与原始影像进行融合叠加。
步骤S5,去除骨骼脑脊液区域,阈值处理该低灌注参数图,得到低灌注参数高于预设阈值的区域,量化该区域的体积,得到患者的梗死区域,完成对血管闭塞后中风患者梗死区域的定位。
本发明提供了一种基于能谱增强CT的图像处理方法,利用能谱CT技术采集注射碘对比剂后患者头部的增强能谱像,并对能谱像进行计算机图像处理,从而可以更敏感的得到碘对比剂分布较低的区域,即患者缺血性卒中血管堵塞后引起的梗死区域,实现对血管闭塞后中风患者的梗死区域的精准定位。
在本说明书的描述中,参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上内容仅仅是对本发明结构所作的举例和说明,所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (1)
1.一种基于能谱增强CT的图像处理方法,其特征在于,具体包括以下步骤:
步骤S1,对血管闭塞后中风患者注射碘对比剂进行动脉期CTA造影后,利用能谱CT技术,延迟3-7秒采集延迟期的能谱像,所述能谱像的每一个体素点均为CT值(HU)-x线能量(KeV)曲线;
步骤S2,截取能谱像任一体素点的CT值-x线能量曲线中25keV至70keV的部分,并求截取部分曲线的一阶导数曲线;
步骤S3,对该一阶导数曲线进行傅里叶变换,将频率等于0时的取值的倒数作为该体素点的低灌注参数;
步骤S4,遍历所有体素点,利用所有体素点的低灌注参数形成低灌注参数图,并与原始影像进行融合叠加;
步骤S5,去除骨骼脑脊液区域,阈值处理该低灌注参数图,得到低灌注参数高于预设阈值的区域,量化该区域的体积,得到患者的梗死区域,完成对血管闭塞后中风患者梗死区域的定位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010295220.7A CN111528885B (zh) | 2020-04-15 | 2020-04-15 | 一种基于能谱增强ct的图像处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010295220.7A CN111528885B (zh) | 2020-04-15 | 2020-04-15 | 一种基于能谱增强ct的图像处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111528885A CN111528885A (zh) | 2020-08-14 |
CN111528885B true CN111528885B (zh) | 2023-09-05 |
Family
ID=71970831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010295220.7A Active CN111528885B (zh) | 2020-04-15 | 2020-04-15 | 一种基于能谱增强ct的图像处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111528885B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010045478A2 (en) * | 2008-10-15 | 2010-04-22 | Brigham And Women's Hospital, Inc. | System and method for differential perfusion analysis of tissues with multiple vascular supplies |
CN101779963A (zh) * | 2010-02-23 | 2010-07-21 | 上海市第八人民医院 | 减影计算机体层灌注功能成像方法 |
CN105997128A (zh) * | 2016-08-03 | 2016-10-12 | 上海联影医疗科技有限公司 | 利用灌注成像识别病灶的方法及*** |
CN108693491A (zh) * | 2017-04-07 | 2018-10-23 | 康奈尔大学 | 稳健的定量磁化率成像***和方法 |
CN109431531A (zh) * | 2018-12-25 | 2019-03-08 | 上海联影医疗科技有限公司 | 基于灌注成像的血管分割方法及装置和计算机装置 |
CN110223285A (zh) * | 2019-06-13 | 2019-09-10 | 中南大学湘雅二医院 | 基于动物迁移和神经网络模型的影像结果预测方法及*** |
CN110288587A (zh) * | 2019-06-28 | 2019-09-27 | 重庆同仁至诚智慧医疗科技股份有限公司 | 一种缺血性脑卒中磁共振影像的病灶识别方法 |
CN110934606A (zh) * | 2019-10-31 | 2020-03-31 | 上海杏脉信息科技有限公司 | 脑卒中早期平扫ct图像评估***及评估方法、可读存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6792302B2 (en) * | 2001-02-21 | 2004-09-14 | Universite De Lausanne | Method and apparatus for determining treatment for stroke |
US7218702B2 (en) * | 2004-05-10 | 2007-05-15 | Wisconsin Alumni Research Foundation | X-ray system for use in image guided procedures |
-
2020
- 2020-04-15 CN CN202010295220.7A patent/CN111528885B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010045478A2 (en) * | 2008-10-15 | 2010-04-22 | Brigham And Women's Hospital, Inc. | System and method for differential perfusion analysis of tissues with multiple vascular supplies |
CN101779963A (zh) * | 2010-02-23 | 2010-07-21 | 上海市第八人民医院 | 减影计算机体层灌注功能成像方法 |
CN105997128A (zh) * | 2016-08-03 | 2016-10-12 | 上海联影医疗科技有限公司 | 利用灌注成像识别病灶的方法及*** |
CN108693491A (zh) * | 2017-04-07 | 2018-10-23 | 康奈尔大学 | 稳健的定量磁化率成像***和方法 |
CN109431531A (zh) * | 2018-12-25 | 2019-03-08 | 上海联影医疗科技有限公司 | 基于灌注成像的血管分割方法及装置和计算机装置 |
CN110223285A (zh) * | 2019-06-13 | 2019-09-10 | 中南大学湘雅二医院 | 基于动物迁移和神经网络模型的影像结果预测方法及*** |
CN110288587A (zh) * | 2019-06-28 | 2019-09-27 | 重庆同仁至诚智慧医疗科技股份有限公司 | 一种缺血性脑卒中磁共振影像的病灶识别方法 |
CN110934606A (zh) * | 2019-10-31 | 2020-03-31 | 上海杏脉信息科技有限公司 | 脑卒中早期平扫ct图像评估***及评估方法、可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111528885A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Scholtz et al. | Advances in cardiac CT contrast injection and acquisition protocols | |
US7853309B2 (en) | X-ray CT apparatus and myocardial perfusion image generating system | |
US9949650B2 (en) | Fractional flow reserve estimation | |
Fleischmann | How to design injection protocols for multiple detector-row CT angiography (MDCTA) | |
Schaible et al. | Quantitative perfusion analysis of hepatocellular carcinoma using dynamic contrast enhanced ultrasound (CEUS) to determine tumor microvascularization | |
Lee et al. | Distal lower extremity arteries: evaluation with two-dimensional MR digital subtraction angiography. | |
JP4857418B2 (ja) | 検体のモーションアーチファクト補正の方法およびシステム | |
Hur et al. | Dual-enhancement cardiac computed tomography for assessing left atrial thrombus and pulmonary veins before radiofrequency catheter ablation for atrial fibrillation | |
RU2604941C2 (ru) | Частота дискретизации динамической регистрации для компьютерной томографической визуализации перфузии (стр) | |
KR20180059510A (ko) | 적응식 스캔 제어를 위한 방법 및 시스템 | |
Goo | Semiautomatic three-dimensional threshold-based cardiac computed tomography ventricular volumetry in repaired tetralogy of Fallot: comparison with cardiac magnetic resonance imaging | |
Mendoza et al. | Viability imaging by cardiac computed tomography | |
CN106469449B (zh) | 一种医学图像中病灶的显示方法及装置 | |
JP2002233525A (ja) | 像形成法を用いた生体の診察方法 | |
JP2012512729A (ja) | 急性脳卒中の組織状態をマッピングする方法及びシステム | |
CN113409328B (zh) | Ct图像的肺动静脉分割方法、装置、介质及设备 | |
Yang et al. | Analysis of elasticity characteristics of ascending aorta, descending aorta and pulmonary artery using 640 slice-volume CT | |
Tridandapani et al. | Coronary computed tomographic angiography: a review of the techniques, protocols, pitfalls, and radiation dose | |
Ippolito et al. | Diagnostic accuracy of 256-row multidetector CT coronary angiography with prospective ECG-gating combined with fourth-generation iterative reconstruction algorithm in the assessment of coronary artery bypass: evaluation of dose reduction and image quality | |
Nasis et al. | Assessment of left ventricular volumes, ejection fraction and regional wall motion with retrospective electrocardiogram triggered 320-detector computed tomography: a comparison with 2D-echocardiography | |
Lu et al. | Evaluation of coronary artery bypass graft patency using three-dimensional reconstruction and flow study on electron beam tomography | |
CN111528885B (zh) | 一种基于能谱增强ct的图像处理方法 | |
Bardo et al. | Cardiac multidetector computed tomography: basic physics of image acquisition and clinical applications | |
Lazoura et al. | Myocardial bridging evaluated with 128-multi detector computed tomography coronary angiography | |
CN108814633A (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 |