GRRRRR: Scipy ndimage.measurement uses 1-D iteration order...

... for assignments
GRR: Fully consider scalar arrays as array indices

This will break matrix[:,np.array(1)] = vals (transposed lhs
during assignment)