Open evas-ssec opened 3 years ago
I tried tracking this down, but am coming up empty on a final solution. If you lower the Y extents about ~300-500k meters south then you don't get the same issue and you even get more longitude lines. With the original code I went into the cartopy Gridliner source code with pycharm's debugger and as far as I can tell there are indeed duplicate Label objects in the list. There are other Labels that are duplicated (they have the same text and are on the same side of the axes) but one of the duplicates is always set to invisible so it isn't a problem. I'm having trouble following the low-level logic that generates these Labels from the list of ticks. And to be clear, the list of ticks to generate is as expected (ex. [15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0]
) and does not include duplicates.
It has to happen somewhere in this logic which leads to the generation of "heads" and "tails" and then from there the generation of the labels:
Given this is the geostationary projection and includes points outside of the lon/lat space, my guess is cartopy is failing to check the bounds/overlaps/intersections between the graticule lines and the map bounding box.
Description
The left and right gridline labels for a geostationary projection are duplicated and misplaced.
Code to reproduce
Traceback
(No error message is generated, only the faulty image output.)
Full environment definition
python 3.6.10 ### Operating system MacOS 10.15.7 ### Cartopy version cartopy 0.18.0 ### conda list ``` # packages in environment at /opt/anaconda3/envs/aitf_ql: # # Name Version Build Channel _anaconda_depends 2020.02 py36_0 aitf-ql CSPP.ABI.Quicklooks.v0.7 dev_0