% 已知圆心和半径画圆:参考:https://blog.csdn.net/ZLK961543260/article/details/70216089
% 对比三种画圆方法,运算时长对比如下
% viscircle:0.22;比较快
% rectangle:0.21;比较快
% function:0.39:慢
[Type Sheet Format]=xlsfinfo('现状OD数据及其他数据.xls');
OD=xlsread('现状OD数据及其他数据.xls',Sheet{1});
center_area=xlsread('现状OD数据及其他数据.xls',Sheet{2});
x_pos=center_area(:,4); %数据第四列为x坐标(米)
y_pos=center_area(:,5); %数据第五列为y坐标(米)
area=center_area(:,3); %数据第3列为面积(平方米):注意与x,y的坐标对应
r_all=sqrt(area/pi); %计算节点覆盖半径,为后文画圆做准备
nanbool=isnan(r_all); %数据1到4行为四个物流园区,通过r为nan,提取出四个物流园区
first_node_num=length(find(nanbool==1)); %统计物流园区个数
r_index=find(nanbool==0,1); %节点开始的行数index
r=r_all(r_index:end); %r存储节点的覆盖半径
jam_coef=center_area(r_index:end,6); %数据第六列