.打开Arcmap2.新建Shapefile文件3.编辑Shapefile属性4.开始编辑5.创建要素并绘制路网6.打断相交线7.保存编辑8.打开图层属性表9.添加字段10.完成字段添加11.计算字段id12.计算点线字段13.选中length字段14.计算length字段15.导出路网数据16.查看最终的路网数据

路网数据转化为图matlab代码

通过构造邻接矩阵在构造的路网中找到一个节点通向另一个节点的最短路径包括DJFloydA*ACOGA算法。

路网与一般的无向图不同点是在一般的无向图中从起点开始不论下一次选择哪个节点最终都可以到达目标节点而路网中存在死路即算法需要回溯才能找到正确的通路

构造路网

右键ArcCatalog中文件夹链接下的目录依次选择“新建”-“Shapefile”

3.编辑Shapefile属性

在弹出的创建新Shapefile对话框中填写文件名称选择要素类型为“折线”点击“编辑”弹出空间参考属性对话框选择XY坐标系:地理坐标系-world-WGS

1984最后点击“确定”坐标将包含M值”和“坐标将包含Z值”千万不要勾选否则后面路网数据处理过程中会报错最后点击“确定”

4.开始编辑

左侧图层中出现了我们刚新增的Shapefile线类型图层我们依次点击“编辑器”-“开始编辑”

5.创建要素并绘制路网

绘制编辑完道路数据后依次点击工具栏中的“编辑器”-“更多编辑工具”-“高级编辑”选中全部要素然后点击高级编辑工具栏中的“打断相交线”

7.保存编辑

注若打断失败可以重启arcmap重新新建文件以上步骤重新来一次

8.打开图层属性表

计算字段id的内容右键Id字段选择“字段计算器”在字段计算器中给Id进行重新计算赋值双击FID字段然后点击再输入数字1即可

12.计算点线字段

计算x1,y1,x2,y2字段选中列x1并右击选择‘计算几何’选择‘线起点的x坐标’点击确定其余的y1x2y2同理

13.选中length字段

计算length字段即计算路径长度选中length列-选择字段计算器

14.计算length字段

在字段计算器中选择解析程序为python输入以下内容计算路径长度点击确定

math.sqrt((

将保存的excel数据加载到matlab中编写并运行以下的matlab代码得到路网与图

无向带权图

1将原始的road数据即保存的excel数据转化为图和邻接矩阵的matlab代码

函数功能

%(1)读取arcmap数据进行点合并绘制路网图保存处理后的数据

%(3)根据邻接矩阵构造图

%uniqueValues:不重复的节点坐标及其对应的唯一id

读取数据

start_pointdata(:,2:3);%记录起点坐标x1

length_listdata(:,6);%记录路径长度length%%

distanceThreshold3;%距离阈值两点距离小于此阈值则合并两个点

numNodes

1:5%通过迭代不断减少距离阈值distanceThresholddistanceThreshold/2;Preprocessed_data[];%记录要替换的数据for

1:numNodesfor

distanceThreshold%两点距离过近Preprocessed_data[Preprocessed_data;start_point(i,

:),start_point(j,

distanceThreshold%两点距离过近Preprocessed_data[Preprocessed_data;end_point(i,

:),end_point(j,

distanceThreshold%两点距离过近Preprocessed_data[Preprocessed_data;start_point(i,

:),end_point(j,

:)];endendendend%for%根据记录处理原始数据for

1:size(data,1)for

data(i,2)Preprocessed_data(j,1)

data(i,3)Preprocessed_data(j,2)data(i,2)Preprocessed_data(j,3);data(i,3)Preprocessed_data(j,4);end

data(i,4)Preprocessed_data(j,1)

data(i,5)Preprocessed_data(j,2)data(i,4)Preprocessed_data(j,3);data(i,5)Preprocessed_data(j,4);end

end

data(:,1);%记录idstart_pointdata(:,2:3);%记录起点坐标end_pointdata(:,4:5);%记录终点坐标length_listdata(:,6);%记录路径长度

end

2)];distancesqrt((start_point(i,

1)-start_point(i,

2),num2str(i),color,r,FontSize,

12);%图上打上名称plot(x,

1:node_lengthuniqueValues(i,3)i;

end

data(i,3)uniqueValues(j,2)new_data(i,1)uniqueValues(j,3);endif

data(i,5)uniqueValues(j,2)new_data(i,2)uniqueValues(j,3);endend

end

new_data(:,3)data(:,6);%载入路径权值即路径长度

save(data.mat,

G.Edges.Weight,EdgeLabelColor,r);

layout(h,