如何解决检查栅格是否翻转
我有一个具有以下属性的栅格
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)
我可以确定翻转的(下部面板)是正确的栅格,因为它跟随 我感兴趣的土地(马来西亚,印度尼西亚)
我的问题是,有什么方法可以检查栅格是否确实被翻转?
解决方法
GDALinfo输出ysign -1
建议不要翻转数据。那就是文件包含的所有信息。
但是,由于您提供的原因,数据显然被翻转了,文件在这方面是错误的。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。