Laboratorio 5

Para esta entrega, tendremos que detectar circulos pero ahora el parametro del radio no está predefinido.

Al igual que la tarea de detección de circulos con un mismo radio, emplearemos estas ecuaciones:

Código:

	#!/usr/bin/python
	from PIL import Image, ImageDraw, ImageFont
	import sys, time, random, math
	def centrocirculo(picx, picy, image, salto = 10, bias=10.0):
	frecuencia = dict()
	circulomat = dict()
	for i in range(image.size[0]):
	for j in range(image.size[1]):
	(Rx, Gx, Bx) = picx[i,j]
	(Ry, Gy, By) = picy[i,j]
	gx = float(Rx+Gx+Bx)/3
	gy = float(Ry+Gy+By)/3
	gradiente = math.sqrt(math.pow(gx, 2) + math.pow(gy, 2))
	if gradiente < -1*bias or gradiente > bias:
	cost = (gx/gradiente)
	sint = (gy/gradiente)
	theta = math.atan2(gy, gx)
	r = 60
	centro = (int( i - r * math.cos(theta+math.radians(90.0))), int( j - r * math.sin(theta+math.radians(90.0))))
	centro = ((centro[0]/salto)*salto, (centro[1]/salto)*salto)
	circulomat[i,j] = centro
	if not centro in frecuencia:
	frecuencia[centro] = 1
	else:
	frecuencia[centro] += 1
	else:
	circulomat[i, j] = None
	return circulomat, frecuencia
	def amarillo():
	return (255, random.randint(100,255), random.randint(0, 50))
	def draw_lines(image, frecuencia, circulomat):
	pic = image.load()
	font = ImageFont.truetype("/usr/share/fonts/truetype/freefont/FreeSansBold.ttf", 18)
	draw = ImageDraw.Draw(image)
	cont = 1
	colores = dict()
	for i in frecuencia.keys():
	colores[i] = amarillo()
	r = 2
	draw.ellipse((i[0]-r, i[1]-r, i[0]+r, i[1]+r), fill=(0, 255, 0))
	draw.text((i[0]+r+3, i[1]), ('C'+str(cont)), fill=(0, 255, 0), font=font)
	print 'C'+str(cont)
	cont += 1
	for i in range(image.size[0]):
	for j in range(image.size[1]):
	if circulomat[i, j] in frecuencia:
	try:
	pic[i,j] = colores[circulomat[i, j]]
	except:
	pass
	return image
	def deteccioncirculo(image_name, output="output.png"):
	f = Filtros()
	imagex = f.aplicar_mascara(image_name, ["sobelx"], [1.0/10.0], cmd="i")
	imagey = f.aplicar_mascara(image_name, ["sobely"], [1.0/10.0], cmd="i")
	image = Image.open(image_name)
	salto = 30
	circulomat, frecuencia = centrocirculo(imagex, imagey, image, salto = salto, bias=10.0)
	for i in frecuencia.keys():
	if frecuencia[i] < salto*8:
	frecuencia.pop(i)
	else:
	print frecuencia[i], i
	draw_lines(image, frecuencia, circulomat)
	image.save(output)
	def main():
	deteccioncirculo(sys.argv[1])
	main()

view raw circulos.py hosted with ❤ by GitHub

El resultado no fue el esperado:

Laboratorio 5

Para el laboratorio 5 se nos pidió realizar un catálogo de los proveedores de electrónica que podríamos utilizar para comprar algún componente, además de revisar que software utilizaremos.

Algunos proveedores que encontré fueron los siguientes:

Adafruit

http://www.adafruit.com/

En esta página encontré este Shield GPS que incluye el Modulo GPS. Me pareció bueno ya que incluye el Shield y el Módulo.

Esta tienda se encuentra en Nueva York, en cuanto a los envíos ofrece diferentes tipos de envíos dependiendo el tiempo en que deseas que llegue, esto se realiza por paquetería UPS.

Ofrece seguros contra daños y/o perdidas, esto tiene un costo pero creo que valdría la pena pagarlo.

Además tiene servicio de devoluciones de 28 días a partir del envío.

Aquí un pdf con algunos productos que ofrece esté proveedor.

https://www.adafruit.com/catalogs/catalog.pdf

5hz-Electronica

http://www.5hz-electronica.com

Esté proveedor se encuentra en Monterrey, la venta de sus productos es mediante pedido.

Tiene garantía y servicio de devolución, también puedes recoger tu pedido personalmente sin necesidad de realizar un deposito.

Cooking Hacks

http://www.cooking-hacks.com/index.php

SeeedStudio

http://www.seeedstudio.com/depot/

http://shop.4dsystems.com.au/

Yo me enfoqué en buscar Modulos, Shield GPS, ya que mi proyecto es la alarma inteligente de automóviles y el GPS es una de las principales funciones que tendrá el proyecto.

En cuanto al software principalmente usaríamos Arduino:

Tarea 4 Detección de círculos

En esta tarea se nos pidió realizar un código que detectara círculos.

Para poder realizar esto hay que calcular para cada pixel el gradiente en “x” y en “y”

Después se calcula Cos y Sin.

Teniendo estos valores se puede calcular los posibles centros del círculo:

Código:

	def conv(img,ancho,altura):
	pixels =img.load()
	angulos = []
	gy = []
	gx = []
	matrizX =([-1,0,1],[-2,0,2],[-1,0,1])
	matrizY =([1,2,1],[0,0,0],[-1,-2,-1])
	for i in range(altura):
	gx.append([])
	gy.append([])
	for j in range(ancho):
	sumx = 0
	sumy = 0
	a=3
	for x in range(a):
	for y in range(a):
	try:
	sumx +=(pixels[j+y-1,i+x-1][0]*matrizX[x][y])
	sumy +=(pixels[j+y-1,i+x-1][0]*matrizY[x][y])
	except:
	pass
	gx[i].append(sumx)
	gy[i].append(sumy)
	grad = math.sqrt(pow(sumx,2)+pow(sumy,2))
	grad = int(grad)
	pixels[j,i] = (grad,grad,grad)
	im= img.save("conv.jpg")
	return gx, gy
	def circ(gx,gy):
	magnitud = []
	angulos = []
	rhos = []
	ancho,altura,pixels,im = cargar(imgesc)
	angulo = 0.0
	ancho2,altura2,pixels2,im2 = cargar(conv)
	coseno = 0
	seno = 0
	votos = list()
	for x in range(altura):
	votos.append([])
	for y in range(ancho):
	votos[x].append(int(0))
	for y in range(altura):
	rhos.append([])
	angulos.append([])
	magnitud.append([])
	votos.append([])
	for x in range(ancho):
	if x > 0 and y > 0 and y < altura -1 and x < ancho -1:
	hor = gx[y][x]
	ver = gy[y][x]
	g = math.sqrt(hor ** 2 + ver ** 2)
	if fabs(g) > 0.0:
	coseno = hor / g
	seno = ver / g
	xc = int(round(x -radio*coseno))
	yc = int(round(y -radio*seno))
	xcm = xc
	ycm = yc
	if xcm > 0 and xcm < ancho -1 and ycm > 0 and ycm < altura:
	votos[ycm][xcm] += 1
	pixels[xcm,ycm] = (255,0,255)

view raw circ.py hosted with ❤ by GitHub

En cuanto al resultado, no logre terminarlo.