Message-ID: <CADvFi5pYBK9rwL4GpmYu+azD7E0UsgZjyWcfq9wABZodfJHkYw@mail.gmail.com>
Date: 2018-11-20T19:15:47Z
From: Hugo Costa
Subject: Doubt: add a boderline map
In-Reply-To: <CALN9TERi67z0QwEzJwF5VhWBHmaFZHRcW=PM5FXPHQBRsQ8anQ@mail.gmail.com>
Hi Lara,
have you tried to run the following?
map <-plot(proj$A_EMmeanByROC_mergedAlgo_mergedRun_mergedData,
main= "", xlab ="x.coords", ylab="y.coords", cex.axis=caxis)
plot(nc.sids, asp=1, main="", add=TRUE)
Hugo
Lara Silva <lara.sfp.silva at gmail.com> escreveu no dia ter?a, 20/11/2018
?(s) 16:04:
> Hello,
>
> I got a map through the plot function.
>
> Code:
>
> map <-plot(proj$A_EMmeanByROC_mergedAlgo_mergedRun_mergedData,
> main= "", xlab ="x.coords", ylab="y.coords", cex.axis=caxis)
>
> However I needed to add a border line in my study area to become more
> perceptible (Figure 1).
>
> I tried the border(), rect () functions but it gave error.
>
> So, I tried to create a border line
>
> Code:
> nc.sids <- readShapePoly("PoligonoSM.shp")
>
> plot(nc.sids, asp=1, main="")
>
> nc.border <- unionSpatialPolygons(nc.sids, rep(1, nrow(nc.sids)))
> plot(nc.border, asp=1, main="")
>
>
> But I am having a lot of difficulties in building the script.
>
> How can cross the information of the 2 plots?
>
> https://www.dropbox.com/sh/y5by1sxjcvdl6rh/AAAIF4srkChW2FxeEfFG4d_ta?dl=0
> Thanks.
>
> [[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]]