code: def _get_atlas_values(vol_info, rr):
      # Transform MRI coordinates (where our surfaces live) to voxels
      rr_vox = apply_trans(vol_info["mri_vox_t"], rr)
      good = (
          (rr_vox >= -0.5)
          & (rr_vox < np.array(vol_info["data"].shape, int) - 0.5)
      ).all(-1)
      idx = np.round(rr_vox[good].T).astype(np.int64)
      values = np.full(rr.shape[0], np.nan)
      values[good] = vol_info["data"][tuple(idx)]
      return values
  
imports: import os.path
  from gzip import GzipFile
  from pathlib import Path
  import numpy
  from ._fiff.constants import FIFF
  from ._fiff.meas_info import read_fiducials
  from .surface import _read_mri_surface, read_surface
  from .transforms import Transform, _ensure_trans, apply_trans, combine_transforms, invert_transform, read_ras_mni_t
  from .utils import _check_fname, _check_option, _import_nibabel, _validate_type, get_subjects_dir, logger, verbose
label_desc: Use image node_id 7 for calling a global function with example usage: _get_atlas_values(vol_info, rr) and returns: values