#include "../FMatrix.hh"
#include "../FMatrixEstimation.hh"
#include "../PMatrixEstimation.hh"
#include "../PMatrixLinear.hh"
#include "../Triangulation.hh"
#include "../../Base/Math/Random.hh"
using namespace BIAS;
using namespace std;
int main()
{
std::vector<FMatrix> vecF;
for (unsigned int i=0; i<1; i++) {
for (unsigned int x=0; x<4; x++)
for (unsigned int y=0; y<3; y++) {
}
vector<HomgPoint3D> points;
vector<HomgPoint2D> p1,p2;
cout << "Randomized correspondences are "<<endl;
for (unsigned int p=0; p<7; p++) {
1.0),
1.0),
1.0)));
p1.push_back(P1 * points[p]); p1[p].Homogenize();
cout << "{ "<<p1[p]<<" ";
p2.push_back(P2 * points[p]); p2[p].Homogenize();
cout << p2[p]<<" } "<<endl;
}
bool NormalizeHartley = true;
vecF.clear();
cout << "Computing epipoles "<<endl;
for (unsigned int s=0; s<vecF.size(); s++) {
vecF[s].GetEpipolesHomogenized(e1, e2);
cout << "True epipoles are "<<theoe1<<" and "<<theoe2<<endl;
cout << "est. epipoles are "<<e1<<" and "<<e2<<endl;
}
}
return 0;
}