Merge remote-tracking branch 'origin/master'

Back/GS/GS.py

@@ -1,2 +1,400 @@
-def main_function(file_path, db_path):
-    print("这是GS")
+import openpyxl
+import xlrd
+import sqlite3
+import os
+from openpyxl import Workbook
+import tkinter as tk
+from tkinter import messagebox
+import math
+import numpy as np
+
+def to_utf8(text):
+    str1 = text.encode('utf-8')
+    return str1
+
+# 弹窗提示用户是否更新数据
+def ask_for_update(file_name):
+    root = tk.Tk()
+    root.withdraw()  # 隐藏主窗口
+    response = messagebox.askyesno("提示", f"检测到数据库中存在相同文件名 '{file_name}'，是否更新数据？")
+    return response
+
+# 读取Excel文件并计算公式结果
+def load_and_calculate_excel(file_path,cell0):
+    # 加载Excel文件并计算单元格的公式结果。
+    workbook = openpyxl.load_workbook(file_path, data_only=True)
+    sheet = workbook.active
+    h1_value = sheet[cell0].value
+    return h1_value
+
+# 数据库插入函数
+def insert_into_database(database,pastname,newname,excelname,listname1,listpastx1,listpasty1,listcgcs1,listnewx1,listnewy1,pjbc,fxzwc,zrbzwc,points,zbs,zfxs,sf,pjbcs,pjfxs):
+    #先把输入的录入数据库
+    db = sqlite3.connect(database)
+    #获取游标
+    cursor = db.cursor()
+    #字符类的都需要转换成字节存进去
+    utf_pan = to_utf8(pastname)
+    utf_nn = to_utf8(newname)
+    utf_tn = to_utf8(excelname)
+    try:
+        # 查询是否已存在相同的记录
+        cursor.execute(
+            "SELECT * FROM GS_Input_Param WHERE TableName = ?",
+            (utf_tn,)
+        )
+        existing_record = cursor.fetchone()
+
+        if existing_record:
+            # 弹窗询问用户是否更新数据
+            if not ask_for_update(excelname):
+                print("用户选择不更新数据")
+                return
+
+            # 更新现有记录GS_Input_Param 
+            cursor.execute(
+                """
+                UPDATE GS_Input_Param 
+                SET Last_ResultName=?,New_ResultName=?,Avg_SL=?,Ms_Dir=?,Ms_WSL=?,Pt_Count=?,SL_Count=?,Dir_Count=?,Scale_Value=?,Avg_MSL=?,Avg_Dir=?
+                WHERE TableName = ?
+                """,
+                (utf_pan,utf_nn,pjbc,fxzwc,zrbzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,utf_tn,)
+            )
+            # 更新现有记录GS_Input_Point 
+            #删了已有的数据，重新插入
+            cursor.execute(
+                """
+                DELETE FROM GS_Input_Point WHERE TableName = ?
+                """,
+                (utf_tn,)
+            )
+            rr = 0
+            while rr < len(listname1):
+                rr1 = listname1[rr]
+                rr2 = listpastx1[rr]
+                rr3 = listpasty1[rr]
+                rr4 = listcgcs1[rr]
+                rr5 = listnewx1[rr]
+                rr6 = listnewy1[rr]
+                utf_rr1 = to_utf8(rr1)
+                cursor.execute(
+                    """
+                    INSERT INTO GS_Input_Point(TableName,Last_ResultName,New_ResultName,Last_X,Last_Y,New_X,New_Y,PointName,cgcs) VALUES (?,?,?,?,?,?,?,?,?)
+                    """,
+                    (utf_tn,utf_pan,utf_nn,rr2,rr3,rr5,rr6,utf_rr1,rr4,)
+                )
+                rr = rr + 1
+            messagebox.showinfo("提示",
+                                f"文件 '{excelname}' 已成功更新到本地数据库中，点击确认以关闭此对话框。对话框关闭后，请点击“计算”以重新计算数据")
+        else:
+            # 插入新记录GS_Input_Param 
+            cursor.execute(
+                """
+                INSERT INTO GS_Input_Param (TableName,Last_ResultName,New_ResultName,Avg_SL,Ms_Dir,Ms_WSL,Pt_Count,SL_Count,Dir_Count,Scale_Value,Avg_MSL,Avg_Dir) 
+                VALUES (?,?,?,?,?,?,?,?,?,?,?,?)
+                """,
+                (utf_tn,utf_pan,utf_nn,pjbc,fxzwc,zrbzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,)
+            )
+            # 插入新记录GS_Input_Point 
+            rr = 0
+            while rr < len(listname1):
+                rr1 = listname1[rr]
+                rr2 = listpastx1[rr]
+                rr3 = listpasty1[rr]
+                rr4 = listcgcs1[rr]
+                rr5 = listnewx1[rr]
+                rr6 = listnewy1[rr]
+                utf_rr1 = to_utf8(rr1)
+                cursor.execute(
+                    """
+                    INSERT INTO GS_Input_Point(TableName,Last_ResultName,New_ResultName,Last_X,Last_Y,New_X,New_Y,PointName,cgcs) VALUES (?,?,?,?,?,?,?,?,?)
+                    """,
+                    (utf_tn,utf_pan,utf_nn,rr2,rr3,rr5,rr6,utf_rr1,rr4,)
+                )
+                rr = rr + 1
+            messagebox.showinfo("提示",
+                                f"文件 '{excelname}' 已成功添加到本地数据库中，点击确认以关闭此对话框。对话框关闭后，请点击“计算”以计算数据")
+
+        db.commit()
+    except sqlite3.Error as e:
+        print(f"插入或更新文件名时发生错误: {e}")
+    except Exception as e:
+        print(f"处理文件时发生错误: {e}")
+    finally:
+        db.close()
+
+#读取xls
+def xlrd_read(excelpath,dbpath):
+    global gszbx
+    global gszby
+    listname1 = []
+    listnewx1 = []
+    listnewy1 = []
+    listpastx1 = []
+    listpasty1 = []
+    listcgcs1 = []
+
+    #excel表名就是表名
+    excelname = os.path.basename(excelpath)
+    # 读取excel
+    xlr = xlrd.open_workbook(excelpath)
+    tabler = xlr.sheets()[0]
+    # 从第3行开始，3A为点号，3B为最新x，3C为最新y，3D为往期x
+    # 3E为往期y
+    # 先看点数
+    rows = tabler.nrows
+    row = 2
+    points = rows - 2
+    # 收集每一列的数据，方便后期计算
+    listname = []
+    listnewx = []
+    listnewy = []
+    listpastx = []
+    listpasty = []
+    listcgcs = []
+    sumnewx = 0
+    sumnewy = 0
+    sumpastx = 0
+    sumpasty = 0
+    # 记录下最早的那几个开头
+    # excelpath是输入表格，outpath是输出路径，pjbc是平均边长
+    # fxzwc是方向中误差，zrbzwc是最弱边边长相对中误差，pjbcs是平均边数，pjfxs是平均方向数，sf是缩放值
+    # 新增总边数zbs，总方向数zfxs
+    # 先计算位移判断值
+    pjbc = float(tabler.cell(0, 6).value)
+    fxzwc = float(tabler.cell(1, 6).value)
+    str1 = str(tabler.cell(2, 6).value)
+    #判断下是哪种情况，转化为数值
+    if '.' in str1:
+        #直接float使用
+        zrbzwc = float(str1)
+    else:
+        #转化为分数
+        # 10的6次方
+        # ten_to_the_six = 10e6
+        fzstr = float(str1.split('/',-1)[0])
+        fmstr = float(str1.split('/',-1)[1])
+        zrbzwc = float(fzstr/fmstr)
+    zbs = float(tabler.cell(3, 6).value)
+    zfxs = float(tabler.cell(4, 6).value)
+    pjbcs = zbs / points
+    pjfxs = zfxs / points
+    sf = float(tabler.cell(5,6).value)
+    # 再拆细点
+    wy1 = pjbc / 1000
+    wy2 = wy1 * zrbzwc * 1e6
+    wypd1 = wy2 * wy2 / pjbcs
+    # 206.265是常数系数
+    wy3 = pjbc * fxzwc / 206.265
+    wypd2 = wy3 * wy3 / pjfxs
+    wypd3 = math.sqrt(wypd1 + wypd2)
+    wypd = round(wypd3 * 3, 4)
+    wypd0 = wypd
+    pastname = tabler.cell(0, 3).value
+    newname = tabler.cell(0, 1).value
+    #如果第二次测量或者首次测量为空，则判断为新建
+    while row < rows:
+        pname = tabler.cell(row, 0).value
+        pcgcs = 1
+        try:
+            pnewx = round(float(tabler.cell(row, 1).value),4)
+        except:
+            pnewx = 0
+            pcgcs = -2
+        try:
+            pnewy = round(float(tabler.cell(row, 2).value),4)
+        except:
+            pnewy = 0
+            pcgcs = -2
+        try:
+            ppastx = round(float(tabler.cell(row, 3).value),4)
+        except:
+            ppastx = 0
+            pcgcs = -2
+        try:
+            ppasty = round(float(tabler.cell(row, 4).value),4)
+        except:
+            ppasty = 0
+            pcgcs = -2
+        sumpastx = sumpastx + ppastx
+        sumnewy = sumnewy + pnewy
+        sumnewx = sumnewx + pnewx
+        sumpasty = sumpasty + ppasty
+        listname1.append(pname)
+        listnewx1.append(pnewx)
+        listnewy1.append(pnewy)
+        listpastx1.append(ppastx)
+        listpasty1.append(ppasty)
+        listcgcs1.append(pcgcs)
+        if pcgcs == -2:
+            points = points - 1
+        row = row + 1
+    #插入数据库
+    insert_into_database(dbpath,pastname,newname,excelname,listname1,listpastx1,listpasty1,listcgcs1,listnewx1,listnewy1,pjbc,fxzwc,zrbzwc,points,zbs,zfxs,sf,pjbcs,pjfxs)
+    #新建点不参与运算
+    rr=0
+    while rr < len(listname1):
+        rr1 = listname1[rr]
+        rr2 = listpastx1[rr]
+        rr3 = listpasty1[rr]
+        rr4 = listcgcs1[rr]
+        rr5 = listnewx1[rr]
+        rr6 = listnewy1[rr]
+        if rr4 == -2:
+            rr = rr + 1
+        else:
+            listname.append(rr1)
+            listpastx.append(rr2)
+            listpasty.append(rr3)
+            listcgcs.append(rr4)
+            listnewx.append(rr5)
+            listnewy.append(rr6)
+            rr = rr + 1
+
+#读取xlsx
+def openpyxl_read(excelpath,dbpath):
+    global gszbx
+    global gszby
+    listname1 = []
+    listnewx1 = []
+    listnewy1 = []
+    listpastx1 = []
+    listpasty1 = []
+    listcgcs1 = []
+
+    #excel表名就是表名
+    excelname = os.path.basename(excelpath)
+    # 读取excel
+    xlr = openpyxl.load_workbook(excelpath)
+    tabler = xlr.worksheets[0]
+    # 从第3行开始，3A为点号，3B为最新x，3C为最新y，3D为往期x
+    # 3E为往期y
+    # 先看点数
+    rows = tabler.max_row
+    row = 3
+    points = rows - 2
+    # 收集每一列的数据，方便后期计算
+    listname = []
+    listnewx = []
+    listnewy = []
+    listpastx = []
+    listpasty = []
+    listcgcs = []
+    sumnewx = 0
+    sumnewy = 0
+    sumpastx = 0
+    sumpasty = 0
+    # 记录下最早的那几个开头
+    # excelpath是输入表格，outpath是输出路径，pjbc是平均边长
+    # fxzwc是方向中误差，zrbzwc是最弱边边长相对中误差，pjbcs是平均边数，pjfxs是平均方向数，sf是缩放值
+    # 新增总边数zbs，总方向数zfxs
+    # 先计算位移判断值
+    pjbc = float(tabler['G1'].value)
+    fxzwc = float(tabler['G2'].value)
+    str1 = str(tabler['G3'].value)
+    #判断下是哪种情况，转化为数值
+    if '.' in str1:
+        #直接float使用
+        zrbzwc = float(str1)
+    else:
+        #转化为分数
+        # 10的6次方
+        # ten_to_the_six = 10e6
+        fzstr = float(str1.split('/',-1)[0])
+        fmstr = float(str1.split('/',-1)[1])
+        zrbzwc = float(fzstr/fmstr)
+    try:
+        zbs = float(tabler['G4'].value)
+    except:
+        zbs = load_and_calculate_excel(excelpath,'G4')
+    try:
+        zfxs = float(tabler['G5'].value)
+    except:
+        zfxs = load_and_calculate_excel(excelpath,'G5')
+    pjbcs = zbs / points
+    pjfxs = zfxs / points
+    sf = float(tabler['G6'].value)
+    # 再拆细点
+    wy1 = pjbc / 1000
+    wy2 = wy1 * zrbzwc * 1e6
+    wypd1 = wy2 * wy2 / pjbcs
+    # 206.265是常数系数
+    wy3 = pjbc * fxzwc / 206.265
+    wypd2 = wy3 * wy3 / pjfxs
+    wypd3 = math.sqrt(wypd1 + wypd2)
+    wypd = round(wypd3 * 3, 4)
+    wypd0 = wypd
+    pastname = tabler['D1'].value
+    newname = tabler['B1'].value
+    #如果第二次测量或者首次测量为空，则判断为新建
+    while row <= rows:
+        pname = tabler.cell(row, 1).value
+        pcgcs = 1
+        try:
+            pnewx = round(float(tabler.cell(row, 2).value),4)
+        except:
+            pnewx = 0
+            pcgcs = -2
+        try:
+            pnewy = round(float(tabler.cell(row, 3).value),4)
+        except:
+            pnewy = 0
+            pcgcs = -2
+        try:
+            ppastx = round(float(tabler.cell(row, 4).value),4)
+        except:
+            ppastx = 0
+            pcgcs = -2
+        try:
+            ppasty = round(float(tabler.cell(row, 5).value),4)
+        except:
+            ppasty = 0
+            pcgcs = -2
+        sumpastx = sumpastx + ppastx
+        sumnewy = sumnewy + pnewy
+        sumnewx = sumnewx + pnewx
+        sumpasty = sumpasty + ppasty
+        listname1.append(pname)
+        listnewx1.append(pnewx)
+        listnewy1.append(pnewy)
+        listpastx1.append(ppastx)
+        listpasty1.append(ppasty)
+        listcgcs1.append(pcgcs)
+        if pcgcs == -2:
+            points = points - 1
+        row = row + 1
+    #插入数据库
+    insert_into_database(dbpath,pastname,newname,excelname,listname1,listpastx1,listpasty1,listcgcs1,listnewx1,listnewy1,pjbc,fxzwc,zrbzwc,points,zbs,zfxs,sf,pjbcs,pjfxs)
+    #新建点不参与运算
+    rr=0
+    while rr < len(listname1):
+        rr1 = listname1[rr]
+        rr2 = listpastx1[rr]
+        rr3 = listpasty1[rr]
+        rr4 = listcgcs1[rr]
+        rr5 = listnewx1[rr]
+        rr6 = listnewy1[rr]
+        if rr4 == -2:
+            rr = rr + 1
+        else:
+            listname.append(rr1)
+            listpastx.append(rr2)
+            listpasty.append(rr3)
+            listcgcs.append(rr4)
+            listnewx.append(rr5)
+            listnewy.append(rr6)
+            rr = rr + 1
+
+# 主函数 读取excel文件
+def main_function(file_path,dbpath):
+    # 调用读取Excel文件的函数
+    file_name = os.path.basename(file_path)
+    #两种加载方式，对于不同的文件
+    if ".xlsx" in file_path:
+        #采用openpyxl
+        openpyxl_read(file_path,dbpath)
+    else:
+        #采用xlrd
+        xlrd_read(file_path,dbpath)
+    
+    

Back/GS/GScompute.py

@@ -1,2 +1,547 @@
-def main_function(file_path, db_path):
-    print("GS")
+import sqlite3
+import os
+import tkinter as tk
+from tkinter import messagebox
+import math
+import numpy as np
+
+#region 各种方法
+
+def to_utf8(text):
+    str1 = text.encode('utf-8')
+    return str1
+
+def sfjs(gsx,listx,gsy,listy):
+    s1 = gsx[len(gsx) - 1] - gsx[0]
+    s2 = s1 * s1
+    s3 = gsy[len(gsy) - 1] - gsy[0]
+    s4 = s3 * s3
+    s5 = s2 + s4
+    s6 = math.sqrt(s5)
+    t1 = listx[len(listx) - 1] - listx[0]
+    t2 = t1 * t1
+    t3 = listy[len(listy) - 1] - listy[0]
+    t4 = t3 * t3
+    t5 = t2 + t4
+    t6 = math.sqrt(t5)
+    num1 = round(s6/t6,8)
+    return num1
+
+def mxjd(gs,nw):
+    ii = 0
+    sum1 = 0
+    while ii < len(gs):
+        s1 = gs[ii] - nw[ii]
+        s2 = s1 * s1
+        sum1 = sum1 + s2
+        ii = ii + 1
+    s3 = sum1 / 2 / len(gs)
+    s4 = math.sqrt(s3)
+    num1 = s4 * 1000
+    return num1
+
+def xzjs(gsx,listx,gsy,listy):
+    s1 = gsx[len(gsx) - 1] - gsx[0]
+    s2 = gsy[len(gsy) - 1] - gsy[0]
+    s3 = s2 / s1
+    s4 = math.atan(s3)
+    t1 = listx[len(listx) - 1] - listx[0]
+    t2 = listy[len(listy) - 1] - listy[0]
+    t3 = t2 / t1
+    t4 = math.atan(t3)
+    n1 = s4 - t4
+    num1 = n1 * 180 / math.pi * 3600
+    return num1
+
+def gsjs(x,y,nlist):
+    listgsx = []
+    xstr1 = nlist[0] * x
+    xstr2 = nlist[1] * y
+    xx = xstr1 + xstr2 + nlist[2]
+    ystr1 = nlist[3] * x
+    ystr2 = nlist[4] * y
+    yy = ystr1 + ystr2 + nlist[5]
+    listgsx.append(xx)
+    listgsx.append(yy)
+    return listgsx
+
+def jzys1(listy, zxzby, listx, zxzbx, sf):
+    rlist = []
+    listlen = len(listy)
+    ll = 0
+    while ll < listlen:
+        # 第一个的值
+        # num1 = round(((listx[ll] - zxzbx) / sf),0)
+        # num2 = round(((listy[ll] - zxzby) / sf),0)
+        num1 = (listx[ll] - zxzbx)
+        num2 = (listy[ll] - zxzby)
+        rlist.append(num1)
+        rlist.append(num2)
+        ll = ll + 1
+    return rlist
+
+def jzys2(listy, zxzby, listx, zxzbx, sf):
+    rlist = []
+    listlen = len(listy)
+    ll = 0
+    while ll < listlen:
+        # 第一个的值
+        num2 = (-1) * (listx[ll] - zxzbx) / sf
+        num1 = (listy[ll] - zxzby) / sf
+        rlist.append(num1)
+        rlist.append(num2)
+        ll = ll + 1
+    return rlist
+
+def jzys3(listy, aa, bb):
+    rlist = []
+    listlen = len(listy)
+    ll = 0
+    while ll < listlen:
+        # 第一个的值
+        num1 = aa
+        num2 = bb
+        rlist.append(num1)
+        rlist.append(num2)
+        ll = ll + 1
+    return rlist
+
+def jzys4(listpasty, zxzbpasty, listpastx, zxzbpastx, listnewy, zxzbnewy, listnewx, zxzbnewx, sf):
+    rlist = []
+    listlen = len(listnewx)
+    ll = 0
+    while ll < listlen:
+        num1 = (listpastx[ll] - listnewx[ll] - zxzbpastx + zxzbnewx) / sf
+        num2 = (listpasty[ll] - listnewy[ll] - zxzbpasty + zxzbnewy) / sf
+        rlist.append(round(num1,6))
+        rlist.append(round(num2,6))
+        ll = ll + 1
+    return rlist
+
+def jzys5(listp, zxzbp, sf):
+    rlist = []
+    listlen = len(listp)
+    ll = 0
+    while ll < listlen:
+        num1 = (listp[ll] - zxzbp) / sf
+        rlist.append(num1)
+        ll = ll + 1
+    return rlist
+
+def jzys6(listp, num):
+    rlist = []
+    listlen = len(listp)
+    ll = 0
+    while ll < listlen:
+        num1 = listp[ll] * num
+        rlist.append(num1)
+        ll = ll + 1
+    return rlist
+
+def cjh(listp):
+    rp = 0
+    listlen = len(listp)
+    ll = 0
+    while ll < listlen:
+        num1 = listp[ll] * listp[ll]
+        rp = rp + num1
+        ll = ll + 1
+    return rp
+
+def cjh2(lista, listb):
+    rp = 0
+    listlen = len(lista)
+    ll = 0
+    while ll < listlen:
+        num1 = lista[ll] * listb[ll]
+        rp = rp + num1
+        ll = ll + 1
+    return rp
+
+def gsys(listnewp, jxlist1, jylist1, lxlist, lylist, arrz, sf, zxzbnewp, zxzbpastp, ii):
+    # 新x+（矩阵[jx1 jy1 lx1 ly1]*矩阵z）*缩放-新重心坐标x+旧重心坐标x
+    rlist = []
+    listlen = len(listnewp)
+    ll = 0
+    while ll < listlen:
+        arr0 = np.array((jxlist1[ll], jylist1[ll],
+                        lxlist[2*ll+ii], lylist[2*ll+ii]))
+        arr1 = np.matmul(arr0, arrz)
+        arr3 = sf * arr1
+        newp = listnewp[ll] + arr3 - zxzbnewp + zxzbpastp
+        rlist.append(newp)
+        ll = ll + 1
+    return rlist
+
+def xcys(newlist, pastlist):
+    listlen = len(newlist)
+    ll = 0
+    rlist = []
+    while ll < listlen:
+        num1 =pastlist[ll] - newlist[ll]
+        num2 = num1 * 1000
+        rlist.append(num2)
+        ll = ll + 1
+    return rlist
+
+def takeFirst(elem):
+        return elem[0]
+
+#endregion
+
+def insert_into_database1(database,utf_tn,listname1,listname,gsx,gsy,listpastx1,listpasty1,listnewx1,listnewy1,nlist,utf_pan,utf_nn,wypd,sfxs,xzj,pycsx,pycsy,mxjdx,mxjdy,zxzbpastx,zxzbpasty,zxzbnewx,zxzbnewy,n1,n2,n3,n4,n5,n6):
+    #将结果输出到数据库
+    db1 = sqlite3.connect(database)
+    #获取游标
+    cursor1 = db1.cursor()
+    #先清除已有数据，用来更新
+    try:
+        sqlstr3 = 'delete from GS_Trans_Param WHERE TableName = ?'
+        sqlstr4 = 'delete from GS_Trans_Point WHERE TableName = ?'
+        sqlstr5 = 'delete from GS_Result_Point WHERE TableName = ?'
+        cursor1.execute(sqlstr3,(utf_tn,))
+        cursor1.execute(sqlstr4,(utf_tn,))
+        cursor1.execute(sqlstr5,(utf_tn,))
+    except:
+        pass
+    cursor1.execute('INSERT INTO GS_Trans_Param(TableName,Last_ResultName,New_ResultName,Dis_RefValue,Pt_CalCount,Scale_Factor,Rotate_Angle,Trans_X,Trans_Y,Model_AccX,Model_AccY,Last_GcX,Last_GcY,New_GcX,New_GcY,Formula_X1,Formula_X2,Formula_X3,Formula_Y1,Formula_Y2,Formula_Y3) VALUES(?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,wypd,len(listname),sfxs,xzj,pycsx,pycsy,mxjdx,mxjdy,zxzbpastx,zxzbpasty,zxzbnewx,zxzbnewy,n1,n2,n3,n4,n5,n6,))
+    #Trans_Point和Result_Point一起存
+    for iname in listname1:
+        if iname in listname:
+            #属于没问题点的
+            #直接使用gsx，gsy
+            utf_pointname = to_utf8(str(iname))
+            wystr = to_utf8('稳定')
+            index1 = listname.index(iname)
+            numm1 = round(gsx[index1],4)
+            numm2 = round(gsy[index1],4)
+            #先找到首次的index
+            index2 = listname1.index(iname)
+            r1 = (listpastx1[index2] - gsx[index1])*1000
+            r2 = (listpasty1[index2] - gsy[index1])*1000
+            r3 = r1 * r1
+            r4 = r2 * r2
+            r5 = r3 + r4
+            r6 = round(math.sqrt(r5),1)
+            numn1 = listpastx1[index2]
+            numn2 = listnewx1[index2]
+            cursor1.execute('INSERT INTO GS_Trans_Point(TableName,Last_ResultName,New_ResultName,PointName,Last_X,Last_Y,New_X,New_Y) VALUES (?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,utf_pointname,listpastx1[index2],listpasty1[index2],listnewx1[index2],listnewy1[index2],))
+            cursor1.execute('INSERT INTO GS_Result_Point(TableName,Last_ResultName,New_ResultName,PointName,Last_X,Last_Y,Result_X,Result_Y,Cal_X,Cal_Y,Last_ResultX,Last_ResultY,Last_ResultP,Last_CalX,Last_CalY,Last_CalP,Dis_Ass) VALUES (?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,utf_pointname,listpastx1[index2],listpasty1[index2],listpastx1[index2],listpasty1[index2],numm1,numm2,0,0,0,r1,r2,r6,wystr,))
+        else:
+            index2 = listname1.index(iname)
+            utf_pointname = to_utf8(str(iname))
+            #首先判断有没有第一个点
+            if listnewx1[index2] != 0:
+                #有问题的点，使用公式跑
+                listxy = gsjs(listnewx1[index2],listnewy1[index2],nlist)   
+                #判断有没有第一次的数据，没有就不写
+                if listpastx1[index2] != 0:
+                    r1 = (listpastx1[index2] - listxy[0])*1000
+                    r2 = (listpasty1[index2] - listxy[1])*1000
+                    r3 = r1 * r1
+                    r4 = r2 * r2
+                    r5 = r3 + r4
+                    r6 = round(math.sqrt(r5),1)
+                    if r6 > 500:
+                        wystr = to_utf8('复建')
+                        cursor1.execute('INSERT INTO GS_Result_Point(TableName,Last_ResultName,New_ResultName,PointName,Last_X,Last_Y,Result_X,Result_Y,Cal_X,Cal_Y,Last_ResultX,Last_ResultY,Last_ResultP,Last_CalX,Last_CalY,Last_CalP,Dis_Ass) VALUES (?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,utf_pointname,listpastx1[index2],listpasty1[index2],listxy[0],listxy[1],listxy[0],listxy[1],r1,r2,r6,r1,r2,r6,wystr,))
+                    else:
+                        wystr = to_utf8('位移')
+                        cursor1.execute('INSERT INTO GS_Result_Point(TableName,Last_ResultName,New_ResultName,PointName,Last_X,Last_Y,Result_X,Result_Y,Cal_X,Cal_Y,Last_ResultX,Last_ResultY,Last_ResultP,Last_CalX,Last_CalY,Last_CalP,Dis_Ass) VALUES (?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,utf_pointname,listpastx1[index2],listpasty1[index2],listxy[0],listxy[1],listxy[0],listxy[1],r1,r2,r6,r1,r2,r6,wystr,))
+                else:
+                    wystr = to_utf8('新建')
+                    cursor1.execute('INSERT INTO GS_Result_Point(TableName,Last_ResultName,New_ResultName,PointName,Last_X,Last_Y,Result_X,Result_Y,Cal_X,Cal_Y,Last_ResultX,Last_ResultY,Last_ResultP,Last_CalX,Last_CalY,Last_CalP,Dis_Ass) VALUES (?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)',(utf_tn,utf_pan,utf_nn,utf_pointname,listpastx1[index2],listpasty1[index2],listnewx1[index2],listnewy1[index2],listpastx1[index2],listpasty1[index2],0,0,0,0,0,0,wystr,))
+    #数据库执行
+    db1.commit()
+    #做完一切后，先关闭游标，再关闭数据库
+    cursor1.close()
+    db1.close()
+
+def tablein(outpath,str1):
+    listcgcs1 = []
+    listname1 = []
+    listpastx1 = []
+    listpasty1 = []
+    listname = []
+    listnewx = []
+    listnewy = []
+    listpastx = []
+    listpasty = []
+    listcgcs = []
+    listnewx1 = []
+    listnewy1 = []
+    points = 0
+    #查询，提取需要的数据
+    db2 = sqlite3.connect(outpath)
+    cursor2 = db2.cursor()
+    utfstr1 = to_utf8(str1)
+    sqlstr1 = 'SELECT * FROM GS_Input_Point WHERE TableName = ?'
+    cursor2.execute(sqlstr1,(utfstr1,))
+    all_da = cursor2.fetchall()
+    for da in all_da:
+        listpastx.append(da[3])
+        listpastx1.append(da[3])
+        listpasty.append(da[4])
+        listpasty1.append(da[4])
+        listnewx.append(da[5])
+        listnewx1.append(da[5])
+        listnewy.append(da[6])
+        listnewy1.append(da[6])
+        listname.append(da[7].decode('utf-8'))
+        listname1.append(da[7].decode('utf-8'))
+        listcgcs.append(da[8])
+        listcgcs1.append(da[8])
+        points = points + 1
+    sqlstr2 = 'SELECT * FROM GS_Input_Param WHERE TableName = ?'
+    cursor2.execute(sqlstr2,(utfstr1,))
+    all_par = cursor2.fetchall()
+    pastname = all_par[0][1].decode('utf-8')
+    newname = all_par[0][2].decode('utf-8')
+    # pjbc是平均边长
+    # fxzwc是方向中误差，zrbzwc是最弱边边长相对中误差，pjbcs是平均边数，pjfxs是平均方向数，sf是缩放值
+    # 新增总边数zbs，总方向数zfxs
+    # 先计算位移判断值
+    pjbc = all_par[0][3]
+    fxzwc = all_par[0][4]
+    zrbzwc = all_par[0][5]
+    zbs = all_par[0][7]
+    zfxs = all_par[0][8]
+    sf = all_par[0][9]
+    pjbcs = zbs / points
+    pjfxs = zfxs / points
+    # 再拆细点
+    wy1 = pjbc / 1000
+    wy2 = wy1 * zrbzwc * 1e6
+    wypd1 = wy2 * wy2 / pjbcs
+    # 206.265是常数系数
+    wy3 = pjbc * fxzwc / 206.265
+    wypd2 = wy3 * wy3 / pjfxs
+    wypd3 = math.sqrt(wypd1 + wypd2)
+    wypd = round(wypd3 * 3, 4)
+    wypd0 = wypd
+    #更新下数据库（参数部分）
+    sqlstr3 = 'update GS_Input_Param set Pt_Count=?,Avg_MSL=?,Avg_Dir=? WHERE TableName = ?'
+    cursor2.execute(sqlstr3,(points,pjbcs,pjfxs,utfstr1,))
+    #数据库执行
+    db2.commit()
+    #做完一切后，先关闭游标，再关闭数据库
+    cursor2.close()
+    db2.close()
+    js = 0
+    return listcgcs1, listname1, listpastx1, listpasty1, wypd0, listname, listnewx, listnewy, listpastx,listpasty, points, listcgcs,sf, newname, pastname, js,listnewx1,listnewy1
+
+def bhjs(dbpath,listcgcs1, listname1, listpastx1, listpasty1, wypd, listname, listnewx, listnewy, listpastx, listpasty, points, listcgcs,  sf, newname, pastname, js,listnewx1,listnewy1,excelname):
+    np.set_printoptions(suppress=False)
+    # pt用于给点数point计数
+    # point会随着回递函数而减少
+    pt = 0
+    sumnewx = 0
+    sumnewy = 0
+    sumpastx = 0
+    sumpasty = 0
+    gszbx = []
+    gszby = []
+    arrz = []
+    # #一个计数，把第一次改算坐标全部记下来
+    # 存的第一次的改算值和起限差
+    jslist = []
+    # 等同于listnewx/y1
+    jslist1 = []
+    # listnewx/y,listpastx/y,listcgcs都是要变的
+    # 带1的是原始的
+    while pt < points:
+        sumnewx = sumnewx + listnewx[pt]
+        sumnewy = sumnewy + listnewy[pt]
+        sumpastx = sumpastx + listpastx[pt]
+        sumpasty = sumpasty + listpasty[pt]
+        pt = pt + 1
+    # 计算4个重心坐标（最新一期和上一期的）
+    zxzbnewx = (sumnewx / points)
+    zxzbnewy = (sumnewy / points)
+    zxzbpastx = (sumpastx / points)
+    zxzbpasty = (sumpasty / points)
+    # 先计算平移参数
+    # pycsx = abs(zxzbnewx - zxzbpastx) * 1000
+    # pycsy = abs(zxzbnewy - zxzbpasty) * 1000
+    pycsx = round((zxzbpastx - zxzbnewx),4)
+    pycsy = round((zxzbpasty - zxzbnewy),4)
+    # 分别得k，sita的值
+    klist = jzys1(listnewy, zxzbnewy, listnewx, zxzbnewx, sf)
+    slist = jzys2(listnewy, zxzbnewy, listnewx, zxzbnewx, sf)
+    # 得dx，dy
+    dxlist = jzys3(listnewy, 1, 0)
+    dylist = jzys3(listnewy, 0, 1)
+    # 矩阵太难用了，还是手动计算
+    # 分别计算k、s、dx、dy的乘积和
+    k2 = cjh(klist)
+    s2 = cjh(slist)
+    dx2 = cjh(dxlist)
+    dy2 = cjh(dylist)
+    # 这里再创建矩阵1
+    arr1 = np.array(((k2, 0, 0, 0), (0, s2, 0, 0),
+                    (0, 0, dx2, 0), (0, 0, 0, dy2)))
+    # 矩阵1的逆矩阵矩阵2
+    arr2 = np.linalg.inv(arr1)
+    # 求矩阵l
+    llist = jzys4(listpasty, zxzbpasty, listpastx, zxzbpastx,
+                    listnewy, zxzbnewy, listnewx, zxzbnewx, sf)
+    # 得到数列e1，e2，e3，e4
+    e1 = cjh2(klist, llist)
+    e2 = cjh2(slist, llist)
+    e3 = cjh2(dxlist, llist)
+    e4 = cjh2(dylist, llist)
+    arrl = np.array(((e1), (e2), (e3), (e4)))
+    # 得到矩阵z
+    arrz = np.matmul(arr2, arrl)
+    # 求转换矩阵
+    jxlist1 = jzys5(listnewx, zxzbnewx, sf)
+    jylist1 = jzys5(listnewy, zxzbnewy, sf)
+    jxlist2 = jylist1
+    # jylist2 = jxlist1 * (-1)
+    jylist2 = jzys6(jxlist1, -1)
+    lxlist = jzys3(jylist1, 1, 0)
+    lylist = jzys3(jylist1, 0, 1)
+    # 求改算坐标
+    gsx = gsys(listnewx, jxlist1, jylist1, lxlist,
+                lylist, arrz, sf, zxzbnewx, zxzbpastx, 0)
+    gsy = gsys(listnewy, jxlist2, jylist2, lxlist,
+                lylist, arrz, sf, zxzbnewy, zxzbpasty, 1)
+    #缩放系数也是计算出来的sfxs
+    sfxs = sfjs(gsx,listnewx,gsy,listnewy)
+    #旋转角，模型精度
+    xzj = xzjs(gsx,listnewx,gsy,listnewy)
+    mxjdx = mxjd(gsx,listpastx)
+    mxjdy = mxjd(gsy,listpasty)
+    # 还是要先求较差
+    xcx = xcys(gsx, listpastx)
+    xcy = xcys(gsy, listpasty)
+    # 这里记录下第一次所有点的改算坐标
+    if js == 0:
+        lengs = len(gsx)
+        xx = 0
+        while xx < lengs:
+            gszbx.append(gsx[xx])
+            gszby.append(gsy[xx])
+            zblist1 = []
+            zblist1.append(gsx[xx])
+            zblist1.append(gsy[xx])
+            zblist1.append(xcx[xx])
+            zblist1.append(xcy[xx])
+            jslist.append(zblist1)
+            zblist2 = []
+            zblist2.append(listnewx[xx])
+            zblist2.append(listnewy[xx])
+            jslist1.append(zblist2)
+            xx = xx + 1
+        # 后面不需要经过这个if
+    # 判定方式是用差值平方和和位移对比
+    gxpd = 0
+    lenxc = len(xcx)
+    yy = 0
+    while yy < lenxc:
+        sumxcx = xcx[yy] * xcx[yy]
+        sumxcy = xcy[yy] * xcy[yy]
+        sumxc = sumxcx + sumxcy
+        xcpd = math.sqrt(sumxc)
+        if xcpd > wypd:
+            gxpd = -1
+            break
+        yy = yy + 1
+    if gxpd == 0:
+        # 输出之前把转换参数都算出来
+        # X=(z1+1)*x+z2*y+(旧重心坐标x-新重心坐标x-新重心坐标x*z1-新重心坐标y*z2+z3)
+        # Y=(-z2)*x+(z1+1)*y+(z1*新重心坐标y+新重心坐标x*z2+z4+旧重心坐标y-新重心坐标y)
+        n1 = float(arrz[0]) + 1
+        n2 = float(arrz[1])
+        s1 = zxzbnewx * float(arrz[0])
+        s2 = zxzbnewy * float(arrz[1])
+        n3 = zxzbpastx - zxzbnewx - s1 - s2 + float(arrz[2])
+        n4 = (-1) * float(arrz[1])
+        n5 = n1
+        s3 = (-1)*float(arrz[0])*zxzbnewy
+        s4 = float(arrz[1])*zxzbnewx
+        n6 = s3 + s4 + float(arrz[3])+zxzbpasty-zxzbnewy
+        sxylist = []
+        sxylist.append(n1)
+        sxylist.append(n2)
+        sxylist.append(n3)
+        sxylist.append(n4)
+        sxylist.append(n5)
+        sxylist.append(n6)
+        # 输出
+        relist1 = []
+        relist1.append(listname)
+        relist1.append(xcx)
+        relist1.append(xcy)
+        relist1.append(listcgcs)
+        relist1.append(gszbx)
+        relist1.append(gszby)
+        relist1.append(sxylist)
+        print(" Finish!")
+        # 输出的那个不是改算坐标，而是有计算的改算坐标
+        newgsx = []
+        newgsy = []
+        lengs1 = len(gszbx)
+        nlist = relist1[6]
+        # return relist
+        #中文再来一次
+        utf_pan = to_utf8(pastname)
+        utf_nn = to_utf8(newname)
+        utf_tn = to_utf8(excelname)
+        insert_into_database1(dbpath,utf_tn,listname1,listname,gsx,gsy,listpastx1,listpasty1,listnewx1,listnewy1,nlist,utf_pan,utf_nn,wypd,sfxs,xzj,pycsx,pycsy,mxjdx,mxjdy,zxzbpastx,zxzbpasty,zxzbnewx,zxzbnewy,n1,n2,n3,n4,n5,n6)
+
+    else:           
+        # 先把所有的合在一起，方便删除
+        lenlist1 = len(xcx)
+        ii = 0
+        relist1 = []
+        while ii < lenlist1:
+            smlist1 = []
+            nn2 = xcx[ii] * xcx[ii]
+            mm2 = xcy[ii] * xcy[ii]
+            nm2 = math.sqrt(nn2 + mm2)
+            smlist1.append(nm2)
+            smlist1.append(listname[ii])
+            smlist1.append(listnewx[ii])
+            smlist1.append(listnewy[ii])
+            smlist1.append(listpastx[ii])
+            smlist1.append(listpasty[ii])
+            relist1.append(smlist1)
+            ii = ii + 1
+        # 直接删除最大的那个
+        relist1.sort(key=takeFirst, reverse=True)
+        num1 = relist1[0]
+        # 获取name
+        num2 = num1[1]
+        ii2 = 0
+        mm2 = 0
+        while ii2 < lenlist1:
+            if listname[ii2] == num2:
+                del listname[ii2]
+                del listnewx[ii2]
+                del listnewy[ii2]
+                del listpastx[ii2]
+                del listpasty[ii2]
+                points = points - 1
+                break
+            else:
+                ii2 = ii2 + 1
+        while mm2 < len(listcgcs):
+            if listname1[mm2] == num2:
+                listcgcs[mm2] = -1
+                break
+            else:
+                mm2 = mm2 + 1
+        print(" Finish!")
+        js1 = 1
+        bhjs(dbpath,listcgcs1, listname1, listpastx1, listpasty1, wypd, listname, listnewx, listnewy, listpastx,
+                listpasty, points, listcgcs, sf, newname, pastname, js1,listnewx1,listnewy1,excelname)
+
+
+# 主函数 计算
+def main_function(dbpath,file_name):
+    #从数据库中调用
+    listcgcs1, listname1, listpastx1, listpasty1, wypd0, listname, listnewx, listnewy, listpastx,listpasty, points, listcgcs,sf, newname, pastname, js,listnewx1,listnewy1 = tablein(dbpath,file_name)
+    #计算
+    bhjs(dbpath,listcgcs1, listname1, listpastx1, listpasty1, wypd0, listname, listnewx, listnewy, listpastx,listpasty, points, listcgcs,sf, newname, pastname, js,listnewx1,listnewy1,file_name)
+    messagebox.showinfo("提示",f"文件 '{file_name}' 计算成功！")

Back/WD/WD.py

@@ -1,2 +1,344 @@
-def main_function(file_path, db_path):
-    print("这是WD")
+import openpyxl
+import xlrd
+import sqlite3
+import os
+from openpyxl import Workbook
+import tkinter as tk
+from tkinter import messagebox
+import math
+import numpy as np
+
+def to_utf8(text):
+    str1 = text.encode('utf-8')
+    return str1
+
+# 弹窗提示用户是否更新数据
+def ask_for_update(file_name):
+    root = tk.Tk()
+    root.withdraw()  # 隐藏主窗口
+    response = messagebox.askyesno("提示", f"检测到数据库中存在相同文件名 '{file_name}'，是否更新数据？")
+    return response
+
+# 读取Excel文件并计算公式结果
+def load_and_calculate_excel(file_path,cell0):
+    # 加载Excel文件并计算单元格的公式结果。
+    workbook = openpyxl.load_workbook(file_path, data_only=True)
+    sheet = workbook.active
+    h1_value = sheet[cell0].value
+    return h1_value
+
+# 数据库插入函数
+def insert_into_database(database,pastname,newname,beforename,excelname,listname1,listpastx1,listpasty1,listcgcs1,pjbc,fxzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,listbex,listbey,listnewx,listnewy):
+    #先把输入的录入数据库
+    db = sqlite3.connect(database)
+    #获取游标
+    cursor = db.cursor()
+    #字符类的都需要转换成字节存进去
+    utf_pan = to_utf8(pastname)
+    utf_nn = to_utf8(newname)
+    utf_bn = to_utf8(beforename)
+    utf_tn = to_utf8(excelname)
+    try:
+        # 查询是否已存在相同的记录
+        cursor.execute(
+            "SELECT * FROM WD_Input_Param WHERE TableName = ?",
+            (utf_tn,)
+        )
+        existing_record = cursor.fetchone()
+
+        if existing_record:
+            # 弹窗询问用户是否更新数据
+            if not ask_for_update(excelname):
+                print("用户选择不更新数据")
+                return
+
+            # 更新现有记录WD_Input_Param 
+            cursor.execute(
+                """
+                UPDATE WD_Input_Param 
+                SET First_ResultName = ?,Last_ResultName=?,New_ResultName=?,Avg_SL=?,Ms_Dir=?,Pt_Count=?,SL_Count=?,Dir_Count=?,Scale_Value=?,Avg_MSL=?,Avg_Dir=?
+                WHERE TableName = ?
+                """,
+                (utf_bn,utf_pan,utf_nn,pjbc,fxzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,utf_tn,)
+            )
+            # 更新现有记录WD_Input_Point 
+            #删了已有的数据，重新插入
+            cursor.execute(
+                """
+                DELETE FROM WD_Input_Point WHERE TableName = ?
+                """,
+                (utf_tn,)
+            )
+            rr = 0
+            while rr < len(listname1):
+                rr1 = listname1[rr]
+                rr2 = listpastx1[rr]
+                rr3 = listpasty1[rr]
+                rr4 = listcgcs1[rr]
+                rr5 = listnewx[rr]
+                rr6 = listnewy[rr]
+                rr7 = listbex[rr]
+                rr8 = listbey[rr]
+                utf_rr1 = to_utf8(rr1)
+                cursor.execute(
+                    """
+                    INSERT INTO WD_Input_Point(TableName,First_ResultName,Last_ResultName,New_ResultName,Last_X,Last_Y,New_X,New_Y,PointName,Wight,First_X,First_Y) VALUES (?,?,?,?,?,?,?,?,?,?,?,?)
+                    """,
+                    (utf_tn,utf_bn,utf_pan,utf_nn,rr2,rr3,rr5,rr6,utf_rr1,rr4,rr7,rr8,)
+                )
+                rr = rr + 1
+            messagebox.showinfo("提示",
+                                f"文件 '{excelname}' 已成功更新到本地数据库中，点击确认以关闭此对话框。对话框关闭后，请点击“计算”以重新计算数据")
+        else:
+            # 插入新记录WD_Input_Param 
+            cursor.execute(
+                """
+                INSERT INTO WD_Input_Param (TableName,First_ResultName,Last_ResultName,New_ResultName,Avg_SL,Ms_Dir,Pt_Count,SL_Count,Dir_Count,Scale_Value,Avg_MSL,Avg_Dir) 
+                VALUES (?,?,?,?,?,?,?,?,?,?,?,?)
+                """,
+                (utf_tn,utf_bn,utf_pan,utf_nn,pjbc,fxzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,)
+            )
+            # 插入新记录WD_Input_Point 
+            rr = 0
+            while rr < len(listname1):
+                rr1 = listname1[rr]
+                rr2 = listpastx1[rr]
+                rr3 = listpasty1[rr]
+                rr4 = listcgcs1[rr]
+                rr5 = listnewx[rr]
+                rr6 = listnewy[rr]
+                rr7 = listbex[rr]
+                rr8 = listbey[rr]
+                utf_rr1 = to_utf8(rr1)
+                cursor.execute(
+                    """
+                    INSERT INTO WD_Input_Point(TableName,First_ResultName,Last_ResultName,New_ResultName,Last_X,Last_Y,New_X,New_Y,PointName,Wight,First_X,First_Y) VALUES (?,?,?,?,?,?,?,?,?,?,?,?)
+                    """,
+                    (utf_tn,utf_bn,utf_pan,utf_nn,rr2,rr3,rr5,rr6,utf_rr1,rr4,rr7,rr8,)
+                )
+                rr = rr + 1
+            messagebox.showinfo("提示",
+                                f"文件 '{excelname}' 已成功添加到本地数据库中，点击确认以关闭此对话框。对话框关闭后，请点击“计算”以计算数据")
+
+        db.commit()
+    except sqlite3.Error as e:
+        print(f"插入或更新文件名时发生错误: {e}")
+    except Exception as e:
+        print(f"处理文件时发生错误: {e}")
+    finally:
+        db.close()
+
+#读取xls
+def xlrd_read(excelpath,dbpath):
+    global gszbx
+    global gszby
+    listname1 = []
+    listpastx1 = []
+    listpasty1 = []
+    listcgcs1 = []
+    #excel表名就是表名
+    excelname = os.path.basename(excelpath)
+    # 读取excel
+    xlr = xlrd.open_workbook(excelpath)
+    tabler = xlr.sheets()[0]
+    # 从第3行开始，3A为点号，3B为最新x，3C为最新y，3D为往期x
+    # 3E为往期y，3F为成果次数，3G为最早x，3H最早y
+    # 先看点数
+    rows = tabler.nrows
+    row = 2
+    points = rows - 2
+    # 收集每一列的数据，方便后期计算
+    listname = []
+    listnewx = []
+    listnewy = []
+    listpastx = []
+    listpasty = []
+    # 成果次数
+    listcgcs = []
+    # 最早的
+    listbex = []
+    listbey = []
+    sumnewx = 0
+    sumnewy = 0
+    sumpastx = 0
+    sumpasty = 0
+    # 记录下最早的那几个开头
+    # excelpath是输入表格，outpath是输出路径，pjbc是平均边长
+    # fxzwc是方向中误差，pjfxs是平均方向数，sf是缩放值
+    # 先计算位移判断值
+    pjbc = float(tabler.cell(0, 9).value)
+    fxzwc = float(tabler.cell(1, 9).value)
+    zbs = float(tabler.cell(2, 9).value)
+    zfxs = float(tabler.cell(3, 9).value)
+    pjbcs = zbs / points
+    pjfxs = zfxs / points
+    sf = float(tabler.cell(4, 9).value)
+    # 再拆细点
+    wy1 = pjbc / 1000
+    wy2 = wy1 + 1
+    wypd1 = wy2 * wy2 / pjbcs
+    # 好像206.265是常数系数
+    wy3 = pjbc * fxzwc / 206.265
+    wypd2 = wy3 * wy3 / pjfxs
+    wypd3 = math.sqrt(wypd1 + wypd2)
+    wypd = round((wypd3 * 2 * math.sqrt(2)), 9)
+    wypd0 = wypd
+    lastname = tabler.cell(0, 3).value
+    finalname = tabler.cell(0, 1).value
+    befname = tabler.cell(0, 6).value
+    utf_ln = to_utf8(lastname)
+    utf_fn = to_utf8(finalname)
+    utf_bn = to_utf8(befname)
+
+    while row < rows:
+        pname = tabler.cell(row, 0).value
+        pnewx = float(tabler.cell(row, 1).value)
+        sumnewx = sumnewx + pnewx
+        pnewy = float(tabler.cell(row, 2).value)
+        sumnewy = sumnewy + pnewy
+        ppastx = float(tabler.cell(row, 3).value)
+        sumpastx = sumpastx + ppastx
+        ppasty = float(tabler.cell(row, 4).value)
+        sumpasty = sumpasty + ppasty
+        pcgcs = float(tabler.cell(row, 5).value)
+        pbex = float(tabler.cell(row, 6).value)
+        pbey = float(tabler.cell(row, 7).value)
+        listname.append(pname)
+        listnewx.append(pnewx)
+        listnewy.append(pnewy)
+        listpastx.append(ppastx)
+        listpasty.append(ppasty)
+        listcgcs.append(pcgcs)
+        listbex.append(pbex)
+        listbey.append(pbey)
+        row = row + 1
+    # 最后返回这里输出
+    rr = 0
+    while rr < len(listname):
+        rr1 = listname[rr]
+        rr2 = listpastx[rr]
+        rr3 = listpasty[rr]
+        rr4 = listcgcs[rr]
+        listname1.append(rr1)
+        listpastx1.append(rr2)
+        listpasty1.append(rr3)
+        listcgcs1.append(rr4)
+        rr = rr + 1
+    #存入数据库
+    insert_into_database(dbpath,lastname,finalname,befname,excelname,listname1,listpastx1,listpasty1,listcgcs1,pjbc,fxzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,listbex,listbey,listnewx,listnewy)
+
+#读取xlsx
+def openpyxl_read(excelpath,dbpath):
+    global gszbx
+    global gszby
+    listname1 = []
+    listpastx1 = []
+    listpasty1 = []
+    listcgcs1 = []
+
+    #excel表名就是表名
+    excelname = os.path.basename(excelpath)
+    # 读取excel
+    xlr = openpyxl.load_workbook(excelpath)
+    tabler = xlr.worksheets[0]
+    # 从第3行开始，3A为点号，3B为最新x，3C为最新y，3D为往期x
+    # 3E为往期y
+    # 先看点数
+    rows = tabler.max_row
+    row = 3
+    points = rows - 2
+    # 收集每一列的数据，方便后期计算
+    listname = []
+    listnewx = []
+    listnewy = []
+    listpastx = []
+    listpasty = []
+    listcgcs = []
+    listbex=[]
+    listbey=[]
+    sumnewx = 0
+    sumnewy = 0
+    sumpastx = 0
+    sumpasty = 0
+    # 记录下最早的那几个开头
+    # excelpath是输入表格，outpath是输出路径，pjbc是平均边长
+    # fxzwc是方向中误差，zrbzwc是最弱边边长相对中误差，pjbcs是平均边数，pjfxs是平均方向数，sf是缩放值
+    # 新增总边数zbs，总方向数zfxs
+    # 先计算位移判断值
+    pjbc = float(tabler['J1'].value)
+    fxzwc = float(tabler['J2'].value)
+    try:
+        zbs = float(tabler['J3'].value)
+    except:
+        zbs = load_and_calculate_excel(excelpath,'J3')
+    try:
+        zfxs = float(tabler['J4'].value)
+    except:
+        zfxs = load_and_calculate_excel(excelpath,'J4')
+    pjbcs = zbs / points
+    pjfxs = zfxs / points
+    sf = float(tabler['J5'].value)
+    # 再拆细点
+    wy1 = pjbc / 1000
+    wy2 = wy1 + 1
+    wypd1 = wy2 * wy2 / pjbcs
+    # 好像206.265是常数系数
+    wy3 = pjbc * fxzwc / 206.265
+    wypd2 = wy3 * wy3 / pjfxs
+    wypd3 = math.sqrt(wypd1 + wypd2)
+    wypd = round((wypd3 * 2 * math.sqrt(2)), 9)
+    wypd0 = wypd
+    befname=tabler['G1'].value
+    lastname = tabler['D1'].value
+    finalname = tabler['B1'].value
+    while row <= rows:
+        try:
+            pname = tabler.cell(row, 1).value
+            pnewx = float(tabler.cell(row, 2).value)
+            sumnewx = sumnewx + pnewx
+            pnewy = float(tabler.cell(row, 3).value)
+            sumnewy = sumnewy + pnewy
+            ppastx = float(tabler.cell(row, 4).value)
+            sumpastx = sumpastx + ppastx
+            ppasty = float(tabler.cell(row, 5).value)
+            sumpasty = sumpasty + ppasty
+            pcgcs = float(tabler.cell(row, 6).value)
+            pbex = float(tabler.cell(row, 7).value)
+            pbey = float(tabler.cell(row, 8).value)
+            listname.append(pname)
+            listnewx.append(pnewx)
+            listnewy.append(pnewy)
+            listpastx.append(ppastx)
+            listpasty.append(ppasty)
+            listcgcs.append(pcgcs)
+            listbex.append(pbex)
+            listbey.append(pbey)
+            row = row + 1
+        except:
+            break
+    # 最后返回这里输出
+    rr = 0
+    while rr < len(listname):
+        rr1 = listname[rr]
+        rr2 = listpastx[rr]
+        rr3 = listpasty[rr]
+        rr4 = listcgcs[rr]
+        listname1.append(rr1)
+        listpastx1.append(rr2)
+        listpasty1.append(rr3)
+        listcgcs1.append(rr4)
+        rr = rr + 1
+    #存入数据库
+    insert_into_database(dbpath,lastname,finalname,befname,excelname,listname1,listpastx1,listpasty1,listcgcs1,pjbc,fxzwc,points,zbs,zfxs,sf,pjbcs,pjfxs,listbex,listbey,listnewx,listnewy)
+
+# 主函数 读取excel文件
+def main_function(file_path,dbpath):
+    # 调用读取Excel文件的函数
+    file_name = os.path.basename(file_path)
+    #两种加载方式，对于不同的文件
+    if ".xlsx" in file_path:
+        #采用openpyxl
+        openpyxl_read(file_path,dbpath)
+    else:
+        #采用xlrd
+        xlrd_read(file_path,dbpath)

Back/WD/WDcompute.py

@@ -1,2 +0,0 @@
-def main_function(file_path, db_path):
-    pass  # 添加计算的方法