请问，如何利用gdal 在python下一次实现多个行政区域数据的统计

[weixin_42683927]

比如，我有全国地区级行政区化地图，有全球范围的某遥感数据，我希望利用有行政区的shp地图去切割该遥感数据(GTiff)。