在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