%% 读入全国省市经纬度坐标
clear;
% 初始化参数
coordinate_file = 'lat_lng.xlsx'; % 经纬度坐标
[num,txt] = xlsread(coordinate_file);
%% 基本信息获取
Lat=num(1:end,1);
Lng=num(1:end,2);
Max_lat=max(Lat);
Min_lat=min(Lat);
Max_lng=max(Lng);
Min_lng=min(Lng);
figure(1);
plot(Lng,Lat,'b.','markersize',6);
%% 聚类分析
[type,center]=kmeans(num,7);
rows=size(center,1);
for i=1:rows
fprintf('聚类号为%1.f的聚类中心坐标为：%f,%f\n',i,center(i,1),center(i,2));
end
h=figure(2);
hold on;
plot(Lng(type==1),Lat(type==1),'ro','markersize',6);
plot(Lng(type==2),Lat(type==2),'bx','markersize',6);
plot(Lng(type==3),Lat(type==3),'g*','markersize',6);
plot(Lng(type==4),Lat(type==4),'m.','markersize',6);
plot(Lng(type==5),Lat(type==5),'c+','markersize',6);
plot(Lng(type==6),Lat(type==6),'m.','markersize',6);
plot(Lng(type==7),Lat(type==7),'c+','markersize',6);
plot(center(:,2),center(:,1),'ro','markersize',12) ;
hold off;
————————————————
版权声明：本文为CSDN博主「LLG_llg_llg」的原创文章，遵循CC 4.0 by-sa版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/LLG_llg_llg/article/details/53192840
图片.png
