雷达数据的运动畸变

在移动机器人的运动过程中,当激光雷达扫描频率较低时,获得的原始扫描数据会发生畸变,并且随着速度的增加变得越来越严重。雷达旋转有顺时针与逆时针2种情况,但导致运动畸变的效果差异并不明显。

  1. 激光点数据不是瞬时获得的,雷达扫描时在旋转,有个扫描频率的问题
  2. 雷达扫描中,机器人也是运动的。 机器人旋转造成的运动畸变比直线运动大。

比如所用激光雷达的扫描频率为10Hz,当移动机器人以角速度为1 rad/s 进行纯旋转时,在激光雷达扫描一圈得到
一帧数据的时间内(0.1s)移动机器人已经旋转了0.1 rad,约为5.7度,这会影响机器人的定位和路径规划。

由于我以前的机器人是低速机器人,最大线速度才0.4m/s,而雷达的扫描频率至少25Hz,所以雷达不必考虑畸变,实际使用时,从rviz上看,也没发现运动畸变。雷达扫描频率最好达到30Hz以上。但是我现在使用的机器人速度最快1.4m/s,雷达频率才10Hz,这就需要考虑运动畸变了。

使用里程计辅助的方法去除运动畸变。里程计的频率一般比较高,move_base里面默认订阅里程计的频率是20Hz,所以里程计至少也得20Hz,我所用的里程计是50Hz。

平时之所以没有这么严重,是因为有定位程序,比如AMCL或者LOAM

去除激光帧畸变原理:把一帧激光的每个激光点坐标变换到不同时刻的机器人里程计上
里程计辅助方法:imu、轮式里程计
(1) IMU
  直接测量获得角速度和线速度,但对于机器人位移和角度需要积分
  更新速度高 1KHz-8KHz
  具有较高的角速度测量精度,但线加速度精度差
2)轮式里程计
  更新速度高 100-200Hz
  (一般机器人轮式里程计200hz,即可认为满足一定的要求,200hz一帧5ms,1.57rad/s下,误差仅为0.45度,可认为机器人没有运动。)

程序流程

1) main-构造函数
订阅scan话题
调用scancallback处理
2) scancallback
scan话题数据提取到vector数组(包括angles,ranges)
vector数组保存到pcl::PointCloud类型的点云中
调用Lidar_Calibration进行畸变矫正
3) Lidar_Calibration
程序对激光雷达数据进行分段线性插值*
调用Lidar_MotionCalibration激光雷达运动畸变去除分段函数*
调用getLaserPose得到激光雷达在odom中的位姿

参考:
基于里程计修复激光雷达数据
激光雷达移动状态下的数据矫正
激光雷达数据效果对比


马氏距离

马氏距离

数据的协方差距离。它是一种有效的计算两个未知样本集的相似度的方法。与欧氏距离不同的是它考虑到各种特性之间的联系(例如:一条关于身高的信息会带来一条关于体重的信息,因为两者是有关联的)并且是尺度无关的(scale-invariant),即独立于测量尺度。

均方根误差

RMSE = Root Mean Square Error=均方根误差。

含义:衡量观测值与真实值之间的偏差。
用途:常用来作为机器学习模型预测结果衡量的标准
RMSE.png


使用 CMake 解决编译出来的so文件过大的问题

CMake加入的一些设置,简化编译结果

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
set(SIMPLIFY "ON")

if(SIMPLIFY STREQUAL "ON")
MESSAGE("catkin_make small libfile.so")
SET(CMAKE_BUILD_TYPE "Release")
# 删除调试符号
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -s")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -s")

# 开启空间优化
if (APPLE)
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -Oz")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -Oz")
else ()
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -Os")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -Os")
endif ()

# 使用gc-section优化
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -ffunction-sections -fdata-sections")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -ffunction-sections -fdata-sections")
set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} -Wl,--gc-sections")
set(CMAKE_MODULE_LINKER_FLAGS "${CMAKE_MODULE_LINKER_FLAGS} -Wl,--gc-sections")
set(CMAKE_SHARED_LINKER_FLAGS "${CMAKE_SHARED_LINKER_FLAGS} -Wl,--gc-sections")
endif()

  1. 编译选项使用-Os,表示以最小化大小为优化方向,假如链接失败,则修改为-O2

在 CMakeLists.txt 文件中加入以下两行:

1
2
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -s")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -s -O4 -Wall")

  1. 去除-g选项,进而去除调试信息

  2. 通过strip裁剪符号及调试信息

  3. 只导出必要符号

Linux默认导出所有符号,并不仅仅导出你开放的接口。

定义如下宏:

#define SYMBOL_EXPORT __attribute__ ((visibility("default")))

在你想导出的符号前加上SYMBOL_EXPORT就好。

同时,需要在脚本中增加如下编译选项: -fvisibility = hidden

譬如:导出符号是int add(int a, int b); 那么添加的结果就是SYMBOL_EXPORT int add(int a, int b);

通过这种方式只导出想导出的符号,也可以减小库大小。

default: 用它定义的符号将被导出,动态库中的函数默认是可见的。

hidden: 用它定义的符号将不被导出,并且不能从其它对象进行使用,动态库中的函数是被隐藏的。

default意味着该方法对其它模块是可见的。而hidden表示该方法符号不会被放到动态符号表里,所以其它模块(可执行文件或者动态库)不可以通过符号表访问该方法。这个选项可能导致编译失败,出现undef reference

还有一个选项 -fvisibility = inlines-hidden

  • 控制C++特性的使用,如无必要,则不使用C++的高级特性。

  • 屏蔽RTTI特性,增加编译选项-fno-rtti。

  • 对性能影响不大时,避免使用C++的inline特性,也就是增加编译选项 -fno-inline

  • 在不影响使用时,控制对STL的组件的使用。尤其iostream的相关模板类。

示例:

1
2
3
4
5
6
7
8
9
10
11
12
# 关闭C++特性set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -fno-inline")
set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -fno-exceptions")
set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -fno-rtti")

set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -fvisibility = inlines-hidden")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -fvisibility = inlines-hidden")
# 删除调试符号
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -s")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -s")
# 开启空间优化
set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -Os")
set(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS_RELEASE} -Os")


tf(四) 常用的程序

使用程序实现父子坐标系的转换

除了直接用static_transform_publisher外,当然可以在程序中实现两个坐标系的转换,以经典的乌龟追乌龟程序为参考,它的tf树是这样的

现在我们写一个这样的程序:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
#include <ros/ros.h>
#include <tf/transform_broadcaster.h>
#include <tf/tf.h>
#include <tf/transform_listener.h>
#include <turtlesim/Pose.h>

void poseCallback(const turtlesim::PosePtr& msg)
{
static tf::TransformBroadcaster br;
tf::Transform transform;
transform.setOrigin( tf::Vector3(0,0, -0.24) ); //子坐标系相对父坐标系的位移
//子坐标系相对父坐标系的旋转,可以用欧拉角和四元数两种形式
tf::Quaternion q;
q.setRPY(0, 0, 0);
//tf::Quaternion q(qx,qy,qz,qw);
transform.setRotation(q);
br.sendTransform(tf::StampedTransform(transform, ros::Time::now(), "turtle1", "robot_base")); //父坐标系在前,子坐标系在后
}

int main(int argc, char** argv)
{
ros::init(argc, argv, "camera_base_tf");
ros::NodeHandle node;
ros::Subscriber sub = node.subscribe("/turtle1/pose", 10, &poseCallback);
ros::spin();
return 0;
}

我们订阅乌龟1的位姿话题/turtle1/pose,在回调函数里定义tf广播和变换的参数,最后发送。这里只能用回调函数的形式,因为在main函数里只能发一次,无法建立父子坐标系的关系。tf::broadcaster类的sendTransform函数中,把父子坐标系交换,只会交换tf树中的位置,坐标关系还是原来的,所以必须使用Transform::inverse ()函数.

使用程序获得坐标系在另一坐标系中的坐标

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#include <ros/ros.h>
#include <tf/tf.h>
#include <tf/transform_listener.h>

int main(int argc, char** argv)
{
ros::init(argc, argv, "turtle_2_1");
ros::NodeHandle n;
tf::TransformListener listener;
ros::Rate rate(10);
while(n.ok())
{
tf::StampedTransform transform;
try
{
listener.waitForTransform("turtle1","turtle2",ros::Time(0),ros::Duration(3));
listener.lookupTransform("turtle1","turtle2",ros::Time(0), transform);
}
catch(tf::TransformException &ex)
{
ROS_ERROR("%s",ex.what());
ros::Duration(1).sleep();
continue;
}
ROS_INFO("x: %f", transform.getOrigin().x());
ROS_INFO("y: %f", transform.getOrigin().y());
ROS_INFO("z: %f\n", transform.getOrigin().z());
rate.sleep();
}
return 0;
}

这个程序获得的是乌龟2在乌龟1坐标系中的坐标

先运行roslaunch turtle_tf turtle_tf_demo.launch,这是ROS自带的一个tf演示程序,然后马上运行上面的程序。我们会看到乌龟1不动,乌龟2向前者移动,如果运行程序太慢,乌龟2就停止了,二者没有相对运动就捕捉不到tf话题上的消息了


同时我们的程序会输出下面内容

1
2
3
4
5
6
7
8
9
10
11
12
[ INFO] [1551530137.681192434]: x: -1.394780
[ INFO] [1551530137.681317271]: y: -3.510608
[ INFO] [1551530137.681352733]: z: 0.000000

[ INFO] [1551530137.781891361]: x: -1.221200
[ INFO] [1551530137.781970030]: y: -3.391552
[ INFO] [1551530137.782007398]: z: 0.000000

[ INFO] [1551530137.982427762]: x: -1.018356
[ INFO] [1551530137.982512261]: y: -3.122349
[ INFO] [1551530137.982559203]: z: 0.000000
...... //忽略

显然从图中能看出乌龟2在乌龟1的坐标系中初始坐标x,y肯定都是负值。随着运动,xy的绝对值越来越小,最终两个乌龟重合,xy都成为0。

程序分析

waitForTransform函数的原型有两个,我们用的这个是:

1
2
3
4
5
6
7
8
9
bool Transformer::waitForTransform
(
const std::string & target_frame,
const std::string & source_frame,
const ros::Time & time, //The time at which to transform
const ros::Duration & timeout, //How long to block before failing
const ros::Duration & polling_sleep_duration = ros::Duration(0.01),
std::string * error_msg = NULL
)const

功能:函数阻塞直到能进行坐标系转换或者时间timeout截止,返回值表示是否成功。函数不会throw,如果给error_msg传变量,它会在函数出错时记录错误信息。

一般用前四个参数就行了,如果开始没有坐标系转换,会报错 “turtle1” passed to lookupTransform argument target_frame does not exist 。如果又有了转换,程序可以正常运行。因为网络问题,如果timeout为0,那会直接报错失败,如果设置太小,也可能捕捉不到转换。使用tf_echo就会发现它要等一点时间才会输出结果,这是一样的道理。

lookupTransform函数的原型也有两个,我们用的是:

1
2
3
4
5
6
7
void Transformer::lookupTransform
(
const std::string & target_frame,
const std::string & source_frame,
const ros::Time & time, //某时间的变换,0表示最近的变换
StampedTransform & transform
)const

功能:获得两坐标系之间的转换,也就是赋值给tf::StampedTransform对象

还有函数canTransform(),这个函数纯粹就是判断用的,也能记录错误信息,但不能指定阻塞时间。

技巧 1

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
std::string tf_error;
// make sure that the transform between base_link and the global frame is available
while (ros::ok()
&& !tf_.waitForTransform(global_frame_, "base_link", ros::Time(), ros::Duration(0.1), ros::Duration(0.01), &tf_error))
{
ros::spinOnce();
if (last_error + ros::Duration(5.0) < ros::Time::now())
{
ROS_WARN("Timed out waiting for transform from %s to %s to become available before running costmap, tf error: %s",
robot_base_frame_.c_str(), global_frame_.c_str(), tf_error.c_str());
last_error = ros::Time::now();
}
// The error string will accumulate and errors will typically be the same, so the last
// will do for the warning above. Reset the string here to avoid accumulation.
tf_error.clear();
}

其他有用的函数

tf::Transformer类的重要函数

  • std::string Transformer::allFramesAsString() const    A way to see what frames have been cached Useful for debugging.

  • bool Transformer::frameExists (const std::string & frame_id_str) const    Check if a frame exists in the tree.

  • void Transformer::getFrameStrings ( std::vector< std::string > & ids ) const    A way to get a std::vector of available frame ids.

  • getParent

1
2
3
4
bool Transformer::getParent ( const std::string &   frame_id,
ros::Time time,
std::string & parent
) const

Fill the parent of a frame.

  • frame_id The frame id of the frame in question
  • parent The reference to the string to fill the parent Returns true unless “NO_PARENT”
  • transformQuaternion()

欧拉角和四元数

概念

欧拉角其实就是表示向三个方向旋转的角度,比较形象直观。四元数是等价的另一种表示方法,不能直观体现,但是具备一些优点。数学理论实在太复杂,只需要记住以下几条:

  • 欧拉角就是绕Z轴、Y轴、X轴的旋转角度,记做ψ、θ、φ,分别为Yaw(偏航)、Pitch(俯仰)、Roll(翻滚)
  • 四元数是w,x,y,z,它们的平方和是1
  • 欧拉角会出现万向节死锁,比如三维空间中有一个平行于 X 轴的向量,将它绕 Y 轴旋转直到它平行于 Z 轴,这时任何绕 Z 轴的旋转都改变不了该向量的方向

三个旋转方向如图所示:

欧拉角到四元数的转换:

四元数到欧拉角的转换:

描述机器人转弯的角就是绕Z轴的旋转,按照右手法则可以知道机器人向左转为正,右转为负。

ROS中的欧拉角 四元数

geometry_msgs/Quaternion消息

1
2
3
4
float64 x
float64 y
float64 z
float64 w

geometry_msgs/Transform消息

1
2
geometry_msgs/Vector3 translation
geometry_msgs/Quaternion rotation

这两个各自有一个加时间戳的消息类型。


除了geometry_msgs/Quaternion,还有一个tf::Quaternion类,重要函数:

1
2
3
4
5
Vector3  getAxis () const    //Return the axis of the rotation represented by this quaternion. 
//Set the quaternion using fixed axis RPY.
void setRPY (const tfScalar &roll, const tfScalar &pitch, const tfScalar &yaw)
//Return an identity quaternion.
tf::Quaternion createIdentityQuaternion()

还有一个有点特殊的函数,直接返回yaw,但是范围是[0, 2Pi],所以尽量不要使用

1
tf::Quaternion::getAngle () const

geometry_msgs::Quaternion —— tf::Quaternion 转换如下
1
2
// [inline, static]
static void tf::quaternionMsgToTF(const geometry_msgs::Quaternion& msg, tf::Quaternion& bt)

四元数——欧拉角

1
2
3
4
tf::Quaternion quat;
double roll, pitch, yaw;
// 赋值给quat,getRPY函数 将quat转为欧拉角
tf::Matrix3x3(quat).getRPY(roll, pitch, yaw);

orientation是表示朝向的消息成员,也就是四元数,类型const geometry_msgs::Quaternion &。 常常配合tf::quaternionMsgToTF函数使用,因为我们一般是先获得geometry_msgs::Quaternion类型

或者用tf2的方式

1
2
3
4
5
6
7
#include “tf2/LinearMath/Matrix3x3.h”

tf2::Matrix3x3 mat(tf2::Quaternion( x, y, z, w) );
mat.getEulerYPR(tf2Scalar& yaw,
tf2Scalar& pitch,
tf2Scalar& roll,
unsigned int solution_number = 1);

欧拉角 —— 四元数

1
2
3
4
5
6
7
8
9
//只通过yaw计算四元数,用于平面小车
static geometry_msgs::Quaternion tf::createQuaternionMsgFromYaw(double yaw);

static tf::Quaternion createQuaternionFromYaw (double yaw)

static tf::Quaternion createQuaternionFromRPY(double roll, double pitch, double yaw)

static geometry_msgs::Quaternion createQuaternionMsgFromRollPitchYaw( double roll,
double pitch, double yaw)

python常见的写法:

1
2
3
pos = Pose()
q = tf.transformations.quaternion_from_euler(0, 0, point.z)
pos.orientation = q

但是经测试,这样写是有问题的,正确的写法如下:
1
2
3
4
5
6
pos = Pose()
q = tf.transformations.quaternion_from_euler(0, 0, point.z)
pos.orientation.x = q[0]
pos.orientation.y = q[1]
pos.orientation.z = q[2]
pos.orientation.w = q[3]

报错积累

报错 MSG to TF: Quaternion Not Properly Normalized , 原因是接受的位姿的四元数不合理。

参考:
四元数与欧拉角(Yaw、Pitch、Roll)的转换


实际导航问题的积累
abstract Welcome to my blog, enter password to read.
Read more
添加新约束:Jerk(加加速度)

加速度的变化率过大会使机器人底盘电机输出的力矩突变而引起机器人受冲击振荡,为了使TEB生成的轨迹 less jerky, 考虑jerk约束。jerky的意思是有突然的停止和行走现象,也就是trajectory不够平滑稳定。

Jerk约束  7元边.png


ROS的global_planner
abstract Welcome to my blog, enter password to read.
Read more
多台车使用gmapping联合建图

三台车使用RRT和gmapping联合建图

  • 集群机器人可以协作进行地图构建吗?
    两个问题:每个机器人的定位;每个机器人和主节点通信。机器人1把扫描的点加到地图,机器人2把扫描的点加到地图,最后形成完整地图。

  • 室内的两个车怎么进行准确定位呢?
    环境不大的话,在天花板上贴二维码,两个机器人都去识别二维码。如果环境很大,初始建图时候,让两个机器人离得近一些,手工量出来二者的相对位置,建图时候在增量定位,这样也可以定位


PCL1.7和1.8版本冲突问题

同事的点云程序必须用PCL-1.8,而ROS默认的是1.7,结果有时出现了版本冲突。

安装了PCL-1.8后,重要的库文件libcostmap_2d.so链接到了pcl-1.8,而不是默认的1.7。此时删除了pcl-1.8之后,再编译会报错:

1
2
3
4
CMake Error at /home/user/catkin_ws/devel/share/costmap_2d/cmake/costmap_2dConfig.cmake:113 (message):
Project 'costmap_2d' specifies '/usr/local/include/pcl-1.8' as an include
dir, which is not found. It does neither exist as an absolute directory
nor in '/home/user/catkin_ws/src/navigation/costmap_2d//usr/local/include/pcl-1.8'.

编译还会找pcl-1.8的路径,这是CMake caches的问题,不是软链接之类的问题,它还会去找旧的路径,所以删除catkin_ws/build下面的catkin_make.cache, Makefile, CMakeCache.txt. 然后重新编译,这样libcostmap_2d.so就不会再去链接1.8了