blog/2023/07/28/gradient-map-fills-r-sf/index

[giscus[bot]]

blog/2023/07/28/gradient-map-fills-r-sf/index

Fix overplotted points on maps by creating bins or filled desntiy gradients using R, {ggplot2}, and {sf}

https://www.andrewheiss.com/blog/2023/07/28/gradient-map-fills-r-sf/index.html

[JoFrhwld]

Great stuff! Geospatial analysis is one of those things that’s gotten so much better in R since I last tried messing around with in about 10 years ago.

[LelViLamp]

Hey Andrew, many thanks for this tutorial! My density plots in option 3 include grid lines and I do not know how to get rid of them. Do you have an idea on how to remove them?

Output of sessionInfo()

R version 4.3.1 (2023-06-16)
Platform: x86_64-apple-darwin20 (64-bit)
Running under: macOS Monterey 12.6.8

Matrix products: default
BLAS:   /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib 
LAPACK: /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/lib/libRlapack.dylib;  LAPACK version 3.11.0

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

time zone: (removed)
tzcode source: internal

attached base packages:
[1] stats     graphics  grDevices datasets  utils     methods   base     

other attached packages:
 [1] stars_0.6-2            abind_1.4-5            patchwork_1.1.2        rnaturalearth_0.3.3   
 [5] tigris_2.0.3           spatstat_3.0-6         spatstat.linnet_3.1-1  spatstat.model_3.2-4  
 [9] rpart_4.1.19           spatstat.explore_3.2-1 nlme_3.1-162           spatstat.random_3.1-5 
[13] spatstat.geom_3.2-4    spatstat.data_3.0-1    sf_1.0-14              lubridate_1.9.2       
[17] forcats_1.0.0          stringr_1.5.0          dplyr_1.1.2            purrr_1.0.1           
[21] readr_2.1.4            tidyr_1.3.0            tibble_3.2.1           ggplot2_3.4.2         
[25] tidyverse_2.0.0        pacman_0.5.1          

loaded via a namespace (and not attached):
 [1] tidyselect_1.2.0      viridisLite_0.4.2     farver_2.1.1          fastmap_1.1.1         digest_0.6.33        
 [6] timechange_0.2.0      lifecycle_1.0.3       magrittr_2.0.3        compiler_4.3.1        rlang_1.1.1          
[11] tools_4.3.1           utf8_1.2.3            yaml_2.3.7            knitr_1.43            sp_2.0-0             
[16] classInt_0.4-9        KernSmooth_2.23-21    withr_2.5.0           grid_4.3.1            polyclip_1.10-4      
[21] fansi_1.0.4           e1071_1.7-13          colorspace_2.1-0      scales_1.2.1          spatstat.utils_3.0-3 
[26] cli_3.6.1             rmarkdown_2.23        generics_0.1.3        rstudioapi_0.15.0     httr_1.4.6           
[31] tzdb_0.4.0            DBI_1.1.3             proxy_0.4-27          splines_4.3.1         parallel_4.3.1       
[36] vctrs_0.6.3           Matrix_1.5-4.1        jsonlite_1.8.7        hms_1.1.3             tensor_1.5           
[41] units_0.8-2           goftest_1.2-3         glue_1.6.2            lwgeom_0.2-13         stringi_1.7.12       
[46] gtable_0.3.3          deldir_1.0-9          munsell_0.5.0         pillar_1.9.0          rappdirs_0.3.3       
[51] htmltools_0.5.5       R6_2.5.1              evaluate_0.21         lattice_0.21-8        renv_0.17.3          
[56] class_7.3-22          Rcpp_1.0.11           uuid_1.1-0            spatstat.sparse_3.0-2 mgcv_1.8-42          
[61] xfun_0.39             pkgconfig_2.0.3  

[joeroe]

Did you leave out color =NA from geom_sf(data = ga_churches_density, aes(fill = v), color = NA)?

[jethro-bodine]

Nice! Maybe congressional districts would be a good geographic unit, since they have, or are supposed to have, pretty much the same population counts.

[ariggs154]

Hey Andrew,
Just discovered your website, it’s amazing and certainly going to be a great tool in the future. I followed your steps and an great density map for my lat/long points. Very happy.

I did have a question; can I carry attributes through when changing to a ppp object? For example, I have annual precipitation attached to my lat/long and I’m attempting to get a density plot of annual precipitation (where its highest and lowest across the state). Any advice would be helpful, thanks in advance!