Skip to content

Instantly share code, notes, and snippets.

@stggh
Created April 20, 2019 01:27
Show Gist options
  • Save stggh/bceb01f1d1ba3a492728046b641dce62 to your computer and use it in GitHub Desktop.
Save stggh/bceb01f1d1ba3a492728046b641dce62 to your computer and use it in GitHub Desktop.
Testing post transform pixel shift for PyAbel onion_bordas
import abel
import matplotlib.pyplot as plt
im = abel.tools.analytical.SampleImage(1001, 'Ominus').image
proj = abel.Transform(im, method='basex', verbose=False,
direction='forward').transform
bas = abel.Transform(proj, method='basex', verbose=False,
angular_integration=True)
onb = abel.Transform(proj, method='onion_bordas',
angular_integration=True)
ims = abel.tools.vmi.angular_integration(im)
plt.plot(*ims, '--', label='Ominus')
plt.plot(*bas.angular_integration, lw=1, label='basex')
plt.plot(*onb.angular_integration, lw=1, label='onion_bordas')
plt.legend()
plt.axis(xmin=330,xmax=350)
plt.savefig('onion-y.png', dpi=75, bb_inches='tight')
plt.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment