def Smooth(self, O):
dn = filters.gaussian_filter(O.data,self.posres,self.order,mode='wrap')
dn = np.ma.masked_array(dn, O.data.mask)
olst = [dn,
O.N,
O.norm,
O.lname,
O.sname,
O.units,
O.colmap,
O.axes,
O.rotation,
self.poshwhm]
nobs = Observable(olst)
self.Nyquist(nobs)
# Keep bounds the same as simulation resolution, for consistency.
nobs.bounds = O.bounds
self.Writer.WriteObservable(nobs)
return nobs
