I started to code in Julia half a year ago and I am quite excited about using GeoMakie. I work on ice sheet dynamics so I mostly want to visualize the poles. Therefore I want to use a stereographic projection and as a first example I tried to adapt the keyword argument in a GeoMakie example of the documentation (Examples · GeoMakie.jl). The problem I face is that I obtain a blank plot:
lons = -180:180;
lats = 0:90; # No need for the whole Earth
field = [exp(cosd(l)) + 3(y/90) for y in lats, l in lons];
G = mat2grid(field, x=lons, y=lats, proj4="+proj=longlat");
imshow(G, proj="+proj=sterea +lat_0=90", coast=true)
Great to see it working on your sides @aramirezreyes and @JanJereczek and sorry for providing only a very general advice. Really happy that GeoMakie as I understand received a significant refresh recently. Hope the development efforts will continue.
Hey, was on vacation for a bit so just saw this. Thanks @aramirezreyes for the solution!
Yeah, in general this is an issue in GeoMakie - if the limits of the plot are greater than the intrinsic limits of the projection, the whole thing goes blank, since the transformed limits become infinity. It’s hard to define a general solution because not all projections have well-defined rectangular boundingboxes.
You will see that the surface is not clipped by the edge of the Axis. My current thought is to just add a white polygon on top to simulate clipping until we get a proper nonlinear clip implemented, but not totally sure how that should work still.