def radialPlot(self, imagedata, pow):
scope = imagedata['scope']
mag = scope['magnification']
tem = scope['tem']
camera = imagedata['camera']
ccd = camera['ccdcamera']
pixelsize = self.calclient.getPixelSize(mag, tem, ccd)
if camera['binning']['x'] != camera['binning']['y']:
self.logger.warning(
'Unequal binning radial average not implemented')
return
shape = pow.shape
if shape[0] != shape[1]:
self.logger.warning(
'Non square image shape radial average not implemented')
return
binned_pixelsize = pixelsize * camera['binning']['x']
rec_pixelsize = 1 / (self.imageshape[0] * 1e10 * binned_pixelsize)
# calculate radial average of the power spectrum
low = 0
high = 0
b = numextension.radialPower(pow, low, high)
# plot
indices = range(0, len(b))
rec_pixels = map((lambda x: x * rec_pixelsize), indices)
self.panelclass.setPlot(self.panel, rec_pixels, b)
def radialPlot(self, imagedata,pow):
scope = imagedata['scope']
mag = scope['magnification']
tem = scope['tem']
camera = imagedata['camera']
ccd = camera['ccdcamera']
pixelsize = self.calclient.getPixelSize(mag,tem,ccd)
if camera['binning']['x'] != camera['binning']['y']:
self.logger.warning('Unequal binning radial average not implemented')
return
shape = pow.shape
if shape[0] != shape[1]:
self.logger.warning('Non square image shape radial average not implemented')
return
binned_pixelsize = pixelsize * camera['binning']['x']
rec_pixelsize = 1/(self.imageshape[0]*1e10*binned_pixelsize)
# calculate radial average of the power spectrum
low = 0
high = 0
b = numextension.radialPower(pow, low, high)
# plot
indices = range(0,len(b))
rec_pixels = map((lambda x: x*rec_pixelsize),indices)
self.panelclass.setPlot(self.panel,rec_pixels,b)
x = args[0]
y = args[1]
fig = self.get_figure()
axes = fig.gca()
axes.plot(x,y)
def clear(self):
fig = self.get_figure()
axes = fig.gca()
axes.cla()
if __name__ == '__main__':
import wx
import numextension
from pyami import mrc
from pyami import imagefun
import sys
im = mrc.read(sys.argv[1])
pow = imagefun.power(im)
b = numextension.radialPower(pow, 5, 20)
app = wx.PySimpleApp()
frame = wx.Frame(None, -1, title='My Plot')
plotpanel = PlotPanel(frame, -1)
#b = [2,3,5,4,3]
plotpanel.plot(b)
frame.Show()
app.MainLoop()
x = args[0]
y = args[1]
fig = self.get_figure()
axes = fig.gca()
axes.plot(x, y)
def clear(self):
fig = self.get_figure()
axes = fig.gca()
axes.cla()
if __name__ == '__main__':
import wx
import numextension
from pyami import mrc
from pyami import imagefun
import sys
im = mrc.read(sys.argv[1])
pow = imagefun.power(im)
b = numextension.radialPower(pow, 5, 20)
app = wx.PySimpleApp()
frame = wx.Frame(None, -1, title='My Plot')
plotpanel = PlotPanel(frame, -1)
#b = [2,3,5,4,3]
plotpanel.plot(b)
frame.Show()
app.MainLoop()
