In [5]:

    

%load_ext autoreload
%autoreload 2

from matplotlib.pyplot import imshow
import matplotlib.pyplot as plt

%matplotlib inline
import numpy as np
import os
import sys
sys.path.append('../')
import gp


    

    




The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [71]:

    

# load CREMI A
input_image = []
input_prob = []
input_gold = []
input_rhoana = []
test_slices = range(0,5)# + range(25+20,50) + range(50+20,75)
for z in test_slices:
    image, prob, gold, rhoana = gp.Util.read_cremi_section(os.path.expanduser('~/data/CREMIGP/TEST/'), z)
    input_image.append(image)
    input_prob.append(255.-prob)
    input_gold.append(gold)
    input_rhoana.append(rhoana)    

original_mean_VI, original_median_VI, original_VI_s = gp.Legacy.VI(input_gold, input_rhoana)


    

In [72]:

    

import matplotlib.gridspec as gridspec
from matplotlib import rcParams
rcParams.update({'figure.autolayout': True})
left, width = .25, .5
bottom, height = .25, .5
right = left + width
top = bottom + height
fig = plt.figure(1, figsize=(30,10))

gs = gridspec.GridSpec(1,2,width_ratios=[1,1])
gs.update(wspace=0.12, hspace=0.05)

BEST95 = 448
DATA = 'cremiA'
clamper=4500

plt.subplot(gs[0])
a,b,c,d = gp.Stats.analyze_users2(['A_autoFP'], ['A_autoGP'], input_gold, input_rhoana, returnplt=True, vilabel=True, 
                        hideYlabels=False, showlegend=False,hideXlabels=False,clabel=True,hline=BEST95,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)
plt.title('Automatic', y=1.02)
plt.subplot(gs[1])
plt.title('Oracle', y=1.02)
a,b,c,d = gp.Stats.analyze_users2(['A_simuserFP'], ['A_simuserGP'], input_gold, input_rhoana, returnplt=True, vilabel=False, 
                        hideYlabels=True, showlegend=True,bestlabel=True, hideXlabels=False,clabel=True,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)

# fig.legend()

plt.savefig('/home/d/PAPERGFX/'+DATA+'_trails.pdf', bbox_inches='tight')


    

In [73]:

    

# load CREMI B
input_image = []
input_prob = []
input_gold = []
input_rhoana = []
test_slices = range(25,30)
for z in test_slices:
    image, prob, gold, rhoana = gp.Util.read_cremi_section(os.path.expanduser('~/data/CREMIGP/TEST/'), z)
    input_image.append(image)
    input_prob.append(255.-prob)
    input_gold.append(gold)
    input_rhoana.append(rhoana)    

original_mean_VI, original_median_VI, original_VI_s = gp.Legacy.VI(input_gold, input_rhoana)


    

In [74]:

    

import matplotlib.gridspec as gridspec
from matplotlib import rcParams
rcParams.update({'figure.autolayout': True})
left, width = .25, .5
bottom, height = .25, .5
right = left + width
top = bottom + height
fig = plt.figure(1, figsize=(30,10))

gs = gridspec.GridSpec(1,2,width_ratios=[1,1])
gs.update(wspace=0.12, hspace=0.05)

BEST95 = 415
DATA = 'cremiB'
clamper=3500

plt.subplot(gs[0])
a,b,c,d = gp.Stats.analyze_users2(['B_autoFP'], ['B_autoGP'], input_gold, input_rhoana, returnplt=True, vilabel=True, 
                        hideYlabels=False, showlegend=False,hideXlabels=False,clabel=True,hline=BEST95,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)
plt.title('Automatic', y=1.02)
plt.subplot(gs[1])
plt.title('Oracle', y=1.02)
a,b,c,d = gp.Stats.analyze_users2(['B_simuserFP'], ['B_simuserGP'], input_gold, input_rhoana, returnplt=True, vilabel=False, 
                        hideYlabels=True, showlegend=True,bestlabel=True, hideXlabels=False,clabel=True,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)

# fig.legend()

plt.savefig('/home/d/PAPERGFX/'+DATA+'_trails.pdf', bbox_inches='tight')


    

In [75]:

    

input_image = []
input_prob = []
input_gold = []
input_rhoana = []
test_slices = range(50,55)
for z in test_slices:
    image, prob, gold, rhoana = gp.Util.read_cremi_section(os.path.expanduser('~/data/CREMIGP/TEST/'), z)
    input_image.append(image)
    input_prob.append(255.-prob)
    input_gold.append(gold)
    input_rhoana.append(rhoana)    

original_mean_VI, original_median_VI, original_VI_s = gp.Legacy.VI(input_gold, input_rhoana)


    

In [77]:

    

import matplotlib.gridspec as gridspec
from matplotlib import rcParams
rcParams.update({'figure.autolayout': True})
left, width = .25, .5
bottom, height = .25, .5
right = left + width
top = bottom + height
fig = plt.figure(1, figsize=(30,10))

gs = gridspec.GridSpec(1,2,width_ratios=[1,1])
gs.update(wspace=0.12, hspace=0.05)

BEST95 = 512
DATA = 'cremiC'
clamper=4500

plt.subplot(gs[0])
a,b,c,d = gp.Stats.analyze_users2(['C_autoFP'], ['C_autoGP'], input_gold, input_rhoana, returnplt=True, vilabel=True, 
                        hideYlabels=False, showlegend=False,hideXlabels=False,clabel=True,hline=BEST95,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)
plt.title('Automatic', y=1.02)
plt.subplot(gs[1])
plt.title('Oracle', y=1.02)
a,b,c,d = gp.Stats.analyze_users2(['C_simuserFP'], ['C_simuserGP'], input_gold, input_rhoana, returnplt=True, vilabel=False, 
                        hideYlabels=True, showlegend=True,bestlabel=True, hideXlabels=False,clabel=True,
                        DATADIR='/home/d/CREMISTUDY/',clampX=clamper,data=DATA)

# fig.legend()

plt.savefig('/home/d/PAPERGFX/'+DATA+'_trails.pdf', bbox_inches='tight')


    

In [ ]: