Sentinel-2 R gdal光栅

时间:2018-04-17 20:54:59

标签: r raster gdal

我想制作一段我更高效的代码。 现在我通过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)

我不知道这个程度意味着什么。

0 个答案:

没有答案