#!/usr/bin/env python
# coding: utf-8

# # Thresholding

# In[1]:


import os
import numpy as np
import skimage
import matplotlib.pyplot as plt
import pandas as pd
import urllib
#import ipyvolume as ipv
import ipywidgets as ipw
if not os.path.isfile('Samples/confocal2.tif'):
    if not os.path.isdir('Samples'):
        os.mkdir('Samples')
    #url = 'https://cildata.crbs.ucsd.edu/media/images/10494/10494.tif'
    url = 'https://cildata.crbs.ucsd.edu/media/images/41694/41694.tif'
    urllib.request.urlretrieve(url, 'Samples/confocal2.tif');

image = skimage.io.imread('Samples/confocal2.tif')


# In[2]:


def thresholding(level, image):
    fig = plt.figure(figsize=(10,10))
    ax1 = plt.subplot(1,2,1)
    plt.imshow(image[:,:,2], cmap = 'gray')
    ax1 = plt.subplot(1,2,2)
    plt.imshow((image[:,:,2]>level).astype(int), cmap = 'gray',vmin = 0, vmax = 1)
    plt.show()


# In[3]:


ipw.interactive(thresholding, level = ipw.IntSlider(min =0, max = 255), 
                image = ipw.fixed(image))


# In[ ]: