利用R下載中國shp格式地圖腳本
利用R下載中國shp格式地圖腳本:
首先,需要安裝和加載一些必要的R包:
復(fù)制代碼
install.packages("maptools")
install.packages("ggplot2")
install.packages("rgeos")
install.packages("RColorBrewer")
library(maptools)
library(ggplot2)
library(rgeos)
library(RColorBrewer)
接著,可以從網(wǎng)上下載中國地圖的shp文件,例如從(gadm.org)這個(gè)鏈接下載中國地圖的shp文件。下載完成后,將文件解壓到本地指定的文件夾中,例如"./china_shp"。
復(fù)制代碼
# 讀取中國地圖的shp文件
china_map <- readShapeSpatial("./china_shp/gadm36_CHN_1.shp")
接下來,可以使用ggplot2包來繪制地圖。首先,將地圖數(shù)據(jù)轉(zhuǎn)換成數(shù)據(jù)框形式。
復(fù)制代碼
# 將地圖數(shù)據(jù)轉(zhuǎn)換成數(shù)據(jù)框形式
china_map_df <- fortify(china_map)
然后,利用ggplot2的geom_polygon函數(shù)來繪制地圖。
復(fù)制代碼
# 繪制地圖
ggplot(china_map_df, aes(x=long, y=lat, group=group)) +
geom_polygon(fill="white", color="black") +
coord_equal() +
theme_void()
最后,可以將地圖導(dǎo)出為shp格式。
復(fù)制代碼
# 導(dǎo)出為shp格式
writeSpatialShape(china_map, "./china_shp/china_map.shp")
注意:導(dǎo)出shp文件需要安裝rgdal包。
本文使用 文章同步助手 同步