Could it be that this should actually be changed from Q[:, 4:] *= 2
to Q[4:, :] *= 2
? As we don't want to multiply by two in the pixel axis.
And shouldn't this be changed from VQ[:,9]
to VQ[9, :]
? As we want to get the eigenvector in VQ.T
to the smallest singular value of Q
.