2010-11-19 5 views
0

私はこのRの座標を使って通常のデータセットをどのようにマッピングできますか?

緯度経度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

私は色のグラデーションとしてmxgYieldをマップしようとしています。私はRでこれをどうやって行うべきですか?

おかげ

+0

可能な複製:http://stackoverflow.com/questions/1896419/plotting-a-3d-surface-plot-contour-map-overlay-using-r –

答えて

1

は、データ量に少し弱いようだが、アキマのinterpの機能の組み合わせを使用して成功し、機能をプロットベースの輪郭:

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

は(上きれいに見えます私の画面が表示されますが、データがある領域の等高線図が表示されます)。

+0

ありがとうございました。しかし、私が望むのは、異なる色のグラデーションによって示される歩留まりを持つ地図です。等高線プロットではありません。 – didimichael

+0

今私はそれを行うためにMAPSパッケージのイメージを使用できると思います。ありがとうたくさんの – didimichael

+0

イメージ(基本関数)もinterpオブジェクトで動作しますが、私にはうれしくありません。 contour(ak.yld、col = 1:10)またはimage(ak.yld、col = terrain.colors(8))を使用して等高線を色付けすることもできます。 –

関連する問題