Question

我正在做这项任务以获得乐趣。

http://groups.csail.mit.edu/graphics/classes/6.837/F04/assignments/assignment0/

如果您想查看它应该如何显示，那么现场有样本输出。它涉及迭代函数系统，其算法根据赋值：

              for "lots" of random points (x0, y0)
                   for k=0 to num_iters 
                       pick a random transform fi
                       (xk+1, yk+1) = fi(xk, yk)
                   display a dot at (xk, yk)

我的实施遇到了麻烦，即：

                void IFS::render(Image& img, int numPoints, int numIterations){

                    Vec3f color(0,1,0);

                    float x,y;
                    float u,v;
                    Vec2f myVector;
                    for(int i = 0; i < numPoints; i++){

                        x = (float)(rand()%img.Width())/img.Width();
                        y = (float)(rand()%img.Height())/img.Height();
                        myVector.Set(x,y);

                        for(int j = 0; j < numIterations;j++){

                            float randomPercent = (float)(rand()%100)/100;

                            for(int k = 0; k < num_transforms; k++){
                                if(randomPercent < range[k]){
                                    matrices[k].Transform(myVector);
                                }
                            }
                        }
                        u = myVector.x()*img.Width();
                        v = myVector.y()*img.Height();  

                        img.SetPixel(u,v,color);

                    }
                }

这就是我从输入矩阵中选择一个随机变换的方法：

                        fscanf(input,"%d",&num_transforms);

                        matrices = new Matrix[num_transforms];  
                        probablility = new float[num_transforms];
                        range = new float[num_transforms+1];

                        for (int i = 0; i < num_transforms; i++) { 
                            fscanf (input,"%f",&probablility[i]);
                            matrices[i].Read3x3(input);

                            if(i == 0) range[i] = probablility[i];
                            else range[i] = probablility[i] + range[i-1];

                        }

我的输出仅显示Sierpinski三角形的开头（1000点，1000次迭代）：

sierpinski

我的龙更好，但仍需要一些工作（1000分，1000次迭代）：

dragon

Answer 1

如果您有RAND_MAX = 4且图片宽度为3，则来自rand（）的均匀分布的序列（如[0,1,2,3,4]）将被映射到[0,1,2,0,1]你的模数代码，即一些数字会更频繁地出现。您需要切断那些高于目标范围最高倍数的数字，这些数字低于RAND_MAX，即高于（（RAND_MAX / 3）* 3）。只需检查此限制并再次调用rand（）。

由于您必须在多个位置修复该错误，请考虑编写实用程序功能。然后，减少变量的范围。 u，v声明使得很难看出这两个代码只用在三行代码中。将它们声明为“unsigned const u = ...”以使其清楚，并且另外让编译器检查您之后是否意外修改它们。

生成不完整的迭代函数系统

1 个答案: