Go go-shp 怎样遍历Shapefile中的几何记录

我第一次把 Shapefile 接进 Go 服务时，最容易写错的不是几何类型，而是把“读取一条记录”误解成只拿一个坐标对象。对 github.com/jonas-p/go-shp 来说，一次读取同时要关注几何记录、记录编号和 DBF 属性行。稳定的写法是：先 Open，再用 Next 驱动循环，在循环内调用 Shape，结束后检查 Err。

官方地址：https://pkg.go.dev/github.com/jonas-p/go-shp

go-shp 遍历 Shapefile 的核心是 for reader.Next()：Shape() 给出当前几何和从 0 开始的记录号，DBF 属性再用这个记录号配合字段下标读取；循环结束后不要把 false 直接当成成功，要继续检查 Err()。

先把 go-shp 的读取对象和文件关系理顺

安装依赖时直接使用模块路径：

# 安装 go-shp 依赖，写入当前 Go module
go get github.com/jonas-p/go-shp

一个可遍历的 Shapefile 通常至少要有主几何文件和索引文件；如果还要读属性，就准备对应的 .dbf。下面的示例假定它们共享 roads 这个基础名，代码只演示读取，不把任何执行结果冒充为图片证据。

这里先记住三个对象关系：Reader 管理文件读取位置，Shape 表示当前几何，Fields 描述 DBF 列。不要把字段定义和字段值混为一谈，前者告诉你列是什么，后者要按当前记录行去取。

用 Next 和 Shape 逐条取得几何记录

Next 每成功推进一次，Shape 就能拿到最近读取的对象和它的零基编号。这个编号非常重要：后面读取 DBF 时，它就是几何与属性对齐的桥。

示例把具体类型处理留在一个小函数里，主循环只负责推进和记录生命周期：

package main

import (
    "fmt"
    "log"

    shp "github.com/jonas-p/go-shp"
)

func main() {
    reader, err := shp.Open("roads.shp")
    if err != nil {
        log.Fatal(err) // 打开失败时不要进入空循环
    }
    defer reader.Close() // Reader 持有的文件句柄由这里统一释放

    for reader.Next() {
        index, geometry := reader.Shape()
        fmt.Printf("record=%d bbox=%+v\n", index, geometry.BBox())
        printGeometryKind(geometry)
    }

    if err := reader.Err(); err != nil {
        log.Fatal(err) // Next 返回 false 后，区分 EOF 与真实读取错误
    }
}

func printGeometryKind(geometry shp.Shape) {
    switch geometry.(type) {
    case *shp.Point:
        fmt.Println("kind=point")
    case *shp.PolyLine:
        fmt.Println("kind=polyline")
    case *shp.Polygon:
        fmt.Println("kind=polygon")
    default:
        fmt.Println("kind=other") // 未覆盖类型不静默当作点处理
    }
}

上面的 BBox 只适合快速确认范围，不等于把所有顶点都取出来。需要分析点、线或面时，再进入具体类型的结构；如果只是统计记录数、做范围筛选，先保留 Shape 接口反而更稳。

根据几何类型选择处理方式

Shapefile 的同一文件通常保持统一的形状类型，但业务代码仍应显式处理类型。点对象适合直接取坐标；折线和面则要继续看它们的 parts、points 等结构。把所有几何都强行转成点，会让空值、边界和后续空间计算变得难以排查。

工具选型上，go-shp 的低层 Reader 适合“逐条读取、自己决定转换”的任务。如果你只想依次拿到几何和属性，可以考虑它提供的顺序读取接口；如果需要投影转换、GeoJSON 或更完整的空间模型，应该把读取层与后续 GIS 库分开，别期待 Reader 自动完成坐标系转换。

用 Fields 和 ReadAttribute 读取 DBF 属性并检查错误

字段和值分别来自两个维度。Fields() 返回列定义，ReadAttribute(row, field) 按零基行号和零基字段号返回字符串。把 Shape() 返回的 index 传给第一个参数，就能维持几何与属性的一一对应。

fields := reader.Fields()
for reader.Next() {
    index, geometry := reader.Shape()
    fmt.Printf("record=%d bbox=%+v\n", index, geometry.BBox())

    for fieldIndex, field := range fields {
        value := reader.ReadAttribute(index, fieldIndex)
        fmt.Printf("%s=%s\n", field.String(), value)
    }
}
if err := reader.Err(); err != nil {
    return fmt.Errorf("read shapefile: %w", err) // 保留底层错误便于定位文件问题
}

实际项目中可以先筛字段，再把字符串转换成业务类型。不要在循环里忽略字段顺序，也不要把字段名当作行号；DBF 的列由 Fields 的位置决定，行由当前几何记录的 index 决定。

常见问题

为什么 Next 返回 false 还要调用 Err？

因为 false 既表示文件已读完，也表示读取过程中遇到错误。调用 Err 才能把正常结束和损坏文件、配套文件不一致等情况分开。

Shape 返回的编号从几开始？

它从 0 开始。读取 DBF 时，使用这个编号作为行参数，再传入字段的零基下标。

没有 DBF 文件还能遍历几何吗？

可以读取几何，但属性字段和值不可用。若业务依赖属性过滤，应在打开前确认配套文件完整，并把缺失属性视为输入条件问题处理。