PostGIS：point列转换为ST_DWithin可用类型时解析错误，求解决方案

我需要统计指定范围内的地点总数，我的places表包含一个类型为geo point default null的geo列。执行以下查询：

select
  count(*)
from
  places
where
  ST_DWithin(
    ST_GeomFromText(concat('POINT', geo::text), 4326),
    ST_GeomFromText('POINT(37.64903402157866, -83.84765625000001)', 4326),
    66 * 1609.34
  );

返回错误：

parse error - invalid geometry

尝试过转为::geometry及其他函数均无效。执行select geo, geo::text, geo::geometry from places where geo is not null limit 10返回的数据示例：

{"x":-84.3871,"y":33.7485}
(-84.3871,33.7485)
01010000000612143FC61855C02B8716D9CEDF4040

疑问：哪里操作有误？是否需要修改表列类型？如何避免丢失现有数据？

1. 错误原因 #

• 冗余转换：无需将geo转成文本再用ST_GeomFromText解析，geo::geometry已经是有效的PostGIS几何对象，直接使用即可。
• 格式错误：即使要转文本，concat('POINT', geo::text)生成的字符串不符合WKT规范——geo::text返回的是(X,Y)格式，拼接后变成POINT(X,Y)，而标准WKT点格式要求坐标用空格分隔而非逗号，这导致了几何解析失败。

2. 修复后的查询 #

直接使用geo::geometry作为空间函数参数，这是最简洁高效的写法：

select
  count(*)
from
  places
where
  ST_DWithin(
    geo::geometry,
    ST_SetSRID(ST_MakePoint(37.64903402157866, -83.84765625000001), 4326),
    66 * 1609.34
  );

注：ST_MakePoint+ST_SetSRID比ST_GeomFromText性能更好，也更不容易出错。

3. 关于列类型优化 #

如果业务频繁涉及空间查询，建议将geo列改为PostGIS标准的geometry(Point, 4326)类型，这样可以：

• 避免每次查询时的类型转换开销
• 支持创建空间索引，大幅提升查询速度

4. 安全修改列类型（无数据丢失） #

修改前务必先备份数据，再执行以下步骤：

1. 备份表数据：

CREATE TABLE places_backup AS SELECT * FROM places;

2. 修改列类型：

ALTER TABLE places ALTER COLUMN geo TYPE geometry(Point, 4326) USING geo::geometry;

3. 创建空间索引（可选，优化查询性能）：

CREATE INDEX idx_places_geo ON places USING GIST (geo);

USING geo::geometry会自动将原point类型转换为对应的PostGIS几何对象，从你的数据示例来看，转换完全兼容，不会丢失数据。

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