使用 OpenCV 在斑点内创建矩形

输入图像：

输出图像：

我的图像中有几个彩色斑点，我试图在每种颜色的最大斑点内创建矩形（或正方形——这似乎更容易）。我发现了如何创建一个包围单个最大斑点的矩形​​的答案，但我不确定如何找到一个适合在斑点内的正方形。它不一定是最大的，它只需要大于某个区域，否则我不会包括它。我还看到了一些关于多边形的工作，但没有看到关于无定形形状的工作。

对于单个 blob，问题可以表述为：找到矩阵中仅包含零的最大矩形.

要找到斑点内最大的轴向矩形，您可以参考该函数findMinRect in 我的另一个答案。该代码是从Python 中的原始代码移植到C++ 中的here.

那么第二个问题是找到所有具有相同颜色的斑点。这有点棘手，因为您的图像是 jpeg，并且压缩会在边界附近产生很多人造颜色。所以我创建了一个 png 图像（如下所示），只是为了表明该算法是有效的。您可以提供没有压缩伪影的图像。

然后，您只需为每种颜色创建一个遮罩，找到该遮罩中每个斑点的连通分量，并计算每个斑点的最小矩形。

初始图像：

在这里，我显示了为每个斑点找到的矩形，按颜色划分。然后，您可以仅采用所需的矩形，无论是每种颜色的最大矩形，还是每种颜色的最大斑点的矩形​​。

Result:

这里是代码：

#include <opencv2/opencv.hpp>
#include <algorithm>
#include <set>
using namespace std;
using namespace cv;

// https://stackoverflow.com/a/30418912/5008845
Rect findMinRect(const Mat1b& src)
{
    Mat1f W(src.rows, src.cols, float(0));
    Mat1f H(src.rows, src.cols, float(0));

    Rect maxRect(0, 0, 0, 0);
    float maxArea = 0.f;

    for (int r = 0; r < src.rows; ++r)
    {
        for (int c = 0; c < src.cols; ++c)
        {
            if (src(r, c) == 0)
            {
                H(r, c) = 1.f + ((r>0) ? H(r - 1, c) : 0);
                W(r, c) = 1.f + ((c>0) ? W(r, c - 1) : 0);
            }

            float minw = W(r, c);
            for (int h = 0; h < H(r, c); ++h)
            {
                minw = min(minw, W(r - h, c));
                float area = (h + 1) * minw;
                if (area > maxArea)
                {
                    maxArea = area;
                    maxRect = Rect(Point(c - minw + 1, r - h), Point(c + 1, r + 1));
                }
            }
        }
    }

    return maxRect;
}


struct lessVec3b
{
    bool operator()(const Vec3b& lhs, const Vec3b& rhs) {
        return (lhs[0] != rhs[0]) ? (lhs[0] < rhs[0]) : ((lhs[1] != rhs[1]) ? (lhs[1] < rhs[1]) : (lhs[2] < rhs[2]));
    }
};

int main()
{
    // Load image
    Mat3b img = imread("path_to_image");

    // Find unique colors
    set<Vec3b, lessVec3b> s(img.begin(), img.end());

    // Divide planes of original image
    vector<Mat1b> planes;
    split(img, planes);
    for (auto color : s)
    {
        // Create a mask with only pixels of the given color
        Mat1b mask(img.rows, img.cols, uchar(255));
        for (int i = 0; i < 3; ++i)
        {
            mask &= (planes[i] == color[i]);
        }

        // Find blobs
        vector<vector<Point>> contours;
        findContours(mask, contours, RETR_EXTERNAL, CHAIN_APPROX_SIMPLE);

        for (int i = 0; i < contours.size(); ++i)
        {
            // Create a mask for each single blob
            Mat1b maskSingleContour(img.rows, img.cols, uchar(0));
            drawContours(maskSingleContour, contours, i, Scalar(255), CV_FILLED);

            // Find minimum rect for each blob
            Rect box = findMinRect(~maskSingleContour);

            // Draw rect
            Scalar rectColor(color[1], color[2], color[0]);
            rectangle(img, box, rectColor, 2);
        }
    }

    imshow("Result", img);
    waitKey();

    return 0;
}