步骤1.栅格化shapefile
创建一个可以确定坐标点是否(x, y)
在该区域中的函数。有关如何将shapefile栅格化为与目标蒙版相同尺寸的数组的更多详细信息,请参见此处。
def point_is_in_mask(mask, point):
# this is just pseudocode
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