通信人家园

标题: 网优业余学习编程：Python小工具--站间距计算工具 [查看完整版帖子] [打印本页]

时间: 2020-2-29 15:50

作者: ufo0373 标题: 网优业余学习编程：Python小工具--站间距计算工具

本帖最后由 ufo0373 于 2020-2-29 15:51 编辑

背景

目的一：解决不同网络站点之间的距离计算，比如要计算全网GSM共站址的LTE站点；

目的二：批量计算不同场景周边的站点，比如一个城市各居民区的覆盖站点清单。

工具源代码

import openpyxl,xlrd
import os,sys
import math,re
import datetime
#-----------------------数据输出到EXCEL文件------------------------
def write_xlsx(dis):
workbook = openpyxl.Workbook()
sheet = workbook.active
sheet.title = '最小站间距'
for i in range(0, len(dis)):
      for j in range(0,len(dis)):
         sheet.cell(row=i+1, column=j+1, value=dis[j])
workbook.save('最小站间距-'+datetime.datetime.now().strftime('%Y%m%d')+'.xlsx')

#-----------------------程序入口------------------------
if __name__ == '__main__':
if os.path.isfile('GSM经纬度.xlsx'):
      GSM_data=xlrd.open_workbook('GSM经纬度.xlsx')
      GSM_table=GSM_data.sheets()[0]
      GSM_Btsname=GSM_table.col_values(0,start_rowx=1,end_rowx=None)
      GSM_lon=GSM_table.col_values(1,start_rowx=1,end_rowx=None)
      GSM_lat=GSM_table.col_values(2,start_rowx=1,end_rowx=None)
if os.path.isfile('LTE经纬度.xlsx'):
      LTE_data=xlrd.open_workbook('LTE经纬度.xlsx')
      LTE_table=LTE_data.sheets()[0]
      LTE_Btsname=LTE_table.col_values(0,start_rowx=1,end_rowx=None)
      LTE_lon=LTE_table.col_values(1,start_rowx=1,end_rowx=None)
      LTE_lat=LTE_table.col_values(2,start_rowx=1,end_rowx=None)

dis = [['网元名称','最小站间距(m)','目标站点']]
print('最小站间距开始计算，请耐心等待:')
GSM_lg = len(GSM_Btsname)

for i in range(0,GSM_lg):
      Bts_dic = 100000000
      tem = [GSM_Btsname]
      name = ''
      LTE_lg = len(LTE_Btsname)
      for j in range(0,LTE_lg):
         lon1=GSM_lon
         lat1=GSM_lat

         lon1, lat1, lon2, lat2 = map(math.radians, [float(lon1), float(lat1), float(LTE_lon[j]), float(LTE_lat[j])]) # 经纬度转换成弧度
         dlon = lon2 - lon1
         dlat = lat2 - lat1
         a = math.sin(dlat / 2) ** 2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon / 2) ** 2
         c = 2 * math.asin(math.sqrt(a))
         r = 6371.137 # 地球平均半径，单位为公里
         dic = (c * r * 1000)
         if dic < Bts_dic:
            Bts_dic = dic
            name = LTE_Btsname[j]
      tem.append(str('%.2f'%Bts_dic))
      tem.append(name)
      dis.append(tem)
      print('+',end='')
      sys.stdout.flush()
write_xlsx(dis)

时间: 2020-2-29 16:20
作者: ran2cn

已阅。

时间: 2020-5-1 19:52
作者: 山东网通

楼主辛苦了，学习了。

时间: 2020-5-12 09:47
作者: 就叫张黑黑吧

谢谢楼主

时间: 2022-2-26 08:33
作者: 浮云游子意

谢谢分享

时间: 2022-2-26 11:24
作者: 涛声依旧79

感谢

通信人家园 (https://www.txrjy.com/) Powered by C114