where L xx , L xy = L y x , and L y y are second-order derivatives of local image patch L , and eig() denotes eigenvalue decomposition. We choose to use SI because of its invariance property with respect to rotation, absolute gray value, and translation.