python - Matplotlib Mollweide/Hammer 投影 : region of interest only

Question

我想知道是否有一种方法可以仅显示基于 basemap (matplotlib)中的 Mollweide/Hammer 投影的绘图中的感兴趣区域。

我正在尝试将绘图边界大致设置为太平洋板 block ，如下面的链接所示。然而set_xlim和set_ylim函数似乎没有任何效果。预先感谢您的任何指导。

http://geology.gsapubs.org/content/29/8/695/F1.large.jpg

Answer 1

根据文档，Hammer 和 Mollweide 投影都不允许这样做，因为它们打印出整个世界地图。下面是一些使用多圆锥投影的代码，但它是由直线包围的。这里的技巧是在创建时定义角点的经度和纬度。

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

my_map = Basemap(projection='poly', lat_0=0, lon_0=-160,
    resolution = 'h', area_thresh = 0.1,
    llcrnrlon=140, llcrnrlat=-60,
    urcrnrlon=-100, urcrnrlat=60)

plt.figure(figsize=(16,12))

my_map.drawcoastlines()
my_map.drawcountries()
my_map.fillcontinents(color='coral', lake_color='aqua')
my_map.drawmapboundary(fill_color='aqua')

my_map.drawmeridians(np.arange(0, 360, 20))
my_map.drawparallels(np.arange(-90, 90, 10))

plt.show()

结果:

请注意，这实际上显示的面积比您提供的图片中的面积要小。