Go go-geom SRID字段为什么不会自动变换坐标

在 go-geom 里把几何对象的 SRID 从 4326 改成 3857，坐标数字仍然不变，这是正常行为。SRID 是对象携带的参考系标识，SetSRID 只负责写入这个标识；它不会根据 EPSG 定义执行投影、单位换算或坐标重算。

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

如果目标是“声明这组坐标属于哪个参考系”，使用 SetSRID；如果目标是“把经纬度换算成另一套坐标数字”，必须显式调用坐标变换逻辑，再把结果标记为目标 SRID。只改 SRID，不能代替投影变换。

为什么改 SRID 后坐标不动

可以把一个 geom.T 理解成两层信息：一层是 FlatCoords() 返回的坐标值，另一层是 SRID() 返回的整数标识。前者决定点、线、面的数字位置，后者告诉下游系统应该怎样解释这些数字。

例如，(116.391, 39.907) 这组数字通常被当作经度和纬度；把 SRID 从 4326 改成 3857，只是把标签换了，并不会自动产生米制投影坐标。继续把原数字当成 3857 使用，反而会造成单位和范围判断错误。

SetSRID 实际只修改哪个字段

go-geom 为具体几何类型提供了 SetSRID 方法，也提供了面向 geom.T 的通用函数。它们的共同语义是设置 SRID，不会遍历并改写平铺坐标。

package main

import (
	"fmt"

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

func main() {
	// 用经度、纬度创建一个二维点，并先标记为 EPSG:4326。
	point := geom.NewPointFlat(geom.XY, []float64{116.391, 39.907}).SetSRID(4326)
	before := append(geom.Coord(nil), point.Coords()[0]...)

	// 这里只替换参考系标识，不执行任何投影计算。
	point.SetSRID(3857)

	fmt.Println("before:", before)
	fmt.Println("after:", point.Coords()[0])
	fmt.Println("srid:", point.SRID())
}

这段代码中，after 仍会打印原来的两个数字，而 srid 变成 3857。这个结果正好说明了职责边界：应用负责确认原始数字确实来自哪个参考系，库负责保存这个声明。

真正变换坐标要显式改写数值

需要把坐标从一个参考系转换到另一个参考系时，应该把“坐标计算”和“SRID 标记”写成两个明确动作。geom.TransformInPlace 会按几何对象的 stride 访问每个坐标，并把回调对坐标切片的修改保留下来；它本身也不是 EPSG 投影数据库。

下面用简单平移模拟外部投影函数的回调位置。平移不是 4326 到 3857 的真实算法，实际项目应在回调中接入已经选定的投影实现，并根据其错误模型处理失败。

package main

import (
	"fmt"

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

func main() {
	// 原始点仍然按 XY 保存，4326 只是它当前的参考系标识。
	point := geom.NewPointFlat(geom.XY, []float64{116.391, 39.907}).SetSRID(4326)

	// 这里用平移演示“回调会改坐标”；这不是实际的 EPSG 投影公式。
	geom.TransformInPlace(point, func(c geom.Coord) {
		c[0] += 1000
		c[1] += 1000
	})

	// 坐标已经被显式改写后，再声明它属于目标参考系。
	point.SetSRID(3857)
	fmt.Println(point.Coords()[0], point.SRID())
}

如果使用的是独立投影库，通常应让它负责计算新坐标，go-geom 只负责承载几何结构。不要为了让输出“看起来像目标坐标”而直接把 SRID 改掉，也不要把示例中的平移常数当成真实投影参数。

EWKB 会保留 SRID，但也不会替你变换

go-geom 的 EWKB 编解码支持携带 SRID。写出非零 SRID 时，编码会把标识放进 EWKB 的类型信息中；读取时再把它设置回几何对象。这个过程解决的是跨进程或数据库传输时的元数据保存，不是坐标系转换。

package main

import (
	"fmt"

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

func main() {
	// 让 EWKB 同时看到坐标值和非零 SRID 元数据。
	point := geom.NewPointFlat(geom.XY, []float64{116.391, 39.907}).SetSRID(4326)
	data, err := ewkb.Marshal(point, ewkb.NDR)
	if err != nil {
		// 编码失败时立即返回，避免把不完整字节交给下游。
		panic(err)
	}

	// 解码会恢复 SRID 和原始坐标，但不会自动投影。
	decoded, err := ewkb.Unmarshal(data)
	if err != nil {
		// 输入来自网络或数据库时，解码错误应交给调用方处理。
		panic(err)
	}
	fmt.Println(decoded.SRID(), decoded.FlatCoords())
}

因此，EWKB 编解码后的 SRID 一致，只能说明元数据被保留下来；判断坐标是否已经转换，仍然要查看坐标数字和负责变换的业务代码。

迁移坐标系时最容易犯的几个错误

工程上可以把转换函数设计成“输入几何、源 SRID、目标 SRID、输出几何”的边界，避免在多个 handler 或存储层隐式改标签。对于原地修改，还要确认调用方没有继续持有同一个对象并把旧坐标当成未转换数据使用。

结论与快速判断

看到 SetSRID 时，把它理解为“给几何对象贴参考系标签”；看到 TransformInPlace 或外部投影函数时，才把它理解为“实际改写坐标值”。前者不会替代后者，EWKB 也只是帮助前者在传输中被保留下来。

• 只想声明输入坐标的参考系：设置并读取 SRID。
• 需要换成另一套坐标数字：先执行显式变换，再设置目标 SRID。
• 需要跨数据库或服务传输：确认编码格式同时携带 SRID，解码后仍要区分标签与坐标。

常见问题

go-geom 会根据 SRID 自动选择投影算法吗？
不会。SRID 是元数据，库不会因为整数从 4326 变成 3857 就自动查找投影公式。

调用 TransformInPlace 后会自动更新 SRID 吗？
不要依赖这种假设。回调负责修改坐标，目标 SRID 应由业务在确认变换成功后显式设置。

EWKB 解码成功是否代表坐标已经转换？
不代表。它说明几何和可编码的 SRID 信息被恢复；是否转换过，要看坐标值和实际变换逻辑。