mdpout_descriptors文件项目过多,虽然可以直接用R进行作图,无奈我不会使用R,所以网上半摸索学习pandas写了一个简单的方法处理文本,如下:

1
2
3
4
5
6
7
8
9
10
11
12
13
#coding=utf-8
# author:kangsgo
# fepocket2 volume 拆分工具
# ref:http://blog.sina.com.cn/s/blog_a73687bc0101eenc.html
import pandas as pd
id = pd.read_table("mdpout_descriptors.txt",sep=" +",engine='python')
name2=id['snapshot']
time=id['pock_volume']
#method1
pd.concat([name2,time],axis=1).to_csv('volume.txt',sep=" ",index=False)
#method2 先全读进来,然后 df.iloc[:,[:3]]
#https://www.v2ex.com/t/306789#reply2