現(xiàn)在，我正在使用以下程序，但我無法獲得所需的圖像時刻、質(zhì)心等。代碼：import cv2import numpyfrom matplotlib.pyplot import imreadfrom numpy import mgrid, sumimage = imread('imagemoment.png')def moments2e(image):    assert len(image.shape) == 2  # only for grayscale images    x, y = mgrid[:image.shape[0], :image.shape[1]]    moments = {}    moments['mean_x'] = sum(x * image) / sum(image)    moments['mean_y'] = sum(y * image) / sum(image)    # raw or spatial moments    moments['m00'] = sum(image)    moments['m01'] = sum(x * image)    moments['m10'] = sum(y * image)    moments['m11'] = sum(y * x * image)    moments['m02'] = sum(x ** 2 * image)    moments['m20'] = sum(y ** 2 * image)    moments['m12'] = sum(x * y ** 2 * image)    moments['m21'] = sum(x ** 2 * y * image)    moments['m03'] = sum(x ** 3 * image)    moments['m30'] = sum(y ** 3 * image)    # central moments    moments['mu01']= sum((y-moments['mean_y'])*image) # should be 0    moments['mu10']= sum((x-moments['mean_x'])*image) # should be 0    moments['mu11'] = sum((x - moments['mean_x']) * (y - moments['mean_y']) * image)    moments['mu02'] = sum((y - moments['mean_y']) ** 2 * image)  # variance    moments['mu20'] = sum((x - moments['mean_x']) ** 2 * image)  # variance    moments['mu12'] = sum((x - moments['mean_x']) * (y - moments['mean_y']) ** 2 * image)    moments['mu21'] = sum((x - moments['mean_x']) ** 2 * (y - moments['mean_y']) * image)    moments['mu03'] = sum((y - moments['mean_y']) ** 3 * image)    moments['mu30'] = sum((x - moments['mean_x']) ** 3 * image)    # opencv versions    # moments['mu02'] = sum(image*(x-m01/m00)**2)    # moments['mu02'] = sum(image*(x-y)**2)    # wiki variations    # moments['mu02'] = m20 - mean_y*m10    # moments['mu20'] = m02 - mean_x*m01你能幫我解決這個問題嗎？非常感謝。