主题：【第十二届原创】一个较难积分的例子

一个较难积分的例子

有时候可能会遇到峰形比较特殊难以积分的情况。例如下面这个例子，前面大峰A一直延申到峰E，呈U型状。A峰和E峰中间的拖尾部分上面还有三个很小的峰b, c, d。根据质谱离子分析，积分切线大致按红线划线处理比较好些。但怎样能够得到这样的积分呢？

图1 待处理的峰的色谱图

*****************************************************************************
1普通简单处理
首先用一般普通的积分参数处理试试。处理结果非常差，误差太大，是无法使用的。

图2 按一般积分参数1处理的效果

***********************************************************************
改变灵敏度，初始阈值由原来的16变成17，积分效果如下，还是不行，达不到要求：

图3 按一般积分参数2处理的效果

*********************************************************************************
再进行撇线处理还是不行。

图4 一般峰拖尾撇线处理效果

****************************************************************************
2 手动积分处理
先对主峰上面的小峰和后面两个较大的峰进行手动积分（大峰先手动积分，小峰就不好手动积分了）。

图5 手动积分处理1

****************************************************************************
然后对大峰进行手动积分。但U型中间的第一个小峰的手动积分已经看不到了。需要用前面的手动积分结果来补充一下和扣除一下。

图6 手动积分处理2

*************************************************************************
用手动积分基本达到了要求。各个峰的面积如下：
A----249203-415614=2079833589
  b----415614
  c----1092597
  d----823251
E----9256764
F----29259394
3 撇线处理，阈值变化，基线回归，峰面积限制和峰面积加和处理
后面采用了“面积加和”与“撇线处理”的方式等积分参数尝试改变积分效果，但只能是部分峰可以得到较好的效果，无法得到所有峰处理的理想结果。后来再进一步尝试撇线处理，阈值变化，基线回归和峰面积加和等参数改变得到如下效果。

图 7 撇线处理和基线回归面积限制处理

********************************************************************************
上面的初始面积限制在50000，过滤了一部分误判的色谱峰；初始峰宽0.02较为合；在13.8分钟进行撇线处理；在14.7分钟到15分钟改变阈值14到17，14可以积分到14.906的峰，过大就无法积分到了；在16.9，17.13和17,4分钟分别让基线归零，使得17.06和17.221两个峰能够得到较好的结果。积分数值如下：
peak R.T. first  max last  PK  peak      corr.
  #  min  scan scan scan  TY  height    area
--- ----- ----- ---- ---- ---  -------  -------
3  13.875  1349 1386 1464 BB  17413149 1675572167
4  14.806  1503 1513 1524 BV 9    9702  280597
5  16.174  1669 1701 1712 PV 4  29585  1257201
6  16.403  1721 1732 1744 PB 4  26523  936022
7  17.061  1809 1823 1832 BBA  2638932 96395299
8  17.221  1832 1845 1870 BBA  763290 28452100
然后从13.6-17.5分钟进行面积加和得到总面积。

图 8 峰面积加和

************************************************************************
峰面积加和得到总面积为2072469456。
peak R.T. first  max last  PK  peak      corr.
  #  min  scan scan scan  TY  height    area
--- ----- ----- ---- ---- --- -------  -------
3 15.553  1348 1616 1883 BB  16984219 2072469456
总面积为2072469456。利用总面积和和b, c, d, E, F的面积就可以计算到特殊峰形A峰的面积。
可能第二种方案得到的结果更为准确一些，但花费的功夫可能要多一些。