Перепроектировать данные MODIS, используя R (приводит к NA или нет пространственного экстента) - PullRequest
Новая
       113

Перепроектировать данные MODIS, используя R (приводит к NA или нет пространственного экстента)

0 голосов
/ 28 ноября 2018

Я использую данные альбедо GLASS, хранящиеся здесь для данных до 2000 года (AVHRR) и здесь для данных после 2000 года (MODIS).Моя конечная цель - создать растровый стек каждого месяца, который будет содержать данные альбедо белого неба за 1982-2015 годы.Проблема, с которой я столкнулся, заключается в том, что данные MODIS и AVHRR находятся в разных системах пространственной привязки, и я не могу перепроектировать их в одну и ту же систему.

Я конвертирую из hdf в tif, используя R какэто: 

fileavhrr <- ".../GLASS02B05.V04.A1990161.2018062.hdf"
filemodis<-".../GLASS02B06.V04.A2013169.2017128.hdf"
gdal_translate(get_subdatasets(filemodis)[10], dst_dataset =
        ".../modis.tif") 
gdal_translate(get_subdatasets(fileavhrr)[8], projwin = c(-180,90,180,50), dst_dataset = ".../avhrr.tif") #ideally I'd only like data north of 50 degrees

avhrr<- raster(".../avhrr.tif")

#class       : RasterLayer 
#dimensions  : 800, 7200, 5760000  (nrow, ncol, ncell)
#resolution  : 0.05, 0.05  (x, y)
#extent      : -180, 180, 50, 90  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=longlat +ellps=clrk66 +no_defs 
#values      : -32768, 32767  (min, max)

modis<- raster(".../modis.tif")

#class       : RasterLayer   
#dimensions  : 3600, 7200, 25920000  (nrow, ncol, ncell) 
#resolution  : 154.4376, 308.8751  (x, y)  
#extent   : -20015109, -18903159, 8895604, 10007555  (xmin, xmax, ymin, ymax)  
#coord. ref. : +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181
    +b=6371007.181 +units=m +no_defs   
#values      : -32768, 32767  (min, max)

Вот что я пробовал:

1.) Воспользуйтесь MODIS Reprojection Tool .По какой-то причине этот инструмент считает, что подмножества файлов MODIS .hdf - это только один фрагмент (самый верхний левый фрагмент, элемент 0,0), а не глобальный набор данных.Насколько я понимаю, данные MODIS являются глобальными (не в виде плиток?), Поэтому я не знаю, почему MRT делает это.

2.) Используйте растровый пакет в R. 

projectedMODIS <- projectRaster(modis,avhrr,method="bilinear")

Возвращает растр со значениями, которые все NA: 

class       : RasterLayer  
dimensions  : 800, 7200, 5760000  (nrow,> ncol, ncell) 
resolution  : 0.05, 0.05  (x, y) 
extent      : -180, 180,> 50, 90  (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +ellps=clrk66 +no_defs  
values      : NA, NA  (min, max)

3.) Используйте пакет gdalUtils в R: 

gdalwarp(srcfile=get_subdatasets(filemodis)[10], dstfile= ".../gdalMODIS_avhrr.tif", s_srs = crs(modis), t_srs =crs(avhrr) )

Возвращает растр, по существу, без пространственногопротяженность. 

gdalMODISavhrr<-raster(".../gdalMODIS_avhrr.tif")
#class       : RasterLayer 
#dimensions  : 357, 12850, 4587450  (nrow, ncol, ncell)
#resolution  : 0.02801551, 0.02801573  (x, y)
#extent      : -180, 179.9993, 79.99838, 90  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=longlat +ellps=clrk66 +no_defs 
#values      : -32768, 32767  (min, max)

Есть идеи, почему перепроектировать эти данные MODIS так сложно?

1 Ответ

0 голосов
/ 29 ноября 2018

Я не пробовал этого, но, посмотрев на пакет gdalUtils, функция gdalwarp() может сделать то, что вам нужно?

...