Go go-geom 空几何与零坐标点有什么区别

在 go-geom 里，空几何和坐标为零的点不是一回事：NewPointEmpty(geom.XY) 表示没有坐标，NewPoint(geom.XY) 则创建一个坐标值全为 0 的有效点，也就是二维坐标 (0,0)。如果把两者都当成“没有位置”，导入数据、计算距离或写回 GeoJSON 时就会出现边界错误。

官方文档：https://pkg.go.dev/github.com/twpayne/go-geom

这篇只讨论 Point、GeometryCollection 的空状态，以及如何在业务代码中做判定，不把空几何扩展成坐标系统或拓扑运算教程。

空几何和零坐标点分别代表什么

NewPointEmpty(layout) 分配的是“没有坐标”的 Point；NewPoint(layout) 分配的是对应布局下值为零的 Point。前者的重点是坐标数量为零，后者虽然数值为零，但仍然有一个完整坐标。

在 geom.XY 布局下，可以用下面的对照理解它们的契约：

用 Empty 和 NumCoords 做第一层判断

判断空状态时，优先看几何对象的 Empty()，再结合 NumCoords() 和布局确认数据形状。不要直接用 X()==0 && Y()==0 代替空判断，因为这会把合法原点误判成缺失。

package main

import (
    "fmt"

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

func describePoint(label string, p *geom.Point) {
    // Empty 表示没有坐标；NumCoords 用来补充检查对象是否有一个坐标。
    fmt.Printf("%s: empty=%v coords=%d layout=%s\n", label, p.Empty(), p.NumCoords(), p.Layout())
    if !p.Empty() {
        // 只有非空点才读取坐标，避免把缺失位置继续传给业务计算。
        fmt.Printf("  coord=(%.1f, %.1f)\n", p.X(), p.Y())
    }
}

func main() {
    emptyPoint := geom.NewPointEmpty(geom.XY)
    origin := geom.NewPoint(geom.XY)

    describePoint("empty", emptyPoint)
    describePoint("origin", origin)
}

这里的关键不是打印格式，而是分支顺序：先用 Empty() 判断有没有坐标，再读取 X()、Y()。这样原点仍会进入后续计算，而真正的空点会走缺失分支。

不要把空坐标的 NaN 当成原点

go-geom 为 Point 空坐标提供了 PointEmptyCoord()，文档说明它是表示空坐标的 NaN 浮点数。它用于表达空坐标的内部约定，不应该被业务代码改写成普通的 0。

如果数据来自平铺坐标切片，可以使用 NewPointFlatMaybeEmpty 让构造函数检查坐标是否全部为 NaN，并据此识别可能的空点。相反，真正的 (0,0) 仍应作为有坐标的点保留。

// 只演示判定边界，不把示例输出当作运行截图。
func isMissingPoint(p *geom.Point) bool {
    // 业务含义由 Empty 决定，而不是由 X/Y 是否等于零决定。
    return p.Empty()
}

func isOrigin(p *geom.Point) bool {
    // 先排除空点，再判断合法原点，避免空坐标参与数值比较。
    return !p.Empty() && p.X() == 0 && p.Y() == 0
}

GeometryCollection 全为空时也要继续看成员

GeometryCollection 的数量和它是否“有业务位置”是两件事。官方文档明确说明：集合里可以包含多个 Geometry，但如果这些对象全部为空，GeometryCollection.Empty() 仍然可以返回 true。

因此，集合非空不等于集合包含有效坐标。需要保留成员索引、判断成员类型或区分“空成员”和“没有成员”时，再读取 NumGeoms()、Geom(i) 或 Geoms()。

func hasUsableGeometry(gc *geom.GeometryCollection) bool {
    // 空集合和“包含但全部为空的成员”都不能直接当作有效位置。
    if gc.Empty() {
        return false
    }
    // 非空集合仍需按成员语义交给上层处理，这里只返回集合级结论。
    return gc.NumGeoms() > 0
}

如果业务要求“至少有一个可计算的成员”，还应遍历 Geoms()，对每个成员使用其具体类型的 Empty()。不要只检查 NumGeoms() > 0。

把判定规则放在导入和计算边界

更稳定的做法是在数据进入距离、范围、索引或序列化逻辑前，先把几何状态归类：

• 缺失位置：对象为空，进入缺失、补采或跳过分支。
• 合法原点：对象非空且坐标为 (0,0)，继续参与计算。
• 普通位置：对象非空且坐标不为零，按正常流程处理。
• 集合对象：先判断集合级空状态，再按需要检查成员级空状态。

这个顺序也能避免把 NaN、零坐标和空切片混在一个“没有数据”的布尔值里。若要写回数据库或 GeoJSON，保留空几何的语义；若要做数值运算，则只把通过非空判断的坐标交给算法。

常见问题

为什么不能用 X 和 Y 都为 0 判断空点？

因为 (0,0) 是有效坐标。空点没有坐标，应该用 Empty() 判断。

GeometryCollection 有成员就一定不为空吗？

不一定。成员可以全部为空，集合的 Empty() 仍可能为 true；要区分成员数量和有效成员，需要继续检查成员对象。

空点能不能直接参与距离计算？

不建议。先在计算边界判断 Empty()，再决定跳过、补值或返回业务错误。

记住一句话：空几何描述“没有坐标”，零坐标点描述“坐标恰好为 (0,0)”。在 go-geom 中把 Empty() 放在数值判断之前，通常就能避开这类最隐蔽的边界问题。