More concise version 

First we compute the r and theta of every pixel, then get all the intensities that are inside that range, then get the radi transformed the same way the intensities were transformed, bin them together using scipy.stats.binned_statistic and put out the resultant x and y arrays.