百叶窗伪影的消除方法、装置、终端设备及可读存储介质制造方法及图纸

本发明专利技术实施例提供了一种百叶窗伪影的消除方法、装置、终端设备及可读存储介质，所述方法包括：获取待处理图像的直流分量，图像包括百叶窗伪影；获取待处理图像的空间域图像的亮度信息；根据直流分量对亮度信息进行修正，得到处理后的图像，处理后的图像的亮度相同。本发明专利技术实施例通过傅里叶变换得到每一幅图像的直流分量，利用这组直流分量对空间域图像的亮度进行修正，获得亮度一致的图像，从而消除百叶窗伪影，解决血管成像的百叶窗伪影问题，快速有效，易于实现。易于实现。易于实现。

Method, device, terminal device and readable storage medium for eliminating louver artifacts

【技术实现步骤摘要】
百叶窗伪影的消除方法、装置、终端设备及可读存储介质


[0001]本专利技术涉及图像处理
，特别是涉及一种百叶窗伪影的消除方法、装置、终端设备及可读存储介质。

技术介绍

[0002]在磁共振血管成像技术中，时间飞跃法(Time of Flight,TOF)这种技术主要是利用了梯度回波序列中流动血液的流入增强效应进行血管成像的，可以叫做流入增强血管成像技术。流入增强血管成像技术最大的特点就是无需使用对比剂可以达到血管成像的目的。该技术主要使用梯度回波序列，通过比较短的TR设置来抑制静止的背景组织信号，从而突显血管，然而该技术也存在一个比较明显的问题，那就是随着血液的流入，这些流动的血液组织由于也会因为受到多个不同层面的射频脉冲激发而产生饱和效应，其结果就是远端的血管显示逐渐变差。
[0003]解决血液饱和现象导致的远端血管显示不清，主要可以采用下列方法。
[0004]首先是逆血流方向采集信号，既然随着血流流入容易产生饱和，那么在信号采集的时候可以逆方向采集，通过这种方法可以有效地减少血液饱和效应。
[0005]其次采用了多叠块重叠扫描组合技术(Multi
‑
Chunk)，该技术主要是在3D的TOF序列中将一个成像容积分为多个组块(Chunk)，这样每个组块内部的饱和效应会减少，然后将重叠的饱和部分剔除进行拼接，这样形成的血管图像则相对比较均匀。
[0006]第三种减少血液饱和的方法是优化射频脉冲的翻转角，这种技术叫做倾斜优化非饱和激励。该技术的原理主要是基于翻转角越大，血液信号越大，但是在相同的TR短时间内其纵向驰豫恢复也越不充分，越容易饱和。TONE技术采用斜坡式的翻转角。具体的角度设置则是一开始血流流入的时候采用较小的翻转角以减少近端血液的饱和效应，随着血液流入，逐渐增大翻转角，这样远端的血液由于翻转角大信号也会增加从而弥补饱和造成的信号下降。
[0007]但是现有技术中的这些处理方式，只能减少百叶窗伪影，并不能消除，因此，如何能够快速消除百叶窗伪影是目前急需解决的问题。

技术实现思路

[0008]鉴于上述问题，提出了本专利技术实施例以便提供一种克服上述问题或者至少部分地解决上述问题的一种百叶窗伪影的消除方法、装置、终端设备及可读存储介质。
[0009]第一个方面，本专利技术实施例提供一种百叶窗伪影的消除方法，所述方法包括：
[0010]获取待处理图像的直流分量，所述图像包括百叶窗伪影；
[0011]获取所述待处理图像的空间域图像的亮度信息；
[0012]根据所述直流分量对所述亮度信息进行修正，得到处理后的图像，所述处理后的图像的亮度相同。
[0013]可选地，所述获取待处理图像的直流分量，包括：
[0014]对所述待处理图像进行傅里叶变换，得到所述待处理图像的多个直流分量。
[0015]可选地，所述根据所述直流分量对所述亮度信息进行修正，得到处理后的图像，包括：
[0016]对所述空间域图像进行傅里叶变换，得到频率域图像信号；
[0017]根据所述频率域图像信号，确定频率最大值；
[0018]根据所述频率最大值，确定频率平均值；
[0019]根据所述频率平均值对所述空间域图像进行修正，得到处理后的图像。
[0020]可选地，所述对所述空间域图像进行傅里叶变换，得到频率域图像信号，包括：
[0021]设定S
k
是待处理图像的空间域图像的信号，对S
k
做傅里叶变换，得到频率域信号f
k
，f
k
通过如下公式得到：
[0022]f
k
＝F(S
k
)，离散形式为：
[0023][0024]k—一组图像中的索引；
[0025]u，v—频率变量；
[0026]x，y—空间域图像变量,即空间位置；
[0027]M，N
‑‑
为图像的长和宽。
[0028]可选地，所述根据所述频率最大值，确定频率平均值，包括：
[0029]所述频率平均值通过如下方式得到：
[0030][0031]其中，avrD为D
k
的平均值；T为图像总幅数。
[0032]可选地，所述根据所述频率平均值对所述空间域图像进行修正，得到处理后的图像，包括：
[0033]逐一像素对S
k
修正，得到图像CS
k
；
[0034][0035]其中，为系数，S
k
(x,y)为待处理图像的每一个像素。
[0036]第二个方面，本专利技术实施例提供一种百叶窗伪影的消除装置，所述装置包括：
[0037]第一获取模块，用于获取待处理图像的直流分量，所述图像包括百叶窗伪影；
[0038]第二获取模块，用于获取所述待处理图像的空间域图像的亮度信息；
[0039]修正模块，用于根据所述直流分量对所述亮度信息进行修正，得到处理后的图像，所述处理后的图像的亮度相同。
[0040]可选地，所述第一获取模块用于：
[0041]对所述待处理图像进行傅里叶变换，得到所述待处理图像的多个直流分量。
[0042]可选地，所述修正模块，用于：
[0043]对所述空间域图像进行傅里叶变换，得到频率域图像信号；
[0044]根据所述频率域图像信号，确定频率最大值；
[0045]根据所述频率最大值，确定频率平均值；
[0046]根据所述频率平均值对所述空间域图像进行修正，得到处理后的图像。
[0047]可选地，所述修正模块，用于：
[0048]设定S
k
是待处理图像的空间域图像的信号，对S
k
做傅里叶变换，得到频率域信号f
k
，f
k
通过如下公式得到：
[0049]f
k
＝F(S
k
)，离散形式为：
[0050][0051]k—一组图像中的索引；
[0052]u，v—频率变量；
[0053]x，y—空间域图像变量,即空间位置；
[0054]M，N
‑‑
为图像的长和宽。
[0055]可选地，所述修正模块用于：
[0056]所述频率平均值通过如下方式得到：
[0057][0058]其中，avrD为D
k
的平均值；T为图像总幅数。
[0059]可选地，所述修正模块用于：
[0060]逐一像素对S
k
修正，得到图像CS
k
；
[0061][0062]其中，为系数，S
k
(x,y)为待处理图像的每一个像素。
[0063]第三个方面，本专利技术实施例提供一种终端设备，包括：至少一个处理器和存储器；
[0064]所述存储器存本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种百叶窗伪影的消除方法，其特征在于，所述方法包括：获取待处理图像的直流分量，所述图像包括百叶窗伪影；获取所述待处理图像的空间域图像的亮度信息；根据所述直流分量对所述亮度信息进行修正，得到处理后的图像，所述处理后的图像的亮度相同。2.根据权利要求1所述的方法，其特征在于，所述获取待处理图像的直流分量，包括：对所述待处理图像进行傅里叶变换，得到所述待处理图像的多个直流分量。3.根据权利要求2所述的方法，其特征在于，所述根据所述直流分量对所述亮度信息进行修正，得到处理后的图像，包括：对所述空间域图像进行傅里叶变换，得到频率域图像信号；根据所述频率域图像信号，确定频率最大值；根据所述频率最大值，确定频率平均值；根据所述频率平均值对所述空间域图像进行修正，得到处理后的图像。4.根据权利要求3所述的方法，其特征在于，所述对所述空间域图像进行傅里叶变换，得到频率域图像信号，包括：设定S
k
是待处理图像的空间域图像的信号，对S
k
做傅里叶变换，得到频率域信号f
k
，f
k
通过如下公式得到：f
k
＝F(S
k
)，离散形式为：k
‑
一组图像中的索引；u，v
‑
频率变量；x，y
‑
空间域图像变量，即空间位置；M，N
‑‑
为图像的长和宽。5.根据权利要求...

【专利技术属性】
技术研发人员：张少斌，
申请(专利权)人：深圳市贝斯达医疗股份有限公司，
类型：发明
国别省市：