[OpenGL] 绘制并且判断凹凸多边形、自相交多边形。

#include <iostream>
#include <ctime>
#include <GL/glut.h>
#include <math.h>
#include <vector>

using namespace std;

struct Pos
{
    int x;
    int y;
};

struct Edge
{
    int x1, x2;
    int y1, y2;
    int vx;
    int vy;
    int a, b, c;
};

struct Poly
{
    // 点集
    int xx[100];
    int yy[100];

    // 边集
    Edge Edges[100];

    int plotNums = 0; //点数量
    int edgeNums = 0; //边数量

    // 记录凹点
    int conv = 0;
};

Poly poly;


// 求交点坐标
Pos CrossPos(int p1,int p2){
    Pos res;
    int A1 = poly.Edges[p1].a;
    int B1 = poly.Edges[p1].b;
    int A2 = poly.Edges[p2].a;
    int B2 = poly.Edges[p2].b;
    int C1 = poly.Edges[p1].c;
    int C2 = poly.Edges[p2].c;

    int m = A1 * B2 - A2 * B1;
    if (m == 0)
        cout <<"第"<<p1<<"边和第"<<p2<<"边"<< "无交点" << endl;
    else
    {
        res.x = (C2*B1 - C1 * B2) / m;
        res.y = (C1*A2 - C2 * A1) / m;
    }
    return res;

}

// 判断点是否在线段内
bool JudgeInLine(Pos pos,Edge edge)
{
    int maxX = edge.x1 >= edge.x2 ? edge.x1 : edge.x2;
    int minX = edge.x1 <= edge.x2 ? edge.x1 : edge.x2;
    int maxY = edge.y1 >= edge.y2 ? edge.y1 : edge.y2;
    int minY = edge.y1 <= edge.y2 ? edge.y1 : edge.y2;
    if (pos.x<=maxX && pos.x>=minX && pos.y<=maxY && pos.y>=minY)
    {
        return true;
    }
    return false;
}

// 求叉积 返回正负值
int CrossProduct(int n)
{
    n = n % poly.edgeNums;
    int np = (n + 1) % poly.edgeNums;
    return (poly.Edges[n].vx*poly.Edges[np].vy - poly.Edges[np].vx*poly.Edges[n].vy) >= 0 ? 1 : -1;
}

// 切割凹多边形
void ChangePoly() {
    int convP = poly.conv; //凹点的下一个点
    Pos interPos;


    for (int i = 0; i < poly.edgeNums; i++)
    {
        if (i<convP-1 || i>convP+1)
        {
            interPos = CrossPos(convP, i);
        }
    }

    glClearColor(0.0, 0.0, 0.0, 0.0);
    glClear(GL_COLOR_BUFFER_BIT);
    glBegin(GL_POLYGON);
    glColor3f(1.0f, 0.0f, 0.0f);
    for (int i = 0; i <= convP; i++)
    {
        glVertex2f(poly.xx[i], poly.yy[i]);
    }
    glVertex2f(interPos.x, interPos.y);
    glEnd();

    glBegin(GL_POLYGON);
    glColor3f(0.0f, 1.0f, 1.0f);
    glVertex2f(interPos.x, interPos.y);
    for (int i=convP+1;i<poly.plotNums;i++)
    {
        glVertex2f(poly.xx[i], poly.yy[i]);
    }
    glEnd();

    glFlush();
}

// 判断是什么多边形
bool Judge()
{
    /*输出边信息*/
    for (int i = 0; i < poly.edgeNums; i++)
    {
        cout << "Vx:" << poly.Edges[i].vx << "  " << "Vy:" << poly.Edges[i].vy << "  " << "A:" << poly.Edges[i].a<< "  " << "B:" << poly.Edges[i].b << "  " << "C:" << poly.Edges[i].c <<endl;
    }

    /*判断自交*/
    Pos interPos;
    if (poly.edgeNums > 3)
        for (int i = 0; i < poly.edgeNums; i++)
        {
            interPos = CrossPos(i, (i + 2) % poly.edgeNums);
            if (JudgeInLine(interPos, poly.Edges[i]) && JudgeInLine(interPos, poly.Edges[(i + 2) % poly.edgeNums]))
            {
                cout << "该多边形为自相交多边形" << endl;
                return false;
            }
        }

    /*判断凹凸*/
    // 判断向量叉积 是否为同一正负
    int judge;
    if (CrossProduct(0) >= 0)
        judge = 1;
    else
        judge = -1;
    //判断每一个角,两边向量乘积是否同符号
    for (int i = 1; i <= poly.edgeNums; i++)
    {
        if (judge*CrossProduct(i) < 0)
        {
            poly.conv = i;
            ChangePoly();
            cout << "该多边形为凹多边形" << endl;
            return false;
        }
    }
    cout << "该多边形为凸多边形" << endl;
    return true;
}

void init(void)
{
    glClearColor(0.0, 0.0, 0.0, 0.0);

    glMatrixMode(GL_PROJECTION);//设置投影矩阵
    gluOrtho2D(0.0, 400.0, 0.0, 300.0);//二维视景区域

    glColor3f(1.0, 0.0, 0.0);
    glPointSize(3.0);//点的大小
}

void plotpoint(GLint x, GLint y)
{
    glBegin(GL_POINTS);
    glVertex2i(x, y);
    glEnd();
}

void displayFcn(void)
{
    glClearColor(0.0, 0.0, 0.0, 0.0);
    glClear(GL_COLOR_BUFFER_BIT);
    for (int i = 0; i < poly.plotNums; i++)
    {
        plotpoint(poly.xx[i], poly.yy[i]);
    }
    glBegin(GL_POLYGON);
    for (int i = 0; i < poly.edgeNums; i++)
    {
        glVertex2f(poly.xx[i], poly.yy[i]);
    }
    glEnd();

    glFlush();
}

void mouse(GLint button, GLint action, GLint x, GLint y)
{
    if (button == GLUT_LEFT_BUTTON && action == GLUT_DOWN)
    {
        poly.xx[poly.plotNums] = x;
        poly.yy[poly.plotNums] = 300 - y;
        cout << "x:" << x << "  " << "y:" << 300 - y << endl;
        poly.plotNums++;
        glutPostRedisplay();//重绘窗口
    }
    if (button == GLUT_RIGHT_BUTTON && action == GLUT_DOWN)
    {
        poly.edgeNums = poly.plotNums;
        if (poly.plotNums > 2)
        {
            for (int i = 1; i <= poly.plotNums; i++)
            {
                poly.Edges[i - 1].x1 = poly.xx[i - 1];
                poly.Edges[i - 1].y1 = poly.yy[i - 1];
                poly.Edges[i - 1].x2 = poly.xx[i%poly.plotNums];
                poly.Edges[i - 1].y2 = poly.yy[i%poly.plotNums];
                poly.Edges[i - 1].vx = poly.Edges[i - 1].x2 - poly.Edges[i - 1].x1;
                poly.Edges[i - 1].vy = poly.Edges[i - 1].y2 - poly.Edges[i - 1].y1;
                poly.Edges[i - 1].a = poly.Edges[i - 1].vy;
                poly.Edges[i - 1].b = -poly.Edges[i - 1].vx;
                poly.Edges[i - 1].c = poly.Edges[i - 1].x2 * poly.Edges[i - 1].y1 - poly.Edges[i - 1].x1 * poly.Edges[i - 1].y2;
            }
            if (Judge())
                glutPostRedisplay();//重绘窗口
        }
    }
}

int 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();

}

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 200,176评论 5 469
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,190评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 147,232评论 0 332
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,953评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,879评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,177评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,626评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,295评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,436评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,365评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,414评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,096评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,685评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,771评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,987评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,438评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,032评论 2 341

推荐阅读更多精彩内容