如何修复从hdf转换为R中合适的栅格的“弯曲”栅格tiff?

时间:2019-07-12 16:31:33

标签: r gdal

我正在使用MOD35之类的MODIS Atmosphere产品,我读到最好将hdf转换为geotiff,并能够在任何软件GIS中使用它。我做了以下事情:

library(raster)
library(rgdal)
library(gdalUtils)

setwd("C:\\Users\\Usuario\\Documents")
filename = "MOD35_L2.A2017003.1520.061.2017311220101.hdf"
#print(filename)

sds = get_subdatasets(filename)
#print(sds)
gdal_translate(sds[7], dst_dataset = "mascara_nubes.tif", b = 1) #sds[7] is cloud_mask sds

gdalwarp("mascara_nubes.tif", "mascara_nubes_projec.tif", 
         s_srs = "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs",
         t_srs = "+proj=longlat +datum=WGS84 +no_defs", r = "near", overwrite = T)

我得到了这个结果:

Result of the code I made

从图片中可以看出,它看起来像弯曲的。

预期结果是这样,它是由ENVI生成的

ENVI result

我想知道我做错了什么或需要添加什么才能在R中获得相同的ENVI结果。 我是这个社区的新手,感谢您的帮助

0 个答案:

没有答案