0

我在使用 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)
4

1 回答 1

1

这不是编码问题,而是错误报告,应在此处报告。这是固定的 terra 1.4-13

于 2021-10-13T23:19:34.363 回答