R语言计算31省份（除港澳台地区）省会城市的球面距离代码

R语言计算31省份（除港澳台地区）省会城市的球面距离代码

require(geosphere)
x=c(126.657717,111.660351,87.564988,125.3222,103.823305,116.395645,112.550864,117.210813,108.939,106.206479,101.767921,117.024967,91.111891,113.649644,
117.282699,114.3162,112.979353,115.893528,102.714601,106.709177,104.0648,106.55,119.330221,113.30765,108.297234,114.522082,118.778074,123.432791,121.4788,
120.219375,110.330802) 
y=c(45.773225,40.828319,43.84038,43.816,36.064226,39.929986,37.890277,39.14393,34.342,38.502621,36.640739,36.682785,29.662557,34.75661,
31.866942,30.581084,28.213478,28.689578,24.882,26.629907,30.57,29.5647,26.047125,23.120049,22.806493,38.048958,32.057236,
41.808645,31.2303,30.259244,20.022071)
muer.lonlat = cbind(x, y)
rownames(muer.lonlat)=c("黑龙江","内蒙古","新疆","吉林","甘肃","北京","山西","天津","陕西","宁夏","青海","山东","西藏","河南","安徽","湖北",
"湖南","江西","云南","贵州","四川","重庆","福建","广东","广西","河北","江苏","辽宁","上海","浙江","海南")
names<-c("北京","天津","河北","山西","内蒙古","辽宁","吉林","黑龙江","上海","江苏","浙江","安徽","福建","江西","山东","河南","湖北","湖南",
"广东","广西","海南","重庆","四川","贵州","云南","西藏","陕西","甘肃","青海","宁夏","新疆")
muer.lonlat2=muer.lonlat[names,]
muer.dists = distm(muer.lonlat2, fun=distVincentyEllipsoid);
rownames(muer.dists) =names
colnames(muer.dists) =names

x和y分别代表经度和维度