图形三维变换

/* Note:Your choice is C IDE */
#include<math.h>
#include<graphics.h>
#include<conio.h>
#include<time.h>
#include<ctype.h>
#define  ZOOM_IN 0.9
#define  ZOOM_OUT 1.1

int turn1[3];
typedef struct{
	float x;
	float y;
	float z;
}point;

typedef struct{
	float x;
	float y;
}point2d;

typedef struct{
	float x;
	float y;
	float h;
	point biao[8];
}fanti;
void make_box(float x,float y,float h,fanti *p)
{ 
	p->x=x;
	p->y=y;
	p->h=h;

	p->biao[0].x=x/2;
	p->biao[0].y=y/2;
	p->biao[0].z=h/2;

	p->biao[1].x=-x/2;
	p->biao[1].y=y/2;
	p->biao[1].z=h/2;

	p->biao[2].x=-x/2;
	p->biao[2].y=-y/2;
	p->biao[2].z=h/2;

	p->biao[3].x=x/2;
	p->biao[3].y=-y/2;
	p->biao[3].z=h/2;

	p->biao[4].x=x/2;
	p->biao[4].y=y/2;
	p->biao[4].z=-h/2;

	p->biao[5].x=-x/2;
	p->biao[5].y=y/2;
	p->biao[5].z=-h/2;

	p->biao[6].x=-x/2;
	p->biao[6].y=-y/2;
	p->biao[6].z=-h/2;

	p->biao[7].x=x/2;
	p->biao[7].y=-y/2;
	p->biao[7].z=-h/2;

}
void trun2d(point *p,point2d *q)
{  
	q->x=p->x+p->z*cos(0.25);
	q->y=p->y+p->z*sin(0.25);
}
void initm(float mat[][4])
{
	int count ;
	for(count =0;count<4;count ++)
	{
		mat[count][0]=0;
		mat[count][1]=0;
		mat[count][2]=0;
		mat[count][3]=0;
		mat[count][count]=1;
	}return;
}
void transfrom(point *p,point *q,float tm[][4])
{
	float xu,yv,zw,h;
	xu=tm[0][0]*p->x+tm[1][0]*p->y+tm[2][0]*p->z+tm[3][0];
	yv=tm[0][1]*p->x+tm[1][1]*p->y+tm[2][1]*p->z+tm[3][1];
	zw=tm[0][2]*p->x+tm[1][2]*p->y+tm[2][2]*p->z+tm[3][2];
	q->x=xu;
	q->y=yv;
	q->z=zw;
	return;
}
void ratationx(point *p,float alfa,float tm[][4])
{ 
	float rad=0.0174532925;
	initm(tm);
	tm[1][1]=cos(rad*alfa);
	tm[1][2]=sin(rad*alfa);
	tm[2][1]=-tm[1][2];
	tm[2][2]=tm[1][1];
	return;
}
void rotationz(point *p,float alfa,float tm[][4])
{
	float rad=0.0174532925;
	initm(tm);
	tm[0][0]=cos(rad*alfa);
	tm[0][1]=sin(rad*alfa);
	tm[1][0]=-tm[0][1];
	tm[1][1]=tm[0][0];
	return;
}
void rotationy(point *p,float alfa,float tm[][4])
{ 
	float rad=0.0174532925;
	initm(tm);
	tm[0][0]=cos(rad*alfa);
	tm[2][0]=sin(rad*alfa);
	tm[0][2]=tm[2][0];
	tm[2][0]=tm[0][0];
	return;
}
void adjust(point *p,point *q)
{
	float t[4][4];
	switch(turn1[0])
	{case 1:
	rotationy(p,2,t);
	transfrom(p,q,t);
	break;
	case -1:
		rotationy(p,-2,t);
		transfrom(p,q,t);
		break;
	default:break;
	}

	switch(turn1[1])
	{case 1:
	rotationz(p,2,t);
	transfrom(p,q,t);
	break;
	case -1:
		rotationz(p,-2,t);
		transfrom(p,q,t);
		break;
	default:break;
	}

	switch(turn1[2])
	{case 1:
	q->x=ZOOM_IN*p->x;
	q->y=ZOOM_IN*p->y;
	q->z=ZOOM_IN*p->z;
	break;
	case -1:
		q->x=ZOOM_IN*p->x;
		q->y=ZOOM_IN*p->y;
		q->z=ZOOM_IN*p->z;
		break;
	default:break;
	}

}
void drawbox(fanti * p)
{
	point2d fan2d[8];
	int i;
	for(i=0;i<=7;i++)
	{
		adjust(&p->biao[i],&p->biao[i]);


		trun2d(&p->biao[i],&fan2d[i]);
		fan2d[i].x+=300;
		fan2d[i].y+=200;
	}
	clearviewport();
	/*outtext("\n -> :right
	\n <-:left
	\n^:up
	\n v:down");*/
	moveto(0,10);
	/*outtext("\n page up:zoom in
	\n space:Redraw
	\n Esc:exit");*/
	for(i=0;i<=3;i++)
	{
		if(i==3)
		{line(fan2d[i].x,fan2d[i].y,fan2d[0].x,fan2d[0].y);

		line(fan2d[i+4].x,fan2d[i+4].y,fan2d[4].x,fan2d[4].y);
		}
		else {
			line(fan2d[i].x,fan2d[i].y,fan2d[i+1].x,fan2d[i+1].y);
			line(fan2d[i+4].x,fan2d[i+4].y,fan2d[i+5].x,fan2d[i+5].y);
		}
		line(fan2d[i].x,fan2d[i].y,fan2d[i+4].x,fan2d[i+4].y);
	}
}

void main()
{
	int gd=DETECT,gm,i,j;
	char key;
	float x,y,h;
	fanti al;
	x=100;y=100;h=100;
	initgraph(&gd,&gm," ");
	make_box(x,y,h,&al);
	drawbox(&al);
	for(;;)
	{
		turn1[0]=0;
		turn1[1]=0;
		turn1[2]=0;
		key=getch();
		switch(key)
		{
		case 77: turn1[0]=1;
			break;
		case 75: turn1[0]=-1;
			break;
		case 72: turn1[1]=1;
			break;
		case 80: turn1[1]=-1;
			break;
		case 73: turn1[2]=1;
			break;
		case 81: turn1[2]=-1;
			break;
		case 32:  make_box(x,y,h,&al);
			break;
		case 27: exit();
			break;
		default:  key=0;
			break;
		}
		if(key!=0)
			drawbox(&al);
	}
	/*closegraph();*/
}

//注：按上下左右键变换，空格键还原，ESC键退出