2013-11-01 33 views
-1

概念ggmap似乎很清楚,我的地圖:獲得與get_map包含位置與2+美國各州(市等)(位置矢量)

  1. 使用get_map以在獲得地圖某個空間縮放處的某個位置。

  2. 使用ggmap() + ggplot()來組合地圖圖像和ggplot圖形。

目前最大的挑戰在於在步驟1中,準確地說,location參數。除了經度/緯度對get_map接受一個字符串,但沒有向量。似乎沒有簡單的方法可以獲得包含兩個或更多州或兩個或更多城市的地圖。因此,要有一個繪製任意數量(可能相鄰)的美國狀態的函數,沒有捷徑,只能通過對每個狀態進行地理編碼並計算最佳位置和縮放(也不知道如何)的詳細過程。

+0

我們很多人都認爲這是你的責任,以提供加載相關的包,並隨後訪問中使用的shape文件,並最終建立了一個測試代碼案件與測試。 –

+0

我修改了這個問題,以避免將其轉向政策/公約/協議討論。雖然我欣賞你的意見。 – topchef

回答

0

你可以在這裏找到 - http://journal.r-project.org/archive/2013-1/kahle-wickham.pdf - 有辦法將shapefile放入ggmap中。

所有你需要做的是:

# read data into R (for example - your states) 
shapefile <- readShapeSpatial(’tr48_d00.shp’,proj4string = CRS("+proj=longlat +datum=WGS84")) 
# convert to a data.frame using fortify function 
data <- fortify(shapefile) 
#and plot your data using qmap or ggmap 
qmap(’texas’, zoom = 6, maptype = ’satellite’) + geom_polygon(aes(x = long, y = lat, group = group), data = data,colour = ’white’, fill = ’black’, alpha = .4, size = .3) 
+0

所以我認爲具體shapefile必須包括那些狀態(或城市,或者我處理的任何多個地理構件),對嗎? – topchef

+0

是的,你是對的;)當然,你也需要稍微改變qmap代碼 - 從'texas'到你的區域,改變縮放等。 –

+0

我想我需要修改我的問題,因爲整個觀點是找到參數化工件列表的方法(可能將其限制在美國境內的州或城市)。 – topchef