[R-sig-Geo] raster: crop error ('nrows > 0 is not TRUE')

[Dominik Schneider]

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.zambrano at 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.zambrano at 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 at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

	[[alternative HTML version deleted]]