ait-survey

normal2psuedo

Jan 31st, 2015
288
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.83 KB | None | 0 0
  1. //*********************************************
  2.  
  3. //Psuedo Color
  4.  
  5.  
  6. //***********************************************//
  7.  
  8.  
  9. #include <stdio.h>
  10. #include <stdlib.h>
  11. #include <math.h>
  12.  
  13. int main(int argc, char *argv[])
  14. {
  15.     //Decralation of variables
  16.     FILE *fpi,*fpo;
  17.     float x,y,z,nx,ny,nz;
  18.     float fnx,fny,fnz;
  19.     float pr,pg,pb;
  20.     unsigned char upr,upg,upb;
  21.     float wr,wg,wb;
  22.     float fr,fg,fb;
  23.     unsigned char ar,ag,ab;
  24.     unsigned char unx,uny,unz,r,g,b;
  25.     char point[4];
  26.  
  27.     //Open CSV file
  28.     if((fpi=fopen(argv[1],"r"))==NULL)
  29.     {
  30.         printf ("The file can't be opened. The program is exit.\n");
  31.  
  32.         return 0;
  33.     }
  34.  
  35.     //Open VTK file for save
  36.     if((fpo=fopen(argv[2],"w"))==NULL)
  37.     {
  38.         printf("The file can't be opened. The program is exit.\n");
  39.  
  40.         return 0;
  41.     }
  42.  
  43.     //Initialization of paramaters
  44.     ar=255;
  45.     ag=255;
  46.     ab=255;
  47.     wr=2*3.14/255;
  48.     wg=3*3.14/255;
  49.     wb=4*3.14/255;
  50.     fr=0;
  51.     fg=0.3*3.14;
  52.     fb=0.6*3.14;
  53.  
  54.     while(fscanf(fpi,"%f %f %f %d %d %d %f %f %f",&x,&y,&z,&r,&g,&b,&nx,&ny,&nz)!=EOF)
  55.     {
  56.             //printf("%f %f %f\n",x,y,z);
  57.         fnx=(nx+1)*255/2;
  58.         fny=(ny+1)*255/2;
  59.         fnz=(nz+1)*255/2;
  60.             unx=(unsigned char)fnx;
  61.             uny=(unsigned char)fny;
  62.             unz=(unsigned char)fnz;
  63.         //printf("wr=%f\n",sin(wr*unx+fr));
  64.         //printf("wg=%f\n",sin(wg*uny+fg));
  65.         //printf("wb=%f\n",sin(wb*unz+fb));
  66.         //printf("wrabs=%f\n",fabs(sin(wr*unx+fr)));
  67.         //printf("wgabs=%f\n",fabs(sin(wg*uny+fg)));
  68.         //printf("wbabs=%f\n",fabs(sin(wb*unz+fb)));
  69.             pr=ar*fabs(sin(wr*unx+fr));
  70.             pg=ag*fabs(sin(wg*uny+fg));
  71.             pb=ab*fabs(sin(wb*unz+fb));
  72.             upr=(char)pr;
  73.             upg=(char)pg;
  74.             upb=(char)pb;
  75.         //printf("%f %f %f %f %f %f\n",x,y,z,pr,pg,pb);
  76.         //printf("%d %d %d\n",upr,upg,upb);
  77.             fprintf(fpo,"%f %f %f %d %d %d\n",x,y,z,upr,upg,upb);
  78.  
  79.     }
  80.  
  81.         fclose(fpi);
  82.         fclose(fpo);
  83.  
  84. }
Advertisement
Add Comment
Please, Sign In to add comment