在多边形中查找点

Question

我有一个数据库，其中一个表（ca_licenses）是商业地址，另一个表是市议会区多边形的公共模式（ca_la_la_council）。

我运行此查询是为了将议会表中的地区号码放在企业地址表中：

UPDATE poursafe.ca_licenses dst
SET prem_council = div_city
FROM public.ca_la_la_council src
WHERE st_within(dst.geom,src.geom) AND dst.geom&&src.geom;

我的问题是我总是得不到任何结果。两个几何列都是几何类型，SRID是4326。

#create poursafe.ca_licenses table

CREATE TABLE ca_licenses (
id INTEGER,
license VARCHAR,
lat DOUBLE PRECISION,
lon DOUBLE PRECISION,
geom GEOMETRY,
prem_council VARCHAR
);

#create public.ca_la_la_council table

CREATE TABLE ca_la_la_council (
gid INTEGER,
div_city VARCHAR,
shape_leng NUMERIC,
shape_area NUMERIC,
geom GEOMETRY
);

示例ca_licenses数据:https://raw . githubusercontent . com/pour safe/data _ samples/master/ca _ licenses . CSV

样本ca_la_la_council数据:https://raw . githubusercontent . com/pour safe/data _ samples/master/ca _ la _ la _ Council . CSV

我需要用div_city填充ca_licenses表列prem_council。为此，我运行上面的查询，查找ca_licenses点位于哪个议会区内。

我为public.ca_la_la_district表创建了一个spatial_ref_sys表，但没有ca_licenses表的spatial_ref_sys表。这张桌子也需要一个吗？我相信两个几何柱的 SRID 都设置为 4326。

我在Centos 7设置上使用postgis 2.4.4和postgres 9.2。

Answer 1

第一

我想推荐你去查一下QGIS

在linux机器上使用wine效果很好

第二

将多边形文件加载到qgis后，代码没有任何问题，这是shapefile结构

然后我加载了点文件

问题是它们之间没有任何交集，原因可能是点的错误构造。

（粉红色的小点是ca_licenses.csv点）

根据我的经验，将多边形与点相交的最佳方法是使用ST_INTERSECT

select * from ca_licenses join ca_la_la_council on st_intersects(ca_licenses.geom,ca_la_la_council.geom)