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