Numpy 什么';s由scipy.interpolate.bisplev给出的B样条偏导数的含义是什么?
我有一个从cloudpoint生成的网格,可以描述为Numpy 什么';s由scipy.interpolate.bisplev给出的B样条偏导数的含义是什么?,numpy,scipy,interpolation,scientific-computing,Numpy,Scipy,Interpolation,Scientific Computing,我有一个从cloudpoint生成的网格,可以描述为z=f(x,y),所以我使用scipy.interpolate.bisplrep和bisplev,效果很好 bisplev可与参数dx=n和/或dy=n一起使用,因此结果在评估点处是顺序n的导数。我计划用它来计算平均曲率和高斯曲率(在Matlab中称为surfature),这应该包括得到survace的二阶偏导数 每次使用一个偏导数(例如dx)的结果非常好,可以清楚地将梯度表示为“着色”效果,如这张来自人类背部的图像所示(代码优先): 到目前
z=f(x,y)
,所以我使用scipy.interpolate.bisplrep
和bisplev
,效果很好
bisplev
可与参数dx=n
和/或dy=n
一起使用,因此结果在评估点处是顺序n
的导数。我计划用它来计算平均曲率和高斯曲率(在Matlab中称为surfature
),这应该包括得到survace的二阶偏导数
每次使用一个偏导数(例如dx
)的结果非常好,可以清楚地将梯度表示为“着色”效果,如这张来自人类背部的图像所示(代码优先):
到目前为止,一切顺利。问题是:当我同时要求两个偏导数时,我无法理解(也无法找到任何解释)结果的含义,因为没有任何明显的数字或视觉意义。例如,如果我使用dx
和dy
:
self.mesh_shadow = inter.bisplev(yy.flat, xx.flat, self.spline, dx=1, dy=1)
我明白了:
所以,我想知道:
(…,dx=1,dy=2)
,函数似乎产生了“有效”的结果,但这有意义吗我必须研究你的观点,它们看起来很有价值,正是我所需要的。d2不连续性的优点,因为我使用的是三阶样条曲线,而且我已经读到它们保证在几阶(二阶?)之前是可微的。谢谢
self.mesh_shadow = inter.bisplev(yy.flat, xx.flat, self.spline, dx=1, dy=1)