查找数字数组曲线的折点

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

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