code: def peak_XYZ(XYZ, Tvalues, labels, label_values):
      """
      Returns (3, n_labels) array of maximum T values coordinates for each label value
      """
      C = np.zeros((3, len(label_values)), int)
      for i in range(len(label_values)):
          I = np.where(labels == label_values[i])[0]
          C[:, i] = XYZ[:, I[np.argmax(Tvalues[I])]]
      return C
  
imports: import numpy
  import warnings
  from nipy.algorithms.graph import wgraph_from_3d_grid
  from nipy.algorithms.graph.field import Field, field_from_graph_and_data
  from ..utils import zscore
  from .onesample import stat
  from .onesample import stat_mfx
  from .twosample import stat
  from .twosample import stat_mfx
label_desc: Use image node_id 11 for calling a global function with example usage: peak_XYZ(XYZ, Tvalues, labels, label_values) and returns: C