如何在曲线上均匀地重新分布点

Question

我有3个维度的任意曲线，由XYZ笛卡尔点列表组成。点数不均匀分布（时间因素）。我怎样才能重建＆＃39;具有给定数量的点的曲线应该构成曲线。我在3D建模程序中看到了这一点，所以我很确定它可能，我只是不知道如何。

基于答案，我在python中需要它，所以我开始努力将interparc转换为python。我得到了线性插值。它可能效率低下并且有冗余，但可能对某人http://pastebin.com/L9NFvJyA

有用

Answer 1

我使用interparc，我的工具就是为了做到这一点。它通过2维或更多维的一般空间曲线拟合样条曲线，然后选择沿该曲线的距离相等的点。在三次样条函数的情况下，解决方案使用一个odesolver来做必须是数值积分的事情，因此它有点慢，但它仍然相当快。在许多情况下，简单的线性插值（我在这里使用）将是完全足够的，并且非常快。

曲线可能是完全一般的，甚至可以跨越自身。我将给出一个三维空间曲线的简单示例：

t = linspace(0,1,500).^3;
x = sin(2*pi*t);
y = sin(pi*t);
z = cos(3*x + y);
plot3(x,y,z,'o')
grid on
box on
view(-9,12)

xyzi = interparc(100,x,y,z,'lin');
plot3(xyzi(:,1),xyzi(:,2),xyzi(:,3),'o')
box on
grid on
view(-9,12)

Answer 2

你的“曲线”是一串连接一堆点的线段。每个线段都有一个长度;曲线的总长度是这些线段长度的总和。

因此，计算d = totalCurveLength / (numberOfPoints - 1)，并将曲线拆分为(numberOfPoints - 1)长度为d的块。

Answer 3

首先，感谢John D'Errico先生的interparc。真是太棒了！

我也遇到了这个问题，但我不熟悉MATLAB引擎API。鉴于此，我尝试将部分interparc Matlab代码转换为Python（仅包括线性插值，因为它足以解决我的问题）。

所以这是我的代码;希望它可以帮助所有pythonics寻求类似的东西：

import numpy as np

def interpcurve(N,pX,pY):
#equally spaced in arclength
N=np.transpose(np.linspace(0,1,N))

#how many points will be uniformly interpolated?
nt=N.size

#number of points on the curve
n=pX.size
pxy=np.array((pX,pY)).T
p1=pxy[0,:]
pend=pxy[-1,:]
last_segment= np.linalg.norm(np.subtract(p1,pend))
epsilon= 10*np.finfo(float).eps

#IF the two end points are not close enough lets close the curve
if last_segment > epsilon*np.linalg.norm(np.amax(abs(pxy),axis=0)):
    pxy=np.vstack((pxy,p1))
    nt = nt + 1
else:
    print('Contour already closed')

pt=np.zeros((nt,2))

#Compute the chordal arclength of each segment.
chordlen = (np.sum(np.diff(pxy,axis=0)**2,axis=1))**(1/2)
#Normalize the arclengths to a unit total
chordlen = chordlen/np.sum(chordlen)
#cumulative arclength
cumarc = np.append(0,np.cumsum(chordlen))

tbins= np.digitize(N,cumarc) # bin index in which each N is in

#catch any problems at the ends
tbins[np.where(tbins<=0 | (N<=0))]=1
tbins[np.where(tbins >= n | (N >= 1))] = n - 1      

s = np.divide((N - cumarc[tbins]),chordlen[tbins-1])
pt = pxy[tbins,:] + np.multiply((pxy[tbins,:] - pxy[tbins-1,:]),(np.vstack([s]*2)).T)

return pt 

Answer 4

不确定我是否关注，但不是存储实际数据，而是可以存储点到点的增量，然后从增量重建曲线，这样就不会出现任何空白点？但是，这会改变曲线的形状。