Read in all the domains and correct flying pixels
      and small holes by reassigning to the closest neighbour
      domains. For the holes this is the domain itself, for
      disconnected pixels it is the domain with greatest area
      of intersection with the dilated pixel. Each corrected
      domain is output to a file with extension _new, additions
      and deletions are output to files with extensions _adds
      an _dels respectively