问题描述
我正在尝试使用'passthrough'包输出一个栅格，其中包含来自两个多波段栅格的组合波段，并在R中带有几个简单的语句。例如：
raster问题在于输入栅格没有波段统计元数据，但是在输出中，均值和StdDev均写为“ 1.＃SNAN”。例如，这是gdalinfo的一些输出：
ras <- stack("five_bands.tif","seven_bands.tif")
writeraster(ras,"output.tif")
我意识到我可以使用gdal_edit重新计算或删除波段统计元数据，但是我希望能够直接从R中直接输出具有准确统计元数据的栅格。Band 6 Block=1696x1 Type=Float32,ColorInterp=Undefined
Min=-110.000 Max=20000.000
Minimum=-110.000,Maximum=20000.000,Mean=1.000,StdDev=1.000
NoData Value=-3.39999999999999996e+38
Metadata:
STATISTICS_MAXIMUM=20000
STATISTICS_MEAN=1.#SNAN
STATISTICS_MINIMUM=-110
STATISTICS_STDDEV=1.#SNAN
函数（来自{{1 }}）似乎没有任何帮助。当writeraster（）调用它来编写GeoTiff时，似乎由rgdal编写了“ STATISTICS _...”，是否可以通过一种方法来获取以元数据形式编写的精确的Mean和StdDev，或者没有编写任何统计数据？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）