# Local geobox in Albers Equal Area projection in GeoMakie

**URL:** <https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975>\
**Category:** Visualization\
**Tags:** geomakie\
**Created:** [May 30, 2024, 4:19pm UTC](https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975 "2024-05-30T16:19:13Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![yum-yab](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/yum-yab/32/208984_2.png) [@yum-yab](https://discourse.julialang.org/u/yum-yab)\
**Post date:** [May 30, 2024, 4:19pm UTC](https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975/1 "2024-05-30T16:19:14Z")

</div>

Hi everyone,

i want to do something very simple: I want to cut out a local map of europe and the atlantic ocaen, but have troubles doing so in GeoMakie and the Documentation is lacking any good examples that do that.

Passing the limits to GeoAxis does not work, my current status is:

```julia
using GeoMakie
using GLMakie

lonlims = (-90, 40)
latlims = extrema(20, 80)

ax = GeoAxis(fig[1, 1], limits=(lonlims, latlims), dest="+proj=aea +lat_1=20 +lat_2=80 +lon_1=-90 +lon_2=40 +lon_0=-25 +lat_0=50")

lines!(ax, GeoMakie.coastlines(); color=:black, transformation=(; translation=(0, 0, 1000)))

```

I only seem to get a cutout by using `proj=merc` and passing the limits, why is it not working? Or is this not possible with GeoMakie (altough I remember that it worked a year back with orthographic projection).

Cheers and thanks for the help.

---

<div class="post-metadata">

**Author:** ![consumere](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/consumere/32/209031_2.png) [@consumere](https://discourse.julialang.org/u/consumere)\
**Post date:** [June 19, 2025, 12:08pm UTC](https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975/2 "2025-06-19T12:08:27Z")

</div>

I’m not exactly sure what you want, but i guess this would be the way. Note that you define the projection within the `GeoAxis` function.

```julia
fig = Figure(fontsize=12)
lonlims = (-90, 40)
latlims = (20, 80)
ax = GeoAxis(fig[1, 1], limits=(lonlims, latlims), 
    dest="+proj=aea +lat_1=20 +lat_2=80 +lon_1=-90 +lon_2=40 +lon_0=-25 +lat_0=50")
    #dest="proj=merc +lat_ts=20 +lon_0=-100 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs")
lines!(ax, GeoMakie.coastlines()) 
fig

```