【GUI界面】基于Python的WSG84三点定位系统(经纬度坐标与平面坐标转换法求解)
文章目录
- 三点定位
- 附录:列表的赋值类型和py打包
- 列表赋值
- BUG复现
- 代码改进
- 优化
- 总结
- py打包
三点定位
方法汇总:
【精准三点定位求解汇总】利用Python或JavaScript高德地图开放平台实现精准三点定位(经纬度坐标与平面坐标转换法求解、几何绘图法求解)
众所周知,如果已知三个点的坐标,到一个未知点的距离,则可以利用以距离为半径画圆的方式来求得未知点坐标。
如果只有两个已知点,则只能得出两个未知点坐标,而第三个圆必定交于其中一个点
如图:
三个圆必定教于一个点
当然,如果第三个绿色圆圆心位于红蓝的圆心连线上,则依然交于两个点,所以在选择对照点时,应尽可能使对照点分布在未知点四周,多取几个点位未尝不是一门好事
主要用到这两个库:
import math
import pyproj
pyproj是用于坐标转换的 这里采用的是utm平面坐标与WGS84经纬度坐标
这里的半径r单位都是米
高德地图用的是GCJ02坐标(还有腾讯) GPS用的是WGS84坐标(谷歌也是) 百度地图用的是BD09坐标
所以在实际计算出来导入地图里面查看时 要么采用WGS84地图 要么就要涉及到坐标转换
首先是坐标转换
在utm坐标中 有个zone值 也就是经度分区 计算公式为int(lon/6+31)
在反坐标转换中 也要输入zone值 这里直接可以输入转换前求得的zone值
def lonlat2utm(lon,lat):
z=int(lon/6+31)
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(lon, lat),z
def utm2lonlat(x,y,z):
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(x, y,inverse=True)
平面上求圆的交点
def insec(p1,r1,p2,r2):
x = p1[0]
y = p1[1]
R = r1
a = p2[0]
b = p2[1]
S = r2
d = math.sqrt((abs(a-x))**2 + (abs(b-y))**2)
if d > (R+S) or d < (abs(R-S)):
# print ("没有公共点")
return
elif d == 0 and R==S :
# print ("两个圆同心")
return
else:
A = (R**2 - S**2 + d**2) / (2 * d)
h = math.sqrt(R**2 - A**2)
x2 = x + A * (a-x)/d
y2 = y + A * (b-y)/d
x3 = x2 - h * (b - y) / d
y3 = y2 + h * (a - x) / d
x4 = x2 + h * (b - y) / d
y4 = y2 - h * (a - x) / d
c1=[x3, y3]
c2=[x4, y4]
return c1,c2
经纬度上求两个圆交点坐标
在反坐标转换中 也要输入zone值 这里直接可以输入转换前求得的两个经纬度坐标的zone平均值
所以 两个要求的区域离得不远 误差就很小 离得太远了 误差就可能大 需要手动去调整 这里是通用函数
def location_trans(p1,r1,p2,r2):
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z=int((z1[1]+z2[1])/2)
C=insec(z1[0],r1,z2[0],r2)
if C:
a=utm2lonlat(C[0][0],C[0][1],z)
b=utm2lonlat(C[1][0],C[1][1],z)
return a,b
else:
return None,None
运行:
a=[[114.304569,30.593354],300000]
b=[[115.857972,28.682976],400000]
c=[[116.378517,39.865246],900000]
print(location_trans(b[0],b[1],c[0],c[1]))
输出:
([114.12482189881902, 31.962881802790577], [117.87031764680636, 31.841927527011755])
在求三个圆的交点时 最多会求出六个点
在六个点中筛选出离另外一个圆最近的点 即可得出三个相近点的坐标
求离得近的那个点的平面坐标
def location_min(p1,p2,p,r):
d1=math.fabs(r-math.sqrt((p[0]-p1[0])**2+(p[1]-p1[1])**2))
d2=math.fabs(r-math.sqrt((p[0]-p2[0])**2+(p[1]-p2[1])**2))
if d1<d2:
return p1
else:
return p2
得到三个点后 求三个点的中心点 即可算出大概位置
中心点的zone值是根据其他三个zone值求平均来确定的
所以 三个要求的区域离得不远 误差就很小 离得太远了 误差就可能大 需要手动去调整 这里是通用函数
def location_judg(p1,r1,p2,r2,p3,r3):
li=[]
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z3=lonlat2utm(p3[0],p3[1])
z12=int((z1[1]+z2[1])/2)
z13=int((z1[1]+z3[1])/2)
z23=int((z2[1]+z3[1])/2)
z=int((z12+z13+z23)/3)
C12=insec(z1[0],r1,z2[0],r2)
C13=insec(z1[0],r1,z3[0],r3)
C23=insec(z2[0],r2,z3[0],r3)
if C12:
m12=location_min(C12[0],C12[1],z3[0],r3)
li.append(utm2lonlat(m12[0],m12[1],z12))
else:
li.append(None)
if C13:
m13=location_min(C13[0],C13[1],z2[0],r2)
li.append(utm2lonlat(m13[0],m13[1],z13))
else:
li.append(None)
if C23:
m23=location_min(C23[0],C23[1],z1[0],r1)
li.append(utm2lonlat(m23[0],m23[1],z23))
else:
li.append(None)
if C12 and C13 and C23:
# print("三个坐标作的圆都有公共点")
m=[(m12[0]+m13[0]+m23[0])/3,(m12[1]+m13[1]+m23[1])/3]
li.append(utm2lonlat(m[0],m[1],z))
return li
elif C12 or C13 or C23:
# print("三个坐标作的圆不全有公共点")
li.append(None)
return li
else:
# print("三个坐标作的圆都没有公共点")
return
最后返回的列表分别是12的最接近坐标 13的最接近坐标 23的最接近坐标 和这三个坐标的中心点坐标 如果不存在 则返回None
运行:
if __name__ == "__main__":
a=[[114.304569,30.593354],300000]
b=[[115.857972,28.682976],400000]
c=[[116.378517,39.865246],900000]
print(location_trans(b[0],b[1],c[0],c[1]))
print(location_judg(a[0],a[1],b[0],b[1],c[0],c[1]))
结果:
[(116.85351953263574, 32.18782636821823), (117.13697531307241, 31.774218803048125), (117.87031764680636, 31.841927527011755), (117.28744847106574, 31.935380071325877)]
整体代码:
# -*- coding: utf-8 -*-
import math
import pyproj
def lonlat2utm(lon,lat):
z=int(lon/6+31)
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(lon, lat),z
def utm2lonlat(x,y,z):
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(x, y,inverse=True)
def insec(p1,r1,p2,r2):
x = p1[0]
y = p1[1]
R = r1
a = p2[0]
b = p2[1]
S = r2
d = math.sqrt((abs(a-x))**2 + (abs(b-y))**2)
if d > (R+S) or d < (abs(R-S)):
# print ("没有公共点")
return
elif d == 0 and R==S :
# print ("两个圆同心")
return
else:
A = (R**2 - S**2 + d**2) / (2 * d)
h = math.sqrt(R**2 - A**2)
x2 = x + A * (a-x)/d
y2 = y + A * (b-y)/d
x3 = x2 - h * (b - y) / d
y3 = y2 + h * (a - x) / d
x4 = x2 + h * (b - y) / d
y4 = y2 - h * (a - x) / d
c1=[x3, y3]
c2=[x4, y4]
return c1,c2
def location_trans(p1,r1,p2,r2):
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z=int((z1[1]+z2[1])/2)
C=insec(z1[0],r1,z2[0],r2)
if C:
a=utm2lonlat(C[0][0],C[0][1],z)
b=utm2lonlat(C[1][0],C[1][1],z)
return a,b
else:
return None,None
def location_min(p1,p2,p,r):
d1=math.fabs(r-math.sqrt((p[0]-p1[0])**2+(p[1]-p1[1])**2))
d2=math.fabs(r-math.sqrt((p[0]-p2[0])**2+(p[1]-p2[1])**2))
if d1<d2:
return p1
else:
return p2
def location_judg(p1,r1,p2,r2,p3,r3):
li=[]
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z3=lonlat2utm(p3[0],p3[1])
z12=int((z1[1]+z2[1])/2)
z13=int((z1[1]+z3[1])/2)
z23=int((z2[1]+z3[1])/2)
z=int((z12+z13+z23)/3)
C12=insec(z1[0],r1,z2[0],r2)
C13=insec(z1[0],r1,z3[0],r3)
C23=insec(z2[0],r2,z3[0],r3)
if C12:
m12=location_min(C12[0],C12[1],z3[0],r3)
li.append(utm2lonlat(m12[0],m12[1],z12))
else:
li.append(None)
if C13:
m13=location_min(C13[0],C13[1],z2[0],r2)
li.append(utm2lonlat(m13[0],m13[1],z13))
else:
li.append(None)
if C23:
m23=location_min(C23[0],C23[1],z1[0],r1)
li.append(utm2lonlat(m23[0],m23[1],z23))
else:
li.append(None)
if C12 and C13 and C23:
# print("三个坐标作的圆都有公共点")
m=[(m12[0]+m13[0]+m23[0])/3,(m12[1]+m13[1]+m23[1])/3]
li.append(utm2lonlat(m[0],m[1],z))
return li
elif C12 or C13 or C23:
# print("三个坐标作的圆不全有公共点")
li.append(None)
return li
else:
# print("三个坐标作的圆都没有公共点")
return
if __name__ == "__main__":
a=[[114.304569,30.593354],300000]
b=[[115.857972,28.682976],400000]
c=[[116.378517,39.865246],900000]
print(location_trans(b[0],b[1],c[0],c[1]))
print(location_judg(a[0],a[1],b[0],b[1],c[0],c[1]))
加上界面:
def central_win(win):
win.resizable(0,0) # 不可缩放
screenwidth = win.winfo_screenwidth() # 获取屏幕分辨率宽
screenheight = win.winfo_screenheight() # 获取屏幕分辨率高
win.update() # 更新窗口
width = win.winfo_width() # 重新赋值
height = win.winfo_height()
size = '+%d+%d' % ((screenwidth - width)/2, (screenheight - height)/2)
# 重新赋值大小 大小为屏幕大小/2
win.geometry(size) # 以新大小定义窗口
def gui_start():
root=tk.Tk()
root.title("WSG84三点定位系统 By 网易独家音乐人Mike Zhou")
mainfram=tk.Frame(root,width=500, height=920)
mainfram.grid_propagate(0)
mainfram.grid()
central_win(root)
labelName=tk.Label(root, text='经度(正E负W)', justify=tk.LEFT)
labelName.place(x=120, y=20, width=180, height=50)
labelName=tk.Label(root, text='纬度(正N负S)', justify=tk.LEFT)
labelName.place(x=300, y=20, width=180, height=50)
labelName=tk.Label(root, text='位置1:', justify=tk.LEFT)
labelName.place(x=20, y=70, width=100, height=50)
labelName=tk.Label(root, text='距离1:', justify=tk.LEFT)
labelName.place(x=20, y=120, width=100, height=50)
labelName=tk.Label(root, text='位置2:', justify=tk.LEFT)
labelName.place(x=20, y=170, width=100, height=50)
labelName=tk.Label(root, text='距离2:', justify=tk.LEFT)
labelName.place(x=20, y=220, width=100, height=50)
labelName=tk.Label(root, text='位置3:', justify=tk.LEFT)
labelName.place(x=20, y=270, width=100, height=50)
labelName=tk.Label(root, text='距离3:', justify=tk.LEFT)
labelName.place(x=20, y=320, width=100, height=50)
labelName=tk.Label(root, text='12推测点:', justify=tk.LEFT)
labelName.place(x=20, y=440, width=100, height=50)
labelName=tk.Label(root, text='13推测点:', justify=tk.LEFT)
labelName.place(x=20, y=490, width=100, height=50)
labelName=tk.Label(root, text='23推测点:', justify=tk.LEFT)
labelName.place(x=20, y=540, width=100, height=50)
labelName=tk.Label(root, text='中心推测点:', justify=tk.LEFT)
labelName.place(x=20, y=590, width=100, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=120, width=20, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=220, width=20, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=320, width=20, height=50)
labelName=tk.Label(root, text='两点定位1:', justify=tk.LEFT)
labelName.place(x=20, y=730, width=100, height=50)
labelName=tk.Label(root, text='两点定位2:', justify=tk.LEFT)
labelName.place(x=20, y=780, width=100, height=50)
e1_lon = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e1_lon.delete(0, tk.END) # 将输入框里面的内容清空
e1_lon.insert(0, 114.304569)
e1_lon.place(x=120, y=70, width=180, height=50)
e1_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e1_r.delete(0, tk.END) # 将输入框里面的内容清空
e1_r.insert(0, 300000)
e1_r.place(x=120, y=120, width=180, height=50)
e2_lon = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e2_lon.delete(0, tk.END) # 将输入框里面的内容清空
e2_lon.insert(0, 115.857972)
e2_lon.place(x=120, y=170, width=180, height=50)
e2_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e2_r.delete(0, tk.END) # 将输入框里面的内容清空
e2_r.insert(0, 400000)
e2_r.place(x=120, y=220, width=180, height=50)
e3_lon = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e3_lon.delete(0, tk.END) # 将输入框里面的内容清空
e3_lon.insert(0, 116.378517)
e3_lon.place(x=120, y=270, width=180, height=50)
e3_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e3_r.delete(0, tk.END) # 将输入框里面的内容清空
e3_r.insert(0, 900000)
e3_r.place(x=120, y=320, width=180, height=50)
e1_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e1_lat.delete(0, tk.END) # 将输入框里面的内容清空
e1_lat.insert(0, 30.593354)
e1_lat.place(x=300, y=70, width=180, height=50)
e2_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e2_lat.delete(0, tk.END) # 将输入框里面的内容清空
e2_lat.insert(0, 28.682976)
e2_lat.place(x=300, y=170, width=180, height=50)
e3_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e3_lat.delete(0, tk.END) # 将输入框里面的内容清空
e3_lat.insert(0, 39.865246)
e3_lat.place(x=300, y=270, width=180, height=50)
ex1 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex1.delete(0, tk.END) # 将输入框里面的内容清空
ex1.insert(0, "")
ex1.place(x=120, y=440, width=360, height=50)
ex1.config(state='readonly')
ex2 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex2.delete(0, tk.END) # 将输入框里面的内容清空
ex2.insert(0, "")
ex2.place(x=120, y=490, width=360, height=50)
ex2.config(state='readonly')
ex3 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex3.delete(0, tk.END) # 将输入框里面的内容清空
ex3.insert(0, "")
ex3.place(x=120, y=540, width=360, height=50)
ex3.config(state='readonly')
ex4 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex4.delete(0, tk.END) # 将输入框里面的内容清空
ex4.insert(0, "")
ex4.place(x=120, y=590, width=360, height=50)
ex4.config(state='readonly')
ey1 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ey1.delete(0, tk.END) # 将输入框里面的内容清空
ey1.insert(0, "")
ey1.place(x=120, y=730, width=360, height=50)
ey1.config(state='readonly')
ey2 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ey2.delete(0, tk.END) # 将输入框里面的内容清空
ey2.insert(0, "")
ey2.place(x=120, y=780, width=360, height=50)
ey2.config(state='readonly')
def input_judg(e):
try:
s=float(str(e.get()))
except:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
return s
if e==e1_lat or e==e2_lat or e==e3_lat:
if s<-90 or s>90:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
elif e==e1_lon or e==e2_lon or e==e3_lon:
if s<-180 or s>180:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
elif e==e1_r or e==e2_r or e==e3_r:
if s<0 or s>400750170:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
return s
def view(e,s):
s=str(s)
e.config(state='normal')
e.delete(0, tk.END)
e.insert(0, s)
e.config(state='readonly')
def event():
lat1 = input_judg(e1_lat)
lon1 = input_judg(e1_lon)
r1 = input_judg(e1_r)
lat2 = input_judg(e2_lat)
lon2 = input_judg(e2_lon)
r2 = input_judg(e2_r)
lat3 = input_judg(e3_lat)
lon3 = input_judg(e3_lon)
r3 = input_judg(e3_r)
result = location_judg([lon1,lat1],r1,[lon2,lat2],r2,[lon3,lat3],r3)
view(ex1,result[0])
view(ex2,result[1])
view(ex3,result[2])
view(ex4,result[3])
def double(n):
lat1 = input_judg(e1_lat)
lon1 = input_judg(e1_lon)
r1 = input_judg(e1_r)
lat2 = input_judg(e2_lat)
lon2 = input_judg(e2_lon)
r2 = input_judg(e2_r)
lat3 = input_judg(e3_lat)
lon3 = input_judg(e3_lon)
r3 = input_judg(e3_r)
if n==1:
result = location_trans([lon1,lat1],r1,[lon2,lat2],r2)
elif n==2:
result = location_trans([lon1,lat1],r1,[lon3,lat3],r3)
elif n==3:
result = location_trans([lon2,lat2],r2,[lon3,lat3],r3)
else:
result = location_trans([lon1,lat1],r1,[lon2,lat2],r2)
view(ey1,result[0])
view(ey2,result[1])
def clean():
view(ex1,"")
view(ex2,"")
view(ex3,"")
view(ex4,"")
view(ey1,"")
view(ey2,"")
e1_lon.delete(0, tk.END)
e1_lon.insert(0, "")
e2_lon.delete(0, tk.END)
e2_lon.insert(0, "")
e3_lon.delete(0, tk.END)
e3_lon.insert(0, "")
e1_lat.delete(0, tk.END)
e1_lat.insert(0, "")
e2_lat.delete(0, tk.END)
e2_lat.insert(0, "")
e3_lat.delete(0, tk.END)
e3_lat.insert(0, "")
e1_r.delete(0, tk.END)
e1_r.insert(0, "")
e2_r.delete(0, tk.END)
e2_r.insert(0, "")
e3_r.delete(0, tk.END)
e3_r.insert(0, "")
def start_event():
t=Thread(target=event)
t.setDaemon(False)
t.start()
def start_clean():
t=Thread(target=clean)
t.setDaemon(False)
t.start()
def start_double(i):
t=Thread(target=double,args=(i,))
t.setDaemon(False)
t.start()
def button_event(s):
if(s.keysym=='C' or s.keysym=='c'):
start_event()
elif(s.keysym=='D' or s.keysym=='d'):
start_clean()
b1=tk.Button(mainfram,width=30,text="计算123三点定位",command=start_event)
b1.bind_all('C',button_event)
b1.bind_all('c',button_event)
b1.place(x=20, y=640, width=460, height=70)
b2=tk.Button(mainfram,width=30,text="清空",command=start_clean)
b2.bind_all('D',button_event)
b2.bind_all('d',button_event)
b2.place(x=20, y=370, width=460, height=50)
b3=tk.Button(mainfram,width=30,text="计算12两点定位",command=lambda:start_double(1))
b3.place(x=20, y=830, width=140, height=70)
b4=tk.Button(mainfram,width=30,text="计算13两点定位",command=lambda:start_double(2))
b4.place(x=180, y=830, width=140, height=70)
b5=tk.Button(mainfram,width=30,text="计算23两点定位",command=lambda:start_double(3))
b5.place(x=340, y=830, width=140, height=70)
start_event()
start_double(1)
root.mainloop()
最后运行:
所有代码:
# -*- coding: utf-8 -*-
import math
import pyproj
import tkinter as tk
from threading import Thread
def lonlat2utm(lon,lat):
z=int(lon/6+31)
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(lon, lat),z
def utm2lonlat(x,y,z):
proj = pyproj.Proj(proj='utm',zone=z,ellps='WGS84')
return proj(x, y,inverse=True)
def insec(p1,r1,p2,r2):
x = p1[0]
y = p1[1]
R = r1
a = p2[0]
b = p2[1]
S = r2
d = math.sqrt((abs(a-x))**2 + (abs(b-y))**2)
if d > (R+S) or d < (abs(R-S)):
# print ("没有公共点")
return
elif d == 0 and R==S :
# print ("两个圆同心")
return
else:
A = (R**2 - S**2 + d**2) / (2 * d)
h = math.sqrt(R**2 - A**2)
x2 = x + A * (a-x)/d
y2 = y + A * (b-y)/d
x3 = x2 - h * (b - y) / d
y3 = y2 + h * (a - x) / d
x4 = x2 + h * (b - y) / d
y4 = y2 - h * (a - x) / d
c1=[x3, y3]
c2=[x4, y4]
return c1,c2
def location_trans(p1,r1,p2,r2):
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z=int((z1[1]+z2[1])/2)
C=insec(z1[0],r1,z2[0],r2)
if C:
a=utm2lonlat(C[0][0],C[0][1],z)
b=utm2lonlat(C[1][0],C[1][1],z)
return a,b
else:
return None,None
def location_min(p1,p2,p,r):
d1=math.fabs(r-math.sqrt((p[0]-p1[0])**2+(p[1]-p1[1])**2))
d2=math.fabs(r-math.sqrt((p[0]-p2[0])**2+(p[1]-p2[1])**2))
if d1<d2:
return p1
else:
return p2
def location_judg(p1,r1,p2,r2,p3,r3):
li=[]
z1=lonlat2utm(p1[0],p1[1])
z2=lonlat2utm(p2[0],p2[1])
z3=lonlat2utm(p3[0],p3[1])
z12=int((z1[1]+z2[1])/2)
z13=int((z1[1]+z3[1])/2)
z23=int((z2[1]+z3[1])/2)
z=int((z12+z13+z23)/3)
C12=insec(z1[0],r1,z2[0],r2)
C13=insec(z1[0],r1,z3[0],r3)
C23=insec(z2[0],r2,z3[0],r3)
if C12:
m12=location_min(C12[0],C12[1],z3[0],r3)
li.append(utm2lonlat(m12[0],m12[1],z12))
else:
li.append(None)
if C13:
m13=location_min(C13[0],C13[1],z2[0],r2)
li.append(utm2lonlat(m13[0],m13[1],z13))
else:
li.append(None)
if C23:
m23=location_min(C23[0],C23[1],z1[0],r1)
li.append(utm2lonlat(m23[0],m23[1],z23))
else:
li.append(None)
if C12 and C13 and C23:
# print("三个坐标作的圆都有公共点")
m=[(m12[0]+m13[0]+m23[0])/3,(m12[1]+m13[1]+m23[1])/3]
li.append(utm2lonlat(m[0],m[1],z))
return li
elif C12 or C13 or C23:
# print("三个坐标作的圆不全有公共点")
li.append(None)
return li
else:
# print("三个坐标作的圆都没有公共点")
li.append(None)
li.append(None)
li.append(None)
li.append(None)
return li
def central_win(win):
win.resizable(0,0) # 不可缩放
screenwidth = win.winfo_screenwidth() # 获取屏幕分辨率宽
screenheight = win.winfo_screenheight() # 获取屏幕分辨率高
win.update() # 更新窗口
width = win.winfo_width() # 重新赋值
height = win.winfo_height()
size = '+%d+%d' % ((screenwidth - width)/2, (screenheight - height)/2)
# 重新赋值大小 大小为屏幕大小/2
win.geometry(size) # 以新大小定义窗口
def gui_start():
root=tk.Tk()
root.title("WSG84三点定位系统 By 网易独家音乐人Mike Zhou")
mainfram=tk.Frame(root,width=500, height=920)
mainfram.grid_propagate(0)
mainfram.grid()
central_win(root)
labelName=tk.Label(root, text='经度(正E负W)', justify=tk.LEFT)
labelName.place(x=120, y=20, width=180, height=50)
labelName=tk.Label(root, text='纬度(正N负S)', justify=tk.LEFT)
labelName.place(x=300, y=20, width=180, height=50)
labelName=tk.Label(root, text='位置1:', justify=tk.LEFT)
labelName.place(x=20, y=70, width=100, height=50)
labelName=tk.Label(root, text='距离1:', justify=tk.LEFT)
labelName.place(x=20, y=120, width=100, height=50)
labelName=tk.Label(root, text='位置2:', justify=tk.LEFT)
labelName.place(x=20, y=170, width=100, height=50)
labelName=tk.Label(root, text='距离2:', justify=tk.LEFT)
labelName.place(x=20, y=220, width=100, height=50)
labelName=tk.Label(root, text='位置3:', justify=tk.LEFT)
labelName.place(x=20, y=270, width=100, height=50)
labelName=tk.Label(root, text='距离3:', justify=tk.LEFT)
labelName.place(x=20, y=320, width=100, height=50)
labelName=tk.Label(root, text='12推测点:', justify=tk.LEFT)
labelName.place(x=20, y=440, width=100, height=50)
labelName=tk.Label(root, text='13推测点:', justify=tk.LEFT)
labelName.place(x=20, y=490, width=100, height=50)
labelName=tk.Label(root, text='23推测点:', justify=tk.LEFT)
labelName.place(x=20, y=540, width=100, height=50)
labelName=tk.Label(root, text='中心推测点:', justify=tk.LEFT)
labelName.place(x=20, y=590, width=100, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=120, width=20, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=220, width=20, height=50)
labelName=tk.Label(root, text='米', justify=tk.LEFT)
labelName.place(x=300, y=320, width=20, height=50)
labelName=tk.Label(root, text='两点定位1:', justify=tk.LEFT)
labelName.place(x=20, y=730, width=100, height=50)
labelName=tk.Label(root, text='两点定位2:', justify=tk.LEFT)
labelName.place(x=20, y=780, width=100, height=50)
e1_lon = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e1_lon.delete(0, tk.END) # 将输入框里面的内容清空
e1_lon.insert(0, 114.304569)
e1_lon.place(x=120, y=70, width=180, height=50)
e1_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e1_r.delete(0, tk.END) # 将输入框里面的内容清空
e1_r.insert(0, 300000)
e1_r.place(x=120, y=120, width=180, height=50)
e2_lon = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e2_lon.delete(0, tk.END) # 将输入框里面的内容清空
e2_lon.insert(0, 115.857972)
e2_lon.place(x=120, y=170, width=180, height=50)
e2_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e2_r.delete(0, tk.END) # 将输入框里面的内容清空
e2_r.insert(0, 400000)
e2_r.place(x=120, y=220, width=180, height=50)
e3_lon = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e3_lon.delete(0, tk.END) # 将输入框里面的内容清空
e3_lon.insert(0, 116.378517)
e3_lon.place(x=120, y=270, width=180, height=50)
e3_r = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
e3_r.delete(0, tk.END) # 将输入框里面的内容清空
e3_r.insert(0, 900000)
e3_r.place(x=120, y=320, width=180, height=50)
e1_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e1_lat.delete(0, tk.END) # 将输入框里面的内容清空
e1_lat.insert(0, 30.593354)
e1_lat.place(x=300, y=70, width=180, height=50)
e2_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e2_lat.delete(0, tk.END) # 将输入框里面的内容清空
e2_lat.insert(0, 28.682976)
e2_lat.place(x=300, y=170, width=180, height=50)
e3_lat = tk.Entry(mainfram)
#e1.grid(ipadx=50,row=0,column=1)
e3_lat.delete(0, tk.END) # 将输入框里面的内容清空
e3_lat.insert(0, 39.865246)
e3_lat.place(x=300, y=270, width=180, height=50)
ex1 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex1.delete(0, tk.END) # 将输入框里面的内容清空
ex1.insert(0, "")
ex1.place(x=120, y=440, width=360, height=50)
ex1.config(state='readonly')
ex2 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex2.delete(0, tk.END) # 将输入框里面的内容清空
ex2.insert(0, "")
ex2.place(x=120, y=490, width=360, height=50)
ex2.config(state='readonly')
ex3 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex3.delete(0, tk.END) # 将输入框里面的内容清空
ex3.insert(0, "")
ex3.place(x=120, y=540, width=360, height=50)
ex3.config(state='readonly')
ex4 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ex4.delete(0, tk.END) # 将输入框里面的内容清空
ex4.insert(0, "")
ex4.place(x=120, y=590, width=360, height=50)
ex4.config(state='readonly')
ey1 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ey1.delete(0, tk.END) # 将输入框里面的内容清空
ey1.insert(0, "")
ey1.place(x=120, y=730, width=360, height=50)
ey1.config(state='readonly')
ey2 = tk.Entry(mainfram)
#e2.grid(ipadx=50,row=1,column=1)
ey2.delete(0, tk.END) # 将输入框里面的内容清空
ey2.insert(0, "")
ey2.place(x=120, y=780, width=360, height=50)
ey2.config(state='readonly')
def input_judg(e):
try:
s=float(str(e.get()))
except:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
return s
if e==e1_lat or e==e2_lat or e==e3_lat:
if s<-90 or s>90:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
elif e==e1_lon or e==e2_lon or e==e3_lon:
if s<-180 or s>180:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
elif e==e1_r or e==e2_r or e==e3_r:
if s<0 or s>400750170:
e.delete(0, tk.END)
e.insert(0, 0)
s=0.0
return s
def view(e,s):
s=str(s)
e.config(state='normal')
e.delete(0, tk.END)
e.insert(0, s)
e.config(state='readonly')
def event():
lat1 = input_judg(e1_lat)
lon1 = input_judg(e1_lon)
r1 = input_judg(e1_r)
lat2 = input_judg(e2_lat)
lon2 = input_judg(e2_lon)
r2 = input_judg(e2_r)
lat3 = input_judg(e3_lat)
lon3 = input_judg(e3_lon)
r3 = input_judg(e3_r)
result = location_judg([lon1,lat1],r1,[lon2,lat2],r2,[lon3,lat3],r3)
view(ex1,result[0])
view(ex2,result[1])
view(ex3,result[2])
view(ex4,result[3])
def double(n):
lat1 = input_judg(e1_lat)
lon1 = input_judg(e1_lon)
r1 = input_judg(e1_r)
lat2 = input_judg(e2_lat)
lon2 = input_judg(e2_lon)
r2 = input_judg(e2_r)
lat3 = input_judg(e3_lat)
lon3 = input_judg(e3_lon)
r3 = input_judg(e3_r)
if n==1:
result = location_trans([lon1,lat1],r1,[lon2,lat2],r2)
elif n==2:
result = location_trans([lon1,lat1],r1,[lon3,lat3],r3)
elif n==3:
result = location_trans([lon2,lat2],r2,[lon3,lat3],r3)
else:
result = location_trans([lon1,lat1],r1,[lon2,lat2],r2)
view(ey1,result[0])
view(ey2,result[1])
def clean():
view(ex1,"")
view(ex2,"")
view(ex3,"")
view(ex4,"")
view(ey1,"")
view(ey2,"")
e1_lon.delete(0, tk.END)
e1_lon.insert(0, "")
e2_lon.delete(0, tk.END)
e2_lon.insert(0, "")
e3_lon.delete(0, tk.END)
e3_lon.insert(0, "")
e1_lat.delete(0, tk.END)
e1_lat.insert(0, "")
e2_lat.delete(0, tk.END)
e2_lat.insert(0, "")
e3_lat.delete(0, tk.END)
e3_lat.insert(0, "")
e1_r.delete(0, tk.END)
e1_r.insert(0, "")
e2_r.delete(0, tk.END)
e2_r.insert(0, "")
e3_r.delete(0, tk.END)
e3_r.insert(0, "")
def start_event():
t=Thread(target=event)
t.setDaemon(False)
t.start()
def start_clean():
t=Thread(target=clean)
t.setDaemon(False)
t.start()
def start_double(i):
t=Thread(target=double,args=(i,))
t.setDaemon(False)
t.start()
def button_event(s):
if(s.keysym=='C' or s.keysym=='c'):
start_event()
elif(s.keysym=='D' or s.keysym=='d'):
start_clean()
b1=tk.Button(mainfram,width=30,text="计算123三点定位",command=start_event)
b1.bind_all('C',button_event)
b1.bind_all('c',button_event)
b1.place(x=20, y=640, width=460, height=70)
b2=tk.Button(mainfram,width=30,text="清空",command=start_clean)
b2.bind_all('D',button_event)
b2.bind_all('d',button_event)
b2.place(x=20, y=370, width=460, height=50)
b3=tk.Button(mainfram,width=30,text="计算12两点定位",command=lambda:start_double(1))
b3.place(x=20, y=830, width=140, height=70)
b4=tk.Button(mainfram,width=30,text="计算13两点定位",command=lambda:start_double(2))
b4.place(x=180, y=830, width=140, height=70)
b5=tk.Button(mainfram,width=30,text="计算23两点定位",command=lambda:start_double(3))
b5.place(x=340, y=830, width=140, height=70)
start_event()
start_double(1)
root.mainloop()
if __name__ == "__main__":
gui_start()
附录:列表的赋值类型和py打包
列表赋值
BUG复现
闲来无事写了个小程序 代码如下:
# -*- coding: utf-8 -*-
"""
Created on Fri Nov 19 19:47:01 2021
@author: 16016
"""
a_list = ['0','1','2','3','4','5','6','7','8','9','10','11','12','13','14','15']
#print(len(a_list))
#b_list = ['','','','','','','','','','','','','','','','']
c_list = [[],[],[],[],[],[],[],[],[],[],[],[],[],[],[],[]]
#for i in range(16):
if len(a_list):
for j in range(16):
a_list[j]=str(a_list[j])+'_'+str(j)
print("序号:",j)
print('a_list:\n',a_list)
c_list[j]=a_list
print('c_list[0]:\n',c_list[0])
print('\n')
# b_list[j]=a_list[7],a_list[8]
# print(b_list[j])
# 写入到Excel:
#print(c_list,'\n')
我在程序中 做了一个16次的for循环 把列表a的每个值后面依次加上"_"和循环序号
比如循环第x次 就是把第x位加上_x 这一位变成x_x 我在输出测试中 列表a的每一次输出也是对的
循环16次后列表a应该变成[‘0_0’, ‘1_1’, ‘2_2’, ‘3_3’, ‘4_4’, ‘5_5’, ‘6_6’, ‘7_7’, ‘8_8’, ‘9_9’, ‘10_10’, ‘11_11’, ‘12_12’, ‘13_13’, ‘14_14’, ‘15_15’] 这也是对的
同时 我将每一次循环时列表a的值 写入到空列表c中 比如第x次循环 就是把更改以后的列表a的值 写入到列表c的第x位
第0次循环后 c[0]的值应该是[‘0_0’, ‘1’, ‘2’, ‘3’, ‘4’, ‘5’, ‘6’, ‘7’, ‘8’, ‘9’, ‘10’, ‘11’, ‘12’, ‘13’, ‘14’, ‘15’] 这也是对的
但是在第1次循环以后 c[0]的值就一直在变 变成了c[x]的值
相当于把c_list[0]变成了c_list[1]…以此类推 最后得出的列表c的值也是每一项完全一样
我不明白这是怎么回事
我的c[0]只在第0次循环时被赋值了 但是后面它的值跟着在改变
如图:
第一次老出bug 赋值以后 每次循环都改变c[0]的值 搞了半天都没搞出来
无论是用appen函数添加 还是用二维数组定义 或者增加第三个空数组来过渡 都无法解决
代码改进
后来在我华科同学的指导下 突然想到赋值可以赋的是个地址 地址里面的值一直变化 导致赋值也一直变化 于是用第二张图的循环套循环深度复制实现了
代码如下:
# -*- coding: utf-8 -*-
"""
Created on Fri Nov 19 19:47:01 2021
@author: 16016
"""
a_list = ['0','1','2','3','4','5','6','7','8','9','10','11','12','13','14','15']
#print(len(a_list))
#b_list = ['','','','','','','','','','','','','','','','']
c_list = [[],[],[],[],[],[],[],[],[],[],[],[],[],[],[],[]]
#for i in range(16):
if len(a_list):
for j in range(16):
a_list[j]=str(a_list[j])+'_'+str(j)
print("序号:",j)
print('a_list:\n',a_list)
for i in range(16):
c_list[j].append(a_list[i])
print('c_list[0]:\n',c_list[0])
print('\n')
# b_list[j]=a_list[7],a_list[8]
# print(b_list[j])
# 写入到Excel:
print(c_list,'\n')
解决了问题
优化
第三次是请教了老师 用copy函数来赋真值
代码如下:
# -*- coding: utf-8 -*-
"""
Created on Fri Nov 19 19:47:01 2021
@author: 16016
"""
a_list = ['0','1','2','3','4','5','6','7','8','9','10','11','12','13','14','15']
#print(len(a_list))
#b_list = ['','','','','','','','','','','','','','','','']
c_list = [[],[],[],[],[],[],[],[],[],[],[],[],[],[],[],[]]
#for i in range(16):
if len(a_list):
for j in range(16):
a_list[j]=str(a_list[j])+'_'+str(j)
print("序号:",j)
print('a_list:\n',a_list)
c_list[j]=a_list.copy()
print('c_list[0]:\n',c_list[0])
print('\n')
# b_list[j]=a_list[7],a_list[8]
# print(b_list[j])
# 写入到Excel:
#print(c_list,'\n')
同样能解决问题
最后得出问题 就是指针惹的祸!
a_list指向的是个地址 而不是值 a_list[i]指向的才是单个的值 copy()函数也是复制值而不是地址
如果这个用C语言来写 就直观一些了 难怪C语言是基础 光学Python不学C 遇到这样的问题就解决不了
C语言yyds Python是什么垃圾弱智语言
总结
由于Python无法单独定义一个值为指针或者独立的值 所以只能用列表来传送
只要赋值是指向一个列表整体的 那么就是指向的一个指针内存地址 解决方法只有一个 那就是将每个值深度复制赋值(子列表内的元素提取出来重新依次连接) 或者用copy函数单独赋值
如图测试:
部分代码:
# -*- coding: utf-8 -*-
"""
Created on Sat Nov 20 16:45:48 2021
@author: 16016
"""
def text1():
A=[1,2,3]
B=[[],[],[]]
for i in range(len(A)):
A[i]=A[i]+i
B[i]=A
print(B)
def text2():
A=[1,2,3]
B=[[],[],[]]
A[0]=A[0]+0
B[0]=A
print(B)
A[1]=A[1]+1
B[1]=A
print(B)
A[2]=A[2]+2
B[2]=A
print(B)
if __name__ == '__main__':
text1()
print('\n')
text2()
py打包
Pyinstaller打包exe(包括打包资源文件 绝不出错版)
依赖包及其对应的版本号
PyQt5 5.10.1
PyQt5-Qt5 5.15.2
PyQt5-sip 12.9.0
pyinstaller 4.5.1
pyinstaller-hooks-contrib 2021.3
Pyinstaller -F setup.py 打包exe
Pyinstaller -F -w setup.py 不带控制台的打包
Pyinstaller -F -i xx.ico setup.py 打包指定exe图标打包
打包exe参数说明:
-F:打包后只生成单个exe格式文件;
-D:默认选项,创建一个目录,包含exe文件以及大量依赖文件;
-c:默认选项,使用控制台(就是类似cmd的黑框);
-w:不使用控制台;
-p:添加搜索路径,让其找到对应的库;
-i:改变生成程序的icon图标。
如果要打包资源文件
则需要对代码中的路径进行转换处理
另外要注意的是 如果要打包资源文件 则py程序里面的路径要从./xxx/yy换成xxx/yy 并且进行路径转换
但如果不打包资源文件的话 最好路径还是用作./xxx/yy 并且不进行路径转换
def get_resource_path(relative_path):
if hasattr(sys, '_MEIPASS'):
return os.path.join(sys._MEIPASS, relative_path)
return os.path.join(os.path.abspath("."), relative_path)
而后再spec文件中的datas部分加入目录
如:
a = Analysis(['cxk.py'],
pathex=['D:\\Python Test\\cxk'],
binaries=[],
datas=[('root','root')],
hiddenimports=[],
hookspath=[],
hooksconfig={},
runtime_hooks=[],
excludes=[],
win_no_prefer_redirects=False,
win_private_assemblies=False,
cipher=block_cipher,
noarchive=False)
而后直接Pyinstaller -F setup.spec即可
如果打包的文件过大则更改spec文件中的excludes 把不需要的库写进去(但是已经在环境中安装了的)就行
这些不要了的库在上一次编译时的shell里面输出
比如:
然后用pyinstaller --clean -F 某某.spec