如何使用C#读取Shape文件存入PostgreSQL并实现重投影

SHP读取写入PostgreSQL实现 #

直接基于NetTopologySuite生态开发即可，不用从零实现SHP格式解析：SHP是一套多文件组合格式，除了主shp文件外还要同步处理索引shx、属性dbf、坐标系prj文件，自己手写解析很容易踩格式兼容、字段类型映射的坑，这套类库已经把所有底层逻辑封装好了。

• 读取环节：用类库自带的ShapefileDataReader加载目标SHP文件，指定匹配的属性编码（国内多数SHP用GBK编码，不指定会出现中文属性乱码），即可逐行遍历拿到所有要素的几何信息、属性字段值；遍历过程中加一层校验，直接过滤掉几何为空、坐标自相交的无效脏数据，避免后续写库报错。
• 入库环节：PostgreSQL提前启用PostGIS扩展，建好对应业务表，几何字段提前指定好目标SRID；用Npgsql连接数据库，搭配NetTopologySuite的PostGIS适配包，直接把读取到的几何对象作为SQL参数传入即可，不需要手动转WKT/WKB格式，类库会自动完成类型映射。如果数据量较大，直接用数据库批量复制接口写入，效率比单条INSERT高数十倍。

核心读取逻辑参考：

using var shpReader = new ShapefileDataReader(
    path: @"D:\data\test.shp",
    geometryFactory: new GeometryFactory(),
    encoding: Encoding.GetEncoding("GBK")
);
var writeBatch = new List<Feature>();
while (shpReader.Read())
{
    var geom = shpReader.Geometry;
    // 跳过无效数据
    if (geom.IsEmpty || !geom.IsValid) continue;
    // 读取属性字段，比如var adminCode = shpReader["adcode"].ToString();
    writeBatch.Add(new Feature { Geom = geom, Props = props });
    // 凑够批次就提交写库
}

空间数据重投影实现 #

完全可以在这个C#工具内直接完成空间数据重投影，不需要依赖QGIS、GDAL等外部工具。

• 实现逻辑：用NetTopologySuite配套的投影转换模块即可，模块内置了绝大多数常用坐标系的EPSG定义，支持不同地理坐标系、投影坐标系之间的双向转换。读取SHP时类库会自动解析prj文件拿到源数据的坐标系，你只需要指定入库要求的目标坐标系，初始化坐标转换实例，对每一条读取到的几何对象执行转换即可。
• 注意事项：投影转换完成后记得再做一次几何有效性校验，个别跨带转换、边界要素转换后可能出现自相交、坐标越界的问题，校验后再写入数据库能避免后续空间查询出错。如果你的业务用的是自定义非标准坐标系，只要拿到坐标系的WKT定义字符串，也可以手动加载创建转换实例，不受内置EPSG库的限制。

实操提醒：批量写库时单批次提交的要素数量控制在1000-5000条即可，批次过大会占用过多内存，批次太小会频繁和数据库交互拖慢写入速度。

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