6
class TestGaussFilter(unittest.TestCase):
9
dirname = os.path.dirname(os.path.abspath(__file__))
10
self.filename = os.path.join(dirname, "test2.mzML")
11
self.exp = pyopenms.MSExperiment()
12
pyopenms.MzMLFile().load(self.filename, self.exp)
15
thisfilter = pyopenms.GaussFilter();
18
thisfilter = pyopenms.GaussFilter();
19
old_firstspec = self.exp[0]
20
thisfilter.filterExperiment(self.exp)
22
self.assertNotEqual(self.exp.size(), 0)
23
self.assertNotEqual(old_firstspec, self.exp[0])
25
# MZ should not change, Intensity should
26
self.assertEqual(old_firstspec[10].getMZ(), self.exp[0][10].getMZ())
27
self.assertNotEqual(old_firstspec[10].getIntensity(), self.exp[0][10].getIntensity())
29
class TestSavitzkyGolayFilter(unittest.TestCase):
32
dirname = os.path.dirname(os.path.abspath(__file__))
33
self.filename = os.path.join(dirname, "test2.mzML")
34
self.exp = pyopenms.MSExperiment()
35
pyopenms.MzMLFile().load(self.filename, self.exp)
38
thisfilter = pyopenms.SavitzkyGolayFilter();
41
thisfilter = pyopenms.SavitzkyGolayFilter();
42
old_firstspec = self.exp[0]
43
thisfilter.filterExperiment(self.exp)
45
self.assertNotEqual(self.exp.size(), 0)
46
self.assertNotEqual(old_firstspec, self.exp[0])
48
# MZ should not change, Intensity should
49
self.assertEqual(old_firstspec[10].getMZ(), self.exp[0][10].getMZ())
50
self.assertNotEqual(old_firstspec[10].getIntensity(), self.exp[0][10].getIntensity())
52
class TestLowessSmoothing(unittest.TestCase):
58
thisfilter = pyopenms.LowessSmoothing();
61
thisfilter = pyopenms.LowessSmoothing();
63
y = [10.0,11.0,12.0,13.0]
65
thisfilter.smoothData(x,y,y_smoothed)
67
self.assertNotEqual( len(y_smoothed), 0)
69
# The smoothed data should be different from the input data
70
self.assertNotEqual( y_smoothed, y)
71
self.assertNotEqual( y_smoothed[0], y[0])
73
if __name__ == '__main__':