PCD文件去除曲率的脚本
在写一个重建算法的时候需要用到点坐标和法向的数据文件,于是向利用pcl中的法向计算模块来生成法向。输出后法向文件中包含曲率信息,但是这是不需要的。于是自己写了一个python小脚本实现格式转换。
#--coding:utf-8-- import time import numpy as np from sys import argv script, input_file = argv input_data = open(input_file,"r") output_data = open("outdata.asc","w") doc = ''' 本脚本用于数据转换,将PCL生成的文件进行转换 输入为 点云数据+法向信息+曲率 PCD文件 输出为 点云数据+法向信息 ASC文件 转换后曲率信息被滤除,并将格式转换为通用的asc文件格式 ''' print(doc) #转换函数 def transform(input_data, output_data): pointList = [] for line in input_data: data = line.split() if is_number(data[0]): output_data.write("%f %f %f %f %f %f\n" %((float)(data[0]), (float)(data[1]), (float)(data[2]), (float)(data[3]), (float)(data[4]), (float)(data[5]))) else: pass #判断字符串是否为数字 def is_number(s): try: float(s) return True except ValueError: pass return False transform(input_data, output_data)
原文件内容为:
# .PCD v0.7 - Point Cloud Data file format VERSION 0.7 FIELDS x y z normal_x normal_y normal_z curvature SIZE 4 4 4 4 4 4 4 TYPE F F F F F F F COUNT 1 1 1 1 1 1 1 WIDTH 26381 HEIGHT 1 VIEWPOINT 0 0 0 1 0 0 0 POINTS 26381 DATA ascii 52.713619 44.306831 2.500001 -0.67111164 -0.08553692 0.73640519 0.00061548286 52.54211 45.728249 2.500001 -0.67243981 -0.090179443 0.73463756 0.00069258711 52.30574 47.160099 2.500001 -0.67485178 -0.10609017 0.73028761 0.00043992131 52.022919 48.57592 2.500001 -0.67386091 -0.11686694 0.72955716 0.00037569425 51.720371 49.97781 2.500001 -0.66632128 -0.11713342 0.73640722 0.00062433031 51.405071 51.394611 2.500001 -0.65537727 -0.11913724 0.74584651 0.00090703077 51.068241 52.847191 2.500001 -0.6404919 -0.1198745 0.75855136 0.001082067 50.717159 54.314861 2.500001 -0.62284762 -0.120147 0.77306247 0.0012497745 50.349918 55.789589 2.500001 -0.60351694 -0.12188542 0.7879793 0.0013898062 49.98164 57.264919 2.500001 -0.56129658 -0.1375796 0.81609923 0.0014789379 49.60405 58.715778 2.500001 -0.53091556 -0.14410463 0.83508235 0.0021452638 49.167542 60.149521 2.500001 -0.45673963 -0.16522408 0.87412238 0.0064957533 48.665119 61.58548 2.500001 -0.27892584 -0.12767911 0.95178694 0.015811486
将该文件命名为transfor.py并在终端执行指令:
$python transfor.py input.pcd
Python中的split函数的默认参数是空格, 即对每行数据安装空格分割, 不管空格的数量是一个还是多个, 全部忽略. 但是一旦添加自定义的参数, 则严格按照定义参数分割. 比如, 当定义的参数为一个空格时, 将按一个空格分割, 两个空格时将按两个空格分割. 该处需要注意.
转换后的文件内容为:
52.713619 44.306831 2.500001 -0.671112 -0.085537 0.736405 52.542110 45.728249 2.500001 -0.672440 -0.090179 0.734638 52.305740 47.160099 2.500001 -0.674852 -0.106090 0.730288 52.022919 48.575920 2.500001 -0.673861 -0.116867 0.729557 51.720371 49.977810 2.500001 -0.666321 -0.117133 0.736407 51.405071 51.394611 2.500001 -0.655377 -0.119137 0.745847 51.068241 52.847191 2.500001 -0.640492 -0.119874 0.758551 50.717159 54.314861 2.500001 -0.622848 -0.120147 0.773062 50.349918 55.789589 2.500001 -0.603517 -0.121885 0.787979 49.981640 57.264919 2.500001 -0.561297 -0.137580 0.816099 49.604050 58.715778 2.500001 -0.530916 -0.144105 0.835082 49.167542 60.149521 2.500001 -0.456740 -0.165224 0.874122 48.665119 61.585480 2.500001 -0.278926 -0.127679 0.951787 48.067780 62.904530 2.500001 -0.263447 -0.142550 0.954083 47.309990 63.906979 2.500001 -0.082274 -0.155248 0.984443 46.385601 64.398064 2.500001 -0.025878 -0.196875 0.980087 45.383808 64.317551 2.500001 0.015420 -0.199873 0.979701 44.408440 63.701160 2.500001 -0.232693 0.125587 -0.964408 43.552479 62.677849 2.500001 -0.300091 0.156567 -0.940974 42.826420 61.414768 2.500001 -0.369369 0.196475 -0.908275 42.207859 60.036179 2.500001 -0.555614 0.201118 -0.806749 41.684139 58.620060 2.500001 -0.607999 0.191823 -0.770416