scipy: इंटरपोलेटिंग प्रक्षेपण

Question

मेरे पास (x, y) जोड़े के अनुक्रम द्वारा गठित एक प्रक्षेपवक्र है। मैं splines का उपयोग कर इस प्रक्षेपवक्र पर अंक interpolate करना चाहता हूँ।

मैं यह कैसे कर सकता हूं? scipy.interpolate.UnivariateSpline का उपयोग करना काम नहीं करता है क्योंकि न तो x और न ही y मोनोटोनिक हैं। मैं प्रक्षेपण के साथ एक parametrization (उदाहरण के लिए लंबाई डी) पेश कर सकते हैं, लेकिन फिर मेरे पास दो निर्भर चर x (डी) और वाई (डी) है।

उदाहरण:

import numpy as np 
import matplotlib.pyplot as plt 
import math 

error = 0.1 
x0 = 1 
y0 = 1 
r0 = 0.5 

alpha = np.linspace(0, 2*math.pi, 40, endpoint=False) 
r = r0 + error * np.random.random(len(alpha)) 
x = x0 + r * np.cos(alpha) 
y = x0 + r * np.sin(alpha) 
plt.scatter(x, y, color='blue', label='given') 

# For this special case, the following code produces the 
# desired results. However, I need something that depends 
# only on x and y: 
from scipy.interpolate import interp1d 
alpha_i = np.linspace(alpha[0], alpha[-1], 100) 
r_i = interp1d(alpha, r, kind=3)(alpha_i) 
x_i = x0 + r_i * np.cos(alpha_i) 
y_i = x0 + r_i * np.sin(alpha_i) 
plt.plot(x_i, y_i, color='green', label='desired') 

plt.legend() 
plt.show() 

Answer 1

splprep का उपयोग करके आप किसी भी ज्यामिति के घटता से अधिक को जोड़ सकते हैं।

from scipy import interpolate 
tck,u=interpolate.splprep([x,y],s=0.0) 
x_i,y_i= interpolate.splev(np.linspace(0,1,100),tck) 

कौन सा दिया की तरह एक साजिश पैदा करता है, लेकिन केवल x और y अंक और नहीं अल्फा और आर पैरामीटर का उपयोग कर।

मेरे मूल उत्तर के बारे में क्षमा करें, मैंने सवाल को गलत तरीके से पढ़ा।