I have a hard time trying to describe the shapefile exactly, but I am trying to generate a grid over the Gulf of St. Lawrence, such that only the gridsquares which lie over water are included. I would like this to extend to the narrow southern strait at (-60, 47)
I'm really new to GIS in R, and below is my attempt at brute forcing it by generating rows of gridsquares and specifying the bounding box for each line individually. I used st_make_grid from the sf library to make the 5 rows. Each gridsquare is 0.25x0.25 degrees.
I feel like there must be a more elegant solution to this... if anyone knows one.