Shapefile | Java library to read point and polygon shape files

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?

There are some other posts out there related to this one, such as these: Post 1, Post 2, Post 3. However, none of them deliver what I am hoping for. What I want is to be able to draw a line segment from a specific point (a sampling location) to the edge of a polygon fully surrounding that point (a lake border) in a specific direction ("due south" aka downward). I then want to measure the length of that line segment in between the sampling point and the polygon edge (really, it's only the distance I want, so if we can get the distance without drawing the line segment, so much the better!). Unfortunately, it doesn't seem like functionality to do this already exists within the sf package: See closed issue here.

I suspect, though, that this is possible through a modification of the solution offered here: See copy-pasted code below, modified by me. However, I am pretty lousy with the tools in sf--I got as far as making line segments that just go from the points themselves to the southern extent of the polygon, intersecting the polygon at some point:

I am working on submitting an R package to CRAN. Right now I am trying to reduce the memory footprint of the package. Because this package deals with spatial data that has a very particular format, I want to include a properly formatted shapefile as an example. If I include the full-size original shapefile, there are no warnings (other than file size) in the R CMD checks. However, if I crop the file and include the cropped version in the package (in "inst/extdata") I get this warning:

I am trying to use pysftp's getfo() to read a shapefile (without downloading it). However the output I get does not seem workable and I'm not sure if its possible do this with a shapefile.

Ideally I would like to read in the file and convert it to a Geopandas GeoDataFrame.

Is it possible to import a .tpk map file either in R or QGIS and use it as a shapefile? I need the coordinates of specific locations in the .tpk map that could be extracted with a shapefile map. I do not have access to ArcGIS.

Grateful for any guidance!

For additional information about .tpk, please see: https://pro.arcgis.com/en/pro-app/latest/help/sharing/overview/tile-package.htm

I'm using this dataset: https://sedac.ciesin.columbia.edu/data/set/gpw-v4-population-density-rev11/data-download (Gridded population density of the world)

With this map: https://data.humdata.org/dataset/uganda-administrative-boundaries-as-of-17-08-2018 (Uganda administrative boundaries shapefile)

I have clipped the uganda map to the region I need, like so:

I am trying to extract countries from NetCDF3 data using the pdsi monthly mean calibrate data from: https://psl.noaa.gov/data/gridded/data.pdsi.html. I am using the following code which performs a spatial merge of coordinates and identifies countries based on a shapefile of the world.

PDSI data format

In previous versions of the spsurvey package, it was possible to draw random points within polygons in a shapefile using a somewhat complicated design specification. (See here for an example).

The newly updated version of spsurvey (5.0.1) appears very user-friendly, except I cannot figure out how to perform a GRTS draw of more than one point within polygons. Below is an example:

Suppose I want to draw 10 random points using GRTS within the states of Montana and Wyoming. The grts() call requires an sf object, so we can get an sf object first.

I have a shapefile with 7 regions. I have an excel file with data about reptiles in these 7 regions. I merged this shapefile with excel.

Using ggplot I tried to generate facet_wrap() from nome_popular, however the rest of the polygon parts were omitted in each facet created.

My tentative code

shapefile: https://drive.google.com/file/d/1I1m9lBX69zjsdGBg2zfpii5H4VFYE1_0/view?usp=sharing

excel:https://docs.google.com/spreadsheets/d/1eKQWWCAalehTTrUuqUlMPQnSTEZxF--g/edit?usp=sharing&ouid=118442515534677263769&rtpof=true&sd=true