(计算机图形学)关于任意直线的对称变换

实验3：关于任意直线的对称变换

实验类型：验证、设计

所需时间：3学时

主要实验内容及要求：

对于任意直线的二维图形对称变化的实验，要求输入的直线是任意直线，直线的端点只能由键盘输入或者鼠标拾取，要做对称变换的图形也是一个任意图形（至少应是一个任意多边形）。

对称变换，先分析如何使用一系列简单变换来构造题目要求的复合变换。本体要实现的变换可以用如下一组变换组合来实现：

①将直线任一点移至与坐标原点重合

②将平移后的直线绕原点旋转至与某一坐标轴重合

③将题目要求的对称变换转为实现已知图形关于上述坐标轴的对称变换

④按逆序求上述①、②变换的逆变换

⑤将上述矩阵依次相乘得到最终的复合变换矩阵

则某一多边形关于任意直线的对称变换就转变为将该多边形的各顶点与上述求得的复合变换进行矩阵乘法，求得变换后的新多边形的各个顶点坐标。

根据上述流程，编程实现，并测试程序功能。

源代码：

#include

#include

using namespace std;

void Initial(void)

{

glClearColor(1.0f,1.0f,1.0f,1.0f);

glMatrixMode(GL_PROJECTION);

gluOrtho2D(0.0,200.0,0.0,150.0);

}

class CPoint

{

public:

int x;

int y;

CPoint(){}

CPoint(int x1,int y1)

{

x=x1;

y=y1;

}

static CPoint ZeroMoveToXY(CPoint p, CPoint XY);//原始坐标向屏幕坐标XY 的平移

static CPoint ToZero(CPoint p);//关于原点对称

static CPoint XYMoveToZero(CPoint p, CPoint XY);//XY坐标向屏幕坐标的平移

static CPoint Charge_AllLine(CPoint p, CPoint line_start, CPoint line_end);//关于Ax+By+C = 0对称

};

CPoint CPoint::ZeroMoveToXY(CPoint p, CPoint XY)

{

//原始坐标向屏幕坐标XY的平移

CPoint result;

int change[3][3] = {{1,0,0},{0,1,0},{-XY.x,-XY.y,1}};

int p1[3] = {0,0,0};

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

{

p1[j]=p.x * change[0][j] + p.y * change[1][j] + change[2][j];

}

result.x = p1[0];

result.y = p1[1];

return result;

}

CPoint CPoint::ToZero(CPoint p)

{

//关于原点对称

CPoint result;

int change[3][3] = {{-1,0,0},{0,-1,0},{0,0,1}};

int p1[3] = {0,0,0};

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

{

p1[j]=p.x * change[0][j] + p.y * change[1][j] + change[2][j];

}

result.x = p1[0];

result.y = p1[1];

return result;

}

CPoint CPoint::XYMoveToZero(CPoint p, CPoint XY)

{

//XY坐标向屏幕坐标的平移

CPoint result;

int change[3][3] = {{1,0,0},{0,1,0},{XY.x,XY.y,1}};

int p1[3] = {0,0,0};

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

{

p1[j]=p.x * change[0][j] + p.y * change[1][j] + change[2][j];

}

result.x = p1[0];

result.y = p1[1];

return result;

}

CPoint CPoint::Charge_AllLine(CPoint p, CPoint line_start, CPoint line_end)

{//关于Ax+By+C = 0对称

double A,B,C;

if(line_start.x == line_end.x)

{

A = 1;

B = 0;

C = -line_start.x;

}

else if(line_start.y == line_end.y)

{

A = 0;

B = 1;

C = -line_start.y;

}

else

{

A = 1000/(line_end.x - line_start.x);

B = -1000/(line_end.y - line_start.y);

C = 1000*(-line_start.x/(line_end.x-line_start.x)+line_start.y/(line_end.y - line_start.y));

}

CPoint result;

float Y_f = (-A*B*p.x + A*A*p.y - B*C)/(B*B+A*A);

float X_f = (-A*B*p.y + B*B*p.x - A*C)/(B*B+A*A);