CN103700083B - 基于自适应模糊神经网络模型ct灌注图像智能融合方法 - Google Patents
基于自适应模糊神经网络模型ct灌注图像智能融合方法 Download PDFInfo
- Publication number
- CN103700083B CN103700083B CN201310733460.0A CN201310733460A CN103700083B CN 103700083 B CN103700083 B CN 103700083B CN 201310733460 A CN201310733460 A CN 201310733460A CN 103700083 B CN103700083 B CN 103700083B
- Authority
- CN
- China
- Prior art keywords
- image
- blood flow
- cerebral blood
- characteristic pattern
- exceptional value
- 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
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于自适应模糊神经网络模型CT灌注图像智能融合方法,涉及人工智能、模式识别与医疗工程的交叉领域,可应用于缺血性脑部的诊断。该方法在CT灌注图像中只提取与缺血性脑血管病相关的信息—脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图,经合理的动态模糊神经网络模型训练学习,通过恰当的图像融合方法处理,得到CT解剖图的融合图像,该融合图像可应用于缺血性脑部的诊断,实现了头颅灌注多参数图像综合分析的智能化和自动化。
Description
技术领域
本发明涉及人工智能、模式识别与医疗工程的交叉领域,更具体涉及CT灌注图像的智能融合方法,可应用于缺血性脑部的诊断。
背景技术
脑血管疾病目前是目前全球范围内致死率、致残率最高的疾病之一,严重危害人类健康,并给社会和患者家庭带来沉重的负担。其中,缺血性脑血管病又是最常见的脑血管疾病类型,占所有脑血管病的65%-80%。不同病因引起的缺血性脑血管病脑血液动力学障碍均可能导致脑组织缺血、梗死而引起患者神经功能障碍,严重的甚至导致患者死亡。因此,进行缺血性脑血管病的研究,发明切实有效地疾病诊断治疗方法,不仅仅有重要的科研价值,同时也具备重大的实际应用价值。
随着计算机成像技术和计算机水平的不断提高,利用模式识别技术进行医学图像的研究逐渐成为图像处理研究的热点。随着研究的不断深入,目前已经具备一定的理论基础和实践经验,也取得了一系列成果。
在治疗缺血性脑病中,计算机体层摄影灌注成像作为一种用于了解人体组织血流灌注状态的功能成像手段现已列入诊断缺血性脑血管病的常规检查方案。据文献报道,CT灌注成像能最早显示出现症状30分钟以内的脑缺血病灶,对诊断早期脑梗死的敏感性和特异性分别为93%、98%。因此,CT灌注成像较之常规解剖学成像对于早期发现缺血性脑病变、对超急性期脑梗死的诊断具有更加重要的临床应用价值。目前,计算机体层摄影和磁共振脑灌注成像是一种用于了解脑组织血流灌注状态的功能成像手段,现已列入诊断缺血性脑血管病的常规检查方法。不同于常规解剖学成像的是,该技术在静脉注射对比剂的同时对选定的感兴趣层面进行连续动态扫描以获得所选层面内各像素的时间-密度曲线,根据指示剂稀释原理和中心法则并通过复杂的数学模型计算得到脑血流量、脑血流容积、对比剂平均通过时间和峰值时间等四个灌注参数(例如机器型号:Toshiba AquilionONE 320排)。对脑组织的灌注状态评价依赖于放射科医生通过对上述血流动力学灌注参数图像上的表现来判断脑组织缺血情况,并需结合解剖学成像所获取的信息来分析脑组织缺血的病变区域和病变程度信息。
神经网络模型是模拟生物神经系统原理而构建的一种新型智能化信息处理系统,简称神经网络。神经网络模型具有很强的自我组织、自我适应的学习训练能力、同时具备高度的鲁棒性和容错能力,被广泛应用于图像处理、信息处理、模式识别、数据挖掘等各个研究领域。本发明中神经网络模型(李清梦等神经网络技术及其在医学图像处理中的应用.中国医学影像学技术,2011,27(6):1291-1294.)构建缺血性脑血管病CT灌注成像参数分析模型,并应用图像融合技术(马国锐等,基于融合和广义高斯模型的遥感影像变化检测,遥感学报,2006,10(6):837-853)将灌注成像诊断信息直接与解剖学图像进行融合实现诊断结果可视化,进一步促使缺血性脑血管病CT灌注成像的诊断过程准确化、智能化。
发明内容
本发明的目的是,提供一种基于自适应模糊神经网络模型CT灌注图像智能融合方法。该方法基于动态模糊神经网络模型以及图像智能融合方法,通过对CT灌注图像中的脑血流量、脑血流容积、对比剂平均通过时间和峰值时间四个异常值图像的融合,得到CT解剖图的融合图像。本发明的优点是:实现了头颅灌注多参数CT图像综合分析的智能化和自动化。
为了达到上述目的,本发明采用如下技术方案。一种基于动态模糊神经网络模型的CT灌注医疗图像智能融合方法,该方法包含下列步骤:
a、在每张CT灌注图像中分别提取脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图;
b、利用神经网络模型分别对脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图样本进行训练,分别获得各CT灌注图像上述四个特征图的异常值图像;
c、采用图像融合方法,将各CT灌注图像上述四个特征图的异常值图像融合到CT解剖图上,得到CT解剖图的融合图像。
本发明基于动态模糊神经网络模型,采用该模型可以处理多类图像和自我归一化更强的能力,训练中在进行融合过程迭代运算时易收敛,对噪声并不敏感等优势,可有效地计算脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的异常值图像。通过分析脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的图像特征关系与CT解剖图的纹理的空间关系,采用图像融合方法中基于代数加权的乘积变化对其进行融合处理,得到CT解剖图的融合图像。
本发明的优点是:在CT灌注图像中只提取与缺血性脑血管病相关的信息—脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图,经合理模型的训练学习,通过恰当的图像融合方法处理,得到CT解剖图的融合图像,该融合图像可应用于缺血性脑部的诊断,实现了头颅灌注多参数图像综合分析的智能化和自动化。
附图说明
图1是在CT灌注图像中提取的脑血流量特征图。
图2是在CT灌注图像中提取的脑血流容积特征图。
图3是在CT灌注图像中提取的对比剂平均通过时间特征图。
图4是在CT灌注图像中提取的峰值时间特征图。
图5是在CT解剖图的融合图像。
具体实施方式
以下结合附图,对本发明作进一步的说明。
实施例一
一种基于自适应模糊神经网络模型CT灌注图像智能融合方法,该方法包含下列步骤:
a、在每张CT灌注图像中分别提取脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图样本;
b、利用神经网络模型分别对脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图样本进行训练,分别获得各CT灌注图像上述四个特征图的异常值图像;
c、采用图像融合方法,将各CT灌注图像上述四个特征图的异常值图像融合到CT解剖图上,得到CT解剖图的融合图像。
上述步骤b中所述的神经网络模型采用下面公式:
其中:当对脑血流量进行训练时,为脑血流量特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是脑血流量特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n纵坐标的最大值;Ii(k)是脑血流量特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m横坐标的最大值;当对脑血流容积进行训练时,为脑血流容积特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是脑血流容积特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n纵坐标的最大值;Ii(k)是脑血流容积特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m横坐标的最大值;当对对比剂平均通过时间进行训练时,为对比剂平均通过时间特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是对比剂平均通过时间特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n纵坐标的最大值;Ii(k)是对比剂平均通过时间特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m横坐标的最大值;当对峰值时间进行训练时,为峰值时间特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是峰值时间特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n纵坐标的最大值;Ii(k)是峰值时间特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m横坐标的最大值。
上述步骤b中所述的训练方法采用神经网络模型自带的学习训练方法。
上述步骤c中所述的图像融合方法采用基于代数加权的乘积变化图像融合方法,公式如下:
F(i,j)=A(i,j)×B(i,j)/maxB(i,j)
A(i,j)=f1(i,j)×f2(i,j)×f3(i,j)×f4(i,j)/4
其中:F(i,j)为将脑血流量、脑血流容积、对比剂平均通过时间和峰值时间特征图的异常值图像融合到CT解剖图上,得到CT解剖图的融合图像;B(i,j)为CT解剖图,maxB(i,j)为CT解剖图上像素的最大值,A(i,j)为脑血流量、脑血流容积、对比剂平均通过时间和峰值时间四个特征图的异常值图像均值,f1(i,j),f2(i,j),f3(i,j)和f4(i,j)分别为脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图的异常值图像。
实施例二
一种基于自适应模糊神经网络模型CT灌注图像智能融合方法,该方法包含下列步骤:
a、在每张CT灌注图像中分别提取脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图样本。图1为脑血流量特征图样本,图2为脑血流容积特征图样本,图3为对比剂平均通过时间特征图样本,图4为峰值时间特征图样本。
b、利用神经网络模型分别对脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图样本进行训练,分别获得各CT灌注图像上述四个特征图的异常值图像;
c、采用图像融合方法,将各CT灌注图像上述四个特征图的异常值图像融合到CT解剖图上,得到CT解剖图的融合图像。图5即为上述四幅特征图样本,经神经网络模型分别训练后,得到的四个特征图的异常值图像通过融合得到的CT解剖图的融合图像。
Claims (1)
1.基于自适应模糊神经网络模型CT灌注图像智能融合方法,其特征在于,该方法包含下列步骤:
a、在每张CT灌注图像中分别提取脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的四个特征图样本;
b、利用神经网络模型分别对脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的四个特征图样本进行训练,神经网络模型采用如下公式:
其中:当对脑血流量进行训练时,为脑血流量特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是脑血流量特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n为纵坐标的最大值;Ii(k)是脑血流量特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m为横坐标的最大值;当对脑血流容积进行训练时,为脑血流容积特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是脑血流容积特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n为纵坐标的最大值;Ii(k)是脑血流容积特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m为横坐标的最大值;当对对比剂平均通过时间进行训练时,为对比剂平均通过时间特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是对比剂平均通过时间特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n纵坐标的最大值;Ii(k)是对比剂平均通过时间特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m为横坐标的最大值;当对峰值时间进行训练时,为峰值时间特征图的异常值图像面积与CT灌注图像面积的比值;Xj(k)是峰值时间特征图的异常值图像纵坐标为j、衰减系数为k时的误差参数,其中:j=1,2,…,n,n为纵坐标的最大值;Ii(k)是峰值时间特征图的异常值图像横坐标为i、衰减系数为k时的收敛参数,其中:i=1,2,…,m,m为横坐标的最大值;
采用神经网络模型自带的学习训练方法进行训练,分别获得各CT灌注图像上述四个特征图的异常值图像;
c、采用图像融合方法,将各CT灌注图像上述四个特征图的异常值图像融合到CT解剖图上,图像融合方法采用基于代数加权的乘积变化图像融合方法,公式如下:
F(i,j)=A(i,j)×B(i,j)/max B(i,j)
A(i,j)=f1(i,j)×f2(i,j)×f3(i,j)×f4(i,j)/4
其中:F(i,j)为将脑血流量、脑血流容积、对比剂平均通过时间和峰值时间特征图的异常值图像融合到CT解剖图上,得到CT解剖图的融合图像;B(i,j)为CT解剖图,max B(i,j)为CT解剖图上像素的最大值,A(i,j)为脑血流量、脑血流容积、对比剂平均通过时间和峰值时间四个特征图的异常值图像均值,f1(i,j),f2(i,j),f3(i,j)和f4(i,j)分别为脑血流量、脑血流容积、对比剂平均通过时间和峰值时间的特征图的异常值图像;得到CT解剖图的融合图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310733460.0A CN103700083B (zh) | 2013-12-26 | 2013-12-26 | 基于自适应模糊神经网络模型ct灌注图像智能融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310733460.0A CN103700083B (zh) | 2013-12-26 | 2013-12-26 | 基于自适应模糊神经网络模型ct灌注图像智能融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103700083A CN103700083A (zh) | 2014-04-02 |
CN103700083B true CN103700083B (zh) | 2017-02-22 |
Family
ID=50361602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310733460.0A Active CN103700083B (zh) | 2013-12-26 | 2013-12-26 | 基于自适应模糊神经网络模型ct灌注图像智能融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103700083B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2698997C1 (ru) * | 2016-09-06 | 2019-09-02 | Электа, Инк. | Нейронная сеть для генерации синтетических медицинских изображений |
CN108606806A (zh) * | 2016-12-09 | 2018-10-02 | 上海西门子医疗器械有限公司 | 确定扫描时延的方法和装置、造影剂诊断扫描方法和装置 |
CN109272024B (zh) * | 2018-08-29 | 2021-08-20 | 昆明理工大学 | 一种基于卷积神经网络的图像融合方法 |
CN109523500B (zh) * | 2018-11-19 | 2020-07-31 | 复旦大学附属华山医院 | 颅内外血流重建术的搭桥血管确定方法及系统 |
CN115040100B (zh) * | 2022-06-14 | 2023-10-27 | 安影科技(北京)有限公司 | 一种视神经血流灌注数值快速采集方法 |
CN116630247B (zh) * | 2023-05-06 | 2023-10-20 | 河北省儿童医院(河北省第五人民医院、河北省儿科研究所) | 脑血流图像处理方法及装置、脑血流监测系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186895A (zh) * | 2013-04-15 | 2013-07-03 | 武汉大学 | 一种基于神经网络模型的ct灌注医疗图像智能融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7998128B2 (en) * | 2006-05-15 | 2011-08-16 | Therataxis, LLC. | Active delivery and flow redirection: novel devices and method of delivery of materials to patients |
-
2013
- 2013-12-26 CN CN201310733460.0A patent/CN103700083B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186895A (zh) * | 2013-04-15 | 2013-07-03 | 武汉大学 | 一种基于神经网络模型的ct灌注医疗图像智能融合方法 |
Non-Patent Citations (1)
Title |
---|
基于神经网络融合模型的头颅灌注图像分析;何洪林等;《医疗卫生装备》;20130131;第34卷(第1期);第19-21,31页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103700083A (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103700083B (zh) | 基于自适应模糊神经网络模型ct灌注图像智能融合方法 | |
Gao et al. | Computer vision in healthcare applications | |
Tian et al. | Multi-path convolutional neural network in fundus segmentation of blood vessels | |
Wu et al. | Gamma challenge: glaucoma grading from multi-modality images | |
Tsao et al. | Stereopsis activates V3A and caudal intraparietal areas in macaques and humans | |
CN109493954A (zh) | 一种基于类别判别定位的sd-oct图像视网膜病变检测系统 | |
CN105913432B (zh) | 基于ct序列图像的主动脉提取方法及装置 | |
CN109325942A (zh) | 基于全卷积神经网络的眼底图像结构分割方法 | |
Pruett Jr et al. | Accurate age classification of 6 and 12 month-old infants based on resting-state functional connectivity magnetic resonance imaging data | |
CN109003299A (zh) | 一种基于深度学习的计算脑出血量的方法 | |
CN111602173B (zh) | 脑部断层扫描数据分析方法 | |
CN104809480A (zh) | 一种基于分类回归树和AdaBoost的眼底图像视网膜血管分割方法 | |
CN110211087A (zh) | 可分享的半自动糖尿病眼底病变标注方法 | |
KR102060895B1 (ko) | 의료 영상 생성 방법 및 디바이스 | |
WO2014036638A1 (en) | Systems and methods for diagnosing strokes | |
CN108665454A (zh) | 一种内镜图像智能分类及不规则病变区域检测方法 | |
Thomas et al. | A novel multiscale and multipath convolutional neural network based age-related macular degeneration detection using OCT images | |
CN110021022A (zh) | 一种基于深度学习的甲状腺核医学图像诊断方法 | |
CN110415216A (zh) | 基于sd-oct和octa视网膜图像的cnv自动检测方法 | |
CN110400301A (zh) | 一种基于眼部特征的脑卒中人工智能筛查方法 | |
Ng et al. | Ultra-high-field neuroimaging reveals fine-scale processing for 3D perception | |
Kihara et al. | Detection of nonexudative macular neovascularization on structural OCT images using vision transformers | |
Liu et al. | Weakly-supervised localization and classification of biomarkers in OCT images with integrated reconstruction and attention | |
Wang et al. | Cycle adaptive multi-target weighting network for automated diabetic retinopathy segmentation | |
Hernanda et al. | Semantic segmentation of venous on deep vein thrombosis (DVT) case using UNet-ResNet |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |