quartz(width = 10, height = 10)
layout(matrix(c(1,2,1,2,3,5,4,5), 4, 2, byrow = TRUE))
lo<- 50; up<-250 
ApCas<-(as.matrix(read.table("JK apatite calib_1.txt", skip=32)))[lo:up,lo:up]
ApPs<-(as.matrix(read.table("JK apatite calib_2.txt", skip=32)))[lo:up,lo:up]
#plot(c(0,2000),c(0,4000), typ='n', main = 'apatite (r), clam (gr)', cex=0.25)

ClCas<-(as.matrix(read.table("JK clam calib_1.txt", skip=32)))[lo:up,lo:up]
ClPs<-(as.matrix(read.table("JK clam calib_2.txt", skip=32)))[lo:up,lo:up]

C3Cas<-(as.matrix(read.table("JK Ca-P calib_1.txt", skip=32)))[lo:up,lo:up]
C3Ps<-(as.matrix(read.table("JK Ca-P calib_4.txt", skip=32)))[lo:up,lo:up]


loy<-31; upy<-80; lox<- 151; upx<-200
CanCas<-(CanCa<-as.matrix(read.table("JK calcite plus exo_1.txt", skip=32)))[loy:upy,lox:upx]
lpx<-length(t(CanCa)[,1])
lpy<-length(t(CanCa)[1,])
CanPs<-(CanP<-as.matrix(read.table("JK calcite plus exo_2.txt", skip=32)))[loy:upy,lox:upx]

image(t(CanP), col = grey(seq(0,1,0.01)))
lines(c(lox,lox)/lpx,c(loy,upy)/lpy, col='purple')
lines(c(upx,upx)/lpx,c(loy,upy)/lpy, col='purple')
lines(c(lox,upx)/lpx,c(loy,loy)/lpy, col='purple')
lines(c(lox,upx)/lpx,c(upy,upy)/lpy, col='purple')
image(t(CanCa), col = grey(seq(0,1,0.01)))
lines(c(lox,lox)/lpx,c(loy,upy)/lpy, col='purple')
lines(c(upx,upx)/lpx,c(loy,upy)/lpy, col='purple')
lines(c(lox,upx)/lpx,c(loy,loy)/lpy, col='purple')
lines(c(lox,upx)/lpx,c(upy,upy)/lpy, col='purple')
image(t(CanCas), col = grey(seq(0,1,0.01)))
image(t(CanPs), col = grey(seq(0,1,0.01)), main=paste("(lx,ux),(ly,uy) = (",lox,",",upx,"), (",loy,",",upy,")",sep=''))

plot(c(0,1600),c(0,3800), typ='n', main = 'apatite (r), clam (gr)', cex=0.25)
points(ApPs,ApCas, typ='p', col='red', cex=0.25)
points(ClPs,ClCas, typ='p', col='green', cex=0.25)
points(C3Ps,C3Cas, typ='p', col='blue', cex=0.25)
points(CanPs,CanCas, typ='p', col='purple', cex=0.25)