python基于北京公交线txt数据画cad

最新推荐文章于 2024-05-29 14:52:57 发布
最新推荐文章于 2024-05-29 14:52:57 发布
阅读量648 收藏 12
点赞数
文章标签： python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43073222/article/details/106795739
版权

这里写自定义目录标题

运用python基于北京公交线网txt数据绘制cad

文章前半部分为学习过程，后半部分为在cad中对北京市两千多条公交线路（单向）的绘制，先上效果图2790个图层的北京公交线网

导入pyautocad

pip install pyautocad

读取cad文件

from pyautocad import Autocad,APoint
import math

acad = Autocad(create_if_not_exists=True)  #必须把cad文件打开

在cad控制台上打字

acad.prompt('hello world')
print(acad.doc.Name)  #输出文件名

读取图层信息

layers_nums=acad.ActiveDocument.Layers.count  #图层总数
print(layers_nums)
layers_names=[acad.ActiveDocument.Layers.Item(i).Name for i in range(layers_nums)]  #图层名称
print(layers_names)

读取多段线顶点坐标

items=[]
for item in acad.iter_objects("PolyLine"):
    print(item.Coordinates)
    items.append(item.Coordinates)
print(len(items))

新建图层置为当前

layerobj=acad.ActiveDocument.Layers.Add('HIT_Layer')
acad.ActiveDocument.ActiveLayer=layerobj

cad上画图

#设置颜色
ClrNum=1
layerobj.color=ClrNum
# ClrNum为颜色索引号，其取值范围为[0,256]；标准颜色的颜色索引号指定如下：：1 红、2 黄、3 绿、4 青、5 蓝、6 洋红、7 白/黑； 0 ByBlock、256 ByLayer；

 #画线l
p1=APoint(0,0)   #所有点必须这样表示
p2=APoint(100,100)
acad.model.AddLine(p1,p2)

#画多段线pl
lon=[[55,66,77],[1,3,44],[755,74,1,20]]
lat=[[55,66,77],[1,3,44],[755,74,1,20]]
pnts=[]
for i in range(len(lon[2])):
    pnts.append(APoint(lon[2][i],lat[2][i]))
pnts = np.array([j for i in pnts for j in i], dtype=np.float)
pnts = Autocad.aDouble(pnts)
pline_obj = acad.model.AddPolyLine(pnts)

#画图形
for i in range(5):   
    text = acad.model.AddText('Hi %s!' % i, p1, 2.5)
    acad.model.AddLine(p1, p2)
    acad.model.AddCircle(p1, 10)
    p1.y += 10
dp = APoint(10, 0)

保存

acad.ActiveDocument.Application.preferences.OpenSave.AutoSaveInterval = 0.51

完整代码

from pyautocad import Autocad,APoint
import math

acad = Autocad(create_if_not_exists=True)       #必须把cad文件打开
# acad.ActiveDocument.Application.Documents.Open("2.dwg")

#在cad控制台上打字
acad.prompt('hello world')
print(acad.doc.Name)  #输出文件名

# # 使用SummaryInfo对象访问图纸的属性（只适用于Windows），也可以使用SummaryInfo访问自动定义属性
# for doc in acad.app.Documents:
#     print(doc.Name)  # 文件名
#     print(doc.SummaryInfo.Author)
#     print(doc.SummaryInfo.Comments)
#     print(doc.SummaryInfo.HyperlinkBase)
#     print(doc.SummaryInfo.Keywords)
#     print(doc.SummaryInfo.LastSavedBy)
#     print(doc.SummaryInfo.RevisionNumber)
#     print(doc.SummaryInfo.Subject)
#     print(doc.SummaryInfo.Title)

# 读取图层信息
layers_nums=acad.ActiveDocument.Layers.count  #图层总数
print(layers_nums)
layers_names=[acad.ActiveDocument.Layers.Item(i).Name for i in range(layers_nums)]  #图层名称
print(layers_names)
print(acad.ActiveDocument.Layers.Item(0))


# 读取多段线顶点坐标
items=[]
for item in acad.iter_objects("PolyLine"):
    print(item.Coordinates)
    items.append(item.Coordinates)
print(len(items))


#新建图层置为当前
layerobj=acad.ActiveDocument.Layers.Add('HIT_Layer')
acad.ActiveDocument.ActiveLayer=layerobj
#设置颜色
ClrNum=1
layerobj.color=ClrNum
# ClrNum为颜色索引号，其取值范围为[0,256]；
 # 标准颜色的颜色索引号指定如下：：1 红、2 黄、3 绿、4 青、5 蓝、6 洋红、7 白/黑；
 # 0 ByBlock、256 ByLayer；


#画线
p1=APoint(0,0)   #所有点必须这样表示
p2=APoint(100,100)
acad.model.AddLine(p1,p2)


#画多段线
lon=[[55,66,77],[1,3,44],[755,74,1,20]]
lat=[[55,66,77],[1,3,44],[755,74,1,20]]
pnts=[]
for i in range(len(lon[2])):
    pnts.append(APoint(lon[2][i],lat[2][i]))
pnts = np.array([j for i in pnts for j in i], dtype=np.float)
pnts = Autocad.aDouble(pnts)

          # 将各点坐标顺序变换为行数据；
          # 添加样条曲线时参数仅支持1行多列的1维数组，即将各点x,y,z坐标顺序排列构成的数组。
          # 如报错可尝试增加该行代码 pnts = aDouble(pnts)
pline_obj = acad.model.AddPolyLine(pnts)


for i in range(5):   #画图形
    text = acad.model.AddText('Hi %s!' % i, p1, 2.5)
    acad.model.AddLine(p1, p2)
    acad.model.AddCircle(p1, 10)
    p1.y += 10
dp = APoint(10, 0)

# 更改标注的位置
for text in acad.iter_objects('Text'):
    print('text: %s at: %s' % (text.TextString, text.InsertionPoint))
    text.InsertionPoint = APoint(text.InsertionPoint) + dp

for obj in acad.iter_objects(['Circle', 'Line']):
    print(obj.ObjectName)

# text = acad.model.AddText("{0}".format(text), p1, 15)             #添加文本，文本字符串，点的位置下显示文本，字体大小
# text.move(p1, p2)                                                 #文本從p1點移動到p2點
# acad.model.AddCircle(p1, 10)                                      #画圆，圆心，半径
# acad.model.AddArc(p1, 500, math.radians(90), math.radians(270))   #画圆弧，圆心，半径，开始弧度，结束弧度

# 保存cad图到指定位置，第一個參數一定是保存的絕對路徑！！，64是保存的一種格式，筆者從1到100都試了一遍，
# 有一些是可以生成图的，都是CAD可以打開的文件類型格式。其實跟着筆者寫64就好！
# acad.doc.SaveAs("D:\\pycharm\\mycode\\this_pc_new\\dwg_read\\test.dwg", 64)
# acad.doc.SaveAs("{0}".format("D:\\pycharm\\mycode\\this_pc_new\\dwg_read\\test.dwg"), 64)
acad.ActiveDocument.Application.preferences.OpenSave.AutoSaveInterval = 0.51  #自动保存

将txt线条打印到cad中

遇到一个需求，统计北京市的所有公交线并置于cad中，目前有txt的文件，大概55.4m，这里放个字段的截图。
北京公交部分线路信息
数据很大，不可能人工画，根据线路代码进行绘制，读取到一个新的线路代码就新建一个cad图层，直接上代码

from pyautocad import Autocad,APoint
import math
from tqdm import tqdm
from numba import jit
import numpy as np

#读取公交线路信息
s = []
code=[]
lon=[]
lat=[]
f= open('bus.txt','r')
for lines in f:
    ls = lines.strip('\n').replace(' ','').replace('、','/').replace('?','').split('/')
    for i in ls:
        s.append(i)
f.close()
del s[0]
for i in s:
    b = i.split(',', 2)
    if b[0] not in code:
        code.append(b[0])
        lon.append([])
        lat.append([])
    lon[-1].append(b[1])
    lat[-1].append(b[2])


# 高德火星坐标系转化为WGS84地理坐标系
def gcj2wgs(jingdu,weidu):
    lon = float(jingdu)
    lat = float(weidu)
    a = 6378245.0  # 克拉索夫斯基椭球参数长半轴a
    ee = 0.00669342162296594323  # 克拉索夫斯基椭球参数第一偏心率平方
    PI = 3.14159265358979324  # 圆周率
    # 以下为转换公式
    x = lon - 105.0
    y = lat - 35.0
    # 经度
    dLon = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * math.sqrt(abs(x))
    dLon += (20.0 * math.sin(6.0 * x * PI) + 20.0 * math.sin(2.0 * x * PI)) * 2.0 / 3.0
    dLon += (20.0 * math.sin(x * PI) + 40.0 * math.sin(x / 3.0 * PI)) * 2.0 / 3.0
    dLon += (150.0 * math.sin(x / 12.0 * PI) + 300.0 * math.sin(x / 30.0 * PI)) * 2.0 / 3.0
    # 纬度
    dLat = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * math.sqrt(abs(x))
    dLat += (20.0 * math.sin(6.0 * x * PI) + 20.0 * math.sin(2.0 * x * PI)) * 2.0 / 3.0
    dLat += (20.0 * math.sin(y * PI) + 40.0 * math.sin(y / 3.0 * PI)) * 2.0 / 3.0
    dLat += (160.0 * math.sin(y / 12.0 * PI) + 320 * math.sin(y * PI / 30.0)) * 2.0 / 3.0
    #转换
    radLat = lat / 180.0 * PI
    magic = math.sin(radLat)
    magic = 1 - ee * magic * magic
    sqrtMagic = math.sqrt(magic)
    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * PI)
    dLon = (dLon * 180.0) / (a / sqrtMagic * math.cos(radLat) * PI)
    wgsLon = lon - dLon
    wgsLat = lat - dLat
    return wgsLon, wgsLat


# 度数坐标系转换成小数坐标系
def deg2flo(jingdu,weidu):
    lon=str(jingdu)
    lat=str(weidu)
    flon=float(lon[:3])+float(lon[3:])/60
    flat = float(lat[:2]) + float(lat[2:])/60
    return flon,flat


# WGS84国际坐标系转GCJ02(火星坐标系)
def wgs2gcj(jingdu, weidu):
    lon = float(jingdu)
    lat = float(weidu)
    a = 6378245.0  # 克拉索夫斯基椭球参数长半轴a
    ee = 0.00669342162296594323  # 克拉索夫斯基椭球参数第一偏心率平方
    PI = 3.14159265358979324  # 圆周率
    # 以下为转换公式
    x = lon - 105.0
    y = lat - 35.0
    # 经度
    dLon = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * math.sqrt(abs(x))
    dLon += (20.0 * math.sin(6.0 * x * PI) + 20.0 * math.sin(2.0 * x * PI)) * 2.0 / 3.0
    dLon += (20.0 * math.sin(x * PI) + 40.0 * math.sin(x / 3.0 * PI)) * 2.0 / 3.0
    dLon += (150.0 * math.sin(x / 12.0 * PI) + 300.0 * math.sin(x / 30.0 * PI)) * 2.0 / 3.0
    # 纬度
    dLat = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * math.sqrt(abs(x))
    dLat += (20.0 * math.sin(6.0 * x * PI) + 20.0 * math.sin(2.0 * x * PI)) * 2.0 / 3.0
    dLat += (20.0 * math.sin(y * PI) + 40.0 * math.sin(y / 3.0 * PI)) * 2.0 / 3.0
    dLat += (160.0 * math.sin(y / 12.0 * PI) + 320 * math.sin(y * PI / 30.0)) * 2.0 / 3.0
    # 转换
    radLat = lat / 180.0 * PI
    magic = math.sin(radLat)
    magic = 1 - ee * magic * magic
    sqrtMagic = math.sqrt(magic)
    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * PI)
    dLon = (dLon * 180.0) / (a / sqrtMagic * math.cos(radLat) * PI)
    wgsLon = lon + dLon
    wgsLat = lat + dLat
    return wgsLon, wgsLat



#画线
acad = Autocad(create_if_not_exists=True)       #必须把cad文件打开
acad.ActiveDocument.Application.Documents.Open("aa.dwg")
for i in tqdm(range(len(code))):     #公交线数量
    layerobj = acad.ActiveDocument.Layers.Add(str(code[i]))  # 新建图层置为当前   PolyLine
    acad.ActiveDocument.ActiveLayer = layerobj
    pnts = []      #多段线坐标转换
    for j in range(len(lon[i])):
        [wgsLon,wgsLat]=wgs2gcj(lon[i][j],lat[i][j])  #坐标转换
        pnts.append(APoint(float(wgsLon)/5.8705*499.29, float(wgsLat)/6.3795*702.64))
        # pnts.append(APoint(float(wgsLon), float(wgsLat)))
    pnts = np.array([m for n in pnts for m in n], dtype=np.float)
    pnts = Autocad.aDouble(pnts)
    pline_obj = acad.model.AddPolyLine(pnts)

acad.ActiveDocument.Application.preferences.OpenSave.AutoSaveInterval = 0.51  #自动保存
weixin_43073222
关注
  • 0
    点赞
  • 踩
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
python cad
12-28
使用Python写的CAD代码，使用OCC库，PyQt 界面，可跨平台使用。
python制作会动的北京市公交线路图（附数据）.rar
09-21
今天教大家用pyecharts制作北京市公交线路动态图今天教大家用pyecharts制作北京市公交线路动态图今天教大家用pyecharts制作北京市公交线路动态图
Python解析.dwg格式文件信息提取
懒惰是科技进步的原始动力
05-29 940
dwg文件提取文本信息，方便通过某零件型号快速定位搜索图纸
Python数据可视化：基于Plotly的城市公交路线地图绘制
NayelyA的博客
08-31 7439
需要达到的目标有爬取相关网站，获得某市的公交车线路数据并对这些数据实现解析，获取每一条线路公交车的具体车站坐标显示在地图上，具体的效果以延安市为例进行展示，如下。 ![在这里插入图片描述](https://img-blog.csdnimg.cn/2019083118251292.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,s...
python3读取autocad图形文件.py实例
09-16
主要介绍了python3读取autocad图形文件.py实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
python读取cad表格_使用Python读取AutoCAD DXF文档
weixin_39588265的博客
12-04 1175
作为《使用Python创建AutoCAD DXF文档》的姊妹篇，本文介绍用Python读取DXF========================================================================================================================================================...
基于python的运动员数据分析源码.zip
最新发布
05-31
基于python的运动员数据分析源码.zip基于python的运动员数据分析源码.zip基于python的运动员数据分析源码.zip基于python的运动员数据分析源码.zip基于python的运动员数据分析源码.zip基于python的运动员数据分析源码...
基于python的NBA球员数据可视化分析
04-07
在本项目"基于python的NBA球员数据可视化分析"中，我们主要探讨了如何使用Python这一强大且广泛使用的编程语言来获取、处理和可视化NBA球员的数据。该项目分为几个关键阶段，包括数据爬取、数据预处理、数据分析以及...
基于python的天气预报系统设计和可视化数据分析
06-18
python或jutyper 都可成功运行，编写的基于爬虫的和TK界面的天气预报系统，可选择多个城市，15天天气，数据绘图处理以及数据保存，代码python和jutyper均可成功运行 (内含py文件和必要文件)
基于Python实现汽车销售数据可视化【500010086】
11-14
详情介绍：...实现基于Python实现汽车销售数据可视化 1、计算销量同比数据 2、汽车年销售量图 3、2023年与2022年同时期汽车销量数据对比图 4、各厂商汽车销量(年) 5、各车型汽车销量(年)
基于Python实现地震数据可视化.zip
01-06
在本项目"基于Python实现地震数据可视化.zip"中，我们将探讨如何使用Python强大的数据处理和可视化库来展示地震数据。这个压缩包包含了名为"visualization_of_seismic_data-master"的文件夹，里面可能包含了源代码、...
python_cad
10-08
Python Extension for Autocad/Intellicad
用python和pythonOCC写的一个简单的CAD程序
09-09
用python和pythonOCC写的一个简单的CAD程序。相当简陋，许多功能尚未实现。可以画简单的几何体：块体，球，圆锥，圆柱，线，圆弧等。
CAD二次开发的Python脚本
08-28
cad二次开发，Python脚本类 CAD二次开发中，CAD需要采用加载DLL集，并通过CommandMethod[diameter]的方法可以自动生成图纸，但这里的diameter能否设置成变量值，并且这个变量值与窗体上的文本框内容一至。 我试了一种方法，那就是在DLL中设置静态变量，然后在窗体中赋值，失败
广州地铁站台线路图（轨迹图）可视化.py
05-10
代码实现了在python下利用百度api获取地铁站点路线（轨迹）的墨卡托坐标数据并转化为WGS-84坐标系下的经纬度，并使用 plotly+mapbox 绘制地铁路线图
如何用python做软件导出cad_<小白学技术>将python脚本导出为exe可执行程序
weixin_39627052的博客
12-02 204
1.简介(为啥需要导出为exe可执行程序)python写完的程序靠命令来执行，显得太专业，不符合python简单的特点(好吧，主要是太low)代码给别人执行，别人没有你的python库也没法用(双击就能运行，方便大家，方便传播)2.pyInstaller安装配置pip install pyinstaller会自动安装pyInstaller和pywin323.打包打包的app里并不包含任何源码，但将...
如何用python做软件导出cad_用Python编写博客导出工具
weixin_39898150的博客
12-03 207
写在前面的话我在 github 上用 octopress 搭建了个人博客，octopress 使用Markdown语法编写博文。之前我在CSDN博客上也写过不少的技术博文，都说自己的孩子再丑也是个宝，所以就起了把CSDN博客里面的文章导出到个人博客上的念头。刚开始想找个工具把CSDN博客导出为xml或文本，然后再把xml或文本转换为Markdown博文。可惜搜了一下现有博客导出工具，大部分要收费才...
Python与CAD系列基础篇（一）操作AutoCAD与读写CAD文件
热门推荐
weixin_43894929的博客
10-10 1万+
后续基础篇将进一步介绍Python与CAD的基础操作方法，如创建点、线、多边形几何对象；创建文本及图案填充；图形移动与复制；图层相关；标注相关等，通过①通过win32com.client或pyautocad连接AutoCAD进行处理（需要AutoCAD软件，可以操作各种CAD格式文件）②通过ezdxf处理dxf格式文件（不需要AutoCAD软件，格式较单一）两种方式来对基础篇功能实现。高级篇将介绍实例应用，如根据范围输出cad文件中文本、颜色、角度等信息；
python autocad 读取excel数据画多段线
05-18
因此，我们可以使用Python编写脚本，实现读取Excel数据，并在AutoCAD中绘制多段线的功能。 具体步骤如下： 1. 使用COM连接访问AutoCAD程序 2. 打开指定的dwg文件 3. 读取Excel中的数据 4. 根据Excel中的数据，在...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值