code: def _generate_icosahedron():
      x = np.array(
          [[0, -1, -phi], [0, -1, +phi], [0, +1, -phi], [0, +1, +phi]]
      )
      return np.concatenate([np.roll(x, i, axis=1) for i in range(3)])
  
imports: import pytest
  import numpy
  from numpy.testing import assert_array_almost_equal
  from scipy.spatial.transform import Rotation
  from scipy.optimize import linear_sum_assignment
  from scipy.spatial.distance import cdist
  from scipy.constants import golden
  from scipy.spatial import cKDTree
label_desc: Use image node_id 4 for calling a global function with example usage: _generate_icosahedron() and returns: np