code: def demo_plot_map_3d():
      map = np.zeros((182, 218, 182))
      # Color a asymmetric rectangle around Broca area:
      x, y, z = -52, 10, 22
      x_map, y_map, z_map = coord_transform(x, y, z, mni_sform_inv)
      map[
          x_map - 5 : x_map + 5,
          y_map - 3 : y_map + 3,
          z_map - 10 : z_map + 10,
      ] = 1
      plot_map_3d(map, mni_sform, cut_coords=(x, y, z))
  
file_path: nipy/nipy/labs/viz_tools/maps_3d.py
imports: import os
  import tempfile
  import numpy
  from scipy import stats
  from .anat_cache import _AnatCache, mni_sform, mni_sform_inv
  from .coord_tools import coord_transform
label_desc: Use image node_id 6 for calling a global function with example usage: demo_plot_map_3d() without return types