test/grilla.r

In profesorpaiche/polilla: Find grids that intersects polygons

# Loading libraries and data
library(ggplot2)
library(dplyr)
library(sf)
library(maps)
library(polilla)
data(shp)
data(gridlon)
data(gridlat)

# Reprojection information
grid_north_pole_latitude = 39.25
grid_north_pole_longitude = -162
lat_center = 90 - grid_north_pole_latitude # displacement from the north
lon_center = grid_north_pole_longitude + 180 # displacement from the east
crs_ortho = paste0("+proj=ortho +lat_0=", lat_center, " +lon_0=", lon_center)

# Testing sfCoords
lonlat_sf = sfCoords(lon = gridlon, lat = gridlat, crs = 4326) |>
    st_transform(crs = crs_ortho)

# Testing extended crop
lonlat_sf = extraCrop(lonlat_sf, shp)

# Testing gridMaker
grid_poly = gridMaker(lonlat_sf, crs = crs_ortho)

# Testing getWeights
intersection = getWeights(grid_poly, shp)
grid_poly = left_join(grid_poly, intersection, by = "id")

# test plot
g = ggplot() +
    geom_sf(data = grid_poly, mapping = aes(fill = weights)) +
    geom_sf(data = shp, fill = NA) +
    geom_sf(data = grid_poly |> st_centroid(), color = "royalblue") +
    geom_sf(data = lonlat_sf, color = "brown")