Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- /*
- * GLUT Shapes Demo
- *
- * Written by Nigel Stewart November 2003
- *
- * This program is test harness for the sphere, cone
- * and torus shapes in GLUT.
- *
- * Spinning wireframe and smooth shaded shapes are
- * displayed until the ESC or q key is pressed. The
- * number of geometry stacks and slices can be adjusted
- * using the + and - keys.
- */
- #include <windows.h>
- #ifdef __APPLE__
- #include <GLUT/glut.h>
- #else
- #include <GL/glut.h>
- #endif
- #include <stdlib.h>
- #include <fstream>
- #include <iostream>
- #include <math.h>
- using namespace std;
- static int slices = 16;
- static int stacks = 16;
- double gamma=1.0;
- const int nx = 1339;
- const int ny = 2000;
- double rOrg[nx][ny];
- double rO0[nx][ny];
- double rNew[nx][ny];
- double gOrg[nx][ny];
- double gO0[nx][ny];
- double gNew[nx][ny];
- double bOrg[nx][ny];
- double bO0[nx][ny];
- double bNew[nx][ny];
- double mUsK[3][3]={{-1,-1,-1},
- {-1,9,-1},
- {-1,-1,-1}};
- double mUs[3][3]={{0,-1,0},
- {0,1,0},
- {0,0,0}};
- double mUs2[3][3]={{0,0,0},
- {-1,1,0},
- {0,0,0}};
- double mUw[3][3]={{-1,1,1},
- {-1,-2,1},
- {-1,1,1}};
- double pS=1;
- int czytanie(const string& nazwaPliku)
- {
- int szer,wys;
- ifstream plik(nazwaPliku.c_str());
- plik>>szer>>wys;
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- plik>>rOrg[i][j];
- plik>>gOrg[i][j];
- plik>>bOrg[i][j];
- rOrg[i][j]/=255.0;
- gOrg[i][j]/=255.0;
- bOrg[i][j]/=255.0;
- // cout<<rOrg[i][j]<<"\t"<<gOrg[i][j]<<"\t"<<bOrg[i][j]<<endl;
- }
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rNew[i][j]=rOrg[i][j];
- gNew[i][j]=gOrg[i][j];
- bNew[i][j]=bOrg[i][j];
- rO0[i][j]=rOrg[i][j];
- gO0[i][j]=gOrg[i][j];
- bO0[i][j]=bOrg[i][j];
- // cout<<rOrg[i][j]<<"\t"<<gOrg[i][j]<<"\t"<<bOrg[i][j]<<endl;
- }
- cout<<"szer="<<szer<<"\twys="<<wys<<endl;
- plik.close();
- return 0;
- }
- /* GLUT callback Handlers */
- static void resize(int width, int height)
- {
- const float ar = (float) width / (float) height;
- glViewport(0, 0, width, height);
- glMatrixMode(GL_PROJECTION);
- glLoadIdentity();
- glOrtho(0,ny,0,nx,-5,5);
- if((double)width/ny>(double)height/nx)
- pS=(double)width/ny;
- else
- pS=(double)height/nx;
- glMatrixMode(GL_MODELVIEW);
- glLoadIdentity() ;
- }
- static void display(void)
- {
- const double t = glutGet(GLUT_ELAPSED_TIME) / 1000.0;
- const double a = t*90.0;
- glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
- glPointSize(pS);
- glBegin(GL_POINTS);
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- glColor3d(rNew[i][j],gNew[i][j],bNew[i][j]);
- glVertex3d(j,i,0);
- }
- glEnd();
- glutSwapBuffers();
- }
- void wypisz()
- {
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- cout<<rOrg[i][j]<<" : ";
- cout<<gOrg[i][j]<<" : ";
- cout<<bOrg[i][j]<<" : ";
- }
- }
- void maxymalny()
- {
- double rmax =0.0,gmax=0.0,bmax=0.0;
- for(int i=1;i<nx-1;i++)
- for(int j=1;j<ny-1;j++)
- {
- rmax =0.0;
- gmax =0.0;
- bmax =0.0;
- for(int k=0;k<3;k++)
- for(int l=0;l<3;l++)
- {
- if(rmax<rOrg[i+k-1][j+l-1])
- rmax=rOrg[i+k-1][j+l-1];
- if(gmax<gOrg[i+k-1][j+l-1])
- gmax=gOrg[i+k-1][j+l-1];
- if(bmax<bOrg[i+k-1][j+l-1])
- bmax=bOrg[i+k-1][j+l-1];
- }
- rNew[i][j]=rmax;
- gNew[i][j]=gmax;
- bNew[i][j]=bmax;
- }
- }
- void minimalny()
- {
- double rmin =0.0,gmin=0.0,bmin=0.0;
- for(int i=1;i<nx-1;i++)
- for(int j=1;j<ny-1;j++)
- {
- rmin =1.0;
- gmin =1.0;
- bmin =1.0;
- for(int k=0;k<3;k++)
- for(int l=0;l<3;l++)
- {
- if(rmin>rOrg[i+k-1][j+l-1])
- rmin=rOrg[i+k-1][j+l-1];
- if(gmin>gOrg[i+k-1][j+l-1])
- gmin=gOrg[i+k-1][j+l-1];
- if(bmin>bOrg[i+k-1][j+l-1])
- bmin=bOrg[i+k-1][j+l-1];
- }
- rNew[i][j]=rmin;
- gNew[i][j]=gmin;
- bNew[i][j]=bmin;
- }
- }
- void funkcja_gamma()
- {
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rNew[i][j]=pow(rOrg[i][j],gamma);
- gNew[i][j]=pow(gOrg[i][j],gamma);
- bNew[i][j]=pow(bOrg[i][j],gamma);
- }
- }
- void bubblesort ()
- {
- }
- static void key(unsigned char key, int x, int y)
- {
- switch (key)
- {
- case 'q':
- exit(0);
- break;
- //reset filtrów
- case '0':
- {
- for(int i=0; i<nx; i++)
- for(int j=0; j<ny; j++)
- {
- rNew[i][j]=rO0[i][j];
- gNew[i][j]=gO0[i][j];
- bNew[i][j]=bO0[i][j];
- rOrg[i][j]=rO0[i][j];
- gOrg[i][j]=rO0[i][j];
- bOrg[i][j]=rO0[i][j];
- }
- gamma=1.0;
- }
- break;
- //Kolor czerwony
- case 'r':
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rNew[i][j]+=0.1;
- }
- break;
- //Kolor zielony
- case 'g':
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- gNew[i][j]+=0.1;
- }
- //Kolor niebieski
- case 'b':
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- bNew[i][j]+=0.1;
- }
- break;
- //rozjaœnienie
- case '1':
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rNew[i][j]+=0.1;
- gNew[i][j]+=0.1;
- bNew[i][j]+=0.1;
- }
- break;
- //korekcja + gamma
- case '3':
- {
- funkcja_gamma();
- gamma+=0.1;
- cout<<"gamma = "<<gamma<<endl;
- }
- break;
- //korekcja - gamma
- case '4':
- {
- funkcja_gamma();
- gamma-=0.1;
- cout<<"gamma = "<<gamma<<endl;
- }
- break;
- //przyciemnienie
- case '2':
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rNew[i][j]-=0.1;
- gNew[i][j]-=0.1;
- bNew[i][j]-=0.1;
- }
- break;
- //filtr krawedziowy (przesuwajacy pionowy)
- case '5':
- {
- double sumaWg,rSum=0.0,gSum=0.0,bSum=0.0;
- for(int i=0; i<3; i++)
- {
- for(int j=0; j<3; j++)
- {
- sumaWg +=mUs[i][j];
- }
- }
- cout<<sumaWg;
- if (sumaWg==0)sumaWg=1;
- for(int i=1; i<nx-1; i++)
- for(int j=1; j<ny-1; j++)
- {
- for(int k=0; k<3; k++)
- for(int l=0; l<3; l++)
- {
- rSum+=(mUs[k][l]*rOrg[i+k-1][j+l-1]);
- gSum+=(mUs[k][l]*gOrg[i+k-1][j+l-1]);
- bSum+=(mUs[k][l]*bOrg[i+k-1][j+l-1]);
- }
- rNew[i][j]=(rSum/sumaWg);
- gNew[i][j]=(gSum/sumaWg);
- bNew[i][j]=(bSum/sumaWg);
- rSum=0.0;
- gSum=0.0;
- bSum=0.0;
- }
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rOrg[i][j]=rNew[i][j];
- gOrg[i][j]=gNew[i][j];
- bOrg[i][j]=bNew[i][j];
- }
- }
- break;
- //filtr krawedziowy (przesuwajacy poziomy)
- case '6':
- {
- double sumaWg,rSum=0.0,gSum=0.0,bSum=0.0;
- for(int i=0;i<3;i++)
- {
- for(int j=0;j<3;j++)
- {
- sumaWg +=mUs2[i][j];
- }
- }
- cout<<sumaWg;
- if (sumaWg==0)sumaWg=1;
- for(int i=1;i<nx-1;i++)
- for(int j=1;j<ny-1;j++)
- {
- for(int k=0; k<3; k++)
- for(int l=0; l<3; l++)
- {
- rSum+=(mUs2[k][l]*rOrg[i+k-1][j+l-1]);
- gSum+=(mUs2[k][l]*gOrg[i+k-1][j+l-1]);
- bSum+=(mUs2[k][l]*bOrg[i+k-1][j+l-1]);
- }
- rNew[i][j]=(rSum/sumaWg);
- gNew[i][j]=(gSum/sumaWg);
- bNew[i][j]=(bSum/sumaWg);
- rSum=0.0;
- gSum=0.0;
- bSum=0.0;
- }
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rOrg[i][j]=rNew[i][j];
- gOrg[i][j]=gNew[i][j];
- bOrg[i][j]=bNew[i][j];
- }
- }
- break;
- //filtr wschód
- case 'w':
- {
- double sumaWg,rSum=0.0,gSum=0.0,bSum=0.0;
- for(int i=0;i<3;i++)
- {
- for(int j=0;j<3;j++)
- {
- sumaWg +=mUw[i][j];
- }
- }
- //cout<<sumaWg;
- if (sumaWg==0)sumaWg=1;
- for(int i=1;i<nx-1;i++)
- for(int j=1;j<ny-1;j++)
- {
- for(int k=0; k<3; k++)
- for(int l=0; l<3; l++)
- {
- rSum+=(mUw[k][l]*rOrg[i+k-1][j+l-1]);
- gSum+=(mUw[k][l]*gOrg[i+k-1][j+l-1]);
- bSum+=(mUw[k][l]*bOrg[i+k-1][j+l-1]);
- }
- rNew[i][j]=(rSum/sumaWg);
- gNew[i][j]=(gSum/sumaWg);
- bNew[i][j]=(bSum/sumaWg);
- rSum=0.0;
- gSum=0.0;
- bSum=0.0;
- }
- for(int i=0;i<nx;i++)
- for(int j=0;j<ny;j++)
- {
- rOrg[i][j]=rNew[i][j];
- gOrg[i][j]=gNew[i][j];
- bOrg[i][j]=bNew[i][j];
- }
- }
- break;
- case 'm':
- {
- maxymalny();
- }
- break;
- case 'n':
- {
- minimalny();
- }
- break;
- }
- glutPostRedisplay();
- }
- static void idle(void)
- {
- glutPostRedisplay();
- }
- const GLfloat light_ambient[] = { 0.0f, 0.0f, 0.0f, 1.0f };
- const GLfloat light_diffuse[] = { 1.0f, 1.0f, 1.0f, 1.0f };
- const GLfloat light_specular[] = { 1.0f, 1.0f, 1.0f, 1.0f };
- const GLfloat light_position[] = { 2.0f, 5.0f, 5.0f, 0.0f };
- const GLfloat mat_ambient[] = { 0.7f, 0.7f, 0.7f, 1.0f };
- const GLfloat mat_diffuse[] = { 0.8f, 0.8f, 0.8f, 1.0f };
- const GLfloat mat_specular[] = { 1.0f, 1.0f, 1.0f, 1.0f };
- const GLfloat high_shininess[] = { 100.0f };
- /* Program entry point */
- int main(int argc, char *argv[])
- {
- cout<<"'q' ->EXIT"<<endl;
- cout<<"'0' ->RESET FILTROW"<<endl;
- cout<<"'r' ->PODNOSZENIE POZIOMU CZERWIENI"<<endl;
- cout<<"'g' ->PODNOSZENIE POZIOMU ZILENI"<<endl;
- cout<<"'b' ->PODNOSZENIE POZIOMU NIEBIESKIEGO"<<endl;
- cout<<"'1' ->ROZJASNIANIE"<<endl;
- cout<<"'2' ->PRZYCIEMNIANIE"<<endl;
- cout<<"'3' ->KORELACJA GAMMA +"<<endl;
- cout<<"'4' ->KORELACJA GAMMA -"<<endl;
- cout<<"'5' ->FILTR PRZESUWAJACY PIONOWY"<<endl;
- cout<<"'6' ->FILTR PRZESUWAJACY POZIOMY"<<endl;
- cout<<"'w' ->FILTR WSCHOD"<<endl;
- cout<<"'m' ->FILTR MAXIMUM"<<endl;
- cout<<"'n' ->FILTR MINIMUM"<<endl;
- czytanie("c:\\tmp\\zd5.txt");
- glutInit(&argc, argv);
- glutInitWindowSize(ny,nx);
- glutInitWindowPosition(10,10);
- glutInitDisplayMode(GLUT_RGB | GLUT_DOUBLE | GLUT_DEPTH);
- glutCreateWindow("GLUT Shapes");
- glutReshapeFunc(resize);
- glutDisplayFunc(display);
- glutKeyboardFunc(key);
- glutIdleFunc(idle);
- glClearColor(1,1,1,1);
- glEnable(GL_CULL_FACE);
- glCullFace(GL_BACK);
- glEnable(GL_DEPTH_TEST);
- glDepthFunc(GL_LESS);
- glEnable(GL_LIGHT0);
- glEnable(GL_NORMALIZE);
- glEnable(GL_COLOR_MATERIAL);
- glEnable(GL_LIGHTING);
- glLightfv(GL_LIGHT0, GL_AMBIENT, light_ambient);
- glLightfv(GL_LIGHT0, GL_DIFFUSE, light_diffuse);
- glLightfv(GL_LIGHT0, GL_SPECULAR, light_specular);
- glLightfv(GL_LIGHT0, GL_POSITION, light_position);
- glMaterialfv(GL_FRONT, GL_AMBIENT, mat_ambient);
- glMaterialfv(GL_FRONT, GL_DIFFUSE, mat_diffuse);
- glMaterialfv(GL_FRONT, GL_SPECULAR, mat_specular);
- glMaterialfv(GL_FRONT, GL_SHININESS, high_shininess);
- glutMainLoop();
- return EXIT_SUCCESS;
- }
Add Comment
Please, Sign In to add comment