高斯卷积核如何生成 C语言实现

    对于学图像专业的人来说,对图像进行高斯滤波应该不会陌生,本质上就是将图像与高斯核进行卷积。

但是高斯核是如何生成的呢。matlab中有函数能自动生成高斯卷积核:

    gs=fspecial('gaussian',3,1)

    gs =

    0.0751    0.1238    0.0751
    0.1238    0.2042    0.1238
    0.0751    0.1238    0.0751

    算法公式:


C语言编程实现如下:

 

  void main()
{
 int i,j;
 double sigma=1;
 const int N=3;
 double gussian[N][N];
 double sum=0.0;

 for(i=0;i<N;i++)
 {
  for(j=0;j<N;j++)
  {
   gussian[i][j]=exp(-((i-N/2)*(i-N/2)+(j-N/2)*(j-N/2))/(2.0*sigma*sigma));
   sum+=gussian[i][j];
  }
 }

 for(i=0;i<N;i++)
 {
  for(j=0;j<N;j++)
  {
   gussian[i][j]=gussian[i][j]/sum;
   printf("%f  ,",gussian[i][j]);
  }
  printf("\n");
 }
}



 


转:http://blog.sina.com.cn/s/blog_7f2d8e4e0102wuxj.html

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在单片机的C语言实现高斯滤波可以通过以下步骤: 1. 定义高斯:根据需要选择合适的高斯大小和标准差。高斯是一个二维矩阵,用于计算滤波器的权重。 2. 加载图像数据:将需要进行滤波的图像数据加载到单片机的内存中。 3. 遍历图像像素:对于每个像素,计算其周围像素的加权平均值,权重由高斯决定。 4. 滤波计算:使用高斯计算像素的加权平均值。可以使用多种方法,如卷运算或逐像素计算。 5. 输出结果:将滤波后的结果保存到内存中或通过外设输出。 以下是一个简单的示例代码,用于在单片机上实现高斯滤波: ```c #include <stdio.h> #include <math.h> // 定义高斯大小和标准差 #define KERNEL_SIZE 3 #define SIGMA 1.0 // 高斯 const float kernel[KERNEL_SIZE][KERNEL_SIZE] = { {0.0625, 0.125, 0.0625}, {0.125, 0.25, 0.125}, {0.0625, 0.125, 0.0625} }; // 图像数据 const unsigned char image[8][8] = { { 10, 20, 30, 40, 50, 60, 70, 80 }, { 15, 25, 35, 45, 55, 65, 75, 85 }, { 20, 30, 40, 50, 60, 70, 80, 90 }, { 25, 35, 45, 55, 65, 75, 85, 95 }, { 30, 40, 50, 60, 70, 80, 90, 100 }, { 35, 45, 55, 65, 75, 85, 95, 105 }, { 40, 50, 60, 70, 80, 90, 100, 110 }, { 45, 55, 65, 75, 85, 95, 105,115 } }; // 高斯滤波函数 void gaussianFilter(const unsigned char input[][8], unsigned char output[][8]) { int i,,m,n; float sum; // 遍历图像像素 for(i=0; i<8; i++) { for(j=0; j<8; j++) { sum = 0.0; // 计算加权平均值 for(m=0; m<KERNEL_SIZE; m++) { for(n=0; n<KERNEL_SIZE; n++) { // 边界处理 int x = i + m - KERNEL_SIZE/2; int y = j + n - KERNEL_SIZE/2; if (x < 0) x = 0; if (x >= 8) x = 7; if (y < 0) y = 0; if (y >= 8) y = 7; sum += input[x][y] * kernel[m][n]; } } // 保存滤波结果 output[i][j] = (unsigned char)sum; } } } int main() { unsigned char filteredImage[8][8]; // 执行高斯滤波 gaussianFilter(image, filteredImage); // 打印滤波结果 int i,j; for(i=0; i<8; i++) { for(j=0; j<8; j++) { printf("%d ", filteredImage[i][j]); } printf("\n"); } return 0; } ``` 注意:上述代码只是一个简单的示例,并未考虑性能和边界处理等问题。实际应用中,可能需要对代码进行优化和改进。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值