All connected vug regions are identified and labeled with a distinct positive number. The rest of the matrix is labeled 0.