问题标签 [pysal]

我正在使用 Pysal 来可视化地理空间数据。我想在个人（网络）之间绘制段，但我不知道如何绘制我的 LineSegment 形状列表（lc）。那么如何在我的地图中显示这些 LineSegment 呢？（下面是绘图代码）

谢谢 ！

我正试图掌握pysal。假设我有一个这样创建的棋盘：

现在我删除最后一行和最后一列以匹配权重矩阵中可用的维度pysal：

现在让我们使用Queen 邻接规则（又名Moore 邻域）来使用 Join Count Statistics（我的值cb是0or ）：1

现在，结果：

我的问题：

1. 为什么我没有得到不同连接的单个值，而是一个数组？此外，数组的每个值似乎都指向一个矩阵单元，但我没有得到总和。
2. 20.0就连接总数而言，是4+4+12。考虑到矩阵的大小和结构，我预计会有更多的连接（更改）。为什么我得到的数字与预期的相差甚远？

我需要创建一个具有一定程度空间自相关的随机 1000x1000 矩阵。我找到了一个 python 库pysal，但是该过程需要创建一个 1000000x1000000 矩阵，这会使程序崩溃。有没有办法解决这个问题？稀疏矩阵有意义吗？任何有替代解决方案的人都将不胜感激。

对于 Uni 的作业，我们必须使用 PySal 来计算 Moran's I。我找到了一个示例代码，并且正在使用该代码，但我收到错误消息“'file' object has no attribute 'by_col'”

任何想法将不胜感激。数据集是纽约的一个文本文件，“T0_MINC”属性是包含不同地区家庭收入中位数的列。脚本附在下面

我正在使用的示例脚本位于此链接中的“Moran”下：http: //pysal.readthedocs.io/en/latest/users/tutorials/autocorrelation.html#moran-si

我正在尝试spgwr使用 PySAL 重新创建一个 R 笔记本。使用 R，可以直接将局部系数提取到 DataFrame 中，如下所示：（CSV 可在此处获得）

R

有没有一种简单的方法可以使用结果来计算这些系数ps.spreg.OLS？

Python

但是，我不确定如何使用其中包含的数据fit来计算局部系数。

我正在尝试查找 1 公里半径内的所有最近邻居。这是我构建树并搜索最近点的脚本，

从我在 pysal 页面中读到的内容，它说 -

kd-tree 建立在 scipy 中的 kd-tree 功能之上。如果使用 scipy 0.12 或更高版本，则使用 scipy.spatial.cKDTree，否则使用 scipy.spatial.KDTree。

就我而言，它应该使用 cKDTree。这对于示例数据集来说工作得很好，但是由于tree.query_ball_point返回索引列表作为结果。每个列表将包含 100 个元素。对于我的数据点（200 万条记录），这会变得越来越大，并且在某个点之后由于内存问题而停止。关于如何解决这个问题的任何想法？

我正在使用 Pysal 工具来分析几个 shapefile（多边形类型），其中包含在 QGIS 中完成的分割结果。

我使用加载 shapefile

segshp = pysal.weights.Queen.from_shapefile(vect)

是对应的 shapefile 的名称。

问题是，当我查看邻居字典时，如果一个多边形共享少于 2 个像素，则它不被视为邻居。让我表现得更好一点：

使用 ArcGIS，我计算了多边形邻居（分析工具 --> Proximity ---> 多边形邻居），例如，多边形 1227 有 4 个邻居，如图所示。

使用 ArcGIS 计算的邻居

但是当使用 Pysal 在 Python 中加载时，多边形 1227 只有 3 个邻居，而缺少的一个是 3307，它只共享 2 个像素的边界。其他仅共享 1 或 2 个像素的多边形也会发生同样的情况。这个 2 像素的“阈值”正是我在分析一些多边形时看到的。

使用 Pysal 计算的邻居

为了考虑所有邻居，是否有改变“阈值”的方法？我在 Pysal 文档中看到过，并认为它可能与命令有关

pysal.weights.user.min_threshold_dist_from_shapefile

或类似的东西，但老实说，我不知道。

我也尝试了 pysal.weights.Rook... 选项，但它识别的邻居更少，发现更多的岛屿。

非常感谢

my_table我在 PostgreSQL 数据库中有一个带有geometry列和其他属性的点表。我有一些my_table如下的示例数据（属性my_table）。

使用 Python PySAL包，我想分析列val1中的值是否是空间自相关的（Moran I）（通过交互绘制它们）。我对交互式空间自相关的预期输出可能类似于（图像源，此处）：

我是 Python 新手。有人可以建议我如何使用 PySAL 做到这一点吗？

我有这个数据框df1，

我在这里构建了搜索树，

现在，我正在尝试搜索数据框 df1 中每个地理点 1KM 范围内的所有地理点。这是我的做法，

这里一切正常，但是我试图并行化这个任务，因为我的 df1 大小是 2M 记录，这个过程运行了 8 个多小时。谁可以帮我这个事？另一件事是，query_ball_point 返回的结果是一个列表，因此当我处理大量记录时它会抛出内存错误。任何方式来处理这个。

编辑：- 内存问题，查看 VIRT 大小。

我正在使用 Fiona 编写一个 shapefile，其中包含美国人口普查国会选区和人口普查区域的几何图形以及一些计算属性。不幸的是，我的文件被输出为空白。我可以看到 QGIS 中的记录，但看不到几何图形。

在我的代码中，我从人口普查 shapefile 中读取形状优美的多边形，计算属性并将具有相关属性的多边形提供给 fiona 文件编写器。

self.w 是在传递给对象后对输出的引用，而 self.subunit 是一个 pysal 多边形几何，我将其转换为 shapely 以便将其写入我的 shapefile。有没有人在使用 fiona 编写这个 shapefile 时遇到问题并知道我的问题可能出在哪里？