[R-sig-Geo] Set color key for stplot

[Frede Aakmann Tøgersen]

Hi Katherine

Yes, it looks like it is lattice plot functions.

I have an example where I wantes to plot the sea as blue and the land according to terrain.colors()

Here is my call to levelplot:

levelplot( Z ~ X + Y, data = xyz, labels = FALSE, aspect = "iso", contour = FALSE,
          cuts = 16, , at = c(-Inf,do.breaks(range(xyz$Z), 15)),
          col.regions = c("#0000FFFF", terrain.colors(15)),
          subset = X > 370000 & X < 380000 & Y > 670000 & Y < 680000)

and a png is attached.

As you can see you need to combine the cuts and col.regions argument to do what you want.



Yours sincerely / Med venlig hilsen


Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling

Technology & Service Solutions
T +45 9730 5135
M +45 2547 6050
frtog at vestas.com
http://www.vestas.com

Company reg. name: Vestas Wind Systems A/S
This e-mail is subject to our e-mail disclaimer statement.
Please refer to www.vestas.com/legal/notice
If you have received this e-mail in error please contact the sender. 


> -----Original Message-----
> From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-
> project.org] On Behalf Of Katharine Walter
> Sent: 14. maj 2014 18:53
> To: r-sig-geo at r-project.org
> Subject: [R-sig-Geo] Set color key for stplot
> 
> Hi,
> Is there a way to use a color ramp when using stplot and simultaneously
> assign a specific color for a single value?  For example, I am trying to
> plot disease cases and would like 0s to be plotted as grey and the rest of
> the cases to be colored by a color ramp.
> 
> #this plots cases including 0s as a color ramp:
> stplot(stfdf[,,"lyme_cases"], t, col.regions = brewer.pal(9, "YlOrRd"),
> cuts=9 )[28]
> cuts=c(1,seq(20,160,by=20), 186)
> 
> I have tried defining cuts so that any polygon with less than 1 case is
> colored grey, but this does not work and instead plots any polygon with <20
> cases as grey.
> 
> stplot(stfdf[,,"lyme_cases"], t, col.regions = c("grey", brewer.pal(9,
> "YlOrRd")), cuts=c(0.1,seq(20,160,by=20), 186) )[28]
> 
> Thank you very much for your help!
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
-------------- next part --------------
A non-text attachment was scrubbed...
Name: seaANDland.png
Type: image/png
Size: 6981 bytes
Desc: seaANDland.png
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20140514/8ce17b12/attachment.png>