|
| 1 | +#include <opencv/cv.h> |
| 2 | +#include <opencv/highgui.h> |
| 3 | +#include <stdio.h> |
| 4 | +#include<iostream> |
| 5 | +using namespace std; |
| 6 | + |
| 7 | +void create_histogram_image(IplImage* bin_img, IplImage* hist_img,int bins) { |
| 8 | + CvHistogram *hist; |
| 9 | + |
| 10 | + int fc = 256/bins; |
| 11 | + int hist_size = 256/fc; |
| 12 | + float range[]={0,256/fc}; |
| 13 | + float* ranges[] = { range }; |
| 14 | + float max_value = 0.0; |
| 15 | + float w_scale = 0.0; |
| 16 | + |
| 17 | + // create array to hold histogram values |
| 18 | + hist = cvCreateHist(1, &hist_size, CV_HIST_ARRAY, ranges, 1); |
| 19 | + |
| 20 | + // calculate histogram values |
| 21 | + cvCalcHist( &bin_img, hist, 0, NULL ); |
| 22 | + |
| 23 | + // Get the minimum and maximum values of the histogram |
| 24 | + cvGetMinMaxHistValue( hist, 0, &max_value, 0, 0 ); |
| 25 | + |
| 26 | + // set height by using maximim value |
| 27 | + cvScale( hist->bins, hist->bins, ((float)hist_img->height)/max_value, 0 ); |
| 28 | + |
| 29 | + // calculate width |
| 30 | + w_scale = ((float)hist_img->width)/hist_size; |
| 31 | + |
| 32 | + // plot the histogram |
| 33 | + for( int i = 0; i < hist_size; i++ ) { |
| 34 | + |
| 35 | + cvRectangle( hist_img, cvPoint((int)i*w_scale , hist_img->height), |
| 36 | + cvPoint((int)(i+1)*w_scale, hist_img->height - cvRound(cvGetReal1D(hist->bins,i))), |
| 37 | + cvScalar(0), -1, 8, 0 ); |
| 38 | + |
| 39 | + } |
| 40 | +} |
| 41 | +IplImage* Pre_Process(IplImage* src,int bins) |
| 42 | +{ |
| 43 | + int c = cvGetSize(src).width; |
| 44 | + int r = cvGetSize(src).height; |
| 45 | + |
| 46 | + IplImage* final_img=cvCreateImage(cvSize(src->width,src->height),src->depth,src->nChannels); |
| 47 | + cvZero(final_img); |
| 48 | + for(int i=0;i<r;i++) |
| 49 | + { |
| 50 | + for(int j=0;j<c;j++) |
| 51 | + { |
| 52 | + CvScalar s1,s; |
| 53 | + s1 = cvGet2D(src,i,j); |
| 54 | + int grval = s1.val[0]; |
| 55 | + if(bins==32) |
| 56 | + grval = grval/8; |
| 57 | + if(bins==64) |
| 58 | + grval = grval/4; |
| 59 | + if(bins==128) |
| 60 | + grval = grval/2; |
| 61 | +// cout<<grval<<" "; |
| 62 | + s.val[0]=grval; |
| 63 | + cvSet2D(final_img,i,j,s); |
| 64 | + } |
| 65 | + } |
| 66 | + return final_img; |
| 67 | +} |
| 68 | + |
| 69 | +int main( int argc, char** argv ) |
| 70 | + |
| 71 | +{ |
| 72 | + IplImage *source = cvLoadImage(argv[1],0); |
| 73 | + char *img_nm,*img_typ; |
| 74 | + img_nm = strtok (argv[2],"."); |
| 75 | + img_typ = strtok (NULL, "."); |
| 76 | + int nbins = atoi(argv[3]); |
| 77 | + IplImage *bin_img = cvCreateImage(cvSize( source->width, source->height ), IPL_DEPTH_8U, 1 ); |
| 78 | + if(nbins!=256) |
| 79 | + bin_img = Pre_Process(source,nbins); |
| 80 | + else |
| 81 | + bin_img = cvCloneImage(source); |
| 82 | + |
| 83 | + int c = cvGetSize(bin_img).width; |
| 84 | + int r = cvGetSize(bin_img).height; |
| 85 | + |
| 86 | + IplImage *bin_img1 = cvCreateImage(cvSize( source->width/2, source->height/2 ), IPL_DEPTH_8U, 1 ); |
| 87 | + int mn4=300,mx4=-1,avg4=0; |
| 88 | + for(int i=0;i<r/2;i++) |
| 89 | + { |
| 90 | + for(int j=0;j<c/2;j++) |
| 91 | + { |
| 92 | + CvScalar s; |
| 93 | + s = cvGet2D(bin_img,i,j); |
| 94 | + if(s.val[0] < mn4) |
| 95 | + mn4 = s.val[0]; |
| 96 | + if(s.val[0] > mx4) |
| 97 | + mx4 = s.val[0]; |
| 98 | + avg4+=s.val[0]; |
| 99 | + cvSet2D(bin_img1,i,j,s); |
| 100 | + } |
| 101 | + } |
| 102 | + avg4 = 4*avg4/(r*c); |
| 103 | + |
| 104 | + IplImage *hist_img = cvCreateImage(cvSize(300,300), 8, 1); |
| 105 | + cvSet( hist_img, cvScalarAll(255), 0 ); |
| 106 | + create_histogram_image(bin_img1, hist_img,nbins); |
| 107 | + char nam[30]; |
| 108 | + int n = sprintf(nam,"%s_1.%s",img_nm,img_typ); |
| 109 | + cvSaveImage(nam,hist_img); |
| 110 | + |
| 111 | + IplImage *bin_img2 = cvCreateImage(cvSize( source->width/2, source->height/2 ), IPL_DEPTH_8U, 1 ); |
| 112 | + int mn1=300,mx1=-1,avg1=0; |
| 113 | + for(int i=0;i<r/2;i++) |
| 114 | + { |
| 115 | + for(int j=c/2;j<c;j++) |
| 116 | + { |
| 117 | + CvScalar s; |
| 118 | + s = cvGet2D(bin_img,i,j); |
| 119 | + if(s.val[0] < mn1) |
| 120 | + mn1 = s.val[0]; |
| 121 | + if(s.val[0] > mx1) |
| 122 | + mx1 = s.val[0]; |
| 123 | + avg1+=s.val[0]; |
| 124 | + cvSet2D(bin_img2,i,j-c/2,s); |
| 125 | + } |
| 126 | + } |
| 127 | + avg1 = 4*avg1/(r*c); |
| 128 | + |
| 129 | + hist_img = cvCreateImage(cvSize(300,300), 8, 1); |
| 130 | + cvSet( hist_img, cvScalarAll(255), 0 ); |
| 131 | + create_histogram_image(bin_img2, hist_img,nbins); |
| 132 | + n = sprintf(nam,"%s_2.%s",img_nm,img_typ); |
| 133 | + cvSaveImage(nam,hist_img); |
| 134 | + |
| 135 | + IplImage *bin_img3 = cvCreateImage(cvSize( source->width/2, source->height/2 ), IPL_DEPTH_8U, 1 ); |
| 136 | + int mn2=300,mx2=-1,avg2=0; |
| 137 | + for(int i=r/2;i<r;i++) |
| 138 | + { |
| 139 | + for(int j=0;j<c/2;j++) |
| 140 | + { |
| 141 | + CvScalar s; |
| 142 | + s = cvGet2D(bin_img,i,j); |
| 143 | + if(s.val[0] < mn2) |
| 144 | + mn2 = s.val[0]; |
| 145 | + if(s.val[0] > mx2) |
| 146 | + mx2 = s.val[0]; |
| 147 | + avg2+=s.val[0]; |
| 148 | + cvSet2D(bin_img3,i-r/2,j,s); |
| 149 | + } |
| 150 | + } |
| 151 | + avg2 = 4*avg2/(r*c); |
| 152 | + |
| 153 | + hist_img = cvCreateImage(cvSize(300,300), 8, 1); |
| 154 | + cvSet( hist_img, cvScalarAll(255), 0 ); |
| 155 | + create_histogram_image(bin_img3, hist_img,nbins); |
| 156 | + n = sprintf(nam,"%s_3.%s",img_nm,img_typ); |
| 157 | + cvSaveImage(nam,hist_img); |
| 158 | + |
| 159 | + |
| 160 | + IplImage *bin_img4 = cvCreateImage(cvSize( source->width/2, source->height/2 ), IPL_DEPTH_8U, 1 ); |
| 161 | + int mn3=300,mx3=-1,avg3=0; |
| 162 | + for(int i=r/2;i<r;i++) |
| 163 | + { |
| 164 | + for(int j=c/2;j<c;j++) |
| 165 | + { |
| 166 | + CvScalar s; |
| 167 | + s = cvGet2D(bin_img,i,j); |
| 168 | + if(s.val[0] < mn3) |
| 169 | + mn3 = s.val[0]; |
| 170 | + if(s.val[0] > mx3) |
| 171 | + mx3 = s.val[0]; |
| 172 | + avg3+=s.val[0]; |
| 173 | + cvSet2D(bin_img4,i-r/2,j-c/2,s); |
| 174 | + } |
| 175 | + } |
| 176 | + avg3 = 4*avg3/(r*c); |
| 177 | + |
| 178 | + hist_img = cvCreateImage(cvSize(300,300), 8, 1); |
| 179 | + cvSet( hist_img, cvScalarAll(255), 0 ); |
| 180 | + create_histogram_image(bin_img4, hist_img,nbins); |
| 181 | + n = sprintf(nam,"%s_4.%s",img_nm,img_typ); |
| 182 | + cvSaveImage(nam,hist_img); |
| 183 | + |
| 184 | + cvReleaseImage( &source ); |
| 185 | + cvReleaseImage( &hist_img ); |
| 186 | + cvReleaseImage( &bin_img ); |
| 187 | + cvReleaseImage( &bin_img1 ); |
| 188 | + cvReleaseImage( &bin_img2 ); |
| 189 | + cvReleaseImage( &bin_img3 ); |
| 190 | + cvReleaseImage( &bin_img4 ); |
| 191 | + |
| 192 | + cout<<mn1<< " "<<mx1<< " "<<avg1<< " "<<mn2<< " "<<mx2<< " "<<avg2<< " "<<mn3<< " "<<mx3<< " "<<avg3<< " "<<mn4<< " "<<mx4<< " "<<avg4<<endl; |
| 193 | + return 0; |
| 194 | +} |
0 commit comments