code: def test_voxel_neighbors():
      """Test finding points above a threshold near a seed location."""
      locs = np.array(np.meshgrid(*[np.linspace(-1, 1, 101)] * 3))
      image = 1 - np.linalg.norm(locs, axis=0)
      true_volume = set(
          [tuple(coord) for coord in np.array(np.where(image > 0.95)).T]
      )
      volume = _voxel_neighbors(
          np.array([-0.3, 0.6, 0.5])
          + (np.array(image.shape[0]) - 1) / 2,
          image,
          thresh=0.95,
          use_relative=False,
      )
      assert volume.difference(true_volume) == set()
      assert true_volume.difference(volume) == set()
  
imports: from pathlib import Path
  import numpy
  import pytest
  from numpy.testing import assert_allclose, assert_array_equal, assert_equal
  from mne import decimate_surface, dig_mri_distances, get_montage_volume_labels, pick_types, read_surface, write_surface
  from mne._fiff.constants import FIFF
  from mne.channels import make_dig_montage
  from mne.datasets import testing
  from mne.io import read_info
  from mne.surface import _compute_nearest, _get_ico_surface, _marching_cubes, _normal_orth, _project_onto_surface, _read_patch, _tessellate_sphere, _voxel_neighbors, fast_cross_3d, get_head_surf, get_meg_helmet_surf, read_curvature
  from mne.transforms import _get_trans
  from mne.utils import _record_warnings, catch_logging, object_diff, requires_freesurfer
label_desc: Use image node_id 13 for calling a global function with example usage: test_voxel_neighbors() without return types