最近用到基于卷积操作的几种二维图像矩阵平滑的方法,简单实现了一下常见的高斯滤波、中值、均值滤波平滑算法
void GaussianFiltering(const Eigen::MatrixXf &src, Eigen::MatrixXf &dst, const int &ksize = 3, const float &sigma = 0.8)
{
//平滑是不处理边缘处
int b = ksize / 2;
//计算高斯核
Eigen::MatrixXf GauMat(ksize, ksize);
//float ksum = 0;
for (int i = 0; i < ksize; i++){
for (int j = 0; j < ksize; j++){
GauMat(i, j) = std::exp(-(std::pow((i - b), 2) + std::pow((j - b), 2)) / (2 * std::pow(sigma, 2)));
//ksum += GauMat(i, j);
}
}
GauMat /= GauMat.sum();
int rows = src.rows(), cols = src.cols();
for (int i = 0; i<rows; ++i)
for (int j = 0; j < cols; ++j) {
if ((i - b >= 0) && (j - b) >= 0 && (i + b)<rows && (j + b)<cols) {
if (src(i, j) >= 0){
Eigen::MatrixXf subGau = Eigen::MatrixXf::Zero(ksize, ksiz