r - 检查栅格是否翻转

Question

我有一个具有以下属性的栅格

  gdalinfo <- try (rgdal::GDALinfo(rasterFile))
  
  # Warning message:
  #  In rgdal::GDALinfo(rasterFile) : statistics not supported by this driver
  
  gdalinfo
  rows        36 
  columns     64 
  bands       1 
  lower left origin.x        91.17 
  lower left origin.y        -14.73 
  res.x       1 
  res.y       1 
  ysign       -1 
  oblique.x   0 
  oblique.y   0 
  driver      GTiff 
  projection  +proj=longlat +datum=WGS84 +no_defs 
  file        C:/earthengine/historical/raw/tmax_1987_ACCESS1-0__0.tif 
  apparent band summary:
    GDType hasNoDataValue NoDataValue blockSize1 blockSize2
  1 Float32          FALSE           0         32         64
  apparent band statistics:
    Bmin       Bmax Bmean Bsd
  1 -4294967295 4294967295    NA  NA
  Metadata:
    AREA_OR_POINT=Area 
  
  

我怀疑这个光栅是反转的，所以我翻转它

  my_ras_cel <- raster(rasterFile)
  
  # reverse it
  my_ras_cel_rev <- flip(my_ras_cel, direction='y')

将两个栅格绘制在一起

  ras_stack <- stack(my_ras_cel, my_ras_cel_rev)
  
  spplot(ras_stack)      

我可以看出翻转（下面板）是正确的栅格，因为它遵循我感兴趣的土地（马来西亚，印度尼西亚）

我的问题是有什么方法可以检查我的栅格是否确实被翻转了？

标签： rrasterrgdal

Answer 1

GDALinfo 输出 ysign -1 表明数据未翻转。这就是文件中有关此的所有信息。

但是由于您给出的原因，数据显然是翻转的，并且文件在这方面是错误的。