非配对多中心磁共振数据和谐化方法、存储介质及应用技术

本发明专利技术涉及一种非配对多中心磁共振数据和谐化方法、存储介质及应用，所述方法包括以下步骤：获取不同厂家的磁共振图像，对所述磁共振图像进行标准化预处理；以循环生成对抗网络的架构为基础，构建初步的风格转换网络；以预处理后的磁共振图像，双向训练所述风格转换网络至收敛，获得训练好的网络参数；指定一厂家的磁共振图像为目标风格，其余为原始风格，加载所述训练好的网络参数，将具有原始风格的第一图像转换具有所述目标风格的第二图像，所述第一图像为预处理后的磁共振图像；进行从所述第一图像到第二图像的直方图匹配。与现有技术相比，本发明专利技术基于非配对图像生成更加接近于目标风格的图像，更好地保留原始解剖结构信息，提高数据可靠性。提高数据可靠性。提高数据可靠性。

【技术实现步骤摘要】
非配对多中心磁共振数据和谐化方法、存储介质及应用


[0001]本专利技术属于计算机辅助医学
，涉及一种医学影像数据处理方法，尤其是涉及一种非配对多中心磁共振数据和谐化方法、存储介质及应用。

技术介绍

[0002]磁共振成像作为一项重要的临床检查，通过丰富的扫描序列，可以对特定临床问题，进行不同对比度的磁共振成像，为临床问题提供丰富的诊断信息。临床医生通常是根据图像的定性特征对图像进行解释，因此不同的医生，甚至是同一医生不同时间给出的诊断结论都不尽相同。特别是当图像存在质量较差、扫描序列不规范等额外因素时，要获得一致的诊断尤为困难。另外人脑的生理结构决定了临床医生对于数据的分析局限于定性特征和宏观特征，无法对定量特征和隐藏在图像之下的复杂特征形式进行识别。因此，诞生了计算机辅助诊断(Computer
‑
Aided Diagnosis，CAD)，可以将其作为“第二名临床医生”，对疾病进行辅助决策，来提高临床诊断的准确性。许多工作都证明了计算机辅助诊断能够帮助临床医生更好进行诊断。
[0003]随着以深度学习为代表的人工智能的迅速发展，越来越多的网络模型被用于计算机辅助诊断。但人工智能在计算机辅助诊断领域面临一个很大的问题，就是诊断模型难以迁移。不同医院采集的数据，很难用同一个模型进行诊断。造成这种模型难以迁移的原因，就在于模型对不同医院扫描的图像之间的差异较为敏感。造成这种差异的原因可以归结为以下两个方面：首先，成像设备不同，即便是同一个厂家的设备，型号也可能不同；另外，扫描参数设置也存在很大差异。另外，计算机辅助诊断中所涉及到的自动检测、自动分割等模型的训练都极大地依赖于医生手工绘制的病灶感兴趣区域(ROI)。因此，利用每家医院的数据重新训练模型也很困难。
[0004]目前，解决这一问题的思路之一是数据标准化，由相关的医院或机构统一制订针对某一具体临床问题诊断的扫描与读片。这一标准化的过程工作量大，需要协调的问题多，并且难以利用大量的、宝贵的回顾性图像资料进行模型训练。除此之外，基于特征数据的和谐化或利用深度学习网络的和谐化也是目前常见的方法，但前者依赖于勾画好的ROI或者优秀的分割结果，后者限制于配对数据集，都存在一定缺陷。

技术实现思路

[0005]本专利技术的目的就是为了克服上述现有技术中存在的工作量大、依赖勾画ROI、限制于配对数据集的问题，而提供一种非配对多中心磁共振数据和谐化方法、设备及存储介质。
[0006]本专利技术的目的可以通过以下技术方案来实现：
[0007]一种非配对多中心磁共振数据和谐化方法，包括以下步骤：
[0008]获取不同厂家的磁共振图像，对所述磁共振图像进行标准化预处理；
[0009]以循环生成对抗网络的架构为基础，构建初步的风格转换网络；
[0010]以预处理后的磁共振图像，双向训练所述风格转换网络至收敛，获得训练好的网
络参数；
[0011]指定一厂家的磁共振图像为目标风格，其余为原始风格，加载所述训练好的网络参数，将具有原始风格的第一图像转换具有所述目标风格的第二图像，所述第一图像为预处理后的磁共振图像；
[0012]进行从所述第一图像到第二图像的直方图匹配。
[0013]进一步地，所述磁共振数据为由3T扫描仪扫描获得的矢状位的T2加权像。
[0014]进一步地，所述标准化预处理包括重采样标准和强度归一化标准。
[0015]进一步地，所述重采样标准具体为：
[0016]基于磁共振数据对应的分辨率统计结果，以一设定分辨率对分辨率超出设定范围的图像进行重采样。
[0017]进一步地，所述设定分辨率基于分辨率统计结果确定。
[0018]进一步地，所述强度归一化标准具体为：
[0019]利用大津算法分割经重采样标准处理后图像的前景和后景区域，对前景区域所包含的像素点按升序进行排列，利用设定的像素点截断阈值进行归一化。
[0020]进一步地，所述像素点截断阈值基于前后截断比例确定。
[0021]进一步地，所述循环生成对抗网络包括采用ResUnet的生成器和采用渐进式生成对抗网络的鉴别器。
[0022]进一步地，所述初步的风格转换网络中，网络参数初始化为均值为0、标准差为0.02的分布。
[0023]进一步地，所述直方图匹配具体为：
[0024]将第一图像和第二图像的图像强度进行归一化处理；
[0025]获取图像的累计直方图分布，计算灰阶映射矩阵，得到灰阶映射表；
[0026]根据灰阶映射表完成从第一图像到第二图像的直方图匹配。
[0027]本专利技术还提供一种计算机可读存储介质，包括供电子设备的一个或多个处理器执行的一个或多个程序，所述一个或多个程序包括用于执行如上所述非配对多中心磁共振数据和谐化方法的指令。
[0028]本专利技术还提供一种如上所述的非配对多中心磁共振数据和谐化方法在磁共振图像分类中的应用，利用和谐化后图像建立图像分类模型。
[0029]与现有技术相比，本专利技术具有以下有益效果：
[0030]1、本专利技术针对不同厂家磁共振采集设备间的差异，设计一个风格转换网络，可以通过不同厂家磁共振设备采集图像之间的风格转换来实现非配对多中心磁共振数据的一致性，从而克服诊断模型难以迁移的缺点，更方便地扩展数据集。
[0031]2、本专利技术在风格转换网络后加入直方图匹配策略，以约束灰阶映射结果，更好地实现解剖结构信息的保真。
[0032]3、本专利技术方法可以基于非配对图像生成更加接近于目标风格的图像，更好地保留原始解剖结构信息，提高数据可靠性，实现更准确的分类诊断效果。
附图说明
[0033]图1为本专利技术方法的流程框图；
[0034]图2为分辨率的统计结果示意图；
[0035]图3为动态归一化的结果示意图，(3a)为重采样后的图像，(3b)为施加大津算法去除背景后的图像区域，(3c)为动态强度归一化后的图像；
[0036]图4为基于cycleGAN的风格转换网络的总体架构示意图，其中，Set A为GE图像，Set B为PHILIPS图像；
[0037]图5为直方图匹配的结果示意图，(5a)为经过预处理的GE图像，(5b)为cycleGAN网络输出的GE图像，(5c)为直方图匹配之后的GE图像，(5d)为(5a)、(5b)和(5c)对应的累计直方图，其中，方框中的不规则实线框是放射科医生勾画的感兴趣病灶区域；
[0038]图6为分类诊断模型的结果示意图，(6a)为和谐化前训练集和测试集的分类AUC，(6b)为用直方图均衡(CLAHE)方法和谐化后的训练集和测试集的分类AUC，(6c)为用本专利技术方法和谐化后的训练集和测试集的分类AUC。
具体实施方式
[0039]下面结合附图和具体实施例对本专利技术进行详细说明。本实施例以本专利技术技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种非配对多中心磁共振数据和谐化方法，其特征在于，包括以下步骤：获取不同厂家的磁共振图像，对所述磁共振图像进行标准化预处理；以循环生成对抗网络的架构为基础，构建初步的风格转换网络；以预处理后的磁共振图像，双向训练所述风格转换网络至收敛，获得训练好的网络参数；指定一厂家的磁共振图像为目标风格，其余为原始风格，加载所述训练好的网络参数，将具有原始风格的第一图像转换具有所述目标风格的第二图像，所述第一图像为预处理后的磁共振图像；进行从所述第一图像到第二图像的直方图匹配。2.根据权利要求1所述的非配对多中心磁共振数据和谐化方法，其特征在于，所述磁共振数据为由3T扫描仪扫描获得的矢状位的T2加权像。3.根据权利要求1所述的非配对多中心磁共振数据和谐化方法，其特征在于，所述标准化预处理包括重采样标准和强度归一化标准。4.根据权利要求3所述的非配对多中心磁共振数据和谐化方法，其特征在于，所述重采样标准具体为：基于磁共振数据对应的分辨率统计结果，以一设定分辨率对分辨率超出设定范围的图像进行重采样。5.根据权利要求3所述的非配对多中心磁共振数据和谐化方法，其特征在于，所述强度归一化标准具体为：利用大津算法分割经重采样标准处理后图像的前...

【专利技术属性】
技术研发人员：蔡昕，常晓，杨光，聂生东，姜小平，郑慧，
申请(专利权)人：上海理工大学，
类型：发明
国别省市：