【OpenCV】计算图像的水平和垂直积分投影

#include <cv.h>  
#include <highgui.h>  
#pragma comment( lib, "cv.lib" )  
#pragma comment( lib, "cxcore.lib" )  
#pragma comment( lib, "highgui.lib" )   
int main()  
{     
    IplImage * src=cvLoadImage("lena.jpg",0);  
    //cvSmooth(src,src,CV_BLUR,3,3,0,0);   
    cvThreshold(src,src,50,255,CV_THRESH_BINARY_INV);  
    IplImage* paintx=cvCreateImage( cvGetSize(src),IPL_DEPTH_8U, 1 );  
    IplImage* painty=cvCreateImage( cvGetSize(src),IPL_DEPTH_8U, 1 );  
    cvZero(paintx);  
    cvZero(painty);  
    int* v=new int[src->width];  
    int* h=new int[src->height];  
    memset(v,0,src->width*4);  
    memset(h,0,src->height*4);  
      
    int x,y;  
    CvScalar s,t;  
    for(x=0;x<src->width;x++)  
    {  
        for(y=0;y<src->height;y++)  
        {  
            s=cvGet2D(src,y,x);           
            if(s.val[0]==0)  
                v[x]++;                   
        }         
    }  
      
    for(x=0;x<src->width;x++)  
    {  
        for(y=0;y<v[x];y++)  
        {         
            t.val[0]=255;  
            cvSet2D(paintx,y,x,t);        
        }         
    }  
      
    for(y=0;y<src->height;y++)  
    {  
        for(x=0;x<src->width;x++)  
        {  
            s=cvGet2D(src,y,x);           
            if(s.val[0]==0)  
                h[y]++;       
        }     
    }  
    for(y=0;y<src->height;y++)  
    {  
        for(x=0;x<h[y];x++)  
        {             
            t.val[0]=255;  
            cvSet2D(painty,y,x,t);            
        }         
    }  
    cvNamedWindow("二值图像",1);  
    cvNamedWindow("垂直积分投影",1);  
    cvNamedWindow("水平积分投影",1);  
    cvShowImage("二值图像",src);  
    cvShowImage("垂直积分投影",paintx);  
    cvShowImage("水平积分投影",painty);  
    cvWaitKey(0);  
    cvDestroyAllWindows();  
    cvReleaseImage(&src);  
    cvReleaseImage(&paintx);  
    cvReleaseImage(&painty);  
    return 0;  
}