###################################################################################################
#  Image processing utilities
###################################################################################################


from ijutils import *
from ch.psi.pshell.imaging.Overlays import *
import ch.psi.pshell.imaging.Pen as Pen

def in_roi(x,y):
    return math.hypot(x-roi_radius, y-roi_radius) < roi_radius


def integrate(ips):    
    roi =  get_roi()
    aux = None
    for i in range(len(ips)):
        if i==0:
            aux = new_image(roi[2], roi[3], image_type="float", title = "sum", fill_color = None)
        op_image(aux, ips[i], "add", float_result=True, in_place=True)    
    return aux

def average (ips):   
    aux = integrate(ips)     
    op_const(aux, "divide", len(ips), in_place=True)    
    return aux

def grab_frames(samples):
    frames = []            
    for i in range(samples):
        aux = get_image()        
        frames.append(aux)        
    return frames
    
def average_frames(samples = 1):
    return average(grab_frames(samples))

def integrate_frames(samples = 1):
    return integrate(grab_frames(samples))


roi_center = (800, 600)
roi_radius = 600

def get_roi():
    return (roi_center[0] - roi_radius, roi_center[1] - roi_radius, 2* roi_radius, 2*roi_radius)


def get_image():
    roi =  get_roi()
    ip = load_image(img.image)                            
    ret = sub_image(ip, roi[0], roi[1], roi[2], roi[3])
    grayscale(ret, do_scaling=True)
    return ret

#def detect_pucks(ip):
#    """
#    """
#    aux = grayscale(ip, in_place=False)
#    threshold(aux,0,50)
#    binary_fill_holes(aux)
#    return analyse_particles(aux, 10000,50000,
#            fill_holes = False, exclude_edges = True,print_table=True,
#            output_image = "outlines", minCirc = 0.4, maxCirc = 1.0)
#
#def detect_samples(ip):      
#    """
#    """  
#    aux = grayscale(ip, in_place=False)
#    invert(aux)
#    subtract_background(aux)
#    auto_threshold(aux)
#    binary_open(aux)
#    return analyse_particles(aux, 250,1000, 
#            fill_holes = False, exclude_edges = True,print_table=True,


r,g,b = [0]*256,[0]*256,[0]*256
b[0]=0xFF
b[1]=0xFF ; g[1] = 0x80; r[1] = 0x80
outline_lut1 = (r,g,b)

r,g,b = [0]*256,[0]*256,[0]*256
g[0]=0x80;r[0]=0x80;
g[1]=0xFF ; r[1] = 0x80; b[1] = 0x80
outline_lut2 = (r,g,b)