code: def test_coverage_union_reduce_axis():
      # shape = (3, 2), all polygons - none of them overlapping
      data = [
          [shapely.box(i, j, i + 1, j + 1) for i in range(2)]
          for j in range(3)
      ]
      actual = shapely.coverage_union_all(data, axis=None)  # default
      assert isinstance(actual, Geometry)
      actual = shapely.coverage_union_all(data, axis=0)
      assert actual.shape == (2,)
      actual = shapely.coverage_union_all(data, axis=1)
      assert actual.shape == (3,)
      actual = shapely.coverage_union_all(data, axis=-1)
      assert actual.shape == (3,)
  
imports: import numpy
  import pytest
  import shapely
  from shapely import Geometry, GeometryCollection, Polygon
  from shapely.errors import UnsupportedGEOSVersionError
  from shapely.testing import assert_geometries_equal
  from shapely.tests.common import all_types, empty, ignore_invalid, point, polygon
label_desc: Use image node_id 23 for calling a global function with example usage: test_coverage_union_reduce_axis() without return types