我正在尝试提取两个重叠栅格集的数据(一个，一堆35个栅格，全部来自同一个源，第二个是高程栅格)，以获取每个栅格的值(值的平均值)的data.frame所有栅格的像素 .

栅格堆栈的描述如下：

> stack_pacifico

class : RasterStack

dimensions : 997, 709, 706873, 35 (nrow, ncol, ncell, nlayers)

resolution : 0.008333333, 0.008333333 (x, y)

extent : -81.62083, -75.7125, 0.3458336, 8.654167 (xmin, xmax, ymin, ymax)

coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0

names : F101992.v//ts.avg_vis, F101993.v//ts.avg_vis, F101994.v//ts.avg_vis, F121994.v//ts.avg_vis, F121995.v//ts.avg_vis, F121996.v//ts.avg_vis, F121997.v//ts.avg_vis, F121998.v//ts.avg_vis, F121999.v//ts.avg_vis, F141997.v//ts.avg_vis, F141998.v//ts.avg_vis, F141999.v//ts.avg_vis, F142000.v//ts.avg_vis, F142001.v//ts.avg_vis, F142002.v//ts.avg_vis, ...

min values : 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...

max values : 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, ...

对于高程栅格：

> elevation_pacifico

class : RasterLayer

dimensions : 997, 709, 706873 (nrow, ncol, ncell)

resolution : 0.008333333, 0.008333333 (x, y)

extent : -81.62083, -75.7125, 0.3458336, 8.654167 (xmin, xmax, ymin, ymax)

coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

data source : in memory

names : COL_msk_alt

values : -16, 5164 (min, max)

这是我第一次使用栅格数据，我希望通过1平方公里(或更小)的网格提取数据 . 我知道两个光栅的分辨率都可以强制适应该区域要求，两个尺寸都相等，因此每个光栅的像素数是相同的 .

我的问题是，我是否只能合并所有栅格(堆栈和高程栅格中的栅格)并提取数据，并确信所有像素重叠(或位于同一位置)？或者我是否必须创建一个主SpatialGrid或SpatialPixel对象，然后将栅格数据提取到这些对象？

提前致谢，

点击此链接即可下载光栅堆栈中的数据(如果要下载所有堆栈，可以使用_2366629中的脚本)：

海拔：

#Download country map and filter by pacific states

colombia_departments

pacific_littoral

pacific_littoral_map

#Download elevation data and filter it for pacific states

elevation

elevation_pacifico

elevation_pacifico