code: def view_thresholdedT(
      design, contrast, threshold, inequality=np.greater
  ):
      """
      A mayavi isosurface view of thresholded t-statistics
  
      Parameters
      ----------
      design : {'block', 'event'}
      contrast : str
      threshold : float
      inequality : {np.greater, np.less}, optional
      """
      maska = np.asarray(MASK)
      tmap = np.array(
          load_image_fiac("group", design, contrast, "t.nii")
      )
      test = inequality(tmap, threshold)
      tval = np.zeros(tmap.shape)
      tval[test] = tmap[test]
  
      # XXX make the array axes agree with mayavi2
      avganata = np.array(AVGANAT)
      avganat_iso = mlab.contour3d(
          avganata * maska,
          opacity=0.3,
          contours=[3600],
          color=(0.8, 0.8, 0.8),
      )
  
      avganat_iso.actor.property.backface_culling = True
      avganat_iso.actor.property.ambient = 0.3
  
      tval_iso = mlab.contour3d(
          tval * MASK, color=(0.8, 0.3, 0.3), contours=[threshold]
      )
      return avganat_iso, tval_iso
  
label_desc: Use image node_id 1 for calling a global function with example usage: view_thresholdedT(design, contrast, threshold, inequality) and returns: avganat_iso,  tval_iso