Arcgis空间自相关出现脚本报错是怎么回事呢

Question

空间自相关出现脚本报错是怎么回事呢？
如果在File "d:arcgis\desktop10.8\ArcToolbox\Scripts\Moransl.py"路径里修改，要怎么修改呢？

Answer 1

引自免费微信小程序：皆我百晓生

看起来你在使用 ArcGIS 进行空间自相关的Moran's I 测试时遇到了一个 MemoryError。这通常是因为处理的数据集太大，超过了系统可用的内存。为了解决这个问题，你可以尝试以下几个步骤：

1. 增加内存：确保你的计算机有足够的RAM来处理数据。如果你的电脑内存不足，考虑升级硬件。

2. 优化工作区设置：在ArcMap中，可以通过管理器 > 工作区 > 属性 > 管理选项卡来调整工作区设置。减少缓存大小可能会有所帮助。

3. 分块分析：将大图层分成较小的部分，分别进行空间自相关分析，然后对结果进行综合。这可以通过使用ArcGIS的「切片工具」（Clipper）或者「分割工具」（Split）来实现。

4. 简化输入数据：如果可能，删除不必要的字段或降低几何精度，以减小数据集的大小。

5. 使用更快的计算方法：在Moran's I脚本中，可以尝试改变距离矩阵的计算方法。在你的例子中，使用的可能是Inverse Distance 或 Euclidean Distance。你可以尝试使用其他方法，比如Nearest Neighbor，它可能会更高效。

由于你提供的代码片段已经包含了错误发生的具体位置，我们可以尝试修改这个部分来提高效率。在文件 MoransI.py 中找到以下代码段：

def processRow(self, row, nhVal, weights):
    wij = self.wij
    for i in range(len(row)):
        for j in range(i + 1, len(row)):
            wij[ij] = wij[ij] + weights[i][j]

这里的问题在于每次迭代都在更新wij数组，这可能导致频繁的内存分配。你可以通过创建一个新的 wij 变量，只在循环结束后再赋值给原 wij 来改进这段代码：

def processRow(self, row, nhVal, weights):
    wij_temp = np.zeros((len(row), len(row)), dtype=np.float32)
    for i in range(len(row)):
        for j in range(i + 1, len(row)):
            wij_temp[i, j] = wij_temp[i, j] + weights[i][j]
    wij = wij_temp.copy()

这里我们引入了NumPy库来创建一个新的 wij_temp 数组，并在循环结束后将其复制到 wij。这样可以减少内存分配次数，提高性能。

请注意，这些步骤并不能保证一定能解决内存问题，因为数据量过大仍然是主要问题。但是它们可以帮助优化算法，减少内存占用。如果问题仍然存在，你可能需要寻求更高级的技术支持，例如使用具有更大内存的服务器，或者使用分布式计算框架（如GDAL/OGR、Python的Dask库等）进行并行处理。