Question

在写入外部文件时会修改栅格值，这是由于数字内部表示的经典问题，我想：

library(terra)
r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
sprintf("%.8f",values(r))
# [1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"
  
t <- writeRaster(r, "test.tif", overwrite=TRUE)

sprintf("%.8f",values(t))
#[1] "5.30000019" "7.09999990" "3.00000000" "1.20000005"

这可能是一个很大的问题。什么是“最佳”适当的解决方法？我测试了整数格式的舍入和写入，并且有效：

r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
sprintf("%.8f",values(r))
#[1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"
r=10*r
 
t <- writeRaster(r, "test.tif", overwrite=TRUE,wopt=list(datatype="INT4S"))
t=t/10
 
sprintf("%.8f",values(t))
#[1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

有更好的解决方法吗？特别是，terra 中是否有任何“精度”属性可以与每个栅格相关联，以便为用户更简单、更有效地管理这些问题？

Answer 1

您可以使用双精度浮点数 (FLT8S)

library(terra)
r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
t <- writeRaster(r, "test.tif", overwrite=TRUE, wopt=list(datatype="FLT8S"))
sprintf("%.8f",values(t))
#[1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

<块引用>

terra 中是否有任何“精度”属性，可以关联到每个栅格

不，没有。如果有，那需要传播规则吗？在写作时设置它可能更容易吗？一个警告是，一些写入临时文件可能会被忽视。要将默认值设置为 FLT8S，您可以使用 terraOptions(datatype="FLT8S")（在开发版本 1.0-2 中）。

通过 terra 管理数字精度：写入文件时修改十进制栅格值

1 个答案: