2010-11-19 15 views
0

Ben buDüzenli veri kümesini R? Koordinatlarını kullanarak nasıl eşleyebilirim?

Lat boylam 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

gibi bir veri kümesi

40.36297989 -123.9729004 -124,058403 1,443492322

38,74369812 -123,1184998 1,938840925

40,64371872 15,00058428 Bir renk degrade olarak mxgYield eşleştirmek için çalışıyorum. Bunu R'de nasıl yapmalıyım?

Teşekkür

+0

Olası kopya: http://stackoverflow.com/questions/1896419/çizme-a-3d-yüzey-arsa-ile-kontur-harita-overlay kullanarak-r –

cevap

1

verilerin hacmine biraz zayıf görünüyor ama Akima en interp fonksiyonu bir arada ve baz kontur komplo işlevini kullanarak başarılı olacaktır:

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

(on güzel görünüyor Ekranım ancak yine de veriyle bölge için bir kontur çizimi görmelisiniz.)

+0

Teşekkürler. Ama istediğim, farklı renklerde bir degrade tarafından gösterilen verim ile bir harita. Kontur planı değil. – didimichael

+0

Şimdi, MAPS paketindeki görüntüyü kullanabileceğimi düşünüyorum. Çok teşekkürler – didimichael

+0

görüntü (temel işlev), aynı zamanda, nesne ile çalışır, ama bana iyi görünmüyor. Kontur çizgilerini de renklendirebilir (ak.yld, col = 1: 10) veya kullanabilirsiniz: image (ak.yld, col = terrain.colors (8)) –

İlgili konular