我有一个Quantum Chemistry程序的输出,我希望从中提取表格数据以输入到我25年前写的FORTRAN程序的Python端口中。
一些输出文件相当长,多达6000行,因此无法使用电子表格进行处理。
典型的表格具有以下形式:
CARTESIAN COORDINATES
1 C 0.011987266 -0.003842185 0.006578784
2 H 1.097152909 -0.003956163 0.013339310
3 H -0.349612312 1.019316731 0.001903075
4 H -0.344276148 -0.517463019 -0.880495291
5 H -0.355315644 -0.513266496 0.891567896
我并不是要有人为我编写Python代码,而是要在可用代码的迷宫中给我一些指导。
我建议您调查np.genfromtxt。以下代码段将从存储在名为的文件中的问题中读取示例数据data.txt
。
import numpy as np
data = np.genfromtxt('data.txt', skip_header=2, dtype=[('id', 'i8'),('label','S1'),('x','f8'),('y','f8'),('z','f8')])
print(data)
输出量
[(1, b'C', 0.01198727, -0.00384219, 0.00657878)
(2, b'H', 1.09715291, -0.00395616, 0.01333931)
(3, b'H', -0.34961231, 1.01931673, 0.00190307)
(4, b'H', -0.34427615, -0.51746302, -0.88049529)
(5, b'H', -0.35531564, -0.5132665 , 0.8915679 )]
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句