实际执行中轴变换？

Question

我已经看到了很多关于该算法所基于的解释，但我似乎找不到中轴变换本身的任何实际代码（伪代码或某种语言）。

对于某些情况（例如具有离散角的多边形），可以将其简化为 delaunay 三角形，但我更关心您有一个斑点（例如正在跟踪的人）并想要计算其骨架的情况。

这个的伪代码是什么样的？

Answer 1

这里你有来自pymorph 库：

基本上是：

def mmskelm(f, B=None, option="binary"):
    from string import upper
    from Numeric import asarray
    if B is None: B = mmsecross()
    assert mmisbinary(f),'Input binary image only'
    option = upper(option)
    k1,k2 = mmlimits(f)
    y = mmgray(mmintersec(f, k1),'uint16')
    iszero = asarray(y)
    nb = mmsesum(B,0)
    for r in range(1,65535):
        ero = mmero( f, nb)
        if mmisequal(ero, iszero): break
        f1 = mmopenth( ero, B)
        nb = mmsedil(nb, B)
        y = mmunion(y, mmgray(f1,'uint16',r))
    if option == 'BINARY':
        y = mmbinary(y)
    return y

当然它使用了另一个函数来自同一个包。

的示例

Mathematica

HTH！