将mzCDF结构转换为峰值列表
[
Peaklist
,时代
) = mzcdf2peaks (mzCDFStruct
)
mzCDFStruct |
MATLAB®结构,其中包含来自netCDF文件的信息,例如由 |
Peaklist |
以下任何一项:
|
时代 |
与液相色谱/质谱(LC/MS)或气相色谱/质谱(GC/MS)数据集相关的保留时间向量的标量。如果 |
[
提取峰信息Peaklist
,时代
) = mzcdf2peaks (mzCDFStruct
)mzCDFStruct
,一个MATLAB结构,包含来自netCDF文件的信息,例如由mzcdfread
功能,并创建Peaklist
,单个矩阵或包含质量/电荷(m/z)值和离子强度值的矩阵单元阵列,以及时代
,与液相色谱/质谱(LC/MS)或气相色谱/质谱(GC/MS)数据集相关的保留时间的标量或矢量。
mzCDFStruct
包含与netCDF文件中的变量和全局属性对应的字段。如果netCDF变量包含本地属性,则创建一个附加字段,该字段的名称是附加的变量名称_属性
.根据质谱仪软件的不同,这些字段的数量和名称也会有所不同,但通常都是有的mass_values
和intensity_values
字段。
在下面的示例中,文件results.cdf
没有提供。
使用mzcdfread
函数将netCDF文件作为结构读入MATLAB软件中。然后从结构中提取峰值信息。
mzcdf_struct = mzcdfread(“results.cdf”);[peaks,time] = mzcdf2peaks(mzcdf_struct) peaks = [7008x2 single] [7008x2 single] [7008x2 single] time = 8.3430 12.6130 16.8830 21.1530
为每个峰值列表(保留时间)创建一个颜色映射。
颜色= hsv(元素个数(山峰));
创建一个峰值的三维图形,并添加标签。
图hold on for I = 1:numel(峰值)t = repmat(时间(I),大小(峰值{I},1),1);plot3(t,峰{i}(:,1),峰{i}(:,2),'color',colors(i,:)结束视图(70,60)xlabel('Time') ylabel(mzcdf_struct.mass_axis_label) zlabel(mzcdf_struct.intensity_axis_label)