Zi 字媒體

版權聲明本文來源於易小蟹說統計本文由數據之王（ID：shujuzhiwang）編輯整理，版權歸原作者所有數據之王（ID：shujuzhiwang）由大中華大數據協會運營，如有邀約合作，首頁回複合作需求，我們將及時回復R語言里的ggmap包是繪製地圖的利器。在確保電腦可以上外網的情況下(為了保證進入google地圖的URL)，ggmap包幫助我們輕鬆下載google map，再加上和ggplot2製圖包的疊加，可以輕鬆、方便地繪製可視化地圖。 ggmap包中最常用的函數是ggmap，該函數用於從網上下載地圖。函數基本構成為： ggmap(get_googlemap(center=c(lon=-95.3632715,lat=29.7632836),zoom=3,maptype='terrain')) 其中，get_googlemap表示從網上獲取google地圖，center表示地圖中心，zoom表示地圖縮放程度，maptype表示地圖種類，可選參數為：』terrain』(地形圖)、』satellite』(衛星圖)、』roadmap』(道路圖)、』hybrid』(混合圖)。配合ggplot包的geom_point描點函數，我們可以在地圖上繪製數據點集。介紹好了基本函數、小蟹君就帶你們一起玩起來啦！~小蟹君從網上獲取了一個世界各地的地震數據集，記錄了地震發生的年份、經緯度、震級等指標。利用這部分數據，可以做出不同年份地震發生的動態地圖。想要繪製動態地圖，還需要在電腦上安裝ImageMagick軟體，該軟體的功能是將數張靜態圖組合成為1張gif動態圖。在R語言中，animation包可以完成ImageMagick軟體與R語言的交互，只需數行代碼便可實現動態圖的繪製。在載入了ggmap包、ggplot包、animation包之後，小蟹君嘗試繪製了美洲地震分布圖：▲美洲地震分佈帶結論：如圖所示的動態地圖展示了70年代美洲部分地區地震發生地點的變化，點的大小反映地震發生的等級大小。不難發現，地震頻發地區為西部臨海地區，內陸地區基本不在地震帶上。另外海里也時常有地震發生。現在，讓我們再把目光轉移到上。▲地震分佈帶(1990-2016)結論：上面這張動態地圖展示了1990年至2016年地震發生的區域。東西部差別很大，西部地震頻發，而東部幾乎沒有地震發生。值得注意的是，台灣成為了。最頻發地震的地區，這也印證了處於各大板塊相互作用的邊界地區是地震帶的主要分布區域。我們截取了一張2008年的地震分布圖，可以看到有一串明顯的紅點分佈在四川汶川附近2、獲取指定地點的經緯度並繪製於地圖ggmap包的強大之處在於它不僅可以獲取諸如地形圖、衛星圖、道路圖等多種形式的地圖，它還支持給定一個地點名稱，計算其經緯度。使用geocode函數，括弧內的參數可以是任意地點名稱，函數為我們返回此地點的經緯度。利用這個信息，我們可以輕鬆地在地圖上標註出我們想要標註的地點。小蟹君以上海各大學的分布圖為例給大家展示一下這個功能~首先獲取上海20所大學的英文名稱，並把它保存在一個數據框中。可以利用循環函數一次性返回20所大學所處的經緯度，再把它們同時標註在上海道路地圖上就好了！下圖是小蟹君得到的結果：▲上海20所高校分布圖地圖中的每個紅點代表一個大學。我們發現，有6所大學分佈在上海郊區的位置，這些大學的同學們生活還真是不太方便呢~此外，位於市區里的學校似乎呈現「扎堆」分佈的現象，我們把地圖放大了看(通過設置zoom參數可以實現地圖的縮放)：▲上海市高校分布圖哇，原來上海高校是集群分佈的，東北方向一群，西南方向一群，同學們之間交流感情很方便呢~3、計算兩地間距離最後，小蟹君給大家介紹ggmap包中的另一個有趣的函數，mapdist。函數的結構如下：mapdist(from, to, mode = c("driving", "walking", "bicycling"),output = c("simple", "all"))這個函數可以計算兩地之間的距離，並且在給定交通方式」mode」下，可以計算從一地出發前往另外一地的所需時間。比如我們想計算復旦大學到上海交通大學的距離和路途時間，則用：mapdist(「Fudan University」,」Shanghai Jiaotong University」,mode=」driving」)，得到結果：於是我們知道從復旦大學到上海交通大學，路程是43公里，開車的話需要48分鐘左右哦。數據之王（ID：shujuzhiwang）——貢獻數據智慧，釋放全球數據價值公眾號運營：大中華大數據協會識別二維碼，關注大數據行業最前瞻性內容