diff --git a/.gitignore b/.gitignore
index 9420924f..e1e5d52a 100644
--- a/.gitignore
+++ b/.gitignore
@@ -16,3 +16,4 @@ tools/
 zstd-1.2.0-linux.tar
 zstd-1.2.0-linux/
 .ipynb_checkpoints/
+examples/Parameters.*.txt
diff --git a/binder/requirements.txt b/binder/requirements.txt
index d0b17b3e..f217d639 100644
--- a/binder/requirements.txt
+++ b/binder/requirements.txt
@@ -1,11 +1,11 @@
-itk-elastix>=0.5.0
+itk-elastix>=0.5.1
 itk>=5.1rc1.post1
 itk-core>=5.1rc1.post1
 itk-numerics>=5.1rc1.post1
 itk-filtering>=5.1rc1.post1
 itk-registration>=5.1rc1.post1
 itk-segmentation>=5.1rc1.post1
-itkwidgets>=0.24.2
+itkwidgets>=0.25.1
 jupyterlab==1.2.5
 imageio
 ipywidgets>=7.5.1
diff --git a/examples/0_HelloRegistrationWorld.ipynb b/examples/0_HelloRegistrationWorld.ipynb
index 67f0ac47..f8ed5d3b 100644
--- a/examples/0_HelloRegistrationWorld.ipynb
+++ b/examples/0_HelloRegistrationWorld.ipynb
@@ -9,12 +9,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [],
    "source": [
     "import itk\n",
-    "from itkwidgets import view, checkerboard"
+    "from itkwidgets import view, compare, checkerboard"
    ]
   },
   {
@@ -30,72 +30,33 @@
    "source": [
     "In image registration, we typically identify the two images as the fixed and moving image. Our goal is to find the spatial transformation that makes the moving image align with the fixed image.\n",
     "\n",
-    "First, let's load our **fixed image**."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 2,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "fixed = itk.imread('data/CT_2D_head_fixed.mha', itk.F)"
+    "First, let's load our **fixed image** and the image we will align to our fixed image, the **moving image**."
    ]
   },
   {
    "cell_type": "code",
    "execution_count": 3,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "e4158972d1a04e098425a4f40d0e7c4b",
-       "version_major": 2,
-       "version_minor": 0
-      },
-      "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "view(fixed, ui_collapsed=True)"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": [
-    "And the image we will align to our fixed image, the **moving image**."
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 4,
-   "metadata": {},
    "outputs": [],
    "source": [
+    "fixed = itk.imread('data/CT_2D_head_fixed.mha', itk.F)\n",
     "moving = itk.imread('data/CT_2D_head_moving.mha', itk.F)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "7a65180dd10d4e2ba85b81ade160e3ab",
+       "model_id": "b3ebeb5e96b64508a72580cd02bf86c3",
        "version_major": 2,
        "version_minor": 0
       },
       "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=False, description='cmap'), Checkbox(v…"
       ]
      },
      "metadata": {},
@@ -103,7 +64,7 @@
     }
    ],
    "source": [
-    "view(moving, ui_collapsed=True)"
+    "compare(fixed, moving, ui_collapsed=True)"
    ]
   },
   {
@@ -115,13 +76,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "8c8f18d0f55344e7b3eb44b56a1f0bb6",
+       "model_id": "9c0aa44dab00479eb61af83594bab01a",
        "version_major": 2,
        "version_minor": 0
       },
@@ -139,7 +100,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -156,18 +117,18 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "ec232d21042e4f159272bec55c8a24a0",
+       "model_id": "30f9e21a18aa4ba3bbd0b07122785301",
        "version_major": 2,
        "version_minor": 0
       },
       "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=False, description='cmap'), Checkbox(v…"
       ]
      },
      "metadata": {},
@@ -175,18 +136,18 @@
     }
    ],
    "source": [
-    "view(registered_moving, ui_collapsed=True)"
+    "compare(fixed, registered_moving, ui_collapsed=True)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "a6e2d7c491ba4d50b45b35e03e85a475",
+       "model_id": "4aea9949acf8440d88d0e2106f41f21e",
        "version_major": 2,
        "version_minor": 0
       },
@@ -211,60 +172,28 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 10,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "fixed = itk.imread('data/CT_3D_lung_fixed.mha', itk.F)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 11,
-   "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "92cdb0e11f354bd49916557e2d6be968",
-       "version_major": 2,
-       "version_minor": 0
-      },
-      "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF3; proxy o…"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "view(fixed)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 12,
+   "execution_count": 9,
    "metadata": {},
    "outputs": [],
    "source": [
+    "fixed = itk.imread('data/CT_3D_lung_fixed.mha', itk.F)\n",
     "moving = itk.imread('data/CT_3D_lung_moving.mha', itk.F)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 13,
+   "execution_count": 12,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "b52c532b5c724faaa5d22cbe4460a504",
+       "model_id": "3435f5147c574e35b581944cb40a3177",
        "version_major": 2,
        "version_minor": 0
       },
       "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF3; proxy o…"
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=True, description='cmap'), Checkbox(va…"
       ]
      },
      "metadata": {},
@@ -272,18 +201,19 @@
     }
    ],
    "source": [
-    "view(moving)"
+    "# A human lung at different time points\n",
+    "compare(fixed, moving, link_cmap=True, link_gradient_opacity=True)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 13,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "103b7484068044fb84ed5f7911aad9d1",
+       "model_id": "b0bba114c9a44411880775350294b876",
        "version_major": 2,
        "version_minor": 0
       },
@@ -305,7 +235,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# With the default, conservative parameters, this can take a few minutes\n# on a low-compute-power system.\n",
+    "# With the default, conservative parameters, this can take a few minutes\n",
+    "# on a low-compute-power system.\n",
     "registered = itk.elastix_registration_method(fixed, moving)"
    ]
   },
diff --git a/examples/1_NumPyArray.ipynb b/examples/1_NumPyArray.ipynb
index 80fec27c..508f813e 100644
--- a/examples/1_NumPyArray.ipynb
+++ b/examples/1_NumPyArray.ipynb
@@ -9,14 +9,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [],
    "source": [
     "import itk\n",
     "import numpy as np\n",
     "import imageio\n",
-    "from itkwidgets import view, checkerboard"
+    "from itkwidgets import view, compare, checkerboard"
    ]
   },
   {
@@ -30,7 +30,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [
     {
@@ -39,7 +39,7 @@
        "numpy.ndarray"
       ]
      },
-     "execution_count": 2,
+     "execution_count": 4,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -52,31 +52,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
-   "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "78f2904bb67d47acba76aa64a819ec1c",
-       "version_major": 2,
-       "version_minor": 0
-      },
-      "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
-   "source": [
-    "view(fixed, ui_collapsed=True)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
@@ -85,7 +61,7 @@
        "numpy.ndarray"
       ]
      },
-     "execution_count": 4,
+     "execution_count": 5,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -98,18 +74,18 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "bf30bacf63324268835c2f7adde53cca",
+       "model_id": "73d7960ef9194b41823abc916f85fad9",
        "version_major": 2,
        "version_minor": 0
       },
       "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=False, description='cmap'), Checkbox(v…"
       ]
      },
      "metadata": {},
@@ -117,7 +93,7 @@
     }
    ],
    "source": [
-    "view(moving, ui_collapsed=True)"
+    "compare(fixed, moving, ui_collapsed=True)"
    ]
   },
   {
@@ -129,13 +105,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "efc57a3976414bddbf411285b209e7af",
+       "model_id": "6885b56476d74c1db2a58cedfe4da230",
        "version_major": 2,
        "version_minor": 0
       },
@@ -153,7 +129,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 11,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
@@ -162,7 +138,7 @@
        "numpy.ndarray"
       ]
      },
-     "execution_count": 11,
+     "execution_count": 8,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -175,18 +151,18 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": 9,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "cfae6e764f134fe2ac0ef124ca5175e8",
+       "model_id": "b743e0733f84453db820d32a3ad9e718",
        "version_major": 2,
        "version_minor": 0
       },
       "text/plain": [
-       "Viewer(geometries=[], gradient_opacity=0.22, point_sets=[], rendered_image=<itkImagePython.itkImageF2; proxy o…"
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=False, description='cmap'), Checkbox(v…"
       ]
      },
      "metadata": {},
@@ -194,7 +170,7 @@
     }
    ],
    "source": [
-    "view(registered, ui_collapsed=True)"
+    "compare(fixed, registered, ui_collapsed=True)"
    ]
   },
   {
@@ -212,7 +188,7 @@
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "b277130b1b2b46aa8336ee94872c9a24",
+       "model_id": "33eb1427d33a4c9293ac5c60ad74c2d3",
        "version_major": 2,
        "version_minor": 0
       },
@@ -227,6 +203,13 @@
    "source": [
     "checkerboard(fixed, registered, pattern=10)"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {
diff --git a/examples/2_RegistrationParameters.ipynb b/examples/2_RegistrationParameters.ipynb
new file mode 100644
index 00000000..3d874bd9
--- /dev/null
+++ b/examples/2_RegistrationParameters.ipynb
@@ -0,0 +1,908 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Examine, Modify, and Save Registration Parameters"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 29,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import itk\n",
+    "from itkwidgets import view, compare, checkerboard"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "[Elastix registration parameters](http://elastix.isi.uu.nl/doxygen/parameter.html) are stored in a *ParameterObject*.\n",
+    "\n",
+    "Let's find out how to examine default registration parameters, modify their values, and serialize them."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 30,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "parameters = itk.ParameterObject.New()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 31,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 24798\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 3)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Get the default parameter map for a 3 resolution translation registration\n",
+    "\n",
+    "# When multiple resolutions are used, a multi-scale image pyramid is generated with downscaled versions of the image.\n",
+    "# Registration results at a lower resolution is used to initialize registration at a higher resolution.\n",
+    "# This improves speed and robustness.\n",
+    "resolutions = 3\n",
+    "\n",
+    "default_translation = parameters.GetDefaultParameterMap(\"translation\", resolutions)\n",
+    "\n",
+    "# Add these parameters to the registration parameters\n",
+    "parameters.AddParameterMap(default_translation)\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Let's use these parameters for registration."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 32,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "fixed = itk.imread('data/CT_3D_lung_fixed.mha', itk.F)\n",
+    "moving = itk.imread('data/CT_3D_lung_moving.mha', itk.F)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 33,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "8d7dd936320443d6826241e66ffdcb0a",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "AppLayout(children=(HBox(children=(Label(value='Link:'), Checkbox(value=True, description='cmap'), Checkbox(va…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "compare(fixed, moving, link_cmap=True, link_gradient_opacity=True)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 34,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "cd073f8997944143afd3e585ce3c3c09",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "VBox(children=(Viewer(annotations=False, interpolation=False, mode='z', rendered_image=<itkImagePython.itkImag…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "checkerboard(fixed, moving, mode='z', ui_collapsed=False)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 35,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "registered = itk.elastix_registration_method(fixed, moving, parameter_object=parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 36,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "29002952ddea4cf9b6e9337c58430da6",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "VBox(children=(Viewer(annotations=False, interpolation=False, mode='z', rendered_image=<itkImagePython.itkImag…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "checkerboard(fixed, registered, mode='z', ui_collapsed=False)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "We can run the registration with multiple stages, using the spatial transformation result from the current stage to initialize registration at the next stage.\n",
+    "\n",
+    "Typically, we start with a simple spatial transformation and advance to a more complex spatial transformation.\n",
+    "\n",
+    "The default, conservative, registration parameters progress from a *translation -> affine -> bspline* transformation, with 4 resolutions for each transformation type.\n",
+    "\n",
+    "A default set of registration parameters are available for\n",
+    "\n",
+    "- *translation*\n",
+    "- *rigid*\n",
+    "- *affine*\n",
+    "- *bspline*\n",
+    "- *spline*\n",
+    "- *groupwise*\n",
+    "\n",
+    "transformations. More information on these transformation can be found in the [elastix manual](http://elastix.isi.uu.nl/download/elastix-5.0.0-manual.pdf)."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 37,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49515\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 3)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 2)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Get the default parameter map for a 3 resolution translation registration\n",
+    "resolutions = 2\n",
+    "default_affine = parameters.GetDefaultParameterMap(\"affine\", resolutions)\n",
+    "\n",
+    "# Add these parameters to the registration parameters\n",
+    "parameters.AddParameterMap(default_affine)\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 38,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49516\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 3)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 2)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 2: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FinalGridSpacingInPhysicalUnits 20)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (GridSpacingSchedule 1)\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\" \"TransformBendingEnergyPenalty\")\n",
+      "  (Metric0Weight 1)\n",
+      "  (Metric1Weight 1)\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 1)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiMetricMultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"BSplineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# For the bspline default parameter map, an extra argument can be specified that define the final bspline grid spacing in physical space. \n",
+    "resolutions = 1\n",
+    "default_bspline = parameters.GetDefaultParameterMap(\"bspline\", resolutions, 20.0)\n",
+    "\n",
+    "# Add these parameters to the registration parameters\n",
+    "parameters.AddParameterMap(default_bspline)\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 41,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49516\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 3)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 2)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 2: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FinalGridSpacingInPhysicalUnits 20)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (GridSpacingSchedule 1)\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\" \"TransformBendingEnergyPenalty\")\n",
+      "  (Metric0Weight 1)\n",
+      "  (Metric1Weight 1)\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 1)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiMetricMultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"BSplineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Change a parameter\n",
+    "\n",
+    "# Note that keys and values are currently both strings\n",
+    "parameters.SetParameter(2, \"NumberOfResolutions\", \"1\")\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 42,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49516\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 2: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FinalGridSpacingInPhysicalUnits 20)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (GridSpacingSchedule 1)\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\" \"TransformBendingEnergyPenalty\")\n",
+      "  (Metric0Weight 1)\n",
+      "  (Metric1Weight 1)\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiMetricMultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (ResultImageFormat \"nii\")\n",
+      "  (Transform \"BSplineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Change a parameter on all parameter maps\n",
+    "parameters.SetParameter(\"NumberOfResolutions\", \"5\")\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 43,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e2865a1e60)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49516\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 2: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FinalGridSpacingInPhysicalUnits 20)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (GridSpacingSchedule 1)\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\" \"TransformBendingEnergyPenalty\")\n",
+      "  (Metric0Weight 1)\n",
+      "  (Metric1Weight 1)\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiMetricMultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"BSplineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Remove a parameter\n",
+    "parameters.RemoveParameter(\"ResultImageFormat\")\n",
+    "\n",
+    "print(parameters)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 44,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "(AutomaticParameterEstimation \"true\")\n",
+      "(AutomaticScalesEstimation \"true\")\n",
+      "(CheckNumberOfSamples \"true\")\n",
+      "(DefaultPixelValue 0.000000)\n",
+      "(FinalBSplineInterpolationOrder 3.000000)\n",
+      "(FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "(ImageSampler \"RandomCoordinate\")\n",
+      "(Interpolator \"LinearInterpolator\")\n",
+      "(MaximumNumberOfIterations 256.000000)\n",
+      "(MaximumNumberOfSamplingAttempts 8.000000)\n",
+      "(Metric \"AdvancedMattesMutualInformation\")\n",
+      "(MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "(NewSamplesEveryIteration \"true\")\n",
+      "(NumberOfResolutions 5.000000)\n",
+      "(NumberOfSamplesForExactGradient 4096.000000)\n",
+      "(NumberOfSpatialSamples 2048.000000)\n",
+      "(Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "(Registration \"MultiResolutionRegistration\")\n",
+      "(ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "(Resampler \"DefaultResampler\")\n",
+      "(Transform \"AffineTransform\")\n",
+      "(WriteIterationInfo \"false\")\n",
+      "(WriteResultImage \"true\")\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Each parameter map can be serialized to a file.\n",
+    "for index in range(parameters.GetNumberOfParameterMaps()):\n",
+    "    parameter_map = parameters.GetParameterMap(index)\n",
+    "    parameters.WriteParameterFile(parameter_map, \"Parameters.{0}.txt\".format(index))\n",
+    "    \n",
+    "%cat Parameters.1.txt"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 46,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Self (0x55e287b14970)\n",
+      "  RTTI typeinfo:   elastix::ParameterObject\n",
+      "  Reference Count: 1\n",
+      "  Modified Time: 49521\n",
+      "  Debug: Off\n",
+      "  Object Name: \n",
+      "  Observers: \n",
+      "    none\n",
+      "ParameterMap 0: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticTransformInitialization \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"TranslationTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 1: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (AutomaticScalesEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\")\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"AffineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "ParameterMap 2: \n",
+      "  (AutomaticParameterEstimation \"true\")\n",
+      "  (CheckNumberOfSamples \"true\")\n",
+      "  (DefaultPixelValue 0)\n",
+      "  (FinalBSplineInterpolationOrder 3)\n",
+      "  (FinalGridSpacingInPhysicalUnits 20)\n",
+      "  (FixedImagePyramid \"FixedSmoothingImagePyramid\")\n",
+      "  (GridSpacingSchedule 1)\n",
+      "  (ImageSampler \"RandomCoordinate\")\n",
+      "  (Interpolator \"LinearInterpolator\")\n",
+      "  (MaximumNumberOfIterations 256)\n",
+      "  (MaximumNumberOfSamplingAttempts 8)\n",
+      "  (Metric \"AdvancedMattesMutualInformation\" \"TransformBendingEnergyPenalty\")\n",
+      "  (Metric0Weight 1)\n",
+      "  (Metric1Weight 1)\n",
+      "  (MovingImagePyramid \"MovingSmoothingImagePyramid\")\n",
+      "  (NewSamplesEveryIteration \"true\")\n",
+      "  (NumberOfResolutions 5)\n",
+      "  (NumberOfSamplesForExactGradient 4096)\n",
+      "  (NumberOfSpatialSamples 2048)\n",
+      "  (Optimizer \"AdaptiveStochasticGradientDescent\")\n",
+      "  (Registration \"MultiMetricMultiResolutionRegistration\")\n",
+      "  (ResampleInterpolator \"FinalBSplineInterpolator\")\n",
+      "  (Resampler \"DefaultResampler\")\n",
+      "  (Transform \"BSplineTransform\")\n",
+      "  (WriteIterationInfo \"false\")\n",
+      "  (WriteResultImage \"true\")\n",
+      "\n"
+     ]
+    }
+   ],
+   "source": [
+    "# To deserialize the parameters\n",
+    "parameter_files = [\"Parameters.{0}.txt\".format(i) for i in range(3)]\n",
+    "restored_parameters = itk.ParameterObject.New()\n",
+    "restored_parameters.ReadParameterFile(parameter_files)\n",
+    "print(restored_parameters)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "For more information on the available parameters and their meaning, see the [Elastix Parameters](http://elastix.isi.uu.nl/doxygen/parameter.html) documentation.\n",
+    "\n",
+    "A database of parameter sets that work well for specific types of data can be found in the [elastix parameter file database](http://elastix.bigr.nl/wiki/index.php/Parameter_file_database)."
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.7.6"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}