Accord.Imaging allows find SURF features and estimate points using RANSAC.
var correlationPoints1 = matches[0];
var correlationPoints2 = matches[1];
var ransac = new RansacHomographyEstimator(0.10, 0.99);
var homography = ransac.Estimate(correlationPoints1, correlationPoints2);
var inliers1 = correlationPoints1.Submatrix(ransac.Inliers);
var inliers2 = correlationPoints2.Submatrix(ransac.Inliers);
var result = new IntPoint[][]
{
inliers1,
inliers2
};