def fitDoppler(self):
"""Read in a Doppler file, remove outliers and then perform a quadratic fit"""
self._wrap()
unw = self._unwrap()
self.pixelIndex = range(len(self.fractionalDoppler))
(self.linear['a'], self.linear['b'], self.linear['stdDev']) = MathModule.linearFit(self.pixelIndex, unw)
(pixels, unw) = self._cullPoints(self.pixelIndex,unw)
(self.linear['a'], self.linear['b'], self.linear['stdDev']) = MathModule.linearFit(pixels, unw)
(pixels, unw) = self._cullPoints(pixels,unw)
(a,b,c) = MathModule.quadraticFit(pixels,unw)
self.quadratic['a'] = a
self.quadratic['b'] = b
self.quadratic['c'] = c
def parameterize(self):
"""
Read in Doppler estimates, remove outliers and then perform a quadratic fit
"""
self.pixelIndex = range(len(self.fractionalDoppler))
(self.linear['a'], self.linear['b'],
self.linear['stdDev']) = MathModule.linearFit(self.pixelIndex,
self.unw)
(pixels, self.unw) = self._cullPoints(self.pixelIndex, self.unw)
(self.linear['a'], self.linear['b'],
self.linear['stdDev']) = MathModule.linearFit(pixels, self.unw)
(pixels, self.unw) = self._cullPoints(pixels, self.unw)
(a, b, c) = MathModule.quadraticFit(pixels, self.unw)
self.quadratic['a'] = a
self.quadratic['b'] = b
self.quadratic['c'] = c