Sierpinski 3D

[1]:
import k3d
import numpy as np

iteration = 4
size = 3**iteration

voxels = np.ones((size, size, size));

def iterate(length, x, y, z):

    nl = length // 3

    if nl < 1:
        return

    margin = (nl-1) // 2

    voxels[z-margin:z+margin+1, y-margin:y+margin+1, :] = 0
    voxels[z-margin:z+margin+1, :, x-margin:x+margin+1] = 0
    voxels[:, y-margin:y+margin+1, x-margin:x+margin+1] = 0

    for ix,iy,iz in np.ndindex((3,3,3)):
        if (1 if ix !=1 else 0) + (1 if iy != 1 else 0) + (1 if iz != 1 else 0) !=2:
            iterate(nl, x + (ix-1) * nl, y + (iy-1) * nl , z + (iz-1) * nl)

iterate(size, size//2, size//2, size//2)

color_map = (0xffff00)
plot = k3d.plot()
plot += k3d.voxels(voxels.astype(np.uint8), color_map)
plot.lighting = 2.0
plot.display()