北航計算機課程-科學(xué)計算可視化算法-大作業(yè)-三維球體體數(shù)據(jù)生成及繪制

北航計算機課程科學(xué)計算可視化算法大作業(yè)masikkk#include<iostream>#include<fstream>#include<vector>#include<math.h>#include<windows.h>#include<gl/gl.h>#include<gl/glu.h>#include<gl/glut.h>#include<gl/glaux.h>#include<math.h>usingnamespacestd;#defineSTEP0.02//采樣間隔//體素結(jié)構(gòu)體typedefstructVoxel{floatx;floaty;floatz;intvalue;}Voxel;vector<Voxel>VolumeData;//體素數(shù)組GLfloath;//視點z軸坐標(biāo)//求兩點間距離floatDistance(floata[3],floatb[3]){returnsqrt((a[0]-b[0])*(a[0]-b[0])+(a[1]-b[1])*(a[1]-b[1])+(a[2]-b[2])*(a[2]-b[2]));}//生成球的體數(shù)據(jù)并保存到文件，參數(shù)是半徑和球心voidGenerateVolumeData(floatRadi,floatCent[3]){floatRadius=Radi;//半徑floatCenter[3];//球心floatBoundBoxStart[3];//包圍盒起點坐標(biāo)，左下角頂點floatBoundBoxEnd[3];//包圍盒終點坐標(biāo)，右上角頂點//初始化for(intm=0;m<3;m++){Center[m]=Cent[m];//球心BoundBoxStart[m]=Center[m]-Radius;//包圍盒起點BoundBoxEnd[m]=Center[m]+Radius+0.02;//包圍盒終點}ofstreamfout_VolumeData;fout_VolumeData.open("VolumeData.txt");//保存體數(shù)據(jù)的文件//遍歷包圍盒內(nèi)的點for(floati=BoundBoxStart[0];i<=BoundBoxEnd[0];i+=STEP){for(floatj=BoundBoxStart[1];j<=BoundBoxEnd[1];j+=STEP){for(floatk=BoundBoxStart[2];k<=BoundBoxEnd[2];k+=STEP*2){floatpoint[3]={i,j,k};floatdist=Distance(point,Center);//點point和球心之間的距離fout_VolumeData.setf(ios::fixed);//不以科學(xué)技術(shù)法顯示fout_VolumeData.precision(5);//設(shè)置精度fout_VolumeData<<i<<""<<j<<""<<k<<"";//輸出坐標(biāo)if(fabs(dist-Radius)<0.01)//點point在球上，輸出0fout_VolumeData<<0;elseif(dist<Radius)//點point在球內(nèi)，輸出-1fout_VolumeData<<-1;else//if(dist>Radius)//點point在球外，輸出1fout_VolumeData<<1;fout_VolumeData<<endl;}}}fout_VolumeData.close();}//窗口初始化和大小改變時，調(diào)用此函數(shù)voidCALLBACKreshape(GLsizeiw,GLsizeih){glMatrixMode(GL_PROJECTION);//設(shè)置當(dāng)前矩陣為投影變換矩陣glLoadIdentity();//初始化當(dāng)前矩陣gluPerspective(20,1,0,3);//設(shè)置透視投影矩陣glMatrixMode(GL_MODELVIEW);//設(shè)置當(dāng)前矩陣為模式變換矩陣glViewport(0,0,w,h);//設(shè)置視區(qū)變換}//顯示函數(shù)voidCALLBACKdisplay(){glClearColor(0,1,1,1);//設(shè)置窗口背景顏色glClear(GL_COLOR_BUFFER_BIT);//清顏色緩沖區(qū)glLoadIdentity();//重置變換矩陣gluLookAt(5,5,h,0,0,0,0,1,0);//設(shè)置視點坐標(biāo)glColor3f(1,0,0);//設(shè)置前景色glBegin(GL_POINTS);//遍歷體素數(shù)組vector<Voxel>::iteratoriter=VolumeData.begin();for(;iter!=VolumeData.end();iter++){if(iter->value==0)glVertex3f(iter->x,iter->y,iter->z);}glEnd();glFlush();//刷新緩沖區(qū)}//上方向鍵voidCALLBACKUp(){h+=0.1;}//下方向鍵voidCALLBACKDown(){h-=0.1;}voidmain(){floatRadius=1;//半徑floatCenter[3]={0,0,0};//球心GenerateVolumeData(Radius,Center);//生成體數(shù)據(jù)并保存到文件ifstreamfin;fin.open("VolumeData.txt");if(!fin)cout<<"無法翻開體數(shù)據(jù)文件"<<endl;//將體數(shù)據(jù)讀入體素數(shù)組while(!fin.eof()){Voxelvox;fin>>vox.x;fin>>vox.y;fin>>vox.z;fin>>vox.value;VolumeData.push_back(vox);//插入數(shù)組}auxInitDisplayMode(AUX_SINGLE|AUX_RGBA);//設(shè)置窗口顯示模式auxInitPosition(0,0,500,500);//設(shè)置窗口位置auxInitWindow("三維球體");//窗口標(biāo)題glShadeModel(GL_FLA