CN108428216A - 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法 - Google Patents

基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法 Download PDF

Info

Publication number
CN108428216A
CN108428216A CN201810038949.9A CN201810038949A CN108428216A CN 108428216 A CN108428216 A CN 108428216A CN 201810038949 A CN201810038949 A CN 201810038949A CN 108428216 A CN108428216 A CN 108428216A
Authority
CN
China
Prior art keywords
formula
indicate
value
image
calculates
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
CN201810038949.9A
Other languages
English (en)
Other versions
CN108428216B (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.)
Liaoning Normal University
Original Assignee
Liaoning Normal University
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 Liaoning Normal University filed Critical Liaoning Normal University
Priority to CN201810038949.9A priority Critical patent/CN108428216B/zh
Publication of CN108428216A publication Critical patent/CN108428216A/zh
Application granted granted Critical
Publication of CN108428216B publication Critical patent/CN108428216B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开一种基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法,属于图像处理领域,首先,通过引入45N度旋转模值不变的算子建立扩散项的控制函数,使得扩散在平滑区域沿着梯度方向和垂直梯度方向进行,在边缘区域只沿着曲率方向进行;其次,通过散布矩阵对图像的特征区域有更好的分辨能力,并以此构造保真项改善了尖角区域和平滑区域的去噪效果,有效缓解了MCM模型在目标尖点处的过平滑效应和PM模型产生的“块状”效应。

Description

基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法
技术领域
本发明涉及图像处理领域,尤其是一种边缘和尖角保持能力强、平滑区域扩散效率高的基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法。
背景技术
遥感图像是利用被观测物体的辐射特性和电磁波的反射性所获取的一类特殊图像,目前已成为对地观测的重要手段,涉及国家政治、经济、军事和社会的众多领域并发挥着重要作用。然而,由于卫星遥感***往往处于复杂的电磁波环境,抑或由于载荷元器件间的电磁干扰或传感器的周期性偏移等因素,遥感图像在获取和传输过程中会不同程度地受到噪声影响,导致目标检测、目标鉴别等后续处理产生一定困难。因此,遥感图像去噪是一切后处理的重要基础之一。
目前,主流的遥感图像去噪方法主要可分为两类:基于非局部均值滤波的去噪方法和基于偏微分方程(Partial Differential Equations, PDE)扩散的去噪方法。尽管非局部均值滤波方法实现简单、无需多次迭代且去噪效果较好,可是该类方法易于在灰度变化缓慢的图像区域产生阶梯效应,而且需要在较大范围内为每个块搜索与其相似的块并进行高斯加权,计算复杂度高。基于PDE扩散的去噪方法主要采用多尺度分析理论设计自适应扩散模型,它又可进一步分为两类:基于各向同性扩散的去噪方法和基于各向异性扩散的去噪方法。考虑到各向同性扩散会随着尺度参数的增大而使图像边缘变得模糊,目前的相关研究和应用不多。而基于各向异性扩散的去噪方法则将图像去噪和边缘检测有机地结合起来,能够在平滑图像的同时较好地保持其边缘特性,取得了较高的主客观质量,典型方法如PM算法等。不过,该类方法对强边缘噪声和诸如椒盐噪声等孤立噪声的平滑效率却不够令人满意,而且众多学者研究发现PM方法是“病态”的,无法保证方程解的存在性和唯一性。为此,Catte等人通过对PM模型进行正则化,解决了方程解的唯一性,却又难以保证方差趋于0时的稳定性。为了更好地保持图像的边缘特征,Alvarez等人提出了基于平均曲率的MCM(Mean Curvature Model)扩散模型,使得扩散总是沿着垂直于梯度的方向进行,但由于该模型的扩散速度和曲率成正比,会在图像的尖角区域产生过度平滑现象。
发明内容
本发明是为了解决现有技术所存在的上述技术问题,提供一种边缘和尖角保持能力强、平滑区域扩散效率高的基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法。
本发明的技术解决方案是:一种基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法,其特征在于按如下步骤进行:
步骤1. 建立遥感图像去噪的二阶偏微分扩散模型,其定义如公式(1)所示:
(1)
所述表示初始图像,表示时间尺度下的平滑图像在坐标处的像素值,表示垂直于梯度矢量方向上的2阶导数,表示平行于梯度方向上的2阶导数,表示用来控制扩散方向的函数,表示用来控制尖角区域和平滑区域的保真度函数,为单调非增的函数且在[0,1]之间取值,为保真项,表示初始图像与时间尺度下的平滑图像在坐标处的像素值之差;
步骤2. 采用中心差分法将公式(1)表示为显式差分格式,其定义由公式(2)给出:
(2)
所述上标均表示迭代次数,下标表示某个像素的空间坐标,表示时间步长;
步骤3. 输入待去噪的初始图像,令
步骤4. 计算第次迭代图像的水平梯度和竖直梯度,其定义由公式(3)和公式(4)给出:
(3)
(4)
所述“*”表示2D卷积操作;
步骤5. 采用公式(5)计算第次迭代图像在每个像素处的梯度模值:
(5)
所述表示梯度算子;
步骤6. 计算函数在第次迭代的值,其定义由公式(6)-公式(7)给出:
(6)
(7)
所述表示取中值的函数;
步骤7. 计算散布矩阵,其定义由公式(8)给出:
(8)
所述表示标准差为的高斯卷积核,“”表示张量积,表示在标准差为的高斯平滑后的结果图像;
步骤8. 计算散布矩阵的特征值,其定义由公式(9)和公式(10)给出:
(9)
(10)
步骤9. 计算函数的值,其定义由公式(11)给出:
(11)
所述是一个很小的正数,其目的是避免分母为0;
步骤10. 根据公式(12)计算
(12)
所述分别表示第次迭代的图像沿着水平方向的2阶偏导数、沿着竖直方向的2阶偏导数和混合偏导数,数值计算采用中心差分格式;
步骤11. 根据公式(13)计算
(13)
步骤12. 将代入公式(2)计算并检查迭代过程是否稳定收敛,若稳定收敛,则停止迭代,输出图像,算法结束;否则,令,转入步骤4。
与现有技术相比,本发明具有以下优点:第一,沿着边缘方向和梯度方向进行不同程度的扩散,提高了在平滑区域的扩散效率,而且能较好地保持图像边缘;第二,扩散项第2项的控制函数选用了更稳定的45N度旋转模值不变的算子模版计算梯度模值,使得对梯度的估计更准确、扩散更稳定;第三,扩散项第3项利用了具有更丰富结构信息的散布矩阵,能够更加准确地分辨平滑区域、边缘区域和尖角区域,实现对尖角区域和平滑区域的保真,一方面避免了MCM模型在去噪同时容易模糊目标尖角的现象,另一方面克服了PM模型所产生的“块状”效应。
附图说明
图1为本发明实施例与其它方法的去噪效果对比图。
具体实施方式
本发明所提供的基于非下采样小波变换的改进四阶偏微分图像去噪方法,按照如下步骤进行:
步骤1. 建立遥感图像去噪的二阶偏微分扩散模型,其定义如公式(1)所示:
(1)
所述表示初始图像,表示时间尺度下的平滑图像在坐标处的像素值,表示垂直于梯度矢量方向上的2阶导数,表示平行于梯度方向上的2阶导数,表示用来控制扩散方向的函数,表示用来控制尖角区域和平滑区域的保真度函数,为单调非增的函数且在[0,1]之间取值,为保真项,表示初始图像与时间尺度下的平滑图像的像素值之差;
步骤2. 采用中心差分法将公式(1)表示为显式差分格式,其定义由公式(2)给出:
(2)
所述上标均表示迭代次数,下标表示某个像素的空间坐标,表示时间步长(本实施例取);
步骤3. 输入待去噪的初始图像,令
步骤4. 计算第次迭代图像的水平梯度和竖直梯度,其定义由公式(3)和公式(4)给出:
(3)
(4)
所述“*”表示2D卷积操作;
步骤5. 采用公式(5)计算第次迭代图像在每个像素处的梯度模值:
(5)
所述表示梯度算子;
步骤6. 计算函数在第次迭代的值,其定义由公式(6)-公式(7)给出:
(6)
(7)
所述表示取中值的函数;
步骤7. 计算散布矩阵,其定义由公式(8)给出:
(8)
所述表示标准差为的高斯卷积核,“”表示张量积,表示在标准差为的高斯平滑后的结果图像;
步骤8. 计算散布矩阵的特征值,其定义由公式(9)和公式(10)给出:
(9)
(10)
步骤9. 计算函数的值,其定义由公式(11)给出:
(11)
所述是一个很小的正数,其目的是避免分母为0;
步骤10. 根据公式(12)计算
(12)
所述分别表示第次迭代的图像沿着水平方向的2阶偏导数、沿着竖直方向的2阶偏导数和混合偏导数,数值计算采用中心差分格式;
步骤11. 根据公式(13)计算
(13)
步骤12. 将代入公式(2)计算并检查迭代过程是否稳定收敛,若稳定收敛,则停止迭代,输出图像,算法结束;否则,令,转入步骤4。
本发明实施例与其它方法的去噪效果比较如图1所示:从左至右分别为(a)原始图像;(b)含噪声图像;(c) PM方法的去噪结果;(d)MCM方法的去噪结果;(e)本发明实施例的去噪结果。
本发明实施例与其它方法的去噪结果(图1)的峰值信噪比对比如表1所示。
图1由上至下分别为图像1~4。
表1
对比结果表明:本发明方法的主客观去噪质量均优于其它方法。

Claims (1)

1.一种基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法,其特征在于按如下步骤进行:
步骤1. 建立遥感图像去噪的二阶偏微分扩散模型,其定义如公式(1)所示:
(1)
所述表示初始图像,表示时间尺度下的平滑图像在坐标处的像素值,表示垂直于梯度矢量方向上的2阶导数,表示平行于梯度方向上的2阶导数,表示用来控制扩散方向的函数,表示用来控制尖角区域和平滑区域的保真度函数,为单调非增的函数且在[0, 1]之间取值,为保真项,表示初始图像与时间尺度下的平滑图像在坐标处的像素值之差;
步骤2. 采用中心差分法将公式(1)表示为显式差分格式,其定义由公式(2)给出:
(2)
所述上标均表示迭代次数,下标表示某个像素的空间坐标,表示时间步长;
步骤3. 输入待去噪的初始图像,令
步骤4. 计算第次迭代图像的水平梯度和竖直梯度,其定义由公式(3)和公式(4)给出:
(3)
(4)
所述“*”表示2D卷积操作;
步骤5. 采用公式(5)计算第次迭代图像在每个像素处的梯度模值:
(5)
所述表示梯度算子;
步骤6. 计算函数在第次迭代的值,其定义由公式(6)和公式(7)给出:
(6)
(7)
所述表示取中值的函数;
步骤7. 计算散布矩阵,其定义由公式(8)给出:
(8)
所述表示标准差为的高斯卷积核,“”表示张量积,表示在标准差为的高斯平滑后的结果图像;
步骤8. 计算散布矩阵的特征值,其定义由公式(9)和公式(10)给出:
(9)
(10)
步骤9. 计算函数的值,其定义由公式(11)给出:
(11)
所述是一个很小的正数,其目的是避免分母为0;
步骤10. 根据公式(12)计算
(12)
所述分别表示第次迭代的图像沿着水平方向的2阶偏导数、沿着竖直方向的2阶偏导数和混合偏导数,数值计算采用中心差分格式;
步骤11. 根据公式(13)计算
(13)
步骤12. 将代入公式(2)计算并检查迭代过程是否稳定收敛,若稳定收敛,则停止迭代,输出图像,算法结束;否则,令,转入步骤4。
CN201810038949.9A 2018-01-16 2018-01-16 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法 Expired - Fee Related CN108428216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810038949.9A CN108428216B (zh) 2018-01-16 2018-01-16 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810038949.9A CN108428216B (zh) 2018-01-16 2018-01-16 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法

Publications (2)

Publication Number Publication Date
CN108428216A true CN108428216A (zh) 2018-08-21
CN108428216B CN108428216B (zh) 2021-09-24

Family

ID=63155903

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810038949.9A Expired - Fee Related CN108428216B (zh) 2018-01-16 2018-01-16 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法

Country Status (1)

Country Link
CN (1) CN108428216B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112862729A (zh) * 2021-03-30 2021-05-28 湖北地信科技集团股份有限公司 一种基于特征曲线指导的遥感图像去噪方法
CN117036962A (zh) * 2023-10-08 2023-11-10 中国科学院空天信息创新研究院 遥感图像变化检测方法、装置、设备及存储介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120106815A1 (en) * 2010-10-28 2012-05-03 Toshiba Medical Systems Corporation Denoising method and system for preserving clinically significant structures in reconstructed images using adaptively weighted anisotropic diffusion filter
CN103559693A (zh) * 2013-11-18 2014-02-05 东南大学 一种基于非连续性指示符的图像局部结构自适应复原方法
CN104346786A (zh) * 2014-11-06 2015-02-11 南京信息工程大学 基于Demons算法的图像去噪算法
CN104463811A (zh) * 2014-12-29 2015-03-25 南京信息工程大学 基于能量泛函的图像平滑与锐化算法
CN104463810A (zh) * 2014-12-25 2015-03-25 南京信息工程大学 基于tv流的自适应扩散滤波图像去噪算法
CN104517266A (zh) * 2014-12-22 2015-04-15 南京信息工程大学 基于边缘检测算子的混合自适应图像去噪方法
CN104899828A (zh) * 2015-04-08 2015-09-09 天津工业大学 一种基于角点保护的偏微分方程图像放大方法
CN105427262A (zh) * 2015-12-15 2016-03-23 南京信息工程大学 基于双向增强扩散滤波的图像去噪方法
CN105894462A (zh) * 2016-03-21 2016-08-24 辽宁师范大学 基于剪切波域隐马尔可夫树模型的遥感图像去噪方法
CN105894472A (zh) * 2016-04-14 2016-08-24 辽宁师范大学 基于非下采样小波变换的改进四阶偏微分图像去噪方法
CN106096517A (zh) * 2016-06-01 2016-11-09 北京联合大学 一种基于低秩矩阵与特征脸的人脸识别方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120106815A1 (en) * 2010-10-28 2012-05-03 Toshiba Medical Systems Corporation Denoising method and system for preserving clinically significant structures in reconstructed images using adaptively weighted anisotropic diffusion filter
CN103559693A (zh) * 2013-11-18 2014-02-05 东南大学 一种基于非连续性指示符的图像局部结构自适应复原方法
CN104346786A (zh) * 2014-11-06 2015-02-11 南京信息工程大学 基于Demons算法的图像去噪算法
CN104517266A (zh) * 2014-12-22 2015-04-15 南京信息工程大学 基于边缘检测算子的混合自适应图像去噪方法
CN104463810A (zh) * 2014-12-25 2015-03-25 南京信息工程大学 基于tv流的自适应扩散滤波图像去噪算法
CN104463811A (zh) * 2014-12-29 2015-03-25 南京信息工程大学 基于能量泛函的图像平滑与锐化算法
CN104899828A (zh) * 2015-04-08 2015-09-09 天津工业大学 一种基于角点保护的偏微分方程图像放大方法
CN105427262A (zh) * 2015-12-15 2016-03-23 南京信息工程大学 基于双向增强扩散滤波的图像去噪方法
CN105894462A (zh) * 2016-03-21 2016-08-24 辽宁师范大学 基于剪切波域隐马尔可夫树模型的遥感图像去噪方法
CN105894472A (zh) * 2016-04-14 2016-08-24 辽宁师范大学 基于非下采样小波变换的改进四阶偏微分图像去噪方法
CN106096517A (zh) * 2016-06-01 2016-11-09 北京联合大学 一种基于低秩矩阵与特征脸的人脸识别方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MICHAEL WENIGER等: "《Spatial Verification Using Wavelet Transforms》", 《QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY》 *
李睿: "《基于45N度旋转梯度等模算子偏微分方程去噪模型》", 《电脑编程技巧与维护》 *
杨平先 等: "《一种基于散布矩阵的自适应非局部均值滤波模型》", 《重庆邮电大学学报( 自然科学版)》 *
王相海 等: "《基于Cauchy分布模型与NSST变换的图像去噪算法》", 《辽宁师范大学学报(自然科学版)》 *
郑满满: "《改进的各向异性扩散图像去噪模型》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112862729A (zh) * 2021-03-30 2021-05-28 湖北地信科技集团股份有限公司 一种基于特征曲线指导的遥感图像去噪方法
CN112862729B (zh) * 2021-03-30 2022-09-06 湖北地信科技集团股份有限公司 一种基于特征曲线指导的遥感图像去噪方法
CN117036962A (zh) * 2023-10-08 2023-11-10 中国科学院空天信息创新研究院 遥感图像变化检测方法、装置、设备及存储介质
CN117036962B (zh) * 2023-10-08 2024-02-06 中国科学院空天信息创新研究院 遥感图像变化检测方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN108428216B (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN109767439B (zh) 一种自适应窗口的多尺度差异与双边滤波的目标检测方法
CN112115953B (zh) 一种基于rgb-d相机结合平面检测与随机抽样一致算法的优化orb算法
Wang et al. Noise detection and image denoising based on fractional calculus
CN109242972B (zh) 一种基于顶点特征的双法向网格模型光顺方法
CN102800063B (zh) 基于各向异性滤波的图像增强与抽象化方法
CN110796616B (zh) 基于范数约束和自适应加权梯度的湍流退化图像恢复方法
CN107967675B (zh) 一种基于自适应投影移动最小二乘的结构化点云去噪方法
CN108010002B (zh) 一种基于自适应隐式移动最小二乘的结构化点云去噪方法
CN108428216A (zh) 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法
CN108932699A (zh) 基于变换域的三维匹配调和滤波图像去噪方法
CN110363775A (zh) 一种基于区域型变分水平集的图像分割方法
Narváez et al. Point cloud denoising using robust principal component analysis
Li et al. Adaptive weighted guided image filtering for depth enhancement in shape-from-focus
Wang et al. Robust nonstationary local slope estimation
CN111880236B (zh) 一种构建多层等效源模型计算化极与数据类型转换的方法
Lott et al. A transition from Kelvin–Helmholtz instabilities to propagating wave instabilities
Leng et al. A novel geometric flow approach for quality improvement of multi-component tetrahedral meshes
CN113421205B (zh) 一种结合红外偏振成像的小目标检测方法
TW202316367A (zh) 決定去模糊化濾波器係數的方法
Zhang et al. Monocular vision SLAM research for parking environment with low light
CN110796609B (zh) 基于尺度感知和细节增强模型的低光图像增强方法
Xiang et al. Fast and simple super resolution for range data
Xiang et al. Real-time spatial and depth upsampling for range data
Hu et al. Denosing 3D ultrasound images by non-local means accelerated by GPU
Liu et al. Research on influence model of fractional calculus derivative Fourier transform on aerobics view image reconstruction

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210924

CF01 Termination of patent right due to non-payment of annual fee