As the voxels in the input data are not always square, the index in the 3D
cube has to be used to calculate the real coordinates based on the voxels'
dimensions.
init starred zeroes
for each column find the first zero if there is no other starred zero in
that row then star the zero, cover the column and row and go onto the
next column