在python中使用shapefile剪切netCDF4 OLR数据？

Question

我获得了全球的netCDF4 OLR数据。现在，我需要印度境内的OLR值。为此，我需要使用INDIA shapefile剪辑OLR数据。我不知道该怎么剪。

有人可以帮我吗？

到目前为止，我已经编写了这段代码，但我不知道会超出此范围。

import netCDF4 as nc

import geopandas as gdp

f = nc.Dataset ('olr_2017.nc','r')

olr = f.variables ['olr'] [:]     #ARRAY SIZE WILL BE 72,144 (2.5*2.5 RESOLUTION)

lat = f.variables ['lat'][:]      # ARRAY SIZE 72

lon = f.variables ['lon'][:]      # ARRAY SIZE 144


shapefile = gdp.read_file('India.shp')

geometry = shapefile['geometry'].head()