使用 ezdxf 库在 DXF 文件中进行 CRS 到 WCS 坐标转换

import ezdxf
from ezdxf import transform
from ezdxf.math import Matrix44

def wcs_to_crs(entities, m: Matrix44):
    """将实体从 WCS 坐标转换为 CRS 坐标。"""
    transform.inplace(entities, m)

def crs_to_wcs(entities, m: Matrix44):
    """将实体从 CRS 坐标转换为 WCS 坐标。"""
    # 转换为 WCS 需要应用逆矩阵
    m_inverse = m.copy()
    m_inverse.inverse()
    transform.inplace(entities, m_inverse)

# 假设您的 DXF 文件名为 "tester.dxf"
doc = ezdxf.readfile("tester.dxf")
msp = doc.modelspace()
geo_data = msp.get_geodata()

# 设置转换方向：True 表示从 CRS 转换为 WCS
CRS_TO_WCS = True 

if geo_data:
    # 从 GEODATA 获取转换矩阵和 EPSG 代码
    m, epsg = geo_data.get_crs_transformation()
    print(f"GEODATA 存在，EPSG: {epsg}")

    if CRS_TO_WCS:
        print("正在将实体从 CRS 转换为 WCS...")
        crs_to_wcs(msp, m)
    else:
        print("正在将实体从 WCS 转换为 CRS...")
        wcs_to_crs(msp, m)
else:
    print("DXF 文件中没有地理参考数据 (GEODATA)。")
    # 如果没有 GEODATA，需要手动提供转换矩阵和 EPSG
    # 这里的 Matrix44() 是一个单位矩阵，表示不进行任何变换
    # epsg = 3395 仅作为标识符，不执行实际投影
    m = Matrix44() 
    epsg = 3395
    print("已使用单位矩阵和默认 EPSG 3395。实体未进行实际坐标转换。")
    # 如果需要实际转换，您必须在这里提供一个有效的转换矩阵
    # 例如：m = Matrix44.scale(factor) 或 Matrix44.translate(x, y, z)
    # 或者从外部源获取一个更复杂的投影矩阵

# 保存修改后的 DXF 文件
output_filename = "tester_transformed.dxf"
doc.saveas(output_filename)
print(f"转换完成，文件已保存为: {output_filename}")

import ezdxf
from ezdxf import transform
from ezdxf.math import Matrix44

# 配置：True 表示从 CRS 转换为 WCS，False 表示从 WCS 转换为 CRS
CRS_TO_WCS = True 

def wcs_to_crs(entities, m: Matrix44):
    """将实体从 WCS 坐标转换为 CRS 坐标。"""
    transform.inplace(entities, m)

def crs_to_wcs(entities, m: Matrix44):
    """将实体从 CRS 坐标转换为 WCS 坐标。"""
    m_inverse = m.copy()
    m_inverse.inverse()
    transform.inplace(entities, m_inverse)

def main():
    input_filename = "tester.dxf"
    output_filename = "tester_transformed.dxf"

    try:
        doc = ezdxf.readfile(input_filename)
        msp = doc.modelspace()
        geo_data = msp.get_geodata()

        if geo_data:
            m, epsg = geo_data.get_crs_transformation()
            print(f"从 '{input_filename}' 读取到 GEODATA，EPSG: {epsg}")

            if CRS_TO_WCS:
                print("正在将模型空间实体从 CRS 转换为 WCS...")
                crs_to_wcs(msp, m)
            else:
                print("正在将模型空间实体从 WCS 转换为 CRS...")
                wcs_to_crs(msp, m)

            doc.saveas(output_filename)
            print(f"转换完成，修改后的文件已保存为: {output_filename}")

        else:
            print(f"'{input_filename}' 中没有地理参考数据 (GEODATA)。")
            print("若要进行坐标转换，您需要手动提供一个有效的转换矩阵。")
            print("当前操作将不进行任何坐标转换，文件将原样保存。")
            # 示例：如果已知外部 CRS 信息并想手动构建矩阵
            # 例如，如果您知道需要将一个特定CRS的坐标平移或缩放
            # manual_transform_matrix = Matrix44.translate(100, 200, 0)
            # crs_to_wcs(msp, manual_transform_matrix) # 假设这是CRS到WCS的矩阵

            # 如果不进行转换，可以只保存原文件或跳过保存
            # doc.saveas("tester_no_geodata_original.dxf") 

    except FileNotFoundError:
        print(f"错误：文件 '{input_filename}' 未找到。请检查文件路径。")
    except Exception as e:
        print(f"发生错误: {e}")

if __name__ == "__main__":
    main()