scipy.interpolate.UnivariateSpline не сглаживается независимо от параметров

Короткий ответ: вам нужно выбрать значение для s более тщательно.

В документации для UnivariateSpline указано, что:

 Positive smoothing factor used to choose the number of knots. Number of knots will be increased until the smoothing condition is satisfied: sum((w[i]*(y[i]-s(x[i])))**2,axis=0) <= s 

Из этого можно вывести, что «разумные» значения для сглаживания, если вы не проходите в явных весах, находятся вокруг s = m * v где m – количество точек данных и v – дисперсия данных. В этом случае s_good ~ 5e7 .

EDIT : разумные значения для s зависят, конечно, от уровня шума в данных. Документы, похоже, рекомендуют выбирать s в диапазоне (m - sqrt(2*m)) * std**2 <= s <= (m + sqrt(2*m)) * std**2 где std является стандартным отклонение, связанное с «шумом», который вы хотите сгладить.

@ Ответ Женя на ручные настройки узлов между точками данных был слишком грубым, чтобы обеспечить хорошие результаты в шумных данных, не будучи выборочным о том, как применяется этот метод. Однако, воодушевленный его предложением, у меня был успех с кластеризацией Mean-Shift из пакета scikit-learn. Он выполняет автоматическое определение количества кластеров и, как представляется, выполняет довольно хорошую работу сглаживания (на самом деле очень гладкая).

 # Imports import numpy import pylab import scipy import sklearn.cluster # Set up original data - note that it's monotonically increasing by X value! data = {} data['original'] = {} data['original']['x'] = [0, 5024.2059124920379, 7933.1645067836089, 7990.4664106277542, 9879.9717114947653, 13738.60563208926, 15113.277958924193] data['original']['y'] = [0.0, 3072.5653360000988, 5477.2689107965398, 5851.6866463790966, 6056.3852496014106, 7895.2332350173638, 9154.2956175610598] # Cluster data, sort it and and save inputNumpy = numpy.array([[data['original']['x'][i], data['original']['y'][i]] for i in range(0, len(data['original']['x']))]) meanShift = sklearn.cluster.MeanShift() meanShift.fit(inputNumpy) clusteredData = [[pair[0], pair[1]] for pair in meanShift.cluster_centers_] clusteredData.sort(lambda pair1, pair2: cmp(pair1[0],pair2[0])) data['clustered'] = {} data['clustered']['x'] = [pair[0] for pair in clusteredData] data['clustered']['y'] = [pair[1] for pair in clusteredData] # Build a spline using the clustered data and predict mySpline = scipy.interpolate.UnivariateSpline(x=data['clustered']['x'], y=data['clustered']['y'], k=1) xi = range(0, round(max(data['original']['x']), -3) + 3000, 20) yi = mySpline(xi) # Plot the datapoints pylab.plot(data['clustered']['x'], data['clustered']['y'], "D", label="Datapoints (%s)" % 'clustered') pylab.plot(xi, yi, label="Predicted (%s)" % 'clustered') pylab.plot(data['original']['x'], data['original']['y'], "o", label="Datapoints (%s)" % 'original') # Show the plot pylab.grid(True) pylab.xticks(rotation=45) pylab.legend( loc="lower right" ) pylab.show() 

Хотя я не знаю какой-либо библиотеки, которая будет делать это для вас, я бы попробовал немного более DIY-подход: я бы начал с создания сплайна с узлами между необработанными точками данных, как в x и в y . В вашем конкретном примере, имея единственный узел между 4-м и 5-м точками, должен сделать трюк, так как он удалит огромную производную около x=8000 .