Ошибка загрузки изображений MODIS с помощью runGdal - PullRequest
0 голосов
/ 07 сентября 2018

все ~

Я безуспешно пытался в R загрузить изображения MODIS моих достопримечательностей, но всегда получаю одну и ту же ошибку. Кто-нибудь может мне помочь?

Заранее спасибо ~

Вот код:

> product <- "MOD13Q1"
> h = 26
> v = 4
> getT <- getTile(tileH = 26, tileV = 4)
> runGdal(job = "NDVI", product = product, begin = dates2$beginDOY, end = dates2$endDOY, extent = getT)
No collection specified, getting the newest for MOD13Q1 
########################
outProj          =  asIn 
pixelSize        =  asIn 
resamplingType   =  near 
Output Directory =  D:/NDVI_Processed//Dali 
########################
Downloading structure on 'LAADS' for: MOD13Q1.006
Downloading structure from 'LAADS'-server! Try: 1 
Downloading structure of '0' from 'LAADS'-server.                        

Error in if (FtpDayDirs[1] == FALSE) { : 
  missing value where TRUE/FALSE needed

До этой ошибки я потратил некоторое время на проблему с rgdal и решил ее, установив OSGeo4W и указав библиотеку rgdal. Вот моя сессияInfo:

> sessionInfo()
R version 3.5.1 (2018-07-02)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 17134)

Matrix products: default

locale:
[1] LC_COLLATE=Chinese (Simplified)_China.936 
[2] LC_CTYPE=Chinese (Simplified)_China.936   
[3] LC_MONETARY=Chinese (Simplified)_China.936
[4] LC_NUMERIC=C                              
[5] LC_TIME=Chinese (Simplified)_China.936    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] MODIS_1.1.2   raster_2.6-7  mapdata_2.3.0 maps_3.3.0    rgdal_1.3-4  
[6] sp_1.3-1     

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.18    compiler_3.5.1  later_0.7.3     nloptr_1.0.4   
 [5] bitops_1.0-6    class_7.3-14    tools_3.5.1     digest_0.6.16  
 [9] memoise_1.1.0   lattice_0.20-35 shiny_1.1.0     DBI_1.0.0      
[13] parallel_3.5.1  spData_0.2.9.3  e1071_1.7-0     withr_2.1.1    
[17] rgeos_0.3-28    devtools_1.13.6 classInt_0.2-3  mapedit_0.3.2  
[21] grid_3.5.1      sf_0.6-3        R6_2.2.2        ptw_1.9-13     
[25] XML_3.98-1.16   foreign_0.8-70  magrittr_1.5    units_0.6-0    
[29] promises_1.0.1  htmltools_0.3.6 maptools_0.9-3  mime_0.5       
[33] xtable_1.8-3    httpuv_1.4.5    RCurl_1.95-4.11

И, мои MODISoptions

> MODISoptions()

STORAGE:
_______________
localArcPath : D:/NDVI/ 
outDirPath   : D:/NDVI_Processed/ 


DOWNLOAD:
_______________
MODISserverOrder : LAADS, LPDAAC 
dlmethod         : auto 
stubbornness     : high 
wait             : 0.5 
quiet            : FALSE 


PROCESSING:
_______________
GDAL           : GDAL 2.2.4, released 2018/03/19 
MRT            : Version 4.1 (March 2011) 
pixelSize      : asIn 
outProj        : asIn 
resamplingType : NN 
dataFormat     : GTiff 

1 Ответ

0 голосов
/ 12 сентября 2018

Не могли бы вы попробовать еще раз, используя MODIS_1.1.3 (или выше), который теперь доступен из CRAN ? Пакет получил исчерпывающее обновление, в котором, в частности, рассматриваются проблемы, связанные с загрузкой данных из LAADS и LP DAAC (для последнего требуется указание учетных данных для входа в Earthdata). С этой версией пакета ваш код работает на моем компьютере без проблем.

library(MODIS)

getT = getTile(tileH = 26, tileV = 4)
tfs = runGdal(job = "NDVI", product = "MOD13Q1", begin = "2000.02.18", end = "2000.03.10", extent = getT)
tfs
# $`MOD13Q1.006`
# $`MOD13Q1.006`$`2000-02-18`
# [1] ".../MODIS_ARC/PROCES~1/NDVI/MOD13Q1.A2000049.250m_16_days_NDVI.tif"                     
# [2] ".../MODIS_ARC/PROCES~1/NDVI/MOD13Q1.A2000049.250m_16_days_EVI.tif"                      
# (etc.)       
# 
# $`MOD13Q1.006`$`2000-03-05`
# [1] ".../MODIS_ARC/PROCES~1/NDVI/MOD13Q1.A2000065.250m_16_days_NDVI.tif"                     
# [2] ".../MODIS_ARC/PROCES~1/NDVI/MOD13Q1.A2000065.250m_16_days_EVI.tif"                      
# (etc.)

Пожалуйста, ознакомьтесь с моими ответами на следующие вопросы

, которые включают дополнительную информацию, которая может иметь отношение к вам (намеренно опущена здесь, чтобы избежать повторения).

...