demo/moebius2.R

In Rgnuplot: R Interface for Gnuplot

# world map on a moebius strip
Gprun("bind \"p\" \"n=strftime('moebiusB%d.%b.%Y-%H.%M.%.3S.png',time(0.0));set terminal pngcairo;set output n;replot;set terminal wxt\"\nset angles degrees\nset parametric\nset isosamples 50\nf(u,v) = cos(u) + v*cos(u/2.0) * cos(u)\ng(u,v) = sin(u)+v*cos(u/2.0)*sin(u)\nh(u,v) = .5*v*sin(u/2.0)\nunset key;unset tics;unset border;unset colorbox\nset palette model RGB file \"earth_day.pal\" u ($1/255):($2/255):($3/255)\nset pm3d corners2color c1\n#set view equal xyz\nset hidden3d\nk=200.0\nsplot \"earth_day_moebiusXYcoords.dat\" u (f(-$1,$2/k)):(g(-$1,$2/k)):(h(-$1,$2/k)):3  w pm3d notit\n", 
    TRUE)