Comment puis-je mapper l'ensemble de données normal en utilisant les coordonnées dans R?

J'ai un ensemble de données comme celui-ciComment puis-je mapper l'ensemble de données normal en utilisant les coordonnées dans R?

Lat Lon mxgYield

38,96152878 -123,5578995 7,925971605

39,24150085 -123,6392975 19,6720449

39,52162933 -123,7214966 19,777747

39,80192947 -123,8044968 12,23218451

40.36297989 -123.9729004 15.00058428

40,64371872 -124,058403 1,443492322

38,74369812 -123,1184998 1,938840925

Je suis en train de la carte mxgYield comme un dégradé de couleurs. Comment dois-je faire cela dans R?

Merci

Source

2010-11-19 didimichael

duplication possible: http://stackoverflow.com/questions/1896419/tracer-un-3d-surface-plot-avec-contour-map-overlay-using-r –

semble un peu faible sur le volume de données, mais réussiront en utilisant une combinaison de la fonction interpolée de Akima et le contour de base fonction de traçage:

yield <- read.table(textConnection("Lat Lon mxgYield 
38.96152878 -123.5578995 7.925971605 
39.24150085 -123.6392975 19.6720449 
39.52162933 -123.7214966 19.777747 
39.80192947 -123.8044968 12.23218451 
40.36297989 -123.9729004 15.00058428 
40.64371872 -124.058403 1.443492322 
38.74369812 -123.1184998 1.938840925"), header=TRUE) 
require(akima) 
ak.yld <- with(yield, interp(x=Lat, y=Lon, z=mxgYield)) 
contour(ak.yld)

alt text

(Attend plus jolie sur mon écran mais vous devriez toujours voir un tracé de contour pour la région avec des données.)

Source

2010-11-19 19:17:19

Merci. Mais ce que je veux, c'est une carte avec un rendement montré par un dégradé de couleurs différentes. Pas le tracé de contour. – didimichael

Maintenant, je pense que je peux utiliser l'image dans le paquet MAPS pour le faire. Merci beaucoup – didimichael

image (la fonction de base) fonctionne également avec l'objet interp, mais ne me semble pas agréable. Peut également colorer les lignes de contour avec: contour (ak.yld, col = 1: 10) ou utiliser: image (ak.yld, col = terrain.colors (8)) –

Comment puis-je mapper l'ensemble de données normal en utilisant les coordonnées dans R?

Répondre

Questions connexes