R 生成坐标向量以过滤特定区域内的数据

R 生成坐标向量以过滤特定区域内的数据,r,ggmap,geocode,data-handling,get-mapping,R,Ggmap,Geocode,Data Handling,Get Mapping,我是这里的新手,非常感谢你们能给我的任何帮助 我的项目要求我在我选择的城市周围创建一个矢量边界框,然后过滤大量数据,这样我就只有相对于该区域的数据。然而,我已经使用R studio好几年了,可以说我对该语言的记忆很少甚至完全没有 我最初使用 geocode(“英国赫里福德”) bbox现在,使用sf软件包可以非常轻松地完成大部分空间工作 下面是一个类似问题的示例代码。评论解释了它的大部分功能 困难的部分可能是理解地图投影(crs)。一些使用单位(米、英尺等),另一些使用纬度/经度。你选择哪一个取

我是这里的新手,非常感谢你们能给我的任何帮助

我的项目要求我在我选择的城市周围创建一个矢量边界框,然后过滤大量数据,这样我就只有相对于该区域的数据。然而,我已经使用R studio好几年了,可以说我对该语言的记忆很少甚至完全没有

我最初使用

geocode(“英国赫里福德”)


bbox现在,使用
sf
软件包可以非常轻松地完成大部分空间工作

下面是一个类似问题的示例代码。评论解释了它的大部分功能

困难的部分可能是理解地图投影(crs)。一些使用单位(米、英尺等),另一些使用纬度/经度。你选择哪一个取决于你在全球的哪个领域工作,以及你想完成什么。大多数网络地图使用crs 4326,但不包括易于使用的距离测量

下面的地图显示,距赫里福德约3英里外的点为红色,而里边的点为深栗色。蓝点用作赫里福德和缓冲区的中心

库(tidyverse)
图书馆(sf)
#>链接到GEOS 3.6.2、GDAL 2.2.3、项目4.9.3
图书馆(地图视图)
种子(4)
#hereford近似位置,ggmap需要api密钥
赫里福德%
st_as_sf(坐标=c('lat','lon'))%>%st_集合(4326)
#ish hereford附近数据点的模拟
随机分数%
st_as_sf(coords=c('lat','lon'))%%>%st_集(4326)%%>%st_变换(27700)
#在赫里福德周围设置约3英里(4800m)的缓冲区
h_缓冲区%st_变换(27700)%>%#将crs更改为以米为单位测量的crs
st_缓冲区(4800)
#仅限~3mi缓冲区内的点
具有6个要素和1个字段的简单要素集合中的点
#>几何图形类型:点
#>尺寸:XY
#>bbox:xmin:346243.2ymin:239070.3xmax:355169.8ymax:243011.4
#>CRS:EPSG:27700
#>点几何
#>1分(353293.1 241673.9)
#>3分(349265.8 239397)
#>4分(349039.5 239217.7)
#>6分(348846.1 243011.4)
#>7分(355169.8 239070.3)
#>10分(346243.2 239690.3)
#如地图视图所示
地图视图(赫里福德,颜色=‘蓝色’)+
地图视图(随机_点,颜色为“红色”,图例为F,列区域为“红色”)+
地图视图(h_缓冲区,图例=F)+
地图视图(其中的点,颜色为“黑色”,图例为F,列区域为“黑色”)


由(v0.3.0)于2020年4月12日创建的大多数空间工作现在可以使用
sf
软件包轻松完成

下面是一个类似问题的示例代码。评论解释了它的大部分功能

困难的部分可能是理解地图投影(crs)。一些使用单位(米、英尺等),另一些使用纬度/经度。你选择哪一个取决于你在全球的哪个领域工作,以及你想完成什么。大多数网络地图使用crs 4326,但不包括易于使用的距离测量

下面的地图显示,距赫里福德约3英里外的点为红色,而里边的点为深栗色。蓝点用作赫里福德和缓冲区的中心

库(tidyverse)
图书馆(sf)
#>链接到GEOS 3.6.2、GDAL 2.2.3、项目4.9.3
图书馆(地图视图)
种子(4)
#hereford近似位置,ggmap需要api密钥
赫里福德%
st_as_sf(坐标=c('lat','lon'))%>%st_集合(4326)
#ish hereford附近数据点的模拟
随机分数%
st_as_sf(coords=c('lat','lon'))%%>%st_集(4326)%%>%st_变换(27700)
#在赫里福德周围设置约3英里(4800m)的缓冲区
h_缓冲区%st_变换(27700)%>%#将crs更改为以米为单位测量的crs
st_缓冲区(4800)
#仅限~3mi缓冲区内的点
具有6个要素和1个字段的简单要素集合中的点
#>几何图形类型:点
#>尺寸:XY
#>bbox:xmin:346243.2ymin:239070.3xmax:355169.8ymax:243011.4
#>CRS:EPSG:27700
#>点几何
#>1分(353293.1 241673.9)
#>3分(349265.8 239397)
#>4分(349039.5 239217.7)
#>6分(348846.1 243011.4)
#>7分(355169.8 239070.3)
#>10分(346243.2 239690.3)
#如地图视图所示
地图视图(赫里福德,颜色=‘蓝色’)+
地图视图(随机_点,颜色为“红色”,图例为F,列区域为“红色”)+
地图视图(h_缓冲区,图例=F)+
地图视图(其中的点,颜色为“黑色”,图例为F,列区域为“黑色”)

由(v0.3.0)于2020年4月12日创建