Я очень новичок в R. Мне удалось создать карту ниже с кодом ниже.окраски полигонов в ggplot2 (R)
Я хочу, чтобы цвет самых правый многоугольника на карте (обозначаемую «WestTerr» в коде) другой цвет, чем остальные три многоугольников (обозначаемого «Серены» в коде). Я пробовал разные вещи, которые на меня набросились. Кто-нибудь знает, как мне помочь?
library(ggmap)
library(RgoogleMaps)
library(foreign)
library(raster)
library(sp)
library(rgeos)
library(spatstat)
library(maptools)
library(rgdal)
library(ggplot2)
CenterOfMap <- geocode("-1.4, 35.08")
MasaiMara <- get_map(c(lon=CenterOfMap$lon, lat=CenterOfMap$lat), zoom = 10,
maptype = "roadmap", source = "google")
MasaiMaraMap <- ggmap(MasaiMara)
setwd("Territories_Jenna/")
SouthTerr <- readOGR(".","SouthMCP")
SouthTerr <- spTransform(SouthTerr, CRS("+proj=longlat +datum=WGS84"))
fortify(SouthTerr)
HZTerr <- readOGR(".","HZMCP")
HZTerr <- spTransform(HZTerr, CRS("+proj=longlat +datum=WGS84"))
fortify(HZTerr)
SandHZ <- gUnion(SouthTerr, HZTerr)
SandHZ <- spTransform(SandHZ, CRS("proj=longlat +datum=WGS84"))
fortify(SandHZ)
MapwithSandHZ <- MasaiMaraMap +
geom_polygon(aes(x=long, y=lat, group=group),
fill = "black", size=.7, color="black",
data=SandHZ, alpha=0)
NorthTerr <- readOGR(".","NorthMCP")
NorthTerr <- spTransform(NorthTerr, CRS("+proj=longlat +datum=WGS84"))
fortify(NorthTerr)
MapwithNorth <- MasaiMaraMap +
geom_polygon(aes(x=long, y=lat, group=group),
fill = "black", size=.7, color="black",
data=NorthTerr, alpha=0)
Serena <- gUnion(SandHZ, NorthTerr)
Serena <- spTransform(Serena, CRS("proj=longlat +datum=WGS84"))
fortify(Serena)
MapwithSerena <- MasaiMaraMap +
geom_polygon(aes(x=long, y=lat, group=group),
fill = "black", size=.7, color="black",
data=Serena, alpha=0)
WestTerr <- readOGR(".","WestMCP")
WestTerr <- spTransform(WestTerr, CRS("+proj=longlat +datum=WGS84"))
fortify(WestTerr)
MapwithWest <- MasaiMaraMap +
geom_polygon(aes(x=long, y=lat, group=group),
fill = "red", size=.7, color="red",
data=WestTerr, alpha=0)
AllTerrs <- gUnion(Serena, WestTerr)
AllTerrs <- spTransform(AllTerrs, CRS("proj=longlat +datum=WGS84"))
fortify(AllTerrs)
MapwithAllTerrs <- MasaiMaraMap +
geom_polygon(aes(x=long, y=lat, group=group),
fill = "red", size=.7, color="red",
data=AllTerrs, alpha=0)
MapwithAllTerrs
А какие ошибки вы получаете? – dg99