熱線電話:13121318867

登錄
首頁精彩閱讀使用R讀取shapefile并畫地圖的過程
使用R讀取shapefile并畫地圖的過程
2018-06-04
收藏

使用R讀取shapefile并畫地圖的過程

下面是一個使用R語言繪制地圖的例子。
    # 首先需要加載需要使用的兩個package:  
    library(maptools)  
    library(rgdal)  
    # 讀取數據并繪圖,下面讀取的是中國行政區域的shapefile到變量x中:  
    x=readShapePoly('bou2_4p.shp')  
    # 現在我們可以直接畫出地圖:  
    plot(x) 

但是上面的繪圖并未使用投影而是直接使用經緯度作為平面坐標進行繪圖,所以是不合適的。下面的代碼首先進行mercator投影然后繪圖。
[plain] view plain copy
    # 首先,我們需要指定x本身的投影信息,下面一行代碼指定x對應的投影為原始的地理坐標(經緯度):  
    proj4string(x) <- CRS("+proj=longlat +ellps=WGS84")  
    # 然后我們指定新的投影方式,并將x投影到其上:  
    projNew <- CRS("+proj=merc +lat_0=45n +lon_0=100e")  
    xProj <- spTransform(x, projNew)  
    # 現在我們可以畫出它了:  
    dev.new()  

plot(xProj) 

數據分析咨詢請掃描二維碼

若不方便掃碼,搜微信號:CDAshujufenxi

數據分析師資訊
更多

OK
客服在線
立即咨詢
日韩人妻系列无码专区视频,先锋高清无码,无码免费视欧非,国精产品一区一区三区无码
客服在線
立即咨詢