Opencv系列9_opencv区域增长算法分割JPG图像并显示保存

﹏ヽ暗。殇╰゛Y 2022-09-03 03:29 200阅读 0赞

#  **实例9：opencv区域增长算法分割JPG图像并显示保存** #

#include "opencv2/highgui/highgui.hpp"
    #include "opencv2/imgproc.hpp"
    #include "opencv2/imgcodecs.hpp"
    #include <iostream>
    #include "math.h"
    
    using namespace cv;
    using namespace std;
    
    Point recent_Point, recent_Point1;
    
    Mat RegionGrow(Mat src, Point2i pt, int th)
    {
        Point2i ptGrowing;                        //待生长点位置
        int nGrowLable = 0;                                //标记是否生长过
        int nSrcValue = 0;                                //生长起点灰度值
        int nCurValue = 0;                                //当前生长点灰度值
        Mat matDst = Mat::zeros(src.size(), CV_8UC1);    //创建一个空白区域，填充为黑色
        //生长方向顺序数据
        int DIR[8][2] = { { -1, -1 }, { 0, -1 }, { 1, -1 }, { 1, 0 }, { 1, 1 }, { 0, 1 }, { -1, 1 }, { -1, 0 } };
        vector<Point2i> vcGrowPt;                        //生长点栈
        vcGrowPt.push_back(pt);                            //将生长点压入栈中
        matDst.at<uchar>(pt.y, pt.x) = 255;                //标记生长点
        nSrcValue = src.at<uchar>(pt.y, pt.x);            //记录生长点的灰度值
    
        while (!vcGrowPt.empty())                        //生长栈不为空则生长
        {
            pt = vcGrowPt.back();                        //取出一个生长点
            vcGrowPt.pop_back();
    
            //分别对八个方向上的点进行生长
            for (int i = 0; i < 9; ++i)
            {
                ptGrowing.x = pt.x + DIR[i][0];
                ptGrowing.y = pt.y + DIR[i][1];
                //检查是否是边缘点
                if (ptGrowing.x < 0 || ptGrowing.y < 0 || ptGrowing.x >(src.cols - 1) || (ptGrowing.y > src.rows - 1))
                    continue;
    
                nGrowLable = matDst.at<uchar>(ptGrowing.y, ptGrowing.x);        //当前待生长点的灰度值
    
                if (nGrowLable == 0)                    //如果标记点还没有被生长
                {
                    nCurValue = src.at<uchar>(ptGrowing.y, ptGrowing.x);
                    if (abs(nSrcValue - nCurValue) < th)                    //在阈值范围内则生长
                    {
                        matDst.at<uchar>(ptGrowing.y, ptGrowing.x) = 255;        //标记为白色
                        vcGrowPt.push_back(ptGrowing);                    //将下一个生长点压入栈中
                    }
                }
            }
        }
        return matDst.clone();
    }
    
    void On_mouse(int event, int x, int y, int flags, void*)//每次点击左键，在相应位置画红点
    {
        if (event == EVENT_LBUTTONDOWN)
        {
            recent_Point = Point(x, y);
            cout << "img_x" << " " << recent_Point.x << " " << "img_y" << " " << recent_Point.y << endl;
            //circle(srcimg, recent_Point, 2, Scalar(0, 0, 255), -1);
            //imshow("srcImg", srcimg);
        }
    }
    
    int main() //区域生长
    {
        Mat binaryimg, greyimg;
        Mat regiongrow, regiongrow1, regiongrow2;
        Mat dst;
        int th = 10;
        Mat src = imread("1.jpg");
        cvtColor(src, greyimg, COLOR_BGR2GRAY);   //转化为灰度图
        Mat temp_regiongrow = Mat::zeros(src.size(), CV_8UC1);    //创建一个空白区域，填充为黑色
        //转化为二值图
        threshold(greyimg, binaryimg, 200, 255, THRESH_BINARY);
        namedWindow("srcImg", 0);
        imshow("srcImg", src);
    
        namedWindow("binaryImg", 0);
        imshow("binaryImg", binaryimg);
        cout << "select one point in binaryImg" << endl;
        setMouseCallback("binaryImg", On_mouse);
    
        for (int i = 0;i < 1;i++) {
            char c = (char)waitKey(0);
            cout << "select one point in binaryImg" << endl;
            setMouseCallback("binaryImg", On_mouse);
            if (c == 'b') {
                regiongrow1 = RegionGrow(binaryimg, recent_Point, th);
                bitwise_or(regiongrow1, temp_regiongrow, regiongrow);    //和前一个分割的图做或运算
                temp_regiongrow = regiongrow1.clone();    //保存前一个分割图
            }
            bitwise_and(greyimg, regiongrow, dst);   //与原图做与运算
            namedWindow("dstimg", 0);
            imshow("dstimg", dst);
            imwrite("region_growing.jpg", dst);
        }
        waitKey(0);
        return 0;
    }

在binaryImag图像中选择区域增长的初始种子点（鼠标左键一下红色区域选取）：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1][]

出现位置坐标后，键盘输入字母‘b'进行区域增长算法运行，dstimg为区域增长分割后图像：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3][]

**[本例程配套素材见源码整理文章下载（点击进入）][Link 1]**

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20210421141043265.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1]: /images/20220829/d6e3086b06724d83b19b3d34398a7068.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2]: https://img-blog.csdnimg.cn/20210421142637145.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3]: /images/20220829/1cc6e624167141da95dea912fe0fad02.png
[Link 1]: https://blog.csdn.net/qq_32809093/article/details/115689942