hist2db<-function(x,y,kugirix,kugiriy){
sx<-kugirix[1:(length(kugirix)-1)]
sy<-kugiriy[1:(length(kugiriy)-1)]
fx<-kugirix[2:(length(kugirix))]
fy<-kugiriy[2:(length(kugiriy))]
boxf<-function(x){
j1<-max(which(sx<x))
}
counts<-0
i1<-0
i2<-0
j1<-0
j2<-0
sur<-0
density<-0
k<-1
mat<-matrix(0,length(sy),length(sx))
matc<-matrix(0,length(sy),length(sx))
x2<-x;
y2<-y;
for(i in 1:length(x)){
ly<-max(which(x[i]>=sy))
lx<-max(which(y[i]>=sx ))
matc[ly,lx]<-matc[ly,lx]+1;
}
mat<-mat/length(matc)
list(data.frame(sx,sy),data.frame(fx,fy),matc,mat)
}
dy<-1.05^(0:484)
dx<-1.05^(0:181)
x<-a[,5]
y<-a[,2]
t<-hist2db(y,x,dy,dx)
image(sort(unique(t[[1]][,3])),sort(unique(t[[1]][,1])), log(t[[2]]) , log="xy" )