目录
.1 文本方式读取
1.2 stl解析
1.3 stl创建
.2 把点转换为.stl
.1 文本方式读取
代码如下
stl_path='/home/pxing/codes/point_improve/data/003_cracker_box/0.stl'points=[]
f = open(stl_path)
lines = f.readlines()
prefix='vertex'
num=3
for line in lines:#print (line)if line.startswith(prefix):values = line.strip().split()#print(values[1:4])if num%3==0:points.append(values[1:4])num=0num+=1#print(type(line))
print(points)
f.close()
该代码可以直接将stl文件中对点直接都输出
改进版本
import numpy as npstl_path='/home/pxing/codes/point_improve/data/003_cracker_box/0.stl'
points=[]
f = open(stl_path)
lines = f.readlines()
prefix='vertex'
num=3
for line in lines:#print (line)if line.startswith(prefix):values = line.strip().split()#print(values[1:4])if num%3==0:points.append(values[1:4])num=0num+=1#print(type(line))
points=np.array(points)
f.close()
print(points.shape)
np.save("/home/pxing/codes/point_improve/feature_get/point_get/test.npy", points)
把输出的点转换为npy并且保存
1.2 stl解析
.stl文件格式如下图所示:
那么读取该文件就需要知道stl的构成
stl是一种表示三角网格的文件格式
STL只能用来表示封闭的面或者体,stl文件有两种:一种是ASCII明码格式,另一种是二进制格式
ASCII明码格式:
STL文件的首行给出了文件路径及文件名
STL三维模型就是由一系列这样的三角面片构成
三角面片的信息单元 facet 是一个带矢量方向的三角面片
每一个facet由7 行数据组成:
facet normal 是三角面片指向实体外部的法矢量坐标
outer loop 代表随后的3行数据分别是三角面片的3个顶点坐标
3顶点沿指向实体外部的法矢量方向逆时针排列
solid [filename] //文件名,可以是任何字符
facet normal [i j k] //面的法线,i、j、k为三个分量,各分量之间用空格隔开,不能用逗号隔开
outor loop
vertex [x y z] //三角面片的第一个点,x、y、z三个坐标之间要用空格隔开
vertex [x y z] //三角面片的第二个点
vertex [x y z] //三角面片的第三个点
endloop
endfacet //完成一个三角面片的定义
、、、、 //其他facet
endsolid [filename] //完成一个stl文件的定义
示例:
在电脑上新建一个txt文件,然后编辑。编辑完了之后保存,并将文件后缀名改为stl,即可完成一个创建了一个stl文件。如下图,利用文本编辑器将一个四面体的几何数据保存到文件中
效果如下图所示:
1.3 stl创建
创建一个三角形面片
solid mystl
facet normal 0 -1 0
outor loop
vertex -100 0 0
vertex 100 0 0
vertex 0 0 -100
endloop
endfacet
end solid mystl
.2 把点转换为.stl
import numpy as np
import sys
sys.path.append('/home/pxing/codes/grasp_ROI_get/')
from gripper_area import area_getarr=np.array(area_get.main())
#print(arr[0,0:3])
with open('stl_process/data/test1.stl','a') as file0:print('solid mystl ',file=file0)print('facet normal' ,0,0,0,file=file0)print('outor loop',file=file0)print('vertex',arr[0,0],arr[0,1],arr[0,2 ],file=file0)print('vertex',arr[1,0],arr[1,1],arr[1,2 ],file=file0)print('vertex',arr[2,0],arr[2,1],arr[2,2 ],file=file0)print('endloop\nendfacet',file=file0)print('end solid mystl',file=file0)