Thank you very much Dominik. You are right, the southern border of the shapefile used as boundary is outside the borders of the RasterStack that I want to crop, and that is the difference with a previous case where the 'crop' command worked without any problem.
However, reading the documentation of 'crop', I think that command should rise a warning and not an error in such a case: "Areas included in ‘y’ but outside the extent of ‘x’ are ignored (see ‘extend’ if you want a larger area)." Thanks for your help. mzb ===================================== "To accomplish great things, we must not only act, but also dream; not only plan, but also believe" (Anatole France) ===================================== Linux user #454569 -- Linux Mint user On 3 October 2016 at 12:34, Dominik Schneider <dosc3...@colorado.edu> wrote: > Isn't the southern border of your boundary shapefile outside the extent of > your raster? hence there are no rows to extract. > > On Mon, Oct 3, 2016 at 7:45 AM, Mauricio Zambrano Bigiarini > <mauricio.zambr...@ufrontera.cl> wrote: >> >> Dear list, >> >> I'm applying the 'crop' command of the raster pacakge to subset a >> RasterStack ('r') with a SpatialPolygonsDataFrame ('boundary'): >> >> raster.crop <- crop(r, boundary, snap="out") >> >> and what I get is: >> >> Error: nrows > 0 is not TRUE >> >> >> The summary of the two previous objects are: >> >> > r >> class : RasterStack >> dimensions : 666, 211, 140526, 8 (nrow, ncol, ncell, nlayers) >> resolution : 0.05, 0.05 (x, y) >> extent : -76.3, -65.75, -50, -16.7 (xmin, xmax, ymin, ymax) >> coord. ref. : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs >> +towgs84=0,0,0 >> names : Year.2003, Year.2004, Year.2005, Year.2006, Year.2007, >> Year.2008, Year.2009, Year.2010 >> min values : 0, 0, 0, 0, 0, >> 0, 0, 0 >> max values : 844.79, 706.53, 616.73, 888.75, 701.53, >> 853.47, 642.07, 986.31 >> >> >> > boundary >> class : SpatialPolygonsDataFrame >> features : 16 >> extent : -75.693, -66.419, -55.914, -17.498 (xmin, xmax, ymin, ymax) >> coord. ref. : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs >> +towgs84=0,0,0 >> variables : 1 >> names : NOM_REG >> min values : 1 >> max values : 9 >> >> the previous 'crop' command worked perfectly with a different >> RasterStack object, with a similar spatial extent, and I don't now >> what is wrong here. >> >> >> The output of my sessionInfo() is: >> >> sessionInfo() >> R version 3.3.1 (2016-06-21) >> Platform: x86_64-pc-linux-gnu (64-bit) >> Running under: Ubuntu 14.04.2 LTS >> >> locale: >> [1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C >> LC_TIME=en_GB.UTF-8 LC_COLLATE=en_GB.UTF-8 >> LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_GB.UTF-8 >> LC_PAPER=en_US.UTF-8 LC_NAME=C >> [9] LC_ADDRESS=C LC_TELEPHONE=C >> LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C >> >> attached base packages: >> [1] stats graphics grDevices utils datasets methods base >> >> other attached packages: >> [1] colorspace_1.2-6 reshape2_1.4.1 rasterVis_0.40 >> latticeExtra_0.6-28 RColorBrewer_1.1-2 lattice_0.20-34 >> rgdal_1.1-10 hydroGOF_0.3-8-7 hydroTSM_0.4-8 xts_0.9-7 >> zoo_1.7-13 >> [12] raster_2.5-8 sp_1.2-3 >> >> loaded via a namespace (and not attached): >> [1] Rcpp_0.12.6 magrittr_1.5 maptools_0.8-39 >> viridisLite_0.1.3 FNN_1.1 stringr_1.0.0 plyr_1.8.4 >> tools_3.3.1 parallel_3.3.1 grid_3.3.1 gstat_1.1-3 >> e1071_1.6-7 >> [13] class_7.3-14 intervals_0.15.1 automap_1.0-14 ncdf4_1.15 >> stringi_1.1.1 spacetime_1.1-5 reshape_0.8.5 >> foreign_0.8-66 hexbin_1.27.1 >> >> >> Any idea about how to solve this issue is highly appreciated. >> >> Thanks in advance, >> >> >> >> Mauricio Zambrano-Bigiarini, PhD >> >> ===================================== >> Department of Civil Engineering >> Faculty of Engineering and Sciences >> Universidad de La Frontera >> PO Box 54-D, Temuco, Chile >> http://hzambran.github.io/ >> ===================================== >> mailto : mauricio.zambr...@ufrontera.cl >> work-phone : +56 45 259 2812 >> ===================================== >> "To accomplish great things, we must not only act, >> but also dream; not only plan, but also believe" >> (Anatole France) >> ===================================== >> Linux user #454569 -- Linux Mint user >> >> -- >> La información contenida en este correo electrónico y cualquier anexo o >> respuesta relacionada, puede contener datos e información confidencial y >> no >> puede ser usada o difundida por personas distintas a su(s) >> destinatario(s). >> Si usted no es el destinatario de esta comunicación, le informamos que >> cualquier divulgación, distribución o copia de esta información constituye >> un delito conforme a la ley chilena. Si lo ha recibido por error, por >> favor >> borre el mensaje y todos sus anexos y notifique al remitente. >> >> _______________________________________________ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > -- La información contenida en este correo electrónico y cualquier anexo o respuesta relacionada, puede contener datos e información confidencial y no puede ser usada o difundida por personas distintas a su(s) destinatario(s). Si usted no es el destinatario de esta comunicación, le informamos que cualquier divulgación, distribución o copia de esta información constituye un delito conforme a la ley chilena. Si lo ha recibido por error, por favor borre el mensaje y todos sus anexos y notifique al remitente. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo