Duplicated factor values

It seems to me that angle number 0 and angle number 39 are the same? Did you use linspace(0, 360, 40) to get the angles? That's definitely wrong. Or did you do this deliberately to get around the rotation symmetry of squares and ellipses?