def plotAnisotropyResolution(self, path):
md = MetaData(path)
y = md.getColumnValues(MDL_COST)
x = md.getColumnValues(MDL_RESOLUTION_SSNR)
xmax = np.amax(x)
xmin = np.amin(x)
ymax = np.amax(y)
ymin = np.amin(y)
resstep = 0.5
anistep = 0.1
from math import floor, ceil
xbin = floor((xmax - xmin) / resstep)
ybin = ceil((ymax - ymin) / 0.05)
print xbin, ybin
from matplotlib.pyplot import contour, contourf
plt.figure()
# plt.plot(x, y,'x')
plt.title('resolution vs anisotropy')
plt.xlabel("Resolution")
plt.ylabel("Anisotropy")
# counts,ybins,xbins,image = plt.hist2d(x,y,bins=100)#)[xbin, ybin])
img, xedges, yedges, imageAx = plt.hist2d(x,
y, (50, 50),
cmap=plt.cm.jet)
plt.colorbar()
plt.show()
xplotter = XmippPlotter(x=1,
y=1,
mainTitle="aaaaa "
"along %s-axis." % self._getAxis())
a = xplotter.createSubPlot("Slice ", '', '')
matrix = img
print matrix
plot = xplotter.plotMatrix(a,
matrix,
0,
200,
cmap=self.getColorMap(),
interpolation="nearest")
xplotter.getColorBar(plot)
print matrix
return [xplotter]
