Se desea que solo se grafique/dibujen los puntos ('+') y parábola dentro de la zona negra que está en el rectángulo.
Lo que se supone que se haciendo es que, se está generando la lista a lo largo del rectángulo para que de ésta forma se delimite un poco el área y que apesar de que haya más zonas negras solo se tome la del centro (la boquita) ya que se "se está delimitando" con el rectángulo.
Ésta primera parte del código según yo no están importante solo es para que tengán el código completo y se pueda visualizar como es que se está desarrollando todo.
Parte uno del código:
from matplotlib import pyplot as plt
import cv2
import numpy as np
import random
image = cv2.imread('radio.jpg',0)
img = cv2.resize(image,(600,300))
hist = cv2.calcHist([img],[0],None,[256],[0,256])
b=sum(hist)
objetivo = b/2
suma=0
numero=0
final=0
for i in hist:
suma=suma+i
numero=numero+len(i)
final=numero+1
#print('Final:',final)
if suma > objetivo:
break
U=final-3
ResFinal=U*2
y=((30*ResFinal)/100)
ret,th1 = cv2.threshold(img,y,255,cv2.THRESH_BINARY)
contours, _ = cv2.findContours(th1, cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE)
lista = []
for c in contours:
area = cv2.contourArea(c)
if area > 3000 and area < 10000:
cv2.drawContours(img, [c], 0, (0, 255, 0), 2, cv2.LINE_AA)
#print('El área es:',area)
lista.append(area)
resultado = sorted(lista, reverse=True)[0:-1]
#print("Áreas encontradas:", resultado)
Parte dos del código:
(x, y, w, h) = cv2.boundingRect(c)
rect= cv2.rectangle(th1, (x, y), (x + w, y + h), (0, 255, 0), 1, cv2.LINE_AA)
#Puntos para la parábola
#Coordenadas de pixeles
vlxU=x
vlxL=x+w
vlyU=y
vlyL=y+h
pnx=list(range(0, len(rect[0])))
pny=list(range(0, len(rect)))
lx = random.choices(pnx, k=70)
ly = random.choices(pny, k=70)
lista_blanca = zip(lx, ly)
lista_oscura = []
#print("Examinar muestra aleatoria")
for x, y in lista_blanca:
pixel = rect[y][x]
#print(f"Pixel({x}, {y}) = {pixel}")
if pixel == 0:
lista_oscura.append((x,y))
#print("Listar lista oscura")
lstx = []
lsty = []
for x, y in lista_oscura:
pixel = rect[y][x]
#print(f"Pixel({x}, {y}) = {pixel}")
#print("Pixel{x} de oscura",[x])
#print("Pixel{y} de oscura",[y])
lstx.append(x)
lsty.append(y)
if len(lstx) and len(lsty) == 30:
print('valores 30 de x ',lstx)
print('valores 30 de y ',lsty)
model = np.poly1d(np.polyfit(listax,listay, 2))
polyline = np.linspace(150, 450, 450)
plt.plot(listax,listay,'+r')
plt.plot(polyline, model(polyline))
plt.axis('off')
plt.legend()
plt.imshow(th1,cmap = 'gray')
plt.title("Funciona")
plt.figure()
cv2.waitKey(0)
cv2.destroyAllWindows()
En la imagen se muestra como se están tomando los puntos('+') en todas las zonas negras e ignora completamente el hecho de que en la parte dos del código se le dice que sea a lo largo de rect para que de ahí es donde se tomen los puntos.