Add a +proj=defmodel transformation for multi-component time-based deformation models

Fixes #1001

Co-authored-by: Chris Crook <ccrook@linz.govt.nz>

17 files changed, 506 insertions, 0 deletions

diff --git a/data/tests/simple_model_degree_3d.json b/data/tests/simple_model_degree_3d.json
new file mode 100644
index 00000000..6a41ce5c
--- /dev/null
+++ b/data/tests/simple_model_degree_3d.json
@@ -0,0 +1,54 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "degree",
+  "horizontal_offset_method": "addition",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "3d",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_degree_3d_grid.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_degree_3d_grid.tif b/data/tests/simple_model_degree_3d_grid.tif
new file mode 100644
index 00000000..3bbff0a6
--- /dev/null
+++ b/data/tests/simple_model_degree_3d_grid.tif
diff --git a/data/tests/simple_model_degree_horizontal.json b/data/tests/simple_model_degree_horizontal.json
new file mode 100644
index 00000000..533c0054
--- /dev/null
+++ b/data/tests/simple_model_degree_horizontal.json
@@ -0,0 +1,53 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "degree",
+  "horizontal_offset_method": "addition",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "horizontal",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_degree_3d_grid.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_metre_3d.json b/data/tests/simple_model_metre_3d.json
new file mode 100644
index 00000000..201aaebb
--- /dev/null
+++ b/data/tests/simple_model_metre_3d.json
@@ -0,0 +1,54 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "metre",
+  "horizontal_offset_method": "addition",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "3d",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_metre_3d_grid.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_metre_3d_geocentric.json b/data/tests/simple_model_metre_3d_geocentric.json
new file mode 100644
index 00000000..1328ff5e
--- /dev/null
+++ b/data/tests/simple_model_metre_3d_geocentric.json
@@ -0,0 +1,54 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "metre",
+  "horizontal_offset_method": "geocentric",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "3d",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_metre_3d_grid.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_metre_3d_grid.tif b/data/tests/simple_model_metre_3d_grid.tif
new file mode 100644
index 00000000..40cf2d70
--- /dev/null
+++ b/data/tests/simple_model_metre_3d_grid.tif
diff --git a/data/tests/simple_model_metre_horizontal.json b/data/tests/simple_model_metre_horizontal.json
new file mode 100644
index 00000000..d0ac477c
--- /dev/null
+++ b/data/tests/simple_model_metre_horizontal.json
@@ -0,0 +1,53 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "metre",
+  "horizontal_offset_method": "addition",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "horizontal",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_metre_3d_grid.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_metre_vertical.json b/data/tests/simple_model_metre_vertical.json
new file mode 100644
index 00000000..70574340
--- /dev/null
+++ b/data/tests/simple_model_metre_vertical.json
@@ -0,0 +1,49 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        -180,
+        -90,
+        180,
+        90
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "vertical",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            -180,
+            -90,
+            180,
+            90
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/simple_model_metre_vertical_grid.tif"
+      },
+      "time_function": {
+        "type": "constant"
+      }
+    }
+  ]
+}
\ No newline at end of file
diff --git a/data/tests/simple_model_metre_vertical_grid.tif b/data/tests/simple_model_metre_vertical_grid.tif
new file mode 100644
index 00000000..058b8081
--- /dev/null
+++ b/data/tests/simple_model_metre_vertical_grid.tif
diff --git a/data/tests/simple_model_polar.json b/data/tests/simple_model_polar.json
new file mode 100644
index 00000000..ef99a0cb
--- /dev/null
+++ b/data/tests/simple_model_polar.json
@@ -0,0 +1,54 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "horizontal_offset_unit": "metre",
+  "horizontal_offset_method": "geocentric",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        0,
+        -90,
+        360,
+        -89
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "3d",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            0,
+            -90,
+            360,
+            -89
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "geocentric_bilinear",
+        "filename": "tests/simple_model_polar.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
diff --git a/data/tests/simple_model_polar.tif b/data/tests/simple_model_polar.tif
new file mode 100644
index 00000000..7371ca1e
--- /dev/null
+++ b/data/tests/simple_model_polar.tif
diff --git a/data/tests/simple_model_projected.json b/data/tests/simple_model_projected.json
new file mode 100644
index 00000000..c97a7c11
--- /dev/null
+++ b/data/tests/simple_model_projected.json
@@ -0,0 +1,51 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:2193",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:2193",
+  "horizontal_offset_unit": "metre",
+  "horizontal_offset_method": "addition",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [
+        1500000.0,
+        5400000.0,
+        1501000.0,
+        5401000.0
+      ]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "3d",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [
+            1500000.0,
+            5400000.0,
+            1501000.0,
+            5401000.0
+          ]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "bilinear",
+        "filename": "tests/test_3d_grid_projected.tif"
+      },
+      "time_function": {
+        "type": "constant"
+      }
+    }
+  ]
+}
diff --git a/data/tests/simple_model_wrap_east.json b/data/tests/simple_model_wrap_east.json
new file mode 100644
index 00000000..2a0a2c0e
--- /dev/null
+++ b/data/tests/simple_model_wrap_east.json
@@ -0,0 +1,42 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [-194.2, -37.5, -193.8, -37.2]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "vertical",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [-194.2, -37.5, -193.8, -37.2]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "geocentric_bilinear",
+        "filename": "tests/simple_model_wrap_east.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
diff --git a/data/tests/simple_model_wrap_east.tif b/data/tests/simple_model_wrap_east.tif
new file mode 100644
index 00000000..816d8a7a
--- /dev/null
+++ b/data/tests/simple_model_wrap_east.tif
diff --git a/data/tests/simple_model_wrap_west.json b/data/tests/simple_model_wrap_west.json
new file mode 100644
index 00000000..54a04bd2
--- /dev/null
+++ b/data/tests/simple_model_wrap_west.json
@@ -0,0 +1,42 @@
+{
+  "file_type": "deformation_model_master_file",
+  "format_version": "1.0",
+  "source_crs": "EPSG:4326",
+  "target_crs": "foo:ignored",
+  "definition_crs": "EPSG:4326",
+  "vertical_offset_unit": "metre",
+  "extent": {
+    "type": "bbox",
+    "parameters": {
+      "bbox": [525.8,-37.5,526.2,-37.2]
+    }
+  },
+  "time_extent": {
+    "first": "1900-01-01T00:00:00Z",
+    "last": "2050-01-01T00:00:00Z"
+  },
+  "components": [
+    {
+      "description": "test",
+      "displacement_type": "vertical",
+      "uncertainty_type": "none",
+      "extent": {
+        "type": "bbox",
+        "parameters": {
+          "bbox": [525.8,-37.5,526.2,-37.2]
+        }
+      },
+      "spatial_model": {
+        "type": "GeoTIFF",
+        "interpolation_method": "geocentric_bilinear",
+        "filename": "tests/simple_model_wrap_west.tif"
+      },
+      "time_function": {
+        "type": "step",
+        "parameters": {
+          "step_epoch": "1900-01-01T00:00:00Z"
+        }
+      }
+    }
+  ]
+}
diff --git a/data/tests/simple_model_wrap_west.tif b/data/tests/simple_model_wrap_west.tif
new file mode 100644
index 00000000..3a8da6f6
--- /dev/null
+++ b/data/tests/simple_model_wrap_west.tif
diff --git a/data/tests/test_3d_grid_projected.tif b/data/tests/test_3d_grid_projected.tif
index e19b1b0e..56138417 100644
--- a/data/tests/test_3d_grid_projected.tif
+++ b/data/tests/test_3d_grid_projected.tif