直接在应用层调用ST_Contains易出错，因需确保SRID一致、坐标系匹配且索引生效；应通过视图封装ST_Intersects与自动坐标转换，并在原始几何列上建GIST索引以保障性能。
为什么直接在应用层调用ST_Contains会出错
多数GIS业务不是“查一个点在不在多边形里”就完事，而是要查“某区域内所有带属性的设施点”，同时还要关联业务表（比如facilities、districts）。如果每次都在应用代码里拼ST_Contains或ST_Intersects，容易漏掉SRID校验、坐标系不一致、索引未命中等问题。PostGIS里ST_Contains(geom1, geom2)要求两个几何对象SRID必须相同，否则返回false或报错Operation on two geometries with different SRIDs。
- 应用层SQL硬编码空间函数，导致同一逻辑在多个服务中重复实现，改坐标系时要全局搜
ST_Transform - 没加
USING INDEX提示或没建GIST索引，10万级面数据+点查询可能秒变全表扫描 -
ST_DWithin写成ST_DWithin(geom, 'POINT(...)', 500)但没指定单位（默认是度），实际想查500米却扫了半个城市
用视图封装ST_Intersects + 坐标系自动对齐
把空间判断和坐标转换收进视图，让业务查询只关心“要什么”，不操心“怎么算”。比如统一用EPSG:4326输入，内部自动转目标投影再计算：
CREATE OR REPLACE VIEW facilities_in_districts AS
SELECT f.id, f.name, f.type,
d.name AS district_name,
d.code AS district_code
FROM facilities f
JOIN districts d
ON ST_Intersects(
ST_Transform(f.geom, 4326),
ST_Transform(d.geom, 4326)
);
注意这里ST_Transform顺序不能反——facilities.geom通常存的是EPSG:3857（Web墨卡托），districts.geom可能是EPSG:4326，强行用ST_Intersects(f.geom, d.geom)大概率不走索引甚至报错。视图里显式转成同一SRID，既安全又可复用。
- 视图定义里别用
*，明确列出字段，避免下游因源表加字段导致视图行为突变 - 如果
facilities表本身已存EPSG:4326，就别套ST_Transform，徒增开销 - PostgreSQL 12+ 支持物化视图，但
ST_Intersects结果无法被物化（非immutable函数），所以只能用普通视图
给空间视图加GIST索引才能真正提速
视图本身不存数据，也没法直接建索引。真正起作用的是底层表的GIST索引。如果只在facilities.geom上建了GIST，但查询条件里用了ST_Transform(f.geom, 4326)，优化器大概率弃用索引——因为函数表达式改变了原始列值。
- 优先在原始存储的
SRID列上建索引：CREATE INDEX idx_facilities_geom ON facilities USING GIST (geom); - 如果业务90%查询都基于
EPSG:4326，考虑在表里加一列geom_4326并建索引，用触发器或写入时同步更新，避免每次查都ST_Transform - 用
EXPLAIN ANALYZE确认是否走了索引：看到Index Scan using idx_... on facilities才算生效；若出现Seq Scan，基本是函数包裹或SRID不匹配导致
视图里用ST_DWithin替代手动算距离
查“某点500米内所有加油站”，别自己写ST_Distance(f.geom, ST_Point(116.4, 39.9)) ，这没法用索引。必须用<code>ST_DWithin，它能触发GIST索引的边界框快速过滤：
CREATE OR REPLACE VIEW nearby_gas_stations AS
SELECT id, name, address,
ST_Distance(geom, ST_SetSRID(ST_Point(116.4, 39.9), 4326)) AS distance_m
FROM facilities
WHERE type = 'gas_station'
AND ST_DWithin(
geom,
ST_SetSRID(ST_Point(116.4, 39.9), 4326),
500
);
关键点：ST_DWithin第三个参数单位取决于前两参数的SRID。如果geom是EPSG:3857，500就是米；如果是EPSG:4326，500就是度（≈55km），结果完全失真。所以务必先确认几何列的实际SRID，再决定是否加ST_Transform。
- 参数顺序不能颠倒：
ST_DWithin(a, b, dist)中a必须是已建GIST索引的列，否则索引失效 - 如果点坐标来自用户输入，用
ST_SetSRID(ST_MakePoint(lon, lat), 4326)比ST_GeomFromText('POINT(...)', 4326)更快，少一次文本解析 - 视图里别直接写死坐标（如
116.4, 39.9），应改为参数化视图（用pg_temp表或应用层传参），否则无法复用
SRID一致性检查和索引是否真被用上——哪怕视图写得再干净，底层没索引或ST_Transform位置错了，查询照样慢十倍。 
