[发明专利]一种TOF-MRA条纹伪影修正方法、装置、电子设备及存储介质

说明书

本发明涉及一种TOF‑MRA条纹伪影修正方法、装置、电子设备及存储介质，属于图像处理的领域；其中，TOF‑MRA条纹伪影修正方法包括：获取TOF‑MRA原始图像，并将TOF‑MRA原始图像作为输入图像；对输入图像进行NLM滤波平滑，获得模版图像；将TOF‑MRA原始图像与模版图像相减，获得残差图像；利用均值滤波器对残差图像进行分离，获得条纹伪影图像；将TOF‑MRA原始图像与条纹伪影图像相减，获得TOF‑MRA修正图像。本发明将NLM滤波平滑和均值滤波器应用于TOF‑MRA的条纹伪影修正中，能够在不损失TOF‑MRA原始图像细节和不降低TOF‑MRA原始图像分辨率的前提下，实现对条纹伪影的修正。

技术领域

本发明涉及图像处理的技术领域，尤其是涉及一种TOF-MRA条纹伪影修正方法、装置、电子设备及存储介质。

背景技术

三维时间飞跃法磁共振血管造影(three-dimensional time-of-flightmagnetic resonance angiography,TOF-MRA)包括2D和3D成像模式。2D采用单板采集方法，板内血流饱和度较轻，有利于显示颅内小血管和矢状窦，但是，2D的分辨率很低。3D通常具有较高的信噪比，可以获得较高的空间分辨率图像，是神经成像的理想选择。多层重叠薄板采集(multiple overlapping thin slab acquisitions,MOTSA)则是将2D采集的对比度优势与3D采集的高信噪比相结合，在保持信噪比的同时进一步降低了血液信号饱和度。但是，MOSTA技术需要多个薄板重叠，才能达到薄板数据采集的目的，导致了成像时间的增加。然而，减小重叠区域可以缩短成像时间，但会导致图像产生条纹伪影，破坏图像的细节信息，影响图像质量。因此，开发一种实用的图像后处理方法，通过减少条纹伪影来提高MOTSA数据的图像性能至关重要。

发明内容

为了在不损失TOF-MRA原始图像细节和不降低TOF-MRA原始图像分辨率的前提下，实现对条纹伪影的修正，第一方面，本发明提供一种TOF-MRA条纹伪影修正方法，采用如下的技术方案：

一种TOF-MRA条纹伪影修正方法，包括：

S11、获取TOF-MRA原始图像，并将所述TOF-MRA原始图像作为输入图像；

S12、对所述输入图像进行NLM滤波平滑，获得模版图像；

S13、将所述TOF-MRA原始图像与所述模版图像相减，获得残差图像；

S14、利用均值滤波器对所述残差图像进行分离，获得条纹伪影图像；

S15、将所述TOF-MRA原始图像与所述条纹伪影图像相减，获得TOF-MRA修正图像。

通过采用上述技术方案，通过NLM滤波平滑获得没有细节和条纹伪影的模版图像，进而获得包含图像细节和条纹伪影的残差图像残差图像，并利用均值滤波器对残差图像进行分离获得条纹伪影图像，以根据条纹伪影图像对TOF-MRA原始图像进行修正，能够在不损失TOF-MRA原始图像细节和不降低TOF-MRA原始图像分辨率的前提下，实现对条纹伪影的修正。

可选的，在所述步骤S12中，NLM滤波的公式是：

其中，p和q是输入图像的像素点，d(B(p),B(q))是以p,q为中心点的邻域块之间的欧式距离，f表示一个单调递减函数，欧式距离越近，值越大，C(p)是归一化函数；

图像u＝(u₁,u₂,u₃)中的像素点p的降噪函数是：

则

其中，B(p,r)表示以p为中心且大小为(2r+1)×(2r+1)个像素的邻域；