library("maptools")<br />
library("akima")<br />
x=readShapePoly('d:/map/jj.shp');<br />
plot(x)<br />
data2 <- read.csv("d:/aa.csv", header=TRUE)<br />
data2.li<-interp(data2$x,data2$y,data2$z)</p>
<p>> data2<br />
x y z<br />
1 29.83092 116.6157 56<br />
2 29.66447 116.2779 157<br />
3 29.61316 115.7791 675<br />
4 29.62643 115.4468 67<br />
5 29.79394 116.1604 1000<br />
6 29.64101 116.0135 578<br />
7 29.35230 116.3277 2483<br />
8 29.26252 115.0096 567<br />
9 29.40059 115.6196 437<br />
10 29.35914 115.9488 786<br />
11 28.99758 114.4326 861<br />
12 29.14067 115.7291 525<br />
以上是我做的代码,但是我不知道怎么用filled.contour()完成在地图上画图的语法,请各位帮忙看下
</p>
