diff --git a/fsl/wrappers/fslmaths.py b/fsl/wrappers/fslmaths.py
index 2bfd7129183739b79ae7af7b150f417c862e9244..21decf6f2a16823f3a31725b6c2dd561e5478f43 100644
--- a/fsl/wrappers/fslmaths.py
+++ b/fsl/wrappers/fslmaths.py
@@ -105,6 +105,11 @@ class fslmaths(object):
             self.__args.append("-dilF")
         return self
 
+    def smooth(self, sigma):
+        """Spatial smoothing - mean filtering using a gauss kernel of sigma mm"""
+        self.__args.extend(("-s", sigma))
+        return self
+
     def add(self, image):
         """Add input to current image."""
         self.__args.extend(("-add", image))
diff --git a/tests/test_wrappers/test_wrappers.py b/tests/test_wrappers/test_wrappers.py
index dbffff1d6958202bade7e104cc211e964041c0e4..e93d494eaffd46221e6efabc4d5dcb18062b69c8 100644
--- a/tests/test_wrappers/test_wrappers.py
+++ b/tests/test_wrappers/test_wrappers.py
@@ -290,14 +290,15 @@ def test_fslmaths():
             .abs().bin().binv().recip().Tmean().Tstd().Tmin().Tmax() \
             .fillh().ero().dilM().dilF().add('addim').sub('subim') \
             .mul('mulim').div('divim').mas('masim').rem('remim')   \
-            .thr('thrim').uthr('uthrim').inm('inmim').bptf(1, 10).run('output')
+            .thr('thrim').uthr('uthrim').inm('inmim').bptf(1, 10) \
+            .smooth(sigma=6).run('output')
 
         expected = [cmd, 'input',
                     '-abs', '-bin', '-binv', '-recip', '-Tmean', '-Tstd',
                     '-Tmin', '-Tmax', '-fillh', '-ero', '-dilM', '-dilF',
                     '-add addim', '-sub subim', '-mul mulim', '-div divim',
                     '-mas masim', '-rem remim', '-thr thrim', '-uthr uthrim',
-                    '-inm inmim', '-bptf 1 10', 'output']
+                    '-inm inmim', '-bptf 1 10', '-s 6', 'output']
         expected = ' '.join(expected)
 
         assert result.stdout[0] == expected