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.
