def test_get_spatial_subset_circular_geometries(self):
"""Test circular geometries. They were causing wrapping errors."""
geoms = TestGeom.get_geometry_dictionaries()
rd = self.test_data.get_rd('cancm4_tas')
ss = SpatialSubsetOperation(rd, wrap=True)
buffered = [element['geom'].buffer(rd.get().spatial.grid.resolution*2) for element in geoms]
for buff in buffered:
record = [{'geom': buff, 'properties': {'UGID': 1}}]
subset_sdim = SpatialDimension.from_records(record)
ret = ss.get_spatial_subset('intersects', subset_sdim)
self.assertTrue(np.all(ret.spatial.grid.extent > 0))
def get_subset_operation(self):
geom = TestGeom.get_geometry_dictionaries()
rd = self.test_data.get_rd('cancm4_tas')
ops = ocgis.OcgOperations(dataset=rd, geom=geom, select_nearest=True)
subset = SubsetOperation(ops)
return subset