Program Listing for File static_agent_layer.cpp

Return to documentation for file (cohan_layers/src/static_agent_layer.cpp)

/*********************************************************************
 *
 * Software License Agreement (BSD License)
 *
 * Copyright (c) 2025 LAAS/CNRS
 * All rights reserved.
 *
 *  Redistribution and use in source and binary forms, with or without
 *  modification, are permitted provided that the following conditions
 *  are met:
 *
 *   * Redistributions of source code must retain the above copyright
 *     notice, this list of conditions and the following disclaimer.
 *   * Redistributions 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.
 *   * Neither the name of the institute nor the names of its
 *     contributors may 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
 *  COPYRIGHT OWNER 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.
 *
 * Author: Phani Teja Singamaneni (email:ptsingaman@laas.fr)
 *********************************************************************/

#include <angles/angles.h>
#include <cohan_layers/static_agent_layer.h>
#include <pluginlib/class_list_macros.h>
#include <tf2_eigen/tf2_eigen.h>

PLUGINLIB_EXPORT_CLASS(cohan_layers::StaticAgentLayer, costmap_2d::Layer)

namespace cohan_layers {
void StaticAgentLayer::onInitialize() {
  AgentLayer::onInitialize();
  ros::NodeHandle private_nh("~/" + name_);
  loadRosParamFromNodeHandle(private_nh);
  server_ = new dynamic_reconfigure::Server<AgentStaticLayerConfig>(private_nh);
  f_ = [this](auto&& PH1, auto&& PH2) { configure(std::forward<decltype(PH1)>(PH1), std::forward<decltype(PH2)>(PH2)); };
  server_->setCallback(f_);
}

void StaticAgentLayer::updateBoundsFromAgents(double* min_x, double* min_y, double* max_x, double* max_y) {
  for (const auto& agent : transformed_agents_) {
    *min_x = std::min(*min_x, agent.pose.position.x - radius_);
    *min_y = std::min(*min_y, agent.pose.position.y - radius_);
    *max_x = std::max(*max_x, agent.pose.position.x + radius_);
    *max_y = std::max(*max_y, agent.pose.position.y + radius_);
  }
}

void StaticAgentLayer::loadRosParamFromNodeHandle(const ros::NodeHandle& private_nh) {
  private_nh.param("robot_radius", robot_radius_, 0.47);
  private_nh.param("agent_radius", agent_radius_, 0.30);
}

void StaticAgentLayer::updateCosts(costmap_2d::Costmap2D& /*master_grid*/, int min_i, int min_j, int max_i, int max_j) {
  boost::recursive_mutex::scoped_lock lock(lock_);
  if (!enabled_) {
    return;
  }

  if (agents_.agents.size() == 0) {
    return;
  }

  costmap_2d::Costmap2D* costmap = layered_costmap_->getCostmap();
  double res = costmap->getResolution();

  for (uint i = 0; i < transformed_agents_.size(); i++) {
    auto agent = transformed_agents_[i];
    auto state = transformed_agents_[i].state;
    auto type = transformed_agents_[i].type;

    // Check the condition to switch the radius
    bool is_human_still = states_.empty() ? (type == 1) : ((state < 2) && (type == 1));

    unsigned int width = std::max(1, static_cast<int>((2 * radius_) / res));
    unsigned int height = std::max(1, static_cast<int>((2 * radius_) / res));

    double cx = agent.pose.position.x;
    double cy = agent.pose.position.y;
    double ox = cx - radius_;
    double oy = cy - radius_;

    int mx;
    int my;
    costmap->worldToMapNoBounds(ox, oy, mx, my);

    int start_x = 0;
    int start_y = 0;
    int end_x = width;
    int end_y = height;

    if (mx < 0) {
      start_x = -mx;
    } else if (mx + width > costmap->getSizeInCellsX()) {
      end_x = std::max(0, static_cast<int>(costmap->getSizeInCellsX()) - mx);
    }

    if ((start_x + mx) < min_i) {
      start_x = min_i - mx;
    }

    if ((end_x + mx) > max_i) {
      end_x = max_i - mx;
    }

    if (my < 0) {
      start_y = -my;
    } else if (my + height > costmap->getSizeInCellsY()) {
      end_y = std::max(0, static_cast<int>(costmap->getSizeInCellsY()) - my);
    }

    if ((start_y + my) < min_j) {
      start_y = min_j - my;
    }

    if ((end_y + my) > max_j) {
      end_y = max_j - my;
    }

    double bx = 0;
    double by = 0;
    double var = 0;
    double rad = 0;

    if (is_human_still) {
      bx = ox + (res / 2);
      by = oy + (res / 2);
      var = radius_;
    }

    else {
      if (type == 1) {
        rad = agent_radius_;
      } else {
        rad = robot_radius_;
      }
    }

    for (int i = start_x; i < end_x; i++) {
      for (int j = start_y; j < end_y; j++) {
        unsigned char old_cost = costmap->getCost(i + mx, j + my);
        if (old_cost == costmap_2d::NO_INFORMATION) {
          continue;
        }

        if (is_human_still) {
          double x = bx + (i * res);
          double y = by + (j * res);
          double val;
          val = Gaussian2D(x, y, cx, cy, amplitude_, var, var);
          double rad = sqrt(-2 * var * log(val / amplitude_));

          if (rad > radius_) {
            continue;
          }
          auto cvalue = static_cast<unsigned char>(val);
          costmap->setCost(i + mx, j + my, std::max(cvalue, old_cost));

        } else {
          double x;
          double y;
          costmap->mapToWorld(i + mx, j + my, x, y);
          x = x - cx, y = y - cy;
          double inrad = sqrt((x * x) + (y * y));

          if (inrad > rad) {
            continue;
          }
          auto cvalue = costmap_2d::LETHAL_OBSTACLE;
          costmap->setCost(i + mx, j + my, std::max(cvalue, old_cost));
        }
      }
    }
  }
}

void StaticAgentLayer::configure(AgentStaticLayerConfig& config, uint32_t /*level*/) {
  amplitude_ = config.amplitude;
  radius_ = config.radius;
  enabled_ = config.enabled;
}
};  // namespace cohan_layers