0

这是我的GitHub 存储库

我使用带有 GRASS 7.0.3 的 QGIS Desktop 2.14.0 将数据添加到 USA shapefile 以创建它，但是当我运行以下代码时

setwd("C:/GIS/US_Governor")

library(ggmap)
library(rgdal)
library(rgeos)
library(maptools)
library(dplyr)
library(tidyr)
library(tmap)



gov <- readOGR(dsn = ".", layer = "Age_of_US_Governor_Mar_2016")
qtm(gov, fill = "Age_of_G_1", fill.title = "Age of US Governor", fill.style="fixed",
    fill.palette = "Blues",
    fill.breaks=c(40,50,60,70)) + 
  tm_layout(legend.position = c("right", "bottom"), bg.color="lightgreen")

我得到一张带有大图例的小美国地图....如何增加地图的大小？谢谢！

在此处输入图像描述

更新

我添加了+ bb(xlim=150)，我得到了错误，Error: could not find function "bb"

我包括 library tmap，因此排除故障....

setwd("C:/GIS/US_Governor")

library(ggmap)
library(rgdal)
library(rgeos)
library(maptools)
library(dplyr)
library(tidyr)
library(tmap)



gov <- readOGR(dsn = ".", layer = "Age_of_US_Governor_Mar_2016")
qtm(gov, fill = "Age_of_G_1", fill.title = "Age of US Governor", fill.style="fixed",
    fill.palette = "Blues",
    fill.breaks=c(40,50,60,70)) + 
  tm_layout(legend.position = c("right", "bottom"), bg.color="lightgreen") +
  bb(xlim=150)
4

2 回答 2

3

@hrbrmstr 非常漂亮的地图！这是我的 tmap 解决方案：

devtools::install_github("mtennekes/tmap/pkg") 
# in the current CRAN version, there was a little problem with mapping
# numeric variables with one feature (in this example, the insets)
library(tmap)
library(grid)

# read shape data from https://github.com/Sonihal/US_Governor_Age
gov <- read_shape("Age_of_US_Governor_Mar_2016.shp")

# download and read simplified shape
f <- tempfile()
download.file("http://www2.census.gov/geo/tiger/GENZ2014/shp/cb_2014_us_state_20m.zip", destfile = f)
unzip(f, exdir = ".")
us <- read_shape("cb_2014_us_state_20m.shp")

# append gov data to this shape
us <- append_data(us, gov@data, key.shp = "NAME", key.data = "NAME_1")

# split US in three: contiguous, Alaska and Hawaii
us_cont <- us[!us$NAME %in% c("Alaska", "Hawaii", "Puerto Rico"), ]
us_AL <- us[us$NAME=="Alaska", ]
us_HI <- us[us$NAME=="Hawaii", ]

# plot contiguous states
tm_shape(us_cont, projection=2163) +
  tm_polygons("Age_of_G_1", title="Age of US Governor", style="fixed",
              palette="Blues",
              breaks=c(40,50,60,70, 80)) +
  tm_layout(frame=FALSE,
            legend.position = c("right", "bottom"), bg.color="lightgreen",
            inner.margins = c(.25,.02,.02,.02))

# create inset map of Alaska
map_AL <- tm_shape(us_AL, projection = 3338) +
  tm_polygons("Age_of_G_1", title="Age of US Governor", style="fixed",
              palette="Blues",
              breaks=c(40,50,60,70, 80),
              legend.show=FALSE) +
tm_layout(frame=FALSE, bg.color="lightgreen")

# create inset map of Hawaii
map_HI <- tm_shape(us_HI, projection = 3759) +
  tm_polygons("Age_of_G_1", title="Age of US Governor", style="fixed",
              palette="Blues",
              breaks=c(40,50,60,70, 80),
              legend.show=FALSE) +
  tm_layout(frame=FALSE, bg.color="lightgreen")

# plot insets
library(grid)
print(map_AL, vp=viewport(x=.15, y=.15, width=.3, height=.3))
print(map_HI, vp=viewport(x=.4, y=.1, width=.2, height=.1))

在此处输入图像描述

代码基于https://github.com/mtennekes/tmap#2-us-choropleth

于 2016-03-30T09:12:27.177 回答
3

如果您愿意使用一些肘部油脂和 ggplot2，以下内容可以帮助您：

  • 更快的绘图（您的 QGIS 多边形对于广泛的专题地图来说实在是太详细了）
  • 重新定位 AK 和 HI 的美国地图投影
  • 更好的休息（真的有4组）
  • 更好的颜色（viridis FTW）

我也尝试对代码进行一些注释。

library(sp)
library(rgeos)
library(rgdal)
library(maptools)
library(ggplot2)
library(ggalt)
library(ggthemes)
library(viridis)

# Read in your map w/data -------------------------------------------------

guvnah <- readOGR("Age_of_US_Governor_Mar_2016.shp",
                  "Age_of_US_Governor_Mar_2016",
                  stringsAsFactors=FALSE)
guvnah_dat <- guvnah@data

# Get a map with simpler polygons -----------------------------------------

URL <- "http://eric.clst.org/wupl/Stuff/gz_2010_us_040_00_20m.json"
fil <- basename(URL)
if (!file.exists(fil)) download.file(URL, fil)

# Handy function ----------------------------------------------------------

make_uniform <- function(SPDF) {
  # bivand & fowler
  pref <- substitute(SPDF)
  newSPDF <- spChFIDs(SPDF,
                      as.character(paste(pref, rownames(as(SPDF,"data.frame")),
                                         sep="_")))
  return(newSPDF)
}

# We'll need these projections later --------------------------------------

albers <- "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs"
latlon <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

# Read in the simpler map -------------------------------------------------

us <- readOGR(fil, "OGRGeoJSON", stringsAsFactors=FALSE)

# I have code that can move AK & HI but the coord are in meters 
# which is fine for albers but not for the generic longlat
# so we switch the lighter map to albers (briefly) and swap it
# back to longlat to prep for use of albers coord_proj in the
# event one wanted to add points or lines down the road.

us_aea <- spTransform(us, CRS(albers))
us_aea <- make_uniform(us_aea)

alaska <- us_aea[us_aea$NAME=="Alaska",]
alaska <- elide(alaska, rotate=-50)
alaska <- elide(alaska, scale=max(apply(bbox(alaska), 1, diff)) / 2.3)
alaska <- elide(alaska, shift=c(-2100000, -2500000))
proj4string(alaska) <- proj4string(us_aea)

hawaii <- us_aea[us_aea$NAME=="Hawaii",]
hawaii <- elide(hawaii, rotate=-35)
hawaii <- elide(hawaii, shift=c(5400000, -1400000))
proj4string(hawaii) <- proj4string(us_aea)

us_aea <- subset(us_aea, !(NAME %in% c("Alaska", "Hawaii", "Puerto Rico")))
us_aea <- make_uniform(us_aea)
us_aea <- rbind(us_aea, alaska, hawaii)

# Back to latlon & make it usable in ggplot2 ------------------------------

us <- spTransform(us_aea, CRS(latlon))
us_map <- fortify(us, region="NAME")

# Get your age data into decent breaks ------------------------------------

guvnah_dat$age_brks <- cut(guvnah_dat$Age_of_G_1, 
                           breaks=c(40, 50, 60, 70, 80), 
                           labels=c("40-50", "50-60", "60-70", "70-80"),
                           include.lowest=TRUE)

# Do the thing ------------------------------------------------------------

gg <- ggplot()
gg <- gg + geom_map(data=us_map, map=us_map,
                    aes(x=long, y=lat, map_id=id),
                    color="white", size=0.1, fill=NA)
gg <- gg + geom_map(data=guvnah_dat, map=us_map,
                    aes(fill=age_brks, map_id=NAME_1),
                    color="white", size=0.1)
gg <- gg + scale_fill_viridis(name="Age of U.S. Governors", discrete=TRUE)
gg <- gg + coord_proj(albers)
gg <- gg + theme_map()
gg <- gg + theme(legend.position="bottom")
gg

在此处输入图像描述

于 2016-03-28T11:15:58.213 回答