Geopandas统计同覆盖小区
defsamefugei_updata(distm,agleabs):#distm:同覆盖距离，单位米；agleabs:同覆盖小区经纬度差
dis_buffer = distm
agle_abs = agleabs
file_yuan ='./原始数据\\工参表.xlsx'
dirout ='./输出结果\\'
p_yuan = pd.read_excel(file_yuan,sheet_name='主小区')
p_mubiao = pd.read_excel(file_yuan,sheet_name='目标小区')
point = gpd.GeoDataFrame(p_yuan,geometry=gpd.points_from_xy(p_yuan.经度, p_yuan.纬度),crs='EPSG:4326')
point = point.to_crs('EPSG:2381')#计算距离到米
point_mubiao = gpd.GeoDataFrame(p_mubiao,
geometry=gpd.points_from_xy(p_mubiao['经度-目标小区'], p_mubiao['纬度-目标小区']),
crs='EPSG:4326')
point_mubiao = point_mubiao.to_crs('EPSG:2381')
point_mubiao['point_2']=point_mubiao['geometry']
point_mubiao = point_mubiao.set_geometry('point_2')
point_buffer = point.buffer(dis_buffer)
plogy_1=point.copy()
plogy_1['point_buffer'] = point_buffer
plogy_1['point'] = plogy_1['geometry']
plogy_1 = plogy_1.set_geometry('point_buffer')
gp_ponitboder = gpd.sjoin(point_mubiao, plogy_1)
print(gp_ponitboder.dtypes)
gp_ponitboder['diatac(米)']=gp_ponitboder['point_2'].distance(gp_ponitboder['point'],align=False)
p_df = pd.DataFrame()
p_df = pd.concat([p_df, gp_ponitboder])
p_df['方位角差_abs'] = np.abs(p_df['方位角'] - p_df['方位角-目标小区'])
p_df['同覆盖'] = np.where(p_df['方位角差_abs'] <= agle_abs,"同覆盖","不同覆盖")
p_df['同小区']=np.where(p_df['eci'] == p_df['eci-目标小区'],"同小区","不同小区")
p_df = p_df[(p_df['同覆盖'] =='同覆盖')&(p_df['同小区']=='不同小区')]
order = ['地市','小区名称','enodbid','cellid','经度','纬度','eci','方位角','小区名称--目标小区',
'enodbid-目标小区','cellid-目标小区','经度-目标小区','纬度-目标小区','eci-目标小区','方位角-目标小区',
'方位角差_abs','同覆盖','diatac(米)']
p_df = p_df[order]
fileout = dirout +'同覆盖小区明细.csv'
p_df.to_csv(fileout,encoding='ANSI',index=False)
