ita <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/penisola.txt",header=FALSE,sep =",") # il grafico plot(ita,asp=1) lines(ita) abline(h=seq(37,47,1),v=seq(7,18,1),lty=3,col="grey") x <- ita$V1; y <- ita$V2 # l'area usando come unità di misura i "quadri" delimitati da # longitudini e latitudini n <- length(x); area <- (y[n]+y[1])*(x[n]-x[1]) for (i in 1:(n-1)) area <- area + (y[i]+y[i+1])*(x[i]-x[i+1]); area <- area/2 area # la sua trasformazione in km^2 tenendo conto che # 1 grado di latitudine è 111 km (valore arrotondato) # 1 grado di longitudine in Italia è mediamente circa 79 km area*111*79; round(area*111*79,-3) # l'area della penisola italiana è effettivamente 252000 km^2