CN107085840A - 基于分数阶偏微分方程的图像去噪方法 - Google Patents

基于分数阶偏微分方程的图像去噪方法 Download PDF

Info

Publication number
CN107085840A
CN107085840A CN201710456403.0A CN201710456403A CN107085840A CN 107085840 A CN107085840 A CN 107085840A CN 201710456403 A CN201710456403 A CN 201710456403A CN 107085840 A CN107085840 A CN 107085840A
Authority
CN
China
Prior art keywords
mrow
msub
msup
msubsup
alpha
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710456403.0A
Other languages
English (en)
Other versions
CN107085840B (zh
Inventor
李远禄
丁亚庆
孟霄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Xinshiyun Science and Technology Co.,Ltd.
Original Assignee
Nanjing University of Information Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201710456403.0A priority Critical patent/CN107085840B/zh
Publication of CN107085840A publication Critical patent/CN107085840A/zh
Application granted granted Critical
Publication of CN107085840B publication Critical patent/CN107085840B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/70—Denoising; Smoothing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

本发明公开了一种基于分数阶偏微分方程的图像去噪方法，本发明根据分数阶导数的非局部性质，在检测边缘时能够减弱噪声的干扰，结合偏微分方程得到一种基于分数阶偏微分方程的图像去噪方法，能够在去噪的同时尽可能地保留原图像的纹理细节；在求解的过程中采用了快速傅立叶变换的方法，避免了复杂的分数阶导数展开运算的同时加快了求解速度；本发明将扩散函数的变量单独设定了分数阶导数，对于不同的图像变化微分阶数可以获得较好的去噪效果，并且收敛速度也较快，所需的迭代次数较少。

Description

基于分数阶偏微分方程的图像去噪方法
技术领域
本发明属于图像处理技术领域，特别是一种基于分数阶偏微分方程的图像去噪方法。
背景技术
随着模式识别、人工智能、遥感科学、计算机视觉等研究领域的不断深入，图像处理问题已经成为当前倍受关注的热点，它是图像分析和图像理解的前提和基础。图像处理主要包括图像去噪、图像增强、图像复原、图像分割、图像压缩、图像融合、图像超分辨率、图像特征提取等问题。在获取和传输过程中，由于传感器元器件的缺陷，传输信道受到干扰，环境的变化如光线和天气的影响，以及后续传输中的人为因素等，会造成信号受到噪声污染。虽然这些噪声的来源和形成因素不同，但它们都会导致原始图像细节模糊，从而进一步影响图像的后续处理，如图像特征提取、图像理解等。因此，对图像进行去噪成为了图像处理的重要步骤之一。
传统的图像去噪方法在去除图像噪声的同时，往往会模糊图像的细节，主要是因为图像的纹理细节和噪声区分不够准确，此外好的去噪方法一般需要复杂的运算消耗大量资源。为了能够达到很好的去噪效果同时又能够提高运算速度和降低资源消耗，就需要引入一些新颖的数学思想和算法，例如分数阶微积分，快速Fourier变换等等，来研究新的图像去噪方法。
图像去噪对图像处理的其他分支也有影响。近年来一些人工智能领域的研究热点例如基于计算机视觉的无人驾驶、即时定位与地图构建、人脸识别、指纹识别等，都和图像去噪中使用的方法有一定的联系。因此，一个好的图像去噪方法往往具有启发性的影响，可以进一步推广到图像处理的其他领域。从这个角度考虑，深入研究图像去噪中一些根本性问题，提出新颖的图像去噪方法，对于促进图像处理的其它分支的发展，也具有重要意义。
偏微分方程应用到图像处理领域最早可以追溯到物理学中的热扩散现象，它描述一个区域内的温度如何随时间的变化，热传导方程用于图像去噪可以表述为随着时间的变化图像的噪声点逐渐趋近于平滑，即达到去噪效果。传统的热传导方程在去噪的同时会导致图像的边缘模糊纹理细节丢失。为了克服这一缺点，Perona和Malik在1990年提出了基于边缘扩散的各向异性扩散模型即著名的P-M扩散模型
该模型利用灰度图像梯度的模|▽u|作为检测图像边缘的工具，同时设计了基于变量|▽u|单调递减的扩散函数g(|▽u|)，对图像的不同区域进行不同程度的扩散。在图像的平坦区域，|▽u|的值一般较小，此时对应的g(|▽u|)较大，因此在平坦区域具有较强的去噪能力同时不影响纹理细节；在图像的边缘区域，|▽u|的值一般较大，此时对应的g(|▽u|)很小，因而在边缘区域去噪能力较弱，因此能够最大程度地保留原图像的纹理细节。但P-M扩散模型的解过度依赖于参数的设定和初始条件，参数设定不当或者特殊的初始条件会导致方程解的不稳定，另外P-M扩散模型在边缘检测时容易出现虚假边缘容易导致去噪处理后的图像存在阶梯效应。
针对求解的问题，Yu-Li You等人分析了P-M扩散方程解的稳定性条件，同时给出了求解稳定的数值计算方法[3]。Catte等人针对P-M扩散模型存在阶梯效应的问题，对该模型进行了一些改进，提出了正则化后的P-M方程。Patrick等人提出了两种改进的四阶偏微分方程模型，重新设计了扩散函数并讨论了参数的影响，一定程度上提升了保留边缘的能力并抑制了阶梯效应，但该方法的去噪效果比较一般；现有技术还未解决在去除噪声的同时尽可能地保留原图像的纹理信息的图像去噪方法，本发明解决这样的问题。
发明内容
为解决现有技术的不足，本发明的目的在于提供一种基于分数阶偏微分方程的各向异性扩散的图像还原方法，本方法根据综合分数阶导数的非局部特性以及各向异性扩散模型的保边缘特性，在去除噪声的同时尽可能地保留原图像的纹理信息，极大地提升了图像还原效果。
为了实现上述目标，本发明采用如下的技术方案：
基于分数阶偏微分方程的图像去噪方法，包括如下步骤：
步骤一，输入一幅被噪声污染的图像u0(x,y)，图像大小为M×N，设定时间间隔Δt＝0.25；
步骤二，将图像进行对称处理，处理后图像的大小为原图像的4倍，计算公式如下：
for x＝1:M
for y＝1:N
u(x,2*N+1-y)＝u0(x,y)
u(2*M+1-x,y)＝u0(x,y)
u(2*M+1-x,2*N+1-y)＝u0(x,y)
end
end；
步骤三，利用MATLAB商业数学软件自带的快速傅立叶变换函数对图像进行傅立叶变换，公式如下：
步骤四，获得图像的分数阶微分；
计算图像的分数阶微分的计算过程如下：
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α，p(mi)表示在频域上对图像的第i行或者第i列的微分算子；
步骤五，获得图像分数阶梯度的模，计算过程如下：
将步骤四中频域上分数阶微分进行傅立叶反变换，得到空间域上的图像的分数阶微分，并将参数α换为扩散函数的参数β，
根据梯度的模的计算公式得到分数阶梯度的模：
步骤六，根据步骤五中求得的分数阶梯度的模计算扩散函数：
g(x)＝1/(1+x2)，
步骤七，计算微分算子的共轭算子，计算过程如下：
p*(m1)＝conj((1-exp(-j2πm1/N))α)，
输出一次迭代求解的结果：
其中：
步骤九，将步骤八中所得结果进行傅立叶反变换，
步骤十，将un+1还原到原图像大小，计算过程如下：
for x＝1:M
for y＝1:N
end
end；
步骤十一，利用下式计算所得结果与原图像进行峰值信噪比计算，判断是否满足终止条件；
其中f为无噪声图像，u为一次迭代还原后的图像，若符合终止条件：ΔPSNR≤ε，其中，ΔPSNR表示两次迭代结果的峰值信噪比的变化绝对值，ε取0.01，则执行步骤十三；若ΔPSNR＞ε，则执行步骤十二；
步骤十二，重新迭代求解，将上一次的迭代结果作为下一次迭代的输入，继续执行步骤四至步骤十一；
步骤十三，输出还原后的结果则为最优的去噪图像。
前述的基于分数阶偏微分方程的图像去噪方法，步骤一，输入一幅被噪声污染的图像u0(x,y)，图像大小为M×N，设定时间间隔Δt≤0.25。
前述的基于分数阶偏微分方程的图像去噪方法，步骤四，获得图像的分数阶微分；
计算图像的分数阶微分的计算过程如下：
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α，p(mi)表示在频域上对图像的第i行或者第i列的微分算子；
图像中，某一点对x方向的一阶偏导数可以写成差分形式：
相应的，根据分数阶导数的定义，某一点对x方向的α阶导数可以表示为：
对上式两边同时进行傅立叶变换可得该点对x方向和y方向的α阶导数的傅立叶变换，公式如下：
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α。
前述的基于分数阶偏微分方程的图像去噪方法，本发明采用的软件是MATLAB。
前述的基于分数阶偏微分方程的图像去噪方法，本发明需要的最低配置的处理器为：I3-2100处理器。
前述的基于分数阶偏微分方程的图像去噪方法，本发明采用的处理器为：Intel(R)Core(TM)i5-6500CPU@3.20GHz处理器。
前述的基于分数阶偏微分方程的图像去噪方法，本发明需要的最低内存是1GBRAM。
前述的基于分数阶偏微分方程的图像去噪方法，本发明采用的内存为16GB RAM。
本发明的有益之处在于：
第一，本发明根据分数阶导数的非局部性质，在检测边缘时能够减弱噪声的干扰，结合偏微分方程得到一种基于分数阶偏微分方程的图像去噪方法，能够在去噪的同时尽可能地保留原图像的纹理细节。
第二，在求解的过程中采用了快速Fourier变换的方法，避免了复杂的分数阶导数展开运算的同时加快了求解速度。
第三，本发明将扩散函数的变量单独设定了分数阶导数，对于不同的图像变化微分阶数可以获得较好的去噪效果，并且收敛速度也较快，所需的迭代次数较少。
附图说明
图1是本发明的一种实施例的流程图；
图2(a)是无噪声的Lena图像；
图2(b)是添加标准差δ＝25高斯白噪声后的Lena图像；
图2(c)使用本发明在实验中得到的去除噪声后的Lena图像；
图3(a)是无噪声的Barbara图像；
图3(b)是添加标准差δ＝25高斯白噪声后的Barbara图像；
图3(c)使用本发明在实验中得到的去除噪声后的Barbara图像；
图4(a)是无噪声的Boat图像；
图4(b)是添加标准差δ＝25高斯白噪声后的Boat图像；
图4(c)使用本发明在实验中得到的去除噪声后的Boat图像。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
基于分数阶偏微分方程的图像去噪方法，包括如下步骤：
步骤一，输入一幅被噪声污染的图像u0(x,y)，图像大小为M×N，设定时间间隔Δt≤0.25；需要说明的是，若时间间隔大于0.25会使得去噪不稳定，所以优选Δt≤0.25；作为一种实施例，如图2所示，图像大小为512*512像素；
步骤二，将图像进行对称处理，处理后图像的大小为原图像的4倍；图像边界区域的点在求导的时候会由于边界区域的有界而产生较大的偏差，将待处理的噪声图像翻转对称为1024*1024像素；计算公式如下：
forx＝1:M
fory＝1:N
u(x,2*N+1-y)＝u0(x,y)
u(2*M+1-x,y)＝u0(x,y)
u(2*M+1-x,2*N+1-y)＝u0(x,y)
end
end；
步骤三，利用MATLAB商业数学软件自带的快速傅立叶变换函数对图像进行傅立叶变换，公式如下：
步骤四，获得图像的分数阶微分；由于数字图像为离散信号，其微分只能用差分来近似，图像中，某一点对x的一阶偏导数可以写成差分形式；
计算图像的分数阶微分的计算过程如下：
由分数阶导数的表达式可以发现，分数阶导数的展开形式比较复杂，因此本发明采用快速Fourier变换的方法来避免繁琐的展开计算，对上式两边同时进行Fourier变换可得该点对x方向和y方向的α阶导数的Fourier变换：
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α，p(mi)表示在频域上对图像的第i行或者第i列的微分算子；
图像中，某一点对x方向的一阶偏导数可以写成差分形式：
相应的，根据分数阶导数的定义，某一点对x方向的α阶导数可以表示为：
对上式两边同时进行傅立叶变换可得该点对x方向和y方向的α阶导数的傅立叶变换，公式如下：
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α；
步骤五，获得图像分数阶梯度的模，计算过程如下：
将步骤四中频域上分数阶微分进行傅立叶反变换，得到空间域上的图像的分数阶微分，并将参数α换为扩散函数的参数β，
根据梯度的模的计算公式得到分数阶梯度的模：
步骤六，根据步骤五中求得的分数阶梯度的模计算扩散函数：
g(x)＝1/(1+x2)，
步骤七，计算微分算子的共轭算子，计算过程如下：
p*(m1)＝conj((1-exp(-j2πm1/N))α)，
输出一次迭代求解的结果：
其中：
步骤九，将步骤八中所得结果进行傅立叶反变换，步骤十，将un+1还原到原图像大小，计算过程如下：
forx＝1:M
fory＝1:N
end
end；
步骤十一，利用下式计算所得结果与原图像进行峰值信噪比计算，判断是否满足终止条件；
其中f为无噪声图像，u为一次迭代还原后的图像，若符合终止条件：ΔPSNR≤ε，其中，ΔPSNR表示两次迭代结果的峰值信噪比的变化绝对值，ε取0.01，则执行步骤十三；若ΔPSNR＞ε，则执行步骤十二；
步骤十二，重新迭代求解，将上一次的迭代结果作为下一次迭代的输入，继续执行步骤四至步骤十一；
步骤十三，输出还原后的结果则为最优的去噪图像。
本发明采用的软件是MATLAB；本发明需要的最低配置的处理器为：I3-2100处理器；本发明需要的最低内存是1GB RAM。
本发明的效果可以通过以下仿真实验进一步说明；
1.实验条件：
本发明的实验仿真环境为：
软件：MATLAB R2015a
处理器：Intel(R)Core(TM)i5-6500CPU@3.20GHz
内存：16GB RAM
本发明的实验所用到的图像来源于标准图像库。
2.实验内容：
本发明的实验具体分为三个实验。
实验一：利用本发明对含有标准差δ＝25的噪声的Lena图像进行去噪，结果如图2所示。
实验二：利用本发明对含有标准差δ＝25的噪声的Barbara图像进行去噪，结果如图3所示。
实验三：利用本发明对含有标准差δ＝25的噪声的Boat图像进行去噪，结果如图4所示。
实验中，使用峰值信噪比PSNR评价指标来评价去噪结果的优劣，其PSNR定义为：
其中，f表示原始无噪声图像，u表示去噪后的图像，图像大小为NxN。
在添加δ＝15,20,25不同噪声的情况下，将本发明与现有的基于偏微分方程的方法进行对比，结果如表1所示：
表1添加不同噪声时本方法与其他方法去噪效果对比结果
3.实验结果分析
从表1的实验数据可以发现，本发明与其他方法相比，PSNR有了明显的提升，主要原因是扩散函数g(x)在去噪过程中起了极为重要的最用，检测边缘区分图像的纹理区域和平坦区域，传统方法中用|▽u|作为边缘检测算子，但由于受噪声的影响会产生很多虚假边缘，而分数阶导数的非局部性可以减弱噪声的影响，即|Dβu|比|▽u|在噪声干扰的情况下具有更好的边缘检测准确度。另外，分数阶微分的非局部性也有效地抑制了阶梯效应，去噪之后图像的视觉效果也有所提升。
本发明提供一种基于分数阶偏微分方程的各向异性扩散的图像还原方法，本发明根据分数阶导数的非局部性质，在检测边缘时能够减弱噪声的干扰，结合偏微分方程得到一种基于分数阶偏微分方程的图像去噪方法，能够在去噪的同时尽可能地保留原图像的纹理细节。在求解的过程中采用了快速Fourier变换的方法，避免了复杂的分数阶导数展开运算的同时加快了求解速度。本发明将扩散函数的变量单独设定了分数阶导数，对于不同的图像变化微分阶数可以获得较好的去噪效果，并且收敛速度也较快，所需的迭代次数较少。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (8)

1.基于分数阶偏微分方程的图像去噪方法，其特征在于，包括如下步骤：
步骤一，输入一幅被噪声污染的图像u0(x,y)，图像大小为M×N，设定时间间隔；
步骤二，将图像进行对称处理，处理后图像的大小为原图像的4倍，计算公式如下：
步骤三，利用MATLAB商业数学软件自带的快速傅立叶变换函数对图像进行傅立叶变换，公式如下：
步骤四，获得图像的分数阶微分；
计算图像的分数阶微分的计算过程如下：
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&amp;ap;</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>2</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>6</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>3</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>...</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>n</mi> <mo>!</mo> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow>
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α，p(mi)表示在频域上对图像的第i行或者第i列的微分算子；
步骤五，获得图像分数阶梯度的模，计算过程如下：
将步骤四中频域上分数阶微分进行傅立叶反变换，得到空间域上的图像的分数阶微分，并将参数α换为扩散函数的参数β
<mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> <mo>=</mo> <mi>i</mi> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>p</mi> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>)</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>(</mo> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> <mo>=</mo> <mi>i</mi> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>p</mi> <mo>(</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>(</mo> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
根据梯度的模的计算公式得到分数阶梯度的模：
<mrow> <mo>|</mo> <msup> <mi>D</mi> <mi>&amp;beta;</mi> </msup> <mi>u</mi> <mo>|</mo> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
步骤六，根据步骤五中求得的分数阶梯度的模计算扩散函数：
g(x)＝1/(1+x2)，
<mrow> <mi>g</mi> <mrow> <mo>(</mo> <mo>|</mo> <msup> <mi>D</mi> <mi>&amp;beta;</mi> </msup> <mi>u</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>/</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤七，计算微分算子的共轭算子，计算过程如下：
p*(m1)＝conj((1-exp(-j2πm1/N))α)，
输出一次迭代求解的结果：
<mrow> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> <mo>=</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>n</mi> </msup> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>g</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>*</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>,</mo> </mrow>
其中：
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>c</mi> <mi>o</mi> <mi>n</mi> <mi>j</mi> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;m</mi> <mn>1</mn> </msub> <mo>/</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>&amp;alpha;</mi> </msup> <mo>)</mo> </mrow> <mo>*</mo> <mn>...</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>*</mo> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>i</mi> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mo>(</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;m</mi> <mn>1</mn> </msub> <mo>/</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>&amp;alpha;</mi> </msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>*</mo> <mn>1</mn> <mo>/</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>c</mi> <mi>o</mi> <mi>n</mi> <mi>j</mi> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;m</mi> <mn>2</mn> </msub> <mo>/</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>&amp;alpha;</mi> </msup> <mo>)</mo> </mrow> <mo>*</mo> <mn>...</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>*</mo> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>i</mi> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mo>(</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;m</mi> <mn>2</mn> </msub> <mo>/</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mi>&amp;alpha;</mi> </msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>*</mo> <mn>1</mn> <mo>/</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;beta;</mi> </msubsup> <mi>u</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
步骤九，将步骤八中所得结果进行傅立叶反变换，
步骤十，将un+1还原到原图像大小，计算过程如下：
步骤十一，利用下式计算所得结果与原图像进行峰值信噪比计算，判断是否满足终止条件；
<mrow> <mi>P</mi> <mi>S</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mn>10</mn> <mi>lg</mi> <mfrac> <mrow> <munder> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mn>0</mn> <mo>&amp;le;</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>&amp;le;</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munder> <mo>|</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> <mrow> <mfrac> <mn>1</mn> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mi>f</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> <mo>-</mo> <mi>u</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>;</mo> </mrow>
其中f为无噪声图像，u为一次迭代还原后的图像，若符合终止条件：ΔPSNR≤ε，其中，ΔPSNR表示两次迭代结果的峰值信噪比的变化绝对值，ε取0.01，则执行步骤十三；若ΔPSNR＞ε，则执行步骤十二；
步骤十二，重新迭代求解，将上一次的迭代结果作为下一次迭代的输入，继续执行步骤四至步骤十一；
步骤十三，输出还原后的结果则为最优的去噪图像。
2.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，步骤一，输入一幅被噪声污染的图像u0(x,y)，图像大小为M×N，设定时间间隔Δt≤0.25。
3.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，步骤四，获得图像的分数阶微分；
计算图像的分数阶微分的计算过程如下：
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&amp;ap;</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>2</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>6</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>3</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>...</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>n</mi> <mo>!</mo> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow>
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α，p(mi)表示在频域上对图像的第i行或者第i列的微分算子；
图像中，某一点对x方向的一阶偏导数可以写成差分形式：
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>u</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
相应的，根据分数阶导数的定义，某一点对x方向的α阶导数可以表示为：
<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&amp;ap;</mo> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>2</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>6</mn> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>3</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>...</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>n</mi> <mo>!</mo> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>u</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
对上式两边同时进行傅立叶变换可得该点对x方向和y方向的α阶导数的傅立叶变换，公式如下：
<mrow> <msup> <mi>u</mi> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> <mo>=</mo> <mi>i</mi> <mi>f</mi> <mi>f</mi> <mi>t</mi> <mn>2</mn> <mrow> <mo>(</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>x</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>D</mi> <mi>y</mi> <mi>&amp;alpha;</mi> </msubsup> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <msub> <mi>m</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>m</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中p(m1)＝(1-exp(-j2πm1/N))α，p(m2)＝(1-exp(-j2πm2/N))α。
4.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，本发明采用的软件是MATLAB。
5.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，本发明需要的最低配置的处理器为：I3-2100处理器。
6.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，本发明采用的处理器为：Intel(R)Core(TM)i5-6500CPU@3.20GHz处理器。
7.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，本发明需要的最低内存是1GB RAM。
8.根据权利要求1所述的基于分数阶偏微分方程的图像去噪方法，其特征在于，本发明采用的内存为16GB RAM。
CN201710456403.0A 2017-06-16 2017-06-16 基于分数阶偏微分方程的图像去噪方法 Active CN107085840B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710456403.0A CN107085840B (zh) 2017-06-16 2017-06-16 基于分数阶偏微分方程的图像去噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710456403.0A CN107085840B (zh) 2017-06-16 2017-06-16 基于分数阶偏微分方程的图像去噪方法

Publications (2)

Publication Number Publication Date
CN107085840A true CN107085840A (zh) 2017-08-22
CN107085840B CN107085840B (zh) 2019-09-27

Family

ID=59605503

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710456403.0A Active CN107085840B (zh) 2017-06-16 2017-06-16 基于分数阶偏微分方程的图像去噪方法

Country Status (1)

Country Link
CN (1) CN107085840B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596859A (zh) * 2018-05-10 2018-09-28 南京邮电大学 一种基于高阶偏微分方程的图像去噪方法
CN109171706A (zh) * 2018-09-30 2019-01-11 南京信息工程大学 基于分类匹配和分数阶扩散的心电信号去噪方法和系统
CN109271913A (zh) * 2018-09-05 2019-01-25 南京信息工程大学 一种基于偏微分方程的maldi质谱峰检测方法
CN111476743A (zh) * 2020-02-18 2020-07-31 烟台大学 一种基于分数阶微分的数字信号滤波与图像处理方法
CN112248886A (zh) * 2020-10-30 2021-01-22 西南交通大学 一种座椅自动调节方法
CN112819701A (zh) * 2019-11-15 2021-05-18 中国科学院长春光学精密机械与物理研究所 一种图像去噪方法、系统及电子设备
CN114549364A (zh) * 2022-03-03 2022-05-27 中国人民解放军国防科技大学 一种图像去噪方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887576A (zh) * 2010-06-04 2010-11-17 西安电子科技大学 基于偏微分方程滤波器的图像去噪方法
CN103337060A (zh) * 2013-07-17 2013-10-02 蒲亦非 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器
CN103427789A (zh) * 2013-07-23 2013-12-04 四川大学 一种基于分数阶计算方程的图书馆图文信息去噪滤波器
CN105930811A (zh) * 2016-04-26 2016-09-07 济南梦田商贸有限责任公司 基于图像处理的手掌纹理特征检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101887576A (zh) * 2010-06-04 2010-11-17 西安电子科技大学 基于偏微分方程滤波器的图像去噪方法
CN103337060A (zh) * 2013-07-17 2013-10-02 蒲亦非 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器
CN103427789A (zh) * 2013-07-23 2013-12-04 四川大学 一种基于分数阶计算方程的图书馆图文信息去噪滤波器
CN105930811A (zh) * 2016-04-26 2016-09-07 济南梦田商贸有限责任公司 基于图像处理的手掌纹理特征检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
杨家红 等: "结合USFFT Curvelet 变换的各向异性扩散图像去噪模型", 《通信学报》 *
杨迎春 等: "基于分数阶导数的自适应各向异性扩散图像去噪模型", 《中北大学学报(自然科学版)》 *
蒋伟: "基于分数阶偏微分方程的图像去噪新模型", 《计算机应用》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596859A (zh) * 2018-05-10 2018-09-28 南京邮电大学 一种基于高阶偏微分方程的图像去噪方法
CN109271913A (zh) * 2018-09-05 2019-01-25 南京信息工程大学 一种基于偏微分方程的maldi质谱峰检测方法
CN109271913B (zh) * 2018-09-05 2022-03-29 南京信息工程大学 一种基于偏微分方程的maldi质谱峰检测方法
CN109171706A (zh) * 2018-09-30 2019-01-11 南京信息工程大学 基于分类匹配和分数阶扩散的心电信号去噪方法和系统
CN109171706B (zh) * 2018-09-30 2020-12-29 南京信息工程大学 基于分类匹配和分数阶扩散的心电信号去噪方法和系统
CN112819701A (zh) * 2019-11-15 2021-05-18 中国科学院长春光学精密机械与物理研究所 一种图像去噪方法、系统及电子设备
CN111476743A (zh) * 2020-02-18 2020-07-31 烟台大学 一种基于分数阶微分的数字信号滤波与图像处理方法
CN111476743B (zh) * 2020-02-18 2023-06-06 烟台大学 一种基于分数阶微分的数字信号滤波与图像处理方法
CN112248886A (zh) * 2020-10-30 2021-01-22 西南交通大学 一种座椅自动调节方法
CN114549364A (zh) * 2022-03-03 2022-05-27 中国人民解放军国防科技大学 一种图像去噪方法及系统
CN114549364B (zh) * 2022-03-03 2024-04-09 中国人民解放军国防科技大学 一种图像去噪方法及系统

Also Published As

Publication number Publication date
CN107085840B (zh) 2019-09-27

Similar Documents

Publication Publication Date Title
CN107085840A (zh) 基于分数阶偏微分方程的图像去噪方法
CN107481264B (zh) 一种自适应尺度的视频目标跟踪方法
Wang et al. Fusion method for infrared and visible images by using non-negative sparse representation
Gottwald et al. Testing for chaos in deterministic systems with noise
CN102496157B (zh) 基于高斯多尺度变换及颜色复杂度的图像检测方法
US11880903B2 (en) Bayesian image denoising method based on distribution constraint of noisy images
Thai et al. Global variational method for fingerprint segmentation by three‐part decomposition
Gilboa Nonlinear eigenproblems in image processing and computer vision
CN104794685A (zh) 一种实现图像去噪的方法及装置
CN103559697A (zh) 基于fft的碎纸片纵切拼接复原算法
Li et al. Signal denoising with random refined orthogonal matching pursuit
Guryanov et al. Fast medical image registration using bidirectional empirical mode decomposition
CN105389839A (zh) 基于流体分析的流体参数估计方法
Tu et al. Logarithmic dyadic wavelet transform with its applications in edge detection and reconstruction
Dai et al. Salient object detection from multi-spectral remote sensing images with deep residual network
CN112381731B (zh) 一种基于图像去噪的单帧条纹图像相位分析方法及系统
Maji et al. Reconstructing an image from its edge representation
Hou et al. A variational formulation for fingerprint orientation modeling
He et al. Beyond ${\chi^ 2} $ Difference: Learning Optimal Metric for Boundary Detection
CN105303538A (zh) 一种基于nsct和pca的高斯噪声方差估计方法
Rew et al. Hybrid Segmentation Scheme for Skin Features Extraction Using Dermoscopy Images.
Martins et al. On the completeness of feature-driven maximally stable extremal regions
Xu et al. Variational total curvature model for multiplicative noise removal
Zhao et al. Small-target leak detection for a closed vessel via infrared image sequences
CN103679201B (zh) 一种用于图像匹配、识别、检索的点集合匹配的校正方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210210

Address after: 210012 19 / F, B1 building, Kechuang City, 34 Dazhou Road, Yuhuatai District, Nanjing City, Jiangsu Province

Patentee after: Jiangsu Xinshiyun Science and Technology Co.,Ltd.

Address before: 210044, No. 219, Ning six road, Pukou District, Jiangsu, Nanjing

Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY

TR01 Transfer of patent right