Hello,
I've been working on an image created with spplot and with help have
gotten it very close to the final product. One final modification I
would like to make involves filling in grid cells with values less
than the minimum value on the legend with the color corresponding to
the bottom of the legend. In the attached figure, the minimum legend
value is 10, but there a few spots with values as low as 0.01. These
cells are currently blank/white. I would like to avoid including this
value on my legend due to the color spread on the plot becoming
homologous, that is, I would like the color spread to be maximized
between 10-1000. Were I to go from 0.01-1000, the range I'm most
interested in (10-1000) would be confined to a more narrow range of
the legend. I have played around with trying to specify a different
number of cuts for the legend vs plot area, but this causes a
disconnect between the legend and the colors displayed for the
different grid cells. The spplot line looks like this (image attached):
pts<-list("sp.points",K.dat,pch=3,col="black",labels=as.character(K.dat$Kh))
cuts<-c(10,20,30,40,50,60,70,80,90,100,200,300,400,500,600,700,800,900,1000)
spplot(lzm.krige.dir["var1.pred"],at=cuts,colorkey=list(at=log10(cuts),labels=as.character(cuts)),scales=list(draw=TRUE),
xlab="Easting",ylab="Northing",key.space="right",cex=1.1,col.regions=terrain.colors(22),main="Hydraulic
Conductivity of Layer 2",sp.layout=list(pts))
Respectfully,
Eric Morway
------------------------------------------------------------------------
------------------------------------------------------------------------
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
