如何将NumPy数组导出为AutoCAD可识别的曲面信息文件

1. 导出带地理参考的GeoTIFF（优先推荐） #

你之前用tifffile.imsave导出的是普通图像TIFF，丢失了原LiDAR影像的地理变换和投影信息，导致AutoCAD无法识别为地理空间数据。正确的做法是用专门处理栅格地理数据的库（比如rasterio），保留原影像的地理元数据后导出GeoTIFF：

示例代码 #

import rasterio
from rasterio.mask import mask
import geopandas as gpd

# 读取原LiDAR影像与边界文件
Existing_Surface = 'FilePath/Existing_Surface.tif'
Fill_Boundary = 'FilePath/Fill_Boundary_polygon.shp'

with rasterio.open(Existing_Surface) as src:
    # 读取裁剪边界
    boundary_gdf = gpd.read_file(Fill_Boundary)
    # 执行裁剪，同时保留裁剪后的地理变换与投影信息
    cropped_array, out_transform = mask(src, boundary_gdf.geometry, crop=True)
    # 复制原影像的元数据模板
    out_meta = src.meta.copy()
    
    # 更新元数据以匹配裁剪后的数据
    out_meta.update({
        "height": cropped_array.shape[1],
        "width": cropped_array.shape[2],
        "transform": out_transform
    })

# 对裁剪后的数组执行你的修改操作（替换为你的Proposed_example_array逻辑）
Proposed_example_array = cropped_array[0]  # 假设是单波段影像
# ... 这里写你修改数组的代码 ...

# 导出带地理参考的GeoTIFF
with rasterio.open("FilePath/Proposed_Surface.tif", "w", **out_meta) as dest:
    dest.write(Proposed_example_array, 1)

导出的GeoTIFF可以直接在AutoCAD/Civil 3D中通过「插入光栅图像参照」并启用地理定位功能，或者用Civil 3D的「导入DTM」工具识别为曲面。

2. 导出为CSV点云格式 #

如果需要更灵活地在AutoCAD中构建曲面，可以将数组转换为带地理坐标的XYZ点云CSV：

示例代码 #

import numpy as np
import csv

# 从原影像获取地理变换（裁剪时已获取的out_transform）
geotransform = out_transform
rows, cols = Proposed_example_array.shape

# 导出CSV
with open("FilePath/Proposed_Surface.csv", "w", newline='') as csv_file:
    writer = csv.writer(csv_file)
    writer.writerow(["X", "Y", "Z"])  # 写入表头
    
    for row_idx in range(rows):
        for col_idx in range(cols):
            z_value = Proposed_example_array[row_idx, col_idx]
            # 跳过无效值（比如你的示例中的0）
            if z_value == 0:
                continue
            # 计算像素对应的真实地理坐标
            x = geotransform[0] + col_idx * geotransform[1] + row_idx * geotransform[2]
            y = geotransform[3] + col_idx * geotransform[4] + row_idx * geotransform[5]
            writer.writerow([x, y, z_value])

在AutoCAD中，你可以用「点」命令批量导入这些点，或者在Civil 3D中使用「从点文件创建曲面」工具直接生成曲面。

3. 导出为Shapefile点集 #

若需要与GIS工具兼容并导入AutoCAD，可导出为带Z值的Shapefile点集：

示例代码 #

import geopandas as gpd
from shapely.geometry import Point

points_list = []
z_values = []

rows, cols = Proposed_example_array.shape
for row_idx in range(rows):
    for col_idx in range(cols):
        z_value = Proposed_example_array[row_idx, col_idx]
        if z_value == 0:
            continue
        # 计算地理坐标
        x = geotransform[0] + col_idx * geotransform[1] + row_idx * geotransform[2]
        y = geotransform[3] + col_idx * geotransform[4] + row_idx * geotransform[5]
        points_list.append(Point(x, y))
        z_values.append(z_value)

# 创建GeoDataFrame并导出Shapefile
points_gdf = gpd.GeoDataFrame(
    data={"Z": z_values},
    geometry=points_list,
    crs=src.crs  # 继承原影像的投影
)
points_gdf.to_file("FilePath/Proposed_Surface_Points.shp")

导入AutoCAD后，可基于这些点生成TIN曲面，Civil 3D支持直接识别Shapefile中的点数据来构建曲面。

内容的提问来源于stack exchange，提问作者Patstro