diff --git a/examples/cartopy/arrows.jl b/examples/cartopy/arrows.jl
index 943b02e2..39e5e2dd 100644
--- a/examples/cartopy/arrows.jl
+++ b/examples/cartopy/arrows.jl
@@ -3,6 +3,8 @@
 
 # The big idea is to plot arrows from one CRS in another!
 
+using GeoMakie, CairoMakie
+
 xs = LinRange(311.9, 391.1, 30)
 ys = LinRange(-23.6, 24.8, 20)
 
diff --git a/examples/cartopy/vesta.jl b/examples/cartopy/vesta.jl
index a0c44189..baa66a67 100644
--- a/examples/cartopy/vesta.jl
+++ b/examples/cartopy/vesta.jl
@@ -22,6 +22,9 @@ vesta_image_matrix = FileIO.load(joinpath(DataDeps.datadep"vesta_image", "PIA170
 # We use Proj's [ellipsoid parameters](https://proj.org/en/9.3/usage/ellipsoids.html#ellipsoid-size-parameters)
 # to define the ellipsoid of Vesta, specifically the semi-major and semi-minor axes (`+a` and `+b` respectively).
 vesta_crs = GeoMakie.GeoFormatTypes.ProjString("+proj=longlat +a=285000 +b=229000 +type=crs")
+# Additionally, since the Vesta CRS is on a non-Earth datum, we have to set this environment variable
+# so that Proj knows that we are aware of this problem:
+ENV["PROJ_IGNORE_CELESTIAL_BODY"] = "yes"
 # Now, we can create a `Raster` from the image.
 vesta_raster = Raster(
     rotr90(vesta_image_matrix);