#Data Electron Microprobe from Nov3, 2005 (1103)
h1<-w1<-9
h2<-w2<-2

paramIn<- function(typ,xo,yo,xl,yl,th){
	     TH<- th*pi/100
        ROT<-matrix(c(cos(TH),sin(TH),-sin(TH),cos(TH)),2,2)
        XY<-matrix(0,4,2)
	     xl<-xl/2; XY[1:2,1]<-xo-xl; XY[3:4,1]<- xo+xl
	     yl<-yl/2; XY[c(1,4),2]<-yo-yl; XY[2:3,2]<-yo+yl
	     mx<-mean(XY[,1]); my<-mean(XY[,2])
	     MXY<- t(matrix(c(mx,my),2,4))
	     XY<-MXY + t(ROT%*%t(XY-MXY))
	     
	     yl<-XY[1:2,2]; xl<-XY[1:2,1]
        yr<-XY[3:4,2]; xr<-XY[3:4,1]
        yu<-XY[2:3,2]; xu<-XY[2:3,1]
        yd<-XY[c(4,1),2]; xd<-XY[c(4,1),1]
	     xleft<-lm(xl~yl)
	     xright<-lm(xr~yr)
	     yup<-lm(yu~xu)
	     ydown<-lm(yd~xd)
	     Out<-c(1, xleft[[1]], xright[[1]], ydown[[1]],yup[[1]])
	     Out	
	     }
TestwParam<- function(p,x,y) {
	 In<-FALSE
	 if(p[1]==1) {
	 	         if(p[2]+p[3]*y <= x) {
	 	         	if(p[4]+p[5]*y >= x) {
	 	         		if(p[6]+p[7]*x <= y){
	 	         			if(p[8]+p[9]*x >= y){
	 	         				In<-TRUE}
	 	         			} } }
	 	 	         }
	  In
	  }
	  
Lines<-read.table("Ascii.listarea", header=TRUE)
attach(Lines)
nlin<- length(Lines[,1])
filcol<-c("red", "green", "blue", "purple", "orange", "yellow", "brown")
for (i in 1:nlin) {
filnam<-	jpat[i]
Pt<-paramIn(1,Xo[i],Yo[i],Xl[i],Yl[i],theta[i])
	
# Calcium plots
quartz(wi=w1, he=h1)
opar<-par(mar=c(0,0,0,0)+0.05, ann=FALSE)
xy<-read.table(paste(filnam,"_2_Ca.txt",sep=''), skip=33)
attach(xy)
YY<-t(as.matrix(xy))
# print(names(xy))
ny<- length(YY[,1])
nx<- length(YY[1,])
image(YY, col=gray(1:255/255),xlab='', ylab='') 
   mtext(paste(filnam, "Ca++", "Xo=", Xo[i], "Yo=",Yo[i]), 
   line=-1.5, col="white" )
detach(xy)
for (j in 1:nx) {for (k in 1:ny) {
		TF<-TestwParam(Pt,j,k)
		if(TF==TRUE) {points(j/nx,k/ny, col=filcol[i], cex=0.1)}
		}
	}  # end for j


# Phosphorous plots
quartz(wi=w1, he=w1)
npar<-par(mar=c(0,0,0,0)+0.05, ann=FALSE)
xy<-read.table(paste(filnam,"_3_P.txt",sep=''), skip=33)
attach(xy)
ZZ<-t(as.matrix(xy))
# print(names(xy))
image(ZZ, col=gray(1:255/255),xlab='', ylab='') 
   mtext(paste(filnam, "PO4", "Xo=", Xo[i], "Yo=",Yo[i]), 
   line=-1.5, col="white" )
detach(xy)

for (j in 1:nx) {for (k in 1:ny) {
		TF<-TestwParam(Pt,j,k)
		if(TF==TRUE) {points(j/nx,k/ny, col=filcol[i], cex=0.1)}
		}
	}# end for j

} 

detach(Lines)
par(opar)