[R-sig-Geo] can spplot add features ?

Edzer Pebesma edzer.pebesma at uni-muenster.de
Tue Jul 7 11:16:29 CEST 2009



Patrick Giraudoux wrote:
> Edzer Pebesma a écrit :
>> Patrick Giraudoux wrote:
>>  
>>> Dear listers,
>>>
>>> I would like to add the border of  a country (SpatialPolygons object
>>> called ChinaBorder) in a spplot built as following:
>>>
>>> mypal<-colorRampPalette(c("white", "cyan",
>>> "blue","green","yellow","red"))
>>> statSpat<-SpatialPointsDataFrame(chinagrid,data.frame(n=z))
>>> gridded(statSpat)<-TRUE
>>> spplot(statSpat,col.regions=mypal)
>>>
>>> I cannot managed to overlap ChinaBorder to statSpat with the 'panel'
>>> argument of spplot.
>>>
>>> Any idea about how to proceed ?
>>>     
>> Patrick, you might want to try
>>
>> spplot(statSpat, col.regions = mypal, sp.layout = list("sp.polygons",
>> ChinaBorder))
>>  
>>> P
>
>
> Exactly what I need (the way to state was unclear to me in the doc),
> except the image is overlapping the border when I would like to get
> the reverse (the contour over the image). Is there a way to make it ?
>
> Cheers,
>
> Patrick
>
Try

spplot(statSpat, col.regions = mypal, sp.layout = list("sp.polygons",
ChinaBorder, first = FALSE))

or maybe TRUE.
--
Edzer

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de/
http://www.springer.com/978-0-387-78170-9 e.pebesma at wwu.de



More information about the R-sig-Geo mailing list