查找数字数组曲线的折点

Question

我有一个图像，我沿着图像对角线提取一条轮廓线，轮廓在对角线的2侧集成了几条线，所以我得到了一条平滑的轮廓线，现在我想找出曲线上表示图像线轮廓的拐点。我尝试过这种方法，但它给了我很多点，而且位置不是我想要的。import numpy as np import matplotlib.pyplot as plt arr = profiles[:,0].astype("int16") d = np.diff(arr) dd = np.diff(d) plt.figure(figsize=(20,15)) plt.plot(np.arange(2047), d) plt.plot(np.arange(2046), dd) plt.plot(arr)for i in range(2046):     if (d[i] != 0) and (d[i] == dd[i]):         plt.plot(i, arr[i], 'o')您知道如何检索曲线上拐点的坐标吗？

Answer 1

据此，您不是在寻找拐点。无论如何，这里有一个代码使用第一个和第二个派生，但在它有一个平滑的via savgol_filter（来自scipy）之前。我认为这是您问题的最佳近似值。

from scipy import signal
smooth = signal.savgol_filter(arr, 55,3)
dxdy =  signal.savgol_filter(np.gradient(smooth),55,3)
dx2dy2 = signal.savgol_filter(np.gradient(dxdy),55,3)
plt.subplot(121)
plt.plot(smooth,'k')
plt.plot(np.where(dx2dy2==dx2dy2.min()),smooth[np.where(dx2dy2==dx2dy2.min())],'ro')
plt.subplot(122)
plt.plot(dxdy,'r',label='first derivate')
plt.plot(dx2dy2,'b',label='second derivate')
plt.legend()

Answer 2

我想如果你能找到一种方法来平滑曲线，你就会得到你想要的东西。