思路:通过鼠标响应函数输入n个顶点(n>=3)显示多边形,鼠标左击加一个顶点,鼠标右击显示顺序相连的多边形。需要建立多边形数据结构(保存顶点表、边信息,及其他信息)
判断方法:为多边形每一个边建立一个向量,凸多边形所有相邻的2边的向量的差积的z分量均为同号。如果出现某些为正,某些为负,则为凹多边形。
/* 数学知识补充:
将向量用坐标表示(三维向量),
若向量a=(a1,b1,c1),向量b=(a2,b2,c2),
则
向量a·向量b=a1a2+b1b2+c1c2 ------------向量点乘
向量a×向量b= ------------向量叉乘
| i j k|
|a1 b1 c1|
|a2 b2 c2|
=(b1c2-b2c1,c1a2-a1c2,a1b2-a2b1)
(i、j、k分别为空间中相互垂直的三条坐标轴的单位向量)。
*/
#include <iostream>
#include <vector>
#include <iterator>
#include "glut.h"
using namespace std;
struct vertex
{
vertex(GLint _x,GLint _y):x(_x),y(_y)
{
}
GLint x;
GLint y;
};
class side
{
public:
side(vertex _begin,vertex _end):begin(_begin),end(_end)
{
x=end.x-begin.x;
y=end.y-begin.y;
z=0;
}
public:
vertex begin;
vertex end;
GLint x;
GLint y;
GLint z;
};
// 一个全局函数,求2个边的差积
GLint getchajiZ(const side& left,const side& right)
{
// (a1,b1,c1) (a2,b2,c2)
// z=a1b2-a2b1
return left.x * right.y - right.x * left.y;
}
vector<vertex> vecVertex; // 因为顶点的数量不固定,由用户左击次数决定,所以采用vector向量容器
vector<side> vecSide;
vector<int> vecZ;
void init(void)
{
glClearColor(0.0,0.0, 1.0, 1.0);
glMatrixMode(GL_PROJECTION);//设置投影矩阵
gluOrtho2D(0.0, 400.0, 0.0, 300.0);//二维视景区域
glColor3f(1.0,0.0,0.0);
glPointSize(13.0);//点的大小
}
void displayFcn(void)
{
glClear(GL_COLOR_BUFFER_BIT);
glBegin(GL_POLYGON); // 绘制模式
glColor3f(200.0f,270.0f,0.0f);
vector<vertex>::iterator it=vecVertex.begin();
int i=0;
for(;it!=vecVertex.end();it++,i++)
{
if (it!=vecVertex.begin()) // 先取得每个边的向量保存到容器vecSide
{
side aSide(vecVertex[i-1],vecVertex[i]);
vecSide.push_back( aSide );
}
glVertex3i((*it).x,(*it).y,0); // 画出一个顶点
if (i == vecVertex.size() -1 )
{
side aSide(vecVertex[i],vecVertex[0]); // 还有一个首尾相连的边
vecSide.push_back( aSide );
}
}
// 现在想办法得到chajiZ ---差积的z
int size=vecSide.size();
for (i=0;i<size;i++)
{
int chajiZ;
chajiZ=getchajiZ(vecSide[i],vecSide[(i+1)%size]);
vecZ.push_back(chajiZ);
cout<<chajiZ<<endl; // 打印出每个差积的z值
}
for (i=0;i<vecZ.size();i++)
{
if (vecZ[i]*vecZ[ (i+1) % vecZ.size() ] <0 ) // 有异号
{
cout<<"凹多边形!!!!!!!"<<endl;
break;
}
if(i==vecZ.size()-1)
{
cout<<"凸多边形!!!"<<endl;
break;
}
}
glEnd();
glFlush();
}
void plotpoint(GLint x, GLint y)
{
// 画出临时的一个点
glBegin(GL_POINTS);
glVertex2i(x,y);
glEnd();
vecVertex.push_back(vertex(x,y)); // 不用new,全局vector容器保存局部对象可以的
}
void mouse(GLint button, GLint action, GLint x,GLint y)
{
if (button==GLUT_LEFT_BUTTON && action==GLUT_DOWN)
{
plotpoint(x,300-y);
}
glFlush();
if (button==GLUT_RIGHT_BUTTON && action==GLUT_DOWN)
{
glutPostRedisplay();//重绘窗口
// vecVertex.clear(); // 发现一旦清除,就不能显示了
// 这边顺便判断下是否为凸多边形(所有向量差积的z都同号)
}
}
void main(int argc, char** argv)
{
glutInit(&argc, argv);
glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB);
glutInitWindowPosition(50, 100);
glutInitWindowSize(400, 300);
glutCreateWindow("mouse");
init();
glutDisplayFunc(displayFcn);
glutMouseFunc(mouse);
glutMainLoop();
}
运行:
凸多边形
一个凹点对应一个和别的周围异号的叉积的z值