Tengo una gran nd.ndarray 3d de datos que representa una variable física muestreada sobre un volumen en una cuadrícula regular (como en el valor en array [0,0,0] representa el valor en coords físico (0,0,0 )).

Me gustaría ir a un espacio de cuadrícula más fino interpolando los datos en la cuadrícula en bruto. En este momento estoy usando la interpolación lineal de scipy griddata pero es bastante lenta (~ 90secs para una matriz de 20x20x20). Está un poco sobre desarrollado para mis propósitos, lo que permite el muestreo aleatorio de los datos de volumen. ¿Hay algo por ahí que pueda aprovechar mis datos espaciados regularmente y el hecho de que solo hay un conjunto limitado de puntos específicos a los que quiero interpolar?