how to clip to area and mask by vector bounds and values

align_clip_and_mask.py

vector = gdal.OpenEx(vector_path, gdal.OF_VECTOR)
layer = vector.GetLayer()
layer.SetAttributeFilter('Country="USA"') # or whatever...

feature = layer.GetNextFeature()
geom = feature.GetGeometryRef()
envelope = geom.GetEnvelope()
target_bb = [envelope[i] for i in [0, 2, 1, 3]])
geom = None
feature = None
layer = None
vector = None

pygeoprocessing.align_and_resize_raster_stack(
    raster_stack, target_stack, ['near']*len(raster_stack),
    target_bb, vector_mask_options={
        'mask_vector_path': vector_path,
        'mask_vector_where_filter': 'Country="USA"',
        })