We present a novel matching method to find the correspondences among different images containing the same object. In the proposed method, by considering each feature point-set as a matrix, two point-sets are projected onto a common subspace using modified projective nonnegative matrix factorization. The core idea of the proposed approach is to jointly factorize of the two feature matrices and the matching operate on embeddings of the two point-sets in the common subspace. Furthermore, it is robust to noise due to the merit of the subspace method. The proposed approach was tested for matching accuracy, and robustness to noise. Its performance on synthetic and real images was compared with state-of-the-art reference algorithms.