根据邻接矩阵计算各点之间的最短路径矩阵(ODM也可以)

最新推荐文章于 2024-06-21 13:32:35 发布

wallezhou

最新推荐文章于 2024-06-21 13:32:35 发布

阅读量6.6k

点赞数 3

文章标签：邻接矩阵最短路径 SP ODM Floyd

本文链接：https://blog.csdn.net/WillChou123/article/details/79056753

版权

本文介绍如何通过Python3和Pandas处理包含邻接矩阵的CSV文件，利用Floyd算法计算任意两点间的最短路径，并输出为新的CSV文件。适用于空间权重矩阵分析和生成ODM。

摘要由CSDN通过智能技术生成

我之前写过一篇将arcgis的swm文件处理成为保存矩阵的文本文件格式的博客，得到的是csv文件。该文件保存的空间权重矩阵。csv文件方便进一步的空间分析。接下分享一下利用存有邻接矩阵的csv文件得到最短路径矩阵(存有任意两个节点之间的最短路径)的csv文件，想获得一般的ODM也同样适用。方法用的是Floyd算法，基于python3，使用pandas和copy模块。

输入文是一个csv文件，存有空间邻接矩阵，1代表相邻，0代表不相邻。比如4和2相邻，4与3不相邻。

import pandas as pd
import copy

A=pd.read_csv(r'E:\atest\NYC\drug\linjie.csv',index_col=0)
#A是邻接矩阵，相邻对象元素为1，不相邻为0

D=copy.deepcopy(A)
#用于储存节点对的最短路径，相邻的为实际权值（本例为1），不相邻设置为很大的数（远大于所有边的权值，本例设置为999）
ilter=[i for i in range(len(A))]
#o代表起始节点ID,d是终点ID,mid是内插节点
for o in ilter:
    for d in ilter:
        if d==o:
            continue
        if D.iloc[o,d]==0:
            D.iloc[o,d]=999
print("得到矩阵D")
#D初始化完毕


#使用Floyd算法计算SP

for mid in il