Skip to content

Commit

Permalink
Fix errors in the examples
Browse files Browse the repository at this point in the history
  • Loading branch information
asinghvi17 committed Jun 8, 2024
1 parent 5a1f60b commit a8483c2
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 0 deletions.
2 changes: 2 additions & 0 deletions examples/cartopy/arrows.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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)

Expand Down
3 changes: 3 additions & 0 deletions examples/cartopy/vesta.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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);
Expand Down

0 comments on commit a8483c2

Please sign in to comment.