I do have two sets of points and I want to find the best transformation between them. In OpenCV, you have the following function:
Mat H = Calib3d.findHomography
OpenCV estimateRigidTransform function is exactly what you need: it returns Translation, Rotation and Scale (use false value for fullAffine flag). And it DOES use RANSAC (see source code to be sure of it).