1. 项目概述：为什么几何图形创建是PostGIS真正的“入门门槛”
刚接触PostGIS的人，常以为装好扩展、连上数据库、跑通SELECT ST_AsText(geom)就算入门了。我带过不少从QGIS或ArcGIS转过来的同事，他们能熟练画点线面，却在第一次用SQL生成一个缓冲区时卡住——不是语法报错，而是根本想不明白“这个多边形坐标串到底该怎么写才合法”“为什么ST_MakePolygon报错说环不闭合”“ST_Buffer出来的结果怎么和界面上画的看起来不一样”。这背后暴露的，是一个被严重低估的事实：PostGIS的几何图形创建，不是简单的函数调用，而是一套空间数据建模的底层思维训练。它直接决定你后续所有空间分析的可靠性、索引效率、跨系统兼容性，甚至影响Web地图服务的渲染性能。我参与过某城市地下管网模拟项目，前期因一个ST_GeomFromText中坐标顺序写反（顺时针 vs 逆时针），导致整个区域的拓扑关系错乱，排查了三天才定位到根源——问题不在分析逻辑，而在最开始那个手工拼接的WKT字符串。所以这篇内容，不讲高深的空间算法，只聚焦“如何稳、准、快地把点、线、面、环、多部件图形，在PostGIS里一砖一瓦地搭出来”。你会看到：WKT/WKB的本质区别不是格式，而是坐标系绑定方式；ST_Point和ST_MakePoint看似一样，实则一个强制校验SRID一个完全不校验；ST_Collect和ST_Union在创建复合图形时，前者保留原始部件拓扑，后者直接融合成新几何体——这些细节，决定了你写的SQL是能跑通，还是能跑得稳、跑得久、跑得对。
2. 核心思路拆解：从“画图思维”到“建模思维”的三重转换
2.1 为什么不能直接照搬GIS软件里的操作逻辑？
在QGIS里画一个多边形，鼠标点完自动闭合，软件内部帮你做了坐标顺序校验、环方向归一化、无效几何修复。但PostGIS没有“界面”，它只认严格符合OGC标准的几何对象。我试过把QGIS导出的GeoJSON直接转WKT再导入PostGIS，结果70%的多边形报错Geometry has Z dimension but column does not——因为QGIS默认导出带Z值（高程），而目标表字段定义是GEOMETRY(POINT,4326)，没声明Z维度。这说明：PostGIS的几何创建，本质是“协议级建模”，不是“绘图级操作”。你必须主动思考三个维度：
- 维度维度：2D（X,Y）、3D（X,Y,Z）、4D（X,Y,Z,M）必须与字段定义完全匹配。
ST_MakePoint(1,2,3)生成的是3D点，若插入2D字段，PostGIS会静默丢弃Z值，但后续ST_3DDistance等函数就失效了； - 环方向维度：多边形外环必须逆时针（CCW），内环（孔洞）必须顺时针（CW）。这是OGC规范硬性要求，不是PostGIS自定规则。
ST_IsValidReason能告诉你“Ring Self-intersection”，但不会告诉你哪个环方向错了； - 拓扑维度：
ST_LineFromText('LINESTRING(0 0,1 1,2 2)')生成的是有效线，但ST_LineFromText('LINESTRING(0 0,1 1,0 0)')会报错“Invalid geometry: Self-intersection”，因为首尾重合构成的是环，不是线。
提示：别依赖
ST_MakeValid兜底。它能把无效几何“修好”，但可能把一个带孔洞的多边形强行合并成实心块，彻底破坏业务语义。正确的做法是：创建时就保证有效。
2.2 函数选型的底层逻辑：何时用WKT，何时用坐标数组，何时用现有几何构造？
PostGIS提供至少5种创建几何的方式，每种适用场景截然不同：
ST_GeomFromText(wkt, srid)：适合静态、已知坐标的简单图形，如测试数据、行政边界代码。优势是可读性强，缺点是字符串拼接易出错，且WKT不携带坐标系元数据（需显式传SRID）；ST_MakePoint(x,y[,z,m])+ST_SetSRID(geom,srid)：适合动态生成单点，如GPS轨迹实时入库。ST_MakePoint不校验SRID，ST_SetSRID只是打标签，不进行坐标变换；ST_PointFromText(wkt,srid)：与ST_GeomFromText功能重叠，但仅限点类型，对非点WKT会报错，属于“强类型防护”；ST_MakeLine(geom1, geom2)/ST_MakeLine(array[geom1,geom2,...])：适合由已有几何对象连接成线，如将一组POI点按时间序连成轨迹。注意：输入点必须共面（同一SRID），否则报错；ST_Polygon(geom, srid)：仅接受闭合线（LinearRing）作为输入，这是最易踩坑的点。很多人误用ST_Polygon(ST_MakeLine(...))，结果报错“Argument to Polygon must be a LinearRing”，因为ST_MakeLine生成的是LineString，不是LinearRing。
我实际项目中总结的选型口诀：“静态查WKT，动态算坐标，组合靠已有，复杂用构造”。比如做城市热力图，网格单元是固定矩形，用WKT最稳妥；做车辆轨迹分析，点来自传感器流，必须用ST_MakePoint动态生成；做行政区划合并，用ST_Union聚合已有几何；做复杂建筑轮廓（含多个孔洞），则必须先用ST_ExteriorRing和ST_InteriorRingN提取环，再用ST_MakePolygon组装。
2.3 SRID绑定策略：为什么“不设SRID”比“设错SRID”更危险？
新手常犯的错误是：CREATE TABLE places (id SERIAL, geom GEOMETRY(POINT))，然后插入ST_GeomFromText('POINT(116.4 39.9)')。表面看没问题，但后果严重：
- 空间索引（GIST）无法生效：
WHERE ST_DWithin(geom, ST_MakePoint(116.4,39.9), 1000)会全表扫描； - 距离计算失真：
ST_Distance返回单位是“度”，不是米； - 跨表关联失败：
JOIN两个不同SRID的表，PostGIS强制要求ST_Transform，否则报错。
正确做法是：建表时强制声明SRID，创建时显式绑定，且全程保持一致。例如北京地理坐标用WGS84（4326），投影坐标用Web Mercator（3857）。我见过最惨的案例：某物流系统用4326存坐标，但前端Leaflet用3857渲染，开发者没做ST_Transform，结果所有配送点偏移3公里——因为4326下1度≈111km，3857下1单位≈1米，单位错位直接放大误差。
3. 核心细节解析与实操要点：手把手拆解7类高频几何创建场景
3.1 单点创建：ST_MakePoint与ST_PointFromText的隐秘差异
单点看似最简单，却是埋雷最多的地方。我们对比两种常用方式：
-- 方式1：ST_MakePoint + ST_SetSRID SELECT ST_AsText(ST_SetSRID(ST_MakePoint(116.4, 39.9), 4326)); -- 返回：POINT(116.4 39.9) -- 方式2：ST_PointFromText SELECT ST_AsText(ST_PointFromText('POINT(116.4 39.9)', 4326)); -- 返回：POINT(116.4 39.9)表面结果一致，但底层行为天差地别：
ST_MakePoint(116.4,39.9)生成的是GEOMETRY类型，但无SRID属性，相当于“裸坐标”。ST_SetSRID只是给这个几何体打上4326标签，不进行任何坐标变换；ST_PointFromText('POINT(116.4 39.9)',4326)在解析WKT时，直接将SRID注入几何元数据，且对WKT格式做严格校验（如空格、括号）。
实操心得：
- 如果坐标来自程序变量（如Python的
x,y = get_gps()），必须用ST_MakePoint，因为WKT需要字符串拼接，易注入风险； - 如果坐标是固定字符串（如配置文件中的
"POINT(116.4 39.9)"），优先用ST_PointFromText，校验更严； - 绝对禁止：
ST_SetSRID(ST_GeomFromText('POINT(116.4 39.9)'),4326)。ST_GeomFromText默认返回GEOMETRY无SRID，ST_SetSRID虽能打标，但若WKT本身含Z值而字段未定义，后续操作仍会失败。
注意：
ST_MakePoint支持Z和M维度，但必须显式传参。ST_MakePoint(116.4,39.9,50)生成3D点，ST_MakePoint(116.4,39.9,50,100)生成4D点。若目标字段是GEOMETRY(POINTZ,4326)，漏传Z值会导致插入失败。
3.2 线串创建：ST_MakeLine的两种模式与环检测陷阱
线串（LineString）创建分“两点连线”和“多点序列”两种模式，行为完全不同：
-- 模式1：两点直接生成线 SELECT ST_AsText(ST_MakeLine( ST_MakePoint(0,0), ST_MakePoint(1,1) )); -- 返回：LINESTRING(0 0,1 1) -- 模式2：点数组生成线（注意：数组元素必须是geometry类型） SELECT ST_AsText(ST_MakeLine(ARRAY[ ST_MakePoint(0,0), ST_MakePoint(1,1), ST_MakePoint(2,0) ])); -- 返回：LINESTRING(0 0,1 1,2 0)关键陷阱在于环检测。当首尾点坐标完全相同时，ST_MakeLine不会自动转为多边形，而是生成一条“自相交”的无效线：
-- 错误示范：试图用ST_MakeLine创建环 SELECT ST_IsValid(ST_MakeLine(ARRAY[ ST_MakePoint(0,0), ST_MakePoint(1,0), ST_MakePoint(1,1), ST_MakePoint(0,1), ST_MakePoint(0,0) -- 首尾重合！ ])); -- 返回：false -- 原因：ST_MakeLine生成的是LineString，首尾重合即构成自相交（Self-intersection） -- 正确做法：先用ST_MakeLine生成线，再用ST_MakePolygon转环 WITH line AS ( SELECT ST_MakeLine(ARRAY[ ST_MakePoint(0,0), ST_MakePoint(1,0), ST_MakePoint(1,1), ST_MakePoint(0,1), ST_MakePoint(0,0) ]) AS geom ) SELECT ST_AsText(ST_MakePolygon(line.geom)) FROM line; -- 返回：POLYGON((0 0,1 0,1 1,0 1,0 0))实操心得：
ST_MakeLine的数组模式，数组长度必须≥2，否则报错“array must have at least two elements”；- 若点来自查询结果（如
SELECT geom FROM points ORDER BY time），必须用array_agg(geom ORDER BY time)聚合，且确保ORDER BY存在，否则数组顺序随机； - 对轨迹数据，务必用
ST_IsSimple检查线是否自相交，ST_Length计算前先验证有效性，避免NaN结果。
3.3 多边形创建：ST_MakePolygon的环顺序铁律与孔洞嵌套
多边形创建是PostGIS最易出错的环节，核心在于环方向（Ring Orientation）。OGC规定：外环必须逆时针（CCW），内环（孔洞）必须顺时针（CW）。违反则ST_IsValid返回false。
-- 正确：外环CCW，内环CW SELECT ST_AsText(ST_MakePolygon( ST_GeomFromText('LINESTRING(0 0,3 0,3 3,0 3,0 0)'), -- CCW ARRAY[ ST_GeomFromText('LINESTRING(1 1,2 1,2 2,1 2,1 1)') -- CW ] )); -- 错误：内环也用CCW（常见于手动画图习惯） SELECT ST_IsValidReason(ST_MakePolygon( ST_GeomFromText('LINESTRING(0 0,3 0,3 3,0 3,0 0)'), ARRAY[ ST_GeomFromText('LINESTRING(1 1,1 2,2 2,2 1,1 1)') -- CCW！错！ ] )); -- 返回：'Hole lies outside shell'验证环方向的技巧：用ST_Orientation函数（PostGIS 3.2+），返回CCW或CW：
SELECT ST_Orientation(ST_GeomFromText('LINESTRING(0 0,3 0,3 3,0 3,0 0)')) as outer_ring, ST_Orientation(ST_GeomFromText('LINESTRING(1 1,2 1,2 2,1 2,1 1)')) as inner_ring; -- outer_ring='CCW', inner_ring='CW'实操心得：
- 手动构造环时，用“左手规则”：沿环行走，区域在左手边即为CCW（外环），在右手边即为CW（内环）；
- 从GIS软件导出环坐标，QGIS默认CCW，但某些CAD软件可能用CW，导入前务必用
ST_ForceRHR（Force Right-Hand Rule）标准化； ST_MakePolygon的孔洞数组必须是非空数组，若无孔洞，第二个参数传NULL，不能传空数组ARRAY[]::geometry[]，否则报错。
3.4 复合几何创建：ST_Collect与ST_Union的本质区别
当需要将多个几何合并为一个集合（GeometryCollection）或单一几何（MultiGeometry）时，ST_Collect和ST_Union常被混用，但语义完全不同：
| 函数 | 输入类型 | 输出类型 | 是否改变拓扑 | 是否去重 | 典型用途 |
|---|---|---|---|---|---|
ST_Collect | 同类型或混合类型 | GEOMETRYCOLLECTION | 否 | 否 | 批量导入时暂存，保留原始部件 |
ST_Union | 同类型 | MULTI*或*（自动简化） | 是 | 是 | 合并相邻地块，消除公共边界 |
-- 示例：两个相邻矩形 WITH polys AS ( SELECT ST_GeomFromText('POLYGON((0 0,2 0,2 1,0 1,0 0))') as geom UNION ALL SELECT ST_GeomFromText('POLYGON((2 0,4 0,4 1,2 1,2 0))') as geom ) SELECT ST_AsText(ST_Collect(geom)) as collect_result, ST_AsText(ST_Union(geom)) as union_result FROM polys; -- collect_result: GEOMETRYCOLLECTION(POLYGON...,POLYGON...) -- union_result: MULTIPOLYGON(((0 0,2 0,2 1,0 1,0 0),(2 0,4 0,4 1,2 1,2 0))) -- 注意：union_result仍是MULTIPOLYGON，因为两矩形仅边界相接，未重叠若两矩形重叠，ST_Union会融合为单个多边形：
-- 重叠矩形 WITH polys AS ( SELECT ST_GeomFromText('POLYGON((0 0,2 0,2 2,0 2,0 0))') as geom UNION ALL SELECT ST_GeomFromText('POLYGON((1 1,3 1,3 3,1 3,1 1))') as geom ) SELECT ST_AsText(ST_Union(geom)) FROM polys; -- 返回：POLYGON((0 0,2 0,2 1,3 1,3 3,1 3,1 2,0 2,0 0))实操心得：
ST_Collect性能极高，O(1)时间复杂度，适合大数据量临时聚合；ST_Union是计算密集型操作，时间复杂度O(n²)，大数据量慎用，建议先用ST_ClusterDBSCAN聚类再ST_Union；ST_Collect不校验几何有效性，ST_Union会自动修复输入中的小瑕疵（如微小缝隙），但不可依赖此特性。
3.5 圆形与缓冲区：ST_Buffer的精度控制与性能权衡
PostGIS没有原生ST_MakeCircle函数，圆形必须通过ST_Buffer生成。其核心参数distance和quad_segs（四分圆段数）直接决定精度与性能：
-- 基础用法：以点为中心，半径1000米的圆 SELECT ST_AsText(ST_Buffer(ST_MakePoint(116.4,39.9), 1000)); -- 精确控制：quad_segs=32（默认值），生成128段圆弧（4*32） SELECT ST_NPoints(ST_Buffer(ST_MakePoint(116.4,39.9), 1000, 'quad_segs=32')); -- 返回129（闭合） -- 高精度：quad_segs=128，生成512段，更接近数学圆 SELECT ST_NPoints(ST_Buffer(ST_MakePoint(116.4,39.9), 1000, 'quad_segs=128')); -- 返回513 -- 低精度：quad_segs=4，生成16段，近似正方形 SELECT ST_AsText(ST_Buffer(ST_MakePoint(116.4,39.9), 1000, 'quad_segs=4')); -- 返回：POLYGON((...)) —— 实际是16边形quad_segs的取值逻辑：每个象限的线段数。quad_segs=1生成菱形（4边），quad_segs=2生成八边形（8边），以此类推。精度提升与顶点数呈线性关系，但渲染和存储开销呈指数增长。
实操心得：
- Web地图显示：
quad_segs=8（32边形）足够，人眼无法分辨与圆的差异； - 精密测量：
quad_segs=32（128边形），误差<0.01%； - 大批量生成（如10万点缓冲区）：
quad_segs=4可提速3倍，但需确认业务能否接受近似； ST_Buffer的join_style参数（round/mitre/bevel）影响角部处理，mitre在尖锐角可能产生极长延伸，生产环境建议round。
3.6 多部件几何：ST_Multi与ST_CollectionExtract的协同使用
当需要将单个几何（如MultiPoint）拆分为独立行，或合并为Multi类型时，ST_Multi和ST_CollectionExtract是黄金搭档：
-- 将多个点聚合为MultiPoint SELECT ST_AsText(ST_Multi(ST_Collect( ST_MakePoint(0,0), ST_MakePoint(1,1), ST_MakePoint(2,2) ))); -- 返回：MULTIPOINT(0 0,1 1,2 2) -- 将MultiPoint拆分为单点行（关键！用于逐点处理） SELECT (ST_Dump(geom)).path as path, ST_AsText((ST_Dump(geom)).geom) as point_wkt FROM ( SELECT ST_Multi(ST_Collect( ST_MakePoint(0,0), ST_MakePoint(1,1), ST_MakePoint(2,2) )) as geom ) t; -- 返回三行：path={1}, {2}, {3}，对应三个点ST_Dump是核心函数，它将集合几何“炸开”为行集，返回path（层级路径）和geom（子几何）。ST_CollectionExtract则用于从GeometryCollection中提取指定类型：
-- 从混合集合中提取所有多边形 SELECT ST_AsText(ST_CollectionExtract( ST_GeomFromText('GEOMETRYCOLLECTION(POINT(0 0),POLYGON((0 0,1 0,1 1,0 1,0 0)),LINESTRING(0 0,1 1))'), 3 -- 3=POLYGON )); -- 返回：POLYGON((0 0,1 0,1 1,0 1,0 0))实操心得：
ST_Dump是处理MultiGeometry的必经之路，几乎所有批量空间分析都需先ST_Dump；ST_CollectionExtract的type参数：1=POINT, 2=LINESTRING, 3=POLYGON, 4=CURVE, 5=POLYHEDRALSURFACE；ST_Multi对单部件几何是“无害包装”，但对已存在的Multi类型会报错，安全写法是ST_CollectionExtract(geom, type)后ST_Multi。
3.7 动态几何生成：基于SQL表达式的实时构造
最高阶的几何创建，是脱离静态坐标，用SQL表达式动态生成。典型场景：根据属性值生成不同尺寸的符号、按时间生成移动轨迹、按统计值生成热力网格。
-- 场景1：按人口密度生成不同半径的圆（假设density字段单位：人/平方公里） SELECT id, name, ST_Buffer( geom, SQRT(density * 10000 / PI()) -- 半径 = sqrt(面积/π)，面积=密度*10000m² ) as symbol_geom FROM cities; -- 场景2：生成网格（1km×1km）覆盖研究区域 WITH bounds AS ( SELECT ST_Extent(geom) as bbox FROM study_area ), grid AS ( SELECT ST_MakeEnvelope( ST_XMin(bbox) + x * 1000, ST_YMin(bbox) + y * 1000, ST_XMin(bbox) + (x+1) * 1000, ST_YMin(bbox) + (y+1) * 1000, ST_SRID(bbox) ) as cell FROM bounds, LATERAL ( SELECT x, y FROM generate_series(0, FLOOR((ST_XMax(bbox)-ST_XMin(bbox))/1000)::int) as x, generate_series(0, FLOOR((ST_YMax(bbox)-ST_YMin(bbox))/1000)::int) as y ) g ) SELECT row_number() OVER() as gid, ST_AsText(cell) as wkt FROM grid;实操心得：
- 动态生成必须配合
ST_SRID确保坐标系一致，ST_Extent返回的bbox无SRID，需显式ST_SetSRID； generate_series是生成网格的核心，但大数据量时（如1000×1000网格），LATERAL JOIN比CROSS JOIN更高效；SQRT等数学函数在ST_Buffer中可直接使用，无需提前计算，PostGIS会优化执行计划。
4. 实操过程与核心环节实现：一个完整城市设施点位管理系统的几何构建流程
4.1 需求分析：从Excel表格到空间数据库的端到端映射
假设我们接手一个城市设施管理系统，原始数据是Excel表格，包含列：id,name,type,lng,lat,altitude,status。目标是构建PostGIS表，支持空间查询、缓冲区分析、热力图渲染。整个流程分五步：
- 建表设计：明确几何类型、维度、SRID；
- 数据清洗：处理空值、异常坐标、重复记录；
- 几何创建：选择最优函数链路；
- 质量校验：有效性、拓扑、索引完整性；
- 性能优化：空间索引、分区、统计信息更新。
注意：跳过任何一步都可能导致后续分析崩溃。我曾见一个项目因未做第2步清洗，
lat列混入文本“N/A”，导致ST_MakePoint报错中断，全量导入失败。
4.2 建表与字段定义：维度、SRID、约束三位一体
-- 创建表：设施点位（含高程Z） CREATE TABLE facilities ( id SERIAL PRIMARY KEY, name VARCHAR(100) NOT NULL, type VARCHAR(50) NOT NULL, altitude NUMERIC(6,2), -- 单位：米 status VARCHAR(20) DEFAULT 'active', geom GEOMETRY(POINTZ,4326), -- 关键：POINTZ + 4326 created_at TIMESTAMP WITH TIME ZONE DEFAULT NOW() ); -- 添加空间约束：强制几何类型和SRID ALTER TABLE facilities ADD CONSTRAINT enforce_dims_geom CHECK (ST_NDims(geom) = 3); -- 必须3D ALTER TABLE facilities ADD CONSTRAINT enforce_srid_geom CHECK (ST_SRID(geom) = 4326); ALTER TABLE facilities ADD CONSTRAINT enforce_geotype_geom CHECK (GeometryType(geom) = 'POINT' OR geom IS NULL); -- 创建空间索引（GIST） CREATE INDEX idx_facilities_geom ON facilities USING GIST(geom);关键点解析：
GEOMETRY(POINTZ,4326)：声明几何类型为带Z值的点，SRID为WGS84。若未来需投影坐标，应另建字段geom_3857，而非修改此字段；- 三重
CHECK约束：ST_NDims确保Z维度存在，ST_SRID确保坐标系正确，GeometryType确保类型为点； idx_facilities_geom索引必须在数据插入前创建，否则全表扫描。
4.3 数据清洗与转换：用SQL完成ETL核心逻辑
假设Excel已导入临时表facilities_raw(id,name,type,lng,lat,altitude,status)。清洗步骤：
-- 步骤1：过滤空坐标和异常值（经纬度范围校验） WITH cleaned AS ( SELECT id, name, type, CASE WHEN lng BETWEEN -180 AND 180 THEN lng ELSE NULL END as lng, CASE WHEN lat BETWEEN -90 AND 90 THEN lat ELSE NULL END as lat, altitude, status FROM facilities_raw WHERE lng IS NOT NULL AND lat IS NOT NULL ) -- 步骤2：生成几何，处理Z值缺失（用0填充） INSERT INTO facilities (id, name, type, altitude, status, geom) SELECT id, name, type, altitude, status, ST_SetSRID( ST_MakePoint( COALESCE(lng, 0), COALESCE(lat, 0), COALESCE(altitude, 0) ), 4326 ) as geom FROM cleaned WHERE lng IS NOT NULL AND lat IS NOT NULL; -- 再次过滤，确保坐标有效实操心得：
COALESCE(lng,0)防止NULL导致ST_MakePoint返回NULL；ST_SetSRID必须在ST_MakePoint之后，顺序不可逆；- 清洗SQL应分步执行，每步后
SELECT COUNT(*)验证数据量，避免错误累积。
4.4 质量校验：五层验证确保几何“零缺陷”
导入后，必须执行以下校验（缺一不可）：
-- 1. 有效性校验：是否存在无效几何？ SELECT COUNT(*) FROM facilities WHERE NOT ST_IsValid(geom); -- 2. 维度校验：是否全部为3D？ SELECT COUNT(*) FROM facilities WHERE ST_NDims(geom) != 3; -- 3. SRID校验：是否全部为4326？ SELECT COUNT(*) FROM facilities WHERE ST_SRID(geom) != 4326; -- 4. 空间索引校验：索引是否生效？ EXPLAIN ANALYZE SELECT * FROM facilities WHERE ST_DWithin(geom, ST_MakePoint(116.4,39.9), 1000); -- 5. 业务逻辑校验：海拔是否在合理范围？（如-400~9000米） SELECT MIN(altitude), MAX(altitude) FROM facilities;若发现无效几何，用ST_IsValidReason定位：
SELECT id, name, ST_IsValidReason(geom) FROM facilities WHERE NOT ST_IsValid(geom); -- 常见原因：'Self-intersection'（自相交）、'Too few points in geometry'（点数不足）4.5 性能优化：空间索引、VACUUM、ANALYZE三剑客
数据就绪后，必须执行：
-- 更新统计信息（让查询规划器知道数据分布） ANALYZE facilities; -- 清理死元组（提高查询速度） VACUUM facilities; -- 重建空间索引（若数据量大，VACUUM后重建更优） DROP INDEX idx_facilities_geom; CREATE INDEX idx_facilities_geom ON facilities USING GIST(geom);提示：
VACUUM不是可选操作。PostgreSQL的MVCC机制会产生死元组，不清理会导致索引膨胀、查询变慢。我维护的一个千万级设施表，未定期VACUUM，空间索引大小从500MB涨到3GB，ST_DWithin查询从20ms升至2s。
5. 常见问题与排查技巧实录：12个真实踩坑场景与解决方案
5.1 “ERROR: Geometry has Z dimension but column does not” —— 维度错配
现象：插入ST_MakePoint(1,2,3)到GEOMETRY(POINT,4326)字段时报错。
根因：目标字段定义为2D，但几何体是3D。
解决：
- 方案1（推荐）：修改字段为
GEOMETRY(POINTZ,4326)； - 方案2：插入时丢弃Z值
ST_Force2D(ST_MakePoint(1,2,3))； - 方案3：用
ST_MakePoint(1,2)生成2D点。
5.2 “ERROR: Argument to Polygon must be a LinearRing” —— 环类型错误
现象：ST_MakePolygon(ST_MakeLine(...))报错。
根因：ST_MakeLine生成LineString，ST_MakePolygon要求LinearRing。
解决：
- 用
ST_AddPoint确保首尾重合，再ST_MakeLine，最后ST_LineToCurve？错！ - 正确：
ST_MakeLine后，用ST_ClosedCurve检查，再ST_ForceClosed（PostGIS 3.2+），或手动追加首坐标。
5.3ST_Buffer返回空几何 —— 距离为0或负数
现象：ST_Buffer(geom, 0)返回NULL或空。
根因：PostGIS 3.0+中，distance=0被视为无效操作。
解决：
- 用
ST_SnapToGrid(geom, 0.000001)替代，或 CASE WHEN distance=0 THEN geom ELSE ST_Buffer(geom,distance) END。
5.4ST_DWithin查询极慢 —— 空间索引未命中
现象：WHERE ST_DWithin(geom, point, 1000)全表扫描。
排查：
EXPLAIN看是否用到idx_facilities_geom；- 检查
point的SRID是否与geom一致（ST_SRID(point)=4326？）； - 检查
point是否为GEOMETRY类型（ST_MakePoint生成的是，字符串WKT需ST_GeomFromText）。
5.5 多边形显示为“空心” —— 孔洞方向错误
现象：QGIS中多边形显示为线框，无填充。
根因：内环方向为CC
