if (FALSE) {
Vector_surveys_NGA_NG <- getVecOcc(country = c("Nigeria", "Niger"))
autoplot(Vector_surveys_NGA_NG)
# Download the predicted distribution of An. dirus species complex Raster and
# vector points for Myanmar and visualise these together on a map.
# Download Myanmar shapefile to use for raster download.
MMR_shp <- getShp(ISO = "MMR", admin_level = "admin0")
# Download An. dirus predicted distribution Raster & plot this
MMR_An_dirus <- getRaster(surface = "Anopheles dirus species complex", shp = MMR_shp)
p <- autoplot(MMR_An_dirus, shp_df = MMR_shp, printed = FALSE)
# Download raw occurrence points & plot these over the top of the raster
species <- getVecOcc(country = "Myanmar", species = "Anopheles dirus")
# p[[1]] +
# geom_point(data = species,
# aes(longitude,
# latitude,
# colour = species))+
# scale_colour_manual(values = "black", name = "Vector suvery locations")+
# scale_fill_distiller(name = "Predicted distribution of An. dirus complex",
# palette = "PuBuGn",
# direction = 1)+
# ggtitle("Vector Survey points\n + The predicted distribution of An. dirus complex")
}
Run the code above in your browser using DataLab