Sentinel-2R gdal光栅

我想让我的代码更高效。 现在我通过open acceshub下载jp2格式的

Sentinel-2

数据。出于某种原因，我下载的jp2文件有一个错误的范围。现在，我用以下方式更正此问题（其中，

file

是jp2的文件名）：

然而，这种方法会写入整个光栅（这需要时间），而我只更改了一小部分细节

有没有一种简洁的方法可以使用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)

---

我不知道这个范围是什么意思。

什么是错误的范围？您能详细说明并可能共享

r

的控制台输出吗？光栅的控制台输出是：光栅层尺寸：183018303348900（nrow，ncol，ncell）分辨率：60，60（x，y）范围：49998060978066902406800040（xmin，xmax，ymin，ymax）坐标。参考：+proj=utm+zone=55+datum=WGS84+units=m+no_defs+ellps=WGS84+towgs84=0,0,0数据源：/home/daniel/R/farmhack/x.jp2名称：x值：0，65535（最小值，最大值），因此范围不在wgs坐标中，很难裁剪。我甚至不知道这个范围是什么意思。

Sentinel-2

瓷砖被投影到UTM中（您的是UTM 55N），它在投影范围内以米为单位显示范围。所以范围没有错，只是不在地理坐标中。谢谢你，这很有帮助。我将看看如何将我的wgs croping窗口转换为UTM than。

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)