使用pointPolygonTest函数获取一个点在轮廓外或轮廓上或轮廓内。
当使用距离测量时,点在轮廓外返回一个负值,越小代表与轮廓距离越远;点在轮廓上时返回0;点在轮廓内将返回一个正值,距离越远值越大。
当不使用距离测量时,负值或正值都是一个恒定值-1或+1
#include "cv.h"
#include "highgui.h"
using namespace cv;
using namespace std;
int main(int argc,char *argv[])
{
Mat src;
src=Mat::zeros(300,300,CV_8UC1);
vector<Point2f> vert(6);
vert[0]=Point(100,75);
vert[1]=Point(200,75);
vert[2]=Point(250,150);
vert[3]=Point(200,225);
vert[4]=Point(100,225);
vert[5]=Point(50,150);
for (int i=0;i<6;i++)
line(src,vert[i],vert[(i+1)%6],Scalar::all(255),1);
imshow("src",src);
vector<vector<Point>> contours;
vector<Vec4i> hierarchy;
findContours(src,contours,hierarchy,RETR_TREE,CHAIN_APPROX_SIMPLE);
Mat raw(src.size(),CV_32FC1);
for (int j=0;j<src.rows;j++)
{
for (int i=0;i<src.cols;i++)
{
raw.at<float>(j,i)=pointPolygonTest(contours[0],Point2f(i,j),true);//true表示支持距离测量
}
}
double minVal;double maxVal;
minMaxLoc(raw,&minVal,&maxVal);
Mat drawing=Mat::zeros(src.size(),CV_8UC3);
for (int j=0;j<src.rows;j++)
{
for (int i=0;i<src.cols;i++)
{
if (raw.at<float>(j,i)<0)
{
drawing.at<Vec3b>(j,i)[0]=255-int(raw.at<float>(j,i)/minVal*255);
}
else if (raw.at<float>(j,i)>0)
{
drawing.at<Vec3b>(j,i)[2]=255-int(raw.at<float>(j,i)/maxVal*255);
}
else
{
drawing.at<Vec3b>(j,i)[1]=255;
}
}
}
imshow("drawing",drawing);
waitKey(0);
destroyAllWindows();
return 0;
}
效果:
不使用距离测量时,结果: