Suite

Comment télécharger les données MODIS à partir de R à l'aide de la bibliothèque MODISTools


Je souhaite télécharger le produit MODIS MOD13Q1 pour utiliser les couches des ensembles de données scientifiques (SDS) de "250m_16_days_EVI", "250m_16_days_pixel_reliability". La bibliothèque MODISTools du logiciel R permet de télécharger les données MODIS à partir du centre de données LP DAAC.

j'ai lu leMODISOutilspdf guide qui fournit la compréhension de base de la façon d'utiliser la bibliothèque pour sous-ensemble les couches de données MOD13Q1 en utilisant l'emplacement des points. La fonction de sous-ensemble suivante que je souhaite utiliser :

library('MODISTools') MODISSoussets(LoadDat = modis.subset, Products = "MOD13Q1", Bands = c("250m_16_days_EVI", "250m_16_days_pixel_reliability"), Size = c(1,1))

IciLoadDat = modis.subsetest le bloc de données des points de latitude et de longitude au format décimal et selonTaille = c(1,1)paramètre qui est utilisé pour extraire les 1 Km à gauche et 1 Km à droite du point focal comme étendue. Mais je veux utiliser l'étendue de la zone d'étude entière comme

lat <- c(23.6947, 37.0894) long <- c(60.8786, 79.3074) period <- data.frame(lat=lat,long=long, start.date=2013,end.date=2014,id=1) MODISSous-ensembles (LoadDat = période, Produits = "MOD13Q1", Bandes = "250m_16_days_EVI", "250m_16_days_pixel_reliability", Taille = c(0,0), TimeSeriesLength = 1)

Mais ici, je suis confus quant à l'utilisation du paramètre Size qui est utilisé pour extraire les valeurs des points. Comment puis-je l'utiliser pour extraire toutes les données disponibles dans cette mesure ? Et deuxièmement, puis-je utiliser ces informations de point dans un fichier ASCII pour afficher un raster dans le logiciel R ? Je ne veux pas télécharger la vignette complète qui contient 250 Mo et dans mon cas, ce sera 250 Go de données si je télécharge toutes les données disponibles à partir de 2000 à 2015.


Il semble que les données de taille seront utilisées pour sélectionner le pixel environnant et autour du pixel central de votre choix.

Size = c(1,1) vous donnera la sortie non seulement du pixel central (ou du point) mais aussi de 8 pixels voisins.

Donc, si vous souhaitez télécharger les données uniquement pour le point concerné, donnez la valeur 'Taille' sous la forme c(0,0). Voir les vignettes pour plus de détails.

Maintenant comme pour télécharger les données, vous pouvez créer une matrice des coordonnées du raster. Ensuite, utilisez-le peut-être pour télécharger toutes les données nécessaires.