矢量裁剪矢量

也不知道为啥，矢量裁剪矢量这么普通的东西这么难找，赶紧放出来让大家用用

import os
import numpy as np
import geopandas as gpd
import warnings
warnings.filterwarnings('ignore', 'GeoSeries.notna', UserWarning)


shp_path = './osm/gis_osm_roads_free_1.shp'    #被裁剪的对象
mask = 'D:/ToolSet/Test_data2/mask.shp'        #用于裁剪的范围

shp_gpd = gpd.read_file(shp_path)
mask_gpd = gpd.read_file(mask)
print("shp_gpd", shp_gpd.crs)

mask_clip = gpd.clip(shp_gpd, mask_gpd)

mask_clip = mask_clip[~mask_clip.is_empty]
mask_clip.to_file("./osm/osm_road_anji.shp")   #存储结果