Estoy intentando entender un código que no he escrito yo.
En la definición de la función calc_R
aparecen cuatro variables:x, y, xc, yc
. Sin embargo, al llamar a esta función tan sólo le pasan dos variables y funciona correctamente.
¿Por que?
A continuación muestro el código:
# -*- coding: utf-8 -*-
# == METHOD 2 ==
# Basic usage of optimize.leastsq
from scipy import optimize as opt
import numpy as np
import f02_graphs as grp
#method_2 = "leastsq"
def calc_R(xc, yc):
""" Calculate the distance of each 2D points from the center (xc, yc) """
return np.sqrt((x-xc)**2 + (y-yc)**2)
def f_2(c):
""" Calculate the algebraic distance between the 2D points and the mean circle centered at c=(xc, yc) """
Ri = calc_R(*c)
return Ri - Ri.mean()
# Datos.
x = [36, 36, 19, 18, 33, 26]
y = [14, 10, 28, 31, 18, 26]
basename = 'arc'
# Coordenadas del baricentro.
x_m = np.mean(x)
y_m = np.mean(y)
barycenter = x_m, y_m
# Aplicamos un algoritmo para aproximar el centro de la circunferencia.
circle_center, ier = opt.leastsq(f_2, barycenter)
xc, yc = circle_center
# Calculamos la distancia al centro de cada par x,y.
R_xy = calc_R(xc, yc)
# Calculamos el promedio de todoas las distancias de R_xy.
R_2 = R_xy.mean()
global
delante. Mírate esta explicación