Matlab中del2()函数学习笔记

最新推荐文章于 2024-09-29 15:36:55 发布

hxgqh

最新推荐文章于 2024-09-29 15:36:55 发布

阅读量1.6w

点赞数 2

分类专栏： Matlab与C/C++混合编程文章标签： matlab dst ui 编程 c

本文链接：https://blog.csdn.net/hxgqh/article/details/6253029

版权

Matlab与C/C++混合编程专栏收录该内容

2 篇文章 0 订阅

订阅专栏

在L=del2(U)表达式中，Matlab帮助文件的解释是：l_i,j=(u_i+1,j+u_i-1,j+u_i,j+1+u_i,j-1)/4 - u_i,j

但帮助文件没有解释边缘点如何处理。根据数学中的原理，在程序中做适当变化处理：

1. 数组左上角角点的公式应为：l_i,j=(u_i+2,j+u_i+2,j+2u_i,j)/4 - (u_i+1,j+u_i,j+1)/2

2. 数组左边边点公式应为：l_i,j=(u_i+1,j+u_i-1,j+u_i,j+2+u_i,j)/4 - (u_i,j+u_i,j+1)/2

依这两个公式类推，可以得到其他角点和其他边点的计算公式，从而得到计算结果。用C编程实现del2时需要注意这些点的计算。

在OpenCV中，cvLaplace的功能能够实现del2。但不知道为什么使用的时候会出现问题，通不过。所以自己写了一个cvDel2。可能是画蛇添足，不过自己写的能用就可以了，写起来也不是太麻烦。

void cvDel2(IplImage *src,IplImage *dst)
{
for(i=0;i<src->height;i++)
{
  for(j=0;j<src->width;j++)
  {
   if(i==0 && j==0)
   {
    cvSetReal1D(dst,i*src->width+j,(2*cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i+2)*src->width+j)+cvGetReal1D(src,i*src->width+j+2))/4-(cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,i*src->width+j+1))/2);
   }
   else if(i==0 && j!=0 && j!=src->width-1)
   {
    cvSetReal1D(dst,i*src->width+j,(cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i+2)*src->width+j)+cvGetReal1D(src,i*src->width+j+1)+cvGetReal1D(src,i*src->width+j-1))/4-(cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,i*src->width+j))/2);
   }
   else if(i==0 && j==src->width-1)
   {
    cvSetReal1D(dst,i*src->width+j,(2*cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i+2)*src->width+j)+cvGetReal1D(src,i*src->width+j-2))/4-(cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,i*src->width+j-1))/2);
   }
   else if(i!=0 && i!=src->height-1 && j==src->width-1)
   {
    cvSetReal1D(dst,i*src->width+j,(cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,i*src->width+j-2)+cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,(i-1)*src->width+j))/4-(cvGetReal1D(src,i*src->width+j-1)+cvGetReal1D(src,i*src->width+j))/2);
   }
   else if(i==src->height-1 && j==src->width-1)
   {
    cvSetReal1D(dst,i*src->width+j,(2*cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i-2)*src->width+j)+cvGetReal1D(src,i*src->width+j-2))/4-(cvGetReal1D(src,(i-1)*src->width+j)+cvGetReal1D(src,i*src->width+j-1))/2);
   }
   else if(i==src->height-1 && j!=0 && j!=src->width-1)
   {
    cvSetReal1D(dst,i*src->width+j,(cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i-2)*src->width+j)+cvGetReal1D(src,i*src->width+j+1)+cvGetReal1D(src,i*src->width+j-1))/4-(cvGetReal1D(src,(i-1)*src->width+j)+cvGetReal1D(src,i*src->width+j))/2);
   }
   else if(i==src->height-1 && j==0)
   {
    cvSetReal1D(dst,i*src->width+j,(2*cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,(i-2)*src->width+j)+cvGetReal1D(src,i*src->width+j+2))/4-(cvGetReal1D(src,(i-1)*src->width+j)+cvGetReal1D(src,i*src->width+j+1))/2);
   }
   else if(i!=src->height-1 && i!=0 && j==0)
   {
    cvSetReal1D(dst,i*src->width+j,(cvGetReal1D(src,(i-1)*src->width+j)+cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,i*src->width+j+2))/4-(cvGetReal1D(src,i*src->width+j)+cvGetReal1D(src,i*src->width+j+1))/2);
   }
   else
   {
    cvSetReal1D(dst,i*src->width+j,(cvGetReal1D(src,(i-1)*src->width+j)+cvGetReal1D(src,(i+1)*src->width+j)+cvGetReal1D(src,i*src->width+j-1)+cvGetReal1D(src,i*src->width+j+1))/4-cvGetReal1D(src,i*src->width+j));
   }
  }
}
}