mzml数据处理及质谱图像绘制

数据处理

from pyteomics import mzml
import numpy as np
import pandas as pd
import matplotlib
import matplotlib.pyplot as plt
import json
import csv
import os
import math
def change_mzML(fd):for spectrum in mzml.read(fd):if(spectrum.get('base peak m/z') >= 114 and spectrum.get('base peak m/z')< 115):dfmz = pd.DataFrame(spectrum.get('m/z array'))dfin = pd.DataFrame(spectrum.get('intensity array'))dfmz.to_csv(fd+"114m_z.csv", mode='a+', index =False, header = True)dfin.to_csv(fd+"114intensity.csv", mode='a+', index =False, header = True)for spectrum in mzml.read(fd):if(spectrum.get('base peak m/z') >= 115 and spectrum.get('base peak m/z')< 116):dfmz = pd.DataFrame(spectrum.get('m/z array'))dfin = pd.DataFrame(spectrum.get('intensity array'))dfmz.to_csv(fd+"115m_z.csv", mode='a+', index =False, header = True)dfin.to_csv(fd+"115intensity.csv", mode='a+', index =False, header = True)for spectrum in mzml.read(fd):if(spectrum.get('base peak m/z') >= 116 and spectrum.get('base peak m/z')< 117):dfmz = pd.DataFrame(spectrum.get('m/z array'))dfin = pd.DataFrame(spectrum.get('intensity array'))dfmz.to_csv(fd+"116m_z.csv", mode='a+', index =False, header = True)dfin.to_csv(fd+"116intensity.csv", mode='a+', index =False, header = True)for spectrum in mzml.read(fd):if(spectrum.get('base peak m/z') >= 117 and spectrum.get('base peak m/z')< 118):dfmz = pd.DataFrame(spectrum.get('m/z array'))dfin = pd.DataFrame(spectrum.get('intensity array'))dfmz.to_csv(fd+"116m_z.csv", mode='a+', index =False, header = True)dfin.to_csv(fd+"116intensity.csv", mode='a+', index =False, header = True)
for fn in os.listdir('C:/hqh/data/breast/'):if fn.endswith('.mzML'):fd = os.path.join('C:/hqh/data/breast/', fn)change_mzML(fd)

质谱图绘制

import numpy as np
import pandas as pd
import matplotlib
import matplotlib.pyplot as plt
import json
import csv
import os
import math
# 获取csv其中一行
def getColumn(filename, column):results = csv.reader(open(filename), delimiter=",")return [result[column] for result in results]
# 转换成浮点数
def getFloat(result):res_float_m = map(float, result)return list(res_float_m)
# 获取数组中的最大值
def getMax(res_float):max_n = 0for res in res_float:if res > max_n:max_n = resreturn max_n
x = getColumn("C:/hqh/1.csv",0)
y = getColumn("C:/hqh/1.csv",1)
x_float = getFloat(x)
y_float = getFloat(y)
max_y = getMax(y_float)
y_float = np.array(y_float)*100/max_y # 将数组中全部值化为相对强度
plt.xlim(0, 2000)
plt.bar(x_float, y_float, width = 5)
plt.savefig("C:/hqh/examples.png")

mzml数据处理及质谱图像绘制相关推荐

【Python】函数图像绘制：二维图像、三维图像、散点图、心形图
[Python]函数图像绘制:二维图像.三维图像.散点图.心形图所有需要用的包二维图像三维图像散点图绘制心形图绘制所有需要用的包 from mpl_toolkits.mplot3d imp ...
【Android 内存优化】自定义组件长图组件 ( 获取图像宽高 | 计算解码区域 | 设置图像解码属性复用像素格式 | 图像绘制 )
文章目录一.获取图像真实宽高二.计算解码区域三.设置解码参数内存复用像素格式四.图像绘制五.执行效果六.源码及资源下载官方文档 API : BitmapRegionDecoder 在 ...
Python：利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等)
Python:利用python编程实现三维图像绘制展示(六面体旋转.三维球柱状体.下雪场景等) 目录利用python编程实现三维图像绘制展示(六面体旋转.三维球柱状体.下雪场景等) 1.实现六面体旋 ...
计算机图形图像电影论文,计算机图形图像绘制技术地现状分析及应用发展.doc...
计算机图形图像绘制技术地现状分析及应用发展计算机图形图像绘制技术的现状分析及应用发展-计算机论文计算机图形图像绘制技术的现状分析及应用发展金晓倩 JIN Xiao-qian (四川建筑职业技术学 ...
计算机图形设计论文真实图形生成技术的发展,绘制技术论文,关于计算机图形图像绘制技术的现状应用相关参考文献资料-免费论文范文...
导读:这是一篇与绘制技术论文范文相关的免费优秀学术论文范文资料,为你的论文写作提供参考. (四川建筑职业技术学院,德阳618000) (Sichuan College of Architectural ...
【OpenCV】OpenCV实战从入门到精通之 -- 基本图像绘制
目录 1.DrawEllipse()函数的写法 2.DrawFilledCircle()函数的写法 3.DrawPolygon()函数的写法 4.DrawLine()函数的写法基本绘制图形代码: m ...
Python符号计算入门及隐函数图像绘制
感谢国防科大刘万伟老师提供了隐函数图像绘制的原始问题以及完美答案,我又补充了一点符号计算的基础知识. >>> from sympy import * #定义符号 >>&g ...
【matplotlib笔记】3D图像绘制
3D图像绘制 3D线图 3D散点图 3D曲面图 z = x^2 + y^2 3D线图 # 3D线图绘制 import matplotlib.pyplot as plt import numpy as ...
Python Imaging Library:ImageDraw Module（图像绘制模块）
Python Imaging Library:ImageDraw Module(图像绘制模块) # 图像绘制库 Mode_NewImg = "RGB" Width_NewImg = ...

mzml数据处理及质谱图像绘制

mzml数据处理及质谱图像绘制相关推荐

最新文章

热门文章