如果需要将坐标转换为其他坐标系统,可以使用st_transform函数。这一步是可选的,根据实际需求来决定是否进行坐标转换。 sf_object<-st_transform(sf_object,crs=3857) 1. crs参数指定了目标坐标系统,这里使用了Web墨卡托投影坐标系统,其对应的EPSG代码为3857。 至此,我们完成了使用R语言的st_as_sf函数实现的流程。
shpdata <- st_read("路径/Shandong16.shp") # 确保tmap处于静态地图模式 tmap_mode("plot") # 绘图 tm_shape(shpdata) + tm_polygons(border.col = "white") + # 添加聊城点 tm_shape(st_as_sf(lc)) + tm_dots(col = "red", size = 0.1) + tm_text("name", just = "left", xmod =...
可以看出,st_read()函数相比于另外两个函数调用格式简单地多,运行速度也更快。 st_read()在sf中还有一个别名函数read_sf(),其使用方式与st_read()一样,区别在于它读取数据后以tibble格式储存,而tibble是tidyverse工具包推出的一种新式的数据框格式,其相比于普通的data.frame在格式上更加规范。 spdata41 <- re...
我们可以使用sf的st_as_sf()函数来做到这一点。 (注意:对空间数据进行操作的sf包函数以st_开头, st_代表“ spatial”和“ temporal”。) [ 通过InfoWorld的机器学习和分析报告时事通讯来了解机器学习,人工智能和大数据分析的最新进展 st_as_sf()接受几个参数。 在下面的代码中,第一个参数是要转换的对象-我的...
解答: 可以使用st_as_sf()函数将地理数据转换为网格数据。 gridded_data <st_as_sf(data, crs = "EPSG:4326") # 将地理数据转换为EPSG:4326坐标系的网格数据 问题2: 如何使用不同的颜色对地理要素进行分类? 解答: 可以在geom_sf()函数中使用aes()函数来指定分类变量和颜色映射。
有时候我们没有.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,...
使用terra包。如果要将光栅转换为sf对象，方法是首先使用terra转换为多边形，然后使用st_as_sf：
3.将csv转为sf对象 b=st_as_sf(a,coords=c("longitude","latitude"),crs=4326) 4.将sf对象导出为shp st_write(b,'F:/test/test.shp',driver = "ESRI Shapefile") 5.batch files<-list.files(pattern=".csv$")dir<-paste0("F:/test/",files)for(iin1:20){sf<-st_as_sf(read.csv(dir[i...
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,...
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...