我已经浏览了R中的等高线图页面(包括stackoverflow上的许多提示)但没有成功.这是我的轮廓数据,包括添加卢旺达地图(数据包括经度,纬度和降雨的14个值为x,y和z):
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)
插值fails.help将不胜感激.
我正试图为我的最终项目绘制一张地图,我正试图在美国制作BLock的犯罪热图.
对于每个街区,我都有Lat,Lon和犯罪率的预测.它遵循以下结构:
Lat / Lon / Prediction
-76.0 / 40.0 / 125
-76.120 / 40.5 / 145
-75.98 / 41.001 / 95
等等.
有没有办法绘制热图,显示预测为填充？
我认为这就是geom_tiles所做的,但是geom不起作用(也许是因为这些点间距不均匀)
任何帮助都非常受欢迎.请!
编辑
这是我到目前为止所尝试的:
-geom_density2d:
ggplot(ny2,aes(x=GEO_CENTROID_LON,y=GEO_CENTROID_LON,fill=prediction))+geom_density2d()
给我错误:"单位错误(tic_pos.c,"mm"):'x'和'units'必须长度> 0"
-geom_tiles:
ggplot(ny2,aes(x=GEO_CENTROID_LON,y=GEO_CENTROID_LON,fill=prediction))+geom_tile()
生成具有适当比例的绘图,但不生成地图上显示的数据.
关于chloropeth,如果我发生整个美国的块级信息,它会工作,但我找不到这样的数据.
可以在此处找到数据的SUBSAMPLE