Mercurial > hg > ltpda

comparison m-toolbox/html_help/help/ug/spikeclean.html @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:f0afece42f48
1 <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
2 "http://www.w3.org/TR/1999/REC-html401-19991224/loose.dtd">
3
4 <html lang="en">
5 <head>
6 <meta name="generator" content=
7 "HTML Tidy for Mac OS X (vers 1st December 2004), see www.w3.org">
8 <meta http-equiv="Content-Type" content=
9 "text/html; charset=us-ascii">
10
11 <title>Spikes reduction in data (LTPDA Toolbox)</title>
12 <link rel="stylesheet" href="docstyle.css" type="text/css">
13 <meta name="generator" content="DocBook XSL Stylesheets V1.52.2">
14 <meta name="description" content=
15 "Presents an overview of the features, system requirements, and starting the toolbox.">
16 </head>
17
18 <body>
19 <a name="top_of_page" id="top_of_page"></a>
20
21 <p style="font-size:1px;">&nbsp;</p>
22
23 <table class="nav" summary="Navigation aid" border="0" width=
24 "100%" cellpadding="0" cellspacing="0">
25 <tr>
26 <td valign="baseline"><b>LTPDA Toolbox</b></td><td><a href="../helptoc.html">contents</a></td>
27
28 <td valign="baseline" align="right"><a href=
29 "interp.html"><img src="b_prev.gif" border="0" align=
30 "bottom" alt="Interpolating data"></a>&nbsp;&nbsp;&nbsp;<a href=
31 "gapfill.html"><img src="b_next.gif" border="0" align=
32 "bottom" alt="Data gap filling"></a></td>
33 </tr>
34 </table>
35
36 <h1 class="title"><a name="f3-12899" id="f3-12899"></a>Spikes reduction in data</h1>
37 <hr>
38
39 <p>
40 <h2>Description</h2>
41 <p>
42 Spikes in data due to different nature can be removed, if desired, from the original data. LTPDA <a href="matlab:doc('ao/spikecleaning')"><tt>spikecleaning</tt></a> detects and replaces
43 spikes of the input AOs. A spike in data is defined as a single sample exceeding a certain value (usually, the
44 floor noise of the data) defined by the user:
45 </p>
46 <br>
47 <p>
48 <div align="center">
49 <IMG src="images/spike10x.png" width="147" height="15" border="0" />
50 </div>
51 </p>
52 <br>
53 <p>
54 where <IMG src="images/spike41x.png" width="50" height="14" border="0" /> is the input data high-pass filtered, <IMG src="images/spike42x.png" width="34" height="13" border="0"/> is a value defined by the user (by default is 3.3) and <IMG src="images/spike43x.png" width="35" height="8" border="0"/> is the standard deviation of <IMG src="images/spike41x.png" width="50" height="14" border="0" /> . In consequence, a spike is defined as the value that exceeds the floor noise of the data by a factor <IMG src="images/spike42x.png" width="34" height="13" border="0"/>, the higher of this parameter the more difficult to "detect" a spike.
55 <br>
56 <br>
57 <h2>Syntax</h2>
58 </p>
59 <div class="fragment"><pre>
60 <br>
61 b = spikecleaning(a, pl)
62 </pre>
63 </div>
64 <h2>Parameters</h2>
65 <p>
66 The following parameters can be set in this method:
67 <ul>
68 <li> <tt>'kspike'</tt> - set the <IMG src="images/spike42x.png" width="34" height="13" border="0"/> value (default is 3.3) </li>
69 <li> <tt>'method'</tt> - method used to replace the "spiky" sample. Three methods are available ---see below for details---:</li>
70 <ul>
71 <li> <tt>'random'</tt> </li>
72 <li> <tt>'mean'</tt> </li>
73 <li> <tt>'previous'</tt> </li>
74 </ul>
75 <li> <tt>'fc'</tt> - frequency cut-off of the high-pass IIR filter (default is 0.025) </li>
76 <li> <tt>'order'</tt> - specifies the order of the IIR filter (default is 2) </li>
77 <li> <tt>'ripple'</tt> - specifies pass/stop-band ripple for bandpass and bandreject filters (default is 0.5) </li>
78 </ul>
79 </p>
80 <p>
81 <h2>Algorithm</h2>
82 </p>
83 <p>
84 <b>Random: </b> this method substitutes the spiky sample by:
85 </p>
86 <p>
87 <br>
88 <div align="center">
89 <IMG src="images/spike20x.png" width="197" height="15" align="center" border="0"/>
90 </div>
91 <br>
92 </p>
93 <p>
94 where <IMG src="images/spike44x.png" width="42" height="14" border="0"/> is a random number of mean zero and standard deviation 1.
95 </p>
96 <p>
97 <b>Mean: </b>this method uses the following equation to replace the spike detected in data.
98 </p>
99 <br>
100 <p>
101 <div align="center">
102 <IMG src="images/spike30x.png" width="158" height="29" align="center" border="0"/>
103 </div>
104 </p>
105 <br>
106 <p>
107 <b>Previous: </b>the spike is substitued by the previous sample, i.e.:
108 </p>
109 <br>
110 <p>
111 <div align="center">
112 <IMG src="images/spike40x.png" width="90" height="14" align="center" border="0"/>
113 </div>
114 </p>
115 <br>
116 <p>
117 <h2>Examples</h2>
118 </p>
119 <p>
120 1. Spike cleaning of a sequence of random data with <tt>kspike = 2</tt>.
121 </p>
122 <p>
123 <div class="fragment"><pre>
124 <br>
125 x = ao(plist( <span class="string">'waveform'</span>, <span class="string">'noise'</span>, <span class="string">'nsecs'</span>,1e4, <span class="string">'fs'</span>,10)); <span class="comment">% create an AO of random data sampled at 1 Hz.</span>
126 pl = plist( <span class="string">'kspike'</span>, 2); <span class="comment">% kspike = 2</span>
127 y = spikecleaning(x, pl); <span class="comment">% spike cleaning function applied to the input AO, x</span>
128 iplot(x, y) <span class="comment">% plot original and "cleaned" data</span>
129 </pre>
130 </div>
131 </p>
132 <p>
133 <div align="center">
134 <img src="images/spike1.png" border="1" >
135 </div>
136 </p>
137 <!-- <p>
138 2. Example of real data: the first image shows data from the real world prior to
139 the application of the spike cleaning tool. It is clear that some spikes are
140 present in data and might be convenient to remove them. The second image shows
141 the same data after the spike samples supression.
142 </p>
143 <p>
144 <img src="images/stat1x.png" border="1" width="400px">
145 </p>
146 <p>
147 <img src="images/stat2x.png" border="1" width="400px">
148 </p>
149 -->
150
151
152
153 </p>
154
155 <br>
156 <br>
157 <table class="nav" summary="Navigation aid" border="0" width=
158 "100%" cellpadding="0" cellspacing="0">
159 <tr valign="top">
160 <td align="left" width="20"><a href="interp.html"><img src=
161 "b_prev.gif" border="0" align="bottom" alt=
162 "Interpolating data"></a>&nbsp;</td>
163
164 <td align="left">Interpolating data</td>
165
166 <td>&nbsp;</td>
167
168 <td align="right">Data gap filling</td>
169
170 <td align="right" width="20"><a href=
171 "gapfill.html"><img src="b_next.gif" border="0" align=
172 "bottom" alt="Data gap filling"></a></td>
173 </tr>
174 </table><br>
175
176 <p class="copy">&copy;LTP Team</p>
177 </body>
178 </html>