1. A method, comprising: integrating an unsupervised feature selection algorithm with a least-square based congealing algorithm; selecting a subset of features from an initial feature representation with the unsupervised feature selection algorithm; and executing the least-square based congealing algorithm to estimate warping parameters for a plurality of images in an ensemble using the subset of features.