cpos = Vector(0, 0, 6)
u = Vector.normalize(0, 0, -1)
c = Camera(cpos, u, scale=0.25)
image = []
l, h = 2**20, 0
for r in c.rays():
if r is not None:
if s.intersects(r):
col = s.intersectionDistance(r)
if col > h:
h = col
if col < l:
l = col
image += [col]
else:
image += [-1]
else:
image += [None]
row = []
imagePixels = []
for p in image:
if p is not None:
if p > 0:
col = int((2**16) * (p-l)/float(h))*5
row += [[col,col,col]]
else: row += [[2**16, 0, 0]]
else:
imagePixels += [row]
row = []
pixmap.save(imagePixels, 'out/sphere.ppm')
def p(x, y, r=10):
for j in range(w):
for i in range(h):
grid[i][j] = grid[i][j] + (lambda x: (r-x) if x <= r else 0)(d(x, y, i, j))
ps = [(0, 0), (20, 20), (500, 500), (750, 750)]
for i in range(len(ps)):
p(*ps[i], r=d(*ps[i] + ps[i-1])) if i > 1 else p(*ps[i])
def n(x, y):
l = (x, y)
cs = (
(x+1, y),
(x+1, y+1),
(x+1, y-1),
(x, y+1),
(x, y-1),
(x-1, y),
(x-1, y+1),
(x-1, y-1)
)
for c in cs:
if grid[c[0]][c[1]] > grid[l[0]][l[1]]:
l = c
print 'c > l', c, l
grid[l[0]][l[1]] = True
grid = [(lambda x:[[int(i*100)] * 3 for i in x])(i) if i is not True else
[int(i*100), 0, 0] for i in grid]
pixmap.save(grid, 'visualization.ppm')
