本发明涉及图像处理技术领域，具体公开了一种脑卒中CT和MRI影像征象定量评估方法，结合CT图像的组织密度影征象和MRI图像的模糊效应特征，通过检测不同急性周期对应的脑部CT图像的豆状核轮廓并计算豆状核区灰质聚合曲线；并根据脑部MRI图像计算得到对应模糊效应拟合曲线，结合豆状核区灰质聚合曲线计算得到脑组织软化指数，能够从海量的脑卒中数据集的大数据中快速的筛选出灰白质信号指标、病灶清晰度和图像分辨率与患者最相似的脑卒中数据图像；可有效避免了由于CT存在骨伪影和MRI分辨率较低导致成像噪声误差干扰医护人员临床诊断结果的不足，作为监测脑卒中患者临床判断的重要指标和评估标准。
1.一种脑卒中CT和MRI影像征象定量评估方法，其特征在于，所述方法包括以下步骤：
S100，获取不同急性周期对应的脑部CT图像和脑部MRI图像；
S200，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线；
S300，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线；
S400，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数；
S500，当脑组织软化指数大于或等于阈值时，从脑卒中数据集中筛选得到和当前的脑部CT图像或脑部MRI图像的相似度最高的脑部CT图像或脑部MRI图像；
其中，在S200中，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线的方法为：S201，将不同急性周期对应的脑部CT图像进行图像配准处理，得到特征点匹配的CT配准图像，根据急性周期对应的时间段按时间顺序构成CT配准图像集合，以i为CT配准图像的序号，记multiCT(x,y,i)为CT配准图像集合中第i个急性周期对应的CT配准图像在像素坐标为(x,y)上的灰度值；其中，设CT配准图像将各像素按矩阵M×N大小排列，x∈[1,M]，y∈[1,N]，i∈[1, n]；
S202，遍历i值的取值范围，叠加所述CT配准图像集合中的所有CT配准图像得到重合CT配准图像，并计算所述CT配准图像集合中所有CT配准图像在各个像素坐标的灰度值对应的算术平均值得到平均CT配准图像；将重合CT配准图像与平均CT配准图像在各个对应的像素坐标上的灰度值分别作差得到CT差值图像；
S203，利用边缘检测算法检测CT差值图像，获得距离图像边界最远或长度大于所有边缘线的一条边缘线，或者，获取各个边缘线所构成的区域中面积最大的区域所对应的边缘线，记为豆状核轮廓边缘线；将豆状核轮廓边缘线所构成的闭合区域记作豆状核区域，提取所述CT配准图像集合中各个CT配准图像对应豆状核区域内的像素点，按照预设顺序存入豆状核区灰质结构矩阵集合，记multiLN(i, j)表示为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的第j个灰度值；其中，j∈[1, Num]，Num为豆状核区域内的像素点个数；
S204，由豆状核区灰质结构矩阵集合计算得到随着急性周期变化的豆状核区灰质聚合曲线，计算豆状核区灰质聚合值LNcurve(i)构成豆状核区灰质聚合曲线，豆状核区灰质聚合值的计算公式为： ；
其中，LNcurve(i)表示为豆状核区灰质聚合曲线上第i个豆状核区灰质聚合值；LNρ(j)为豆状核区灰质结构矩阵集合内对应的第j个灰质结构倾向度，meanLN(i)为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的平均灰度值；max[LNρ(j)]表示为遍历j取值范围计算LNρ(j)的绝对值中的最大值；
其中，在S300中，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线的方法为：S301，遍历i值取值范围，将第i个急性周期对应的T1加权成像、T2加权成像、磁共振弥散加权成像和液体反转恢复成像分别记作T1WI(i)、T2WI(i)、DWI(i) 和FLAIR(i)；
S302，遍历T1WI(i)、T2WI(i)、DWI(i)和FLAIR(i)的各个像素点在i值取值范围内进行累加求均值，并由各个像素点对应的均值构成平均图像，分别记作aveT1WI、aveT2WI、aveDWI和aveFLAIR；分别对aveT1WI、aveT2WI、aveDWI和aveFLAIR进行归一化处理，并利用OTSU阈值分割算法分别获得对应的非脑组织区域阈值，分别将aveT1WI、aveT2WI、aveDWI和aveFLAIR上数值大于或等于对应的非脑组织区域阈值的各个像素点标记为脑组织点，其余的像素点标记为非脑组织点；并根据脑组织点利用边界膨胀处理、边缘检测算法提取出对应的脑组织轮廓线，其脑组织轮廓线所构成的闭合区域记作对应的脑组织区域；
S303，分别对T1WI、T2WI、DWI和FLAIR对应的脑组织区域内进行图像增强处理，并利用模糊C均值聚类算法将图像划分得到对应的灰白质子区域以及脑脊液子区域；
S304，遍历i值取值范围，分别计算T1WI(i)、T2WI(i)、DWI(i)以及FLAIR(i)对应的灰白质子区域内所有像素点对应的灰度值的算术平均值，分别记作T1(i)、T2(i)、Dwg(i)、以及Fwg(i)；利用最小二乘法分别对T1(i)、T2(i)、Dwg(i)、以及Fwg(i)进行拟合，分别获得T1WI、T2WI、DWI和FLAIR对应的模糊效应拟合曲线，分别记作T1fitting(i)、T2fitting(i)、Dwgfitting(i)、以及Fwgfitting(i)；
其中，在S400中，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数的方法为：S401，根据豆状核区灰质聚合曲线和不同的模糊效应拟合曲线计算得到对应的脑组织模糊液化曲线，其计算公式为：；
其中BlurT1(i)表示为T1WI对应的随着i值变化的脑组织模糊液化曲线；BlurT2(i)表示为T2WI对应的随着i值变化的脑组织模糊液化曲线；
S402，由脑组织模糊液化曲线计算得到对应的脑组织软化指数，其计算公式为：
其中SoftenT1(i)表示为第i个急性周期时T1WI对应的脑组织软化指数；SoftenT2(i)表示为第i个急性周期时T2WI对应的脑组织软化指数。
2.根据权利要求1所述的一种脑卒中CT和MRI影像征象定量评估方法，其特征在于，在S100中，利用计算机断层扫描仪获取不同急性周期对应的脑部CT图像和核磁共振仪获得不同急性周期对应的脑部MRI图像；其中脑部MRI图像包括T1加权成像记作T1WI、T2加权成像记作T2WI、磁共振弥散加权成像记作DWI以及液体反转恢复成像记作FLAIR；其中，不同急性周期是根据脑卒中发病时间之后的时间段划分的周期，i值表示为急性周期序号，以6h作为一个急性周期的时间间隔，按i值从小到大顺序排列划分得到n个急性周期，i∈[1, n]，n∈[56, 120]。
3.根据权利要求1所述的一种脑卒中CT和MRI影像征象定量评估方法，其特征在于，在S500中，所述相似度的计算方法包括：结构相似性度量、余弦相似度算法、感知哈希算法、向量相似度算法、PSNR峰值信噪比法。
4.一种脑卒中CT和MRI影像征象定量评估系统，其特征在于，所述一种脑卒中CT和MRI影像征象定量评估系统包括：处理器、存储器及存储在所述存储器中并在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1至3所述任一项的一种脑卒中CT和MRI影像征象定量评估方法中的步骤，所述一种脑卒中CT和MRI影像征象定量评估系统还包括：脑卒中病变程度评估单元，用于根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度；
其中，根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度具体包括以下步骤：S501，分别计算SoftenT1(i)和SoftenT2(i)在i值取值范围内的算术平均值，分别记作meanST1和meanST2，以及分别计算Dwgfitting(i)和Fwgfitting(i) 在i值取值范围内的算术平均值，分别记作meanDwg和meanFwg；遍历i值取值范围，进一步由步骤S5031 S5034评估~脑卒中病变程度：
S5031，当SoftenT1(i)≤meanST1且SoftenT2(i)≤meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于超急性早期；否则跳转至S5034；
S5032，当SoftenT1(i)≤meanST1且SoftenT2(i)＞meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)＞meanFwg条件，则评估脑卒中在第i个急性周期处于超急性晚期；
S5033，当SoftenT1(i)＞meanST1时，进一步判断是否满足SoftenT2(i)＜meanST2条件，是则若满足Dwgfitting(i) ≤meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于亚急性早期；否则均评估脑卒中在第i个急性周期处于亚急性晚期；
S5034，判断是否满足i≤(n/2)值条件，是则评估脑卒中在第i个急性周期处于急性期；
所述一种脑卒中CT和MRI影像征象定量评估系统运行于桌上型计算机、笔记本电脑、掌上电脑及云端数据中心的计算设备中。
一种脑卒中CT和MRI影像征象定量评估方法及系统
技术领域
[0001] 本发明涉及图像处理技术领域，具体涉及一种脑卒中CT和MRI影像征象定量评估方法及系统。
背景技术
[0002] 随着影像学技术如计算机断层扫描(CT)与磁共振扫描(MRI)等在临床的广泛应用，CT技术能清晰呈现脑卒中病灶情况，区分出血性和缺血性类型；而MRI术能多方位直接成像，显示动脉瘤等血管异常。因此，通过分析脑卒中患者的特定典型的CT和MRI影像征象，如脑肿胀、脑萎缩、组织密度影、脑梗死范围、皮髓界限模糊以及颅内动脉钙化征等，能为临床动态观察、治疗诊断和判断预后提供技术指标。
[0003] 但是在基于临床影像的脑卒中诊疗过程中，由于影像数据涉及多模式多序列，医生难以短时间内依据大量影像数据准确找出脑卒中病灶并分析成因进行救治；而急救救治需要高时效性，会由于临床诊断效率与准确性的不足错过最佳时间窗造成患者不可逆性的损伤。而基于影像的疾病诊断应用主要依赖于医生临床经验用肉眼观察做出诊断，但医学影像存在噪声多、边缘模糊、层数多、数据量大等困难，且患者存在个体化差异因素强，常规的辅助诊断工具不足以保证高准确率。因此，目前亟需一种高效精准的基于影像学分析的脑卒中定量评估方法，能有效辅助临床诊断且评估预后，降低脑卒中患者的致死率、致残率和复发率。
[0004] 而CT和MRI影像在临床应用上仍存在一些局限性，CT成像无法避免颅骨等非脑组织的遮挡影响，不能清晰呈现小脑、脑干等病灶位置容易出现误判漏诊；MRI能显示脑部结构无死角无骨伪影，获得清晰的脑组织成像，但分辨率较低且成像速度较慢。近年来伴随影像学技术的不断革新，开展了应用多模影像对脑卒中病灶的研究，但目前仍未有将CT和MRI影像两者优势结合的脑卒中定量评估标准及应用，并且，在已经公开的海量的脑卒中数据集的大数据中，尽管已经收录有无数种脑卒中数据图像，医生难以根据灰白质信号指标、病灶清晰度和图像分辨率与患者最相似的脑卒中数据图像。
发明内容
[0005] 鉴于以上所述现有方法的局限，本发明的目的在于提供一种脑卒中CT和MRI影像征象定量评估方法及系统，结合CT图像的组织密度影征象和MRI图像的模糊效应特征，通过检测不同急性周期对应的脑部CT图像的豆状核轮廓并计算豆状核区灰质聚合曲线；并根据脑部MRI图像计算得到对应模糊效应拟合曲线，结合豆状核区灰质聚合曲线计算得到脑组织软化指数，进一步评估脑卒中病变程度，可有效避免了由于CT存在骨伪影和MRI分辨率较低导致成像噪声误差干扰医护人员临床诊断结果的不足，作为监测脑卒中患者临床判断的重要指标和评估标准，并且，在已经公开的海量的脑卒中数据集的大数据中，尽管已经收录有无数种脑卒中数据图像，医生难以根据灰白质信号指标、病灶清晰度和图像分辨率与患者最相似的脑卒中数据图像。
[0006] 为了实现上述目的，根据本公开的一方面，提供一种脑卒中CT和MRI影像征象定量评估方法，所述方法包括以下步骤：
[0007] S100，获取不同急性周期对应的脑部CT图像和脑部MRI图像；
[0008] S200，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线；
[0009] S300，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线；
[0010] S400，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数；
[0011] S500，当脑组织软化指数大于或等于阈值时，从脑卒中数据集中筛选得到和当前的脑部CT图像或脑部MRI图像的相似度最高的脑部CT图像或脑部MRI图像(能够从海量的脑卒中数据集的大数据中快速的筛选出灰白质信号指标、病灶清晰度和图像分辨率与患者最相似的脑卒中数据图像)。
[0012] 进一步地，在S100中，利用计算机断层扫描仪获取脑部CT图像和核磁共振仪获得不同急性周期对应的脑部MRI图像；其中脑部MRI图像包括T1加权成像记作T1WI、T2加权成像记作T2WI、磁共振弥散加权成像记作DWI以及液体反转恢复成像记作FLAIR；优选地，不同急性周期可根据脑卒中发病时间之后的时间段划分周期，i值表示为急性周期序号，以6h作为一个急性周期的时间间隔，按i值从小到大顺序排列划分得到n个急性周期，i∈[1, n]，n∈[56, 120]。
[0013] 进一步地，在S200中，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线的方法为：
[0014] S201，将不同急性周期对应的脑部CT图像进行图像配准处理，得到特征点匹配的CT配准图像，根据急性周期对应的时间段按时间顺序构成CT配准图像集合，以i为CT配准图像的序号，记multiCT(x,y,i)为CT配准图像集合中第i个急性周期对应的CT配准图像在像素坐标为(x,y)上的灰度值；其中，设CT配准图像将各像素按矩阵M×N大小排列，x∈[1,M]，y∈[1,N]，i∈[1, n]；
[0015] 优选地，图像配准处理可为Harris‑Laplacian提取算法，Susan角点算法、Harris角点算法等其中一种，特征点提取算子包括Canny、Moravec、Forstner、罗伯特算子等其中一种；
[0016] S202，遍历i值的取值范围，叠加所述CT配准图像集合中的所有CT配准图像得到重合CT配准图像，并计算所述CT配准图像集合中所有CT配准图像在各个像素坐标的灰度值对应的算术平均值得到平均CT配准图像；将重合CT配准图像与平均CT配准图像在各个对应的像素坐标上的灰度值分别作差得到CT差值图像；
[0017] S203，利用边缘检测算法检测CT差值图像，获得距离图像边界最远或长度大于所有边缘线的一条边缘线，或者，获取各个边缘线所构成的区域中面积最大的区域所对应的边缘线，记为豆状核轮廓边缘线；将豆状核轮廓边缘线所构成的闭合区域记作豆状核区域，提取所述CT配准图像集合中各个CT配准图像对应豆状核区域内的像素点，按照预设顺序存入豆状核区灰质结构矩阵集合，记multiLN(i, j)表示为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的第j个灰度值；其中，j∈[1, Num]，Num为豆状核区域内的像素点个数；
[0018] S204，由豆状核区灰质结构矩阵集合计算得到随着急性周期变化的豆状核区灰质聚合曲线，计算豆状核区灰质聚合值LNcurve(i)构成豆状核区灰质聚合曲线，豆状核区灰质聚合值的计算公式为：
[0019] ；
[0020] ；
[0021] ；
[0022] 其中，LNcurve(i)表示为豆状核区灰质聚合曲线上第i个急性周期对应的豆状核区灰质聚合值；LNρ(j)为豆状核区灰质结构矩阵集合内对应的第j个灰质结构倾向度，meanLN(i)为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的平均灰度值；max[LNρ(j)]表示为遍历j取值范围计算LNρ(j)的绝对值中的最大值。
[0023] 进一步地，在S300中，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线的方法为：
[0024] S301，遍历i值取值范围，将第i个急性周期对应的T1加权成像、T2加权成像、磁共振弥散加权成像和液体反转恢复成像分别记作T1WI(i)、T2WI(i)、DWI(i) 和FLAIR(i)；
[0025] S302，遍历T1WI(i)、T2WI(i)、DWI(i)和FLAIR(i)的各个像素点在i值取值范围内进行累加求均值，并由各个像素点对应的均值构成平均图像，分别记作aveT1WI、aveT2WI、aveDWI和aveFLAIR；分别对aveT1WI、aveT2WI、aveDWI和aveFLAIR进行归一化处理，并利用OTSU阈值分割算法分别获得对应的非脑组织区域阈值，分别将aveT1WI、aveT2WI、aveDWI和aveFLAIR上数值大于或等于对应的非脑组织区域阈值的各个像素点标记为脑组织点，其余的像素点标记为非脑组织点；并根据脑组织点利用边界膨胀处理、边缘检测算法提取出对应的脑组织轮廓线，其脑组织轮廓线所构成的闭合区域记作对应的脑组织区域；
[0026] S303，分别对T1WI、T2WI、DWI和FLAIR对应的脑组织区域内进行图像增强处理，并利用模糊C均值聚类算法将图像划分得到对应的灰白质子区域以及脑脊液子区域；
[0027] S304，遍历i值取值范围，分别计算T1WI(i)、T2WI(i)、DWI(i)以及FLAIR(i)对应的灰白质子区域内所有像素点对应的灰度值的算术平均值，分别记作T1(i)、T2(i)、Dwg(i)、以及Fwg(i)；利用最小二乘法分别对T1(i)、T2(i)、Dwg(i)、以及Fwg(i)进行拟合，分别获得T1WI、T2WI、DWI和FLAIR对应的模糊效应拟合曲线，分别记作T1fitting(i)、T2fitting(i)、Dwgfitting(i)、以及Fwgfitting(i)。
[0028] 进一步地，在S400中，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数的方法为：
[0029] S401，根据豆状核区灰质聚合曲线和不同的模糊效应拟合曲线计算得到对应的脑组织模糊液化曲线，其计算公式为：
[0030] ；
[0031] ；
[0032] 其中BlurT1(i)表示为T1WI对应的随着i值变化的脑组织模糊液化曲线；BlurT2(i)表示为T2WI对应的随着i值变化的脑组织模糊液化曲线；
[0033] S402，由脑组织模糊液化曲线计算得到对应的脑组织软化指数，其计算公式为：
[0034] ；
[0035] ；
[0036] 其中SoftenT1(i)表示为第i个急性周期时T1WI对应的脑组织软化指数；SoftenT2(i)表示为第i个急性周期时T2WI对应的脑组织软化指数。
[0037] 进一步地，在S500中，阈值为人工设置的值，或者，阈值为通过S100～S400计算5~10个脑卒中患者得到的脑组织软化指数的平均值，或者，阈值为脑卒中数据集中随机抽取5
100组脑部CT图像和脑部MRI图像数据通过S100～S400计算得到的脑组织软化指数的平均~
值。
[0038] 其中，脑卒中数据集包括： BraTS2018数据集、CQ500数据集、ISLE数据集、MRBrainS数据集, Hippocampus(MRI, 1759) 数据集、 Brain Hemarhage数据集。
[0039] 其中，相似度的计算方法包括：结构相似性度量、余弦相似度算法、感知哈希算法、向量相似度算法、PSNR峰值信噪比法。
[0040] 进一步地，在S500中，将相似度最高的脑部CT图像或脑部MRI图像推送到客户端或者移动设备中输出，以作为监测脑卒中患者临床判断的重要指标和评估标准。
[0041] 本发明还提供了一种脑卒中CT和MRI影像征象定量评估系统，所述系统包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序；所述处理器执行所述计算机程序时实现一种脑卒中CT和MRI影像征象定量评估方法中的步骤，[0042] 所述一种脑卒中CT和MRI影像征象定量评估系统还包括：
[0043] 豆状核区灰质结构处理单元，用于对不同急性周期对应的脑部CT图像检测得到豆状核轮廓边缘线，并计算对应的豆状核区灰质聚合曲线；
[0044] MRI模糊效应拟合处理单元，用于根据不同急性周期对应的脑部MRI图像计算得到T1WI、T2WI、DWI和FLAIR对应的模糊效应拟合曲线；
[0045] 脑组织软化指数计算单元，用于根据豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数；
[0046] 脑卒中病变程度评估单元，用于根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度。
[0047] 其中，根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度具体包括以下步骤：
[0048] S501，分别计算SoftenT1(i)和SoftenT2(i)在i值取值范围内的算术平均值，分别记作meanST1和meanST2，以及分别计算Dwgfitting(i)和Fwgfitting(i) 在i值取值范围内的算术平均值，分别记作meanDwg和meanFwg；遍历i值取值范围，进一步由步骤S5031 S5034~评估脑卒中病变程度：
[0049] S5031，当SoftenT1(i)≤meanST1且SoftenT2(i)≤meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于超急性早期；否则跳转至S5034；
[0050] S5032，当SoftenT1(i)≤meanST1且SoftenT2(i)＞meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)＞meanFwg条件，则评估脑卒中在第i个急性周期处于超急性晚期；
[0051] S5033，当SoftenT1(i)＞meanST1时，进一步判断是否满足SoftenT2(i)＜meanST2条件，是则若满足Dwgfitting(i) ≤meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于亚急性早期；否则均评估脑卒中在第i个急性周期处于亚急性晚期；
[0052] S5034，判断是否满足i≤(n/2)值条件，是则评估脑卒中在第i个急性周期处于急性期；否则评估脑卒中在第i个急性周期处于慢性期；
[0053] 所述一种脑卒中CT和MRI影像征象定量评估系统可以运行于桌上型计算机、笔记本、掌上电脑及云端数据中心等计算设备中。
[0054] 如上所述，本发明所述的一种脑卒中CT和MRI影像征象定量评估方法及系统，具有以下有益效果：（1）通过叠加抵消能较好地豆状核区域的模糊轮廓，识别豆状核区域，解决了CT成像由于骨伪影难以对脑灰白质的病灶清晰成像的局限性；（2）识别不同MRI图像的灰白质子区域计算得到对应的模糊效应拟合曲线，能提供更准确的灰白质信号指标，在一定程度上解决MRI成像分辨率较低的问题；（3）能够从海量的脑卒中数据集的大数据中快速的筛选出灰白质信号指标、病灶清晰度和图像分辨率与患者最相似的脑卒中数据图像；（4）结合CT图像的组织密度影特征和MRI图像的模糊效应特征，可有效避免由于CT存在骨伪影和MRI分辨率较低导致成像噪声误差干扰医护人员临床诊断结果的不足，计算得到脑组织软化指数进一步评估脑卒中病变程度，能作为监测脑卒中患者临床判断的重要指标和评估标准。
附图说明
[0055] 通过对结合附图所示出的实施方式进行详细说明，本公开的上述以及其他特征将更加明显，本公开附图中相同的参考标号表示相同或相似的元素，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，在附图中：
[0056] 图1所示为一种脑卒中CT和MRI影像征象定量评估方法于一实施例中的流程图；
[0057] 图2所示为一种脑卒中CT和MRI影像征象定量评估系统于一实施例中的系统结构图。
具体实施方式
[0058] 以下将结合实施例和附图对本公开的构思、具体结构及产生的技术效果进行清楚、完整的描述，以充分地理解本公开的目的、方案和效果。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
[0059] 如图1所示为根据本发明的一种脑卒中CT和MRI影像征象定量评估方法的流程图，下面结合图1来阐述根据本发明的实施方式的一种脑卒中CT和MRI影像征象定量评估方法。本公开提出一种脑卒中CT和MRI影像征象定量评估方法，所述方法具体包括以下步骤：
[0060] S100，获取不同急性周期对应的脑部CT图像和脑部MRI图像；
[0061] S200，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线；
[0062] S300，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线；
[0063] S400，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数；
[0064] S500，当脑组织软化指数大于或等于阈值时，从脑卒中数据集中筛选得到和当前的脑部CT图像或脑部MRI图像的相似度最高的脑部CT图像或脑部MRI图像。
[0065] 进一步地，在S100中，利用计算机断层扫描仪获取脑部CT图像和核磁共振仪获得不同急性周期对应的脑部MRI图像；其中脑部MRI图像包括T1加权成像记作T1WI、T2加权成像记作T2WI、磁共振弥散加权成像记作DWI以及液体反转恢复成像记作FLAIR；优选地，不同急性周期可根据脑卒中发病时间之后的时间段划分周期，i值表示为急性周期序号，以6h（6小时）作为一个急性周期的时间间隔，按i值从小到大顺序排列划分（或者也可以为：脑卒中发病时间到当前时刻的时间段）得到n个急性周期，i∈[1, n]，n∈[56, 120]；
[0066] 优选地，所述计算机断层扫描仪采用双排螺旋CT断层扫描，以om线为扫描基线，层厚5.0 mm，螺距1.5 mm，扫描参数为160 kV，50 mAs；所述核磁共振仪的视场角为250×220 mm，扫描尺寸为256×256，扫描参数分别为FSE‑T1WI（TE/TR= 9/2000 ms），FSE‑T2WI（TE/TR= 105/6000 ms），DWI（TE/TR= 102/5500 ms），FLAIR（TE/TR= 96/8000 ms）。
[0067] 进一步地，在S200中，获得不同急性周期对应的脑部CT图像的豆状核轮廓，并计算豆状核区灰质聚合曲线的方法为：
[0068] S201，将不同急性周期对应的脑部CT图像进行图像配准处理，得到特征点匹配的CT配准图像，根据急性周期对应的时间段按时间顺序构成CT配准图像集合，以i为CT配准图像的序号，记multiCT(x,y,i)为CT配准图像集合中第i个急性周期对应的CT配准图像在像素坐标为(x,y)上的灰度值；其中，设CT配准图像将各像素按矩阵M×N大小排列，x∈[1,M]，y∈[1,N]，i∈[1, n]；
[0069] 优选地，图像配准处理可为Harris‑Laplacian提取算法，Susan角点算法、Harris角点算法等其中一种，特征点提取算子包括Canny、Moravec、Forstner、罗伯特算子等其中一种；
[0070] S202，遍历i值的取值范围，叠加所述CT配准图像集合中的所有CT配准图像得到重合CT配准图像，并计算所述CT配准图像集合中所有CT配准图像在各个像素坐标的灰度值对应的算术平均值得到平均CT配准图像；将重合CT配准图像与平均CT配准图像在各个对应的像素坐标上的灰度值分别作差得到CT差值图像；
[0071] S203，利用边缘检测算法检测CT差值图像，获得距离图像边界最远或长度大于所有边缘线的一条边缘线，或者，获取各个边缘线所构成的区域中面积最大的区域所对应的边缘线，记为豆状核轮廓边缘线；将豆状核轮廓边缘线所构成的闭合区域记作豆状核区域，提取所述CT配准图像集合中各个CT配准图像对应豆状核区域内的像素点，按照预设顺序存入豆状核区灰质结构矩阵集合，记multiLN(i, j)表示为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的第j个灰度值；其中，预设顺序为将豆状核区域内像素坐标为(x,y)对应的(x+y)值按从小到大顺序排列，其排列中第j个像素坐标(x,y)对应的灰度值对应为multiLN(i, j)中的第j个灰度值，j∈[1, Num]，Num为豆状核区域内的像素点个数；
[0072] S204，由豆状核区灰质结构矩阵集合计算得到随着急性周期变化的豆状核区灰质聚合曲线，计算豆状核区灰质聚合值LNcurve(i)构成豆状核区灰质聚合曲线，豆状核区灰质聚合值的计算公式为：
[0073] ；
[0074] ；
[0075] ；
[0076] 其中，LNcurve(i)表示为豆状核区灰质聚合曲线上第i个急性周期对应的豆状核区灰质聚合值；LNρ(j)为豆状核区灰质结构矩阵集合内对应的第j个灰质结构倾向度，meanLN(i)为豆状核区灰质结构矩阵集合中第i个急性周期对应的CT配准图像在豆状核区域内的平均灰度值；max[LNρ(j)]表示为遍历j取值范围计算LNρ(j)的绝对值中的最大值。（其原理为：当脑灰质的血流灌注减少会导致CT图像征象上出现豆状核征、脑岛带征及皮髓界限模糊，而豆状核征是由于大脑中动脉供血区细胞水肿导致的，随着脑卒中的急性周期增加，对应的豆状核轮廓会逐渐模糊或密度逐渐降低，对应灰白质界限消失，脑实质密度降低呈片状或斑片状的低密度影，处于慢性期时会发生“模糊效应”，即受损脑组织逐渐出现软化灶，在CT图像上表现为脑脊液样低密度影，但由于骨伪影难以对脑灰白质的病灶清晰成像，无法提供完整的灰白质疾病的诊断依据；而步骤S200通过叠加抵消能较好地豆状核区域的模糊轮廓，识别豆状核区域并计算豆状核区灰质聚合曲线，可在无清晰CT成像的条件下反映豆状核区域内脑实质的X线吸收系数随着急性周期变化的密度差异，有效量化低密度影的变化程度）。
[0077] 进一步地，在S300中，根据不同急性周期对应的脑部MRI图像计算得到模糊效应拟合曲线的方法为：
[0078] S301，遍历i值取值范围，将第i个急性周期对应的T1加权成像、T2加权成像、磁共振弥散加权成像和液体反转恢复成像分别记作T1WI(i)、T2WI(i)、DWI(i) 和FLAIR(i)；
[0079] S302，遍历T1WI(i)、T2WI(i)、DWI(i)和FLAIR(i)的各个像素点在i值取值范围内进行累加求均值，并由各个像素点对应的均值构成平均图像，分别记作aveT1WI、aveT2WI、aveDWI和aveFLAIR；分别对aveT1WI、aveT2WI、aveDWI和aveFLAIR进行归一化处理，并利用OTSU阈值分割算法分别获得对应的非脑组织区域阈值，分别将aveT1WI、aveT2WI、aveDWI和aveFLAIR上数值大于或等于对应的非脑组织区域阈值的各个像素点标记为脑组织点，其余的像素点标记为非脑组织点；并根据脑组织点利用边界膨胀处理、边缘检测算法提取出对应的脑组织轮廓线，其脑组织轮廓线所构成的闭合区域记作对应的脑组织区域；
[0080] S303，分别对T1WI、T2WI、DWI和FLAIR对应的脑组织区域内进行图像增强处理，并利用模糊C均值聚类算法将图像划分得到对应的灰白质子区域以及脑脊液子区域；
[0081] S304，遍历i值取值范围，分别计算T1WI(i)、T2WI(i)、DWI(i)以及FLAIR(i)对应的灰白质子区域内所有像素点对应的灰度值的算术平均值，分别记作T1(i)、T2(i)、Dwg(i)、以及Fwg(i)；利用最小二乘法分别对T1(i)、T2(i)、Dwg(i)、以及Fwg(i)进行拟合，分别获得T1WI、T2WI、DWI和FLAIR对应的模糊效应拟合曲线，分别记作T1fitting(i)、T2fitting(i)、Dwgfitting(i)、以及Fwgfitting(i)。（其原理为：由于脑部MRI图像由脑组织和非脑组织两部分组成，其中非脑组织包括头皮、头骨等，脑组织包括脑白质、脑灰质、脑脊液等；CT成像会存在骨伪影干扰，而MRI能获得清晰的脑组织成像，通过步骤S302可有效去除各个MRI图像中非脑组织部分，克服了CT成像由于颅骨等非脑组织的遮挡易出现误判漏诊的局限；而步骤S303 S304识别不同MRI图像的灰白质子区域计算得到对应的模糊效应拟合曲线，能提~供更准确的灰白质信号指标，在一定程度上解决MRI成像分辨率较低的问题）。
[0082] 进一步地，在S400中，结合豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数的方法为：
[0083] S401，根据豆状核区灰质聚合曲线和不同的模糊效应拟合曲线计算得到对应的脑组织模糊液化曲线，其计算公式为：
[0084] ；
[0085] ；
[0086] 其中BlurT1(i)表示为T1WI对应的随着i值变化的脑组织模糊液化曲线；BlurT2(i)表示为T2WI对应的随着i值变化的脑组织模糊液化曲线；
[0087] S402，由脑组织模糊液化曲线计算得到对应的脑组织软化指数，其计算公式为：
[0088] ；
[0089] ；
[0090] 其中SoftenT1(i)表示为第i个急性周期时T1WI对应的脑组织软化指数；SoftenT2(i)表示为第i个急性周期时T2WI对应的脑组织软化指数。（有益效果为：步骤S400可结合CT图像的组织密度影特征和MRI图像的模糊效应特征，可有效避免了由于CT存在骨伪影和MRI分辨率较低导致成像噪声误差干扰医护人员临床诊断结果的不足，能获得更高检测精度的灰白质病灶较为完整的诊断依据，通过计算不同急性周期对应的脑组织软化指数并结合MRI图像对应的模糊效应拟合曲线，能够准确的从脑卒中数据集中筛选出脑组织软化和患者最接近的图像数据集，使医护人员能够通过筛选出的图像数据集准确的进一步评估脑卒中病变程度，从而能准确的有效辅助临床诊断且评估预后）。
[0091] 进一步地，在S500中，阈值为人工设置的值，或者，阈值为通过S100～S400计算5~10个脑卒中患者得到的脑组织软化指数的平均值，或者，阈值为脑卒中数据集中随机抽取5
100组脑部CT图像和脑部MRI图像数据通过S100～S400计算得到的脑组织软化指数的平均~
值。
[0092] 其中，脑卒中数据集包括： BraTS2018数据集、CQ500数据集、ISLE数据集、MRBrainS数据集, Hippocampus(MRI, 1759) 数据集、 Brain Hemarhage数据集。
[0093] 其中，相似度的计算方法包括：结构相似性度量、余弦相似度算法、感知哈希算法、向量相似度算法、PSNR峰值信噪比法。
[0094] 如图2所示为本公开的一种脑卒中CT和MRI影像征象定量评估系统结构图，所述处理器执行所述计算机程序运行在以下系统的单元中：
[0095] 脑部图像获取单元，用于获取不同急性周期对应的脑部CT图像和脑部MRI图像；
[0096] 豆状核区灰质结构处理单元，用于对不同急性周期对应的脑部CT图像检测得到豆状核轮廓边缘线，并计算对应的豆状核区灰质聚合曲线；
[0097] MRI模糊效应拟合处理单元，用于根据不同急性周期对应的脑部MRI图像计算得到T1WI、T2WI、DWI和FLAIR对应的模糊效应拟合曲线；
[0098] 脑组织软化指数计算单元，用于根据豆状核区灰质聚合曲线和模糊效应拟合曲线计算得到脑组织软化指数；
[0099] 脑卒中病变程度评估单元，用于根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度；
[0100] 其中，根据不同急性周期对应的脑组织软化指数和MRI对应的模糊效应拟合曲线量化评估脑卒中病变程度具体包括以下步骤：
[0101] S501，分别计算SoftenT1(i)和SoftenT2(i)在i值取值范围内的算术平均值，分别记作meanST1和meanST2，以及分别计算Dwgfitting(i)和Fwgfitting(i) 在i值取值范围内的算术平均值，分别记作meanDwg和meanFwg；遍历i值取值范围，进一步由步骤S5031 S5034~评估脑卒中病变程度：
[0102] S5031，当SoftenT1(i)≤meanST1且SoftenT2(i)≤meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于超急性早期；否则跳转至S5034；
[0103] S5032，当SoftenT1(i)≤meanST1且SoftenT2(i)＞meanST2时，若满足Dwgfitting(i)≥meanDwg且Fwgfitting(i)＞meanFwg条件，则评估脑卒中在第i个急性周期处于超急性晚期；
[0104] S5033，当SoftenT1(i)＞meanST1时，进一步判断是否满足SoftenT2(i)＜meanST2条件，是则若满足Dwgfitting(i) ≤meanDwg且Fwgfitting(i)≤meanFwg条件，则评估脑卒中在第i个急性周期处于亚急性早期；否则均评估脑卒中在第i个急性周期处于亚急性晚期；
[0105] S5034，判断是否满足i≤(n/2)值条件，是则评估脑卒中在第i个急性周期处于急性期；否则评估脑卒中在第i个急性周期处于慢性期；
[0106] 优选地，在本具体实施例中，结合MRI信号进一步评估得到随着急性周期递增分别对应着不同的病变阶段；其中，
[0107] （1）超急性期是红细胞内氧合血红蛋白阶段，氧合血红蛋白抗磁性，主要分为早期和晚期，超急性早期对应的SoftenT1(i)和SoftenT2(i)未见异常，但和Dwgfitting(i)呈现高信号；而超急性晚期血肿含水较多，SoftenT2(i)和Fwgfitting(i)均呈现高信号；
[0108] （2）急性期是红细胞内脱氧血红蛋白顺磁阶段，但脱氧血红蛋白结构所限，SoftenT1(i)、SoftenT2(i)、Dwgfitting(i)和Fwgfitting(i)均呈现低信号；
[0109] （3）亚急性期也分为早期和晚期，亚急性早期是红细胞内正铁血红蛋白顺磁阶段，SoftenT2(i)呈现低信号，但由于正铁血红蛋白（和脱氧比铁原子位置变了）转正使SoftenT1(i)呈现高信号，Dwgfitting(i)和Fwgfitting(i)保持低信号；而亚急性晚期时红细胞外正铁血红蛋白顺磁，且血肿里面水分增加，SoftenT1(i)、SoftenT2(i)、Dwgfitting(i)和Fwgfitting(i)均呈现高信号；
[0110] （4）而慢性期时是红细胞外含铁血黄素超顺磁阶段，SoftenT1(i)和SoftenT2(i) 均呈现低信号，Dwgfitting(i)和Fwgfitting(i)保持低信号。
[0111] 所述一种脑卒中CT和MRI影像征象定量评估系统可以运行于桌上型计算机、笔记本、掌上电脑及云端数据中心等计算设备中。所述一种脑卒中CT和MRI影像征象定量评估系统包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，所述例子仅仅是一种脑卒中CT和MRI影像征象定量评估方法及系统的示例，并不构成对一种脑卒中CT和MRI影像征象定量评估方法及系统的限定，可以包括比例子更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述一种脑卒中CT和MRI影像征象定量评估系统还可以包括输入输出设备、网络接入设备、总线等。
[0112] 所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器 (Digital Signal Processor，DSP)、专用集成电路 (Application Specific Integrated Circuit，ASIC)、现场可编程门阵列 (Field‑ProgrammHCO3‑le Gate Array，FPGA) 或者其他可编程逻辑器件、分立元器件门电路或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述一种脑卒中CT和MRI影像征象定量评估系统的控制中心，利用各种接口和线路连接整个一种脑卒中CT和MRI影像征象定量评估系统的各个分区域。
[0113] 所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述一种脑卒中CT和MRI影像征象定量评估方法及系统的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡（Smart Media Card, SMC），安全数字（Secure Digital, SD）卡，闪存卡（Flash Card）、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
[0114] 尽管本公开的描述已经相当详尽且特别对几个所述实施例进行了描述，但其并非旨在局限于任何这些细节或实施例或任何特殊实施例，从而有效地涵盖本公开的预定范围。此外，上文以发明人可预见的实施例对本公开进行描述，其目的是为了提供有用的描述，而那些目前尚未预见的对本公开的非实质性改动仍可代表本公开的等效改动。
