code: def demo_plot_map(do3d=False, **kwargs):
      """Demo activation map plotting."""
      map = np.zeros((182, 218, 182))
      # Color a asymmetric rectangle around Broca area:
      x, y, z = -52, 10, 22
      mapped = coord_transform(x, y, z, mni_sform_inv)
      x_map, y_map, z_map = (int(v) for v in mapped)
      # Compare to values obtained using fslview. We need to add one as
      # voxels do not start at 0 in fslview.
      assert x_map == 142
      assert y_map + 1 == 137
      assert z_map + 1 == 95
      map[
          x_map - 5 : x_map + 5,
          y_map - 3 : y_map + 3,
          z_map - 10 : z_map + 10,
      ] = 1
      return plot_map(
          map,
          mni_sform,
          threshold="auto",
          title="Broca's area",
          do3d=do3d,
          **kwargs,
      )
  
file_path: nipy/nipy/labs/viz_tools/activation_maps.py
imports: import numbers
  import warnings
  import matplotlib.pyplot
  import numpy
  from nipy.utils import is_numlike
  from .anat_cache import _AnatCache, mni_sform, mni_sform_inv
  from .coord_tools import coord_transform, find_maxsep_cut_coords
  from .edge_detect import _fast_abs_percentile
  from .slicers import SLICERS, _xyz_order
label_desc: Use image node_id 4 for calling a global function with example usage: demo_plot_map(do3d) and returns: plot_map