void image_print3(){ int r,c; for(r=0;r<pic3.width;r++) { for(c=0;c<pic3.height;c++) { glBegin(GL_POINTS); float rr=(float)(pic3.pixel(r,c,0)); float gg=(float)(pic3.pixel(r,c,1)); float bb=(float)(pic3.pixel(r,c,2)); glColor3f(rr/255,gg/255.0,bb/255.0); glVertex2d(r,c); glEnd(); } } glTranslatef(0,-250,0); for(r=0;r<pic3.width;r++) { for(c=0;c<pic3.height;c++) { glBegin(GL_POINTS); double magn=(F_real2[r][c]*F_real2[r][c]); magn+=(F_img2[r][c]*F_img2[r][c]); magn=sqrt(magn); float rr=(float)magn; rr/=255.0; glColor3f(rr,rr,rr); glVertex2d(r,c); glEnd(); } } glTranslatef(0,250,0); }
void fourier_transform4(){ int i,j,n,m,k,l; // take input . . . f(i,j) n=pic4.height; m=pic4.width; for(i=0;i<pic4.height;i++) { for(j=0;j<pic4.width;j++) { float rr=(float)(pic4.pixel(i,j,0)); f_real[i][j]=rr; f_img[i][j]=0; if((i+j)%2){ f_real[i][j]=-rr; f_img[i][j]=-0; } } } // calculating . . . P(i,j); for(j=0;j<m;j++){ int k; for (k = 0; k < n; k++) { /* For each output element */ double sumreal = 0; double sumimag = 0; int t; for (t = 0; t < n; t++) { /* For each input element */ double angle = 2 * M_PI * t * k / (double)n; sumreal += f_real[t][j] * cos(angle) + f_img[t][j] * sin(angle); sumimag += -f_real[t][j] * sin(angle) + f_img[t][j] * cos(angle); } P_real[k][j] = sumreal; P_img[k][j] = sumimag; } } // taking Input . . . F(i,j); for(i=0;i<n;i++){ int k; for (k = 0; k < m; k++) { /* For each output element */ double sumreal = 0; double sumimag = 0; int t; for (t = 0; t < m; t++) { /* For each input element */ double angle = 2 * M_PI * t * k / (double)n; sumreal += P_real[i][t] * cos(angle) + P_img[i][t] * sin(angle); sumimag += -P_real[i][t] * sin(angle) + P_img[i][t] * cos(angle); } F_real3[i][k] = sumreal/(double)m; F_img3[i][k] = sumimag/(double)m; } } puts("FINISH"); return ; }
BMPError BMPLoad(std::string fname,BMPClass& bmp) { FILE * f=fopen(fname.c_str(),"rb"); //open for reading in binary mode if(!f) return BMPNOOPEN; char header[54]; fread(header,54,1,f); //read the 54bit main header if(header[0]!='B' || header[1]!='M') { fclose(f); return BMPNOTABITMAP; //all bitmaps should start "BM" } if(sizeof(int)==4) //Annoyingly I can't just assume this { bmp.width=*(int*)(header+18); bmp.height=*(int*)(header+22); } else { bmp.width=header[18]+256*header[19]+256*256*header[20]+256*256*256*header[21]; bmp.height=header[22]+256*header[23]+256*256*header[24]+256*256*256*header[25]; } bmp.allocateMem(); int bits=int(header[28]); //colourdepth int x,y,c; BYTE cols[256*4]; switch(bits) { case(24): fread(bmp.bytes,bmp.width*bmp.height*3,1,f); //24bit is easy for(x=0;x<bmp.width*bmp.height*3;x+=3) //except the format is BGR, grr { BYTE temp=bmp.bytes[x]; bmp.bytes[x]=bmp.bytes[x+2]; bmp.bytes[x+2]=temp; } break; case(8): fread(cols,256*4,1,f); //read colortable for(y=0;y<bmp.height;++y) //(Notice 4bytes/col for some reason) for(x=0;x<bmp.width;++x) { BYTE byte; fread(&byte,1,1,f); //just read byte for(int c=0;c<3;++c) bmp.pixel(x,y,c)=cols[byte*4+2-c]; //and look up in the table } break; case(4): fread(cols,16*4,1,f); for(y=0;y<256;++y) for(x=0;x<256;x+=2) { BYTE byte; fread(&byte,1,1,f); //as above, but need to exract two for(c=0;c<3;++c) //pixels from each byte bmp.pixel(x,y,c)=cols[byte/16*4+2-c]; for(c=0;c<3;++c) bmp.pixel(x+1,y,c)=cols[byte%16*4+2-c]; } break; //NOT CURRENTLY WORKING PROPERLY //case(1): // fread(cols,8,1,f); //colourtable // for(y=0;y<256;++y) // for(x=0;x<256;x+=8) // { // BYTE byte; // fread(&byte,1,1,f); // //The idea here is that every byte is eight pixels // //so I'm shifting the byte to the relevant position, then masking out // //all but the lowest bit in order to get the index into the colourtable. // for(int x2=0;x2<8;++x2) // { // for(int c=0;c<3;++c) // bmp.pixel(x+x2,y,c)=cols[((byte>>(7-x2))&1)*4+2-c]; // } // } // break; default: fclose(f); return BMPUNKNOWNFORMAT; } fclose(f); if(ferror(f)) return BMPFILEERROR; return BMPNOERROR; }
BMPError BMPLoad(std::string fname,BMPClass& bmp) { if(sizeof(int)!=4) return BMPBADINT; FILE* f=fopen(fname.c_str(),"rb"); //open for reading in binary mode if(!f) return BMPNOOPEN; char header[54]; fread(header,54,1,f); //read the 54bit main header if(header[0]!='B' || header[1]!='M') { fclose(f); return BMPNOTABITMAP; //all bitmaps should start "BM" } //it seems gimp sometimes makes its headers small, so we have to do this. hence all the fseeks int offset=*(unsigned int*)(header+10); bmp.width=*(int*)(header+18); bmp.height=*(int*)(header+22); //now the bitmap knows how big it is it can allocate its memory bmp.allocateMem(); int bits=int(header[28]); //colourdepth int x,y,c; BYTE cols[256*4]; //colourtable switch(bits) { case(24): fseek(f,offset,SEEK_SET); fread(bmp.bytes,bmp.width*bmp.height*3,1,f); //24bit is easy for(x=0;x<bmp.width*bmp.height*3;x+=3) //except the format is BGR, grr { BYTE temp=bmp.bytes[x]; bmp.bytes[x]=bmp.bytes[x+2]; bmp.bytes[x+2]=temp; } break; case(8): fread(cols,256*4,1,f); //read colortable fseek(f,offset,SEEK_SET); for(y=0;y<bmp.height;++y) //(Notice 4bytes/col for some reason) for(x=0;x<bmp.width;++x) { BYTE byte; fread(&byte,1,1,f); //just read byte for(int c=0;c<3;++c) bmp.pixel(x,y,c)=cols[byte*4+2-c]; //and look up in the table } break; case(4): fread(cols,16*4,1,f); fseek(f,offset,SEEK_SET); for(y=0;y<256;++y) for(x=0;x<256;x+=2) { BYTE byte; fread(&byte,1,1,f); //as above, but need to exract two for(c=0;c<3;++c) //pixels from each byte bmp.pixel(x,y,c)=cols[byte/16*4+2-c]; for(c=0;c<3;++c) bmp.pixel(x+1,y,c)=cols[byte%16*4+2-c]; } break; case(1): fread(cols,8,1,f); fseek(f,offset,SEEK_SET); for(y=0;y<bmp.height;++y) for(x=0;x<bmp.width;x+=8) { BYTE byte; fread(&byte,1,1,f); //Every byte is eight pixels //so I'm shifting the byte to the relevant position, then masking out //all but the lowest bit in order to get the index into the colourtable. for(int x2=0;x2<8;++x2) for(int c=0;c<3;++c) bmp.pixel(x+x2,y,c)=cols[((byte>>(7-x2))&1)*4+2-c]; } break; default: fclose(f); return BMPUNKNOWNFORMAT; } if(ferror(f)) { fclose(f); return BMPFILEERROR; } fclose(f); return BMPNOERROR; }
void mainmenu(){ /// have to change the color (sajid) //puts("MainMenu"); glBegin( GL_POINTS ); glColor3f(1,0,0); for(int i=0;i<=100;i++){ for(int j=0;j<=100;j++){ glVertex2d( i,j); } } glEnd(); glColor3f(0,0,0); glPushMatrix(); glLoadIdentity(); glTranslatef(20,100,0); //texture maping of left field int r,c,r1; double p1,c1; for(r=0;r<menu5.width;r++) { for(c=0;c<menu5.height;c++) { glBegin(GL_POINTS); float rr=(float)(menu5.pixel(r,c,0)); float gg=(float)(menu5.pixel(r,c,1)); float bb=(float)(menu5.pixel(r,c,2)); glColor3f(rr/255,gg/255.0,bb/255.0); if(!print_flag && j<30) printf("%3.0f,%3.0f,%3.0f, ",rr,gg,bb); glVertex2d(r,c); glEnd(); } if(!print_flag)puts(""); } if(!print_flag) printf("%d %d\n",r,c); print_flag=1; glPopMatrix(); glPushMatrix(); glLoadIdentity(); glTranslatef(250,100,0); for(r=0,r1=0;r<=menu5.width;r++,r1+=1) { for(c=0,c1=0;c<=menu5.height;c++,c1+=1.55) { glBegin(GL_POINTS); float rr=(float)(menu5.pixel(r,c,0)); float gg=(float)(menu5.pixel(r,c,1)); float bb=(float)(menu5.pixel(r,c,2)); glColor3f(rr/255,gg/255.0,bb/255.0); for(p1=0;p1<=2;p1+=.25){ glVertex2d(r1+p1,c1+p1); } glEnd(); } } glPopMatrix(); }