Add python scripts to run registration on the cluster

scripts/extension/registration/apply_registration.py

@@ -98,9 +98,16 @@ def apply_for_file(input_path, output_path,
                    transformation_file, fiji_executable,
                    elastix_directory, tmp_folder, n_threads):
 
-    # grrr design anti-pattern
+    assert os.path.exists(elastix_directory)
+    assert os.path.exists(tmp_folder)
+    assert os.path.exists(input_path)
+    assert os.path.exists(transformation_file)
+    assert os.path.exists(os.path.split(output_path)[0])
+
+    # transformix arguments need to be passed as one string,
+    # with individual arguments comma separated
     # the argument to transformaix needs to be one large comma separated string
-    transformix_argument = ["elastix_directory=\'%s\'" % elastix_directory,
+    transformix_argument = ["elastixDirectory=\'%s\'" % elastix_directory,
                             "workingDirectory=\'%s\'" % tmp_folder,
                             "inputImageFile=\'%s\'" % input_path,
                             "transformationFile=\'%s\'" % transformation_file,
@@ -116,15 +123,26 @@ def apply_for_file(input_path, output_path,
     # "elastixDirectory='/g/almf/software/elastix_v4.8', workingDirectory='$TMPDIR',
     # inputImageFile='$INPUT_IMAGE',transformationFile='/g/cba/exchange/platy-trafos/linear/TransformParameters.BSpline10-3Channels.0.txt
     # outputFile='$OUTPUT_IMAGE',outputModality='Save as BigDataViewer .xml/.h5',numThreads='1'"
-    # NOTE: I ommit --run here, because fiji throws a warning that it does not recognise the argument
-    cmd = [fiji_executable, "--ij2", "--headless", "\"Transformix\"", transformix_argument]
+    cmd = [fiji_executable, "--ij2", "--headless", "--run", "\"Transformix\"", transformix_argument]
 
     cmd_str = " ".join(cmd)
     fu.log("Calling the following command:")
     fu.log(cmd_str)
 
+    # the elastix wrapper only works properly if we set these as environment variables as well, see
+    # TODO make issue about this
+    # os.environ['TMPDIR'] = tmp_folder
+    os.environ['TRAFO'] = transformation_file
+
     try:
-        check_output(cmd)
+        # check_output(cmd)
+        # the CLI parser is very awkward (to put it nicely).
+        # I could only get it to work by passing the whole command string
+        # and setting shell to True.
+        # otherwise, it would parse something wrong, and do nothing but
+        # throwing a warning:
+        # [WARNING] Ignoring invalid argument: --run
+        check_output([cmd_str], shell=True)
     except CalledProcessError as e:
         raise RuntimeError(e.output)