我一直在尝试计算光栅文件中的波段数。虽然下面的代码大部分工作,但当我在投影的栅格上尝试代码时,它似乎抛出了一个异常。<init>(GeoTiffReader.java:266) at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:76)
我有一个很大的RasterStack (114个geotiff图像),我已经成功地对其进行了地理处理(屏蔽等)。在R中,但是我很难让它对每个栅格应用一个简单的条件语句(所有栅格层都是相同的范围,分辨率,并且是共同注册的)。我希望将所有小于每个栅格最大值的95%的像素值设置为NA。下面是我的代码: #Get max value from each raster layer
#Set all values < 95% of m