R語言ggmap空間可視化機(jī)動車交通事故地圖|附代碼數(shù)據(jù)
原文鏈接:http://tecdat.cn/?p=12350
最近我們被客戶要求撰寫關(guān)于空間可視化的研究報告,包括一些圖形和統(tǒng)計輸出。
在本文中,我使用ggmap可視化紐約市的交通事故
數(shù)據(jù)來自紐約市開放數(shù)據(jù)。我的數(shù)據(jù)范圍是2012年至2015年。該數(shù)據(jù)跟蹤車輛的類型,發(fā)生事故的街道的名稱以及事故的經(jīng)度和緯度坐標(biāo)。兩個坐標(biāo)都保存為單個字符變量,稱為“ LOCATION”。
在下面,我加載數(shù)據(jù),刪除所有沒有位置坐標(biāo)的事故,并解析LOCATION變量以獲取經(jīng)度和緯度坐標(biāo)。我還解析日期變量以創(chuàng)建年份變量,并使用該變量創(chuàng)建兩個數(shù)據(jù)集:一個數(shù)據(jù)集包含2013年的所有車輛事故,另一個數(shù)據(jù)集包含2014年的所有車輛事故。
d_clean=d[which(regexpr(',',d$LOCATION)!=-1),]?comm=regexpr(',',d_clean$LOCATION)d_clean$loc=as.character(d_clean$LOCATION)d_clean$lat=as.numeric(substr(d_clean$loc,2,comm-1))d_clean$long=as.numeric(substr(d_clean$loc,comm+1,nchar(d_clean$loc)-1))d_clean$year=substr(d_clean$DATE,7,10)?d_2013=d_clean[which(d_clean$year=='2013'),c('long','lat')]d_2014=d_clean[which(d_clean$year=='2014'),c('long','lat')]
?
?
接下來,我使用get_map()查詢Google Maps并獲取紐約市的地圖。我使用stat_density2d()向該地圖添加了一個二維密度層。我對2013年和2014年的數(shù)據(jù)都執(zhí)行此操作,并使用gridExtra的grid.arrange()并排放置地圖。
ny_plot=ggmap(get_map('New York, New York',zoom=12, maptype='terrain'))?grid.arrange(plot1, plot2,nrow=1,ncol=2)
?
?

?
點(diǎn)擊標(biāo)題查閱往期內(nèi)容

R語言在地圖上繪制散點(diǎn)餅圖可視化

左右滑動查看更多

01

02

03

04

接下來,我按行政區(qū)域繪制了2013年的事故密度。我編寫了一個函數(shù)boro(),該函數(shù)會刪除所有缺少街道名稱的觀測值,并基于自治市鎮(zhèn)(2013)子集。我將它們堆疊在一起,同時添加自治市鎮(zhèn)名稱,然后在向量col_vals中為每個自治市鎮(zhèn)分配一種顏色。我使用stat_density2d的group參數(shù)以不同的顏色分別繪制每個自治市鎮(zhèn)的密度層。
plot4=ny_plot+??stat_density2d(data=full, geom='polygon',bins = 10, aes(x=full$long,y=full$lat,fill = Borough, alpha=..level..))+??scale_fill_manual(values=col_vals)+??#guides(fill = guide_colorbar(barwidth = 1, barheight = 12)) +??scale_alpha(guide = FALSE)+??xlab(' ')+ylab(' ')+??ggtitle('NYC Vehicle Accident Density by Borough, 2013')plot4
?

?

點(diǎn)擊文末?“閱讀原文”
獲取全文完整代碼數(shù)據(jù)資料。
本文選自《R語言ggmap空間可視化機(jī)動車交通事故地圖》。
點(diǎn)擊標(biāo)題查閱往期內(nèi)容
【視頻】關(guān)聯(lián)規(guī)則模型、Apriori算法及R語言挖掘商店交易數(shù)據(jù)與交互可視化|數(shù)據(jù)分享
R語言空氣污染數(shù)據(jù)的地理空間可視化和分析:顆粒物2.5(PM2.5)和空氣質(zhì)量指數(shù)(AQI)
上海無印良品地理空間分布特征與選址策略可視化研究
R語言空間可視化:繪制英國脫歐投票地圖
R語言在地圖上繪制散點(diǎn)餅圖可視化?r語言空間可視化繪制道路交通安全事故地圖
在GIS中用ggmap地理空間數(shù)據(jù)分析
tableau的騎行路線地理數(shù)據(jù)可視化
R語言推特twitter轉(zhuǎn)發(fā)可視化分析
618電商大數(shù)據(jù)分析可視化報告
用RSHINY DASHBOARD可視化美國投票記錄
python主題LDA建模和t-SNE可視化
R語言高維數(shù)據(jù)的主成分pca、 t-SNE算法降維與可視化分析案例報告
R語言動態(tài)圖可視化:如何、創(chuàng)建具有精美動畫的圖
Tableau 數(shù)據(jù)可視化:探索性圖形分析新生兒死亡率數(shù)據(jù)
R語言動態(tài)可視化:制作歷史全球平均溫度的累積動態(tài)折線圖動畫gif視頻圖