Go go-geom 怎样用平铺坐标构造二维折线

如果手里已经有一串按 X、Y、X、Y 排列的 float64，Go 的 go-geom 可以直接用 geom.NewLineStringFlat 构造二维折线。关键不是把切片再套一层，而是让 layout 和切片的分组方式一致：geom.XY 的步长是 2，所以每两个数才代表一个控制点。

官方地址：https://pkg.go.dev/github.com/twpayne/go-geom

二维场景优先使用 geom.XY，把平铺切片按每两个元素分组，再传给 NewLineStringFlat。构造后用 NumCoords、Coords 和 Length 读取几何信息；如果数据来自会继续复用的缓冲区，要先处理切片别名问题。

先把“平铺坐标”对齐到 XY 布局

go-geom 用 Layout 描述每个坐标包含哪些维度。XY 只有横坐标和纵坐标，Stride() 返回 2；XYZ 则返回 3。也就是说，下面这串数据不是六个独立点，而是三个二维控制点：

flat := []float64{
    10, 20, // 第一个点：X=10，Y=20
    14, 23, // 第二个点：X=14，Y=23
    19, 21, // 第三个点：X=19，Y=21
}

当布局是 geom.XY 时，索引 0、1 属于第一个点，索引 2、3 属于第二个点，以此类推。布局和数据不能各说各话：把三维数据交给 XY 会改变分组含义，把二维数据交给 XYZ 则无法按预期还原控制点。

用 NewLineStringFlat 构造二维折线

NewLineStringFlat 接收布局和一维浮点切片，返回一个 *geom.LineString。对于二维折线，最小写法如下：

package main

import (
    "fmt"

    "github.com/twpayne/go-geom"
)

func main() {
    flat := []float64{
        10, 20, // 第一个二维控制点
        14, 23, // 第二个二维控制点
        19, 21, // 第三个二维控制点
    }

    // XY 决定每两个 float64 组成一个坐标。
    line := geom.NewLineStringFlat(geom.XY, flat)

    // NumCoords 返回控制点数量，而不是 flat 切片的元素数量。
    fmt.Println("layout:", line.Layout())
    fmt.Println("stride:", line.Stride())
    fmt.Println("points:", line.NumCoords())
    fmt.Println("length:", line.Length())
}

这里的 points 应按坐标组理解：平铺切片长度为 6，geom.XY 的步长为 2，因此控制点数量是 3。LineString 表示一条由控制点线性连接而成的连续折线，不等同于一个闭合面；首尾点是否相同，应由输入数据和业务语义决定。

构造完成后怎样读取坐标和长度

构造后常见的三类读取需求分别是：按结构读取坐标、按原始形态读取平铺数据、计算折线长度。Coords() 适合逐点处理，FlatCoords() 适合和编码器或已有缓冲区衔接，Length() 则按照当前布局的步长计算折线长度。

func summarize(line *geom.LineString) {
    // Coords 把平铺数据还原成逐点的 Coord，便于业务层遍历。
    for i, coord := range line.Coords() {
        fmt.Printf("point[%d]=(%v, %v)\n", i, coord[0], coord[1])
    }

    // FlatCoords 保留 X、Y、X、Y 的原始排列，适合交给编码层。
    fmt.Println("flat:", line.FlatCoords())

    // Length 使用当前几何的坐标布局和步长计算线段总长。
    fmt.Println("length:", line.Length())
}

如果只关心点数，优先用 NumCoords()，不要用 len(line.FlatCoords()) 代替。后者得到的是标量数量，在 XYZ、XYM 或 XYZM 场景下尤其容易把维度数误当成点数。

什么时候改用 SetCoords

如果输入本来就是一组 geom.Coord，直接用 NewLineString 加 SetCoords 更容易表达意图。它把“每个点是一条 Coord”保留在调用侧，不需要调用方自己维护平铺索引。

func buildFromCoords() (*geom.LineString, error) {
    // 结构化坐标适合来自 JSON、表单或逐点计算的输入。
    coords := []geom.Coord{
        {10, 20},
        {14, 23},
        {19, 21},
    }

    line := geom.NewLineString(geom.XY)
    // SetCoords 会按 line 的布局整理这些坐标，错误交给调用方处理。
    if _, err := line.SetCoords(coords); err != nil {
        return nil, fmt.Errorf("设置折线坐标失败: %w", err)
    }
    return line, nil
}

可以这样选择：数据已经是连续的 float64 缓冲区时使用 NewLineStringFlat；数据以逐点对象存在，或者需要在写入前逐点检查维度时使用 SetCoords。两条路径最后都得到 LineString，区别主要在输入数据的边界和转换成本。

两个容易踩到的边界

平铺切片长度必须匹配布局步长

使用 geom.XY 时，平铺切片应该能按 2 个元素完整分组。不要把一个额外的 Z 值塞在某个点后面，也不要把点数组和布局常量拆开由不同层随意决定。需要三维坐标时，明确改用 geom.XYZ，并让每组三个值保持一致。

注意底层切片的复用和修改

NewLineStringFlat 接收的是现成的 flatCoords 切片。若这块内存随后被复用、清空或填充下一条记录，几何对象看到的数据也可能随之变化。需要独立生命周期时，先复制一份再交给构造器；只读且生命周期明确的缓冲区则可以直接传入。

func cloneFlatCoords(src []float64) *geom.LineString {
    // 复制切片，让 LineString 不依赖上游缓冲区后续的写入。
    flat := append([]float64(nil), src...)
    return geom.NewLineStringFlat(geom.XY, flat)
}

另外，SetCoords 和 MustSetCoords 的错误语义不同：前者返回错误，适合处理外部输入；后者在坐标不符合布局时会触发 panic，更适合调用方已经保证数据形状的内部路径。

一张表记住选择

常见问题

为什么 flat 切片长度是 6，但折线只有 3 个点？

因为 geom.XY 的步长是 2，6 个标量会分成 3 组二维坐标。点数应使用 NumCoords() 获取。

平铺坐标可以直接传给 NewLineStringFlat 吗？

可以，但要保证布局和每个坐标包含的维度一致。构造器适合已有平铺缓冲区的场景，数据来自逐点输入时可改用 SetCoords。

修改原来的 flat 切片会不会影响 LineString？

可能会。若切片被复用或后续会写入，先复制数据再构造独立的折线对象。

LineString 会自动闭合吗？

不会。它表达的是连续折线；如果业务需要闭合环，应明确使用线性环或让首尾坐标相同，并按目标几何类型处理。

总结起来，Go go-geom 构造二维折线的关键只有三件事：用 geom.XY 表达二维布局，按两个浮点数一组组织平铺坐标，再用 NewLineStringFlat 创建 LineString。随后用 NumCoords、Coords、FlatCoords 和 Length 读取结果；如果输入缓冲区会变化，就在构造前明确复制边界。