estimateRigidTransform returns non sensical values

asked 2013-04-03 16:23:59 -0600

HIy
11 ●2 ●3

updated 2013-04-03 16:43:31 -0600

Guanta
6736 ●6 ●25 ●79

Hi all,

I am trying to estimate an affine transform from two point sets. Here is the sample code that I am running to verify estimateRigidTransform.

std::vector< Point2f > obj1;
std::vector< Point2f > scene1;
for (int i=0;i<10;i++)
{
    obj1.push_back(Point2f(i,i+10));
    scene1.push_back(Point2f(i+10,i+20));
}
Mat H1 = estimateRigidTransform(obj1,scene1,false);
cout << H1 << endl;

The output I get is:

[-4.836697827223e-026, -4.836697827223e-026, -4.836697827223e-026;
  -4.836697827223e-026, -4.836697827223e-026, -4.836697827223e-026]

What am I doing wrong?

Also, I would like to estimate the transformation robustly. Is there any sample code there for using RANSAC/or other robust methods on top of estimateRigidTransform?

Thanks HIy

edit retag flag offensive close merge delete

add a comment

2 answers

Sort by » oldest newest most voted

answered 2013-08-14 04:25:07 -0600

Neptilo
21 ●1 ●2

Your code looks right. However, to estimate an affine transformation you need at least 3 points, but if these points are all aligned, there is an infinity of solutions, which seems to lead the estimation to result in a degenerate solution.

Although that might no longer be true as you set fullAffine to false, perhaps the function is implemented in such a way that it still needs points that are not aligned. But that is a supposition.

What is the result if you add another point? For example:

obj1.push_back(Point2f(0, 0));
scene1.push_back(Point2f(10, 10));

edit flag offensive delete link

add a comment

answered 2013-04-19 10:33:37 -0600

Guanta
6736 ●6 ●25 ●79

I guess, this: http://docs.opencv.org/modules/calib3d/doc/camera_calibration_and_3d_reconstruction.html#findhomography is more what you are looking for, since you can specify also RANSAC there, or you write RANSAC your self and use http://docs.opencv.org/modules/imgproc/doc/geometric_transformations.html?highlight=affinetransform#getaffinetransform . Good luck!

edit flag offensive delete link

Comments

What if we really want to use estimateRigidTransform? What is the right way to use it?

Neptilo ( 2013-08-14 03:25:34 -0600 )edit

Haven't used estimateRigidTransform so far, but think @HIy actually calls it correctly. I also just checked the code (it is located in /modules/video/src/lkpyramid.cpp) and it seems that the transformation is actually computed via RANSAC but maybe its internal used thresholds are not very good chosen (e.g. RANSAC_GOOD_RATIO=0.5 seems pretty low). However RANSAC is really not difficult to compute on your own, so I'd use getAffineTransform with a self-implemented RANSAC-version.

Guanta ( 2013-08-14 07:23:52 -0600 )edit

I just downloaded the source code lkpyramid.cpp and set RANSAC's local parameters as arguments so I can modify whatever I want, haha!

Neptilo ( 2013-08-14 08:19:44 -0600 )edit

@Neptilo: Hi, I have the same problem as you had so I also took the source from lkpyramid.cpp. Have you gotten it to work? Could you post a simple working example in that case? I suck at debugging :(

Tobbe ( 2014-02-26 04:10:34 -0600 )edit

add a comment

estimateRigidTransform returns non sensical values

2 answers

Comments

Links

Question Tools

Stats

Related questions

estimateRigidTransform returns non sensical values edit

2 answers

Comments

Links

Question Tools

Stats

Related questions

estimateRigidTransform returns non sensical values