Python大疆导出csv文件转化大地2000的dxf文件

大疆导出三维模型里面有个models\pc\0\terra_grid\csv\terra_grid.csv文件，里面记录所有点的坐标和高程，但坐标是经纬度坐标，需要转化为大地2000坐标。

我参照了：经纬度坐标转换为CGCS2000大地坐标系对应XY值（PYTHON实现）_经纬度转换大地2000坐标工具-CSDN博客和百度ai提供的代码：

import ezdxf
doc = ezdxf.new()  # 创建新 DXF 文档
msp = doc.modelspace()  # 获取模型空间# 添加点
msp.add_point((100, 200, 0), dxfattribs={'layer': 'Roads'})# 添加线
points = [(150, 300), (200, 400)]
msp.add_lwpolyline(points, dxfattribs={'layer': 'Buildings'})doc.saveas('output.dxf')

改写了一个能从csv文件转化为dxf文件的程序：

# -*- coding: utf-8 -*-
"""
Created on Wed Mar 19 10:38:48 2025@author: YBK
"""import math
import ezdxf
import csv
doc = ezdxf.new()  # 创建新 DXF 文档
msp = doc.modelspace()  # 获取模型空间def GaussToBLToGauss(longitude, latitude):output = [0.0, 0.0]longitude1, latitude1, longitude0, X0, Y0, xval, yval = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0# ConstantsiPI = 0.0174532925199433  # 3.1415926535898/180.0;a = 6378137.0f = 1 / 298.257222101  # CGCS2000坐标系参数# a=6378137.0; f=1/298.2572236; #wgs84坐标系参数longitude0 = 117.0  # 中央子午线 根据实际进行配置longitude0 = longitude0 * iPI  # 中央子午线转换为弧度longitude1 = longitude * iPI  # 经度转换为弧度latitude1 = latitude * iPI  # 纬度转换为弧度e2 = 2 * f - f * fee = e2 * (1.0 - e2)NN = a / math.sqrt(1.0 - e2 * math.sin(latitude1) * math.sin(latitude1))T = math.tan(latitude1) * math.tan(latitude1)C = ee * math.cos(latitude1) * math.cos(latitude1)A = (longitude1 - longitude0) * math.cos(latitude1)M = a * ((1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256) * latitude1 - (3 * e2 / 8 + 3 * e2 * e2 / 32 + 45 * e2 * e2* e2 / 1024) * math.sin(2 * latitude1)+ (15 * e2 * e2 / 256 + 45 * e2 * e2 * e2 / 1024) * math.sin(4 * latitude1) - (35 * e2 * e2 * e2 / 3072) * math.sin(6 * latitude1))xval = NN * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 72 * C - 58 * ee) * A * A * A * A * A / 120)yval = M + NN * math.tan(latitude1) * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24+ (61 - 58 * T + T * T + 600 * C - 330 * ee) * A * A * A * A * A * A / 720)X0 = 500000Y0 = 0xval = xval + X0yval = yval + Y0# 转换为投影output[0] = xvaloutput[1] = yvalreturn outputtry:with open('terra_grid.csv', 'r', encoding='utf-8') as file:reader = csv.reader(file)for row in reader:# Example usage:longitude_input = float(row[0])latitude_input = float(row[1])hight_imput = float(row[2]) - 14.81798172 # 高程校对 - 14.81798172result = GaussToBLToGauss(longitude_input, latitude_input)msp.add_point((result[0], result[1], hight_imput), dxfattribs={'layer': 'dji_gcd'})print(f"Converted Coordinates: X = {result[0]}, Y = {result[1]}")except FileNotFoundError:print("错误：文件不存在！")
except UnicodeDecodeError:print("错误：文件编码非 UTF-8，请尝试 'encoding='gbk'' 或 'latin-1'")
except Exception as e:print(f"未知错误: {e}")doc.saveas('output0319.dxf')

用时记得修改中央子午线值和高程校对值。