2013-12-09 66 views
0

這是我第一次在這裏,我真的很難過這個。使用ggplot在ggmaps上分層連續數據的shapefile(熱圖)

因此,我有灣區房價(這裏稱爲SAMPLE3.csv)和海灣地區郵政編碼的形狀文件的文件。我能夠匹配這兩個,並使用強化來融化郵政編碼,並且能夠在地圖上繪製數據,但無法正確調整配色方案以製作熱圖(例如,使用heat.colors) 。這裏是我的代碼:

houseprices <- read.csv("~/DigitalFlyer/GIS/bayarea/SAMPLE3.csv", header=T) 
bayarea <- USA[USA$city == 'bayarea',] 
sf = fortify(bayarea, region="ZCTA5CE10") 
sf$score <- engagement[match(sf$id, houseprices$zip, nomatch=NA), 2] 

我可以用

qmap('oakland, ca', zoom=10,maptype='hybrid')+ 
geom_polygon(aes(x=long, y=lat, group=group), data=sf, colour="white", fill=sf$score, alpha=.5, size=.3) 

得到一個不錯的谷歌地圖爲基礎層,我可以使用

qplot(long, lat, data=sd, group=group , fill= score, geom="polygon"), 

但隨後得到與QMAP一個不錯的熱圖無法在谷歌地圖上分層。

這裏是我所得到的每個: 1)http://i.stack.imgur.com/vJOHe.jpg 2)http://i.stack.imgur.com/bCeG4.jpg

感謝您的幫助提前!

回答

0

我不能對此進行測試,因爲我沒有你的數據幀sf但你的第一個代碼問題是fill=應該是geom_polygon()aes()內的地方,以保證填充值是根據score值的數據幀。

qmap('oakland, ca', zoom=10,maptype='hybrid') + 
    geom_polygon(aes(x=long, y=lat, group=group, fill=score), 
        data=sf, colour="white",alpha=.5, size=.3) 
+0

嗨Didzis。非常感謝你的迴應,它的效果非常好!有沒有辦法上傳數據幀,以便將來我可以向用戶提供所有信息? – user3082069

+0

閱讀[這個問題](http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)和關於可重現問題的答案 –