概述
离散点求曲率
第一步:离散点拟合函数
第二步:求一阶及二阶导,并带入求曲率公式。
import numpy as np
# x = 离散点x-axis坐标
# y = 离散点y-axis坐标
# x_target = 为拟合的函数采样点(这里是采样点的x-axis值)
# 第一步拟合函数
poly = np.polyfit(x, y, 3) # 拟合n次函数(这里为3次)
poly = np.poly1d(poly)
y_target = poly(x_target)
# 第二步求导
x_t = np.gradient(x_target) # x' 一阶导
xx_t = np.gradient(x_t) # x'' 二阶导
y_t = np.gradient(y_target) # y' 一阶导
yy_t = np.gradient(y_t)
# y'' 二阶导
curvature_val = np.abs(xx_t * y_t - x_t * yy_t) / (x_t * x_t + y_t * y_t) ** 1.5 # 求曲率公式
最后
以上就是可耐小白菜为你收集整理的离散点求曲率的全部内容,希望文章能够帮你解决离散点求曲率所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复