36

R のコンター プロットのページとページ (stackoverflow に関する多くのヒントを含む) を成功せずに調べました。ルワンダの地図を追加するなど、等高線を作成するデータは次のとおりです (データは、経度、緯度、雨量の x、y、z の 14 の値で構成されています)。

Lon Lat Rain
28.92   -2.47   83.4
29.02   -2.68   144
29.25   -1.67   134.7
29.42   -2.07   174.9
29.55   -1.58   151.5
29.57   -2.48   224.1
29.6    -1.5    254.3
29.72   -2.18   173.9
30.03   -1.95   154.8
30.05   -1.6    152.2
30.13   -1.97   126.2
30.33   -1.3    98.5
30.45   -1.81   145.5
30.5    -2.15   151.3

これが私がstackoverflowから試したコードです:

datr <- read.table("Apr0130precip.txt",header=TRUE,sep=",")
x <- datr$x
y <- datr$y
z <- datr$z

require(akima)

fld <- interp(x,y,z)

par(mar=c(5,5,1,1))
filled.contour(fld)

補間が失敗します。助けていただければ幸いです。

4

1 に答える 1