{"action":"created","issue":{"url":"https://api.github.com/repos/SciTools/cartopy/issues/525","labels_url":"https://api.github.com/repos/SciTools/cartopy/issues/525/labels{/name}","comments_url":"https://api.github.com/repos/SciTools/cartopy/issues/525/comments","events_url":"https://api.github.com/repos/SciTools/cartopy/issues/525/events","html_url":"https://github.com/SciTools/cartopy/issues/525","id":53055389,"number":525,"title":"Comparing basemap and cartopy","user":{"login":"guziy","id":900941,"avatar_url":"https://avatars.githubusercontent.com/u/900941?v=3","gravatar_id":"","url":"https://api.github.com/users/guziy","html_url":"https://github.com/guziy","followers_url":"https://api.github.com/users/guziy/followers","following_url":"https://api.github.com/users/guziy/following{/other_user}","gists_url":"https://api.github.com/users/guziy/gists{/gist_id}","starred_url":"https://api.github.com/users/guziy/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/guziy/subscriptions","organizations_url":"https://api.github.com/users/guziy/orgs","repos_url":"https://api.github.com/users/guziy/repos","events_url":"https://api.github.com/users/guziy/events{/privacy}","received_events_url":"https://api.github.com/users/guziy/received_events","type":"User","site_admin":false},"labels":[],"state":"closed","locked":false,"assignee":null,"milestone":null,"comments":1,"created_at":"2014-12-29T17:16:21Z","updated_at":"2015-01-06T08:59:16Z","closed_at":"2015-01-06T08:59:16Z","body":"Hi:\r\n\r\nI was checking out cartopy since basemap author told that this is future....\r\n\r\nCan someone tell me if I can do a similar plot with cartopy? \r\n![image](https://cloud.githubusercontent.com/assets/900941/5570609/6d9016c4-8f52-11e4-9c61-572bd307f1fb.png)\r\n\r\nBelow is my effort to do it with cartopy, which does not seem to work very well.\r\n![image](https://cloud.githubusercontent.com/assets/900941/5570601/4143b774-8f52-11e4-8f3e-fb13ff8980e3.png)\r\n\r\nThe code for the images can be found here:\r\nhttp://nbviewer.ipython.org/github/guziy/PyNotebooks/blob/master/cartopy_compare_with_basemap.ipynb\r\n\r\nThanks"},"comment":{"url":"https://api.github.com/repos/SciTools/cartopy/issues/comments/68840970","html_url":"https://github.com/SciTools/cartopy/issues/525#issuecomment-68840970","issue_url":"https://api.github.com/repos/SciTools/cartopy/issues/525","id":68840970,"user":{"login":"pelson","id":810663,"avatar_url":"https://avatars.githubusercontent.com/u/810663?v=3","gravatar_id":"","url":"https://api.github.com/users/pelson","html_url":"https://github.com/pelson","followers_url":"https://api.github.com/users/pelson/followers","following_url":"https://api.github.com/users/pelson/following{/other_user}","gists_url":"https://api.github.com/users/pelson/gists{/gist_id}","starred_url":"https://api.github.com/users/pelson/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/pelson/subscriptions","organizations_url":"https://api.github.com/users/pelson/orgs","repos_url":"https://api.github.com/users/pelson/repos","events_url":"https://api.github.com/users/pelson/events{/privacy}","received_events_url":"https://api.github.com/users/pelson/received_events","type":"User","site_admin":false},"created_at":"2015-01-06T08:59:16Z","updated_at":"2015-01-06T08:59:16Z","body":"You might be interested in the docstring (not yet in the latest release) at https://github.com/SciTools/cartopy/blob/master/lib/cartopy/crs.py#L1036.\r\n\r\nI think the CRS you want is:\r\n\r\n```\r\nrpole_crs = ccrs.RotatedPole(\r\n                             pole_longitude=bmap_params[\"lon_0\"] + 180, \r\n                             pole_latitude=bmap_params[\"o_lat_p\"],\r\n                             )\r\n```\r\n\r\nNotice that the o_lon_p is not being used here. The reason is that the RotatedPole in cartopy has been implemented to work well with meteorological type projections where a pole rotation is used to provide an approximately equidistant limited area domain. The meteorological definition of rotated pole (as documented by the WMO) specifies the latitude longitude rotation in a different order to the way proj4 has implemented it, therefore cartopy moved to ratify the problem, hence the confusion here (esp. if you are used to the proj4 definition). For more information see the thread entitled \"proj4 specifications for a rotated latlon grid on a sphere?\" [here](http://lists.maptools.org/pipermail/proj/2012-August/006391.html).\r\n\r\nHTH\r\n\r\n"}}