本文将提供一段完整地生成一幅分形图像文件的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
%d %d
255
", 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);
}
相关文章: