polygons | playable post on how harmless choices | Frontend Framework library

How can I add a mouse 'click' event to each polygon plotted on a leaflet map? I want to do this so I can then filter a separate widget based on data from the features (in this case WD21CD).

I am working on a project where I am using a shape file to make a choropleth map of the United States. To do this, I downloaded the standard shape file here from the US Census Bureau. After a little bit of cleaning up (there were some extraneous island territories which I removed by changing the plot's axis limits), I was able to get the contiguous states to fit neatly within the bounds of the matplotlib figure. For reference, please see Edit 4 below.

Edit 1: I am using the cb_2018_us_state_500k.zip [3.2 MB] shape file.

The only problem now is that by setting axis limits I now am no longer able to view Alaska and Hawaii (as these are obviously cut out by restricting the axis limits). I would now like to add both of these polygons back in my map but now towards the lower part of the plot figure (the treatment that is given by most other maps of this type) despite its geographical inaccuracy.

To put this more concretely, I am interested in selecting the polygon shapes representing Alaska and Hawaii and moving them to the lower left hand side of my figure. Is this something that would be possible?

I can create a Boolean mask using:

I have a problem with drawing the legends in a ggplot. I have already searched the net for hours for a solution but have not found anything yet.

I am trying to create a ggplot with different line and polygon shapefiles. I draw the shapefiles with the function geom_sf(). The drawing is not a problem. My problem is related to the legend creation.

The goal is to create a plot that has a single legend and has all the different shapefile types including their colors and or fill astetic. That means, a line-shape should be represented in the legend as a simple line and a polygon-shape as a simple polygon. I have managed to do all that.

The main difficulty is that one of the polygonshape (for better readability) not only uses the "fill" argument as astetic but also the "color" argument.

If I create a plot where all polygonshapes always use the fill argument as astetic and all lineshapes always use the color argument as astetic everything works fine. But as soon as I include a polygonshape in the plot which uses both arguments at the same time, it destroys the whole legend structure. In this case, for example, the symbols for the line shapefiles suddenly have a frame and/or background although they should not have one.

Is it not possible to build something like this with ggplot or am I just being clumsy?

Summary:

Analogous to the minimal example created below, I would like to create a single legend in this plot where:

• A = a polygon symbol is without border (as it can be seen in the
  plot)
• B = a polygon symbol is with border (as seen in the plot) and
• C = a simple line without border or background (as it can be seen in the plot).

How do I get this to work?

I have a question concerning rasterization of polygons by maximum overlap, i.e assign the value of the polygon that has the highst area overlap with the raster cell.

The real world exercise is to rasterize polygons of soil-IDs in R, in order to produce relatively low resolution maps of soil properties as model inputs.

The problem is that the rasterize() function of the terra package (and similar stars' st_rasterize()) assigns the cell value from the polygon that contains the cell midpoint. If a raster cell contains multiple polygons, I would rather like to select the value of the polygon (soil-ID), which has the highest aerea cover in a raster cell.

Here is a small self-contained example that visualizes my problem, using terra.

Here is a simple shiny app:

I have a large polygon file, small polygon file and points file. What I do here is loop through large polygons to find which small polygons intersect. Then calculate the area of each small polygon within the large one. And then I loop through the small polygons to find points statistics in each of them.

I have found number_of_somethin value in each small polygon. And the question would be how to can I sum all number_of_somethin small polygons values within the large polygon and store the results in original large_polygon file as a new column, let's say large_polygon['smth_sum']?

With df_res_2.loc[idx, 'smth'] = number_of_somethin I get number_of_somethin values in each small polygon inside the large ones. Now I need to sum them in large_polygon['smth_sum']

Note: FID is the id for large polygons and ID is the id for small polygons

EDIT : Adding code that generates the polygons. df is the base data loaded from csv and processed into clusters :

In my very simple case I would like to display the heatmap of the points in the points GeoJSON file but not on the geographic density (lat, long). In the points file each point has a confidence property (a value from 0 to 1), how to display the heatmap on this parameter? weight=points.confidence don't seem to work.

for exemple:

I have a shape file of towns in the north of Spain that I have to join into groups (municipalities or comarcas in Spanish). I've used st_union from the sf package to join them successfully (and each one is their own SpatialPolygonsDataFrame object with a single polygon). I plot each of the municipalities individually and they look fine.

However, once I want to combine the municipalities into a single SpatialPolygonsDataFrame object with multiple polygons, I can't for the life of me manage to do it. I've tried three approaches mostly based on this answer: https://gis.stackexchange.com/questions/155328/merging-multiple-spatialpolygondataframes-into-1-spdf-in-r and this one https://gis.stackexchange.com/questions/141469/how-to-convert-a-spatialpolygon-to-a-spatialpolygonsdataframe-and-add-a-column-t – If I use raster::union it throws out the error Error in .rowNamesDF<-(x, value = value) : invalid 'row.names' length – If I use a simple rbind it throws out the error Error in SpatialPolygonsDataFrame(pl, df, match.ID = FALSE) : Object length mismatch: pl has 7 Polygons objects, but df has 4 rows Or something similar for 6/11 of the municipalities. – If I try a lapply approach (more convoluted) it seems to work but one I plot it using leaflet the municipalities that gave the error when trying to raster::union or rbind don't look as they should/don't look as they do when I plot them individually.

** Municipalities 1 and 2 work fine. 3 and 4 for example do not. **

Here's a link to the two files needed to reproduce my code below: – Link to shape files: https://www.dropbox.com/sh/z9632hworbbchn5/AAAiyq3f_52azB4oFeU46D5Qa?dl=0 – Link to xls file that contains the mapping from towns to municipalities: https://www.dropbox.com/s/4w3fx6neo4t1l3d/listado-comarcas-gipuzkoa.xls?dl=0

And my code: