我正在将一个用R编写的脚本移植到Python上。在R中，我使用smooth.spline，在Python中，我使用SciPy UnivariateSpline。它们不会产生相同的结果(即使它们都基于三次样条线方法)。有没有一种方法或替代UnivariateSpline的方法，使Python样条返回与R相同的样条线？
我是个数学家。我理解样条的一般概念。但没有详细介绍它们在Python或R中的实现。
下面是用R编写的代码，然后是Python。两者的输入数据是相同的。
以下是输入数据：
x = 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8
浏览 1提问于2019-06-19得票数 2
我有一个问题，只有在使用stdin输入从控制台(或ProcessBuilder )运行python脚本时，才会出现在Windows 10上。
当我在Ubuntu上运行这段代码时：
import sys
print(sys.stdin.read())
使用此命令：
python3 Execute.py < json.txt
我在我的控制台上获得了json.txt的内容："Hello“
当我在Windows上使用相同的代码运行相同的命令(没有"<")时，进程的反应就像一个无限循环。(我的开发环境有什么问题？)
在这两种环境中，我的python版本都是3.8.10
浏览 1提问于2022-06-01得票数 0
回答已采纳
我有一些数据点，我喜欢用三次B样条拟合来近似。在另一个程序中，我喜欢只使用拟合的结和系数来插值一些点。这意味着我需要编写代码来计算给定节点和系数的点。使用，我已经能够正确地评估各个点，其结果是
knots,coeff,n=scipy.interpolate.splrep(x,y)
奇怪的是，我需要删除数组中的第一个和最后三个条目才能让代码正常工作(注意:前4个和最后4个节条目是相同的，因此会产生0除法错误)。有人知道为什么会有这些重复的条目吗？除了删除这些条目外，wiki页面上的等式也能工作。从结果来看
knots=[ 15., 15., 15., 15., 75.
浏览 3提问于2014-05-07得票数 1
2回答
我有一堆坐标，它们是2D平面上固定的均匀三次B样条曲线的控制点。我想使用Cairo调用(在Python中，使用Cairo的Python绑定)来绘制这条曲线，但据我所知，Cairo只支持贝塞尔曲线。我也知道两个控制点之间的B样条曲线段可以用Bézier曲线绘制，但我在任何地方都找不到确切的公式。给定控制点的坐标，我如何导出相应的贝塞尔曲线的控制点？有没有什么有效的算法来解决这个问题？
浏览 2提问于2010-03-29得票数 3
回答已采纳
1回答
我正在使用numpy将c代码转换为python。
什么是numpy等价于下面的matlab函数？
1. mxGetNumberOfDimensions
2. mxIsDouble
3. mxGetM
4. mxCreateDoubleMatrix
5. mxGetPr
6. lookupspline
浏览 1提问于2014-01-21得票数 1
回答已采纳
我正在开发一个程序，其中一个要求是以DXF作为输入。输入仅限于2D情况。程序本身在C++/Qt中，但是要测试它，我需要一些示例DXF输入。样条导入已经实现，下一步是增加样条拟合点或控制点的折线。我决定使用Python/ezdxf生成这样的polyline，因为我没有Autocad。
我的第一种方法是利用add_spline_control_frame从fit点创建样条，然后将其转换为polyline。问题是没有从样条转换到polyline (虽然我想我在文档中看到了它，但是找不到它)。
目前的方法是用add_polyline2d(points)制作多边形，使每个点与DXF标志字段相等，8 (样
浏览 7提问于2020-09-01得票数 1
回答已采纳
我在使用python工具"statsmodels.api.GLM“时遇到了一个问题，我无法理解。我来这里寻求帮助。 我正在做这个页面https://www.analyticsvidhya.com/blog/2018/03/introduction-regression-splines-python-codes/ (页面或here中包含数据链接)上的“三次和自然三次样条”的示例(请参阅)。 问题是。在对数据进行拟合后，我尝试预测x的给定位置的值(例如，以下代码中的xp00和xp01 )。然后我发现，一旦请求的位置与训练x集(即xp)具有不同的min和max (即xp01)，结果就变成了
浏览 85提问于2021-06-24得票数 0
1回答
我很难将样条插值的数学概念与样条滤波器在python中的应用联系起来。我对样条插值的基本理解是，它以分段的方式拟合数据，而拟合的分段多项式称为样条。但是它在图像处理中的应用涉及到对图像进行预处理，然后进行插值，这一点我很难理解。
为了给出一个例子，我想使用scipy.ndimage.map_coordinates(input, coordinates, prefilter=True)和关键字prefilter根据文档对图像进行插值：
确定输入数组在插值前是否使用spline_filter预先过滤。
scipy.ndimage.interpolation.spline_filter的文档简单地
浏览 2提问于2019-10-08得票数 2
4回答
大数据集，插值
、、、
我正在寻找一种“方法”，以获得一个公式，公式来自拟合一组数据(3000点)。我使用的是Legendre多项式，但对于> 20的点，它给出的值并不准确。我可以写chi2测试，但是算法需要大量的时间来计算N个参数，而且一开始我不知道函数是什么样子的，所以它需要时间。我在想样条线。也许..。
所以输入是: 3000品脱
输出: f(x) = ...某物
我想要一个来自fit的公式。在python中做这件事的最好方法是什么？
让原力与我们同在！尼康
浏览 0提问于2011-02-17得票数 1
回答已采纳
我试图训练一个预测模型，并希望使用样条样条插值来表示一些函数，构成模型的一部分。然而，这并不是一个简单的情况，在某个区域上拟合某些x，y数据以求样条系数，而是用样条近似的函数构成一个非线性生成模型的一部分。为了找到系数，我需要使用非线性最小化算法来对训练数据集进行优化。这意味着我需要能够直接指定一组系数，而不是使用scipy.interpolate中的拟合方法(如scipy.interpolate.UnivariateSpline)。
是否有一种方法可以指定样条系数，然后将得到的对象作为模型中的函数使用？如果这在this中是不可能的，那么是否还有另一个支持此功能的python库？
浏览 0提问于2014-01-15得票数 0
回答已采纳
我正在尝试使用Python的SciPy包通过一些数据生成样条。关于的问题，我认为在使用interpolate.splrep的开始和结束时存在重复结的原因是因为它正在制造一个夹紧样条。我想让这些结通过数据(s=0)。分割器的输出重复第一个结和最后一个结各4倍(如上面的问题)，它包括内部数据点作为节，除了第二点和下一点到最后一点。这是一个MWE：
import math
import numpy as np
import scipy.interpolate
import matplotlib.pyplot as plt
x_dec = np.array([ 3.0289979374768700e
浏览 5提问于2015-12-28得票数 0
回答已采纳
我正在使用Python3.6，并尝试运行来自的以下代码
from __future__ import division
import numpy as np
import matplotlib.pyplot as plt
import scipy.interpolate
x = np.array([ 2., 1., 1., 2., 2., 4., 4., 3.])
y = np.array([ 1., 2., 3., 4., 2., 3., 2., 1.])
plt.plot(x,y, label='poly')
t = np.arange(x.
浏览 66提问于2020-03-04得票数 0
1回答
我正在使用python读取netcdf文件，并需要使用python中的matplotlib库绘制图表。netcdf文件包含3个变量：u、v和w组件。我必须在垂直的尺度上画出这三个组件。由于这些数据将被用于天气预报雷达，所以我需要在拐角处用光滑的曲线来绘制这些图表。现在情节看上去
代码是.
from netCDF4 import Dataset
from pylab import *
import numpy
from scipy import interpolate
from scipy.interpolate import spline
#passing the filename
浏览 2提问于2013-02-05得票数 13
回答已采纳
2回答
给定一组2D点，我想用给定数量的内部节点来拟合这个数据的最优样条。
我已经看到，我们可以使用来指定结的数目和位置，但是它不允许我们只指定结的数目。
从文档来看，这似乎意味着它们有一种方法可以用给定的节数来拟合样条，因为平滑因子s的文档状态(强调s)：
正平滑因子用于选择节点数。在满足平滑条件之前，结数会增加。
因此，虽然我可以用一种反向的方式，通过平滑因子搜索，直到它产生一个样条，与所需的节数，这似乎是一个相当荒谬的方式，从计算效率的角度来处理。两个额外的搜索步骤正在发生，仅仅是为了相互取消并获得一个已经在开始时直接计算的结果。
我已经搜索过了，但是还没有找到一个函数来直接访问这个给定
浏览 0提问于2019-12-12得票数 1
在我的电火花代码im读取测试csv文件，过滤它，并编写。在控制台中，我可以看到以LoggingLineageDispatcher格式的json格式的所有操作，但我希望找到一种方法，在我的python代码中正确地获取这些数据。找不到任何选择。
我的pyspark代码：
session = create_spark_session()
test_df: DataFrame = session.read.csv(
"test.csv",
sep =',',
header = True
)
mc_df = test_df.filter(col
浏览 8提问于2022-05-30得票数 0
1回答
我目前正在进行一个项目，使用opencv和python来测量通常是弯曲的对象，例如，下面所示的箭头尽可能准确。
我认为一种策略可能是使用枕叶Voronoi函数来获得箭头中心脊柱的点，但现在有麻烦了。这是我的代码：
img = cv2.imread('example_rubystreak_2.PNG')
img.shape
gray = cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
ret, thresh = cv2.threshold(gray,190,255,cv2.THRESH_BINARY)
countimage, contours,
浏览 10提问于2015-08-07得票数 2
回答已采纳
1回答
我的日期是对象格式的，因此我很难将它转换为浮点数或整数。
我在python中使用样条函数得到错误。
如何从样条图中提取数学方程？谢谢
浏览 4提问于2022-06-14得票数 2
在Matlab中，利用样条插值函数spapi生成样条曲线，使曲线及其前两个导数的给定值相匹配。例如：
spapi([0 0 0 0 1 2 2 2 2],[0 1 1 1 2],[2 0 1 2 -1])
这就产生了区间0..2上唯一的三次样条f，在1处正好有一个内结，满足以下五个条件
f(0)=2，f(1)=0，f'(1)=1，f''(1)=2，f(2)=-1
它们包括1处的3倍匹配，即与函数及其前两个导数的指定值匹配。
然而，在Python中，样条插值函数scipy.interpolate.splrep(x,y)不能工作，因为x应该是单调递增的。
那么，我怎
浏览 3提问于2018-09-02得票数 1
我试图求多项式的曲率。X和Y分别是X和Y坐标的python列表。我使用scipy.interpolate是因为我能够在我的图像中看到更好的曲线。但是，一旦我找到一个2D多项式的系数并将它们重新绘制回图像，重新绘制的曲线看起来就太离谱了。
如何求多项式曲线的精确系数？
插值= interpolate.interp1d(X，Y) z= np.polyfit(X，插值(X)，2)#系数poly_y = [z*x*x + z1*x + z2 for x in X] #重新计算Y坐标plt.plot(X，poly_y)
浏览 3提问于2018-02-28得票数 0
我正在使用Python3中优秀的svgpathtools库来处理SVG文件中的一些路径，该文件是在一个矢量绘图应用程序中创建的。
我想为SVG中包含的每个路径创建详细的点数组，其中点在路径上是等距的。下面的代码就是这样做的，但是如果超过几千个样本，就会变得非常慢。
SAMPLES_PER_PX = 1
fname = "/path/to/file.svg"
paths, attributes = svg2paths(fname)
myPaths = {}
for path,attr in zip(paths, attributes):
myPathList = []
浏览 4提问于2017-08-22得票数 3
我试图通过定义b样条的阶数，基函数，节点和求值范围来生成一个基样条函数。请向我推荐一个合适的Python函数，它可以帮助我。 我当前的实现是使用johntfoster/b样条方法。它不允许我定义基函数的数量，并且结果与MATLAB的结果不相似。https://github.com/johntfoster/bspline scipy.interpolate.BSpline.basis_element函数不允许我定义样条的阶数、基函数、结点 Matlab实现： nbreaks = 20;
nbasis = nbreaks + norder - 2;
breaks = linspace(0,tau
浏览 39提问于2019-05-17得票数 2
回答已采纳
1回答
R中的单调样条
、、、
我试图用R中的来构建一个单调样条。
我在评估自变量的新值的模型时遇到了困难。一般来说，我很难掌握“预测”的R处理及其与spline2的关系，以及如何使用产生的bs对象。
我试着跟踪使用的spline1。我的数据在一个名为BRIyII的数据文件中，自变量t和因变量P因此：
plot(BRIyII$t,BRIyII$P)
产量：
所以我想：
knots=c(9)
myMat=mSpline(BRIyII$t, knots = knots, degree = 3, intercept = TRUE)
mylm=lm(BRIyII$P~myMat)
现在，如果我们：
pr = predict(m
浏览 4提问于2017-07-26得票数 0
回答已采纳
我使用对大量数据进行平滑插值。效果很好。我得到一个对象，它的行为类似于一个函数。
现在我想保存样条点以备以后使用，并在Matlab中使用它们(还有Python，但这不是那么紧急)，而不需要原始数据。我该怎么做呢？
在scipy中，我没有线索；UnivariateSpline似乎没有提供具有先前计算的结点和系数的构造函数。
在MATLAB中，我尝试了Matlab函数和pchip()，虽然这两个函数都很接近，但它们在端点附近都有类似于的错误。
以下是我拥有的Matlab格式的样本数据集：
splinedata = struct('coeffs',[-0.04127391809552
浏览 0提问于2013-06-05得票数 7
1回答
维基百科上的内容和图片表明B样条没有通过控制点，而我在其他地方找到的信息与维基百科上的内容相同。然而，样条插值可以在scipy.interpolate中进行，例如interp1d 'cubic'或splrep splev(这两个操作应该是类似的)。通过绘制，我发现曲线通过了所有控制点。但文档和库代码表明，这些操作的基础是B样条。
我现在对此深感困惑。我知道B样条是一种样条，除了B样条插值外，还有多种样条插值，如自然三次样条插值和三次Hermitian插值。
我现在的问题是：
样条曲线与样条插值的关系。样条插值总是通过所有控制点吗？
B样条和B样条插值和三次样条插值是一样
浏览 18提问于2020-11-22得票数 1
回答已采纳
1回答
我正在编写函数，使用scipy.interpolate函数在python中计算一维插值。利用文档的帮助，我为三次样条插值和三次样条插值编写了两个不同的函数。
# calculate cubic interpolation
def linear_interpolation(x):
linear = interpolate.interp1d(support_x, support_y, 'cubic')
return linear(x)
# calculate cubic spline interpolation
def cubic_spline_interpol
浏览 3提问于2017-10-09得票数 3
回答已采纳
2回答
我有一个数据集，其中我有三年的观测值(例如，2000、2005和2010)，并且需要使用R对其间年份的值进行插值。我尝试使用某种类型的样条来实现这一点，但是，插值的值超出了原始范围。在下面的情况下，它们甚至变成了负值。 years <- c(2000, 2005, 2010)
outcome_values <- c(1, 10, 90)
plot(spline(years, outcome_values, xout = seq(min(years), max(years))))
points(years, outcome_values, pch = 16) ? 有人用Py
浏览 20提问于2020-06-19得票数 0
1回答
我有一个catmull-rom样条线，有6个CV，我想沿着该样条线在任意点剪切它。
如何计算两条新样条线的端点，以使它们保持被一分为二之前的形状？
因此，假设这是具有6个点的样条线：
p0 p2 p4
\ / \ / \
\ / \C / \
\ / \ / \
p1 p3 p5
如果我想在点C处切割这条样条曲线，那么我将有两条新的样条曲线，我想每条都有5个CV？
Spline #1: p0, p1, p2, X, C
Spline #2: C, X, p3, p4, p5
其中X是为了保
浏览 2提问于2018-09-20得票数 0
1回答
我一直致力于让Catmull-Rom splines为一个辅助项目工作，但我很难让它做我需要的事情。我尝试了以下两个实现，但都不适用于我，并且我无法跟踪我的代码中相对于他们的代码的任何错误(我必须假设这些错误已经过测试)。我称他们的为"ABC“解决方案：
然后我实现了以下解决方案(我称之为“矩阵”解决方案)，并使用编辑后的版本3帖子进行了工作：
然而，此matrix解决方案仅实现了Catmull-Rom，并在矩阵中内置了0.5 'a‘值。我想让Chordal工作，因此我需要'a‘== 1。
考虑到我的ABC版本的解决方案造成了问题，我尝试使用这里的矩阵()来传入我自
浏览 1提问于2015-06-22得票数 0
1回答
我正在开发一个3D图形应用程序，用户可以在其中绘制曲线。我记录由用户绘制的曲线，我想从记录的一组点创建一个平滑的nurb。我试着使用openNurbs库，但是我找不到一种方法来使用这个库来进行装配。我怎样才能把一组点配给护士呢？
浏览 5提问于2014-12-19得票数 2
回答已采纳
AutoCAD允许将样条实体存储在仅由拟合点定义的DXF文件中，问题是样条定义具有无限的数值正确解，Autodesk不提供从给定的拟合点计算所需参数所需的信息。
tl;dr -缺失的信息是估计的开始和结束切线的方向和大小输入切线的全局B样条插值的末端导数，谁能帮助计算这个值吗？
我使用BricsCAD进行测试，但是"Trueview 2020“显示了相同的结果。
1.设想情况
只给出拟合点，使用不受任何约束的全局曲线插值，得到由控制顶点定义的样条：
# First spline defined by control vertices interpolated from given f
浏览 3提问于2020-06-19得票数 7
library(ISLR)
fit=lm(wage~bs(age ,knots =c(25 ,40 ,60),degree = 5,),data=Wage)
fit=lm(wage~ns(age ,knots =c(25 ,40 ,60),degree = 5,),data=Wage)
我能够构建5次多项式回归样条，但是如何构建5次自然样条，因为ns()函数缺少次数参数。
我只能用ns()生成一个三次自然样条线。还有没有其他函数可以用来生成二次自然样条，等等？
浏览 2提问于2018-11-16得票数 1
1回答
有没有可能在Python中找到基于De Boor方法的平滑样条？用于数据近似。
之前我在matlab中使用了平滑样条，我需要在python中使用完全相同的算法。
浏览 2提问于2017-04-06得票数 1
除了这一种情况之外，我使用了来自UnivariateSpline模块的data.It来适应几乎所有的情况，这会导致Process finished with exit code -1073741819 (0xC0000005)错误。如果我将平滑因子s改为0，它也能工作。任何解决这个问题的建议都会有帮助。
Update1
我的工作环境是：
python 3.7scipy 1.3.2numpy 1.17.4
import matplotlib.pyplot as plt
import numpy as np
from scipy.interpolate import UnivariateSpl
浏览 2提问于2019-11-13得票数 0
1回答
样本点是在固定采样率(e，g 1s=25帧)下从飞行器上提取的，它包含x，y，z位置。
其要求是导入点列表并编辑曲线，然后将其导出到新的点列表。我想使用三次Bezier曲线来显示点，问题是：
飞行轨迹是任意的，可以用贝塞尔曲线来画吗？
如何确定顶点和控制点？
顶点和控制点的原因是为了编辑，我想拖动控制点来修改曲线。
这是一个新的领域，任何建议或样例代码都将不胜感激。:)
浏览 1提问于2012-08-02得票数 3
回答已采纳
我正在使用monogame创建一个游戏，它接受用户的触摸输入，并沿着他们绘制的路径移动一个精灵。我希望能够从屏幕上的某一点(比如一个正方形)画一条线，然后在屏幕上的任何方向上绘制这条线，让它平滑(贝塞尔曲线)，然后让一个字符跟随这条线。
有人知道monogame是否支持xna 4.0中的曲线类吗？我真的不知道从哪里开始使用上面提到的用户触摸输入的技术来绘制和平滑线，因为我的数学不是很好。
以前有没有人这样做过，不介意给我指个正确的方向？任何您认为可能有用的片段或文章，请将其发布。我找到的所有文章都使用预先确定的路径，这些路径是使用beziers平滑的，而不是来自用户的触摸输入。
干杯
浏览 6提问于2012-08-19得票数 0
回答已采纳
我已经找到了像这样的图像的边界：
现在我需要填充这个区域的内部(它是一个乳房组织)。opencv或python中的其他模块(如skimage )中有没有什么函数可以让我做到这一点？比如MATLAB中的imfill。
浏览 2提问于2016-06-02得票数 2
我正在尝试将一个重复的对象映射到一个环形结的路径上，有人熟悉这个过程吗？
我知道如何使用“路径约束”和“快照”沿路径映射对象，但我似乎不知道如何使用像圆环结这样的基元来实现这一点。
浏览 2提问于2015-01-29得票数 0
我试图从GAM模型中提取节点的位置，以便将我的预测变量划分为另一个模型的类别。我的数据包含一个二进制响应变量(使用)和一个连续预测器(打开)。
data <- data.frame(Used = rep(c(1,0,0,0),1250),
Open = round(runif(5000,0,50), 0))
我适合这样的GAM：
mod <- gam(Used ~ s(Open), binomial, data = data)
我可以得到预测的值，以及模型矩阵等，在type=c("response", "lpmatrix
浏览 0提问于2016-05-24得票数 2
回答已采纳
3回答
我在三维空间中有N个点。我要用一条线加入他们。然而，如果我用一条简单的线条去做，它是不光滑的，而且看起来很难看。
我目前的方法是使用Bezier曲线，对4点使用DeCasteljau算法，并对数据集中的每组4点运行该曲线。然而，问题在于，由于我分别在1-4、5-8、9-12等点上运行，所以这条线在4-5、8-9等之间并不平坦。
我还寻找其他方法；具体而言，我找到了关于Catmull样条的，它似乎更适合我的目的，因为曲线通过所有控制点，不像Bezier曲线。所以我几乎开始实现它，但是后来，我在那个站点上看到这个公式是"assuming uniform spacing of control
浏览 4提问于2011-06-28得票数 3
回答已采纳
1回答
有人能解释在画布上画一条样条线的逻辑吗？我理解样条背后的基本几何逻辑，但我很难理解如何以编程方式实现这个逻辑。我看了几个例子(和)，但是没有一个在源代码中有足够的注释来确切地了解需要做什么。
我不想继续修改一个工作样条示例，因为我希望理解它背后的逻辑。
浏览 3提问于2013-09-13得票数 0
回答已采纳
1回答
我们有一个应用程序，需要绘制样条曲线(三次曲线、bezier曲线、b样条曲线等).我们目前在C中有工作算法，把这些曲线的控制点划成直线条。
我们要解决的问题是，需要不断地重新绘制基于缩放的曲线，以及曲线的可见程度。我猜这将是更实际和更有效率的GPU。
处理这个问题的最标准的方法是什么？我听说过几何着色器，但也听说它们效率低下，所以我不确定是否有更好的方法。如果不是，这样抚摸曲线的一般过程是什么？
一般情况下，冲刺算法采用曲线的控制点，确定一条线条的一组新的点集(基于从缩放水平计算出的公差，该公差定义了拍距)。
如果有什么不同的话，我们正在使用Vulkan。谢谢!
浏览 0提问于2021-10-16得票数 2
当我有非均匀的间隔输入时，我试图用Python获得一个很好的采样器。任何建议都会有帮助。我试过很多实习医生的功能。下面是一个例子：
from scipy.interpolate import InterpolatedUnivariateSpline
from numpy import linspace, arange, append
from matplotlib.pyplot import plot
F=[0, 1000,1500,2000,2500,3000,3500,4000,4500,5000,5500,22050]
M=[0.,2.85,2.49,1.65,1.55,1.81,1.3
浏览 14提问于2014-12-06得票数 2
回答已采纳
1回答
我正在尝试使用python自动将文本插入到dxf轮廓中。我有大量用于激光切割的dxf文件。我们经常想在钣金零件上刻上零件编号。我的尝试是制作一个正方形方框，其中的长度和with等于文本的高度和宽度。在编程之后，我找到了一个位于内部轮廓和外部轮廓之外的地方，我想用文本填充这个框。我尝试用多边形提取轮廓，并开始检查它是否适合，哪种类型的工作正常。还没有完全完成。
我想知道是否有某种库/tool有这个功能，因为在把更多的钱投入到程序中之前，计算时间是相当高的，或者任何人有比我更容易的方法。
浏览 4提问于2021-03-09得票数 0
我有三个3D矩阵X、Y和Z，它们在一些3D空间上定义了相同大小的矩阵V。矩阵的间距是规则的。现在，我正在尝试执行插值，并计算V的空间偏导数，即对于每个像素，计算V如何随x，y和z变化。我读到过，用样条插值和计算导数会得到很好的结果。例如，我以前使用过splinefit和ppdiff ()
如何对我拥有的数据集使用样条线？有没有更好的MATLAB代码(Python和C语言也可以)来执行这些类型的计算？
假设我只想要由X、Y和Z定义的采样位置的导数，我可以对每个维度进行一维样条近似，并以这种方式计算偏导数吗？也许这应该是数学交换的一个问题。这可能需要一段时间，但它应该可以工作，对吗？
谢谢你的帮忙
浏览 1提问于2013-01-29得票数 0
3回答
假设我在python中有两个数组，我希望得到(并实际使用)这两个点之间的三次样条插值。(即:我希望集成该函数)。我更喜欢使用numpy scipy的方式。
我知道的事。然而，这只允许我评估这些点，比如说非常简单的函数：
现在我可以做一些简单的事情，比如：
import numpy as np
import scipy.interpolate
import matplotlib.pyplot as plt
y = np.array([0,2,3,4,8,10,12,12,12,10,9,8,7,1,0,0,1,2])
x = np.array(range(len(y)))
xvals = np
浏览 2提问于2015-04-05得票数 1
3回答
一点背景。我有一个模拟，使用立方样条的一维轨迹。在这种情况下，三次样条将对象的位置、速度、加速度和加速度作为时间的函数指定。
如果你有：
位置、速度、加速度和时间的初始值和最终值
关于最大和最小速度、加速度和挺举的常数值约束
然后是一个独特的样条。如果你没有指定最后的时间，而是想要最小的时间轨迹，那么还有一个唯一的样条。
实际上，找到这些样条可能是一种巨大的痛苦。在指定时间的情况下，样条将包含最多7个多项式，并且节点(多项式之间的过渡点)不知道提前。
这不是将样条拟合到一组数据的通常情况，它是从边界条件和一些附加约束中创建样条。我读过一些文件，人们使用过类似的安排，也有类似的需
浏览 1提问于2010-12-06得票数 13
1回答
所以我有一个由N个数组组成的大集合，每个数组的长度为N_pix。我想在Python语言中的所有数组之间进行插值(即执行N_pix插值，每个数组有N个数据点)。
问题是每个数组都非常大(每个数组可以高达100sMB，其中有100个数组)，所以我希望避免将它们全部读取到内存中。
我最初的想法是使用生成器表达式逐个读取大文件，并将生成器传递给插值例程：
def read_data_from_disc(paths):
# Cycle through paths, and yield each large array
for path in paths:
yield
浏览 0提问于2017-10-23得票数 0
1回答
我想用B样条来计算numpy中闭曲线的曲率。为了得到一个平滑的结果，我想对样条表示上的导数而不是数据进行评估。但是，下面的代码返回一个错误：
Traceback (most recent call last):
File "<stdin>", line 16, in <module>
File "/Users/jfl/anaconda/lib/python2.7/site-packages/scipy/interpolate/fitpack.py", line 1212, in splder
浏览 0提问于2017-07-14得票数 2
回答已采纳
我有一组离散的点(x_n，y_n)，我想将它们近似/表示为B样条基函数的线性组合。我需要能够手动更改该方法使用的B样条基函数的数量，我正在尝试使用scipy在python中实现这一点。具体地说，下面是我使用的一些代码：
import scipy
spl = scipy.interpolate.splrep(x, y)
但是，除非我误解或遗漏了文档中的某些内容，否则我似乎无法更改scipy使用的B样条基函数的数量。这似乎是由x和y的大小设置的。所以，我的具体问题是：
我可以在上面使用的"splrep“函数中更改scipy使用的B-spline基函数吗？一旦我执行了上面代码中所示的变换，我
浏览 252提问于2018-07-06得票数 0
2回答
如何确保样条通过点
、、、
我试图用样条来计算图像一部分的导数。问题是，我似乎不知道如何使样条经过一个特定的点。我使用scipy.interpolate库和UnivariateSpline方法。
浏览 1提问于2015-07-13得票数 1
回答已采纳
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频
腾讯云开发者
