这是一个查找轮廓算法,他的原理是从图像的0行,0列开始扫描,找到当前像素,与之前一个像素不同的像素并改变他的颜色,以作为标记。其中要是用旋转算法,对图像进行顺时针和逆时针旋转,在进行下一轮查找,已完成所有轮廓的区分。
#include <opencv2/opencv.hpp>
using namespace std;
using namespace cv;
void find(Mat m) {
int h = m.rows;
int w = m.cols;
for (int i = 0; i < h; i++)
{
for (int j = 1; j < w; j++)
{
if (m.at<Vec3b>(i, j-1)[0] != m.at<Vec3b>(i, j)[0]
&& m.at<Vec3b>(i, j-1)[1] != m.at<Vec3b>(i, j)[1]
&& m.at<Vec3b>(i, j-1)[2] != m.at<Vec3b>(i, j)[2]) {
m.at<Vec3b>(i, j)[0] = 255;
m.at<Vec3b>(i, j)[1] = 0;
m.at<Vec3b>(i, j)[2] = 255;
}
}
}
}
Mat imgRotate(Mat matSrc, float angle, bool direction)
{
float theta = angle * CV_PI / 180.0;
int nRowsSrc = matSrc.rows;
int nColsSrc = matSrc.cols;
// 如果是顺时针旋转
if (!direction)
theta = 2 * CV_PI - theta;
// 全部以逆时针旋转来计算
// 逆时针旋转矩阵
float matRotate[3][3]{
{std::cos(theta), -std::sin(theta), 0},
{std::sin(theta), std::cos(theta), 0 },
{0, 0, 1}
};
float pt[3][2]{
{ 0, nRowsSrc },
{nColsSrc, nRowsSrc},
{nColsSrc, 0}
};
for (int i = 0; i < 3; i++)
{
float x = pt[i][0] * matRotate[0][0] + pt[i][1] * matRotate[1][0];
float y = pt[i][0] * matRotate[0][1] + pt[i][1] * matRotate[1][1];
pt[i][0] = x;
pt[i][1] = y;
}
// 计算出旋转后图像的极值点和尺寸
float fMin_x = min(min(min(pt[0][0], pt[1][0]), pt[2][0]), (float)0.0);
float fMin_y = min(min(min(pt[0][1], pt[1][1]), pt[2][1]), (float)0.0);
float fMax_x = max(max(max(pt[0][0], pt[1][0]), pt[2][0]), (float)0.0);
float fMax_y = max(max(max(pt[0][1], pt[1][1]), pt[2][1]), (float)0.0);
int nRows = cvRound(fMax_y - fMin_y + 0.5) + 1;
int nCols = cvRound(fMax_x - fMin_x + 0.5) + 1;
int nMin_x = cvRound(fMin_x + 0.5);
int nMin_y = cvRound(fMin_y + 0.5);
// 拷贝输出图像
Mat matRet(nRows, nCols, matSrc.type(), Scalar(0));
for (int j = 0; j < nRows; j++)
{
for (int i = 0; i < nCols; i++)
{
// 计算出输出图像在原图像中的对应点的坐标,然后复制该坐标的灰度值
// 因为是逆时针转换,所以这里映射到原图像的时候可以看成是,输出图像
// 到顺时针旋转到原图像的,而顺时针旋转矩阵刚好是逆时针旋转矩阵的转置
// 同时还要考虑到要把旋转后的图像的左上角移动到坐标原点。
int x = (i + nMin_x) * matRotate[0][0] + (j + nMin_y) * matRotate[0][1];
int y = (i + nMin_x) * matRotate[1][0] + (j + nMin_y) * matRotate[1][1];
if (x >= 0 && x < nColsSrc && y >= 0 && y < nRowsSrc)
{
matRet.at<Vec3b>(j, i) = matSrc.at<Vec3b>(y, x);
}
}
}
return matRet;
}
int main() {
Mat m = imread("D:/jj.png");//载入图像
float angle = 90.0;//旋转图像的度数
find(m);//查找轮廓
Mat matret = imgRotate(m, angle, true);//逆时针旋转图像
find(matret);//查找轮廓
matret = imgRotate(matret, angle, false);//顺时针旋转图像
imshow("jj.png", m);//显示图像
imshow("rotate", matret);//显示图像
waitKey(0);
return 0;
}
查找之前的图像
查找之后的图像,可以看到黑色的物体已经被圈起来了。
但是有个小问题,这个旋转图像的算法在进行旋转的时候会出现黑边,不知道为什么,希望有人可以指出来。