1.读取shp,并按照属性提取特定的region代码:
#read region shp
path = './shp/'
lake_gdf = gpd.read_file(os.path.join(path,'region.shp'))
region = lake_gdf[lake_gdf.name == 'Great Slave Lake']
2.total_bounds:限定图像的画幅范围,其返回依次记录了整列矢量数据所在最小矩形区域左下角x、左下角y、右上角x以及右上角y的numpy数组。
#get the outer envelop rectangular coordinates.
#minx, miny, maxx, maxy
region_lat = [float(region.total_bounds[1]), float(region.total_bounds[3])]
region_lon = [float(region.total_bounds[0]), float(region.total_bounds[2])]
print(region_lat,region_lon)