步骤1.栅格化shapefile
创建一个可以确定坐标点是否
(x,y)在该区域中的函数。有关如何将shapefile栅格化为与目标蒙版相同尺寸的数组的更多详细信息，请参见此处。
def point_is_in_mask(mask, point): # this is just pseudopre return mask.contains(point)
步骤2.创建您的面具
mask = np.zeros((height, width))value = np.zeros((height, width))for y in range(height): for x in range(width): if not point_is_in_mask(mask, (x, y)): value[y][x] = np.nan
