correctly pass args in recon script

mu_map/recon/recon.py

@@ -69,16 +69,16 @@ OSMAPOSLParameters :=
   ; warning: do not normally use together with a prior
   ;inter-iteration filter subiteration interval := 4
   ;inter-iteration filter type := Separable Gaussian
-  post-filter type := Separable Gaussian
-  separable gaussian filter parameters :=
-    x-dir filter fwhm (in mm) := 6
-    y-dir filter fwhm (in mm) := 6
-    z-dir filter fwhm (in mm) := 6
-    x-dir maximum kernel size := 129
-    y-dir maximum kernel size := 129
-    z-dir maximum kernel size := 31
-    Normalise filter to 1 := 1
-  end separable gaussian filter parameters :=
+  ;post-filter type := Separable Gaussian
+  ;separable gaussian filter parameters :=
+    ;x-dir filter fwhm (in mm) := 6
+    ;y-dir filter fwhm (in mm) := 6
+    ;z-dir filter fwhm (in mm) := 6
+    ;x-dir maximum kernel size := 129
+    ;y-dir maximum kernel size := 129
+    ;z-dir maximum kernel size := 31
+    ;Normalise filter to 1 := 1
+  ;end separable gaussian filter parameters :=
 
 END :=
 """
@@ -112,6 +112,7 @@ def reconstruct(
     init: Optional[Tuple[Dict[str, str], np.ndarray]] = None,
     n_subsets: Optional[int] = 4,
     n_iterations: Optional[int] = 10,
+    **kwargs,
 ):
     # sanitize parameters
     n_subiterations = n_subsets * n_iterations
@@ -222,7 +223,11 @@ if __name__ == "__main__":
     else:
         projection = load_as_interfile(args.projection)
 
-    header, image = reconstruct(projection, mu_map=mu_map)
+    kwargs = vars(args)
+    del kwargs["mu_map"]
+    del kwargs["projection"]
+    print(kwargs)
+    header, image = reconstruct(projection, mu_map=mu_map, **kwargs)
     image = image[:, :, ::-1]
 
     write_interfile(args.out, header, image)