如何解决通过 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 中)。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。