(GIS开发) Python dxf 文件解析

GIS程序媛—椰子

已于 2023-03-06 15:27:47 修改

阅读量1.2k

点赞数 2

CC 4.0 BY-SA版权

文章标签： python

于 2022-10-18 17:47:04 首次发布

本文链接：https://blog.youkuaiyun.com/OTTOkongbai/article/details/127392750

该代码实现从DXF文件中读取图层和块，并将其转换为SHP地理空间文件。使用geopandas库进行转换，dxfgrabber库读取图层和块信息。每个图层的数据被单独导出为SHP文件，同时对列名进行截断以符合长度限制。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

dxf2shp

import os

import geopandas as gpd


def dxf2shp(file):
    gdf = gpd.read_file(file)
    gdf.crs = {'init': 'epsg:3414'}
    shape = gdf.to_crs(epsg=4326)
    shape['Layer'].value_counts()
    arr = shape['Layer'].unique()

    for i in shape.columns:
        if (len(i) > 10):
            new_name = i[:9]
            shape.rename(columns={i: new_name}, inplace=True)

    for i in arr:
        export = shape[(shape.Layer == i)]
        if (export.shape[0] > 0):
            try:
                export.to_file("shp/" + i + ".shp")
            except Exception as e:
                print("### dxf2shp(): error: ", e)
                pass
        else:
            print("### dxf2shp(): export.shape = 0 : ", i)
        continue
    print("### dxf2shp(): done")


if __name__ == "__main__":
    try:
        if not os.path.exists('shp'):
            os.mkdir('shp')
        dxf2shp('XX.dxf')
        print("### main(): done")
    except Exception as e:
        print("!!! Error found: ", e)
        pass

读取 layers 和 blocks

import dxfgrabber


def printDxf(file):
    dxf = dxfgrabber.readfile(file)
    for layer in dxf.layers:
        print("### layer: ", layer.name)
    for block in dxf.blocks:
        print("### block: ", block.name)


if __name__ == "__main__":
    try:
        printDxf('XX.dxf')
        print("### main(): done")
    except Exception as e:
        print("!!! Error found: ", e)
        pass