第一种:Zhang并行算法(个人表示来搞笑的)
p1 p2 p3
p8 p1 p4
p7 p6 p5
一次细化算法份两次迭代
第一次删去满足下列三个条件的点
2<= p2+p3+p4+p5+p6+p7+p8+p9<=6(0孤立点1端点78内部点)
p2->p9的排列顺序中,01模式的数量(这里假设二值图非零值为1)为1
P2*p4*p6 = 0且p4*p6*p8 = 0(移去东南角点)
第二次删去满足下列三个条件的点
2<= p2+p3+p4+p5+p6+p7+p8+p9<=6(0孤立点1端点78内部点)
p2->p9的排列顺序中,01模式的数量(这里假设二值图非零值为1)为1
P2*p4*p8 = 0且p2*p6*p8 = 0(移去西北角点)
#include<opencv2/opencv.hpp>
using namespace cv;
void cvThin(cv::Mat& src, cv::Mat& dst, int intera){
if(src.type()!=CV_8UC1){
printf("只能处理二值或灰度图像\n");
return;
}
//非原地操作时候,copy src到dst
if(dst.data!=src.data)src.copyTo(dst);
int i, j, n,width, height;
width = src.cols -1;//之所以减1,是方便处理8邻域,防止越界
height = src.rows -1;
int step = src.step,p2,p3,p4,p5,p6,p7,p8,p9;
uchar* img;
bool ifEnd;
int A1;
cv::Mat tmpimg;//n表示迭代次数
for(n = 0; n<intera; n++){
dst.copyTo(tmpimg);
ifEnd = false;
img = tmpimg.data;
for(i = 1; i < height; i++){
img += step;
for(j =1; j<width; j++){
uchar* p = img + j;
A1 = 0;
if( p[0] > 0){
if(p[-step]==0&&p[-step+1]>0)A1++; //p2,p3 01模式
if(p[-step+1]==0&&p[1]>0)A1++; //p3,p4 01模式
if(p[1]==0&&p[step+1]>0)A1++; //p4,p5 01模式
if(p[step+1]==0&&p[step]>0)A1++; //p5,p6 01模式
if(p[step]==0&&p[step-1]>0)A1++; //p6,p7 01模式
if(p[step-1]==0&&p[-1]>0)A1++; //p7,p8 01模式
if(p[-1]==0&&p[-step-1]>0)A1++; //p8,p9 01模式
if(p[-step-1]==0&&p[-step]>0)A1++; //p9,p2 01模式
p2 = p[-step]>0?1:0;
p3 = p[-step+1]>0?1:0;
p4 = p[1]>0?1:0;
p5 = p[step+1]>0?1:0;
p6 = p[step]>0?1:0;
p7 = p[step-1]>0?1:0;
p8 = p[-1]>0?1:0;
p9 = p[-step-1]>0?1:0;
if((p2+p3+p4+p5+p6+p7+p8+p9)>1 && (p2+p3+p4+p5+p6+p7+p8+p9)<7 && A1==1){
if((p2==0||p4==0||p6==0)&&(p4==0||p6==0||p8==0)){ //p2*p4*p6=0 && p4*p6*p8==0
dst.at<uchar>(i,j) = 0; //满足删除条件,设置当前像素为0
ifEnd = true;
}
}
}
}
}
dst.copyTo(tmpimg);
img = tmpimg.data;
for(i = 1; i < height; i++){
img += step;
for(j =1; j<width; j++){
A1 = 0;
uchar* p = img + j;
if( p[0] > 0){
if(p[-step]==0&&p[-step+1]>0)A1++; //p2,p3 01模式
if(p[-step+1]==0&&p[1]>0)A1++; //p3,p4 01模式
if(p[1]==0&&p[step+1]>0)A1++; //p4,p5 01模式
if(p[step+1]==0&&p[step]>0)A1++; //p5,p6 01模式
if(p[step]==0&&p[step-1]>0)A1++; //p6,p7 01模式
if(p[step-1]==0&&p[-1]>0)A1++; //p7,p8 01模式
if(p[-1]==0&&p[-step-1]>0)A1++; //p8,p9 01模式
if(p[-step-1]==0&&p[-step]>0)A1++; //p9,p2 01模式
p2 = p[-step]>0?1:0;
p3 = p[-step+1]>