mirror of
https://github.com/opencv/opencv.git
synced 2024-11-25 11:40:44 +08:00
223 lines
8.5 KiB
C++
223 lines
8.5 KiB
C++
/*M///////////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
|
|
//
|
|
// By downloading, copying, installing or using the software you agree to this license.
|
|
// If you do not agree to this license, do not download, install,
|
|
// copy or use the software.
|
|
//
|
|
//
|
|
// Intel License Agreement
|
|
// For Open Source Computer Vision Library
|
|
//
|
|
// Copyright (C) 2000, Intel Corporation, all rights reserved.
|
|
// Third party copyrights are property of their respective owners.
|
|
//
|
|
// Redistribution and use in source and binary forms, with or without modification,
|
|
// are permitted provided that the following conditions are met:
|
|
//
|
|
// * Redistribution's of source code must retain the above copyright notice,
|
|
// this list of conditions and the following disclaimer.
|
|
//
|
|
// * Redistribution's in binary form must reproduce the above copyright notice,
|
|
// this list of conditions and the following disclaimer in the documentation
|
|
// and/or other materials provided with the distribution.
|
|
//
|
|
// * The name of Intel Corporation may not be used to endorse or promote products
|
|
// derived from this software without specific prior written permission.
|
|
//
|
|
// This software is provided by the copyright holders and contributors "as is" and
|
|
// any express or implied warranties, including, but not limited to, the implied
|
|
// warranties of merchantability and fitness for a particular purpose are disclaimed.
|
|
// In no event shall the Intel Corporation or contributors be liable for any direct,
|
|
// indirect, incidental, special, exemplary, or consequential damages
|
|
// (including, but not limited to, procurement of substitute goods or services;
|
|
// loss of use, data, or profits; or business interruption) however caused
|
|
// and on any theory of liability, whether in contract, strict liability,
|
|
// or tort (including negligence or otherwise) arising in any way out of
|
|
// the use of this software, even if advised of the possibility of such damage.
|
|
//
|
|
//M*/
|
|
|
|
#include "test_precomp.hpp"
|
|
#include "opencv2/calib3d/calib3d_c.h"
|
|
|
|
using namespace cv;
|
|
using namespace std;
|
|
|
|
class CV_POSITTest : public cvtest::BaseTest
|
|
{
|
|
public:
|
|
CV_POSITTest();
|
|
protected:
|
|
void run(int);
|
|
};
|
|
|
|
|
|
CV_POSITTest::CV_POSITTest()
|
|
{
|
|
test_case_count = 20;
|
|
}
|
|
|
|
void CV_POSITTest::run( int start_from )
|
|
{
|
|
int code = cvtest::TS::OK;
|
|
|
|
/* fixed parameters output */
|
|
/*float rot[3][3]={ 0.49010f, 0.85057f, 0.19063f,
|
|
-0.56948f, 0.14671f, 0.80880f,
|
|
0.65997f, -0.50495f, 0.55629f };
|
|
|
|
float trans[3] = { 0.0f, 0.0f, 40.02637f };
|
|
*/
|
|
|
|
/* Some variables */
|
|
int i, counter;
|
|
|
|
CvTermCriteria criteria;
|
|
CvPoint3D32f* obj_points;
|
|
CvPoint2D32f* img_points;
|
|
CvPOSITObject* object;
|
|
|
|
float angleX, angleY, angleZ;
|
|
RNG& rng = ts->get_rng();
|
|
int progress = 0;
|
|
|
|
CvMat* true_rotationX = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* true_rotationY = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* true_rotationZ = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* tmp_matrix = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* true_rotation = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* rotation = cvCreateMat( 3, 3, CV_32F );
|
|
CvMat* translation = cvCreateMat( 3, 1, CV_32F );
|
|
CvMat* true_translation = cvCreateMat( 3, 1, CV_32F );
|
|
|
|
const float flFocalLength = 760.f;
|
|
const float flEpsilon = 0.5f;
|
|
|
|
/* Initilization */
|
|
criteria.type = CV_TERMCRIT_EPS|CV_TERMCRIT_ITER;
|
|
criteria.epsilon = flEpsilon;
|
|
criteria.max_iter = 10000;
|
|
|
|
/* Allocating source arrays; */
|
|
obj_points = (CvPoint3D32f*)cvAlloc( 8 * sizeof(CvPoint3D32f) );
|
|
img_points = (CvPoint2D32f*)cvAlloc( 8 * sizeof(CvPoint2D32f) );
|
|
|
|
/* Fill points arrays with values */
|
|
|
|
/* cube model with edge size 10 */
|
|
obj_points[0].x = 0; obj_points[0].y = 0; obj_points[0].z = 0;
|
|
obj_points[1].x = 10; obj_points[1].y = 0; obj_points[1].z = 0;
|
|
obj_points[2].x = 10; obj_points[2].y = 10; obj_points[2].z = 0;
|
|
obj_points[3].x = 0; obj_points[3].y = 10; obj_points[3].z = 0;
|
|
obj_points[4].x = 0; obj_points[4].y = 0; obj_points[4].z = 10;
|
|
obj_points[5].x = 10; obj_points[5].y = 0; obj_points[5].z = 10;
|
|
obj_points[6].x = 10; obj_points[6].y = 10; obj_points[6].z = 10;
|
|
obj_points[7].x = 0; obj_points[7].y = 10; obj_points[7].z = 10;
|
|
|
|
/* Loop for test some random object positions */
|
|
for( counter = start_from; counter < test_case_count; counter++ )
|
|
{
|
|
ts->update_context( this, counter, true );
|
|
progress = update_progress( progress, counter, test_case_count, 0 );
|
|
|
|
/* set all rotation matrix to zero */
|
|
cvZero( true_rotationX );
|
|
cvZero( true_rotationY );
|
|
cvZero( true_rotationZ );
|
|
|
|
/* fill random rotation matrix */
|
|
angleX = (float)(cvtest::randReal(rng)*2*CV_PI);
|
|
angleY = (float)(cvtest::randReal(rng)*2*CV_PI);
|
|
angleZ = (float)(cvtest::randReal(rng)*2*CV_PI);
|
|
|
|
true_rotationX->data.fl[0 *3+ 0] = 1;
|
|
true_rotationX->data.fl[1 *3+ 1] = (float)cos(angleX);
|
|
true_rotationX->data.fl[2 *3+ 2] = true_rotationX->data.fl[1 *3+ 1];
|
|
true_rotationX->data.fl[1 *3+ 2] = -(float)sin(angleX);
|
|
true_rotationX->data.fl[2 *3+ 1] = -true_rotationX->data.fl[1 *3+ 2];
|
|
|
|
true_rotationY->data.fl[1 *3+ 1] = 1;
|
|
true_rotationY->data.fl[0 *3+ 0] = (float)cos(angleY);
|
|
true_rotationY->data.fl[2 *3+ 2] = true_rotationY->data.fl[0 *3+ 0];
|
|
true_rotationY->data.fl[0 *3+ 2] = -(float)sin(angleY);
|
|
true_rotationY->data.fl[2 *3+ 0] = -true_rotationY->data.fl[0 *3+ 2];
|
|
|
|
true_rotationZ->data.fl[2 *3+ 2] = 1;
|
|
true_rotationZ->data.fl[0 *3+ 0] = (float)cos(angleZ);
|
|
true_rotationZ->data.fl[1 *3+ 1] = true_rotationZ->data.fl[0 *3+ 0];
|
|
true_rotationZ->data.fl[0 *3+ 1] = -(float)sin(angleZ);
|
|
true_rotationZ->data.fl[1 *3+ 0] = -true_rotationZ->data.fl[0 *3+ 1];
|
|
|
|
cvMatMul( true_rotationX, true_rotationY, tmp_matrix);
|
|
cvMatMul( tmp_matrix, true_rotationZ, true_rotation);
|
|
|
|
/* fill translation vector */
|
|
true_translation->data.fl[2] = (float)(cvtest::randReal(rng)*(2*flFocalLength-40) + 60);
|
|
true_translation->data.fl[0] = (float)((cvtest::randReal(rng)*2-1)*true_translation->data.fl[2]);
|
|
true_translation->data.fl[1] = (float)((cvtest::randReal(rng)*2-1)*true_translation->data.fl[2]);
|
|
|
|
/* calculate perspective projection */
|
|
for ( i = 0; i < 8; i++ )
|
|
{
|
|
float vec[3];
|
|
CvMat Vec = cvMat( 3, 1, CV_32F, vec );
|
|
CvMat Obj_point = cvMat( 3, 1, CV_32F, &obj_points[i].x );
|
|
|
|
cvMatMul( true_rotation, &Obj_point, &Vec );
|
|
|
|
vec[0] += true_translation->data.fl[0];
|
|
vec[1] += true_translation->data.fl[1];
|
|
vec[2] += true_translation->data.fl[2];
|
|
|
|
img_points[i].x = flFocalLength * vec[0] / vec[2];
|
|
img_points[i].y = flFocalLength * vec[1] / vec[2];
|
|
}
|
|
|
|
/*img_points[0].x = 0 ; img_points[0].y = 0;
|
|
img_points[1].x = 80; img_points[1].y = -93;
|
|
img_points[2].x = 245;img_points[2].y = -77;
|
|
img_points[3].x = 185;img_points[3].y = 32;
|
|
img_points[4].x = 32; img_points[4].y = 135;
|
|
img_points[5].x = 99; img_points[5].y = 35;
|
|
img_points[6].x = 247; img_points[6].y = 62;
|
|
img_points[7].x = 195; img_points[7].y = 179;
|
|
*/
|
|
|
|
object = cvCreatePOSITObject( obj_points, 8 );
|
|
cvPOSIT( object, img_points, flFocalLength, criteria,
|
|
rotation->data.fl, translation->data.fl );
|
|
cvReleasePOSITObject( &object );
|
|
|
|
Mat _rotation = cvarrToMat(rotation), _true_rotation = cvarrToMat(true_rotation);
|
|
Mat _translation = cvarrToMat(translation), _true_translation = cvarrToMat(true_translation);
|
|
code = cvtest::cmpEps2( ts, _rotation, _true_rotation, flEpsilon, false, "rotation matrix" );
|
|
if( code < 0 )
|
|
break;
|
|
|
|
code = cvtest::cmpEps2( ts, _translation, _true_translation, flEpsilon, false, "translation vector" );
|
|
if( code < 0 )
|
|
break;
|
|
}
|
|
|
|
cvFree( &obj_points );
|
|
cvFree( &img_points );
|
|
|
|
cvReleaseMat( &true_rotationX );
|
|
cvReleaseMat( &true_rotationY );
|
|
cvReleaseMat( &true_rotationZ );
|
|
cvReleaseMat( &tmp_matrix );
|
|
cvReleaseMat( &true_rotation );
|
|
cvReleaseMat( &rotation );
|
|
cvReleaseMat( &translation );
|
|
cvReleaseMat( &true_translation );
|
|
|
|
if( code < 0 )
|
|
ts->set_failed_test_info( code );
|
|
}
|
|
|
|
TEST(Calib3d_POSIT, accuracy) { CV_POSITTest test; test.safe_run(); }
|
|
|
|
/* End of file. */
|