raise errors when trying to request local grids on wrong grid types

polytope/datacube/transformations/datacube_mappers/mapper_types/healpix.py

@@ -6,14 +6,16 @@
 
 class HealpixGridMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
         self._resolution = resolution
         self._axis_reversed = {mapped_axes[0]: True, mapped_axes[1]: False}
         self._first_axis_vals = self.first_axis_vals()
         self.compressed_grid_axes = [self._mapped_axes[1]]
 
+        if local_area != []:
+            raise NotImplementedError("Local area grid not implemented for healpix grids")
+
     def first_axis_vals(self):
         rad2deg = 180 / math.pi
         vals = [0] * (4 * self._resolution - 1)

polytope/datacube/transformations/datacube_mappers/mapper_types/healpix_nested.py

@@ -6,7 +6,6 @@
 
 class NestedHealpixGridMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
         self._resolution = resolution
@@ -18,6 +17,9 @@ def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
         self.Npix = 12 * self.Nside * self.Nside
         self.Ncap = (self.Nside * (self.Nside - 1)) << 1
 
+        if local_area != []:
+            raise NotImplementedError("Local area grid not implemented for healpix grids")
+
     def first_axis_vals(self):
         rad2deg = 180 / math.pi
         vals = [0] * (4 * self._resolution - 1)

polytope/datacube/transformations/datacube_mappers/mapper_types/local_regular.py

@@ -5,9 +5,12 @@
 
 class LocalRegularGridMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
+
+        if local_area == [] or len(local_area) != 4:
+            raise TypeError("Local area grid region not or wrongly specified")
+
         self._first_axis_min = local_area[0]
         self._first_axis_max = local_area[1]
         self._second_axis_min = local_area[2]

polytope/datacube/transformations/datacube_mappers/mapper_types/octahedral.py

@@ -6,7 +6,6 @@
 
 class OctahedralGridMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
         self._resolution = resolution
@@ -16,6 +15,9 @@ def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
         self._axis_reversed = {mapped_axes[0]: True, mapped_axes[1]: False}
         self.compressed_grid_axes = [self._mapped_axes[1]]
 
+        if local_area != []:
+            raise NotImplementedError("Local area grid not implemented for octahedral grids")
+
     def gauss_first_guess(self):
         i = 0
         gvals = [

polytope/datacube/transformations/datacube_mappers/mapper_types/reduced_ll.py

@@ -5,14 +5,16 @@
 
 class ReducedLatLonMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
         self._resolution = resolution
         self._axis_reversed = {mapped_axes[0]: False, mapped_axes[1]: False}
         self._first_axis_vals = self.first_axis_vals()
         self.compressed_grid_axes = [self._mapped_axes[1]]
 
+        if local_area != []:
+            raise NotImplementedError("Local area grid not implemented for reduce lat-lon grids")
+
     def first_axis_vals(self):
         resolution = 180 / (self._resolution - 1)
         vals = [-90 + i * resolution for i in range(self._resolution)]

polytope/datacube/transformations/datacube_mappers/mapper_types/regular.py

@@ -5,7 +5,6 @@
 
 class RegularGridMapper(DatacubeMapper):
     def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
-        # TODO: if local area is not empty list, raise NotImplemented
         self._mapped_axes = mapped_axes
         self._base_axis = base_axis
         self._resolution = resolution
@@ -14,6 +13,9 @@ def __init__(self, base_axis, mapped_axes, resolution, local_area=[]):
         self._first_axis_vals = self.first_axis_vals()
         self.compressed_grid_axes = [self._mapped_axes[1]]
 
+        if local_area != []:
+            raise TypeError("Use local_regular grid type for local area regular lat-lon grids")
+
     def first_axis_vals(self):
         first_ax_vals = [90 - i * self.deg_increment for i in range(2 * self._resolution)]
         return first_ax_vals

polytope/options.py

@@ -66,7 +66,6 @@ class Config(ConfigModel):
 class PolytopeOptions(ABC):
     @staticmethod
     def get_polytope_options(options):
-
         parser = argparse.ArgumentParser(allow_abbrev=False)
         conflator = Conflator(app_name="polytope", model=Config, cli=False, argparser=parser, **options)
         config_options = conflator.load()