raster: crop error ('nrows > 0 is not TRUE')

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

raster: crop error ('nrows > 0 is not TRUE')

MAURICIO ZAMBRANO BIGIARINI
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     : [hidden email]
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
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: 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 <
[hidden email]> 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     : [hidden email]
> 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
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: raster: crop error ('nrows > 0 is not TRUE')

MAURICIO ZAMBRANO BIGIARINI
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 <[hidden email]> 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
> <[hidden email]> 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     : [hidden email]
>> 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
>> [hidden email]
>> 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
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo