Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034 #ifndef DEPTH_IMAGE_PROC_DEPTH_CONVERSIONS
00035 #define DEPTH_IMAGE_PROC_DEPTH_CONVERSIONS
00036
00037 #include <sensor_msgs/Image.h>
00038 #include <pcl/point_types.h>
00039 #include <image_geometry/pinhole_camera_model.h>
00040 #include <depth_image_proc/depth_traits.h>
00041
00042 #include <limits>
00043
00044 namespace depth_image_proc {
00045
00046
00047 template<typename T>
00048 void convert(
00049 const sensor_msgs::ImageConstPtr& depth_msg,
00050 pcl::PointCloud<pcl::PointXYZ>::Ptr& cloud_msg,
00051 const image_geometry::PinholeCameraModel& model,
00052 double range_max = 0.0)
00053 {
00054
00055 float center_x = model.cx();
00056 float center_y = model.cy();
00057
00058
00059 double unit_scaling = DepthTraits<T>::toMeters( T(1) );
00060 float constant_x = unit_scaling / model.fx();
00061 float constant_y = unit_scaling / model.fy();
00062 float bad_point = std::numeric_limits<float>::quiet_NaN();
00063
00064 pcl::PointCloud<pcl::PointXYZ>::iterator pt_iter = cloud_msg->begin();
00065 const T* depth_row = reinterpret_cast<const T*>(&depth_msg->data[0]);
00066 int row_step = depth_msg->step / sizeof(T);
00067 for (int v = 0; v < (int)cloud_msg->height; ++v, depth_row += row_step)
00068 {
00069 for (int u = 0; u < (int)cloud_msg->width; ++u)
00070 {
00071 pcl::PointXYZ& pt = *pt_iter++;
00072 T depth = depth_row[u];
00073
00074
00075 if (!DepthTraits<T>::valid(depth))
00076 {
00077 if (range_max != 0.0)
00078 {
00079 depth = DepthTraits<T>::fromMeters(range_max);
00080 }
00081 else
00082 {
00083 pt.x = pt.y = pt.z = bad_point;
00084 continue;
00085 }
00086 }
00087
00088
00089 pt.x = (u - center_x) * depth * constant_x;
00090 pt.y = (v - center_y) * depth * constant_y;
00091 pt.z = DepthTraits<T>::toMeters(depth);
00092 }
00093 }
00094 }
00095
00096 }
00097
00098 #endif