Plotting median values globally

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

Plotting median values globally

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

Re: Plotting median values globally

Vijay Lulla
Did you try `levelplot(raster(RCP1pctCO2median, 61))`?

On Tue, Oct 29, 2019 at 12:17 AM rain1290--- via R-sig-Geo <
[hidden email]> 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
> [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: Plotting median values globally

R-sig-geo mailing list
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 <[hidden email]>
To: rain1290 <[hidden email]>
Cc: R-sig-geo Mailing List <[hidden email]>
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 <[hidden email]> 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
[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: Plotting median values globally

Bede-Fazekas Ákos
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 <[hidden email]>
> To: rain1290 <[hidden email]>
> Cc: R-sig-geo Mailing List <[hidden email]>
> 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 <[hidden email]> 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
> [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
>

_______________________________________________
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: Plotting median values globally

R-sig-geo mailing list
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 <[hidden email]>
To: r-sig-geo <[hidden email]>
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 <[hidden email]>
> To: rain1290 <[hidden email]>
> Cc: R-sig-geo Mailing List <[hidden email]>
> 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 <[hidden email]> 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
> [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
>

_______________________________________________
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