def roundtrip_test():
imshape = (100,100)
rmap = r.make_radmap(imshape, (25.0, 51.2), np.arange(0, 100), np.linspace(0, 2*np.pi, 17))
marr1 = rmap['map']
inds = r.maparr2indices(marr1)
marr2 = r.indices2maparr(imshape, inds)
for i in xrange(np.size(marr1)):
assert(marr1.flat[i] == marr2.flat[i])
def plot_radind(radind):
"""Display a dict containing radial indices.
"""
marr = r.indices2maparr(radind['imshape'], radind['indices'])
plt.imshow(marr2rgba(marr))
