lspia's Introduction

LSPIA

对LSPIA方法对实现。

该方法详细论述请看论文《Progressive and iterative approximation for least squares B-spline curve and surface fitting》.

注意事项

使用LSPIA拟合曲面及曲线时，需要数据点数量足够大。在初始测试时，我使用10个点拟合曲线，无法得到理想效果；使用20x20个点拟合曲面，误差十分大；
关于曲面拟合，论文中没有给出miu的定义，miu影响曲面拟合的速度，可自行根据数据进行调整；
初始点选择，对于初始点定义，可以按照文章中的方式间隔采样取点，也可以随机取点。我在测试时将所有数据点设为相同，仍能得到理想结果。初始点的定义只影响迭代的次数。

实验结果

LSIPIA拟合曲面

LSPIA拟合曲面

lspia's People

Contributors

Stargazers

Watchers

lspia's Issues

How to change number of control points P_h / P_l?

How to set number of control points P_h ? When I change to P_h = int(row / 20), it'll get error:

C:\Users\junew\Miniconda3\envs\tf_1.14.0\lib\site-packages\numpy\core\fromnumeric.py:86: RuntimeWarning: overflow encountered in reduce
  return ufunc.reduce(obj, axis, dtype, out, **passkwargs)
iteration  72
C:\Users\junew\Desktop\LSPIA-master\LSPIA-master\surface_fitting_error.py:24: RuntimeWarning: overflow encountered in square
  error = error + np.sum(np.square(D_dim - np.dot(np.dot(Nik_u, P_dim), np.transpose(Nik_v))))
LSPIA.py:295: RuntimeWarning: invalid value encountered in double_scalars
  while (abs(e[-1] - e[-2]) >= 1e-7):
inf

Also, what if the input data isn't gridded?

Recommend Projects