【无标题】

oo372

于 2024-07-15 11:24:35 发布

阅读量69

点赞数 1

文章标签： opencv 人工智能计算机视觉

本文链接：https://blog.csdn.net/weixin_69483988/article/details/140430369

版权

方框滤波

#include<opencv2/opencv.hpp>
#include<quickopencv.h>
#include<iostream>
#include<math.h>
using namespace cv;
using namespace std;


int main(int argc, char** argv) {
	Mat equalLena = imread("D:/images/dog.png",IMREAD_ANYDEPTH);//输入图片
	if (equalLena.empty()) {
		printf("could not load image...\n");
		return -1;
	}
	imshow("原图", equalLena);

	//验证方框滤波算法的数据矩阵
	float points[25] = { 1,2,3,4,5,
						6,7,8,9,10,
						11,12,13,14,15,
						16,17,18,19,20,
						21,22,23,24,25 };
	Mat data(5, 5, CV_32FC1, points);
	//将CV_8U转变成CV_32F
	Mat equalLena_32f;
	equalLena.convertTo(equalLena_32f, CV_32F, 1.0 / 255);
	Mat resultNorm, result, dataSqrNorm, dataSqr, equalLena_32FSqr;
	//方框滤波函数boxFilter(),sqrBoxFilter()
	boxFilter(equalLena,result,-1,Size(3,3),Point(-1,-1),true);		//不进行归一化
	boxFilter(equalLena, resultNorm, -1, Size(3, 3), Point(-1, -1), false);		//进行归一化
	sqrBoxFilter(data,dataSqrNorm,-1,Size(3,3),Point(-1,-1),true,BORDER_CONSTANT);//进行归一化
	sqrBoxFilter(data, dataSqr, -1, Size(3, 3), Point(-1, -1), false, BORDER_CONSTANT);//不进行归一化
	sqrBoxFilter(equalLena_32f, equalLena_32FSqr, -1, Size(3, 3), Point(-1, -1), true, BORDER_CONSTANT);
	//显示处理结果
	imshow("resultNorm",resultNorm);
	imshow("result", result);
	imshow("equalLena_32FSqr", equalLena_32FSqr);

	waitKey(0);//此时图片显示时间为一直停留。（x）为x毫秒
	destroyAllWindows();
	return 0;
}

可分离滤波

#include<opencv2/opencv.hpp>
#include<quickopencv.h>
#include<iostream>
#include<math.h>
using namespace cv;
using namespace std;


int main(int argc, char** argv) {

	system("color F0");
	float points[25] = { 1,2,3,4,5,
						6,7,8,9,10,
						11,12,13,14,15,
						16,17,18,19,20,
						21,22,23,24,25 };
	Mat data(5, 5, CV_32FC1, points);
	//x方向、y方向和联合滤波器的构建
	Mat a = (Mat_<float>(3, 1) << -1, 3, -1);
	Mat b = a.reshape(1, 1);
	Mat ab = a * b;
	//验证高斯滤波的可分离性
	Mat gaussX = getGaussianKernel(3, 1);
	Mat gaussData, gaussDataXY;
	GaussianBlur(data, gaussData, Size(3, 3), 1, 1, BORDER_CONSTANT);
	sepFilter2D(data, gaussDataXY, -1, gaussX, gaussX, Point(-1, -1), 0, BORDER_CONSTANT);
	//输入高斯滤波的计算结果
	cout << "gaussData" << endl << gaussData << endl;
	cout << "gaussDataXY" << endl << gaussDataXY << endl;
	//线性滤波的可分离性
	Mat dataYX, dataY, dataXY, dataXY_sep;
	filter2D(data, dataY, -1, a, Point(-1, -1), 0, BORDER_CONSTANT);
	filter2D(dataY, dataYX, -1, b, Point(-1, -1), 0, BORDER_CONSTANT);
	filter2D(data, dataXY, -1, ab, Point(-1, -1), 0, BORDER_CONSTANT);
	sepFilter2D(data, dataXY_sep, -1, b, b, Point(-1, -1), 0, BORDER_CONSTANT);
	//输出可分离滤波和联合滤波的计算结果
	cout<<"dataY=" << endl << dataY << endl;
	cout << "dataYX=" << endl << dataYX << endl;
	cout << "dataXY=" << endl << dataXY << endl;
	cout << "dataXY_sep=" << endl << dataXY_sep << endl;

	Mat src = imread("D:/images/dog.png");//输入图片
	if (src.empty()) {
		printf("could not load image...\n");
		return -1;
	}
	//namedWindow("输入窗口", WINDOW_FREERATIO);//调整图片显示窗口大小；后面参数可调整-输出灰度图等
	imshow("原图", src);
	Mat srcY, srcXY, srcYX;
	filter2D(src, srcY, -1, a, Point(-1, -1), 0, BORDER_CONSTANT);
	filter2D(srcY, srcYX, -1, b, Point(-1, -1), 0, BORDER_CONSTANT);
	filter2D(src, srcXY, -1, ab, Point(-1, -1), 0, BORDER_CONSTANT);
	imshow("srcY", srcY);
	imshow("srcYX", srcYX);
	imshow("srcXY", srcXY);
	waitKey(0);//此时图片显示时间为一直停留。（x）为x毫秒
	destroyAllWindows();
	return 0;
}

疑问：代码中，a是3*1矩阵，b和ab是1*1的矩阵。理论上不应该是a是3*1的矩阵，b是1*3矩阵，ab是3*3矩阵吗？？麻烦知道的童鞋评论区留言，谢谢！