关闭

[算法]柏林噪声 Perlin Noise

标签: C++图形学算法游戏开发OpenGL
6043人阅读 评论(0) 收藏 举报
分类:

时间:2017/1/26 大三上学年寒假

       每次寒假都挺闲的,恰好最近家里网又不太稳定,不能LOL!于是想趁这个机会,利用OpenGL,尝试着写一个仿造的《Minecraft》。但是,思来想去发现Minecraft的地形生成算法实在难以下手,因为地形需要具有:随机但又连续、复杂但又自然的效果。

       凑巧在github上偶尔看到一个项目,一个大神写的minecraft,做的非常非常好。然后在他的说明文档里,注意到这个算法PerlinNoise。于是,决定写几篇博客,来总结一下关于PerlinNoise的知识。

       这里也附上我使用Qt+OpenGL实现的minecraft的结果。

相关文章:

  [1]博客 - [数学][转载]柏林噪声 :http://www.cnblogs.com/Memo/archive/2008/09/08/1286963.html

  [2]gitbub - A simple Minecraft clone written in C using modern OpenGL (shaders). https://github.com/fogleman/Craft

  [3] 我的另一篇博客《我的世界minecraft》:http://blog.csdn.net/mahabharata_/article/details/55505403


我实现的《myMinecraft》的截图:











我实现的Minecraft演示程序的下载链接:

         http://download.csdn.net/detail/mahabharata_/9756981

下面这个PerlinNoise演示程序的下载链接,大家可以自己调整参数来设计自己的柏林噪声~~~

         http://download.csdn.net/detail/mahabharata_/9800578


柏林噪声Perlin Noise:

       关于柏林噪声的原理和算法在参考文献[1]的博客中,已经描述的非常详细。我在阅读Perlin的论文和其它博客后,先写了一个一维PerlinNoise的演示程序~。在写完这个演示程序之后,发现原来不懂的东西,现在变得直观多了!!下面写一下在演示程序中,总结的一些柏林噪声的规律。

首先上一下PerlinNoise演示程序的截图~(很简陋,一晚上写的),大家可以自己下载下来这个程序去简单模拟一下~

         

1. 持续度Persistence [ float ]:我设定的范围是(0,1]。

        持续度越大,频率frequency越大,振幅amplitude也越大(反应在图像上,直观的表现为:波动更剧烈,更尖锐)。

 

2. 倍频Octaves [ int ]:我设定的范围是[0,99]

            倍频越大,频率和振幅也有所增加(图像波动变大)


3.  噪声点间距Delta [ float ],在文章中没说到的一个值,简单的说就相当于我们数学求导中的dx,在这里相邻噪声数据的间隔。我第一次设定为1,显然生成的噪声波动很剧烈。之后较多采用的一个值是0.001,产生的数据就规律多了~上图说话!

      


上面Perlin Noise演示程序采用的柏林噪声的2D版([x,y]->z):

float persistence = 0.50;
int Number_Of_Octaves = 4;

double Noise(int x,int y)    // 根据(x,y)获取一个初步噪声值
{
    int n = x + y * 57;  
    n = (n<<13) ^ n;
    return ( 1.0 - ( (n * (n * n * 15731 + 789221) + 1376312589) & 0x7fffffff) / 1073741824.0);
}

double SmoothedNoise(int x, int y)   //光滑噪声
{
    double corners = ( Noise(x-1, y-1)+Noise(x+1, y-1)+Noise(x-1, y+1)+Noise(x+1, y+1) ) / 16;
    double sides = ( Noise(x-1, y) +Noise(x+1, y) +Noise(x, y-1) +Noise(x, y+1) ) / 8;
    double center = Noise(x, y) / 4;
    return corners + sides + center;
}
double Cosine_Interpolate(double a,double b, double x)  // 余弦插值
{
    double ft = x * 3.1415927;
    double f = (1 - cos(ft)) * 0.5;
    return a*(1-f) + b*f;
}

double InterpolatedNoise(float x,float y)   // 获取插值噪声
{
    int integer_X = int(x);
    float  fractional_X = x - integer_X;
    int integer_Y = int(y);
    float fractional_Y = y - integer_Y;
    double v1 = SmoothedNoise(integer_X, integer_Y);
    double v2 = SmoothedNoise(integer_X + 1, integer_Y);
    double v3 = SmoothedNoise(integer_X, integer_Y + 1);
    double v4 = SmoothedNoise(integer_X + 1, integer_Y + 1);
    double i1 = Cosine_Interpolate(v1, v2, fractional_X);
    double i2 = Cosine_Interpolate(v3, v4, fractional_X);
    return Cosine_Interpolate(i1, i2, fractional_Y);
}

double PerlinNoise(float x,float y)    // 最终调用:根据(x,y)获得其对应的PerlinNoise值
{
    double total = 0;
    double p = persistence;
    int n = Number_Of_Octaves;
    for(int i=0; i<n; i++)
    {
        double frequency = pow(2,i);
        double amplitude = pow(p,i);
        total = total + InterpolatedNoise(x * frequency, y * frequency) * amplitude;
    }

    return total;
}


3
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:69132次
    • 积分:1310
    • 等级:
    • 排名:千里之外
    • 原创:47篇
    • 转载:0篇
    • 译文:3篇
    • 评论:115条
    最新评论