[英]How to compute the area or number of pixels for selected region in R?
我正在栅格文件上绘制一个区域,但我需要知道此框中覆盖了多少像素(通过面积表示):栅格文件为1440 pixels*720 lines``(25km*25km)
。
例:
saf <- stack(system.file("external/rlogo.grd", package="raster"))
plotRGB( saf )
e <- drawExtent()
因此,在此之后我将e
绘制为一个框,但是面积为多少像素/多少? 谢谢
尝试使用raster::crop
...
crop(saf , e )
#class : RasterBrick
#dimensions : 40, 50, 2000, 3 (nrow, ncol, ncell, nlayers)
#resolution : 1, 1 (x, y)
#extent : 23, 73, 26, 66 (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=merc
#data source : in memory
#names : red, green, blue
#min values : 0, 0, 0
#max values : 255, 255, 255
如果您只想要细胞数量...
ncell( crop(saf , e ) )
#[1] 2000
并消除NA。
x <- crop( saf , e )
ncell( ! is.na(x[]) )
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.