R语言,raster包area函数,细胞大小 # R 版本:R x64 4.0.2# cgh163email@163.com#个人笔记不负责任# —— 拎了个梨🍐rm(list=ls());gc().rs.restartR()require(raster)r<-raster(nrow=18,ncol=36)a<-area(r)plot(a)# Tue Sep 08 15:38:04 2020 -require(rgdal)require(rgeos)p<-shapefile(...
如果要转换栅格数据的坐标系统,使用的是projectRaster()函数: elev.r.gcs <- projectRaster(elev.r, crs="+proj=longlat +datum=WGS84") crs(elev.r.gcs) ## CRS arguments: +proj=longlat +datum=WGS84 +no_defs 如果要使用EPSG代码,需要采用+init=EPSG:...语句: elev.r.gcs <- projectRaster(elev...
par(mfrow=c(1,2))# value= ID_2shape_r=rasterize(shape,r,"ID_2")plot(shape_r)plot(shape,add=T)title(main="value=ID_2")shape_r# value= AREAshape_r=rasterize(shape,r,"AREA")plot(shape_r)plot(shape,add=T)title(main="value=AREA")shape_r image.png NA处理 有时候生成的raster里面...
# value=AREAshape_r=rasterize(shape,r,"AREA")plot(shape_r)plot(shape,add=T)title(main="value=AREA")shape_r image.png NA处理 有时候生成的raster里面有NA数据,那么如何替换掉呢,(reclassify)[http://search.r-project.org/library/raster/html/reclassify.html]可以实现该过程。主要参数cbind(0,a,b...
1library(raster) 2library(sp) 3library(dplyr) 4library(magrittr) 5 6# 定义一个函数,将rasterLayer栅格数据转化为data.frame 7# 将rasterLayer栅格数据转化为dataframe数据。 8rasterL_to_DF <- function(climate_mask) { 9climate_mask_df <- as.data.frame(cbind(coordinates(climate_mask),# 合并坐...
("E:/01RawData/test2/",ndvilist)shp=vect("E:/02ResearchArea/test1.shp")for(iin1:length(ndvidir)){ndvi=rast(ndvidir[i])ndvi_shp=trim(mask(ndvi,shp))writeRaster(ndvi_shp,filename=paste0("E:/03预处理/test3/",ndvilist[i]))print(paste("finish...
(Income)) %>% .$medinc#计算面积ME.inc$Area<-gArea(ME.inc,byid=TRUE)/1000000#intersect()截取相交部分,rasterclp1<-intersect(s1,s2)#按照输入数据顺序,来决定输出数据类型,clp2为线空间数据clp2<-intersect(l1,s1)#clp3为面空间数据clp3<-intersect(s1,l1)#union()相交融合,raster包un1<-union(s1...
area <- st_area(intersection) 最后,可以打印出阴影多边形的面积: 代码语言:txt 复制 print(area) 需要注意的是,以上步骤中的x1, y1, x2, y2, ...是阴影多边形的坐标点,可以根据实际情况进行替换。另外,还可以根据具体需求对地图数据进行处理,例如裁剪、缩放等操作。
circos.lines(x, y, col, area, baseline, border) 例如 n <- 90 df <- data.frame( sectors = rep(letters[1:9], each = 10), x = rep(1:10, 9), y = runif(n) ) circos.par("track.height" = 0.3) circos.initialize(df$sectors, x = df$x) ...
# AREA_OR_POINT=Area 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. 27. 28. 加载栅格 #load raster into R DSM_HARV <- raster(paste0(wd,"HARV/DSM/HARV_dsmCrop.tif")) ...