% create random 3D data data = randn(20, 20, 20); % create random coordinate to sample from coords = 1 + 19 .* rand(1000, 3); % sample sampled = flexinterpn(data, coords);