曲折的递归分形山

Troubling Recursive Fractal Mountain

本文关键字:分形山 递归      更新时间:2023-10-16

我正在写一些OpenGl代码来创建一个分形山。我正在使用递归,并且遇到了困难。

我遇到的问题是三角形的撕裂和重叠:https://i.stack.imgur.com/ha7gx.jpg

问题是,我以为我已经解决了这个问题,这是我的代码。

我的三角类,vec2只是一个通用的二维向量

struct Triangle {
vec2 A, B, C;
Triangle(vec2 a, vec2 b, vec2 c) {
    A = a;
    B = b;
    C = c;
};

这是我的递归函数,它最初是用默认查看体积中的顶点(-1,-1)、(0,1)和(1,-1)的三角形调用的。

 void divide_triangle(const Triangle& baseTriangle, int count) {
if (count > 0) {
    std::vector<Triangle> subTriangles = subdivideTriangle(baseTriangle);
    for (int i = 0; i < subTriangles.size(); i++) {
        divide_triangle(subTriangles[i], count - 1);
    }
}
else {
    triangle(baseTriangle); // Draw triangle at the end of recursion
}

这是我的subDivideTriangles函数,它简单地将经过的三角形划分为四个独立的三角形

std::vector<Triangle> subdivideTriangle(const Triangle& triangle){
vec4 aPrime = getMidpoint(triangle.B, triangle.C);
vec4 bPrime = getMidpoint(triangle.C, triangle.A);
vec4 cPrime = getMidpoint(triangle.A, triangle.B);
return{
    { Triangle(triangle.A, cPrime, bPrime) }, // Top triangle
    { Triangle(cPrime, triangle.B, aPrime) }, // Bottom left triangle
    { Triangle(bPrime, aPrime, triangle.C) }, // Bottom right triangle
    { Triangle(cPrime, aPrime, bPrime) }      // Center triangle
};

这是一个我认为可以停止撕裂的函数,我跟踪一对向量对应的中点值,如果我找到对应的向量对,我会返回原始中点,因此不会移动中点两次。

然而,这种方法似乎并不像我想的那样奏效。

vec4 getMidpoint(const vec4& a, const vec4& b) {
auto AB = std::make_pair(a, b), BA = std::make_pair(b, a);
auto foundResult1 = memo.find(AB);
auto foundResult2 = memo.find(BA);
if (foundResult1 != memo.end())
    return memo[AB];
if (foundResult2 != memo.end())
    return memo[BA];
auto result = (a + b) / 2; 
result += randVector() * getLength(a, b) * randomness; 
memo.insert(std::make_pair(AB, result)); 
memo.insert(std::make_pair(BA, result));
return result;

你们有什么见解吗,我的中点函数有问题吗?

您的问题来自相邻的边。

假设你有两个三角形ABC和DBC,BC是公共边,你想通过添加一些有噪声的中点来修改它们。

根据三角形的不同,可以使用两个不同的点作为BC的中间。

  • 当你在ABC上修改它时,让E是BC的中间
  • 当您对DBC进行修改时,让F是BC的中间

E和F必须相同。否则就会出现很大的差距。它在你的照片上很容易看到。

在我看来,您应该使用索引顶点,并使用一个表(或任何方法)来查找给定边的现有中间。即:当在DBC中计算BC的F中间时,我们检查并发现它已经作为E.存在