United States missing using geostationary satellite view projection only about rnaturalearth HOT 5 OPEN

library(tidyverse)
library(rnaturalearth)
library(sf)
#> Linking to GEOS 3.12.1, GDAL 3.8.3, PROJ 9.1.1; sf_use_s2() is TRUE
library(s2)

proj <- "+proj=geos +h=35785831.0 +lat_0=45 +lon_0=-90.476194 +sweep=y"

world <- ne_countries(
  continent = c("north america", "south america"),
  scale = "medium",
  returnclass = "sf"
)

pts <- world |>
  # st_bbox() |>
  st_sample(size = 2e3, method = "regular") |>
  st_transform(proj)

ch <- pts[!st_is_empty(pts), drop = FALSE] |>
  s2_convex_hull() |>
  st_as_sf() |>
  st_transform(st_crs(world))

world2 <- world |>
  st_crop(ch)
#> Warning: attribute variables are assumed to be spatially constant throughout
#> all geometries

worldGeos <- st_transform(world2, crs = proj)

ggplot(data = worldGeos) +
  geom_sf(color = "lightgray", fill = "lightgray") +
  ggtitle("Geostationary Satellite View")

library(ggplot2)
library(rnaturalearth)
library(sf)
#> Linking to GEOS 3.12.1, GDAL 3.8.0, PROJ 9.2.0; sf_use_s2() is TRUE

world <- ne_countries(
  continent = c("north america", "south america"),
  scale = "medium",
  returnclass = "sf"
)

world <- st_cast(world, "POLYGON")
#> Warning in st_cast.sf(world, "POLYGON"): repeating attributes for all
#> sub-geometries for which they may not be constant

worldGeos <- st_transform(
  world,
  crs = "+proj=geos +h=35785831.0 +lat_0=45 +lon_0=-90.476194 +sweep=y"
)
ggplot(data = worldGeos) +
  geom_sf(color = "lightgray", fill = "lightgray") +
  ggtitle("Geostationary Satellite View")