jueves, 7 de marzo de 2013

Detección de círculos

Para esta tarea se hicieron unas pequeñas modificaciones al código hecho en la entrada pasada, lo primero que hice fue establecer una especie de rango de radios que voy a buscar en la imagen, pero esto no es suficiente ya que genera mucho ruido, y aparece círculos donde no hay. Se ve algo así:

Entonces lo que se me ocurrió hacer es primero en los radios que hice el "for" hacer como saltos, yo los hice de 2 en 2 me pareció que eso daba un mejor resultado.
Lo segundo fue sacar nuevamente los bordes de la imagen para esto yo use otra máscara la laplaciana y con esto si era posible obtener completamente el contorno del círculo de manera más fácil, entonces para saber cuales eran los bordes los cargue en una lista, esta lista es algo así [(x1,y1), (x2, y2), ...,(xN, yN)].
Una vez teniendo esto y todos los posibles centros de los círculos que son demasiados), hice otra lista que contuviera radio y posible centro, con esto me fue posible con un "for" recorrer toda esa lista y empezar a hacer cálculos de puntos de la circunferencia, con las siguientes formula:
 xc = int(round(x - radio * cosTheta))
 yc = int(round(y - radio * sinTheta))
Entonces estos eran puntos de la circunferencia, para la medida del angulo hice otra lista con 50 ángulos aleatorios.
Entonces para cada radio con centro posibles hice el calculo anterior, y sume en un contador cada que había más de 30 puntos que coincidían con los bordes obtenidos del laplaciano.
De esta manera los que si tenían al menos 30 puntos iguales a los bordes obtenidos anteriormente si pinta el círculo para comprobar en la imagen. Y es todo :) para sacar los círculos.
Algunos ejemplos:





      


Por último lo que falta es imprimir un listado que indica para cada ID el diámetro el círculo (como porcentaje del diagonal máxima de la imagen).



Otro ejemplo:





Código relevante



 Repositorio

1 comentario: