“查出仓库两公里内的客户”看似只需一个距离函数。如果经纬度被调换、投影坐标被标成角度，或把度数直接与 2000 比较，查询仍可能返回一张像样的列表。本文面向维护位置数据与智能问数接口的工程团队，限定于地表点到点的二维距离，不讨论路线规划或测绘级三维计算。
先定义距离问题，再接收坐标
首先写清业务问题中的起点和终点。仓库坐标代表建筑中心、装卸入口还是园区大门，客户坐标来自实地定位还是地址匹配，这些选择会影响两公里附近的记录。计算函数不能纠正坐标代表了错误业务对象的问题。本文建议在数据契约中保存对象 ID、位置含义、来源、采集时间与精度说明。
把服务半径写成地表直线距离还是道路行驶距离，也必须先决定。河流两岸的点可能相距很近，车辆却需绕行；直线距离适合候选筛选，不足以承诺配送时长。若业务需要道路距离，应接入有明确数据版本和道路限制的路线计算，不能给直线结果换一个字段名就交付。
对模型或厂区局部坐标，可参照三维模型的单位与坐标变换验收理解单位、原点和轴向。但局部米制坐标并不自动拥有地理参考；若缺少可靠控制点或已确认的转换关系，应保持独立，不能直接赋为 WGS 84 后参与全国范围查询。本文新增的是地理参考系与距离运算的边界，而非重复静态模型的矩阵变换。
GeoJSON 的顺序明确，但数值范围不是来源证明
RFC 7946 使用 WGS 84 地理坐标，位置数组的前两项按经度、纬度排列，单位为十进制度。[1] 因此接口契约应明确 lon、lat 或数组顺序，不接受一个未说明意义的 x、y 列表再凭经验猜测。其他接口是否采用相同顺序，需要分别核对，不能由文件中出现“坐标”二字推断。
导入时先检查类型、有限数值、缺失状态和业务覆盖范围。经度超出约定范围、纬度超出正负 90 度等明显错误应进入隔离区。更难发现的是两项都落在合法区间：构造点 (12, 45) 被交换为 (45, 12)，两者数值都合法，但位置完全不同。范围检查只能筛错，不能证明轴顺序正确。
建议用来源方确认的至少几个地理分散控制点，核对其地名、经纬顺序与地图位置；不要全部使用靠近原点或对称位置的样本。原始字段应保留，清洗后的 lon、lat 单独存储，并记录转换规则版本。发现整批偏移时先核对来源声明，不要逐条拖动点位去迎合底图。
缺失位置不能填成 (0, 0)，因为零坐标是有效数值。位置过旧、地址匹配到城市中心、只有低精度行政区位置，都应有单独状态。若没有精度依据，至少标记未知；增加小数位只改变显示，不增加定位可信度。
声明 SRID 与转换坐标是两种不同操作
PostGIS 的 ST_SetSRID 只设置参考系标识，不改变坐标数值。[2] 在已确认输入是 WGS 84 经度、纬度时，ST_SetSRID(ST_MakePoint(lon, lat), 4326) 可以建立带声明的点。关键前提是来源已确认；这个函数本身并未检测数据是否真的属于该参考系。
如果源数据是已知投影坐标，要先保留正确源 SRID，再用 ST_Transform 转到目标参考系；该函数实际改变坐标，并依赖可用的转换支持。[3] 不要将以米计的投影数值直接标为 4326，也不要把来源未知数据反复尝试不同 SRID，直到图上看起来合理。
建议将处理链分成四个持久字段或记录：原始值及来源、已确认源参考系、转换后的目标值、转换软件与规则版本。转换失败应保留错误和原始值，不自动退回原坐标却继续贴目标标签。涉及高精度或特殊区域时，还应由具备相应能力的人员确认适用转换及所需资源。
双向转换能检查部分实现问题，却不能单独证明参考系正确：错误地声明一个源参考系后，再按同一套错误声明来回转换，也可能回到原数值。因此需要独立控制点与目标环境位置核对，不能仅凭可逆就通过。
同一个距离函数，数据类型决定单位
PostGIS 中，ST_Distance 对 geometry 计算二维平面距离，结果采用该参考系单位；对 geography 默认计算椭球面测地距离并返回米。[4] 所以两个 SRID 4326 的 geometry 点相减，不能把结果当米。写成 ST_Distance(a, b) < 2000 而没有说明类型和单位，是一个应在评审中阻止的模糊接口。
对已确认的 WGS 84 地表点，可以选择 geography 表达米制半径；对受控区域的工程平面分析，也可以选择适合该区域的投影 geometry。二者不是速度与准确性的一句口号比较：需要结合数据覆盖、目标误差、可用运算和索引进行测试。仅因某投影单位写作米，不代表它在所有纬度都保持真实地面距离。
构造示例：仓库 A=(0, 0)，客户 B=(0.01, 0)，客户 C=(1, 0)，均假设为已确认的 WGS 84 经度、纬度。geometry 的平面角度距离分别为 0.01 和 1；与 2000 作比较，两者都会落入阈值，但这不是两公里判断。赤道附近上述经度差对应的地表距离分别约为 1.11 公里和 111 公里，因此两公里筛选应包含 B、排除 C。近似数用于解释量级，不作为实测结果。
可供目标环境验证的查询表达式是 ST_DWithin(a::geography, b::geography, 2000)。这里 a、b 已是正确构造的 SRID 4326 点，而不是任意原始列；2000 的单位为米。ST_DWithin 的 geography 版本使用米制阈值，并具备利用可用索引进行边界框筛选的机制。[5] 是否实际采用索引仍需检查执行计划，不能由函数名保证。
从原始导入到可查询数据的实施步骤
第一步建立暂存表。保留外部对象 ID、原始坐标文本、来源参考系、轴顺序、位置时间和批次号；解析错误与元数据缺失分别计数。不要在导入函数里同时猜单位、交换轴、转换位置并覆盖原值，否则故障后难以还原是哪一步造成偏移。
第二步按来源分组验证。范围检查通过后，用控制点与业务覆盖区检查位置，再按已确认规则生成标准地理点。未知参考系进入待核查状态，不参与默认距离统计。合格数、隔离数与输入总数应能对账；少返回客户时，查询结果应说明排除了多少未知位置。
第三步选择服务端统一的数据类型。若主要需求是跨地区的点半径查询，可以物化 geography 列并建立对应 GiST 索引，避免每次查询都在临时表达式上进行未规划的转换。若保留 geometry 并在查询中转换或强制类型转换，应按实际表达式评估索引方案。用代表性数据量验证计划与耗时，不在没有测量时承诺加速倍数。
第四步封装参数边界。要求半径非负、单位固定、中心点来源明确，并为允许的最大范围设置业务约束。返回对象 ID、距离值及单位、位置版本和计算口径。若结果需要分页，用距离与稳定 ID 共同排序，并保持一致的数据版本，避免更新中的位置使同一对象跨页跳动。
第五步发布前留存对照。保存原始批次、标准化记录、控制点结果、查询参数和参考结果。转换规则改变后生成新候选版本，比较对象数量、区域覆盖和边界名单，再切换默认查询。若失败，回到上一版本；不要直接用新版坐标覆盖唯一可恢复的数据。
用边界与反例验收，不只看地图上有没有点
至少覆盖六类用例：同一点距离为零；上述 B、C 的两公里分类；轴顺序交换但范围合法；缺失坐标不变成零点；源 SRID 错误或缺失被隔离；距离阈值附近的内外点。分别记录预期状态和允许误差，不能只看查询是否返回成功。
阈值附近的控制点应由独立参考或已验证的计算链构造，并测试略小于、等于和略大于约定半径的情况。先明确边界是否包含，再区分数值误差和定位误差。若设备定位误差已达几十米，显示六位小数并不能使两公里边界具有厘米级可信度。
对接近边界、精度不足的记录，本文建议输出“确定在内、确定在外、需复核”的业务状态，容差带由定位资料和业务损失共同确定，不能一律固定为某个米数。这是应用层建议，不是 PostGIS 自动提供的置信区间；没有误差模型时不要虚构概率。
最后加入跨越日期变更线的样本以及目标业务覆盖的高纬度样本，检查是否错误地用经度差直接判断远近。本文流程不证明道路可达、楼层接近或设备真实在场；也未运行 PostGIS 性能实验。官方函数文档查阅于 2026 年 11 月 8 日，部署验收应冻结 PostGIS、PROJ 与数据版本。