Menger Sponge

画了一个 Menger 海绵。用povray 运行代码即可。

1. #include "colors.inc"
   
2. 
   
3. camera {
   
4.   location <150, 450, -300>
   
5.   look_at  <150, 150, 150>
   
6.   right x*image_width/image_height
   
7.   up y
   
8. }
   
9. 
   
10. light_source {
    
11.   <150, 800, -500>
    
12.   color rgb <1, 1, 1>
    
13. }
    
14. 
    
15. #declare T_Cube =
    
16.   texture {
    
17.     pigment {
    
18.       color rgb <1, 0.2, 0>
    
19.     }
    
20.   }
    
21. 
    
22. #declare T_Menger =
    
23.   texture {
    
24.     pigment {
    
25.       bozo
    
26.       color_map {
    
27.         [0 color White*0.7]
    
28.         [1 color White*1.1]
    
29.       }
    
30.     }
    
31.    
    
32.     normal {
    
33.       agate 0.3
    
34.       scale <1, 1, 1>*1.6
    
35.     }
    
36.    
    
37.     finish {
    
38.       ambient Black
    
39.       diffuse 0.8
    
40.       specular 0.1
    
41.       roughness 0.1
    
42.       brilliance 1
    
43.       metallic
    
44.     }
    
45.    
    
46.     scale <1, 1, 1>*0.4
    
47.   }
    
48. 
    
49. #declare prod1 = -0.1;
    
50. #declare prod2 = 1.1;
    
51. 
    
52. #macro RemoveCubes(level)
    
53.   #local cubelen = 1.0 / pow(3, level+1);
    
54.   #local dz = cubelen;
    
55.   
    
56.   #while (dz <= 1)
    
57.     #local dy = 1.0 / pow(3, level+1);
    
58.     #while (dy <= 1)
    
59.       box{ <prod1, dy, dz>, <prod2, dy+cubelen, dz+cubelen> texture { T_Cube }}
    
60.       box{ <dy, dz, prod1>, <dy+cubelen, dz+cubelen, prod2> texture { T_Cube }}
    
61.       box{ <dy, prod1, dz>, <dy+cubelen, prod2, dz+cubelen> texture { T_Cube }}
    
62.       
    
63.       #if (level<2)
    
64.         light_source { <(prod1+prod2)/2, dy, dz> color rgb <1, 1, 1>}
    
65.         light_source { <dy, dz, (prod1+prod2)/2> color rgb <1, 1, 1>}
    
66.         light_source { <dy, (prod1+prod2)/2, dz> color rgb <1, 1, 1>}
    
67.       #end
    
68.       
    
69.       #local dy = dy + 1.0 / pow(3, level);
    
70.     #end
    
71.     #local dz = dz + 1.0 / pow(3, level);
    
72.   #end
    
73. #end
    
74. 
    
75. 
    
76. #declare menger =
    
77.   difference {
    
78.     box { <0, 0, 0>, <1, 1, 1> }
    
79.     RemoveCubes(0)
    
80.     RemoveCubes(1)
    
81.     RemoveCubes(2)
    
82.     RemoveCubes(3)
    
83. }
    
84. 
    
85. object {
    
86.   menger
    
87.   texture{ T_Menger }
    
88.   scale 300
    
89. }

复制代码

是否为中空的呢？
Asymptote 的主页那个是我当年画的，是用最少的曲目片构成，因此里面是空的，人可以无障碍在里面行走，是一个连通的建筑，但我代码那时候写得很乱，也不知如何以中空的形式达到你的这么高阶。

herbertfederer 发表于 2018-2-1 20:37
是否为中空的呢？
Asymptote 的主页那个是我当年画的，是用最少的曲目片构成，因此里面是空的，人可以无障 ...

这个图大概不是你指的那种中空。这个就是普通的CSG，在一个单位正方体里面反复挖去小的立方体，没有挖去的部分仍然是实心的。

对光线追踪技术来说，看不到曲面是不会去计算其颜色的。