问题标签 [pysal]

我正在寻找一种有效的方法来找到给定多边形的一阶邻居。我的数据是shapefile格式。

我的第一个想法是计算多边形质心的 x 和 y 坐标，以便找到邻居的质心。

这种方法可能对常规网格有效，但就我而言，我需要找到一个“更通用”的解决方案。该图显示了问题。考虑黄色多边形有裁判。邻居的多边形是灰色多边形。使用 centroids-neighbors 方法，清晰的蓝色多边形被认为是邻居，但它与黄色多边形没有共同的边。

最近从有效地找到 200k 多边形的一阶邻居修改的解决方案如下：

使用此解决方案，我有一个以顶点坐标为键的字典和一个块组 ID 列表，这些 ID 在该坐标处具有一个顶点作为值。

Pysal 的新手，刚刚开始处理加载、阅读等。我有一些十进制度 (WGS84) 的测试点数据，并且想要使用 get_point_dist 函数。运行它后，我得到了正确的响应，但我想要米/公里等的距离而不是 DD。

任何想法如何做到这一点，或者答案是“投影输入数据”？

我用这个statsmodels包来估计我的 OLS 回归。现在我想要Breusch Pagan test。我使用该pysal包进行此测试，但此函数返回错误：

返回错误：

AttributeError：“OLSResults”对象没有属性“u”

我应该怎么办？

我想要给定多边形中相交点的结果，但出现错误。

我的代码是：

但我收到错误消息：

“PurePyShpWrapper”对象没有属性“bounding_box”

我正在尝试将.dbf文件转换为pandas.DataFrame. 但是，我从以下代码行中收到错误消息：

文件“/Users/UserName/anaconda/lib/python2.7/site-packages/pysal/core/IOHandlers/pyDbfIO.py”，第 77 行，在init assert terminator == '\r' AssertionError

这是我的代码：

我有一个文本文件，其中包含以下形式的坐标：

我希望能够遍历每个坐标以检查它在哪个多边形中（shapefile 中的英国县），但是我不确定如何标记数字，以便我可以有一个代码...

目前它们是字符串，但我希望每条线由两个点组成，以便程序可以尝试将它们定位在多边形中。

我正在尝试绘制我拥有的一组坐标所在的县。我希望代码遍历每个坐标，检查它所在的县，然后打印相应的县。但是，当前在空闲状态下运行时，光标只是闪烁，似乎没有任何反应。任何有关我的代码可能出错的地方的帮助将不胜感激。我不确定代码是否正确检查县或是否可以从 shapefile 中获取县的名称？

我有一个二维值数组，我正在尝试分析空间相关性。为了在 python 中计算像Moran's I这样的 2D 自相关，pysal 提供了一个实现。1）如何将我的 2D 数据转换为 pysal 预期的 1D 数组？2）如何构造一个基于距离的权重数组 w（输入点数组在核距离函数中的含义是什么？）？

我正在尝试使用 PySAL（具有 HH/HL/LH/LL 象限的那个）创建一个 Moran 散点图，并且认为我已经到了那里，但想检查我的理解/解释/代码。下面的代码使用内置的 North Carolina SIDS 数据集和行标准化。

这产生了一些看似合理的东西，但我认为我已经认为自己陷入了困境，因为我还没有实际计算出 Moran 的 I （通过ps.Moran_Local(...)），这被称为 Moran 散点图......

首先让我说，除了 stata .do 文件和 R 中的一些 MCMC 之外，自从 VB 中的 AOLwarez 组应用程序以来，我就没有编码过；所以我很抱歉在这里和你们一起混日子。

我正在写一篇环境正义论文，使用人口普查区块组级别的人口统计和暴露数据。因为街区组的规模相对较小，污染源和居住在一个街区组中的人们很容易影响至少一阶邻居。我本来是个笨蛋，只是通过 FIPS 代码汇总备份，但这只是糟糕的数学。

得到了我想要的 ACS 年份的形状文件，首先尝试了 arcGIS，但我没有得到任何结果。然后我阅读了 Pysal 并安装了它

导入的形状文件对所有 216,000 个块组运行（12 小时）皇后邻居分析

在[52]中：w.histogram Out[52]：

[(0, 87), (1, 709), (2, 3634), (3, 16627), (4, 48736), (5, 56952), (6, 42848), (7, 24878), ( 8, 12646), (9, 6294), (10, 3040), (11, 1515), (12, 759), (13, 432), (14, 233), (15, 128), (16, 85), (17, 44), (18, 34), (19, 20), (20, 21), (21, 13), (22, 8), (23, 7), (24, 6) , (25, 1), (26, 3), (27, 1), (28, 2), (29, 1), (30, 2), (31, 1), (32, 0), ( 33, 2), (34, 0), (35, 1), (36, 1), (37, 1), (38, 0), (39, 0), (40, 0), (41, 0), (42, 0), (43, 0), (44, 0), (45, 0), (46, 1), (47, 0), (48, 0), (49, 0) , (50, 0), (51, 0), (52, 0), (53, 0), (54, 0), (55, 0), (56, 0), (57, 0), ( 58, 0), (59, 0), (60, 0), (61, 1)]

我需要的是一个 .csv（或者老实说，如果我在某处复制/粘贴它，任何事情都会做），它通过 FIPS 枚举每个块组（这应该是 ACS shapefile 用于 ID 的内容）及其邻居列表。

如果我能拿到清单，我可以把它移到我更舒服的环境中。昨晚我坐在那里玩了几个小时，并且可以在 numpy.savetext 上得到一些破解来工作，但它只有一列，并且数字以科学计数法存储，因为 FIPS 代码是 12 位数字。有一次它告诉我元组超出范围，我认为这是我得到的最接近的

我事先相当广泛地搜索了数据本身，或者我保证我不会在这里浪费你的时间。

谢谢你，戴夫