1

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

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 中是否有任何“精度”属性，可以与每个栅格相关联，从而允许用户更简单、更有效地管理这些问题？

4

1 回答 1

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 中）。

于 2020-12-19T22:28:58.473 回答