We propose the linear Laplacian discrimination (LLD) algorithm/or discriminant feature extraction, which is an extension of linear discriminant analysis (LDA). Our motivation is to address the issue that LDA cannot work well in cases where sample spaces are non-Euclidean.