opencv实践::对象提取与测量

问题描述
照片是来自太空望远镜的星云图像，科学家想知道它的面 积与周长。 

解决思路
方法一： 通过二值分割+图像形态学+轮廓提取

#include <opencv2/opencv.hpp>
#include <iostream>
#include <math.h>

using namespace cv;
using namespace std;

int main(int argc, char** argv) {
    Mat src = imread("D:/case7.png");
    if (src.empty()) {
        printf("could not load image...\n");
        return -;
    }
    namedWindow("input image", CV_WINDOW_AUTOSIZE);
    imshow("input image", src);

    Mat blurImage;
    //高斯模糊
    GaussianBlur(src, blurImage, Size(, ), , );
    //imshow("blur", blurImage);

    Mat gray_src, binary;
    cvtColor(blurImage, gray_src, COLOR_BGR2GRAY);
    threshold(gray_src, binary, , , THRESH_BINARY | THRESH_TRIANGLE);
    //imshow("binary", binary);

    // 形态学操作
    Mat morphImage;
    Mat kernel = getStructuringElement(MORPH_RECT, Size(, ), Point(-, -));
    morphologyEx(binary, morphImage, MORPH_CLOSE, kernel, Point(-, -), );
    //imshow("morphology", morphImage);

    // 获取最大轮廓
    vector<vector<Point>> contours;
    vector<Vec4i> hireachy;
    findContours(morphImage, contours, hireachy, CV_RETR_EXTERNAL, CHAIN_APPROX_SIMPLE, Point());
    Mat connImage = Mat::zeros(src.size(), CV_8UC3);
    for (size_t t = ; t < contours.size(); t++) {
        Rect rect = boundingRect(contours[t]);
        if (rect.width < src.cols / ) continue;
        if (rect.width > (src.cols - )) continue;
        //计算面积
        double area = contourArea(contours[t]);
        //计算长度
        double len = arcLength(contours[t], true);
        drawContours(connImage, contours, static_cast<int>(t), Scalar(, , ), , , hireachy);
        printf("area  of star could : %f\n", area);
        printf("length  of star could : %f\n", len);
    }
    imshow("result", connImage);

    waitKey();
    return ;
}