autoplot_MAPraster: Quickly visualise Rasters downloaded from MAP

In malariaAtlas: An R Interface to Open-Access Malaria Data, Hosted by the 'Malaria Atlas Project'

Quickly visualise Rasters downloaded from MAP

Description

autoplot_MAPraster is a wrapper for autoplot.MAPraster that calls as.MAPraster to allow automatic map creation for RasterLayer/RasterStack objects downloaded from MAP.

Usage

autoplot_MAPraster(object, ...)

Arguments

object	RasterLayer/RasterStack object to be visualised.
...	other optional arguments to autoplot.MAPraster (e.g. shp_df, legend_title, page_title...)

Value

autoplot_MAPraster returns a list of plots (gg objects) for each supplied raster.

See Also

getRaster:

to download rasters directly from MAP.

as.MAPraster:

to convert RasterLayer/RasterStack objects into a 'MAPraster' object (data.frame) for easy plotting with ggplot.

autoplot.MAPraster:

to quickly visualise MAPraster objects created using as.MAPraster.

Examples

## Not run: 
#Download PfPR2-10 Raster (Bhatt et al 2015) and raw survey points for Madagascar in
#  2013 and visualise these together on a map.

# Download madagascar shapefile to use for raster download.
MDG_shp <- getShp(ISO = "MDG", admin_level = "admin0")

# Download PfPR2-10 Raster for 2013 & plot this
MDG_PfPR2_10 <- getRaster(surface = "Plasmodium falciparum PR2-10", shp = MDG_shp, year = 2013)
# p <- autoplot_MAPraster(MDG_PfPR2_10)

# Download raw PfPR survey points & plot these over the top of the raster
pr <- getPR(country = c("Madagascar"), species = "Pf")
# p[[1]] +
# geom_point(data = pr[pr$year_start==2013,],
#            aes(longitude, latitude, fill = positive / examined, size = examined), shape = 21) +
#   scale_size_continuous(name = "Survey Size") +
#   scale_fill_distiller(name = "PfPR", palette = "RdYlBu") +
#   ggtitle("Raw PfPR Survey points\n + Modelled PfPR 2-10 in Madagascar in 2013")
## End(Not run)


# Download global raster of G6PD deficiency (Howes et al 2012) and visualise this on a map.
## Not run: 
G6PDd_global <- getRaster(surface = "G6PD Deficiency Allele Frequency")
#autoplot_MAPraster(G6PDd_global)

## End(Not run)

malariaAtlas documentation built on Oct. 27, 2023, 9:07 a.m.