R中提取栅格区域的有效方法
我尝试使用30米(1弧秒)分辨率NED(国家海拔数据集)按州为5个州的区域提取5000英尺以下的区域。我有一个国家的shapefile,以及一个NED的大(约4GB)栅格。我正在使用velox软件包,该软件包可以大大加快速度,但是我觉得必须有一种更快的方法来完成(这仍然需要很长时间)。
使用velox软件包,我将栅格变成了velox对象,然后使用velox的提取功能提取像元。由于我知道每个像元的分辨率,因此我对每种状态下5000英尺以下的像元数进行计数,然后乘以每个像元的面积。
DEM_velox <- v...