search

mdtanker / RIS_basement_sediment

Ross Ice Shelf basement depths and sediment thickness
MIT License
9 stars 2 forks source link

some wrong #1

Open Bohr-niels opened 3 hours ago

Bohr-niels commented 3 hours ago

Hello, @mdtanker, there are some errors occur when I run your code, I think the error is caused by the add_map_features() function, the code of the function is following:

微信图片_20241117210008

And the error message is: image One purpose of the function is to add geographical graticules on the figure which uses the linear projection, as we can see, to achieve this purpose, you use the basemap module of pygmt, and the projection is set to Polar Stereographic projection (EPSG3031), but the region you set is the region which used in linear projection, and the error message seems to be due to this. Is the error caused by the different pygmt version? And if we can't add geographical graticules by this method, how can I add them to the linear projection figure? Thanks a lot.😁😁😁

mdtanker commented 2 hours ago

Thanks for pointing this out! Since I made this repository I have switched all my plotting to use a Python package I created called PolarToolkit. It's basically a wrapper of PyGMT. Check out this function for adding lat long lines to a figure: https://polartoolkit.readthedocs.io/en/latest/_modules/polartoolkit/maps.html#add_gridlines

Let me know if that helps!