diff --git a/visualize.R b/visualize.R new file mode 100644 index 0000000..e0c5dfb --- /dev/null +++ b/visualize.R @@ -0,0 +1,23 @@ +#!/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]]) + } +} + +jpeg(paste0(world, "_iso.jpg"), width=720, height=720) +persp(1:100, 1:100, as.matrix(w), axes=F, theta=110, phi=50, + expand=0.4, col=cols, box=F, shade=0.3) +dev.off()