Create a gist now

Instantly share code, notes, and snippets.

#!/usr/bin/python
##importamos las librerias necesarias
import pygame ##interfaz
from pygame.locals import * ##Para los botones y presion del raton
import Image ##para cargar imagenes, matriz de los pixeles
from math import * ##operaciones
import sys ##cosas de sistema como para terminar el programa
import numpy ##arreglos
import time
import math
from math import *
import random
from PIL import Image, ImageDraw,ImageFont
Inombre = raw_input("Dame el nombre de la imagen con extencio: ")
nnormal = 'nimagen.png'
ngrises = 'filtro.png'
nfiltro = 'difusion.png'
nconvolucion = 'convolu.png'
numbra = 'umbral.png'
nnorma = 'normalizacion.png'
convx = "convox.png"
convy = "convoy.png"
min = 110.0 ##parametro definido
max = 190.0 ##parametro definido
boton = pygame.Surface((100,25)) ##para crear el rectangulo del boton
boton.fill((100,100,255)) ##para el color del boton
matrii = ([0,0,0],[-1,0,0],[0,0,0])
def cargar(imagen):
im = Image.open(imagen)
ancho, altura = im.size
pixels = im.load()
return ancho,altura,pixels,im
def sii(visitados,ancho,altura):
for i in range(ancho):
for j in range(altura):
if visitados[j][i] == 0:
return j,i
return 0,0
def bordes(imagen):
colores = ([130,130,130],[0,255,0],[0,0,255],[255,255,0],[0,255,255],[255,0,255],[250,0,0]) ##colores
cola = []
ancho,altura,pixels,im = cargar(imagen)
visitados = []
contcolor = 0
for y in range(altura):
visitados.append([])
for x in range(ancho):
if pixels[x,y][1] == 0:
visitados[y].append(1)
#if len(cola) == 0:
# cola.append((y,x))
else:
if len(cola) == 0:
cola.append((y,x))
visitados[y].append(0)
while len(cola) > 0:
tam = len(cola) ##extraemos datos de la colas
a = int(cola[tam-1][0]) ##obtenemos direccion a y
b = int(cola[tam-1][1]) ##b = x
visitados[a][b] = 1 ##marcamos como visitados
cola.pop(tam-1) ##borramos de la cola
try:
if pixels[b,a+1][0] == pixels[b,a][0] and visitados[a][b+1] == 0:
#print "llego a pasar por aqui chngado"
cola.append((a,b+1))
except:
pass
try:
if pixels[b,a-1][0] == pixels[b,a][0] and visitados[a][b-1] == 0:
#print "llego a pasar por aqui chngado"
cola.append((a,b-1))
except:
pass
try:
if pixels[b+1,a][0] == pixels[b,a][0] and visitados[a+1][b] == 0:
#print "llego a pasar por aqui chngado"
cola.append((a+1,b))
except:
pass
try:
if pixels[b-1,a][0] == pixels[b,a][0] and visitados[a-1][b] == 0:
# print "llego a pasar por aqui chngado"
cola.append((a-1,b))
except:
pass
r =colores[contcolor][0]
g =colores[contcolor][1]
h =colores[contcolor][2]
pixels[b,a] = (r,g,h)
#if puntos > 0:
# plista.append((a,b))
if len(cola) == 0:
#print "se termino la figura :P lol"
#print contcolor
j,i = sii(visitados,ancho,altura)
if i==0 and j ==0:
break
else:
cola.append((j,i))
contcolor = contcolor + 1
if contcolor == 7:
contcolor = 1
im.save(pin)
print "finaliso"
return pygame.image.load(pin)
###funcion de filtro
def filtro(imagen):
prom = 0
ancho,altura,pixels,im = cargar(imagen)
for x in range(altura): ###
for i in range(ancho):
suma = 0
cont = 0
for j in range(-1,2): ##altura
for k in range(-1,2): ##ancho
if i+k >= 0 and i+k < ancho and x+j >= 0 and x+j < altura:
suma += pixels[i+k,x+j][1]
cont += 1
prom = int(suma/cont)
pixels[i,x] = (prom,prom,prom)
im.save(nfiltro)
return pygame.image.load(nfiltro)
###funcion de la doctora
def gvotos(frecs, ancho,altura):
frec = list()
for valor in histo:
if valor is None:
continue
frecuencia = frecs[valor]
acepta = False
if len(frec) <= cantidad:
acepta = True
if not acepta:
for (v, f) in frec:
if frecuencia > f:
acepta = True
break
if acepta:
frec.append((valor, frecuencia))
frec = sorted(frec, key = lambda tupla: tupla[1])
if len(frec) > cantidad:
frec.pop(0)
incluidos = list()
for (valor, frecuencia) in frec:
incluidos.append(valor)
return incluidos
##funcion de la doctora
def pints(mayor,rmin):
##ancho 1
##verticaul 0
cola = []
ancho,altura,pixels,im = cargar(nconvolucion)
cont = 0
aceptado = 0
mod = []
for k in range(len(mayor)):
x = mayor[k][1]
y = mayor[k][0]
aceptado = 0
for c in range(-rmin,rmin):
if not y == rmin:
if x+c >0 and x+c < ancho-1 and y+c >0 and y+c< altura-1:
if pixels[x+c,y][1] == 255 or pixels[x,y+c][1]==255:
cont = cont + 1
if cont > 5:
aceptado = 1
print "valor de x aceptado: ",x
print "valor de y acpetado: ",y
mod.append((y,x))
return aceptado,mod
def sinradio(rmin,rmax,gx,gy,mxy):
print "Buscando radios"
r = []
maxa = 0
rmaxa = 0
ancho,altura,pixels,im = cargar(nconvolucion)
ancho1,altura1,pixels1,im1 = cargar(ngrises)
while rmin <= rmax:
frecuencias,cord,mayor = circulos(gx,gy,rmin)
if mayor > maxa:
maxa = mayor
print "mayor",rmin
if len(cord) > 0:
print "pasado"
#maxa = mayor
#rmaxa = rmin
##print maxa
#print "rado: ",rmin
aceptado,mod = pints(cord,rmin)
if aceptado == 1:
#print "radio aceptado: ",rmin
for xx in range(len(mod)-1):
for angulo in range(360):
x = mod[xx][1] + rmin*cos(angulo)
y = mod[xx][0] + rmin*sin(angulo)
if x > 0 and y > 0 and y < altura -1 and x < ancho -1:
pixels1[x,y] = (0,0,255)
rmin = rmin + 1
im1.save(ngrises)
print "Termine"
return pygame.image.load(ngrises)
#print "valor de puntos:",maxa
#print "valor del radio:",rmaxa
def gvotos(votos,ancho,altura):
dim = int(ancho*altura)
for rango in range (1,2):
agregado = True
while agregado:
agregado = False
for y in range(altura):
for x in range(ancho):
v = votos[y][x]
if v > 0:
for dy in xrange(-rango, rango):
for dx in xrange(-rango, rango):
if not (dx == 0 and dy == 0):
if y + dy > 0 and y + dy < altura and x + dx > 0 and x + dx < ancho:
w = votos[y + dy][x + dx]
if w > 0:
if v - rango >= w:
votos[y][x] = v + w
votos[y + dy][x + dx] = 0
agregado = True
#if rango >= 200:
# print rango
return votos
def circulos(gx,gy,radio):
tiempo = time.time()
CERO = 0.0001 ##para comparaciones de los angulos
magnitud = [] ##guardamos el gradiente
angulos = [] ##guaramos angulos
rhos = [] ##guardamos rho
ancho,altura,pixels,im = cargar(ngrises)
#ancho2,altura2,pixels2,im2 = cargar(
angulo = 0.0 ##iniciamos rho como cero inicia
#normalisacion()
ancho2,altura2,pixels2,im2 = cargar(nconvolucion) ##ca
#print "inicio de rho: "
cose = 0
sen = 0
votos = list()
#print "poniendoles todos cero"
for x in range(altura):
votos.append([])
for y in range(ancho):
votos[x].append(int(0))
#print "cchekando las frecuencias"
for y in range(altura):
rhos.append([])
angulos.append([])##
magnitud.append([])##agre
votos.append([])
#xm = altura/2 - x
for x in range(ancho):
# ym = y - ancho/2
if x > 0 and y > 0 and y < altura -1 and x < ancho -1:
if pixels2[x,y][1] == 255:
hor = gx[y][x] ##dat
ver = gy[y][x] ##dato de
#mag = mxy[y][x] ##dato de gradiente combinacion
g = math.sqrt(hor ** 2 + ver ** 2)
if fabs(g) > 0.0:
#print "paso aqui si acao"
#angulo = atan(ver/hor)
cose = hor / g
sen = ver / g
xc = int(round(x -radio*cose))
yc = int(round(y -radio*sen))
xcm = xc
ycm = yc
#xcm = xc + ancho/2
#ycm = altura/2 - yc
#print "si paso aqui"
if xcm > 0 and xcm < ancho -1 and ycm > 0 and ycm < altura:
votos[ycm][xcm] += 1
pixels[xcm,ycm] = (0,255,0)
#print "chekando vecinos"
votos = gvotos(votos,ancho,altura)
maximo = 0
suma = 0.0
#print "sumando"
for x in range(altura):
for y in range(ancho):
v = votos[x][y]
suma += v
if v > maximo:
maximo = v
promedio = suma / (ancho * altura)
umbral = (promedio+maximo)/2.0
#umbral = maximo-promedio
#print "detectando"
puntosx = []
puntosy = []
mayor = 0
cord = []
mayorx = 0
mayory = 0
#print "lalalala"
for x in range(altura):
for y in range(ancho):
v = votos[x][y]
if v > umbral:
#print "pase"
#print "contador: "
#cord.append((x,y))
if v > mayor:
mayory = y
mayorx = x
mayor = v
cord.append((mayorx,mayory))
#mayor = v
#print mayor
#print "punto detectado en x:",x,"en y: ",y
#puntosx.append(y)
#puntosy.append(x)
#pixels[y,x] = (255,0,0)
#print mayor
#cord.append((mayorx,mayory))
return votos,cord,mayor
#print "valor de ancho",y,"valor de la altura",x
#for cont in range(len(puntosx)):
# print cont
# for angulo in range(360):
# x = puntosx[cont] + radio*cos(angulo)
# y = puntosy[cont] + radio*sin(angulo)
# pixels[x,y] = (0,0,255)
#im.save(ngrises)
#print "termono"
#return pygame.image.load(ngrises)
def frecuentes(histo, cantidad):
frec = list()
for valor in histo:
if valor is None:
continue
frecuencia = histo[valor]
acepta = False
if len(frec) <= cantidad:
acepta = True
if not acepta:
for (v, f) in frec:
if frecuencia > f:
acepta = True
break
if acepta:
frec.append((valor, frecuencia))
frec = sorted(frec, key = lambda tupla: tupla[1])
if len(frec) > cantidad:
frec.pop(0)
incluidos = list()
for (valor, frecuencia) in frec:
incluidos.append(valor)
print frecuencia
return incluidos
def lineas(mxy,gx,gy):
tiempoi = time.time()
CERO = 0.0001 ##para comparaciones de los angulos
magnitud = [] ##guardamos el gradiente
angulos = [] ##guaramos angulos
rhos = [] ##guardamos rho
ancho,altura,pixels,im = cargar(ngrises) ##cargamos imagen a escalade grises
#ancho2,altura2,pixels2,im2 = cargar(nconvolucion) ##cargamos la imagen de bordes
angulo = 0.0 ##iniciamos rho como cero inicia
# normalisacion()
ancho2,altura2,pixels2,im2 = cargar(nconvolucion) ##cargamos la imagen de bordes
#gy = normalisacion(gy)
print "inicio de rho: "
for x in range(altura):
rhos.append([]) ##agregamos fila
angulos.append([])##agregamos fila
magnitud.append([])##agregamos fila
for y in range(ancho):
if x > 0 and y > 0 and x < altura -1 and y < ancho -1:
hor = gx[x][y] ##dato de gradientes gx
ver = gy[x][y] ##dato de gradiente gy
mag = mxy[x][y] ##dato de gradiente combinacion
if fabs(hor) > 0.0:
angulo = atan(ver/hor)
else:
if fabs(hor) + fabs(ver) <= 0.0:#nada en ninguna direccion
angulo = None
#elif hor == 0 and ver == 255: ##para saber si es vertical
#angulo = pi/2
#print "angulo de la mitad: ",angulo
elif fabs(ver - hor) < CERO: # casi iguales diagonales
angulo = atan(1.0)
elif hor * ver > 0.0: # mismo signo horizontal
angulo = pi
else: # negativo -> -pi##horizontal
angulo = 0.0
if angulo is not None:
#cosa = '%0.2f'%angulo
# asegurar que este entre cero y pi (0 y 180) y comparaciones mas adelante
while angulo < CERO:
angulo += pi
while angulo > pi:
angulo -= pi
cosa = float('%0.2f'%angulo)
#if cosa == float('%0.2f'%pi/2):
# print "y"
#if cosa == 0.79:
# print "x"
rho = fabs( y * cos(angulo) + x * sin(angulo))
magnitud[x].append(mag)
angulos[x].append(cosa)
rhos[x].append(float('%0.2f'%rho))
#print cosa
if angulo == None:
magnitud[x].append(mag)
angulos[x].append(angulo)
rhos[x].append(None)
#angulos[x].append(int(180*(angulo/pi))/18)
# rhos[x].append(rho)
##aqu contamos las frecuencia de cada uno de las rho
grupos = dict() ##diccionario
cont = 0 ##contador para ciertas cosas como contar los vecinos si tiene el mismo angulo
#print len(magnitud)
#print altura
for x in range(altura):
for y in range(ancho):
if x > 0 and y > 0 and x < altura -1 and y < ancho -1:
# print "rhos: ",len(rhos[x])
#x print "ancho: ",ancho
# ancho2,altura2,pixels2,im2 = cargar(nconvolucion) ##cargamos la imagen de bordes
try:
if rhos[x][y] != None:
dato = ((rhos[x][y]),(angulos[x][y]))
if dato in grupos: ##decimas me falla esa part
grupos[dato] += 1 ##aqui chekamos si el
else: ##se encuentra en el d
grupos[dato] = 1 ##si no lo agrega y pone
except:
pass
print "frecuencias"
frec = frecuentes(grupos, int(ceil(len(grupos) * 50))) #
print "termino"
for x in range(altura):
for y in range(ancho):
#if pixels2[y,x][1] == 255:
if x > 0 and x< altura-1 and y > 0 and y < ancho-1:
#if pixels2[y,x][0] == 0:
try:
if pixels2[y,x][1] == 255:
#if gx[x][y] == 255 or gy[x][y] == 255:
rho, ang = rhos[x][y],angulos[x][y]
if (rho, ang) in frec:
#print "valor de theta: ",theta
if ang == 0.79:
pixels[y,x] = (0,0,255)
if ang == 3.14 or ang == 0:
pixels[y,x] = (255,255,0)
if 1.50 < ang and ang < 1.60 :
pixels[y,x] = (255,0,0)
except:
pass
print "finalizo "
im.save(ngrises) ##guardamos puntos
tiempof = time.time()
print "Tiempo estimado: ",tiempof-tiempoi
return pygame.image.load(ngrises)
####bfss para ellipse
def bfs(imagen,earray,centro):
ancho,altura,pixels,im = cargar(imagen)
pro = []
#pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 0
total = []
cord = []
print "iniciando bfs"
print len(centro)
print len(earray)
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
area = altura * ancho
for q in range(altura):
visitados.append([])
for w in range(ancho):
#total += 1
if pixels[w,q] == (0,0,0) or pixels[w,q] == (255,255,255):
visitados[q].append(0)
#print "hahahahah"
else:
visitados[q].append(1)
if len(centro) > 0:
for k in range (len(centro)):
cx,cy = centro[k]
pro.append((cx,cy))
while len(pro) > 0 :
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
if i == 0 and j == 0:
pass
else:
if pixels[x+j,y+i] != (0,0,255):
#if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
#if visitados[y+i][x+j] == 0:
# if i == 0 and j == 0:
# pass
# else:
# if (x+j,y+i) is not earray[k]:
# pro.append((x+j,y+i))
# visitados[y+i][x+j] = 1
#earray[k]
pixels[x,y] = (a,b,c)
visitados[y][x]= 1
cont += 1
if len(pro) == 0:
print "finalize figura :)"
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
total.append(cont)
cont = 0
draw = ImageDraw.Draw(im)
fuente = ImageFont.truetype('/usr/share/fonts/truetype/ubuntu-font-family/Ubuntu-C.ttf',16)
for i in range (len(centro)):
draw.text((centro[i][0], centro[i][1]), 'ellipse: ' + str(i), fill=(0,0,0), font=fuente)
print "Ellipse : ",i,"Tiene un porcentaje en la imagen de :",(total[i]*100)/area,"%."
for j in range (-10,10):
for k in range (-10,10):
pixels[centro[i][0]+j,centro[i][1]+k] = (0,0,0)
print "termino"
im.save(imagen)
return pygame.image.load(imagen)
###bfs para ellipse
### ellipse metodo
def ellipse(imagen,gx,gy,g):
ancho,altura,pixels,im = cargar(imagen)
#ancho1,altura1,pixels1,im1 = cargar(Inombre)
pro = []
#pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
cord = []
k = 0
bcenterarray = []
ccenterarray = []
centerarray = []
print "iniciando bfs"
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (255,255,255):
visitados[q].append(0)
if k == 0:
cord.append((w,q))
#print gy[q][w]
pro.append((w,q))
k = 1
else:
visitados[q].append(1)
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
cord.append((x,y))
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
#if pixels[x+j,y+i] == (0,0,0):
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
#pixels[x,y] = (a,b,c)
#pixels[x,y] = (255,255,255)
visitados[y][x] = 1
cont += 1
##w es el anchocord((w,q)) , gy[q][w]
if len(pro) == 0:
y,x = chekar(visitados,ancho,altura)
prom = (cont*100)/(total)
#print cord
#teta = arctan(gy /gx
#m=tgteta
man = len(cord)
while i < (len(cord)):
hx,hy = cord[i]
ver = gy[hy][hx]
hor = gx[hy][hx]
#if fabs(hor) < 20 or fabs(hy) <70 or (fabs(hor) < 20 and fabs(hy) < 70):
#if fabs()
if fabs(ver) > 100 or fabs(hor) < 30:
cord.pop(i)
else:
i = i + 1
tam = len(cord)
menorx =((cord[tam-1][0],cord[tam-1][1]))
mayorx = ((0,0))
menory =((cord[tam-1][0],cord[tam-1][1]))
mayory = ((0,0))
i = 0
for i in range(len(cord)):
if menorx[0] > cord[i][0]:
menorx = ((cord[i][0],cord[i][1]))
if mayorx[0] < cord[i][0]:
mayorx = ((cord[i][0],cord[i][1]))
if menory[1] > cord[i][1]:
menory = ((cord[i][0],cord[i][1]))
if mayory[1] < cord[i][1]:
mayory = ((cord[i][0],cord[i][1]))
for i in range (len(cord)):
hx,hy = cord[i]
pixels[hx,hy] = (255,0,255)
acept = 0
pos = 0
dicci = dict()
condicionr = 0
for k in range (700):
cordxx = []
cordyy = []
acept,pos = 0,0
while acept == 0 or pos == 0:
al1 = random.randint(0, tam-1)
al2 = random.randint(0, tam-1)
hx,hy = cord[al1]
hx1,hy1 = cord[al2]
if gx[hy][hx] != 0 and gx[hy1][hx1] != 0:
m = gy[hy][hx]/gx[hy][hx]
m2 = gy[hy1][hx1]/gx[hy1][hx1]
if m != 0 and acept == 0 and m < -0.01:
#if gx[hy][hx] != 0:
acept = 1
#print "posicion de x :",hx,"posicion de y : ",hy
hy2 = hy
hx2 = hx
m11 = gy[hy][hx]/gx[hy][hx]
#condicionr = 0
if m2 >= 1 and pos == 0 and acept != 0:
pos = 1
#print "posicion de x :",hx1,"posicion de y : ",hy1
m22 = gy[hy1][hx1]/gx[hy1][hx1]
break
condicionr = condicionr + 1
if condicionr == 100:
acept = 1
pos = 1
m11 = 1
m22 = 1
break
#print "primera pendiente : ",m11
condicionr = 0
for i in range (-50,50):
my = m11*(hx2+i - hx2) + hy2
if hx2+i >= 0 and my >=0 and hx2+i <ancho and my <altura:
#pixels[hx2+i,my] = (0,255,255)
cordxx.append((hx2+i,my))
#print "segunda pendiente: ",m22
for i in range (-50,50):
my = m22*(hx1+i - hx1) + hy1
if hx1+i >= 0 and my >=0 and hx1+i <ancho and my <altura:
#pixels[hx1+i,my] = (255,255,0) #my = hx+j
cordyy.append((hx1+i,my))
xxx = 0
yyy = 0
direh = 0 ##derecha ->
direv = 0 ##abajo
control = 0
if (mayory[1] + menory[1])/2 > hy2:
direv = 1 #va hacia arriba
else:
direv = 0 #va hacia abajp
for i in range(len(cordxx)):
xx,yy = cordxx[i]
if yy != 0:
if (xx,yy) in cordyy:
xxx = xx
yyy = yy
if xxx != 0:
#k = 0
kk = 0
while True:
if direv == 0:
if hx2 >= 0 and hy2-kk >=0 and hx2 < ancho and hy2-kk <altura:
if pixels[hx2,hy2-(kk)] != (255,255,255):
#if hx2 >= 0 and hy2-kk >=0 and hx2 < ancho and hy2-kk <altura:
#pixels[hx2,hy2-kk] = (255,0,0)
pass
else:
break
else:
break
else:
if hx2 >= 0 and hy2+kk >=0 and hx2 <ancho and hy2+kk <altura:
if pixels[hx2,hy2+(kk)] != (255,255,255):
#if hx2 >= 0 and hy2+kk >=0 and hx2 <ancho and hy2+kk <altura:
#pixels[hx2,hy2+kk] = (255,0,0)
pass
else:
break
else:
break
kk = kk + 1
for i in range(-200,200):
if xxx+i >= 0 and yyy >=0 and xxx+i <ancho and yyy <altura:
#pixels[xxx+i,yyy] = (255,0,255)
if xxx+i > hx2:
##acemos la votacion XD
if xxx+i in dicci:
dicci[xxx+i,yyy] += 1
else:
dicci[xxx+i,yyy] = 1
cont = 0
cord = []
if x == 0 and y == 0:
cosa = 1
else:
pro.append((x,y))
#print "valor de x: ",x,"valor de y : ",y
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
#cord = []
#print dicci
dic = dicci.items()
if len(dic) > 0:
mayor = dic[0][1]
pcx = dic[0][0]
for i in range(len(dic)):
if mayor < dic[i][1]:
mayor = dic[i][1]
pcx = dic[i][0]
#if mayor < 25:
# print "no es ellipse u circulo"
centerarray.append(pcx)
#print "punto centr: ",centerarray
#dicci = dict()
#print dicci(len(dicci)-1)
print "termino"
#im.save("bfs.png")
#im.save("centro.png")
ancho1,altura1,pixels1,im1 = cargar("convolu.png")
for i in range(len(centerarray)):
x,y = centerarray[i]
print "x",x,"y",y
suma1 = 0
suma2 = 0
k = 0
a = 0
while a == 0:
if x+k >= 0 and y >=0 and x+k <ancho and y <altura and x-k >= 0 and x-k < ancho:
if pixels1[x+k,y] == (255,255,255) or pixels1[x-k,y] == (255,255,255):
a = 1
#print "bye bye"
break
else:
suma1 = suma1 + 1
suma2 = suma2 + 1
else:
break
k = k + 1
#if k > ancho/4:
# break
k = 0
foco1 = x + suma1/2 ##derecho
foco2 = x - suma1/2
radio = x - foco2
#print "foco1 : ",foco1,"foco2: ",foco2, "valor de y :",y
umbrall = 0
for angulo in range (0,360):
puntox = x + (2*radio)*cos(angulo)
puntoy = y + radio*sin(angulo)
for pumbrall in range (-10,10):
if puntox+pumbrall >= 0 and puntoy+pumbrall >=0 and puntox+pumbrall <ancho and puntoy+pumbrall <altura:
if pixels1[puntox+pumbrall,puntoy+pumbrall] == (255,255,255):
umbrall = umbrall + 1
#for umbrally in range (radio,radio):
if umbrall > 400:
print "si es ellipse creemos "
bcenterarray.append([])
ccenterarray.append((x,y))
btam = len(bcenterarray) - 1
for angulo in range (0,360):
#print umbrall
puntox = x + (2*radio)*cos(angulo)
puntoy = y + radio*sin(angulo)
if puntox >= 0 and puntoy >=0 and puntox < ancho and puntoy <altura:
bcenterarray[btam].append((puntox,puntoy))
for h in range (-3,3):
for o in range(-3,3):
if puntox+h >= 0 and puntoy+o >=0 and puntox+h < ancho and puntoy+o <altura:
pixels1[puntox+h,puntoy+o] = (0,0,255)
for j in range (-10,10):
for g in range (-10,10):
if foco1+j >= 0 and y+g >=0 and foco1+j < ancho and y+g < altura:
pass
#pixels1[foco1+j,y+g] = (255,0,0)
if foco2+j >= 0 and y+g >=0 and foco2+j < ancho and y+g < altura:
pass
#pixels1[foco2+j,y+g] = (255,0,0)
if x+j >= 0 and y+g >=0 and x+j < ancho and y+g < altura:
pass
#pixels1[x+j,y+g] = (0,255,0)
else:
print umbrall
umbrall = 0
print "temrine el rollo por fin :)"
im1.save("centro.png")
return pygame.image.load("centro.png"),bcenterarray,ccenterarray
#return pygame.image.load("bfs.png")
### bfs metodo
def chekar(visitados,ancho,altura):
for y in range(altura):
for x in range(ancho):
if visitados[y][x] == 0:
return y,x
return 0,0
def convex_hull(imagen):
ancho,altura,pixels,im = cargar(nconvolucion)
hulls = []
mx = 0
mex = 0
my = 0
mey = 0
alle = []
print "inicio del mugrero"
for i in range(altura):
for j in range(ancho):
if pixels[j, i][1] == 255:
#print pixels[j,i]
cord = bfss(j,i)
#print "i",i
for k in range(len(cord)):
#print "eee"
cx = cord[k][0]
cy = cord[k][1]
#print "cx,",cx,"cy",cy
pixels[cx , cy] = (255,0,0)
#print "valor de qui: ",pixels[cx,cy]
#print "aquii"
mex,mx,mey,my = cord[0][0],cord[0][0],cord[0][1],cord[0][1]
xxx = 0
yyy = 0
xx = 0
yy = 0
for x in range(len(cord)):
if cord[x][0] > mx:
mx = cord[x][0]
xxx = cord[x][1]
if cord[x][0] < mex:
mex = cord[x][0]
xx = cord[x][1]
if cord[x][1] >my:
my = cord[x][1]
yyy = cord[x][0]
if cord[x][1] <mey:
mey = cord[x][1]
yy = cord[x][0]
puntos = [(mx,xxx),(mex,xx),(my,yyy),(mey,yy)]
#puntos = proceso(mx,mex,my,mey,cord)
#for o in range(len(puntos)):
# alle.append((puntos[o][0],puntos[o][1]))
print "termino"
im.save("HULLS.png", 'png')
print len(alle)
return pygame.image.load("HULLS.png"),puntos ##regresemos la imagen cargada
def proceso(mx,mex,my,mey,cord):
puntos = []
for x in range(len(cord)):
cx,cy = cord[x][0],cord[x][1]
if cx == mx:
puntos.append((cx,cy))
if cx == mex:
puntos.append((cx,cy))
if cy == my:
puntos.append((cx,cy))
if cy == mey:
puntos.append((cx,cy))
return puntos
def bfsss(j,i):
ancho,altura,pixels,im = cargar(nconvolucion)
pro = []
pro.append((j,i))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
cord = []
#print "iniciando bfs"
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (0,0,0):
visitados[q].append(0)
else:
visitados[q].append(1)
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
#if pixels[x+j,y+i] == (0,0,0):
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
cord.append((x,y))
pixels[x,y] = (a,b,c)
visitados[y][x] = 1
cont += 1
return cord
## convolucion
def convolucion(imagen):
tiempoi = time.time()
print "Iniciando convolucion: "
ancho,altura,pixels,im = cargar(imagen)
matrix = ([-1,0,1],[-2,0,2],[-1,0,1])
matriy = ([1,2,1],[0,0,0],[-1,-2,-1])
nueva = Image.new("RGB", (ancho, altura))
npixels = nueva.load()
gx = [] ##gradientes de x
gy = [] ##gradientes de y
mxy = [] ##la convinacon de ambos
for i in range(altura):
gx.append([])
gy.append([])
mxy.append([])
for j in range(ancho):
sumax = 0
sumay = 0
for y in range(-1,2):
for x in range(-1,2):
if j+x >= 0 and j+x < ancho and i+y >=0 and i+y < altura:
sumax += matrix[y+1][x+1] * pixels[j+x,i+y][1]
sumay += matriy[y+1][x+1] * pixels[j+x,i+y][1]
r = int(math.sqrt(sumax**2+sumay**2)) ##calculamos gradiente mayor
gx[i].append(sumax) ##guardamos los gradientes en x
gy[i].append(sumay) ##guardamos los gradientes en y
mxy[i].append(r) ##guardamos el resultado de los gradientes en x e y
#r = sumay
if r < 0:
r = 0
if r > 255:
r = 255
npixels[j, i] = (r,r,r)
nueva.save(nconvolucion)
tiempof = time.time()
print "Se tardo: ",tiempof-tiempoi,"segundos"
return pygame.image.load(nconvolucion),gx,gy,mxy
##funcion para la escala de grises
def escala(imagen):
ancho,altura,pixels,im = cargar(imagen)
for i in range(ancho):
for j in range(altura):
##sacamos promedio de (r,g,b) de los pixeles
(a,b,c) = pixels[i,j]
suma = a+b+c
prom = int(suma/3)
a = prom ##igualamos
b = prom ##igualamos
c = prom ##igualamos
pixels[i,j] = (a,b,c) ##igualamos
im.save(ngrises) ##guardamos la imagen nueva
return pygame.image.load(ngrises) ##enviamos la imagen cargada
##funcion para el umbral
def mumbra(imagen):
ancho,altura,pixels,im = cargar(imagen)
for i in range(ancho):
for j in range(altura):
(a,b,c) = pixels[i,j] ##obtenemos valore
prom = a+b+c/3 ##promedio
if min >= prom: #condiciones para poner blanco o negro
prom = 0
if max <= prom:
prom = 255
a = prom ##asignamos el promedio a los pixeles
b = prom
c = prom
pixels[i,j]=(a,b,c) ##asiganmos
im.save(numbra) ##guardamos
return pygame.image.load(numbra) ##regresemos la imagen cargada
##funcio para la diferenciacion entre la original y la imagen filtrada
def normalisacion(imagen):
print "inicia normalizacion"
tiempoi = time.time()
ancho,altura,pixels,im = cargar(imagen)
prom = 0
maxi = pixels[0,0][1]
mini = pixels[0,0][1]
for j in range(altura):
for i in range(ancho):
if maxi < pixels[i,j][1]:
maxi = pixels[i,j][1]
if mini > pixels[i,j][1]:
mini = pixels[i,j][1]
div = 256.0/(maxi-mini)
for j in range(altura):
for i in range(ancho):
prom = 0
if i > 0 and j > 0 and i < ancho-1 and j < altura-1:
prom = int(math.floor((pixels[i,j][1] - mini)*div))
pixels[i,j] = (prom,prom,prom)
if pixels[i,j][1] > 60:
pixels[i,j] = (255,255,255)
else:
pixels[i,j] = (0,0,0)
else:
pixels[i,j] = (0,0,0)
print "terminino normalizacion"
print "Binarisacion"
im.save(nconvolucion) ##guardamos imagenes
return pygame.image.load(nconvolucion)
def diferencia(imagen1,imagen2,min,max):
ancho,altura,pixels,im = cargar(imagen1)
ancho,altura,pixels,im = cargar(imagen2)
for i in range (ancho -2):
for j in range (altura -2):
(a,b,c) = pixels[i,j] ##guardamos pixeles de ambas
(a1,b1,c1) = pixels1[i,j] ##imagenes
prom = a+b+c/3 ##sacamos propedio de ambas
prom1 = a1+b1+c1/3
prom2 = prom - prom1 ##restamos
v = prom2*255/500
print v
if v > 25: ##aqui asemos binarisasion
prom2 = 255
else:
prom2 = 0
a = prom2 ##reasignamos la diferencia
b = prom2
c = prom2
pixels[i,j] = (a,b,c) ##cambiamos pixeles
im.save(nimagen) ##guardamos imagenes
return pygame.image.load(nimagen) ##enviamos imagenes cargadas
#def normalizacion():
def ollos(imagen):
ancho,altura,pixels,im = cargar(imagen)
ancho1,altura1,pixels1,im1 = cargar(imagen)
histogramah = []
histogramav = []
suma = 0
promh = 0
promv = 0
##pixels [ancho,largo]
# pixels[i,j]
##horizontal
for i in range(altura):
suma = 0
for j in range(ancho):
pix = pixels[j,i][1]
suma += pix
histogramah.append(suma)
##horizontal
##vartical
for i in range (ancho):
suma = 0
for j in range (altura):
pix = pixels[i,j][1]
suma += pix
histogramav.append(suma)
##vertical
suma = 0
for i in range (len(histogramah)):
suma += histogramah[i]
promh = suma/(len(histogramah))
suma = 0
for i in range (len(histogramav)):
suma += histogramav[i]
promv = suma/(len(histogramav))
print promh
print promv
#controlador = 0
anterior = histogramah[0]
for i in range (len(histogramah)-1):
sigue = histogramah[i+1]
if anterior > histogramah[i] < sigue and histogramah[i] < promh-10:
for j in range (ancho):
pixels1[j,i] = (0,0,255)
anterior = histogramah[i]
anterior = histogramav[0]
for i in range (len(histogramav)-1):
sigue = histogramav[i+1]
if anterior > histogramav[i] < sigue and histogramav[i] < promv-10:
for j in range (altura):
pixels1[i,j] = (255,0,255)
anterior = histogramav[i]
#print histogramah
im1.save("lineas.png")
pos = []
anterior = histogramah[0]
anterior1 = histogramav[0]
for i in range (len(histogramah)-1):
sigue = histogramah[i+1]
if anterior > histogramah[i] < sigue and histogramah[i] < promh-20:
for j in range (len(histogramav)-1):
sigue1 = histogramav[j+1]
if anterior1 > histogramav[j] < sigue1 and histogramav[j] < promv-20:
pixels[j,i] = (255,0,255)
pos.append((j,i))
anterior1 = histogramav[j]
anterior = histogramah[i]
archivoh=open('datosh.txt','w')
archivov=open('datosv.txt','w')
for i in range (len(histogramah)):
archivoh.write(''+ str(i) +' '+str(histogramah[i])+'\n')
for i in range (len(histogramav)):
archivov.write(''+ str(i) +' '+str(histogramav[i])+'\n')
archivoh.close()
archivov.close()
im.save("agujeros.png")
return pygame.image.load("agujeros.png"),pos
###bfs para los ollos
def pintoyos(imagen,pos):
ancho,altura,pixels,im = cargar(imagen)
pro = []
pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
cord = []
cx = 0
cy = 0
actma = 0
lisc = []
lisc.append((0,0))
controlador = 0
lugar = 0
#draw = ImageDraw.Draw(im)
#fuente = ImageFont.truetype('/usr/share/fonts/truetype/ubuntu-font-family/Ubuntu-C.ttf',40)
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
print "iniciando bfs"
print "iniciando el otro rollo"
visitados = []
pro = []
pro.append((lisc[lugar][0],lisc[lugar][1]))
print "iniciando bfs"
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (0,0,0):
visitados[q].append(0)
else:
visitados[q].append(1)
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
#if pixels[x+j,y+i] == (0,0,0):
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
pixels[x,y] = (255,0,255)
visitados[y][x]= 1
print "termino el otro rollo"
#visitados[y][x]= 1
#cont += 1
##################################################################################################################
pro = []
visitados = []
print "Contando posibles Agujeros"
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (0,0,0):
visitados[q].append(0)
else:
visitados[q].append(1)
#pro.append((pos[0][0],pos[0][1]))
print "valor de llitos: ",len(pos)
cambio = 0
ancho1,altura1,pixels1,im1 = cargar("1.JPG")
while len(pos) > 0:
pro.append((pos[0][0],pos[0][1]))
pos.pop(0)
controlador = 0
contador = 0
while len(pro) > 0:
x,y = pro[0][0],pro[0][1]
pro.pop(0)
for i in range(-2,3):
for j in range(-2,3):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
visitados[y+i][x+j] = 1
pro.append((x+j,y+i))
acept,borror = buscar(pos, (x+j,y+i))
if acept == 1:
controlador = 0
pos.pop(borror)
pass
else:
controlador += 1
if controlador == 200:
for k in range(len(pro)):
x,y = pro[k][0],pro[k][1]
if x >= 0 and y >=0 and x <ancho and y <altura:
visitados[y][x] = 1
pixels[x,y] = (255,0,255)
pro = []
constancia = 1
controlador = 0
#print "no es y ya"
break
pixels[x,y] = (255,0,255)
contador += 1
#cambio += 1
#print "cambio ",cambio
im.save("ollos.png")
###parte final#####################################################
colorear = []
ancho,altura,pixels,im = cargar("ollos.png")
pro = []
visitados = []
tollos = 1
contblanco= 0
contnegro = 0
contadorlinea = 0
aollo = 0
print "checando Agujeros"
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (255,0,255):
visitados[q].append(0)
else:
visitados[q].append(1)
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
pro.append((0,0))
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if pixels[x+j,y+i] == (255,255,255):
contblanco += 1
contadorlinea += 1
#else:
if pixels[x+j,y+i] == (0,0,0):
contnegro += 1
contadorlinea += 1
#contadorlinea += 1
if visitados[y+i][x+j] == 0:
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
#pixels[x,y] = (a,b,c)
colorear.append((x,y))
visitados[y][x]= 1
cont += 1
if len(pro) == 0:
promedio = float(contblanco) / float(contadorlinea)
#print "aaa: ",promedio
#if promedio > 0.30:
if promedio > 0.60 and promedio < 0.90:
#print "promedio ",promedio
#print "blanco: ",contblanco
#print "linea: ",contadorlinea
mayorx = colorear[len(colorear)-1]
menorx = colorear[0]
mayory = colorear[len(colorear)-1]
menory = colorear[0]
cxm = 0
cym = 0
cyn = 0
cxn = 0
#draw.text((centro[i][0], centro[i][1]), 'ellipse: ' + str(i), fill=(0,0,0), font=fuente)
for i in range(len(colorear)):
x,y = colorear[i]
if mayorx[0] < x:
mayorx = colorear[i]
if menorx[0] > x:
menorx = colorear[i]
if mayory[1] < y:
mayory = colorear[i]
if menory[1] > y:
menory = colorear[i]
cx = (mayorx[0] + menorx[0])/2
cy = (mayory[1] + menory[1])/2
#pixels[x,y] = (a,b,c)
#pixels1[x,y] = (255,0,255)
aollo += 1
area = altura*ancho
porcentaje = (float(aollo)/float(area))*100.0
if porcentaje > 1:
#x,y = colorear[i]
for p in range(len(colorear)):
x,y = colorear[p]
pixels1[x,y] = (255,0,255)
aollo += 1
#print "Agujero ",tollos,"porcentaje del agujero",(float(aollo)/float(area))*100.0
draw = ImageDraw.Draw(im1)
fuente = ImageFont.truetype('/usr/share/fonts/truetype/ubuntu-font-family/Ubuntu-C.ttf',20)
draw.text((cx, cy), 'Agujero: ' + str(tollos), fill=(255,0,0), font=fuente)
for l in range(-2,2):
for u in range (-2,2):
if cx+u >= 0 and cy+l >=0 and cx+u <ancho and cy+l <altura:
pixels1[cx+u,cy+l] = (255,255,0)
area = ancho*altura
print "Agujero ",tollos,"porcentaje del agujero",(float(aollo)/float(area))*100.0
tollos += 1
#draw.text((colorear[0][0], colorear[0][1]), 'ollo: ', fill=(255,0,0), font=fuente)
aollo = 0
contblanco= 0
contnegro = 0
contadorlinea = 0
colorear = []
promedio = 0
y,x = chekar(visitados,ancho,altura)
if x == 0 and y == 0:
cosa = 1
else:
pro.append((x,y))
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
contblanco= 0
contnegro = 0
contadorlinea = 0
im1.save("final.png")
im.save("ollos.png")
print "termino"
return pygame.image.load("final.png")
###bfs para los ollos
###solo para buscar elementos
def buscar(lista, elemento):
for i in range(len(lista)):
if lista[i] == elemento:
return 1,i
return 0,0
## solo para buscar elemtos
###
def comprobacion(imagen,gx,gy):
ancho,altura,pixels,im = cargar(imagen)
pro = []
pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
cord = []
cond = 0
print "iniciando bfs"
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (255,255,255):
visitados[q].append(0)
if cond == 0:
cond = 1
pro.append((w,q))
else:
visitados[q].append(1)
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
pixels[x,y] = (a,b,c)
visitados[y][x]= 1
cont += 1
cord.append((x,y))
if len(pro) == 0:
formas(imagen,gx,gy,cord)
y,x = chekar(visitados,ancho,altura)
#prom = (cont*100)/(total)
#print "Porcentaje del objeto: ",prom,"%"
cont = 0
cord = []
if x == 0 and y == 0:
cosa = 1
else:
pro.append((x,y))
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
#im.save("bfs.png")
return pygame.image.load(imagen)
#############################################################################################################################
###########################################################deteccion de formas #################################################
def formas(imagen,gx,gy,cord):
ancho,altura,pixels,im = cargar(imagen)
pro = []
#pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
cordenadas = []
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
cond = 0
print "iniciando verificacion:"
yep = 0
print len(cord)
draw = ImageDraw.Draw(im)
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
yep = 0
try:
l = cord.index((w,q))
visitados[q].append(0)
if cond == 0:
cond = 1
pro.append((w,q))
except:
visitados[q].append(1)
#if pixels[w,q] == (255,255,255):
# visitados[q].append(0)
# if cond == 0:
# cond = 1
# pro.append((w,q))
#else:
# visitados[q].append(1)
tras = []
kk = 0
print "todo bien"
# for i in range (len(cord)):
# for j in range(len(visitados)):
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
tras.append([])
###nueva pendiente
if gx[y][x] != 0:
m = gy[y][x]/gx[y][x]
else:
m = gy[y][x]/1
##nueva pendiente
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
#m = gy[hy][hx]/gx[hy][hx]
if i == 0 and j == 0: ##para no tomar el mismo valor
cosa = 0
else:
if gx[y+i][x+j] != 0:
m2 = gy[y+i][x+j]/gx[y+i][x+j]
else:
m2 = gy[y+i][x+j]/1
if pixels[x,y][1] == pixels[x+j,y+i][1]: ##siga siendo borde
if m-3 <= m2 <= m+3:
#if m == m2:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
#pixels[x,y] = (a,b,c)
visitados[y][x]= 1
cont += 1
cordenadas.append((x,y))
tras[kk].append((x,y))
#yep += 1
#kk += 1
if len(pro) == 0:
kk = kk + 1
cont = 0
cond = 0
k = 0
y,x = chekar(visitados,ancho,altura)
if x == 0 and y == 0:
cosa = 1
cordenadas = []
else:
cordenadas = []
pro.append((x,y))
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
#ancho1,altura1,pixels1,im1 = cargar("ev.jpeg")
##chekamos linea por linea
cosa = []
medio = []
for i in range (len(tras)-1):
print len(tras[i])
if len(tras[i]) > 100:
#print len(tras[i])
puntosmedios = []
#medio = []
menorx,menory = tras[i][0]
mayorx,mayory = tras[i][len(tras[i])-1]
lines = []
#x,y = tras[i][j]
#cosa.append((x,y))
##calculamos puntos medios
for j in range (len(tras[i])):
yep += 1 ## eso es para le porcentaje usados mas adelante
x,y = tras[i][j]
cosa.append((x,y))
if menorx > x:
menorx,menory = tras[i][j]
if mayorx < x:
mayorx,mayory = tras[i][j]
#cosa.append((x,y))
mediox,medioy = ((mayorx+menorx)/2,(mayory+menory)/2) ##calculamos medio
puntosmedios.append((mediox,medioy)) ##los agregamos a una lista para uso posterior
lines.append((menorx,menory)) ##estos es para juntar y dibujar
lines.append((mayorx,mayory)) ##esto es para juntar y dibujar
medio.append((mayorx,mayory,menorx,menory)) ##esto es para calcular el centro
#draw = ImageDraw.Draw(im)
draw.line(lines, fill=500) ##dibujamos lineas
##aqui dibujamos laas lienas
##dibujamos puntos medios
for j in range (-10,11):
for k in range (-10,11):
pixels[mediox+k,medioy+j] = (255,255,0)
##aqui calculamos centro de enmedio
menorx,menory = cosa[0]
mayorx,mayory = cosa[len(cosa)-1]
for i in range(len(cosa)):
x,y = cosa[i]
if x > mayorx:
mayorx = x
if x < menorx:
menorx = x
if y > mayory:
mayory = y
if y < menory:
menory = y
centrox,centroy = ((mayorx+menorx)/2,(mayory+menory)/2)
for i in range (-10,10):
for j in range (-10,10):
pixels[centrox+i,centroy+j] = (255,0,255)
###movemos todo el mugrero hajshajhsjahsjahjsa de aqui es para mover las lineas
hor = 0
ver = 0
print "total de lineas: ",len(medio)
for i in range (len(medio)):
mx,my,mex,mey = medio[i]
mediox,medioy = ((mx+mex)/2,(my+mey)/2)
hor = 0
ver = 0
lines = []
diferenciax = centrox - mediox
diferenciay = centroy - medioy
mx = mx + diferenciax
my = my + diferenciay
mex = mex + diferenciax
mey = mey + diferenciay
print "cordenadas = ",(mx,my,mex,mey)
#lines = []
lines.append((mex,mey))
lines.append((mx,my))
draw.line(lines, fill=100)
porcentaje = int(float(yep*100)/float(len(cord)))
print "porcentaje hacer figura: ",porcentaje
fuente = ImageFont.truetype('/usr/share/fonts/truetype/ubuntu-font-family/Ubuntu-C.ttf',20)
draw.text((centrox, centroy), 'Porcen Poligono: ' + str(porcentaje)+'%', fill=(255,0,0), font=fuente)
##movemos todo el mugrero ###
#im.save("bfs.png")
im.save(imagen)
#im1.save("ev.jpeg")
print "termino"
return
#return pygame.image.load("bfs.png")
#a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
###########################################################deteccion de formas ##################################################
##################################################---esquinas ---##################################
def esquinas(imagen):
ancho,altura,pixels,im = cargar(imagen)
for i in range(ancho):
for j in range(altura):
lista = []
for k in range(-1,2):
for l in range(-1,2):
if i+k >= 0 and j+l >= 0 and i+k < ancho and j+l < altura:
lista.append(pixels[i+k,j+l][1])
lista.sort()
rec = len(lista)/2
if len(lista) % 2 == 0:
c = (lista[rec - 1] + lista[rec])/2
else:
c = lista[rec]
pixels[i,j] = (c,c,c)
im.save("filtrom.png")
ancho1,altura1,pixels1,im1 = cargar(imagen)
cordenadas = []
for i in range(ancho):
for j in range(altura):
diferencia = pixels[i,j][1] - pixels1[i,j][1]
if diferencia < 0:
diferencia = 0
if diferencia > 255:
diferencia = 255
if diferencia > 25:
pixels1[i,j] = (diferencia,diferencia,diferencia)
cordenadas.append((i,j))
else:
pixels1[i,j] = (0,0,0)
im1.save("esquina.png")
return pygame.image.load("esquina.png"),cordenadas
###################################################--esquinas--###########################################################
##################################################--uniones--#############################################################
def union(imagen,cord):
ancho,altura,pixels,im = cargar(imagen)
ancho1,altura1,pixels1,im1 = cargar("esquina.png")
pro = []
draw = ImageDraw.Draw(im1)
#pro.append((0,0))
visitados = []
pintar = []
visitados = []
cont = 1
total = 0
#cord = []
print "iniciando bfs"
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
for q in range(altura):
visitados.append([])
for w in range(ancho):
total += 1
if pixels[w,q] == (255,255,255):
visitados[q].append(0)
if cont == 1:
pro.append((w,q))
cont = 0
print "seleccionando uno :P"
else:
visitados[q].append(1)
#identp = 0
union = []
print "valores de las esquinas",cord
while len(pro) > 0:
x = pro[len(pro)-1][0]
y = pro[len(pro)-1][1]
pro.pop(len(pro)-1)
######################
try:
l = cord.index((x,y))
identp = 1
union.append((x,y))
print "esquina encontrada"
except:
pass
for i in range(-1,2):
for j in range(-1,2):
if x+j >= 0 and y+i >=0 and x+j <ancho and y+i <altura:
if visitados[y+i][x+j] == 0:
if i == 0 and j == 0:
cosa = 0
else:
if pixels[x,y][1] == pixels[x+j,y+i][1]:
pro.append((x+j,y+i))
visitados[y+i][x+j] = 1
pixels[x,y] = (a,b,c)
visitados[y][x]= 1
cont += 1
if len(pro) == 0:
y,x = chekar(visitados,ancho,altura)
prom = (cont*100)/(total)
print "Porcentaje del objeto: ",prom,"%"
cont = 0
#print "union :",union
#union.append(union[0])
print "union :",union
draw.line(union, fill=500)
union = []
if x == 0 and y == 0:
cosa = 1
else:
pro.append((x,y))
a,b,c = random.randint(1, 255),random.randint(1, 255),random.randint(1, 255)
im1.save("union.png")
im.save("bfs.png")
print "termino"
return pygame.image.load("union.png")
##################################################--uniones--#############################################################
##funcion principla del programa
def main(nombreI):
pygame.init() ##inicialisamos python
ancho,altura,pixels,im = cargar(nombreI)
pantalla = pygame.display.set_mode((ancho + 150,altura)) ##inicialisamos pantalla
pygame.display.set_caption("Rutinas para imagenes") ##ponemos nombre
imagen = pygame.image.load(nombreI) ##cargamos imagen principal a mostrar
fuente = pygame.font.Font(None, 20) ##para la fuente de las letras
grisle = 'Escala de grices' ##variable para guardar las letras a mostrar y cambiar posteriormente
umbral = fuente.render('Umbral',1,(255,255,255)) ##para poner la fuente a variable o predefinido
normal = fuente.render('Normal',1,(255,255,255)) ##igual que la de arriba
normali = fuente.render('Conve xhull',1,(255,255,255))
cont = 0 ##para diferentes acciones de los botones psoteriores
##ciclo principal
#lista = []
while True:
esgris = fuente.render(grisle,1,(255,255,255)) ##poner fuente a la letras grisle anteriormente declarada arriba
for event in pygame.event.get(): ##verificar si existe algun evento
if event.type == MOUSEBUTTONDOWN: ##chekar si se presiono el boton del raton
cordx, cordy = pygame.mouse.get_pos() ##tomar cordenadas del puntero del raton despues de ser presionado
##primer evento para el primer boton
if 0 < cordx < 100 and 0 < cordy < 30:
##se usa condiciones para cambiar la letra
## adenmas de hacer escala de grises y posteriormente el mismo boton aplica filtro
##la tercera condicion es para aplicar la diferencia para que sala borrosa la imagen
if cont == 0:
imagen = escala(nombreI) ##hace llamar a la funcion de escala y garda el resultado
grisle = 'probando'
#imagen,pos = ollos(ngrises)
#imagen = filtro(ngrises)
imagen,cordenadas = esquinas(ngrises)
imagen = normalisacion("esquina.png")
#imagen = normalisacion(ngrises)
imagen = filtro(ngrises)
imagen,gx,gy,mxy = convolucion(nfiltro)
imagen = normalisacion(nconvolucion)
imagen = union(nconvolucion,cordenadas)
#imagen = comprobacion(nconvolucion,gx,gy)
#imagen = formas(nconvolucion,gx,gy,mxy)
#imagen = pintoyos(nconvolucion,pos)
if cont == 1:
imagen = filtro(ngrises) ##lo mismo que arriba pero para filtro
#imagen,pos = ollos(ngrises)
grisle = 'convolucion'
if cont == 2:
imagen,gx,gy,mxy = convolucion(nfiltro) ##lo mismo de arriba para difucion
grisle = 'normalizacion'
if cont == 3:
imagen = normalisacion(nconvolucion,mxy)
grisle = 'ellipse'
if cont == 4:
imagen,earray,centro = ellipse(nconvolucion,gx,gy,mxy) ##para los contornos
#imagen = lineas(mxy,gx,gy) ##lo mismo de arriba para$
# grisle = 'cosas de mas'
grisle = 'bfs'
# cont = cont + 1
if cont == 5:
cosai = "centro.png"
imagen = bfs(cosai,earray,centro)
if cont == 6:
#imagen = lineas(mxy,gx,gy)
imagen = circulos(gx,gy,45)
if cont == 7:
imagen = pygame.image.load(nombreI)
cont = 0
cont = cont + 1
##accion para el boton de umbral
if 0 < cordx < 100 and 30 < cordy < 60:
imagen = mumbra(nombreI) ##llama al metodo umbral
if 0 < cordx < 100 and 60 < cordy < 90:
imagen = pygame.image.load(nombreI) ##carga la imagen normal u original
grisle = 'Escala de grices'
cont = 0
if 0 < cordx < 100 and 90 < cordy < 120:
imagen,alle = convex_hull(nconvolucion)
## esto es para quitar la ventan y salir del programa
if event.type == pygame.QUIT:
sys.exit()
pantalla.fill((0,0,0)) ##borrar todo lo grafico
pantalla.blit(boton,(0,0)) ##colocar boton 1
pantalla.blit(boton,(0,30)) ##colocar boton 2
pantalla.blit(boton,(0,90)) ##boton 4
pantalla.blit(umbral,(0,30)) ##coloar letra umbral
pantalla.blit(boton,(0,60)) ##colocar boton 3
pantalla.blit(normal,(0,60)) ##colocar letra normal
pantalla.blit(esgris,(0,0)) ##colocar letra gris
pantalla.blit(normali,(0,90))
pantalla.blit(imagen, (100,0)) ##ppner la imagen actual
#try:
# pygame.draw.lines(pantalla, (255,0,0), False, alle)
# tama = len(puntos)-1
# pygame.draw.line(pantalla,(0,0,0),(puntos[0][0],puntos[0][1]),(puntos[tama][0],puntos[tam][1]))
#except:
# pass
pygame.display.update() ##refrescamos la pantalla con los nuevos elemntos
main(Inombre) ##llamamos la funcion principa lol
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment