-
Notifications
You must be signed in to change notification settings - Fork 27
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Global extent GDAL mosaic problem #54
Comments
Hi I am just guessing now, but as you do global coverage, try without
extent at all.
|
With no extent given it opens an interactive map to pick tiles. |
This seems to be related to tfs1 = runGdal(product = "MOD13A3.006",
extent = raster::extent(-180, 180, -90, 90),
begin = "2016.01.01", end = "2016.01.01",
SDSstring = "01", job = "global-test-ndvi",
outProj = "+init=epsg:4326", overwrite = TRUE)
raster(unlist(tfs1))
# class : RasterLayer
# dimensions : 21600, 43210, 933336000 (nrow, ncol, ncell)
# resolution : 0.008331493, 0.008333333 (x, y)
# extent : -180, 180.0038, -90, 90 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
# data source : /.../MODIS_ARC/PROCESSED/global-test-ndvi/MOD13A3.A2016001.1_km_monthly_EVI.tif
# names : MOD13A3.A2016001.1_km_monthly_EVI
# values : -32768, 32767 (min, max) However, I wonder why |
I'm using v1.1.3 with GDAL 2.1.3, I don't have MRT. Trying to produce a global MODIS landcover raster using the following code:
All tiles are downloaded correctly and everything appears to work fine, but when I look at the raster it's the correct height by only 1 cell wide:
The text was updated successfully, but these errors were encountered: