如果需要将坐标转换为其他坐标系统,可以使用st_transform函数。这一步是可选的,根据实际需求来决定是否进行坐标转换。 sf_object<-st_transform(sf_object,crs=3857) 1. crs参数指定了目标坐标系统,这里使用了Web墨卡托投影坐标系统,其对应的EPSG代码为3857。 至此,我们完成了使用R语言的st_as_sf函数实现的流程。
可以看出,st_read()函数相比于另外两个函数调用格式简单地多,运行速度也更快。 st_read()在sf中还有一个别名函数read_sf(),其使用方式与st_read()一样,区别在于它读取数据后以tibble格式储存,而tibble是tidyverse工具包推出的一种新式的数据框格式,其相比于普通的data.frame在格式上更加规范。 spdata41 <- re...
在sf::st_as_sf()这种形式的函数调用方式中,::左边为包名称,::右边为该包中的函数名称,例如这里我们指定“使用sf包中的st_as_sf()函数”。这种调用方式支持包的动态加载,即以上代码即使在为调用library(sf)命令的环境中也可以执行。当然,用户必须已经安装了sf包到当前的R环境中。更重要的是,这种调用方式可以...
我们可以使用sf的st_as_sf()函数来做到这一点。 (注意:对空间数据进行操作的sf包函数以st_开头, st_代表“ spatial”和“ temporal”。) [ 通过InfoWorld的机器学习和分析报告时事通讯来了解机器学习,人工智能和大数据分析的最新进展 st_as_sf()接受几个参数。 在下面的代码中,第一个参数是要转换的对象-我的...
有时候我们没有.shp文件,只有点的信息数据。这时候可以利用st_as_sf函数进行创建空间sf类。 df<-data.frame(lon=c(-68.783, -69.6458, -69.7653),lat = c(44.8109,44.5521,44.3235),Name= c("Bangor","Waterville","Augusta"))#coords=设置坐标,crs=设置坐标系,4326代表WGS84坐标系point<-st_as_sf(df,...
grid(x = seq(-1,1,.1),y = seq(-1,1,.1)) %>% st_as_sf(coords = c("x","y")) ### 分组创建 multipoints 对象 data %>% st_as_sf(coords = c("x", "y")) %>% group_by(label) %>% summarise() %>% st_cast('MULTIPOINT') ### 根据标签分组识别点集的轮廓 hll_group...
scatter_df_tro<-st_as_sf(scatter_df,coords=c("lon","lat"),crs=4326) 结果如下: 注意红色框中的内容 接下来将数据在地图中显示,具体代码如下: 代码语言:javascript 复制 ggplot()+geom_sf(data=china,fill="NA",size=.5,color="black")+geom_sf(data=nine_line)+geom_sf(data=scatter_df_tro,...
sites.sf <- st_as_sf(sites, coords = c("lon", "lat"), crs=4326)#25831 sites.sf.robin <- st_transform(sites.sf, proj) 导入底图 # land = read_sf('10m_physical/ne_10m_land.shp') %>% # # subset(CONTINENT != 'Antarctica') %>% ...
simple featrue对象定义了两个序列化标准:常见文本表现形式(WKT)和二进制文本表现形式(WKB)。常见文本表现形式是日常打印时默认的输出格式,sfc列可以利用st_as_sfc函数从WKT格式的字符串向量中直接读取。 代码语言:javascript 复制 >library(sf)Linking toGEOS3.5.1,GDAL2.1.2,proj.44.9.3>(pt<-st_point(c(0,...
R语言的空间包(sp)可以帮助我们处理空间数据,如st_as_sf()函数可以将点数据转换为空间数据框,方便后续分析。基于流行病学研究设计的R语言实践应用涵盖了描述性统计分析、相关性分析、生存分析、多元线性回归分析等多个方面。掌握这些技能将有助于我们更有效地开展流行病学研究,为公共卫生政策制定提供有力支持。1....