cad（2000坐标系）转kml

拼凑来的，找不到原文章了

import win32com.client
import pythoncom
import simplekml
from tkinter import messagebox,filedialog
import mathdef XY2LatLon(X, Y, L0):iPI = 0.0174532925199433a = 6378137.0f= 0.00335281006247ZoneWide = 3 #按3度带进行投影ProjNo = int(X / 1000000)L0 = L0 * iPIX0 = ProjNo * 1000000 + 500000Y0 = 0xval = X - X0yval = Y - Y0e2 = 2 * f - f * f #第一偏心率平方e1 = (1.0 - math.sqrt(1 - e2)) / (1.0 + math.sqrt(1 - e2))ee = e2 / (1 - e2) #第二偏心率平方M = yvalu = M / (a * (1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256))fai = u + (3 * e1 / 2 - 27 * e1 * e1 * e1 / 32) * math.sin(2 * u) + (21 * e1 * e1 / 16 - 55 * e1 * e1 * e1 * e1 / 32) * math.sin(4 * u) + (151 * e1 * e1 * e1 / 96) * math.sin(6 * u)+ (1097 * e1 * e1 * e1 * e1 / 512) * math.sin(8 * u)C = ee * math.cos(fai) * math.cos(fai)T = math.tan(fai) * math.tan(fai)NN = a / math.sqrt(1.0 - e2 * math.sin(fai) * math.sin(fai))R = a * (1 - e2) / math.sqrt((1 - e2 * math.sin(fai) * math.sin(fai)) * (1 - e2 * math.sin(fai) * math.sin(fai)) * (1 - e2 * math.sin(fai) * math.sin(fai)))D = xval / NN#计算经纬度（弧度单位的经纬度）longitude1 = L0 + (D - (1 + 2 * T + C) * D * D * D / 6 + (5 - 2 * C + 28 * T - 3 * C * C + 8 * ee + 24 * T * T) * D * D * D * D * D / 120) / math.cos(fai)latitude1 = fai - (NN * math.tan(fai) / R) * (D * D / 2 - (5 + 3 * T + 10 * C - 4 * C * C - 9 * ee) * D * D * D * D / 24 + (61 + 90 * T + 298 * C + 45 * T * T - 256 * ee - 3 * C * C) * D * D * D * D * D * D / 720)#换换为deglongitude = longitude1 / iPIlatitude = latitude1 / iPIreturn  longitude,latitude
acad = win32com.client.Dispatch("AutoCAD.Application")
doc = acad.ActiveDocument
doc.Utility.Prompt("\n醉后不知天在水\n满船清梦压星河\n")
mp = doc.ModelSpace  # 模型空间
messagebox.showinfo('提示', "请在屏幕拾取图元，以Enter键结束")
try:doc.SelectionSets.Item("SS1").Delete()
except:messagebox.showinfo('警告', "Delete selection failed")
slt = doc.SelectionSets.Add("SS1")
slt.SelectOnScreen()
kml = simplekml.Kml()for x in slt:if x.ObjectName == "AcDbPolyline":pts = x.Coordinatespoints = [(XY2LatLon(pts[x],pts[x+1],120)) for x in range(0,len(pts),2)]pol = kml.newpolygon(name="czh", outerboundaryis=points)pol.style.polystyle.color = '550000FF' #aabbggrrpol.style.polystyle.outline = 1pnt = kml.newpoint(name="czh", coords=[points[0]])pnt.description = "HPPs in Nepal"pnt.style.labelstyle.color = 'FFFFFF'  # html color code
savepath = filedialog.asksaveasfilename(filetypes=[("kml", "kml")], title="保存kml")+".kml"
kml.save(savepath)

cad（2000坐标系）转kml相关推荐

Global Mapper20生成可以进CAD的国家2000坐标系的DXF等高线
一.加载苏州市12米分辨率地形数据拖拽加载: 二.查看.设置投影信息 2.1查看投影信息双击图层,或右键->[options] 从图新GIS获取的12米分辨率地形数据默认是wgs84,经纬度 ...
快速推进国家2000坐标系的空间数据成果坐标转换解决方案
原文发布时间:2017/6/12 14:36:32 国土资源部.国家测绘局近期发布通知要求各地在2018年7月1日起全面启用国家2000坐标系,停止使用其他坐标系,并明确提出具体任务:1.完成存量数 ...
2000坐标系xy坐标几位_2000坐标(国家2000坐标xy各几位)
2000坐标系,我国当前最新的国家大地坐标系,英文名称为China Geodetic Coordinate System 2000,英文缩写为CGCS2000 2000.80.54这几个坐标系,都是我 ...
arcgis 坐标系 2000坐标系_干货 | 关于投影和坐标系那些事，你是否都已还给了老师？...
很多新接触GIS的人员对地图投影以及坐标系统很难理解,甚至做GIS开发做了好几年的人也有这方面的疑惑,地球仪式的地图是如何变成纸上的平面地图的?平面的二维地图是如何在三维GIS里面进行展示的,因为三维 ...
中添加2000坐标系_ArcGIS API for JavaScript 4.16在三维场景中以天地图为底图加载2000坐标系的倾斜摄影数据...
2021年目标每日更新一篇公众号文章!!! 需求描述 WebGIS项目中如果有三维需求的话,开发中我们经常会遇到倾斜摄影数据,有了倾斜摄影数据之后,我们可以进行生成真实的三维场景等操作. 本文从前端 ...
ENVI：如何自定义beijing54、西安80、国家2000坐标系
这仅仅是一条目录目录这仅仅是一条目录 01 为什么要自定义坐标系? 02 自定义坐标系的思路 03 自定义坐标系实操 (以Beijing1954坐标系.西安80.国家2000坐标系为例) 1. 添 ...
mapbox 支持国家2000 坐标系的数据
最近有个项目使用mapbox作为地图展示,但是mapbox 只支持web 墨卡托(3857)坐标系的数据.而客户想接入其他单位发布的共享数据,共享数据都是基于arcgisserver 发布的国家200 ...
什么是CAD模型坐标系？能够修改吗？
CAD模型坐标系是三维CAD系统中建立的过程中带入的坐标系,在测量软件中导入CAD模型时,程序也同时将模型坐标系导入.当需要在测量软件中对CAD模型移动时,一般需要通过建立CAD模型临时坐标系,将CA ...
【胡学长带你学 Global Mapper 】如何将影像地图转换为国家2000坐标系
多部门联合开展工作坐标数据不匹配怎么办?Hi,我是胡学长!今天教大家如何利用Global Mapper将WGS84影像数据转换为国家2000坐标系影像数据 : (一)EPSG是什么? EPSP的英文全 ...
2000坐标系高程与85高程转换_科普 | 如何在大疆智图中设置坐标系
点击上方蓝字关注我们在大疆智图"输出坐标系设置"中可设置坐标系,例如选择WGS 84/ UTM zone 49N或CGCS2000(2000国家大地坐标系)/Gauss-Kr ...

cad（2000坐标系）转kml

cad（2000坐标系）转kml相关推荐

最新文章

热门文章