ヘッダー文字列と数値配列の読み書き

Question

私はPythonには本当に新しく、とても簡単なことはできません。文字列と数値を扱う場合の例を見ていますが、数行の文字列ヘッダとそれに続く数値の配列を扱う方法ではありません。私は部分的な解決策を見つけましたが、データの最初の行が部分的に切り刻まれています。これは私がやったことです：

import numpy as np 
infil="infile.txt" 
outfil="outfile.txt" 
fi = open(infil,"r") 
fo = open(outfil,"w") 
lines = fi.readlines() 
fo.writelines(lines[0:3]) 
fi.close() 
x1, x2, x3, x4, x5, x6, x7 = np.loadtxt(infil, skiprows=3, unpack=True) 
# Some computations with those columns of numbers occurs here. 
data = np.array([x1,x2,x3,x4,x5,x6,x7] 
data = data.T 
np.savetxt(outfil, data, fmt='%f') 
fo.close()  

上記もほぼ同じです。最初の行のデータの最初の半分が欠落しているだけです。

ご協力いただきますようお願い申し上げます。

おかげで、 ビル

Answer 1

私は多分、よりエレガントな解決策を追加したい： np.savetxtは前のデータをファイルに書き込まれたヘッダ文字列を受け付けます。我々は、データがどのように見えるか見ることができるように

np.savetxt(filename, data, header="my header\nand a second line")

Answer 2

はい！それを解決！その解決策は面白いですが。問題は、3行のヘッダーの後に適切にフォーマットされたデータの多くの行を含む出力ファイルを持つことでした。最初の1.5行のデータは上記のアルゴリズムで切断されました。私が以下に示す解決策は機能しますが、一時ファイルを使用して書式設定されたデータ行を最終出力ファイルに書き込む文字列に渡すため、エレガントではありません。このアルゴリズムを改善し、改良するための提案は感謝しています。それまでは、これはうまく動作します。

import numpy as np    # Load useful routines for arrays. 

infil="infile.txt"    # Input file name. 
tfil ="tmp.txt"     # Temporary file name. 
outfil="outfile.txt"   # Output file name. 

fi = open(infil,"r")   # Open input file. 
ft = open(tfil,"w")    # Open temporary file for writing into. 
fo = open(outfil,"w")   # Open output file. 

lines = fi.readlines()   # Read in the input file. 

fi.close()      # Close input file. 




x1, x2, x3, x4, x5, x6, x7 = np.loadtxt(infil, skiprows=3, unpack=True) 
          # Read in and unpack the numerical arrays from the input file. 

# Computation with those numerical arrays will occur here when I'm ready.        


data = np.array([x1,x2,x3,x4,x5,x6,x7]) 
          # Pack the data arrays together. 

data = data.T     # Transpose the data into columns. 


np.savetxt(tfil, data, fmt='%f\t')  # Write to temporary file. 


ft.close()      # Close temporary file. 

ft = open(tfil,"r")    # Now open temporary file to read in the formatted data lines. 

dlines = ft.readlines()   # Read in the temporary file. 

fo.writelines(lines[0:3])  # Write the header to the output file. 

fo.writelines(dlines)   # Write the formatted data lines to the output file. 

fo.close()      # Close output file.