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]
lox<- 295; upx<-310; loy<- 130; upy<-290
CanCas<-(CanCa<-as.matrix(read.table("JK canal plus exo plus endo_1.txt", skip=32)))[lox:upx,loy:upy]
CanPs<-(CanP<-as.matrix(read.table("JK canal plus exo plus endo_2.txt", skip=32)))[lox:upx,loy:upy]
image(t(CanP), col = grey(seq(0,1,0.01)))
image(t(CanCa), col = grey(seq(0,1,0.01)))
image(t(CanCas), col = grey(seq(0,1,0.01)))
image(t(CanPs), col = grey(seq(0,1,0.01)))

plot(c(0,2000),c(0,4000), 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)