Skip to content
Back to top
View this page

# Polyconic Projection

    from mpl_toolkits.basemap import Basemap
    import numpy as np
    import matplotlib.pyplot as plt
    # setup polyconic basemap
    # by specifying lat/lon corners and central point.
    # area_thresh=1000 means don't plot coastline features less
    # than 1000 km^2 in area.
    m = Basemap(llcrnrlon=-35.,llcrnrlat=-30,urcrnrlon=80.,urcrnrlat=50.,\
                resolution='l',area_thresh=1000.,projection='poly',\
                lat_0=0.,lon_0=20.)
    m.drawcoastlines()
    m.fillcontinents(color='coral',lake_color='aqua')
    # draw parallels and meridians.
    m.drawparallels(np.arange(-80.,81.,20.))
    m.drawmeridians(np.arange(-180.,181.,20.))
    m.drawmapboundary(fill_color='aqua')
    plt.title("Polyconic Projection")
    plt.show()

(`Source code`)