Open yangqingnian opened 4 months ago
I'm unable to reproduce your issue.
I added some printings in tutorial-apriltag-detector-live-rgbd-realsense.cpp
. Can you copy/past and try again.
//! \example tutorial-apriltag-detector-live-rgbd-realsense.cpp
#include <visp3/core/vpConfig.h>
#ifdef VISP_HAVE_MODULE_SENSOR
#include <visp3/sensor/vpRealSense2.h>
#endif
//! [Include]
#include <visp3/detection/vpDetectorAprilTag.h>
//! [Include]
#include <visp3/core/vpImageConvert.h>
#include <visp3/gui/vpDisplayGDI.h>
#include <visp3/gui/vpDisplayOpenCV.h>
#include <visp3/gui/vpDisplayX.h>
#include <visp3/vision/vpPose.h>
int main(int argc, const char **argv)
{
//! [Macro defined]
#if defined(VISP_HAVE_APRILTAG) && defined(VISP_HAVE_REALSENSE2)
//! [Macro defined]
vpDetectorAprilTag::vpAprilTagFamily tagFamily = vpDetectorAprilTag::TAG_36h11;
vpDetectorAprilTag::vpPoseEstimationMethod poseEstimationMethod = vpDetectorAprilTag::HOMOGRAPHY_VIRTUAL_VS;
double tagSize = 0.053;
float quad_decimate = 1.0;
int nThreads = 1;
bool display_tag = false;
int color_id = -1;
unsigned int thickness = 2;
bool align_frame = false;
bool opt_verbose = true;
#if !(defined(VISP_HAVE_X11) || defined(VISP_HAVE_GDI) || defined(VISP_HAVE_OPENCV))
bool display_off = true;
std::cout << "Warning: There is no 3rd party (X11, GDI or openCV) to dislay images..." << std::endl;
#else
bool display_off = false;
#endif
for (int i = 1; i < argc; i++) {
if (std::string(argv[i]) == "--pose_method" && i + 1 < argc) {
poseEstimationMethod = (vpDetectorAprilTag::vpPoseEstimationMethod)atoi(argv[i + 1]);
}
else if (std::string(argv[i]) == "--tag_size" && i + 1 < argc) {
tagSize = atof(argv[i + 1]);
}
else if (std::string(argv[i]) == "--quad_decimate" && i + 1 < argc) {
quad_decimate = (float)atof(argv[i + 1]);
}
else if (std::string(argv[i]) == "--nthreads" && i + 1 < argc) {
nThreads = atoi(argv[i + 1]);
}
else if (std::string(argv[i]) == "--display_tag") {
display_tag = true;
}
else if (std::string(argv[i]) == "--display_off") {
display_off = true;
}
else if (std::string(argv[i]) == "--color" && i + 1 < argc) {
color_id = atoi(argv[i + 1]);
}
else if (std::string(argv[i]) == "--thickness" && i + 1 < argc) {
thickness = (unsigned int)atoi(argv[i + 1]);
}
else if (std::string(argv[i]) == "--tag_family" && i + 1 < argc) {
tagFamily = (vpDetectorAprilTag::vpAprilTagFamily)atoi(argv[i + 1]);
}
else if (std::string(argv[i]) == "--z_aligned") {
align_frame = true;
}
else if (std::string(argv[i]) == "--help" || std::string(argv[i]) == "-h") {
std::cout << "Usage: " << argv[0]
<< " [--tag_size <tag_size in m> (default: 0.053)]"
" [--quad_decimate <quad_decimate> (default: 1)]"
" [--nthreads <nb> (default: 1)]"
" [--pose_method <method> (0: HOMOGRAPHY, 1: HOMOGRAPHY_VIRTUAL_VS, "
" 2: DEMENTHON_VIRTUAL_VS, 3: LAGRANGE_VIRTUAL_VS, "
" 4: BEST_RESIDUAL_VIRTUAL_VS, 5: HOMOGRAPHY_ORTHOGONAL_ITERATION) (default: 0)]"
" [--tag_family <family> (0: TAG_36h11, 1: TAG_36h10 (DEPRECATED), 2: TAG_36ARTOOLKIT (DEPRECATED),"
" 3: TAG_25h9, 4: TAG_25h7 (DEPRECATED), 5: TAG_16h5, 6: TAG_CIRCLE21h7, 7: TAG_CIRCLE49h12,"
" 8: TAG_CUSTOM48h12, 9: TAG_STANDARD41h12, 10: TAG_STANDARD52h13) (default: 0)]"
" [--display_tag] [--z_aligned]";
#if (defined(VISP_HAVE_X11) || defined(VISP_HAVE_GDI) || defined(VISP_HAVE_OPENCV))
std::cout << " [--display_off] [--color <color id>] [--thickness <line thickness>]";
#endif
std::cout << " [--help]" << std::endl;
return EXIT_SUCCESS;
}
}
try {
//! [Construct grabber]
std::cout << "Use Realsense 2 grabber" << std::endl;
vpRealSense2 g;
rs2::config config;
unsigned int width = 640, height = 480;
config.enable_stream(RS2_STREAM_COLOR, static_cast<int>(width), static_cast<int>(height), RS2_FORMAT_RGBA8, 30);
config.enable_stream(RS2_STREAM_DEPTH, static_cast<int>(width), static_cast<int>(height), RS2_FORMAT_Z16, 30);
config.enable_stream(RS2_STREAM_INFRARED, static_cast<int>(width), static_cast<int>(height), RS2_FORMAT_Y8, 30);
vpImage<unsigned char> I;
vpImage<vpRGBa> I_color(height, width);
vpImage<uint16_t> I_depth_raw(height, width);
vpImage<vpRGBa> I_depth;
g.open(config);
const float depth_scale = g.getDepthScale();
std::cout << "I_color: " << I_color.getWidth() << " " << I_color.getHeight() << std::endl;
std::cout << "I_depth_raw: " << I_depth_raw.getWidth() << " " << I_depth_raw.getHeight() << std::endl;
rs2::align align_to_color = RS2_STREAM_COLOR;
g.acquire(reinterpret_cast<unsigned char *>(I_color.bitmap), reinterpret_cast<unsigned char *>(I_depth_raw.bitmap),
nullptr, nullptr, &align_to_color);
std::cout << "Read camera parameters from Realsense device" << std::endl;
vpCameraParameters cam;
cam = g.getCameraParameters(RS2_STREAM_COLOR, vpCameraParameters::perspectiveProjWithoutDistortion);
//! [Construct grabber]
std::cout << cam << std::endl;
std::cout << "poseEstimationMethod: " << poseEstimationMethod << std::endl;
std::cout << "tagFamily: " << tagFamily << std::endl;
std::cout << "nThreads : " << nThreads << std::endl;
std::cout << "Z aligned: " << align_frame << std::endl;
vpImage<vpRGBa> I_color2 = I_color;
vpImage<float> depthMap;
vpImageConvert::createDepthHistogram(I_depth_raw, I_depth);
vpDisplay *d1 = nullptr;
vpDisplay *d2 = nullptr;
vpDisplay *d3 = nullptr;
if (!display_off) {
#ifdef VISP_HAVE_X11
d1 = new vpDisplayX(I_color, 100, 30, "Pose from Homography");
d2 = new vpDisplayX(I_color2, I_color.getWidth() + 120, 30, "Pose from RGBD fusion");
d3 = new vpDisplayX(I_depth, 100, I_color.getHeight() + 70, "Depth");
#elif defined(VISP_HAVE_GDI)
d1 = new vpDisplayGDI(I_color, 100, 30, "Pose from Homography");
d2 = new vpDisplayGDI(I_color2, I_color.getWidth() + 120, 30, "Pose from RGBD fusion");
d3 = new vpDisplayGDI(I_depth, 100, I_color.getHeight() + 70, "Depth");
#elif defined(HAVE_OPENCV_HIGHGUI)
d1 = new vpDisplayOpenCV(I_color, 100, 30, "Pose from Homography");
d2 = new vpDisplayOpenCV(I_color2, I_color.getWidth() + 120, 30, "Pose from RGBD fusion");
d3 = new vpDisplayOpenCV(I_depth, 100, I_color.getHeight() + 70, "Depth");
#endif
}
//! [Create AprilTag detector]
vpDetectorAprilTag detector(tagFamily);
//! [Create AprilTag detector]
//! [AprilTag detector settings]
detector.setAprilTagQuadDecimate(quad_decimate);
detector.setAprilTagPoseEstimationMethod(poseEstimationMethod);
detector.setAprilTagNbThreads(nThreads);
detector.setDisplayTag(display_tag, color_id < 0 ? vpColor::none : vpColor::getColor(color_id), thickness);
detector.setZAlignedWithCameraAxis(align_frame);
//! [AprilTag detector settings]
std::vector<double> time_vec;
for (;;) {
double t = vpTime::measureTimeMs();
//! [Acquisition]
g.acquire(reinterpret_cast<unsigned char *>(I_color.bitmap),
reinterpret_cast<unsigned char *>(I_depth_raw.bitmap), nullptr, nullptr, &align_to_color);
//! [Acquisition]
I_color2 = I_color;
vpImageConvert::convert(I_color, I);
vpImageConvert::createDepthHistogram(I_depth_raw, I_depth);
depthMap.resize(I_depth_raw.getHeight(), I_depth_raw.getWidth());
#ifdef VISP_HAVE_OPENMP
#pragma omp parallel for
#endif
for (int i = 0; i < static_cast<int>(I_depth_raw.getHeight()); i++) {
for (int j = 0; j < static_cast<int>(I_depth_raw.getWidth()); j++) {
if (I_depth_raw[i][j]) {
float Z = I_depth_raw[i][j] * depth_scale;
depthMap[i][j] = Z;
}
else {
depthMap[i][j] = 0;
}
}
}
vpDisplay::display(I_color);
vpDisplay::display(I_color2);
vpDisplay::display(I_depth);
std::vector<vpHomogeneousMatrix> cMo_vec;
detector.detect(I, tagSize, cam, cMo_vec);
// Display camera pose for each tag
for (size_t i = 0; i < cMo_vec.size(); i++) {
vpDisplay::displayFrame(I_color, cMo_vec[i], cam, tagSize / 2, vpColor::none, 3);
}
//! [Pose from depth map]
std::vector<std::vector<vpImagePoint> > tags_corners = detector.getPolygon();
std::vector<int> tags_id = detector.getTagsId();
std::map<int, double> tags_size;
tags_size[-1] = tagSize; // Default tag size
std::vector<std::vector<vpPoint> > tags_points3d = detector.getTagsPoints3D(tags_id, tags_size);
for (size_t i = 0; i < tags_corners.size(); i++) {
vpHomogeneousMatrix cMo;
double confidence_index;
if (vpPose::computePlanarObjectPoseFromRGBD(depthMap, tags_corners[i], cam, tags_points3d[i], cMo,
&confidence_index)) {
if (confidence_index > 0.5) {
vpDisplay::displayFrame(I_color2, cMo, cam, tagSize / 2, vpColor::none, 3);
}
else if (confidence_index > 0.25) {
vpDisplay::displayFrame(I_color2, cMo, cam, tagSize / 2, vpColor::orange, 3);
}
else {
vpDisplay::displayFrame(I_color2, cMo, cam, tagSize / 2, vpColor::red, 3);
}
std::stringstream ss;
ss << "Tag id " << tags_id[i] << " confidence: " << confidence_index;
vpDisplay::displayText(I_color2, 35 + i * 15, 20, ss.str(), vpColor::red);
if (opt_verbose) {
std::cout << "cMo[" << i << "] orig: \n" << cMo_vec[i] << std::endl;
std::cout << "cMo[" << i << "] from RGBD: \n" << cMo << std::endl;
}
}
}
//! [Pose from depth map]
vpDisplay::displayText(I_color, 20, 20, "Pose from homography + VVS", vpColor::red);
vpDisplay::displayText(I_color2, 20, 20, "Pose from RGBD fusion", vpColor::red);
vpDisplay::displayText(I_color, 35, 20, "Click to quit.", vpColor::red);
t = vpTime::measureTimeMs() - t;
time_vec.push_back(t);
std::stringstream ss;
ss << "Detection time: " << t << " ms for " << detector.getNbObjects() << " tags";
vpDisplay::displayText(I_color, 50, 20, ss.str(), vpColor::red);
if (vpDisplay::getClick(I_color, false))
break;
vpDisplay::flush(I_color);
vpDisplay::flush(I_color2);
vpDisplay::flush(I_depth);
}
std::cout << "Benchmark loop processing time" << std::endl;
std::cout << "Mean / Median / Std: " << vpMath::getMean(time_vec) << " ms"
<< " ; " << vpMath::getMedian(time_vec) << " ms"
<< " ; " << vpMath::getStdev(time_vec) << " ms" << std::endl;
if (!display_off) {
delete d1;
delete d2;
delete d3;
}
}
catch (const vpException &e) {
std::cerr << "Catch an exception: " << e.getMessage() << std::endl;
}
return EXIT_SUCCESS;
#else
(void)argc;
(void)argv;
#ifndef VISP_HAVE_APRILTAG
std::cout << "Enable Apriltag support, configure and build ViSP to run this tutorial" << std::endl;
#else
std::cout << "Install librealsense 3rd party, configure and build ViSP again to use this example" << std::endl;
#endif
#endif
return EXIT_SUCCESS;
}
On my side I have printings like:
$ ./tutorial-apriltag-detector-live-rgbd-realsense
Use Realsense 2 grabber
I_color: 640 480
I_depth_raw: 640 480
Read camera parameters from Realsense device
Camera parameters for perspective projection without distortion:
px = 608.263 py = 608.282
u0 = 324.869 v0 = 246.946
poseEstimationMethod: HOMOGRAPHY_VIRTUAL_VS
tagFamily: 36h11
nThreads : 1
Z aligned: 0
cMo[0] orig:
0.9133185316 -0.1806541258 0.3649840362 0.5132143997
-0.1242392292 -0.9771005989 -0.1727397854 0.2027256606
0.3878322752 0.1124211118 -0.9148484136 1.215087718
0 0 0 1
cMo[0] from RGBD:
0.9133181784 -0.1806540897 0.3649849379 0.5132139584
-0.1242391623 -0.977100645 -0.1727395728 0.2027255573
0.3878331285 0.1124207691 -0.914848094 1.215086773
0 0 0 1
cMo[1] orig:
0.7510274501 -0.1860919048 0.6335042006 -0.001361306844
-0.06464711833 -0.9755756599 -0.2099354231 0.1222851822
0.6570985613 0.1167130445 -0.7447144057 0.6585298159
0 0 0 1
cMo[1] from RGBD:
0.7510284957 -0.186091707 0.6335030192 -0.00136135469
-0.06464784282 -0.9755757926 -0.2099345835 0.1222852265
0.657097295 0.1167122508 -0.7447156473 0.658530116
0 0 0 1
What do you have on your side ?
@yangqingnian
Could be some linear algebra issue? You should record all the necessary data + the intrinsic matrix in order to be able to reproduce your issue.
thanks,i try your code, and get the follow result cMo[0] orig: 0.7869443847 0.2849374236 0.5472926092 0.05809407295 0.3143618831 -0.9483850557 0.04174197735 -0.1897006697 0.5309379831 0.1391993205 -0.8358997591 2.258655847 0 0 0 1 cMo[0] from RGBD: 0.7024851961 0.3033318084 0.6438201327 1.269809057e+62 0.240923481 -0.9525698558 0.1859208063 -4.156847072e+62 0.6696793453 0.02450477338 -0.7422459771 4.948832852e+63 0 0 0 1 the result from rgbd maybe also wrong
when i run the program, it also show [ INFO:0@1.343] global registry_parallel.impl.hpp:96 cv::parallel::ParallelBackendRegistry::ParallelBackendRegistry core(parallel): Enabled backends(3, sorted by priority): ONETBB(1000); TBB(990); OPENMP(980) [ INFO:0@1.343] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load F:\visp-ws\3rdparty\opencv-4.8.0\build\x64\vc17\bin\opencv_core_parallel_onetbb480_64d.dll => FAILED [ INFO:0@1.346] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load opencv_core_parallel_onetbb480_64d.dll => FAILED [ INFO:0@1.347] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load F:\visp-ws\3rdparty\opencv-4.8.0\build\x64\vc17\bin\opencv_core_parallel_tbb480_64d.dll => FAILED [ INFO:0@1.351] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load opencv_core_parallel_tbb480_64d.dll => FAILED [ INFO:0@1.351] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load F:\visp-ws\3rdparty\opencv-4.8.0\build\x64\vc17\bin\opencv_core_parallel_openmp480_64d.dll => FAILED [ INFO:0@1.359] global plugin_loader.impl.hpp:67 cv::plugin::impl::DynamicLib::libraryLoad load opencv_core_parallel_openmp480_64d.dll => FAILED
why the matrix from orig and rgb is different
Because the underlying algorithm uses depth info in the second case.
What is the full content of your ViSP-third-party.txt
file
ViSP-third-party.txt this is my visp-third-party.txt
Maybe you can try to record the data when you have the issue with:
when i use tutorial-apriltag-detector-live-rgbd-realsense.cpp to detect Apriltag,vpPose::computePlanarObjectPoseFromRGBD() maybe get a wrong result, like -0.01018921118 -0.9994735211 -0.03080358246 -1.138677748e+63 -0.9994234887 0.009181315827 0.03268629311 7.796915965e+62 -0.03238626704 0.03111887139 -0.9989908636 1.84748416e+64 0 0 0 1