setwd("~/Programming/stunts/grcarread")
verx<-read.table(file="STANSX_r01_vex.dat")
poly<-read.table(file="STANSX_r01_pol.dat",fill=T)
ax1<-3
ax2<-2
plot(verx[,ax1],verx[,ax2],type="p",col="NA",asp=1)
for (i in 1:length(poly[,1])){
  polvec<-as.vector(poly[i,])
  polvec<-polvec[!is.na(polvec)]
#  print(polvec)
  polcoordx<-vector(length=length(polvec))
  for (j in 1:length(polvec)){
    polcoordx[j]<-verx[polvec[j],ax1]
    print(polcoordx)
  } 
  polcoordy<-vector(length=length(polvec))
  for (j in 1:length(polvec)){
    polcoordy[j]<-verx[polvec[j],ax2]
  } 
  polygon(polcoordx,polcoordy,border="red")
  text(mean(polcoordx),mean(polcoordy),format.hexmode(i),col="blue")
}
for (i in 1:length(verx[,1])){
  text(verx[i,ax1],verx[i,ax2],format.hexmode(i-1))
}
