如何解决在R中将矩形像素转换为正方形像素
我正在使用R中的一个asc文件,并尝试在其他人编写的脚本上通过某些GIS处理来运行它。我正在处理的当前代码是
asc <- (raster(inputRas))
asc.extent <- extent(asc)
head <- read.ascii.grid.header(inputRas)
writeraster(asc,"outputRas",format="SAGA",overwrite=TRUE)
,直到最后一行。这里失败了,因为我的栅格是在直角坐标中的,而SAGA不喜欢。我花了很多时间在Google搜索和尝试不同的功能上,但是没有任何事情能够解决其他问题。我考虑使用projectRaster()一段时间是因为可以在其中更改分辨率,但实际上我不想更改投影。任何帮助将不胜感激,谢谢。
解决方法
您可以使用let arr = [1.0,-42.0,0.0,-5.0,42.0,7.0];
let max = arr.iter().copied().fold(f64::NAN,f64::max) // 42.0
代替raster::projectRaster
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。