搞研究需要解析质谱的mzml文件,看了好多内容也不知道应该怎么弄,今天突然就行了。太激动了纪念下。
from pyteomics import mzml
import numpy as np
import pandas as pd
# 指定mzML文件路径
file_path = ‘你的路径’
dict={}
# 使用mzml.read()函数读取mzML文件
with mzml.read(file_path) as spectra:
for spectrum in spectra:
# 获取m/z数组和强度数组
mz_array = spectrum['m/z array']
intensity_array = spectrum['intensity array']
dict['y']=mz_array
dict['z']=intensity_array
df = pd.DataFrame(dict)
# 指定CSV文件的保存路径
csv_file_path = ‘你的路径’
# 将DataFrame写入CSV文件
df.to_csv(csv_file_path, index=False, encoding='utf-8')
是不是特别简单!!!本来pyteomics安装不了,结果整体升级了下就装上了。文件多写个循环就OK了,哈哈哈哈哈。