如何通过gdal或Python中的其他模块使用shapefile剪辑.hdf？

Question

有一个文件（* .hdf，或者我们可以调用一个栅格），我可以通过gdal ReadAsArray()在Python中将其作为数组打开。然后，我需要用shapefile剪辑它并保留其余部分（B - A ∩ B）。我随处可见，它可能很有用：clip a raster with shapefile。

但是有很多问题，我只是不知道我需要什么模块或功能。