Hi Akos and everyone, Thank you so much, as well! Yes, this worked very well - I cannot believe it was just a bracket that was misplaced that caused that! Many thanks, once again, for your help - it was extremely valuable and, thus, much appreciated!
-----Original Message----- From: Bede-Fazekas Ákos <bfalevl...@gmail.com> To: r-sig-geo <r-sig-geo@r-project.org> Sent: Tue, Oct 29, 2019 12:58 pm Subject: Re: [R-sig-Geo] Plotting median values globally Hello, It is because the wrong place of the closing bracket. This might work, since margin, at, cuts etc. are the parameters of levelplot(9 not raster(): levelplot(raster(RCP1pctCO2Median,61), margin=F,at=cutpoints,cuts=11,pretty=TRUE,par.settings=mapTheme, main="Median values") HTH, Ákos 2019.10.29. 17:38 keltezéssel, rain1290--- via R-sig-Geo írta: > Hi Vijay and others, > Thank you so, so much for this reply!! Yes, this worked, just as you > suggested. Just one little detail that I wanted to know - I was trying to > change the color range of my plot, like this: > mapTheme<-rasterTheme(region=rev(brewer.pal(10,"Spectral"))) > cutpoints<-c(200,175,150,125,100,75,50,25,0) >> levelplot(raster(RCP1pctCO2Median,61, >> margin=F,at=cutpoints,cuts=11,pretty=TRUE,par.settings=mapTheme, >> main="Median values")) > > > > But I receive this error: > Error in .local(x, ...) : > unused arguments (margin = F, at = cutpoints, cuts = 11, pretty = TRUE, >par.settings = mapTheme, main = "Median values") > What could be the reason behind this error? > Thank you, once again!! > ~Trav.~ > -----Original Message----- > From: Vijay Lulla <vijaylu...@gmail.com> > To: rain1290 <rain1...@aim.com> > Cc: R-sig-geo Mailing List <r-sig-geo@r-project.org> > Sent: Tue, Oct 29, 2019 9:02 am > Subject: Re: [R-sig-Geo] Plotting median values globally > > Did you try `levelplot(raster(RCP1pctCO2median, 61))`? > > On Tue, Oct 29, 2019 at 12:17 AM rain1290--- via R-sig-Geo > <r-sig-geo@r-project.org> wrote: > > Hi there, > I am trying to create a global plot with median precipitation values. What I > would like to do is isolate the "median" values and plot those values on a > global raster map. I already first created a raster stack, called > "RCP1pctCO2median", which contains the median values to be plotted. However, > what I would like to do is simply plot only those values at Year 61 for each > grid cell. There are 138 layers, with each layer representing one year. The > idea would be to plot only those values for Year/layer 61 globally (i,e, for > each grid cell), but I am unsure how to isolate that specific year/layer? > "RCP1pctCO2median" has the following attributes: > > class : RasterStack > dimensions : 64, 128, 8192, 138 (nrow, ncol, ncell, nlayers) > resolution : 2.8125, 2.789327 (x, y) > extent : -181.4062, 178.5938, -89.25846, 89.25846 (xmin, xmax, ymin, > ymax) > coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 > names : layer.1, layer.2, layer.3, layer.4, layer.5, > layer.6, layer.7, layer.8, layer.9, layer.10, layer.11, > layer.12, layer.13, layer.14, layer.15, ... > min values : 0.42964514, 0.43375653, 0.51749371, 0.50838983, 0.45366730, > 0.53099146, 0.49757186, 0.45697752, 0.41382199, 0.46082401, 0.45516687, > 0.51857087, 0.41005131, 0.45956529, 0.47497867, ... > max values : 96.30350, 104.08584, 88.92751, 97.49373, 89.57201, > 90.58570, 86.67651, 88.33519, 96.94720, 101.58247, 96.07792, > 93.21948, 99.59785, 94.26218, 90.62138, ... > > > And this is how I would like to create my global raster map. However, I am > not sure how to plot only one value (i.e. the median value, which occurs at > Year 61) for each grid cell. This is what I have so far, but all that this > does is plot 138 mini plots, and I just want one plot showing the median > values (i.e. Year 61): > > > mapTheme<-rasterTheme(region=rev(brewer.pal(10,"Spectral"))) > cutpoints<-c(200,175,150,125,100,75,50,25,0) > levelplot(RCP1pctCO2median,margin=F,at=cutpoints,cuts=11,pretty=TRUE,par.settings=mapTheme, > main="Median values") > > > > Any help would be greatly appreciated! > I look forward to your response! > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo