我正在使用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)
我得到了这个结果:
从图片中可以看出,它看起来像弯曲的。
预期结果是这样,它是由ENVI生成的
我想知道我做错了什么或需要添加什么才能在R中获得相同的ENVI结果。 我是这个社区的新手,感谢您的帮助