Adaptive thresholding using the integral image

大家好，又见面了，我是你们的朋友全栈君。#include <opencv2/opencv.hpp>  

#include <opencv2/legacy/compat.hpp> 

#include <opencv2/imgproc/types_c.h>

#include <fstream>

using namespace std;  

void vvAdaptiveThreshold( IplImage* inImg, IplImage* outImg)

{


int S = inImg->width >> 5;


int T = 10;

char *input, *bin;


input = inImg->imageData;


bin = outImg->imageData;

int width = inImg->width;


int height = inImg->height;

unsigned long* integralImg = 0;


int i, j;


long sum=0;


int count=0;


int index;


int x1, y1, x2, y2;


int s2 = S/2;

//bin = new unsigned char[width*height];


// create the integral image


integralImg = (unsigned long*)malloc(width*height*sizeof(unsigned long*));


for (i=0; i<width; i++)


{


// reset this column sum


sum = 0;


for (j=0; j<height; j++)


{


index = j*width+i;


sum += input[index];


if (i==0)


integralImg[index] = sum;


else


integralImg[index] = integralImg[index-1] + sum;


}


}


// perform thresholding


for (i=0; i<width; i++)


{


for (j=0; j<height; j++)


{


index = j*width+i;


// set the SxS region


x1=i-s2; x2=i+s2;


y1=j-s2; y2=j+s2;


// check the border


if (x1 < 0) x1 = 0;


if (x2 >= width) x2 = width-1;


if (y1 < 0) y1 = 0;


if (y2 >= height) y2 = height-1;


count = (x2-x1)*(y2-y1);


// I(x,y)=s(x2,y2)-s(x1,y2)-s(x2,y1)+s(x1,x1)


sum = integralImg[y2*width+x2] –


integralImg[y1*width+x2] –


integralImg[y2*width+x1] +


integralImg[y1*width+x1];


if ((long)(input[index]*count) < (long)(sum*(100-T)/100))


bin[index] = 255;


else


bin[index] = 0;


}


}


free (integralImg);

}

int main()  //欢迎大家加入图像识别技术交流群：271891601

{

// 从文件中加载原图  

//IplImage *pSrcImage1 = cvLoadImage(“adaptive.jpg”, CV_LOAD_IMAGE_UNCHANGED);  

IplImage *pSrcImage1 = cvLoadImage(“ding9.jpeg”, CV_LOAD_IMAGE_UNCHANGED); 

IplImage *pSrcImage = cvCreateImage(cvGetSize(pSrcImage1), IPL_DEPTH_8U,1);

cvCvtColor(pSrcImage1,pSrcImage,CV_BGR2GRAY);//cvCvtColor(src,des,CV_BGR2GRAY)  

//创建输出的图像

IplImage *pOutImage = cvCreateImage(cvGetSize(pSrcImage), IPL_DEPTH_8U,1);

vvAdaptiveThreshold(pSrcImage,pOutImage);

const char *pstrWindowsATitle = “Adptive Thresholding using the Integral Image”;  

cvNamedWindow(pstrWindowsATitle, CV_WINDOW_AUTOSIZE); 

cvShowImage(pstrWindowsATitle, pOutImage);   

const char *pstrWindowsATitle1 = “origial Image”;  

cvNamedWindow(pstrWindowsATitle1, CV_WINDOW_AUTOSIZE); 

cvShowImage(pstrWindowsATitle1, pSrcImage1);   

cvWaitKey(0); 

//这里记得释放掉不用的资源哦，由于这是测试程序，就不写了，正式的程序一定要写上哦

return 0;

}