写一个GIS中计算曲率的python算法

下面是一个简单的计算曲率的python算法：

import numpy asnp

def curvature(x, y):
    dx = np.gradient(x)
    dy = np.gradient(y)
    ddx = np.gradient(dx)
    ddy = np.gradient(dy)
    num = np.abs(ddx * dy - dx * ddy)
    den = np.power(dx * dx + dy * dy, 1.5)
    k = num / den
    return k

上面的算法使用了Numpy库