Landmasks for each basin

[geanders]

Landmasks for each basin

[geanders]

In this case, we would need a function to ID the correct basin and then use that to select correct landmask from basin-specific ones.

[geanders]

It looks like this is in Liz's PR (~ 14 MB for landmask file).

We also have a landmask for part of the NW Pacific basin (in the file named "landmask_0p2deg_15to55N_70to140E.csv"). This file is around 900 KB.

landmask_nwpacific <- read_csv("data-raw/landmask_0p2deg_15to55N_70to140E.csv", 
               col_names = c("lon", "lat", "land")) %>% 
  mutate(land = as.factor(land), 
         land = fct_recode(land, land = "1", water = "0"), 
         land = fct_relevel(land, "land"))
ggplot(data = landmask_nwpacific , aes(x = lon, y = lat, color = land)) + 
  geom_point()

[geanders]

Here's the new landmask from Liz's PR:

load("data/landmask.rda")
ggplot(data = landmask, aes(x = longitude, y = latitude, color = land)) + 
  geom_point() + 
  geom_vline(aes(xintercept = 360)) + 
  coord_map() + 
  scale_color_manual(values = c("green", "blue"))

Here it is zoomed in for North America:

ggplot(data = landmask, aes(x = longitude, y = latitude, color = land)) + 
  geom_point() + 
  geom_vline(aes(xintercept = 360)) + 
  coord_map() + 
  scale_color_manual(values = c("green", "blue")) +
  xlim(c(260, 300)) + 
  ylim(c(15, 50))

[geanders]

I don't think we have documentation on the source of the global landmask. Can try to get that from Liz.

[geanders]

We will need to update the documentation (in "R/data.R") if we use this landmask.

[geanders]

As a comparison, here's the North Atlantic landmask we currently have:

ggplot(data = landmask_nwpacific , aes(x = lon, y = lat, color = land)) + 
     geom_point()

[geanders]

For global landmask, maybe change so 0/360 is prime meridian and then western Europe in low range? (currently continues past 360 moving past Europe, up to 540 in Pacific)

[geanders]

Illustration of previous comment:

[geanders]

Here are the lat/lon conventions for IBTrACS: