# 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)