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 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116
| #include<opencv2\opencv.hpp> #include<iostream>
using namespace std; using namespace cv; Mat src, dst; int threshold_value = 100; int threshold_max = 255; RNG rng; string input_title = "input img"; string output_title = "output img";
int main(int argc, char * argv) { src = imread("C:\\Users\\Administrator\\Pictures\\lab\\红桃.png"); namedWindow(input_title, CV_WINDOW_AUTOSIZE); namedWindow(output_title, CV_WINDOW_AUTOSIZE); for (int row = 0; row < src.rows; row++) { for (int col = 0; col < src.cols; col++) { if (src.at<Vec3b>(row,col) == Vec3b(255,255,255)) { src.at<Vec3b>(row, col)[0] = 0; src.at<Vec3b>(row, col)[1] = 0; src.at<Vec3b>(row, col)[2] = 0; } } } Mat structureElement = getStructuringElement(MORPH_RECT, Size(5, 5), Point(-1, -1)); dilate(src, src, structureElement, Point(-1, -1), 1);
Mat kernel = (Mat_<float>(3, 3) << 1, 1, 1, 1, -8, 1, 1, 1, 1); Mat imgLaplance; Mat shapenImg = src; filter2D(src, imgLaplance, CV_32F, kernel, Point(-1, -1), 0, BORDER_DEFAULT); src.convertTo(shapenImg, CV_32F); Mat resultImg = shapenImg - imgLaplance; resultImg.convertTo(resultImg, CV_8UC3); imgLaplance.convertTo(imgLaplance, CV_8UC3);
Mat binaryImg; cvtColor(resultImg, resultImg, CV_BGR2GRAY); threshold(resultImg, binaryImg, 40, 255, CV_THRESH_BINARY | THRESH_OTSU);
Mat distImg; distanceTransform(binaryImg, distImg, DIST_L1, 3, 5); normalize(distImg, distImg, 0, 1, NORM_MINMAX);
threshold(distImg, distImg, 0.4, 1, THRESH_BINARY); Mat kernel1 = Mat::zeros(13, 13, CV_8UC1); erode(distImg, distImg, kernel1);
Mat dist_8u; distImg.convertTo(dist_8u, CV_8U); vector<vector<Point>> contours; findContours(dist_8u, contours, RETR_TREE, CHAIN_APPROX_SIMPLE, Point(0, 0));
Mat markers = Mat::zeros(src.size(), CV_32SC1); for (size_t i = 0; i < contours.size(); i++) { drawContours(markers, contours, static_cast<int>(i), Scalar::all (static_cast<int>(i) + 1),-1); } circle(markers, Point(5, 5), 3, Scalar(255, 255, 255), -1);
watershed(src, markers); Mat res = Mat::zeros(markers.size(), CV_8UC1); markers.convertTo(res, CV_8UC1); bitwise_not(res, res, Mat());
vector<Vec3b> colors; for (size_t i = 0; i < contours.size(); i++) { int r = rng.uniform(0, 255); int g = rng.uniform(0, 255); int b = rng.uniform(0, 255);
colors.push_back(Vec3b((uchar)b, (uchar)g, (uchar)r)); }
dst = Mat::zeros(markers.size(), CV_8UC3); for (int row = 0; row < src.rows; row++) { for (int col = 0; col < src.cols; col++) { int index = markers.at<int>(row, col); if (index>0 && index<=static_cast<int>(contours.size())) { dst.at<Vec3b>(row, col) = colors[index - 1]; } else { dst.at<Vec3b>(row, col) = Vec3b(0, 0, 0); } } }
imshow(input_title, src); imshow(output_title, dst);
waitKey(0); return 0; }
|