code: def test_gaussian_filter1d_higher_dims():
      a = np.zeros((120, 120))
      a[:60] = 2
      a1 = mh.gaussian_filter1d(a, 1.0, axis=0)
  
      for a3 in ([a], [a, a], [a, a, a]):
          a1_3d = mh.gaussian_filter1d(np.array(a3), 1.0, axis=1)
          a1_3dm = a1_3d.max(0)
          assert np.allclose(a1, a1_3dm)
          assert np.abs(a1 - a1_3dm).max() < 0.01
  
      a1_4d = mh.gaussian_filter1d(
          np.array([[a], [a], [a]]), 1.0, axis=2
      )
      a1_4dm = a1_4d.max(0).max(0)
      assert np.allclose(a1, a1_4dm)
      assert np.abs(a1 - a1_4dm).max() < 0.01
  
imports: import numpy
  import mahotas
  import mahotas.convolve
  import mahotas
  from mahotas.convolve import convolve1d, gaussian_filter
  import mahotas._filters
  from mahotas.tests.utils import luispedro_jpg
  import pytest
label_desc: Use image node_id 26 for calling a global function with example usage: test_gaussian_filter1d_higher_dims() without return types