python两个shp面矢量之间求交集（类似arcgis中相交）

需要特别注意geopandas包的安装，可以参考博客：手把手教学如何安装geopandas

import geopandas as gpd 
def getIntersec(shpfile1,shpfile2,outshp):
    shp1 = gpd.read_file(shpfile1)
    shp2 = gpd.read_file(shpfile2)
    intersec = gpd.overlay(shp1,shp2,how="intersection")
    intersec.to_file(outshp,encoding="utf-8")
if __name__ == "__main__":
	shpfile1=r"D:/1.shp"
	shpfile2=r"D:/2.shp"
	outshp=r"D:/jiaoji.shp"
	getIntersec(shpfile1,shpfile2,outshp)