对于一个不规则的shape文件,如何提取该shape文件中的坐标点/网格点?
主要思路:利用geopandas
库读取shape文件,判断点是否在区域内。
import numpy as np
import geopandas as gpd
from shapely.geometry import Point, MultiPoint, Polygon
def get_xy_in_region(shp_path, origin_xy):
'''
input
shp_path: 不规则shp文件的路径, 要求为单个多边形Polygon
orig