poster
2019-12-14, 20:13
可以使用bwtraceboundary (http://www.mathworks.com/help/toolbox/images/ref/bwtraceboundary.html)函数在Matlab中提取对象的参数边界。它返回一个Q×2矩阵B,其中Q是对象的边界像素数,第一和第二列分别存储边界像素的行和列坐标。
我想做的是对Q个元素的边界采样N个点,这些点将原始边界分成等长的弓段。
我认为直截了当的解决方案包括通过将所有两个连续边界像素的距离相加来计算边界的长度L。这些距离为1或sqrt(2)。然后我将L除以N,以找到所需的弧长。最后,我再次对边界进行迭代,求和所有两个连续边界像素的距离。当总和大于或等于期望的弧长时,将当前边界像素选为构成采样边界的N个像素之一。
那是一个好的解决方案吗?有没有更有效/更简单的解决方案?
回答:
多年来,我看过这个问题的次数似乎很多。所以我写了一个小工具可以做到这一点。以一般尺寸采样分段线性或曲线(样条曲线)弧,以使连续点沿该弧均匀或指定距离。
在仅使用分段线性弧的情况下,这相当容易。您可以对曲线的总弧长求和,然后对弧长进行插值,但是由于已知分段线性,因此只需要沿着该长度进行线性插值即可,这是累积弧长的函数。
在弧形的情况下,最容易将其作为常微分方程组的解来完成,并观察沿途的事件。 ODE45很好地做到了这一点。
您可以使用在MATLAB Central File Exchange上找到的interparc (http://www.mathworks.com/matlabcentral/fileexchange/27096-interparc)为您执行此操作,或者,如果您想学习简单的分段线性情况,请仔细阅读代码的第一部分,其中我进行了分段线性弧长插补。一件好事是线性情况是以完全矢量化的形式完成的,因此不需要显式循环。
更多&回答... (https://stackoverflow.com/questions/4690215)
我想做的是对Q个元素的边界采样N个点,这些点将原始边界分成等长的弓段。
我认为直截了当的解决方案包括通过将所有两个连续边界像素的距离相加来计算边界的长度L。这些距离为1或sqrt(2)。然后我将L除以N,以找到所需的弧长。最后,我再次对边界进行迭代,求和所有两个连续边界像素的距离。当总和大于或等于期望的弧长时,将当前边界像素选为构成采样边界的N个像素之一。
那是一个好的解决方案吗?有没有更有效/更简单的解决方案?
回答:
多年来,我看过这个问题的次数似乎很多。所以我写了一个小工具可以做到这一点。以一般尺寸采样分段线性或曲线(样条曲线)弧,以使连续点沿该弧均匀或指定距离。
在仅使用分段线性弧的情况下,这相当容易。您可以对曲线的总弧长求和,然后对弧长进行插值,但是由于已知分段线性,因此只需要沿着该长度进行线性插值即可,这是累积弧长的函数。
在弧形的情况下,最容易将其作为常微分方程组的解来完成,并观察沿途的事件。 ODE45很好地做到了这一点。
您可以使用在MATLAB Central File Exchange上找到的interparc (http://www.mathworks.com/matlabcentral/fileexchange/27096-interparc)为您执行此操作,或者,如果您想学习简单的分段线性情况,请仔细阅读代码的第一部分,其中我进行了分段线性弧长插补。一件好事是线性情况是以完全矢量化的形式完成的,因此不需要显式循环。
更多&回答... (https://stackoverflow.com/questions/4690215)