【发布时间】:2021-06-01 20:54:57
【问题描述】:
我以这种方式创建了这个列表:
#Alotting Black Holes at z=6
bhs=[0]*1000
for i in tqdm(range(0,1000),position=0, leave=True):
if len(mass_array1[i])!=0:
bhs[i]=np.zeros(len(mass_array1[i]))
else:
bhs[i]=np.zeros(1)
for j in range (len(mass_array1[i])):
bhs[i][j]=np.random.lognormal(np.log(MbhthShimasaku(mass_array1[i],6)[j]),np.log(5))
我需要将结果保存在文本文件中。我试过 numpy.savetxt、pickle.dump 和 open():
open()
with open("bhs.txt", 'w') as file:
for row in bhs:
s = " ".join(map(str, row))
file.write(s+'\n')
#Result .txt file:
0.0
0.0
0.0
0.0
1937651.7861915156 246221.20328840986 226756.87389065413
0.0
0.0
numpy.savetxt()
bhs=np.array(bhs)
np.savetxt('bhs.txt',bhs,fmt='%s')
#Result .txt file:
[0.]
[0.]
[0.]
[0.]
[26447480.89508711 1097038.92200952 971383.67441455]
[0.]
[0.]
[0.]
[0.]
[0.]
泡菜
bhs.append(bhs)
tupleA=tuple(bhs)
filename = 'bhs.p'
with open(filename, 'wb') as filehandler:
pickle.dump(tupleA, filehandler)
#Result .p file
array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([1937651.78619152, 246221.20328841, 226756.87389065])
我无法从所有这些保存的文件中取回原始数组/列表。当我尝试使用这些加载的列表中的任何一个时,都会出现某种错误:
np.loadtxt
could not convert string to float: '[0.]'
open()
my_file = open("bhs.txt", "r")
content = my_file.read()
content_list = content.split(",")
my_file.close()
print(content_list)
[0.]\n[0.]\n[26447480.89508711 1097038.92200952 971383.67441455]\n[0.]\n[0.]\n[0.]\n[0.]\n[0.]\n[0.]\n[0.]\n[0.]\n
列表中的 bhs 示例
array([1461403.98258597]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([0.]), array([26447480.89508711, 1097038.92200952, 971383.67441455]),
我怎样才能说出我的多维列表,以便我可以准确地回到我开始时的内容?
额外:mass_array1 文件
https://drive.google.com/file/d/1Kdmv1fcbDelEzGmi4BOE4HjUbM7Cg23b/view?usp=sharing
这就是我将它导入 python 的方式:
您需要先将文件解压到一个文件夹中。
dirlist=["bh2e10"]
import time
mass_array1=[0]*1000
#print(mass_array)
#read all the files
for i,X in enumerate(dirlist):
exec('filelist=glob.glob("%s/test*.dat")'%(X))
#exec("mass_array%s=[]"%X)
initial_mass=[]
for j,Y in tqdm(enumerate(filelist),position=0, leave=True, total=1000):
Y=Y.replace(os.sep, '/')
#Z=int(Y[10:13])
Z=int(re.findall("\d+", Y)[2])
#print(Z)
mass_array1[Z]=[]
#print('i=',Z,end="\r")
#print('i=',Z,end="\r")
exec("initial_partial=np.loadtxt('%s',max_rows=1)"%(Y))
exec("initial_mass=np.append(initial_mass,initial_partial)")
exec("mass_partial=np.loadtxt('%s',skiprows=1)"%(Y))
mass_array1[Z]=np.append(mass_partial,mass_array1[Z])
#mass_array1[Z]=mass_partial
【问题讨论】:
-
你能用h5py吗?
-
JSON 是一种非常实用的方式。
-
您能否更新您的帖子并将 bhs 示例粘贴为列表。
-
真的,每一行都应该有相同数量的元素。
-
为什么你为什么要使用
exec?这里非常危险,完全没有必要。