scipy BSpline Einbau in python

Dies ist mein erstes mal mit BSpline, und ich will fit einer Kurve auf meine Daten Punkte. Ich habe versucht, mit Hilfe von Univariaten Spline und versucht, auf splev und splrep aber ich würde wirklich gerne lernen, wie dies zu tun mit BSpline.

Es sieht aus wie mein fitting ist wirklich abgehackt und die Linie ist nicht selbst gehen durch die Punkte. 

arraymagU = linspace(U_timeband.min(),U_timeband.max(),300) #array for my x data points
UfunctionBS = BSpline(U_timeband,U_magband,k=4,extrapolate=False)
arraymagU2 = UfunctionBS(arraymagU)

plt.plot(arraymagU,arraymagU2)

U_timeband ist mein x-Koordinaten und U_magband ist einfach mein y. k=4 ich denke, dass zeigt eine kubische passen? Ich habe gespielt, um mit diesem Wert und es scheint nicht, um es besser zu machen.

Es produziert dies:

scipy BSpline Einbau in python

Wie kann ich das verbessern, konsistent?
Ich glaube, ich kann haben, zu definieren, Haltepunkte, aber ich bin mir nicht sicher, wie zu tun, dass entweder.

InformationsquelleAutor Sophia Medallon | 2017-07-19
  1. 7

    splrep gibt ein Tupel (t,c,k) mit der vector der Knoten der B-spline-Koeffizienten und den Grad des spline. Diese können gefüttert werden, um interpolate.BSpline zum erstellen einer BSpline-Objekt:

    import numpy as np
import scipy.interpolate as interpolate
import matplotlib.pyplot as plt

x = np.array([ 0. ,  1.2,  1.9,  3.2,  4. ,  6.5])
y = np.array([ 0. ,  2.3,  3. ,  4.3,  2.9,  3.1])

t, c, k = interpolate.splrep(x, y, s=0, k=4)
print('''\
t: {}
c: {}
k: {}
'''.format(t, c, k))
N = 100
xmin, xmax = x.min(), x.max()
xx = np.linspace(xmin, xmax, N)
spline = interpolate.BSpline(t, c, k, extrapolate=False)

plt.plot(x, y, 'bo', label='Original points')
plt.plot(xx, spline(xx), 'r', label='BSpline')
plt.grid()
plt.legend(loc='best')
plt.show()

    scipy BSpline Einbau in python

    • Danke. Ich habe UnivariateSpline bereits und führte es in eine bessere Passform, aber mein Berater ist, glaube ich, das beharren auf die Verwendung BSpline für Licht-Kurven (dies ist Photometrie-Daten von einer supernova). Wie konnte ich wissen, das best-fit-spline dann? Mit der spline aus UnivariateSpline oder vielleicht ein interp1d?
    InformationsquelleAutor unutbu
Schreibe einen Kommentar