#Data Electron Microprobe from Nov3, 2005 (1103)
quartz(wi=9, he=9)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
filnam<-"K1-2b"
xy<-read.table(paste(filnam,"_bse.size",sep=''), skip=2)
attach(xy)
xi<-V1; yi<-V2
detach(xy)
xy<-read.table(paste(filnam,"_bse.pgm",sep=''), skip=4)
attach(xy)
XX<-V1;
detach(xy)
XX<-matrix(XX,xi,yi)
XX[1025:1080,1:790]<-XX[1025:1080,790:1]
XX[1:1024,790:770]<-XX[1:1024,770:790]
image(XX, col=gray(1:255/255),xlab='', ylab='')

quartz(wi=9, he=9)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_2_Ca.txt",sep=''), skip=33)
attach(xy)
YY<-t(as.matrix(xy))
# print(names(xy))
image(YY, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=1.5, he=1.8)
lx<-126; ux<-175; ly<-146; uy<-205
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(YY[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
yy<-YY[lx:ux,ly:uy]
detach(xy)



quartz(wi=9, he=9)
npar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_4_P.txt",sep=''), skip=33)
attach(xy)
ZZ<-t(as.matrix(xy))
# print(names(xy))
image(ZZ, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=1.5, he=1.8)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(ZZ[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
zz<-ZZ[lx:ux,ly:uy]
detach(xy)

quartz(wi=9, he=9)
npar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
plot(zz[1:3000], yy[1:3000])

quartz(wi=9, he=9)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_1_Na.txt",sep=''), skip=33)
attach(xy)
WW<-t(as.matrix(xy))
# print(names(xy))
image(WW, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=2, he=2)
lx<-126; ux<-200; ly<-150; uy<-225
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(WW[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
ww<-WW[lx:ux,ly:uy]
detach(xy)

quartz(wi=9, he=9)
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
xy<-read.table(paste(filnam,"_3_K.txt",sep=''), skip=33)
attach(xy)
UU<-t(as.matrix(xy))
# print(names(xy))
image(UU, col=gray(1:255/255),xlab='', ylab='')
quartz(wi=2, he=2)
lx<-126; ux<-200; ly<-150; uy<-225
opar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
image(UU[lx:ux,ly:uy], col=gray(1:255/255),xlab='', ylab='')
uu<-UU[lx:ux,ly:uy]
detach(xy)

quartz(wi=9, he=9)
npar<-par(mar=c(0,0,0,0)+0.1, ann=FALSE)
plot(uu[1:5625], ww[1:5625])

quartz(wi=9, he=9)
dens<-matrix(0,80,80)
for (i in 1:5625) {dens[uu[i]+1,ww[i]+1]<- dens[uu[i]+1,ww[i]+1] +1}
contour(dens^0.25)


par(opar)