图像的高斯滤波

最新推荐文章于 2023-08-08 20:22:20 发布

darkerJ

最新推荐文章于 2023-08-08 20:22:20 发布

阅读量537

点赞数 1

分类专栏：图像处理文章标签：图像处理

本文链接：https://blog.csdn.net/jsf921942722/article/details/51526987

版权

图像处理专栏收录该内容

7 篇文章 0 订阅

订阅专栏

MakeGaussianData.h MakeGaussianData.c 是产生高斯滤波的数据的头文件、函数
MakeGaussianData.h

#ifndef _MAKEGAUSSIANDATA_H_
#define _MAKEGAUSSIANDATA_H_

#ifdef __cplusplus
extern "C"
{
#endif
    int makeGaussianData(float ** data, int width, int height, float filterRatio);
#ifdef __cplusplus
}
#endif

#endif

MakeGaussianData.c

#include <stdio.h>
#include <math.h>
#include <memory.h>
#include <malloc.h>
#define ENumber 2.71828183f
#define PI 3.1415926f
int makeGaussianData(float ** data, int width, int height, float filterRatio){
    int i, j;
    *data = (float*)malloc(sizeof(float) * width * height);
    for(i = 0; i < height; i ++){
        for(j = 0; j < width; j ++){
            float temp = 0.0f;
            temp = -1.0f * ((j - width / 2) * (j - width / 2) + (i - height / 2) * (i - height / 2)) / (2.0f * filterRatio * filterRatio);
            temp = powf(ENumber, temp);
            (*data)[i * width + j] = (1.0f / (2.0f * PI * filterRatio * filterRatio)) * temp;
        }
    }
    return 0;
}

2.GaussianFlur.h GaussianFlur.c是高斯滤波的头文件、核心函数
GaussianFlur.h

#ifndef _GAUSSIANFLUR_H_
#define _GAUSSIANFLUR_H_

#ifdef __cplusplus
extern "C"
{
#endif
    int gaussianFlurFunction(float * inData, float * outData, int dataWidth, int dataHeight, int gaussianWidth, int gaussianHeight, float filterRatio, int channels);
#ifdef __cplusplus
}
#endif

#endif

GaussianFlur.c

#include <stdio.h>
#include <stdlib.h>
#include "MakeGaussianData.h"
#include "GaussianFlur.h"
int gaussianFlurFunction(float * inData, float * outData, int dataWidth, int dataHeight, int gaussianWidth, int gaussianHeight, float filterRatio, int channels){
    float * data = NULL;
    int i, j, z;
    int m, n;
    int widthTemp, heightTemp;
    float sumNumber = 0.0f;
    if(channels > 32){
        printf("channels is larger.\n");
        return -1;
    }

    if(gaussianWidth % 2 == 0){
        widthTemp = gaussianWidth / 2 - 1;
    }else{
        widthTemp = gaussianWidth / 2;
    }

    if(gaussianHeight % 2 == 0){
        heightTemp = gaussianHeight / 2 - 1;
    }else{
        heightTemp = gaussianHeight / 2;
    }

    makeGaussianData(&data, gaussianWidth, gaussianHeight, filterRatio);
    for(z = 0; z < channels; z ++){
        for(i = 0; i < dataHeight; i ++){
            for(j = 0; j < dataWidth; j ++){
                sumNumber = 0.0f;

                for(m = -heightTemp; m < gaussianHeight - heightTemp; m ++){
                    for(n = -widthTemp; n < gaussianWidth - widthTemp; n ++){
                        if(i + m >= 0 && j + n >= 0 && i + m < dataHeight && j + n < dataWidth){
                            sumNumber += inData[((i + m) * dataWidth + (j + n)) * channels + z];
                        }
                    }
                }
                outData[(i * dataWidth + j) * channels + z] = sumNumber;
            }
        }
    }
    free(data);
    return 0;
}

3.主函数

#include "stdafx.h"
#include "MakeGaussianData.h"
#include "Tools.h"
#include "GaussianFlur.h"
#include <opencv2\opencv.hpp>
using namespace cv;
int main()
{
    int imageWidth, imageHeight;
    int channels = 1;
    Mat mImage = imread("D:\\workSpace\\VSWorkSpace\\ImageProcess\\ImageProcess\\gaussianImage.jpg", 1);
    imshow("srcImage", mImage);
    float * imageData = NULL;
    float * imageOutData = NULL;
    imageWidth = mImage.cols;
    imageHeight = mImage.rows;
    channels = mImage.channels();
    imageData = (float*)malloc(sizeof(float) * imageWidth * imageHeight * channels);
    imageOutData = (float*)malloc(sizeof(float) * imageWidth * imageHeight * channels);
    uCharDataToFloatData(mImage.data, imageData, imageWidth, imageHeight, channels);

    gaussianFlurFunction(imageData, imageOutData, imageWidth, imageHeight, 5, 5, 0.1f, channels);

    floatDataToUCharData(imageOutData, mImage.data, imageWidth, imageHeight, channels, 1);
    imshow("gaussianImage", mImage);
    waitKey(0);
    free(imageData);
    free(imageOutData);
    return 0;
}

uCharDataToFloatData();
floatDataToUCharData();两个函数见链接

http://blog.csdn.net/jsf921942722/article/details/51526673
4.效果图
*********************************3通道高斯滤波彩色图原图***************************

*****************************3通道彩色高斯滤波后效果图******************

*********************************1通道高斯滤波灰度图原图***************************

*****************************1通道灰度高斯滤波后效果图******************