Nope. Only what custom axes allow
To finish this (for now). I can’t make it work on Linux until GDAL fixes the bug (opened an issue). I thought that I could cheat it by running twice saving on disk and reloading but no way. Not even that worked.
On Win we can now do (function names are not definitive)
julia> G, cl = worldrectangular("@earth_relief_30m_g", pm=30, latlim=(80,-70), coast=:low);
julia> grid = worldrectgrid(proj="+proj=vandg");
julia> grdimage(G, shade=true, plot=(data=cl,), cmap=:geo, B=:none)
julia> GMT.plotgrid!(G, grid);
(@Andreas your attempts at are perfectly doable with some changes in the functions that do this, and not only coastlines but grids too.)
2 Likes
For reference: the issue @Joaquim opened can be found over at GDAL under #7644.
As I unfortunately can’t experiment with your work – is it still possible to extend beyond the +/-180° longitude domain as in your earlier work above? To get a chart that spans 400° longitude at the Equator?
That example goes from [-540 540], so it should be a matter of calculating what are the projected coords at +/- 400 and trim the plot at those XX.
@KristofKoch, you saw Evan’s answer in the GDAL issue. Only way for now is to have on-the-edge GDAL+PROJ build and follow my recipe above to use a local GMT build.
I have now a minimal working code, though are still many issues that need future work. Namely, the default title offset is not enough to not overlap the custom annotations and those, easily clutter near the prime meridian.
G, cl = worldrectangular("@earth_relief_30m_p", proj="vandg", latlim=(-70,85), coast=true);
grid = worldrectgrid(G);
grdimage(G, shade=true, plot=(data=cl,), cmap=:geo, B=:none, title="Van der Grinten")
plotgrid!(G, grid, figname="vandg.png", show=true)
1 Like
G, cl = worldrectangular("@earth_relief_30m_p", proj="aitoff", coast=true, latlim=(-80,89));
grid = worldrectgrid(G, annot_x=[-180,-150,0,150,180]);
grdimage(G, shade=true, plot=(data=cl,), cmap=:geo, B=:none, title="Aitoff")
plotgrid!(G, grid, figname="aitoff.png", show=true)
1 Like
And there is something strange with Winkel Tripel. It takes much longer to compute and the file is > 3 times larger
G, cl = worldrectangular("@earth_relief_30m_p", proj="wintri", coast=true, latlim=(-80,89));
grid = worldrectgrid(G);
grdimage(G, shade=true, plot=(data=cl,), cmap=:geo, B=:none, title="Winkel Tripel")
plotgrid!(G, grid, figname="wintri.jpg", show=true)
2 Likes
That is very impressive @Joaquim - thank you! Now on to installing the bleeding edge GDAL+PROJ versions and then follow your recipe … that day job takes up way too much of my time …
@KristofKoch, how is this project going?
@Andreas its on the backburner at the moment. Not dead but unfortunately not high priority either.
Well, the problem is solved in GMT.jl
@Joaquim you are right - you kindly implemented all the required technology to make it happen. A big thank you for that.
Unfortunately I have not enough time at the moment to make good use of it. I simply need more time and as soon as I produced a proper world map (3-4m wide, your implementation of Leland Browns Texture Shading, OSM Borders & Coastlines, etc) I’ll post it here.
1 Like