最简单的分形图像生成算法

      本文将提供一段完整地生成一幅分形图像文件的C语言代码,并且极为简单.我相信这应该是最简单的分形图像生成算法.大部分的分形图像代码也都很短,但一有递归迭代就难以理解了.而这段代码则很好懂,并且其生成的图像会使人意想不到.

#include <iostream>
#include <cmath>
#include <cstdlib>

#define DIM 1000

void pixel_write(int,int);
FILE *fp;

int main()
{
    fp = fopen("image.ppm","wb");
    if (!fp)
    {
        return -1;
    }

    fprintf(fp, "P6\n%d %d\n255\n", DIM, DIM);
    for(int j=0;j<DIM;j++)
    {
        for(int i=0;i<DIM;i++)
        {
            pixel_write(i,j);
        }
    }
    fclose(fp);

    return 0;
}

void pixel_write(int i, int j)
{
    static unsigned char color[3];
    float t = j + i*0.001f;
    memcpy(color, &t, 3);
    fwrite(color, 1, 3, fp);

    // 其实更简单粗爆的方式是
     //fwrite(&t, 1, 3, fp);
}

代码运行完后,会生成一个PPM格式的图像文件,为:

最简单的分形图形应该是康托(Cantor)三分集,只是不停地将线段一分为二.如下图所示:

该算法生成的图像很康托三分集有些相似.尤其是它的绿色通道:

将代码修改下,其分形的效果会更明显:

void pixel_write(int i, int j)
{
    static unsigned char color[3];
    float t = j/4 + i*0.001f;
    memcpy(color, &t, 3);
    fwrite(color, 1, 3, fp);
}

相关文章:

C语言中将0到1000的浮点数用强制指针类型转换的方式生成一幅图像