程序不打印结果

Question

现在，我正在使用以下程序，但我无法获得所需的图像时刻、质心等。

代码：

import cv2

import numpy

from matplotlib.pyplot import imread

from numpy import mgrid, sum

image = 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

你能帮我解决这个问题吗？非常感谢。

Answer 1

乍一看，我看到的一个问题是，您已经定义了函数moments2e(image)但没有调用它。您需要moments2e(image)在函数定义之外调用该函数。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]]    .    .    .    return momentsmoments = moments2e(image)print(moments)