OpenCV边缘提取(VS2019 & C++)

#include<opencv2/opencv.hpp>
#include<iostream>
#include<opencv2/highgui/highgui_c.h>
#include<vector>

using namespace cv;
using namespace std;


#define max(a,b) a>b?a:b
int main()
{
Mat pupil_image = cv::imread("D:\\opencv\\4.PNG",0);   //imread按单通道的方式读入图像，即灰度图像
imshow("yuantu", pupil_image);
Mat binary_img, spec_mask, kernel;
Mat histogram;
int histSize;
histSize = 256;
float range[] = { 0, 256 };
const float* histRange = { range };
int max_intensity = 40, lowest_spike_index = 255, highest_spike_index = 0;
calcHist(&pupil_image, 1, 0, Mat(), histogram, 1, &histSize, &histRange, true, false);  //计算一系列数组的直方图
for (int i = 0; i < histogram.rows; i++) {
	const float intensity = histogram.at<float>(i, 0);
	if (intensity > 40) {
		max_intensity = max(intensity, max_intensity);
		lowest_spike_index = min(lowest_spike_index, i);
		highest_spike_index = max(highest_spike_index, i);
	}
}
inRange(pupil_image, Scalar(0), Scalar(lowest_spike_index + 20), binary_img);  //指定阈值实现二值化功能

kernel = getStructuringElement(MORPH_ELLIPSE, Size(5, 5));    //构造一个椭圆形结构元素

dilate(binary_img, binary_img, kernel, Point(-1, -1), 2);     //膨胀

inRange(pupil_image, Scalar(0), Scalar(highest_spike_index - 20), spec_mask);

erode(spec_mask, spec_mask, kernel);                           //腐蚀

kernel = getStructuringElement(MORPH_ELLIPSE, Size(9, 9));
//下面开始处理原图
morphologyEx(pupil_image, pupil_image, MORPH_OPEN, kernel);    //形态学开运算
medianBlur(pupil_image, pupil_image, 3);                       //中值滤波来平滑(模糊)处理图像

Mat edges;

Canny(pupil_image, edges, 15, 40, 3);                          //边缘检测
(min)(edges, spec_mask, edges);
(min)(edges, binary_img, edges);
imshow("out", edges);
imwrite("D://out.JPG", edges);
 waitKey(0);
 return 0;
}

检测结果如下：