diff --git a/skimage/measure/_moments.pyx b/skimage/measure/_moments.pyx index 8cb5ae47..f5bb63af 100644 --- a/skimage/measure/_moments.pyx +++ b/skimage/measure/_moments.pyx @@ -39,7 +39,16 @@ def moments(double[:, :] image, Py_ssize_t order=3): .. [4] http://en.wikipedia.org/wiki/Image_moment """ - return moments_central(image, 0, 0, order) + cdef Py_ssize_t p, q, r, c + cdef double[:, ::1] mu = np.zeros((order + 1, order + 1), dtype=np.double) + cdef double val + for r in range(image.shape[0]): + for c in range(image.shape[1]): + val = image[r, c] + for p in range(order + 1): + for q in range(order + 1): + mu[p, q] += val * r ** q * c ** p + return np.asarray(mu) def moments_central(double[:, :] image, double cr, double cc,