library(akima)

n<-seq(0.05,1,by=0.05)
x<-rep(n,each=20)
y<-rep(n,20)

cnt<-read.table(file="grid_counts_7460.txt",header=TRUE,row.names=1)
z<-cnt$CNT

# NO SCALE
a<-interp(y,x,z)

nrz <- nrow(a$z)
ncz <- ncol(a$z)
mycolors <- colorRampPalette( c("aquamarine", "darkblue", "yellow", "darkred") ) 
nbcol <- nrz*ncz
# Add white as the first colour (to clearly see zeros)
color <- c("white",mycolors(nbcol-1))
# Compute the z-value at the facet centres
zafacet <- a$z[-1, -1] + a$z[-1, -ncz] + a$z[-nrz, -1] + a$z[-nrz, -ncz]
# Recode facet z-values into color indices
facetcola <- cut(zafacet, nbcol)

pdf("my_landscape.pdf")
persp(a, shade=0.3, theta=-70, phi=10, col=color[facetcola], ticktype="detailed", nticks=10, xlim=c(0.05,1), ylim=c(0.05,1), cex.axis=0.6, xlab="Single-copyness", ylab="Universality", zlab="Group counts")
dev.off()


# SCALE
z2<-z/(z/(z^0.5))
z2[z2=="NaN"]<-0

a<-interp(y,x,z2)

nrz <- nrow(a$z)
ncz <- ncol(a$z)
mycolors1 <- colorRampPalette( c("cyan", "blue", "darkblue","yellow", "red") ) 
mycolors2 <- colorRampPalette( c("green", "darkgreen") ) 
mycolors3 <- colorRampPalette( c("purple", "pink") ) 
nbcol <- nrz*ncz
color <- c("white",mycolors1(nbcol-((nbcol/4)*3)),mycolors2(nbcol/4),mycolors3((nbcol/4)+500))
zafacet <- a$z[-1, -1] + a$z[-1, -ncz] + a$z[-nrz, -1] + a$z[-nrz, -ncz]
facetcola <- cut(zafacet, nbcol)

pdf("my_landscape2.pdf")
persp(a, shade=0.5, theta=-50, phi=20, col=color[facetcola], ticktype="detailed", nticks=10, xlim=c(0.05,1), ylim=c(0.05,1), cex.axis=0.6, xlab="Single-copyness", ylab="Universality", zlab="Scaled group counts")
dev.off()