北航计算机课程-科学计算可视化算法-大作业-三维球体体数据生成及绘制.docxVIP

北航计算机课程

科学计算可视化算法

大作业

masikkk

#includeiostream

#includefstream

#includevector

#includemath.h

#includewindows.h

#includegl/gl.h

#includegl/glu.h

#includegl/glut.h

#includegl/glaux.h

#includemath.h

usingnamespacestd;

#defineSTEP0.02//采样间隔

//体素结构体

typedefstructVoxel

{

floatx;

floaty;

floatz;

intvalue;

}Voxel;

vectorVoxelVolumeData;//体素数组

GLfloath;//视点z轴坐标

//求两点间距离

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]));

}

//生成球的体数据并保存到文件，参数是半径和球心

voidGenerateVolumeData(floatRadi,floatCent[3])

{

floatRadius=Radi;//半径

floatCenter[3];//球心

floatBoundBoxStart[3];//包围盒起点坐标，左下角顶点

floatBoundBoxEnd[3];//包围盒终点坐标，右上角顶点

//初始化

for(intm=0;m3;m++)

{

Center[m]=Cent[m];//球心

BoundBoxStart[m]=Center[m]-Radius;//包围盒起点

BoundBoxEnd[m]=Center[m]+Radius+0.02;//包围盒终点

}

ofstreamfout_VolumeData;

fout_VolumeData.open(VolumeData.txt);//保存体数据的文件

//遍历包围盒内的点

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);//不以科学技术法显示

fout_VolumeData.precision(5);//设置精度

fout_VolumeDataijk;//输出坐标

if(fabs(dist-Radius)0.01)//点point在球上，输出0

fout_VolumeData0;

elseif(distRadius)//点point在球内，输出-1

fout_VolumeData-1;

else//if(distRadius)//点point在球外，输出1

fout_VolumeData1;

fout_VolumeDataendl;

}

}

}

fout_VolumeData.close();

}

//窗口初始化和大小改变时，调用此函数

voidCALLBACKreshape(GLsizeiw,GLsizeih)

{

glMatrixMode(GL_PROJECTION);//设置当前矩阵为投影变换矩阵

glLoadIdentity();//初始化当前矩阵

gluPerspective(20,1,0,3);//设置透视投影矩阵

glMatrixMode(GL_MODELVIEW);//设置当前矩阵为模式变换矩阵

glViewport(0,0,w,h);//设置视区变换

}

//显示函数

voidCALLBACKdisplay()

{

glClearColor(0,1,1,1);//设置窗口背景颜色

glClear(GL_COLOR_BUFFER_BIT);//清颜色缓冲区

glLoadIdentity();//重置变换矩阵

gluLookAt(5,5,h,0,0,0,0,1,0);//设置视点坐标

glColor3