图象增强--梯度锐化

作者:Sender  来源:----  发布日期:2000-02-13  最后修改日期:2008-02-12

1、微分法

  在图象中,边缘是由灰度级和相邻域点不同的象素点构成的。因而,若想增强边缘,就应该突出相邻点间的灰度级的变化。微分运算可用来求信号的变化率,因而具有加强高频分量的作用。如果将其应用在图象上,可使图象的轮廓清晰。由于我们常常无法事先确定轮廓的取向,因而挑选用于轮廓增强的微分算子时,必须选择那些不具备空间方向性的和具有旋转不变的线形微分算子。
  图象处理中最常用的微分方法是求梯度。对于图象f(x, y), 它在点(x,y)处的梯度是一个矢量。

  微分运算一般用差分来代替。常用的差分形式有两种:

1) GM(x,y)=|f(x,y)-f(x+1,y)|+|f(x,y)-f(x,y+1)|
2) GM(x,y)=|f(x,y)-f(x+1,y+1)|+|f(x+1,y)-f(x,y+1)|

  利用差分运算时,图象的最后一行和最后一列的象素的梯度无法求得,一般用前一行或前一列的梯度值近似代替。

  算出梯度后让梯度图象的灰度值g(x,y)等于该点的梯度幅度,即g(x,y)=GM(x,y)。这是常用的方法。

  还有就是:
  1)
      function 1
  2)
      function 2  Lg为一指定的灰度值。

  3)
      fucntion 3  Lb为一对背景指定的灰度值。

  4)
      function 4  Lg 和 Lb 的意义同上。

2、卷积

  一般可使用如下高通滤波矩阵:

0 -1 0
-1 5 -1
0 -1 0
-1 -1 -1
-1 9 -1
-1 -1 -1
1 -2 1
-2 5 -2
1 -2 1
-1 -2 -1
-2 19 -2
-1 -2 -1
-2 1 -2
1 6 1
-2 1 -2

下面是微分法第一种方法的程序,适用于灰阶图象。

/* Contents a_gradient Sharpen images with differential method */ #include <stdio.h> #include <vicdefs.h> #include <vicfcts.h> #include <vicerror.h> #include <string.h> #include <math.h> #include <float.h> extern int _cdecl checkrange_(imgdes *); /* Sharpen images. Returns NO_ERROR, BAD_RANGE, BAD_FAC, NO_EMM, EMM_ERR, NO_XMM, or XMM_ERR */ int _cdcel a_gradient(imgdes * srcimg, imgdes * desimg, int kind_method) { int pixel_gray_1, pixel_gray_2, pixel_gray_3, pixel_gray_4, j, k, i, l; int rcode=NO_ERROR; int sx, sy, ex, ey; int gradient; int kind; kind = kind_method; /* Check range of start, end position */ if (checkrange_(srcimg)) return (BAD_RANGE); if (kind != 1 && kind != 2) return (BAD_FAC); sx = srcimg->stx; sy = srcimg->sty; ex = srcimg->endx; ey = srcimg->endy; copyimgdes (srcimg, desimg); rcode = copyimage(srcimg, desimg); if (rcode != NO_ERROR) return (rcode); for (k= sx; k<= ex-1; k++) { for (j=sy; j<= ey-1; j++) { gradient = 0; if (kind ==1) { pixel_gray_1 = getpixelgray (srcimg, k, j); pixel_gray_2 = getpixelgray (srcimg, k+1, j); pixel_gray_3 = getpixelgray (srcimg, k, j); pixel_gray_4 = getpixelgray (srcimg, k, j+1); } else { pixel_gray_1 = getpixelgray (srcimg, k, j); pixel_gray_2 = getpixelgray (srcimg, k+1, j+1); pixel_gray_3 = getpixelgray (srcimg, k+1, j); pixel_gray_4 = getpixelgray (srcimg, k, j+1); } if (pixel_gray_1 < 0) return (pixel_gray_1); if (pixel_gray_2 < 0) return (pixel_gray_2); if (pixel_gray_3 < 0) return (pixel_gray_3); if (pixel_gray_4 < 0) return (pixel_gray_4); gradient = (int)(abs(pixel_gray_1 - pixel_gray_2)) + (int)(abs(pixel_gray_3 - pixel_gray_4)); if (gradient > 255) gradient = 255; rcode = setpixelgray(desimg, k, j, (UCHAR)gradient); if (rcode!=NO_ERROR) return (rcode); } } /* for last column and row */ for (j=sy; j<= ey-1; j++) { rcode=getpixelgray(desimg, ex-1, j); if (rcode!=NO_ERROR) return (rcode); rcode=setpixelgray(desimg, ex, j, rcode); if (rcode!=NO_ERROR) return (rcode); } for (k=sx; k<=sx-1; k++) { rcode = getpixelgray(desimg, k, ey-1); if (rcode!= NO_ERROR) return (rcode); rcode = setpixelgray (desimg, k, ey, rcode); if (rcode != NO_ERROR) return (rcode); } rcode = getpixelgray (desimg, ex-1, ey-1); if (rcode!= NO_ERROR) return (rcode); rcode = setpixelgray (desimg, ex, ey, rcode); if (rcode != NO_ERROR) return (rcode); }

  程序中用参数KIND来选择近似计算梯度的两种方法。

本栏目相关
  •  2000-01-20 亮度专业词语
  •  2000-02-13 图象增强--梯度锐化
  •  1999-10-23 Alpha通道的实现方法
  • 微信订阅号二维码

    本页网址二维码:

    本栏目热门内容
  • Acrobat虚拟PDF打印机执行打印时挂起,解决办法竟然...
  • LINKSYS交换机登录WEB界面显示不正确的解决方法
  • 又一次RAID 5阵列故障记录
  • 解决VMware vSphere ESXi 5.0 Update 1 中虚机不能...
  • 修改CentOS发行信息以绕过Dell服务器BIOS更新和DSET...
  • 解决虚拟化运行的 Windows Server 2003 标准版出现...
  • Windows Server 2008 重命名域和域控制器
  • Intel Nehalem CPU Errata 导致 VMWare ESXi(vSpher...
  • 一次很精神的电脑组装过程记录(但不是自己的电脑)...
  • 解决MySQL Cluster 备份总是失败,提示文件已存在的...
  • MegaCli安装及使用杂记
  • 解决WSUS显示客户端不全的问题
  • 解决 VMWare vSphere 6 客户端无法修改用户密码
  • 解决Windows Server 2008 R2域控制器显示无法连接到...
  • 本站服务器RAID 5阵列双硬盘失效挽救记录
  • 网站数据库从MySQL 5.0升级到5.6的记录
  • 解决MariaDB使用Percona XtraBackup增量备份出错
  • DELL PowerEdge 820 报CPU3 INTERNAL ERROR 的解决...
  • Linux 下的分区调整工具GParted实战
  • 修改arpwatch使通知邮件主题显示IP地址
  • 程序员漫画:如何用8种不同的编程语言去解救公主
  • 解决很好用的多合一即时通信软件pidgin的崩溃问题
  • 解决Samba WINS服务的错误解释问题
  • 使用 GParted 进行虚拟机硬盘分区调整操作
  • 解决Squid代理HTTP时在浏览器出现Content Encoding ...
  • 用Delphi编写使用到ADO的DLL的一些问题
  • 网站简单改版
  • 索尼系列手提电脑备份失败,出现700错误的解决办法
  • Dell R900服务器 BMC firmware incompatible with C...
  • 更多...