ray-tracing.py

for v in range(height):
  for u in range(width):
    # Get end of ray and round to sdf array coordinate...
    pos = start[v,u,:] + travel[v,u] * norm[v,u,:]
    pos += 0.5 * numpy.array(sdf.shape) # Center of sdf to 3D origin
    sdf_coordinate = (pos+0.5).astype(int)

    # Loop the 3 dimensions and check each is not out of bounds...
    for d in range(3): 
      if sdf_coordinate[d]<0:
        sdf_coordinate[d] = 0
      if sdf_coordinate[d]>=sdf.shape[d]:
        sdf_coordinate[d] = sdf.shape[d] - 1

    # Evaluate function...
    ret[v,u] = sdf[sdf_coordinate]