PostGIS ST_DWithin查询不符合预期问题排查求助

我在Station模型中添加了location_point字段，定义如下：

location_point: {
    type: Sequelize.GEOMETRY('POINT'),
},

插入Station新记录时，我会将经度和纬度转换为POINT类型存储到location_point字段中：

var location_point = {
    type: "Point",
    coordinates: [
         location_longitude,
         location_latitude,
    ],
};

现在我需要查找指定点fixedPoint（定义如下）50米范围内的所有Station，使用Sequelize针对PostgreSQL编写了如下查询：

const fixedPoint = {
      type: "Point",
      coordinates: [end_location_longitude, end_location_latitude], 
};
const distanceInMeters = 50;

const stationsWithinRadius = await Station.findAll({
            where: db.where(
                db.fn(
                    "ST_DWithin",
                    db.col("location_point"),
                    db.fn(
                        "ST_SetSRID",
                        db.fn(
                            "ST_GeomFromText",
                            `POINT(${fixedPoint.coordinates[0]} ${fixedPoint.coordinates[1]})`
                        ),
                        4326
                    ),
                    distanceInMeters,
                    true // Use meters as the unit of distance
                ),
                true
            ),
        });

我已在stations表中添加了多个距离fixedPoint不足10米的记录，但执行上述查询时，却返回了超出50米范围的记录。当前使用的SRID为4326，技术栈包含PostgreSQL、PostGIS和Sequelize，请问可能是什么问题？

核心问题是SRID 4326属于经纬度地理坐标系，单位是度而非米。你直接传入50作为距离参数时，PostGIS会把它当成50度来计算——这个范围覆盖数千公里，自然会返回大量超出预期的结果。

解决步骤如下：

• 转换坐标系到支持米单位的投影坐标系：比如全球通用的EPSG:3857（Web墨卡托），或者针对特定区域的高精度投影坐标系（如国内常用的EPSG:4523）。
• 修改查询逻辑，用ST_Transform将数据库中的坐标和查询点统一转换为目标投影坐标系后，再执行ST_DWithin判断。

修改后的Sequelize查询示例：

const fixedPoint = {
      type: "Point",
      coordinates: [end_location_longitude, end_location_latitude], 
};
const distanceInMeters = 50;

const stationsWithinRadius = await Station.findAll({
    where: db.where(
        db.fn(
            "ST_DWithin",
            db.fn("ST_Transform", db.col("location_point"), 3857),
            db.fn(
                "ST_Transform",
                db.fn(
                    "ST_SetSRID",
                    db.fn(
                        "ST_GeomFromText",
                        `POINT(${fixedPoint.coordinates[0]} ${fixedPoint.coordinates[1]})`
                    ),
                    4326
                ),
                3857
            ),
            distanceInMeters
        ),
        true
    ),
});

额外注意事项：

• 如果数据集中在特定区域，优先使用该区域对应的本地投影坐标系，距离计算精度比3857更高。
• 可通过PostGIS的ST_SRID(location_point)函数，验证数据库中location_point字段是否正确设置了SRID 4326。

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