R實例：根據經緯度座標批量返回行政區域信息

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&region=",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   廣東省 廣州市   天河區   獵德街道