#!/usr/bin/python
from SimpleCV import Camera, Display, Image
import numpy as np
import scipy.misc
img = Image("img4.png")
img = img.transpose()
(imgred, imggreen, imgblue) = img.splitChannels(False)
rhor = np.array([0.6650,0.7699,0.8573,0.4735,0.5482,0.6104,0.4339,0.5023,0.5593])
epsilon = 0.42
Edges = np.zeros((480,640))
x = [0,0,0,0,0,0,0,0]
for i in range (1,479):
for j in range (1,639):
[r1,g1,b1] = imgred.getPixel(i-1,j)
[r2,g2,b2] = imgred.getPixel(i+1,j)
[r3,g3,b3] = imggreen.getPixel(i-1,j)
[r4,g4,b4] = imggreen.getPixel(i+1,j)
[r5,g5,b5] = imgblue.getPixel(i-1,j)
[r6,g6,b6] = imgblue.getPixel(i+1,j)
kappahoriz1 = np.array([r1/r2,r1/g4,r1/b6,g3/r2,g3/g4,g3/b6,b5/r2,b5/g4,b5/b6])
kappahoriz2 = np.array([r2/r1,r2/g3,r2/b5,g4/r1,g4/g3,g4/b5,b6/r1,b6/g3,b6/b5])
difhor1 = kappahoriz1 - rhor
difhor2 = kappahoriz2 - rhor
for k in range (0,7):
if difhor1[k] < 0:
