SYNOPSIS:
FIND_DG_NEIGHBOURS inputset outfd
DESCRIPTION:
Find all voxels from inputset, that have neighbours (also in
inputset) which differ in depth or grey value.
Write all such pairs to outfd.