Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Plotting error with moll projection for combined Atlantic and Southern Ocean #53

Closed
pwolfram opened this issue Sep 26, 2016 · 6 comments
Assignees
Labels

Comments

@pwolfram
Copy link
Contributor

Using the file at https://gist.github.com/43b8c23a27a84da35eadd9304a1b3753, there is a plotting error for the moll projection (which does not show up for moll2):

screenshot 2016-09-26 16 18 41

Script to reproduce error from repository root:

mkdir test_feature_collection
cd test_feature_collection/
ln -s ../ocean/region/North_Atlantic_Ocean/region.geojson NA_region.geojson
ln -s ../ocean/region/South_Atlantic_Ocean/region.geojson SA_region.geojson
ln -s ../ocean/region/Southern_Ocean/region.geojson SO_region.geojson
cd ..
./merge_features.py -d test_feature_collection/ -o AO_SO.geojson
./plot_features.py -f AO_SO.geojson 

Note, found in the process of testing #52. I believe this is a plotting error but it could be an error with preliminary changes in #52.

@pwolfram
Copy link
Contributor Author

Going further:

./difference_features.py -f AO_SO.geojson -m landice/region/Antarctica_IMBIE1/region.geojson -o AO_SO_masked.geojson
./plot_features.py -f AO_SO_masked.geojson 

moll is even more visibly messed up (geojson file at https://gist.github.com/79c3cc90ef72ad8320e799c4a54eb203)...

screenshot 2016-09-26 16 39 28

@xylar
Copy link
Collaborator

xylar commented Sep 29, 2016

@pwolfram, I can reproduce these bugs and will look into them when I have a chance.

@xylar
Copy link
Collaborator

xylar commented Sep 29, 2016

@pwolfram, I strongly suspect this is a bug in how shapes including the south (or presumably also the north) pole are handled in cartopy for this projection. It should be pretty easy to set up an example where we plot a box similar to the SouthernOcean in a simplified script to indicate the problem. Note that there are glitches in the Antarctic land coverage in both moll and moll2 projections as well.
Anyway, I don't think this is an issue for this repo. Please close if you agree and open an issue with cartopy instead.

@xylar
Copy link
Collaborator

xylar commented Sep 29, 2016

@pwolfram, here is a simpler script demonstrating that the problem is with cartopy, not with our plotting routine:
https://gist.github.com/xylar/b5e1835f031d122d98196ba52cacf027

@pwolfram
Copy link
Contributor Author

@xylar, I agree. I'll submit an issue with them, will pass on your gist, and reference this issue.

@pwolfram
Copy link
Contributor Author

The issue is SciTools/cartopy#802, I'm closing this PR for now.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants