Most general operations like this can be performed with e2proc2d.py (for 2-D images), e2proc3d.py (for 3-D volumes) or e2filtertool.py (for dynamic tuning of parameters).
e2proc3d.py --help
will give a list of all of the program options, but many things you can do is via "processors"
e2help.py processors
What you are trying to do would be
e2proc3d.py <input map> <out text file> --calcradial <0 or 1>
if you give it 0, it will compute radial amplitudes, which is probably what you want (1 will compute radial "intensity", which mainly makes sense in Fourier space)
You can use any plotting program you like, of course, but:
e2display.py --plot <text file>
then middle-click for the control panel. This will allow you to select which column to use for the x and y axes easily. The text file will contain
radius
mean min
max sigma
for each radial shell