returns the Ids of neighboring grains
% neighbouring relationships within grains pairs = neighbors(grains) pairs = neighbors(grains1,grains2)
% all neighbouring relationships pairs = neighbors(grains,'full')