我在使用 terra 包时遇到了问题。当我对 weights=TRUE 的多边形数据集执行栅格提取时,有时会收到许多警告消息。从对警告的解释来看,似乎 terra 正在尝试写入临时文件来执行此操作(?)并且没有写入权限。因此,我将 terra 包的 tempdir 重新指定为我创建的具有读写权限的文件夹,但我仍然收到这些警告。
是否有其他人收到来自 terra 的类似警告并且知道如何解决问题,或者这些警告是否可能产生任何后果(即它们会影响提取值的可靠性)?
谢谢阅读。
在下面的代码中, bisondist 是一个栅格图层,US.County 是一个带有多边形的 sf 对象(10 是包含几何图形的列)
bison.mean.us <- terra::extract(bisondist,vect(US.County[,10]), fun=mean,na.rm=TRUE,weights=TRUE)
warnings()
Warning messages:
1: Attempt to create new tiff file `/spat_j1J6nrrhKYxy9MQ.tif' failed: Read-only file system (GDAL error 4)
2: Pointer 'hDS' is NULL in 'GDALGetRasterBand'.
(GDAL error 10)
3: Pointer 'hObject' is NULL in 'GDALSetDescription'.
(GDAL error 10)
4: Pointer 'hBand' is NULL in 'GDALSetRasterNoDataValue'.
(GDAL error 10)
5: Pointer 'hBand' is NULL in 'GDALFillRaster'.
(GDAL error 10)
6: Pointer 'hDS' is NULL in 'GDALSetGeoTransform'.
(GDAL error 10)
7: Pointer 'hDS' is NULL in 'GDALSetProjection'.
(GDAL error 10)
8: Pointer 'hDS' is NULL in 'GDALGetRasterBand'.
(GDAL error 10)
9: Pointer 'hBand' is NULL in 'GDALGetRasterDataType'.
(GDAL error 10)
10: Pointer 'hBand' is NULL in 'GDALGetRasterNoDataValue'.
(GDAL error 10)
11: Pointer 'hDS' is NULL in 'GDALRasterizeGeometries'.
(GDAL error 10)