idl_subset_bbox.py

import ocgis
from shapely.geometry import box
from shapely.geometry.multipolygon import MultiPolygon


path = '/media/benkoziol/Extra Drive 1/data/ocgis_test_data/CanCM4/tas_day_CanCM4_decadal2000_r2i1p1_20010101-20101231.nc'
rd = ocgis.RequestDataset(path)
wrapped_state = rd.get().spatial.wrapped_state

# If the data is on a 0 to 360 coordinate system, "unwrap" the bounding box.
if wrapped_state == ocgis.constants.WrappedState.UNWRAPPED:
    geom = [160, -20, 340, 20]
# If the data is on a -180 to 180 coordinate system, break the bounding box into two elements and pass to the subsetter.
else:
    geom1 = box(*[160, -20, 180, 20])
    geom2 = box(*[-180, -20, -20, 20])
    geom = MultiPolygon([geom1, geom2])

ops = ocgis.OcgOperations(dataset=rd, geom=geom)
ret = ops.execute()