如何反转OpenGL代码输出？科赫曲线反转修改及顶点公式咨询

看起来你在OpenGL绘制科赫曲线时遇到了输出反转的问题，同时想了解更多计算等边三角形顶点的方法——我来帮你一步步解决这些问题：

一、搞定输出反转的两种方法 #

OpenGL默认的窗口坐标系是左下角为原点(0,0)，右上角对应窗口的宽高；而很多时候我们会习惯左上角作为原点的坐标系（比如常见的UI框架），这就是输出反转的核心原因。这里有两种靠谱的修复方式：

1. 调整投影矩阵（推荐） #

在你的初始化函数（比如myinit）或者显示函数里，修改正交投影的参数，直接翻转y轴的范围：

void myinit() {
    glClearColor(0.0,1.0,0.0,1.0);
    glColor3f(0.0,0.0,1.0);
    glPointSize(3.0);
    
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    // 把y轴范围翻转，让左上角成为原点
    int winWidth = glutGet(GLUT_WINDOW_WIDTH);
    int winHeight = glutGet(GLUT_WINDOW_HEIGHT);
    glOrtho(0, winWidth, winHeight, 0, -1, 1);
    glMatrixMode(GL_MODELVIEW);
}

这样整个坐标系就会和你习惯的一致，科赫曲线就能正常显示了。

2. 手动翻转顶点y坐标 #

如果你不想改动投影矩阵，也可以在计算每个顶点的时候，用窗口高度减去原来的y值：

// 在绘制前获取窗口高度
int winHeight = glutGet(GLUT_WINDOW_HEIGHT);
// 翻转y坐标
int flippedY1 = winHeight - y1;
int flippedY2 = winHeight - y2;
// 之后用flippedY1、flippedY2来绘制线段

二、科赫曲线等边三角形顶点的其他计算方法 #

你当前用预定义的COS和SIN常量是一种简化方式，还有几种更灵活、可读性更强的方法：

1. 向量旋转法（最通用） #

科赫曲线的核心是把线段分成三段，在中间段向外（上/下）生成等边三角形。用向量旋转可以清晰计算顶点：
假设线段两端点是A(x1,y1)和B(x2,y2)：

• 计算向量AB = (dx, dy) = (x2-x1, y2-y1)
• 线段1/3长度：segmentLen = sqrt(dx*dx + dy*dy)/3
• 中间分点C：(x1 + dx/3, y1 + dy/3)，D：(x1 + 2*dx/3, y1 + 2*dy/3)
• 把向量CD（也就是AB/3）逆时针旋转60度，得到三角形的顶点E：

const float angle60 = M_PI / 3; // 60度转弧度
float cos60 = cos(angle60);
float sin60 = sin(angle60);

// 旋转向量AB/3
float rotatedX = (dx/3)*cos60 - (dy/3)*sin60;
float rotatedY = (dx/3)*sin60 + (dy/3)*cos60;

// E点坐标
float ex = x1 + dx/3 + rotatedX;
float ey = y1 + dy/3 + rotatedY;

如果想让三角形向下生成，就把旋转方向改成顺时针，只需要把sin60的符号反过来：

float rotatedX = (dx/3)*cos60 + (dy/3)*sin60;
float rotatedY = -(dx/3)*sin60 + (dy/3)*cos60;

2. 几何公式法 #

利用等边三角形的高的性质直接计算：
等边三角形的高是(边长 * √3)/2，这里的边长是原线段的1/3，所以高是(原线段长度 * √3)/6。结合线段中点，垂直于线段方向移动这个高度就能得到顶点：

float dx = x2 - x1;
float dy = y2 - y1;
// 计算垂直方向的偏移量
float offsetX = -dy * (sqrt(3)/6);
float offsetY = dx * (sqrt(3)/6);

// E点坐标（基于线段AB的1/3处点C）
float ex = x1 + dx/3 + offsetX;
float ey = y1 + dy/3 + offsetY;

这个方法不用计算旋转，直接利用几何关系，代码更简洁。

3. 动态计算三角函数值 #

如果你不想用预定义的0.5和0.866，可以直接用标准库函数计算60度的三角函数值，可读性更好：

// 60度转弧度
const float rad60 = M_PI / 3;
#define COS cos(rad60)
#define SIN sin(rad60)

这样如果以后要调整分形的角度（比如改成45度），只需要修改rad60的值即可，非常灵活。

内容的提问来源于stack exchange，提问作者Niraj Wagh