位置: IT常识 - 正文

KITTI数据集可视化(一):点云多种视图的可视化实现

编辑:rootadmin
KITTI数据集可视化(一):点云多种视图的可视化实现

推荐整理分享KITTI数据集可视化(一):点云多种视图的可视化实现,希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:,内容如对您有帮助,希望把文章链接给更多的朋友!

如有错误,恳请指出。

在本地上,可以安装一些软件,比如:Meshlab,CloudCompare等3D查看工具来对点云进行可视化。而这篇博客是将介绍一些代码工具将KITTI数据集进行可视化操作,包括点云鸟瞰图,FOV图,以及标注信息在图像+点云上的显示。

文章目录1. 数据集准备2. 环境准备3. KITTI数据集可视化4. 点云可视化5. 鸟瞰图可视化1. 数据集准备

KITTI数据集作为自动驾驶领域的经典数据集之一,比较适合我这样的新手入门。以下资料是为了实现对KITTI数据集的可视化操作。首先在官网下载对应的数据:http://www.cvlibs.net/datasets/kitti/eval_object.php?obj_benchmark=3d,下载后数据的目录文件结构如下所示:

├── dataset│ ├── KITTI│ │ ├── object│ │ │ ├──KITTI│ │ │ ├──ImageSets│ │ │ ├──training│ │ │ ├──calib & velodyne & label_2 & image_22. 环境准备

这里使用了一个kitti数据集可视化的开源代码:https://github.com/kuixu/kitti_object_vis,按照以下操作新建一个虚拟环境,并安装所需的工具包。其中千万不要安装python3.7以上的版本,因为vtk不支持。

# 新建python=3.7的虚拟环境conda create -n kitti_vis python=3.7 # vtk does not support python 3.8conda activate kitti_vis# 安装opencv, pillow, scipy, matplotlib工具包pip install opencv-python pillow scipy matplotlib# 安装3D可视化工具包(以下指令会自动安转所需的vtk与pyqt5)conda install mayavi -c conda-forge# 测试python kitti_object.py --show_lidar_with_depth --img_fov --const_box --vis3. KITTI数据集可视化

下面依次展示 KITTI 数据集可视化结果,这里通过设置 data_idx=10 来展示编号为000010的数据,代码中dataset需要修改为数据集实际路径。(最后会贴上完整代码)

def visualization(): import mayavi.mlab as mlab dataset = kitti_object(os.path.join(ROOT_DIR, '../dataset/KITTI/object')) # determine data_idx data_idx = 100 # Load data from dataset objects = dataset.get_label_objects(data_idx) print("There are %d objects.", len(objects)) img = dataset.get_image(data_idx) img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB) img_height, img_width, img_channel = img.shape pc_velo = dataset.get_lidar(data_idx)[:,0:3] calib = dataset.get_calibration(data_idx)

代码来源于参考资料,在后面会贴上我自己修改的测试代码。以下包含9种可视化的操作:

1. 图像显示def show_image(self): Image.fromarray(self.img).show() cv2.waitKey(0)

结果展示:

2. 图片上绘制2D bbox def show_image_with_2d_boxes(self): show_image_with_boxes(self.img, self.objects, self.calib, show3d=False) cv2.waitKey(0)

结果展示:

3. 图片上绘制3D bbox def show_image_with_3d_boxes(self): show_image_with_boxes(self.img, self.objects, self.calib, show3d=True) cv2.waitKey(0)

结果展示:

4. 图片上绘制Lidar投影 def show_image_with_lidar(self): show_lidar_on_image(self.pc_velo, self.img, self.calib, self.img_width, self.img_height) mlab.show()

结果展示:

5. Lidar绘制3D bbox def show_lidar_with_3d_boxes(self): show_lidar_with_boxes(self.pc_velo, self.objects, self.calib, True, self.img_width, self.img_height) mlab.show()

结果展示:

6. Lidar绘制FOV图 def show_lidar_with_fov(self): imgfov_pc_velo, pts_2d, fov_inds = get_lidar_in_image_fov(self.pc_velo, self.calib, 0, 0, self.img_width, self.img_height, True) draw_lidar(imgfov_pc_velo) mlab.show()

结果展示:

KITTI数据集可视化(一):点云多种视图的可视化实现

7. Lidar绘制3D图 def show_lidar_with_3dview(self): draw_lidar(self.pc_velo) mlab.show()

结果展示:

8. Lidar绘制BEV图

BEV图的显示与其他视图不一样,这里的代码需要有点改动,因为这里需要lidar点云的其他维度信息,所以输入不仅仅是xyz三个维度。改动代码:

# 初始pc_velo = dataset.get_lidar(data_idx)[:, 0:3]# 改为(要增加其他维度才可以查看BEV视图)pc_velo = dataset.get_lidar(data_idx)[:, 0:4]

测试代码:

def show_lidar_with_bev(self): from kitti_util import draw_top_image, lidar_to_top top_view = lidar_to_top(self.pc_velo) top_image = draw_top_image(top_view) cv2.imshow("top_image", top_image) cv2.waitKey(0)

结果展示:

9. Lidar绘制BEV图+2D bbox

同样,这里的代码改动与3.8节一样,需要点云的其他维度信息

def show_lidar_with_bev_2d_bbox(self): show_lidar_topview_with_boxes(self.pc_velo, self.objects, self.calib) mlab.show()

结果展示:

完整测试代码

参考代码:

import mayavi.mlab as mlabfrom kitti_object import kitti_object, show_image_with_boxes, show_lidar_on_image, \ show_lidar_with_boxes, show_lidar_topview_with_boxes, get_lidar_in_image_fov, \ show_lidar_with_depthfrom viz_util import draw_lidarimport cv2from PIL import Imageimport timeclass visualization: # data_idx: determine data_idx def __init__(self, root_dir=r'E:\Study\Machine Learning\Dataset3d\kitti', data_idx=100): dataset = kitti_object(root_dir=root_dir) # Load data from dataset objects = dataset.get_label_objects(data_idx) print("There are {} objects.".format(len(objects))) img = dataset.get_image(data_idx) img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB) img_height, img_width, img_channel = img.shape pc_velo = dataset.get_lidar(data_idx)[:, 0:3] # 显示bev视图需要改动为[:, 0:4] calib = dataset.get_calibration(data_idx) # init the params self.objects = objects self.img = img self.img_height = img_height self.img_width = img_width self.img_channel = img_channel self.pc_velo = pc_velo self.calib = calib # 1. 图像显示 def show_image(self): Image.fromarray(self.img).show() cv2.waitKey(0) # 2. 图片上绘制2D bbox def show_image_with_2d_boxes(self): show_image_with_boxes(self.img, self.objects, self.calib, show3d=False) cv2.waitKey(0) # 3. 图片上绘制3D bbox def show_image_with_3d_boxes(self): show_image_with_boxes(self.img, self.objects, self.calib, show3d=True) cv2.waitKey(0) # 4. 图片上绘制Lidar投影 def show_image_with_lidar(self): show_lidar_on_image(self.pc_velo, self.img, self.calib, self.img_width, self.img_height) mlab.show() # 5. Lidar绘制3D bbox def show_lidar_with_3d_boxes(self): show_lidar_with_boxes(self.pc_velo, self.objects, self.calib, True, self.img_width, self.img_height) mlab.show() # 6. Lidar绘制FOV图 def show_lidar_with_fov(self): imgfov_pc_velo, pts_2d, fov_inds = get_lidar_in_image_fov(self.pc_velo, self.calib, 0, 0, self.img_width, self.img_height, True) draw_lidar(imgfov_pc_velo) mlab.show() # 7. Lidar绘制3D图 def show_lidar_with_3dview(self): draw_lidar(self.pc_velo) mlab.show() # 8. Lidar绘制BEV图 def show_lidar_with_bev(self): from kitti_util import draw_top_image, lidar_to_top top_view = lidar_to_top(self.pc_velo) top_image = draw_top_image(top_view) cv2.imshow("top_image", top_image) cv2.waitKey(0) # 9. Lidar绘制BEV图+2D bbox def show_lidar_with_bev_2d_bbox(self): show_lidar_topview_with_boxes(self.pc_velo, self.objects, self.calib) mlab.show()if __name__ == '__main__': kitti_vis = visualization() # kitti_vis.show_image() # kitti_vis.show_image_with_2d_boxes() # kitti_vis.show_image_with_3d_boxes() # kitti_vis.show_image_with_lidar() # kitti_vis.show_lidar_with_3d_boxes() # kitti_vis.show_lidar_with_fov() # kitti_vis.show_lidar_with_3dview() # kitti_vis.show_lidar_with_bev() kitti_vis.show_lidar_with_bev_2d_bbox() # print('...') # cv2.waitKey(0)

此外,下面再提供两份可视化代码。

4. 点云可视化

这里的同样使用的是上述的图例,且直接输入的KITTI数据集的.bin文件,即可显示点云图像。

参考代码:import numpy as npimport mayavi.mlabimport os# 000010.bin这里需要填写文件的位置# bin_file = '../data/object/training/velodyne/000000.bin'# assert os.path.exists(bin_file), "{} is not exists".format(bin_file)kitti_file = r'E:\Study\Machine Learning\Dataset3d\kitti\training\velodyne\000100.bin'pointcloud = np.fromfile(file=kitti_file, dtype=np.float32, count=-1).reshape([-1, 4])# pointcloud = np.fromfile(str("000010.bin"), dtype=np.float32, count=-1).reshape([-1, 4])print(pointcloud.shape)x = pointcloud[:, 0] # x position of pointy = pointcloud[:, 1] # y position of pointz = pointcloud[:, 2] # z position of pointr = pointcloud[:, 3] # reflectance value of pointd = np.sqrt(x ** 2 + y ** 2) # Map Distance from sensorvals = 'height'if vals == "height": col = zelse: col = dfig = mayavi.mlab.figure(bgcolor=(0, 0, 0), size=(640, 500))mayavi.mlab.points3d(x, y, z, col, # Values used for Color mode="point", colormap='spectral', # 'bone', 'copper', 'gnuplot' # color=(0, 1, 0), # Used a fixed (r,g,b) instead figure=fig, )x = np.linspace(5, 5, 50)y = np.linspace(0, 0, 50)z = np.linspace(0, 5, 50)mayavi.mlab.plot3d(x, y, z)mayavi.mlab.show()输出结果:

ps:这里的输出点云结果相比上面的点云输出结果更加的完善,而且参考的中心坐标点也不一样。

5. 鸟瞰图可视化

代码中的鸟瞰图范围可以自行设置。同样,输入的也只需要是.bin文件即可展示其鸟瞰图。

参考代码:import numpy as npfrom PIL import Imageimport matplotlib.pyplot as plt# 点云读取:000010.bin这里需要填写文件的位置kitti_file = r'E:\Study\Machine Learning\Dataset3d\kitti\training\velodyne\000100.bin'pointcloud = np.fromfile(file=kitti_file, dtype=np.float32, count=-1).reshape([-1, 4])# 设置鸟瞰图范围side_range = (-40, 40) # 左右距离# fwd_range = (0, 70.4) # 后前距离fwd_range = (-70.4, 70.4)x_points = pointcloud[:, 0]y_points = pointcloud[:, 1]z_points = pointcloud[:, 2]# 获得区域内的点f_filt = np.logical_and(x_points > fwd_range[0], x_points < fwd_range[1])s_filt = np.logical_and(y_points > side_range[0], y_points < side_range[1])filter = np.logical_and(f_filt, s_filt)indices = np.argwhere(filter).flatten()x_points = x_points[indices]y_points = y_points[indices]z_points = z_points[indices]res = 0.1 # 分辨率0.05mx_img = (-y_points / res).astype(np.int32)y_img = (-x_points / res).astype(np.int32)# 调整坐标原点x_img -= int(np.floor(side_range[0]) / res)y_img += int(np.floor(fwd_range[1]) / res)print(x_img.min(), x_img.max(), y_img.min(), x_img.max())# 填充像素值height_range = (-2, 0.5)pixel_value = np.clip(a=z_points, a_max=height_range[1], a_min=height_range[0])def scale_to_255(a, min, max, dtype=np.uint8): return ((a - min) / float(max - min) * 255).astype(dtype)pixel_value = scale_to_255(pixel_value, height_range[0], height_range[1])# 创建图像数组x_max = 1 + int((side_range[1] - side_range[0]) / res)y_max = 1 + int((fwd_range[1] - fwd_range[0]) / res)im = np.zeros([y_max, x_max], dtype=np.uint8)im[y_img, x_img] = pixel_value# imshow (灰度)im2 = Image.fromarray(im)im2.show()# imshow (彩色)# plt.imshow(im, cmap="nipy_spectral", vmin=0, vmax=255)# plt.show()结果展示:

后续的工作会加深对点云数据的理解,整个可视化项目的工程见:KITTI数据集的可视化项目,有需要的朋友可以自行下载。

参考资料:

1. KITTI自动驾驶数据集可视化教程

2. kitti数据集在3D目标检测中的入门

3. kitti数据集在3D目标检测中的入门(二)可视化详解

4. kitti_object_vis项目

本文链接地址:https://www.jiuchutong.com/zhishi/300214.html 转载请保留说明!

上一篇:js算法 字母大小写转换(如何在js中给字母排序)

下一篇:【vue2】使用elementUI进行表单验证实操(附源码)(vue el-)

  • 苹果11如何分屏多窗口(苹果11如何分屏微信)

    苹果11如何分屏多窗口(苹果11如何分屏微信)

  • 抖音短视频怎么去水印保存(抖音短视频怎么剪辑制作)

    抖音短视频怎么去水印保存(抖音短视频怎么剪辑制作)

  • 投屏网速慢怎么办(投屏网速不稳定)

    投屏网速慢怎么办(投屏网速不稳定)

  • 华为如何自由复制文字(华为复制功能怎么用)

    华为如何自由复制文字(华为复制功能怎么用)

  • iphone11指纹支付咋办(iphone11pro指纹支付咋办)

    iphone11指纹支付咋办(iphone11pro指纹支付咋办)

  • 登录中国移动收不到验证码怎么办(登录中国移动收不到验证码怎么办多久能收到)

    登录中国移动收不到验证码怎么办(登录中国移动收不到验证码怎么办多久能收到)

  • 苹果6s plus什么时候在中国上市的(苹果6splus什么屏幕材质)

    苹果6s plus什么时候在中国上市的(苹果6splus什么屏幕材质)

  • 序列号ff开头是哪里的(序列号ff开头是官换机吗)

    序列号ff开头是哪里的(序列号ff开头是官换机吗)

  • qq下拉小程序没有了(qq下拉小程序没了怎么办)

    qq下拉小程序没有了(qq下拉小程序没了怎么办)

  • 淘宝极速退款条件申请(淘宝极速退款条件是用菜鸟裹裹才能退吗?)

    淘宝极速退款条件申请(淘宝极速退款条件是用菜鸟裹裹才能退吗?)

  • 为什么手机会自动反向充电(为什么手机会自动重启)

    为什么手机会自动反向充电(为什么手机会自动重启)

  • 苹果手机home键按下去咔咔响声(苹果手机home键失灵了怎么开锁)

    苹果手机home键按下去咔咔响声(苹果手机home键失灵了怎么开锁)

  • word上标和下标怎么打(word上标和下标怎么设置)

    word上标和下标怎么打(word上标和下标怎么设置)

  • 手机相册怎么不显示照片(手机相册怎么不显示视频)

    手机相册怎么不显示照片(手机相册怎么不显示视频)

  • 抖音如何关掉商品橱窗功能(抖音如何关掉商城功能)

    抖音如何关掉商品橱窗功能(抖音如何关掉商城功能)

  • 刷宝提现不了(刷宝提现越来越麻烦)

    刷宝提现不了(刷宝提现越来越麻烦)

  • 退群以后发的消息还在吗(退群以后发的消息能看吗)

    退群以后发的消息还在吗(退群以后发的消息能看吗)

  • 小红书怎么搜不到了(小红书怎么搜不到通讯录好友了)

    小红书怎么搜不到了(小红书怎么搜不到通讯录好友了)

  • 怎样在手机上修改银行预留手机号(怎样在手机上修改医保卡密码)

    怎样在手机上修改银行预留手机号(怎样在手机上修改医保卡密码)

  • 海底捞微信怎么排号(海底捞微信怎么预约美甲)

    海底捞微信怎么排号(海底捞微信怎么预约美甲)

  • 如何使用淘宝优惠券(淘宝优选入口在哪)

    如何使用淘宝优惠券(淘宝优选入口在哪)

  • 路由器怎么重置?(路由器怎么重置网络)

    路由器怎么重置?(路由器怎么重置网络)

  • Vue3通透教程【十三】TS简单类型详解

    Vue3通透教程【十三】TS简单类型详解

  • 自己房子出租发布信息
  • 增值税是价内税,消费者是税款的最终负担者( ) A对 B错
  • 利润表的期末余额
  • 个体户单位性质是什么
  • 如何确定固定资产的原始价值
  • 职工福利费允许结转以后年度扣除吗
  • 小微企业印花税减免最新政策
  • ATM可以存款到对公账户吗?
  • 外经证预交几个点
  • 非居民企业直接投资居民企业取得股息
  • 不得抵扣的进项税额转出会计分录
  • 资产的计量属性主要包括哪些
  • 一般纳税人收到专票怎么做分录
  • 三栏式明细分类账的账页格式适用于
  • 建筑工程居间服务协议
  • 长期股权投资成本法初始计量
  • 领备用金时会计怎么做账
  • 行政单位以前年度会计错账调整
  • 离婚后房产过户需要多少钱
  • 实际成本法的会计分录怎么写
  • 公司室内装修
  • 员工劳务报酬分录
  • 系统安装的步骤
  • 购买一台电脑会计分录
  • 关于工程材料的图书有哪些
  • Apache+php+mysql在windows下的安装与配置图解(最新版)
  • PHP:stream_set_blocking()的用法_Stream函数
  • 民办非企业的税金及附加放在哪个科目
  • 企业所得税的纳税人包括哪些
  • 在建工程的核算内容和范围是什么
  • 分次预缴税额怎么算
  • phpcms如何配置404
  • 免征和抵扣税额的区别
  • 行政事业性收据丢了怎么办
  • 将织梦dedecms转换到wordpress
  • 企业并购的主要特征是a扩大企业规模
  • mongodb数据库操作题
  • 在建工程转固定资产摘要怎么写
  • mysql常见错误
  • android 矩阵运算
  • sql删除表中的某一行
  • mysql主从配置详解
  • 改良生产线发生的职工薪酬计入期间费用
  • 什么是异地企业
  • 公司购买汽车如何抵成本
  • 六税一费减免
  • 增值税加计扣除最新政策2023
  • 暂估入库的商品怎么做账
  • 佣金支出和手续费一样吗
  • 内含增长率的推导
  • 收到法人的借款怎样写摘要
  • 来料加工和进料加工的相同之处
  • 农产品收购单怎么做会计分录
  • 售后回租含手续费吗
  • 增值税普票遗失可以补开吗
  • 购买的垃圾袋计入什么科目
  • 如何获取数据库的sid
  • sqlserver数据库怎么导出
  • Windows 2003作中转VPN服务器多路由共享上网的方法
  • solaris版本查询
  • Win10 Mobile RS2预览版14951升级遭遇卡在0%的解决办法
  • deepin-win
  • xp系统怎么改文件类型
  • windows 10 mobile--移动版
  • 如何卸载win8系统
  • linux无法使用ls
  • shader教程
  • chrome heat
  • js querySelector() 使用方法
  • “android”
  • 健壮的什么
  • python datetime.now
  • python怎么制作
  • javascript parseInt 函数分析(转)
  • 手机截取电脑屏幕
  • 税务举报电话号码
  • 浙江银行上班时间查询
  • 庐山坐缆车上去后还要买门票吗?
  • 什么是集体约谈制
  • 车辆完税证明有啥用
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设