opencv3编程入门（第八章8.3-8.6）

最新推荐文章于 2024-03-02 11:58:36 发布

flyer飞亚

最新推荐文章于 2024-03-02 11:58:36 发布

阅读量105

点赞数

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

本文链接：https://blog.csdn.net/weixin_43544406/article/details/107563124

版权

图像处理专栏收录该内容

15 篇文章 1 订阅

订阅专栏

8.4图像的矩：

#include<opencv2/highgui/highgui.hpp>
#include<opencv2/imgproc/imgproc.hpp>
#include<opencv2/opencv.hpp>
#include<iostream>
using namespace cv;
using namespace std;
#define WINDOW_NAME1 "原始图窗口"
#define WINDOW_NAME2 "效果图窗口"
Mat g_srcImage, g_grayImage;
int g_nThresh = 100, g_MaxThresh = 255;
RNG g_rng(12345);
Mat g_cannyMat_output;
vector<vector<Point>>g_vContours;
vector<Vec4i>g_vHierarchy;
static void ShowHelpText();
void on_ThreshChange(int, void*);
int main(int argc,char** aegv)
{
	system("color 1E");
	g_srcImage = imread("hehua.jpg", 1);
	if (!g_srcImage.data) { printf("读取图片错误\n");return false; }
	cvtColor(g_srcImage, g_grayImage, COLOR_BGR2GRAY);
	blur(g_grayImage, g_grayImage, Size(3, 3));
	namedWindow(WINDOW_NAME1, WINDOW_AUTOSIZE);
	imshow(WINDOW_NAME1, g_srcImage);
	//创建滚动条
	createTrackbar("阈值：", WINDOW_NAME1, &g_nThresh, g_MaxThresh, on_ThreshChange);
	on_ThreshChange(0, 0);//调用一次回调函数进行初始化
	waitKey(0);
	return 0;
}
void on_ThreshChange(int, void*)
{
	Canny(g_grayImage, g_cannyMat_output, g_nThresh, g_nThresh * 2, 3);
	findContours(g_cannyMat_output, g_vContours, g_vHierarchy, RETR_TREE,
		CHAIN_APPROX_SIMPLE, Point(0, 0));
	//计算矩
	vector<Moments>mu(g_vContours.size());
	for(unsigned int i=0;i<g_vContours.size();i++)
	{
		mu[i] = moments(g_vContours[i], false);
	}
	//计算中心距
	vector<Point>mc(g_vContours.size());
	for (unsigned int i = 0;i < g_vContours.size();i++)
	{
		mc[i] = Point2f(static_cast<float>(mu[i].m10 / mu[i].m00),
			static_cast<float>(mu[i].m01 / mu[i].m00));
	}
	//绘制轮廓
	Mat drawing = Mat::zeros(g_cannyMat_output.size(), CV_8UC3);
	for (unsigned int i = 0;i < g_vContours.size();i++)
	{
		Scalar color = Scalar(g_rng.uniform(0, 255), g_rng.uniform(0, 255), g_rng.uniform(0, 255));
		drawContours(drawing, g_vContours, i, color, 2, 8, g_vHierarchy, 0, Point());//绘制外层和内层轮廓
		circle(drawing, mc[i], 4, color, -1, 8, 0);//绘制圆
	}
	namedWindow(WINDOW_NAME2, WINDOW_AUTOSIZE);
	imshow(WINDOW_NAME2, drawing);
	//通过m00计算轮廓面积并且和Opencv函数比较
	printf("\t输出内容：面积和轮廓长度\n");
	for (unsigned int i = 0;i < g_vContours.size();i++)
	{
		printf(">通过m00计算出轮廓[%d]的面积：(M_00)=%.2f\nOpencv函数算出的面积 = %.2f,长度：%.2f\n\n",
			i,mu[i].m00,
			contourArea(g_vContours[i]), arcLength(g_vContours[i], true));
		Scalar color = Scalar(g_rng.uniform(0, 255), g_rng.uniform(0, 255), g_rng.uniform(0, 255));
		drawContours(drawing, g_vContours, i, color, 2, 8, g_vHierarchy, 0, Point());
		circle(drawing, mc[i], 4, color, -1, 8, 0);
	}
}

程序执行结果如图：
在这里插入图片描述
8.5分水岭算法：

#include<opencv2/highgui/highgui.hpp>
#include<opencv2/imgproc/imgproc.hpp>
#include<opencv2/opencv.hpp>
#include<iostream>
using namespace cv;
using namespace std;
#define WINDOW_NAME "窗口1"
Mat g_srcImage, g_maskImage;
Point prevPt(-1, -1);
static void ShowHelpText();
static void on_Mouse(int event, int x, int y, int flags, void*);
int main(int argc, char** aegv)
{
	g_srcImage = imread("huashan.jpg", 1);
	if (!g_srcImage.data) { printf("读取图片错误\n");return false; }
	namedWindow(WINDOW_NAME, WINDOW_AUTOSIZE);
	imshow(WINDOW_NAME, g_srcImage);
	Mat srcImage, grayImage;
	g_srcImage.copyTo(srcImage);
	cvtColor(g_srcImage, g_maskImage, COLOR_BGR2GRAY);
	cvtColor(g_maskImage, grayImage, COLOR_GRAY2BGR);
	g_maskImage = Scalar::all(0);
	setMouseCallback(WINDOW_NAME, on_Mouse, 0);
	while (1)//轮询按键进行处理
	{
		int c = waitKey(0);
		if ((char)c == 27)
			break;
		if ((char)c == '2')
		{
			g_maskImage = Scalar::all(0);
			srcImage.copyTo(g_srcImage);
			imshow("image", g_srcImage);
		}
		if ((char)c == '1' || (char)c == ' ')
		{
			int i, j, compCount = 0;
			vector<vector<Point>>contours;
			vector<Vec4i>hierarchy;
			findContours(g_maskImage, contours, hierarchy,CHAIN_APPROX_SIMPLE, RETR_TREE, Point(0, 0));
			if (contours.empty())
				continue;
			Mat maskImage(g_maskImage.size(), CV_32S);//复制掩膜
			maskImage = Scalar::all(0);
			for (int index = 0;index >= 0;index = hierarchy[index][0], compCount++)
				drawContours(maskImage, contours, index, Scalar::all(compCount + 1),
					-1, 8, hierarchy, INT_MAX);
			if (compCount == 0)
				continue;
			vector<Vec3b>colorTab;
			for (i = 0;i < compCount;i++)
			{
				int b = theRNG().uniform(0, 255);
				int g = theRNG().uniform(0, 255);
				int r = theRNG().uniform(0, 255);
				colorTab.push_back(Vec3b((uchar)b, (uchar)g, (uchar)r));
			}
			double dTime = (double)getTickCount();
			watershed(srcImage, maskImage);
			dTime = (double)getTickCount() - dTime;
			printf("\t处理时间=%gms\n", dTime*1000. / getTickFrequency());
			//双层循环，将分水岭图像遍历存入watershedImage中
			Mat watershedImage(maskImage.size(), CV_8UC3);
			for (i = 0;i < maskImage.rows;i++)
				for (j = 0;j < maskImage.cols;j++)
				{
					int index = maskImage.at<int>(i, j);
					if (index == -1)
						watershedImage.at<Vec3b>(i, j) = Vec3b(255, 255, 255);
					else if (index <= 0 || index > compCount)
						watershedImage.at<Vec3b>(i, j) = Vec3b(0, 0, 0);
					else
						watershedImage.at<Vec3b>(i, j) = colorTab[index - 1];
				}
			//混合灰度图和分水岭效果图并显示最终的窗口
			watershedImage = watershedImage*0.5 + grayImage*0.5;
			imshow("watershed transform", watershedImage);
		}
	}
	return 0;
}
static void on_Mouse(int event, int x, int y, int flags, void*)
{
	if (x < 0 || x >= g_srcImage.cols || y < 0 || y >= g_srcImage.rows)
		return;
	if (event == EVENT_LBUTTONUP || !(flags&EVENT_FLAG_LBUTTON))
		prevPt = Point(-1, -1);
	else if(event==EVENT_LBUTTONDOWN)
		prevPt = Point(x, y);
	else if (event == EVENT_MOUSEMOVE && (flags&EVENT_FLAG_LBUTTON))
	{
		Point pt(x, y);
		if (prevPt.x < 0)
			prevPt = pt;
		line(g_maskImage, prevPt, pt, Scalar::all(255), 5, 8, 0);
		line(g_srcImage, prevPt, pt, Scalar::all(255), 5, 8, 0);
		prevPt = pt;
		imshow(WINDOW_NAME, g_srcImage);

	}
}

程序执行结果如图：
在这里插入图片描述
8.6图像修补：

#include<opencv2/highgui/highgui.hpp>
#include<opencv2/imgproc/imgproc.hpp>
#include<opencv2/opencv.hpp>
#include<iostream>
using namespace cv;
using namespace std;
#define WINDOW_NAME1 "窗口1"
#define WINDOW_NAME2 "修补后的效果图"
Mat srcImage1, inpaintMask;
Point previousPoint(-1, -1);
static void on_Mouse(int event, int x, int y, int flags, void*)
{
	if (event== EVENT_LBUTTONUP || !(flags&EVENT_FLAG_LBUTTON))
		previousPoint=Point(-1,-1);
	else if (event == EVENT_LBUTTONDOWN)
		previousPoint = Point(x, y);
	else if (event == EVENT_MOUSEMOVE && (flags&EVENT_FLAG_LBUTTON))
	{
		Point pt(x, y);
		if (previousPoint.x < 0)
			previousPoint = pt;
		line(inpaintMask, previousPoint, pt, Scalar::all(255), 5, 8, 0);
		line(srcImage1, previousPoint, pt, Scalar::all(255), 5, 8, 0);
		previousPoint = pt;
		imshow(WINDOW_NAME1, srcImage1);

	}
}
int main(int argc, char** aegv)
{
	Mat srcImage = imread("huashan.jpg", -1);
	if (!srcImage.data) { printf("读取图片错误\n");return false; }
	srcImage1 = srcImage.clone();
	inpaintMask = Mat::zeros(srcImage1.size(), CV_8U);
	imshow(WINDOW_NAME1, srcImage1);
	setMouseCallback(WINDOW_NAME1, on_Mouse, 0);
	while (1)//轮询按键进行处理
	{
		int c =(char)waitKey();
		if ((char)c == 27)
			break;
		if ((char)c == '2')
		{
			inpaintMask= Scalar::all(0);
			srcImage.copyTo(srcImage1);
			/*copyTo有两种用法1.image.copyTo(imageROI)，作用是把image的内容粘贴到imageROI
			2.image.copyTo(imageROI，mask),作用是把mask和image重叠后，把mask中像素值为0的
			点变为image的对应点的像素，而保留其他点。*/
			imshow("WINDOW_NAME1", srcImage1);
		}
		if ((char)c == '1' || (char)c == ' ')
		{
			Mat inpaintedImage;
			inpaint(srcImage1, inpaintMask, inpaintedImage, 3, INPAINT_TELEA);
			//inpaintMask为修复掩膜，inpaintedImage为函数调用后的运算结果
			imshow(WINDOW_NAME2, inpaintedImage);
		}
	}
	return 0;
}

程序执行后的结果如图：
在这里插入图片描述

flyer飞亚

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
opencv3编程入门（第八章8.3-8.6）

8.4图像的矩：#include<opencv2/highgui/highgui.hpp>#include<opencv2/imgproc/imgproc.hpp>#include<opencv2/opencv.hpp>#include<iostream>using namespace cv;using namespace std;#define WINDOW_NAME1 "原始图窗口"#define WINDOW_NAME2 "效果图窗口"
复制链接

扫一扫