微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

在R中将矩形像素转换为正方形像素

如何解决在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 举报,一经查实,本站将立刻删除。