[R-sig-Geo] Question with spplot for publication

[Ziyu Ma]

Thanks Jon and Roger!

The panel function worked effortlessly!
And I learned a lot from your codes.

Cheers,
Ma

Ziyu Ma
PhD Student
Ecoinformatics & Biodiversity,
Department of Bioscience, Aarhus University
Ny Munkegade 114, DK-8000 Aarhus C, Denmark



On Fri, Aug 22, 2014 at 9:45 AM, Jon Skoien <jon.skoien at jrc.ec.europa.eu>
wrote:

> Hi Ma,
>
> You should be able to plot the polygon on top of your raster with a panel
> function. Here is an example with the meuse data:
>
> r <- raster(system.file("external/test.grd", package="raster"))
> s <- stack(r, r*2)
> names(s) <- c('meuse', 'meuse x 2')
> spplot(s)
> data(meuse.riv)
>
> meuse.sr = SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"
> meuse.riv")))
>
> spplot(s,
>   panel = function(x, y, ...) {
>     panel.gridplot(x, y, ...)
>     sp.polygons(meuse.sr)
>   }
> )
>
> Is this more or less how you want it?
>
> Cheers,
> Jon
>
>
>
> On 8/22/2014 9:15 AM, Ziyu Ma wrote:
>
>> Dear all,
>>
>> I miss you guys from GEOSTAT Bergen!
>>
>> Here's my question: I want to make a map directly from R for my
>> manuscript,
>> using a raster layer overlaid with spatial polygons as outlines.
>>
>> Please see the attachment for what I mean: http://1drv.ms/1sYYqPg
>>
>> I did the upper plot in the attachement using these codes:
>>
>> NRI = stack(NRI.A,NRI.G)  # NRI.A and NRI.G are raster layer names
>> names(NRI) = c("Angiosperms NRI", "Gymnosperms NRI")
>> spplot(NRI)
>>
>>
>> I need to put two rasters together in spplot because I want to use the
>> same
>> legend scale.
>>
>> I also want to add the outlines, so that it can look like the Photoshopped
>> plot in the attachment.
>>
>> Normally, if not using spplot, I can do this:
>>
>> plot(NRI.G) # NRI.G is a raster layer
>> plot(Outlines)  # Outlines is a SpatialPolygons object
>>
>>
>> But how can I do it for both panels in spplot? I think the panels and
>> titles in spplot are very neat, so it is preferred if there is a straight
>> forward way to tweak the settings.
>>
>> Or should I use other functions than spplot? (e.g. try coding the color
>> scheme of the 2 plots and place them side by side?)
>>
>> Thank you for the help. I really appreciate those who made R able to
>> produce publication quality plots.
>>
>> Cheers,
>> Ma​
>>
>> Ziyu Ma
>> PhD Student
>> Ecoinformatics & Biodiversity,
>> Department of Bioscience, Aarhus University
>> Ny Munkegade 114, DK-8000 Aarhus C, Denmark
>>
>>         [[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
>>
>
>
>

	[[alternative HTML version deleted]]