rggplot2r-sfggspatial
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)
- Calculate transition probabilities
- How to subset R dataframe based on specific values in several columns?
- Change colon to period in figure caption in Quarto
- Applying different functions to different groups using 2 data frames
- Disable button for 5 seconds with different label and change after it
- How to group columns with an odd index, based on columns with an even index?
- Different plots of marginal effects with interaction term and offset when using different packages
- Animate a MWE GIF in R
- populating a matrix by rows
- How to relevel in ROC curve in R/plotROC?
- R - glmer different results on different machines (non-deterministic)
- Dynamically operate on column (inside :=)
- Random change in obj_addr() output when including the objects into a list and vectorizing over them
- Add to curly-curly argument
- h2o predict error: Test/Validation dataset has no columns in common with the training set
- Is is possible to convert a dataframe object to a tribble constructor?
- encoding issues in R
- How to Scrape NBA stats page using rvest
- Conversion of Decimal Minutes Seconds into Decimal Degrees when cells within the Columns Contains the Factor 'Missed' in R
- Difficulty using case_when() to add column that, conditionally, pastes value from another column
- Flip the matrix
- Creating test for date columns incorrectly throwing error
- Exists() and getSymbols() can't find symbol listed in yahoo
- spearman's rho plot
- Use OpenMP on M2 Mac with R and data.table
- How to find the connected instances from a minimum spanning trees model in R
- Hide a column in shiny datatable but keep it searchable
- How to extract data from an edited datatable in shiny?
- Render Data Table with double header
- filter dataTables in Shiny Dashboad based on selectInput Values