检查光栅是否翻转

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)

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

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

Answer 1

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

但是数据明显是颠倒了，你给的原因，这方面文件是错误的。

检查光栅是否翻转

check if a raster is flipped

r

raster

rgdal