我有一张healpy全天图，需要选择一些以给定像素为中心的天空补丁。 hp.query_disc找到给定光盘中的像素，这很不错。但我希望能够在不同的天空补丁上加/减，因此不同的光盘不能相同的尺寸会导致问题。
我也用过
hp.gnomview(return_projected_map=True)
查看我的天空补丁，发现每个返回的2D阵列都有相同的大小。但是，由于我可能会在许多子图上执行操作，因此我不希望每次都绘制地图。是否有可能只获得2D阵列？如果是这样，在添加/减去这些数组后使用plt.imshow是否可以？
我觉得应该有更优雅的方式来实现我想要的东西 - 是否存在我忽略的任何现有例程？谢谢！