File: TOPP_example_signalprocessing.html

package info (click to toggle)
openms 1.11.1-5
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 436,688 kB
  • ctags: 150,907
  • sloc: cpp: 387,126; xml: 71,547; python: 7,764; ansic: 2,626; php: 2,499; sql: 737; ruby: 342; sh: 325; makefile: 128
file content (79 lines) | stat: -rw-r--r-- 5,156 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
<HTML>
<HEAD>
<TITLE>Profile data processing</TITLE>
<LINK HREF="doxygen.css" REL="stylesheet" TYPE="text/css">
<LINK HREF="style_ini.css" REL="stylesheet" TYPE="text/css">
</HEAD>
<BODY BGCOLOR="#FFFFFF">
<A href="index.html">Home</A> &nbsp;&middot;
<A href="classes.html">Classes</A> &nbsp;&middot;
<A href="annotated.html">Annotated Classes</A> &nbsp;&middot;
<A href="modules.html">Modules</A> &nbsp;&middot;
<A href="functions_func.html">Members</A> &nbsp;&middot;
<A href="namespaces.html">Namespaces</A> &nbsp;&middot;
<A href="pages.html">Related Pages</A>
<HR style="height:1px; border:none; border-top:1px solid #c0c0c0;">
<!-- Generated by Doxygen 1.8.5 -->
</div><!-- top -->
<div class="header">
  <div class="headertitle">
<div class="title">Profile data processing </div>  </div>
</div><!--header-->
<div class="contents">
<div class="textblock"><p><b>Goal:</b> You want to find all peaks in your profile data.</p>
<p>The first step shown here is the elimination of noise using a <b>NoiseFilter</b>. The now smoothed profile data can be further processed by subtracting the baseline with the <b>BaselineFilter</b>. Then use one of the <b>PeakPickers</b> to find all peaks in the baseline-reduced profile data.</p>
<div class="image">
<img src="TOPP_raw_data.png" alt="TOPP_raw_data.png"/>
</div>
 <p>We offer two different smoothing filters: NoiseFilterGaussian and NoiseFilterSGolay. If you want to use the Savitzky Golay filter, or our <b>BaselineFilter</b> with non equally spaced profile data, e.g. TOF data, you have to generate equally spaced data using the <b>Resampler</b> tool.</p>
<h1><a class="anchor" id="TOPP_example_signalprocessing_peakpicker"></a>
Picking peaks with a PeakPicker</h1>
<p>The <b>PeakPicker</b> tools allow for picking peaks in profile data. Currently, there are two different TOPP tools available, PeakPickerWavelet and PeakPickerHiRes.</p>
<table  border="1">
<tr>
<td><b>PeakPickerWavelet</b> &#160;&#160; <b>Input data:</b> profile data (low/medium resolution)   </td></tr>
<tr>
<td><p class="starttd"><b>Description:</b> <br/>
 This peak picking algorithm uses the continuous wavelet transform of a raw data signal to detect mass peaks. Afterwards a given asymmetric peak function is fitted to the raw data and important peak parameters (e.g. fwhm) are extracted. In an optional step these parameters can be optimized using a non-linear opimization method. <br/>
</p>
<p>The algorithm is descripted in detail in Lange et al. (2006) Proc. PSB-06.</p>
<p><b>Application:</b> <br/>
 This algorithm was designed for low and medium resolution data. It can also be applied to high-resolution data, but can be slow on large datasets.<br/>
</p>
<p class="endtd">See the PeakPickerCWT class documentation for a parameter list.   </p>
</td></tr>
</table>
<table  border="1">
<tr>
<td><b>PeakPickerHiRes</b> &#160;&#160; <b>Input data:</b> profile data (high resolution)   </td></tr>
<tr>
<td><p class="starttd"><b>Description:</b> <br/>
 This peak-picking algorithm detects ion signals in raw data and reconstructs the corresponding peak shape by cubic spline interpolation. Signal detection depends on the signal-to-noise ratio which is adjustable by the user (see parameter <em>signal_to_noise</em>). A picked peak's m/z and intensity value is given by the maximum of the underlying peak spline. Please notice that this method is still <b>experimental</b> since it has not been tested thoroughly yet.<br/>
</p>
<p><b>Application:</b> <br/>
The algorithm is best suited for high-resolution MS data (FT-ICR-MS, Orbitrap). In high-resolution data, the signals of ions with similar mass-to-charge ratios (m/z) exhibit little or no overlapping and therefore allow for a clear separation. Furthermore, ion signals tend to show well-defined peak shapes with narrow peak width. These properties faciliate a fast computation of picked peaks so that even large data sets can be processed very quickly.</p>
<p class="endtd">See the PeakPickerHiRes class documentation for a parameter list.   </p>
</td></tr>
</table>
<h1><a class="anchor" id="TOPP_example_signalprocessing_parameters"></a>
Finding the right parameters for the</h1>
<p>NoiseFilters, the BaselineFilter and the PeakPickers</p>
<p>Finding the right parameters is not trivial. The default parameters will not work on most datasets. In order to find good parameters, we propose the following procedure:</p>
<ol type="1">
<li>Load the data in TOPPView</li>
<li>Extract a single scan from the middle of the HPLC gradient (Right click on scan)</li>
<li>Experiment with the parameters until you have found the proper settings</li>
</ol>
<ul>
<li>You can find the <b>NoiseFilters</b>, the <b>BaselineFilter</b>, and the <b>PeakPickers</b> in <b>TOPPView</b> in the menu 'Layer' - 'Apply TOPP tool' </li>
</ul>
</div></div><!-- contents -->
<HR style="height:1px; border:none; border-top:1px solid #c0c0c0;">
<TABLE width="100%" border="0">
<TR>
<TD><font color="#c0c0c0">OpenMS / TOPP release 1.11.1</font></TD>
<TD align="right"><font color="#c0c0c0">Documentation generated on Thu Nov 14 2013 11:19:25 using doxygen 1.8.5</font></TD>
</TR>
</TABLE>
</BODY>
</HTML>