#%%
plt.imshow(gc + 0.5*bsel, cmap='bwr')
plt.tight_layout()
#%%
plt.imshow(bsel + gsel)
#%%
plt.imshow(bsel + rsel)
# coef = ogre.calc_coloc_fromlabs(blabs, bc, gc, gmask)
#%%
from scipy.stats import pearsonr
plive = pearsonr(bcslice.flatten(), gcslice.flatten())[0]
pdead = pearsonr(bcslice.flatten(), rcslice.flatten())[0]
print 'plive = {0}, pdead = {1}'.format(plive, pdead)
plive1 = ogre.pearsonr(bcslice.flatten(), gcslice.flatten())[0]
pdead1 = ogre.pearsonr(bcslice.flatten(), rcslice.flatten())[0]
print 'plive1 = {0}, pdead1 = {1}'.format(plive1, pdead1)
plive2 = ogre.pearsonr2(bcslice, gcslice)
pdead2 = ogre.pearsonr2(bcslice, rcslice)
print 'plive2 = {0}, pdead2 = {1}'.format(plive2, pdead2)
mlive = ogre.mandersr(bcslice, gcslice)
mdead = ogre.mandersr(bcslice, rcslice)
print 'mlive = {0}, mdead = {1}'.format(mlive, mdead)
klive = ogre.overlapk(bcslice, gcslice)
kdead = ogre.overlapk(bcslice, rcslice)
print 'klive = {0}, kdead = {1}'.format(klive, kdead)
