rm(list=ls())
filroot<- "LobOct02"
require(pixmap)
quartz(width=10.3, height=8)
PP<-read.pnm(paste(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)
XX<-read.csv(paste(filroot,'.csv',sep=''))
attach(XX)
Yi<- max(Ypos[!is.na(Ypos)])-Ypos
lenxx<- length(id)
pid<- 1:lenxx
breaks<- pid[is.na(Xpos)]
foci<-1:pid[breaks[1]-1]
outl<-pid[breaks[1]+1]:lenxx

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

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

browser()
detach(XX)