sf

问题 I would like to create a leaflet map where you can select multiple polygons and this will update the selectizeInput() in a shiny app. This would including removing a selected polygon, when it is removed in the selectizeInput() . I have slightly changed/updated the code from the answer here (use of sf instead of sp and more dplyr where I could work out what the base R was). The polygons could probably be updated with an observeEvent tied in with input$clicked_locations , but not sure exactly

问题 I was wondering if there is a way to create linestring from two points given in the same row in a dataframe in a new geometry column. In other words longitudes and latitudes of the two points are given in a dataframe like the following: df <- data.frame(id = c("a", "b"), lon1 = c(1,2), lat1 = c(3,4), lon2 = c(5,6), lat2 = c(7,8)) where lon1 and lat1 represent the coordinates of the first point and lon2 and lat2 are the coordinates of the second point. The desired dataframe would have two rows

问题 I'm trying to fill a combined categorical and continuous variable on a map. So, for instance, in my minimally reproducible example below, say I want to display the number of KrispyKreme Donut shops in each county, which is generally a continuous variable I want to fill on a gradient. But I also have counties that forbid KrispyKremes indicated by a "-1" and those that have them under construction "-2". I want to display these in a different color not mapped on the gradient. I also have NA in

问题 I would like to make pie graph with ggplot and sf package. I have very simple data and un think there is simple method to do. All my data are pourcent : data <- data.frame(Territory = c(1, 2, 3, 4, 5), Pins = c(25, 45, 45, 60, 75), oak = c(45, 50, 45, 20, 15), land = c(30, 5, 10, 20, 10), sum = c(100, 100, 100, 100, 100)) And my graph code : read_sf("territories.shp") %>% left_join(data, by = "Territory") %>% ggplot() + geom_sf(aes(fill = Pins), color = "black") + theme_bw() + xlab("") + ylab

问题 The required libraries are in non-standard locations and I can install rgdal by the command below: install.packages('rgdal', type = "source", configure.args=c('--with-gdal-config=/home/programs/anaconda3/bin/gdal-config', '--with-proj-include=/home/programs/anaconda3/include', '--with-proj-lib=/home/programs/anaconda3/lib')) But I can't install sf package. It keeps saying configure: error: proj_api.h not found in standard or given locations. . The installation code and output is below. Do you

问题 If you have spatial point and polygon time series data, how do you do a spatial join/merge and "normal" merge of a non-spatial variable together? Point data over years that I want to merge into yearly polygons and then summarise ( xvar ) by year: #spatial point data by year library(sf) set.seed(10) df_point <- data.frame(id = c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7), year = c(2016, 2017, 2018, 2019, 2016, 2017, 2018, 2019, 2016, 2017, 2016, 2017, 2016, 2017, 2016, 2017, 2016,

问题 If you have spatial point and polygon time series data, how do you do a spatial join/merge and "normal" merge of a non-spatial variable together? Point data over years that I want to merge into yearly polygons and then summarise ( xvar ) by year: #spatial point data by year library(sf) set.seed(10) df_point <- data.frame(id = c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7), year = c(2016, 2017, 2018, 2019, 2016, 2017, 2018, 2019, 2016, 2017, 2016, 2017, 2016, 2017, 2016, 2017, 2016,

问题 I am creating a map of Brazil with its state boundaries, which is straight forward to achieve using ggplot2 and geom_sf . However, this time around, instead of color filling each state with data, I want to fill each state's boundaries with an external image (png), similar to this example of the largest employer in each state. I have tried a few settings of geom_image and even geom_flag to no avail (hence why I am not posting any code here). Are there any suggestions on how to approach this

问题 Using the tutorial here: https://www.r-spatial.org/r/2018/10/25/ggplot2-sf-2.html (in the section called "States (polygon data)"), I'm trying to convert polygons to an sf object. library("maps") states <- st_as_sf(map("state", plot = FALSE, fill = TRUE)) head(states) I get the error message: Error in as_mapper(.f, ...) : argument ".f" is missing, with no default Thank you in advance! 回答1: Looks like it's confusing map() from the purrr package with map() from the maps package. Try: states <-

问题 Calculate the total number of sampling points within each grid cell of a spatial grid. I would like to make a grid and calculate the total count of sampling points within each grid cell. I created a randomly generated data and grid, and tried to calculate the number of records within a grid cells using both the sf and raster packages, using previous similar SO questions, but wthout success. I have also looked into the extract function. Im fairly new to spatial analysis. library(sf) library