code: def test_moments_coords():
      image = np.zeros((20, 20), dtype=np.float64)
      image[13:17, 13:17] = 1
      mu_image = moments(image)
  
      coords = np.array(
          [[r, c] for r in range(13, 17) for c in range(13, 17)],
          dtype=np.float64,
      )
      mu_coords = moments_coords(coords)
      assert_almost_equal(mu_coords, mu_image)
  
imports: import itertools
  import numpy
  import pytest
  from scipy import ndimage
  from skimage import draw
  from skimage._shared import testing
  from skimage._shared.testing import assert_allclose, assert_almost_equal, assert_equal
  from skimage._shared.utils import _supported_float_type
  from skimage.measure import centroid, inertia_tensor, inertia_tensor_eigvals, moments, moments_central, moments_coords, moments_coords_central, moments_hu, moments_normalized
label_desc: Use image node_id 4 for calling a global function with example usage: test_moments_coords() without return types