R实例:根据经纬度坐标批量返回行政区域信息

2019独角兽企业重金招聘Python工程师标准>>> hot3.png

library(RCurl)
data<-list()
site<-read.table("C:/Users/steph/Desktop/visual/site.csv",sep=",",header=TRUE)
len<-nrow(site)
head(site)
       lon      lat
1 113.1135 23.02531
2 112.5148 22.80004
3 113.0218 22.95940
4 113.1344 23.01153
5 112.8725 23.04658
6 113.0884 22.96158
for(n in 1:len)
	{
	lon<-site$lon[n]
	lat<-site$lat[n]
	out<-getURL(paste("https://geoindex.amap.com/geoIndex/area?callback=siptea_200823435167697&adGrade=4®ion=",lon,"%2C",lat,sep=""))
	res<-strsplit(out,'\"data\"')[[1]][2]
	res<-strsplit(res,'\",\"')[[1]][c(3,5,7,9)]
	res<-gsub('\":\"',"",gsub("[a-z]","",res))
	region<-data.frame(t(res))
	data[[n]]<-region
	next
	}
mydata<-do.call(rbind,data)
names(mydata)<-c("province","city","district","town")
print(mydata)
   province   city district       town
1    广东省 佛山市   禅城区   祖庙街道
2    广东省 佛山市   高明区     更合镇
3    广东省 佛山市   禅城区     南庄镇
4    广东省 佛山市   禅城区 石湾镇街道
5    广东省 佛山市   三水区     白坭镇
6    广东省 佛山市   顺德区     乐从镇
7    广东省 佛山市   高明区   荷城街道
8    广东省 佛山市   南海区   桂城街道
9    广东省 佛山市   顺德区     乐从镇
10   广东省 佛山市   南海区     大沥镇
11   广东省 佛山市   禅城区 石湾镇街道
12   广东省 佛山市   顺德区   容桂街道
13   广东省 佛山市   南海区     狮山镇
14   广东省 佛山市   南海区     大沥镇
15   广东省 佛山市   顺德区   大良街道
16   广东省 佛山市   顺德区   大良街道
17   广东省 佛山市   顺德区     乐从镇
18   广东省 佛山市   三水区     乐平镇
19   广东省 佛山市   禅城区   祖庙街道
20   广东省 广州市   南沙区   南沙街道
21   广东省 广州市   天河区   五山街道
22   广东省 广州市   南沙区     东涌镇
23   广东省 广州市   番禺区   桥南街道
24   广东省 广州市   海珠区   凤阳街道
25   广东省 广州市   天河区   石牌街道
26   广东省 广州市   天河区   猎德街道

转载于:https://my.oschina.net/u/3093769/blog/1142709

你可能感兴趣的:(r语言,python)