Question

编辑：我不确定这是不是一个新问题所以我现在只是更新这个问题。除非我更改原始坐标，否则现在正在生成雪花。例如，如果我的原始三角形像图片1，则5次迭代后的结果是图片2：

enter image description here

但是，如果我的原始三角形有任何不同，例如图片3，则结果会偏斜：

enter image description here

我再次认为问题出在我的身上，但我真的迷路了。我已经尝试了几个小时来找出正确的公式来做到这一点，但我并没有真正取得任何进展。由于似乎很多时候问题是新三角形被倒置了，我怀疑atan在给它一个正值时给了我一个负值。有没有数学方法可以解决这个问题？非常感谢你的帮助！

我的代码（减去openGL部分，因为我不认为它们是问题）是：

const int NumTimesToSubdivide = 5;
const int NumSegments = 3072;  // 3x4^n lines generated {3, 12, 48, 192, 768, 3072..}
const int NumVertices = NumSegments * 2; // 2 vertices for each segment

vec2 arrayA[NumVertices];
vec2 arrayB[NumVertices];

//----------------------------------------------------------------------------
void koch( const vec2 &a, const vec2 &b, vec2 * verts) {        

/* Calculate new segments.
 *         v1      
 *         /\
 *        /  \
 *       /    \
 * a----/      \----b
 *     v0      v2
 */
vec2 v0;
vec2 v1;
vec2 v2;

GLfloat distance(sqrt(pow(b.x - a.x, 2) + pow(b.y - a.y,2)));
GLfloat deltaX = abs(b.x - a.x);
GLfloat deltaY = abs(b.y - a.y);
GLfloat normalX((b.x - a.x) / deltaX);
GLfloat normalY((b.y - a.y) / deltaY);
GLfloat theta = atan2(b.y - a.y, b.x - a.x) - M_PI / 3.0;
cout << " theta = " << theta << endl;

/*************************
 * Find trisection points
 *************************/
// horizontal line _____________
if(a.y == b.y) {
    vec2 temp0(a.x + (deltaX * normalX / 3) , a.y);
    vec2 temp2(a.x + (deltaX * normalX * 2/3) , a.y);
    vec2 temp1((a.x + b.x) / 2, a.y + distance * sin(theta) / 3);
    v0 = temp0;
    v2 = temp2;
    v1 = temp1;
}

//                 |
// vertical line   |
//                 |
else if(a.x == b.x){
    vec2 temp0(a.x , (a.y + (deltaY * normalY / 3)));
    vec2 temp2(a.x , (a.y + (deltaY * normalY * 2/3)));
    vec2 temp1(a.x + distance * cos(theta) / 3 , (a.y + b.y) / 2);
    v0 = temp0;
    v2 = temp2;
    v1 = temp1;
}

// slope != 0 && slope != 1
else {
    vec2 temp0(a.x + (deltaX * normalX / 3), a.y + (deltaY * normalY / 3));
    vec2 temp2(a.x + (deltaX * normalX * 2/3), a.y + (deltaY * normalY * 2/3));
    // Andrew is the greatest!
    vec2 temp1(temp0.x + distance * cos(theta) / 3, 
               temp0.y + distance * sin(theta) / 3);        
    v0 = temp0;
    v2 = temp2;
    v1 = temp1;
}

verts[0] = a;
verts[1] = v0;
verts[2] = v0;
verts[3] = v1;
verts[4] = v1;
verts[5] = v2;
verts[6] = v2;
verts[7] = b;

}

//----------------------------------------------------------------------------

void divide_line( const vec2& a, const vec2& b, const vec2& c, int n )
{ 

// arrayA = {a, b, b, c, c, a} i.e., the sides of the initial triangle
arrayA[0] = a;
arrayA[1] = b;
arrayA[2] = b;
arrayA[3] = c;
arrayA[4] = c;
arrayA[5] = a;

// If at least one iteration:
if(n > 0) {
    // The current iteration, starting with 0
    int currentIteration = 1;
    // Do for every iteration from 0 - n:
    while (currentIteration <= n) {
        int i;
        int j = 0;
        int size = 3 * 2 * (pow(4,currentIteration-1));
        // Call koch() for each pair of vertices in arrayA
        for(i = 0; i < size; i = i+2) {
            vec2 verts[8];
            koch(arrayA[i], arrayA[i+1], verts);
            // Store each vertex in arrayB
            int k;
            for(k = 0; k <= 7; k++)
                arrayB[j++] = verts[k];
        }
        // Copy arrayB to arrayA for next iteration.
        size = 3 * 2 * pow(4, currentIteration);
        for(i = 0; i < NumVertices; i++) {
            arrayA[i] = arrayB[i];
        }
        // Increase count of currentIteration.
        currentIteration++;
    }
} else
    printf("The number of iterations must be >= 0.\n");
}

我目前正在尝试用c ++实现Koch曲线。我几乎正常工作。我正在解决的基本上有三种不同的情况：我需要分成四个段的线段是水平，垂直或其他。

问题在于，当程序计算线段的新顶点时，有两种可能性：三角形可以面向“向上”或三角形可以面向“向下”。我试图通过规范化向量来考虑这一点，但我要么做错了，要么稍微有点别的。下面的图片是两个例子;两者都有3次迭代的分形。

分割非水平和非垂直线段的代码如下所示，因为三角形看起来像是垂直/水平线段的正确方式：

if(a.x == b.x) {
  ...
}
else if (a.y == b.y) {
  ...
}
// slope != 0 && slope != 1
else {
    GLfloat deltaX = abs(b.x - a.x);
    GLfloat deltaY = abs(b.y - a.y);
    vec2 temp0(a.x + (deltaX * normalX / 3), a.y + (deltaY * normalY / 3));
    vec2 temp2(a.x + (deltaX * normalX * 2/3), a.y + (deltaY * normalY * 2/3));
    GLfloat dist(sqrt(pow(temp2.x - temp0.x, 2) + pow(temp2.y - temp0.y,2)));
    GLfloat theta = (a.x - b.x)/ (b.y - a.y);
    vec2 temp1((a.x + b.x) / 2 + dist * cos(atan(theta)) ,
               (a.y + b.y) / 2 + dist * sin(atan(theta)));
    v0 = temp0;
    v2 = temp2;
    v1 = temp1;
}

a和b是该段的向量。 normalX和normalY是：

GLfloat normalX((b.x - a.x) / (abs(b.x - a.x)));
GLfloat normalY((b.y - a.y) / (abs(b.y - a.y)));

我可以采取哪些措施来解决这个问题？

Answer 1

您的问题可能是atan(theta)。它的域太小，无法用于通用角度确定。您应该使用atan2(y,x)：

GLfloat theta = atan2(b.y - a.y, a.x - b.x);
vec2 temp1((a.x + b.x) / 2 + dist * cos(theta) ,
           (a.y + b.y) / 2 + dist * sin(theta));

（我没有测试过这段代码，因此可能会出现一些符号错误。）

了解详情：Atan2 (Wikipedia)

编辑：

我的猜测是你以错误的顺序给出坐标。如果从面向终点的起点的角度考虑一个段，则曲线将始终向右挤出。如果以相反的顺序给出坐标，则曲线将被镜像。

Answer 2

以逆时针顺序提供所有初始顶点，它应该按照您期望的方式工作。

不能旋转科赫雪花

2 个答案: