计算机图形学 图形几何变换的实现

计算机图形学图形几何变换的实现

————————————————————————————————作者：————————————————————————————————日期：

实验五图形几何变换的实现

班级08信计2 学号89姓名徐阳分数

一、实验目的和要求：

1、掌握理解二维、三维变换的数学知识、变换原理、变换种类、变换方法；进一步理解采用齐次坐标进行二维、三维变换的必要性；利用Turboc实现二维、三维图形的基本变换和复合变换。

二、实验内容：

1、理解采用齐次坐标进行图形变换的必要性——变换的连续性，使复合变换得以实现。

2、掌握二维、三维图形基本变换（平移、缩放、对称、旋转、错切）的原理及数学公式。

3、利用Turboc实现二维、三维图形的基本变换、复合变换，在屏幕上显示变换过程或变换结果。

三、实验结果分析：

程序代码如下：

/*二维图形（直线）平移变换*/

#include

#include

#include

main()

{int x0,y0,x1,y1,i,j;

int a[3][3];

char key;

for(i=0;i<3;i++)

for(j=0;j<3;j++)

a[i][j]=0;

for(i=0;i<3;i++)

a[i][i]=1;

int graphdriver=DETECT;

int graphmode=0;

initgraph(&graphdriver,&graphmode," ");

cleardevice();

x0=250;y0=120;x1=350;y1=220;

line(x0,y0,x1,y1);

for( ; ;)

{outtextxy(100,400,"<-:left->:right^:up v:down Esc->exit");

key=getch();

switch(key)

{

case 75: a[2][0]=-10;

break;

case 77: a[2][0]=10;

break;

case 72: a[2][1]=-10;

break;

case 80: a[2][1]=10;

break;

case 27: exit(0);

break;

}

x0=x0*a[0][0]+y0*a[1][0]+a[2][0];

y0=x0*a[0][1]+y0*a[1][1]+a[2][1];

x1=x1*a[0][0]+y1*a[1][0]+a[2][0];

y1=x1*a[0][1]+y1*a[1][1]+a[2][1];

cleardevice();

line(x0,y0,x1,y1);

}

closegraph();

}

2.三维图形（立方体）旋转和比例变换的算法：

#include

#include

#include

#include

#include

#include

#define ZOOM_IN 0.9

#define ZOOM_OUT 1.1

int turn1[3]; /* [0] rx,[1] ry,[3] zoom*/

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.;