读取经纬度数据:使用raster、sp或sf包中的函数来读取已有的经纬度数据。例如,可以使用raster()函数来读取栅格数据,readOGR()函数来读取矢量数据。 设定栅格参数:栅格化需要设置栅格的分辨率和范围。使用raster或sp包中的函数来设定栅格参数,如res()和extent()函数。 进行栅格化:使用raster或sp包中的函数来对经纬度数...
zion_b = raster::raster("C:/Users/Administrator/Google云端硬盘/caohai_s2.tif",band = 3) print(zion_r) zion_rbg = sqrt(raster::stack(zion_r, zion_g, zion_b))#拉伸 raster::plotRGB(zion_rbg,scale=100)#显示 4.统一投影坐标系(采用双线性内插) zion_elevation_utm = raster::projectRas...
r <- raster("D:/R-3.6.3/library/raster/external/test.grd") r的数据结构如下: 可以看出,raster工具包的栅格对象与sp工具包中的矢量对象的数据结构是类似的。 r ## class : RasterLayer ## dimensions : 115, 80, 9200 (nrow, ncol, ncell) ## resolution : 40, 40 (x, y) ## extent : 1784...
同时,边界匹配界限参数还可以用于裁剪栅格,只显示感兴趣区域的数据,提高数据可视化和分析的效果。 在R中,绘制栅格可以使用多个包和函数,例如raster包中的raster函数和plot函数,以及ggplot2包中的geom_raster函数等。具体使用哪个函数取决于数据的类型和需求。 以下是一些常见的边界匹配界限参数及其作用: xmin、xmax、...
plot(r)## Here is where the graticule routine startscrs.longlat<-CRS("+init=epsg:4326")prj<-CRS(projection(r))extLL<-projectExtent(r,crs=crs.longlat)lons<-pretty(c(xmin(extLL),xmax(extLL)))lats<-pretty(c(ymin(extLL),ymax(extLL)))## optionally, specify the extents of the mer...
在shape转栅格之前,首先需要建议一个新的空白的栅格,并指定控制栅格分辨率的行列,用extent制定空间范围 blank_raster<-raster(nrow=100,ncol=100,extent(bound)) 接下来给栅格赋值 values(blank_raster)<-1 plot(blank_raster) 因为给栅格的赋值都为1,因此上图显示的也只有一个颜色,你也可以给栅格赋值多个,即nrow...
R plot.raster 绘制光栅图像R语言 plot.raster 位于graphics 包(package)。 说明 此函数实现了光栅图像的 plot 方法。 用法 ## S3 method for class 'raster' plot(x, y, xlim = c(0, ncol(x)), ylim = c(0, nrow(x)), xaxs = "i", yaxs = "i", asp = 1, add = FALSE, ...) ...
extent(r) animate对Raster *对象的图层进行动画处理 # Tue Sep 08 13:47:47 2020 --animate对Raster *对象的图层进行动画处理---b<-brick(system.file("external/rlogo.grd",package="raster"))# 创建一个RasterBrick对象plot(b)animate(b,n=2)#基础# 多参数:animate(b,n=2,# 循环2次pause=1# 层...
下面介绍R的实现方式,主要是scale_fill_stepsn/scale_fill_binned的函数设置。 scale_fill_stepsn(name = "Priority ranks", # 改图例名称 labels = rev(c("2.5%", "5%", "10%", "17%","30%","50%","75%","100%")), # 改图例标签 ...
# class : RasterLayer # dimensions : 1367, 1697, 2319799 (nrow, ncol, ncell) # resolution : 1, 1 (x, y) # extent : 731453, 733150, 4712471, 4713838 (xmin, xmax, ymin, ymax) # crs : +proj=utm +zone=18 +datum=WGS84 +units=m +no_defs ...