data-raw/plottas.r

library(codeutils)


tas <- read.csv(file="C:/Users/had06a/Dropbox/abalone/tas_100k.csv",header=TRUE)

plotprep(width=7,height=8,newdev = FALSE)

plot(c(145,148),c(-40,-44),type="n",xlim=c(143.5,149),xaxs="i",
     ylim=c(-44,-39),yaxs="i",cex=1.0,font=7,panel.first=grid())

ans <- c(802,2340,3732,6292,6302,6363,6686,7092,8128,8358,8488,8501)
ans <- unique(tas$poly)
for (i in ans) {
  pick <- which(tas$poly == i)
  polygon(tas$Long[pick],tas$Lat[pick],col=0)
}
haddonm/abspatial documentation built on June 7, 2019, 9:54 a.m.