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

问题描述

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

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,values(r))
#[1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"
r=10*r
 
t <- writeraster(r,overwrite=TRUE,wopt=list(datatype="INT4S"))
t=t/10
 
sprintf("%.8f",values(t))
#[1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

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

解决方法

您可以使用双精度浮点数 (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 中）。