def p(longitude,latitude,angle,distance):
ang=angle/180*PI #角BAC
lon0=longitude/180*PI
lat0=latitude/180*PI
nAB=distance/R #AB对应圆心角
ab=R*m.sin(nAB/2)*2
ac=ab*m.cos(ang)
bc=ab*m.sin(ang)
nAC=m.asin(ac/R)*2 #圆心角
nBC=m.asin(bc/R)*2
lat1=lat0+nAC
r=R*m.cos(lat1) #纬度圈半径
nbc=m.asin(bc/(2*r))*2 #纬度圈经度偏移量
longitude1=longitude+nbc*180/PI
latitude1=lat1*180/PI
#------角度修正
if(angle>=0 and angle<=90):
latitude=latitude
elif(angle>=90 and angle<=180):
latitude1-=2*(latitude1-latitude)
elif(angle>=180 and angle<=270):
longitude1-=2*(longitude1-longitude)
latitude1-=2*(lat