yo tengoun catálogo de datos y quiero usarlo en miMCMC código. Lo que es crucial es la velocidad de implementación, para evitar frenar mi muestreo de Monte Carlo de la cadena Markov.El problema: En el catálogo, tengo en la primera y segunda columna dos parámetros llamadosra ydec cuales soncoordenadas del cielo:

data=np.loadtxt('Final.Cluster.Shear.NegligibleShotNoise.Redshift.cat')
ra=data[:,0]
dec=data[:,1]

luego en las siete y ocho columnasX yY posiciones, es decir, las coordenadas de la cuadrícula, son puntos en un espacio de cuadrícula

Xpos=data[:,6]
Ypos=data[:,7]

En la función que he escrito y es necesario llamarla como un millón de veces, le daré unaXcenter yYcenter posiciones (por ejemplo, Xcenter = 200.6, Ycenter = 310.9) como entradas a la función y quiero encontrar los puntos de correspondencia en elra ydec columnas Sin embargo, puede ocurrir que las entradas no tengan ninguna correspondencia real en elra ydec. Así que quiero hacer una interpolación en caso de que no haya entradas similares paraX yY yra ydec datos en el catálogo y obtener las coordenadas interpoladas basadas en realra ydec entradas en el catálogo.