Blue marble doesn't map to regional ortho projection correctly

Hello,

I could easily be doing something wrong. I posted this issue about two
years ago and never saw a response, so I'll try again. It's more of an
academic question - I can work around the issue. In the following code, I'm
trying to plot a regional ortho projection over central Europe (I've
commented out the global projection), and then trying to use a bluemarble
background. As you can see from the attached image, global bluemarble
background is plotted, not a regional one. This also happens with
shadedrelief and etopo.

<http://matplotlib.1069221.n5.nabble.com/file/n46944/ortho.png>

import mpl_toolkits.basemap as bm
import matplotlib.pyplot as plt
import numpy as np

"""
theMap = bm.Basemap(projection='ortho',
                    lat_0=48.2, lon_0=16.4)
"""

theMap = bm.Basemap(projection='ortho',
                    lat_0=48.2, lon_0=16.4,
                    llcrnry=-1000000.0, llcrnrx=-1000000.0,
                    urcrnry=1000000.0, urcrnrx=1000000.0,
                    resolution='l', area_thresh=1000)

theMap.drawcoastlines()
theMap.bluemarble()

plt.show()

···

--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Blue-marble-doesn-t-map-to-regional-ortho-projection-correctly-tp46944.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Indeed, that is obviously wrong. Could you file a bug report up on
basemap's github page? Basemap is now under new management, so I will see
what I can do to address the problem for you.

Cheers!
Ben Root

···

On Fri, Apr 8, 2016 at 7:49 PM, Don Morton <don.morton at borealscicomp.com> wrote:

Hello,

I could easily be doing something wrong. I posted this issue about two
years ago and never saw a response, so I'll try again. It's more of an
academic question - I can work around the issue. In the following code,
I'm
trying to plot a regional ortho projection over central Europe (I've
commented out the global projection), and then trying to use a bluemarble
background. As you can see from the attached image, global bluemarble
background is plotted, not a regional one. This also happens with
shadedrelief and etopo.

<http://matplotlib.1069221.n5.nabble.com/file/n46944/ortho.png&gt;

import mpl_toolkits.basemap as bm
import matplotlib.pyplot as plt
import numpy as np

"""
theMap = bm.Basemap(projection='ortho',
                    lat_0=48.2, lon_0=16.4)
"""

theMap = bm.Basemap(projection='ortho',
                    lat_0=48.2, lon_0=16.4,
                    llcrnry=-1000000.0, llcrnrx=-1000000.0,
                    urcrnry=1000000.0, urcrnrx=1000000.0,
                    resolution='l', area_thresh=1000)

theMap.drawcoastlines()
theMap.bluemarble()

plt.show()

--
View this message in context:
http://matplotlib.1069221.n5.nabble.com/Blue-marble-doesn-t-map-to-regional-ortho-projection-correctly-tp46944.html
Sent from the matplotlib - users mailing list archive at Nabble.com.
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users at python.org
Matplotlib-users Info Page

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/matplotlib-users/attachments/20160418/f7ee9eb2/attachment.html&gt;