sparkR-根据GPS坐标判断行政区域（二）

最新推荐文章于 2022-04-08 20:38:52 发布

VIP文章 tututototo

最新推荐文章于 2022-04-08 20:38:52 发布

阅读量1.6k

点赞数

分类专栏： R 文章标签： sparkR 点区域判定

本文链接：https://blog.csdn.net/tututototo/article/details/53392595

版权

# sparkR-根据GPS坐标判断行政区域（二）
上篇主要使用高德API接口调用，判断其行政区域，本次采用R里的sp包对点进行区域判定，这样做的主要目的是在先建立好模型后，直接调用模型，能否减少时间的开销。
边界坐标需要上高德API下载，其中城市POI
城市编码

1.sp简介

library(RCurl)
library(sp) 
library(stringr)
library(leaflet)
##抓取北京市边界
url_head1 <- "http://restapi.amap.com/v3/config/district?keywords="
url_tail1 <- "&subdistrict=0&extensions=all&output=xml&key=自己的key"
city <- "110000"
url1 <-  paste0(url_head1, city, url_tail1)
str1 <- getURL(url1)
rule <- gregexpr("<polyline>(.*?)</polyline>", str1)
value <- as.data.frame(unlist(regmatches(str1, rule)))
value <- gsub(".*<polyline>(.*?)</polyline>.*", '\\1', value$`unlist(regmatches(str1, rule))`)
value
#转换为data.frane
value_data <- data.frame(lonlat = unlist(strsplit(value,";")),stringsAsFactors = F)
value_data <- data.frame(str_split_fixed(value_data$lonlat, ",", 2),stringsAsFactors = F)
names(value_data) <- c("lon","lat")
value_data$lon <- as.double(value_data$lon)
value_data$lat <- as.double(value_data$lat)
#转换为Sp格式
r1 = cbind(value_data$lon,value_data$lat)
sr1=Polygons(list(Polygon(r1)),"r1")
sr=SpatialPolygons(list(sr1))
#测试数据
test = data.frame(lon=c(115.5660583464,116.0255749621,116.3999056987,116.4603790196,116.7411326333,117.3849387344,115.1680472056),
               lat=c(40.1109818799,39.8600058894,39.9031383143,39.5424745878,39.5418774826,40.6195326408,39.6489488846)