用 PySal 计算 Moran's I？答案

【问题标题】：Calculating Moran's I with PySal?用 PySal 计算 Moran's I？
【发布时间】：2017-11-16 15:29:48
【问题描述】：

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

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

import os
import pysal
import numpy as np

os.chdir = ("C:\\Users\\hsv0336\\Desktop")
wd = os.getcwd()

w = pysal.open(wd + "\\Rook.gal")
f = pysal.open(wd + "\\NewYorkData.txt")
y = np.array(f.by_col["T0_MINC"])

mi = pysal.Moran(y, w, two_tailed=False)

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

【问题讨论】：

标签： python spatial pysal

【解决方案1】：

pysal.open 尝试根据扩展名和检查文件来确定文件类型。 pysal.open 返回的类型并不总是清楚的，应该使用内置的类型函数进行检查。在您的情况下，它返回一个普通的 python 文件对象，这意味着 pysal 无法为您解析它。 “by_col”方法仅在 PySAL 将您的文件识别为 DataTable 时可用。

支持的 DataTable 类型包括带有特殊“GeoDa”标头的 .csv、.dbf 或 .txt 文件。 GeoDa 表头在第一行包含数据行数和列数，在第二行包含列名。

例如示例 STL_HOM.txt 数据集的前几行。其中有 78 行数据和 4 列。

78,4
"FIPSNO","HR8488","HR8893","HC8488"
17107,1.290722,1.624458,2

您可以将文件重新格式化为一种受支持的 DataTable 类型，或者编写自己的代码来解析您的文本文件。

【讨论】：