#CanalMap_area2ab_URL.R
rm(list=ls())
url<-"http://www.bio.umass.edu/biology/kunkel/pub/lobster/Rjoe/R8Oct03/"
filroot<- "area2b"
csvroot<- "area2ab"
require(pixmap)
quartz(width=9.5, height=8)
PP<-read.pnm(paste(url,filroot,".ppm",sep=''))
Att<-attributes(PP)[[1]]
print(names(PP))
AspRat<-0.70
Asp<-(Att[1]/Att[2])/AspRat
plot(PP, asp=Asp, axes=TRUE, main=filroot, cex.main=2)
XX<-read.csv(paste(url,csvroot,'.csv',sep=''))
attach(XX)
#Yi<- max(Ypos[!is.na(Ypos)])-Ypos
Yi<- Att[1]-YposB
lenxx<- length(id)
pid<- 1:lenxx
breaks<- pid[is.na(XposB)]
foci<-1:pid[breaks[1]-1]
# outl<-pid[breaks[1]+1]:lenxx

points(XposB[foci], Yi[foci], typ='p', pch=1, col='red', cex=Ptype[foci])

#lines(Xpos[outl], Yi[outl], col='white', lwd=2)
text(XposB[foci], Yi[foci],paste(pid[foci]) , col='blue', pos=4)

browser()
detach(XX)