Postgresql Postgis SQL复杂联接（但是不一定与GIS相关）

Question

抱歉，标题欠佳，但很难描述...

假设我有两个表（cad和cad_polygon）...

cad和cad_polygon共享与另一个相关的同一列（cad_pid）...

cad具有以下列：cad_pid，jrsdctn_id

而cad_polygon具有以下列：cad_pid，ogc_fid，wkb_geometry

现在，我正在工作的以下查询（半天尝试）基于一个较长的经纬度坐标从小得多的多边形子集中根据度数选择一个宗地，它查找距给定的较长的经纬度坐标以米为单位的距离-ord，然后仅显示质心在长而纬坐标500m以内的多边形。

SELECT SUBQUERY.cad_pid, SUBQUERY.ogc_fid, SUBQUERY.dist_meters,
    SUBQUERY.wkb_geometry FROM (
    SELECT cad_pid, ogc_fid,
    CAST(ST_Distance_Sphere(
        ST_Centroid(wkb_geometry),
            ST_GeomFromText(
                'POINT(00.0000 -00.0000)',
            900914)
        ) AS numeric
    ) AS dist_meters, wkb_geometry
    FROM cad_polygon
    WHERE ST_DWithin(
        ST_Centroid(wkb_geometry),
        ST_GeomFromText(
            'POINT(00.0000 -00.0000)',
        900914),
    0.01)
    ORDER BY dist_meters ASC
) AS SUBQUERY
WHERE SUBQUERY.dist_meters < 500;

我想添加到此列表，并使用该查询弹出的列表，加入其他表（cad），以便为每个结果行提供额外的列“ jrsdctn_id” ...即：

样本数据为：cad：

cad_pid | jrsdctn_id
0001    | abc123
0002    | def456
0003    | dhk778
0004    | dsk730


cad_polygon:

cad_pid | ogc_fid | wkb_geometry
0001    | ht0101  | 67686687601010000200063D7987FF15ASD1518541DAW
0002    | hz4561  | 435453457601010000200063D7987FF15ASDFW4GF8DE4
0003    | yv0301  | 2626WD687601010000200063D7987FF15ASD1WE851D4D
0004    | vt9701  | D484DW4D8441D8W1C684V63D7987FF15ASD1D7DW4848D

预期成绩：

cad_pid | jtsdctn_id | ogc_fid | dist_meters | wkb_geometry
0002    | def456     | hz4561  | 192.769     | 43545...
0004    | dsk730     | vt9701  | 342.548     | D484D...

如果有一些sql向导可以提供帮助，那就太好了！

Answer 1

有点晚了，当然，关于JOIN的公认答案是绝对正确的，但这实际上与GIS密切相关，并且它的意识最终使您更轻松了：

您似乎正在使用自定义CRS或ogr2ogr （或任何GDAL / OGR函数）在PostGIS的space_ref_sys表中找不到匹配的SRID /投影； 任何使用LonLat作为球面/类弧形代数输入的PostGIS函数，都将始终假定您的坐标在EPSG：4326（WGS84）中。

如果您的LonLat与WGS84不完全匹配，结果将不正确 ！

现在，PostGIS也具有地理类型，它再次假设EPSG：4326坐标，如果与这些功能一起使用，将隐式使用米作为单位，如果与默认的ST_Distance参数一起使用，将在WGS84椭球体上进行计算（更精确，但比use_spheroid := false稍慢，后者将根据球体来计算距离）。

考虑到这一点，您的查询可以表示为：

WITH
  pt AS (
    SELECT ST_Transform(ST_SetSRID(ST_MakePoint(0, 0), 900914), 4326)::geography AS geog
  ),

  ctr AS (
    SELECT *,
           ST_Transform(ST_Centroid(wkt_geometry), 4326)::geography AS geog
    FROM cad_polygon
  )

SELECT ctr.cad_pid,
       cad.jtsdctn_id,
       ctr.ogr_fid,
       ST_Distance(ctr.geog, pt.geog) AS distance_meter,
       ctr.wkt_geometry
FROM ctr
JOIN cad
  ON ctr.cad_pid = cad.cad_pid
WHERE ST_DWithin(ctr.geog, pt.geog, 500)
ORDER BY distance_meter ASC;

请注意，使用CTE可以避免为每个已处理的行进行转换/广播，并使其变得更加结构化。

我只是不能放手...

Answer 2

您可以使用JOIN

    SELECT SUBQUERY.cad_pid, SUBQUERY.ogc_fid, SUBQUERY.dist_meters,
        SUBQUERY.wkb_geometry, SUBQUERY.jrsdctn_id FROM (
        SELECT cad_pid, ogc_fid,
        CAST(ST_Distance_Sphere(
            ST_Centroid(wkb_geometry),
                ST_GeomFromText(
                    'POINT(00.0000 -00.0000)',
                900914)
            ) AS numeric
        ) AS dist_meters, wkb_geometry, cad.jrsdctn_id
        FROM cad_polygon
        INNER JOIN cad on cad.cad_pid = cad_polygon.cad_pid

        WHERE ST_DWithin(
            ST_Centroid(wkb_geometry),
            ST_GeomFromText(
                'POINT(00.0000 -00.0000)',
            900914),
        0.01)
        ORDER BY dist_meters ASC
    ) AS SUBQUERY
    WHERE SUBQUERY.dist_meters < 500;

或者最好在外部添加连接

SELECT SUBQUERY.cad_pid
        , SUBQUERY.ogc_fid
        , SUBQUERY.dist_meters,
        SUBQUERY.wkb_geometry, cad.jrsdctn_id 
    FROM (
        SELECT cad_polygon.cad_pid, ogc_fid,
        CAST(ST_Distance_Sphere(
            ST_Centroid(wkb_geometry),
                ST_GeomFromText(
                    'POINT(00.0000 -00.0000)',
                900914)
            ) AS numeric
        ) AS dist_meters
        , wkb_geometry
        FROM cad_polygon
        WHERE ST_DWithin(
            ST_Centroid(wkb_geometry),
            ST_GeomFromText(
                'POINT(00.0000 -00.0000)',
            900914),
        0.01)
        ORDER BY dist_meters ASC
    ) AS SUBQUERY
    INNER JOIN cad on cad.cad_pid = SUBQUERY.cad_pid
    WHERE SUBQUERY.dist_meters < 500;