OccupancyGridMap.cpp

Go to the documentation of this file.

/*+-------------------------------------------------------------------------+
  |                       MultiVehicle simulator (libmvsim)                 |
  |                                                                         |
  | Copyright (C) 2014  Jose Luis Blanco Claraco (University of Almeria)    |
  | Copyright (C) 2017  Borys Tymchenko (Odessa Polytechnic University)     |
  | Distributed under GNU General Public License version 3                  |
  |   See <http://www.gnu.org/licenses/>                                    |
  +-------------------------------------------------------------------------+ */

#include <mvsim/WorldElements/OccupancyGridMap.h>
#include <mvsim/World.h>
#include "xml_utils.h"

#include <mrpt/opengl/COpenGLScene.h>
#include <mrpt/poses/CPose2D.h>
#include <mrpt/utils/CFileGZInputStream.h>
#include <mrpt/system/filesystem.h>
#include <mrpt/utils/CObject.h>
#include <rapidxml.hpp>

#include <mrpt/version.h>
#if MRPT_VERSION >= 0x130
#include <mrpt/maps/CSimplePointsMap.h>
using mrpt::maps::CSimplePointsMap;
#else
#include <mrpt/slam/CSimplePointsMap.h>
using mrpt::slam::CSimplePointsMap;
#endif

using namespace rapidxml;
using namespace mvsim;
using namespace std;

OccupancyGridMap::OccupancyGridMap(
        World* parent, const rapidxml::xml_node<char>* root)
        : WorldElementBase(parent),
          m_gui_uptodate(false),
          m_show_grid_collision_points(true),
          m_restitution(0.01),
          m_lateral_friction(0.5)
{
        loadConfigFrom(root);
}

OccupancyGridMap::~OccupancyGridMap() {}
void OccupancyGridMap::loadConfigFrom(const rapidxml::xml_node<char>* root)
{
        m_gui_uptodate = false;

        // <file>FILENAME.{png,gridmap}</file>
        xml_node<>* xml_file = root->first_node("file");
        if (!xml_file || !xml_file->value())
                throw std::runtime_error(
                        "Error: <file></file> XML entry not found inside gridmap node!");

        const string sFile = m_world->resolvePath(xml_file->value());
        const string sFileExt =
                mrpt::system::extractFileExtension(sFile, true /*ignore gz*/);

        // MRPT gridmaps format:
        if (sFileExt == "gridmap")
        {
                mrpt::utils::CFileGZInputStream f(sFile);
                f >> m_grid;
        }
        else
        // Assume it's an image:
        {
                std::map<std::string, TParamEntry> other_params;
                double xcenterpixel = -1, ycenterpixel = -1;
                double resolution = 0.10;

                other_params["resolution"] = TParamEntry("%lf", &resolution);
                other_params["centerpixel_x"] = TParamEntry("%lf", &xcenterpixel);
                other_params["centerpixel_y"] = TParamEntry("%lf", &ycenterpixel);

                parse_xmlnode_children_as_param(*root, other_params);

                if (!m_grid.loadFromBitmapFile(
                                sFile, resolution, xcenterpixel, ycenterpixel))
                        throw std::runtime_error(
                                mrpt::format(
                                        "[OccupancyGridMap] ERROR: File not found '%s'",
                                        sFile.c_str()));
        }

        {
                // Other general params:
                std::map<std::string, TParamEntry> ps;
                ps["show_collisions"] =
                        TParamEntry("%bool", &m_show_grid_collision_points);
                ps["restitution"] = TParamEntry("%lf", &m_restitution);
                ps["lateral_friction"] = TParamEntry("%lf", &m_lateral_friction);

                parse_xmlnode_children_as_param(*root, ps);
        }
}

void OccupancyGridMap::gui_update(mrpt::opengl::COpenGLScene& scene)
{
        using namespace mrpt::math;

        // 1st time call?? -> Create objects
        if (!m_gl_grid)
        {
                m_gl_grid = mrpt::make_aligned_shared<mrpt::opengl::CSetOfObjects>();
                SCENE_INSERT_Z_ORDER(scene, 0, m_gl_grid);
        }
        if (m_gl_obs_clouds.size() != m_obstacles_for_each_obj.size())
        {
                m_gl_obs_clouds.resize(m_obstacles_for_each_obj.size());
        }

        // 1st call OR gridmap changed?
        if (!m_gui_uptodate)
        {
                m_grid.getAs3DObject(m_gl_grid);
                m_gui_uptodate = true;
        }

        // Update obstacles:
        {
                std::lock_guard<std::mutex> csl(m_gl_obs_clouds_buffer_cs);
                for (size_t i = 0; i < m_gl_obs_clouds.size(); i++)
                {
                        mrpt::opengl::CSetOfObjects::Ptr& gl_objs = m_gl_obs_clouds[i];
                        if (!gl_objs)
                        {
                                gl_objs =
                                        mrpt::make_aligned_shared<mrpt::opengl::CSetOfObjects>();
                                MRPT_TODO("Add a name, and remove old ones in scene, etc.")
                                SCENE_INSERT_Z_ORDER(scene, 1, gl_objs);
                        }

                        // Now that we are in a safe thread (with the OpenGL scene lock
                        // adquired from the caller)
                        // proceed to replace the old with the new point cloud:
                        gl_objs->clear();
                        if (m_gl_obs_clouds_buffer.size() > i)
                                gl_objs->insert(m_gl_obs_clouds_buffer[i]);
                }

                m_gl_obs_clouds_buffer.clear();
        }  // end lock
}

void OccupancyGridMap::simul_pre_timestep(const TSimulContext& context)
{
        // Make a list of objects subject to collide with the occupancy grid:
        // - Vehicles
        // - Blocks
        {
                const World::TListVehicles& lstVehs =
                        this->m_world->getListOfVehicles();
                const World::TListBlocks& lstBlocks = this->m_world->getListOfBlocks();
                const size_t nObjs = lstVehs.size() + lstBlocks.size();

                m_obstacles_for_each_obj.resize(nObjs);
                size_t obj_idx = 0;
                for (World::TListVehicles::const_iterator itVeh = lstVehs.begin();
                         itVeh != lstVehs.end(); ++itVeh, ++obj_idx)
                {
                        TInfoPerCollidableobj& ipv = m_obstacles_for_each_obj[obj_idx];
                        ipv.max_obstacles_ranges =
                                itVeh->second->getMaxVehicleRadius() * 1.50f;
                        const mrpt::math::TPose3D& pose = itVeh->second->getPose();
                        ipv.pose = mrpt::poses::CPose2D(
                                pose.x, pose.y,
                                0 /* angle=0, no need to rotate everything to later rotate back again! */);
                }
                for (World::TListBlocks::const_iterator it = lstBlocks.begin();
                         it != lstBlocks.end(); ++it, ++obj_idx)
                {
                        TInfoPerCollidableobj& ipv = m_obstacles_for_each_obj[obj_idx];
                        ipv.max_obstacles_ranges = it->second->getMaxBlockRadius() * 1.50f;
                        const mrpt::math::TPose3D& pose = it->second->getPose();
                        ipv.pose = mrpt::poses::CPose2D(
                                pose.x, pose.y,
                                0 /* angle=0, no need to rotate everything to later rotate back again! */);
                }
        }

        // For each object, create a ground body with fixtures at each scanned point
        // around the vehicle, so it can collide with the environment:
        // Upon first usage, reserve mem:
        {  // lock
                std::lock_guard<std::mutex> csl(m_gl_obs_clouds_buffer_cs);
                const size_t nObjs = m_obstacles_for_each_obj.size();
                m_gl_obs_clouds_buffer.resize(nObjs);

                for (size_t obj_idx = 0; obj_idx < nObjs; obj_idx++)
                {
                        // 1) Simulate scan to get obstacles around the vehicle:
                        TInfoPerCollidableobj& ipv = m_obstacles_for_each_obj[obj_idx];
                        CObservation2DRangeScan::Ptr& scan = ipv.scan;
                        // Upon first time, reserve mem:
                        if (!scan)
                                scan = mrpt::make_aligned_shared<CObservation2DRangeScan>();

                        const float veh_max_obstacles_ranges = ipv.max_obstacles_ranges;
                        const float occup_threshold = 0.5f;
                        const size_t nRays = 50;

                        scan->aperture = 2.0 * M_PI;  // 360 field of view
                        scan->maxRange = veh_max_obstacles_ranges;

                        m_grid.laserScanSimulator(
                                *scan, ipv.pose, occup_threshold, nRays, 0.0f /*noise*/);

                        // Since we'll dilate obstacle points, let's give a bit more space
                        // as compensation:
                        const float range_enlarge = 0.25f * m_grid.getResolution();
                        for (size_t k = 0; k < scan->scan.size(); k++)
                        {
#if MRPT_VERSION >= 0x150
                                scan->setScanRange(k, scan->getScanRange(k) + range_enlarge);
#else
                                scan->scan[k] += range_enlarge;
#endif
                        }
                        // 2) Create a Box2D "ground body" with square "fixtures" so the
                        // vehicle can collide with the occ. grid:
                        b2World* b2world = m_world->getBox2DWorld();

                        // Create Box2D objects upon first usage:
                        if (!ipv.collide_body)
                        {
                                b2BodyDef bdef;
                                ipv.collide_body = b2world->CreateBody(&bdef);
                                ASSERT_(ipv.collide_body)
                        }
                        // Force move the body to the vehicle origins (to use obstacles in
                        // local coords):
                        ipv.collide_body->SetTransform(
                                b2Vec2(ipv.pose.x(), ipv.pose.y()), .0f);

                        // GL:
                        // 1st usage?
                        mrpt::opengl::CPointCloud::Ptr& gl_pts =
                                m_gl_obs_clouds_buffer[obj_idx];
                        if (m_show_grid_collision_points)
                        {
                                gl_pts = mrpt::make_aligned_shared<mrpt::opengl::CPointCloud>();
                                gl_pts->setPointSize(4.0f);
                                gl_pts->setColor(0, 0, 1);

                                gl_pts->setVisibility(m_show_grid_collision_points);
                                gl_pts->setPose(mrpt::poses::CPose2D(ipv.pose));
                                gl_pts->clear();
                        }

                        // Physical properties of each "occupied cell":
                        const float occCellSemiWidth = m_grid.getResolution() * 0.4f;
                        b2PolygonShape sqrPoly;
                        sqrPoly.SetAsBox(occCellSemiWidth, occCellSemiWidth);
                        sqrPoly.m_radius = 1e-3;  // The "skin" depth of the body
                        b2FixtureDef fixtureDef;
                        fixtureDef.shape = &sqrPoly;
                        fixtureDef.restitution = m_restitution;
                        fixtureDef.density = 0;  // Fixed (inf. mass)
                        fixtureDef.friction = m_lateral_friction;  // 0.5f;

                        // Create fixtures at their place (or disable it if no obstacle has
                        // been sensed):
                        const CSinCosLookUpTableFor2DScans::TSinCosValues& sincos_tab =
                                m_sincos_lut.getSinCosForScan(*scan);
                        ipv.collide_fixtures.resize(nRays);
                        for (size_t k = 0; k < nRays; k++)
                        {
                                if (!ipv.collide_fixtures[k].fixture)
                                        ipv.collide_fixtures[k].fixture =
                                                ipv.collide_body->CreateFixture(&fixtureDef);

                                if (!scan->validRange[k])
                                {
                                        ipv.collide_fixtures[k].fixture->SetSensor(
                                                true);  // Box2D's way of saying: don't collide with
                                                                // this!
                                        ipv.collide_fixtures[k].fixture->SetUserData(
                                                INVISIBLE_FIXTURE_USER_DATA);
                                }
                                else
                                {
                                        ipv.collide_fixtures[k].fixture->SetSensor(
                                                false);  // Box2D's way of saying: don't collide with
                                                                 // this!
                                        ipv.collide_fixtures[k].fixture->SetUserData(NULL);

                                        b2PolygonShape* poly = dynamic_cast<b2PolygonShape*>(
                                                ipv.collide_fixtures[k].fixture->GetShape());
                                        ASSERT_(poly != NULL)

                                        const float llx = sincos_tab.ccos[k] * scan->scan[k];
                                        const float lly = sincos_tab.csin[k] * scan->scan[k];

                                        const float ggx = ipv.pose.x() + llx;
                                        const float ggy = ipv.pose.y() + lly;

                                        const float gx = m_grid.idx2x(m_grid.x2idx(ggx));
                                        const float gy = m_grid.idx2y(m_grid.y2idx(ggy));

                                        const float lx = gx - ipv.pose.x();
                                        const float ly = gy - ipv.pose.y();

                                        poly->SetAsBox(
                                                occCellSemiWidth, occCellSemiWidth, b2Vec2(lx, ly),
                                                .0f /*angle*/);

                                        if (gl_pts && m_show_grid_collision_points)
                                        {
                                                gl_pts->mrpt::opengl::CPointCloud::insertPoint(
                                                        lx, ly, .0);
                                        }
                                }
                        }
                }  // end for obj_idx

        }  // end lock
}