磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法技术

本发明专利技术公开了一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法。它采集被测大脑的结构像，对结构像进行分割产生灰质概率图GM和白质概率图WM，灰质概率图GM和白质概率图WM组成参考图像pmCBF，利用参考图像pmCBF去除被测大脑在n个时间点的脑血流图CBF1～CBFn中的伪迹图。本发明专利技术能够有效去除采集到的脑血流图像序列中的伪迹图，降低误差，提高图像质量。

【技术实现步骤摘要】

本专利技术涉及磁共振动脉自旋标记灌注成像
，尤其涉及一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法。
技术介绍
脑灌注又称脑血流(cerebral blood flow，CBF)是一项重要的生理指标。测量脑血流能给临床诊断器质性或功能性脑病变提供定量依据。因为脑功能激活会导致局部脑血流变化，测量脑血流也成了脑功能研究的一个重要手段。临床上测量脑血流基本上都要注射外源性示踪剂或造影剂。但这些示踪剂或造影剂有不同程度的副作用比如有放射性或肾毒性，所以适用面常有限，也不能在短时间内多次重复测量。基于磁共振的动脉血自旋标记灌注成像(Arterial Spin Labeling(ASL)perfusion MRI)是基于磁共振成像(MRI)的一种完全无创无毒副作用的血流测量技术，在临床和脑科学研究中有非常广阔的应用前景。由于总体标记的时间有限和标记后自旋的传输过程中的信号衰减，ASL MRI的信噪比通常很小，信号只占背景的1％左右。在ASL MRI中灌注信号的提取是从带标记和没有标记的图像相减的差值中获取的。这个过程使得ASL MRI对于运动非常敏感。低信噪比又进一步恶化这个过程，导致ASL CBF图像序列中会出现大量的伪迹图。这些伪迹图对ASL CBF质量影响很大必须去除。以往的去除伪迹图技术仅仅考虑到运动和标准差，没有考虑CBF本身的特性，伪迹图去除不到位，无法降低误差，有时甚至会加大误差。
技术实现思路
本专利技术的目的是克服磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法去除伪迹图不到位的技术问题，提供了一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，其能够有效去除采集到的脑血流图像序列中的伪迹图，降低误差，提高图像质量。为了解决上述问题，本专利技术采用以下技术方案予以实现：本专利技术的磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，磁共振动脉自旋标记脑灌注成像数据包括被测大脑在n个时间点的脑血流图CBF1～CBFn，包括：采集被测大脑的结构像，对结构像进行分割产生灰质概率图GM和白质概率图WM，灰质概率图GM和白质概率图WM组成参考图像pmCBF，利用参考图像pmCBF去除被测大脑在n个时间点的脑血流图CBF1～CBFn中的伪迹图，去除伪迹图的方法包括以下步骤：S1：计算当前剩余的每一幅脑血流图与参考图像pmCBF的相关性系数CC；S2：从相关性系数CC最小的脑血流图开始，按照相关性系数CC从小到大的顺序依次检测对应脑血流图是否为伪迹图，当检测到某幅脑血流图是伪迹图时，删除该幅脑血流图，计算当前剩余的其他脑血流图的均值，得到平均图像meanCBF，计算参考图像pmCBF与平均图像meanCBF的均值，得到平均图像averageCBF，将得到的平均图像averageCBF作为新的参考图像pmCBF，重新执行步骤S1；判断被检测的脑血流图是否为伪迹图的方法包括以下步骤：计算当前剩余的所有脑血流图灰质区域的标准差std0(GM)及白质区域的标准差std0(WM)，计算除了被检测的脑血流图之外当前剩余的其他脑血流图灰质区域的标准差std1(GM)及白质区域的标准差std1(WM)，判断std1(GM)+std1(WM)是否小于std0(GM)+std0(WM)，如果是，则判断被检测的脑血流图是伪迹图，否则判断被检测的脑血流图不是伪迹图。在本技术方案中，先计算当前剩余的每一幅脑血流图与参考图像pmCBF的相关性系数CC，找出相关性系数CC第一小的脑血流图，检测该脑血流图是否为伪迹图，如果是，则将其删除并计算出新的参考图像pmCBF，从头开始重新检测剩余的脑血流图，如果不是，则检测相关性系数CC第二小的脑血流图，如果相关性系数CC第二小的脑血流图是伪迹图，则将其删除并计算出新的参考图像pmCBF，从头开始重新检测剩余的脑血流图，如果相关性系数CC第二小的脑血流图不是伪迹图，则检测相关性系数CC第三小的脑血流图，如果相关性系数CC第三小的脑血流图是伪迹图，则将其删除并计算出新的参考图像pmCBF，从头开始重新检测剩余的脑血流图，如果相关性系数CC第三小的脑血流图不是伪迹图，则检测相关性系数CC第四小的脑血流图，如此循环，直到检测到相关性系数CC最大的脑血流图，如果相关性系数CC最大的脑血流图是伪迹图，则将其删除并计算出新的参考图像pmCBF，从头开始重新检测剩余的脑血流图，如果相关性系数CC最大的脑血流图不是伪迹图，则结束检测，当前剩余的脑血流图都不是伪迹图。脑血流图和大脑质子密度成像或T1加权结构图像高度相似，灰质和白质区域内脑血流的变化比较平滑。现有的去除伪迹图方法都是基于图像的均值，当大量伪迹图存在时，均值本身就会受到影响，导致质量好的图像帧被抛弃，真正伪迹图被保留，且去除伪迹的方法与脑血流图像特性没有关系。本方法引入了由被测大脑的结构像分割产生的灰质概率图GM和白质概率图WM组成参考图像pmCBF，避免大量伪迹图存在时，质量好的图像被误判为伪迹图去除，充分利用脑血流图像特性，降低误差。作为优选，由结构像分割产生的灰质概率图GM和白质概率图WM组成的参考图像pmCBF＝GM*β+WM，β＞1，β是一个常量，表示灰质和白质的对比度。为了让参考图像pmCBF的灰质部分信号比白质部分高。作为优选，计算一幅脑血流图与参考图像pmCBF的相关性系数CC的方法包括以下步骤：将该幅脑血流图的像素按列依次拼接成一维数据阵列X＝[x1，x2，x3，…，xN]T，将参考图像pmCBF的像素按列依次拼接成一维数据阵列Y＝[y1，y2，y3，…，yN]T，x、y表示像素，T表示转置变换，其中，本专利技术的磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，磁共振动脉自旋标记脑灌注成像数据包括被测大脑在n个时间点的脑血流图CBF1～CBFn，包括：采集被测大脑的结构像，结构像和脑血流图CBF1～CBFn都是由M层单层扫描图组成的三维立体图，将脑血流图CBF1～CBFn的同一层单层脑血流扫描图组成一组，分别对每组图像进行处理，去除每组图像中的伪迹图；去除某一组图像中的伪迹图的方法包括以下步骤：N1：对结构像中与该组图像对应的一层单层结构像扫描图进行分割产生灰质概率图GM和白质概率图WM，灰质概率图GM和白质概率图WM组成参考图像pmCBF；N2：计算该组图像中当前剩余的每一幅单层脑血流扫描图与参考图像pmCBF的相关性系数CC；N3：从相关性系数CC最小的单层脑血流扫描图开始，按照相关性系数CC从小到大的顺序依次检测对应单层脑血流扫描图是否为伪迹图，当检测到某幅单层脑血流扫描图是伪迹图时，删除该幅单层脑血流扫描图，计算当前剩余的其他单层脑血流扫描图的均值，得到平均图像meanCBF，计算参考图像pmCBF与平均图像meanCBF的均值，得到平均图像averageCBF，将得到的平均图像averageCBF作为新的参考图像pmCBF，重新执行步骤N2；判断被检测的单层脑血流扫描图是否为伪迹图的方法包括以下步骤：取出该组图像中当前剩余的每一幅单层脑血流扫描图，将每一幅单层脑血流扫描图与其在所在的三维脑血流图中的相邻层脑血流扫描图取平均得到一幅平均脑血流扫描图，计算所有得到的平均脑血流扫描图灰质区域的标准差std2(GM)及白质区域的标准差s本文档来自技高网...

【技术保护点】
一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，磁共振动脉自旋标记脑灌注成像数据包括被测大脑在n个时间点的脑血流图CBF1～CBFn，其特征在于，包括：采集被测大脑的结构像，对结构像进行分割产生灰质概率图GM和白质概率图WM，灰质概率图GM和白质概率图WM组成参考图像pmCBF，利用参考图像pmCBF去除被测大脑在n个时间点的脑血流图CBF1～CBFn中的伪迹图，去除伪迹图的方法包括以下步骤：S1：计算当前剩余的每一幅脑血流图与参考图像pmCBF的相关性系数CC；S2：从相关性系数CC最小的脑血流图开始，按照相关性系数CC从小到大的顺序依次检测对应脑血流图是否为伪迹图，当检测到某幅脑血流图是伪迹图时，删除该幅脑血流图，计算当前剩余的其他脑血流图的均值，得到平均图像meanCBF，计算参考图像pmCBF与平均图像meanCBF的均值，得到平均图像averageCBF，将得到的平均图像averageCBF作为新的参考图像pmCBF，重新执行步骤S1；判断被检测的脑血流图是否为伪迹图的方法包括以下步骤：计算当前剩余的所有脑血流图灰质区域的标准差std0(GM)及白质区域的标准差std0(WM)，计算除了被检测的脑血流图之外当前剩余的其他脑血流图灰质区域的标准差std1(GM)及白质区域的标准差std1(WM)，判断std1(GM)+std1(WM)是否小于std0(GM)+std0(WM)，如果是，则判断被检测的脑血流图是伪迹图，否则判断被检测的脑血流图不是伪迹图。...

【技术特征摘要】
1.一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，磁共振动脉自旋标记脑灌注成像数据包括被测大脑在n个时间点的脑血流图CBF1～CBFn，其特征在于，包括：采集被测大脑的结构像，对结构像进行分割产生灰质概率图GM和白质概率图WM，灰质概率图GM和白质概率图WM组成参考图像pmCBF，利用参考图像pmCBF去除被测大脑在n个时间点的脑血流图CBF1～CBFn中的伪迹图，去除伪迹图的方法包括以下步骤：S1：计算当前剩余的每一幅脑血流图与参考图像pmCBF的相关性系数CC；S2：从相关性系数CC最小的脑血流图开始，按照相关性系数CC从小到大的顺序依次检测对应脑血流图是否为伪迹图，当检测到某幅脑血流图是伪迹图时，删除该幅脑血流图，计算当前剩余的其他脑血流图的均值，得到平均图像meanCBF，计算参考图像pmCBF与平均图像meanCBF的均值，得到平均图像averageCBF，将得到的平均图像averageCBF作为新的参考图像pmCBF，重新执行步骤S1；判断被检测的脑血流图是否为伪迹图的方法包括以下步骤：计算当前剩余的所有脑血流图灰质区域的标准差std0(GM)及白质区域的标准差std0(WM)，计算除了被检测的脑血流图之外当前剩余的其他脑血流图灰质区域的标准差std1(GM)及白质区域的标准差std1(WM)，判断std1(GM)+std1(WM)是否小于std0(GM)+std0(WM)，如果是，则判断被检测的脑血流图是伪迹图，否则判断被检测的脑血流图不是伪迹图。2.根据权利要求1所述的磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，其特征在于，由结构像分割产生的灰质概率图GM和白质概率图WM组成的参考图像pmCBF＝GM*β+WM，β＞1，β是一个常量，表示灰质和白质的对比度。3.根据权利要求1所述的磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，其特征在于，计算一幅脑血流图与参考图像pmCBF的相关性系数CC的方法包括以下步骤：将该幅脑血流图的像素按列依次拼接成一维数据阵列X＝[x1，x2，x3，...，xN]T，将参考图像pmCBF的像素按列依次拼接成一维数据阵列Y＝[y1，y2，y3，...，yN]T，x、y表示像素，T表示转置变换，其中，4.一种磁共振动脉自旋标记脑灌注成像数据伪迹图去除方法，磁共振动脉自旋标记脑灌注成像数据包括被测大脑在n个时间点的脑血流图CBF1～CBFn，其特征在于，包括：采集被测大脑的结构像，结构像和脑血流图CBF1～CBFn都是由M层单层扫描图组成的三...

【专利技术属性】
技术研发人员：王泽，
申请(专利权)人：杭州师范大学，
类型：发明
国别省市：浙江;33