#!/usr/bin/Rscript
### Visualize the topographical csv files output by terranostra.lisp
### Daniel Vedder, 12/9/2018

world = "world_large"
w = read.csv(paste0(world, ".csv"), header=F)

jpeg(paste0(world, "_topo.jpg"))
image(as.matrix(w), col=terrain.colors(range(w)[2]-range(w)[1]))
dev.off()

cols=c()
palette = terrain.colors(range(w)[2]-range(w)[1])
for (y in 1:length(w[1,])) {
    for (x in 2:length(w[,1])) {
        cols = c(cols, palette[w[x,y]-range(w)[1]])
    }
}

l = 720*(length(w)/100)
jpeg(paste0(world, "_iso.jpg"), width=l, height=l)
persp(1:(length(w[,1])), 1:(length(w[1,])), as.matrix(w), axes=F,
      theta=110, phi=50, expand=0.4, col=cols, box=F, shade=0.5, border=NA)
dev.off()