python 在测绘作业中的一些小应用（与cad交互）-1

虽然笔者已经基本上告别了本科的测绘工程专业，但是笔者的本科同学他们在实际作业中难免会遇到一些批量化、重复性劳动问题，如果会编程，写上一个小脚本，无疑会提高工作效率。

下面是笔者本科同学处理测量数据时遇到的问题，笔者通过python与CAD的交互完美的解决了！

问题1：提取断面上地性线上的高程

问题2：场平设计高程，十字丝交点的坐标（x,y,z)

首先，python 作为一个胶水语言，与各行各业的接口是无与伦比，与cad的接口库有：autocad、dxfgrabber、dxfwrite。笔者目前就知道这三个，孤陋寡闻。如有兴趣，可以关注freecad，一个开源的python 与cad软件。

笔者程序使用的是dxfgrabber库，读取dxf文件。

对于第一个问题：笔者写了一个GUI

from PyQt5 import  QtWidgets
from cass import Ui_MainWindow
from PyQt5.QtWidgets import QFileDialog
import dxfgrabber
class height(QtWidgets.QMainWindow,Ui_MainWindow):
def __init__(self):
super(height,self).__init__()
self.setupUi(self)
self.openfile.triggered.connect(self.o1)
self.extract_h.triggered.connect(self.e1)
self.save_txt.triggered.connect(self.s1)
def o1(self):
self.file,ok=QFileDialog.getOpenFileName(self,"打开",None,"*.dxf")
def e1(self):
dxf=dxfgrabber.readfile(self.file)
h1=[]
c=[]
allGCD=[h for h in dxf.entities if h.layer=='GCD']
alldmx=[line for line in dxf.entities if line.layer=='dmx']
l1=[line for line in alldmx if line.dxftype=='LINE']
for heigh in allGCD:
h1.append(heigh.insert)
for i in l1:
c.append([i.start,i.end])
#计算每个线的斜率
k1=[]
for i in range(len(c)):
k1.append((c[i][1][1]-c[i][0][1])/(c[i][1][0]-c[i][0]0]))
#计算高程点与直线的斜率,#，以所有线的起点为基准，判断每个高程点和起点连线的斜率与所有线的斜率，如果相同，则高程点再地性线上
k2=[]
for i in range(len(c)):
for j in range(len(h1)):
k3=(h1[j][1]-c[i][0][1])/(h1[j][0]-c[i][0][0])
if abs(k3)-abs(k1[i])>-0.2 and abs(k3)-abs(k1[i])<0.2:
k2.append(j)
self.allresult_h=[]
for i in k2:
self.allresult_h.append(h1[i])
def s1(self):
path,ok=QFileDialog.getSaveFileName(self,"保存结果",None,"*.txt")
with open(path,'w') as file:
for line in self.allresult_h:
for i in line:
file.writelines(str(i)+',')
file.writelines('\n')
if __name__=='__main__':
import sys
app=QtWidgets.QApplication(sys.argv)
myapp=height()
myapp.show()
app.exec_()

“`