Lab 5- Detección de Círculos con radios desconocidos

Para este laboratorio se encargo detectar círculos con diferentes tamaños.

Intente varias pruebas para detectar los círculos de diferentes, pensé en hacer que desde el centro de mas hacia la frontera del circulo y sacar la distancia , pero me tope con muchos problemas y no pude sacarlo a tiempo. 

Se detectan los círculos  aquí dejo el código.

	import pygame
	from pygame.locals import *
	from PIL import Image, ImageDraw,ImageFont
	import sys, random
	from math import sqrt,fabs,sin,cos,atan2, pow
	def convolucion(image,gResultant):
	ancho, altura = image.size
	pixeles = image.load()
	imageN = Image.new('RGB', (ancho,altura))
	pixels = imageN.load()
	for x in range(ancho):
	for y in range(altura):
	sumatoria = 0
	for i in range(x-1, x+2):
	for j in range(y-1, y+2):
	try:
	sumatoria += gResultant[i - (x-1)][j - (y-1)] * pixeles[i,j][1]
	except:
	pass
	pixels[x,y] = (sumatoria,sumatoria,sumatoria)
	return imageN
	def detectarcentros(imgx, imgy, image, s,um):
	ancho,altura = image.size
	frecuencia = dict()
	circulos = dict()
	pixelesx = imgx.load()
	pixelesy = imgy.load()
	for x in range(ancho):
	for y in range(altura):
	(rx,gx,bx) = pixelesx[x,y]
	(ry,gy,by) = pixelesy[x,y]
	valor=10
	gx = float(rx+gx+bx)/3
	gy = float(ry+gy+by)/3
	gxalcuadrado = pow(gx, 2)
	gyalcuadrado = pow(gy, 2)
	gResultant = sqrt((gxalcuadrado+gyalcuadrado))
	if fabs(gResultant) > 0:
	cosTheta = gx/gResultant
	sinTheta = gy/gResultant
	theta = atan2(gy,gx)
	radio=50
	centro = (int( x - radio * cos(theta)), int( y - radio * sin(theta)))
	centro = ((centro[0]/valor)*valor, (centro[1]/valor)*valor)
	circulos[x,y]=centro
	if not centro in frecuencia:
	frecuencia[centro] = 1
	else:
	frecuencia[centro] += 1
	else:
	circulos[x,y]=None
	return circulos,frecuencia
	def amarillo():
	return (255, random.randint(100,255), random.randint(0,50))
	def nombraracirc(image, frecuencia, circulos):
	pixeles = image.load()
	ancho, altura = image.size
	fuente = ImageFont.truetype('/usr/share/fonts/truetype/ubuntu-font-family/Ubuntu-LI.ttf',18)
	imagen = ImageDraw.Draw(image)
	cont = 1
	color = dict()
	aux = 2
	radio = 50
	for i in frecuencia.keys():
	color[i] = amarillo()
	imagen.ellipse((i[0]-aux, i[1]-aux, i[0]+aux, i[1]+aux), fill=(0,255,0))
	imagen.text((i[0]+aux+3, i[1]), ('Circulo '+str(cont)), fill=(0,255,0), font=fuente)
	cont +=1
	for c in color.keys():
	i ,j = c
	for angulo in range(360):
	x = i + radio*cos(angulo)
	y = j + radio*sin(angulo)
	try:
	pixeles[x,y] = color[c]
	except:
	pass
	return image
	if __name__ =="__main__":
	pygame.init()
	pantalla = pygame.display.set_mode((275,250))
	imagens = pygame.image.load("sal.png")
	image = Image.open("sal.png")
	pixeles = image.load()
	ancho, altura =image.size
	mascarax = [[-1, -1, -1], [2, 2, 2], [-1, -1, -1]]
	mascaray = [[-1, 2, -1], [-1, 2, -1], [-1, 2, -1]]
	#mx = ([-1, 0, 1], [-2, 0, 2], [-1, 0, 1]) #Para gradiente de x.
	# my = ([1, 2, 1], [0, 0, 0], [-1, -2, -1]) #Para gradiente de y.
	imgx = convolucion(image, mascarax)
	imgy = convolucion(image, mascaray)
	s=10
	um=10.0
	circulos,frecuencia = detectarcentros(imgx, imgy, image, s,um)
	maximo = 0
	suma = 0.0
	for i in frecuencia.keys():
	if frecuencia[i] < s*8:
	frecuencia.pop(i)
	else:
	print frecuencia[i], i
	nombraracirc(image,frecuencia,circulos)
	image.save('circulos.png')
	while True:
	for eventos in pygame.event.get():
	if eventos.type == pygame.QUIT:
	exit()
	pantalla.blit(imagens,(0,0))
	pygame.display.update()

view raw gistfile1.py hosted with ❤ by GitHub

X

Aquí esta la imagen de prueba y sus resultados

Lab 5 - Listado de proveedores

En esta entrada mostrare un listado de proveedores que podrían ser de mucha ayuda.

El dispositivo móvil es parte fundamental de muchos proyectos , así como lo es en nuestro proyecto. Para realizar pruebas se contara con un Samsung Galaxy Ace. 

Pero si se desea adquirir algún smartphone nuevo o diferente en la pagina de TELCEL  en el apartado equipos pueden ver que celular podría ser útil ya que se muestran las características de cada uno.Para lograr verlos equipos no hace falta nada mas que añadir nuestro estado y se mostrara esta pagina.

Aquí ya si se busca algún modelo en especifico se selecciona la marca deseada , la casilla tecnología te da opción de seleccionar 3G,4GLTE o GSM. Y en características puedes seleccionar la que estés buscando Wi-Fi. Radio etc.

Una vez que viste cual es el que deseas vas a la opción Tienda que es la tienda en línea de Telcel , seleccionas el aparato deseado y lo agregas al carrito. 

También en nuestro proyecto y en muchos de los demás equipos se necesitan LEDs, cables,resistencias , protoboards, micrófonos, cámaras, webcams etc. Y Steren podría resultar en una opción rápida de conseguir muchos dispositivos electrónicos.  Steren se especializa en Electrónica y la forma mas fácil de conseguir los productos es ir a alguna sucursal, pero en línea también se puede hacer.  

Se entra a la pagina y aparecerá una pantalla donde se debe seleccionar la opción que dice Catálogo que los mandara a esta pantalla , donde aquí están agrupados todos los dispositivos por grupos.

Basta con dar clic en la opción que se este buscando, por ejemplo algunos LEDs, pues entonces se da clic en la palabra LEDs y se mostraran los existentes ahora se le da clic en la opción comprar y se añadirá al carrito, posteriormente en la parte de arriba en la pagina dice productos en tu carrito , se le da clic y nos manda a una pantalla que solamente se palomea la opción de conformidad con las políticas de privacidad y se acepta. 

Ya simplemente se ingresan los datos personales para efectuar la compra.

Ahora otra tienda que se encuentra aquí en Monterrey es 5hz ,igual la forma mas segura de conseguir lo que se busca es acudir a la tienda, también puede comprarse en línea.

Esta tienda es proveedora de resistencias,leds, protoboards, miniprotoboards, pantallas LCD, aunque el que podemos destacar en esta tienda son los Arduinos que no maneja Steren y manejan varios tipos como el mas usado Arduino Uno.

Uno ingresa a su pagina que por cierto cabe mencionar que no es muy bonita ni con un buen diseño.

En la parte izquierda verán un listado de todas las categorías que manejan y un buscador, si se sabe el nombre del producto que se desea comprar basta con escribirlo y dar enter, en este caso mostrare el ejemplo con un Arduino uno.

Escribí Arduino Uno en el buscador , y aparecen muchas dispositivos relacionados con este , simplemente lo buscas en los resultados y lo seleccionas , aparecerá una leyenda con un texto que dice su precio , su estado ( si esta agotado o existente) y un link que dice add toa new shopping list. Aquí abajo dejo como aparece.

Precio: MX$415.00

Agotado [envíanos un correo]

Art. #: DEV5H-11021

Add to a new shopping list

Email this page to a friend

Se le da clic en Add to a new shopping list para que se agrege a nuestro carrito y al final todo lo que se desee comprar se agrega aquí  ya nada mas vana ver el carrito y dan clic en pagar , y seleccionan la forma de pago.

Otra tienda aquí en monterrey es AG Electrónica que como su nombre lo dice es especialista en electrónica , cuenta con una gran variedad de dispositivos como  adaptadores USB, microcontroladores, cautines etc.

Ingresas a la pagina AG y cuentan con un buscador se escribe lo que desea comprar en este ejemplo escribí microcontrolador ,y aparecerá una lista delo que tienen, ahora se debe hacer login , después se  escribe la cantidad que se desee agregar al carrito de compra . Posteriormente se selecciona la forma de pago. Igual es recomendable ir a la tienda personalmente.

Ahora una tienda que no existe en Nuevo León , es Sparkfun , aquí muchos hemos realizado pedidos para nuestros proyectos y esta ves quizás no sea la excepción nuevamente. Esta tienda vende arduinos, sensores, kits de robótica , dispositivos GPS etc.

Se entra a su pagina de Internet Sparkfun  y se debe hacer un registro para poder realizar la compra, de el lado izquierdo nos muestra las categorías con las que cuenta , ahora seleccionas alguna en mi caso seleccione GPS y un producto al azar, Ahora ya se agrega el producto al carrito "Add to Cart" ya que se tienen todos los productos que se pedirán se debe realizar el deposito en un banco y los productos son enviados. 

Tarea 4 - Detección de Circulos

Para esta tarea se encargo detectar círculos , utilizando un proceso similar al de las tarea anteriores

donde se detectaban lineas , pero en esta ocasión se calculara un posible centro para cada punto.

Para identificar el centro con dos ecuaciones  :

a= x-radioCosθ

b=y-radioSinθ

Se calcula el angulo y se toma ese punto suponiendo que es parte del círculo y se calcula un supuesto centro y  después se hace un calculo para ver que punto en la imagen es el que tiene mas frecuencias y de ahí se identifica el centro del círculo.

El centro del círculo  que se corresponde con el centro de la circunferencia, del cual equidistan todos los puntos de esta.

Se cálculo el gradiente para cada punto , con los gradientes "X" "Y"  como en las tareas pasadas y se aplica posteriormente estas dos formulas que dejo aquí abajo.

Este gradiente para cada punto se hace de esta manera :

Ya después  se aplican estas dos formulas:

Cosθ = gx/GradienteCalculada

Sinθ = gy/GradienteCalculada

Y después se aplican las dos formulas mencionadas arriba para calcular el posible centro para cada pixel:

a= x-radioCosθ

b=y-radioSinθ

Este fue el circuló utilizado para las pruebas 

  Aquí esta  el código:

 

	import pygame
	from pygame.locals import *
	import Image
	from math import floor, atan, fabs, pi, cos, sin, ceil
	import math
	pygame.init()
	pantalla = pygame.display.set_mode((190,250))
	# FOTO NORMAL
	imagen = pygame.image.load("salidaconvolubinariza.png")
	#CONVOLUCION A ESCALA DE GRISES Y FILTRADA
	image = Image.open("rayitas.png")
	pixeles = image.load()
	pixelesY= image.load()
	ancho, altura =image.size
	rango=100
	msobelX = ([-1, 0, 1], [-2, 0, 2], [-1, 0, 1]) #Para gradiente de x.
	msobelY = ([1, 2, 1], [0, 0, 0], [-1, -2, -1]) #Para gradiente de y.
	#prewittX=([-1, 0, 1], [-1, 0, 1], [-1, 0, 1])#EJE X PREWITT
	#prewittY=([1, 1, 1], [0, 0, 0], [-1,-1,-1])#EJE Y PREWITT
	tamanomatriz=3
	sumatoriaX = 0
	sumariaY = 0
	seleccion=raw_input("INGRESA 1 PARA SOBEL y DOS PARA PREWITT ")
	matrizagarrada=int(seleccion)
	radio=raw_input("INGRESA El radio ")
	radioselec=int(radio)
	if matrizagarrada==1:
	for x in range(altura):
	for y in range(ancho):
	sumatoriaX = 0
	sumatoriaY = 0
	if x != 0 and y != 0 and y != ancho and x != altura:
	for i in range(tamanomatriz):
	for j in range(tamanomatriz):
	try:
	gx = msobelX[i][j]*pixeles[y+j, x+i][1]
	gy = msobelY[i][j]*pixelesY[y+j, x+i][1]
	except:
	productosGX = 0
	productosGY = 0
	sumatoriaX = gx+sumatoriaX
	sumatoriaY = gy+sumatoriaY
	gxalcuadrado = pow(sumatoriaX, 2)
	gyalcuadrado = pow(sumatoriaY, 2)
	gradienteResultante = int(math.sqrt(gxalcuadrado+gyalcuadrado))
	#pixelNuevo=gradienteResultante
	#if pixelNuevo> 255:
	#pixelNuevo = 255
	#if pixelNuevo < 0:
	#pixelNuevo = 0
	#pixeles[y,x] = ( pixelNuevo, pixelNuevo,pixelNuevo)
	#imagen.save("salidaconvcirc.png")
	#pixeles[y,x] = ( gxalcuadrado, gxalcuadrado,gxalcuadrado)
	#image.save("imagenhorizo.png",'png')
	#pixelesY[y,x] = ( gyalcuadrado, gyalcuadrado,gyalcuadrado)
	#image.save("imagenverti.png",'png')
	#image.show()
	if abs(gradienteResultante)>0:
	cos=gx/gradienteResultante#coseno
	sin=gy/gradienteResultante#seno
	a = int(round(x - radioselec * cos))
	b = int(round(y - radioselec * sin))
	xcm = a
	ycm = b
	if xcm >= 0 and xcm < ancho and ycm >= 0 and ycm < altura:
	votos[ycm][xcm] += 1#aumenta uno a los votos
	#se comparan los puntos
	Maximo = 0
	sumatoria = 0.0
	#si el punto recibe muchos votos es un posible centro
	for x in range(altura):
	for y in range(ancho):
	v = votos[x][y]
	sumatoria += v
	if v > Maximo:
	Maximo = v
	#
	promedio = sumatoria / (ancho * altura)
	umbral = (Maximo+promedio)/2.0
	pixeless = []
	for x in range(altura):
	for y in range(ancho):
	v = votos[x][y]
	if v > umbral:
	pixeless.append((x,y))
	while True:
	for eventos in pygame.event.get():
	if eventos.type == pygame.QUIT:
	exit()
	pantalla.blit(imagen,(0,0))
	pygame.display.update()

view raw gistfile1.py hosted with ❤ by GitHub

El programa no logre terminarlo .