【GUI界面】基于Python的WSG84三点定位系统(经纬度坐标与平面坐标转换法求解)

news2025/2/26 19:36:32

【GUI界面】基于Python的WSG84三点定位系统(经纬度坐标与平面坐标转换法求解)

方法汇总:
blog.csdn.net/weixin_53403301/article/details/128441789
【精准三点定位求解汇总】利用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 = round(x2 - h * (b - y) / d,2)
        y3 = round(y2 + h * (a - x) / d,2)
        x4 = round(x2 + h * (b - y) / d,2)
        y4 = round(y2 - h * (a - x) / d,2)
        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)
        c1=[a[0], a[1]]
        c2=[b[0], b[1]]
        return c1,c2
    else:
        return

运行:

 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 = round(x2 - h * (b - y) / d,2)
        y3 = round(y2 + h * (a - x) / d,2)
        x4 = round(x2 + h * (b - y) / d,2)
        y4 = round(y2 - h * (a - x) / d,2)
        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)
        c1=[a[0], a[1]]
        c2=[b[0], b[1]]
        return c1,c2
    else:
        return

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=750)
    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)
    
    e1_lat = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e1_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e1_lat.insert(0, 114.304569)  
    e1_lat.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_lat = tk.Entry(mainfram)
    #e2.grid(ipadx=50,row=1,column=1)
    e2_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e2_lat.insert(0, 115.857972)
    e2_lat.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_lat = tk.Entry(mainfram)
    #e2.grid(ipadx=50,row=1,column=1)
    e3_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e3_lat.insert(0, 116.378517)
    e3_lat.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_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e1_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e1_lon.insert(0, 30.593354)  
    e1_lon.place(x=300, y=70, width=180, height=50)
    
    e2_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e2_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e2_lon.insert(0, 28.682976)  
    e2_lon.place(x=300, y=170, width=180, height=50)
    
    e3_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e3_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e3_lon.insert(0, 39.865246)  
    e3_lon.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')
    
    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_lon or e==e1_lon or e==e1_lon:
            if s<-90 or s>90:
                e.delete(0, tk.END)
                e.insert(0, 0)
                s=0.0
        elif e==e1_lat or e==e1_lat or e==e1_lat:
            if s<-180 or s>180:
                e.delete(0, tk.END)
                e.insert(0, 0)
                s=0.0
        elif e==e1_r or e==e1_r or e==e1_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([lat1,lon1],r1,[lat2,lon2],r2,[lat3,lon3],r3)
        
        view(ex1,result[0])
        view(ex2,result[1])
        view(ex3,result[2])
        view(ex4,result[3])
    
    def clean():
        view(ex1,"")
        view(ex2,"")
        view(ex3,"")
        view(ex4,"")
        
        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():
        Thread(target=event).setDaemon(False)
        Thread(target=event).start()
    
    def start_clean():
        Thread(target=clean).setDaemon(False)
        Thread(target=clean).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="计算三点定位",command=start_event)
    b1.bind_all('C',button_event)
    b1.bind_all('c',button_event)
    b1.place(x=20, y=660, 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)
    
    start_event()
    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 = round(x2 - h * (b - y) / d,2)
        y3 = round(y2 + h * (a - x) / d,2)
        x4 = round(x2 + h * (b - y) / d,2)
        y4 = round(y2 - h * (a - x) / d,2)
        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)
        c1=[a[0], a[1]]
        c2=[b[0], b[1]]
        return c1,c2
    else:
        return

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=750)
    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)
    
    e1_lat = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e1_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e1_lat.insert(0, 114.304569)  
    e1_lat.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_lat = tk.Entry(mainfram)
    #e2.grid(ipadx=50,row=1,column=1)
    e2_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e2_lat.insert(0, 115.857972)
    e2_lat.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_lat = tk.Entry(mainfram)
    #e2.grid(ipadx=50,row=1,column=1)
    e3_lat.delete(0, tk.END)  # 将输入框里面的内容清空
    e3_lat.insert(0, 116.378517)
    e3_lat.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_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e1_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e1_lon.insert(0, 30.593354)  
    e1_lon.place(x=300, y=70, width=180, height=50)
    
    e2_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e2_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e2_lon.insert(0, 28.682976)  
    e2_lon.place(x=300, y=170, width=180, height=50)
    
    e3_lon = tk.Entry(mainfram)
    #e1.grid(ipadx=50,row=0,column=1)
    e3_lon.delete(0, tk.END)  # 将输入框里面的内容清空
    e3_lon.insert(0, 39.865246)  
    e3_lon.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')
    
    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_lon or e==e1_lon or e==e1_lon:
            if s<-90 or s>90:
                e.delete(0, tk.END)
                e.insert(0, 0)
                s=0.0
        elif e==e1_lat or e==e1_lat or e==e1_lat:
            if s<-180 or s>180:
                e.delete(0, tk.END)
                e.insert(0, 0)
                s=0.0
        elif e==e1_r or e==e1_r or e==e1_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([lat1,lon1],r1,[lat2,lon2],r2,[lat3,lon3],r3)
        
        view(ex1,result[0])
        view(ex2,result[1])
        view(ex3,result[2])
        view(ex4,result[3])
    
    def clean():
        view(ex1,"")
        view(ex2,"")
        view(ex3,"")
        view(ex4,"")
        
        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():
        Thread(target=event).setDaemon(False)
        Thread(target=event).start()
    
    def start_clean():
        Thread(target=clean).setDaemon(False)
        Thread(target=clean).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="计算三点定位",command=start_event)
    b1.bind_all('C',button_event)
    b1.bind_all('c',button_event)
    b1.place(x=20, y=660, 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)
    
    start_event()
    root.mainloop()
    
if __name__ == "__main__":
    gui_start()
    


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/118381.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

[开源工具]使用Fiddler简单计算QPS[新手开箱可用]

使用Fiddler简单计算QPS1.什么是QPS?2.怎么计算QPS?3.如何使用Fiddler得到一个API接口的QPS?3.1配置&#xff1a;打开Fiddler文件夹&#xff0c;点击Fiddler.exe运行fiddler进行配置4.如何得到本机的核心数?5.根据公式计算QPS?6.扩展计算单机可支撑PV(理论值)?1.什么是QP…

springboot中controller层接收参数,servers层调用mapper层,一条sql搞定排序

前言 很多小伙伴们在公司不管是测试C端产品还是B端产品&#xff0c;都会测到排序的业务需求&#xff1b;那么我们就会好奇排序是如何实现的呢&#xff1f;下面我们开始介绍代码的实现 数据库建表 我们需要创建一个书籍book表结构&#xff0c;如下图所示 CREATE TABLE book ( …

嵌入式C语言面向对象编程 --- 总结

什么是 C 语言面向对象? 在开始嵌入式 C 语言设计模式系列文章之前,先通过三篇文章讲述了如何使用 C 语言实现面向对象的三大特性,封装,继承,多态。 图片来源公众号:码农翻身 对于“面向对象”这个词语,相信很多软件工程师都不会感觉到陌生,并且很多软件工程师在刚开…

ABAP: 定义关键字的区别

问题&#xff1a;TYPE、LIKE、LIKE TABLE OF、LIKE LINE OF、TYPE TABLE OF 的区别&#xff1f; 1、TYPE 用于变量的类型定义&#xff0c;可以是表中预定义好的字段&#xff0c;也可以是C(字符)&#xff0c;F(浮点型)&#xff0c;I(整型)等。 例如&#xff1a; DATA: NAME TY…

Spring Cloud简介

一、什么是SpringCloud&#xff1f; Spring Cloud 是一系列框架的有序集合。 Spring Cloud 并没有重复制造轮子&#xff0c;它只是将目前各家公司开发的比较成熟、经得起实际考验的服务框架组合起来。 通过 Spring Boot 风格进行再封装屏蔽掉了复杂的配置和实现原理&#xff…

Python 图像边缘检测 | 利用 opencv 和 skimage 的 Canny 算法

文章目录一、简介二、opencv 实践三、skimage 实践CSDN 叶庭云&#xff1a;https://yetingyun.blog.csdn.net/ 一、简介 提取图片的边缘信息是底层数字图像处理的基本任务之一。边缘信息对进一步提取高层语义信息有很大的影响。大部分边缘检测算法都是上个世纪的了&#xff0c…

【Java基础】day11

day11 一、BIO、NIO、AIO 三种 IO 模型分别是什么&#xff1f; BIO &#xff08;Blocking I/O&#xff09;同步阻塞的 I/O 、NIO&#xff08;New/Non-blocking I/O&#xff09; 同步非阻塞的 I/O 、AIO&#xff08;Asynchronous I/O&#xff09; 异步非阻塞的 I/O 。这三种 I…

前端例程20221227:下雪动画

演示 动图太大了不好上传&#xff0c;这里就放个静态图吧&#xff0c;实际上这里是雪花从上到下飘落的效果。 代码 <!DOCTYPE html> <html lang"en"><head><meta charset"UTF-8"><meta name"viewport" content&quo…

【Spring【AOP】】——20、搭建一个AOP测试环境?

&#x1f4eb;作者简介&#xff1a;zhz小白 公众号&#xff1a;小白的Java进阶之路 专业技能&#xff1a; 1、Java基础&#xff0c;并精通多线程的开发&#xff0c;熟悉JVM原理 2、熟悉Java基础&#xff0c;并精通多线程的开发&#xff0c;熟悉JVM原理&#xff0c;具备⼀定的线…

Linux系统运行时参数命令--CPU性能监控

目录 2 CPU性能监控 2.1 平均负载和CPU使用率 1 平均负载基础 2 使用uptime命令分析平均负载 3 平均负载与 CPU 使用率 4 CPU使用率监测命令 ps查找进程信息 top命令查询进程的cpu、内存信息 mpstat pidstat 场景一&#xff1a;CPU 密集型进程 场景二&#xff1a;I…

嵌入式开发学习之--DMA(上)

提示&#xff1a;本篇文章主要了解学习一下DMA 文章目录前言一、DMA是什么&#xff1f;二、DMA使用步骤1、外设通道选择2.2、仲裁器2.3、FIFO2.4、存储器端口 外设端口2.5、编程端口三、DMA结构体初始化详解总结前言 上一篇文章学习了串口通信&#xff0c;这一篇我们来学习一下…

虹科案例 | 实现了新的核磁共振应用!(下)

摘要 光纤传感器已成为推动MRI最新功能套件升级和新MRI设备设计背后的关键技术。在这篇文章中&#xff0c;我们介绍了三个基于MRI的运动控制应用&#xff0c;展示了最近开发的、可在市场上买到的基于MRI的安全光纤反馈传感器的操作和使用。 案例研究 案例研究#1用于验证MRI技…

【聆思CSK6 视觉AI开发套件试用】体验AI功能

本篇文章来自极术社区与聆思科技组织的CSK6 视觉AI开发套件活动&#xff0c;更多开发板试用活动请关注极术社区网站。作者&#xff1a;YiCheng 很高兴极术社区和聆思科技举办的开发板申请试用活动&#xff0c;放张开箱图&#xff1a; 开发板介绍NanoKit开发板# NanoKit开发板是…

【算法】前缀和与差分

文章目录1.一维前缀和输入格式输出格式数据范围2.二维前缀和输入格式输出格式数据范围3.一维差分输入格式输出格式数据范围4.二维差分输入格式输出格式数据范围前缀和差分是一对逆运算1.一维前缀和 有一个长度为n的数组an:a1,a2…an; 对于前缀和&#xff1a;Si a1a2…ai 如何…

刷爆力扣之二进制求和

刷爆力扣之二进制求和 HELLO&#xff0c;各位看官大大好&#xff0c;我是阿呆 &#x1f648;&#x1f648;&#x1f648; 今天阿呆继续记录下力扣刷题过程&#xff0c;收录在专栏算法中 &#x1f61c;&#x1f61c;&#x1f61c; 该专栏按照不同类别标签进行刷题&#xff0c;…

Spark零基础入门实战(八)Spark运行时架构

Spark有多种运行模式,可以运行在一台机器上,称为本地(单机)模式,也可以以YARN或Mesos作为底层资源调度系统以分布式的方式在集群中运行,称为Spark On YARN模式,还可以使用Spark自带的资源调度系统,称为Spark Standalone模式。 本地模式通过多线程模拟分布式计算,通常…

如何使用Python创建一个虚拟助理

使用Python创建一个虚拟助手 虚拟助理是基于人工智能的程序。它们是一种智能计算机程序&#xff0c;可以通过语音命令或文本理解人类的自然语言&#xff0c;并为用户执行任务。本文将指导你通过使用Python库来创建你自己的语音助手的过程。 在本教程中&#xff0c;我将在Jetb…

2022已经到了尾声,做一个漂亮的倒计时页面迎接2023年的到来

个人主页&#xff1a;天寒雨落的博客_CSDN博客-初学者入门C语言,python,数据库领域博主 目录 前言 效果图 网页直链 完整代码 倒计时的实现原理 前言 2022即将到来&#xff0c;快来一起做一个元旦倒计时分享给自己喜欢的人吧~ 效果图 网页直链 下面是可以直接访问的网页…

离散数学实践三判别图的连通性【JAVA实现】

文章目录思路重要算法&Warshall 算法代码全部代码效果展示思路 对于给定的邻接矩阵 A&#xff0c;可以用可达矩阵 Warshall 算法求出 A 所表示的图的可达矩阵 P。 对于可达矩阵 P 来说&#xff0c;如果 P 的所有元素均为 111&#xff0c;则所给的有向图是强连通的&#xf…

Python实现的一个简洁轻快的后台管理框架.支持拥有多用户组的RBAC管理后台,不用配置各种运行环境

Mini Admin 完整代码下载地址&#xff1a;Python实现的一个简洁轻快的后台管理框架.支持拥有多用户组的RBAC管理后台 Mini Admin,一个简洁轻快的后台管理框架.支持拥有多用户组的RBAC管理后台 &#x1f680; 应用场景&#xff1a;2-5人的管理团队&#xff0c;需要管理的资源数…