使用缓冲区和变换投影，PostGIS

Question

我正在创建一个缓冲区并尝试将我的 SRID 转换为 geog，然后再转换回原始投影和 geom。 我的点层有我的 WKT 约定：

POINT (-73.73891083685496 40.59593727361931)

CREATE TABLE intrsct_buff AS ( 
  SELECT
    ST_Transform(ST_Buffer(ST_Transform(pt.geom,2263), 60),4326) as geom, pt.count
  FROM public.intrsct_pts as pt
);

但是当我使用这段代码时，我的新缓冲表有这个 WKT：

POLYGON ((-77.51963927750339 40.11249512168173..

我是否为地理使用了错误的 SRID？ WGS84 = 4326。它需要的投影是 2263，因为这是我的点表和数据库中所有其他表的 SRID。

Answer 1

对于这种情况，请拿一张纸和一支笔，写下您需要（输入）和您期望（输出）的每个步骤。

几乎每一个假设都有一个错误。

1) Geography 与 Geometry 不同。虽然两者都使用 WGS84 4326 CRS，但前者以米为单位，而后者以度为单位。前者使用大圆弧，而后者使用直线连接两点（距离计算，相交等在大距离上影响很大）

2) 您说您的数据在 CRS 2263（纽约长岛 (ftUS)）中，具有以下界限：

WGS84 边界：-74.2700、40.4700、-71.7500、41.3100
预计范围：909126.0155、110626.2880、1610215.3590、424498.0529

您显示的坐标不在此 CRS 中，但更有可能在 WGS84 中。

让我们验证这个假设：如果我们（错误地）假设您提供的坐标是 2263，如果我们将它们投影到 4326 会发生什么？

select st_asText(st_transform(st_geomFromEWKT('SRID=2263;POINT (-73.73891083685496 40.59593727361931)'),4326));
                st_astext
------------------------------------------
 POINT(-77.5198536012811 40.112488506058)

Bingo，这就是你得到的结果。
==> 首先，要么修复你正在使用的 CRS，要么修复坐标。其中一个是错误的。

但即使您确实有正确的 CRS，查询也不会按预期工作。

您说您想要 1) 将点更改为地理，2) 将其缓冲 60 米，3) 转换回原始 CRS。

您所做的是 1) 将点更改为 2263（它已经在 CRS 中，它没有效果）2) 缓冲 60 英尺 3) 转换为 WGS84（几何）

为了实现你所描述的，宁愿是：

SELECT ST_Transform(
     ST_Buffer(
      ST_Transform(pt.geom,4326)::geography, -- change to 4326 and cast to geography
      60), -- buffer 60 meters
    2263) as geom -- change back to original CRS