我想制作一段我更高效的代码。
现在我通过open acceshub以jp2格式下载Sentinel-2
数据。我下载的jp2文件由于某种原因而有错误的程度。现在我通过以下方式纠正这个问题(其中file
是jp2的文件名):
r = raster(file)
extent(r) = new_extent
writeRaster(r, file)
然而,这种方法会写入整个栅格(需要很长时间),而我只更改了一个小细节。
使用 gdal 或栅格包有更好的方法来更有效地执行此操作吗?
如果我打印光栅,我看到:
class : RasterLayer
dimensions : 1830, 1830, 3348900 (nrow, ncol, ncell)
resolution : 60, 60 (x, y)
extent : 499980, 609780, 6690240, 6800040 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=55 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : /home/daniel/R/farmhack/x.jp2
names : x
values : 0, 65535 (min, max)
我不知道这个程度意味着什么。