如何使用范围信息获取栅格的单元格编号？

Question

我有一个光栅，我正在使用光栅包。

class       : RasterLayer 
dimensions  : 103, 118, 12154  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -83.075, -82.09167, 34.95833, 35.81667  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\fb.tif 
names       : fdr_fb 
values      : 1, 128  (min, max)

我知道如何分配和所有。但是如何通过使用Lat-Long值找到cellnumber（首选）或cellvalue？

例如，我可以使用lat / long找到单元格值：

extract(ras,SpatialPoints(cbind(-82.8,35.2))) 

但我想找到对应于（Say）Long = -82.1和Lat = 35.0的单元格编号（row，col）

栅格：{{3}}

Answer 1

要从某个点获取单元格编号，您可以执行以下操作：

cellFromXY(ras, cbind(-82.8, 35.2))

如果你有一个Extent对象e，你可以这样做：

cellsFromExtent(ras, e)