OpenCV mat类实现水平投影和垂直投影

淡淡的烟草味﹌ 2022-07-13 14:43 96阅读 0赞

图像经过灰度化和otsu阈值分割，分别绘制水平和垂直投影

#include<iostream>  
    #include <cv.h>  
    #include <highgui.h>  
    #include <cvaux.h>  
    #include <tchar.h> 
    
    using namespace std ;
    using namespace cv ;
    
    int main()
    {
    	//原图
    	Mat srcImage=imread("test.png");  
        imshow("原图",srcImage);  
    	//灰度化
        cvtColor(srcImage,srcImage,CV_RGB2GRAY);  
    	//otsu阈值分割
        threshold(srcImage,srcImage,0,255,CV_THRESH_OTSU+CV_THRESH_BINARY);   
        imshow("二值图",srcImage);  
    	//计算垂直投影
        int *colheight =new int[srcImage.cols];    
    	 //数组必须赋初值为零，否则出错。无法遍历数组。 
        memset(colheight,0,srcImage.cols*4);  
        int value;   
        for(int i=0;i<srcImage.rows;i++)   
            for(int j=0;j<srcImage.cols;j++)    
            {    
                value=srcImage.at<uchar>(i,j);  
                if(value==255)    
                {    
                    colheight[j]++; 
                }    
      
            }    
    
            Mat histogramImage(srcImage.rows,srcImage.cols,CV_8UC1);  
            for(int i=0;i<srcImage.rows;i++)    
                for(int j=0;j<srcImage.cols;j++)    
                {    
                    value=0;  //设置为黑色。  
                    histogramImage.at<uchar>(i,j)=value;    
                }    
              
                for(int i=0;i<srcImage.cols;i++)    
                    for(int j=0;j<colheight[i];j++)    
                    {    
                        value=255;  //设置为白色  
                        histogramImage.at<uchar>(j,i)=value;  
                    }    
                    imshow("垂直投影",histogramImage); 
                    //计算水平投影
    				int *colheighttwo =new int[srcImage.rows];    
    				//数组必须赋初值为零，否则出错。无法遍历数组。 
    				memset(colheighttwo,0,srcImage.rows*4);  
    				int valuetwo; 
    				for(int i=0;i<srcImage.rows;i++)    
    					for(int j=0;j<srcImage.cols;j++)    
    					{    
    						valuetwo=srcImage.at<uchar>(i,j);  
    						if(valuetwo==255)    
    						{    
    							colheighttwo[i]++; 
    						}    
    					}    
    
    					Mat plantImage(srcImage.rows,srcImage.cols,CV_8UC1);  //创建一个新的mat型
    					//把这个图全部画成黑色
    					for(int i=0;i<srcImage.rows;i++)    
    						for(int j=0;j<srcImage.cols;j++)    
    						{    
    							valuetwo=0;  //设置为黑色。  
    							plantImage.at<uchar>(i,j)=valuetwo;    
    						}    
    
    						for(int i=0;i<srcImage.rows;i++)    
    							for(int j=0;j<colheighttwo[i];j++)    
    							{    
    								valuetwo=255;  //设置为白色  
    								plantImage.at<uchar>(i,j)=valuetwo;  
    							}    
    							imshow("水平投影",plantImage); 
    			             	waitKey(0);
    			            	return 0;
    }

效果：  
![这里写图片描述][487b8be8d7e49fadc7095984bfcde99a.png]

![这里写图片描述][483d097a3ffde4bd16b9f39b8d3541e2.png]

![这里写图片描述][844b13f97fca59dbb35b391bdd561cbf.png]

![这里写图片描述][ae190872034f1a6027cf3b0b8697283b.png]

[487b8be8d7e49fadc7095984bfcde99a.png]: /images/20220713/27c6d343bc8f4425b9e089f035903a62.png
[483d097a3ffde4bd16b9f39b8d3541e2.png]: /images/20220713/130bc295f2e84d6ea4211b3f370fec47.png
[844b13f97fca59dbb35b391bdd561cbf.png]: /images/20220713/7aff8a8b57ed4f56afe1cc9d50e33d4f.png
[ae190872034f1a6027cf3b0b8697283b.png]: /images/20220713/e5f25912689e4915893447307735d5a4.png