Hide area outside of area of interest in ggspatial
I made a map using ggplot2, ggspatial and sf. I used ggspatial so that I could download open source base maps without relying on Google and needing an API key. The ggspatial annotation_map_tile call downloads a set of tiles that provide the base map. I am trying to white out the map outside of the state I am interested in so that the base map only shows within the state.
Here's my code. The dput for the sf object is at the bottom. When you run it, you see a map of New Mexico with a black line outlining the state. I'd like to make the map white outside of the state.
library(ggplot)
library(ggspatial)
library(sf)
ggplot() +
annotation_map_tile(zoom = 7, type = "hikebike") +
geom_sf(data = sf, fill = NA, color = "black", lwd = 2)
I tried two things. First, I though about somehow inverting the fill argument in geom_sf, but don't think there's a way to do that. My other attempt was to include all states in an extra geom_sf call and fill those white. That might have worked but it required downloading over 4,000 map tiles, which I'd like to avoid, even though map tiles can be saved using the cachedir argument.
sf <- structure(list(GEOID = "35", geometry = structure(list(structure(list(
list(structure(c(-109.050044, -109.050173, -109.049843, -109.049813,
-109.049112, -109.049195, -109.048763, -109.049106, -109.048769,
-109.04859, -109.048599, -109.048731, -109.048296, -109.048286,
-109.047612, -109.047653, -109.047653, -109.047645, -109.047638,
-109.047117, -109.04748, -109.047453, -109.046905, -109.047013,
-109.047117, -109.047116, -109.047237, -109.04747, -109.046827,
-109.046909, -109.047045, -109.04687, -109.046564, -109.047298,
-109.046662, -109.047145, -109.046426, -109.047006, -109.046182,
-109.046182, -109.046156, -109.046086, -109.045363, -109.046104,
-109.045624, -109.046072, -109.045851, -109.046084, -109.046796,
-109.046481, -109.046509, -109.046296, -109.046295, -109.046024,
-109.046055, -109.046054, -109.046011, -109.045973, -109.045729,
-109.046183, -109.045431, -109.045433, -109.045407, -109.045272,
-109.045244, -109.045223, -108.958868, -108.954404, -108.620309,
-108.619689, -108.320721, -108.320464, -108.2884, -108.288086,
-108.250635, -108.249358, -108.000623, -107.420913, -106.877292,
-106.869796, -106.750591, -106.675626, -106.661344, -106.628733,
-106.628652, -106.617125, -106.617159, -106.500589, -106.343139,
-106.293279, -106.248675, -106.247705, -106.201469, -106.006634,
-105.997472, -105.996159, -105.716471, -105.66472, -105.62747,
-105.533922, -105.512485, -105.508836, -105.465182, -105.447255,
-105.442459, -105.41931, -105.251296, -105.220613, -105.1208,
-105.029228, -105.000554, -104.73212, -104.732031, -104.645029,
-104.625545, -104.624556, -104.519257, -104.338833, -104.250536,
-104.007855, -103.734364, -103.733247, -103.155922, -103.002199,
-103.002247, -103.001964, -103.002198, -103.002518, -103.002252,
-103.002188, -103.002565, -103.002434, -103.041924, -103.041745,
-103.041674, -103.040824, -103.041305, -103.042186, -103.041716,
-103.041917, -103.041146, -103.041554, -103.042366, -103.042775,
-103.042497, -103.042377, -103.042366, -103.042339, -103.042395,
-103.042568, -103.042711, -103.0426, -103.04252, -103.043261,
-103.042642, -103.042521, -103.042781, -103.04277, -103.042769,
-103.042827, -103.043286, -103.043072, -103.043594, -103.043589,
-103.043588, -103.043582, -103.043538, -103.043583, -103.043611,
-103.043585, -103.043613, -103.043613, -103.043614, -103.04363,
-103.043693, -103.043919, -103.043944, -103.043946, -103.043979,
-103.043936, -103.043719, -103.043644, -103.043569, -103.043516,
-103.043686, -103.043744, -103.043767, -103.043721, -103.043771,
-103.043746, -103.043555, -103.043531, -103.043617, -103.04395,
-103.044893, -103.045698, -103.045644, -103.046907, -103.047346,
-103.049096, -103.049608, -103.050148, -103.050532, -103.051087,
-103.051535, -103.051363, -103.051664, -103.05261, -103.056655,
-103.057487, -103.057856, -103.059242, -103.05972, -103.060103,
-103.063905, -103.06398, -103.064452, -103.064625, -103.064679,
-103.064657, -103.064569, -103.064701, -103.064862, -103.064807,
-103.064916, -103.064889, -103.064672, -103.064699, -103.064711,
-103.064698, -103.064807, -103.064827, -103.064799, -103.064798,
-103.064864, -103.064633, -103.064815, -103.064761, -103.064788,
-103.064761, -103.064696, -103.064422, -103.064348, -103.064344,
-103.064423, -103.085876, -103.088698, -103.215641, -103.267633,
-103.267708, -103.270383, -103.278521, -103.326501, -103.748317,
-103.875476, -103.980179, -104.024521, -104.531756, -104.531937,
-104.640918, -104.643526, -104.918272, -105.077046, -105.078605,
-105.11804, -105.131377, -105.132916, -105.14824, -105.15031,
-105.153994, -105.390396, -105.427049, -105.428582, -105.429281,
-105.731362, -105.750527, -105.854061, -105.886159, -105.9006,
-105.998003, -106.125534, -106.18184, -106.200699, -106.205915,
-106.313307, -106.376861, -106.394298, -106.411075, -106.565142,
-106.566056, -106.587972, -106.595333, -106.598639, -106.599096,
-106.618486, -106.619448, -106.623568, -106.631182, -106.636492,
-106.639529, -106.638186, -106.630114, -106.626466, -106.623216,
-106.621873, -106.619569, -106.618745, -106.619371, -106.620454,
-106.624299, -106.625535, -106.625123, -106.622819, -106.617708,
-106.614702, -106.616136, -106.623659, -106.622377, -106.622117,
-106.622529, -106.625322, -106.629747, -106.628663, -106.623933,
-106.611846, -106.614346, -106.623445, -106.625947, -106.633668,
-106.64084, -106.645479, -106.645646, -106.645296, -106.6429,
-106.638154, -106.633927, -106.630692, -106.629197, -106.630799,
-106.634873, -106.63588, -106.635926, -106.627808, -106.625763,
-106.621857, -106.614637, -106.605845, -106.605245, -106.602045,
-106.601945, -106.605267, -106.602727, -106.593826, -106.589045,
-106.588045, -106.582144, -106.581344, -106.577244, -106.570944,
-106.566844, -106.563444, -106.562945, -106.558444, -106.547144,
-106.545344, -106.544714, -106.542144, -106.542097, -106.535843,
-106.535343, -106.535154, -106.534743, -106.533043, -106.533,
-106.53248, -106.530515, -106.527943, -106.527738, -106.527623,
-106.527997, -106.528543, -106.528543, -106.750547, -106.750547,
-106.993544, -106.998235, -107.00056, -107.00056, -107.296824,
-107.422246, -107.422495, -108.208394, -108.208087, -108.208521,
-108.208572, -108.208573, -108.707657, -108.788711, -108.851105,
-108.861028, -109.050044, 31.332502, 31.480004, 31.499515,
31.499528, 31.636598, 31.796551, 31.810776, 31.843715, 31.861383,
31.870791, 32.013651, 32.028174, 32.084093, 32.089114, 32.426377,
32.681379, 32.686327, 32.689988, 32.693439, 32.77757, 33.06842,
33.069427, 33.091931, 33.092917, 33.137559, 33.137995, 33.208965,
33.250063, 33.365272, 33.36557, 33.36928, 33.372654, 33.37506,
33.409783, 33.625055, 33.74001, 33.875052, 34.00005, 34.522393,
34.522553, 34.579291, 34.771016, 34.785406, 34.799981, 34.814226,
34.828566, 34.959718, 35.250025, 35.363606, 35.546326, 35.54644,
35.614251, 35.616517, 35.8798, 35.888721, 35.92586, 35.925896,
36.002338, 36.117028, 36.181751, 36.500001, 36.874589, 36.874998,
36.968871, 36.969489, 36.999084, 36.998913, 36.998906, 36.999287,
36.999249, 36.99951, 36.999499, 36.99952, 36.999555, 36.999561,
36.999015, 37.000001, 37.000005, 37.000139, 36.992426, 36.992461,
36.993123, 36.993243, 36.993161, 36.993175, 36.993004, 36.992967,
36.993768, 36.99423, 36.99389, 36.994288, 36.994266, 36.994122,
36.995343, 36.995417, 36.995418, 36.995849, 36.995874, 36.995679,
36.995875, 36.995777, 36.995895, 36.995991, 36.996017, 36.995994,
36.995856, 36.995605, 36.995169, 36.995428, 36.992729, 36.993264,
36.993484, 36.993447, 36.993378, 36.993599, 36.994377, 36.993766,
36.993535, 36.994644, 36.996239, 36.998041, 36.998016, 37.000232,
37.000104, 36.911587, 36.909573, 36.719427, 36.675186, 36.61718,
36.602716, 36.526588, 36.500397, 36.500439, 36.318267, 36.317534,
36.055231, 35.837694, 35.825217, 35.814072, 35.796441, 35.791583,
35.622487, 35.250056, 35.241237, 35.211862, 35.183149, 35.182786,
35.181922, 35.178573, 35.159318, 35.144735, 35.142766, 35.135596,
35.125058, 35.109913, 34.899546, 34.850243, 34.792224, 34.747361,
34.671188, 34.653099, 34.619782, 34.46266, 34.459774, 34.459662,
34.455657, 34.405463, 34.400678, 34.397105, 34.393716, 34.390442,
34.388679, 34.384969, 34.38469, 34.383578, 34.380916, 34.37966,
34.379555, 34.312764, 34.302585, 34.289441, 34.256903, 34.087947,
34.079382, 34.063078, 34.049986, 34.043545, 34.04232, 34.041538,
34.037294, 34.032714, 34.018014, 34.003633, 33.974629, 33.945617,
33.906299, 33.901537, 33.8503, 33.824675, 33.74627, 33.737766,
33.701971, 33.672408, 33.658186, 33.650487, 33.64195, 33.629489,
33.570599, 33.388438, 33.329477, 33.315234, 33.260371, 33.256262,
33.219225, 33.042055, 33.038693, 33.01029, 32.999899, 32.964373,
32.959097, 32.900014, 32.879355, 32.868346, 32.857696, 32.85726,
32.849359, 32.82847, 32.827531, 32.784593, 32.783602, 32.777303,
32.726628, 32.708694, 32.690761, 32.682647, 32.64642, 32.624537,
32.601863, 32.600397, 32.587983, 32.522193, 32.145006, 32.123041,
32.087051, 32.000518, 32.000465, 32.000453, 32.000513, 32.000475,
32.000324, 32.000326, 32.000419, 32.00037, 32.000198, 32.000554,
32.000125, 32.00001, 32.000117, 32.000311, 32.000396, 32.000443,
32.000496, 32.000579, 32.000533, 32.000485, 32.000524, 32.000518,
32.000485, 32.000497, 32.000497, 32.000607, 32.000638, 32.0006,
32.000577, 32.001564, 32.002206, 32.00235, 32.00197, 32.0021,
32.002328, 32.002533, 32.00205, 32.001785, 32.001762, 32.001512,
32.001172, 32.001484, 32.001334, 32.000736, 32.000759, 32.000749,
32.000778, 32.000754, 32.000731, 32.000495, 31.994733, 31.990999,
31.989809, 31.985719, 31.980348, 31.97682, 31.971258, 31.97069,
31.97291, 31.972933, 31.971578, 31.966955, 31.964777, 31.963403,
31.961054, 31.957476, 31.954531, 31.952891, 31.956008, 31.956,
31.948439, 31.94551, 31.940863, 31.936621, 31.934863, 31.930053,
31.92657, 31.923614, 31.925335, 31.920003, 31.918003, 31.914034,
31.912227, 31.90979, 31.904598, 31.89867, 31.895649, 31.894859,
31.892933, 31.891663, 31.889184, 31.886411, 31.883717, 31.879697,
31.874478, 31.871514, 31.866235, 31.860593, 31.856276, 31.852854,
31.84649, 31.846305, 31.845905, 31.844405, 31.839605, 31.827912,
31.825024, 31.824901, 31.822706, 31.822106, 31.815506, 31.813906,
31.810406, 31.810206, 31.813306, 31.812606, 31.811104, 31.810406,
31.807305, 31.805007, 31.804287, 31.802107, 31.802146, 31.798607,
31.797507, 31.797089, 31.796107, 31.791907, 31.791829, 31.791914,
31.792103, 31.790507, 31.789761, 31.789119, 31.786945, 31.784407,
31.783907, 31.783706, 31.783898, 31.783689, 31.783671, 31.783679,
31.783513, 31.783762, 31.783599, 31.783599, 31.783599, 31.613489,
31.499798, 31.499742, 31.333395, 31.333191, 31.332365, 31.332301,
31.332315, 31.332502), .Dim = c(401L, 2L)))), class = c("XY",
"MULTIPOLYGON", "sfg"))), class = c("sfc_MULTIPOLYGON", "sfc"
), precision = 0, bbox = structure(c(xmin = -109.050173, ymin = 31.332301,
xmax = -103.001964, ymax = 37.000232), class = "bbox"), crs = structure(list(
input = "EPSG:4326", wkt = "GEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n CS[ellipsoidal,2],\n AXIS[\"geodetic latitude (Lat)\",north,\n ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n AXIS[\"geodetic longitude (Lon)\",east,\n ORDER[2],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n USAGE[\n SCOPE[\"Horizontal component of 3D system.\"],\n AREA[\"World.\"],\n BBOX[-90,-180,90,180]],\n ID[\"EPSG\",4326]]"), class = "crs"), n_empty = 0L)), row.names = c(NA,
-1L), sf_column = "geometry", agr = structure(c(GEOID = NA_integer_), .Label = c("constant",
"aggregate", "identity"), class = "factor"), class = c("sf",
"tbl_df", "tbl", "data.frame"))
Solution
You can create a big white mask around NM state, and overlay it over the annotation. It will clip the outside area.
You can get the NM state outline via {tigris} - highly recommended for any kind US vector data - and create the mask via a combination of sf::st_buffer() and sf::st_difference().
As the mask reaches, by necessity, outside of the plot area you will need to set the x and y limits of your plot by including a coord_sf() call.
Consider this code:
nm <- tigris::states() %>%
filter(STUSPS == "NM")
big_white <- st_buffer(nm, units::as_units(150, "km")) %>%
st_difference(nm)
ggplot() +
annotation_map_tile(zoom = 7, type = "hikebike") +
geom_sf(data = sf, fill = NA, color = "black", lwd = 2) +
geom_sf(data = big_white, fill = "white", color = NA) +
coord_sf(xlim = c(-102.5, -109), ylim = c(31.5, 37), crs = 4326)
- Retrieving expected data.frame for testthat expectation
- How to use an anonymous function with pipe operator?
- R passing RSelenium driver environment as function argument
- Keep file names merging list with lapply?
- replace strings that contains one another
- Stop parsing out zeros after decimals in ggplot2's annotate
- Plot data.tree coloring and labelling by level
- Extracting the date and time out of a date and time format
- How to reshape a dataframe into wide format with specified column pairs
- Show partial segments removed after limit scales for `geom_pointrange`
- Control ggplot2 legend look without affecting the plot
- Concatenate strings over all rows in a single column of a data frame
- Is there some way to keep variable names from.SD+.SDcols together with non .SD variable names in data.table?
- How to create mutually dependent checkbox columns in a DT table?
- A problem with an stacked, grouped and faceted graph y gg plot
- Installing R on Linux: configure: error: libcurl >= 7.28.0 library and headers are required with support for https
- Create a column conditional on multiple other character and numerical columns
- data.frame to array (2 columns)
- Changing options() in a function environment without changing options() in global environment in R?
- pivot(), group_by() and summarise() with nested data
- Why legend is not showing up in plot layout using patchwork in R?
- in dplyr::mutate, dplyr::starts_with works for .before but not .after?
- Finding optimal cut-off points for dividing a variable into equal groups considering ties in which quantiles do not work
- How to add a legend to hline?
- R Shiny Suppress slider handle for sliderInput until click
- brmultinom brglm2 vifs and dredge dont work
- R, how to get current time miliseconds, without extra information ? (Sys.time(), format("%X"))
- Installing multiple versions of R
- Command to see 'R' path that RStudio is using
- Fast pairwise longest common substring from start