测绘距离交会解算

chenqiu521

已知A,B两点坐标及到未知点P的距离，求解P点坐标，最好是批量处理

wp231957

真有刚需，就转战python吧   excel不行的

1. from scipy.optimize import fsolve
   
2. 
   
3. def equations(vars):
   
4.     x, y = vars
   
5.     eq1 = (x - 4391518.06)**2 + (y - 595013.803)**2 - 50.31**2
   
6.     eq2 = (x - 4391536.648)**2 + (y - 595072.09)**2 - 39.79**2
   
7.     return [eq1, eq2]
   
8. 
   
9. # 初始猜想值
   
10. initial_guess = [4391518, 595013]
    
11. 
    
12. # 使用fsolve求解
    
13. solution = fsolve(equations, initial_guess)
    
14. x, y = solution
    
15. 
    
16. print(f'解：x = {x}, y = {y}')
    
17. 
    
18. 
    
19. '''
    
20. PS D:\wp> py -3.8-64 tmp.py
    
21. 解：x = 4391560.747204542, y = 595040.4281530951
    
22. '''

复制代码