opencv/opencv_contrib
1
0
Fork 0
mirror of https://github.com/opencv/opencv_contrib.git synced 2025-10-15 20:37:07 +08:00
Code Issues Releases Wiki Activity

fix for 1-camera and stereocalib

Browse Source
This commit is contained in:
Rostislav Vasilikhin
2023-01-26 03:03:33 +01:00
parent b01f2f90a9
commit 6b15ef3f62
1 changed files with 8 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
modules/ccalib/src/omnidir.cpp
View File

@@ -1782,13 +1782,15 @@ void cv::omnidir::internal::estimateUncertainties(InputArrayOfArrays objectPoint
Mat sigma_x; Mat sigma_x;
meanStdDev(reprojError.reshape(1,1), noArray(), sigma_x); meanStdDev(reprojError.reshape(1,1), noArray(), sigma_x);
sigma_x *= sqrt(2.0*(double)reprojError.total()/(2.0*(double)reprojError.total() - 1.0));
double s = sigma_x.at<double>(0);
Mat _JTJ_inv, _JTE; Mat _JTJ_inv, _JTE;
computeJacobian(objectPoints, imagePoints, parameters, _JTJ_inv, _JTE, flags, 0.0); computeJacobian(objectPoints, imagePoints, parameters, _JTJ_inv, _JTE, flags, 0.0);
sqrt(_JTJ_inv, _JTJ_inv); sqrt(_JTJ_inv, _JTJ_inv);
int nParams = _JTJ_inv.rows;
sigma_x *= sqrt(2.0*(double)reprojError.total()/(2.0*(double)reprojError.total() - nParams));
double s = sigma_x.at<double>(0);
errors = 3 * s * _JTJ_inv.diag(); errors = 3 * s * _JTJ_inv.diag();
rms = 0; rms = 0;
@@ -1868,13 +1870,15 @@ void cv::omnidir::internal::estimateUncertaintiesStereo(InputArrayOfArrays objec
Mat sigma_x; Mat sigma_x;
meanStdDev(reprojErrorAll.reshape(1,1), noArray(), sigma_x); meanStdDev(reprojErrorAll.reshape(1,1), noArray(), sigma_x);
sigma_x *= sqrt(2.0*(double)reprojErrorAll.total()/(2.0*(double)reprojErrorAll.total() - 1.0));
double s = sigma_x.at<double>(0);
Mat _JTJ_inv, _JTE; Mat _JTJ_inv, _JTE;
computeJacobianStereo(objectPoints, imagePoints1, imagePoints2, _parameters, _JTJ_inv, _JTE, flags, 0.0); computeJacobianStereo(objectPoints, imagePoints1, imagePoints2, _parameters, _JTJ_inv, _JTE, flags, 0.0);
cv::sqrt(_JTJ_inv, _JTJ_inv); cv::sqrt(_JTJ_inv, _JTJ_inv);
int nParams = _JTJ_inv.rows;
sigma_x *= sqrt(2.0*(double)reprojErrorAll.total()/(2.0*(double)reprojErrorAll.total() - nParams));
double s = sigma_x.at<double>(0);
errors = 3 * s * _JTJ_inv.diag(); errors = 3 * s * _JTJ_inv.diag();
rms = 0; rms = 0;