# La terra e la rotta  longitudine = latitudine
# Riduco i margini dell'immagine (lo "spazio bianco" intorno): cosa non essenziale
par( mai = c(0.2,0.2,0.2,0.2) )
# Una ripartizione dell'intervallo [0,2*pi]
t <- seq(0,2*pi,len=200)
z0 <- c(-1,1); u <- rep(z0[1],4)
z <- array(u,dim=c(2,2)); x <- c(-1,1); y <- c(-1,1)
# In x,y,z,z0 ho definito il box. Ora scelgo il p. di vista e definisco la figura
traccia <- function(x,y,z,z0,th,ph,di) {
 F <- persp(x,y,z,theta=th,phi=ph,scale=TRUE,zlim=z0,xlim=x,ylim=y,d=di,border="white",box=FALSE)
 # equatore
 lines(trans3d(cos(t),sin(t),0,pmat=F),col="blue",lwd=2)
 # merdiano di long. 0
 lines(trans3d(cos(t),0,sin(t),pmat=F),col="brown",lwd=2)
 # paralleli
 for (ang in seq(-90,90,len=10)*pi/180)
 lines(trans3d(cos(ang)*cos(t),cos(ang)*sin(t),sin(ang),pmat=F),col="red",lty=3)
 # merdiani
 for (ang in seq(0,170,len=10)*pi/180)
 lines(trans3d(cos(ang)*cos(t),sin(ang)*cos(t),sin(t),pmat=F),col="red",lty=3)
 # la rotta longitudine=latitudine
 ang <- seq(-90,90,len=200)*pi/180
 lines(trans3d(cos(ang)*cos(ang),cos(ang)*sin(ang),sin(ang),pmat=F),col="black",lwd=2)
}
# Traccio la figura vista con theta=80°,phi=20°,distanza=1
traccia(x,y,z,z0,80,20,1)
#

# vediamo da un altro punto di vista
traccia(x,y,z,z0,80,90,20)
# vediamo da un altro punto di vista
traccia(x,y,z,z0,90,0,20)