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 ;
}
Esempio n. 3
0
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;
}
Esempio n. 4
0
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;
}
Esempio n. 5
0
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();

}