[R-sig-Geo] rotate-ing a raster brick?

Andy Hartley andrew.hartley at metoffice.gov.uk
Thu Jan 12 10:47:35 CET 2012


Hi Robert,
 
Thanks for replying. Unfortunately I can't change my version of R (we
have a centrally installed version that gets rolled out every so often).
I think raster is up to date though. Here's the info you asked for - the
traceback doesn't make much sense to me!
 
I've already done a work around (looping through each band, and stacking
at the end), but it would be nice to understand why this is happening.
 
Cheers,
Andy
 
 
> sessionInfo()
R version 2.13.1 (2011-07-08)
Platform: x86_64-unknown-linux-gnu (64-bit)
 
locale:
[1] C
 
attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

 
other attached packages:
[1] ncdf_1.6.6    raster_1.9-52 sp_0.9-93    
 
loaded via a namespace (and not attached):
[1] grid_2.13.1     lattice_0.19-30 tools_2.13.1   
> rbb <- rotate(rb)
Error in rowSums(is.na(x)) : 
  'x' must be an array of at least two dimensions
> traceback()
10: stop("'x' must be an array of at least two dimensions")
9: rowSums(is.na(x))
8: .local(x, y, ...)
7: merge(x, tolerance = tolerance, filename = filename, format = format,

       datatype = datatype, overwrite = overwrite, progress = progress, 
       ext = ext, overlap = overlap, check = FALSE)
6: merge(x, tolerance = tolerance, filename = filename, format = format,

       datatype = datatype, overwrite = overwrite, progress = progress, 
       ext = ext, overlap = overlap, check = FALSE)
5: .local(x, y, ...)
4: merge(r1, r2, ...)
3: merge(r1, r2, ...)
2: rotate(rb)
1: rotate(rb)

 

-- 
Andrew Hartley Climate Impacts Risk Analyst
Met Office FitzRoy Road Exeter EX1 3PB United Kingdom
Tel: +44 (0)1392 885720 Fax: +44 (0)1392 885681
E-mail: andrew.hartley at metoffice.gov.uk http://www.metoffice.gov.uk
<blocked::http://www.metoffice.gov.uk/> 

Join us in celebrating 150 years of forecasting for the nation and enter
our photo competition
http://www.metoffice.gov.uk/about-us/who/our-history
<blocked::http://www.metoffice.gov.uk/about-us/who/our-history>  

 

________________________________

From: Robert Hijmans [via R-sig-geo]
[mailto:ml-node+s2731867n7179261h67 at n2.nabble.com] 
Sent: 12 January 2012 09:18
To: Hartley, Andrew
Subject: Re: rotate-ing a raster brick?


Andrew, 
Strange, the same thing (for the same file) works for me. Could you send

your sessionInfo() and a traceback() after the error occurs ( that
should 
reveal where the call to rowSums(is.na(x)) occurs (I have no idea) that
the 
error message complains about ) 
Robert 

On Tue, Jan 10, 2012 at 4:32 AM, Hartley, Andrew < 
[hidden email]> wrote: 


> Hi list, 
> 
> First, a big thank you to the people behind the development of the 
> "raster" package - I find it very useful, especially with the addition

> of the Z dimension. Now, a question - I've imported climate data into
R 
> from the CMIP5 project. All the data is in NetCDF format, which is
much 
> easier to handle now (big thanks for that!). I find the 'rotate' 
> function especially useful because most of these datasets seem to be 
> between 0 and 360 degrees. 
> 
> I've used rotate successfully on one band of the imported data, but it

> throws up an error with bricks or stacks. See the code below. Is this
a 
> known error, or am I doing something wrong? 
> 
> Cheers, 
> Andy 
> 
> Example file downloaded from 
>
http://pcmdi3.llnl.gov/esgcet/dataset/cmip5.output1.CCCma.CanESM2.esmHis
> torical.mon.atmos.Amon.r1i1p1.html 
> 
>
f="CMIP5/CCCma/CanESM2/historical/mon/atmos/Amon/r1i1p1/v20111028/pr/pr_

> Amon_CanESM2_historical_r1i1p1_185001-200512.nc" 
> 
> rb <- brick(f, crs="+proj=longlat +datum=WGS84") 
> rb 
> class       : RasterBrick 
> dimensions  : 64, 128, 8192, 1872  (nrow, ncol, ncell, nlayers) 
> resolution  : 2.8125, 2.789327  (x, y) 
> extent      : -1.40625, 358.5938, -89.25846, 89.25846  (xmin, xmax, 
> ymin, ymax) 
> coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
> values      : /media/WD 
>
MYBOOK/CMIP5/CCCma/CanESM2/historical/mon/atmos/Amon/r1i1p1/v20111028/pr

> /pr_Amon_CanESM2_historical_r1i1p1_185001-200512.nc 
> zvar        : pr 
> rbr <- rotate(rb) 
> Error in rowSums(is.na(x)) : 
>  'x' must be an array of at least two dimensions 
> 
> r <- subset(rb, 1) 
> rr <- rotate(r) 
> rr 
> class       : RasterLayer 
> dimensions  : 64, 128, 8192  (nrow, ncol, ncell) 
> resolution  : 2.8125, 2.789327  (x, y) 
> extent      : -181.4062, 178.5938, -89.25846, 89.25846  (xmin, xmax, 
> ymin, ymax) 
> coord. ref. : NA 
> values      : in memory 
> min value   : 0 
> max value   : 0.0003839434 
> -- 
> Andrew Hartley Climate Impacts Risk Analyst 
> Met Office FitzRoy Road Exeter EX1 3PB United Kingdom 
> Tel: +44 (0)1392 885720 Fax: +44 (0)1392 885681 
> E-mail: [hidden email] http://www.metoffice.gov.uk
> <blocked::http://www.metoffice.gov.uk/> 
> 
> Join us in celebrating 150 years of forecasting for the nation and
enter 
> our photo competition 
> http://www.metoffice.gov.uk/about-us/who/our-history
> <blocked::http://www.metoffice.gov.uk/about-us/who/our-history> 
> 
> 
>        [[alternative HTML version deleted]] 
> 
> _______________________________________________ 
> 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



________________________________

If you reply to this email, your message will be added to the discussion
below:
http://r-sig-geo.2731867.n2.nabble.com/rotate-ing-a-raster-brick-tp71718
58p7179261.html 
To start a new topic under R-sig-geo, email
ml-node+s2731867n2731867h53 at n2.nabble.com 
To unsubscribe from R-sig-geo, click here
<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=u
nsubscribe_by_code&node=2731867&code=YW5kcmV3LmhhcnRsZXlAbWV0b2ZmaWNlLmd
vdi51a3wyNzMxODY3fC0xNTc1OTEyMTQw> .
NAML
<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=m
acro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.names
paces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.nam
l.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nab
ble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamesp
ace-nabble.view.web.template.InstantMailNamespace&breadcrumbs=instant+em
ails%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_ins
tant_email%21nabble%3Aemail.naml>  


--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/rotate-ing-a-raster-brick-tp7171858p7179691.html
Sent from the R-sig-geo mailing list archive at Nabble.com.



More information about the R-sig-Geo mailing list