空间连接问题

Question

我是新来的sql，并试图使用它来加快空间分析的~120万旅行从一个csv，包含lat和lon的拾取和下降点。我想用简单的英语说：

• 将从感兴趣区域开始的所有行程（作为shapefile加载到我的数据库中）选择到一个表中

我正在使用的代码如下（如果有帮助，NTA或邻域制表区域，是我想在此操作结束时显示在我的表中的邻域名称）：

    --Select all trips that end in project area
SELECT *
INTO end_PA
FROM trips, projarea
WHERE ST_Intersects(trips.dropoff, projarea.geom);

--for trips that end in project area - index by NTA of pick up point

ALTER TABLE end_PA ADD COLUMN GID SERIAL;

CREATE TABLE points_ct_end AS
SELECT nyct2010.ntacode as ct_nta, end_PA.gid as point_id 
from nyct2010, end_PA WHERE ST_Intersects(nyct2010.geom , end_PA.pickup);


--Count most common NTA
--return count for each NAT as a csv
copy(
select count(ct_nta) from points_ct_end
group by ct_nta
order by count desc)
to 'C://TaxiData//Analysis//Trips_Arriving_LM.csv' DELIMITER ',' CSV HEADER;

但是，我从一开始就遇到了问题-ST_交叉点不会返回感兴趣区域内的任何点！

到目前为止，我尝试过的故障排除解决方案：

我的第一个想法是这些点不在正确的SRID中。当我创建下拉点时，我将SRID设置为4326。我尝试使用ST_SetSRID将两个数据集的投影更改为4326，并在ArcMap中手动将形状文件重新投影为4326，但都不起作用。

我从Arc Map中的旅行数据集中绘制了一个小样本点，以确保它们被正确地投影并与ProjArea形状文件重叠。他们是。

我将这个创建的多点shapefile导入到我的geo数据库，以测试它是否适用于ST_Intersects。没有。

我试着在内部使用ST_。这引发了错误消息：

错误：中的函数st_（字符变化、几何体）不存在
…
提示：没有与给定名称和参数类型匹配的函数。您可能需要添加显式类型转换。

我正在使用大SQL和postgres谢谢！！

Answer 1

我的第一个想法是，这些点不在正确的SRID中。创建“衰减”点时，我将SRID设置为4326。我尝试使用ST_SetGrid将两个数据集的投影都更改为4326，并在ArcMap中手动将ShapeFile重新投影到4326，但都不起作用。

ST_SetGrid不会更改投影（重投影）。它只是改变了内部表示。如果以前的SRID与输入数据匹配，这可能会完全搞糟一切。您可能想要ST\u Transform（）。

这里没有足够的信息来解决这个问题。然而，我们可以回答这个...

ERROR: function st_within(character varying, geometry) does not exist

这仅仅意味着第一个参数不是几何。当然，我们对此无能为力，因为我们没有您尝试过的ST_Within（）查询。

您的ST_Intersects（）语法看起来是正确的。但是，根本没有足够的信息来提供帮助。显示一些模式和示例数据。