diff --git a/fsl/transform/nonlinear.py b/fsl/transform/nonlinear.py
index 55fd48ec84264552436223dbbc8f7621c05a1017..3010cda86f246e1e8697a4cfaee412a000fc91b5 100644
--- a/fsl/transform/nonlinear.py
+++ b/fsl/transform/nonlinear.py
@@ -750,10 +750,9 @@ def applyDeformation(image, field, ref=None, order=1, mode=None, cval=None):
     if not field.sameSpace(ref):
         field = resample.resampleToReference(field,
                                              ref,
-                                             order=1,
+                                             order=order,
                                              mode='constant',
                                              cval=-1)[0]
-
     else:
         field = field.data