code: def moments_hu(nu):
      hu = np.zeros((7,), dtype=np.double)
      t0 = nu[3, 0] + nu[1, 2]
      t1 = nu[2, 1] + nu[0, 3]
      q0 = t0 * t0
      q1 = t1 * t1
      n4 = 4 * nu[1, 1]
      s = nu[2, 0] + nu[0, 2]
      d = nu[2, 0] - nu[0, 2]
      hu[0] = s
      hu[1] = d * d + n4 * nu[1, 1]
      hu[3] = q0 + q1
      hu[5] = d * (q0 - q1) + n4 * t0 * t1
      t0 *= q0 - 3 * q1
      t1 *= 3 * q0 - q1
      q0 = nu[3, 0] - 3 * nu[1, 2]
      q1 = 3 * nu[2, 1] - nu[0, 3]
      hu[2] = q0 * q0 + q1 * q1
      hu[4] = q0 * t0 + q1 * t1
      hu[6] = q1 * t0 - q0 * t1
      return np.asarray(hu)
  
label_desc: Use image node_id 1 for calling a global function with example usage: moments_hu(nu) and returns: np