基于距离场融合的MR图像海马体分割方法技术

一种基于距离场融合的MR图像海马体分割方法。包括：（1）对初始待分割的MR测试图像归一化处理，去除颅骨和偏移场，得到归一化处理后的MR测试图像；（2）将MR训练集图像和训练集标号图像配准到MR测试图像；（3）将配准后的标号图像经过距离变换得到距离场DF；（4）对MR测试图像中的点取图像块并转化成一个列向量；（5）定义搜索窗选取MR字典和DF字典；（6）用局部线性表达测试样本，求字典权值系数向量；（7）得到测试样本的DF预测图像块向量，并把转化为图像块；（8）重复步骤（4）-（7），获得每个点的DF值；（9）得到测试图像每个点所对应的标号。本发明专利技术方法能够从MR图像中准确地分割目标。

【技术实现步骤摘要】
基于距离场融合的MR图像海马体分割方法
本专利技术涉及医学图像分析
，具体涉及一种基于距离场融合的MR图像海马体分割方法。
技术介绍
海马体属于中脑的灰质结构，海马体萎缩是阿尔兹海默病、精神分裂症、抑郁症等精神疾病的一种病理表现，临床上可以通过分析海马体的体积和形态对这些精神疾病做出诊断。MRI可以提供高分辨率解剖结构、清晰对比度和广泛的图像序列，这使得MRI技术在临床诊断方面广受欢迎。医生可以通过病人的脑部MR图像，对病人的脑部结构进行分析，从而对疾病做出诊断。对MR脑部图像的海马体进行分割，是海马体体积大小的测量、形态分析的前提。由于海马体的手动分割非常耗时、枯燥，而且不同医生的分割结果也会有所差异，海马体的自动分割在临床应用上显得尤为重要。由于MR图像中海马体的体积很小、形状复杂、边界模糊、部分容积效应明显，海马体的自动分割非常具有挑战性。迄今为止，MR图像海马体分割方法主要有如下两类：第一类是基于活动轮廓模型(ACM)的方法。这类方法是根据图像灰度统计学信息或者图像梯度对曲线(曲面)进行演化，最终曲线(曲面)收敛于目标边界处。ACM方法对图像噪声比较敏感。第二类是基于多图谱的方法。这类方法把图谱配准到测试图像，对形变后的训练集标号图像，通过某种方法实现标号融合，得到最终的分割结果。标号融合并没有利用到待分割目标的形状先验信息，分割精度有所限制。针对现有技术的不足，本专利技术方法提出一种基于距离场融合的MR图像海马体分割方法。
技术实现思路
本专利技术方法提供一种基于距离场融合的MR图像海马体分割方法，本专利技术方法成功地应用在MR脑部图像的海马体分割，能准确地分割MR图像中的海马体。本专利技术的上述目的通过如下技术手段实现。一种基于距离场融合的MR图像海马体分割方法，该方法基于两种假设：Ⅰ，MR图像块和DF图像块位于两个非线性流形上，任意一个MR图像块可以由其所在流形的局部空间内的近邻样本线性表达；Ⅱ，在局部约束条件下，MR流形到DF流形的映射近似于一个微分同胚映射；该基于距离场融合的MR图像海马体分割方法通过如下步骤进行：(1)对初始待分割的MR测试图像归一化处理，去除颅骨和偏移场，得到归一化处理后的MR测试图像；(2)将预先准备的训练集中的MR训练集图像和与MR训练集图像对应的训练集标号图像配准到MR测试图像，使得MR训练集图像、训练集标号图像与MR测试图像在空间位置上对齐；(3)对步骤(2)配准后的训练集标号图像进行距离变换，得到训练集标号图像的DF训练集距离场，点x所对应的DF训练集距离场的值为：其中C表示分割目标的边界，B表示距离x最近的点，并且B∈C，dist(x，B)表示点x和点B之间的欧氏距离；(4)对MR测试图像中的点x，以x为中心取一个图像块xMR，把xMR转化成一个列向量用来作为点x的特征，得到测试样本其中m代表向量的维数；(5)在训练集MR图像、训练集距离场DF图像上分别以点x为中心定义一个搜索窗选取图像块和构建MR字典和DF字典n为字典的大小；(6)基于假设Ⅰ，用MR字典DMR局部线性表示MR测试样本字典DMR权值系数向量用LAE方法求解，具体表达式如下：表示测试样本在字典DMR中的k个近邻；(7)用步骤(6)所得的字典DMR权值系数向量线性组合DF字典xDF中的样本，得到测试样本的DF预测图像块向量并把转化为图像块xDF；具体是：基于假设Ⅰ，可得：基于假设Ⅱ，可得：由于f是局部线性的，故：把转化为图像块xDF；(8)对MR测试图像中的每个点重复按照步骤(4)-(7)操作，获得每个点的DF训练集距离场的值；用P(x)表示以点x为中心、大小与xDF相等的图像块，对于P(x)中的任一点u，其权重为：点x的DF训练集距离场的值为：表示以点u为中心的图像块所预测的点x的DF训练集距离场的值；(9)对所预测出来的DF预测图像进行阈值处理，得到测试图像每个点所对应的标号，具体是：由式(1)所定义的距离场可知，点x的标号可以由式(8)导出：标号为1表示该点属于分割目标，标号为0表示背景，根据分割目标获得海马体分割图像。优选的,上述步骤(1)具体是采用灰度归一化方法对MR测试图像进行归一化处理，用BET算法去除颅骨，用N4算法去除偏移场。优选的,上述步骤(2)具体是使用DRAMMS软件工具将训练集中的MR训练集图像和训练集标号图像配准到MR测试图像。优选的,上述步骤(4)具体是提取MR测试图像中的每个点的特征，用来作为测试样本。优选的,上述步骤(5)具体是在MR训练集图像和DF训练集距离场中使用局部搜索窗构建MR字典和DF字典。优选的,上述步骤(8)具体是对重叠的DF预测图像块加权平均处理，最终获得每个点的DF训练集距离场的预测值。本专利技术方法提供一种基于距离场融合的MR图像海马体分割方法，本专利技术方法能准确地分割MR图像中的海马体，可用于阿尔兹海默病等精神疾病的诊断。附图说明利用附图对本专利技术作进一步的说明，但附图中的内容不构成对本专利技术的任何限制。图1为本专利技术方法的流程示意图；图2为本专利技术方法的框架图；图3为本专利技术方法的细节流程示意图；图4为本专利技术方法用不同大小的搜索窗对15个测试图像的分割结果箱线图；其中，图4(a)为针对右海马体的分割结果箱线图，图4(b)为针对左海马体的分割结果箱线图；图5为用标号融合方法和距离场融合方法对15个测试图像的分割结果对比图，其中，图5(a)为针对右海马体的分割结果对比图，图5(b)为针对左海马体的分割结果对比图；图6为采用距离场融合方法和标号融合方法在四个测试图像上分割结果的冠状面示意图，其中，第一行表示手动分割结果，第二行是采用距离场融合方法的分割结果，第三行是采用标号融合方法的分割结果，每一列表示一个测试图像。具体实施方式结合以下实施例对本专利技术作进一步描述。实施例1。一种基于距离场融合的MR图像海马体分割方法，该方法基于两种假设：Ⅰ，MR图像块和DF图像块位于两个非线性流形上，任意一个MR图像块可以由其所在流形的局部空间内的近邻样本线性表达；Ⅱ，在局部约束条件下，MR流形到DF流形的映射近似于一个微分同胚映射。该基于距离场融合的MR图像海马体分割方法通过如下步骤进行：(1)对初始待分割的MR测试图像归一化处理，去除颅骨和偏移场，得到归一化处理后的MR测试图像。具体是采用灰度归一化方法对MR测试图像进行归一化处理，用BET算法去除颅骨，用N4算法去除偏移场。(2)将预先准备的训练集中的MR训练集图像和与MR训练集图像对应的训练集标号图像配准到MR测试图像，使得MR训练集图像、训练集标号图像与MR测试图像在空间位置上对齐。步骤(2)中具体是使用DRAMMS软件工具将训练集中的MR训练集图像和训练集标号图像配准到MR测试图像。(3)对步骤(2)配准后的训练集标号图像进行距离变换，得到训练集标号图像的训练集距离场DF，点x所对应的训练集距离场DF的值为：其中C表示分割目标的边界，B表示距离x最近的点，并且B∈C，dist(x，B)表示点x和点B之间的欧氏距离。(4)对MR测试图像中的点x，以x为中心取一个图像块xMR，把xMR转化成一个列向量用来作为点x的特征，其中m代表向量的维数。(5)在训练集MR图像、训练集距离场DF图像上分别以点x为中心定义一个搜索窗选本文档来自技高网...

【技术保护点】
一种基于距离场融合的MR图像海马体分割方法，其特征在于：该基于距离场融合的MR图像海马体分割方法基于两种假设：Ⅰ，MR图像块和DF图像块位于两个非线性流形上，任意一个MR图像块可以由其所在流形的局部空间内的近邻样本线性表达；Ⅱ，在局部约束条件下，MR流形到DF流形的映射近似于一个微分同胚映射；该基于距离场融合的MR图像海马体分割方法通过如下步骤进行：（1）对初始待分割的MR测试图像归一化处理，去除颅骨和偏移场，得到归一化处理后的MR测试图像；（2）将预先准备的训练集中的MR训练集图像和与MR训练集图像对应的训练集标号图像配准到MR测试图像，使得MR训练集图像、训练集标号图像与MR测试图像在空间位置上对齐；（3）对步骤（2）配准后的训练集标号图像进行距离变换，得到训练集标号图像的训练集距离场DF，点所对应的训练集距离场DF的值为：……  （1）；其中C表示分割目标的边界，B表示距离x最近的点，并且，表示点x和点B之间的欧氏距离；（4）对MR测试图像中的点，以为中心取一个图像块，把转化成一个列向量，用来作为点的特征，其中m代表向量的维数；（5）在训练集MR图像、训练集距离场DF图像上分别以点为中心定义一个搜索窗选取图像块和，构建MR字典和DF字典，n为字典的大小；（6）基于假设Ⅰ，用MR字典局部线性表示MR测试样本，字典权值系数向量用LAE方法求解，具体表达式如下：，……（2）；；表示测试样本在字典中的k个近邻；（7）用步骤（6）所得的字典权值系数向量线性组合DF字典中的样本，得到测试样本的DF预测图像块向量，并把转化为图像块；具体是：基于假设Ⅰ，可得：……（3）；基于假设Ⅱ，可得：……（4）；由于f是局部线性的，故：……（5）；把转化为图像块即可得到预测的距离场图像块；（8）对MR测试图像中的每个点重复按照步骤（4）‑（7）操作，获得每个点的DF值；用表示以点x为中心、大小与相等的图像块，对于中的任一点u，其权重为：……（6）；点x的DF值为：……（7）；表示以点u为中心的图像块所预测的点x的DF值；（9）对所预测出来的DF预测图像进行阈值处理，得到测试图像每个点所对应的标号，具体是：由式（1）所定义的距离场可知，点x的标号可以由式（8）导出：…… （8）；标号为1表示该点属于分割目标，标号为0表示背景，根据分割目标获得海马体分割图像。...

【技术特征摘要】
1.一种基于距离场融合的MR图像海马体分割方法，其特征在于：该基于距离场融合的MR图像海马体分割方法基于两种假设：I，MR图像块和DF图像块位于两个非线性流形上，任意一个MR图像块由其所在流形的局部空间内的近邻样本线性表达；II，在局部约束条件下，MR流形到DF流形的映射近似于一个微分同胚映射；该基于距离场融合的MR图像海马体分割方法通过如下步骤进行：(1)对初始待分割的MR测试图像归一化处理，去除颅骨和偏移场，得到归一化处理后的MR测试图像；(2)将预先准备的训练集中的MR训练集图像和与MR训练集图像对应的训练集标号图像配准到MR测试图像，使得MR训练集图像、训练集标号图像与MR测试图像在空间位置上对齐；(3)对步骤(2)配准后的训练集标号图像进行距离变换，得到训练集标号图像的DF训练集距离场，点x所对应的DF训练集距离场的值为：其中C表示分割目标的边界，B表示距离x最近的点，并且B∈C，dist(x，B)表示点x和点B之间的欧氏距离；(4)对MR测试图像中的点x，以x为中心取一个图像块xMR，把xMR转化成一个列向量用来作为点x的特征，得到测试样本其中m代表向量的维数；(5)在MR训练集图像、DF训练集距离场图像上分别以点x为中心定义一个搜索窗选取图像块和构建MR字典和DF字典n为字典的大小；(6)基于假设I，用MR字典DMR局部线性表示MR测试样本字典DMR权值系数向量用LAE方法求解，具体表达式如下：表示测试样本在字典DMR中的k个近邻；(7)用步骤(6)所得的字典DMR权值系数向量线性组合DF字典DDF中的样本，得到测试样本的DF预测图像块...

【专利技术属性】
技术研发人员：冯前进，庞树茂，阳维，卢振泰，
申请(专利权)人：南方医科大学，
类型：发明
国别省市：广东;44