生成高斯核

// makeGause.cpp : 定义控制台应用程序的入口点。 // #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; float sum = 0.0f; *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; sum += (*data)[i * width + j]; } } printf("sum = %f\n", sum); for(i = 0; i < height; i ++){ for(j = 0; j < width; j ++){ (*data)[i * width + j] = (*data)[i * width + j] / sum; } } return 0; } int makeOneDimGaussKernel(float ** data, int num, float sigma) { int i = 0; *data = (float*)malloc(sizeof(float) * num); float sum = 0.0f; for(i = 0; i < num; i ++) { float temp = 0.0f; temp = -1.0f * (i - num / 2) * (i - num / 2) / (2 * sigma * sigma); temp = powf(ENumber, temp); (*data)[i] = (1.0f / (2.0f * PI * sigma * sigma)) * temp; sum += (*data)[i]; } printf("onesum = %f\n", sum); for(i = 0; i < num; i ++) { (*data)[i] = (*data)[i] / sum; } return 0; } int main(int argc, char* argv[]) { float * gaussData = NULL; int width = 5; int height = 5; float sigma = 0.8f; makeGaussianData(&gaussData, width, height, sigma); for(int i = 0; i < height; i ++) { for(int j = 0; j < width; j ++) { printf("%f ", gaussData[width * i + j]); } printf("\n"); } float * dataOne = NULL; makeOneDimGaussKernel(&dataOne, width, sigma); printf("\n"); for(int i = 0; i < width; i ++) { printf("%f ", dataOne[i]); } printf("\n"); free(dataOne); free(gaussData); return 0; }